摘要：我國當前電力系統(tǒng)已經(jīng)開始向大規(guī)模和高集成化方向發(fā)展，Matlab軟件在研究過程中可以對電力系統(tǒng)發(fā)揮模擬作用，但是對于電力潮流來說，需要在獲得數(shù)據(jù)以及具體節(jié)點之后，通過程序處理獲得參數(shù)。在電力系統(tǒng)潮流計算階段需要選用最科學的算法，由于Gauss-Seidel迭代法可以降低計算誤差，所以在提高計算精度方面發(fā)揮著重要作用。基于此，介紹了Gauss-Seidel迭代算法基本原理以及Gauss-Seidel迭代算法用于潮流計算的基本原理，進一步分析了基于Matlab的Gauss-Seidel迭代算法電力系統(tǒng)潮流計算方法。

關(guān)鍵詞：Matlab軟件;Gauss-Seidel迭代法;電力系統(tǒng)潮流計算

0 引言

在迭代法的使用過程中，無論何種處理方法都會形成一定誤差，經(jīng)過多次迭代后，當采用同一數(shù)值參與計算時，則該誤差會持續(xù)增加，甚至在一定情況下產(chǎn)生的誤差會高于原始數(shù)值。對于本文選用的Gauss-Seidel迭代法，其優(yōu)勢是能夠在形成迭代結(jié)果后，通過對該結(jié)果與下一步處理結(jié)果進行橫向?qū)Ρ炔⒅苯永?，從而形成精準的計算結(jié)果。

1 Gauss-Seidel迭代算法基本原理

1.1 ? ?算法使用流程

該算法的使用流程是根據(jù)數(shù)學原理得到某一參數(shù)，計算方法可以表現(xiàn)為y=f（x），則迭代過程可以表示為：

y（1）=f（y）y（2）=f（y（1））y（3）=f（y（2））y（n）=f（y（n-1））

將不同流程得到的外值代入到同一個計算方程，以獲得第二步結(jié)果，經(jīng)過多次運算之后，每一步得到的結(jié)果都會產(chǎn)生一定差異。為了確定最終取得的運算參數(shù)，需要設(shè)定對比參數(shù)，計算方式為：

|y（n）-y（n-1）|≤ε

其中方程右側(cè)為計算對象在多次計算中允許的最大誤差值，需要注意的是，對于選用的算法中，函數(shù)加工廠“f（x）”中可以有多個自變量，但是因變量只能有一個，即：

y=f（x1，x2，x3，…，xn）

要根據(jù)實際的計算對象，合理確定自變量的數(shù)量和數(shù)值。

1.2 ? ?相關(guān)參數(shù)的處理

在參數(shù)處理過程中，一個參數(shù)為自變量，另一個參數(shù)為因變量，無論是何種算法都會設(shè)定一個標準值，而該標準值通常用于說明計算過程中產(chǎn)生的各階段的誤差程度[1]。通常情況下，該結(jié)果需要按照行業(yè)以及學科知識的要求設(shè)置，而具體計算工作由專業(yè)算法和軟件自主完成，并且只有在確定多次計算結(jié)果之后，才可以得到最精確的指標。

1.3 ? ?最終結(jié)果的獲取

事實上，在每一次迭代過程中都會產(chǎn)生一定的固定結(jié)果，但是不同的計算階段數(shù)據(jù)生成的實際誤差存在差異，只有在誤差能夠保持學科要求時，獲得的結(jié)果才可實際使用。另外在每一步的計算階段，都需要將產(chǎn)生的誤差值加入到持續(xù)性的迭代算法內(nèi)，而該迭代算法在融合了誤差值后具有更高的精度，利用該方法可以使實際獲得的結(jié)果具有可靠性。

2 Gauss-Seidel迭代算法用于潮流計算的基本原理

2.1 ? ?計算參數(shù)的確定

計算參數(shù)的確定要根據(jù)電力系統(tǒng)的自身運行標準和研究對象確定，比如對于電力潮流來說，由于主要的評價標準是某節(jié)點中的功率和其余節(jié)點間的功率比值，尤其是對于線路的上游和下游區(qū)域，當下游的實際功率高于上游時，則會出現(xiàn)逆潮流現(xiàn)象，對電力系統(tǒng)造成破壞，所以無論是實際檢測，還是模擬運算，首先確定的參數(shù)是當前電力系統(tǒng)中的各個節(jié)點，其次通過對電壓電流的分析，確定這兩個節(jié)點的功率數(shù)值，并將獲得的參數(shù)納入計算方程中，可以在多次迭代之后和設(shè)定的標準值對比，當發(fā)現(xiàn)某參數(shù)經(jīng)過多次迭代之后和標準值之間的誤差小于設(shè)定標準，可確定當前整個潮流系統(tǒng)處于穩(wěn)定運行狀態(tài)。最后是實際分析，制定線路系統(tǒng)的優(yōu)化方案。

2.2 ? ?計算公式的確定

在計算公式的建設(shè)階段，要確定潮流的變化方法和計算流程，并且將各時間點和工作步驟中產(chǎn)生的實際參數(shù)納入后續(xù)的計算范疇。通常情況下，由于迭代算法中計算公式不出現(xiàn)變化，所以可按照該方式研究最終的計算成果，讓獲得的結(jié)果可以更好地描述實際情況[2]。考慮到功率的計算方式為電壓和電流的乘積，且對電流的計算和檢測難度較小，所以主要分析當前不同節(jié)點的電流數(shù)值，此外根據(jù)不同節(jié)點的有功和無功功率數(shù)值，可獲得第2個節(jié)點的電壓以及第n個節(jié)點的電壓數(shù)值，將其納入實際的潮流計算公式內(nèi)并可以得到實際的計算結(jié)果，以對平衡節(jié)點的電力參數(shù)進行討論。由于已經(jīng)在計算中獲得了電壓值，將該參數(shù)直接代入潮流迭代公式中進行多次迭代計算即可。由于向某節(jié)點注入的有功功率和電壓大小已經(jīng)給定，所以可以將其使用到迭代計算公式內(nèi)，并對該節(jié)點的無功功率預先計算，根據(jù)功率和電流的表達式，并且加入需要考慮的電壓值參數(shù)，則可使獲得的結(jié)果得以精確計算。

2.3 ? ?計算節(jié)點的確定

在計算節(jié)點的確定過程中，由于不同的計算方法具有一定的誤差，比如某節(jié)點為發(fā)電機節(jié)點而另一節(jié)點為輸電端節(jié)點，這兩個節(jié)點的計算公式必然不同。基于Matlab的數(shù)據(jù)分析過程要求，要根據(jù)已經(jīng)建成的電力系統(tǒng)模型確定不同的節(jié)點方式，同時將獲得的電壓參數(shù)和電流參數(shù)代入迭代算法，獲得的不同電壓值和電流值都可以納入該算法體系，以得到最終結(jié)果，并通過比較確定潮流情況。

3 基于Matlab的Gauss-Seidel迭代算法電力系統(tǒng)潮流計算方法

3.1 ? ?編程原理的確定

編程原理的確定階段要根據(jù)Matlab的專用編程語言，在實際研究過程確定某節(jié)點類型，分為發(fā)電機節(jié)點、配電側(cè)節(jié)點等，由于計算參數(shù)和計算方式已知，所以在編程階段可以將工作數(shù)值納入實際的編程體系中，并且按照時間變化設(shè)定相關(guān)參數(shù)。另外對于共軛電壓值，只需要做出相關(guān)處理即可，采用迭代算法可以得到節(jié)點電壓和運行相位角，在計算過程需要建立不同節(jié)點的導納矩陣，對其初始化后設(shè)定初始運行狀態(tài)值，并分析在經(jīng)過一段時間以及多步驟分析后可以取得的作用效果，把獲得的電壓值輸入原有的迭代算法，同時分析在不同時間節(jié)點中該節(jié)點的實際運行參數(shù)。

3.2 ? ?軟件編程階段

軟件編程階段要容納電力學的各類參數(shù)、電力系統(tǒng)各節(jié)點的運行位置、不同節(jié)點電力學參數(shù)的變化方程等，通過所有數(shù)據(jù)的分析才可以將其加入實際的潮流計算階段。此外在后續(xù)的運行過程需要對相關(guān)方程和編程結(jié)果進行預處理，該過程中一方面會設(shè)定初始運行值，另一方面要研究該數(shù)值的后續(xù)運行情況，以形成實際的運行指標，當確定該指標符合對比值要求時，則可確定目前的節(jié)點潮流可以符合系統(tǒng)的正常運行指標。

3.3 ? ?數(shù)據(jù)導入工作

在數(shù)據(jù)導入階段，首先要設(shè)置整個供配電網(wǎng)絡(luò)的初始運行值，包括電壓參數(shù)、電流參數(shù)以及節(jié)點參數(shù)，并設(shè)定不同節(jié)點的電壓變化情況，按照實際運行階段的電力參數(shù)運行模式設(shè)定[2]。其次是設(shè)定專用的描述方程，包括電壓的變化方程、電流的變化方程以及實際的對比方程等，并且將實際的迭代算法與電力學的相關(guān)變化方程進行對比，當發(fā)現(xiàn)實際的作用階段，數(shù)據(jù)在經(jīng)過了長時間運行之后，可以取得在不同運行時間點上的節(jié)點功率參數(shù)，經(jīng)過橫向?qū)Ρ?，可研究當前該參?shù)的差值是否超出了限定標準，若發(fā)現(xiàn)存在這一問題時，可確定當前的作用形勢下該電力系統(tǒng)無法保持穩(wěn)定運行狀態(tài)。最后是對電力系統(tǒng)的優(yōu)化，以解決目前存在的潮流問題。

3.4 ? ?計算結(jié)果的獲取

在計算結(jié)果的取得和研究階段，要研究的參數(shù)包括電力輸出節(jié)點的電壓初始值、發(fā)電機節(jié)點的相位角初始值以及經(jīng)過一段運行時間之后的電力參數(shù)，并讓該參數(shù)參與對比，尤其是對于發(fā)電機節(jié)點和功率輸出節(jié)點，要分析不同發(fā)電機節(jié)點的運行參數(shù)。通過該形式一方面可以研究當前發(fā)電機的發(fā)電功率是否可以符合用電指標，另一方面可研究不同的發(fā)電機節(jié)點運行過程是否會出現(xiàn)逆潮流問題。比如，發(fā)現(xiàn)上游的發(fā)電機節(jié)點功率參數(shù)為U1，下游的發(fā)電機節(jié)點功率參數(shù)是U2，其中前個數(shù)值小于后個數(shù)值，則可確定該供電系統(tǒng)會出現(xiàn)逆潮流問題。

4 結(jié)語

綜上所述，基于Matlab的迭代算法使用階段，由于迭代算法可以對不同運行過程的操作進行全方位考慮，所以在實際運行過程中，可通過當前電力系統(tǒng)的相關(guān)參數(shù)確定運行狀態(tài)，找到實際的運行指標。Matlab軟件經(jīng)過編程之后，可以在數(shù)據(jù)優(yōu)化和基礎(chǔ)數(shù)據(jù)獲取階段，經(jīng)過多次計算得到成果，以便研究該電力系統(tǒng)是否存在故障以及逆潮流問題。

[參考文獻]

[1] 韋雅，文耀寬，田珂，等.基于時間序列線性大數(shù)據(jù)分析的電力系統(tǒng)潮流計算方法[J].機械與電子，2020，38（1）：32-35.

[2] 安軍，宋俊達，葛維春.大規(guī)模電力系統(tǒng)潮流計算收斂性診斷與調(diào)整方法[J].電力自動化設(shè)備，2020（2）：1-6.

收稿日期：2020-01-16

作者簡介：宋陽（1979—），男，四川樂山人，碩士研究生，高級工程師，研究方向：電力系統(tǒng)分析、電網(wǎng)穩(wěn)定分析與潮流計算等。