多维阶数控制的多步Taylor级数暂态稳定分析方法与流程

本发明是多维阶数控制的多步taylor级数暂态稳定分析方法，属于电力系统暂态稳定性分析领域。

背景技术：

暂态稳定分析一直是电力系统领域的重要内容。在暂态稳定分析中，最基本的工作是求解表征动态元件的高维微分代数方程组，寻找一种高效、数值稳定的微分代数方程求解方法至关重要。在数学解法层次上，快速的高阶taylor级数法作为典型，其利用了状态变量的高阶导数而允许采用较大的积分步长，提高了计算速度；显式多步taylor级数法，用多步低阶导数值代替了单步高阶导数值，避免求解高阶导数，节省了计算量。

现有技术中对于taylor级数法暂态稳定分析而言，主要的计算量都集中在网络状态变量高阶导数的递推求解上。通常对于一个特定的系统，taylor级数法每一步的计算量是导数最高阶数p的幂函数，随着阶数p的增大，计算量将显著增加。常规的多步taylor级数法对整个积分过程的每一状态变量都采用统一固定的展开阶数，而采用统一固定阶数就会使得仿真过程有时计算冗余，有时精度达不到要求。

技术实现要素：

本发明的目的在于根据计算精度，对时间常数不同的各个发电机、同一发电机的不同时刻，计算不同阶数的状态变量导数，实现时间和空间上的动态多维阶数控制，消除固定阶数taylor级数法暂态稳定分析的计算冗余。从而能够有效的实现动态多维阶数控制，避免冗余计算量从而提高计算效率。

为实现上述目的，本发明采用了以下技术方案：一种多维阶数控制的多步taylor级数暂态稳定计算方法，该方法包括以下步骤：

(1)确定系统误差精度确定误差精度ε和展开步长h，设置精度阶数pi的初值为1；

(2)依次计算第i台发电机的pi+1阶误差常数和功角的pi+1阶导数

(3)分别判断第i台发电机精度阶数为pi时的截断误差是否满足精度要求；

(4)对于满足精度要求的发电机，此时对应的导数阶数即为状态变量δi应该采用的最高阶数，忽略其更高阶导数；

(5)对于没有达到精度条件的发电机功角，令pi＝pi+1，可以继续向其高阶导数递推，重复步骤(2)～(5)。

本发明的推导过程：

电力系统的暂态稳定分析需求解如下形式的微分方程初值问题，其中η为待求变量的初值，f(y(t),t)为所求微分方程的导数。

多步taylor级数法求解处置问题的计算通式为：当且时，上述多步高阶暂态稳定计算方法称为k步s阶taylor级数法，其中h代表积分步长，y^(j)n-i代表状态变量第i步j阶导数，αij代表第i步j阶导数的系数。该方法用前k步的导数值可方便地求得k+1步的状态变量数值解。

发电机采用暂态电抗x'd后的暂态电势e'恒定的经典模型，忽略励磁系统、原动机及调速器动态，网络线性，恒定阻抗，设负荷阻抗已经归入节点导纳矩阵，并将网络节点导纳矩阵扩展到包含发电机内节点。则基于降阶后的发电机的微分方程为cij＝eiejbij,dij＝eiejgij形式。其中δi和ωi分别表示第i台机组的转角和转速，pei代表机组i的电磁功率，pmi为机组i的机械功率，mi为机组i的转动惯量，ei和ej分别为发电机i和j的电压，gij和bij为收缩后导纳矩阵相应元素的实部和虚部，peij＝cijsinδij+dijcosδij，qeij＝dijsinδij-cijcosδij分别为第i台发电机和第j台发电机之间的有功功率和无功功率。

根据以上推导可以建立发电机功角导数从低阶向高阶的递推计算体系

式中：代表第i台机组转角在tk的j阶导数。

若计算精度阶数为p阶，则k步s阶taylor级数公式的截断误差为

η(p+1)＝cp+1y^(p+1)(t)h^p+1+o(h^p+2)。其中，o(h^p+2)为步长h的p+2次迭代结果的无穷小项；

其中s代表步数，αis代表第i步s阶导数的系数。

根据多步taylor级数式，求解暂态稳定问题的taylor级数展开式为式中：h是状态变量第k步的积分步长，pi是状态变量i的展开阶数。

对于给定误差精度ε和一定的展开步长h，截断误差满足下式时，此时δi对应的pi便是当前步状态变量i应该采用的展开阶数。

根据上述截断误差的判别式来确定第s台发电机功角最高阶数为m，即前述的导数最高阶数pi＝m，忽略第s台发电机功角m阶以上导数，由[0015]中功角导数公式可知不会影响到其他发电机的m+3阶以下功角导数。

而对其他发电机的m+3阶功角导数的影响，以发电机i为例，忽略第s台发电机功角m阶以上导数，不会对多项式p产生影响，只会对m+3阶导数产生影响，具体分析如下。

采用多步taylor级数法动态多维阶数控制后，忽略第s台发电机功角的m阶以上导数，第i台发电机功角计算与m+3阶导数相关的展开项为：

其中qeis＝dissinδis-ciscosδis为第i台发电机和第s台发电机之间的无功功率。

由式|η(pi+1)|＜ε可知令则第i台发电机的m+3阶功角导数展开式可近似写成

在实际计算中接近10^-2，因此e是一个非常小的数，一般情况下其绝对值在10^-4以下。可见eεi对δi的计算精度影响很小，完全可以忽略。也就是说，忽略第s台发电机功角的m阶以上导数，多步taylor级数法对其他发电机功角的m+2阶以上导数的影响很小，不会影响计算精度。

与现有技术相比，本发明具有以下优点：

由以上过程看出，对于阶数大于pi的各项，因为它们在整个展开式中所占比重甚小，故可把阶数大于pi的各项导数直接取作零，停止高阶导数的计算。对于没有达到精度条件的发电机功角，可以继续向其高阶导数递推，直至达到精度条件，这样便得到各个发电机功角的展开阶数pi。

最高阶导数pi的确定是在获得发电机转角导数后再结合误差常数依据截断误差判别式进行直接判断，并不增加关于pi的额外复杂计算，因此计算工作量小，实用性强。

通过分析常微分方程组的数学解法，根据计算精度对不同积分时刻的状态变量、同一时刻不同状态变量采用不同的taylor展开阶数，实现暂态稳定分析状态变量时间和空间上的动态多维阶数控制，消除固定阶数的显式多步taylor级数法暂态稳定分析的计算冗余。

附图说明

为了更清楚地说明本发明所提出算法的有效性和实用性，下面采用newengland10机39节点系统为例，假定故障为17-27线路靠近17节点三相短路，0.1秒时切除故障，仿真时间取为2秒。

下面将对本案例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例仿真结果图，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明提供的多维阶数控制的多步taylor级数暂态稳定计算方法实施的流程图；

图2是采用单步12阶taylor级数法得到的转角曲线图；

图3是采取本发明提出的计算方法，动态多维变阶taylor级数法的仿真转角曲线图；

图4是对图1、图2中的摆幅最大的34号机(图中加点标记曲线)取0.1秒为间隔进行逐点对比34号机转角曲线对比结果图。

图5是newengland系统中，采用动态多维阶数控制时，每一时间点发电机的展开阶数对比图。

具体实施方式

下面将结合案例中的附图，对本发明所提出的计算方法进行清楚、完整地描述。

如图1所示，该图为本发明提供的分层调节电网储能的方法实施流程图。该方法包括以下步骤：

(1)确定系统误差精度确定误差精度ε和展开步长h，设置精度阶数pi的初值为1；

(2)利用公式依次计算第i台发电机的pi+1阶误差常数

(3)利用公式依次计算第i台发电机的功角导数

(4)分别判断第i台发电机精度阶数为pi时的截断误差是否满足精度要求

(5)对于满足精度要求的发电机，此时对应的导数阶数即为状态变量δi应该采用的最高阶数，忽略其更高阶导数；

(6)对于没有达到精度条件的发电机功角，令pi＝pi+1，可以继续向其高阶导数递推，重复步骤(2)～(6)。

图2是采用单步12阶taylor级数法得到的转角曲线图；图3是采取本发明提出的计算方法，动态多维变阶taylor级数法的仿真转角曲线图。

对比图2、图3可以看出本文的动态多维阶数控制方法可以很好进行暂态稳定分析，满足工程精度10度的要求。

为了便于比较结果，对图2、图3中的摆幅最大的34号机(图中加点标记曲线)取0.1秒为间隔进行逐点对比，对比结果如图4。

从图4看出，采用统一变阶多步taylor级数法和动态多维阶数控制方法均能得到满足工程精度要求的转角曲线。但动态多维阶数控制方法表现更为出色，保证计算精度的同时，有效地降低了计算量。与统一变阶taylor级数法相比，提高了计算效率(0.039s-0.035s)/0.035s＝10％左右。

newengland系统中，采用动态多维阶数控制时，每一时间点发电机的展开阶数如图5所示。

计算时仿真步长取为0.05秒，表中仅给出了间隔0.1秒时发电机的展开阶数结果。

从图5看出，采用易于变阶的显式多步taylor级数法时，每一积分时间点的所有发电机转角导数取为相同阶数；采用动态多维阶数控制方法时，在某些时间点避免了统一阶数控制方法中部分高阶导数的求取，各台发电机可以采用不同的展开阶数，从而节省计算量。