井下超声流量计的误差分析

林龙贵

【摘 要】井下超声流量计是贵州航天控制技术有限公司在国内首次研制成功并开辟了广阔市场的产品，它是油田实施二次、三次采油工艺的主要测试装备，也是国家重大攻关项目之—。井下超声流量计的开发，曾先后被列入贵州省科技攻关计划，国家级重点新产品试制鉴定计划和国家级火炬计划等。

【Abstract】The underground ultrasonic flowmeter is the first product developed by Guizhou Aerospace Control Technology Co. Ltd. in China and has opened up a vast market. It is the main test equipment for the secondary and tertiary oil recovery process in the oilfield， and it is also a major national key project. The development of underground ultrasonic flowmeters has been included in the Guizhou Province scientific and technological research plan， the national key new product trial production and certification program and the national-level torch plan.

【关键词】超声流量计;相位法;误差

【Keywords】ultrasonic flowmeter; phase method; error

【中图分类号】TE931                                            【文獻标志码】A                                【文章编号】1673-1069（2019）07-0190-03

1 引言

在贵州航天控制技术有限公司的井下超声流量计开发成功之前，油田使用的流量计主要由浮子式、涡轮式、电磁式等几种原理制造的产品，由于受井内原油、泥沙和液体黏度的影响，这些流量计的测试精度和可靠性都较差。井下超声流量计从原理上解决了粘度法测定聚合物流量的测试问题。

2 随机误差分析

随机误差是产品在相同条件下，多次重复测量同一物理量时，误差值呈无规律变化的一种误差，随机误差不能用实验方法消除，也不能修正，只能依靠统计规律了解其分布特性[1]。

2.1 流量微分方程

2.2 随机误差计算

对测试的数据进行统计后，得出了在同一条件下的流量测试的原始数据，其传播时间“T”随机变化dT=1～2μs，相位“Δ？覫”随机变化dΔ？覫=10～20个A/D采集位，由于“T”和“Δ？覫”的随机变化，导致了在同一条件下的流量Q的随机变化。

在以下的误差计算中，设T=134μs（多数仪器测试的时间“T”为130μs-140μs=134μs），dT=1μs，dΔ？覫=10（针对多数仪器的dT、dΔ？覫的变化范围），带入式（2）得：

3 “相位法”误差分析

3.1 原理误差分析

3.2 原理误差估算

在误差估算中，取时间T2=1000（这是相位法计算中给定的值），声速c从1480～1555m/s，即水温从20～70℃时声速的变化（没有考虑压力对声速c的影响）。则：

在式（8）中，Δt取2μs（即2μs所产生的相位的A/Dv采集值为4095），c取1500m/s，式（5）中的“-”号表示当声速“c”变大时，相位的A/D采集值变小，即实测流量变小。

4 “相位-时间法”误差分析

根据式（1），考虑到Δ？覫和“T”的综合影响，求全微分方程为：

此式与式（2）的区别在于没有绝对值符号，这给两部分误差的相互抵消带来了可能。

当声速“c”从1480m/s变到1555m/s时，dΔ？覫A/D（max）的表达式和式（5）一样，即：dΔ？覫A/D=-35.2v=-9.68×10-2（A/Dv-A/Dv=0）;

dΔ？覫=T2-T1=L/c2-L/c1=0.2/1555-0.2/1480=-6.518（μs）

在式（8）中，取T=134μs，结合上两式得：

dQ=（1.1138×10-3a1+2.48×10-6a2Δ？覫+4.146×10-9a3Δ？覫2）×dΔ？覫A/D+（1.66×10-5a1Δ？覫+3.703×10-8a2Δ？覫2+6.19×10-11

a3Δ？覫3）×6.518 （9）

式（9）即为“相位-时间法”在声速“c”从1480m/s变到1555m/s时的误差估算公式。

5 误差的验证

5.1 随机误差估算的验证

统计出04A057号仪器的相位Δ？覫max=800，带入式（3）得：

dQ100=1.1138×10-3a1+1.985×10-2a2+2.2653×10-2a3+

1.33×10-2a1+2.37×10-2a2+3.168×10-2a3（10）

对04A057号仪器的标定（标定精度为1.74%）得：

a1=126.58903，a2=8.71896，a3=-5.56843代入式（10）得：dQ100=0.166+1.71=1.876（m3/d），即随机相对误差=1.88%。

在随机误差估算中，流量误差的产生由两部分组成：即由“dΔ？覫”产生的相位测试误差和由“dT”产生的时间测试误差之和。从A-100m3/d、A-200m3/d、A-300m3/d的标定数据看：相位测试误差的绝对值比较小（均在1m3/d以下）、而时间测试误差的绝对值随量程的增大而增大，相对误差之和随量程的增大而减小，但最大不大于3%（F·S）。

根据误差分配原则+，在温度、压力变化不大的情况下，采用相位法（x法）克服了“dT”的影响，即排除了时间测试误差，所以更容易标检合格。

以上随机误差的大小依赖于仪器的测试精度和换能器的稳定度。

根据误差理论，该随机误差遵循正态分布，并且可针对具体的流量台阶的多次测试数据计算出特定流量条件下的标准偏差“σ”及精密度指数“h”。由于随机误差的概率密度分布函数计算需要依靠大量的统计数据才能保证其准确性，统计和计算的工作量大，故忽略此环节。

5.2 相位法误差估算的验证

以下验证中的方程形式為：Q=a3x3+a2x2+a1x+a0，其中x=

Δ？覫/1000，Δ？覫为相位的“A/D”采集值（即Δ？覫=A/Dv）。

对03A052号仪器进行相位标定得实际拟合方程：

Q=4.58049x3-13.47665x2+151.74806x-12.20589，

取A/Dv=794（实测最大值），A/Dv=0=80（实测零位），则由式（5）得：dΔ？覫A/D（max）=-9.68×10-2×714=69.12，将dΔ？覫A/D（max）代入式（7）得：

从上看出，当声速“c”从1480m/s变到1555m/s（即水温从20℃变到70℃）时，相同的实际流量，测试时可差10%左右，反映到实际生产过程中，超差的仪器大部分出现在水温变化较大的时候。

6 结论

依据对流量计存在的随机误差、“相位法”误差和“相位-时间法”误差的分析结果，结合流量计目前的技术水平，可以认为：

①随机误差是暂时无法消除的（误差一般在2%～4%左右），它主要由换能器决定，这也是造成流量标定精度较高、而检定精度较低;或流量标定精度较低、而检定精度较高的主要原因;

②“相位法”误差属于原理性误差。当工作环境变化较大时会产生较大的误差，如果考虑压力的影响，则误差将会更大;

③“相位-时间法”误差也属于原理性误差。主要由电路产生的相位零位引起，因为相位零位是流速=0m/s时的相位，由此看见，从原理上就把相位分成了两部分：固定相位，和流速相位，相位本来随声速“”变化，样这一来，为常数的“固定相位”和“流速相位”一起变成了自变量。目前，生产的流量计的设计零相位均在70～150个A/D值左右，所以由于零相位带来的误差均在2m3/d左右。

消除“相位-时间法”误差的办法是在拟合方程时减去零相位，这一点在以后的计算中可以进一步完善。

综上所述，可以得出两条结论：

第一，在今后的流量标检中，只能采用“相位-时间法”;第二，随机误差遵循正态分布，且主要受测试时间“”的影响，所以，在仪器标检超差（超差幅度为：4%以内）的情况下，可以将上下换能器对换后重新标检，但超差大于4%时要另找其他的原因。

【参考文献】

【1】杨正一，贾振安.误差理论与测量不确定度[M].北京：石油工业出版社，2000.