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

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

工程投資單機(jī)容量

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

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

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

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