电力系统数字仿真中的功率计算问题-天津大学研究生e-Learning平台.DOC

电力系统数字仿真中的功率计算问题 李霞林 （天津大学电气与自动化工程学院，1011203002） 1 引言 数字仿真作为电力系统分析中一个重要研究方向，其用途十分广泛，包括：安全性稳定性分析，小扰动分析，电力系统谐波分析，继电保护装置整定和校验等等。数字仿真的实质是求解微分代数方程组的时域解，目前商业软件PSCAD采用EMTP法，MATLAB则通过选择不同的常微分方程组解法来计算。两种方法的计算结果是电压和电流的离散时域量，在应用于实际电力系统问题过程中，往往还需要求解出电路各个节点流过的功率。通常来说，电力系统是基于50Hz的正弦波，但是在实际仿真中由于元件自身的电磁过程、电网频率的微小便宜以及仿真本身的误差，如何快速准确地求解功率成为了一个实际工程问题。 2 电磁仿真基本原理 2.1 基于EMTP原理的数值计算方法 EMTP方法是上世纪70年代由DOMMEL教授提出的一种电磁仿真方法[1]。这种方法具有较高的计算速度和计算精度，因此广泛应用于电力系统的商业软件中，其基本原理以电阻电感串联电路为例进行说明。 图1 电阻电感串联电路 Fig. 1 Series circuit of R-L 图1中所示的是最简单的电路元件，由一个电阻和电感串联而成，其数学方程如公式（1）所示。 （1） （2） 其中： EMTP法即对公式（1）进行隐式梯形化差分，形成公式（2），并将此数学模型等值成为一个电阻和电流源的电路模型如图1所示，从而将时域的电路化简成为在每一个仿真时刻的直流电路模型进行求解。 2.2 电路功率求解原理 电磁仿真的结果是电压和电流的瞬时时域量，通常情况下，电力系统分析需要知道流过电路的功率大小。实际电网为交流电路，交流电路的功率由有功功率P和无功功率Q组成。实际工程实践主要关心的是基波（50hz）的电路功率值，由此产生出几种不同的功率求解方法，本文即通过电磁仿真的电路结果对功率采用不同的算法进行求解。 3 基于时域频域变换思想的功率求解方法 根据电路原理[2]，有功功率和无功功率的表达式为公式（3）所示： (3) 其中：U，I均为电压和电流在复平面下的幅值，为电压和电流的夹角。为了求解功率，需要得出幅频域下电压和电流的数值。因此，需要对电磁仿真的数值解进行时域—频域变换，本文采用了三种不同的变换方式进行求解，从中选出最合适的一种方法。 3.1 基于快速傅里叶变换（FFT）的求解方法 离散傅里叶变换的思想是将一组离散的数字量变换到频域，从而对其进行频谱分析，其基础是基于三角函数形式的傅里叶级数形式。当给定离散点集（电磁仿真的离散时间数据）Xj已知时，可以求出函数f(x)的一组三角差值函数： (4) 电力系统所发出的电压电流信号是正弦波，因此，通过傅里叶变换，能够很好地分析出电网信号的频域特性。其中，a0对应电磁仿真中电压电流相量的直流分量，a1 、b1对应电磁仿真中电压电流相量的基频分量（即基频向量值），a2 、b2对应电磁仿真中电压电流相量的二倍基频分量。其他高次谐波则可以忽略不计。 本文采用开源软件fftw[3]作为仿真的算法包嵌入到已有的电磁仿真软件中，通过求解电压电流信号，获取功率值。 3.2 基于QR分解的最小二乘法的求解方法 按照给定公式（4）求解离散量的频域值还可以利用最小二乘法(lsm)进行求解。并且假定电磁仿真的信号中不存在3次以上的谐波分量，从而将公式（4）进一步化简为公式（5）： （5） 对公式（5）进行变形，形成公式（6）： （6） 其中： 对应最小二乘法的方程组如公式（7）所示： （7） 为了保证仿真过程中，最小二乘法的稳定性，需要避免构造正规方程组，在这种情况下，本文选用基于反施密特正交化方法（RGS）的QR分解法进行求解： 对于F可以分解为：F=QR，其中Q为m*n;R为n*n的上三角矩阵。 采用RGS(反施密特正交化的方法进行求解)，其中： 具体算法流程则如下所示：： 输入电磁仿真数据X矩阵，建立F矩阵； 依次采用算法计算R矩阵（二维矩阵）和Q矩阵（按列存储为n行）； 按照QR分解法的方法计算C阵； 反求电压、电流的幅值相角。 3.3 基于改进最小二乘法（CFA）的求解方法 采用QR分解法的最小二乘法需要求解出五维的未知量，相对计算量较大，在考虑到计算速度的场合下，上述方法显得不够快速。因此，本文采用了更为简化的求解方法，即改进最小二乘法（CFA）[4]进行求解。其基本思想是认为电磁仿真的信号仅仅由基波组成，即进一步忽略直流分量和高次谐波。从而公式（5）化简成为公式（8）： （8） 对于公式（8）仍然采用最小二乘曲线拟合的方法，但是由于方程仅仅为二阶方程组，其可以不需要采用QR分解