導語：供水管網(wǎng)改應用論文一文來源于網(wǎng)友上傳，不代表本站觀點，若需要原創(chuàng)文章可咨詢客服老師，歡迎參考。

摘要：本文將字典序線性目標規(guī)劃基本模型應用于大型供水管網(wǎng)改、擴建優(yōu)化決策，建立了大型供水管網(wǎng)改、擴建工程優(yōu)化決策字典序線性目標規(guī)劃模型，并用字典序線性目標規(guī)劃多階段算法進行計算，編制了相應的計算程序，且根據(jù)其計算特點，編制了字典序線性目標規(guī)劃靈敏度分析程序，對所建立的模型進行靈敏度分析。計算實例表明：字典序線性目標規(guī)劃法的理論和方法非常適合大型供水管網(wǎng)改、擴建工程優(yōu)化決策計算，結(jié)果令人滿意。所編制程序具有準確性和通用性，對生產(chǎn)具有科學和實際的指導意義。

關(guān)鍵詞：字典序目標規(guī)劃供水管網(wǎng)靈敏度

TheLexicographicLinearGoalprogrammingisappliedtothe

OptimalDecision-makingofthelargeWaterNetworks

Rebuliding&Enlargingproject

Abstract：ThebasicmodelofthelexicographiclinearGoalprogrammingisappliedtocalculatingtheoptinaldecision-makingofthelargewaternetworksrebuliding&enlargingproject.Amodelusedforcalculatingtheoptimaldecision-makingofthelargewaternetworksrebuliding&enlargingprojectisconstructed.ThelexicographiclinearGoalprogrammingmutiplestagealgorithmisusedtocalculate.A

correlativecalculatingprogramisprogrammed.Accordingtoitscalculatingcharacteristics,theprocedure

ofthelexicographiclinearGoalprogrammingsensitivityanalysisismadetoanalysetheconstructed

model.Thecalculatingexampleindicatesthat，thetheoryandthemethodofthelexicographiclinearGoal

programmingisverysuitableforthecalculationoftheoptimaldecision—makingofthelargewaternetworksrebuliding&enlargingproject.Theresultprovessatisfactory.Theprogramsarepreciseandstable.Theyhavescientificandpracticalguidingvaluetoproduction.

keywords：LexiGographicGoalprogrmmingwaternetworkssensitivity

1字典序線性目標規(guī)劃基本模型[1][2]

1.1字典序線性目標規(guī)劃方法的概念

字典序線性目標規(guī)劃(LexicopraphicLinearGoalprogramming)簡稱LLGP是美國阿波羅計劃工程師J.P伊格尼齊奧于1976年在研究目標規(guī)劃法的基礎(chǔ)上提出的一種用目標偏差量有序向量衡量多目標模型目標達成度的字典序最小方法。

1.2數(shù)學模型

LLGP方法的數(shù)學模型為：

求x（x1,x2,…,xn）

使Lexminα={g1(η,δ),g2(η,δ),…,gk(η,δ)}（1-1-1）

s.t∑nj=1=cI,jxj+ηi-δi=bi，i=1,2,…,m（1-1-2）

x,η,δ≧0（1-1-3）

其中：x—n維決策變量；η—負偏差變量；δ—正偏差變量；gk(η,δ)—為第K個優(yōu)先級別中目標偏差變量的線性系數(shù)：α—尋求字典序極小化達成系數(shù)；cI,j—是對應于第i個目標或約束中第j個變量的系數(shù)；bi—為目標或約束分端常數(shù)：（1-1-1）為LLGP問題的目標；（1-1-2）為目標的約束條件，可以有硬約束和軟約束；（1-1-3）為非負約束，即x、η與δ均為非負變量；k—為優(yōu)先次序數(shù)。m—為全部目標和約束數(shù)。

1.3字典序線性目標規(guī)劃模型的求解方法：

字典序線性目標規(guī)劃模型最常用的求解方法是多階段單純形算法，它的目標行里有多個目標存在，進行單純形迭代時，按確定的優(yōu)先次序分層次的解決問題，它得到多個目標滿足約束條件下的滿意解。

2字典序線性目標規(guī)劃模型靈敏度分析方法[3]

2.1靈敏度分析

靈敏度分析是考查模型結(jié)構(gòu)中，某元素變化對問題解的影響。

2.2字典序線性目標規(guī)劃靈敏度分析的三種方法：

(1)結(jié)構(gòu)性離散變化：①重捧優(yōu)先級②增加新變量⑧增加變目標；

(2)非結(jié)構(gòu)性離散變化：有vh+毗k，yl+及h的離散變化：

(3)在一個區(qū)域連續(xù)取值的參數(shù)分析：即參數(shù)線性目標規(guī)劃。

3大型供水管網(wǎng)改、擴建工程優(yōu)化決策字典序線性目標規(guī)劃計算模型

大型供水管網(wǎng)改、擴建工程優(yōu)化決策問題的目標是多個的，既有總水量目標，又有各輸送管道及轉(zhuǎn)壓站輸送水量目標，還有控制費用目標，而且這些目標的重要程度是不同的：首先，須滿足總水量目標，其次是控制費用及各輸送輸送管道及轉(zhuǎn)壓站輸送水量目標，除滿足上述目標外，還須滿足在每種輸送方式下各種輸送管道及增壓站的輸送水量小于或等于該種輸送方式下的輸送水量。下述表明：大型供水管網(wǎng)改、擴建工程優(yōu)化決策問題，實際上是個多目標規(guī)劃問題，其約束和目標方程都是線性的，因而完全適合字典序線性目標規(guī)劃法求解條件。

3.1已知條件：

(1)所考察的各類水量；S1，S2，……，Sp；

(2)每種輸送方式費用單價：C11，C12，…，C1q1；…，…，…，…；

Cp1，Cp2，…，Cpqp；

(3)各條輸送管道及各轉(zhuǎn)壓站的輸送水量；A1，A2，……，Ar

(4)所有改.擴建管道及轉(zhuǎn)壓站項目的控制費用：W

3.2待求變量：

(1)每類水量適宜的輸送方式輸送該類水量的量：

X11，X12，…，X1q1；…，…，…，…；xp1，xp2，…，xpqp；

(2)每種輸送方式下?lián)Q輸管道和轉(zhuǎn)壓站所輸送水量：

x111，x112，…，x11r11；…，…，…，…；xpqp1，xpqp2，…，xpqprpqp；

(3)各種輸送方式下各種輸送管道和轉(zhuǎn)壓站的輸送水量：；

x11A1，x11A2，…，x11Ar；…，…，…，…；xpqpA1，xpqpA2，…，xpqpAr；

3.3計算模型：

3.3.1目標：

(1)每類水量各輸送方式所輸送水量之和與該類水量相等，即：

∑q1j1=1x1j1=S1；∑q2j2=1x2j2=S2；…；∑qpjp=1xpjp=Sp；（3-1-1）

(2)輸送各類水量各輸送方式管道及轉(zhuǎn)壓站改、擴建費用之和等于總控制費用：

∑q1j1=1C1j1x1j1+∑q2j2=1C2j2x2j2+…∑qpjp=1Cpjpxpjp=w；（3-1-2）

(3)輸送各類水量同一種管道求轉(zhuǎn)壓站所輸送水量應該等于該管道求轉(zhuǎn)壓站所輸送的水量：

∑q1j1=1x1j1A1+∑q2j2=1x2j2A1+…∑qpjp=1xpjpA1；…；∑q1j1=1x1j1Ar+∑q2j2=1x2j2Ar+…∑qpjp=1xpjpAr=Ar

3.3.2約束條件：

(1)水量約束：

每種輸送方式下?lián)Q輸管道和轉(zhuǎn)壓站所送水量小于等于該輸送方式所輸送的水量：

(2)非負約束：X≧0；

3.3.3優(yōu)先級：一級為各類水量目標；二級為控制費用目標；三級為各管道和轉(zhuǎn)壓站輸送水量日標；

4大型供水管網(wǎng)改、擴建工程優(yōu)化決策字典序線性目標規(guī)劃模型靈敏度分析

4.1某些問題的處理：

(1)考慮到四舍五入帶來的影響，為方便編程計算，特作如下規(guī)定：凡本模型靈敏度分析中間值的絕對值小于或等于0.00001的元素值均設(shè)為0.0。

(2)為方便編程計算：在本模型靈敏度分析中，用1.0e+35代替+∞，-(1.0e+35)代替-∞。

4.2單元素變化及多元素參數(shù)分析：

均針對第t個方案，在保持目前最優(yōu)基的條件下，求滿足下列不等式：B-1(b+ubt*)≧0的u的取值范圍(因篇幅所限，公式推導及結(jié)果省略).

4.3離散變化的靈敏度分析：

從單元素變化及多元素參數(shù)分析所確定的各取一u值，計算XTb*B-1(b+ubt*)及atk*=∑m+ns=1Wk,s*XTb,s*+∑mI=1uI,k*XTb,I*

其中：t—單元素變化或多元素參數(shù)分析的第t個方案；XTb,s*—問題變量中的負偏差變量值；XTb,I*—問題變量中的正偏差變量值。

5大型供水管網(wǎng)改擴建工程優(yōu)化決策字典序線性目標規(guī)劃模型多階段單純形算法程序框圖及靈敏度分析程序框圖[4]

5.1多階段單純形算法程序框圖

5.2靈敏度分析程序框圖

6例計算和分析

6.1實例

1.武漢武昌地區(qū)“湖改江”工程，含東湖水廠“湖改江”和團山水廠“湖改江”，根據(jù)武昌地區(qū)供水規(guī)劃，實現(xiàn)東湖水廠“湖改江”須補充水量12萬噸／日，實現(xiàn)團山水廠“湖改江”須增加水量18萬噸／日。

2.實現(xiàn)東湖水廠"湖改江"有以下六種可能輸送方式：

序號前輸轉(zhuǎn)輸改、擴建單位費用

(元／m3)

1南湖北路轉(zhuǎn)壓站南湖北路DN1200管600

2關(guān)山路轉(zhuǎn)壓站南湖北路DN1200管450

3關(guān)山路轉(zhuǎn)壓站南湖南路DN800管300

4關(guān)山路轉(zhuǎn)壓站關(guān)山一路延長線DNl000管300

5南湖北路DNl200管關(guān)山一路延長線DN1000管450

6湖邊泵站南湖北路DN1200600

(B)實現(xiàn)團山水廠“湖改江”有以下十種可能輸送方式：

序號前輸轉(zhuǎn)輸改、擴建單位費用

(元／m3)

1三環(huán)路DN400管三環(huán)路轉(zhuǎn)壓站600

2三環(huán)路DN400管天京路DN300管600

3三環(huán)路DN400管三環(huán)路DN100管900

4南湖北路轉(zhuǎn)壓站魯巷廣場600

5關(guān)山路轉(zhuǎn)壓站南湖南路300

6關(guān)山路轉(zhuǎn)壓站關(guān)山路延長線DN1000管300

7關(guān)山路轉(zhuǎn)壓站民院路DN400管600

8關(guān)山路轉(zhuǎn)壓站天京路ON300管900

9湖邊泵站環(huán)湖路DN300管450

10三環(huán)路轉(zhuǎn)壓站關(guān)山一路延長線DN1000管300

3.根據(jù)管網(wǎng)水力計算，各轉(zhuǎn)壓站及管道增輸水量(萬噸)情況如下：

南湖北路轉(zhuǎn)壓站：12；南湖北路DNl200管：12；湖邊泵站：3；三環(huán)路轉(zhuǎn)壓站：6；南湖南路DN800管；4.5；關(guān)山路轉(zhuǎn)壓站；7；關(guān)山一路延長線DNl000管；7：環(huán)湖路DN300；0.5；天泉路DN300：0.5；民院路DN400管：0.5；三環(huán)路DN400管：3；魯巷廣場DN800管：4。

4.根據(jù)控制費用：實現(xiàn)武昌地區(qū)"湖改江"管網(wǎng)及轉(zhuǎn)壓站改、擴建控制費用為：1.5億元。

要求選擇適宜的輸送方式并求各適宜輸送方式下各輸送管道及轉(zhuǎn)壓站的轉(zhuǎn)輸水量。

6.2本問題是一個典型的大型供水管網(wǎng)，改、擴建工程優(yōu)化決策問題。

它須滿足總水量及單水量目標，還須滿足控制費用目標，硬目標為是總水量目標，單水量及控制費用為軟目標，因此，該問題屬字典序線性目標規(guī)劃問題，故可用本文的字典序線性目標規(guī)劃模型求解程序求解，并進行模型的靈敏度分析。

6.2.1字典序線性目標規(guī)劃分階段算法程序求解結(jié)果及分析

(1)東湖水廠"湖改江"

序號前輸水量(萬噸)轉(zhuǎn)輸水量(萬噸)

1南湖北路轉(zhuǎn)壓站8南湖北路DN1200管8

2關(guān)山路轉(zhuǎn)壓站1南湖北路DN1200管1

3關(guān)山路轉(zhuǎn)壓站0南湖北路DN800管0

4關(guān)山路轉(zhuǎn)壓站0關(guān)山一路延長線DN1000管0

5南湖北路DN1200管0.5關(guān)山一路延長線DN1000管0.5

6湖邊泵站2.5南湖北路DN12002.5

(2)團山水廠"湖改江"

序號前輸水量(萬噸)轉(zhuǎn)輸水量(萬噸)

1三環(huán)路DN400管0三環(huán)路轉(zhuǎn)壓站0

2三環(huán)路DN400管0天京路DN300管0

3三環(huán)路DN400管1.5三環(huán)路DN100管1.5

4南湖北路轉(zhuǎn)壓站4.0魯巷廣場4.0

5關(guān)山路轉(zhuǎn)壓站4.5南湖南路4.5

6關(guān)山路轉(zhuǎn)壓站0.5關(guān)山路延長線DN1000管0.5

7關(guān)山路轉(zhuǎn)壓站0.5民院路DN400管0.5

8關(guān)山路轉(zhuǎn)壓站0.5天京路ON300管0.5

9湖邊泵站0.5環(huán)湖路DN300管0.5

10三環(huán)路轉(zhuǎn)壓站6關(guān)山一路延長線DN1000管6

由上表可知：東湖水廠"湖改江"中改、擴建單位費較低的兩種輸送方式被淘汰，而在團山水廠"湖改江"中改、擴建單位費較高的輸送方式被淘汰，這表明所確定的輸送方式并不存在一致性偏向問題。

6.2.2單元素參數(shù)靈敏度分析結(jié)果(因篇幅所限，多元素參數(shù)分析及離散變化分析結(jié)果省略).

在保持原模型實例所求最優(yōu)基仍為最優(yōu)的條件下，各輸送管道及轉(zhuǎn)壓站輸送水量允許變化范圍如下：南湖北路DN1200管，三環(huán)路轉(zhuǎn)壓站，天泉路DN300管：0—+∞；南湖北路轉(zhuǎn)壓站，湖邊泵站，三環(huán)路DN400管，南湖南路DN800管，關(guān)山路轉(zhuǎn)壓站，關(guān)山路延長線DN1000管，蚌潮路DN300管，民院路DN400管，魯巷廣場DN800管分別為：0—1.3，0—1.3，0—1，0—1，0—6，0—1，0—2.5，0—1，0—2.5。

參考文獻

[1]IGNIZIO,J.P“Goalprogramming：ApoolformultiobjectiveAnalysis”，JournalofoperationalResearch，Vol，29，Ⅱ，1978．

[2]J．P．Ignizio，閩仲求《單目標和多目標系統(tǒng)線性規(guī)劃》上海：同濟大學出版杜，1986．

[3]韓旭昆《線性覿劃靈敏度分析》,系統(tǒng)工程，1992．(1)．

[4]潭浩強，《C程序設(shè)計》，北京：清華大學出版社，1991.