超声波水表流量系数处理方法
【专利摘要】本发明涉及超声波流量检测【技术领域】，特别涉及一种超声波水表流量系数处理方法，包括如下步骤：从内存中取出本次检测得到的时差数值，记为t；根据程序内置时差分界数组A[m]得出时差分区编号i；从内存中取出本次检测得到的水温数值，即为w；根据程序内置水温分界数组B[n]得出水温分区编号j；根据时差分区编号i、水温分区编号j从程序内置二维系数矩阵K[m,n]中取出四个元素ki,j、ki,j+1、ki+1,j、ki+1,j+1并进行处理得到拟合系数k0＝S1＋(S2－S1)×δB；超声波声速系数W＝c4×w4＋c3×w3＋c2×w2＋c1×w；流体体积Q＝t×k0×W。通过一系列的处理，使得超声波水表对流量的计算更为精确，避免了各种因素对处理结果的影响，同时简化处理过程，大幅降低水表制造成本。
【专利说明】超声波水表流量系数处理方法
【技术领域】
[0001]本发明涉及超声波流量检测【技术领域】，特别涉及一种超声波水表流量系数处理方法。
【背景技术】
[0002]超声波水表是采用超声波时差原理，采用工业级电子元器件制造而成的全电子水表，与机械式水表相比较具有精度高，可靠性好，量程比宽，使用寿命长，无任何活动部件，无需设置参数，任意角度安装等优点，随着社会的进步，超声波水表的使用越来越广。
[0003]超声波水表在测量时，利用超声波脉冲在流体中顺流传播时间和逆流传播时间的差值来测量流体流速，再根据管道截面面积计算出管道内流体的流量。在测量时，由于超声波在水中传播速度很快，无论是顺流传播时间还是逆流传播时间都很小，故对时间的测量精度要求很高。现有的测量都是使用一个高频计数时钟，超声波脉冲在流体顺流传播过程中，即从发射到接收的过程中计数器对高频计数时钟进行计数，其顺流传播时间t就是计数值L除以高频计数时钟的频率f ;逆流传播时间同理，也是通过计数器的计数值计算得出的。现有技术中，对流体的流量的处理，是通过超声波的传播时间、超声波的声速得出的，其并没有根据各方面因素对结果进行修正，使得处理结果存在误差，使得流体流量计算的不够准确。
【发明内容】
[0004]本发明的目的在于提供一种超声波水表流量系数处理方法，保证超声波水表流量计算的准确性。
[0005]为实现以上目的，本发明采用的技术方案为:一种超声波水表流量系数处理方法，包括如下步骤:(A)从内存中取出本次检测得到的时差数值，记为t ; (B)根据程序内置时差分界数组A[m]得出时差分区编号i ; (C)从内存中取出本次检测得到的水温数值，即为w ;
(D)根据程序内置水温分界数组B[n]得出水温分区编号j ; (E)根据时差分区编号1、水温分区编号j从程序内置二维系数矩阵K[m，η]中取出四个元素<」、1^」+1、1^1，」、1^1，」+1并进行处理得到拟合系数kQ = S1 + (S2 - S1) X δΒ ;(F)超声波声速系数W = C4Xw4 + C3Xw3 +C2Xw2 + C1 Xw ； (G)流体体积 Q = t Xk0Xff0
[0006]与现有技术相比，本发明存在以下技术效果:通过一系列的处理，使得超声波水表对流量的计算更为精确，避免了环境温度、流速等因素对处理结果造成影响；同时，这种处理方法简单，大幅 降低电路制造成本，使得超声波水表能够得到广泛应用。
【专利附图】
【附图说明】
[0007]图1是本发明的流程示意图。
【具体实施方式】[0008]下面结合图1，对本发明做进一步详细叙述。
[0009]一种超声波水表流量系数处理方法，包括如下步骤:(A)从内存中取出本次检测得到的时差数值，记为t;(B)根据程序内置时差分界数组A[m] = {&1，％，％，...，&π}得出时差分区编号，其中超声波水表的时差检测范围为[tmin，tmax]且ai < tmin < a2〈…< am_i< tmax < am,若t < tmin或t > tmax,超声波水表报错并停止流量检测,若ai≤t < ai+1,时差分区编号即为i ; (C)从内存中取出本次检测得到的水温数值，即为w ; (D)根据程序内置水温分界数组B[n] = Ibpb2^3，…，bn}得出水温分区编号，其中超声波水表的温度检测范围为[Wniin, WmaJ 且 Id1 < Wmin < b2 < …< V1 < wmax < bn，若 w < Wmin 或 w > Wmax,超声波水表报错并停止流量检测，若bj ^w< bA，水温分区编号即为j。以上步骤也可以按照步骤C、D、A、B的顺序进行。(E)根据时差分区编号1、水温分区编号j从程序内置二维系数矩阵K[m, n] = {{k1;1, k12,…，kl n}，{k2;1, k2，2,…，k2，J，…，{km；1, km；2,..., km；n}}中取出四个元素Κ」+1、1--+1，」、1^+1，」+1,按作如下处理:δ Α = (t — Ci1)/ (ai+1 — Ci1) , δ B = (w — bj)/ (bJ+1 —bj)，S1 = (k1；1 + k1；2 H l.k1； j) + (k1； j + k2，j H l.ki；J) + ki+1；JX δ A, S2 = (k1；1 +k1>2 H l.k1；J+1) + (k1； J+1 + k2，j+1 H l.ki+1；J+1) + ki+1；J+1X δΑ，拟合系数 kQ = S1 + (S2 —S1) X δ B ； (F)超声波声速系数 W = C4Xw4 + C3Xw3 + C2Xw2 + C1Xw ； (G)流体体积 Q =tXk0Xffo通过这一系列的处理，得出的流量更为准确。
[0010]作为本发明的优选方案，所述的二维系数矩阵K[m，η]按如下步骤确定:(a)依据超声波水表管端的流态流体力学数据建模确定二维曲面r = r (t, w) ； (b)将二维曲面r在t=a1； a2, a3,..., am 以及 w = Id1, b2, b3,...，bn共 mXn 个坐标点上离散化得到 = Hai, bj)；
(c)令ki;j = ιj将ki;j存储至二维系数矩阵K[m, η]。
[0011]更进一步地,所述的步骤c中，选取参数α和β，令I^j= (ri；J — β ) / a使得k,j为一个较小的数值便于存储和计算，将ku存储至二维系数矩阵K[m，η];所述的步骤E中，根据拟合系数h计算流量系数r(ι = a Xk0+ β，其中α，β即为前面所选取的参数；所述的步骤G中，流体体积Q = t X r(ι X W。通过这样处理，减少了处理过程中数据的计算和
存储量。
[0012]本方法可在保证数据处理精度要求的前提下极大简化流量数据计算过程，提高处理效率，降低处理器及存储器的硬件资源要求，使用IMHz至4MHz主频的处理器，存储空间在256Byte至1024Byte之间，即可正常完成流量计算处理过程。比使用本方案前与本方案后的硬件资源要求如下:
【权利要求】
1.一种超声波水表流量系数处理方法，包括如下步骤: (A)从内存中取出本次检测得到的时差数值，记为t； (B)根据程序内置时差分界数组A[m]= {a^w.’aj得出时差分区编号，其中超声波水表的时差检测范围为[tmin, tmax]且S1 < tmin < a2〈…< < tmax < am,若t < tmin或t > tmax，超声波水表报错并停止流量检测，若ai ( t < ai+1，时差分区编号记为i ; (C)从内存中取出本次检测得到的水温数值，即为w； (D)根据程序内置水温分界数组B[n]= Kb2, bf，bn}得出水温分区编号，其中超声波水表的温度检测范围为[wmin, WmaJ且Id1 < Wmin < b2〈…< Iv1 < Wmax < bm,若w < Wmin或w > Wmax，超声波水表报错并停止流量检测，若bj≤w < bj+1，水温分区编号记为j ； (E)根据时差分区编号1、水温分区编号j从程序内置二维系数矩阵K[m，n]={{k1；1, k1；2,..., k1；n}, {k2;1, k2，2,…，k2，J,…，{km；1, km；2,..., km；n}}中取出四个元素 ki;』、ki;J+1、ki+i，j、ki+1，j+1,按作如下处理:
δ a = (t — aj/ (ai+1 — at), δ B = (w — b」)/ (bJ+1 — b」),
51= (k1；1 + k1；2 H h k1；J) + (k1； j + k2，j H l.ki;j) + ki+1； j X δ A,
52= (ki,i 十 ki’2 H ^ k1；J+1) + (k1；J+1 + k2’j+1 H l.ki+1；J+1) + ki+1’j+1X δ A,
拟合系数 k0 = S1 + (S2 — S1) X δ B ; (F)超声波声速系数W = C4Xw4 + C3Xw3 + C2Xw2 + C1 Xw ； (G)流体体积Q= tXk0Xff0
2.如权利要求1所述的超声波水表流量系数处理方法，其特征在于:所述的二维系数矩阵K[m, η]按如下步骤确定: (a)依据超声波水表管端的流态流体力学数据建模确定二维曲面r= r (t, w)； (b)将二维曲面r在t= a” a2, a3,…，am以及w = b” b2, b3,...，bn共mXn个坐标点上离散化得到r” j = r (au b」)； (c)令ki;j= A，将ki;j存储至二维系数矩阵K[m, n]。
3.如权利要求2所述的超声波水表流量系数处理方法，其特征在于:所述的步骤c中，选取参数α和β，令ku= (ri；J - β)/α使得为一个较小的数值便于存储和计算，将ky存储至二维系数矩阵K[m，n]； 所述的步骤E中，根据拟合系数1?计算流量系数r0(ι= a Xk0+ β，其中α，β即为前面所选取的参数； 所述的步骤G中,流体体积Q = tXi^XW。
【文档编号】G01F15/02GK103913203SQ201410148902
【公开日】2014年7月9日 申请日期:2014年4月14日 优先权日:2014年4月14日
【发明者】姜跃炜, 薛明星 申请人:姜跃炜