基于動(dòng)態(tài)燃料成本雙層優(yōu)化的燃煤發(fā)電企業(yè)管控模型

廖志偉 張文錦 陳琳韜 黃杰棟 莊競(jìng)

摘 ? 要：為使燃煤發(fā)電企業(yè)降低燃料成本，實(shí)現(xiàn)燃料科學(xué)管控，提出一種綜合考慮發(fā)電可靠性和生產(chǎn)經(jīng)濟(jì)性的燃料成本雙層優(yōu)化模型，上層為電煤采購(gòu)庫(kù)存混合整數(shù)規(guī)劃模型，下層為混煤摻燒非線性優(yōu)化模型，雙層模型交替迭代協(xié)調(diào)優(yōu)化，實(shí)現(xiàn)燃料管控動(dòng)態(tài)決策. 針對(duì)優(yōu)化模型高維多約束情況，提出一種混沌映射初始化和高斯變異結(jié)合的改進(jìn)灰狼優(yōu)化算法. 以某燃煤電廠模擬運(yùn)營(yíng)進(jìn)行算例驗(yàn)證分析，結(jié)果表明，所提模型可以降低燃煤電廠7.80%的燃煤成本，證明了所提模型與求解算法的有效性和可行性.

關(guān)鍵詞：燃煤電廠;動(dòng)態(tài)規(guī)劃;電煤采購(gòu);混煤摻燒;灰狼優(yōu)化算法

中圖分類號(hào)：TM621 ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

Management and Control Model of Coal-fired Power Generation

Enterprise Based on Dynamic Fuel Cost Bi-level Optimization

LIAO Zhiwei?，ZHANG Wenjin，CHEN Lintao，HUANG Jiedong，ZHAUNG Jing

（School of Electric Power Engineering，South China University of Technology，Guangzhou 510641，China）

Abstract：In order to reduce fuel costs and achieve scientific fuel control in coal-fired power generation enterprises， a bi-level fuel cost optimization model that comprehensively considers the reliability of power generation and the economy of production is proposed in this paper. The upper layer is a mixed integer programming model for power coal purchasing and inventory， while the lower layer is a nonlinear optimization model for mixed coal combustion. The bi-level model is the alternating iterative coordination optimization to realize the dynamic decision-making of fuel control. Then， an improved gray wolf optimization algorithm combining chaos mapping initialization and Gaussian mutation is proposed for the high dimensional and multiple constraints of the optimization model. Finally， a simulated operation of a coal-fired power plant is used to verify the results. The results show that the proposed model can reduce the coal cost of coal fired power plant by 7.80%， which proves the effectiveness and feasibility of the proposed model and solution algorithm.

Key words：coal fired power plants;dynamic programming;electric coal purchase;mixed coal blending;grey wolf optimization algorithm

隨著電煤價(jià)格雙軌制取消以及供給側(cè)結(jié)構(gòu)性改革引起電煤市場(chǎng)波動(dòng)性增強(qiáng)，電煤價(jià)格連年走高[1]，直接導(dǎo)致燃煤發(fā)電企業(yè)生產(chǎn)成本提高;與此同時(shí)，風(fēng)電、光電等低邊際成本可再生能源參與電力市場(chǎng)競(jìng)爭(zhēng)，導(dǎo)致市場(chǎng)出清價(jià)格接近或低于燃煤機(jī)組發(fā)電成本，并使得燃煤電廠中標(biāo)發(fā)電量進(jìn)一步降低[2]. 以燃煤電廠為代表的傳統(tǒng)火電企業(yè)陷入電力市場(chǎng)運(yùn)營(yíng)困境，甚至出現(xiàn)虧損情況[3]. 因此，火電企業(yè)亟須優(yōu)化現(xiàn)行運(yùn)營(yíng)策略以降低發(fā)電成本從而增強(qiáng)電力市場(chǎng)競(jìng)爭(zhēng)力.

燃煤電廠的燃料成本占據(jù)發(fā)電成本的70%～80%，分別受電煤采購(gòu)入廠環(huán)節(jié)、電煤庫(kù)存管理環(huán)節(jié)及電煤入爐摻燒環(huán)節(jié)影響[4]，國(guó)內(nèi)外不少學(xué)者對(duì)燃煤電廠燃料管控策略展開(kāi)研究. 文獻(xiàn)[5-7]將電廠電煤庫(kù)存管理與電煤采購(gòu)決策相聯(lián)系，對(duì)電煤庫(kù)存管理問(wèn)題建立混合整數(shù)線性規(guī)劃模型，在不同煤炭供應(yīng)商中選擇最優(yōu)采購(gòu)組合方案. 文獻(xiàn)[8-9]討論不同風(fēng)險(xiǎn)類型電廠決策者考慮現(xiàn)貨價(jià)格波動(dòng)情況下的不同燃煤采購(gòu)策略，在不同煤炭供應(yīng)商中求解最優(yōu)采購(gòu)組合方案. 文獻(xiàn)[10]結(jié)合動(dòng)力煤期貨和現(xiàn)貨市場(chǎng)，優(yōu)化電廠庫(kù)存策略，減少電煤庫(kù)存對(duì)電廠流動(dòng)資金的占用. 文獻(xiàn)[11]以電廠年度收益條件風(fēng)險(xiǎn)價(jià)值最小為目標(biāo)，考慮煤價(jià)和上網(wǎng)電價(jià)不確定性，運(yùn)用線性規(guī)劃對(duì)電廠燃煤庫(kù)存優(yōu)化模型求解. 然而，上述文獻(xiàn)研究燃煤電廠庫(kù)存管理策略和采購(gòu)方案是基于電煤發(fā)電消耗固定的前提，根據(jù)電煤價(jià)格波動(dòng)，盡可能選擇購(gòu)買低價(jià)煤，忽視電廠內(nèi)鍋爐混煤摻燒比例對(duì)發(fā)電燃料成本的影響. 為此，也有不少學(xué)者研究鍋爐最佳混煤摻燒方案. 文獻(xiàn)[12]通過(guò)調(diào)整混煤摻燒方案使其燃燒特性接近鍋爐設(shè)計(jì)目標(biāo)煤質(zhì)，并優(yōu)化煤場(chǎng)庫(kù)存存放位置減少入爐混煤調(diào)運(yùn)時(shí)間. 文獻(xiàn)[13-15]研究不同混煤摻燒方案下對(duì)鍋爐氣體排放和碳化合物排放的影響，優(yōu)化鍋爐燃燒對(duì)環(huán)境污染的排放. 文獻(xiàn)[16]建立入爐混煤標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)價(jià)差模型計(jì)算燃煤電廠燃料成本，證明直接摻燒便宜劣煤不能有效降低燃煤成本. 文獻(xiàn)[17-18]兼顧燃煤電廠鍋爐運(yùn)行的安全性和經(jīng)濟(jì)性，通過(guò)摻燒試驗(yàn)得出保證鍋爐安全運(yùn)行下經(jīng)濟(jì)性較好的混煤摻燒方案.

綜上可知，傳統(tǒng)火力電廠由于信息系統(tǒng)孤島及燃料部門(mén)和采購(gòu)部門(mén)對(duì)發(fā)電可靠性和生產(chǎn)經(jīng)濟(jì)性的側(cè)重點(diǎn)不同，混煤摻燒方案難以與電煤采購(gòu)方案實(shí)現(xiàn)協(xié)同全局優(yōu)化燃料管控.

隨著計(jì)算機(jī)、數(shù)字化和大數(shù)據(jù)技術(shù)的發(fā)展，使燃料管控整體動(dòng)態(tài)全局尋優(yōu)成為可能，本文提出一種綜合考慮發(fā)電可靠性和生產(chǎn)經(jīng)濟(jì)性的燃料成本雙層優(yōu)化模型，通過(guò)多個(gè)環(huán)節(jié)動(dòng)態(tài)耦合，將部分環(huán)節(jié)局部靜態(tài)優(yōu)化轉(zhuǎn)換成全局動(dòng)態(tài)優(yōu)化燃料管控方案. 首先，梳理燃煤電廠傳統(tǒng)燃料管控流程工作邏輯，分析燃料管控對(duì)發(fā)電成本的影響;其次，分別構(gòu)建上層電煤采購(gòu)庫(kù)存混合整數(shù)規(guī)劃模型和下層鍋爐混煤摻燒非線性優(yōu)化模型，其中，上層模型以采購(gòu)成本和庫(kù)存成本最低為目標(biāo)函數(shù)，確定采購(gòu)方案，向下層模型傳遞可用電煤;下層模型以鍋爐燃燒混煤平均標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)最低為優(yōu)化目標(biāo)，確定摻燒方案，并向上層模型傳遞電煤需求計(jì)劃. 然后，將上下層模型交替迭代協(xié)調(diào)優(yōu)化，并使用混沌映射初始化和高斯變異結(jié)合的改進(jìn)灰狼優(yōu)化算法求解燃料成本雙層優(yōu)化模型;最后，以某燃煤電廠為算例驗(yàn)證雙層優(yōu)化模型的可行性.

1 ? 燃煤電廠運(yùn)營(yíng)機(jī)制及雙層優(yōu)化模型

1.1 ? 燃煤電廠運(yùn)營(yíng)機(jī)制成本模型

燃煤電廠的燃料成本與其燃料管控水平息息相關(guān)，電煤從采購(gòu)到生產(chǎn)的全過(guò)程如圖1所示.

影響燃料成本的過(guò)程主要為廠外采購(gòu)環(huán)節(jié)、電煤庫(kù)存管理環(huán)節(jié)和鍋爐混煤摻燒環(huán)節(jié). 在不同的環(huán)節(jié)，發(fā)電燃料成本有不同的表現(xiàn)形式，為方便本文描述，做如下定義：

入廠原煤?jiǎn)蝺r(jià)pi，m，反映電廠采購(gòu)水平：

pi，m = ? ? ? ? ? ? （1）

入爐標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)p*

j，m，體現(xiàn)庫(kù)存管理水平：

p*

j，m={+}

（2）

入爐混煤標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)p*

m，決定發(fā)電燃料成本：

p*

m = εj，m p*

j，m ? ? ? ? （3）

式中：pi，m為第m周第i種電煤入廠原煤?jiǎn)蝺r(jià);pc

i，m為港口煤?jiǎn)蝺r(jià);pt

m為入廠運(yùn)費(fèi)單價(jià);ph

m為入廠印花稅; pz

m為入廠雜散費(fèi);ηm為運(yùn)輸損耗;p*

j，m為第m周第j類電煤的入爐標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià);Pj，m為當(dāng)周該類煤采購(gòu)量;Sj，m-1為上周該類煤囤積量;φ為電煤庫(kù)存管理費(fèi)用單價(jià);θ為電煤囤積導(dǎo)致的資金積壓利息月利率; Q*

net、Qnet，j分別代表標(biāo)準(zhǔn)煤和對(duì)應(yīng)電煤低位發(fā)熱量; Nj為燃煤電廠庫(kù)存電煤分類種數(shù);p*

m為第m周混煤入爐標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià);εj，m為第j類電煤在第m周的摻燒比例.

1.2 ? 燃料成本雙層優(yōu)化模型

由1.1節(jié)可知，采購(gòu)環(huán)節(jié)與庫(kù)存管理直接相關(guān)，本文針對(duì)對(duì)應(yīng)環(huán)節(jié)分別建立電煤采購(gòu)庫(kù)存模型和混煤摻燒優(yōu)化模型. 考慮到采購(gòu)環(huán)節(jié)和配煤環(huán)節(jié)是相互影響、相互制約的關(guān)系，將二者構(gòu)建為如圖2所示的燃料成本雙層優(yōu)化模型.

上層電煤采購(gòu)庫(kù)存模型中，以電煤采購(gòu)總成本最低為目標(biāo)函數(shù)，采購(gòu)周期內(nèi)各類電煤采購(gòu)量作為決策變量，考慮電廠庫(kù)存和煤炭供應(yīng)滿足最小發(fā)電計(jì)劃需求從而保證發(fā)電可靠性，求解采購(gòu)庫(kù)存方案并向下層混煤摻燒優(yōu)化模型傳遞電煤入廠標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)和可用電煤.

下層混煤摻燒優(yōu)化模型中，以入爐混煤平均標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)最低為目標(biāo)函數(shù)，決策變量為各類電煤摻燒比例，約束條件以供熱約束、庫(kù)存約束和鍋爐約束的形式保證鍋爐燃燒生產(chǎn)的可靠性，通過(guò)求解混煤摻燒方案向上層優(yōu)化模型傳遞電煤消耗情況.

通過(guò)迭代上下層模型協(xié)同優(yōu)化目標(biāo)函數(shù)，避免單獨(dú)求解模型時(shí)陷入燃料成本局部最優(yōu)情況，探尋全局最優(yōu)燃料管控方案.

1.2.1 ? 電煤采購(gòu)庫(kù)存模型

1）目標(biāo)函數(shù). 本文在文獻(xiàn)[4]基礎(chǔ)上將電煤采購(gòu)總成本C分為廠外購(gòu)入成本C1和廠內(nèi)倉(cāng)儲(chǔ)成本C2，廠外購(gòu)入成本C1集中體現(xiàn)為入廠原煤?jiǎn)蝺r(jià)和電煤采購(gòu)量，廠內(nèi)倉(cāng)儲(chǔ)成本C2包括電煤庫(kù)存管理費(fèi)用和采購(gòu)資金積壓利息：

min C = C1 + C2 ? ? ? ? ? （4）

C1 = xi，m Ui pi，m ? ? ? ? （5）

C2 = （xi，m Ui φ + yi，m pi，m θ） ? ? ? ? （6）

式中：Nm為采購(gòu)周期內(nèi)周數(shù);Ns為可選電煤供應(yīng)商數(shù);整數(shù)型決策變量xi，m代表電廠第m周向第i個(gè)煤炭供應(yīng)商的購(gòu)買批次;Ui代表向第i個(gè)煤炭供應(yīng)商單一批次電煤購(gòu)買量;yi，m為跨月提前采購(gòu)電煤數(shù)量.

2）約束條件分為供貨約束，電廠庫(kù)存約束，分類庫(kù)存約束.

①供貨約束. 根據(jù)有無(wú)月前預(yù)訂合同，燃煤電廠向電煤供應(yīng)商訂購(gòu)的動(dòng)力煤可分為長(zhǎng)協(xié)煤和現(xiàn)貨煤. 長(zhǎng)協(xié)煤?jiǎn)卧虏少?gòu)量應(yīng)在合同的預(yù)定范圍內(nèi)，對(duì)應(yīng)供貨約束如下：

Smin，i≤xi，m Ui ≤Smax，i ? ? ? （7）

式中：xi，m Ui 代表電廠向第i個(gè)煤炭供應(yīng)商單月采購(gòu)的電煤量;Smin，i和Smax，i分別為煤炭供應(yīng)商在合同中簽訂的最低采購(gòu)量和最高采購(gòu)量.

②電廠庫(kù)存約束. 由于電廠存放煤炭場(chǎng)地有限，同時(shí)要保證最小連續(xù)發(fā)電天數(shù)，電煤庫(kù)存數(shù)量介于電廠庫(kù)存上下限之間，即

Pm =xi，m Ui ? ? ? ? （8）

St

m = St

m-1 + Pm - Ct

m ? ? ? （9）

St

min ≤ St

m ≤ St

max ? ? ? （10）

式中：Pm為電廠第m周的電煤采購(gòu)總量;Ct

m為電廠第m周的電煤消耗總量;St

m為電廠第m周的電煤庫(kù)存總量，由式（8）（9）計(jì)算得到;St

min和St

max分別為電廠總庫(kù)存下限和上限.

③分類庫(kù)存約束. 為保證電廠生產(chǎn)安全連續(xù)進(jìn)行，各類電煤庫(kù)存數(shù)量應(yīng)保持在安全范圍內(nèi)，即

Pj，m = xi，m Ui ? ? ? ?（11）

Sc

j，m = Sc

j，m-1+ Pj，m - Cc

j，m ? ? ?（12）

Sc

min≤ Sc

j，m≤ Sc

max ? ? ? （13）

式中：Cc

j，m為電廠第m周第j類電煤的消耗量;Sc

j，m為第m周第j類電煤的庫(kù)存量，由式（11）（12）計(jì)算得到;Sc

min和Sc

max分別代表電廠分類庫(kù)存下限和上限.

1.2.2 ? 混煤摻燒優(yōu)化模型

1）目標(biāo)函數(shù). 當(dāng)發(fā)電計(jì)劃確定時(shí)，機(jī)組發(fā)電所需的標(biāo)準(zhǔn)煤耗量確定，發(fā)電成本取決于入爐混煤平均標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)，因此混煤摻燒優(yōu)化模型以最低入爐混煤平均標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)為目標(biāo)函數(shù)[16]：

min p*

m = εj，m p*

j，m ? ? ? ? （14）

2）約束條件分為供熱約束，庫(kù)存約束，鍋爐約束.

①供熱約束. 鍋爐燃燒提供的熱量須滿足根據(jù)反平衡法計(jì)算機(jī)組發(fā)電的熱量需求，即

Qnet，m =εj，m Qnet，j ? ? ? ? （15）

Ct*

m = Gm gm /（1 - ηg） ? ? ? （16）

Ct

mQnet，m ηQ ≥Ct*

mQ*

net ? ? ? ? ? （17）

式中：Qnet，m為混煤平均低位發(fā)熱量，滿足混煤摻燒線性可加性;Ct*

m為第m周的發(fā)電標(biāo)準(zhǔn)煤耗量;Gm為第m周的發(fā)電量;gm為第m周的機(jī)組發(fā)電單位標(biāo)準(zhǔn)煤耗量;ηg為廠用電率;ηQ為鍋爐熱效率，受混煤低位發(fā)熱量影響，由混煤試燒試驗(yàn)得出.

②庫(kù)存約束. 機(jī)組發(fā)電使用電煤耗量受電廠當(dāng)周可用電煤限制，即

0≤Ct

m≤St

m-1+Pm ? ? ? ? （18）

Cc

j，m= Ct

mεj，m ? ? ?   （19）

0≤Cc

j，m≤Sc

j，m-1+ Pj，m ? ? ? ? （20）

③鍋爐約束. 為保證鍋爐安全生產(chǎn)，混煤煤質(zhì)應(yīng)在鍋爐燃燒設(shè)計(jì)范圍內(nèi)，具體如下：

Qnet，min≤Qnet，m ≤Qnet，max ? ? ? （21）

Sm = εj，m Sj ? ? ? ? （22）

Sm≤Smax ? ? ? （23）

Vdaf，m = εj，m Vdaf，j ? ? ? ? （24）

Vdaf，min≤Vdaf，m≤Vdaf，max ? ? ? ?（25）

Tamp，min≤Tamp，m≤Tamp，max ? ? ? ?（26）

式中：Sm和Vdaf，m分別為第m周的混煤全硫分和混煤干燥無(wú)灰基揮發(fā)份，均滿足混煤摻燒煤質(zhì)線性可加性;Tamp，m為第m周的混煤灰熔點(diǎn)溫度，由混煤試燒試驗(yàn)得出.

2 ? 灰狼優(yōu)化算法

灰狼優(yōu)化算法 （Grey Wolf Optimizer，GWO） 是Mirjalili等[19]在2014年受野外狼群的等級(jí)機(jī)制和狩獵方式啟迪而提出的群體智能優(yōu)化算法，在低維尋優(yōu)問(wèn)題中表現(xiàn)出較好的搜索能力和收斂速度，對(duì)混煤摻燒非凸優(yōu)化模型求解有不錯(cuò)的適應(yīng)性和魯棒性，但是面對(duì)電煤采購(gòu)庫(kù)存模型中高維的可行解空間和復(fù)雜時(shí)空連續(xù)性約束情況，GWO算法出現(xiàn)種群個(gè)體質(zhì)量低和早熟收斂問(wèn)題. 為此，本文從種群初始化和個(gè)體變異兩方面出發(fā)，分別使用混沌映射種群初始化算子和高斯變異算子改進(jìn)灰狼優(yōu)化算法（Improved Grey Wolf Optimizer，IGWO）.

2.1 ? 改進(jìn)灰狼優(yōu)化算法

2.1.1 ? 混沌映射種群初始化算子

電煤采購(gòu)庫(kù)存模型的求解質(zhì)量依賴于種群初始化的采購(gòu)庫(kù)存方案分布情況，優(yōu)秀的初始化種群采購(gòu)庫(kù)存方案盡可能在可行解空間呈現(xiàn)多樣化，對(duì)算法的全局搜索能力和收斂速度都有不錯(cuò)的提升. 相較于計(jì)算機(jī)偽隨機(jī)數(shù)生成，非線性的混沌映射具有初值敏感性和空間遍歷性的特點(diǎn)，可以提高初始化方案的多樣性.

本文采用Tent映射產(chǎn)生[0，1]隨機(jī)混沌序列，再通過(guò)等比放大映射為初始化電煤采購(gòu)庫(kù)存方案種群，Tent映射的表達(dá)式如下：

zk+1 = zk/β，zk∈（0，β]

（1 - zk）/（1 - β），zk∈（β，1] ? ? ?（27）

Xi，m = round（Xmin

i，m + zk（Xmax

i，m - Xmin

i，m）） ? ?（28）

式中：zk為初始隨機(jī)數(shù);β為混沌因子控制參數(shù)，取β = 0.3;Xmin

i，m、Xmax

i，m分別為第i種電煤第m周可購(gòu)訂單最大和最小值，取整后為電煤采購(gòu)模型初始整數(shù)決策變量Xi，m .

2.1.2 ? 高斯變異算子

針對(duì)GWO算法在電煤采購(gòu)庫(kù)存模型求解中早熟收斂的情況，本文受遺傳算法的啟發(fā)，在灰狼種群更新之后增加高斯變異算子以提高局部搜索能力. 高斯變異算子計(jì)算公式如下：

X*

i，m = round（Xi，m（1 + g）） ? ?（29）

式中：g為滿足均值0，方差1的高斯分布隨機(jī)數(shù);X*

i，m為高斯變異后的第m周第i個(gè)煤炭供應(yīng)商的采購(gòu)決策變量.

2.2 ? 算法流程

IGWO算法針對(duì)電煤采購(gòu)庫(kù)存模型的求解步驟為：

步驟1 ? 設(shè)置狼群中采購(gòu)庫(kù)存方案總數(shù)，最大迭代次數(shù)，根據(jù)約束條件確定電煤采購(gòu)庫(kù)存模型可行域;

步驟2 ? 根據(jù)式（27）（28）混沌映射初始化種群電煤采購(gòu)方案;

步驟3 ? 隨機(jī)選擇部分灰狼個(gè)體方案按式（29）進(jìn)行高斯變異，檢驗(yàn)種群內(nèi)灰狼個(gè)體方案是否滿足約束條件，對(duì)于越限采購(gòu)庫(kù)存方案映射到約束可行域內(nèi);

步驟4 ? 計(jì)算各方案采購(gòu)總成本并選擇最優(yōu)的3個(gè)采購(gòu)庫(kù)存方案作為灰狼領(lǐng)袖α、β和δ狼，其余灰狼ω按領(lǐng)袖狼信息更新位置;

步驟5 ? 重復(fù)步驟3、4，直至達(dá)到最大迭代次數(shù);

步驟6 ? 算法結(jié)束，輸出α狼對(duì)應(yīng)的電煤采購(gòu)庫(kù)存方案的采購(gòu)總成本并向混煤摻燒優(yōu)化模型傳遞可用電煤.

求解燃料成本雙層優(yōu)化模型算法流程如圖3所示，其中取相鄰迭代入爐平均標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)最大差值為雙層優(yōu)化模型收斂判據(jù)：

max（[p]*（k）

1 ? ? ? - [p]*（k-1）

1 ? ? ? ? ?，…，[p]*（k）

m ? ? ? - [p]*（k-1）

m ? ? ? ，

…，[p]*（k）

Nm ? ? ? - [p]*（k-1）

Nm ? ? ? ） < ε ? ? ? （30）

式中：[p]*（k）

m ? ? ?、[p]*（k-1）

m ? ? ? 分別為相鄰迭代的第m周入爐混煤平均標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià);ε為收斂判據(jù)，本文取ε = 0.001.

3 ? 算例分析

本文選擇某2 × 600 MW燃煤電廠某年6月至7月期間連續(xù)8周的燃煤管控決策算例分析. 選取6月至7月作為采購(gòu)周期是因?yàn)?月為迎夏用電高峰，全國(guó)發(fā)電量顯著增加，動(dòng)力煤價(jià)格同期上漲，是燃煤電廠重點(diǎn)關(guān)注的調(diào)整庫(kù)存時(shí)段. 該電廠庫(kù)存信息、鍋爐約束和發(fā)電計(jì)劃分別如表1、表2和表3所示;以中國(guó)2.3 × 104 kJ煤炭指數(shù)（Chinese Coal Index 5500，CCI5500）作為動(dòng)力煤價(jià)格參照，如圖4所示;具體的電煤價(jià)格和煤質(zhì)信息見(jiàn)表4.

輸入電廠歷史混煤摻燒方案，設(shè)定最大迭代次數(shù)為50，通過(guò)求解電煤采購(gòu)庫(kù)存模型獲得如圖5所示的電廠庫(kù)存曲線. 在采購(gòu)方案周期內(nèi)，燃煤電廠在動(dòng)力煤市場(chǎng)價(jià)格上漲時(shí)期以購(gòu)買長(zhǎng)協(xié)煤為主，穩(wěn)定發(fā)電成本;在現(xiàn)貨煤價(jià)格較低時(shí)，增加現(xiàn)貨煤采購(gòu)量以降低燃料成本. 另外，燃煤電廠通過(guò)提前采購(gòu)下周發(fā)電用煤，降低采購(gòu)成本，實(shí)現(xiàn)“價(jià)低拉庫(kù)存，價(jià)高保發(fā)電”的采購(gòu)優(yōu)化目標(biāo).

圖6為遺傳算法（Genetic Algorithm，GA），粒子群優(yōu)化算法（Particle Swarm Optimization，PSO）以及原始GWO算法與IGWO算法求解性能對(duì)比結(jié)果. 由圖6可知，4種求解算法均能在最大迭代次數(shù)內(nèi)收斂，相較于GA算法和PSO算法，2種GWO算法的采購(gòu)成本均能迅速下降到最優(yōu)解附近，然而原始GWO算法陷入局部最優(yōu)，收斂的解質(zhì)量最差. 對(duì)比4種算法的收斂速度和解的質(zhì)量，說(shuō)明IGWO算法在電煤采購(gòu)庫(kù)存模型求解擁有更優(yōu)的全局搜索能力和求解精度.

將上述使用電廠歷史混煤配比方案進(jìn)行電煤采購(gòu)庫(kù)存模型求解設(shè)置為算例1;將經(jīng)過(guò)一次混煤摻燒優(yōu)化后再對(duì)電煤采購(gòu)庫(kù)存模型求解設(shè)置為算例2，模擬單獨(dú)對(duì)電廠燃料管控環(huán)節(jié)優(yōu)化的情形;使用本文提出的燃料成本雙層優(yōu)化模型模擬電廠整體優(yōu)化燃料成本方案設(shè)置為算例3.

各算例各周入爐標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)求解結(jié)果如表5所示. 算例2平均單價(jià)較算例1降低35.57元/t，算例3較算例1降低44.27元/t，這是因?yàn)樗憷?受限于電廠可用電煤的數(shù)量，混煤配比優(yōu)化空間有限，而算例3使用燃料成本雙層優(yōu)化模型，電煤采購(gòu)方案與混煤配比方案協(xié)同優(yōu)化，可以實(shí)現(xiàn)鍋爐混煤配比方案主動(dòng)優(yōu)化效果. 其次，得益于提前購(gòu)入的低價(jià)長(zhǎng)協(xié)煤，第5周的平均入爐標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)成為各算例中7月單價(jià)最低周. 另外，算例3的第4周和第7周入爐標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)在算例1單價(jià)相同的情況下分別較第3周和第6周降低2.3元/t、1.2元/t，這是由于燃煤電廠事先簽訂的供應(yīng)合同要求必須使用一定數(shù)量的高價(jià)電煤，當(dāng)發(fā)電量升高時(shí)，耗煤量增高，燃煤電廠可以使用更多廉價(jià)動(dòng)力煤，混煤平均入爐標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)隨之降低.

算例1、算例2和算例3電煤采購(gòu)總成本分別為25 358.55萬(wàn)元、23 611.74萬(wàn)元和23 379.78萬(wàn)元，算例2和算例3分別比算例1降低 1 746.81萬(wàn)元、6.89%，1 978.77萬(wàn)元、7.80%，說(shuō)明對(duì)燃煤電廠燃料管控環(huán)節(jié)優(yōu)化可以減少電煤采購(gòu)成本. 本文提出的燃料成本雙層優(yōu)化模型對(duì)電煤采購(gòu)庫(kù)存管理和混煤配比的協(xié)同優(yōu)化處理可以更為明顯地降低電煤采購(gòu)成本，有效提高燃煤電廠的市場(chǎng)競(jìng)爭(zhēng)力.

4 ? 結(jié) ? 論

本文構(gòu)建一種燃煤電廠的燃料成本雙層優(yōu)化模型，并使用改進(jìn)灰狼優(yōu)化算法對(duì)其求解. 算例結(jié)果分析表明：

1）對(duì)于電煤采購(gòu)庫(kù)存模型中可行解空間復(fù)雜且間斷的情況，基于混沌初始化與高斯變異的IGWO算法的收斂速度和求解質(zhì)量?jī)?yōu)于GA、PSO與GWO，證明了IGWO算法在該模型中具有良好的收斂性和精度.

2）對(duì)燃煤電廠燃料管控環(huán)節(jié)優(yōu)化可以減少電煤采購(gòu)成本，而燃料成本雙層優(yōu)化模型兼顧電煤采購(gòu)方案和混煤配比方案的協(xié)同優(yōu)化，對(duì)比傳統(tǒng)燃料管控方案與單獨(dú)管控環(huán)節(jié)優(yōu)化方案，在降低采購(gòu)成本和入爐標(biāo)準(zhǔn)煤?jiǎn)蝺r(jià)方面更具優(yōu)越性.

參考文獻(xiàn)

[1] ? ?張艷芹，陳權(quán)寶. 國(guó)內(nèi)外煤炭市場(chǎng)價(jià)格動(dòng)態(tài)關(guān)系研究：基于電煤價(jià)格市場(chǎng)化改革前后變化的對(duì)比分析[J]. 價(jià)格理論與實(shí)踐，2015（1）：85—87.

ZHANG Y Q，CHEN Q B. Study on the dynamic relationship of coal market price at home and abroad[J]. Price：Theory & Practice，2015（1）：85—87. （In Chinese）

[2] ? ?廣東電力交易中心. 廣東電力市場(chǎng) 2019 年年度報(bào)告[R]. 廣州：廣東電力交易中心，2020.

Guangdong Power Exchange Center. Annual report of Guangdong power market in 2019[R]. Guangzhou：Guangdong Power Exchange Center，2020. （In Chinese）

[3] ? ?宋永華，包銘磊，丁一，等. 新電改下我國(guó)電力現(xiàn)貨市場(chǎng)建設(shè)關(guān)鍵要點(diǎn)綜述及相關(guān)建議[J]. 中國(guó)電機(jī)工程學(xué)報(bào)，2020，40（10）：3172—3187.

SONG Y H，BAO M L，DING Y，et al. Review of Chinese electricity spot market key issues and its suggestions under the new round of Chinese power system reform[J]. Proceedings of the CSEE，2020，40（10）：3172—3187. （In Chinese）

[4] ? ?黃雷. 火力發(fā)電廠燃煤全過(guò)程成本分析與優(yōu)化研究[D]. 北京：華北電力大學(xué)（北京），2018：16—24.

HUANG L. Cost analysis and optimization of the whole process coal burning thermal power plant[D]. Beijing：North China Electric Power University，2018：16—24. （In Chinese）

[5] ? ?馬姍姍. 發(fā)電企業(yè)燃煤庫(kù)存管理優(yōu)化研究[D]. 北京：華北電力大學(xué)（北京），2014：18—26.

MA S S. Research on the optimization of coal inventory of power plants[D]. Beijing：North China Electric Power University，2014：18—26.（In Chinese）

[6] ? ?張曦木. 火電廠燃煤配比與采購(gòu)計(jì)劃優(yōu)化方法與系統(tǒng)開(kāi)發(fā)[D]. 沈陽(yáng)：東北大學(xué)，2014：17—34.

ZHANG X M. Research on optimization and system design for coal blending and purchase planning in thermal power plants[D]. Shenyang：Northeastern University，2014：17—34. （In Chinese）

[7] ? ?張廣宏. 數(shù)字燃煤環(huán)境下發(fā)電企業(yè)燃煤庫(kù)存管理優(yōu)化研究[D]. 北京：華北電力大學(xué)（北京），2016：18—23.

ZHANG G H. Research on optimization of coal inventory management of power plants under digital coal[D]. Beijing：North China Electric Power University，2016：18—23. （In Chinese）

[8] ? ?吳橋. 現(xiàn)貨價(jià)格波動(dòng)下原材料最優(yōu)采購(gòu)決策研究[D]. 杭州：浙江大學(xué)，2012：33—44.

WU Q. Optimal procurement strategies of raw materials under spot price uncertainty[D]. Hangzhou：Zhejiang University，2012：33—44.（In Chinese）

[9] ? ?謝瑛. 煤炭?jī)r(jià)格波動(dòng)對(duì)火力發(fā)電廠經(jīng)營(yíng)利潤(rùn)影響風(fēng)險(xiǎn)分析[D]. 北京：華北電力大學(xué)（北京），2010：33—55.

XIE Y. The analysis of the efect of coal price fluctuation on the economic benerit of power plant[D]. Beijing：North China Electric Power University ，2010：33—55. （In Chinese）

[10] ?石佳星. 考慮動(dòng)力煤期貨的電煤庫(kù)存優(yōu)化模型研究[D]. 北京：華北電力大學(xué)（北京），2017：19—38.

SHI J X. Research on electricity-coal inventory optimization model considering steam coal futures[D]. Beijing：North China Electric Power University，2017：19—38. （In Chinese）

[11] ?楊甲甲，何洋，鄒波，等. 電力市場(chǎng)環(huán)境下燃煤電廠電煤庫(kù)存優(yōu)化的CVaR模型[J]. 電力系統(tǒng)自動(dòng)化，2014，38（4）：51—59.

YANG J J，HE Y，ZOU B，et al. A CVaR-based coal inventory optimization model for coal-fired power plants in electricity market environment[J]. Automation of Electric Power Systems，2014，38（4）：51—59. （In Chinese）

[12] ?黃海東，劉福國(guó)，國(guó)欽光，等. 配煤摻燒最小煤質(zhì)偏差模型及煤場(chǎng)存放優(yōu)化[J]. 熱能動(dòng)力工程，2020，35（5）：97—104.

HUANG H D，LIU F G，GUO Q G，et al. Modeling of minimum coal quality deviation for coal blending and optimization of coal yard storage[J]. Journal of Engineering for Thermal Energy and Power，2020，35（5）：97—104. （In Chinese）

[13] ?BAEK S H，PARK H Y，KO S H. The effect of the coal blending method in a coal fired boiler on carbon in ash and NOx emission[J]. Fuel，2014，128：62—70.

[14] ?GUERRAS L S，MARTíN M. Optimal gas treatment and coal blending for reduced emissions in power plants：a case study in Northwest Spain[J]. Energy，2019，169：739—749.

[15] ?NAWAZ Z，ALI U. Techno-economic evaluation of different operating scenarios for indigenous and imported coal blends and biomass co-firing on supercritical coal fired power plant performance[J]. Energy，2020，212：118721.

[16] ?李金晶，焦開(kāi)明，趙振寧，等. 火力發(fā)電廠配煤摻燒的燃煤成本模型[J]. 中國(guó)電力，2018，51（9）：15—19.

LI J J，JIAO K M，ZHAO Z N，et al. Model of coal cost for blending coal combustion in thermal power plant[J]. Electric Power，2018，51（9）：15—19. （In Chinese）

[17] ?姚偉，郝兵，劉家利，等. 配煤摻燒方式主要特點(diǎn)及燃煤適應(yīng)性分析[J]. 中國(guó)電力，2018，51（9）：20—27.

YAO W，HAO B，LIU J L，et al. Main characteristics of coal blending method and adaptability analysis for blended coal[J]. Electric Power，2018，51（9）：20—27. （In Chinese）

[18] ?劉貴鋒，韓雪松. 煙煤鍋爐摻燒褐煤對(duì)機(jī)組安全經(jīng)濟(jì)性影響的研究[J]. 中國(guó)電力，2018，51（9）：46—52.

LIU G F，HAN X S. Study on the impacts of blending lignite in bituminous coal boiler on unit security and economy[J]. Electric Power，2018，51（9）：46—52.（In Chinese）

[19] ?MIRJALILI S，MIRJALILI S M，LEWIS A. Grey wolf optimizer[J]. Advances in Engineering Software，2014，69（3）：46—61.

收稿日期：2021-06-23

基金項(xiàng)目：國(guó)家自然科學(xué)基金資助項(xiàng)目（51937005），National Natural Science Foundation of China（51937005）

作者簡(jiǎn)介：廖志偉（1973—），男，廣西平樂(lè)人，華南理工大學(xué)副教授

通信聯(lián)系人，E-mail：epliao@scut.edu.cn