1.外网等值网络边界电压无功支撑充裕性计算方法，其特征在于：所述方法的具体步骤包括以下内容；
1)确定m时刻数
依据内外网相连的端口数，通过公式(1)确定求解等值网络参数所需的PMU时刻数；
2mn＞n2+6n (1)
式中，n为内外网相连的端口数，m为待求的PMU采样时刻数；
2)输入测量数据
输入m个时刻各边界节点处PMU的电压量测值和m个时刻各边界节点向内网注入的等效电流量测值其中i＝1,2,...,n,n为内外网相连的端口数，t＝t1,t2,...,tm,t为PMU采样时刻，m为总的PMU采样时刻数；
3)建立等值网络测量方程
根据第2)步输入的电压量测量与等效量测量，建立等值网络的量测方程和
f 1 i t ( x ) = α [ E i , Re - U i , Re t - Z i , Re ( I i , Re - B i U i , Im t ) + Z i , Im ( I i , Im + B i U i , Re t ) ] - β [ E i , Im - U i , Im t - Z i , Re ( I i , Im + B i U i , Re t ) - Z i , Im ( I i , Re - B i U i , Im t ) ] - Σ j = 1 , j ≠ i n [ η ( U i , Re t - U j , Re t ) - μ ( U i , Im t - U j , Im t ) ] = 0 - - - ( 2 ) ]]>
f 2 i t ( x ) = α [ E i , Im - U i , Im t - Z i , Re ( I i , Im + B i U i , Re t ) - Z i , Im ( I i , Re - B i U i , Im t ) ] + β [ E i , Re - U i , Re t - Z i , Re ( I i , Re - B i U i , Im t ) + Z i , Im ( I i , Im + B i U i , Re t ) ] - Σ j = 1 , j ≠ i n [ η ( U i , Im t - U j , Im t ) + μ ( U i , Re t - U j , Re t ) ] = 0 - - - ( 3 ) ]]>
与分别表示t时刻第i个边界节点处量测方程的实部与虚部，x＝[Ei,Re,Ei,Im,Zi,Re,Zi,Im,Zij,Re,Zij,Im,SLi,Re,SLi,Im,Bi]T,i＝1,2,...,n，j＝1,2,...,n,j≠i，x为待求的等值网络参数矩阵，i代表第i个边界节点，j代表第j个边界节点，n为内外网相连的端口数；其中Ei,Re和Ei,Im分别为第i个边界节点处所对应等值虚拟发电机节点电压实部与虚部；Zi,Re和Zi,Im分别为第i个边界节点处所对应等值支路的电阻和电抗；Zij,Re和Zij,Im分别为第i个边界节点与第j个边界节点之间等值支路的电阻和电抗；SLi,Re和SLi,Im分别为第i个边界节点处的负荷对应的等值电流的实部和虚部；Bi为第i个边界节点处对应的对地支路；和分别为第i个边界节点的电压实部和虚部；和分别为第j个边界节点的电压实部和虚部，t表示第t个时刻；Ii,Re和Ii,Im分别为第i个边界节点处负荷等值电流与边界节点向内网注入的等效电流之和的实部和虚部；α和β分别代表了的实部和虚部，η和μ分别代表了的实部和虚部，其中Zij为边界节点i与边界节点j之间的等值阻抗，Zik为边界节点i与边界节点k之间的等值阻抗，k＝1,2,...,n，k≠i,j；
4)建立等值参数约束方程
公式(4)表示为输入的最大运行方式下等值网络的等值参数向量xmaxo；
xmaxo＝[Ei,maxo,Ri,maxo,Xi,maxo,Rij,maxo,Xij,maxo,PLi,maxo,QLi,maxo,Bi,maxo]T (4)
公式(5)表示为输入的最小运行方式下的等值网络的等值参数向量xmino；
xmino＝[Ei,mino,Ri,mino,Xi,mino,Rij,mino,Xij,mino,PLi,mino,QLi,mino,Bi,mino]T (5)
依据最大最小运行方式等值参数建立公式(6)-(13)的等值参数约束方程；
Bi,mino≤Bi'≤Bi,maxo (6)
Ri,maxo≤Zi,Re'≤Ri,mino (7)
Xi,maxo≤Zi,Im'≤Xi,mino (8)
Rij,maxo≤Zij,Re'≤Rij,mino (9)
Xij,maxo≤Zij,Im'≤Xij,mino (10)
PLi,mino≤SLi,Re'≤PLi,maxo (11)
QLi,mino≤SLi,Im'≤QLi,maxo (12)
E i , m i n ≤ E i , Re ′ 2 + E i , Im ′ 2 ≤ E i , m a x - - - ( 13 ) ]]>
其中Bi,mino和Bi,maxo分别为最小运行方式下和最大运行方式下第i个边界节点所对应的等值对地支路；Bi'为当前待求的第i个边界节点处对应的对地支路；Ri,mino和Ri,maxo分别为最小运行方式下和最大运行方式下第i个边界节点所对应等值支路的电阻；Zi,Re'为当前待求的第i个边界节点处对应等值支路的电阻；Xi,mino和Xi,maxo分别为最小运行方式下和最大运行方式下第i个边界节点所对应等值支路的电抗；Zi,Im'为当前待求的第i个边界节点处对应等值支路的电抗；Rij,mino和Rij,maxo分别为最小运行方式下和最大运行方式下第i个边界节点与第j个边界节点之间等值支路的电阻；Zij,Re'为当前待求的第i个边界节点与第j个边界节点之间等值支路的电阻；Xij,mino和Xij,maxo分别为最小运行方式下和最大运行方式下第i个边界节点与第j个边界节点之间等值支路的电抗；Zij,Im'为当前待求的第i个边界节点与第j个边界节点之间等值支路的电抗；Ei,min和Ei,max分别为最小运行方式下和最大运行方式下第i个边界节点处所对应等值虚拟发电机节点电压；Ei,Re'和Ei,Im'分别为待求的第i个边界节点处所对应等值虚拟发电机节点电压实部和虚部；PLi,mino和PLi,maxo分别为最小运行方式下和最大运行方式下第i个边界节点处的负荷有功功率；QLi,mino和QLi,maxo分别为最小运行方式下和最大运行方式下第i个边界节点处的负荷无功功率；SLi,Re'和SLi,Im'分别为待求的第i个边界节点处的负荷对应的等值负荷的有功功率和无功功率；
5)建立等值网络的优化模型
利用公式(14)建立等值网络优化模型的目标函数；
m i n ( Σ t = 1 m Σ i = 1 n [ ( f 1 i t ( x ) ) 2 + ( f 2 i t ( x ) ) 2 ] ) - - - ( 14 ) ]]>
式中m为PMU采样时刻数，n为内外网相连的端口数；与分别表示t时刻第i个边界节点处量测方程的实部与虚部；
6)求解外网等值参数
通过内点法直接求得外网等值参数x；
7)等值外网的无功优化模型
根据第6)步中的外网等值参数建立等值外网的无功优化模型；
建立外网无功损耗为最小的目标函数；
minfE(QG,Qb) (15)
式中，QG为等值发电机无功功率向量，Qb为等值对地支路无功功率向量；
建立无功优化模型的功率平衡等式方程；
P a - U a Σ b = 1 N U b ( G a , b cosδ a , b + B a , b sinδ a , b ) = 0 - - - ( 16 ) ]]>
Q a - U a Σ b = 1 N U b ( G a , b sinδ a , b - B a , b cosδ a , b ) = 0 - - - ( 17 ) ]]>
式中，a、b为等值外网节点，a＝1,2,...,N,b＝1,2,…,N，N为等值外网节点总数，Pa和Qa分别为节点a的注入有功功率与无功功率；Ua、δa分别为节点a处的电压幅值和相角，Ub、δb分别为节点b处的电压幅值和相角，式中δa,b＝δa-δb；Ga,b、Ba,b分别为节点导纳矩阵的第a行、第b列项元素的实部和虚部，sin与cos分别为正弦函数与余弦函数；
建立无功优化模型的不等式约束方程；
Q G e min ≤ Q G e ≤ Q G e max , e = 1 , 2 , ... , N B - - - ( 18 ) ]]>
U f min ≤ U f ≤ U f max , f = 1 , 2 , ... , N - - - ( 19 ) ]]>
式中，为节点e处等值发电机发出的无功功率，e＝1,2,...,NB，NB为等值发电机总数，和为节点e处等值发电机发出的无功功率的上下限；Uf为节点f处的电压幅值，f＝1,2,...,N，N为等值外网节点总数，和为节点f处的电压幅值上下限；用内点法求解上述非线性优化问题可求得Qb，其中Qbi∈Qb；Qbi为节点处等值对地支路无功功率；
8)求取无功功率上下限
通过外网最大最小运行方式下的等值对地支路可求取等值对地支路无功功率上下限；
Qbi,min＝Ui2Bi,mino (20)
Qbi,max＝Ui2Bi,maxo (21)
式(20)与(21)中，Bi,maxo和Bi,mino分别为第i个边界节点处最大和最小运行方式下的等值对地支路；Ui为第i个边界节点处的电压幅值；Qbi,max和Qbi,min分别为第i个边界节点等值对地支路无功功率上下限；
9)等值方法的选择
如果Qbi＜Qbi,max，则边界节点处有充足的电压无功支撑能力，外网用简单的PV等值方法等值；如果Qbi＞Qbi,max，则边界节点处没有充足的电压无功支撑能力，外网只能采用考虑元件全面性的外网等值方法。
