火電廠有關(guān)煤的調(diào)運(yùn)問題編程

1、表 7 2011年湖南省火力發(fā)電年產(chǎn)量統(tǒng)計2011年湖南省火力發(fā)電年產(chǎn)量統(tǒng)計月份本月/億千萬時本月止累計/億千萬時累計增長%本月比同月增長%一月76.676.6-二月54.9131.5-1.1-4.2三月75.5207.3-0.70.3四月69.8277.18.247.3五月623391562六月52.5391-七月814732758八月86559523九月826422539十月73.3715.523.424.6十一月73791220十二月888792219 表8 省內(nèi)電煤點(diǎn)-湖南電廠的編號及所在位置編號供煤點(diǎn)及電廠東經(jīng)北緯編號供煤點(diǎn)及電廠東經(jīng)北緯1長沙市113.9328.2212湘西土家族苗

2、族自治州109.7328.302株洲市113.1327.8213大唐耒陽電廠112.8526.403湘潭市112.9328.8214大唐華銀金竹山電廠110.4327.684衡陽市112.5726.8815華潤漣源電廠110.6627.685邵陽市111.4727.2316韶能辰溪電廠110.1828.006常德市111.6829.0217大唐華銀金竹山電廠110.4327.687張家界市110.4729.1018華潤漣源電廠110.6627.688郴州市113.0025.7719華能岳陽電廠112.9727.879永州市113.6026.4220湘潭電廠113.1529.4310懷化市110

3、.0027.5521國電益陽電廠112.3728.5711婁底市113.9827.7022株洲電廠113.1527.784衡陽市112.5726.8815華潤漣源電廠110.6627.685邵陽市111.4727.2316韶能辰溪電廠110.1828.006常德市111.6829.0217大唐針對湖南省火力發(fā)電量情況，結(jié)合煤炭價格：（如下圖）2011年全國煙煤優(yōu)混平均價格走勢日期價格（單位：元/噸）2011年1月上半月6462011年1月下半月6492011年2月上半月6502011年2月下半月6522011年3月上半月6522011年3月下半月6532011年4月上半月6542011年4月下

4、半月6542011年5月上半月6582011年5月下半月6602011年6月上半月6612011年6月下半月6642011年7月上半月6662011年7月下半月6672011年8月上半月6662011年8月下半月6672011年9月上半月6682011年9月下半月6722011年10月上半月6772011年10月下半月6812011年11月上半月6832011年11月下半月6832011年12月上半月6832011年12月下半月6812012年1月上半月6812012年1月下半月681由以上數(shù)據(jù)表與曲線圖，可建立湖南省內(nèi)煤炭調(diào)度模型：對省內(nèi)主要火電廠及各省市進(jìn)行編號，如下圖：（表一）編號供煤點(diǎn)

5、及電廠東經(jīng)北緯編號供煤點(diǎn)及電廠東經(jīng)北緯1長沙市113.9328.2212湘西土家族苗族自治州109.7328.302株洲市113.1327.8213大唐耒陽電廠112.8526.403湘潭市112.9328.8214大唐華銀金竹山電廠110.4327.68華銀金竹山電廠110.4327.687張家界市110.4729.1018華潤漣源電廠110.6627.688郴州市113.0025.7719華能岳陽電廠112.9727.879永州市113.6026.4220湘潭電廠113.1529.4310懷化市110.0027.5521國電益陽電廠112.3728.5711婁底市113.9827.7022

6、株洲電廠113.1527.78模型通過表8所示的各電煤點(diǎn)-電廠的所在位置，可計算出i,j兩點(diǎn)間的直線距離D(i,j)D(i,j)= （1）設(shè)dk為實(shí)際距離與經(jīng)緯度下的距離的比值,realD為點(diǎn)點(diǎn)實(shí)際距離。由realD=dk*D （2） 可折算出實(shí)際距離。其中，由“地球緯度每相差1度實(shí)際距離相差111.13千米”， “地球經(jīng)圈長度由赤道依次向兩極遞減,赤道處長度為40075.7千米,赤道每相差1度實(shí)際距離相差111.32千米”可得：dk=157.3（km）設(shè)fd為每噸電煤每千米下的運(yùn)費(fèi)，danjia為每噸電煤的點(diǎn)點(diǎn)運(yùn)費(fèi)矩陣。則可得： danjia=fd*realD （3）其中，由“整車機(jī)保運(yùn)行基

7、價為0.0675元/噸公里”（數(shù)據(jù)來源： 鐵道部鐵路貨物運(yùn)價率表） 可得：fd=0.0675。截取矩陣danjia的第13列至第22列，第1行至第12行即可得電煤點(diǎn)電廠每噸電煤運(yùn)輸價格矩陣dan。根據(jù)長沙市、株洲市、湘潭市、衡陽市、邵陽市、常德市、張家界市、郴州市、永州市、懷化市、婁底市及湘西土家族苗族自治州的每噸電煤的購入價，列出煤炭點(diǎn)電廠每噸電煤出礦價格矩陣coalcost。表9 湖南省各縣市年產(chǎn)煤炭量統(tǒng)計表湖南省各縣市年產(chǎn)煤炭量統(tǒng)計表地級市縣級年生產(chǎn)能力（萬噸）煤炭均價（元/噸）長沙市瀏陽市18920寧鄉(xiāng)縣108920株洲市天元區(qū)6860株洲縣6860攸縣697860茶陵縣46860湘潭

8、市雨湖區(qū)6860湘潭縣142860衡陽市耒陽市701830常寧市46830邵陽市雙清區(qū)21890邵東縣85890新邵縣18890邵陽縣124890洞口縣8890新寧縣9890武岡市44890常德市澧縣69830石門縣120830張家界市永定區(qū)6850慈利縣26850桑植縣100850郴州市蘇仙區(qū)25860永興縣100860嘉苗縣27860資興市173860永州市冷水灘區(qū)96860祁陽市74860懷化市辰溪縣189.8860溆浦縣48860芷江侗族自治縣6860婁底市雙峰縣24860新化縣199860冷水江市429860漣源縣30860湘西土家族苗族自治州保靖縣20850龍山縣53850總計38

9、99.8- 設(shè)cost為電煤點(diǎn)電廠每噸電煤價格矩陣。這里認(rèn)為煤炭點(diǎn)電廠每噸電煤出礦價格與電煤點(diǎn)電廠每噸電煤運(yùn)輸價格之和即為電煤點(diǎn)電廠每噸電煤價格，也即cost=coalcost+dan （4） 表10 2011年湖南省各主要火電廠火力發(fā)電及需煤量單位:yuan/ton電廠火力發(fā)電量（億千萬時）單位耗煤量（克/千瓦/時）需煤量（萬噸）大唐耒陽電廠91.5936345.09314.8530大唐石門電廠69.2404351.22243.1861華潤鯉魚江電廠44.5701364.6162.5026韶能辰溪電廠36.3921372.8135.6697大唐華銀金竹山電廠165.0593390.2644.

10、0614華潤漣源電廠79.3260383.5304.2175華能岳陽電廠172.4195316.31545.3801湘潭電廠124.3056420.3522.4564國電益陽電廠104.1332375.1390.6036株洲電廠58.8816411.4242.2389總計945.922-3505.1695LB為格式輸出控制矩陣。 LB=0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

11、 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ； 最后，調(diào)用線性規(guī)劃函數(shù)wx,fval=linprogwcost,waeq,wbeq,wLB 其中，x代表調(diào)用線性規(guī)劃函數(shù)分別計算出每個電煤點(diǎn)-電廠的調(diào)煤噸數(shù),fval代表總的最短調(diào)煤噸數(shù)fval。編程如下：x=28.22 27.82 28.82 26.88 27.23 29.02 29.10 25.77 26.42 27.55 27.70 28.30 26.40 29.58

12、26.12 28.00 27.68 27.68 27.87 29.43 28.57 27.78; %電煤點(diǎn)、電廠地理坐標(biāo)輸入：北緯y=113.93 113.13 112.93 112.57 111.47 111.68 110.47 113.00 113.60 110.00 113.98 109.73 112.85 111.37 113.23 110.18 111.43 111.66 112.97 113.15 112.37 113.15; %電煤點(diǎn)、電廠地理坐標(biāo)輸入：東經(jīng)for i=1:22,for j=1:22,D(i,j)=sqrt(x(i)- x(j)2+(y(i)- y(j)2); %

13、計算距離矩陣Dend,End>> realD=157.3*D; %將經(jīng)緯度上的距離差折算為實(shí)際距離（單位：千米）>> danjia=0.0675*realD; %根據(jù)煤炭運(yùn)輸費(fèi)用將實(shí)際距離矩陣轉(zhuǎn)變?yōu)槊繃嶋娒旱狞c(diǎn)點(diǎn)運(yùn)費(fèi)矩陣>> dan=danjia(13:22,1:12); %截取點(diǎn)點(diǎn)間的單位煤炭運(yùn)費(fèi)矩陣，使之成為電煤點(diǎn)電廠每噸電煤矩陣>> coalcost=920860860830890830850860860850860850;920860860830890830850860860850860850;9208608608308908308508

14、60860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850; %列出煤炭點(diǎn)電廠的每噸電煤出口價格矩陣 >&

15、gt; cost=coalcost+dan %得到電煤點(diǎn)電廠的每頓電煤價格，再列出煤炭點(diǎn)電廠的每噸電煤出口價格矩陣與不等約束條件的右側(cè)結(jié)構(gòu)矩陣，調(diào)用線性規(guī)劃函數(shù)分別計算出每個電煤點(diǎn)-電廠的調(diào)煤噸數(shù)wx,并返回總的最短調(diào)煤噸數(shù)fval（注：MATLAB編程最終結(jié)果附最后）>>%列出煤炭點(diǎn)電廠的每噸電煤出口價格矩陣>> beq=40 251.7 49.3 249 103 63 44 108.3 56.7 81.3 227.3 24.3 128.2907 99.1 66.2 55.3 262.4 124.0 222.2 212.88 159.2 98.7； %列出不等約束條件

16、的右側(cè)結(jié)構(gòu)矩陣>> LB=0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ；>> x,fval=linprog(cost,aeq,beq,LB)

17、 %調(diào)用線性規(guī)劃函數(shù)分別計算出每個電煤點(diǎn)-電廠的調(diào)煤噸數(shù)wx,并返回總的最短調(diào)煤噸數(shù)fvalx = 0.6604 1.1589 1.0213 0.7949 0.9641 0.7727 0.2507 1.5647 0.8268 1.3089 114.3114 4.6558 0.7189 0.6896 1.6493 1.3666 51.0542 0.9586 0.3094 38.9261 0.7391 1.1430 0.6274 0.9177 0.6267 1.9336 3.3334 1.2090 8.5900 1.764236.3519 2.4966 2.0127 0.6896 3.4349 3

18、.7575 2.6433 7.5459 8.5946 30.2825 1.0899 0.6115 0.4142 1.0962 0.7511 0.8510 0.8497 0.5701 0.7237 0.9589 2.7291 205.0415 1.6549 4.1883 2.8467 31.7494 1.9536 3.2044 4.0448 3.30470.6262 107.9895 0.8807 2.1863 4.1265 1.6142 0.4808 1.6422 1.0246 0.91451.7739 0.7404 30.5320 1.9947 5.7654 2.5196 22.3085 1.0026 0.0003 16.719246.4178 1.3505 91