CN101877482B - 基于慢模式特征值对于线路参数灵敏度的弱连接识别方法 - Google Patents
基于慢模式特征值对于线路参数灵敏度的弱连接识别方法 Download PDFInfo
- Publication number
- CN101877482B CN101877482B CN2009102383341A CN200910238334A CN101877482B CN 101877482 B CN101877482 B CN 101877482B CN 2009102383341 A CN2009102383341 A CN 2009102383341A CN 200910238334 A CN200910238334 A CN 200910238334A CN 101877482 B CN101877482 B CN 101877482B
- Authority
- CN
- China
- Prior art keywords
- generator
- node
- sensitivity
- circuit
- line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明公开了属于电力系统分析和控制技术领域的一种基于慢模式相关特征值对于线路参数灵敏度的弱连接识别方法,其特征在于:首先计算系统的慢模式相关特征值组对各条线路参数变化灵敏度,再选取灵敏度数值比较大的线路作为系统的弱连接。主要按照以下几个步骤进行:计算机初始化,输入包括系统网络参数、节点电压和功率信息的初始化数据;运用牛顿拉夫逊法计算系统潮流;形成系统线性化的动态方程;运用QR方法求取系统特征值,并根据大小进行排序,选取系统慢模式相关特征值组;计算慢模式相关特征值对于系统各条线路参数变化的灵敏度,并按大小进行排序;选取系统弱连接。本发明的有益效果是相比于以前首先处理节点的方法,能够直接处理线路。
Description
技术领域
本发明属于电力系统分析和控制技术领域,特别涉及一种基于慢模式特征值对于线路参数灵敏度的弱连接识别方法。
背景技术:
解列是电力系统受到大干扰,在系统完整性无法保持的情况下将系统分成几个孤岛的控制。以在线计算、全局控制、协调动作为特点的主动解列控制是一个重要的研究方向。基于电力系统慢同调的主动解列思路认为,如果解列操作发生在同调机群之间的弱连接上对于系统的冲击最小,解列后孤岛更容易存活。基于这种思路,弱连接识别在基于电力系统慢同调的主动解列控制中占据着举足轻重的地位。对于弱连接识别,以前的做法是首先对节点进行同调分群,各个同调群之间的连接就是弱连接,但是这种做法有以下缺陷:
1)负荷节点处理会遇到困难。尽管引入了广义特征值概念,能够将负荷节点加以考虑,但是仍然面临着发电机节点和负荷节点对系统慢模式特征值的灵敏度在数值上差别很大以及阀值难以选取的困难;
2)这种基于节点灵敏度的弱连接识别方法直接对节点进行操作,并没有直接处理线路,“弱连接”的概念并不能直接地显现出来。但是,解列最后操作的对象是线路,这种通过首先处理节点再确定线路的做法并不能直观反映出所挑选线路的特性。
本发明的理论基础就是电力系统慢同调理论,电力系统中同调的定义有多种,本发明采用的是一种模式不可观性定义。
同调定义:对于系统线性x=[x1 x2...xn]T为其的状态变量组, (i=1,...,n)为系统第i个状态变量对时间的一阶导数,即σs={λ1,λ2,..λs}为系统的一组特征值模式,如果x中的两个状态变量xi(t),xj(t)满足yij(t)=xi(t)-xj(t)不含σr中的任何一个模式分量,那么就称xi(t),xj(t)关于模式组σr是同调的。如果系统的状态变量x=[x1 x2...xn]T可以分成S组,每一组内的任何两个变量之间都是关于模式组σs同调的,那么就称这个系统是关于模式组σr同调的。同时如果σs={λ1,λ2,..λs}是系统的一组绝对值最小的特征值,那么就说系统是严格慢同调的。
这种同调定义直接应用到电力系统还有障碍,发电机取经典二阶模型同时忽略阻尼的情况下,电力系统的线性化动态方程的形式为:
Δδ=[Δδ1 Δδ2...ΔδnG]T为系统所有发电机功角在平衡点的偏差列向量,nG为系统发电机的台数;为系统所有发电机角速度在平衡点的偏差列向量; 下文中,某个列向量的上标有一个黑点“·”都表示此列向量中每一个元素对时间一阶导数组成的列向量。
可以进一步忽略角速度项:
其中:A=M-1K,M=diag{M1,M2,...,Mn}是由所有发电机惯性时间常数组成的对角矩阵,Mi=2Hi(i=1,2,..,n),K是系统的网络矩阵;下文中某个列向量上标有两个黑点“··”都表示此列向量中每一个元素对于时间的两阶导数组成的列向量。
这是一个二阶微分方程的形式,在一阶微分方程形式下的同调定义应用到二阶微分方程的形式要做相应的修改。同样经过推导可以得到如下结论:A的特征值组σs={λ1,λ2,..λs}并不是系统真正的机电振荡模式,而对应的真正的机电振荡模式是为了与模式特征值相区别,称σs={λ1,λ2,..λs}中的每个特征值为模式相关特征值。进一步,如果σs={λ1,λ2,..λs}是A的绝对值最小的一组特征值,那么称σs={λ1,λ2,..λs}为慢模式相关特征值组。如果说A关于慢模式相关特征值组σs={λ1,λ2,..λs}是慢同调的,意思就是系统关于模式特征值组是慢同调的。
进一步得到在形式下系统的同调定义:对于系统线性X=[x1 x2...xn]T为其的状态变量, (i=1,...,n)为系统第i个状态变量对时间的二阶导数,即σs={λ1 λ2...λs}为系统的一组特征值,如果x中的两个状态变量xi(t),xj(t)满足yij(t)=xi(t)-xj(t)不含中的任何一个模式分量,那么就称xi(t),xj(t)关于特征值组σs是同调的。如果系统的状态变量x=[x1 x2...xn]T可以分成s组,每一组内的任何两个变量之间都是关于模式组σr同调的,那么就称这个系统是关于模式组σs同调的。同时如果σs={λ1 λ2...λs}是系统的一组绝对值最小的特征值,那么就说系统是严格慢同调的。
对于A的一组慢特征值σs={λ1,λ2,..λs},如果系统能够关于σr是严格慢同调的,那么A的对应于σs={λ1,λ2,..λs}的特征相量矩阵V中同调机组对应的行的各个元素都相同,这是一个很重要的结论。
由系统的严格慢同调可以推出一些有益的结论。在系统严格慢同调的条件下,假设系统的s个同调机群已知,每个同调机群选择一个参考状态变量组成X1=[x1 x2...xr]T,称X1为参考变量相量,那么剩下的n-s个状态变量组成X2=[xr+1 xr+2...xn]T,X2为剩余变量相量。那么系统的状态变量相量就按照如下的顺序排列:
下面定义一个分群矩阵Lg,Lg具有以下特征:
1)Lg是一个(n-s)×s维的矩阵,每一行代表每一个待分群的状态(X2=[xs+1 xs+2...xn]T中的某个状态),每一列代表着以某个状态(X1=[x1 x2...xs]T)为参考状态的机群。
2)Lg的第i(i<=n-s)元素按照如下的方式定义:如果状态xr+i属于以xj(j<=s)为参考变量的机群,那么Lg(i,j)=1,否则Lg(i,j)=0。
其中: 而 Ir为s×s维的单位矩阵,M2=diag{Ms+1,Ms+2,...,Mn}。
根据系统严格慢同调条件,经过严格的推导可以得到如下结论:
1)能够实现系统的区域分解和聚合。
其中:
继续上面的推导,K作为网络矩阵具有下列的性质:
也就是说K的非正元素只有对角元素kii,并且它与本行的非对角元素之和的相反数相等。可以知道K的特征值中有0特征值,且对应的特征向量为:u=[1 1,1]T。在已知系统分群的情况下对于系统的状态变量x重新排序,以使得在排序后的状态向量Px中,同一区域的状态量在一起。假设系统分成s个区域,重新排序后的系统矩阵具有如下的形式:
Kαα表示区域内部的联结,而Kαβ(β≠α)表示区域α,β之间的连接。重写KP以使内部和外部连接区分开来:KP=KI+KE。注意:这里面的KI,KE不是简单的取对角块和非对角块的结果。而是KI中减去了对于与区域之间连接的元素,这部分元素加到了KE上。也就是说KE在原来KI对应的对角块位置上的元素并不全是0。
这样得到的KI有如下特征:
这样:
而能够推得:
又由于:
所以,就有:
进而:
又知λ(Aa)=σs={λ1,λ2,..λs},所以可以得到了一个很重要的结论:系统的模式特征值只与系统的同调区域之间的弱连接KE有关,而内连接KI对它没有任何影响。
可以把系统区域之间的弱连接定义成KE所对应的线路,但是我们知道以上的系统模型并不是结构保留模型,而是收缩到发电机节点的网络,所以KE并不对应着实际的物理线路,这样就无法直接通过KE来识别弱连接了。但是,上面的结论还是能够提供很有益的启示:区域之间的弱连接线路对于系统的慢模式特征值影响比较大,而区域内部的线路对于系统的慢模式影响比较小。
发明内容
本发明的目的是提出了一种慢模式相关特征值对于线路参数灵敏度的弱连接识别方法,其特征在于,具体算法步骤如下:
步骤(1):计算机初始化;
1.1)输入以下参数:
电力系统的网络参数,包括:输电线路的串联电阻、串联电抗、并联电导和并联电纳;变压器的变比和阻抗;并联在输电线路上的电容器和电抗器的阻抗;发电机的动态参数包括发电机的惯性时间常数;发电机的暂态电抗;
1.2)设定以下参数:各母线节点的电压和有功功率和无功功率的实时量测值;
步骤(2):根据步骤(1)的参数运用牛顿拉夫逊方法计算系统潮流,并记录至少包括系统各个节点电压的幅值和相位当前潮流断面信息;
形成按以下思路进行:
为了计算和推导的方便,本发明采用了简化的电力系统模型,具体模型如下:研究电力系统同调时有一个基本的假设:发电机模型的复杂程度不改变同调机群之间基本的振荡模态;所以本发明中选取很简单的电力系统模型来设计算法,其中发电机取经典二阶模型表同时忽略阻尼;而负荷模型则采用恒阻抗负荷模型;
发电机取经典二阶模型并忽略阻尼的情况下,电力系统的机电振荡模型如下:
第一个等式中:δ为所有发电机功角组成的列相量,nG为系统中发电机的个数;H为所有发电机惯性时间常数组成的对角矩阵, Pm为所有发电机的机械输入功率组成的列向量, 为常数相量;Pe为所有发电机的电磁输出功率组成的列向量,
第二个等式中:δi为发电机i的功角,Edi为发电机i的暂态电抗后电势,x′di为发电机i的暂态电抗,Ugi为发电机i的机端电压;Ugi为发电机i的机端电压,θUgi为发电机i的机端电压的相角;所以f(Ed,UG,x′d,δ,θUG)是一个矩阵表达式,其形式如下:
每台发电机的暂态电抗后电势及当前功角按照如下思路确定:
对于系统网络的处理,对于发电机和负荷分别进行简化处理;对于发电机,经典二阶模型下的发电机相当于有内阻的电压源内阻是x′di,而电压是Edi,那么经过诺顿等值将这个电压源转换成一个电流源:阻抗仍为x′di,即x′di为发电机i的暂态电抗; 对于负荷节点k进行如下处理:假设当前潮流下,节点的电压模值为VLoad_k,而负荷有功和无功分别为Pk和Qk,那么等效负荷阻抗为 相应得负荷等效导纳为
把负荷节点的等效阻抗以及发电机的暂态电抗加以考虑之后,系统的节点导纳矩阵变为:Y=GY+jBY;其中:GY、BY分别表示为节点导纳矩阵Y的实部和虚部;Y、GY、BY都是n×n维的矩阵;令表示各个节点电压组成的列向量,每个元素为各个节点的电压矢量;代表各个节点的注入电流向量,后面某个列向量上标是一个星号(*)都表示此列向量每个元素都是复数矢量;其中:Ux为各个节点电压实部组成的向量,Uy为各个节点电压虚部组成的向量,Ix为各个节点注入电流实部组成的向量,Iy为个节点注入电流虚部组成的向量;
经过上述简化之后,只有发电机节点具有注入电流,其他的节点没有电流注入;
令系统所有节点的电压模值向量为U=[U1 U2...UN],相角向量为θU={θu1 θu2...θuN}T,N为系统所有节点的个数,那么:
Ux=[U1cosθu1 U2cosθu2...UncosθuN]T,Uy=[U1sinθu1 U2sinθu2...UnsinθuN]T
Ix,Iy都只有发电机节点对应元素非零,其他元素均为零;令:
那么Ix,Iy可以用以下的形式表示:Ix=T×IGx,Iy=T×IGy
其中:IGi表示发电机i节点进行诺顿等值后的注入电流模值;T是一个N×nG维的矩阵;T的每一行对应系统的一个节点(编号为k),如果本节点没有发电机接入,即非发电机节点,那么本行元素全为0,即T(k,:)=[0...0];如果本节点有发电机接入,即为发电机节点,假设此发电机序号为i,那么T(k,i)=1,本行其他元素均为0,即T(k,:)=[0...1...0];
对(2)式两边线性化得:
其中:
ΔU,ΔθU,Δδ分别为U,θU,δ的偏列差量,N和nG分别为系统所有节点的个数以及发电机节点的个数,表达式中的各个矩阵的形式如下:
cosθU=diag{cosθU1 cosθU2...cosθUN};sinθU=diag{sinθU1 sinθU2...sinθUN};
UcosθU=diag{U1cosθU1 U2cosθU2...UNcosθUN};
UsinθU=diag{U1sinθU1 U2sinθU2...UNsinθUN};
令: E是2N×nG维的,那么:
对于式(1)所示的电力系统机电振荡模型在平衡点处进行线性化,可得:
其中:ΔUG,ΔθUG,Δδ为UG,θUG,δ的偏差量, 有:ΔUG=T2×ΔU,ΔθUG=T2×ΔθU,代入(5)式有:
其中: T2=TT
带入(4)式就得到了多机系统线性化模型为:
步骤(4):求取矩阵A的特征值并选择慢模式相关特征值组;
4.1)采用QR算法求取A的全部特征值;
4.2)把A的特征值按照绝对值由小到大的顺序进行排序;
假设排序后的特征值顺序如下: 其中nG为发电机的台数,也正是A的维数;值得注意的是λ1=0;
4.3)选取慢模式相关特征值组,关键是确定选取慢模式相关特征值的个数s;
按照如下的方法选择s:比较所有的λi/λi+1,i=2,3,...,nG-1,令:这样就确定了慢模式相关特征值的个数s;这种选择的依据是电力系统多时间尺度特性原理,λs/λs+1的值越小表示系统的多时间尺度特性越明显,系统的区域性也越明显,系统也越适合解列;经过这一步骤,就确定了慢模式相关特征值组σs={λ1,λ2,..λs};
步骤(5):求取慢模式相关特征值组对于每条线路参数变化灵敏度;
线路就是指电力系统中实际的输电线路,首先对于所有节点进行编号:1,2,...,N,再对线路进行编号:1,2,...,Nl;一条线路主要有如下元素特征:起始节点编号、终了节点编号、线路电阻r、线路电抗x,线路电导g和线路电纳b,值得注意的是在选取线路起止节点时,哪一端点作为起始节点是无所谓的,所以在确定起始节点时随便指定其中一个端点作为起始节点即可。
某条线路的编号是L,其起始节点编号为S,终了节点编号为J,线路的电阻和电抗分别为r和x,电导和电纳分别为g和b,求取慢模式相关特征值组对于线路L参数变化的灵敏度,可按以下两步进行:
特征值λi对于线路L参数变化的灵敏度可按照以下思路进行:
λi为A的一个特征值,λi对应的左右特征向量分别为φi(nG×1)和ψi(1×nG);
那么λi对于线路L的导纳b变化的灵敏度表达式如下:
其中: M=[...1...-1...]T,其中1的位置是S,而-1的位置是J;
λi对于线路L的导纳g变化的灵敏度表达式如下:
其中:
5.2)确定慢模式相关特征值组σs={λ1,λ2,..λs}对于线路L参数变化的灵敏度,按下式确定:
经过这一步骤就能够得到所有线路的慢模式相关特征值组对于系统所有线路的参数的灵敏度信息;
步骤(6):根据步骤5得到的慢模式相关特征值组对于线路参数变化的灵敏度信息来确定弱连接;
6.1)把慢模式相关特征值组对于系统所有的Nl条线路的灵敏度按照绝对值从大到小的顺序进行排序,排序后的顺序为
6.2)选取弱连接线路,基本思路就是选取灵敏度比较大的一组线路作为系统弱连接线路;至于弱连接线路的条数wl由用户自己按照需要设定,或自动生成,所述自动生成wl的方法:对于一个多时间尺度特性比较明显的系统,弱连接线路与非弱连接线路的灵敏度数值差异会比较大;这种方法首先对于系统所有线路的灵敏度进行从大到小的排序,然后按照上述原理来选择弱连接线路的条数wl,选取wl的公式表达形式如下:
其中,符号abs表示取一个实数的绝对值。
本发明的有益效果是该弱连接识别方法能够直接处理线路,通过计算系统的一组慢模式对应的特征值对各条线路参数变化灵敏度,选取灵敏度数值比较大的线路作为系统的弱连接。
附图说明
图1算法流程图。
图2发电机诺顿等值示意图
图3网络处理效果图
图4Kundurll节点测试系统拓扑图及相应的弱连接示意图
具体实施方式:
本发明提出了一种基于慢模式相关特征值对于线路参数灵敏度的弱连接识别方法。本节利用Kundurll节点系统来测试本发明的选取系统弱连接的方法。Kundurll节点系统包括11条母线,4台发电机,12条线路(包含变压器支路),系统拓扑图见图4。
参照附图1所示的流程步骤识别此系统的弱连接:
步骤(1):开始,计算机初始化;
1.1)输入以下参数
电力系统的网络参数,包括:输电线路的串联电阻、串联电抗、并联电导和并联电纳;变压器的变比和阻抗;并联在输电线路上的电容器和电抗器的阻抗;发电机的动态参数包括发电机的惯性时间常数;发电机的暂态电抗;
1.2)设定以下参数:各母线节点的电压和有功功率和无功功率的实时量测值;
步骤(2):根据步骤(1)的参数运用牛顿拉夫逊方法计算系统潮流,并记录至少包括系统各个节点电压的幅值和相位当前潮流断面信息;
步骤(3):本发明中采用系统简化的模型来进行计算,首先需要对网络进行处理,对于发电机和负荷分别进行简化处理;对于发电机,经典二阶模型下的发电机相当于有内阻的电压源(内阻是x′di,而电压是Edi),那么经过诺顿等值将这个电压源转换成一个电流源:阻抗仍为x′di, 诺顿等值示意图参见图2;对于负荷节点k进行如下处理:假设当前潮流下,节点的电压模值为VLoad_k,而负荷有功和无功分别为Pk和Qk,那么等效负荷阻抗为 相应得负荷等效导纳为 此网络处理效果的示意图参见图3;
经过以上的网络处理后,根据式(7) 计算矩阵A=BT′2D-1ZE+C其中各个矩阵的含义参见发明内容中的详细推导过程。由此得到:
步骤(4):求取矩阵A的特征值并选择慢模式相关特征值组;
4.1)采用QR算法求取A的全部特征值及对应的左右特征向量并存储;
4.2)把A的特征值按照绝对值由小到大的顺序进行排序,所得结果如表1;
表1系统特征值
4.3)选取慢模式相关特征值组,关键是确定选取慢模式相关特征值的个数r;
按照如下的方法选择s:比较所有的λi/λi+1,i=2,3,...,nG-1,令:这样就确定了慢模式相关特征值的个数s;这种选择的依据是电力系统多时间尺度特性原理,λs/λs+1的值越小表示系统的多时间尺度特性越明显,系统的区域性也越明显,系统也越适合解列;经过这一步骤,就确定了慢模式相关特征值组σs={λ1,λ2,..λs};
对于本实施例:
{λ2/λ3,λ3/λ4}={-0.0265/-0.1386 -0.1386/-0.1461}={0.1912 0.9487}
即:s=2,所以σs={λ1,λ2}={0 -0.0265}
步骤(5):求取慢模式相关特征值组对于每条线路参数变化灵敏度;
5.1)根据式(8)(9)(10)求取慢模式相关特征值组σs={λ1,λ2,..λs}对于各条线路参数变化的灵敏度;
5.2)根据式(11)确定慢模式相关特征值组σs={λ1,λ2,..λs}对于各条线路线路参数变化的灵敏度;
本实施例中,模式特征值组中各个特征值以及此特征值组整体对于各条线路参数变化灵敏度参见表2:
表2慢模式相关特征值对于线路参数变化的灵敏度
步骤(6):根据步骤(5)得到的慢模式相关特征值组对于线路参数变化的灵敏度信息来确定弱连接;
6.1)把慢模式相关特征值组对于系统所有的Nl(本实施例中Nl=12)条线路的灵敏度按照绝对值从大到小的顺序进行排序;排序后的结果参见表3:
表3慢模式相关特征值组对线路参数灵敏度排序后的结果
6.2)按照式(12)选取弱连接线路的条数wl:
即:wl=4,那么就得到此实施例中的弱连接线路如表4:
表4系统的弱连接线路
此实施例中系统的弱连接线路在图4中用虚线标识,参见图4。
Claims (1)
1.一种慢模式相关特征值对于线路参数灵敏度的弱连接识别方法,其特征在于,具体算法步骤如下:
步骤(1):计算机初始化;
1.1)输入以下参数:
电力系统的网络参数,包括:输电线路的串联电阻、串联电抗、并联电导和并联电纳;变压器的变比和阻抗;并联在输电线路上的电容器和电抗器的阻抗;发电机的动态参数包括发电机的惯性时间常数;发电机的暂态电抗;
1.2)设定以下参数:各母线节点的电压和有功功率和无功功率的实时量测值;
步骤(2):根据步骤(1)的参数运用牛顿拉夫逊方法计算系统潮流,并记录至少包括系统各个节点电压的幅值和相位当前潮流断面信息;
形成按以下思路进行:为了计算和推导的方便,采用了简化的电力系统模型,具体模型如下:研究电力系统同调时有一个基本的假设:发电机模型的复杂程度不改变同调机群之间基本的振荡模态;所以选取很简单的电力系统模型来设计算法,其中发电机取经典二阶模型表同时忽略阻尼;而负荷模型则采用恒阻抗负荷模型;
发电机取经典二阶模型并忽略阻尼的情况下,电力系统的机电振荡模型如下:
第一个等式中:δ为所有发电机功角组成的列相量,nG为系统中发电机的个数;H为所有发电机惯性时间常数组成的对角矩阵, Pm为所有发电机的机械输入功率组成的列向量, 为常数相量;Pe为所有发电机的电磁输出功率组成的列向量,
第二个等式中:δi为发电机i的功角,Edi为发电机i的暂态电抗后电势,x′di为发电机i的暂态电抗,Ugi为发电机i的机端电压;θUgi为发电机i的机端电压的相角;所以f(Ed,UG,x′d,δ,θUG)是一个矩阵表达式,其形式如下:
每台发电机的暂态电抗后电势及当前功角按照如下思路确定:
对于系统网络的处理,对于发电机和负荷分别进行简化处理;对于发电机,经典二阶模型下的发电机相当于有内阻的电压源内阻是x′di,而电压是Edi,那么经过诺顿等值将这个电压源转换成一个电流源:阻抗仍为x′di,即x′di为发电机i的暂态电抗; 对于负荷节点k进行如下处理:假设当前潮流下,节点的电压模值为VLoad_k,而负荷有功和无功分别为Pk和Qk,那么等效负荷阻抗为 相应得负荷等效导纳为
把负荷节点的等效阻抗以及发电机的暂态电抗加以考虑之后,系统的节点导纳矩阵变为:Y=GY+jBY;其中:GY、BY分别表示为节点导纳矩阵Y的实部和虚部;Y、GY、BY都是n×n维的矩阵;令表示各个节点电压组成的列向量,每个元素为各个节点的电压矢量;代表各个节点的注入电流向量,后面某个列向量上标是一个星号(*)都表示此列向量每个元素都是复数矢量;其中:Ux为各个节点电压实部组成的向量,Uy为各个节点电压虚部组成的向量,Ix为各个节点注入电流实部组成的向量,Iy为个节点注入电流虚部组成的向量;
经过上述简化之后,只有发电机节点具有注入电流,其他的节点没有电流注入;
令系统所有节点的电压模值向量为U=[U1 U2 … UN],相角向量为θU={θu1 θu2 … θuN}T,N为系统所有节点的个数,那么:
Ux=[U1cosθu1 U2cosθu2 … UncosθuN]T,Uy=[U1sinθu1 U2sinθu2 … UnsinθuN]T,
Ix,Iy都只有发电机节点对应元素非零,其他元素均为零;令:
那么Ix,Iy可以用以下的形式表示:Ix=T×IGx,Iy=T×IGy
其中:T是一个N×nG维的矩阵;T的每一行对应系统的一个节点(编号为k),如果本节点没有发电机接入,即非发电机节点,那么本行元素全为0,即T(k,:)=[0 … 0];如果本节点有发电机接入,即为发电机节点,假设此发电机序号为i,那么T(k,i)=1,本行其他元素均为0,即T(k,:)=[0 … 1… 0];
对(2)式两边线性化得:
其中:
ΔU,ΔθU,Δδ分别为U,θU,δ的偏列差量,N和nG分别为系统所有节点的个 数以及发电机节点的个数,表达式中的各个矩阵的形式如下:
cosθU=diag{cosθU1 cosθu2 … cosθUN};sinθU=diag{sinθU1 sinθU2 … sinθUN};
UcosθU=diag{U1cosθU1 U2cosθU2 … UNcosθUN};
UsinθU=diag{U1sinθU1 U2sinθU2 … UNsinθUN};
对于式(1)所示的电力系统机电振荡模型在平衡点处进行线性化,可得:
步骤(4):求取矩阵A的特征值并选择慢模式相关特征值组;
4.1)采用QR算法求取A的全部特征值;
4.2)把A的特征值按照绝对值由小到大的顺序进行排序;
4.3)选取慢模式相关特征值组,关键是确定选取慢模式相关特征值的个数s;
按照如下的方法选择s:比较所有的λi/λi+1,i=2,3,...,nG-1,令: 这样就确定了慢模式相关特征值的个数s;这种选择的依据是电力系统多时间尺度特性原理,λs/λs+1的值越小表示系统的多时间尺度特性越明显,系统的区域性也越明显,系统也越适合解列;经过这一步骤,就确定了慢模式相关特征值组σs={λ1,λ2,..λs};
步骤(5):求取慢模式相关特征值组对于每条线路参数变化灵敏度;
线路就是指电力系统中实际的输电线路,首先对于所有节点进行编号:1,2,...,N,再对线路进行编号:1,2,...,Nl;一条线路主要有如下元素特征:起始节点编号、终了节点编号、线路电阻r、线路电抗x,线路电导g和线路电纳b,值得注意的是在选取线路起止节点时,哪一端点作为起始节点是无所谓的,所以在确定起始节点时随便指定其中一个端点作为起始节点即可;
某条线路的编号是L,其起始节点编号为S,终了节点编号为J,线路的电阻 和电抗分别为r和x,电导和电纳分别为g和b,求取慢模式相关特征值组对于线路L参数变化的灵敏度,可按以下两步进行:
特征值λi对于线路L参数变化的灵敏度可按照以下思路进行:
λi为A的一个特征值,λi对应的左右特征向量分别为Φi(nG×1)和ψi(1×nG);
那么λi对于线路L的导纳b变化的灵敏度表达式如下:
λi对于线路L的导纳g变化的灵敏度表达式如下:
5.2)确定慢模式相关特征值组σs={λ1,λ2,..λs}对于线路L参数变化的灵敏 度,按下式确定:
即为慢模式相关特征值组σs={λ1,λ2,..λs}对于线路L参数变化的灵敏度,用这个量来衡量线路参数变化对于慢模式组的影响;
经过这一步骤就能够得到所有线路的慢模式相关特征值组对于系统所有线路的参数的灵敏度信息;
步骤(6):根据步骤(5)得到的慢模式相关特征值组对于线路参数变化的灵敏度信息来确定弱连接;
6.2)选取弱连接线路,基本思路就是选取灵敏度比较大的一组线路作为系统弱连接线路;至于弱连接线路的条数wl由用户自己按照需要设定,或自动生成,所述自动生成wl的方法:对于一个多时间尺度特性比较明显的系统,弱连接线路与非弱连接线路的灵敏度数值差异会比较大;这种方法首先对于系统所有线路的灵敏度进行从大到小的排序,然后按照上述原理来选择弱连接线路的条数wl,选取wl的公式表达形式如下:
其中,符号abs表示取一个实数的绝对值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102383341A CN101877482B (zh) | 2009-11-25 | 2009-11-25 | 基于慢模式特征值对于线路参数灵敏度的弱连接识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102383341A CN101877482B (zh) | 2009-11-25 | 2009-11-25 | 基于慢模式特征值对于线路参数灵敏度的弱连接识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101877482A CN101877482A (zh) | 2010-11-03 |
CN101877482B true CN101877482B (zh) | 2013-05-22 |
Family
ID=43019970
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102383341A Active CN101877482B (zh) | 2009-11-25 | 2009-11-25 | 基于慢模式特征值对于线路参数灵敏度的弱连接识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101877482B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103578040A (zh) * | 2013-09-05 | 2014-02-12 | 国家电网公司 | 薄弱环节跟踪系统 |
CN108964010A (zh) * | 2017-05-19 | 2018-12-07 | 国网安徽省电力公司 | 一种确定电网设备对电网安全指标的灵敏度的方法和系统 |
CN107196337B (zh) * | 2017-06-29 | 2019-08-20 | 广东电网有限责任公司电力调度控制中心 | 一种孤岛静态工作点校验方法及系统 |
CN112034251B (zh) * | 2020-08-28 | 2021-11-23 | 武汉大学 | 评价诺顿型间谐波源接入点激励并联谐振能力的方法 |
CN113131473B (zh) * | 2021-04-22 | 2023-04-07 | 云南电网有限责任公司电力科学研究院 | 一种基于参数灵敏度的输电系统可用传输容量计算方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1305196C (zh) * | 2003-05-10 | 2007-03-14 | 清华大学 | 灾变事故下避免电力系统崩溃的解列方法 |
CN101136549A (zh) * | 2007-10-12 | 2008-03-05 | 清华大学 | 电力系统解列决策空间筛选方法 |
CN101521407A (zh) * | 2008-11-21 | 2009-09-02 | 四川省电力公司 | 一种三段式失步自适应解列的方法 |
-
2009
- 2009-11-25 CN CN2009102383341A patent/CN101877482B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1305196C (zh) * | 2003-05-10 | 2007-03-14 | 清华大学 | 灾变事故下避免电力系统崩溃的解列方法 |
CN101136549A (zh) * | 2007-10-12 | 2008-03-05 | 清华大学 | 电力系统解列决策空间筛选方法 |
CN101521407A (zh) * | 2008-11-21 | 2009-09-02 | 四川省电力公司 | 一种三段式失步自适应解列的方法 |
Non-Patent Citations (3)
Title |
---|
JP特开平10-258170A 1998.09.29 |
乔颖等.《大电网解列决策空间筛选及快速搜索方法》.《中国电机工程学报》.2008,第28卷(第22期),23-28. * |
乔颖等.《大规模电网解列控制可行性判断》.《中国电机工程学报》.2008,第28卷(第25期),50-55. * |
Also Published As
Publication number | Publication date |
---|---|
CN101877482A (zh) | 2010-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101877482B (zh) | 基于慢模式特征值对于线路参数灵敏度的弱连接识别方法 | |
CN101505061B (zh) | 一种基于时域仿真的戴维南等值参数跟踪的计算方法 | |
Wu et al. | The SOC estimation of power Li-Ion battery based on ANFIS model | |
CN109800520B (zh) | 一种基于神经网络的电动汽车充电站谐波建模方法 | |
CN105938578A (zh) | 一种基于聚类分析的大规模光伏电站等值建模方法 | |
CN107505392A (zh) | 基于纹理表面接触加速度触觉信息的材质分析方法及装置 | |
CN105846448A (zh) | 一种基于随机矩阵理论的配电网无功补偿容量确定方法 | |
CN105391059A (zh) | 一种基于电流量测变换的分布式发电系统状态估计方法 | |
CN107204617A (zh) | 基于线性规划的直角坐标形式的区间潮流计算方法 | |
CN106208099A (zh) | 一种基于二层规划的电力系统无功优化方法及其应用 | |
CN106921161A (zh) | 电力系统薄弱点识别方法及薄弱点电压分布特性分析方法 | |
CN106372724A (zh) | 一种人工神经网络算法 | |
CN107036808A (zh) | 基于支持向量机概率估计的风电机组齿轮箱复合故障诊断方法 | |
CN105956252A (zh) | 基于生成型深信度网络的超短期风速多尺度预报建模方法 | |
CN102545205A (zh) | 基于自动微分技术的分布式电力系统状态估计方法 | |
CN104156504B (zh) | 一种发电机励磁系统参数可辨识性判定方法 | |
De Tuglie et al. | A coherency recognition based on structural decomposition procedure | |
CN109579967A (zh) | 智能动态称重方法与系统 | |
Joo et al. | Enhancement of coherency identification techniques for power system dynamic equivalents | |
CN107230982A (zh) | 一种微电网线性化潮流计算方法 | |
Chao et al. | A maximum power point tracker with automatic step size tuning scheme for photovoltaic systems | |
CN113010844B (zh) | 一种基于子空间动态模式分解的参与因子计算方法 | |
Xiao et al. | A novel impedance measurement technique for power electronic systems | |
CN104600700B (zh) | 基于发电机诺顿等值和广义逆的网损微增率计算方法 | |
CN102707161B (zh) | 一种基于动态转移阻抗的短路容量在线监测方法及其装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |