1.一种抑制多机电力系统暂态频率下降率的三相重合时序整定方法，其特征在于：当
多机电力系统的交流输电线路发生三相短路故障时，通过修改系统正常状态下的导纳矩
阵，得到故障线路分别由线路首、末端首先三相重合时的导纳矩阵；将导纳矩阵写成发电机节点和负荷节点构成的电导分块矩阵和电纳分块矩阵，从而计算得到故障线路首、末端分别重合时的雅克比矩阵和同步功率系数矩阵；据此计算出两种重合时序下多机电力系统的暂态频率下降率，比较两种重合时序下计算出的暂态频率下降率，选取暂态频率下降率较小者对应的三相重合时序为优化的重合方案。
2.根据权利要求1所述的抑制多机电力系统暂态频率下降率的三相重合时序整定方
法，其特征在于具体步骤为：
(1)构成多机系统的导纳矩阵YB：
根据多机系统的网络接线结构，采用节点导纳法形成系统的节点导纳矩阵：
若多机系统的网络节点数为m，描述它的导纳矩阵为m·m阶矩阵YB；
导纳矩阵YB各对角元素，即各节点的自导纳Yii等于相应节点所连支路的导纳之和：
其中yij为节点i与节点j间所连支路的导纳；
导纳矩阵YB各非对角元素的互导纳Yij就等于任意两个节点之间所连接导纳的负值：Yij
＝-yij；
(2)模拟线路首端投入三相重合闸：
1)修改导纳矩阵YB中故障线路对应元素，模拟线路首端投入三相重合闸，即故障线路首
端断线消失，得到的导纳矩阵记作Y，具体方法如下：
若故障线路首端断路器母线侧和线路侧的节点编号分别为p、q，则导纳矩阵中p、q点对
应自导和互导纳按式(1)修改：

式中，ypp、yqq为节点p、q的自导纳，ypq为节点p、q间的互导纳，y0为模拟p、q点间的故障消失而采用的导纳，其标么值取9999，单位pu，导纳矩阵Y如式(2)所示：
Y＝G+jB   (2)
式中，电导G为导纳矩阵Y的实部，为电导矩阵；电纳B为导纳矩阵Y的虚部，为电纳矩阵；
分别将电导矩阵G、电纳矩阵B按照发电机节点和负荷节点列写成分块矩阵的形式，如式(3)所示：

式中，GGG为各发电机节点的自电导矩阵，GGD为发电机节点与负荷节点之间的互电导矩
阵，GDG为负荷节点与发电机节点之间的互电导矩阵，GDD为各负荷节点的自电导矩阵；BGG为各发电机节点的自电纳矩阵，BGD为发电机节点和负荷节点之间的互电纳矩阵，BDG为负荷节点和发电机节点之间的互电纳矩阵，BDD为负荷节点的自电纳矩阵；
2)计算偏导数矩阵
将式(3)所示的电导分块矩阵中的GGD、GDD矩阵和电纳分块矩阵中的BGD、BDD矩阵代入式(4)和式(5)中，分别计算有功功率对转子角的偏导数矩阵HGD和有功功率对电压相角的偏导数矩阵HDD；


式中，PG为各发电机的有功功率PGi组成的矩阵，PGi表示第i台发电机的有功功率，δD为各负荷节点的电压相量相角δDi组成的矩阵，δDi表示第i个负荷节点的电压相量相角；PD为各负荷节点的有功功率PDi组成的矩阵，PDi表示第i个负荷节点的有功功率，θD为各负荷节点的电压相量相角θDi组成的矩阵，θDi表示第i个负荷节点电压相量相角，UG为各发电机节点的电压幅值UGi组成的矩阵，UGi表示第i台发电机的电压幅值；UD为各负荷节点的电压幅值UDi组成的矩阵，UDi表示第i个负荷节点的电压幅值，δGD为发电机节点电压相量和负荷节点电压相量之间的相角差组成的矩阵；θDD为任意两个负荷节点电压相量之间的相角差组成的矩阵；
3)计算同步功率系数矩阵L：
将式(4)和(5)计算得到的HGD、HDD代入式(6)中，计算得到发电机内节点与负荷节点之间的同步功率系数矩阵L：

4)计算多机电力系统暂态角频率下降率：
由式(6)中计算得到的矩阵L代入式(7)中，计算得到多机系统暂态角频率下降率Δω：

式中，τ为所有发电机的转动惯量之和，即 Mi为第i台发电机的转动惯量；1n为
n维的单位列向量， 为1n的转置，ΔΡm0为原动机功率扰动量，ΔΡD0为负荷功率扰动量；
5)计算多机电力系统的暂态频率下降率：
由于角频率ω和频率f之间满足式(8)，将计算得到的多机电力系统暂态角频率Δω代
入式(9)中，计算首端三相重合时多机电力系统的暂态频率下降率Δf：
f＝ω/2π   (8)
即：
(3)模拟线路末端投入三相重合闸
1)修改导纳矩阵YB中故障线路对应元素，模拟线路末端投入三相重合闸，即故障线路末
端断线状况消失，得到导纳矩阵Y′，具体方法如下：
若故障线路末端断路器对应母线侧和线路侧的节点编号为p′、q′，则导纳矩阵中p′、q′点对应自导纳和互导纳按式(10)进行修改：

式中，yp′p′、yq′q′为节点p′、q′的自导纳，yp′q′为节点p′、q′间的互导纳，y0为模拟p′、q′点间的故障消失而采用的导纳，其标么值取9999，单位pu，导纳矩阵Y′如式(11)所示：
Y′＝G′+jB′   (11)
式中，电导G′为导纳矩阵Y′的实部，为电导矩阵；电纳B′为导纳矩阵Y′的虚部，为电纳矩阵，分别将电导矩阵G′、电纳矩阵B′按照发电机节点和负荷节点列写成分块矩阵的形式，如式(12)所示：

式中，G′GG为各发电机节点的自电导矩阵，G′GD为发电机节点与负荷节点之间的互电导矩阵，G′GG为负荷节点与发电机节点之间的互电导矩阵，G′GD为各负荷节点的自电导矩阵；
B′GG为各发电机节点的自电纳矩阵，B′GD为发电机节点和负荷节点之间的互电纳矩阵，B′DG为负荷节点和发电机节点之间的互电纳矩阵，B′DD为负荷节点的自电纳矩阵；
2)重复步骤(2)中的2)—5)，将式(12)所示的电导分块矩阵中的G′GD、G′DD矩阵和电纳分块矩阵中的B′GD、B′DD矩阵代入式(4)和式(5)中，分别计算有功功率对转子角的偏导数矩阵H′GD和有功功率对电压相角的偏导数矩阵H′DD；将H′GD、H′DD代入式(6)中计算得到发电机内节点与负荷节点之间的同步功率系数矩阵L′；将L′代入式(7)中，计算得到多机系统暂态角频率下降率Δω′；将Δω′代入式(9)中，计算得到末端三相重合时多机电力系统暂态频率下降率Δf′；
(4)得出结论：
比较线路首端三相首先重合和末端三相首先重合时多机电力系统暂态频率下降Δf、
Δf′，选取较小者对应的三相重合时序作为输出的优化重合时序方案；
若Δf＞Δf′，则由线路末端首先三相重合；
反之，由线路首端三相首先重合。