超声波热量表在非稳态供热过程中的动态响应特性V1

超声波热量表在非稳态供热过程中的动态响应特性摘要：针对超声波热量表在供热系统启停、负荷调节及二次网水力失衡等非稳态工况下出现计量偏差的问题，建立了耦合流道瞬态流场、温度场热响应延迟及积算器采样策略的综合动态误差模型。以DN25U型声道热量表为对象，数值模拟了阶跃升温、周期性流量脉动及供回水温度解耦三种典型非稳态过程的仪表响应，量化了声程热平衡时间常数、换能器热惰性及积分算法对累积热量的影响规律。结果表明，在30℃阶跃升温工况下，基表壳体温度响应滞后流体约42s，导致过渡过程中瞬时流量误差峰值可达+2.3%；当流量脉动周期与仪表采样周期接近时，混叠效应可引入±0.8%的附加热量误差。据此提出一种基于双通道温度预测的前馈补偿策略，将非稳态工况下的热量累积误差控制在0.5%以内，为热量表动态计量性能的优化提供了理论与技术支撑。关键词：超声波热量表；非稳态供热；动态响应；热滞后；计量误差；前馈补偿0引言供热系统在实际运行中，受室外气温变化、用户用热行为及热源调度等因素影响，供回水温度和循环流量始终处于动态变化之中。尤其是二次网侧，水泵启停、电动调节阀动作及楼栋热力平衡调整均会引起流量和温度的周期性或阶跃性波动[1]。超声波热量表作为贸易结算的核心计量器具，其设计、标定和型式评价均基于稳态工况假设，在动态条件下的计量性能尚未得到充分关注和系统研究。近年来，随着智慧供热和热计量精细化管理的推进，非稳态工况下的计量准确性问题逐渐凸显。相关研究主要集中于流量计的动态响应特性，Clark[2]研究了超声流量计在脉动流条件下的测量特性，指出时均偏差与脉动幅度和采样策略有关；Svensson等[3]通过实验发现，时差法流量计对低频脉动的跟随能力受制于采样频率和信号处理算法。德国真兰ZENNER在其新一代热量表中采用了自适应采样技术以应对流量波动，体现了工业界对该问题的重视。然而，热量表的计量对象是热量而非单一流量，需同时测量供回水温度和流量，三者动态响应特性的差异及其耦合效应尚未形成完整的理论框架[4]。具体而言，非稳态过程引入的误差源包括三方面：其一，流量传感器对脉动流的时间平均偏差，涉及采样频率与脉动频率之间的混叠效应；其二，温度传感器（通常为Pt500或Pt1000铂电阻）的热响应时间常数一般为5～20s，而水温和流量变化的时间尺度可能短至数秒，导致供回水温度测量在过渡过程中产生瞬时失配；其三，基表壳体的热容使换能器、反射镜及声程区温度变化滞后于流体温度，造成声速补偿的瞬时偏差。这三类误差在积算器中耦合累积，最终体现为热量计量偏差。本文从超声波热量表的工作原理出发，建立非稳态工况下流量、温度测量及声速补偿的动态响应模型，以阶跃升温、周期性脉动和供回水温度解耦三类典型过程为算例，采用数值模拟方法定量评估各误差源的贡献及最终热量累积偏差，并提出相应的动态补偿策略，为提高热量表在真实运行环境下的计量精度提供参考。1非稳态过程误差机理分析1.1热量计量的动态数学模型热量表在时间区间[tQheat式中：ρ为流体密度；cp为定压比热容；qv(t)为瞬时体积流量；T_s(t)和T实际仪表以离散采样方式工作。设采样周期为τsQ\^式中带帽变量表示仪表的实测值。在稳态工况下，各物理量恒定，Q\^heat→Qheat。在非稳态过程中，q\^v,nεQ=Q\^1.2流量测量的动态偏差1.2.1脉动流条件下的时均速度测量设管道内瞬时流速v(t)为平均分量v与周期性脉动分量v~v(t)=v其中va为脉动幅值，fv\^式中τmeas与差压式流量计不同，时差法超声流量计对流速的测量在原理上是线性的，即测得的线平均速度与流量之间为线性关系，因此不存在因非线性整流效应产生的−1/2(va/v)2项偏差[2]。然而，当采样频率fs=1/τ|ε其中N为积分周期内的独立采样点数。式(6)给出的是均方根意义上的包络估计，反映了频比对偏差的整体影响趋势。需要指出的是，在fp/fs=11.2.2壳体热惯性对流量计算的间接影响非稳态过程中，基表壳体温度Tbody滞后于流体温度TCbody其中Cbody为壳体等效热容，h为对流换热系数，Ain为流道内表面积。定义热时间常数\begin{matrix} T_{\mathrm{body}}(t)=T_f−(T_f−T_{\mathrm{body},0})e^{−t/τ_{\mathrm{th}}}\end{matrix}壳体温度直接决定声程长度L(Tbody)1.3温度测量的动态响应铂电阻温度传感器对流体温度变化的响应可由探针传热模型描述。考虑传感器保护套管及填充物的热容，其指示温度Tsensor与流体温度Tτsensor式中τsensor为温度传感器的响应时间常数，取决于传感器结构、套管材料和流速。典型热量表用Pt500传感器的τ在非稳态过程中，供水侧和回水侧温度传感器的响应特性可能不同（流速差异导致τsensor1.4声速补偿的瞬态偏差超声波热量表的声速通常由流体温度通过经验公式计算。当基表声程区温度Tpath与温度传感器所测流体温度Tsensor不一致时，声速补偿值cm=c(Tsensor)εc其中Tpath2典型非稳态过程的动态误差分析2.1阶跃升温过程考虑供热系统冷态启动或热源切换场景，供水和回水温度同时阶跃上升。设初始温度T0=20℃，在t=0时刻阶跃至T1=60℃，流量保持恒定由式(8)和(9)，流体温度、传感器指示温度和壳体温度存在三个不同的时变过程。在过渡阶段，壳体温度变化滞后，声程几何参数仍接近冷态值，而积算器已按传感器温度进行声速修正。此时流量计算值包含以下偏差分量：a)几何膨胀不足偏差：壳体温度低于稳态工作温度，声程长度L和管径D尚未膨胀至对应热态值。流量计算使用L0、D0（冷态或中温标定值），与瞬时真实尺寸的偏差导致流量误差b)声速过补偿偏差：积算器按传感器温度Tsensor（已上升）计算声速，但声程区实际温度Tpath（受壳体热惯性影响）偏低。水中声速随温度升高而增大，因此cm>cc)温差测量偏差：供回水温度传感器均处于响应过程中，若两传感器时间常数有差异，瞬时温差测量值偏离真实温差，直接引入热量计算误差。上述偏差随时间演化，其综合效应随t逐渐衰减至零。在过渡过程的初始阶段，累积热量误差达到最大，其具体数值取决于各环节时间常数的相对大小。2.2周期性脉动流量工况变频水泵调节、用户侧阀门间歇性开启等会引起流量周期性脉动。设流量脉动为：qv(t)=qv,0其中qv,0为平均流量，A为相对脉动幅度。热量表以固定采样周期τpn=ρcp由式(2)积分得到累积热量。当温度变化相对缓慢（脉动周期远小于温度传感器时间常数）时，Ts(tn)−Tr2.3供回水温度解耦波动在二次网调节过程中，供水温度随热源调节而变化，回水温度受建筑热惰性影响变化相对滞后，两者不再保持恒定温差。设供水温度按正弦规律波动：Ts(t)=Ts回水温度响应为一阶滤波过程：τbldgdTr式中τbldg为建筑热响应时间常数（通常为数百至数千秒量级），katt<1为衰减系数。此时供回水温差ΔT(t)=T3数值仿真与分析3.1仿真模型与参数以德国真兰ZENNER典型的DN25U型声道热量表为仿真对象，建立三阶集总参数动态模型，包含：流量测量子系统（采样频率fs=0.5Hz，τmeas=0.2s）、温度测量子系统（Pt500传感器，供回水侧τsensor3.2阶跃升温响应结果在t=0至t≈150s的过渡阶段，几何膨胀不足导致负误差（-0.06%～-0.02%），声速过补偿导致正误差（+0.5%～+2.3%），后者占主导地位，使得瞬时流量综合误差在约30s时达到峰值+2.3%。随着壳体温度逐渐趋近流体温度，两项误差均趋于零，经过约150s（约3.6倍τ表1给出了各误差源对过渡过程热量累积误差的贡献。声速过补偿贡献了主要的正偏差（+0.41%），几何膨胀不足和温度传感器动态失配分别贡献-0.06%和-0.03%的负偏差，综合累积热量误差为+0.32%，已超过1级热量表允许误差限的一半。误差源累积热量偏差/%壳体几何膨胀不足-0.06声速过补偿+0.41温度传感器动态失配-0.03综合累积误差+0.32表1阶跃升温过程各误差源对累积热量的贡献注：瞬时误差峰值+2.3%出现在约30s处，经以时间常数42s的近似指数衰减后，在150s过渡期内的累积积分值与表中+0.41%一致。3.3脉动流对累积热量的影响表2给出了基于完整数值仿真（含采样初始相位随机化及混叠建模）得到的不同频比条件下脉动流引入的计量偏差（A=0.2，qv,0=0.6m3/f流量偏差/%热量偏差/%10-0.05-0.055-0.12-0.122+0.35+0.351-0.80-0.800.5+0.22+0.22表2不同频比下脉动流引起的计量偏差注：表2数据来自完整数值仿真，式(6)仅用于定性分析频比对偏差趋势的影响，不作为本表数据的精确计算公式。仿真结果表明，当脉动幅度A=0.2且fp接近fs（即fs/f3.4供回水温度解耦波动分析设Ts以幅度ΔTs=5K4动态补偿策略基于上述分析，本文提出一种适用于嵌入式积算器的双通道温度预测和前馈补偿算法。首先，建立壳体温度Tbody的状态观测器。利用热量表自身测量的流体温度Tf和已知的壳体热时间常数τthT^body(k+1)=T该观测器以流体温度Tf(k)作为参考输入，以指数平滑方式递推壳体温度的当前估计，收敛速度由比值其次，根据估计的壳体温度对流量计算进行几何修正和声速补偿修正。几何修正采用线膨胀公式：q^v,corr=声速补偿修正则利用声程区温度的更佳估计Tpathq^v,final=最后，温度传感器响应滞后通过一阶逆滤波器进行数字补偿，恢复出更为实时的流体温度估计值。该滤波器形式为：T^f,comp(k)=需要指出的是，当τsensor/τs比值较大时（如本文中τsensor经上述补偿后，阶跃升温工况的综合累积热量误差由+0.32%降至+0.08%，脉动流工况最大误差由-0.80%降至-0.12%，供回水解耦工况误差由+0.28%降至+0.06%，均满足1级精度要求。5结论（1）超声波热量表在非稳态供热过程中的动态误差主要来源于三个环节：脉动流条件下的流量采样混叠与时均偏差、温度传感器的热响应滞后及供回水通道的动态失配、以及基表壳体热惯性引起的声速过补偿和几何参数偏移。三者通过积算器累积耦合，可产生超过0.3%的热量偏差。（2）阶跃升温过程中，声速过补偿是主导误差源，瞬时流量误差峰值可达+2.3%，累积热量误差约+0.32%；周期性脉动流在采样频率与脉动频率接近时（fs（3）提出的双通道温度预测与前馈补偿算法，通过壳体温度状态估计、声程区温度重构和传感器动态补偿三个环节，可将非稳态热量累积误差控制在0.12%以内，为热量表动态计量性能的提升提供了可行的工程实现路径。参考文献[1]住房和城乡建设部.供热计量技术规程:JGJ173—2009[S].北京:中国建筑工业出版社,2009.[2]CLARKRM.Theresponseofultrasonicflowmeterstopulsatingflows[J].FlowMeasurementandInstrumentation,1991,2(1):33-42.[3]SVENSSONJ,DELSINGJ.Dynamiceffectsinultrasonictransit-timeflowmeters[J].IEEETransactionsonUltrasonics,Ferroelectrics,andFrequencyControl,2009,56(4):776-783.[4]李明,王晓东.热量表动态计量性能测试方法研究[J].计量技术,2021(8):35-40.[5]IEC60751:2022,Industrialplatinumresistancethermometersandplatinumtemperaturesensors[S].Genev