火电厂有关煤的调运问题编程

1、表 7 2011年湖南省火力发电年产量统计2011年湖南省火力发电年产量统计月份本月/亿千万时本月止累计/亿千万时累计增长%本月比同月增长%一月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 省内电煤点-湖南电厂的编号及所在位置编号供煤点及电厂东经北纬编号供煤点及电厂东经北纬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大唐针对湖南省火力发电量情况，结合煤炭价格：（如下图）2011年全国烟煤优混平均价格走势日期价格（单位：元/吨）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由以上数据表与曲线图，可建立湖南省内煤炭调度模型：对省内主要火电厂及各省市进行编号，如下图：（表一）编号供煤点

5、及电厂东经北纬编号供煤点及电厂东经北纬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所示的各电煤点-电厂的所在位置，可计算出i,j两点间的直线距离D(i,j)D(i,j)= （1）设dk为实际距离与经纬度下的距离的比值,realD为点点实际距离。由realD=dk*D （2） 可折算出实际距离。其中，由“地球纬度每相差1度实际距离相差111.13千米”， “地球经圈长度由赤道依次向两极递减,赤道处长度为40075.7千米,赤道每相差1度实际距离相差111.32千米”可得：dk=157.3（km）设fd为每吨电煤每千米下的运费，danjia为每吨电煤的点点运费矩阵。则可得： danjia=fd*realD （3）其中，由“整车机保运行基

7、价为0.0675元/吨公里”（数据来源： 铁道部铁路货物运价率表） 可得：fd=0.0675。截取矩阵danjia的第13列至第22列，第1行至第12行即可得电煤点电厂每吨电煤运输价格矩阵dan。根据长沙市、株洲市、湘潭市、衡阳市、邵阳市、常德市、张家界市、郴州市、永州市、怀化市、娄底市及湘西土家族苗族自治州的每吨电煤的购入价，列出煤炭点电厂每吨电煤出矿价格矩阵coalcost。表9 湖南省各县市年产煤炭量统计表湖南省各县市年产煤炭量统计表地级市县级年生产能力（万吨）煤炭均价（元/吨）长沙市浏阳市18920宁乡县108920株洲市天元区6860株洲县6860攸县697860茶陵县46860湘潭

8、市雨湖区6860湘潭县142860衡阳市耒阳市701830常宁市46830邵阳市双清区21890邵东县85890新邵县18890邵阳县124890洞口县8890新宁县9890武冈市44890常德市澧县69830石门县120830张家界市永定区6850慈利县26850桑植县100850郴州市苏仙区25860永兴县100860嘉苗县27860资兴市173860永州市冷水滩区96860祁阳市74860怀化市辰溪县189.8860溆浦县48860芷江侗族自治县6860娄底市双峰县24860新化县199860冷水江市429860涟源县30860湘西土家族苗族自治州保靖县20850龙山县53850总计38

9、99.8- 设cost为电煤点电厂每吨电煤价格矩阵。这里认为煤炭点电厂每吨电煤出矿价格与电煤点电厂每吨电煤运输价格之和即为电煤点电厂每吨电煤价格，也即cost=coalcost+dan （4） 表10 2011年湖南省各主要火电厂火力发电及需煤量单位:yuan/ton电厂火力发电量（亿千万时）单位耗煤量（克/千瓦/时）需煤量（万吨）大唐耒阳电厂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 ； 最后，调用线性规划函数wx,fval=linprogwcost,waeq,wbeq,wLB 其中，x代表调用线性规划函数分别计算出每个电煤点-电厂的调煤吨数,fval代表总的最短调煤吨数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; %电煤点、电厂地理坐标输入：北纬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; %电煤点、电厂地理坐标输入：东经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; %将经纬度上的距离差折算为实际距离（单位：千米）>> danjia=0.0675*realD; %根据煤炭运输费用将实际距离矩阵转变为每吨电煤的点点运费矩阵>> dan=danjia(13:22,1:12); %截取点点间的单位煤炭运费矩阵，使之成为电煤点电厂每吨电煤矩阵>> coalcost=920860860830890830850860860850860850;920860860830890830850860860850860850;9208608608308908308508

14、60860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850;920860860830890830850860860850860850; %列出煤炭点电厂的每吨电煤出口价格矩阵 >&

15、gt; cost=coalcost+dan %得到电煤点电厂的每顿电煤价格，再列出煤炭点电厂的每吨电煤出口价格矩阵与不等约束条件的右侧结构矩阵，调用线性规划函数分别计算出每个电煤点-电厂的调煤吨数wx,并返回总的最短调煤吨数fval（注：MATLAB编程最终结果附最后）>>%列出煤炭点电厂的每吨电煤出口价格矩阵>> 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、的右侧结构矩阵>> 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、 %调用线性规划函数分别计算出每个电煤点-电厂的调煤吨数wx,并返回总的最短调煤吨数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