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

1、某地區(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è)方案安排為：在第一

2、階段，也即擴(kuò)建原有的舊火電站階段，我們預(yù)備在第一年擴(kuò)建兩臺單機(jī)容量為10千瓦的發(fā)電機(jī)組；而在第二年和第三年各擴(kuò)建一臺和兩臺；第二階段，也即新建水電站階段，從第四年初到第七年末的四年間，每年建設(shè)一臺單機(jī)容量為25千瓦的發(fā)電機(jī)組臺數(shù)；第三階段，也即新建火電站階段，在第八年和第九年初，各建設(shè)并投入使用一臺單機(jī)容量為30千瓦的發(fā)電機(jī)組；至此，全部建設(shè)完成并投入正常運(yùn)作。由此，我們算出最終的投入成本為23.3866億。【關(guān)鍵詞】經(jīng)濟(jì)最優(yōu)化、成本最小化、靜態(tài)模型、MATLAB一、 問題的重述 某地區(qū)在制定十年電力發(fā)展規(guī)劃時(shí)遇到一個(gè)問題：根據(jù)電力需求預(yù)測得知，該地區(qū)在十年后發(fā)電裝機(jī)容量需要增加180萬千瓦，

3、那時(shí)的年發(fā)電量需要增加100億度。根據(jù)調(diào)查和討論，電力規(guī)劃的備選技術(shù)方案有三個(gè)：1） 擴(kuò)建原有的火電站，但最多只能再安裝五臺10萬千瓦的發(fā)電機(jī)組；2） 新建水電站，但最多只能安裝四臺25萬千瓦的發(fā)電機(jī)組；3） 再新建一個(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ò)建舊火電站 2110501034110662新建水電站504702540057822

4、8043新建火電站24065304010336507 表 一表中負(fù)荷因子為全年滿功率運(yùn)行天數(shù)與全年總天數(shù)之比。根據(jù)該地實(shí)際調(diào)查原有火電站平均全年滿功率運(yùn)行天數(shù)為241天，水電站和新建的火電站應(yīng)分別為146天和255天，而全年365天，故折算得表中數(shù)據(jù)。表中資本回收因子是由如下數(shù)據(jù)所確定的，火電站的回收年限取15年，年利息率為006；水電站的回收年限取30年，年利息率為004，即得表中所列數(shù)值。我們的任務(wù)是：在滿足上述技術(shù)要求的前提下，選取經(jīng)濟(jì)效果最優(yōu)的建設(shè)方案.即：1. 原來的火電站應(yīng)如何擴(kuò)建？2. 新建的火電站應(yīng)如何確定單機(jī)容量為25萬千瓦的發(fā)電機(jī)組的數(shù)量？3. 新建的火力發(fā)電站應(yīng)如何確定單

5、機(jī)容量為30萬千瓦的發(fā)電機(jī)組的數(shù)量？二 、 問題的分析與假設(shè)因?yàn)槊颗_發(fā)電機(jī)在開始工作之后，每年必然會有所盈利，且由于在上限的限制上我們只有允許最大裝機(jī)臺數(shù)的限制，若如此，必然是在全部應(yīng)用時(shí)是最優(yōu)的，所以我們在確定這個(gè)模型的經(jīng)濟(jì)最優(yōu)解時(shí)，我們的目標(biāo)函數(shù)不應(yīng)是收入與支出的順差的最大值，而是將目標(biāo)函數(shù)定為滿足技術(shù)要求前提下，成本付出達(dá)到最小值。為此，我們做出了以下假設(shè)：1、發(fā)電機(jī)都是在每年的年初投入使用的，且每年的發(fā)電機(jī)都是滿負(fù)荷工作的，即負(fù)荷因子不變；2 、不計(jì)每年的發(fā)電剩余量的價(jià)值，即除去本地區(qū)的用電需求之外，無跨地區(qū)的輸電；3、十年中所有的發(fā)電機(jī)都是正常運(yùn)作，且不會出現(xiàn)拉閘限電現(xiàn)象，每年所發(fā)電

6、量滿足每年各月的用電需求量；4 、現(xiàn)階段的發(fā)電量剛好滿足生產(chǎn)生活所需的用電量；5、把職工的薪金，稅收等額外付費(fèi)每年當(dāng)作一固定開銷支付，故其對發(fā)電所得收益無波動(dòng)影響；三 、符號說明： 1，2，3分別表示擴(kuò)建的舊火電站，新建的水電站，以及新建的火電站；： 時(shí)間，單位：年；： 火電站和水電站最終所裝機(jī)的臺數(shù)；： 各新增的發(fā)電站的工程投資資金，單位：萬元；： 各發(fā)電站的年運(yùn)行成本，單位：萬元；： 各發(fā)電站的年運(yùn)行總成本，單位：萬元；： 該地區(qū)的用電量的年平均增長率；： 各發(fā)電站的貸款年利息率；： 火電站或水電站的負(fù)荷因子；： 水電站或火電站的回收期限，單位：年；： 各發(fā)電站的每年等額支付現(xiàn)金，單位：萬

7、元；： 現(xiàn)階段需求電量，也即現(xiàn)階段發(fā)電能力，單位：萬度；：發(fā)電站新增發(fā)電量所需支出，單位：萬元；： 經(jīng)資本回收折算后投資的每年應(yīng)支出量，單位：萬元；： 各新增的發(fā)電站的單機(jī)容量，單位：萬千瓦；： 各新增的發(fā)電站的裝機(jī)臺數(shù)；： 各新增的發(fā)電站的資本回收因子，單位：百萬元/億度；： 某年用電需求量，單位：萬度；： 擴(kuò)建的或增加的發(fā)電站的年發(fā)電量，單位：萬度；： 在年電的總供應(yīng)量；單位：萬度；： 擴(kuò)建的或增加的發(fā)電站的裝機(jī)容量，單位：萬度；四 、模型的建立與求解4.1 建立一次性投入模型在該模型的建立中，我們首先考慮最簡單的情況：第一年的投資就滿足十年以后的用電需求量。因?yàn)槊磕甑闹С銮闆r是一個(gè)定值，

8、我們就可以考慮轉(zhuǎn)化成每年的支出情況。得到折算后投資的年支出量根據(jù)經(jīng)濟(jì)學(xué)知識可知，現(xiàn)值在回收期限為年，年利率為情況下，其資本回收因子有如下公式：則等額系列償還資金有如下計(jì)算公式：類似的，對于對應(yīng)的擴(kuò)建的舊火電站和新建的水電站、火電站也有以下公式表示它們的等額償還資金： 發(fā)電站投入使用之后，各發(fā)電站的年運(yùn)行總成本經(jīng)分析，年運(yùn)行總成本是與該發(fā)電站的年發(fā)電量成正比關(guān)系的，而年發(fā)電量則取決于火電站和水電站的裝機(jī)臺數(shù)、單機(jī)容量和負(fù)荷因子。有如下關(guān)系式：新的水電站和火電站的投資都有個(gè)前期工程投資，但是因?yàn)槭莻€(gè)固定的常數(shù)，所以它所引起的等額償還資金也是個(gè)常數(shù)，我們這里在考慮成本最小化時(shí)可以忽略它的影響。所以，

9、根據(jù)題目的技術(shù)要求，我們可以將目標(biāo)函數(shù)表示如下：但是我們是有所約束的，為非限制性約束：用LINDO解得各個(gè)相應(yīng)的值為：具體實(shí)現(xiàn)程序見附錄一。4.2 考慮該地區(qū)用電需求量以遞增的趨勢增長的模型 確定用電年增長量和年需求增長率因?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)階段用電需求量，代表原來的電力需求量，代表年平均

10、需求增長率，代表年限，即第幾年。我們利用網(wǎng)絡(luò)資源找到了一組全國的從1981年至2000年的歷年電力裝機(jī)和發(fā)電量的構(gòu)成比的數(shù)據(jù)，如下表：   歷年電力裝機(jī)和發(fā)電量的構(gòu)成比（19812000）      年 份 裝機(jī)容量（萬千瓦） 比 重（%） 發(fā)電量（億千瓦時(shí)） 比 重（%） 水 電 火 電 水 電 火 電 水 電 火 電 水 電 火 電 1981 2193 4720 31.7 68.3 656 2437 21.2 78.8 1982 2296 4940 31.7 68.3 744 2533 22.7 77.3 1983 24

11、16 5228 31.6 68.4 864 2651 24.6 75.4 1984 2560 5452 32.0 68.0 868 2902 23.0 77.0 1985 2641 6064 30.3 69.7 924 3183 22.5 77.5 1986 2754 6628 29.4 70.6 945 3551 21.0 79.0 1987 3019 7271 29.3 70.7 1002 3971 20.1 79.9 1988 3270 8280 28.3 71.7 1092 4359 20.0 80.0 1989 3458 9206 27.3 72.7 1185 4662 20.3 7

12、9.7 1990 3605 10184 26.1 73.9 1263 4950 20.3 79.7 1991 3788 11359 25.0 75.0 1248 5527 18.4 81.6 1992 4068 12585 24.4 75.6 1315 6227 17.4 82.6 1993 4489 13802 24.5 75.5 1516 6868 18.1 81.9 1994 4906 14874 24.5 74.4 1668 7470 18.0 80.5 1995 5218 16294 24.0 75.0 1868 8074 18.6 80.2 1996 5558 17886 23.5

13、 75.6 1869 8781 17.3 81.3 1997 5973 19241 23.5 75.6 1946 9252 17.2 81.6 1998 6507 20988 23.5 75.7 2043 9388 17.6 81.1 1999 7297 22343 24.4 74.8 2129 10047 17.3 81.5 2000 7935 23754 24.9 74.4 2431 11079 17.8 81.0 表 二我們假定上面這組數(shù)據(jù)中的電的需求和供應(yīng)量始終是平衡的，這樣我們就可以從標(biāo)了標(biāo)記的兩組數(shù)據(jù)的和的變化去求得這個(gè)年增長率。在對數(shù)據(jù)進(jìn)行了前期取對數(shù)之后，我們利用MATLAB軟

14、件對數(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è)。解出水電站和火電站的建造順序和建造時(shí)間正如我們的模型的分析中所述，每臺發(fā)電機(jī)投入生產(chǎn)之后，其每年都在盈利，經(jīng)計(jì)算比較，

15、其金額比發(fā)電站的年支出總資金大一個(gè)數(shù)量級，故直接考慮順差的差值大小來決定先建造哪一個(gè)發(fā)電站是不合理的，我們對于建造順序的安排也取決于如何使所花成本最小。由于每年的支付金額都是等額的，而且建成發(fā)電站之后全部在負(fù)荷因子下工作，所以各個(gè)發(fā)電站的年運(yùn)行成本也是一個(gè)固定的數(shù)值，所以我們可以考慮用一年的支出總金額來代替十年支出總金額。在一年的費(fèi)用支出中，主要有兩方面的因素：1） 每年運(yùn)行時(shí)的年運(yùn)行成本；2） 購買機(jī)器時(shí)的工程投入所引起的等額償還資金。 在上一模型中，我們已經(jīng)知道某一發(fā)電站的年運(yùn)行總成本是與該發(fā)電站的年發(fā)電量成正比關(guān)系的，年發(fā)電量又取決于火電站和水電站的裝機(jī)臺數(shù)、單機(jī)容量和負(fù)荷因子。工程投資

16、所引起的等額償還資金也可由上一模型中的公式可以算得。故水電站或火電站的年費(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ò)建舊火電站之后，在

17、第3年開始建造并使用新建的水電站，第七年建造使用新的火電站。利用MATLAB作圖如下（具體程序?qū)崿F(xiàn)見附錄三）： 圖 二4.3 考慮裝機(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.

18、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年份67

19、8910年需求電量131.0953142.1073154.0443166.9840181.0107年裝機(jī)臺數(shù)年初完成11110年發(fā)電量136144.76163.12181.48181.48表 四現(xiàn)利用表四用MATLAB作出其圖形如下： 圖 三4.4 求出最終的成本由于每年的發(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)：1） 本模型在考慮最經(jīng)濟(jì)的方案時(shí)，并沒有直接考慮其收益的最大值，

20、而是取其成本的最小值；2） 本模型將復(fù)雜的動(dòng)態(tài)規(guī)劃問題通過合理的假設(shè)簡化成靜態(tài)的規(guī)劃問題；3） 本模型對問題以點(diǎn)連線，步步深入。但并沒有什么晦澀難懂的地方，易于理解；模型的缺點(diǎn)：1） 本模型在考慮發(fā)電能力時(shí)，沒有考慮電力生產(chǎn)彈性系數(shù)；2） 模型對電站建設(shè)進(jìn)行了理想化的簡化假設(shè)，在每年初投入使用，且一有需要就能馬上提供，沒有考慮施工進(jìn)程的限制不太符合實(shí)際，在實(shí)際的發(fā)電站建設(shè)中應(yīng)投入更大。改進(jìn)方向：1） 現(xiàn)在企業(yè)都是向環(huán)保型企業(yè)發(fā)展了，社會基本設(shè)施也不能對環(huán)境擁有特權(quán)，考慮到火電站的發(fā)電原料主要是煤，里面含有硫物質(zhì)，會對大氣環(huán)境產(chǎn)生比較大的污染，如果不加以脫硫處理的話；而且煤是我國最重要的燃料資源

21、之一，煤的儲存量正在減少。而水電站是完全綠色型的，且是可再生資源。故在發(fā)電站的新建上應(yīng)優(yōu)先考慮水電站；2） 我們把用電需求量只是簡單地處理成一個(gè)指數(shù)函數(shù)，實(shí)際情況遠(yuǎn)非如此簡單，從長期時(shí)間段來看，應(yīng)該近似滿足Logestic曲線，所以我們可以向著中長期方向發(fā)展，以求有更廣的適用范圍，更好的預(yù)測效果；3） 我們的假設(shè)中提出無跨地區(qū)的輸電，在現(xiàn)實(shí)生活中卻存在著“西電東送”的現(xiàn)象，我們將考慮跨地區(qū)的輸電以減少浪費(fèi)值，可以彌補(bǔ)成本；4） 由于現(xiàn)實(shí)生活中可以有拉閘限電的情況存在，我們可以在裝機(jī)過程中使那條供應(yīng)的折線圖略小于需求曲線，從生活用電量的節(jié)儉來降低成本。六、 參考文獻(xiàn)1、吳添組 工程經(jīng)濟(jì)（第二版）

22、 高等教育出版社 20042、蔡鎖章 數(shù)學(xué)建模原理及方法 海洋出版社 2000 3、姜啟源 數(shù)學(xué)模型（第三版） 高等教育出版社 20034、網(wǎng)址： 七、 附錄附錄一：用LINDO解線性規(guī)劃問題：min 5957.52N1+5063.28N2+18910.8N3+2377.22N1+1997.28N2+6701.4N3st57840N1+87600N2+183600N3>=100000010N1+25N2+30N3>=180N1<=5N2<=4N3<=4endgin 3附錄二：X中的值表示年需求量增率的對數(shù)>>x=3.490379923.5154764413.5459253293.576341353.6135247033.6528263033.696