基于实时积分算法的非线性元件电磁暂态仿真方法和装置制造方法及图纸

本发明专利技术涉及一种基于实时积分算法的非线性元件电磁暂态仿真方法和装置，方法包括：对电力网络中的非线性元件进行建模，建立非线性常微分方程组，并改写为齐次项加非齐次项组合的形式，得到统一形式的数学模型；对数学模型初始化；将数学模型转化为同解的积分方程，对方程的齐次项采用积分算法计算，对方程的非齐次积分项通过Calahan算法计算，并结合Richardson外推法提高Calahan算法的精度，最终得到下一时刻的模型数值解；对模型数值解进行每一时刻的计算，直至达到预设的时间阈值，最终输出电磁暂态数值仿真结果。与现有技术相比，本发明专利技术新算法可以避免求解Jacobian矩阵和求逆计算，计算效率相当于是一种显式的方法，大大提高了Rosenbrock积分算法的计算效率。大大提高了Rosenbrock积分算法的计算效率。大大提高了Rosenbrock积分算法的计算效率。

【技术实现步骤摘要】
基于实时积分算法的非线性元件电磁暂态仿真方法和装置


[0001]本专利技术涉及电力系统电磁暂态仿真
，尤其是涉及基于实时积分算法的非线性元件电磁暂态仿真方法和装置。

技术介绍

[0002]电力系统电磁暂态仿真程序已经成为现代电力系统规划和运行中的一种基础仿真工具，而非线性元件的电磁暂态仿真是一个特殊的挑战，因为它需要精确的表示和有效的解决方案。电力系统中含有大量的非线性元件，例如变压器、电抗器、避雷器、气体放电管、压敏电阻等都是非线性元件或时变元件。对于这些元件许多暂态问题的研究都是用常微分方程来描述其数学模型，这导致非线性电路的复杂性增加，对于非线性常微分方程数值求解中，刚性问题比较常见，从而对数值积分算法的稳定域提出了更高的要求。
[0003]在电力系统电磁暂态仿真中常用的数值积分方法有欧拉法、梯形法、Gear法、Runge
‑
Kutta法(RK算法)等。目前，电力系统电磁暂态仿真程序(electromagnetic transients program，EMTP)所采用数值算法主要为隐式梯形积分法。然而，该方法对于含有较多非线性元件的大型电力网络系统中，尤其是当非线性元件是三角形连接的电路时，无法运行程序进行计算，具有一定的局限性。在数值稳定性上，隐式梯形积分法是A
‑
稳定但不是L
‑
稳定的，不具有数值阻尼特性，因此当网络中状态变量突变时将出现数值振荡现象。并且由于非线性元件产生的非线性隐式方程必须采用各种类型的Newton迭代法来计算，这就使得隐式方法较显式方法而言计算量与计算难度大大提高了。
[0004]为解决传统电磁暂态仿真程序中数值算法存在的问题，研究人员提出了一种2级3阶的Rosenbrock积分方法，即Calahan算法，求解时每向前积分一步只需计算一个Jacobian矩阵，提高了计算速度。上述基于Rosenbrock积分方法的Calahan算法常常被称为半隐式(或半显式)RK算法。对于非线性方程组求解时，不需要迭代，可以直接递推，且具有绝对稳定性，即A
‑
稳定性。但是这类半隐格式算法在求解公式中为了避免与外加非线性函数有关的数值振荡问题，必须计算右函数的解析偏导数，即Jacobian矩阵，求解比较复杂；且求解时需要求解逆矩阵，随着系统方程阶数的提高会导致求逆时间大大增加。

技术实现思路

[0005]本专利技术的目的就是为了克服上述现有技术存在的缺陷而提供一种基于实时积分算法的非线性元件电磁暂态仿真方法和装置，计算效率相当于是一种显式的方法，大大提高了Rosenbrock积分算法的计算效率。
[0006]本专利技术的目的可以通过以下技术方案来实现：
[0007]一种基于实时积分算法的非线性元件电磁暂态仿真方法，包括以下步骤：
[0008]对电力网络中的非线性元件进行建模，并结合电力网络原始参数，建立电力系统m维非线性常微分方程组；
[0009]将所述电力系统m维非线性常微分方程组改写为齐次项加非齐次项组合的形式，
得到统一形式的数学模型；
[0010]对所述数学模型的电磁暂态数值进行初始化；
[0011]将所述数学模型转化为同解的积分方程，对所述积分方程的齐次项采用积分算法计算下一时刻的值，对所述积分方程的非齐次积分项通过Calahan算法计算下一时刻值，并结合Richardson外推法提高Calahan算法的精度，最终得到下一时刻的模型数值解；
[0012]若当前模型数值解所对应的时刻小于预设的时间阈值，则继续计算下一时刻的模型数值解，并判断计算过程中待输出的状态变量是否发生数值振荡，若发生数值振荡则采用滤波器进行仿真结果波形平滑；若不小于预设的时间阈值，则完成电磁暂态仿真，输出电磁暂态数值仿真结果。
[0013]进一步地，对所述非线性元件进行建模的过程包括：
[0014]根据非线性元件的特性曲线，通过分段线性法、函数拟合法或分段线性法与函数拟合法相结合的方法进行曲线拟合，得到对应的描述函数。
[0015]进一步地，所述分段线性法与函数拟合法相结合的方法具体为：
[0016]对于所述特性曲线中的分段线性段和拐弯段采用函数拟合法进行拟合，对于所述特性曲线中的饱和段采用线性拟合法进行拟合。
[0017]进一步地，所述统一形式的数学模型的表达式为：
[0018][0019]式中，y为状态变量，为状态变量的一阶导数，H为常系数矩阵，r(y,t)为非线性广义项，t为时间。
[0020]进一步地，在任意积分区间[t
k
,t
k+1
]内，所述数学模型转化为的同解积分方程的表达式为：
[0021][0022]式中，y
k+1
为k+1时刻的状态变量值，e为自然常数，h为t
k
与t
k+1
之间的积分步长，y
k
为k时刻的状态变量值，ξ为积分变量；
[0023]所述积分方程的齐次项为e
hH
y
k
；
[0024]所述积分方程的非齐次项为
[0025]进一步地，令：
[0026][0027]对所述积分方程的非齐次积分项通过Calahan算法的计算表达式为：
[0028][0029][0030]y
k+1
＝y
k
+0.75k1+0.25k2[0031]式中，k1和k2均为中间量，y
k
为t
k
时刻的Calahan算法计算结果，y
k+1
为t
k+1
时刻的Calahan算法计算结果，γ为算法系数，α
21
为算法系数，J为函数f的Jacobian矩阵，为时间t
k
和状态变量y
k
的一维函数，为F(t
k
+α
21
h)的函数关系t
k
+α
21
h和状态变量y
k
+α
21
k1的一维函数。
[0032]进一步地，所述Calahan算法结合Richardson外推法的计算过程包括：
[0033]对于一个p阶的方法，把以步长h计算到t
k+1
点的数值解记为y
k+1
(h)，以h/2为步长计算到相同点t
k+1
的数值解记为y
k+1
(h/2)，微分方程在t
k+1
点的精确解为y
k+1
，则有：
[0034][0035]以ε(t
k+1
)为主误差函数，则：
[0036][0037][0038]式中，y
k+1
作为新的数值解。
[0039]进一步地，所述滤波器为一阶数字低通滤波器。
[0040]进一步地，所述一阶数字本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于实时积分算法的非线性元件电磁暂态仿真方法，其特征在于，包括以下步骤：对电力网络中的非线性元件进行建模，并结合电力网络原始参数，建立电力系统m维非线性常微分方程组；将所述电力系统m维非线性常微分方程组改写为齐次项加非齐次项组合的形式，得到统一形式的数学模型；对所述数学模型的电磁暂态数值进行初始化；将所述数学模型转化为同解的积分方程，对所述积分方程的齐次项采用积分算法计算下一时刻的值，对所述积分方程的非齐次积分项通过Calahan算法计算下一时刻值，并结合Richardson外推法提高Calahan算法的精度，最终得到下一时刻的模型数值解；若当前模型数值解所对应的时刻小于预设的时间阈值，则继续计算下一时刻的模型数值解，并判断计算过程中待输出的状态变量是否发生数值振荡，若发生数值振荡则采用滤波器进行仿真结果波形平滑；若不小于预设的时间阈值，则完成电磁暂态仿真，输出电磁暂态数值仿真结果。2.根据权利要求1所述的一种基于实时积分算法的非线性元件电磁暂态仿真方法，其特征在于，对所述非线性元件进行建模的过程包括：根据非线性元件的特性曲线，通过分段线性法、函数拟合法或分段线性法与函数拟合法相结合的方法进行曲线拟合，得到对应的描述函数。3.根据权利要求2所述的一种基于实时积分算法的非线性元件电磁暂态仿真方法，其特征在于，所述分段线性法与函数拟合法相结合的方法具体为：对于所述特性曲线中的分段线性段和拐弯段采用函数拟合法进行拟合，对于所述特性曲线中的饱和段采用线性拟合法进行拟合。4.根据权利要求1所述的一种基于实时积分算法的非线性元件电磁暂态仿真方法，其特征在于，所述统一形式的数学模型的表达式为：式中，y为状态变量，为状态变量的一阶导数，H为常系数矩阵，r(y,t)为非线性广义项，t为时间。5.根据权利要求4所述的一种基于实时积分算法的非线性元件电磁暂态仿真方法，其特征在于，在任意积分区间[t
k
,t
k+1
]内，所述数学模型转化为的同解积分方程的表达式为：式中，y
k+1
为k+1时刻的状态变量值，e为自然常数，h为t
k
与t
k+1
之间的积分步长，y
k
为k时刻的状态变量值，ξ为积分变量；所述积分方程的齐次项为e
hH
y
k
；所述积分方程的非齐次项为6.根据权利要求5所述的一种基于实时积分算法的非线性元件电磁...

【专利技术属性】
技术研发人员：王永，叶婧，谢继豪，张磊，吴秀海，许杨，吕金伟，桂传林，马伟，王少杰，黎慧，崔律，
申请(专利权)人：三峡大学，
类型：发明
国别省市：