为解决系统小扰动情况下的同步发电机参数辨识问题,应用改进遗传算法进行计算。针对轮盘赌选择操作的随机性,通过改变选择方式,有效地规避了搜索后期易陷入局部最优的缺陷。以PMU实测数据为基础,选择改进Euler法求解微分方程,根据最小二乘原理构建目标函数,分别采用改进的遗传算法和粒子群算法进行辨识,比较了两种算法的优劣,辨识结果表明改进的遗传算法具有更高辨识精度。 A genetic algorithm for the identification of synchronous generator parameters is proposed in the case of system minor disturbance. In view of the randomness of the operation of roulette, the selection method is improved, and the problem of local optimum is avoided effectively. In order to solve the differential equation, the improved Euler method is selected, and the objective function is built according to the application of least square principle, then the improved genetic algorithm and particle swarm optimization (PSO) used for the identification of synchronous generator parameters are identified which based on PMU measured data. The identification results which compare the ad-vantages and disadvantages of two kinds of algorithm show that the improved genetic algorithm has higher identification accuracy.
朱文豪
青岛大学自动化与电气工程学院,山东 青岛
收稿日期:2018年3月25日;录用日期:2018年4月11日;发布日期:2018年4月18日
为解决系统小扰动情况下的同步发电机参数辨识问题,应用改进遗传算法进行计算。针对轮盘赌选择操作的随机性,通过改变选择方式,有效地规避了搜索后期易陷入局部最优的缺陷。以PMU实测数据为基础,选择改进Euler法求解微分方程,根据最小二乘原理构建目标函数,分别采用改进的遗传算法和粒子群算法进行辨识,比较了两种算法的优劣,辨识结果表明改进的遗传算法具有更高辨识精度。
关键词 :同步发电机,六阶模型,参数辨识,改进Euler法,遗传算法
Copyright © 2018 by author and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
随着电网的规模和智能化水平的不断提高,电力系统的安全维护愈加成为广大电力从事者所关注的焦点。同步发电机是电力系统的核心元件,若同步发电机参数偏离实际,将严重影响电力系统的稳定性,因而,同步发电机的分析、控制、保护、运行均对其参数的准确性有着很高的要求。进而言之,同步发电机参数的准确辨识对电力系统的稳定运行分析有着极其重要的意义 [
同步发电机参数辨识的方法主要分为频域响应法和时域响应法 [
同步发电机参数的辨识需要首先建立其数学模型,选取合适的电压测量值作为输入,通过改进的遗传算法来寻找一组参数,使得经过模型计算得到的输出电流量和实测电流量的误差最小,此时对应的参数即为辨识出的参数。因此,选取合适的同步发电机数学模型是同步发电机参数辨识的重要前提。
同步发电机的实用参数数学模型在不同场合下有不同的选取原则,本文选用利于描述实心转子的六阶模型作为同步发电机的数学模型,此与实际运行发电机相一致。现对同步发电机六阶模型作出如下假设 [
1) 忽略定子磁暂态过程;
2) 计及转子d轴、q轴的暂态过程及次暂态过程。
根据上述假设,建立了以标幺值(时间常数除外)表示的同步发电机电气模型:
d轴模型
{ T ′ d 0 ⋅ d E ′ q d t = E f − E ′ q − x d − x ′ d x ′ d − x ″ d ( E ′ q − E ″ q ) T ″ d 0 ⋅ d E ″ q d t = E ′ q − E ″ q − ( x ′ d − x ″ d ) i d + T ″ d 0 ⋅ d E ′ q d t u q = E ″ q − x ″ d i d E f = k ⋅ u f d (1)
q轴模型
{ T ′ q 0 ⋅ d E ′ d d t = − E ′ d − x q − x ′ q x ′ q − x ″ q ( E ′ d − E ″ d ) T ″ q 0 ⋅ d E ″ d d t = E ′ d − E ″ d + ( x ′ q − x ″ q ) i q + T ″ q 0 ⋅ d E ′ d d t u d = E ″ d + x ″ q i q (2)
式(1)、(2)中参数以d轴为例进行说明,其中k为电机内电势Ef与励磁电压ufd的比例系数;ud和id为d轴电压和电流;xd为d轴同步电抗; x ′ d 、 x ″ d 、 T ′ d0 、 T ″ d0 、 E ′ d 和 E ″ d 是d轴暂态电抗、次暂态电抗、暂态开路时间常数、次暂态开路时间常数、暂态电势和次暂态电势。
由同步发电机的电气模型可知,电机的d轴、q轴解耦,所以可将d轴和q轴参数分开辨识。辨识时先计算稳态参数,再辨识暂态参数,流程见图1。
在稳态参数计算时,由d轴电压Ud和q轴电流Iq可计算q轴的同步电抗,即
x q = U d / I q (3)
根据隐极式的同步发电机特点,d轴同步电抗与q轴同步电抗近似相等,可取
x d = x q (4)
在参数辨识过程中需要寻优,这就先要构造一个目标函数。由(1)、(2)可知,根据 u q = E ″ q − x ″ d i d 和 u d = E ″ d − x ″ q i q 可以得到d轴和q轴电流的设定值,通过辨识计算可得到d轴和q轴电流的计算值,应用最小二乘原理以电流实际值和辨识值的差值加权最小作为目标函数,d轴、q轴的目标函数建立如下:
J d = ∑ t = 1 n ( i dc − i dm ) 2 (5)
J q = ∑ t = 1 n ( i qc − i qm ) 2 (6)
其中,Jd和Jq分别为d轴和q轴模型的目标函数,idc和iqc分别为根据d轴、q轴的辨识参数计算出的电流值,idm和iqm是根据PMU数据通过Park变换处理后的d轴、q轴电流实测值。
本文在用遗传算法寻优时,并不采用上述的目标函数,而是对目标函数进行变换。这里采用了遗传
算法的适应度函数Fd和Fq,现以d轴为例说明:目标函数Jd之值越接近于零,则辨识效果越好,而适应度函数Fd的构建则以选择较大的适应度值为原则,以便于参数寻优。因此,参数辨识中为了将辨识原理与优化方法联系起来,可取目标函数倒数当作适应度函数,即Fd = 1/Jd (Jd不为0)。
以PMU数据处理后的实测d轴电压ud、q轴电压uq、励磁电压ufd作为输入量,求解d轴电气模型(1)、q轴电气模型(2)中的微分方程中的4个状态变量,最终得到d轴电流id、q轴电流iq。
解方程时选用具有四阶精度的改进Euler法,先利用Euler公式求得一个初步的近似值 y ¯ n + 1 ,称之为预测值,预测值 y ¯ n + 1 的精度相对较差,再用梯形公式将其校正一次,得到 y n + 1 ,这个结果称为校正值,这样建立的预测——校正方法通常称为改进的Euler法,改进Euler法的示意图如图2所示。
y ¯ n + 1 = y n + h f ( x n , y n ) (7)
y n + 1 = y n + h 2 [ f ( x n , y n ) + f ( x n + 1 , y ¯ n + 1 ) ] (8)
遗传算法(GA)是一种借鉴生物界自然选择可实现全局并行优化的搜索工具 [
遗传算法不同于传统的进化类法,主要步骤包括:选择、交叉和变异,其中交叉和变异是遗传算法的主要部分,即待寻优参数通过按照一定的交叉概率和变异概率产生新个体。本文设置算法的交叉概率为0.6,设置变异概率为0.02 [
图1. 同步发电机参数辨识流程
图2. 改进Euler法
稳态数据计算得出稳态参数,然后将稳态参数代入暂态辨识过程。暂态辨识中待辨识的暂态参数有: x ′ d ( p .u ) 、 x ″ d ( p .u ) 、 T ′ d 0 ( s ) 和 T ″ d 0 ( s ) ,辨识过程以d轴电压ud和励磁电压ufd作为模型的输入值,以d轴电流id为输出值。
考虑到待辨识参数的计算量问题,决定对待辨识参数的初值进行浮点数编码。为了克服算法在搜索后期效率低和易于收敛的缺点,因此本文对该算法的选择算子进行了改进:放弃轮盘赌选择方法,先对待辨识参数初值在有意义的辨识区域内进行随机设置,先计算个体的适应度,然后对个体适应度的大小进行排序,以平均适应度为阈值,选择适应度值大于平均适应度值的个体,并以最高适应度值为模板,选择不同模板的个体组成群体,直至到达种群规模。寻优时搜索到的每一组参数,都要代入d轴电气模型(1)中,并使用改进的Euler法求解d轴电流的辨识值,通过设置合理的迭代次数,即可得到最终的最优值。程序验证时发现:迭代次数过大并不会使辨识结果的精度更高,反而会增加辨识的计算时间,所以迭代次数设置在300上下即可得到为理想的辨识效果。改进遗传算法的流程图如图3所示。
以PMU在线记录的小扰动数据为基础,使用同步发电机六阶数学模型,选用具有四阶精度的改进Euler法求解电机微分方程,应用最小二乘原理构建模型的目标函数。之后分别采用粒子群算法(PSO)和改进的遗传算法(AGA)进行同步发电机参数辨识,以比较粒子群算法和改进遗传算法的优劣性。
本文基于华北电网PMU所测得的实际数据。选取内蒙古托克托电厂一号机组2010年4月14日上午10点47分的实际数据,我们截取扰动发生的42秒至46秒这个时间段,机端正序电压、正序电流和励磁电压的波形如图3所示。
之后分别采用粒子群算法和改进的遗传算法对PMU实测数据进行辨识验证,得到的辨识结果如下:其中表1和表2为两种优化算法的d轴、q轴参数辨识的结果比较。图4~图7为使用两种优化算法进行
图3. 改进遗传算法流程简示
图4. 实测数据波形
图5. d轴电流拟合曲线图–粒子群算法
图6. d轴电流拟合曲线图–改进遗传算法
图7. q轴电流拟合曲线图–粒子群算法
参数辨识得到的d轴、q轴的电流拟合曲线。
从图5~图8的拟合曲线的拟合效果可以看出,改进的遗传算法对电流曲线的拟合效果相比于粒子群
图8. q轴电流拟合曲线图–改进遗传算法
参数 | 目标值 | Xd'(p.u) | Xd"(p.u) | Td0'(s) | Td0"(s) |
---|---|---|---|---|---|
设定值 | 0.220 | 0.190 | 8.200 | 0.037 | |
PSO | 4.87 × 10−4 | 0.299 | 0.1230 | 8.934 | 0.071 |
AGA | 5.28 × 10−5 | 0.257 | 0.1320 | 8.234 | 0.041 |
PSO误差% | 35.90 | −35.26 | 8.95 | 91.89 | |
AGA误差% | 16.8 | −30.5 | 0.41 | 10.81 |
表1. 辨识结果比较-d轴
参数 | 目标值 | Xq'(p.u) | Xq"(p.u) | Tq0'(s) | Tq0 "(s) |
---|---|---|---|---|---|
设定值 | 0.380 | 0.190 | 1.000 | 0.045 | |
PSO | 5.24 × 10−4 | 0.297 | 0.252 | 1.382 | 0.137 |
AGA | 3.24 × 10−5 | 0.351 | 0.202 | 1.192 | 0.075 |
PSO误差% | −21.84 | 32.63 | 38.20 | 204.04 | |
AGA误差% | −7.63 | 6.31 | 19.20 | 66.67 |
表2. 辨识结果比较-q轴
算法要更加理想;表1和表2中的参数辨识结果表明,改进遗传算法由于全局寻优能力强,且规避了陷入局部最优的问题,使得搜寻出的同步发电机的参数与实际值更加接近,辨识误差更小;且计算得到的目标函数在数量级上远小于粒子群算法,辨识过程的精度更高。
本文以提出更好的同步发电机参数辨识的优化方法为目的,对于传统的遗传算法,针对选择操作的随机性,放弃了轮盘赌选择方法,并改进了选择方法,有效地规避了算法在搜索后期易陷于局部最优的问题。以PMU量测数据为基础,分别用粒子群算法和改进的遗传算法进行同步发电机参数辨识,并对两种优化算法的结果进行对比,参数辨识结果显示了改进的遗传算法更为有效且具更高的精度。改进的遗传算法不仅可以应用于电力系统参数辨识领域,同样可应用于其他的优化问题。
朱文豪. 基于改进遗传算法的同步发电机参数辨识 Parameter Identification of Synchronous Generator Based on Improved Genetic Algorithm[J]. 智能电网, 2018, 08(02): 135-142. https://doi.org/10.12677/SG.2018.82015