某地區(qū)未來十年電力發(fā)展規(guī)劃模型

某地區(qū)未來十年電力發(fā)展規(guī)劃模型【摘要】本文建立了一個電力規(guī)劃中經(jīng)濟最優(yōu)化靜態(tài)模型。本文認為該問題中的經(jīng)濟最優(yōu)化不是經(jīng)濟利益達到最大化，而是在滿足技術要求下的成本最小化。該模型首先利用數(shù)學軟件LINDO解出滿足技術要求下的各方案的分配份額：擴建舊火電站:新建水電站:新建火電站=5:4:2。其次，本文考慮了各項建設的時間、順序上的安排對于節(jié)省成本的影響，用MATLAB軟件通過畫圖確定了在第一年年初即隊原舊火電站進行擴建；在第四年投入使用新的水電站；最后在第八年年初起投入使用新的火電站。再次，本文考慮了從裝機臺數(shù)的角度降低成本，以期在最小的成本投入下完成各項建設。我們最后得出的建設方案安排為：在第一階段，也即擴建原有的舊火電站階段，我們預備在第一年擴建兩臺單機容量為10千瓦的發(fā)電機組；而在第二年和第三年各擴建一臺和兩臺；第二階段，也即新建水電站階段，從第四年初到第七年末的四年間，每年建設一臺單機容量為25千瓦的發(fā)電機組臺數(shù)；第三階段，也即新建火電站階段，在第八年和第九年初，各建設并投入使用一臺單機容量為30千瓦的發(fā)電機組；至此，全部建設完成并投入正常運作。由此，我們算出最終的投入成本為23.3866億。【關鍵詞】經(jīng)濟最優(yōu)化、成本最小化、靜態(tài)模型、MATLAB問題的重述某地區(qū)在制定十年電力發(fā)展規(guī)劃時遇到一個問題：根據(jù)電力需求預測得知，該地區(qū)在十年后發(fā)電裝機容量需要增加180萬千瓦，那時的年發(fā)電量需要增加100億度。根據(jù)調(diào)查和討論，電力規(guī)劃的備選技術方案有三個：擴建原有的火電站，但最多只能再安裝五臺10萬千瓦的發(fā)電機組；新建水電站，但最多只能安裝四臺25萬千瓦的發(fā)電機組；再新建一個火電站，最多只能安裝四臺30萬千瓦的發(fā)電機組。通過調(diào)研和計算，獲得有關的參數(shù)如表1所示：表1:備選技術方案參數(shù)表

工程投資單機容量

（萬千瓦）允許裝機臺數(shù)資本回收因子年運行成本（百萬元/億度）負荷因子備選方案工程特點前期工程投資（百萬元）單機設備投資（百萬元）1擴建舊火電站

211050．1034．110．662新建水電站504702540．05782．280．43新建火電站240653040．1033．650．7表一表中負荷因子為全年滿功率運行天數(shù)與全年總天數(shù)之比。根據(jù)該地實際調(diào)查原有火電站平均全年滿功率運行天數(shù)為241天，水電站和新建的火電站應分別為146天和255天，而全年365天，故折算得表中數(shù)據(jù)。表中資本回收因子是由如下數(shù)據(jù)所確定的，火電站的回收年限取15年，年利息率為0．06；水電站的回收年限取30年，年利息率為0．04，即得表中所列數(shù)值。我們的任務是：在滿足上述技術要求的前提下，選取經(jīng)濟效果最優(yōu)的建設方案.即：原來的火電站應如何擴建？新建的火電站應如何確定單機容量為25萬千瓦的發(fā)電機組的數(shù)量？新建的火力發(fā)電站應如何確定單機容量為30萬千瓦的發(fā)電機組的數(shù)量？年發(fā)電量則取決于火電站和水電站的裝機臺數(shù)、單機容量和負荷因子。有如下關系式：新的水電站和火電站的投資都有個前期工程投資，但是因為是個固定的常數(shù)，所以它所引起的等額償還資金也是個常數(shù)，我們這里在考慮成本最小化時可以忽略它的影響。所以，根據(jù)題目的技術要求，我們可以將目標函數(shù)表示如下：但是我們是有所約束的，為非限制性約束：用LINDO解得各個相應的值為：具體實現(xiàn)程序見附錄一?？紤]該地區(qū)用電需求量以遞增的趨勢增長的模型4.2.1確定用電年增長量和年需求增長率因為實際生活中，人們的生產(chǎn)生活對用電需求量會因為經(jīng)濟水平的提高，恩格爾系數(shù)的降低而不斷的增長，且我們可預見，隨著時間的不斷后延，我們的科技水平也在不斷的進步，對用電的需求增量也會不斷增大，在短時期內(nèi)，用電增長量應該滿足一個規(guī)律，故大膽假設：這個用電需求量的增長曲線應該是一個指數(shù)增長型的，有如下公式：將公式作如下變形：其中，代表現(xiàn)階段用電需求量，代表原來的電力需求量，代表年平均需求增長率，代表年限，即第幾年。我們利用網(wǎng)絡資源找到了一組全國的從1981年至2000年的歷年電力裝機和發(fā)電量的構成比的數(shù)據(jù)，如下表：

歷年電力裝機和發(fā)電量的構成比（1981～2000）

年份裝機容量（萬千瓦）比重（%）發(fā)電量（億千瓦時）比重（%）水電火電水電火電水電火電水電火電19812193472031.768.3656243721.278.819822296494031.768.3744253322.777.319832416522831.668.4864265124.675.419842560545232.068.0868290223.077.019852641606430.369.7924318322.577.519862754662829.470.6945355121.079.019873019727129.370.71002397120.179.919883270828028.371.71092435920.080.019893458920627.372.71185466220.379.7199036051018426.173.91263495020.379.7199137881135925.075.01248552718.481.6199240681258524.475.61315622717.482.6199344891380224.575.51516686818.181.9199449061487424.574.41668747018.080.5199552181629424.075.01868807418.680.2199655581788623.575.61869878117.381.3199759731924123.575.61946925217.281.6199865072098823.575.72043938817.681.1199972972234324.474.821291004717.381.5200079352375424.974.424311107917.881.0表二我們假定上面這組數(shù)據(jù)中的電的需求和供應量始終是平衡的，這樣我們就可以從標了標記的兩組數(shù)據(jù)的和的變化去求得這個年增長率。在對數(shù)據(jù)進行了前期取對數(shù)之后，我們利用MATLAB軟件對數(shù)據(jù)進行了線性擬合，得到我們所需要的數(shù)據(jù)。并用MATLAB對這20個數(shù)據(jù)進行了數(shù)據(jù)線性擬合（具體實現(xiàn)過程見附錄二）。圖一如上圖所示，數(shù)據(jù)的線性情況較好，故我們認為假設合理，于是我們有：由于現(xiàn)階段電量的供求剛好達到平衡，且我們的目標函數(shù)是求所花成本投入最小，顯而易見，擴建舊火電站是必須先做的，這樣才有可能達到成本最小值。所以我們接下來要解決的就是新的水電站和新的火電站的建造順序和建造時間問題，決定水電站和火電站的建造時間和順序，以求得在最小的最終建設成本條件下完成建設。4.2.2解出水電站和火電站的建造順序和建造時間正如我們的模型的分析中所述，每臺發(fā)電機投入生產(chǎn)之后，其每年都在盈利，經(jīng)計算比較，其金額比發(fā)電站的年支出總資金大一個數(shù)量級，故直接考慮順差的差值大小來決定先建造哪一個發(fā)電站是不合理的，我們對于建造順序的安排也取決于如何使所花成本最小。由于每年的支付金額都是等額的，而且建成發(fā)電站之后全部在負荷因子下工作，所以各個發(fā)電站的年運行成本也是一個固定的數(shù)值，所以我們可以考慮用一年的支出總金額來代替十年支出總金額。在一年的費用支出中，主要有兩方面的因素：每年運行時的年運行成本；購買機器時的工程投入所引起的等額償還資金。在上一模型中，我們已經(jīng)知道某一發(fā)電站的年運行總成本是與該發(fā)電站的年發(fā)電量成正比關系的，年發(fā)電量又取決于火電站和水電站的裝機臺數(shù)、單機容量和負荷因子。工程投資所引起的等額償還資金也可由上一模型中的公式可以算得。故水電站或火電站的年費用支出資金應為：經(jīng)計算，水電站的年費用支出資金要小于火電站的費用支出，所以我們判定三個發(fā)電站的最終建造順序方案為：先擴建舊火電站，再新建水電站，最后新建火電站。在模型中，我們可以預見用戶對于用電的年需求量是一個指數(shù)上升的曲線，而電的供應量是一個分段的函數(shù)模式。根據(jù)最小成本原則，我們可以解出在什么年份時原有的電力供應會出現(xiàn)供不應求的現(xiàn)象，然后在這一年到來之際，增加投入建新的發(fā)電站。顯然，這樣的投資將會使成本最小，故我們可依據(jù)以下平衡關系式解出相應的年份：再由上面模型的一些公式和結論，解得：即在當年擴建舊火電站之后，在第3年開始建造并使用新建的水電站，第七年建造使用新的火電站。利用MATLAB作圖如下（具體程序?qū)崿F(xiàn)見附錄三）：圖二考慮裝機臺數(shù)的時間模型在上圖中，我們看到在需求曲線上方的分段函數(shù)中，有部分電是處于浪費狀態(tài)的，這必將使投入的資金不能產(chǎn)生最大的經(jīng)濟效益，也間接的體現(xiàn)出成本還未達到最小，還有更為節(jié)省的投資方案。我們是否可以考慮將機器的投入生產(chǎn)采取分段投入的方案呢？這個方案就是使電供應的分段函數(shù)更為靠近需求曲線，也即兩個函數(shù)與坐標軸圍住的面積達到最小，我們提出的解決方案如下：算出每年的需求量，在分段函數(shù)上取盡可能小的函數(shù)曲線與之對應（程序見附錄四），列表如下：年份12345年需求電量87.587294.9445102.9199111.5651120.9366年份678910年需求電量131.0953142.1073154.0443166.9840181.0107表三我們考慮在滿足用電年需求量的技術前提下，使供應量達到最小，使成本最小化。由上面的圖4與表3的數(shù)據(jù)及其題目給我們的各發(fā)電站的單機年發(fā)電量，我們可以進而求得在年時滿足需求量下的最小裝機臺數(shù)（程序見附錄五），列表如下：年份12345年需求電量87.587294.9445102.9199111.5651120.9366年裝機臺數(shù)年初完成21211年發(fā)電量92.36898.152109.72118.48127.24年份678910年需求電量131.0953142.1073154.0443166.9840181.0107年裝機臺數(shù)年初完成11110年發(fā)電量136144.76163.12181.48181.48表四現(xiàn)利用表四用MATLAB作出其圖形如下：圖三求出最終的成本由于每年的發(fā)電機投入都是不同的，故應該計算十年的總資本投入。因為每臺發(fā)電機的投入資金量對另外的發(fā)電機沒有影響，所以可以獨立地進行計算?？偝杀居蓛刹糠謽嫵桑汗こ掏度雰斶€資金和年運行成本投入。經(jīng)計算，最終成本為23.386634億元。五、模型的優(yōu)缺點及其改進方向模型的優(yōu)點：本模型在考慮最經(jīng)濟的方案時，并沒有直接考慮其收益的最大值，而是取其成本的最小值；本模型將復雜的動態(tài)規(guī)劃問題通過合理的假設簡化成靜態(tài)的規(guī)劃問題；本模型對問題以點連線，步步深入。但并沒有什么晦澀難懂的地方，易于理解；模型的缺點：本模型在考慮發(fā)電能力時，沒有考慮電力生產(chǎn)彈性系數(shù)；模型對電站建設進行了理想化的簡化假設，在每年初投入使用，且一有需要就能馬上提供，沒有考慮施工進程的限制不太符合實際，在實際的發(fā)電站建設中應投入更大。改進方向：現(xiàn)在企業(yè)都是向環(huán)保型企業(yè)發(fā)展了，社會基本設施也不能對環(huán)境擁有特權，考慮到火電站的發(fā)電原料主要是煤，里面含有硫物質(zhì)，會對大氣環(huán)境產(chǎn)生比較大的污染，如果不加以脫硫處理的話；而且煤是我國最重要的燃料資源之一，煤的儲存量正在減少。而水電站是完全綠色型的，且是可再生資源。故在發(fā)電站的新建上應優(yōu)先考慮水電站；我們把用電需求量只是簡單地處理成一個指數(shù)函數(shù)，實際情況遠非如此簡單，從長期時間段來看，應該近似滿足Logestic曲線，所以我們可以向著中長期方向發(fā)展，以求有更廣的適用范圍，更好的預測效果；我們的假設中提出無跨地區(qū)的輸電，在現(xiàn)實生活中卻存在著“西電東送”的現(xiàn)象，我們將考慮跨地區(qū)的輸電以減少浪費值，可以彌補成本；由于現(xiàn)實生活中可以有拉閘限電的情況存在，我們可以在裝機過程中使那條供應的折線圖略小于需求曲線，從生活用電量的節(jié)儉來降低成本。六、參考文獻[1]、吳添組工程經(jīng)濟（第二版）高等教育出版社2004[2]、蔡鎖章數(shù)學建模原理及方法海洋出版社2000[3]、姜啟源數(shù)學模型（第三版）高等教育出版社2003[4]、網(wǎng)址：/newsp/zgdl/dltj/gndltj26.htm/newsp/zgdl/dltj/gndltj23.htm/newsp/zgdl/dltj/gndltj27.htm七、附錄附錄一：用LINDO解線性規(guī)劃問題：min5957.52N1+5063.28N2+18910.8N3+2377.22N1+1997.28N2+6701.4N3st57840N1+87600N2+183600N3>=100000010N1+25N2+30N3>=180N1<=5N2<=4N3<=4endgin3附錄二：X中的值表示年需求量增率的對數(shù)>>x=[3.490379923.5154764413.5459253293.576341353.6135247033.6528263033.6966184593