基于増维精细积分的高压输电线路电磁暂态数值计算方法
【技术领域】
[0001] 本发明涉及一种电力系统电磁暂态数值计算方法，具体是涉及一种基于増维精细 积分方法的高压输电线路电磁暂态数值计算方法。
【背景技术】
[0002] 高压输电线路是电力系统的重要元件之一。利用数值计算方法分析高压输电线路 的电磁暂态过程，可解决电力系统诸多领域的问题:可对输电线路故障进行仿真，得到短路 电流以及谐波分量等;可用于输电线路故障测距算法的研究以及输电线路空载合闸过电压 计算研究等。
[0003] 概括起来，高压输电线路的电磁暂态数值计算方法大致包括频域法和时域法这2 大类方法。频域类方法通常是先利用快速傅立叶变换或拉普拉斯变换求解输电线的频率响 应，然后利用相应的反变换以得到其时域解。频域类方法概念比较简单，但涉及的计算极为 繁琐，而且此类方法的求解过程主要是基于线性叠加原理，因而很难处理非线性问题。 [0004]时域法的基本思路是采用时间离散或空间离散方法将偏微分方程转化为常微分 方程，并在此基础上利用各类数值积分方法求解各状态变量的值。目前使用较多的时域类 方法主要包括:时域有限差分法、基于隐式梯形积分方法和隐式欧拉法相结合的临界阻尼 调整法（Critical Damping Adjustment，CDA)等。
[0005] 时域有限差分法的主要思想是把高压输电线路的电报方程进行时间、空间同时离 散，用差分方程代替一阶偏微分方程，再通过求解线性方程组以计算出各状态变量的值。时 域有限差分法算法简单，但该方法中时间步长的选取受稳定性条件的限制，当空间步长很 小时亦须选取很小的时间步长，从而导致计算效率低下。
[0006] 目前，商用化的EMTP(Electromagnetic Transient Program，EMTP)程序中主要是 采用临界阻尼调整法(CDA法)来求解高压输电线路的电磁暂态计算问题。CDA方法计算过程 简单，从理论上讲可以有效避免传统的电磁暂态计算方法中的数值振荡问题，但对一些无 法准确判断的突变现象，使用该方法进行计算仍将产生"虚假的"数值振荡现象。
[0007] 除上述数值计算方法以外，还有研究人员将精细积分方法应用到高压输电线路的 电磁暂态计算中。精细积分方法是钟万勰院士针对结构动力学方程的求解问题所提出的一 类新的数值计算方法，这种方法在求解齐次常微分方程时具有明显的优势：具有极高的计 算精度;计算过程简单、计算效率高。然而，对于高压输电线路，离散后的常微分方程组通常 是非齐次的，若直接采用精细积分法进行求解则将涉及矩阵求逆问题，而这种矩阵求逆运 算不仅计算量大，有时甚至是不稳定的。
【发明内容】
[0008] 本发明所要解决的技术问题，就是提供一种基于増维精细积分的高压输电线路电 磁暂态数值计算方法，其能根据高压输电线路的输入或激励是正弦信号这一特点，将空间 离散后所形成的非齐次常微分方程组通过适当增维转换成齐次常微分方程组，并采用精细 积分法来求解齐次常微分方程组，可以精确、高效地求解出高压输电线路的电磁暂态响应。
[0009] 解决上述技术问题，本发明采取如下的技术方案：
[0010] -种基于増维精细积分的高压输电线路电磁暂态数值计算方法，其特征在于包括 以下步骤:将高压输电线路的输入或激励正弦信号改写为状态变量，通过状态变量增维，将 描述高压输电线路电磁暂态计算的非齐次常微分方程转换为齐次形式的常微分方程，在此 基础上采用精细积分法对齐次常微分方程组进行数值积分计算。
[0011] 具体步骤包括：
[0012] 1)将高压输电线路进行空间离散化，建立相应的电磁暂态数值计算的数学模型， 即：
[0013] i = Ax + w(t) (1)；
[0014] 式(1)中，x是高压输电线路电磁暂态计算中的状态变量;A是与高压输电线路参数 [0015]相关的定常系数矩阵;w(t)即是电磁暂态计算的输入信号或激励信号，也就是所 谓的
[0016] 非齐次项，它是一个稀疏列向量；
[0017] 2)电磁暂态数值计算初始化
[0018] 置t = 〇.〇s，积分步数n = 0;
[0019] 确定数值积分步长h、电磁暂态数值计算时程T;
[0020] 确定各状态变量的初值，即x(t=0)=x〇;
[0021] 输入电磁暂态数值计算的故障或操作；
[0022] 3)故障或操作判断
[0023]依据时刻t判断系统有无故障或操作：
[0024]若无故障或操作，则直接转至步骤4);
[0025] 若有故障或操作，则修改系数矩阵A以及激励信号w(t)并重新形成方程(1);
[0026] 同时，若涉及状态变量的突变，则还需依据具体情况修改相应的状态变量值Xn (t)；
[0027] 4)常微分方程齐次化 [0028]依据时刻t判断w(t)是否为0:
[0029] 若w(t)=0,则无需齐次化，直接转至步骤5);
[0030] 若w(t)矣0,则利用状态变量増维方法将非齐次常微分方程即式（1)转换成以下形 式的齐次常微分方程：
[0031] i = m (2)：
[0032] 上式(2)中:Η是与原系数矩阵A等有关的一个新的常系数矩阵；
[0033]
(3)：；
[0034]上式(3)中，ΔΧ即是増维状态变量；
[0035] 5)数值积分
[0036] 采用经典的精细积分方法对齐次常微分方程(2)进行求解，计算出所有状态变量 的值；
[0037] 6)t = t+h；n = n+l ；
[0038] 7)数值计算是否终止判断
[0039] 若t<T，则返回步骤3):
[0040] 若t 2 T，则转至步骤8);
[0041] 8)数值计算结果输出。
[0042] 所述步骤4)常微分方程齐次化是本发明的核心部分，其具体的实施步骤如下： [0043]第一步:依据w(t)为正弦信号这一特点，引入新的变量并对状态变量进行増维
[0044] 具体过程可叙述如下：
[0045] 设非齐次项w(t) = [Wl(t)]T中含有Μ个非零元素，将这些不为零的元素依据其先后 顺序记为Wj(t)，je(l，M);
[0046] 由于高压输电线路电磁暂态数值计算的输入信号为正弦信号，因此，用以下通用 表达式来描述Wj (t):
[0047]
(4);
[0048] 上式(4)中：ω = 2πΧ50;α」为常量；g为相位角，是已知量。
[0049] 对应上述非零元素，引入如下变量：
[0054] (7>；
[0050] (5)；
[0051]
[0052] (6); 12345678910
[0053] 因此，对应w(t)中每一个非零元素，定义以下増维状态变量：
(8)； 2 第二步:非齐次常微分方程增维齐次化 3
[0057] 具体过程可叙述如下： 4
[0058]设原状态变量X的维数为s， 5 定义 6
[0060]
(9); 7 D = diag(dj)，je(l，M);DER2MX2M (10)； 8
[0062]
(11); 9 上式（11)中：若w(t)中第i个元素为零，贝ljei = 0;若w(t)中第i个元素不为零，则ei 是一个2M维的稀疏行向量，该向量的第2j_l个元素为1，其余元素全部为0; 10 基于上述定义，将非齐次方程即式(1)转换为形如式(2)的齐次方程，即：
[0065] (12)；
[0066]显然，Η是一个常系数矩阵；
[0067]第三步:计算增维状态变量的初值
[0068]
(13)；
[0069]第四步:置齐次化结果标志 [0070]完成上述齐次化后，置w(t)=0。
[0071 ]所述步骤5)数值积分中，利用精细积分法求解各状态变量的过程如下：
[0072]第一步:计算指数矩阵T
[0073] 利用精细积分法求解齐次常微分方程，需求解如下的指数矩阵：
[0074] T = eHh (14)；
[0075] 指数矩阵的求解过程叙述如下。
[0076] 将积分步长进一步细分为21^段，则每一小段的时间间隔为：
[0077] At = h/(2m) (15)；
[0078] 依据泰勒级数展开，首先计算如下初始矩阵：
[0079]
(16)；
[0080]上式（16)中:K为泰勒多项式的阶数;在上述基础上，循环求解：
[0081 ]
(17) ；
[0082]则式（14)中的指数矩阵为：
[0083] Γ~/ + ΓΛ(,,,) (1-8);
[0084] 第二步:利用指数矩阵Τ计算各状态变量在t = tn+1时刻的值
[0085] 利用上式(18)的结果，求解各状态变量在t = tn+1时刻的值：
[0086]
(19)〇
[0087] 上述精细积分方法具有极高的计算精度，因而能有效避免"虚假的"数值振荡现 象。此外，指数矩阵T只与Η和步长h有关:若系数矩阵Η和步长h不变，则只需一次性计算出指 数矩阵T即可;换言之，只有当发生故障或操作后才需要重新计算指数矩阵T。因此，上述精 细积分方法的计算效率高于传统的数值积分方法。
[0088] 本发明的重点就在于:利用高压输电线路的输入信号或激励信号为正弦信号这一 特点，通过状态变量的适当增维，将非齐次常微分方程组转化为齐次常微分方程组，并采用 精细积分法来求解齐次常微分方程组，可以精确、高效地完成高压输电线路的电磁暂态计 算。
[0089] 本发明的优点和积极效果是:基于増维精细积分方法的高压输电线路电磁暂态数 值计算方法，继承或保持了经典精细积分方法所具有的高精度这一优点，因而可以有效避 免传统电磁暂态数值计算方法所存在的数值振荡问题;与已有的电磁暂态数值积分方法相 比较，基于増维精细积分方法的高压输电线路电磁暂态数值计算方法具有更高的计算效 率。
【附图说明】
[0090] 图1:高压均匀长输电线路空载合闸示意图；
[0091] 图2:高压均匀长输电线路离散化等值电路示意图；
[0092] 图3:高压均匀长输电线路空载合闸的电磁暂态计算结果(初相角p =
[0093] 图4:高压均匀长输电线路空载合闸的电磁暂态计算结果(初相角p = 90 );
[0094] 图5:高压输电线路电磁暂态计算结果三维曲线图（初相角炉= 90°);
[0095] 图6:CDA方法计算高压输电线路空载合闸的电磁暂态结果(初相角，=0°);
[0096] 图7:CDA方法计算高压输电线路空载合闸的电磁暂态结果(初相角供= 90°);
[0097] 图8:相同步长下两种不同方法计算结果的对比曲线(初相角0 =
[0098] 图9:不同步长下两种不同方法计算结果的对比曲线(初相角炉=〇°);
【具体实施方式】
[0099]以下结合附图对本发明的实施方式进行概要说明。
[0100] 以如图1所示的高压工频均匀长输电线路在不同初相角情况下进行空载合闸的电 磁暂态计算为例，本发明的具体实施步骤可以概述如下：
[0101] 1)将高压输电线路进行空间离散化，建立相应的电磁暂态数值计算的数学模型
[0102] 输入输电线单位长度参数办、1^、0)，线路长度1^，电压源函数6(〇及其内阻1^
[0103] 空间离散化:取Ν=20,将L = 300km的输电线路均匀分为20段，表示为如图2所示的 空间离散化形式；离散后每一段线路的电阻、电感、电容分别为：
c = 。由此，可以建立各元件的微分方程如下：
[0104] (20)；
[0105] (21);
[0106] 上述方程共含2N+1个待求状态变量。将上述方程写成矩阵形式，即是：
[0107]
(22)；
[0108] 上式（22)中