CN106895890B - 一种多声道超声波气体流量计声道权系数计算方法 - Google Patents

Abstract

本发明公开了一种多声道超声波气体流量计声道权系数计算方法。在流量计标定过程中，获取超声波气体流量计各个声道的顺、逆流渡越时间，通过温度和压力测量装置读取流量计工作温度和压力，结合标准表测量结果，使用粒子群算法优化参数的支持向量机算法，得出多声道超声波气体流量计声道权系数。支持向量机算法的使用能够有效降低流量计***误差。使用粒子群算法优化支持向量机算法的参数，能够有效降低人为设定参数带来的支持向量机算法计算结果的偏差。该方法可以统一应用于低速区和非低速区的测量，不必按照雷诺数大小划分流速区来分别进行流速校正和拟合。该方法可适应不同声道布置方式和位置，有效降低计量误差，实用性更强。

Description

一种多声道超声波气体流量计声道权系数计算方法

技术领域

本发明涉及一种多声道超声波气体流量计声道权系数计算方法，特别涉及一种基于粒子群算法优化参数的支持向量机算法的多声道超声波气体流量计声道权系数计算方法。

背景技术

超声波气体流量仪表是近年来广泛应用于天然气管道的流量测量装置。其具有的双向性、无阻流元件、重复性高、精度高和不带来管道压损等优点，使多声道超声波气体流量计取代传统的孔板流量计和涡轮流量计等，成为流量测量的主要装置。超声波气体流量计使用的原理是超声波在流体的影响下，顺、逆流传播时声波速度与流速叠加或者相抵，导致顺逆流传播时间不一致，并与流体流速相关，从而可通过对顺、逆流时间的测量反推流体流速。多声道超声波气体流量计是在单声道超声波气体仪表的基础上加入多个声道，声道所在的平面是平行于管道纵向的过直径横截面或者过弦横截面。通过测量管道不同位置的流体速度分布，按照相应的声道权系数方案进行计算后，得出流体流速。常用的多声道超声波气体流量计权系数方案有Gauss-Legendre方案、Tchebychev方案，Tailored方案和OWICS方案。与单声道超声波气体流量计相比，多声道超声波流量计的优点在于多声道的布局可以抵消一部分由于安装误差和电路时延等带来的***误差，对于速度分布的多位置测量能够更加准确地反应流动的真实状况，利于精确测量和流态监测。

在国际上，多声道超声波气体流量计是天然气贸易结算的标准器具。目前制造商主要为欧洲和美国的大型仪表制造企业，中国国内具有自主知识产权的精密超声波气体流量计制造公司几乎没有，仪表使用来源依赖进口，价格高、维修不便。多声道超声波气体流量计的研制是一个包含了软件和硬件的设计、结合和优化的过程。工作主要集中在在硬件方面，超声波换能器的制造、性能提高、优化电路以提高信号信噪比等，在软件方面，对于超声波在流体中传播方式进行深入探讨以改进流量计的测量原理、流量计校准方法的优化等。在多声道超声波气体流量计声道系数方案的研究上，Gauss-Legendre方案、Tchebychev方案，Tailored方案和OWICS方案。这些方案都建立在特定的声道位置布局基础之上，缺点在于在仪表实际制造过程中，无法保证声道位置完全按照设计位置排布，从而带来***误差，同时由于这些设计方案都基于充分发展的管道流体速度分布半经验模型，并不完全适合实际测量中气体流体的流动状态。从而导致多声道超声波气体流量计的测量准确度降低。

与本专利相关的文献和专利中，Yeh，T等人发表在2001年IEEE Instrumentation&Measurement Technology会议的论文“An intelligence ultrasonic flow meter forimproved flow measurement and flow calibration facility”中，提出了使用人工神经网络来计算多声道流量计的声道权系数，但并没有通过实验数据来证明这种方法的有效性。专利《一种采用高斯-雅可比多项式确定声道位置的设置新方法》(申请公布号：201610117702.7)使用高斯-雅可比多项式来确定声道的位置，要求声道位置必须依照计算方案固定安装，在实际仪表制造中难以保证，灵活性差。

发明内容

本发明目的在于克服上述现有研究和技术存在的问题和缺陷，提出一种多声道超声波气体流量计声道权系数计算方法。

本发明的目的通过以下的技术方案实现：一种多声道超声波气体流量计声道权系数计算方法，包括以下步骤：

1)将多声道超声波气体流量计安装在流量测试管道中，测试段上游需存在至少20倍管道直径长度的直管段或者相应整流设施，以保证超声波气体流量计测试管道内流体流动充分发展，测试管道包含的测量装置包括标准流量计、温度测量装置、压力测量装置；

2)根据超声波流量计测量范围和雷诺数计算公式，将超声波气体流量计测量范围划分为低速区和非低速区，最低流速到雷诺数4000对应的流速范围为低速区，雷诺数4000对应的流速至最高流速范围为非低速区；

3)在低速区选取均匀分布的若干流速点，进行流量计量测试；

4)在非低速区选取均匀分布的若干流速点，进行流量计量测试；

5)将流量计量测试中测得的数据作为支持向量机的输入；

6)选取支持向量机的惩罚因子C和核函数参数σ，给定惩罚因子C和核函数参数σ的上下限作为搜索域，随机生成若干组C和σ的初始值，将这若干组初始值记为(C,σ)，作为每个粒子在搜索域内的位置，将其作为支持向量机算法的初始参数设置，同时随机设定每个粒子的速度初始值；

7)使用支持向量机进行多维流量曲线拟合，得出流量v和温度T、压力p、各声道顺流渡越时间t_down、逆流渡越时间t_up的关系式v＝f(t_up,t_down,T,p)；

8)将关系式v＝f(t_up,t_down,T,p)的系数矩阵W作为该流量计声道权系数，并带入输入的各流速及其对应的温度T、压力p、各声道顺流渡越时间t_down、逆流渡越时间t_up数据，计算对应流速，将计算出的流速值与真实流速数据进行比对，并计算均方误差，作为当前粒子的适应度；

9)计算当前每个粒子的适应度，即采用了当前的惩罚因子C和核函数参数σ设定方案的支持向量机算法得出的流速计算值和真实值之间的均方误差值，并求出每个粒子的个体最优适应度和所有粒子的全局最优适应度；

10)根据粒子群算法中设定的粒子速度和位置进化规则，对每个粒子的位置和速度进行进化；

11)判断是否达到粒子群算法设定的最大代数，若达到设定最大代数，停止进行参数优化，选取粒子群算法的全局最优适应度，作为支持向量机的参数，否则返回步骤7)；

12)在流速范围内随机取一系列流速点进行流量计量测试，并用关系式v＝f(t_up,t_down,T,p)计算出流量值作为该流量计的测量流量值；

13)将测量流量值与标准流量计的计量值进行比对，计算相对误差、量程误差和精度等级，并重复步骤12)，检验流量计测量的可重复性；

14)判断量程误差和精度等级等指标是否达到技术指标要求，输出关系式v＝f(t_up,t_down,T,p)及其系数矩阵W，分别作为该台多声道超声波气体流量计的仪表特性曲线和声道权系数，否则返回步骤5)。

进一步地，该方法与传统声道布局权系数相比，优点在于不要求各个声道位置布局按照某种方案固定排布，权系数由流量仪表的流量测量数据作为输入。

进一步地，适应度通过均方误差来体现，其计算公式为：

其中：MSE是均方误差，N是输入支持向量机的测量数据组数，V_i是第i组测量数据的流量测量值，V_{M_i}是第i组测量数据由支持向量机得出的流速测量模型(v＝f(t_up,t_down,T,p))计算出的流速计算值。

与现有技术相比，本发明具有以下优点：

1、算法先进。支持向量机算法的使用解决了多声道超声波气体流量计测量原理中流量和顺、逆流渡越时间非线性关系和流量计生产制造中无法避免的硬件误差和电路时延不易测量等问题，能够有效降低流量计***误差。

2、精确度高。该方法使用粒子群优化算法来选择支持向量机的参数，这样的参数选择方式能够有效避免人为给支持向量机算法选择参数带来的支持向量机计算模型不精确，导致流量计测量误差大的结果，同时更易实现仪表的智能化和自动化。

3、实用性强。该方法可以统一应用于低速区和非低速区的测量，不必按照雷诺数大小划分流速区来分别进行流速校正和拟合。与传统的声道权系数计算方法相比，可适应不同声道布置方式和位置，有效降低计量误差，实用性更强。

4、扩展性好。该方法不要求声道布局必须在某些固定位置上，适用于按照Gauss-Legendre方案、Tchebychev方案，Tailored方案和OWICS方案等布局声道的多声道流量计，适用性广。同时，该研究方法和成果可以有效推广拓展至类似的问题求解中，解决相应问题。

附图说明

图1是本发明的多声道超声波气体流量计声道权系数计算方法的流程图。

图2是六声道超声波气体流量计结构示意图。

图3是由本发明的优化方法应用于实施例与传统计算方法应用于实施例的流量计算结果误差。

具体实施方式

以下结合附图对本发明的实施作如下详述：

实施例

本实施例对一台口径为100毫米的高斯-勒让德方案布局的六声道超声波气体流量计进行各个声道的权系数计算。该气体流量计测量对象为空气，工况为常温常压，声道布局结构如图2。使用一台安装在同一管路上、精度等级为0.5的罗茨流量计作为标准流量计。同时测量温度和压强数据。确定该六声道超声波气体流量计的实施步骤具体如下：

1)将多声道超声波气体流量计安装在流量测试管道中，测试段上游需存在至少20倍管道直径长度的直管段或者相应整流设施，以保证超声波气体流量计测试管道内流体流动充分发展，测试管道包含的测量装置包括标准流量计、温度测量装置、压力测量装置；

2)根据超声波流量计测量范围和雷诺数计算公式，将超声波气体流量计测量范围划分为低速区和非低速区，最低流速到雷诺数4000对应的流速范围为低速区，雷诺数4000对应的流速至最高流速范围为非低速区；

3)在低速区选取均匀分布的若干流速点，进行流量计量测试；

4)在非低速区选取均匀分布的若干流速点，进行流量计量测试；

5)将流量计量测试中测得的数据作为支持向量机的输入；

6)选取支持向量机的惩罚因子C和核函数参数σ，给定惩罚因子C和核函数参数σ的上下限作为搜索域，随机生成若干组C和σ的初始值，将这若干组初始值记为(C,σ)，作为每个粒子在搜索域内的位置，将其作为支持向量机算法的初始参数设置，同时随机设定每个粒子的速度初始值；

7)使用支持向量机进行多维流量曲线拟合，得出流量v和温度T、压力p、各声道顺流渡越时间t_down、逆流渡越时间t_up的关系式v＝f(t_up,t_down,T,p)；

8)将关系式v＝f(t_up,t_down,T,p)的系数矩阵W作为该流量计声道权系数，并带入输入的各流速及其对应的温度T、压力p、各声道顺流渡越时间t_down、逆流渡越时间t_up数据，计算对应流速，将计算出的流速值与真实流速数据进行比对，并计算均方误差，作为当前粒子的适应度；

9)计算当前每个粒子的适应度，即采用了当前的惩罚因子C和核函数参数σ设定方案的支持向量机算法得出的流速计算值和真实值之间的均方误差值，并求出每个粒子的个体最优适应度和所有粒子的全局最优适应度；

10)根据粒子群算法中设定的粒子速度和位置进化规则，对每个粒子的位置和速度进行进化；

11)判断是否达到粒子群算法设定的最大代数，若达到设定最大代数，停止进行参数优化，选取粒子群算法的全局最优适应度，作为支持向量机的参数，否则返回步骤7)；

12)在流速范围内随机取一系列流速点进行流量计量测试，并用关系式v＝f(t_up,t_down,T,p)计算出流量值作为该流量计的测量流量值；

13)将测量流量值与标准流量计的计量值进行比对，计算相对误差、量程误差和精度等级，并重复步骤12)，检验流量计测量的可重复性；

14)判断量程误差和精度等级等指标是否达到技术指标要求，输出关系式v＝f(t_up,t_down,T,p)及其系数矩阵W，分别作为该台多声道超声波气体流量计的仪表特性曲线和声道权系数，否则返回步骤5)。

支持向量机的适应度通过均方误差来体现，其计算公式为：

其中：MSE是均方误差，N是输入支持向量机的测量数据组数，V_i是第i组测量数据的流量测量值，V_{M_i}是第i组测量数据由支持向量机得出的流速测量模型计算出的流速计算值。

本发明基于粒子群算法优化参数的支持向量机算法的多声道超声波气体流量计声道权系数计算方法主要包括粒子群算法优化参数、计算当前支持向量机模型的适应度函数值、进一步更新粒子群算法中的粒子速度和位置、优化支持向量机参数以得到更小的适应度函数值等环节组成。图1为基于粒子群算法优化参数的支持向量机算法的多声道超声波气体流量计声道权系数计算方法具体流程。整个实施例是按照图1中所示的流程，进行多声道超声波气体流量计声道权系数计算。图2是实施例中使用本方法进行声道权系数计算的对象，一台按照高斯-勒让德方法布局的六声道超声波气体流量计。图3是由本发明的优化方法应用于实施例的计算结果的相对误差，其中圆形实心点为使用该方法计算出来的声道权系数进行计算得出流量的相对误差，空心三角形为使用传统的高斯-勒让德积分方法计算的流量相对误差。

在全量程范围内，使用该方法计算出来的声道权系数进行计算得出流量的相对误差在±1％以内，使用传统的高斯-勒让德积分方法计算的流量相对误差则在±3％以内。使用基于粒子群算法优化参数的支持向量机算法的多声道超声波气体流量计声道权系数计算方法计算结果显示，该方法使用支持向量机算法的解决了多声道超声波气体流量计测量原理中流量和顺、逆流渡越时间非线性关系和流量计生产制造中无法避免的硬件误差和电路时延不易测量等问题，能够有效降低流量计***误差。同时，该方法具有较高的精确度，使用粒子群优化算法来选择支持向量机的参数，这样的参数选择方式能够有效避免人为给支持向量机算法选择参数带来的支持向量机计算模型不精确，导致流量计测量误差大的结果，同时更易实现仪表的智能化和自动化。同时，方法可以统一应用于低速区和非低速区的测量，不必按照雷诺数大小划分流速区来分别进行流速校正和拟合。与传统的声道权系数计算方法相比，可适应不同声道布置方式和位置，有效降低计量误差，实用性更强。此外，该方法不要求声道布局必须在某些固定位置上，适用于按照Gauss-Legendre方案、Tchebychev方案，Tailored方案和OWICS方案等布局声道的多声道流量计，适用性广。同时，该研究方法和成果可以有效推广拓展至类似的问题求解中，解决相应问题。

Claims (2)

1.一种多声道超声波气体流量计声道权系数计算方法，其特征在于，包括以下步骤：

1)将多声道超声波气体流量计安装在流量测试管道中，测试段上游需存在至少20倍管道直径长度的直管段或者相应整流设施，以保证超声波气体流量计测试管道内流体流动充分发展，测试管道包含的测量装置包括标准流量计、温度测量装置、压力测量装置；

2)根据超声波流量计测量范围和雷诺数计算公式，将超声波气体流量计测量范围划分为低速区和非低速区，最低流速到雷诺数4000对应的流速范围为低速区，雷诺数4000对应的流速至最高流速范围为非低速区；

3)在低速区选取均匀分布的若干流速点，进行流量计量测试；

4)在非低速区选取均匀分布的若干流速点，进行流量计量测试；

5)将流量计量测试中测得的数据作为支持向量机的输入；

6)选取支持向量机的惩罚因子C和核函数参数σ，给定惩罚因子C和核函数参数σ的上下限作为搜索域，随机生成若干组C和σ的初始值，将这若干组初始值记为(C,σ)，作为每个粒子在搜索域内的位置，将其作为支持向量机算法的初始参数设置，同时随机设定每个粒子的速度初始值；

7)使用支持向量机进行多维流量曲线拟合，得出流量v和温度T、压力p、各声道顺流渡越时间t_down、逆流渡越时间t_up的关系式v＝f(t_up,t_down,T,p)；

8)将关系式v＝f(t_up,t_down,T,p)的系数矩阵W作为该流量计声道权系数，并带入输入的各流速及其对应的温度T、压力p、各声道顺流渡越时间t_down、逆流渡越时间t_up数据，计算对应流速，将计算出的流速值与真实流速数据进行比对，并计算均方误差，作为当前粒子的适应度；

9)计算当前每个粒子的适应度，即采用了当前的惩罚因子C和核函数参数σ设定方案的支持向量机算法得出的流速计算值和真实值之间的均方误差值，并求出每个粒子的个体最优适应度和所有粒子的全局最优适应度；

10)根据粒子群算法中设定的粒子速度和位置进化规则，对每个粒子的位置和速度进行进化；

11)判断是否达到粒子群算法设定的最大代数，若达到设定最大代数，停止进行参数优化，选取粒子群算法的全局最优适应度，作为支持向量机的参数，否则返回步骤7)；

12)在流速范围内随机取一系列流速点进行流量计量测试，并用关系式v＝f(t_up,t_down,T,p)计算出流量值作为该流量计的测量流量值；

13)将测量流量值与标准流量计的计量值进行比对，计算相对误差、量程误差和精度等级，并重复步骤12)，检验流量计测量的可重复性；

14)判断量程误差和精度等级等指标是否达到技术指标要求，输出关系式v＝f(t_up,t_down,T,p)及其系数矩阵W，分别作为该台多声道超声波气体流量计的仪表特性曲线和声道权系数，否则返回步骤5)。

2.根据权利要求1所述的一种多声道超声波气体流量计声道权系数计算方法，其特征在于，使用粒子群算法优化参数的支持向量机算法进行模型训练的过程中，使用均方误差最小来作为适应度函数，其计算公式为：