CN103810646A - 一种基于改进投影积分算法的有源配电***动态仿真方法 - Google Patents
一种基于改进投影积分算法的有源配电***动态仿真方法 Download PDFInfo
- Publication number
- CN103810646A CN103810646A CN201410020279.XA CN201410020279A CN103810646A CN 103810646 A CN103810646 A CN 103810646A CN 201410020279 A CN201410020279 A CN 201410020279A CN 103810646 A CN103810646 A CN 103810646A
- Authority
- CN
- China
- Prior art keywords
- algorithm
- simulation
- distribution system
- integral
- dynamic
- 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.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种基于改进投影积分算法的有源配电***动态仿真方法:输入有源配电***算例潮流数据、动态元件数据以及故障及操作设置信息,设置仿真参数;根据***算例参数和动态数据形成全***节点导纳矩阵;仿真初始化计算;设置仿真时间t=0;设置积分步数s=1;作tn~tn+1时段计算;判断***是否发生故障或操作;判断仿真时间t是否大于设置时间T;判断改进投影积分算法内部积分步数s是否大于设置步数k+1;判断未来Mh时间内是否存在故障或操作;设置仿真时间t=t+Mh;执行改进投影积分算法的外部积分器校正步;判断仿真时间t是否达到仿真设置时间T。本发明的方法数值精度高、数值稳定性好、计算效率高、适用于具有刚性特征的有源配电***动态仿真。
Description
技术领域
本发明涉及一种有源配电***动态仿真方法。特别是涉及一种基于改进投影积分算法的有源配电***动态仿真方法。
背景技术
分布式电源(DG)的大规模广泛接入以及需求侧响应技术实施后***负荷特性的改变,无不对配电***的规划与运行带来新的挑战。不含DG的配电网是“被动的”,接入用户所使用的电能由上一级输电网提供,当配电网接入DG产生双向潮流时,称该***为“有源配电***”。有源配电***是具备组合控制各种分布式能源(DER,如DG、可控负荷、储能等)能力的复杂配电***。
在未来的有源配电***中,DG的接入容量可以轻易超过(至少在特定时间段内)配电***中的负荷总量,此时有源配电***作为外部电源向外部输电网络输送电能。即使DG总容量不超过负荷总量,高渗透率下DG仍会造成配电网络动态响应特性发生变化进而影响整个电力***的动态特性,特别是故障特性。在***层面,相关问题的分析与研究往往无法直接在实际***上进行试验,因此必须采用有效的数字仿真工具作为重要的研究手段。
传统电力***时域仿真针对***动态过程的不同时间尺度分别发展出电磁暂态仿真、机电暂态仿真和中长期动态仿真三种电力***数字仿真方法,三者从元件数学模型到仿真计算方法都具有不同的特征。电力***电磁暂态仿真侧重于***中电场与磁场相互影响产生的电压电流的变化过程;机电暂态仿***要研究电力***在大扰动下(如故障、切机、切负荷、重合闸操作等情况)的动态行为和保持同步稳定运行的能力,即暂态稳定性,所关注的时间范围通常为几秒至几十秒,因而也称为暂态稳定性仿真;中长期动态过程仿真是电力***受到扰动后较长过程的动态仿真,即通常的电力***长过程动态稳定计算。
电力***暂态稳定性仿真除关注传统电力***的暂态稳定运行能力外,近年来还着重于分析含各种分布式电源及储能装置的有源配电***运行时其工频电气量在***扰动下(开关操作、故障、分布式电源及负荷波动等)的动态响应特性,此时可称为有源配电***动态仿真。有源配电***接入了种类繁多的分布式电源和大量的电力电子装置,包括旋转电机和各种静态直流型分布式电源,具有明显的多时间尺度特征。有源配电***动态仿真的多时间尺度特征在数学上可以归结为刚性问题,其对所采用的数值算法的精度和稳定性要求更高。
有源配电***动态仿真本质上可归结为对动力学***时域响应的求取,分为数学建模和模型求解两部分。动态仿真首先根据元件间的拓扑关系将有源配电***各元件模型构成全***模型,形成一组联立的微分-代数方程组,然后以稳态工况或潮流解为初值,求解扰动下的数值解,即逐步求得***状态量和代数量随时间的变化曲线。一般可将有源配电***数学模型通过一个高维非线性且连续自治的微分-代数方程组来描述,如式(1)所示。
式中,为微分方程,为代数方程,为***状态变量,代表电机转子转速、电力电子器件控制***和负荷动态参数等,为代数变量,表征母线电压幅值与相角。数学模型的求解一般通过特定的数值算法和相应的仿真程序来实现。因此,有源配电***动态仿真在数学上可以归结为求解一个微分-代数方程组的初值问题。
有源配电***动态仿真算法按照对于式(1)中微分方程和代数方程解算形式的不同可以分为交替求解法和联立求解法两大类。交替求解法首先采用特定的数值积分算法,根据初始化计算结果求解微分方程,得到本时步状态变量的值,然后将其代入到代数方程中求解,得到该时步代数变量的值,最后再将代数变量代入微分方程进行下一时步状态变量求解,以此类推实现微分-代数方程组的交替求解;联立求解法则是将微分方程差分化之后,和代数方程联立成一个完整的代数方程组,同时求解状态变量和代数变量。对于式(1)中的微分方程,除少数可得到解析解以外,大多数只能采用数值解法进行求解,其中,单步法在有源配电***动态仿真中应用广泛。根据求解过程的不同,单步法可分为显式积分方法和隐式积分方法,显式积分方法可根据当前时刻状态变量直接计算下一时刻状态变量,而隐式积分方法则需要对含有当前时刻和下一时刻状态变量的方程进行求解才能求得下一时刻状态变量。常见的显式积分方法包括欧拉法、改进欧拉法和龙格-库塔法,而隐式积分方法主要有后向欧拉法和隐式梯形法。
有源配电***的刚性特征使得传统积分方法的应用存在诸多缺陷。对于显式积分方法,其每一时步内的运算量较小,但由于显式方法的数值稳定性较差,因此针对刚性问题求解往往只能采取较小的仿真步长,仿真速度受到极大限制。隐式积分方法虽然数值稳定性较好,可以在刚性问题的求解过程中保证数值稳定性,但其每一时步都需要迭代求解方程组,相比显式积分算法其计算与编程工作复杂,限制了其在有源配电***动态仿真中的应用。
可见,提出一种数值精度高、数值稳定性好、计算效率高、适于具有刚性特征的有源配电***的动态仿真方法十分重要。
发明内容
本发明所要解决的技术问题是,提供一种数值精度高、数值稳定性好、计算效率高且适于刚性有源配电***的基于改进投影积分算法的有源配电***动态仿真方法。
本发明所采用的技术方案是:一种基于改进投影积分算法的有源配电***动态仿真方法,包括如下步骤:
1)输入有源配电***算例潮流数据、动态元件数据以及故障及操作设置信息,设置仿真参数;
2)根据***算例参数和动态数据形成全***节点导纳矩阵,其中包括***节点导纳矩阵和并入动态元件的节点导纳矩阵,并重新进行潮流校核;
3)对***中的动态元件进行仿真初始化计算;
4)设置仿真时间t=0;
5)设置当前改进投影积分算法内部积分器的积分步数s=1;
6)作tn~tn+1时段计算,设置仿真时间t=t+h,h为仿真步长,采用改进投影积分算法内部积分器对有源配电***模型积分一个步长得到tn+1时刻***状态变量和代数变量,然后设置内部积分器的积分步数s=s+1,内部积分器采用显式交替求解方法对微分-代数方程组进行求解,对其中的微分方程求解则采用显式四阶龙格-库塔法;
7)根据步骤1)输入的故障及操作设置信息判断***是否发生故障或操作,若发生,则返回步骤5),否则进入下一步骤;
8)判断仿真时间t是否大于设置时间T,若仿真时间t大于设置时间T,则仿真结束,否则进入下一步骤;
9)判断改进投影积分算法内部积分步数s是否大于设置步数k+1,若不大于,则返回步骤6),否则进入下一步骤;
10)根据步骤1)输入的故障及操作设置信息判断未来Mh时间内是否存在故障或操作,若存在,则返回步骤5),否则进入下一步骤;
11)设置仿真时间t=t+Mh,执行改进投影积分算法的外部积分器预测步,根据内部积分器得到的x(tn+k+1)的值,利用前向欧拉法得到tn+k+1+M时刻***状态变量的预测值x(0)(tn+k+1+M),步长为内部积分器积分步长的M倍,迭代求解代数方程g(x(0)(tn+k+1+M),y(0)(tn+k+1+M))=0得到代数变量y(0)(tn+k+1+M);
12)执行改进投影积分算法的外部积分器校正步,对预测值x(0)(tn+k+1+M)进行校正得到tn+k+1+M时刻的***状态变量x(tn+k+1+M),迭代求解代数方程g(x(tn+k+1+M),y(tn+k+1+M))=0得到代数变量y(tn+k+1+M),改进投影积分算法外部积分器计算结束;
13)判断仿真时间t是否达到仿真设置时间T,若达到,则仿真结束,否则返回步骤5),依此反复进行直至仿真结束。
步骤1)中所述的仿真参数,包括仿真设置时间T、仿真步长h和投影积分算法参数,其中,所述的投影积分算法参数包括内部积分器的积分步数k和外部积分器步长相对内部积分器步长的倍数M。
步骤6)中所述的内部积分器采用显式交替求解方法对微分-代数方程组进行求解,对其中的微分方程求解则采用显式四阶龙格-库塔法,具体步骤是:
(1)采用显式四阶龙格-库塔法积分k步,步长为h,时间从tn到tn+k,由x(tn)得到x(tn+k),递推公式如下,同时迭代求解代数方程由代数变量y(tn)得到y(tn+k);
其中
k1=hf(x(tn),tn)
k4=hf(x(tn)+k3,tn+h)
(2)再次采用显式四阶龙格-库塔法积分一步,步长仍为h,得到x(tn+k+1),同时迭代求解代数方程得到代数变量y(tn+k+1)。
步骤11)所述的外部积分器预测步的算法是:
根据内部积分器得到的x(tn+k+1)的值,利用前向欧拉法得到tn+k+1+M时刻***状态变量的预测值x(0)(tn+k+1+M),步长为内部积分器积分步长的M倍,即Mh,如下式所示:
步骤12)所述的外部积分器校正步的算法是:
采用下式对预测值x(0)(tn+k+1+M)进行校正,得到tn+k+1+M时刻的***状态变量x(tn+k+1+M):
上式也可写为
本发明的一种基于改进投影积分算法的有源配电***动态仿真方法,考虑了有源配电***的刚性特征,采用显式交替求解方法对描述有源配电***模型的微分-代数方程组进行交替求解,对其中的微分方程利用改进的投影积分算法进行差分求解。本发明的方法为2阶精度算法,适用于具有刚性特征的有源配电***动态仿真问题求解,相比传统数值积分方法大大提升了仿真计算速度,为高效、可靠的有源配电***仿真程序的开发奠定了良好的基础。本发明的方法数值精度高、数值稳定性好、计算效率高、适用于具有刚性特征的有源配电***动态仿真。
附图说明
图1是本发明方法的整体流程图;
图2是低压有源配电***算例结构图;
图中1:燃料电池;2:第二蓄电池;3:第二光伏电池;4:第一光伏电池;5:第一蓄电池;M1:中压母线;L1~L19:低压母线;Load1~Load7:负荷;
图3是L16母线电压仿真结果及局部放大图;
图4是L17母线电压仿真结果及局部放大图;
图5是L17母线频率仿真结果及局部放大图;
图6是L17母线电压改进投影积分算法仿真结果较DIgSILENT相对误差。
具体实施方式
下面结合实施例和附图对本发明的一种基于改进投影积分算法的有源配电***动态仿真方法做出详细说明。
本发明的一种基于改进投影积分算法的有源配电***动态仿真方法,属于显式积分方法。有源配电***接入了种类繁多的分布式电源和大量的电力电子装置,包括旋转电机和各种静态直流型分布式电源,具有明显的多时间尺度特征。因此,有源配电***动态仿真的多时间尺度特征在数学上可以归结为刚性问题,其对所采用的数值算法的数值精度和数值稳定性要求更高。本发明的方法考虑了有源配电***的多时间尺度特征,采用显式交替求解方法对有源配电***模型进行交替求解,对其中的微分方程利用改进的投影积分算法进行差分求解。本发明的方法为2阶精度算法,适用于具有刚性特征的有源配电***动态仿真问题求解,相比传统数值积分算法大大提升了仿真计算速度,为高效、可靠的有源配电***仿真程序的开发奠定了良好的基础。
本发明采用显式交替求解算法对通过微分-代数方程组描述的有源配电***数学模型进行求解,对其中的微分方程求解则采用改进的投影积分算法:首先进行若干小步长的积分计算,对应***的快动态过程;而后根据小步长的积分计算结果,先进行一步大步长的预测步计算得到下一时刻状态变量的预测值,然后再进行一步校正步对该预测值进行校正,该过程与***的慢动态过程对应。其中,小步长积分计算过程称为内部积分器,采用数值稳定性较好且精度较高的显式四阶龙格-库塔法(explicit four-order Runge-Kutta method)以提高算法的稳定性和数值精度;大步长的预测-校正积分计算过程称为外部积分器,通过预测-校正投影积分计算实现在保证数值精度的同时,提升算法的计算效率。
如图1所示,本发明的一种基于改进投影积分算法的有源配电***动态仿真方法,包括如下步骤:
1)输入有源配电***算例潮流数据、动态元件数据以及故障及操作设置信息,设置仿真参数;所述的仿真参数包括仿真设置时间T、仿真步长h和投影积分算法参数,其中,所述的投影积分算法参数包括内部积分器的积分步数k和外部积分器步长相对内部积分器步长的倍数M,其中k和M均为正整数。
2)根据***算例参数和动态数据形成全***节点导纳矩阵,其中包括***节点导纳矩阵和并入动态元件的节点导纳矩阵,并重新进行潮流校核;
3)对***中的动态元件进行仿真初始化计算;
4)设置仿真时间t=0;
5)设置当前改进投影积分算法内部积分器的积分步数s=1;
6)作tn~tn+1时段计算,设置仿真时间t=t+h,h为仿真步长,采用改进投影积分算法内部积分器对有源配电***模型积分一个步长得到tn+1时刻***状态变量和代数变量,然后设置内部积分器的积分步数s=s+1,内部积分器采用显式交替求解方法对微分-代数方程组进行求解,对其中的微分方程求解则采用显式四阶龙格-库塔法;
所述的有源配电***模型是下式所示的一个高维非线性且连续自治的微分-代数方程组:
式中,为微分方程,用于描述有源配电***设备动态特征;为代数方程,用于描述设备之间的电气联系;为***状态变量,代表电机转子转速、电力电子器件控制***和负荷动态参数等;为代数变量,表征母线电压幅值与相角。
所述的内部积分器采用显式交替求解方法对微分-代数方程组所描述的有源配电***模型进行求解,对其中的微分方程求解则采用显式四阶龙格-库塔法,具体步骤是:
(1)采用显式四阶龙格-库塔法积分k步,步长为h,时间从tn到tn+k,由x(tn)得到x(tn+k),递推公式如下,同时迭代求解代数方程由代数变量y(tn)得到y(tn+k);
其中
k1=hf(x(tn),tn)
k4=hf(x(tn)+k3,tn+h)
(2)再次采用显式四阶龙格-库塔法积分一步,步长仍为h,得到x(tn+k+1),同时迭代求解代数方程得到代数变量y(tn+k+1)。
采用显式交替求解方法对微分-代数方程组所描述的有源配电***模型进行求解,包括如下步骤:
(1)对有源配电***模型显式差分化后得到下式:
(2)根据tn时刻的状态变量x(tn)和代数变量y(tn),计算向量k1=hf(x(tn),y(tn));
(6)计算tn+1时刻然后求解代数方程g(x(tn+1),y(tn+1))=0得出y(tn+1)。
7)根据步骤1)输入的故障及操作设置信息判断***是否发生故障或操作,若发生,则返回步骤5),否则进入下一步骤;
8)判断仿真时间t是否大于设置时间T,若仿真时间t大于设置时间T,则仿真结束,否则进入下一步骤;
9)判断改进投影积分算法内部积分步数s是否大于设置步数k+1,若不大于,则返回步骤6),否则进入下一步骤;
10)根据步骤1)输入的故障及操作设置信息判断未来Mh时间内是否存在故障或操作,若存在,则返回步骤5),否则进入下一步骤;
11)设置仿真时间t=t+Mh,执行改进投影积分算法的外部积分器预测步,根据内部积分器得到的x(tn+k+1)的值,利用前向欧拉法得到tn+k+1+M时刻***状态变量的预测值x(0)(tn+k+1+M),步长为内部积分器积分步长的M倍,迭代求解代数方程g(x(0)(tn+k+1+M),y(0)(tn+k+1+M))=0得到代数变量y(0)(tn+k+1+M);
所述的外部积分器预测步的算法是:根据内部积分器得到的x(tn+k+1)的值,利用前向欧拉法得到tn+k+1+M时刻***状态变量的预测值x(0)(tn+k+1+M),步长为内部积分器积分步长的M倍,即Mh,如下式所示:
12)执行改进投影积分算法的外部积分器校正步,对预测值x(0)(tn+k+1+M)进行校正得到tn+k+1+M时刻的***状态变量x(tn+k+1+M),迭代求解代数方程g(x(tn+k+1+M),y(tn+k+1+M))=0得到代数变量y(tn+k+1+M),改进投影积分算法外部积分器计算结束;
所述的外部积分器校正步的算法是:采用下式对预测值x(0)(tn+k+1+M)进行校正,得到tn+k+1+M时刻的***状态变量x(tn+k+1+M):
上式也可写为
使得外部积分器预测步计算完成后,不必存储x(tn+k+1)的值,有利于节省内存。
13)判断仿真时间t是否达到仿真设置时间T,若达到,则仿真结束,否则返回步骤5),依此反复进行直至仿真结束。
下面给出具体实例:
数字仿真和电网计算程序(DIgSILENT PowerFactory)是德国DIgSLENTGmbH公司开发的一款商业电力***仿真软件。本实例以C++编程语言环境为基础,实现了适于有源配电***的改进投影积分动态仿真方法,通过含分布式电源的低压有源配电***算例(附图2)对本发明的方法进行测试验证,并与DIgSILENT PowerFactory进行比较。
低压有源配电***算例电压等级为400V,主馈线通过0.4/10kV变压器接至中压母线M1处,变压器采用常用的DYn11联结方式,低压侧设有无功补偿电容,主馈线节点间距为50m。考虑到配电***的实际运行特点,测试算例中既有三相对称和不对称负荷,也有单相负荷。另外,算例中接入了多种类型的分布式电源,包括:具备最大功率跟踪控制的光伏发电***,燃料电池发电***和蓄电池储能***,各分布式电源控制方式、接入容量及输出功率如表1所示。
表1 分布式电源控制方式、接入容量及输出功率
采用本发明提出的一种基于改进投影积分算法的有源配电***动态仿真方法对测试算例进行动态仿真计算,设置仿真时间为5s,仿真步长为0.3ms,2.0s时刻L1母线发生三相短路故障,2.1s时刻故障清除。关于算法参数的选取本发明以k=6,M=2,k=4,M=2和k=4,M=7为例对该算法进行测试,在具体实现时,算法参数可根据实际应用情况,在满足数值稳定性和数值精度的条件下进行任意取值,本发明的实施对此不做限制。执行仿真计算的硬件平台为Intel(R)Core(TM)i5-3470CPU3.20GHz,4GB RAM的PC机;软件环境为32位Windows7操作***。
将本发明的一种基于改进投影积分算法的有源配电***动态仿真方法与步长取0.1ms时DIgSILENT PowerFactory的仿真结果进行比较,仿真结果如附图3至附图5所示,附图6为不同参数投影算法较DIgSILENT仿真结果的相对误差。比较分析结果表明改进投影积分算法具有良好的精度,其仿真结果与DIgSILENT吻合,而显式交替求解方法由于存在一定的交接误差,导致计算结果存在一定的误差,且随着k值的减小或M值的增大,投影积分算法与DIgSILENT仿真结果的偏差逐渐增大,但最大误差维持在0.15%以内,显示出该算法具有良好的数值精度。
本发明的一种基于改进投影积分算法的有源配电***动态仿真方法,可以实现对其内部积分算法,即显式四阶龙格-库塔法仿真速度的提升,将选取不同参数的投影积分算法的仿真计算时间与同步长显式四阶龙格-库塔法以及商业软件DIgSILENT定步长仿真进行比较,如表2所示。
表2 算法性能比较
表2给出了选取不同算法参数时改进的投影积分算法与同步长显式四阶龙格-库塔法以及商业软件DIgSILENT定步长仿真计算时间的对比。可以看出,改进投影积分算法的仿真计算时间明显小于传统显式四阶龙格-库塔法和DIgSILENT,且随着k值的减小或M值的增大,算法仿真速度的提升越明显,相比显式四阶龙格-库塔法其最大加速比可达142.4%。由此可见,本发明提出的一种基于改进投影积分算法的有源配电***动态仿真方法可以实现仿真计算效率的大幅度提升。
Claims (5)
1.一种基于改进投影积分算法的有源配电***动态仿真方法,其特征在于,包括如下步骤:
1)输入有源配电***算例潮流数据、动态元件数据以及故障及操作设置信息,设置仿真参数;
2)根据***算例参数和动态数据形成全***节点导纳矩阵,其中包括***节点导纳矩阵和并入动态元件的节点导纳矩阵,并重新进行潮流校核;
3)对***中的动态元件进行仿真初始化计算;
4)设置仿真时间t=0;
5)设置当前改进投影积分算法内部积分器的积分步数s=1;
6)作tn~tn+1时段计算,设置仿真时间t=t+h,h为仿真步长,采用改进投影积分算法内部积分器对有源配电***模型积分一个步长得到tn+1时刻***状态变量和代数变量,然后设置内部积分器的积分步数s=s+1,内部积分器采用显式交替求解方法对微分-代数方程组进行求解,对其中的微分方程求解则采用显式四阶龙格-库塔法;
7)根据步骤1)输入的故障及操作设置信息判断***是否发生故障或操作,若发生,则返回步骤5),否则进入下一步骤;
8)判断仿真时间t是否大于设置时间T,若仿真时间t大于设置时间T,则仿真结束,否则进入下一步骤;
9)判断改进投影积分算法内部积分步数s是否大于设置步数k+1,若不大于,则返回步骤6),否则进入下一步骤;
10)根据步骤1)输入的故障及操作设置信息判断未来Mh时间内是否存在故障或操作,若存在,则返回步骤5),否则进入下一步骤;
11)设置仿真时间t=t+Mh,执行改进投影积分算法的外部积分器预测步,根据内部积分器得到的x(tn+k+1)的值,利用前向欧拉法得到tn+k+1+M时刻***状态变量的预测值x(0)(tn+k+1+M),步长为内部积分器积分步长的M倍,迭代求解代数方程g(x(0)(tn+k+1+M),y(0)(tn+k+1+M))=0得到代数变量y(0)(tn+k+1+M);
12)执行改进投影积分算法的外部积分器校正步,对预测值x(0)(tn+k+1+M)进行校正得到tn+k+1+M时刻的***状态变量x(tn+k+1+M),迭代求解代数方程g(x(tn+k+1+M),y(tn+k+1+M))=0得到代数变量y(tn+k+1+M),改进投影积分算法外部积分器计算结束;
13)判断仿真时间t是否达到仿真设置时间T,若达到,则仿真结束,否则返回步骤5),依此反复进行直至仿真结束。
2.根据权利要求1所述的一种基于改进投影积分算法的有源配电***动态仿真方法,其特征在于,步骤1)中所述的仿真参数,包括仿真设置时间T、仿真步长h和投影积分算法参数,其中,所述的投影积分算法参数包括内部积分器的积分步数k和外部积分器步长相对内部积分器步长的倍数M。
3.根据权利要求1所述的一种基于改进投影积分算法的有源配电***动态仿真方法,其特征在于,步骤6)中所述的内部积分器采用显式交替求解方法对微分-代数方程组进行求解,对其中的微分方程求解则采用显式四阶龙格-库塔法,具体步骤是:
(1)采用显式四阶龙格-库塔法积分k步,步长为h,时间从tn到tn+k,由x(tn)得到x(tn+k),递推公式如下,同时迭代求解代数方程由代数变量y(tn)得到y(tn+k);
其中
k1=hf(x(tn),tn)
k4=hf(x(tn)+k3,tn+h)
(2)再次采用显式四阶龙格-库塔法积分一步,步长仍为h,得到x(tn+k+1),同时迭代求解代数方程得到代数变量y(tn+k+1)。
4.根据权利要求1所述的一种基于改进投影积分算法的有源配电***动态仿真方法,其特征在于,步骤11)所述的外部积分器预测步的算法是:
根据内部积分器得到的x(tn+k+1)的值,利用前向欧拉法得到tn+k+1+M时刻***状态变量的预测值x(0)(tn+k+1+M),步长为内部积分器积分步长的M倍,即Mh,如下式所示:
5.根据权利要求1所述的一种基于改进投影积分算法的有源配电***动态仿真方法,其特征在于,步骤12)所述的外部积分器校正步的算法是:
采用下式对预测值x(0)(tn+k+1+M)进行校正,得到tn+k+1+M时刻的***状态变量x(tn+k+1+M):
上式也可写为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410020279.XA CN103810646B (zh) | 2014-01-16 | 2014-01-16 | 一种基于改进投影积分算法的有源配电***动态仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410020279.XA CN103810646B (zh) | 2014-01-16 | 2014-01-16 | 一种基于改进投影积分算法的有源配电***动态仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103810646A true CN103810646A (zh) | 2014-05-21 |
CN103810646B CN103810646B (zh) | 2017-02-15 |
Family
ID=50707372
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410020279.XA Active CN103810646B (zh) | 2014-01-16 | 2014-01-16 | 一种基于改进投影积分算法的有源配电***动态仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103810646B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156542A (zh) * | 2014-08-26 | 2014-11-19 | 天津大学 | 一种基于隐式投影的有源配电***稳定性仿真方法 |
CN105808825A (zh) * | 2016-02-29 | 2016-07-27 | 天津大学 | 一种适于有源配电网动态仿真的滚动式随机投影积分方法 |
CN109245151A (zh) * | 2018-08-31 | 2019-01-18 | 中国电力科学研究院有限公司 | 一种面向大规模分布式发电的仿真方法及装置 |
CN110705066A (zh) * | 2019-09-20 | 2020-01-17 | 天津大学 | 基于投影积分的气电耦合园区综合能源***动态仿真方法 |
CN111260333A (zh) * | 2020-02-11 | 2020-06-09 | 长江师范学院 | 一种基于***动力学仿真的高校基建项目管理方法 |
CN111900718A (zh) * | 2020-06-17 | 2020-11-06 | 东南大学 | 基于多级优化追赶变分迭代法的有源配电网动态仿真方法 |
CN111984046A (zh) * | 2020-08-19 | 2020-11-24 | 国网山西省电力公司 | 一种基于仿真算法环境主动预警及调节*** |
CN113036801A (zh) * | 2021-02-25 | 2021-06-25 | 东北大学 | 一种考虑时间特性的分布式多能***及其控制方法 |
CN114114958A (zh) * | 2021-11-17 | 2022-03-01 | 天津大学 | 一种考虑源荷温度控制特性的供热***状态离散仿真方法 |
CN115688447A (zh) * | 2022-11-08 | 2023-02-03 | 南方电网数字电网研究院有限公司 | 一种高性能电力***云仿真***架构 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103400207A (zh) * | 2013-08-01 | 2013-11-20 | 天津大学 | 一种含可调度分布式电源的配电网运行优化方法 |
CN103440374A (zh) * | 2013-08-23 | 2013-12-11 | 中国南方电网有限责任公司电网技术研究中心 | 基于状态变量节点分析混合框架的电磁暂态仿真建模方法 |
-
2014
- 2014-01-16 CN CN201410020279.XA patent/CN103810646B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103400207A (zh) * | 2013-08-01 | 2013-11-20 | 天津大学 | 一种含可调度分布式电源的配电网运行优化方法 |
CN103440374A (zh) * | 2013-08-23 | 2013-12-11 | 中国南方电网有限责任公司电网技术研究中心 | 基于状态变量节点分析混合框架的电磁暂态仿真建模方法 |
Non-Patent Citations (2)
Title |
---|
C.W.GEAR 等: "projective methods for stiff differential equations:problems with gaps in their eigenvalue spectrum", 《SOCIETY FOR INDUSTRIAL AND APPLIED MATHEMATICS》 * |
王成山: "分布式发电***电力电子控制器通用建模方法", 《电力***自动化》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156542A (zh) * | 2014-08-26 | 2014-11-19 | 天津大学 | 一种基于隐式投影的有源配电***稳定性仿真方法 |
CN104156542B (zh) * | 2014-08-26 | 2018-01-02 | 天津大学 | 一种基于隐式投影的有源配电***稳定性仿真方法 |
CN105808825A (zh) * | 2016-02-29 | 2016-07-27 | 天津大学 | 一种适于有源配电网动态仿真的滚动式随机投影积分方法 |
CN105808825B (zh) * | 2016-02-29 | 2018-09-04 | 天津大学 | 一种适于有源配电网动态仿真的滚动式随机投影积分方法 |
CN109245151B (zh) * | 2018-08-31 | 2022-09-02 | 中国电力科学研究院有限公司 | 一种面向大规模分布式发电的仿真方法及装置 |
CN109245151A (zh) * | 2018-08-31 | 2019-01-18 | 中国电力科学研究院有限公司 | 一种面向大规模分布式发电的仿真方法及装置 |
CN110705066A (zh) * | 2019-09-20 | 2020-01-17 | 天津大学 | 基于投影积分的气电耦合园区综合能源***动态仿真方法 |
CN111260333A (zh) * | 2020-02-11 | 2020-06-09 | 长江师范学院 | 一种基于***动力学仿真的高校基建项目管理方法 |
CN111900718A (zh) * | 2020-06-17 | 2020-11-06 | 东南大学 | 基于多级优化追赶变分迭代法的有源配电网动态仿真方法 |
CN111900718B (zh) * | 2020-06-17 | 2022-08-09 | 东南大学 | 基于多级优化追赶变分迭代法的有源配电网动态仿真方法 |
CN111984046A (zh) * | 2020-08-19 | 2020-11-24 | 国网山西省电力公司 | 一种基于仿真算法环境主动预警及调节*** |
CN113036801A (zh) * | 2021-02-25 | 2021-06-25 | 东北大学 | 一种考虑时间特性的分布式多能***及其控制方法 |
CN114114958A (zh) * | 2021-11-17 | 2022-03-01 | 天津大学 | 一种考虑源荷温度控制特性的供热***状态离散仿真方法 |
CN114114958B (zh) * | 2021-11-17 | 2023-09-29 | 天津大学 | 一种考虑源荷温度控制特性的供热***状态离散仿真方法 |
CN115688447A (zh) * | 2022-11-08 | 2023-02-03 | 南方电网数字电网研究院有限公司 | 一种高性能电力***云仿真***架构 |
CN115688447B (zh) * | 2022-11-08 | 2023-10-20 | 南方电网数字电网研究院有限公司 | 一种高性能电力***云仿真***架构 |
Also Published As
Publication number | Publication date |
---|---|
CN103810646B (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103810646B (zh) | 一种基于改进投影积分算法的有源配电***动态仿真方法 | |
CN103646152B (zh) | 一种基于矩阵指数的电力***电磁暂态仿真方法 | |
CN103700036B (zh) | 一种适于电力***多时间尺度的暂态稳定性投影积分方法 | |
CN104298809B (zh) | 一种基于矩阵指数电磁暂态仿真的非线性建模求解方法 | |
CN103077268B (zh) | 面向电力***电磁暂态仿真的状态空间自动建模方法 | |
CN104156542A (zh) | 一种基于隐式投影的有源配电***稳定性仿真方法 | |
CN106329523A (zh) | 考虑不确定性的有源配电网智能软开关鲁棒优化建模方法 | |
CN103049617B (zh) | 保留无源性的大规模配电网络电磁暂态仿真模型化简方法 | |
CN103336866B (zh) | 一种电磁暂态仿真中含负电阻支路的处理方法 | |
Shu et al. | A multirate EMT co-simulation of large AC and MMC-based MTDC systems | |
CN104104081B (zh) | 一种基于优化方法的非迭代不确定潮流分析方法 | |
CN104794541A (zh) | 一种基于模拟退火和锥优化的配电网运行优化方法 | |
CN105574809A (zh) | 基于矩阵指数的电磁暂态仿真图形处理器并行计算方法 | |
CN103455948B (zh) | 一种配电***多维多分辨率建模与分析方法 | |
CN103440374A (zh) | 基于状态变量节点分析混合框架的电磁暂态仿真建模方法 | |
CN105184027A (zh) | 一种基于交互式多模型算法的电力负荷建模方法 | |
CN102855382B (zh) | 一种电力***三相短路故障临界切除时间的在线求取方法 | |
Liu et al. | Fast power system dynamic simulation using continued fractions | |
CN105224728A (zh) | 一种含详细发电机模型的电网暂态稳定能量函数分析方法及*** | |
CN104979840A (zh) | 一种主动配电网三相无功优化方法 | |
Dubyago | Mathematical description of interaction of the elements in the power network system | |
CN103714212B (zh) | 一种面向暂态仿真的配电***模型化简误差控制方法 | |
CN106126871A (zh) | 一种pscad电磁暂态仿真用的调速器模型建模方法 | |
CN103226655B (zh) | 适用pwm变流器平均模型的改进emtp算法 | |
CN102609576B (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 |