本发明专利技术涉及一种油水两相流含水率自适应估计方法，包括：将来自电导传感器上、下游的测量数据进行时间配准；基于自适应更新的状态估计；判定两相流是否处于水连续的状态；对传感器的含率测量值进行卡尔曼状态估计和更新。本发明专利技术能够提升含水率的估计精度。
【技术实现步骤摘要】
本专利技术属于流体测量
，利用电导传感器与自适应卡尔曼算法测量含水率信息，用于提高两相流含水率估计精度。技术背景两相流广泛存在于日常生活与工业生产过程中，如食品加工，生物工程，化工产业，冶金工业和石油产业等行业。与一般单相流相比，两相流的流动状态更为复杂，建模更加困难，因此两相流的在线过程参数检测一直是学术界和工业界关注的难点和重点。特别在石油工业中，油水两相流含水率对油井监测、生产预测和开采计划制定具有重要意义。目前两相流含水率的测量手段包括快关阀测量方法、射线方法、微波方法、Wire-mesh方法、超声方法和电阻抗方法等。其中，电学方法具有测量速度快，敏感性强，传感器结构简单，易于安装，成本低廉等优点。两相流是一种随机的动态流动过程，流动状态复杂多变，传感器的测量信号有随机误差，且只能反映系统的外部特性。为获得多相流含水率的准确信息，需要根据实测数据对流动状态进行估计。状态估计对认知和分析测量系统具有重要的应用价值，常用的估计理论包括最小二乘估计、线性最小方差估计和卡尔曼估计等。不同的估计方法性能指标不同，使用的信息不同，其适用范围也不同。最小二乘估计通常应用于随机向量和常值向量的估计，算法结构简单但精度有限；线性最小方差估计以均方误差最小作为优化参数，适用于平稳过程，其测量信息越丰富，估计精度越高，计算量越大，而且会随着估计过程的累积，造成计算负担较重，不适用于复杂的长时间的估计过程。卡尔曼估计采用递推的形式，通过一个状态模型和一个观测模型的迭代，利用上一时刻的估计值和这一时刻的观测值，对这一刻的状态进行估计。
技术实现思路
本专利技术提供一种油水两相流含水率自适应估计方法，该方法利用阵列式电导传感器上下游电极测量数据，实时更新卡尔曼估计迭代过程的相关参数，进行自适应状态转换，提高含水率的估计精度。本专利技术的技术方法如下：一种油水两相流含水率自适应估计方法，利用阵列式电导传感器上下游电极测量数据进行估计，包括如下步骤：1)将来自电导传感器上、下游的测量数据进行时间配准：将取自上游和下游测量单元的含率测量值z1(n)和z2(n)看作来自各态历经的平稳随机过程的两个样本函数，其中n＝1,2,3,…N为测量时间内的数据点，两传感器间距为L，设渡越时间为τ0，对应离散系统中上、下游测量数据平移组数n0；计算z1(n)和z2(n)在离散条件下的互相关函数公式如下：的最大值出现在n0处，即互相关函数峰值对应流体从上游测量单元到下游测量单元的渡越时间，通过传感器上游测量单元获取的前、后时刻的实测系统状态向量确定状态转移矩阵，反应上游测量所得的系统状态矩阵的转化情况：其中，z1(k-1)和z1(k)是取自上游测量单元，分别为测量区域内前后时刻的含率测量值；利用n0结合电导传感器采样频率，对上下游测量数据进行时间配准，如公式(3)所示：z2(k)＝z'2(k+n0)(3)其中，z'2为下游传感器单元采集数据，z2是配准后的数据；2)基于自适应更新的状态估计：将由传感器上游测量单元确定的状态转移矩阵用于时间配准后的下游含水率状态更新中：其中，F2(k)为传感器下游测量单元状态转移矩阵，F1(k)为上游测量单元的状态转移矩阵；传感器测量数据时间配准后，依据传感器上下游具体数据对F(k)进行实时的自适应更新，如公式(5)所示：3)判定两相流是否处于水连续的状态：将电导传感器下游测量单元k时刻含率测量值z2(k)带入公式(6)，计算k时刻的量测噪声w(k)，将量测噪声w(k)带入公式(8)，计算流态判别系数q(k)；z2(k)-H(k)x(k|k-1)＝w(k)(6)S＝H(k)p(k|k-1)+R(k)(7)q(k)＝wT(k)S-1w(k)(8)公式(6)中x(k|k-1)是k时刻系统状态预测值，估计方式如公式(9)所示；H(k)是量测矩阵，取值为1；公式(7)中，R(k)是k时刻量测噪声的方差，可取传感器含率测量数据的标准差作为恒值代入；p(k|k-1)是传感器k-1时刻对k时刻估计的协方差，计算方式如公式(10)所示；流态判别系数q(k)在卡尔曼估计过程中针对具体对象衡量环境信息的门限值，当q过大，代表流动状况处于非水相连续的状态，此处，q的限值取0.118。4)对传感器的含率测量值进行卡尔曼状态估计和更新：x(k|k-1)＝F(k)x(k-1|k-1)(9)p(k|k-1)＝F(k-1)p(k-1|k-1)FT(k-1)+Q(k-1)(10)根据量测信息z(k)修正估计值：x(k|k)＝x(k|k-1)+K(k)×[z(k)-H(k)x(k|k-1)](11)p(k|k)＝[I-K(k)H(k)]p(k|k-1)(12)其中，K(k)是传感器系统在k时刻的卡尔曼增益：K(k)＝p(k|k-1)HT(k)×[H(k)p(k|k-1)HT(k)+R(k)]-1(13)其中，p(k-1|k-1)是k-1时刻的协方差，其迭代初值取含水率测量值的方差；z(k)是传感器在k时刻对含水率的测量值，x(k|k-1)是传感器在k时刻对状态x(k)的预测值，x(k|k)是传感器k时刻对状态x(k)的最小均方误差估计，p(k|k)是对k时刻估计x(k|k)的协方差，Q(k-1)为k-1系统过程噪声的方差，可取恒值为0.01，I为单位矩阵，其维数等同于状态向量。本专利技术的实质性特点是：针对多相流流动过程，采用电导传感器上下游测量所得数据，在测量精度难以提高的情况下，利用可以对被测对象估计过程中的状态转移矩阵进行实时更新和变化的特点，改进卡尔曼估计算法的更新和迭代过程，进行自适应卡尔曼估计。利用“不确定度补偿”提高估计精度。本专利技术的有益效果及优点如下：1、该方法结合传感器具体结构，充分利用传感器上游和下游测量数据；2、该方法相对于传感器直接测量，提高了估计精度；3、该方法在计算量较小，可以实现实时在线估计。附图说明以下附图描述了本专利技术所选择的实施例，均为示例性附图而非穷举或限制性，以电导传感器测量油水两相流为例，其中：图1本专利技术的测量方法电导传感器测量结构示意图，(a)为正视图，(b)为侧视图，(c)为立体图；图2本专利技术的在线自适应卡尔曼估计算法计算步骤。具体实施方式当信号频率在10kHz到100kHz之间时，油水两相流的负载特性体现为电阻性，即电极对之间的电压与分布其间的流体阻抗成正比。在油水两相流流体可以形成回路时，油水两相流的阻抗与相分布有对应关系，因此可以通过测量电极对之间的电压来测量其相含率。下面结合说明书附图说明本专利技术的计算方法。图1为传感器测量空间示意图。电导传感器实物图如图1所示，详细流程如图2所示，由六个沿轴向排列的环形电极组成，其中1、6为激励电极，2、3、4、5为测量电极。采用电流激励电压测量的方式，通过向电极1中注入恒定电流并将电极6接地，可在电极1和6中间形成稳定的电学测量场，分别测量电极对2-3、电极对3-4以及电极对4-5之间的电势差V可获得其间的两相流平均电阻率。电导传感器为6环结构电导环传感器阵列以轴向排列的6个环状金属电极嵌入管路中的形式构成，管径D为50mm，相邻2个环状电极的距离分别为：1D、0.4D、1.2D、0.4D、1D。在此，依照图2所示流动方向，选取电极2和电极3之间的区域为上本文档来自技高网...
【技术保护点】
一种油水两相流含水率自适应估计方法，利用阵列式电导传感器上下游电极测量数据进行估计，包括如下步骤：1)将来自电导传感器上、下游的测量数据进行时间配准：将取自上游和下游测量单元的含率测量值z1(n)和z2(n)看作来自各态历经的平稳随机过程的两个样本函数，其中n＝1,2,3,…N为测量时间内的数据点，两传感器间距为L，设渡越时间为τ0，对应离散系统中上、下游测量数据平移组数n0；计算z1(n)和z2(n)在离散条件下的互相关函数公式如下：Rz1z2(n)=1NΣm=0N-1z1(m)z2(m+n)---(1)]]>Rz1z2(n)的最大值出现在n0处，即互相关函数峰值对应流体从上游测量单元到下游测量单元的渡越时间，通过传感器上游测量单元获取的前、后时刻的实测系统状态向量确定状态转移矩阵，反应上游测量所得的系统状态矩阵的转化情况：F1(k)=z1(k)z1(k-1)---(2)]]>其中，z1(k‑1)和z1(k)是取自上游测量单元，分别为测量区域内前后时刻的含率测量值；利用n0结合电导传感器采样频率，对上下游测量数据进行时间配准，如公式(3)所示：z2(k)＝z'2(k+n0) (3)其中，z'2为下游传感器单元采集数据，z2是配准后的数据；2)基于自适应更新的状态估计：将由传感器上游测量单元确定的状态转移矩阵用于时间配准后的下游含水率状态更新中：F2(k)=F1(k)=z1(k)z1(k-1)---(4)]]>其中，F2(k)为传感器下游测量单元状态转移矩阵，F1(k)为上游测量单元的状态转移矩阵；传感器测量数据时间配准后，依据传感器上下游具体数据对F(k)进行实时的自适应更新，如公式(5)所示：F2(k+n0)=z1(k)z1(k-1)---(5)]]>3)判定两相流是否处于水连续的状态：将电导传感器下游测量单元k时刻含率测量值z2(k)带入公式(6)，计算k时刻的量测噪声w(k)，将量测噪声w(k)带入公式(8)，计算流态判别系数q(k)；z2(k)‑H(k)x(kk‑1)＝w(k) (6)S＝H(k)p(k|k‑1)+R(k) (7)q(k)＝wT(k)S‑1w(k) (8)公式(6)中x(k|k‑1)是k时刻系统状态预测值，估计方式如公式(9)所示；H(k)是量测矩阵，取值为1；公式(7)中，R(k)是k时刻量测噪声的方差，可取传感器含率测量数据的标准差作为恒值代入；p(k|k‑1)是传感器k‑1时刻对k时刻估计的协方差，计算方式如公式(10)所示；流态判别系数q(k)在卡尔曼估计过程中针对具体对象衡量环境信息的门限值，当q过大，代表流动状况处于非水相连续的状态，此处，q的限值取0.118。4)对传感器的含率测量值进行卡尔曼状态估计和更新：x(k|k‑1)＝F(k)x(k‑1|k‑1) (9)p(k|k‑1)＝F(k‑1)p(k‑1|k‑1)FT(k‑1)+Q(k‑1) (10)根据量测信息z(k)修正估计值：x(k|k)＝x(k|k‑1)+K(k)×[z(k)‑H(k)x(k|k‑1)] (11)p(k|k)＝[I‑K(k)H(k)]p(k|k‑1) (12)其中，K(k)是传感器系统在k时刻的卡尔曼增益：K(k)＝p(k|k‑1)HT(k)×[H(k)p(k|k‑1)HT(k)+R(k)]‑1 (13)其中，p(k‑1|k‑1)是k‑1时刻的协方差，其迭代初值取含水率测量值的方差；z(k)是传感器在k时刻对含水率的测量值，x(k|k‑1)是传感器在k时刻对状态x(k)的预测值，x(k|k)是传感器k时刻对状态x(k)的最小均方误差估计，p(k|k)是对k时刻估计x(k|k)的协方差，Q(k‑1)为k‑1系统过程噪声的方差，可取恒值为0.01，I为单位矩阵，其维数等同于状态向量。...
【技术特征摘要】
1.一种油水两相流含水率自适应估计方法，利用阵列式电导传感器上下游电极测量数据进行估计，包括如下步骤：1)将来自电导传感器上、下游的测量数据进行时间配准：将取自上游和下游测量单元的含率测量值z1(n)和z2(n)看作来自各态历经的平稳随机过程的两个样本函数，其中n＝1,2,3,…N为测量时间内的数据点，两传感器间距为L，设渡越时间为τ0，对应离散系统中上、下游测量数据平移组数n0；计算z1(n)和z2(n)在离散条件下的互相关函数公式如下：Rz1z2(n)=1NΣm=0N-1z1(m)z2(m+n)---(1)]]>Rz1z2(n)的最大值出现在n0处，即互相关函数峰值对应流体从上游测量单元到下游测量单元的渡越时间，通过传感器上游测量单元获取的前、后时刻的实测系统状态向量确定状态转移矩阵，反应上游测量所得的系统状态矩阵的转化情况：F1(k)=z1(k)z1(k-1)---(2)]]>其中，z1(k-1)和z1(k)是取自上游测量单元，分别为测量区域内前后时刻的含率测量值；利用n0结合电导传感器采样频率，对上下游测量数据进行时间配准，如公式(3)所示：z2(k)＝z'2(k+n0)(3)其中，z'2为下游传感器单元采集数据，z2是配准后的数据；2)基于自适应更新的状态估计：将由传感器上游测量单元确定的状态转移矩阵用于时间配准后的下游含水率状态更新中：F2(k)=F1(k)=z1(k)z1(k-1)---(4)]]>其中，F2(k)为传感器下游测量单元状态转移矩阵，F1(k)为上游测量单元的状态转移矩阵；传感器测量数据时间配准后，依据传感器上下游具体数据对F(k)进行实时的自适应更新，如公式(5)所示：F2(k+n0)=z1(k)z1(k-1)---(5)]]>3)判定两相流是否处于水连续的状态：将电导传感器下游测量单元k时刻含率测量值z...
【专利技术属性】
技术研发人员：谭超，付广智，董峰，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津;12
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。