0 引言

    電力系統(tǒng)隨機生產(chǎn)模擬是在考慮機組的隨機故障及電力負荷隨機性的情況下,，通過模擬和優(yōu)化發(fā)電機組的生產(chǎn)運行,，計算系統(tǒng)在最優(yōu)運行方式下不同發(fā)電機組的發(fā)電量、可靠性指標以及生產(chǎn)成本指標等的算法^[1],。當前,，隨著資源稟賦,、環(huán)境污染和氣候問題的日益嚴峻，以及清潔能源的大規(guī)模發(fā)展,，提高了對時序性,、互動性、隨機性電力系統(tǒng)生產(chǎn)模擬技術的需求,。目前,，已有文獻對其進行了研究，文獻[2]圍繞電力系統(tǒng)生產(chǎn)模擬中可再生能源的不確定性問題,，通過多場景隨機規(guī)劃理論反映可再生能源的不確定性,，建立了隨機規(guī)劃數(shù)學模型模擬系統(tǒng)運行，并采用Benders分解法進行求解,。文獻[3],、文獻[4]基于序列運算，探索了適用于綜合資源規(guī)劃和電力市場的隨機生產(chǎn)模擬方法,，提出了基于序列運算理論的隨機生產(chǎn)模擬,。文獻[5]將需求側(cè)的能效電廠納入生產(chǎn)模擬的研究范圍，為電力系統(tǒng)雙向互動的生產(chǎn)模擬提供了方向,。文獻[6]根據(jù)智能工程理論,，用Agent模擬的方法描述機組和能效電廠，實現(xiàn)了能效電廠運行的電力系統(tǒng)生產(chǎn)Agent模擬,，并證明了有效性,。

    本文將在以往研究的基礎上，將電力系統(tǒng)分為發(fā)-輸-配-用幾個環(huán)節(jié),，基于復雜適應系統(tǒng)原理,，通過供需雙方資源的調(diào)度過程反映電力系統(tǒng)運行，以可靠性為約束條件,，考慮發(fā)電元件的隨機性,，在可計及多態(tài)機組的情況下，建立包含傳統(tǒng)能源和可再生能源的電力系統(tǒng)時序生產(chǎn)模擬模型,，并對模型進行求解,。通過模型計算結(jié)果與電力發(fā)展的實際數(shù)據(jù)進行對比，證明其在電力發(fā)展規(guī)劃和分析研究方面的支撐性作用,。

1 復雜適應系統(tǒng)理論

    復雜適應系統(tǒng)（Complex Adaptive System,，CAS）是指由兩個以上按照一定規(guī)則或模式進行相互作用的行為主體所組成的復雜動態(tài)系統(tǒng)^[7]，具有以下特點：(1)復雜適應系統(tǒng)強調(diào)主體(Agent)的主動性,，即Agent有各自的知識和目標等屬性,，由于屬性不同，在相同環(huán)境下不同Agent可能表現(xiàn)出不同的行為方式,；(2)復雜適應系統(tǒng)認為Agent之間,，及Agent與環(huán)境之間存在交互,、影響、反饋等作用,，這些作用的綜合結(jié)果是系統(tǒng)進化,、演變的動力和原因；(3)不同主體通過特定規(guī)則相互聯(lián)系形成系統(tǒng),，而系統(tǒng)又在更宏觀的層面作為“Agent”與其他“Agent”交互成新系統(tǒng)^[8]。一個基于智能工程的簡單復雜適應系統(tǒng)模型中包含外部環(huán)境,、上層系統(tǒng),、協(xié)調(diào)系統(tǒng)和下層系統(tǒng)，這里不同的系統(tǒng)層級可以是單一Agent,，也可以由多個Agent組成的系統(tǒng),。圖1展示了一個各部分都由單一Agent組成的復雜適應系統(tǒng)。

    圖1中,，上層系統(tǒng)和下層系統(tǒng)依據(jù)自身的目標,，按照知識庫和規(guī)則根據(jù)外部環(huán)境的變化進行調(diào)整；協(xié)調(diào)層感知上下兩層的狀態(tài),，若上下兩層中的單個目標超出允許的閾值范圍,，協(xié)調(diào)層感知狀態(tài)的變化并協(xié)調(diào)上下兩層進行調(diào)整，并將協(xié)調(diào)后的狀態(tài)返回給上下兩層,，使其再次進行調(diào)整直至整體達到最優(yōu)狀態(tài)并將結(jié)果返回給外部環(huán)境,，保證整體的向目標或最優(yōu)方向進化。

2 復雜適應系統(tǒng)生產(chǎn)模擬模型

2.1 系統(tǒng)模型

    電力系統(tǒng)運行的過程實際上是供給側(cè)電源系統(tǒng)與需求側(cè)電力負荷的供需匹配過程,。要實現(xiàn)這一過程需要發(fā)電企業(yè),、電網(wǎng)企業(yè)、電力市場和用戶四個環(huán)節(jié)的協(xié)調(diào)配合,。因此,，基于復雜適應系統(tǒng)的電力系統(tǒng)生產(chǎn)模擬模型由4個模塊模擬電力系統(tǒng)功能：發(fā)電企業(yè)模塊、電網(wǎng)企業(yè)模塊,、電力市場模塊和用戶需求模塊,，其結(jié)構(gòu)如圖2。

2.1.1 用戶需求模塊

    用戶需求模塊作為其他模塊需要適應的外部環(huán)境功能存在,，模擬了電力用戶的用電特征和規(guī)律,。其核心功能是根據(jù)生產(chǎn)模擬需要為其他模塊提供電力負荷需求信息，用戶需求模塊的變化是系統(tǒng)發(fā)展的動力以及其他模塊產(chǎn)生生產(chǎn)行為的目的,。

    本文的生產(chǎn)模擬以小時為時間尺度,，模擬目標年全年8 760小時為模擬周期。根據(jù)模型的功能要求,，電力負荷變化的特征計算公式如下：

其中,，E為預測電量,；P為每個時間點的負荷，P=α×P_max,，P_max為系統(tǒng)最大負荷,，α為負荷的標幺值。將最大負荷和負荷特征分開,，可以方便對不同周期或時間跨度的電力系統(tǒng)進行分析,。

2.1.2 電力生產(chǎn)模塊

    電力生產(chǎn)模塊模擬不同的發(fā)電機組，包括：火電,、水電,、天然氣、風電,、光伏和核電,，其主要功能為計算和調(diào)整不同機組的發(fā)電成本并將結(jié)果發(fā)送至其他模塊；根據(jù)可再生能源特征信息生成其各小時的出力狀態(tài),；最后根據(jù)發(fā)電計劃進行生產(chǎn)模擬,，計算污染物排放量、可再生能源棄電率等指標,。具體的機組報價計算和調(diào)整的原則如下：

    將發(fā)電機組i每小時的成本函數(shù)表示為二次函數(shù),，即：

    (1)風電

    風電的出力呈現(xiàn)隨機性，本文采用雙參數(shù)的威布爾分布來表述風速的概率密度,，而風機的出力主要由風輪機輪轂高度處的風速來決定,，其具體關系為：

2.1.3 電網(wǎng)調(diào)度模塊

    電網(wǎng)調(diào)度模塊以最小碳排放或最小運行成本為目標，根據(jù)電力市場模塊計算的出清價格以及電力生產(chǎn)模塊共享的機組裝機規(guī)模,、報價,，根據(jù)負荷信息，統(tǒng)計電力供需情況,，安排發(fā)電計劃,；最后根據(jù)電網(wǎng)約束條件，在時間周期T內(nèi)計算電力不足概率（Loss of Load Probability,，LOLP）和電力不足期望值（Expected Energy Not Served,，EENS），計算公式如下：

2.1.4 電力市場模塊

    電力市場模塊的功能為計算電力出清價格及模擬電力輔助服務市場生成機組停運概率矩陣,。在不考慮網(wǎng)絡約束的平衡交易市場中,，用戶以購電成本C最小化為目標，即

求導可得拉格朗日乘子λ為用戶購電成本的最小化條件為將其與功率平衡方程聯(lián)立求解可以得到系統(tǒng)的市場出清價格以及個發(fā)電廠商的發(fā)電出力,。

    電力輔助服務市場方面,，考慮常規(guī)機組的強迫停運，在模擬時先對其的運行狀態(tài)進行抽樣,，然后通過多次模擬系統(tǒng)運行來對隨機強迫停運進行近似模擬,，其目標函數(shù)如下：

    上述模型的計算流程如圖3,。

3 算例分析

3.1 數(shù)據(jù)和相關參數(shù)設定

    模型通過Java和Python編程實現(xiàn)。為了簡化計算,，生產(chǎn)模擬輸入的數(shù)據(jù)中,，最大負荷量、負荷需求,、可再生能源出力特征,、不同機組的初始規(guī)模等需要設定初始值，成本參數(shù),、補貼參數(shù),、電力生產(chǎn)的排放系數(shù)、機組運行小時數(shù)等均考慮為常量,。

3.1.1 風光出力數(shù)據(jù)

    本文從統(tǒng)計學的角度出發(fā)，收集不同地區(qū)風,、光能統(tǒng)計和歷史數(shù)據(jù),，建立風光出力模型。根據(jù)數(shù)據(jù)統(tǒng)計結(jié)果,，以坐標軸的橫坐標為可再生能源的出力水平,，縱坐標為出力水平相應概率，得到的風能/太陽能模擬月度出力特性如圖4,。

3.1.2 生產(chǎn)模擬數(shù)據(jù)

    模型考慮的電源機組類型包括：水電,、煤電、天然氣,、核電,、風電和太陽能。根據(jù)2016年中國的電力結(jié)構(gòu),，其凈發(fā)電量為59 747億千瓦時,，總發(fā)電裝機容量為165 051萬千瓦，其中,，煤電94 624萬千瓦,、氣電7 011萬千瓦、水電33 207萬千瓦,、風電14 747萬千瓦,、核電3 364萬千瓦、太陽能7 631萬千瓦,。

3.2 計算結(jié)果

    在2016年裝機結(jié)構(gòu)下,，不同類型機組的年發(fā)電量、年利用小時數(shù)模擬與實際值的對比結(jié)果,，以及年運行費用結(jié)果如表1,。生產(chǎn)模擬計算的失負荷小時數(shù),、失負荷率、棄電量,、棄電率等可靠性指標與實際情況的對比結(jié)果如表2,。 

    2016年的電力裝機結(jié)構(gòu)下進行的8 760小時的生產(chǎn)模擬并在結(jié)果中分別選取冬季和夏天的典型周模擬結(jié)果顯示如圖5。

3.3 結(jié)果分析

    模擬結(jié)果與實際年總發(fā)電量的相對誤差約-5.05%,，火電,、水電、天然氣,、風電,、太陽能和核電的相對誤差分別為-16.47%、-17.36%,、46.23%%,、41.64%、65.51%,、-4.14%,。同時，根據(jù)《中國電力行業(yè)年度發(fā)展報告2017》來看,，2016年全國戶均停電時間為17.11小時,，供電可靠率為99.81%；而生產(chǎn)模擬結(jié)果中的年度停電時間為20.23小時,，供電可靠率為99.78%,。

    從生產(chǎn)模擬結(jié)果可知,，煤電機組在當前的系統(tǒng)中仍然占據(jù)主導地位,，煤電出力水平占整體的50%以上,。可再生能源棄電很大程度上是由于在負荷低谷時段可再生能源出力激增,，而其他常規(guī)電源出力已經(jīng)壓低至最小技術出力,，無法給可再生能源提供足夠的消納空間。而失負荷的原因是由于在負荷高峰時,，可再生能源出力不足,，即使全部常規(guī)電源出力均提高至最大依然無法滿足電力需求。

    模擬結(jié)果反映了當前電源結(jié)構(gòu)的幾個特征：(1)煤電是當前電力結(jié)構(gòu)中的主要電源出力形式,；(2)電力需求的峰谷差較大,，系統(tǒng)靈活性不足，難以單純的依靠調(diào)整電源出力規(guī)模滿足電力需求,；(3)傳統(tǒng)電源固定容量占比較大,，在負荷需求較低時會擠占可再生能源的消納空間，導致可再生能源棄電。

4 結(jié)論

    本文在復雜適應系統(tǒng)理論的基礎上建立了電力系統(tǒng)生產(chǎn)模擬模型,，并在2016年電源結(jié)構(gòu)的條件下驗證模型的電力系統(tǒng)發(fā)電量與實際發(fā)電量的相對誤差約4.71%,。根據(jù)模型求解結(jié)果顯示，其可以在給定電源結(jié)構(gòu)的情況下對滿足電力需求的不同機組類型發(fā)電量,、污染物排放量,、電力系統(tǒng)經(jīng)濟性、供電可靠性,、可再生能源利用效率等進行分析和預測,。

    另外，本研究仍需深入展開如下工作：細化用戶需求和電網(wǎng)企業(yè)模塊,，考慮電力負荷多樣性的特征,，引入電網(wǎng)安全運行特征分析；深化電力市場模塊,，豐富電力市場交易機制,；而電力企業(yè)模塊中不同機組的策略選擇單一，多場景方式簡化模擬可再生能源發(fā)電的合理性也需要進一步論證,。

參考文獻

[1] 王錫凡,，王秀麗.隨機生產(chǎn)模擬及其應用[J].電力系統(tǒng)自動化，2003,，27(8)：10-15.

[2] 姚力，王秀麗,，肖漢,，等.基于多場景隨機規(guī)劃的電力系統(tǒng)生產(chǎn)模擬[J].電力建設，2016,，37(12)：74-81.

[3] 康重慶,，夏清，相年德,，等.隨機生產(chǎn)模擬的序列化分析[J].中國電機工程學報,，2002，22(4)：8-12.

[4] 夏清,，相年德,，白利超，等.序列運算理論及其應用[J].電力系統(tǒng)自動化,，2002,，26(17)：6-12.

[5] 周景宏，胡兆光,，田建偉,，等.含能效電廠的力系統(tǒng)生產(chǎn)模擬[J].電力系統(tǒng)自動化，2010(18)：27-31.

[6] 段煒,，胡兆光,，姚明濤,，等.考慮需方響應資源的電力系統(tǒng)生產(chǎn)模擬[J].電網(wǎng)技術，2014,，38(6)：1523-1529.

[7] HOLLAND H J.Complex adaptive system[M].Addison Wesley,，1995.

[8] YIN R.Discovering the future of the case study method in evaluation research[J].Evaluation Practice，1994,，15(3)：283-290.

[9] 劉東冉,，陳樹勇，馬敏,，等.光伏發(fā)電系統(tǒng)模型綜述[J].電網(wǎng)技術,，2011，35(8)：47-52.

作者信息:

黃  何,，梁大鵬

（哈爾濱工業(yè)大學 經(jīng)濟與管理學院,，黑龍江 哈爾濱150001）