发明名称 一种钢铁企业发电成本优化的方法
摘要 一种钢铁企业发电成本优化的方法,属于钢铁企业综合节电技术领域。步骤包括:确定需要进行发电成本优化的调度周期T,获取机组发电成本优化参数;获取用于日前发电计划编制的网络断面;确定周期内考虑煤气掺烧情况下的发电机组燃料成本;根据企业机组能耗,外购电成本及外送电收益实际情况建立以全系统发电成本最小为目标的优化调度模型;依据优化求解获得企业自发电机组出力,外购电量,外送电量的结果等。优点在于,可满足生产负荷需要、发电机组出力、线路传输容量等任务,同时保证企业电网的安全稳定运行,结合负荷预测曲线,生产和检修计划,进一步优化调整自发电出力,以及自发电和外购电的比例关系,使企业发电成本最小。
申请公布号 CN103617552B 申请公布日期 2016.11.16
申请号 CN201310598477.X 申请日期 2013.11.22
申请人 冶金自动化研究设计院 发明人 刘庆贺;于立业;徐化岩;曾玉娇;贾天云;苏胜石;赵博;朱寅;盛刚
分类号 G06Q10/04(2012.01)I;G06Q50/06(2012.01)I;G06Q50/04(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 北京华谊知识产权代理有限公司 11207 代理人 刘月娥
主权项 一种钢铁企业发电成本优化的方法,其特征在于,包括以下步骤:步骤一,确定需要进行发电成本优化的调度周期T,获取机组发电成本优化参数,这些参数包括:调度周期内企业负荷预测曲线、生产和检修计划、外购送电计划、机组的启停状态、可用状态、固定的发电机组出力计划,调整机组出力计划;步骤二,获取用于日前发电计划编制的网络断面,网络断面是指在电网安全或电量交易上相互关联的一组输电线路,同时结合企业生产检修计划,利用电网模型自动生成各时段的网络拓扑,并计算各时段的潮流分布及监视元件影响系数;步骤三,确定周期内考虑煤气掺烧情况下的发电机组燃料成本:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>C</mi><mi>zf</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>[</mo><mi>F</mi><msub><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&times;</mo><msubsup><mi>C</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow><mi>f</mi></msubsup><mo>]</mo></mrow>]]></math><img file="FDA0000420360380000011.GIF" wi="650" he="143" /></maths>F(P)<sub>i,t</sub>=aP<sub>i,t</sub><sup>2</sup>+bP<sub>i,t</sub>+c<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>C</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow><mi>f</mi></msubsup><mo>=</mo><msup><mi>C</mi><mi>coal</mi></msup><mo>&times;</mo><mfrac><msup><mi>Q</mi><mi>B</mi></msup><msup><mi>Q</mi><mi>coal</mi></msup></mfrac><mo>&times;</mo><msubsup><mi>Y</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow><mi>coal</mi></msubsup><mo>+</mo><msup><mi>C</mi><mi>bfg</mi></msup><mo>&times;</mo><mfrac><msup><mi>Q</mi><mi>B</mi></msup><msup><mi>Q</mi><mi>bfg</mi></msup></mfrac><mo>&times;</mo><msubsup><mi>Y</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow><mi>bfg</mi></msubsup><mo>+</mo><msup><mi>C</mi><mi>cog</mi></msup><mo>&times;</mo><mfrac><msup><mi>Q</mi><mi>B</mi></msup><msup><mi>Q</mi><mi>cog</mi></msup></mfrac><mo>&times;</mo><msubsup><mi>Y</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow><mi>cog</mi></msubsup></mrow>]]></math><img file="FDA0000420360380000012.GIF" wi="1211" he="139" /></maths>其中,C<sub>zf</sub>为发电机组燃料成本,F(P)<sub>i,t</sub>为第i台机组t时段的标煤耗,P<sub>i,t</sub>为第i台机组t时段的有功出力,<img file="FDA0000420360380000013.GIF" wi="77" he="83" />为第i台机组t时段煤与煤气混合燃烧的燃料单价,C<sup>coal</sup>、C<sup>bfg</sup>、C<sup>cog</sup>分别为外购煤、高炉煤气,焦炉煤气的价格;Q<sup>coal</sup>、Q<sup>bfg</sup>、Q<sup>cog</sup>分别为煤、高炉煤气,焦炉煤气的热值;其中Q<sup>B</sup>为标煤热值,其值恒定为29308KJ/Kg;<img file="FDA0000420360380000014.GIF" wi="248" he="83" /><img file="FDA0000420360380000015.GIF" wi="85" he="75" />分别为第i台机组t时段煤、高炉煤气,焦炉煤气占总燃料热量的比例;考虑煤气掺烧的情况,从企业现有的系统中读取煤,高炉煤气,焦炉煤气的量以及其对应的机组的有功出力,将煤,高炉煤气,焦炉煤气折合成标煤量,然后拟合得到一个周期T内的机组能耗模型F(P)<sub>i,t</sub>=aP<sub>i,t</sub><sup>2</sup>+bP<sub>i,t</sub>+c其中a,b,c为拟合系数;步骤四,根据企业机组能耗,外购电成本及外送电收益实际情况建立以全系统发电成本最小为目标的优化调度模型;以15分钟作为一个优化的逻辑时段,以调度周期内负荷预测曲线作为研究对象,结合生产和检修计划,优化企业内部调度发电机组的出力,优化目标是企业自发电成本最小;钢铁企业发电成本优化模型为,目标函数:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>min</mi><msub><mi>C</mi><mrow><mi>cos</mi><mi>t</mi></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mi>t</mi></munder><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>C</mi><mrow><mi>zf</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow></msub><mo>+</mo><msub><mi>C</mi><mi>wg</mi></msub><mo>-</mo><msub><mi>C</mi><mi>ws</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000420360380000021.GIF" wi="935" he="140" /></maths>约束条件:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>C</mi><mi>wg</mi></msub><mo>=</mo><munder><mi>&Sigma;</mi><mi>t</mi></munder><msub><mi>P</mi><mrow><mi>w</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&times;</mo><msub><mi>C</mi><mi>t</mi></msub><mo>&times;</mo><msub><mi>&delta;</mi><mi>t</mi></msub></mrow>]]></math><img file="FDA0000420360380000022.GIF" wi="431" he="113" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>C</mi><mi>ws</mi></msub><mo>=</mo><munder><mi>&Sigma;</mi><mi>t</mi></munder><msub><mi>P</mi><mrow><mi>ws</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&times;</mo><msub><mi>C</mi><mi>s</mi></msub><mo>&times;</mo><msub><mi>&sigma;</mi><mi>t</mi></msub></mrow>]]></math><img file="FDA0000420360380000023.GIF" wi="451" he="113" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>w</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>load</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow>]]></math><img file="FDA0000420360380000024.GIF" wi="459" he="139" /></maths>P<sub>i,min</sub>≤P<sub>i,t</sub>≤P<sub>i,max</sub>‑UR<sub>i</sub>≤P<sub>i,t</sub>‑P<sub>i,t‑1</sub>≤UR<sub>i</sub>,i=1,2...K<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>ij</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>&Element;</mo><mi>M</mi></mrow></munder><mo>[</mo><mi>P</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>]</mo><msub><mi>&omega;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow>]]></math><img file="FDA0000420360380000025.GIF" wi="632" he="115" /></maths>P<sub>ij,min</sub>≤P<sub>ij,t</sub>≤P<sub>ij,max</sub><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>w</mi><mo>,</mo><mi>min</mi></mrow></msub><mo>&le;</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mn>24</mn></munderover><msub><mi>P</mi><mrow><mi>w</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mrow><mi>w</mi><mo>,</mo><mi>max</mi></mrow></msub></mrow>]]></math><img file="FDA0000420360380000026.GIF" wi="457" he="139" /></maths>其中,K为企业内可调度的自备电厂发电机组数,N为系统调度周期T内所包含的时段数,C<sub>zf(i,t)</sub>为第i台机组t时段的自发电成本,C<sub>wg</sub>外购电成本,P<sub>w,t</sub>外购电量,C<sub>t</sub>为t时段的外购电价,δ<sub>t</sub>为0,1量表征有无外供电,C<sub>ws</sub>外送电收益,P<sub>ws,t</sub>外送电量,C<sub>s</sub>为外送电价,σ<sub>t</sub>为0,1量表征无外送电,P<sub>i,t</sub>为i台机组t时段的有功出力,P<sub>load,t</sub>为t时段企业负荷预测值,P<sub>i,min</sub>,P<sub>i,max</sub>分别为机组i出力及其上下限,‑UR<sub>i</sub>,UR<sub>i</sub>为机组i每时段可加减负荷的最大值,P<sub>ij,t</sub>为支路ij在t时间段上的潮流,P(i,t)为节点i在t时间段的注入功率,M为节点i下包括的监视元件数,S<sub>i,t</sub>为节点负荷功率,ω<sub>i,j,t</sub>为节点注入功率对支路ij在t时间段上的影响系数,P<sub>ij,min</sub>,P<sub>ij,max</sub>表示支路ij潮流的上下限,P<sub>w,min</sub>,P<sub>w,max</sub>别表示发电机i在时间T内由合同确定的供购电量上下限;步骤五,依据优化求解获得企业自发电机组出力,外购电量,外送电量的结果,考虑全部网络监视元件,对各时段进行安全校核,若没有新增监视元件电压或电流越限,则进入步骤六,否则计算新增越限监视元件的影响系数信息,然后进入步骤二;步骤六,迭代结束,得出钢铁企业自发电计划,优化结束,给出成本最小的自发电方案。
地址 100071 北京市丰台区西四环南路72号