皮禮仕

（中國市政工程中南設(shè)計研究總院有限公司，湖北 武漢 430010）

0 引言

截至2017年底，我國已建成陜京線、西氣東輸、川氣東送、中緬線等骨干管網(wǎng)系統(tǒng)，支線管道正實現(xiàn)全國覆蓋，總里程約為7.7×104km[1-2]。面對輸氣規(guī)模逐年增大，運行復(fù)雜度逐漸增加的天然氣管網(wǎng)，對其運行參數(shù)進(jìn)行有效優(yōu)化，從而實現(xiàn)輸氣管道的經(jīng)濟運行具有重要意義[3]。輸氣管道運行優(yōu)化研究是對管道系統(tǒng)運行方案進(jìn)行優(yōu)化分析的過程。在目標(biāo)函數(shù)確定的前提下，根據(jù)水力、熱力，壓縮機性能曲線等約束條件，求解出管道系統(tǒng)在一定輸送條件下的最省運行方案。在上世紀(jì)60年代，歐美等國已開展了該課題的研究，并開發(fā)了適用于管道穩(wěn)態(tài)工況的優(yōu)化算法[4]，以壓氣站總能耗最低為目標(biāo)函數(shù)，計算確定各站最優(yōu)開機方案。1968年，Peter首次將動態(tài)規(guī)劃法應(yīng)用于單條輸氣管道總能耗最低的優(yōu)化[5]，為該方法的推廣奠定了基礎(chǔ)，到60年代末期，動態(tài)規(guī)劃法已經(jīng)成為解決輸氣管道運行優(yōu)化的主流算法。它能較好處理多階段決策中的非線性約束條件，優(yōu)化結(jié)果具備全局最優(yōu)性。

1 模型的建立與求解

1.1 數(shù)學(xué)模型的建立

壓縮機是壓氣站在運行中的主要能耗裝置，其經(jīng)營成本占輸氣管道運營成本的40%～50%，研究壓氣站能耗最優(yōu)具有重要意義[6-7]。假設(shè)某管道有n座壓氣站，配備壓縮機m臺，同型號的機組并聯(lián)運行下流量均分時最優(yōu)[8]，以各壓氣站的總能耗費用之和最小建立目標(biāo)函數(shù)，并根據(jù)管道水力、熱力及壓縮機等約束條件建立具體數(shù)學(xué)模型[9-10]如下：

式中，C為n座壓氣站的能耗總和，MW；Ci為第i座壓氣站能耗，MW；fg為管道水力、熱力仿真計算函數(shù)，計算各站進(jìn)站壓力和溫度；fz為壓氣站仿真優(yōu)化計算函數(shù)，計算各站出站溫度、開機方案及能耗費用；為第i站的進(jìn)站壓力，MPa；為第i站的出站壓力，MPa；TS,(i)為第i站的進(jìn)站溫度，K；Td,(i)為第i站的出站溫度，K；Qi為第i站的進(jìn)站流量，m3/s；Qi，j為第i站第j臺壓縮機的流量，m3/s；Di，j為第i站j編號壓縮機的可行域約束，無因次。

1.2 動態(tài)規(guī)劃優(yōu)化算法

單條輸氣管道的穩(wěn)態(tài)運行優(yōu)化，其決策變量為各站開機方案，利用動態(tài)規(guī)劃法分層級求解該問題，可可將其分為兩級優(yōu)化過程，依此為管線級和站級。

1）管線級優(yōu)化求解過程

對沿途設(shè)置n座壓氣站的管道，通過動態(tài)規(guī)劃法，將該管道系統(tǒng)的優(yōu)化過程分為n個階段，將第k-1站的出站到第k的出站視為第k階段，第k-1站的出站狀態(tài)變量用xk表示，那么第1階段為首站的開機方案決策，如圖1所示。將第k階段壓氣站的出站壓力Pd，（k）和第k站壓比dk的確定看成一個決策過程，則全線所有站的出站壓力確定可以看成一個多階段決策問題，應(yīng)用動態(tài)規(guī)劃法求解很方便，建立輸氣管道各站壓比及出站壓力優(yōu)化動態(tài)規(guī)劃模型如下[11-13]：

由于氣源的壓力p0，初始邊界條件為：

式中，pd,(0)為第一個壓氣站的入口壓力，MPa；根據(jù)順序遞推法計算出全線最小能耗為Fn（xn+1），MW。

圖1 輸氣管道動態(tài)規(guī)劃示意圖

對上述動態(tài)規(guī)劃模型求解步驟如下：①確定各站出站壓力上下限，并對其進(jìn)行離散化處理，得到k-1站 有 限 個 出 站 壓 力 狀 態(tài) 點｛Pd，（k-1），1，Pd，（k-1），2，…，Pd，（k-1），n｝，計作狀態(tài)參數(shù)xk；第k-1站的最小出站壓力可以由k站的最小要求進(jìn)站壓力Ps，（k），min反算求得，因此利用終點最低要求用氣壓力Pz，min依次反算可求得各站最低要求出站壓力。②利用管道水力、熱力仿真計算下一站進(jìn)站壓力，得到k站有限個進(jìn)站壓力狀態(tài)點｛Ps，（k），1，Ps，（k），2，…，Ps，（k），n｝，此時即可篩選掉小于Ps，（k），min的狀態(tài)點，減少k站站內(nèi)遞推匹配次數(shù)。③站內(nèi)遞推，對于k站某一確定出站狀態(tài)Pd，（k），i，將該站 進(jìn) 站狀態(tài)點｛Ps，（k），1，Ps，（k），2，…，Ps，（k），n｝與之依次匹配，首先判斷從Ps，（k），1到pd，（k），i的壓比di，1是否在壓縮機可行域內(nèi)，若不在，則令該路徑費用為一極大值；若在可行域內(nèi)，則調(diào)用站級優(yōu)化模塊，計算該路徑階段費用，將該費用與Ps，（k），1攜帶的 從 氣 源 到 該 點 最 低 指 標(biāo) 費 用Fk-1（xk，1）相 加 得 到Fk（xk+1，i），并記錄相應(yīng)開機方案。同理依次將其他進(jìn)站狀態(tài)點Ps，（k），j與Pd，（k），i匹配，并計算階段費用，若得到Fk（xk+1，i）的值更低，則記錄該值，并更新到達(dá)Pd，（k），i的決策路徑di，j（i代表k站第i個出站壓力，j代表對應(yīng)的最優(yōu)進(jìn)站壓力編號）。依此類推，找出第k個壓氣站所有出站壓力從氣源到該點最低指標(biāo)費Fk（xk+1）及對應(yīng)決策壓比dk。④在輸氣末端，將各方案對應(yīng)的指標(biāo)費用從小到達(dá)依次排序，記錄能耗費用值最低的那組可行方案，并從末站到首站依次逆推各站相應(yīng)的開機方案。

2）站級優(yōu)化求解過程

站級優(yōu)化是輸氣管道穩(wěn)態(tài)運行優(yōu)化中最低層級的優(yōu)化，該優(yōu)化模塊在整個優(yōu)化過程中的調(diào)用次數(shù)最多，每當(dāng)壓氣站有一組確定的進(jìn)出站狀態(tài)參數(shù)時，就會調(diào)用該模塊進(jìn)行開機方案優(yōu)化計算。采用壓縮機流量均分方式，如圖2所示。另外，對離心壓縮機壓頭—流量曲線、效率—流量曲線、喘振流量—轉(zhuǎn)速曲線、滯止流量—轉(zhuǎn)速曲線的求解均采用最小二乘法。

圖2 站級優(yōu)化算法圖

2 某管道運行優(yōu)化案例

輸氣管道仿真計算是穩(wěn)態(tài)運行優(yōu)化的基礎(chǔ)，仿真計算是對各站進(jìn)行一次確定性方案的模擬運算，而運行優(yōu)化則需進(jìn)行有限種方案的篩選，并比較出最優(yōu)方案。二者都需調(diào)用輸氣管道水力、熱力仿真模塊和壓氣站仿真模塊，只不過運行優(yōu)化調(diào)用基礎(chǔ)模塊的次數(shù)要更頻繁。

2.1 基礎(chǔ)數(shù)據(jù)

管線全長為1 100 km，管徑為1 422 mm，設(shè)計壓力為12 MPa，全線設(shè)置壓氣站4座，沿途分輸點4處，設(shè)計輸量為3.60×1010m3/a，各站壓縮機組型號一樣，每站設(shè)置3臺壓縮機，均為電力驅(qū)動。利用最小二乘法對4個壓氣站的壓縮機壓頭—流量曲線進(jìn)行二項式擬合得擬合系數(shù)A，同理對效率—流量曲線、喘振流量—轉(zhuǎn)速曲線、滯止流量—轉(zhuǎn)速曲線進(jìn)行二項式擬合[14]，擬合系數(shù)分別為B、C、D，如表1所示。

2.2 優(yōu)化計算

由于管道在達(dá)到設(shè)計輸量后基本沒有優(yōu)化空間[15]，研究采用對運行流量為設(shè)計輸量85%的工況對管道全線進(jìn)行優(yōu)化計算，并利用壓氣站優(yōu)化模塊與水力計算模塊設(shè)置一組自定義開機方案，與優(yōu)化方案作對比。同時將該自定義方案用國外管網(wǎng)仿真軟件TGNET進(jìn)行建模，TGNET軟件曾用于西氣東輸、川氣東送等大型管道和川西管網(wǎng)的模擬計算[16-17]。對比計算結(jié)果，以驗證本次優(yōu)化仿真程序的準(zhǔn)確性。優(yōu)化仿真結(jié)果見表2。

表1 壓縮機曲線方程二項式擬合系數(shù)表

表2 動態(tài)規(guī)劃法優(yōu)化仿真模擬結(jié)果統(tǒng)計表

由表2模擬結(jié)果可知，TGNET自定義方案與優(yōu)化仿真自定義方案計算結(jié)果較為接近，兩方案相比可知：壓氣站1的開機方案完全相同，功率計算的相對誤差僅為0.4%；壓氣站4的出站壓力設(shè)置相同，末站的終點壓力計算相對誤差為3%。綜上可知優(yōu)化仿真程序模擬精度完全符合工程要求。

優(yōu)化方案下的開機策略總能耗較自定義方案降低了8.8%，對比開機方案可看出，雖然自定義方案開機總臺數(shù)減少了，但不合理的越站方案導(dǎo)致了更高的總能耗，驗證了所研究算法的優(yōu)化性。

3 結(jié)論

1）對輸氣管道穩(wěn)態(tài)運行優(yōu)化進(jìn)行研究能為已投運管道優(yōu)化方案的制定提供有效指導(dǎo)，研究提出了以總能耗最低為目標(biāo)函數(shù)的輸氣管道運行優(yōu)化模型，并利用動態(tài)規(guī)劃法進(jìn)行管線級優(yōu)化和站級優(yōu)化求解，已成功應(yīng)用于某管道，求解出其在85%設(shè)計輸量運行下的各站最優(yōu)開機方案。

2）根據(jù)自定義方案與TGNET軟件自定義方案計算結(jié)果對比表明了優(yōu)化仿真程序的模擬精度符合工程實際要求，根據(jù)自定義方案與優(yōu)化方案計算對比結(jié)果，可使全線壓縮機總功率有效降低。

3）研究建立的輸氣管道穩(wěn)態(tài)運行優(yōu)化模型不僅適用于已投產(chǎn)的管道，也適用于前期設(shè)計階段，利用設(shè)計階段推薦方案的典型運行工況，結(jié)合廠家匹配出合適型號的壓縮機，就可提前對該管道進(jìn)行多種典型工況的方案優(yōu)化，將設(shè)計與優(yōu)化相結(jié)合，提高設(shè)計方案的合理性。