基于最小二乘法的注水管网摩阻因数反演.pdfVIP

第 6期 王玉学等 ：基于最小二乘法的注水管网摩 阻因数反演 方程 (2)的矩阵向量形式为 Aq— Q． (3) 记管道铺设时的摩阻为经验摩阻C0，利用已知节点压力、压降方程和经验摩阻 可以求 出管道经验 流量 q。[{；]． 寻找能够满足方程(2)，又能使 llq--q。l1达到最小的q，求出q后，利用已知压力和压降方程可以得 到管道摩阻因数．以此摩阻因数计算节点压力与已知压力相等，且能使得计算摩阻与经验摩阻在某种意义 下最为接近．因此，该问题的数学模型为 minll口一口。lI2， q (4) S．t．Aq— Q． 又由于方程 (3)有解 ，且A的秩小于未知量 的个数 ，因此其解与其最小二乘解一致 ，均有无穷多解 ，且 任一解都能使 lIQ—Aq 一O．因此 ，问题的数学模型可以改为求q，使得 rainlIq—q0ll2， q ‘ (5) S．t．rainlJQ—Aqll；． 为了方便利用亏秩最小二乘法求解 ，需对所求变量进行变形．令 q一口一qo，Q一Q--Aq。，则问题的数 学模型变成求 q，使得 rainllq，ll2， … S．t．rainllQ一Aq，ll；． 若能求出模型(6)的解 q，则可得 q=q+q。，利用 已知点压力和压 降方程可 以求出管道摩阻因数 C_ 2 模 型求解 数学模型(6)属于标准的求最小二乘最小范数解的问题 ，文献E9]~EN这类最小二乘最小范数解是存 在并且唯一的． 将矩阵A做奇异值分解A—u[：IVT,则UTAV=[吾：]，所以有 llQ一Aq，ll一 IlU (Q一Aq)ll；一 llUQ一(uAV)(Vq))II一 lI(：)一(：)()==：lI_『c 一c·一y“+--czi“． c7 式中：(：)一，【Q，()一V口· 只要)．，一三c，任取，所对应的q一vf1必为方程组Aq一Q的最小二乘解．又因为II口II一 Il(，为了使口既是最小二乘解，并且范数最小，需取z为o．所求的最小--~@1／、范数解s为 一 v【uJ1． (8) 利用奇异值分解求最小二乘最小范数解的方法 ，不仅能够处理秩亏损的最小二乘问题 ，而且对于方程 的个数小于未知量 的个数 、方程的个数大于未知量的个数 、方程的个数等于未知量的个数 (包括系数矩阵 可逆和不可逆的情况)等能够求其最 ／]~--乘最小范数解． 利用单工况数据反演管道摩阻因数步骤 ： (1)利用已知数据确定系数矩 阵A． (2)利用 已知压力和经验摩 阻因数确定各管道经验流萤 q。． 大 庆 石 油 学 院 学 报 第 35卷 2011年 (3)求满足方程 (6)的解 q；此处利用亏秩最小二 乘法求解． (4)令q—q+g。． (5)由5—lH 一H，l／Iql 求得 5． (6)由 s一 求得各管道摩阻因数c ，其 中l为该管段长度 ；d 为该管段直径． 3 实例计算 某一理想注水管网由2座注水站、l6个节点、24条 管道、9个环构成 ，节点编号 2和 15为泵站所在位置 ， 18 见图 1． 首先利用实际摩阻因数