CN105096209A - 一种针对电网断面热稳定极限的快速评估方法 - Google Patents

Abstract

本发明提供了一种针对电网断面热稳定极限的快速评估方法，包括以下步骤：(1)计算电网的初始潮流；(2)对断面进行N-1开断故障仿真；(3)判定断面是否运行与热稳定极限，根据判断结果转至步骤(4)或(9)；(4)生成断面支路有功潮流变化率矩阵λ；(5)生成断面运行于热稳极限状态下的断面有功潮流矩阵P_T；(6)根据得到有功潮流矩阵P_T对断面的热稳定极限进行评估；(7)根据评估结果调整电网运行方式；(8)对调整后的运行方式进行仿真，根据判断结果转至步骤(4)或(9)；(9)输出断面热稳定极限的评估结果。通过该评估方法，可实现电网断面热稳定极限的快速评估，并准确定位限制故障的类型，为电网调度人员提供一套可靠的热稳定极限评估方法、流程，为软件开发人员提供相关的理论依据。

Description

一种针对电网断面热稳定极限的快速评估方法

技术领域

本发明属于电力***领域，具体涉及一种针对电网断面的热稳定极限进行快速评估的方法。

背景技术

随着我国电网的不断发展，主网架结构日益增强。为提升运行的可靠性和灵活性，以特高压交直流线路构成主网架的特点更加突出。在电网结构发展变化的过程中，电网稳定问题的主要矛盾和焦点由发展初期的暂稳问题，逐步过渡到受限于电力设备自身容量的热稳定问题。因此，研究电网的热稳定问题及其解决方案，尤其是断面的热稳极限问题已经成为电网生产实际中不可回避的重要课题。

关于断面热稳定极限的评估方法，国内外开展了不少研究，包括：线性或二次规划法、连续潮流法、遗传算法、蒙特卡罗模拟法等。采用这类方法往往需要完善的后期软件维护工作，由于电网发展变化较快，而相关软件的实时更新无法尽如人意。因此在实际电网调度中，涉及到断面热稳定极限的求取仍需要调度人员进行手动计算，而这类方法与调度人员实际采用的方式调整方法差别往往较大，不利于在实际电网中的推广应用。

发明内容

为了克服现有技术的上述缺陷，本发明提出了一种能够快速基于断面潮流按阻抗分配原理的热稳定极限快速评估方法。

为实现上述发明目的，本发明是通过如下技术方案实现的：

一种电网断面热稳定极限的快速评估方法，包括以下步骤：

(1).计算构成电网断面各条线路的初始有功潮流；

(2).对电网断面进行N-1开断故障仿真，以得到电网断面各条线路开断后在运线路的有功潮流；

(3).判定电网断面各条线路是否运行于热稳定极限，根据判断结果选择执行步骤(4)或步骤(9)；

(4).生成电网断面支路有功潮流变化率矩阵λ；

(5).生成运行于热稳极限状态下的电网断面有功潮流矩阵P_T；

(6).根据电网断面有功潮流矩阵P_T对电网断面的热稳定极限有功功率进行评估；

(7).根据评估结果调整电网运行方式；

(8).对调整后的电网运行方式进行N-1开断故障仿真，根据判断结果选择执行步骤(4)或步骤(9)；

(9).输出断面热稳定极限的评估结果。

优选地，所述步骤(1)中，通过PSD-BPA潮流程序计算构成电网断面各条线路的初始有功潮流。

优选地，所述步骤(2)中，通过PSD-BPA暂态稳定程序对所研究的电网断面进行N-1开断故障仿真，得到电网断面各条线路开断后运行线路的有功潮流；进行所述开断故障仿真的元件为：组成电网断面的N条线路，N为正整数。

优选地，所述步骤(3)的方法为：

首先，通过下式将电网断面各条线路的热稳定极限电流折算为热稳定极限有功功率：

式中，I_max为电网断面中线路的热稳定极限电流，I_max＝KI_s；K为温度校正系数；I_s为环境温度为25℃时线路的安全电流；为功率因数，一般取值为0.95；P_t为电网断面中线路的热稳定极限有功功率；

其次，判断电网断面各条线路的有功潮流是否达到了相应线路的热稳定极限有功功率：

当断面中n条线路的有功潮流先于其它(N-n)条线路达到了其热稳定极限有功功率P_t值，说明所研究断面的有功潮流达到了热稳定极限有功功率，则执行步骤(9)；反之，执行步骤(4)；其中，n为整数且1≤n≤N。

优选地，所述步骤(4)包括：电网断面线路N-1开断后，根据运行线路有功潮流的变化生成有功潮流变化率矩阵λ如下所示：

其中，将矩阵λ中第k行、第j列的参数定义为线路j开断在线路k上产生的有功潮流变化率λ_kj，η_kj为线路j开断后在线路k上的有功潮流转移比，k≠j，x_k、x_j分别为线路k、j的电抗。

优选地，所述步骤(5)中，运行于热稳定极限状态下的电网断面有功潮流矩阵P_T如下所示：

其中，矩阵P_T各元素的含义为：对应于矩阵P_T的第j列，即为线路j开断前电网断面中各条线路的有功潮流；矩阵的对角线元素P_j通过下式求得：

$P_j=\underset{\underset{k\≠j}{k=1}}{\overset{n}{\mathrm{\Σ}}}\mathrm{\Δ}{P}_{\mathrm{kj}}=\underset{\underset{k\≠j}{k=1}}{\overset{n}{\mathrm{\Σ}}}(P_{\mathrm{kt}}-P_{\mathrm{kj}})$

式中，P_kt为电网断面中线路k的热稳定极限有功功率，I_kmax为电网断面中线路k的热稳定极限电流，I_kmax＝KI_ks，K为温度校正系数，I_ks为环境温度为25℃时线路k的安全电流，为功率因数；P_kj为线路k在线路j未开断前的有功潮流，k≠j。矩阵P_T的反映并给出了断面各元件开断前的潮流分布情况。

优选地，所述步骤(6)中，对矩阵P_T各列分别求和，得到电网断面在线路j开断前的运行功率，所有求和值中的最小值，即为电网断面的实际热稳定极限有功功率P_Σt，所述P_Σt如下所示：

$P_{\mathrm{\Σt}}=\min (P_j+\underset{\underset{k\≠j}{k=1}}{\overset{N}{\mathrm{\Σ}}}{P}_{\mathrm{kj}})j=1,\·\·\·,N$

式中，对应于P_Σt取最小值时的j，即为电网断面热稳定极限有功功率的支路开断约束故障；N为组成电网断面的线路总数，P_j为线路j开断前电网断面各线路的有功潮流，P_kj为线路k在线路j未开断前的有功潮流。因此，采用本方法即可以快速地对断面热稳极限进行评估，同时也定位了断面的约束故障。

优选地，所述步骤(7)中，根据电网断面热稳定极限有功功率的评估结果，得到电网断面当前有功潮流与热稳定极限有功功率间的差额，根据此差额调整电网断面两端的开机量。

优选地，所述步骤(8)中，通过PSD-BPA暂态稳定程序对调整后的电网断面进行N-1开断故障仿真，仿真所选择的故障为步骤(6)中确定的断面约束故障。若电网断面达到热稳定极限有功功率，则执行步骤(9)；否则，执行步骤(4)。

优选地，所述步骤(9)中，输出电网断面热稳定极限评估结果包括：电网断面的热稳定极限有功功率和对应的约束故障。

与现有技术相比，本发明的有益效果在于：

基于断面潮流按阻抗分配原理的热稳定极限快速评估方法，可以实现电网断面热稳定极限的快速评估，并能准确定位断面热稳定极限的约束故障。本方法最大的优点在于可实现热稳定极限的快速手动求解，为调度人员进行热稳定极限计算提供了一种标准方法、流程，为软件开发人员提供了可供参考的理论依据。

附图说明

图1经由断面互联的两等值***示意图；

图2断面热稳定极限的评估流程示意图；

图3豫北断面一次接线示意图。

具体实施方式

下面结合附图对本发明的评估方法做进一步的详细说明。

本例所提出的一种针对电网断面热稳定极限的快速评估方法，用于调度生产中对电网断面热稳定极限的快速评估，包括以下步骤：

(1).计算电网的初始潮流；

(2).对断面进行N-1开断故障仿真；

(3).判定断面是否运行与热稳定极限，根据判断结果转至步骤(4)或(9)；

(4).生成断面支路有功潮流变化率矩阵λ；

(5).生成断面运行于热稳极限状态下的断面有功潮流矩阵P_T；

(6).根据得到有功潮流矩阵对断面的热稳定极限进行评估；

(7).根据评估结果调整电网运行方式；

(8).对调整后的运行方式进行N-1开断仿真，根据判断结果转至步骤(4)或(9)；

(9).输出断面热稳定极限的评估结果。

在步骤(1)中，对电网进行一次初始潮流计算，确定所研究断面的初始运行状态，得到构成电网断面的各条线路的有功潮流。

在步骤(2)中，对研究的电网断面进行N-1开断故障仿真，对组成断面的N条线路进行一轮开断。

在步骤(3)中，断面各线路的热稳定极限折算为有功潮流，对应的公式如下：

I_max＝KI_s

上式中，I_max为线路的热稳定极限电流，I_s为环境温度为25℃时线路的安全电流，K为温度校正系数，功率因数一般取0.95，P_t为线路的热稳定极限有功功率。当断面中n条线路的有功潮流先于其它(N-n)条线路达到了其热稳定极限有功功率P_t值，则所研究断面的有功潮流达到了热稳定极限，其中n的取值范围是：1≤n≤N。

在步骤(4)中，断面线路N-1开断后，根据在运线路有功潮流的变化生成有功潮流变化率矩阵λ。其具体定义可结合附图1来予以解释。附图1所示的两等值***间断面有功潮流可按式(1)进行计算。

$P_i=\frac{V_s{V}_r\sin {\mathrm{\δ}}_{\mathrm{sr}}}{x_i}$

$P_{\mathrm{\Σ}}=V_s{V}_r\sin {\mathrm{\δ}}_{\mathrm{sr}}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\frac{1}{x_i}---\left(1\right)$

i＝1…N

由式(1)可以看出，断面的有功潮流按照各条线路的电抗成比例分配，即P_i1/P_i2＝x_i2/x_i1。

若线路l_j开断，由于断面两侧***有功的供需没有发生变化，断面总的有功潮流将保持不变，线路l_j开断前有功潮流P_j将分配至断面其它线路。其中，线路l_k在l_j开断后的有功潮流如式(2)所示。

${\stackrel{\‾}{P}}_{\mathrm{kj}}=\frac{V_s{V}_r}{x_k}\sin {\stackrel{\‾}{\mathrm{\δ}}}_{\mathrm{srj}}---\left(2\right)$

因此，可以得到线路l_k有功潮流的增量ΔP_kj如式(3)所示。

$\mathrm{\Δ}{P}_{\mathrm{kj}}={\stackrel{\‾}{P}}_{\mathrm{kj}}-P_k=\frac{V_s{V}_r}{x_k}(\sin {\stackrel{\‾}{\mathrm{\δ}}}_{\mathrm{srj}}-\sin {\mathrm{\δ}}_{\mathrm{sr}})---\left(3\right)$

P_j即为在运各条线路有功潮流增量之和，如式(4)所示。

$P_j=\underset{k=1}{\overset{N}{\mathrm{\Σ}}}\mathrm{\Δ}{P}_{\mathrm{kj}}=V_s{V}_r(\sin {\stackrel{\‾}{\mathrm{\δ}}}_{\mathrm{srj}}-\sin {\mathrm{\δ}}_{\mathrm{sr}})\underset{k=1}{\overset{N}{\mathrm{\Σ}}}\frac{1}{x_k}---\left(4\right)$

其中，k≠j。根据式(3)～(4)，可以得到线路l_j开断后在线路l_k上的有功潮流转移比η_kj，如式(5)所示。

${\mathrm{\η}}_{\mathrm{kj}}=\frac{\mathrm{\Δ}{P}_{\mathrm{kj}}}{P_j}=\frac{\underset{k=1}{\overset{N}{\mathrm{\Σ}}}\frac{1}{x_k}}{x_k}---\left(5\right)$

k≠j

可以看出，P_j将按照电抗比例关系分配至剩余的各条线路。至此，可以给出线路l_j开断在线路l_k上产生的有功潮流变化率λ_kj定义，其具体含义如式(6)所示。

${\mathrm{\λ}}_{\mathrm{kj}}=\frac{\mathrm{\Δ}{P}_{\mathrm{kj}}}{P_k}=\frac{{\mathrm{\η}}_{\mathrm{kj}}{P}_j}{P_k}={\mathrm{\η}}_{\mathrm{kj}}\frac{x_k}{x_j}---\left(6\right)$

λ_kj反映了线路l_j的开断在线路l_k上造成的有功潮流变化，由式(6)可以看出，λ_kj是与断面线路参数相关的确定值。根据式(6)，若k＝j，则有ΔP_kk＝P_k，因此λ_kk＝1.0。

基于λ_kj的含义，基于初始潮流，对附图1所示断面进行一轮开断模拟之后，可以得到一系列λ_kj的值，生成矩阵λ，如下所示：

在步骤(5)中，附图1所示***断面若已运行在热稳定极限状态下，则必将会有某条线路l_j的开断，导致至少一条线路l_k达到其热稳极限的情况出现。因此，可以计算得到线路l_k在线路l_j未开断前的有功潮流如下所示：

${\mathrm{\λ}}_{\mathrm{kj}}=\frac{P_{\mathrm{kt}}-P_{\mathrm{kj}}}{P_{\mathrm{kj}}}$ (7)

$P_{\mathrm{kj}}=\frac{P_{\mathrm{kt}}}{1+{\mathrm{\λ}}_{\mathrm{kj}}},k\≠j$

式(7)中，P_kt为线路l_k的热稳限额。结合λ_kj的定义并根据根据式(7)，可以得到运行于热稳极限状态下的断面有功潮流矩阵P_T如式(8)所示：

其中，矩阵的对角线元素P_j根据式(4)求出，具体如式(9)所示：

$P_j=\underset{\underset{k\≠j}{k=1}}{\overset{n}{\mathrm{\Σ}}}\mathrm{\Δ}{P}_{\mathrm{kj}}=\underset{\underset{k\≠j}{k=1}}{\overset{n}{\mathrm{\Σ}}}(P_{\mathrm{kt}}-P_{\mathrm{kj}})---\left(9\right)$

至此，矩阵P_T中的所有元素均已求出。对应于矩阵P_T的第j列，即为线路j开断前断面各线路的有功潮流。矩阵P_T的反映并给出了断面各元件开断前的潮流分布情况。

在步骤(6)中，对矩阵P_T各列分别求和，可以得到断面在线路j开断前的运行功率，这些和中的最小值，即为断面的实际热稳定极限P_Σt，如下所示：

$P_{\mathrm{\Σt}}=\min (P_j+\underset{\underset{k\≠j}{k=1}}{\overset{N}{\mathrm{\Σ}}}{P}_{\mathrm{kj}})j=1,\·\·\·,N$

对应于P_Σt取最小值时的j，即为断面热稳定极限的支路开断约束故障。因此，采用本方法即可以快速地对断面热稳极限进行评估，同时也定位了断面的约束故障。

步骤(7)中，根据断面热稳定极限的评估结果，可以得到断面当前有功潮流与热稳极限值之间的差额，根据此差额调整断面两端的开机量。

步骤(8)中，对调整后的电网断面进行N-1开断仿真，所选择的故障为步骤(6)中确定断面约束故障。若断面达到热稳定极限，则转至步骤(9)，否则转至步骤(4)。

步骤(9)中，输出断面热稳定极限评估结果，包括：断面的热稳定极限值和对应的约束故障。

【实施例】

下面以河南电网2014年夏季大负荷方式为例，采用本文提出的方法对豫北断面的热稳定极限进行评估。豫北断面的接线构成如附图3所示，其由豫HH双回线和豫TC单回线组成。豫HH双回线每条线路的热稳定极限功率为2056MW，豫TC单回线的热稳定极限功率为1650MW。

(1)计算电网的初始潮流

采用中国版BPA潮流计算程序(详见中国电力科学研究院于2005年9月公开的PSD-BPA潮流程序用户手册)对河南电网2014年夏季大负荷方式进行计算，得到豫北断面的初始潮流，断面线路的初始有功潮流计算结果如下表1所示：

表1豫北断面各线路的初始潮流与热稳限额

	豫HH I回线	豫HH II回线	豫TC线
				初始潮流	595.8	595.8	644.9
热稳限额	2056	2056	1655

(2)对断面进行N-1开断故障仿真

采用中国版BPA暂态稳定计算程序(详见中国电力科学研究院于2006年6月公开的PSD-BPA暂态稳定程序用户手册)对豫北断面进行N-1开断的故障仿真，得到断面各条线路一轮开断后在运线路的有功潮流值。

(3)判定断面是否运行与热稳定极限，根据判断结果转至步骤(4)或(9)

根据附表2，豫北断面线路第一轮N-1开断后，在运线路豫HH单回有功潮流为990MW，豫TC单回有功潮流为850MW，均未达到热稳定极限。因此转入步骤(4)。

表2豫北断面的热稳定极限的评估结果

(4)生成断面支路有功潮流变化率矩阵λ

根据豫北断面线路N-1开断后在运线路有功潮流的变化情况，结合有功潮流变化率矩阵λ定义，可生成豫北断面的λ矩阵。如下所示：

$\mathrm{\λ}=\left[\begin{array}{ccc}1& 0.662& 0.541\\ 0.662& 1& 0.541\\ 0.318& 0.318& 1\end{array}\right]$

(5)生成断面运行于热稳极限状态下的断面有功潮流矩阵P_T

基于豫北断面线路第一轮开断得到的有功潮流变化矩阵λ，以及各线路的热稳定极限值，可以生成豫北断面运行于热稳定极限状态下的断面有功潮流矩阵P_T，结果如下：

$P_T=\left[\begin{array}{ccc}1218& 1237.3& 1237.3\\ 1237.3& 1218& 1334.3\\ 1255.7& 1255.7& 1443.2\end{array}\right]$

(6)根据得到有功潮流矩阵P_T对断面的热稳定极限进行评估

对步骤(5)得到的有功潮流矩阵P_T各列求和，并取最小值，以定位断面热稳定极限的约束故障，第一轮断面线路N-1开断得到的P_T矩阵各列之和中，最小的为前两列，其值为3711MW。因此经过第一轮断面N-1开断后，豫北断面的热稳定极限评估为3711MW，约束故障为豫HH双回线N-1故障。

(7)根据评估结果调整电网运行方式

根据断面热稳定极限的评估结果，可以得到断面当前有功潮流与热稳极限值之间的差额，根据此差额调整断面两端的开机量。

(8)对调整后的运行方式进行N-1开断仿真，根据判断结果转至步骤(4)或(9)

对调整后的电网断面进行N-1开断仿真，所选择的故障为步骤(6)中确定断面约束故障。若断面达到热稳定极限，则转至步骤(9)，否则转至步骤(4)。

(9)输出断面热稳定极限的评估结果

对豫北断面进行3轮的N-1计算和方式调整后，得到了豫北断面的热稳定极限为3534MW，约束故障为豫HH双回线N-1故障。第一轮评估的结果误差为5％。因此，在对断面热稳定极限进行粗略评估时，进行一轮计算，也即计算进行到第(6)步，本方法也可以保证一定的评估结果精度。

该快速评估方法可以实现电网断面热稳定极限的快速评估，并能准确定位断面热稳定极限的约束故障，为调度人员进行热稳定极限计算提供了一种标准方法、流程，为软件开发人员提供了可供参考的理论依据。

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。

Claims (10)

1.一种针对电网断面热稳定极限的快速评估方法，其特征在于，所述方法包括以下步骤：

(1).计算构成电网断面各条线路的初始有功潮流；

(2).对电网断面进行N-1开断故障仿真，以得到电网断面各条线路开断后运行线路的有功潮流；

(3).判定电网断面各条线路是否运行于热稳定极限，根据判断结果选择执行步骤(4)或步骤(9)；

(4).生成电网断面支路有功潮流变化率矩阵λ；

(5).生成运行于热稳极限状态下的电网断面有功潮流矩阵P_T；

(6).根据电网断面有功潮流矩阵P_T对电网断面的热稳定极限有功功率进行评估；

(7).根据评估结果调整电网运行方式；

(8).对调整后的电网运行方式进行N-1开断故障仿真，根据判断结果选择执行步骤(4)或步骤(9)；

(9).输出断面热稳定极限的评估结果。

2.根据权利要求1所述的方法，其特征在于：所述步骤(1)中，通过PSD-BPA潮流程序计算构成电网断面各条线路的初始有功潮流。

3.根据权利要求1所述的方法，其特征在于，所述步骤(2)中，通过PSD-BPA暂态稳定程序对所研究的电网断面进行N-1开断故障仿真，得到电网断面各条线路开断后运行线路的有功潮流；进行所述开断故障仿真的元件为：组成电网断面的N条线路，N为正整数。

4.根据权利要求1所述的方法，其特征在于，所述步骤(3)的方法为：

首先，通过下式将电网断面各条线路的热稳定极限电流折算为热稳定极限有功功率：

式中，I_max为电网断面中线路的热稳定极限电流，I_max＝KI_s；K为温度校正系数；I_s为环境温度为25℃时线路的安全电流；为功率因数；P_t为电网断面中线路的热稳定极限有功功率；

其次，判断电网断面各条线路的有功潮流是否达到了相应线路的热稳定极限有功功率：

当断面中n条线路的有功潮流先于其它(N-n)条线路达到了其热稳定极限有功功率P_t值，说明所研究断面的有功潮流达到了热稳定极限有功功率，则执行步骤(9)；反之，执行步骤(4)；其中，n为整数且1≤n≤N。

5.根据权利要求1所述的方法，其特征在于，所述步骤(4)包括：电网断面线路N-1开断后，根据运行线路有功潮流的变化生成有功潮流变化率矩阵λ如下所示：

其中，将矩阵λ中第k行、第j列的参数定义为线路j开断在线路k上产生的有功潮流变化率λ_kj，η_kj为线路j开断后在线路k上的有功潮流转移比，k≠j，x_k、x_j分别为线路k、j的电抗。

6.根据权利要求1所述的方法，其特征在于：所述步骤(S)中，运行于热稳定极限状态下的电网断面有功潮流矩阵P_T如下所示：

其中，矩阵P_T各元素的含义为：对应于矩阵P_T的第j列，即为线路j开断前电网断面中各条线路的有功潮流；矩阵的对角线元素P_j通过下式求得：

$P_j=\underset{\underset{k\≠j}{k=1}}{\overset{n}{\mathrm{\Σ}}}\mathrm{\Δ}{P}_{\mathrm{kj}}=\underset{\underset{k\≠j}{k=1}}{\overset{n}{\mathrm{\Σ}}}(P_{\mathrm{kt}}-P_{\mathrm{kj}})$

式中，P_kt为电网断面中线路k的热稳定极限有功功率，I_kmax为电网断面中线路k的热稳定极限电流，I_kmax＝KI_ks，K为温度校正系数，I_ks为环境温度为25℃时线路k的安全电流，为功率因数；P_kj为线路k在线路j未开断前的有功潮流， $P_{\mathrm{kj}}=\frac{P_{\mathrm{kt}}}{1+{\mathrm{\λ}}_{\mathrm{kj}}},k\≠j.$

7.根据权利要求1所述的方法，其特征在于：所述步骤(6)中，对矩阵P_T各列分别求和，得到电网断面在线路j开断前的运行功率，所有求和值中的最小值，即为电网断面的实际热稳定极限有功功率P_Σt，所述P_Σt如下所示：

$P_{\mathrm{\Σt}}=\min (P_j+\underset{\underset{k\≠j}{k=1}}{\overset{N}{\mathrm{\Σ}}}{P}_{\mathrm{kj}})j=1,\·\·\·,N$

式中，对应于P_Σt取最小值时的j，即为电网断面热稳定极限有功功率的支路开断约束故障；N为组成电网断面的线路总数，P_j为线路j开断前电网断面各线路的有功潮流，P_kj为线路k在线路j未开断前的有功潮流。

8.根据权利要求1所述的方法，其特征在于：所述步骤(7)中，根据电网断面热稳定极限有功功率的评估结果，得到电网断面当前有功潮流与热稳定极限有功功率间的差额，根据此差额调整电网断面两端的开机量。

9.根据权利要求1所述的方法，其特征在于：所述步骤(8)中，通过PSD-BPA暂态稳定程序对调整后的电网断面进行N-1开断故障仿真，仿真所选择的故障为电网断面热稳定极限有功功率的支路开断约束故障；若电网断面达到热稳定极限有功功率，则执行步骤(9)；否则，执行步骤(4)。

10.根据权利要求1所述的方式，其特征在于：所述步骤(9)中，输出电网断面热稳定极限评估结果包括：电网断面的热稳定极限有功功率和对应的约束故障。