CN110096767B - 一种交直流混联电网连锁故障仿真方法 - Google Patents

Abstract

本发明公开一种交直流混联电网连锁故障仿真方法，包括：S1、确定电网的参数；S2、对电网进行循环仿真，相对于前一次仿真，按照一定比例提高发电机容量和负荷水平，按照一定比例增加前一次仿真中发生过故障的线路的容量；S3、每次仿真所有线路按照一定概率断开；S4、对断线后的电网进行潮流优化，判断是否有线路达到重载，若有线路达到重载，则进入S5，否则，进入S6；S5、对重载线路按照一定概率断开，如果有新线路断开则重复S4，否则，进入S6；S6、判断仿真次数是否达到预设次数，如果达到，结束仿真，否则，重复S2。本发明提供的连锁故障仿真方法，能够将直流线路模型和直流保护措施考虑进来，适用于交直流混联电网的连锁故障仿真。

Description

一种交直流混联电网连锁故障仿真方法

技术领域

本发明涉及电气工程领域，更具体地，涉及一种交直流混联电网连锁故障仿真方法。

背景技术

在当前我国大力建设交直流特高压互联电网的背景下，随着电网规模日益庞大以及特高压交直流输电工程的建设，***传输功率逐渐增长，***动态日益复杂，交直流混联电网的运行安全受到极大威胁。另一方面，***中的不确定性因素逐渐增多，由于故障、天气等因素导致的局部电网功率不平衡、失稳、振荡等问题可能扩散，从而导致连锁故障，甚至诱发大停电事故。研究表明，电力***大停电事故一般是由于***中某一元件故障导致一系列其他元件停运的连锁反应所造成，线路有功潮流的大规模转移导致输电线路连锁过载和保护的拒动/误动是连锁故障发生的主要因素。因此建立面向交直流混联电网的连锁故障模型，揭示***可能发生的连锁故障模式及其风险，为***连锁故障提供指导，将有效降低电网运行风险，产生巨大的经济和社会效益。而传统的连锁故障OPA仿真方法采用直流潮流优化模型，仅能应用于纯交流***，不能进行交直流混联电力***的连锁故障分析；因此，研究交直流混联电网的连锁故障仿真方法具有重要意义。

发明内容

针对现有技术的缺陷，本发明的目的在于解决传统的连锁故障OPA仿真方法仅能应用于纯交流***，不能进行交直流混联电力***的连锁故障分析的技术问题。

为实现上述目的，本发明提供一种交直流混联电网连锁故障仿真方法，包括以下步骤：

步骤S1、确定交直流混联电网的交流线路、直流线路、联络线、发电机、换流站以及负荷；发电机容量和负荷水平可用于计算电网的潮流，所述潮流包括各条线路的有功功率、无功功率、有功损耗、无功损耗以及电网各母线节点的电压幅值和相角；

步骤S2、对交直流混联电网进行预设次数的循环仿真，每次仿真相对于前一次仿真，按照第一预设比例提高发电机容量和负荷水平，按照第二预设比例增加前一次仿真中发生过故障的线路的容量；

步骤S3、每次仿真中，所有线路按照第一预设概率断开；

步骤S4、若交直流混联电网有线路断开，则对断线后的交直流混联电网进行潮流优化，将断开线路的功率分配到未断线的线路上，优化后判断是否有线路达到重载，若有线路达到重载，则进入步骤S5，否则，进入步骤S6；

步骤S5、对重载线路按照第二预设概率断开，然后判断有无新线路断开，如果有，则重复步骤S4，否则，进入步骤S6；

步骤S6、判断仿真次数是否达到预设次数，如果达到，则结束连锁故障仿真过程，否则，重复步骤S2。

可选地，所述步骤S4具体包括如下步骤：

步骤S401、在每次仿真中，所有线路都有可能随机开断，如果有线路开断则进入快动态过程；

步骤S402、对断线后的新的网络采用支路潮流模型进行一次潮流优化，所述支路潮流模型包括：交流节点有功功率约束、交流节点无功功率约束、交流线路有功损耗约束、交流线路无功损耗约束、交流线路两端电压关系约束、交流网络相角差关系约束、换流节点功率平衡约束、直流线路两端电压关系约束、换流站处功率耦合关系约束以及换流站传输容量约束；

步骤S403、优化调度完成后判断是否有线路达到重载，若有线路达到重载，进入步骤S5；否则，进入步骤S6。

可选地，所述交流节点的有功功率约束为：

所述交流节点的无功功率约束为：

其中，P_Gi和Q_Gi分别是第i个发电机发出的有功功率和无功功率，N_lAC是交流线路总条数，P_cvi和Q_cvi分别是第i个换流站注入的有功功率和无功功率，P_rACj和Q_rACj分别是第j条线路末端接收到的有功功率和无功功率，P_lsACj和Q_lsACj分别是第j条线路上损失的有功功率和无功功率，B_ii是第i个交流节点的对地总电纳，W_ACi是第i个交流节点处的电压的平方，P_cuti和Q_cuti分别是因功率不平衡所需要的切除的有功负荷量和无功负荷量，M_PQAC(i,j)是根据节点i和支路j的关系而确定的线路末端接收到的有功功率和无功功率的系数取值，M_lAC(i,j)是根据节点i和支路j的关系而确定的线路有功损耗和无功损耗的系数取值，P_Di是第i个交流节点处的有功负荷，Q_Di是第i个交流节点处的无功负荷。

可选地，所述交流线路有功损耗约束为：

所述交流线路无功损耗约束：

其中，W_rACj代表线路j末端电压的平方，R_lACjj、X_ljj分别是第j条交流线路的电阻和电抗。

可选地，所述交流线路两端电压关系约束为：

所述交流网络相角差关系约束为：

其中，矩阵M_WAC为M_PQAC的转置矩阵，为网络的节点关联矩阵，C_kj是基本回路矩阵，_nC是基本回路数。

可选地，所述换流节点功率平衡约束为：

其中，矩阵M_PDC表示根据节点i和支路j的关系而确定的直流线路末端接收到的有功功率的系数取值，矩阵M_lDC表示根据节点i和支路j的关系而确定的直流线路有功损耗的系数取值，P_rDCj为第j条直流线路末端输送的功率；P_lsDCj为第j条直流线路上的有功损耗，P_cvm为换流站注入有功功率，N_lDC为直流线路总条数。

可选地，所述直流线路两端电压关系约束为：

其中，矩阵M_WDC为M_PDC的转置矩阵，R_DCjj是第j条直流线路的电阻，W_DCi为第i个直流节点处的电压的平方。

可选地，所述换流站处有无功功率耦合关系约束为：

P_cvm＝|Q_cvm tanα|

其中，Q_cvm为换流站吸收无功功率，角度α为一固定值，由实际运行方式决定。

可选地，所述换流站传输容量约束为：

其中，S_{max cvm}表示第m个换流站的传输容量上限。

可选地，所述步骤S402、对断线后的新的网络采用支路潮流模型进行一次潮流优化，具体包括如下步骤：

若断开的是交流线路，则固定直流线路传输功率不变，将直流线路视为交流***的固定负荷，确定交流线路断开后交直流混联电网的潮流；若断开的是直流线路，则相应换流站注入交流***的功率为0，换流站处的无功功率注入仍存在，固定其他直流线路传输功率不变，将其他直流线路视为交流***的固定负荷，确定直流线路断开后交直流混联电网的潮流；如果没有新的线路开断，则退出快动态过程。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有以下有益效果：

(1)本发明提供的方法所采用的交直流支路潮流模型能够计及网络中的无功功率以及有功和无功损耗，且并未固定电压，适用于交流***且能够拓展到交直流混联***。

(2)本发明提供的方法能够考虑直流线路的控制，在优化调度的过程中，认为直流线路传输功率可以调节，同时认为直流线路可能因为重载断开，进一步引发故障；而且直流线路过载断开后，在下一次仿真时，其容量可以得到更新。

(3)本发明提供的方法考虑换相失败保护，认为直流保护是在调度之后响应而断开线路或闭锁换流器。保护动作的判据为换流器交流测电压低于给定阈值，低于该阈值则换流器闭锁，且不考虑故障时间，即认为调度在整定时间内完成。

附图说明

图1是本发明提供的交直流混联电网连锁故障仿真方法的整体流程图；

图2是本发明实施例中修改后的IEEE30节点***的拓扑图；

图3是本发明实施例提供的修改后的IEEE30节点***停电规模分布曲线图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

本发明提供了一种交直流混联电网的连锁故障仿真方法，其目的在于通过基于交直流潮流模型的潮流优化计算，得到交直流混联电网的故障次数、连锁故障步数、连锁故障的规模和风险大小，由此为交直流电力***的连锁故障风险分析和故障应对以及运行规划提供参考依据，有助于减少交直流混联电力***的故障风险，提高电力***运行的安全可靠性。

如图1所示，本发明提供的交直流混联电网连锁故障仿真方法包括以下流程：

步骤S1、收集所在交直流混联电网的所有技术参数，计算电网初始潮流。

S101、交直流混联电网的所有技术参数包括有：容量基准值、节点数据、支路数据、发电机数据、联络线数据、直流线路和换流站数据。

S102、容量基准值baseMVA设置为100MVA。

S103、节点数据包括有交流***节点个数N_bAC，节点类型type，第i个交流节点处的有功负荷P_Di和无功负荷Q_Di，第i个交流节点处的对地电导g_i和对地电纳b_i，第i个交流节点处电压的幅值V_g和相角V_a以及电压上限V_{ACi max}和电压下限V_{ACi min}，直流***节点个数N_bDC，第i个直流节点处的电压上限V_{DCi max}和电压下限V_{DCi min}。

S104、支路数据包括有交流***线路总条数N_lAC,第j条交流线路的电阻R_lACj、电抗X_lj和电纳b_lj，第j条交流线路的额定容量f_lACj；直流***线路总条数N_lDC,第j条直流线路的电阻R_lDCj，第j条直流线路的容量f_lDCj。

S105、发电机数据包括有发电机有功出力P_g和无功出力Q_g，每台发电机的最大无功出力Q_max和最小无功出力Q_min，最大有功出力P_max和最小有功出力P_min。

S106、联络线数据包括各联络线注入有功功率P_inj和无功功率Q_inj。

S107、换流站数据包括各换流站母线节点编号bus_m、换流站额定容量S_cvm、换流站有功无功耦合系数α、换流站注入有功功率P_cvm。

S108、采用MATpower软件包进行潮流计算，输入运行程序得到电网初始潮流状态。

S2、给定仿真天数k，每次仿真相对于前一次仿真，按照一定比例提高负荷水平和发电机容量，按一定比例增加上次仿真发生过故障的线路的容量。

具体地，按比例给定负荷与发电机出力的增长因子r₁、r₂，对于在上一次仿真中因为过载或重载开断的线路进行容量更新，选取线路容量增长因子u以对线路进行改造增加线路的容量。

S201、在一个示例中，给定仿真天数k＝2000天，每次仿真代表一天，随着电力***的发展，***发电能力和负荷水平不断提高，取负荷与发电机出力增长因子r₁＝r₂＝1.00041。

S202、对于上一次仿真时因过载而开断的线路，一般会被电力规划部门认为需要进行改造，按一定比例增加其传输容量，取线路容量增长因子u＝1.005。

S3、每次仿真中，所有线路均可能按照一定概率断开进入快动态过程。

具体地，在每次仿真中，所有线路都有可能因天气等原因而随机开断，可以取随机开断概率t＝0.0007，如果有线路开断则进入快动态过程。

S4、对断线后新的网络进行一次优化，调度完成后判断是否有线路达到重载，如果有则进入步骤S5，否则，进入步骤S6；

S401、对断线后的新的网络采用支路潮流模型进行一次潮流优化。

具体地，在每次仿真中，所有线路均有可能因为天气等原因按照一定的概率随机开断，给定线路随机开断概率t，如果有线路开断则进入快动态过程。快动态过程中首先需要对新的网络进行一次优化，由于直流潮流模型未考虑线路无功功率因此不能应用于交直流混联电网，在此选用基于支路潮流模型的优化方法进行优化。应用支路潮流模型能够计及网络中的无功和损耗，且并未固定电压，适用于交流***且能够拓展到直流***，适合本项目的研究。所述模型具体包括以下约束条件：

(a)交流节点有功功率和无功功率约束：

其中，P_Gi和Q_Gi分别是第i个发电机发出有功和无功功率，P_cvi和Q_cvi分别是第i个换流站注入的有功和无功功率，P_rACj和Q_rACj分别是第j条线路末端接收到的有功和无功功率，矩阵M_PQAC表示根据节点i和支路j的关系而确定的直流线路末端接收到的有功功率的系数取值，矩阵M_lAC表示根据节点i和支路j的关系而确定的直流线路有功损耗的系数取值，P_lsACj和Q_lsACj分别是第j条线路上损失的有功和无功功率，B_ii是第i个交流节点的对地总电纳，W_ACi是第i个交流节点处的电压的平方，P_cuti和Q_cuti分别是因功率不平衡所需要的切除的有功和无功负荷量。

(b)交流线路有功损耗与无功损耗的关系约束：

其中，W_rACj代表线路j末端电压的平方，R_lACjj、X_ljj分别是第j条交流线路的电阻和电抗。现有技术将W_rACj视为不变值，本发明能够考虑W_rACj的变化以及线路的有功损耗和无功损耗。

(c)交流线路两端电压幅值的关系约束和交流网络中每个独立回路的相角差关系约束：

其中，矩阵M_WAC为M_PQAC的转置矩阵，为网络的节点关联矩阵。

其中，独立基本回路的个数N_C＝N_lAC-N_bAC+1；矩阵C为基本回路矩阵，C_kj表示支路j和回路k的关系。

(d)换流节点功率平衡约束：

其中，矩阵M_PDC和M_lDC表示根据节点i和支路j的关系而确定的直流线路末端接收到的有功功率的系数取值，矩阵表示根据节点i和支路j的关系而确定的直流线路有功损耗的系数取值，P_rDCj为第j条直流线路末端输送的功率；P_lsDCj为第j条直流线路上的有功损耗。

(e)直流线路两端电压关系约束：

其中，矩阵M_WDC为M_PDC的转置矩阵，R_DCjj是第j条直流线路的电阻，W_DCi为第i个直流节点处的电压的平方。

(f)换流站处有无功功率耦合关系约束：

P_cvm＝|Q_cvm tanα|

其中，式中绝对值可依据换流器运行方式来消去，角度α为一固定值，由实际运行方式决定。

(g)换流站传输容量约束：

其中，S_{max cvm}表示第m个换流站的传输容量上限。

(h)模型中其他变量的上下限约束：

P_{Gk min}≤P_Gk≤P_{Gk max},Q_{Gk min}≤Q_Gk≤Q_{Gk max}

-f_LACj≤P_rACj≤f_LACj

-f_LDCj≤P_rDCj≤f_LDCj

其中，P_{Gk min}是第k台发电机发出的有功功率P_Gk的最小值；P_{Gk max}是P_Gk的最大值；Q_{Gk min}表示第k台发电机发出的无功功率Q_Gk的最小值；Q_{Gk max}表示Q_Gk的最大值；f_LACj和f_LDCj分别是第j条交流线路和直流线路的线路容量上限，V² _{ACi max}和V² _{ACi min}分别是第i个交流节点的电压最大值平方和最小值平方；V² _{DCi max}和V² _{DCi min}分别是第i个交流节点的电压最大值平方和最小值平方。

S402、优化调度完成后判断是否有线路达到重载，若有线路达到重载，进入步骤S5；否则，进入步骤S6。

具体地，优化调度完成后给定线路重载阈值a判断是否有线路达到重载。在一个示例中，线路重载阈值取为a＝0.99。

S5、对于重载线路按照一定的概率断开，然后判断有无新线路断开，如果有，则形成新的网络重复步骤S4过程，否则，进入步骤S6；

具体地，对于达到重载阈值的线路按照给定的概率b进行开断，然后判断如果有新线路断开，需要划分孤岛并进行孤岛功率平衡优化以形成新的网架然后继续重复步骤S4。若断开的是交流线路，则固定直流线路传输功率为初始值不变，视为交流***的固定负荷，由电网初始潮流状态计算交流线路断开后交直流***的潮流；若断开的是直流线路，则相应换流站注入交流***的功率为0，换流站处的无功补偿装置提供的无功功率注入仍存在，固定其他直流线路传输功率为初始值不变，视为交流***的固定负荷，由电网初始潮流状态计算直流线路断开后交直流***的潮流。如果没有新的线路开断，则退出快动态过程。

在一个示例中，取重载线路开断概率b＝0.999。

S6、判断仿真天数k是否达到设定最大值，如果达到，则结束连锁故障仿真过程，否则，重复步骤S2过程。

具体地，如果达到最大值，则退出连锁故障仿真过程。此时，可根据连锁故障仿真结果统计故障次数、最多断线次数、最多连锁故障步数、不同置信度下的故障规模等。

其中，在断线后需要求解潮流优化问题，本发明采用基于支路潮流模型的优化模型进行潮流优化求解。当交直流***网架参数已知后，考虑到其交流线路的容量限制、换流站交流侧的电压限制等安全约束，以及功率平衡、发电机有无功出力范围、换流站注入功率范围等运行约束，能够对断线后的交直流电力***进行潮流优化，然后形成新的网架。

以下公开详细的示范实施例，此处公开的具体结构和功能细节仅仅是处于描述示范实施例的目的。然后应当理解，本发明不局限于公开的具体示范实施例，而是覆盖落入本公开范围内的所有修改、等同物和替换物。在对全部附图的描述中，相同的附图标记表示相同的元件。

图2为本发明所选取的IEEE30节点测试***的***图，以此为例进行实例说明。实例中，在节点2和节点8上分别增加一台换流器，在其间增加一条直流线路，形成交直流互联电力***。收集所在交直流互联电网常规技术参数，同时给出慢动态过程和快动态过程的各仿真参数值。

如图2所示的交直流互联电力***包含交流节点个数N_bAC＝30，交流线路条数N_lAC＝41，换流站个数N_cv＝2，直流***节点个数N_bDC＝2，直流***线路总条数N_lDC＝1。

具体地，交流节点参数如表1：

表1交流***节点参数

交流线路参数如表2：

表2交流***线路参数

发电机参数如表3：

表3发电机参数

换流器参数见表4：

表4换流器参数

交流节点	直流节点	容量(MW)
			2	1	100
8	2	100

直流节点参数见表5：

表5直流***节点参数

编号	电压上限(p.u.)	电压下限(p.u.)	基准电压(kV)
				1	1.06	0.94	320
2	1.06	0.94	320

直流线路参数见表6：

表6直流线路参数

首端节点	末端节点	线路电阻(p.u.)	线路容量(MVA)
				1	2	0.01	100

此外，考虑实际运行情况，取换流站的有功无功耦合系数α＝0.8。取直流线路输送的功率为27MW，通过MATpower计算可得初始的交直流***潮流状态，逆变站所在的交流节点8的电压为0.9538p.u。

给定仿真天数k为2000天，取线路容量增长因子u＝1.005，负荷与发电机增长因子r₁＝r₂＝1.00041，线路重载阈值a＝0.99，重载线路开断概率b＝0.999，因天气等原因线路随机开断概率t＝0.0007。

得到连锁故障的计算总体结果如下表7所示：

表7连锁故障仿真总体结果

由表7连锁故障仿真总体结果能够得到交直流混联***在所有规定仿真天数中的随机故障次数和重载故障次数，某个置信度下的最大切负荷率以及最多断线条数和连锁故障步数。

故障次数最多的前十条线路如下表8所示：

表8故障次数排序表

不同置信度下的最大故障规模如表9所示：

表9不同置信度下最大故障规模

表8和表9分别给出了交直流混联***在规定仿真天数中故障次数最多的前十条线路和不同置信度下的最大故障规模。表7～表9所给出的数据能够为交直流电力***的连锁故障风险和规划运行提供参考依据，促使我们能够及时采取措施对过载和重载线路进行扩容，降低连锁故障的发生风险。

图3给出了该交直流混联***连锁故障仿真后的停电规模分布曲线图，其中横轴是以被切除负荷量x为量度的停电规模，纵轴是停电规模分布。从图3可知，改进的交直流***连锁故障模型的停电规模分布近似呈现幂律特性，符合实际情况，能够间接证明该连锁故障仿真方法的正确性。

本发明公开一种交直流混联电网连锁故障仿真方法，属于电气工程领域。该方法分为快动态和慢动态两个迭代过程，快动态过程用于模拟电力***连锁故障，而慢动态过程则用于模拟电力***的演化。快动态过程的思想是：当某条线路因天气等原因随机断开后，其余线路可能会因重载或过载而相继开断从而导致电网连锁故障；慢动态过程的思想是：随着电力***发展，***发电能力和负荷水平不断提高，线路潮流相应增加，当其达到线路传输容量极限时会导致线路开断，而这些因过载而开断的线路，一般会被电力规划部门认为需要进行改造，如增加其传输容量，以提高输电***的安全性。

通过本发明的连锁故障仿真方法，能够将直流线路模型和直流保护措施考虑进来，适用于交直流混联电网的连锁故障仿真，能够为交直流电力***的故障风险和运行规划提供参考依据，从而能够采取措施降低连锁故障发生的风险，提高电力***运行的安全可靠性。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种交直流混联电网连锁故障仿真方法，其特征在于，包括以下步骤：

步骤S1、确定交直流混联电网的交流线路、直流线路、联络线、发电机、换流站以及负荷；发电机容量和负荷水平可用于计算电网的潮流，所述潮流包括各条线路的有功功率、无功功率、有功损耗、无功损耗以及电网各母线节点的电压幅值和相角；

步骤S2、对交直流混联电网进行预设次数的循环仿真，每次仿真相对于前一次仿真，按照第一预设比例提高发电机容量和负荷水平，按照第二预设比例增加前一次仿真中发生过故障的线路的容量；

步骤S3、每次仿真中，所有线路按照第一预设概率断开；

步骤S4、若交直流混联电网有线路断开，则对断线后的交直流混联电网采用支路潮流模型进行潮流优化，将断开线路的功率分配到未断线的线路上，优化后判断是否有线路达到重载，若有线路达到重载，则进入步骤S5，否则，进入步骤S6；

步骤S5、对重载线路按照第二预设概率断开，然后判断有无新线路断开，如果有，则重复步骤S4，否则，进入步骤S6；

步骤S6、判断仿真次数是否达到预设次数，如果达到，则结束连锁故障仿真过程，否则，重复步骤S2；

其中，所述步骤S4具体包括如下步骤：

步骤S401、在每次仿真中，所有线路都有可能随机开断，如果有线路开断则进入快动态过程；

步骤S402、对断线后的新的网络采用支路潮流模型进行一次潮流优化，所述支路潮流模型包括：交流节点有功功率约束、交流节点无功功率约束、交流线路有功损耗约束、交流线路无功损耗约束、交流线路两端电压关系约束、交流网络相角差关系约束、换流节点功率平衡约束、直流线路两端电压关系约束、换流站处功率耦合关系约束以及换流站传输容量约束；

步骤S403、优化调度完成后判断是否有线路达到重载，若有线路达到重载，进入步骤S5；否则，进入步骤S6。

2.根据权利要求1所述的交直流混联电网连锁故障仿真方法，其特征在于，所述交流节点的有功功率约束为：

所述交流节点的无功功率约束为：

其中，P_Gi和Q_Gi分别是第i个发电机发出的有功功率和无功功率，N_lAC是交流线路总条数，P_cvi和Q_cvi分别是第i个换流站注入的有功功率和无功功率，P_rACj和Q_rACj分别是第j条线路末端接收到的有功功率和无功功率，P_lsACj和Q_lsACj分别是第j条线路上损失的有功功率和无功功率，B_ii是第i个交流节点的对地总电纳，W_ACi是第i个交流节点处的电压的平方，P_cuti和Q_cuti分别是因功率不平衡所需要的切除的有功负荷量和无功负荷量，M_PQAC(i,j)是根据节点i和支路j的关系而确定的线路末端接收到的有功功率和无功功率的系数取值，M_lAC(i,j)是根据节点i和支路j的关系而确定的线路有功损耗和无功损耗的系数取值，P_Di是第i个交流节点处的有功负荷，Q_Di是第i个交流节点处的无功负荷。

3.根据权利要求2所述的交直流混联电网连锁故障仿真方法，其特征在于，所述交流线路有功损耗约束为：

所述交流线路无功损耗约束：

其中，W_rACj代表线路j末端电压的平方，R_lACjj、X_ljj分别是第j条交流线路的电阻和电抗。

4.根据权利要求3所述的交直流混联电网连锁故障仿真方法，其特征在于，所述交流线路两端电压关系约束为：

所述交流网络相角差关系约束为：

其中，矩阵M_WAC为M_PQAC的转置矩阵，为网络的节点关联矩阵，C_kj是基本回路矩阵，_nC是基本回路数，N_bAC为交流***节点个数。

5.根据权利要求3所述的交直流混联电网连锁故障仿真方法，其特征在于，所述换流节点功率平衡约束为：

其中，矩阵M_PDC表示根据节点i和支路j的关系而确定的直流线路末端接收到的有功功率的系数取值，矩阵M_lDC表示根据节点i和支路j的关系而确定的直流线路有功损耗的系数取值，P_rDCj为第j条直流线路末端输送的功率；P_lsDCj为第j条直流线路上的有功损耗，P_cvm为换流站注入有功功率，N_lDC为直流线路总条数。

6.根据权利要求5所述的交直流混联电网连锁故障仿真方法，其特征在于，所述直流线路两端电压关系约束为：

其中，矩阵M_WDC为M_PDC的转置矩阵，R_DCjj是第j条直流线路的电阻，W_DCi为第i个直流节点处的电压的平方,N_bDC为直流***节点个数。

7.根据权利要求5所述的交直流混联电网连锁故障仿真方法，其特征在于，所述换流站处有无功功率耦合关系约束为：

P_cvm＝|Q_cvm tanα|

其中，Q_cvm为换流站吸收无功功率，角度α为一固定值，由实际运行方式决定。

8.根据权利要求7所述的交直流混联电网连锁故障仿真方法，其特征在于，所述换流站传输容量约束为：

其中，S_maxcvm表示第m个换流站的传输容量上限。

9.根据权利要求1至8任一项所述的交直流混联电网连锁故障仿真方法，其特征在于，所述步骤S402、对断线后的新的网络采用支路潮流模型进行一次潮流优化，具体包括如下步骤：

若断开的是交流线路，则固定直流线路传输功率不变，将直流线路视为交流***的固定负荷，确定交流线路断开后交直流混联电网的潮流；若断开的是直流线路，则相应换流站注入交流***的功率为0，换流站处的无功功率注入仍存在，固定其他直流线路传输功率不变，将其他直流线路视为交流***的固定负荷，确定直流线路断开后交直流混联电网的潮流；如果没有新的线路开断，则退出快动态过程。

Images