圖24所示，具體步驟如下：1）輸入各能源參數(shù)及耦合單元參數(shù)，并根據(jù)原理方程計算能量流關(guān)系式；2）基于N-R法求解電力子系統(tǒng)的潮流計算；3）更新電力系統(tǒng)的負荷參數(shù)，求得電力系統(tǒng)的功率需求；4）基于N-R法求解天然氣系統(tǒng)的能流計算，求得天然氣系統(tǒng)中電壓縮機等用電設(shè)備的功率需求；5）判斷電壓縮機的功率前后兩次迭代結(jié)果的差值是否滿足精度要求，若滿足則輸出計算結(jié)果，若不滿足則重復(fù)步驟2）至5），直至輸出結(jié)果。圖STYLEREF1\s2SEQ圖表\*ARABIC\s14電-氣區(qū)域綜合能源系統(tǒng)潮流計算流程圖1.5算例分析1.5.1算例系統(tǒng)在配置為i5處理器、8GBRAM的計算機上，采用軟件MatlabR2018a，編制電-氣區(qū)域綜合能源系統(tǒng)潮流求解程序。電-氣區(qū)域綜合能源算例系統(tǒng)的拓撲結(jié)構(gòu)如REF_Ref68013447\h圖25。其中，配電網(wǎng)的基準容量為100MVA，基準頻率為50Hz。配電網(wǎng)總負荷有功為3.7150p.u.，負荷無功為1.3000p.u.；2~32節(jié)點均為PQ節(jié)點，給定DG1、DG2、DG3的注入功率為1.25+j0.5p.u.；IEEE33節(jié)點配電網(wǎng)負荷參數(shù)如表2-1所示，線路參數(shù)如表2-2所示。天然氣網(wǎng)絡(luò)中，天然氣節(jié)點G1為氣源節(jié)點壓力為60bar，即平衡節(jié)點；天然氣節(jié)點G8為恒出口壓力壓縮機壓力為60bar，壓縮比大于1.2，小于1.8，效率為0.9；設(shè)置標(biāo)準狀態(tài)下天然氣流動指數(shù)m為2、多變指數(shù)α為1.27。天然氣網(wǎng)絡(luò)管徑參數(shù)如表2-3所示，EH位置參數(shù)如表2-4所示。圖STYLEREF1\s2SEQ圖表\*ARABIC\s15電-氣區(qū)域綜合能源系統(tǒng)拓撲結(jié)構(gòu)表2-1負荷參數(shù)表節(jié)點號負荷/p.u.節(jié)點號負荷/p.u.20.1+j0.06180.09+j0.0430.09+j0.04190.09+j0.0440.12+j0.08200.09+j0.0450.06+j0.03210.09+j0.0460.06+j0.02220.09+j0.0470.2+j0.1230.09+j0.0580.2+j0.1240.42+j0.290.06+j0.02250.42+j0.2100.06+j0.02260.06+j0.025110.045+j0.03270.06+j0.025120.06+j0.035280.06+j0.02130.06+j0.035290.12+j0.07140.12+j0.08300.2+j0.6150.06+j0.01310.15+j0.07160.06+j0.02320.21+j0.1170.06+j0.02330.06+j0.04表2-2線路參數(shù)表編號首-末節(jié)點阻抗/p.u.線路首-末節(jié)點阻抗/p.u.11-20.000574+j0.00021717-180.004558+j0.010722-30.00307+j0.00156182-190.001021+j0.000933-40.002279+j0.00111919-200.009366+j0.008444-50.002373+j0.00122020-210.00255+j0.0029755-60.0051+j0.0044022121-220.004414+j0.005866-70.001166+j0.0038223-230.002809+j0.001977-80.00443+j0.001462323-240.005592+j0.004488-90.006413+j0.00462424-250.005579+j0.004399-100.006501+j0.0046256-260.001264+j0.00061010-110.001224+j0.00042626-270.00177+j0.000901111-120.002332+j0.00072727-280.006594+j0.00581212-130.009141+j0.00712828-290.005007+j0.00431313-140.003372+j0.00442929-300.00316+j0.001611414-150.00368+j0.003273030-310.006067+j0.00591515-160.004647+j0.00333131-320.001933+j0.00221616-170.008026+j0.01073232-330.002123+j0.0033表2-3天然氣管徑參數(shù)管道編號管徑/mm管道長度/m管道編號管徑/mm管道長度/m1（1-2）40020008（8-9）30010002（2-3）20010009（8-10）22518003（2-4）200100010（8-11）22518004（2-5）225180011（8-12）22518005（2-6）17550012（12-13）20010006（2-7）200100013（12-14）20010007（2-8）3002000表2-4EH位置參數(shù)EH編號電力節(jié)點天然氣節(jié)點編號EH編號電力節(jié)點天然氣節(jié)點編號EH13G3EH620G9EH28G4EH723G10EH312G5EH829G11EH413G6EH932G13EH515G71.5.2潮流計算結(jié)果及分析采用電-氣區(qū)域綜合能源系統(tǒng)的分解求解法求解系統(tǒng)潮流可得，EH的參數(shù)計算結(jié)果如表2-5所示，電力網(wǎng)絡(luò)潮流計算結(jié)果如REF_Ref68291478\h圖26~REF_Ref68817861\h圖28所示，天然氣網(wǎng)絡(luò)潮流計算結(jié)果如REF_Ref68291525\h圖29~REF_Ref68291531\h圖210所示。從計算結(jié)果可以發(fā)現(xiàn)，電力網(wǎng)絡(luò)中電壓幅值均處于1.005(p.u.)~1.035(p.u.)之間，電壓未發(fā)生越限；節(jié)點25、節(jié)點17以及節(jié)點33處和相連，因此，節(jié)點的電壓幅值高于周圍其余節(jié)點；符合系統(tǒng)要求。天然氣網(wǎng)絡(luò)中，天然氣節(jié)點的壓力維持在30bar以上，能夠?qū)崿F(xiàn)系統(tǒng)的供需平衡。

表2-5EH參數(shù)計算結(jié)果EH編號電能輸入量/p.u.天然氣輸入量/m3h-1電能輸出量/p.u.10.4530277.770.192320.6325333.3330.230730.2991333.3330.230740.8291221.2220.153850.3675388.8890.269260.7179221.2220.153870.7607166.6670.115380.5641277.7780.192390.4957221.2220.1538圖STYLEREF1\s2SEQ圖表\*ARABIC\s16電力網(wǎng)絡(luò)節(jié)點電壓圖STYLEREF1\s2SEQ圖表\*ARABIC\s17電力網(wǎng)絡(luò)節(jié)點注入功率圖STYLE