電力-天然氣綜合能源系統(tǒng)建模與仿真

常曉勇，王玉婷，陳福鋒，李玉平，徐程驥

（1.國(guó)電南京自動(dòng)化股份有限公司，南京 210032；2.南京國(guó)電南自電網(wǎng)自動(dòng)化有限公司，南京 211153）

0 引言

能源問(wèn)題和環(huán)境問(wèn)題是人類(lèi)21世紀(jì)面臨的兩大難題，制約著人類(lèi)社會(huì)的可持續(xù)發(fā)展。近年來(lái)，雖然我國(guó)的智能電網(wǎng)建設(shè)不斷取得突破性進(jìn)展，但是人們逐漸意識(shí)到，單獨(dú)依靠智能電網(wǎng)仍然無(wú)法很好地解決大規(guī)模可再生能源的消納問(wèn)題[1]。此外，隨著社會(huì)能源消耗總量的不斷增長(zhǎng)，采用傳統(tǒng)的各能源系統(tǒng)單獨(dú)規(guī)劃、建設(shè)、運(yùn)行的發(fā)展模式無(wú)法實(shí)現(xiàn)不同能源形式之間的優(yōu)勢(shì)互補(bǔ)，影響能源的綜合利用效率，不利于降低能源整體成本[2]。為此，提出了綜合能源系統(tǒng)的概念。

綜合能源系統(tǒng)是指在規(guī)劃、建設(shè)和運(yùn)行等過(guò)程中，通過(guò)對(duì)能源的產(chǎn)生、傳輸與分配（能源供應(yīng)網(wǎng)絡(luò)）、轉(zhuǎn)換、存儲(chǔ)和消費(fèi)等環(huán)節(jié)實(shí)施有機(jī)協(xié)調(diào)與優(yōu)化，進(jìn)而形成的能源產(chǎn)供銷(xiāo)一體化系統(tǒng)，是能源互聯(lián)網(wǎng)的重要載體[3]。綜合能源系統(tǒng)主要由供電、供氣、供冷/熱、電氣化交通等不同能源系統(tǒng)組成，利用先進(jìn)的信息網(wǎng)絡(luò)和通信技術(shù)，實(shí)現(xiàn)多種能源系統(tǒng)的深度融合和緊密互動(dòng)，促進(jìn)能源及設(shè)備的高效利用，提高能源系統(tǒng)安全性和可靠性。通常，按照地理范圍可把綜合能源系統(tǒng)分為跨區(qū)綜合能源系統(tǒng)、區(qū)域綜合能源系統(tǒng)和終端能源系統(tǒng)[4]。

目前，國(guó)內(nèi)外對(duì)于綜合能源系統(tǒng)的研究尚處于起步階段，當(dāng)前的關(guān)注焦點(diǎn)主要集中在多能源系統(tǒng)的協(xié)同規(guī)劃、穩(wěn)態(tài)分析、潮流計(jì)算、經(jīng)濟(jì)調(diào)度和運(yùn)行優(yōu)化等方面。文獻(xiàn)[5]研究了基于冷/熱/電三聯(lián)供系統(tǒng)的多區(qū)域綜合能源系統(tǒng)的規(guī)劃問(wèn)題，實(shí)現(xiàn)了系統(tǒng)中各類(lèi)能源的優(yōu)化利用。文獻(xiàn)[6]對(duì)電力-天然氣區(qū)域綜合能源系統(tǒng)進(jìn)行了穩(wěn)態(tài)建模和求解，研究了天然氣系統(tǒng)網(wǎng)絡(luò)狀態(tài)改變對(duì)區(qū)域綜合能源系統(tǒng)的影響。文獻(xiàn)[7]提出了一種適用于電、熱、氣的綜合能源系統(tǒng)擴(kuò)展Newton-Raphson多能潮流計(jì)算方法，可以在系統(tǒng)不同運(yùn)行場(chǎng)景下對(duì)多能潮流進(jìn)行計(jì)算和分析。文獻(xiàn)[8]提出了一種計(jì)及風(fēng)、光出力和電、熱負(fù)荷不確定性的綜合能源系統(tǒng)經(jīng)濟(jì)調(diào)度模型，利用CVaR理論研究了不確定性因素對(duì)系統(tǒng)調(diào)度產(chǎn)生的影響。文獻(xiàn)[9]以綜合能源系統(tǒng)的供暖體系為研究對(duì)象，利用供應(yīng)側(cè)和需求側(cè)的經(jīng)濟(jì)性模型，實(shí)現(xiàn)了供應(yīng)側(cè)電廠運(yùn)行效率和需求側(cè)居民用電成本的雙目標(biāo)優(yōu)化。

除了上述研究方向以外，綜合能源系統(tǒng)中各能源網(wǎng)絡(luò)相互之間的動(dòng)態(tài)影響問(wèn)題也是一個(gè)需要研究的重要課題。然而，目前國(guó)內(nèi)外針對(duì)該課題的研究成果還不是很多。文獻(xiàn)[10]通過(guò)建立天然氣網(wǎng)絡(luò)的動(dòng)態(tài)仿真模型，研究了天然氣網(wǎng)絡(luò)和電力網(wǎng)絡(luò)之間的相互影響，然而文中主要側(cè)重天然氣網(wǎng)絡(luò)的研究，對(duì)電力網(wǎng)絡(luò)進(jìn)行了簡(jiǎn)化，并且只考慮了能量通過(guò)燃?xì)廨啓C(jī)發(fā)電系統(tǒng)從天然氣網(wǎng)絡(luò)向電力網(wǎng)絡(luò)流動(dòng)的情況，并未考慮能量從電力網(wǎng)絡(luò)向天然氣網(wǎng)絡(luò)流動(dòng)的情況。

隨著燃?xì)廨啓C(jī)發(fā)電技術(shù)和P2G（電轉(zhuǎn)氣）技術(shù)的快速發(fā)展和應(yīng)用，天然氣系統(tǒng)在能源系統(tǒng)中的重要性越來(lái)越顯著，天然氣網(wǎng)絡(luò)與電力網(wǎng)絡(luò)之間的耦合也越來(lái)越緊密。因此，本文選取電-氣互聯(lián)的終端能源系統(tǒng)為研究對(duì)象，基于MATLAB/Simulink仿真軟件搭建了包含電力網(wǎng)絡(luò)、天然氣網(wǎng)絡(luò)和耦合環(huán)節(jié)（微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)和P2G系統(tǒng)）的終端能源系統(tǒng)仿真平臺(tái)，在此基礎(chǔ)上研究了微電網(wǎng)動(dòng)態(tài)過(guò)程對(duì)天然氣供應(yīng)網(wǎng)絡(luò)的影響。

1 電-氣互聯(lián)終端能源系統(tǒng)建模

以一個(gè)簡(jiǎn)化的終端能源系統(tǒng)為研究對(duì)象，系統(tǒng)拓?fù)淙鐖D1所示。該系統(tǒng)主要由一個(gè)交流微電網(wǎng)和一個(gè)天然氣供應(yīng)網(wǎng)絡(luò)組成，二者之間通過(guò)一個(gè)微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)和一個(gè)P2G系統(tǒng)相互耦合。其中，微電網(wǎng)由光伏發(fā)電系統(tǒng)、儲(chǔ)能系統(tǒng)和負(fù)載等部分組成，通過(guò)PCC（公共連接點(diǎn)）的升壓變壓器接入上一級(jí)配電網(wǎng)，具有孤島運(yùn)行和并網(wǎng)運(yùn)行2種不同模式；天然氣供應(yīng)網(wǎng)絡(luò)為一個(gè)簡(jiǎn)單的低壓供氣網(wǎng)絡(luò)，通過(guò)一個(gè)調(diào)壓閥接入上一級(jí)天然氣網(wǎng)，調(diào)壓閥同時(shí)負(fù)責(zé)對(duì)網(wǎng)絡(luò)內(nèi)的管道壓力進(jìn)行動(dòng)態(tài)調(diào)節(jié)，維持管道壓力在規(guī)定范圍內(nèi)。

圖1 電-氣互聯(lián)終端能源系統(tǒng)拓?fù)?/p>

圖1所示系統(tǒng)中，微型燃?xì)廨啓C(jī)作為天然氣網(wǎng)絡(luò)中的負(fù)荷，利用天然氣燃燒做功，驅(qū)動(dòng)同步發(fā)電機(jī)發(fā)電，同時(shí)作為微電網(wǎng)中的電源為電力負(fù)載供電；P2G系統(tǒng)作為微電網(wǎng)中的負(fù)載，根據(jù)電化學(xué)反應(yīng)和催化反應(yīng)，將電能轉(zhuǎn)化為氫氣和甲烷，同時(shí)作為天然氣網(wǎng)絡(luò)中的氣源，將生成的氣體注入天然氣網(wǎng)絡(luò)進(jìn)行存儲(chǔ)或供其他負(fù)載使用。

1.1 微電網(wǎng)建模

針對(duì)微電網(wǎng)中的光伏發(fā)電系統(tǒng)，首先采用文獻(xiàn)[11]中的方法建立光伏組件的數(shù)學(xué)模型，其次將光伏組件進(jìn)行串聯(lián)和并聯(lián)后組成光伏陣列，最后通過(guò)電壓源型DC-AC變流器接入微電網(wǎng)。通過(guò)對(duì)DC-AC變流器直流側(cè)母線電壓進(jìn)行閉環(huán)調(diào)節(jié)，實(shí)現(xiàn)光伏發(fā)電的MPPT（最大功率點(diǎn)跟蹤）控

制[12]。

針對(duì)微電網(wǎng)中的電池儲(chǔ)能系統(tǒng)，采用等效電路的建模方法搭建了電池單體模型[13]。電池陣列通過(guò)雙向DC-AC變流器接入微電網(wǎng)。當(dāng)微電網(wǎng)運(yùn)行于并網(wǎng)模式時(shí)，電池儲(chǔ)能系統(tǒng)采用PQ控制，根據(jù)微電網(wǎng)中央控制單元下發(fā)的功率指令，進(jìn)行有功功率和無(wú)功功率調(diào)節(jié)；當(dāng)微電網(wǎng)運(yùn)行于孤島模式時(shí)，電池儲(chǔ)能系統(tǒng)配合微電網(wǎng)中的微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)和光伏發(fā)電系統(tǒng)，參與微電網(wǎng)的電壓和頻率調(diào)節(jié)。

1.2 天然氣網(wǎng)絡(luò)建模

采用文獻(xiàn)[10]中的建模方法，將天然氣流假設(shè)為一維流體，利用質(zhì)量守恒定律和牛頓第二運(yùn)動(dòng)定律，得到表示動(dòng)態(tài)天然氣流的連續(xù)性方程和動(dòng)量方程。采用有限元近似的思想，將連續(xù)性方程和動(dòng)量方程由偏微分方程轉(zhuǎn)化為常微分方程：

式中：p為壓力；A為管道的橫截面積；c為天然氣中聲音的傳播速度；M為以質(zhì)量計(jì)量的天然氣流量；D為管道直徑；f為摩擦系數(shù)；L為管道的長(zhǎng)度。

天然氣網(wǎng)絡(luò)中調(diào)壓閥的功能是維持其出口壓力在設(shè)定值，為了模擬調(diào)壓閥的動(dòng)態(tài)特性，采用簡(jiǎn)化模型為：

1.3 耦合環(huán)節(jié)建模

微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)中主要包括微型燃?xì)廨啓C(jī)、發(fā)電機(jī)和變流器。微型燃?xì)廨啓C(jī)主要由壓氣機(jī)、燃燒室和透平組成，采用單軸型結(jié)構(gòu)[14]。根據(jù)微型燃?xì)廨啓C(jī)各部分的數(shù)學(xué)模型建立微型燃?xì)廨啓C(jī)的整體模型，微型燃?xì)廨啓C(jī)驅(qū)動(dòng)永磁同步發(fā)電機(jī)發(fā)電，然后經(jīng)過(guò)整流器、逆變器及濾波器并入微電網(wǎng)。微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)拓?fù)淙鐖D2所示。其中永磁同步發(fā)電機(jī)與微型燃?xì)廨啓C(jī)同軸，由于其轉(zhuǎn)速較高，需要首先經(jīng)過(guò)AC-DC變流器將高頻交流轉(zhuǎn)換為直流，再經(jīng)過(guò)DC-AC變流器將直流轉(zhuǎn)換為工頻交流。

如圖2所示的控制器中，AC-DC變流器采用轉(zhuǎn)速控制：ωg為測(cè)量到的發(fā)電機(jī)轉(zhuǎn)速，與額定轉(zhuǎn)速ωn比較，對(duì)誤差進(jìn)行PI控制，得到q軸參考信號(hào)isqref，d軸參考信號(hào)isdref=0；根據(jù)永磁同步發(fā)電機(jī)的三相電壓usabc、三相電流isabc和相角θs，進(jìn)行dq變換得到d軸電流isd和q軸電流isq，分別與isdref和isqref比較，對(duì)誤差進(jìn)行PI控制，輸出AC-DC變流器的SPWM控制信號(hào)。DC-AC變流器采用功率控制：根據(jù)DC-AC變流器輸出三相電壓uabc、三相電流iabc和相角θ，進(jìn)行dq變換得到電壓的d軸分量ud和q軸分量uq，以及電流的d軸分量 id和 q 軸分量 iq； 根據(jù) ud， uq， id， iq計(jì)算得到DC-AC變流器輸出有功功率Pgrid和無(wú)功功率Qgrid，分別與參考值進(jìn)行比較，經(jīng)過(guò)PI控制，得到d軸參考信號(hào)idref和q軸參考信號(hào) iqref；idref和iqref分別與id和iq進(jìn)行比較，經(jīng)過(guò)PI控制，輸出DC-AC變流器的SPWM控制信號(hào)。

圖2 微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)拓?fù)?/p>

P2G技術(shù)是指利用電能將H2O和CO2轉(zhuǎn)化為H2或CH4的過(guò)程。P2G的實(shí)現(xiàn)過(guò)程主要包括電解水和氫氣甲烷化2個(gè)步驟，如圖3所示。其中，微電網(wǎng)中的電能首先需要通過(guò)AC-DC轉(zhuǎn)換為直流，其次輸入P2G系統(tǒng)。

圖3 P2G系統(tǒng)拓?fù)?/p>

對(duì)于水電解環(huán)節(jié)，針對(duì)目前小容量裝置中普遍采用的PEM（聚合物電解質(zhì)膜）水電解技術(shù)，采用一種基于等效模型的建模方法，根據(jù)實(shí)驗(yàn)數(shù)據(jù)得到裝置中單體的I-V外特性表達(dá)式和H2產(chǎn)生速率表達(dá)式，從而得到其數(shù)學(xué)模型[15]。模型單體的外特性和制氫速率如圖4所示。電解水裝置根據(jù)其容量大小由若干個(gè)單體共同組成。

圖4 聚合物電解質(zhì)電解水模型單體特性曲線

目前關(guān)于甲烷化反應(yīng)裝置動(dòng)態(tài)運(yùn)行過(guò)程的研究成果較少，因此在建模過(guò)程中忽略了這部分的動(dòng)態(tài)過(guò)程。

2 電-氣互聯(lián)終端能源系統(tǒng)仿真

為了對(duì)文中所搭建的仿真模型的動(dòng)態(tài)特性進(jìn)行驗(yàn)證，在MATLAB/Simulink仿真平臺(tái)中按照?qǐng)D1所示系統(tǒng)搭建了仿真模型，設(shè)計(jì)了不同算例，分別進(jìn)行仿真和分析。

2.1 算例1

當(dāng)微電網(wǎng)中分布式光伏發(fā)電系統(tǒng)的輸出功率不足時(shí)，微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)投入運(yùn)行，采用功率控制模式進(jìn)行發(fā)電，此時(shí)P2G系統(tǒng)處于退出運(yùn)行狀態(tài)。研究了微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)輸出有功功率波動(dòng)時(shí)對(duì)天然氣網(wǎng)絡(luò)的影響，仿真波形如圖5所示。其中，當(dāng)t=20 s時(shí)，微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)輸出有功功率由額定值下降50%。根據(jù)圖5中的仿真波形可知，當(dāng)微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)輸出的有功功率發(fā)生波動(dòng)時(shí)，其需要的天然氣量也隨之波動(dòng)，由于天然氣動(dòng)態(tài)流動(dòng)是一個(gè)具有很大慣性時(shí)間常數(shù)的長(zhǎng)時(shí)間過(guò)程，因此天然氣網(wǎng)絡(luò)中各節(jié)點(diǎn)的壓力值會(huì)受到微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)進(jìn)氣口流量波動(dòng)的影響，產(chǎn)生相應(yīng)的波動(dòng)，并且這種波動(dòng)現(xiàn)象會(huì)由節(jié)點(diǎn)2開(kāi)始，由近及遠(yuǎn)依次進(jìn)行傳播。

圖5 算例1仿真波形

2.2 算例2

當(dāng)微電網(wǎng)中分布式光伏發(fā)電系統(tǒng)輸出有功充足時(shí)，微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)退出，P2G系統(tǒng)投入運(yùn)行。研究了P2G系統(tǒng)的動(dòng)態(tài)過(guò)程對(duì)天然氣網(wǎng)絡(luò)的影響，仿真波形如圖6所示。其中，當(dāng)t=20 s時(shí)，P2G系統(tǒng)的輸出由額定值下降50%。由圖6的仿真結(jié)果可知，與微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)相似，P2G系統(tǒng)的動(dòng)態(tài)過(guò)程同樣會(huì)對(duì)天然氣網(wǎng)絡(luò)的氣體壓力產(chǎn)生影響。圖6中，由于P2G系統(tǒng)輸出的氣體流量突然減少，導(dǎo)致其接入節(jié)點(diǎn)（節(jié)點(diǎn)3）的氣壓突然降低，由于天然氣動(dòng)態(tài)流動(dòng)的大慣性時(shí)間常數(shù)特性，與節(jié)點(diǎn)3相鄰的節(jié)點(diǎn)2和節(jié)點(diǎn)4的氣壓在經(jīng)過(guò)一段延時(shí)后也會(huì)發(fā)生相應(yīng)的波動(dòng)，并且天然氣網(wǎng)絡(luò)各節(jié)點(diǎn)的氣壓波動(dòng)需要一個(gè)較長(zhǎng)的過(guò)渡過(guò)程才能逐漸恢復(fù)穩(wěn)定。

圖6 算例2仿真波形

3 結(jié)語(yǔ)

針對(duì)電-氣互聯(lián)綜合能源系統(tǒng)中電力網(wǎng)絡(luò)與天然氣網(wǎng)絡(luò)在動(dòng)態(tài)運(yùn)行過(guò)程中的相互影響進(jìn)行了初步研究。選取一個(gè)微電網(wǎng)和天然氣網(wǎng)絡(luò)相互耦合的終端能源系統(tǒng)，研究了系統(tǒng)中各組成單元的數(shù)學(xué)模型，并在MATLAB/Simulink仿真平臺(tái)上完成了系統(tǒng)仿真模型的搭建。在此基礎(chǔ)上，通過(guò)仿真研究了微電網(wǎng)的動(dòng)態(tài)過(guò)程對(duì)天然氣網(wǎng)絡(luò)的影響，分析了微電網(wǎng)和天然氣網(wǎng)絡(luò)之間的耦合環(huán)節(jié)。分析了當(dāng)微型燃?xì)廨啓C(jī)發(fā)電系統(tǒng)和P2G系統(tǒng)的輸出發(fā)生波動(dòng)時(shí)天然氣管道中氣體壓力的動(dòng)態(tài)變化過(guò)程。研究成果為進(jìn)一步分析綜合能源系統(tǒng)的動(dòng)態(tài)過(guò)程提供了理論基礎(chǔ)，為深入研究綜合能源系統(tǒng)的運(yùn)行控制和故障特征提供了試驗(yàn)平臺(tái)。