CN103700036B - 一种适于电力系统多时间尺度的暂态稳定性投影积分方法 - Google Patents

一种适于电力系统多时间尺度的暂态稳定性投影积分方法 Download PDF

Info

Publication number
CN103700036B
CN103700036B CN201310730821.6A CN201310730821A CN103700036B CN 103700036 B CN103700036 B CN 103700036B CN 201310730821 A CN201310730821 A CN 201310730821A CN 103700036 B CN103700036 B CN 103700036B
Authority
CN
China
Prior art keywords
time
power
simulation
length
integral
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.)
Expired - Fee Related
Application number
CN201310730821.6A
Other languages
English (en)
Other versions
CN103700036A (zh
Inventor
王成山
原凯
李鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
State Grid Corp of China SGCC
Shanghai Municipal Electric Power Co
Original Assignee
Tianjin University
State Grid Corp of China SGCC
Shanghai Municipal Electric Power Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University, State Grid Corp of China SGCC, Shanghai Municipal Electric Power Co filed Critical Tianjin University
Priority to CN201310730821.6A priority Critical patent/CN103700036B/zh
Publication of CN103700036A publication Critical patent/CN103700036A/zh
Application granted granted Critical
Publication of CN103700036B publication Critical patent/CN103700036B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明的一种适于电力系统多时间尺度的暂态稳定性投影积分方法,采用显式交替求解算法对电力系统数学模型进行求解,关于电力系统数学模型中微分方程的求解则采用投影积分算法:首先进行若干小步长的积分运算,对应系统的快动态过程;而后根据小步长积分计算结果,进行一个投影步长的计算,与系统的慢动态过程对应。其中,小步长积分计算过程称为内部积分器,采用数值稳定性较好的显式四阶龙格‑库塔法以提高算法的稳定性;大步长积分计算过程称为投影积分器。本发明不仅适用于传统电力系统暂态稳定性仿真,也适用于智能配电网动态仿真问题,在满足仿真数值稳定性和数值精度要求的同时,提升了仿真计算速度,为高效、可靠的电力系统仿真程序的开发奠定了良好的基础。

Description

一种适于电力系统多时间尺度的暂态稳定性投影积分方法
技术领域
[0001] 本发明涉及一种电力系统暂态稳定性仿真方法。特别是涉及一种适于电力系统多 时间尺度的暂态稳定性投影积分方法。
背景技术
[0002] 电力系统是由发电、变电、输电、配电、用电等设备和相应的辅助系统组成的复杂 大系统。在电力系统层面,相关问题的分析与研究往往无法直接在实际系统上进行试验,因 此必须采用有效的数字仿真作为重要的研究手段。
[0003] 在电力系统时域仿真的研究中,针对系统动态过程的不同时间尺度分别发展出电 磁暂态仿真、机电暂态仿真和中长期动态仿真三种电力系统数字仿真方法,三者从元件数 学模型到仿真计算方法都具有不同的特征。电力系统电磁暂态仿真侧重于系统中电场与磁 场相互影响产生的电压电流的变化过程;机电暂态仿真主要研究电力系统在大扰动下(如 故障、切机、切负荷、重合闸操作等情况)的动态行为和保持同步稳定运行的能力,即暂态稳 定性,所关注的时间范围通常为几秒至几十秒,因而也称为暂态稳定性仿真;中长期动态过 程仿真是电力系统受到扰动后较长过程的动态仿真,即通常的电力系统长过程动态稳定计 算。
[0004] 电力系统暂态稳定性仿真除关注传统电力系统的暂态稳定运行能力外,还侧重于 分析含各种分布式电源及储能装置的智能配电系统运行时其工频电气量在系统扰动下(开 关操作、故障、分布式电源及负荷波动等)的动态响应特性,此时也可称为智能配电系统动 态仿真。智能配电系统接入了种类繁多的分布式电源和大量的电力电子装置,包括旋转电 机和各种静态直流型分布式电源,具有明显的多时间尺度特征。电力系统暂态稳定性仿真 的多时间尺度问题在数学上可以归结为刚性问题,其对所采用的数值算法的精度和稳定性 要求更高。
[0005] 电力系统暂态稳定性仿真本质上可归结为对动力学系统时域响应的求取,分为数 学建模和模型求解两部分。暂态稳定性仿真首先根据元件间拓扑关系将电力系统各元件模 型构成全系统模型,形成一组联立的微分-代数方程组,然后以稳态工况或潮流解为初值, 求解扰动下的数值解,即逐步求得系统状态量和代数量随时间的变化曲线。
[0006] 电力系统数字仿真建模的过程则是根据系统仿真关注的时间尺度范围,由物理原 型抽象出数学模型。电力系统暂态稳定性仿真中的数学模型包括两部分:描述设备动态特 征的微分方程和描述设备之间电气联系的代数方程。
[0007] 动态设备之间的电气连接关系在运行过程中可能改变,如负荷的投切、机组的启 停、线路开断和重合闸等操作,若计及继电保护装置,还应包含大量连续和(或)离散的逻辑 时变参数。一般可将电力系统数学模型通过一个高维非线性且连续自治的微分-代数方程 组来描述,如式(1)所示。
[O^g^x.y) (1)
[0009] 式中,/(/:心xM~4m)为微分方程,发bM~)为代数方程,eR~) 为系统状态变量,代表电机转子转速、电力电子器件控制系统和负荷动态参数等,:Kj e R"' > 为代数变量,表征母线电压幅值与相角。数学模型的求解一般通过特定的数值算法和相应 的仿真程序来实现。因此,电力系统暂态稳定性仿真在数学上可以归结为求解一个微分-代 数方程组的初值问题。
[0010] 暂态稳定性仿真算法按照对于式(1)中微分方程和代数方程解算形式的不同可以 分为交替求解法和联立求解法两大类。交替求解法首先采用特定的数值积分算法,根据初 始化计算结果求解微分方程,得到本时步状态变量的值,然后将其代入到代数方程中求解, 得到该时步代数变量的值,最后再将代数变量代入微分方程进行下一时步状态变量求解, 以此类推实现微分-代数方程组的交替求解;联立求解法则是将微分方程差分化之后,和代 数方程联立成一个完整的代数方程组,同时求解状态变量和代数变量。
[0011] 对于式(1)中的微分方程,除少数可得到解析解以外,大多数只能采用数值解法进 行求解。微分方程数值解求解有两种基本方法:一种是把近似解表示成有限个独立函数之 和;另一种是差分法,该方法常用于电力系统数字仿真。差分法是寻求在一系列离散点上X 的近似值的一类方法的统称,两相邻离散点之间的距离为步长,记为h。根据状态变量 x(t) 在tn+i时刻的数值解x(tn+i)能否由tn时刻的值直接求出,差分法又可分为单步法(one step method)和线性多步法(linear multistep method),其中,单步法在电力系统暂态稳定性 仿真中应用广泛。
[0012]根据求解过程的不同,单步法可分为显式积分方法和隐式积分方法,显式积分方 法可根据当前时刻状态变量直接计算下一时刻状态变量,而隐式积分方法则需要对含有当 前时刻和下一时刻状态变量的方程进行求解才能求得下一时刻状态变量。常见的显式积分 方法包括欧拉法、改进欧拉法和龙格_库塔法,而隐式积分方法主要有后向欧拉法和隐式梯 形法。电力系统动态过程时间尺度差异较大,具有较明显的刚性特征。显式积分方法每一时 步内的运算量较小,但由于其数值稳定性较差,即在某一时步计算中产生的误差在以后逐 步积分过程中不断累积而可能导致计算结果数值不收敛,因此针对刚性问题求解往往只能 采取较小的仿真步长,仿真速度受到极大限制。隐式积分方法虽然数值稳定性较好,可以在 刚性问题的求解过程中保证数值稳定性,但其每一时步都需要迭代求解方程组,相比显式 积分算法其计算与编程工作复杂,极大地限制了其在大规模电力系统暂态稳定性仿真中的 应用。
[0013] 可见,开发一种计算简单、数值稳定性好、适用于刚性问题的电力系统暂态稳定性 仿真方法十分重要。
发明内容
[0014] 本发明所要解决的技术问题是,提供一种不仅适用于传统电力系统暂态稳定性仿 真,也适用于智能配电网动态仿真问题的适于电力系统多时间尺度的暂态稳定性投影积分 方法。
[0015] 本发明所采用的技术方案是:一种适于电力系统多时间尺度的暂态稳定性投影积 分方法,包括如下步骤:
[0016] 1)输入电力系统潮流数据、动态元件数据以及故障及操作设置信息,设置仿真参 数;
[0017] 2)根据系统算例参数和动态数据形成全系统节点导纳矩阵,其中包括系统节点导 纳矩阵和并入动态元件的节点导纳矩阵,并重新进行潮流校核;
[0018] 3)对系统算例的动态元件进行仿真初始化计算;
[0019] 4)设置仿真计算时间t=0;
[0020] 5)设置当前投影积分算法内部积分器的积分步数s=l;
[0021] 6)采用投影积分算法内部积分器对电力系统模型积分一个步长,步长为h,然后设 置仿真计算时间t=t+h,内部积分器的积分步数s=s+l,内部积分器采用显式交替求解方法 对微分_代数方程组进行求解,对其中的微分方程求解则采用显式四阶龙格-库塔法;
[0022] 7)根据步骤1)输入的故障及操作设置信息判断系统是否发生故障或操作,若发 生,则返回步骤5),否则进入下一步骤;
[0023] 8)判断仿真计算时间t是否大于设置时间T,若计算时间t大于设置时间T,则仿真 结束,否则进入下一步骤;
[0024] 9)判断投影积分算法内部积分步数s是否大于设置步数k+1,若不大于,则返回步 骤6),否则进入下一步骤;
[0025] 10)根据步骤1)输入的故障及操作设置信息判断未来Mh时间内是否存在故障或操 作,若存在,则返回步骤5),否则进入下一步骤;
[0026] 11)判断未来Mh时间内仿真计算时间t是否会超过设置时间T,若超过,则返回步骤 5),否则进入下一步骤;
[0027] 12 )根据内部积分器得到的x (tn+k)和x (tn+k+i)的值进行一步外部积分步,步长为内 部积分器积分步长的]^倍,利用式1(1:11+1 {+1+")=(]\/[+1)1(1:11+1{+1)-]\^(1:11+1〇得到1(1:11+1 {+1+"),然后 求解代数方程g ( x ( tn+k+1+M),y ( tn+k+1+M) ) =0得到y (tn+k+1+M),设置t=t+Mh ;
[0028] 13)判断此时仿真计算时间t是否达到设置的仿真时间T,若达到,则仿真结束,否 则返回步骤5),依此反复进行直至仿真结束。
[0029] 步骤1)中所述的仿真参数,包括仿真时间T、仿真步长h和投影积分算法参数,其 中,所述的投影积分算法参数包括内部积分器的积分步数k和外部积分器步长相对内部积 分器步长的倍数M。
[0030] 所述的内部积分器算法是:
[0031] (1)采用显式四阶龙格-库塔法积分k步,步长为h,时间从trjlju+k,由X(tn)得到X (t n+k),递推公式如下;
Figure CN103700036BD00061
[0037] k4=hf (x (tn) +k3, tn+h);
[0038] (2)再次采用显式四阶龙格-库塔法积分一步,步长仍为h,得到x(tn+k+1)。
[0039]所述的外部积分器算法是:
[0040]根据内部积分器得到的x(tn+k)和x(tn+k+1)的值,利用下式进行一步外部积分步,步 长为内部积分步长的M倍,时间从tn+k+l到tn+k+1+M,得到X( tn+k+1+M):
[0041 ] X(tn+k+l+M) = (M+l)X(tn+k+l)-Mx(tn+k)。
[0042]本发明的一种适于电力系统多时间尺度的暂态稳定性投影积分方法,考虑了电力 系统的多时间尺度特征,采用显式交替求解方法对描述电力系统模型的微分-代数方程组 进行交替求解,对其中的微分方程利用投影积分算法进行差分求解。本发明的方法不仅适 用于传统电力系统暂态稳定性仿真,也适用于智能配电网动态仿真问题,具有一般性,在满 足仿真数值稳定性和数值精度要求的同时,提升了仿真计算速度,为高效、可靠的电力系统 仿真程序的开发奠定了良好的基础。本发明的方法计算简单、数值稳定性好、适用于具有刚 性问题的电力系统暂态稳定性仿真。
附图说明
[0043]图1是本发明方法的整体流程图;
[0044]图2是低压配电系统算例结构图;
[0045] 图中1:燃料电池;2:第二蓄电池;3:第二光伏电池;4:第一光光伏电池;
[0046] 5:第一蓄电池;Ml:中压母线;L1~L19:低压母线;Loadl~Load7:负荷;
[0047] 图3是L16母线电压仿真结果及局部放大图;
[0048] 图4是L17母线电压仿真结果及局部放大图;
[0049] 图5是L17母线频率仿真结果及局部放大图;
[0050] 图6是L17母线电压投影积分算法仿真结果较DIgSILENT相对误差。
具体实施方式
[0051] 下面结合实施例和附图对本发明的一种适于电力系统多时间尺度的暂态稳定性 投影积分方法做出详细说明。
[0052]本发明的一种适于电力系统多时间尺度的暂态稳定性投影积分方法,属于显式积 分方法,考虑了电力系统的多时间尺度特征,采用显式交替求解方法对电力系统模型进行 交替求解,对其中的微分方程利用投影积分算法进行差分求解。不仅适用于传统电力系统 暂态稳定性仿真,也同样适用于含各种分布式电源及储能装置的智能配电系统动态仿真。 智能配电系统接入了种类繁多的分布式电源和大量的电力电子装置,包括旋转电机和各种 静态直流型分布式电源,具有明显的多时间尺度特征。因此,电力系统暂态稳定性仿真的多 时间尺度问题在数学上可以归结为刚性问题,其对所采用的数值算法的精度和稳定性要求 更高。本发明的方法具有一般性,不仅适用于传统电力系统暂态稳定性仿真,也适用于含分 布式电源及储能的智能配电网动态仿真。本发明在满足仿真数值稳定性和数值精度要求的 同时,提升了仿真计算速度,为高效、可靠的电力系统仿真程序的开发奠定了良好的基础。 [0053]本发明采用显式交替求解算法对电力系统数学模型进行求解,关于电力系统数学 模型中微分方程的求解则采用投影积分算法:首先进行若干小步长的积分运算,对应系统 的快动态过程;而后根据小步长积分计算结果,进行一个投影步长的计算,与系统的慢动态 过程对应。其中,小步长积分计算过程称为内部积分器,采用数值稳定性较好的显式四阶龙 格-库塔法(explicit four-order Runge-Kutta method)以提高算法的稳定性;大步长积 分计算过程称为外部积分器,通过大步长投影计算以提升算法效率。
[0054]如图1所示,本发明的一种适于电力系统多时间尺度的暂态稳定性投影积分方法, 包括如下步骤:
[0055] 1)输入电力系统潮流数据、动态元件数据以及故障及操作设置信息,设置仿真参 数;
[0056] 所述的仿真参数,包括仿真时间T、仿真步长h和投影积分算法参数,其中,所述的 投影积分算法参数包括内部积分器的积分步数k和外部积分器步长相对内部积分器步长的 倍数M 〇
[0057]所述的内部积分器算法是:
[0058] (1)米用显式四阶龙格-库塔法(explicit four-order Runge-Kutta method)积 分k步,步长为h,时间从tn到tn+k,由x(tn)得到x(t n+k),递推公式如下;
Figure CN103700036BD00081
[0060] 其中
Figure CN103700036BD00082
[0064] k4=hf (x (tn) +k3, tn+h)
[0065] (2)再次采用显式四阶龙格-库塔法积分一步,步长仍为h,得到x(tn+k+1)。
[0066]所述的外部积分器算法是:
[0067] 根据内部积分器得到的x(tn+k)和x(tn+k+1)的值,利用下式进行一步外部积分步,步 长为内部积分步长的M倍,时间从tn+k+l到tn+k+1+M,得到X( tn+k+1+M):
[0068] X ( tn+k+1+M) = (M+1 ) X ( tn+k+1) ~Mx ( tn+k)
[0069] 2)根据系统算例参数和动态数据形成全系统节点导纳矩阵,其中包括系统节点导 纳矩阵和并入动态元件的节点导纳矩阵,并重新进行潮流校核;
[0070] 3)对系统算例的动态元件进行仿真初始化计算;
[0071] 4)设置仿真计算时间t=0;
[0072] 5)设置当前投影积分算法内部积分器的积分步数s=l;
[0073] 6)采用投影积分算法内部积分器对电力系统模型积分一个步长,步长为h,然后设 置仿真计算时间t=t+h,内部积分器的积分步数s=s+l,内部积分器采用显式交替求解方法 对微分_代数方程组进行求解,对其中的微分方程求解则采用显式四阶龙格-库塔法;
[0074]所述的电力系统模型是下式所示的一个高维非线性且连续自治的微分-代数方程 组:
[0076] 式中,/(/ : f h^R">)为微分方程,用于描述电力系统设备动态特征; f(f M为代数方程,用于描述电力系统设备之间的电气联系;坤尤識~)为系 统状态变量,代表同步电机转子转速、电力电子器件控制系统和负荷动态参数等;e K~) 为代数变量,表征母线电压幅值与相角。
[0077]对电力系统模型显式差分化后得到下式:
Figure CN103700036BD00091
[0079]所述的显式交替求解方法对显式差分化后的电力系统模型进行求解包括如下步 骤:
[0080] (1)根据1时刻电力系统的状态变量奴1)和代数变量7(1),计算向量1^=1^(1 (tn) ,y(tn));
Figure CN103700036BD00092
[0081] (2)计算向量 ,然后求解代数方程g(xi,yi)=0得出yi,计算向量k2=hf (xi,yi);
[0082] (3)计算向量, :,然后求解代数方程g(X2,y2)=〇得出y2,计算向量k3=hf (X2,y2);
[0083] (4)计算向量_ ,然后求解代数方程g(X3,y3)=〇得出y3,计算向量k4=hf (X3,y3);
[0084] (5)计算tn+1时刻
Figure CN103700036BD00093
,然后求解代数方程g(x (tn+i),y(tn+i))=0得出y(tn+i)。
[0085] 7)根据步骤1)输入的故障及操作设置信息判断系统是否发生故障或操作,若发 生,则返回步骤5),否则进入下一步骤;
[0086] 8)判断仿真计算时间t是否大于设置时间T,若计算时间t大于设置时间T,则仿真 结束,否则进入下一步骤;
[0087] 9)判断投影积分算法内部积分步数s是否大于设置步数k+1,若不大于,则返回步 骤6),否则进入下一步骤;
[0088] 10)根据步骤1)输入的故障及操作设置信息判断未来Mh时间内是否存在故障或操 作,若存在,则返回步骤5),否则进入下一步骤;
[0089] 11)判断未来Mh时间内仿真计算时间t是否会超过设置时间T,若超过,则返回步骤 5),否则进入下一步骤;
[0090 ] 12 )根据内部积分器得到的x (tn+k)和x (tn+k+1)的值进行一步外部积分步,步长为内 部积分器积分步长的]^倍,利用式1(1:11+1{+1+")=(]\/[+1)1(1:11+1 {+1)-]\^(1:11+1〇得到1(1:11+1{+1+"),然后 求解代数方程g ( x ( tn+k+l+M),y ( tn+k+l+M) ) =0得到y (tn+k+l+M),设置t=t+Mh ;
[0091] 13)判断此时仿真计算时间t是否达到设置的仿真时间T,若达到,则仿真结束,否 则返回步骤5),依此反复进行直至仿真结束。
[0092]下面给出具体实例:
[0093] 数字仿真和电网计算程序(DIgSILENT PowerFactory)是德国DIgSLENTGmbH公司 开发的一款商业电力系统仿真软件。本实例以C++编程语言环境为基础,实现了面向电力系 统暂态稳定性仿真的投影积分仿真方法,通过含分布式电源的低压配电系统算例(附图2) 对本发明的方法进行测试验证,并与商业软件DIgSILENT PowerFactory进行比较。
[0094] 低压配电系统算例电压等级为400V,主馈线通过0.4/10kV变压器接至中压母线Ml 处,变压器采用常用的DYnll联结方式,低压侧设有电容无功补偿系统,主馈线节点间距为 50m。考虑到配电系统的实际运行特点,测试算例中既有三相对称和不对称负荷,也有单相 负荷。另外,算例中接入了多种类型的分布式电源,包括:具备最大功率跟踪控制的光伏发 电系统,燃料电池发电系统和蓄电池储能系统,各分布式电源控制方式、接入容量及输出功 率如表1所示。
[0095]表1分布式电源控制方式、接入容量及输出功率
Figure CN103700036BD00101
[0097] 采用本发明提出的投影积分算法对测试系统算例进行动态仿真计算,设置仿真时 间为5s,仿真步长为0.3ms,2. Os时刻L1母线发生三相短路故障,2. Is时刻故障清除。关于算 法参数的选取本发明以k=6,M=2,k=8,M=2和k=8,M=3为例对该算法进行测试,在具体实现 时,算法参数可根据实际应用情况,在满足数值稳定性和数值精度的条件下进行任意取值, 本发明的实施对此不做限制。执行仿真计算的硬件平台为Intel(R)C 〇re(TM)i5-3470CPU@ 3 • 20GHz,4GB RAM的PC机;软件环境为32位Windows7操作系统。
[0098] 将本发明提出的一种适于电力系统多时间尺度的暂态稳定性投影积分方法与步 长取0. lms时DIgSILENT PowerFactory的仿真结果进行比较,仿真结果如附图3至附图5所 示,附图6为不同参数投影算法较DIgSILENT仿真结果的相对误差。比较分析结果表明投影 积分算法具有较好的精度,其仿真结果与DIgSILENT吻合,而显式交替求解方法由于存在一 定的交接误差,导致计算结果存在一定的误差,且随着k值的减小或M值的增大,投影积分算 法与DIgSILENT仿真结果的偏差逐渐增大,但最大误差维持在0.15%以内,显示出该算法具 有良好的数值精度。
[0099] 本发明提出的一种适于电力系统多时间尺度的暂态稳定性投影积分方法可以实 现对其内部积分算法,即显式四阶龙格-库塔法仿真速度的提升,将选取不同参数的投影积 分算法的仿真计算时间与同步长显式四阶龙格-库塔法以及商业软件DIgSILENT定步长仿 真进行比较,如表2所示。
[0100] 表2算法性能比较
Figure CN103700036BD00111
[0102]表2给出了选取不同参数时投影积分算法与同步长显式四阶龙格-库塔法以及商 业软件DIgSILENT定步长仿真计算时间的对比。可以看出,投影积分算法的仿真计算时间明 显小于传统显式四阶龙格-库塔法和DIgSILENT,且随着k值的减小或M值的增大,投影算法 仿真速度的提升越明显,相比显式四阶龙格-库塔法其最大加速比可达47.9%。由此可见,本 发明提出的一种适于电力系统多时间尺度的暂态稳定性投影积分方法可以实现仿真计算 效率的提升。

Claims (4)

1. 一种适于电力系统多时间尺度的暂态稳定性投影积分方法,其特征在于,包括如下 步骤: 1) 输入电力系统潮流数据、动态元件数据以及故障及操作设置信息,设置仿真参数; 2) 根据系统算例参数和动态数据形成全系统节点导纳矩阵,其中包括系统节点导纳矩 阵和并入动态元件的节点导纳矩阵,并重新进行潮流校核; 3) 对系统算例的动态元件进行仿真初始化计算; 4) 设置仿真计算时间t = 0; 5) 设置当前投影积分算法内部积分器的积分步数s = 1; 6) 采用投影积分算法内部积分器对电力系统模型积分一个步长,步长为h,然后设置仿 真计算时间t = t+h,内部积分器的积分步数s = s+l,内部积分器采用显式交替求解方法对 微分-代数方程组进行求解,对其中的微分方程求解则采用显式四阶龙格-库塔法; 7) 根据步骤1)输入的故障及操作设置信息判断系统是否发生故障或操作,若发生,则 返回步骤5),否则进入下一步骤; 8) 判断仿真计算时间t是否大于设置时间T,若计算时间t大于设置时间T,则仿真结束, 否则进入下一步骤; 9) 判断投影积分算法内部积分步数s是否大于设置步数k+Ι,若不大于,则返回步骤6), 否则进入下一步骤; 10) 根据步骤1)输入的故障及操作设置信息判断未来Mh时间内是否存在故障或操作, 若存在,则返回步骤5),否则进入下一步骤; 11) 判断未来Mh时间内仿真计算时间t是否会超过设置时间T,若超过,则返回步骤5), 否则进入下一步骤; 12) 根据内部积分器得到的x(tn+k)和x(tn+k+1)的值进行一步外部积分步,步长为内部积 分器积分步长的Μ倍,利用式X(t n+k+l+M) = (M+l)X(tn+k+l)-Mx(tn+k)得到X(tn+k+l+M),然后求解 代数方程g(x(tn+k+i+M),y(t n+k+i+M)) =0得到y(tn+k+i+M),设置t = t+Mh; 13) 判断此时仿真计算时间t是否达到设置的仿真时间T,若达到,则仿真结束,否则返 回步骤5),依此反复进行直至仿真结束。
2. 根据权利要求1所述的一种适于电力系统多时间尺度的暂态稳定性投影积分方法, 其特征在于,步骤1)中所述的仿真参数,包括仿真时间Τ、仿真步长h和投影积分算法参数, 其中,所述的投影积分算法参数包括内部积分器的积分步数k和外部积分器步长相对内部 积分器步长的倍数Μ。
3. 根据权利要求2所述的一种适于电力系统多时间尺度的暂态稳定性投影积分方法, 其特征在于,内部积分器算法是: (1 )采用显式四阶龙格-库塔法积分k步,步长为h,时间从tn到tn+k,由X ( tn)得到X( tn+k), 递推公式如下;
Figure CN103700036BC00021
其中 kl = hf (x(tn) ,tn)
Figure CN103700036BC00031
k4 = hf (x(tn)+k3,tn+h); (2)再次采用显式四阶龙格-库塔法积分一步,步长仍为h,得到x(tn+k+1)。
4.根据权利要求2所述的一种适于电力系统多时间尺度的暂态稳定性投影积分方法, 其特征在于,外部积分器算法是: 根据内部积分器得到的x(tn+k)和x(tn+k+1)的值,利用下式进行一步外部积分步,步长为 内部积分步长的Μ倍,时间从tn+k+l到tn+k+1+M,得到X(tn+k+l+M): x(tn+k+l+M) = (M+l)x(tn+k+l)-Mx(tn+k) 〇
CN201310730821.6A 2013-12-19 2013-12-19 一种适于电力系统多时间尺度的暂态稳定性投影积分方法 Expired - Fee Related CN103700036B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310730821.6A CN103700036B (zh) 2013-12-19 2013-12-19 一种适于电力系统多时间尺度的暂态稳定性投影积分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310730821.6A CN103700036B (zh) 2013-12-19 2013-12-19 一种适于电力系统多时间尺度的暂态稳定性投影积分方法

Publications (2)

Publication Number Publication Date
CN103700036A CN103700036A (zh) 2014-04-02
CN103700036B true CN103700036B (zh) 2016-09-07

Family

ID=50361555

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310730821.6A Expired - Fee Related CN103700036B (zh) 2013-12-19 2013-12-19 一种适于电力系统多时间尺度的暂态稳定性投影积分方法

Country Status (1)

Country Link
CN (1) CN103700036B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104156542B (zh) * 2014-08-26 2018-01-02 天津大学 一种基于隐式投影的有源配电系统稳定性仿真方法
CN104376158B (zh) * 2014-11-05 2017-05-10 天津大学 一种用于矩阵指数的暂态仿真多时间尺度输出方法
CN105808825B (zh) * 2016-02-29 2018-09-04 天津大学 一种适于有源配电网动态仿真的滚动式随机投影积分方法
CN106405330A (zh) * 2016-08-31 2017-02-15 广西电网有限责任公司电力科学研究院 一种基于暂态电流投影分量极性比较的故障选线方法
CN106443330A (zh) * 2016-09-07 2017-02-22 中国石油大学(华东) 一种基于暂态投影分量幅值与极性关系比较的小电流接地系统高阻接地故障选线方法
CN111984046A (zh) * 2020-08-19 2020-11-24 国网山西省电力公司 一种基于仿真算法环境主动预警及调节系统
CN112464604B (zh) * 2021-01-28 2021-08-20 通号(长沙)轨道交通控制技术有限公司 一种多变流器集中并网系统高效仿真方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Projective methods for stiff differential equations:problems with gaps in their eigenvalue spectrum;Gear C W,et al.;《SIAM Journal on Scientific Computing》;20030131;第24卷(第4期);全文 *
分布式发电、微网与智能配电网的发展与挑战;王成山等;《电力系统自动化》;20100125;第34卷(第2期);全文 *
计及保护的电力系统暂态稳定性评价的研究;杨静;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑(月刊)》;20090915;全文 *

Also Published As

Publication number Publication date
CN103700036A (zh) 2014-04-02

Similar Documents

Publication Publication Date Title
CN103700036B (zh) 一种适于电力系统多时间尺度的暂态稳定性投影积分方法
CN103942372B (zh) 基于fpga的有源配电网暂态实时仿真多速率接口方法
CN101882167B (zh) 一种大规模风电集中接入电网的风电场等值建模方法
CN102819641B (zh) 适于电磁暂态仿真的大规模配电网络整体模型化简方法
CN102012956B (zh) 一种考虑风电场输入风速和风向随机波动的风电场等值方法
CN105932666A (zh) 复杂配电网多时间尺度数模混合仿真系统及其仿真方法
CN102436534B (zh) 一种电力系统仿真方法
CN103646152B (zh) 一种基于矩阵指数的电力系统电磁暂态仿真方法
CN103049617B (zh) 保留无源性的大规模配电网络电磁暂态仿真模型化简方法
CN102346442B (zh) 适于节点分析框架的微网暂态仿真并行计算方法
CN103077268B (zh) 面向电力系统电磁暂态仿真的状态空间自动建模方法
CN204314716U (zh) 一种用于风电场控制系统硬件在环测试平台
Jalili-Marandi et al. A real-time dynamic simulation tool for transmission and distribution power systems
CN104732033A (zh) 一种电力系统机电暂态过程的仿真分析方法
CN104298809B (zh) 一种基于矩阵指数电磁暂态仿真的非线性建模求解方法
CN104732459A (zh) 大规模电力系统病态潮流分析系统
CN104156542B (zh) 一种基于隐式投影的有源配电系统稳定性仿真方法
CN103810646B (zh) 一种基于改进投影积分算法的有源配电系统动态仿真方法
CN105574809A (zh) 基于矩阵指数的电磁暂态仿真图形处理器并行计算方法
CN103440374A (zh) 基于状态变量节点分析混合框架的电磁暂态仿真建模方法
CN108448631A (zh) 基于模型降阶的含分布式电源接入配电网的动态相量建模方法
CN203102283U (zh) 一种基于实时数字仿真仪rtds的风电场建模仿真系统
CN106126871A (zh) 一种pscad电磁暂态仿真用的调速器模型建模方法
CN104361170B (zh) 适于大规模风电并网暂态稳定分析的风电场等值建模方法
CN105787604A (zh) 一种输配协同的配网态势快速感知方法

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160907

Termination date: 20191219

CF01 Termination of patent right due to non-payment of annual fee