两相流含水率自适应串联估计方法技术

本发明专利技术涉及一种两相流含水率自适应串联估计方法，利用上、下游不同传感器的测量数据进行估计，包括如下步骤：对上游传感器和下游传感器测量数据进行时间配准；状态估计算初始化；对上游传感器含率测量值进行线性卡尔曼状态估计和更新，获得对状态向量的预测和协方差矩阵的预测；3)利用下游传感器测量数据，计算两相流过程状态转换矩阵：对下游含率测量值进行基于无迹卡尔曼滤波非线性状态估计。

Adaptive cascade estimation of water flow in two-phase flow

The invention relates to a method for the estimation of series rate adaptive water two-phase flow were estimated by the upper and lower reaches of different sensor measurement data, which comprises the following steps: time of registration data upstream and downstream sensor sensor; state estimation is initialized; on upstream sensor holdup measurement value of linear state estimation and prediction of Calman update the state vector and covariance matrix forecasts; 3) using the downstream sensor measurement data, calculation of two-phase flow process state transition matrix of the downstream holdup measurement value of unscented filter for nonlinear state estimation based on Calman.

【技术实现步骤摘要】
两相流含水率自适应串联估计方法
本专利技术属于流体测量
，利用自适应无迹卡尔曼算法提升相含率传感器对两相流含水率估计准确度。技术背景两相流广泛存在于日常生活与工业生产过程中，如食品加工，生物工程，化工产业，冶金工业和石油产业等行业。与一般单相流相比，两相流的流动状态更为复杂，建模更加困难，因此两相流的在线过程参数检测一直是学术界和工业界关注的难点和重点。特别在石油工业中，油水两相流含水率对油井监测、生产预测和开采计划制定具有重要意义。目前两相流含水率的测量手段包括快关阀测量方法、射线方法、微波方法、Wire-mesh方法、超声方法和电阻抗方法等。其中，电学方法具有测量速度快，敏感性强，传感器结构简单，易于安装，成本低廉等优点。两相流是一种随机的动态流动过程，流动状态复杂多变，传感器的测量信号有随机误差，且只能反映系统的外部特性。为获得多相流含水率的准确信息，需要根据实测数据对流动状态进行估计，以状态估计算法提升相含率传感器的测量精度。状态估计对认知和分析测量系统具有重要的应用价值，常用的估计理论包括最小二乘估计、线性最小方差估计和卡尔曼估计等。不同的估计方法性能指标不同，使用的信息不同，其适用范围也不同。最小二乘估计通常应用于随机向量和常值向量的估计，算法结构简单但精度有限；线性最小方差估计以均方误差最小作为优化参数，适用于平稳过程，其测量信息越丰富，估计精度越高，计算量越大，而且会随着估计过程的累积，造成计算负担较重，不适用于复杂的长时间的估计过程。卡尔曼估计采用递推的形式，通过一个状态模型和一个观测模型的迭代，利用上一时刻的估计值和这一时刻的观测值，对这一刻的状态进行估计。而基于不同结构的卡尔曼算法，卡尔曼估计又衍生出基本的线性卡尔曼估计，扩展卡尔曼估计，基于无迹卡尔曼滤波的非线性卡尔曼估计等估计方式，分别适用于不同对象。
技术实现思路
本专利技术提供一种两相流含水率自适应估计方法，该方法协同利用管道上、下游相含率传感器测量数据，上、下游相含率传感器可为同类传感器或数据采样特征匹配的异类传感器。对上游传感器测量数据进行线性卡尔曼滤波，计算其实时的状态转移矩阵，将该矩阵作为下游传感器状态转移矩阵的先验信息，代入下游传感器测试数据的无迹卡尔曼估计迭代过程中，实现含水率估计的自适应状态转换，同时调整无迹卡尔曼算法中Sigma点集结构，提高含水率的估计精度。本专利技术的技术方法如下：一种两相流含水率自适应串联估计方法，利用上、下游不同传感器的测量数据进行估计，包括如下步骤：1)对上游传感器和下游传感器测量数据进行时间配准：获取上游传感器和下游传感器所获相含率时间序列z1(j)和z2(j)，j＝1,2,3,…N为测量时间内的数据点，下标1和2分别代表上游传感器和下游传感器，计算渡越时间，利用渡越时间结合数据采样频率，对上游传感器和下游传感器测量数据进行时间配准；2)状态估计算初始化：其中，是上游状态向量的估计初值，是上游协方差的迭代初值，z1,k是上游传感器在k时刻的测量值，且是下游状态向量的估计初值；3)对上游传感器含率测量值进行线性卡尔曼状态估计和更新，获得对状态向量的预测和协方差矩阵的预测；在估计过程中量测矩阵H取值均为1，状态向量x1代表含水率，k时刻的系统模型为：x1,k+1＝F1,kx1,k+v1,k(4)z1,k＝H1,kx1,k+w1,k(5)其中，F1,k是k时刻的状态转移矩阵，x1,k+1和x1,k分别是上游在k+1时刻和k时刻的状态向量，H1,k是k时刻的量测矩阵，z1,k是量测向量，v1,k和w1,k分别表示系统在k时刻的过程噪声和量测噪声；则状态向量的预测的预测和协方差矩阵的预测为：根据上游量测信息z1,k修正估计值：其中，K1,k是上游传感器系统在k时刻的卡尔曼增益：其中，H1,k是k时刻上游传感器量测矩阵，是上游传感器在k时刻对状态的预测值，是上游传感器k时刻对状态x1,k的最小均方误差估计，是k时刻估计的协方差，R1,k为上游滤波过程在k时刻系统量测噪声的协方差，可取值为0.01至0.2；Q1,k-1为上游滤波过程在k-1时刻系统过程噪声的协方差，可取值为0.01至0.09，I为单位矩阵，其维数等同于状态向量的维数；3)利用下游传感器测量数据，计算两相流过程状态转换矩阵：将传感器上游测量单元确定的状态转移过程加以记录，用于时间配准后下游含水率状态估计中：其中，F2,k为下游测量单元的状态转移矩阵，k＝n0,n0+1,…,N；4)对下游含率测量值进行基于无迹卡尔曼滤波非线性状态估计：构造Sigma点集，其基本构造方式如式(12)所示：其中，L是状态向量的维度，是下游估计过程k-1时刻的状态估计，p2,k-1是下游估计过程k-1时刻的全局协方差，γ与λ可由下式计算：其中，λ为规模因子，α是决定Sigma点集中点的分散程度扩展因子，取值为0到1之间；κ是缩放系数，设为0，计算状态预测权系数Wi(m)和协方差预测权系数Wi(c)：其中，β先验分布因子，通常设置为2，i为矩阵在Sigma点集中对应的元素序号，下游状态估计过程的状态预测方式如下：估计过程中的协方差矩阵分为全局协方差矩阵和互协方差矩阵，更新方式分别如下：其中，p2,kk-1是对全局协方差的预测，Q2,k-1是下游过程噪声的协方差，H2,k是下游估计过程在k时刻的量测矩阵；其中，是下游传感器测量值在k时刻的协方差，是k时刻互协方差矩阵，R2,k为下游估计过程中的系统量测噪声的协方差，取值为0.01至0.09；K2,k是k时刻下游估计过程的卡尔曼增益，按照式(19)方式进行更新；基于无迹卡尔曼的非线性状态估计方式为：其中，H2,k是下游估计过程在k时刻的量测矩阵，是对含水率的最终估计结果，其均值即为测量时间段内平均相含率的最优估计。本专利技术的实质性特点是：针对多相流流动过程，采用上、下游相含率传感器测量所得相含率进行估计。利用测量对象的上、下游含率状态并行测量的特点，将估计过程分为一次线性滤波和一次非线性滤波，对上游传感器测量数据进行线性卡尔曼滤波，计算其实时的状态转移矩阵，将该矩阵作为下游传感器状态转移矩阵的先验信息，对下游传感器测试数据的进行基于无迹卡尔曼滤波的非线性状态估计，随估计过程改变Sigma点集而进行自适应卡尔曼估计。本专利技术的有益效果及优点如下：1、充分挖掘两相流上、下游相含率传感器测量数据之间的内在联系，实现同类传感器测试信息的有机融合；2、增加了对于两相流对象状态转换的先验信息；3、利用不同传感器测量的互补性，提高了估计精度；4、计算量较小，精度较好，可以实现在线估计。附图说明以下附图描述了本专利技术所选择的实施例，均为示例性附图而非穷举或限制性，以电导传感器测量油水两相流为例，其中：图1本专利技术实施例采用的电导传感器结构示意图。图2本专利技术的在线自适应卡尔曼估计算法计算步骤。具体实施方式以电导传感器为具体实施例，不局限于电导传感器，其它类型的传感器同样适用。当信号频率在10kHz到100kHz之间时，油水两相流的负载特性体现为电阻性，即电极对之间的电压与分布其间的流体阻抗成正比。在油水两相流流体可以形成回路时，油水两相流的阻抗与相分布有对应关系，因此可以通过测量电极对之间的电压来测量其相含率。如图1所示，由六本文档来自技高网...

【技术保护点】
一种两相流含水率自适应串联估计方法，利用上、下游不同传感器的测量数据进行估计，包括如下步骤：1)对上游传感器和下游传感器测量数据进行时间配准：获取上游传感器和下游传感器所获相含率时间序列z

【技术特征摘要】
1.一种两相流含水率自适应串联估计方法，利用上、下游不同传感器的测量数据进行估计，包括如下步骤：1)对上游传感器和下游传感器测量数据进行时间配准：获取上游传感器和下游传感器所获相含率时间序列z1(j)和z2(j)，j＝1,2,3,…N为测量时间内的数据点，下标1和2分别代表上游传感器和下游传感器，计算渡越时间，利用渡越时间结合数据采样频率，对上游传感器和下游传感器测量数据进行时间配准；2)状态估计算初始化：其中，是上游状态向量的估计初值，是上游协方差的迭代初值，z1,k是上游传感器在k时刻的测量值，且是下游状态向量的估计初值；3)对上游传感器含率测量值进行线性卡尔曼状态估计和更新，获得对状态向量的预测和协方差矩阵的预测；在估计过程中量测矩阵H取值均为1，状态向量x1代表含水率，k时刻的系统模型为：x1,k+1＝F1,kx1,k+v1,k(4)z1,k＝H1,kx1,k+w1,k(5)其中，F1,k是k时刻的状态转移矩阵，x1,k+1和x1,k分别是上游在k+1时刻和k时刻的状态向量，H1,k是k时刻的量测矩阵，z1,k是量测向量，v1,k和w1,k分别表示系统在k时刻的过程噪声和量测噪声；则状态向量的预测的预测和协方差矩阵的预测为：根据上游量测信息z1,k修正估计值：其中，K1,k是上游传感器系统在k时刻的卡尔曼增益：其中，H1,k是k时刻上游传感器量测矩阵，是上游传感器在k时刻对状态的预测值，是上游传感器k时刻对状态x1,k的最小均方误差估计，是k时刻估计的协方差，R1,k为上游滤波过程在k时刻系统量测噪声的协方差，可取值为0.01至0.2；Q1,k-1为上游滤波过程在k-1时刻系统过程噪声的协方差，可取值为0.01至0.09，I为单位矩阵，其维数等同于状态向量的维数；3)利用下游传感器测量数据，计算两相流过程状态转换矩阵：将传感器上游测量单元确定的状态转移过程加以记录，用于时间配准后下游含水率状态估计中：其中，F2,k为下游测量单元的状态转移矩阵，k＝n0,n0+1,…,N；4)对下游含率测量值...

【专利技术属性】
技术研发人员：谭超，付广智，董峰，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津,12