油田開發(fā)規(guī)劃的合理編制

1、a題：油田開發(fā)規(guī)劃的合理編制 油田開發(fā)規(guī)劃的合理編制摘 要 本文討論了油田開發(fā)規(guī)劃的合理編制問題。首先對油田規(guī)劃期內(nèi)三項指標（產(chǎn)油量、產(chǎn)水量、費用）進行了預測，然后結合預測值建立多目標規(guī)劃模型，得出規(guī)劃期內(nèi)新增措施的具體實施方案，并對各種不同方案進行了評價。對于問題一，經(jīng)過對數(shù)據(jù)的統(tǒng)計、作圖分析，發(fā)現(xiàn)三項指標的歷史數(shù)據(jù)存在一定的規(guī)律，針對不同的變化規(guī)律，選擇不同的函數(shù)進行曲線擬合，可得到三項指標在“十一五”期間的預測值。為確定新增方案，建立多目標規(guī)劃模型。在滿足油田最大生產(chǎn)能力的前提下，以產(chǎn)油量上限、產(chǎn)水量下限、費用上限以及各項增產(chǎn)措施工作量上限作為約束條件，分別突出不同因素的影響，對目標函數(shù)

2、權重系數(shù)取不同值，得出不同方案采取各種增產(chǎn)措施的數(shù)量。對于問題二，針對問題一得出的不同方案，從產(chǎn)油量、產(chǎn)水量、費用三方面進行分析，指出不同方案的優(yōu)、缺點以及適用性，并給出合理化建議。但以上三種方案的各項增產(chǎn)措施均衡度不高，改變權重系數(shù)與正負偏差變量，并規(guī)定規(guī)劃期內(nèi)每年新增措施總數(shù)之差不超過30；每項新增措施在每的年數(shù)目之差不超過30。最終綜合考慮總產(chǎn)油量、總產(chǎn)水量、總費用、投入產(chǎn)出比、綜合含水率的影響，均衡方案一更具實用性，更具推廣性。均衡方案一如下表： 年份措施2006年2007年2008年2009年2010年總計110010086707042629666966666390總產(chǎn)油量（萬噸）1

3、521.51524.81415.81486.11396.07344.2總產(chǎn)水量（萬噸）5583.35789.66044.56400.76852.530670.8總費用（萬元）58948.495864.2161433.3196135.3213177.9725559.1注：詳見論文12-14頁。關鍵詞：曲線擬合；預測曲線；多目標規(guī)劃一、問題重述編制油田開發(fā)規(guī)劃是油田開發(fā)的核心問題，它是確定在一個時期內(nèi)，油田開發(fā)生產(chǎn)的戰(zhàn)略決策和具體部署，直接影響到油田的開發(fā)效果和開發(fā)效益的好壞，這就要求所編制的油田開發(fā)規(guī)劃要具有科學性、合理性和可行性。油田為了實現(xiàn)規(guī)劃期（五年）的各項規(guī)劃指標（主要包括產(chǎn)油量、綜合含

4、水及費用三項指標），首先要計算在規(guī)劃初期老井（規(guī)劃期之前實施的增產(chǎn)措施）的各項指標在規(guī)劃期的五年內(nèi)的預測值。與此同時，要想完成規(guī)劃指標，就要采取大量的增產(chǎn)措施。編制油田開發(fā)規(guī)劃方案，就是在滿足油田最大生產(chǎn)能力的前提下，制定出采取各種增產(chǎn)措施的數(shù)量，使得盡可能達到規(guī)劃指標。另外，在制定規(guī)劃方案時要盡可能做到均衡安排各項增產(chǎn)措施，實現(xiàn)科學開采、持續(xù)發(fā)展。二、問題分析2.1問題一的分析問題一應為一多目標規(guī)劃問題。問題的關鍵在于確定規(guī)劃期內(nèi)各種增產(chǎn)措施的數(shù)量，使得產(chǎn)油量、產(chǎn)水量及費用三項指標盡可能達到規(guī)劃指標。油區(qū)內(nèi)的井網(wǎng)分為四類：“七五”井、“八五”井、“九五”井及“十五”井。題目給出了各類井網(wǎng)近幾

5、年的產(chǎn)油量、產(chǎn)水量，經(jīng)過對數(shù)據(jù)進行統(tǒng)計、作圖、分析得知，產(chǎn)油量和產(chǎn)水量隨時間的變化存在一定的規(guī)律，因此，可根據(jù)數(shù)據(jù)的走勢圖，選擇適當?shù)暮瘮?shù)進行曲線擬合，預測出在規(guī)劃初期，老井的各項指標在規(guī)劃期的五年內(nèi)的預測值。產(chǎn)油量、產(chǎn)水量和費用均包含兩部分內(nèi)容：各類井網(wǎng)在規(guī)劃期內(nèi)的各項指標預測值和新增措施的實際值。在滿足各項約束條件的前提下，以產(chǎn)油量最大、費用最低、產(chǎn)水量最低為目標函數(shù)，建立多目標規(guī)劃模型，可求得各項增產(chǎn)措施的安排方案。2.2問題二的分析 對于多目標規(guī)劃問題，最終結果會因權重系數(shù)的不同而不同，由此可以突出不同因素的影響，得到不同方案。針對不同方案，根據(jù)突出因素的不同，給出優(yōu)、缺點及實用性。三

6、、模型假設（1）假設歷史數(shù)據(jù)具有較高的準確性； （2）假設不同類油井的產(chǎn)油量、產(chǎn)水量隨時間變化趨勢相同；（3）假設“十一五”期間不會發(fā)生重大自然災害影響到各項指標。四、符號說明第年采取第項措施的井口數(shù)（；，以下相同）第年的產(chǎn)油量（萬噸）四類井網(wǎng)在第年產(chǎn)油量預測值總和（萬噸）第年的產(chǎn)水量（萬噸）四類井網(wǎng)在第年產(chǎn)水量總和（萬噸）第年的總費用（萬元）四類井網(wǎng)第年總費用預測值（萬元）第年第項措施的單井產(chǎn)油量（萬噸） 第年第項措施的單井產(chǎn)水量（萬噸） 第年第項措施的單井年費用（萬元） 五、模型建立及求解5.1問題一5.1.1模型準備5.1.1.1產(chǎn)油量的預測題目中給出了該油區(qū)各類井網(wǎng)在近幾年的產(chǎn)油量、產(chǎn)

7、水量，但“十一五”期間的產(chǎn)油量、產(chǎn)水量卻未知，通過作圖可得到近幾年產(chǎn)油量、產(chǎn)水量隨時間變化的走勢圖（如圖1）。由圖1并結合實際情況分析可知，產(chǎn)油量隨時間變化的趨勢應為先增長后降低的過程，但數(shù)據(jù)量過少，因此只能選則有限點進行曲線擬合，得到預測結果。在求解“七五”井、“八五”井、“九五”井在規(guī)劃期內(nèi)的預測值時，根據(jù)走勢，“七五”、“八五”井選擇反比例函數(shù)進行擬合，“九五”井只選取2001-2005年的五個數(shù)據(jù)點進行一次函數(shù)擬合。“十五”井則選取二次函數(shù)進行擬合。圖1 各類井產(chǎn)量隨時間變化走勢圖由此得到四類井網(wǎng)產(chǎn)油量隨時間變化的函數(shù)方程式如下：“七五”井： “八五”井： “九五”井： “十五”井：

8、由此可以求得四類井網(wǎng)在“十一五”期間每年產(chǎn)油量的預測值：表1 四類井網(wǎng)在“十一五”期間每年產(chǎn)油量預測值年份產(chǎn)量（萬噸）2006年2007年2008年2009年2010年“七五”井335.4320.9309.1299.1290.6“八五”井219.3208.3199.2191.7185.4“九五”井414.5380.4346.3312.2278.1“十五”井514.45512.3479.0414.5318.8總計1483.71421.91333.61217.51072.95.1.1.2產(chǎn)水量的預測對于產(chǎn)水量的預測，首先得到了產(chǎn)水量隨時間變化的走勢圖（圖2）。圖2 產(chǎn)水量隨時間變化的走勢圖從圖上走

9、勢可知，“七五”井產(chǎn)水量隨時間逐年遞減，選擇指數(shù)函數(shù)進行擬合?！鞍宋濉边x擇二次數(shù)函數(shù)進行擬合，“九五”井選擇泊松公式進行擬合，“十五”井產(chǎn)水量隨時間逐年遞增，選擇一次函數(shù)進行擬合。由此得到四類井網(wǎng)產(chǎn)水量隨時間變化的表達式：“七五”井： “八五”井： “九五”井： “十五”井： 由此可以求得四類井網(wǎng)在“十一五”期間每年產(chǎn)水量的預測值：表2 四類井網(wǎng)在“十一五”期間每年產(chǎn)水量預測值年份產(chǎn)量（萬噸）井網(wǎng)2006年2007年2008年2009年2010年“七五”井2263.9213420201919.91832.2“八五”井661.3633.4590.5532.5459.4“九五”井1538.2155

10、8.51613.31615.31597.6續(xù)表2“十五”井10851292.51499.91707.41914.8總計5548.45648.45723.75775.158045.1.1.3對產(chǎn)油量、產(chǎn)水量、費用的分析由題目知，在規(guī)劃期內(nèi)，每種增產(chǎn)措施的單井年產(chǎn)油量、產(chǎn)水量及費用是隨著年份的不同而變換的， 因此，每年三項指標的構成均不同，但三項指標的組成形式是類似的。以產(chǎn)油量為例，2006年的產(chǎn)油量為該年各類井網(wǎng)的總產(chǎn)油量與該年新增措施總產(chǎn)油量之和，2007年的產(chǎn)油量為該年各類井網(wǎng)的總產(chǎn)油量與該年新增措施總產(chǎn)油量之和，再加上2006年的新增措施在2007年的總產(chǎn)油量，之后依此類推，即可得到每年產(chǎn)

11、油量的目標函數(shù)表達式。同理可得出產(chǎn)水量、費用目標函數(shù)表達式。5.1.2多目標規(guī)劃模型的建立在得到四類井網(wǎng)在“十一五”期間的產(chǎn)油量、產(chǎn)水量預測值之后，就可以得到每年產(chǎn)油量、產(chǎn)水量、費用目標函數(shù)：產(chǎn)油量： 產(chǎn)水量： 費用：根據(jù)題中要求，可以得到如下約束條件：1） 規(guī)劃期內(nèi)增產(chǎn)措施的工作量上限： 2）規(guī)劃期內(nèi)增產(chǎn)措施的產(chǎn)油量下限： 3）規(guī)劃期內(nèi)增產(chǎn)措施的產(chǎn)水量上限：4）規(guī)劃期內(nèi)增產(chǎn)措施的費用上限： 5.1.3多目標規(guī)劃模型的求解本問題為一多目標規(guī)劃問題，問題中要考慮具有多個相互矛盾的目標的優(yōu)化。因為愿望與實際結果之間有一定的差異（稱為目標偏差），一般問題都要計算目標偏差，因此不管是哪種目標表達式，都

12、應加上一個負偏差量（），并減去一個正偏差量（），使它們轉(zhuǎn)換成目標規(guī)劃的格式。然后建立能反映和衡量各種解的目標達成程度的關系式（達成函數(shù)），為了滿足目標（或約束），要求有關偏差量取極小，在達成函數(shù)中反映這種要求。對解的評價，一般采用使加權的目標偏差之和達到極小的方法。因此要對約束條件1)、2)、3)、4）的每一個表達式，加上正負偏差變量（），并且得到如下達成函數(shù)：其中； ；.其中、為權重系數(shù)，權重系數(shù)越高，表明約束條件越強、優(yōu)先級越高。1）方案一 當權重系數(shù)取，時，即規(guī)劃期內(nèi)新增措施的工作量上限為硬約束，優(yōu)先級最高，所以使值很大，表示產(chǎn)油量為第二約束條件，將產(chǎn)水量、費用作為第三約束條件，突出產(chǎn)油

13、量的影響，追求產(chǎn)油量最大。由lingo軟件進行求解得到：表3 方案一各項新增措施數(shù)量 年份措施2006年2007年2008年2009年2010年總計137516344252119271390319319839430續(xù)表343743667405410146480545347476476總計68747610724393353009由此計算出方案一2006-2010年各年三項指標的總量：表4 2006-2010年各年三項指標值年份措施2006年2007年2008年2009年2010年總計總產(chǎn)油量（萬噸）下限預測下限預測下限預測下限預測下限預測下限預測15001535.514901552.314851

14、465.314501504.514231423.173487480.9總產(chǎn)水量（萬噸）上限預測上限預測上限預測上限預測上限預測上限預測57065358.758435459.861525610.963135853.067206214.33073432026.0 總費用（萬元）上限預測上限預測上限預測上限預測上限預測上限預測8658086563.7119100101623.017198.0171933.4196132195993.6199321199285.5773113755399.22）方案二當權重系數(shù)取，時，即追求產(chǎn)水量最小。由lingo軟件進行求解得到：表5 方案二各項新增措施數(shù)量 年份措

15、施2006年2007年2008年2009年2010年總計142542523183183392374294363363515429444682452534744531476總計5228977864663182989由此計算出方案二2006-2010年各年三項指標的總量：表6 2006-2010年各年三項指標值年 份措施2006年2007年2008年2009年2010年總計總產(chǎn)油量（萬噸）下限預測下限預測下限預測下限預測下限預測下限預測15001523.714901567.314851458.014501502.714231423.073487474.8總產(chǎn)水量（萬噸）上限預測上限預測上限預測上限

16、預測上限預測上限預測57065348.058435454.261525606.663135802.967206178.83073430514.3 總費用（萬元）上限預測上限預測上限預測上限預測上限預測上限預測8658075200.8119100108359.8171980169751.1196132192541.0199321198827.5773113744680.23）方案三：當權重系數(shù)取，時，即追求費用最低。由lingo軟件進行求解得到：表7 方案三各項新增措施數(shù)量 年份措施2006年2007年2008年2009年2010年總計139035142621421572883803430430

17、4374374544544565345347106370476總計9245857665022883065由此計算出方案三2006-2010年各年三項指標的總量：表8 2006-2010年各年三項指標值年 份措施2006年2007年2008年2009年2010年總計總產(chǎn)油量（萬噸）下限預測下限預測下限預測下限預測下限預測下限預測15001546.114901586.514851468.114501529.414231423.073487553.3總產(chǎn)水量（萬噸）上限預測上限預測上限預測上限預測上限預測上限預測57065364.958435486.261525633.163135906.16720

18、6279.03073430773.2 總費用（萬元）上限預測上限預測上限預測上限預測上限預測上限預測8658086524.6119100117120.4171980169377.3196132196110.5199321199112.6773113768245.55.2問題二5.2.1方案一的評價方案一針對產(chǎn)油量因素比較突出所制定的方案,此因素突出是由油價上漲,儲油豐富等原因所造成的,優(yōu)點是針對性很強, 但作用時間很短,具有短期性的特點,而且需要根據(jù)情況做靈活調(diào)整,適用于油價上漲,儲油豐富等產(chǎn)油效益較突出的情況,有一定局限性。5.2.2方案二的評價方案二是針對產(chǎn)水量因素比較突出所制定的方案,此

19、因素突出是由產(chǎn)水突出造成的,優(yōu)點是針對性很強,但作用時間短,需要根據(jù)情況做靈活調(diào)整,適用于水不好處理,水量多、大的情況,適用度不廣。5.2.3方案三的評價方案三針對資金因素較突出時,選用這個方案,具有很好的效果,能維持油田繼續(xù)運行而不會造成因資金短缺而停產(chǎn)的嚴重后果,但這只是短期調(diào)整,當情況變化時要及時調(diào)整,否則效益會受到很大影響。5.2.4三種方案的比較設投入產(chǎn)出比為費用與產(chǎn)油量的比值（表示單位費用的產(chǎn)油量）綜合含水率為產(chǎn)水量與產(chǎn)水、產(chǎn)油總量之和的比值，得到表9：表9 三種方案的投入產(chǎn)出比與綜合含水率 方案指標方案一方案二方案三投入產(chǎn)出比（萬噸/萬元）101.8897.4297.42綜合含水

20、率79.1%79.07%79.15%從實際情況出發(fā)，油田應追求投入產(chǎn)出比最大，綜合含水率最低的方案，綜合兩因素分析方案一相對較優(yōu)。5.2.5均衡方案 以上三種方案是簡單的線性規(guī)劃問題的解，分別突出了三項指標的影響，從結果看，三種方案的均衡度并不高，有的年份甚至沒有實施新增措施，因此需要考慮均衡方案。為達使新增措施均勻分布在規(guī)劃方案中，增加另外兩約束條件：規(guī)劃期內(nèi)每年新增措施總數(shù)之差不超過30；每項新增措施在每年數(shù)目之差不超過30。由此得到均衡后的方案（如表10）：1）對方案一的均衡：表10 方案一均衡后各項新增措施數(shù)量 年份措施2006年2007年2008年2009年2010年總計110010

21、086707042629666966666390398 989868684304878687575737457777771071074456929310512212253477710778107107476總計6276276275975973075表11 方案一均衡后2006-2010年各年三項指標值年 份措施2006年2007年2008年2009年2010年總計總產(chǎn)油量（萬噸）下限預測下限預測下限預測下限預測下限預測下限預測15001521.514901524.814151415.814501486.114231396.073487344.2總產(chǎn)水量（萬噸）上限預測上限預測上限預測上限預測上

22、限預測上限預測57065583.358435789.661526044.563136400.767206852.53073430670.8 總費用（萬元）上限預測上限預測上限預測上限預測上限預測上限預測8658058948.411910095864.2171980161433.3196132196135.3199321213177.9773113725559.12） 對方案二的均衡：表12 方案二均衡后各項新增措施數(shù)量年份措施2006年2007年2008年2009年2010年總計1989795686842629090906060390398989868684304926592636237456

23、687699596413694989412412453471415144444131總計5525505525225222698表13 方案二均衡后2006-2010年各年三項指標值年 份措施2006年2007年2008年2009年2010年總計總產(chǎn)油量（萬噸）下限預測下限預測下限預測下限預測下限預測下限預測15001518.014901515.714851398.314501462.614231367.773487262.5總產(chǎn)水量（萬噸）上限預測上限預測上限預測上限預測上限預測上限預測57065578.558435770.661526004.163136312.567206707.33073

24、430373.1 總費用（萬元）上限預測上限預測上限預測上限預測上限預測上限預測8658054673.411910088026.0171980152095.9196132184380.2199321199229.9773113678405.53） 對方案三的均衡：表14 方案三均衡后各項新增措施數(shù)量 年份措施2006年2007年2008年2009年2010年總計1979797676742529090906060390310186101717143049363926363374續(xù)表14562 926292924006921059312212253471719174747147總計55255255

25、25225222700表15 方案三均衡后2006-2010年各年三項指標值年 份措施2006年2007年2008年2009年2010年總計總產(chǎn)油量（萬噸）下限預測下限預測下限預測下限預測下限預測下限預測15001521.514901524.814851415.814501486.114231396.073487262.3總產(chǎn)水量（萬噸）上限預測上限預測上限預測上限預測上限預測上限預測57065348.258435442.361525569.563135781.9672060728 總費用（萬元）上限預測上限預測上限預測上限預測上限預測上限預測8658058948.4

26、11910095864.2171980161433.3196132196135.3199321213177.9773113679705.7在得到新的均衡方案后與前三種方案進行對比：表16 均衡前后的對比 方案指標方案一均衡一方案二均衡二方案三均衡三總產(chǎn)油量（萬噸）7437.57344.27474.87262.57553.37262.3總產(chǎn)水量（萬噸）32026.030670.830514.330373.130773.230375.8總費用（萬元）724572.2725559.1744680.2678405.5768245.5679705.7投入產(chǎn)出比100.9798.7999.6293.241

27、01.7098.79綜合含水率79.20%79.34%79.15%79.37%79.14%79.34%由上述表中數(shù)據(jù)可以看到，在對三種方案進行均衡后，新方案的各新增措施的數(shù)量在每年分布比較均勻。跟前三種方案比，雖然總產(chǎn)油量有所下降，但產(chǎn)水量、費用也有所下降，因此，從總體上看，均衡方案一更具推廣性。六、模型改進與推廣對于產(chǎn)油量、產(chǎn)水量的預測，為保證精度，對不同的走勢圖，選擇了一次函數(shù)、二次函數(shù)、反比例函數(shù)分別進行曲線擬合，沒能選擇統(tǒng)一的函數(shù)進行預測，實際操作比較繁瑣，要根據(jù)數(shù)據(jù)的走勢來判斷函數(shù)類別進行擬合，可推廣度不高。因此可選擇統(tǒng)一函數(shù)，如指數(shù)函數(shù)、對數(shù)函數(shù)進行擬合，提高計算效率。對于多目標規(guī)

28、劃模型，每一種方案只突出了一種指標，綜合考慮不足，可改變權重系數(shù)，得到新的規(guī)劃方案，對比已經(jīng)求得的方案，選出綜合度最高的方案。本問題的預測模型可應用于城市供水量預測、生物種群數(shù)量的預測、商品銷售額預測等預測問題中。對于多目標規(guī)劃模型，可應用于管理、工業(yè)生產(chǎn)等實際問題中。七、模型評價本文主要有以下優(yōu)點：（1）預測模型的數(shù)據(jù)按照不同的走勢分開來處理，具體問題具體分析，這樣可以減少個別數(shù)據(jù)發(fā)生劇烈變動，方便我們把握規(guī)律；（2）多目標規(guī)劃模型，著重考慮了產(chǎn)油量、產(chǎn)水量、費用三個指標，使問題得到簡化。本文主要有以下缺點：（1）由于歷史數(shù)據(jù)有限，我們在預測的過程中，預測精度受到了限制；（2）多目標規(guī)劃的解

29、應該有很多，所求解并不一定是最優(yōu)解，要進行多次取權重并對不同結果進行分析，得到相對較優(yōu)值，因此加大了工作量。參考文獻1 薛定宇等，高等應用數(shù)學問題的matlab求解，北京：清華大學出版社，2004年5月第一版。2 陳杰，matlab寶典，北京：電子工業(yè)出版社，2007年。3 袁新生，邵大宏，lingo和excel在數(shù)學建模中的應用，北京：科學出版社，2007年1月第一版。4 謝金星，薛毅，優(yōu)化建模與lindo/lingo軟件，北京：清華大學出版社，2005年7月第一版。附 錄附錄1：反比例函數(shù)擬合程序function f=sq(a,t7)f=a(1)./(t7+a(2)+a(3);x7=500

30、.6 442.4 428.6 370.1 343.1;t7=2001:2005;a0=1,1,1;a=lsqcurvefit(sq,a0,t7,x7);z=sq(a,t7);plot(t7,x7,*,t7,z)xlabel(年份)ylabel(產(chǎn)油量（萬噸）)yuce7=sq(a,2006:2010)附錄2：二次函數(shù)擬合求解預測值程序function f=nihe1(c2,t2)f=-15.6.*t2.2+c2(1).*t2+c2(2);t1=1997:2005;x9=197.4 297.7 412.8 547.0 579.8 547.5 527.0 492.3 437.0;c1=polyfi

31、t(t1,x9,2);yc9=polyval(c1,t1);plot(t1,x9,-*,t1,yc9,-+); t2=2001:2005;x10= 72.3 218.2 297.1 416.1 508.7;c20=c1(2),c1(3);c2=lsqcurvefit(nihe1,c20,t2,x10);yuce=-15.6.*t2.2+c2(1).*t2+c2(2);plot(t2,x10,-+,t2,yuce,-*) t=2001:2010;hunhe10=-15.6.*t.2+c2(1).*t+c2(2);plot(t2,x10,-+,t,hunhe10,-*)附錄3：多目標模型求解程序及

32、結果model:sets:nf/1,2.5/:y,s,f,yx,ss,fs,q,w,m,k;cs/1,2.7/:js;links(cs,nf):x,yk,sk,fk; hh/1,2,3/:h;endsetsdata:y=1483.7 1421.9 1333.6 1217.5 1072.9;s=5548.4 5648.4 5723.7 5775.1 5804;f=16861 20477 60851 80114 85014;yx=1500 1490 1485 1450 1423;ss=5706 5843 6152 6313 6720;fs=86580 119100 171980 196132 19

33、9321;js=426 390 430 374 445 534 476;yk=0.0820 0.2160 0.2679 0.23040.19810.0720 0.0960 0.1479 0.10040.08810.0650 0.0860 0.1670 0.13440.09760.0540 0.0862 0.1509 0.12280.10060.0380 0.0760 0.0979 0.07040.05810.0570 0.0820 0.1009 0.0828 0.06860.0450 0.0750 0.1074 0.08930.0551;sk=0.0720 0.2120 0.3750 0.56

34、660.75850.0630 0.1880 0.2050 0.40860.67850.0570 0.1280 0.3120 0.57660.78460.0500 0.1760 0.3788 0.49860.64380.0320 0.1790 0.2131 0.42530.66530.0442 0.1262 0.2045 0.35840.47440.0675 0.1775 0.3196 0.56350.7874;fk=128.1000 107.0000 72.2500 56.5500 44.700078.1000 60.0000 52.2500 39.5500 30.700064.1000 51

35、.0000 42.2500 33.5500 21.700043.0000 31.0500 22.2000 17.0500 13.150058.1000 49.0000 37.2500 30.5500 24.700036.9000 30.7500 23.7000 17.3000 12.000050.5000 43.8000 35.7000 27.3500 22.0500;enddata!每年產(chǎn)油總量;q(1)=y(1)+sum(cs(i):yk(i,1)*x(i,1);q(2)=y(2)+sum(cs(i):(yk(i,1)*x(i,2)+yk(i,2)*x(i,1);q(3)=y(4)+sum

36、(cs(i):(yk(i,1)*x(i,3)+yk(i,2)*x(i,2)+yk(i,3)*x(i,1);q(4)=y(4)+sum(cs(i):(yk(i,1)*x(i,4)+yk(i,2)*x(i,3)+yk(i,3)*x(i,2)+yk(i,4)*x(i,1);q(5)=y(5)+sum(cs(i):(yk(i,1)*x(i,5)+yk(i,2)*x(i,4)+yk(i,3)*x(i,3)+yk(i,4)*x(i,2)+yk(i,5)*x(i,1);!每年產(chǎn)水總量;w(1)=s(1)+sum(cs(i):(sk(i,1)*x(i,1);w(2)=s(2)+sum(cs(i):(sk(i,

37、1)*x(i,2)+sk(i,2)*x(i,1);w(3)=s(3)+sum(cs(i):(sk(i,1)*x(i,3)+sk(i,2)*x(i,2)+sk(i,3)*x(i,1);w(4)=s(4)+sum(cs(i):(sk(i,1)*x(i,4)+sk(i,2)*x(i,3)+sk(i,3)*x(i,2)+sk(i,4)*x(i,1);w(5)=s(5)+sum(cs(i):(sk(i,1)*x(i,5)+sk(i,2)*x(i,4)+sk(i,3)*x(i,3)+sk(i,4)*x(i,2)+sk(i,5)*x(i,1);!每年總費用;m(1)=f(1)+sum(cs(i):(fk(i

38、,1)*x(i,1);m(2)=f(2)+sum(cs(i):(fk(i,1)*x(i,2)+fk(i,2)*x(i,1);m(3)=f(3)+sum(cs(i):(fk(i,1)*x(i,3)+fk(i,2)*x(i,2)+fk(i,3)*x(i,1);m(4)=f(4)+sum(cs(i):(fk(i,1)*x(i,4)+fk(i,2)*x(i,3)+fk(i,3)*x(i,2)+fk(i,4)*x(i,1);m(5)=f(5)+sum(cs(i):(fk(i,1)*x(i,5)+fk(i,2)*x(i,4)+fk(i,3)*x(i,3)+fk(i,4)*x(i,2)+fk(i,5)*x(

39、i,1);!約束條件一：規(guī)劃期內(nèi)增產(chǎn)措施的工作量上限;sum(nf(i):x(1,i)+d1_-d1=426;sum(nf(i):x(2,i)+d2_-d2=390;sum(nf(i):x(3,i)+d3_-d3=430;sum(nf(i):x(4,i)+d4_-d4=374;sum(nf(i):x(5,i)+d5_-d5=445;sum(nf(i):x(6,i)+d6_-d6=534;sum(nf(i):x(7,i)+d7_-d7=476;!約束條件二：規(guī)劃期內(nèi)增產(chǎn)措施的產(chǎn)油量下限;q(1)+d8_-d8=1500;q(2)+d9_-d9=1490;q(3)+d10_-d10=1458;q(

40、4)+d11_-d11=1450;q(5)+d12_-d12=1423;!約束條件三：規(guī)劃期內(nèi)增產(chǎn)措施的產(chǎn)水量上限;w(1)+d13_-d13=5760;w(2)+d14_-d14=5843;w(3)+d15_-d15=6152;w(4)+d16_-d16=6313;w(5)+d17_-d17=6720;!約束條件四：規(guī)劃期內(nèi)增產(chǎn)措施的費用上限;m(1)+d18_-d18=86580;m(2)+d19_-d19=119100;m(3)+d20_-d20=171980;m(4)+d21_-d21=196132;m(5)+d22_-d22=199321;!規(guī)劃期內(nèi)總產(chǎn)油量;h(1)=sum(nf(

41、i):q(i);!規(guī)劃期內(nèi)總產(chǎn)水量;h(2)=sum(nf(i):w(i);!規(guī)劃期內(nèi)總費用;h(3)=sum(nf(i):m(i);for(links(i,j):gin(x(i,j);min=100000*(d1+d2+d3+d4+d5+d6+d7)+1000*(d8_+d9_+d10_+d11_+d12_)+(d13+d14+d15+d16+d17)+(d18+d19+d20+d21+d22);end加入均衡條件后的程序段：!規(guī)劃期內(nèi)每年新增措施總數(shù);k(1)=sum(cs(i):x(i,1);k(2)=sum(cs(i):x(i,2);k(3)=sum(cs(i):x(i,3);k(4)=sum(cs(i):x