基于多变量拟合的发电机调速系统参数优化建模方法
基于多变量拟合的发电机调速系统参数优化建模方法
【技术领域】
[0001] 本发明属于电力系统分析与控制技术领域，具体设及一种基于多变量拟合的发电 机调速系统参数优化建模方法。
【背景技术】
[0002] -次调频性能的研究重点在于控制系统的瞬态性能指标，其应达到超调量尽可能 小，上升时间和调节时间尽可能快，控制系统稳、准、快的要求。一次调频性能指标受调节系 统参数的影响，但目前并没有一次调频性能指标与控制参数之间的定量关系。传统调节系 统参数的选择通常根据经验，缺乏具体的理论计算依据，往往得不到最优控制参数。
【发明内容】
[0003] 本发明的目的在于克服已有技术存在的不足，提供一种基于多变量拟合的发电机 调速系统参数优化建模方法。
[0004] 本发明的目的通过W下技术方案实现。
[0005] 基于多变量拟合的发电机调速系统参数优化建模方法，包括W下步骤：
[0006] (1)系统超调量由四个线性控制参数，单位调节功率、油动机时间常数、汽轮机蒸 汽时间常数和发电机惯性时间常数共同决定，上升时间只与发电机惯性时间常数有关，调 节时间与汽轮机蒸汽时间常数和发电机惯性时间常数有关，故将系统超调量01%、系统第 二个周波超调量〇2 %及调节时间ts进行多变量拟合，对上升时间tr进行单变量拟合。
[0007] 多变量拟合的表达式形式需根据单变量拟合的结果来确定，故首先假设：目标因 变量T因单变量X变化时的多项式拟合形式为式(1)如下：
[000引
[0009]式中，T为超调量、上升时间或调节时间；X为单位调节功率K、油动机时间常数Ts、 汽轮机蒸汽时间常数Tch或发电机惯性时间常数化。
[0010]根据单变量拟合的结果，假设T与K、Ts、Tch、Th拟合出的多项式的次数分别为i，j， k, 1，指数变化范围分别为m = 0, i ;n = 0, j ;p = 0，k;q = 0,1则多变量拟合因变量Τ如（2)所 示：
[0011]
[0012]将上述方法运用到调速系统性能指标的多变量拟合中去，于是对于系统超调量 01%的拟合过程就有:首先进行单变量拟合，得其随1(、了3、1'。4、1'4均成线性关系变化，根据此 结果，可知式(2)中i = j = k=l = l，于是可得如式(3)所示的〇1%的多变量拟合式；
[0016]同理可得系统第二个周波超调量02%、调节时间ts、上升时间tr的拟合表达式如式 (4) -(6)，其中，02%随1(、了3、1'。4、化均成线性关系变化，*3随了3成线性关系变化，随化成平方 关系变化，tr随调速系统线性参数化变化。
[0022] 其中，K为单位调节功率，Ts为油动机时间常数，Tc功汽轮机蒸汽时间常数，Th为发 电机惯性时间常数，aia = l，2,...，n)为每一项的系数，通过仿真获得大量离散数据，将运 些数据带入各拟合表达式中，将所得到的方程联立，就可W求得ai的值；；
[0023] (2)调速系统稳定运行时，给定调速系统控制参数取值范围及其他约束条件，建立 调节系统参数的不等式约束方程：
[0024]
[002引其中，8<01/02。0为对调速系统中衰减率的限制；50| ΔΡ(1+0ι%化I <0.2为对 系统中频率最大波动值50 ± 0.2Ηζ的限制条件；
[0026] (3)在步骤（1)和(2)的基础上，建立W系统一次调频的性能指标最为优化的目标 函数，即式(8)，W下两式表示要使超调量、上升时间和调节时间取到最小的值，使系统的频 率响应更为优化，
[0027]
[0028] (4)基于W上条件，应用目标多约束非线性优化方法进行寻优求解过程：
[0029] 1)取初始惩罚因子'^>〇，允许误差6>〇;
[0030] 2)选取迭代初值，取值如式(9)所示；
[0031] x=(〇i,tr,ts,K,Ts,Tch,Th) = (0.25,20,23,0.11,0.2,7,0.05) (9)
[0032] 3)构造惩罚函数
如式（10)所示，
[0033]
[0034] 其中，rW是正数序列，逐个递减，称惩罚因子：
産罚因子的取值为: =1.0,0.1，0.01，0.001，1(^(°)>'(1)>'口)>1>卢)>卢")>1>0，邑。口)为式(7)中的不等 式约束条件，各惩罚函数为并列关系；
[0035] 4)从点出发用无约束优化方法求解惩罚函数的的极值点
[0036] 5)设定相对误差默认值为10-6作为检查迭代是否终止的准则，若误差在默认值范 围内，则停止迭代计算，并为原目标函数f(X)的约束最优解，否则迭代次数加1， 转入步骤3)，继续求解。
[0037] 进一步，上述步骤建立了一次调频性能指标与控制参数的数量关系，在此基础上 W调频性能指标最优为优化目标，并满足相应的不等式约束条件，通过多目标多约束非线 性优化方法，将一次调频指标与控制参数之间关系量化，将调节系统参数优化问题转化为 求解式(8)最小值问题，即可得到最优的调速系统参数。
[0038] 上述数学模型是一个非线性规划模型，通过多目标多约束非线性优化方法可求得 其最优解，既满足有约束的条件，又能使目标函数取得最小值。
[0039] 与现有技术相比，本发明具有W下优点和有益效果：
[0040] 本发明所述方法W超调量、上升时间、调节时间作为一次调频的性能评价指标，考 虑调速系统四个线性参数(单位调节功率、油动机时间常数、汽轮机蒸汽时间常数或发电机 惯性时间常数)对指标的影响，采用多变量拟合的方法，利用仿真数据拟合控制参数和性能 评价指标之间的定量关系，得出它们之间的表达式。在此基础上，通过插值拟合建立一次调 频指标与控制参数之间的定量关系，将调节系统参数的优化问题转化为求解非线性规划的 最优解问题，改善了传统靠经验确定控制参数不理想的问题。通过本申请方法可得到一次 调频性能最优的一组调速系统控制参数，提高调速器一次调频性能，使超调量尽可能小，上 升时间和调节时间尽可能快，达到控制系统稳、准、快的要求。
【附图说明】
[0041] 图1是某地区大规模电力外送系统结构图。
[0042] 图2是超调量指标的单变量拟合图。
[0043] 图3是系统第二个周波超调量02%的单变量拟合图。
[0044] 图4是上升时间指标的单变量拟合图。
[0045] 图5是调节时间的单变量拟合图。
[0046] 图6是系统发生单相接地后的发电机转速偏差曲线图；
[0047] 图7是系统发生两相短路后的发电机转速偏差曲线图；
[0048] 图8是系统切除国能一台机后的其他发电机转速偏差曲线图。
【具体实施方式】
[0049] 下面结合附图对本发明进一步详细说明。本发明提出的基于多变量拟合的发电机 调速系统参数优化建模方法说明如下：
[0050] 某地区大规模电力外送系统结构如图1所示。±800kV特高压直流配套火电电源共 6600MW，其中神华国能电厂4X660MW、国投电厂、国电电厂和瑞虹哈密电厂均为2X660MW。 10台火电机组通过500kV线路接入新哈换500kV交流变电站，新哈换500kV交流母线连接± 800kV特高压直流整流侧换流变压器，连接750/500kV联络变压器，通过750/500kV联络变压 器将500kV送端电网接入750kV地区主网。系统的运行方式为直流孤岛运行方式，即切除该 地区送端交流网络，仅投直流配套电源。
[0051] (1)依据电厂试验报告所提供的参数和多变量拟合方法得到拟合出的系统超调 量、系统第二个周波超调量、调节时间和上升时间的多变量拟合式。
[0052] 系统超调量σι%随调速 系统线性参数1(^3、1'。4、化变化，对系统超调量进行单变量 拟合，拟合结果如图2所示。
[0053] 由单变量拟合的结果可见:系统超调量01%随1(、了3、1'。4、化均成线性关系变化，根据 此结果，写出01%的多变量拟合式如式(3)。
[0054] 通过仿真获得大量离散数据，将运些数据带入各拟合表达式中，将所得到的方程 联立，就可W求得各系数日1的值。
[0055] σι %的多变量拟合表达式的各系数如表1所示。
[0056] 表1超调量的多变量拟合表达式系数 [0化7]
[0059] 系统第二个周波超调量02%的单变量拟合图像如图3所示，由拟合结果可见，系统 第二个周波超调量02%也随1(、了3、1'。4、扣均成线性关系变化，根据系统超调量〇1%的拟合式 书写过程，可W得到系统第二个周波超调量02%的拟合式如式(4)。同理，我们可W通过仿 真得到大量离散数据求得拟合式中各系数的值，系数值如表2所示。
[0060] 表2第二周期超调量的多变量拟合表达式系数
[0061]
[0062] 上升时间tr随调速系统线性参数化变化。对上升时间进行单变量拟合，拟合结果如 图4所示。可见上升时间tr随调速系统线性参数化呈线性关系变化。
[0063] 则上升时间tr的拟合表达式如式(5)所示。将仿真所得离散数据代入式(5)，联立 方程求解各系数，得拟合出的上升时间拟合式系数如表3所示。
[0064] 表3上升时间的多变量拟合表达式系数
[00 化]
[0066] 调节时间ts随调速系统线性参数Ts、Th变化。利用MTLAB对调节时间进行多变量拟 合。图5为ts随Ts成线性关系变化，随化成平方关系变化。
[0067] 根据ts的单变量拟合结果，得到调节时间ts的拟合表达式为式(6)所示。将仿真所 得离散数据代入式(6)，联立方程求解各系数，得拟合出的上升时间拟合式系数如表4所示。
[0068] 表4调节时间的多变量拟合表达式系数
[0069]
[0070] (2)得到一次调频调速系统参数寻优目标函数规划，其中等式约束为步骤(1)中得 到的一次调频的多变量拟合式，不等式约束为控制参数取值范围及其他条件；
[0071] (3)步骤（1)和(2)的基础上，建立W系统一次调频的性能指标最为优化的目标函 数如下式：
[0072]
[0073] (4)基于步骤(2)的约束及(3)的目标函数，应用多目标多约束非线性优化方法进 行寻优求解过程：
[0074] 1)取初始惩罚因子rW>〇,允许误差ε>〇;
[0075] 2)选取迭代初值，取值如式(9)所示；
[0076] x=(〇i,tr,ts,K,Ts,Tch,Th) = (0.25,20,23,0.11,0.2,7,0.05) (9)
[0077] 3)构造惩罚函数口 (义，卢1)如式（10)所示；
[007引
[0079] 其中，是正数序列，逐个递减，称惩罚因子，
惩罚因子的取值为: = 1.0,0.1，0.01，0.001，K，rW>r(l)>r(リ>L>r(k)>r(k")>L>0，gu(X)为式(7)中的不等 式约束条件，各惩罚函数为并列关系；
[0080] 4)从xfk-il点出发用无约束优化方法求解惩罚函数的极值点Χ*0-?);
[0081] 5)设定相对误差默认值为10-6作为检查迭代是否终止的准则，若误差在默认值范 围内，则停止迭代计算，并为原目标函数f(X)的约束最优解，否则迭代次数加1， 转入步骤3)，继续求解。求解满足要求后，得到的控制参数取值结果如表5所示：
[0082] 表5内点法寻优的参数取值结果
[0083]
[0084] (5)将得到的最优控制参数及系统内电厂采用的控制参数分别用于该系统的调速 系统模型进行仿真，得到系统在直流孤岛运行情况下的仿真结果对比如图6-图8所示。
[0085] 由仿真结果来看，在直流孤岛运行方式下，系统中发生单相接地短路、两相相间短 路及切除国能一台机时，寻优所得参数的优势体现为频响超调量小、稳定时间快，一次调频 的性能优势体现明显。
【主权项】
1.基于多变量拟合的发电机调速系统参数优化建模方法，其特征在于，包括以下步骤： (1) 将系统超调量、系统第二个周波超调量σ2%及调节时间ts进行多变量拟合，对 上升时间tr进行单变量拟合，建立一次调频性能指标与控制参数之间的数量关系式（1)~ (4)：其中，K为单位调节功率，Ts为油动机时间常数，Td1为汽轮机蒸汽时间常数，Th为发电机 惯性时间常数，ai(i = l，2,. . .，n)为每一项的系数，通过大量离散实验数插值拟合求得； (2) 调速系统稳定运行时，给定调速系统控制参数取值范围及其他约束条件，建立调节 系统参数的不等式约束方程：其中，^ 10为对调速系统中衰减率的限制；50 I ΔΡ(1+σι%)Κ I <0.2为对系统 中频率最大波动值50 ± 0.2Hz的限制条件； (3) 在步骤（1)和(2)的基础上，建立以系统一次调频的性能指标最为优化的目标函数， 即式(6)，(4) 应用目标多约束非线性优化方法迭代求解函数最优解，具体算法如下： 1) 取初始惩罚因子r(())>0,允许误差ε>0; 2) 选取迭代初值，取值如式(7)所示： X = (〇i,tr,ts,K, Ts,Tch, Th) = (0.25,20,23,0.11,0.2,7,0.05) (7) 3) ?? ? I8)^,其中，r(k)是正数序列，逐个递减，称惩罚因子，惩罚因子的取值为:r(k) = 1.0,0.1,0.01,0.001，K，r(0)>r(1)>r(2)>L>r( k)>r(k+1)>L>0，gu(X)为式(5)中的不等式 约束条件，各惩罚函数为并列关系； 4) 从点出发用无约束优化方法求解惩罚函数...................|的极值点浐(r(k));5) 设定相对误差默认值为HT6作为检查迭代是否终止的准则，若误差在默认值范围内， 则停止迭代计算，并以浐(r (k))为原目标函数f (X)的约束最优解，否则迭代次数加1，转入步 骤3)，继续求解。2.根据权利1所述的基于多变量拟合的发电机调速系统参数优化建模方法，其特征在 于，建立一次调频性能指标与控制参数的数量关系，在此基础上以调频性能指标最优为优 化目标，并满足相应的不等式约束条件，通过多目标多约束非线性优化方法，将一次调频指 标与控制参数之间关系量化，将调节系统参数优化问题转化为求解式(6)最小值问题，即可 得到最优的调速系统参数。
【专利摘要】本发明公开了一种基于多变量拟合的发电机调速系统参数优化建模方法。该方法首先采用多变量拟合方法，拟合发电机一次调频三个性能指标超调量、上升时间、调节时间与调速系统四个控制参数单位调节功率、油动机时间常数、汽轮机蒸汽时间常数和发电机惯性时间常数之间的关系，然后将得到的一次调频的多变量拟合式作为等式约束，将控制参数取值范围及其他条件作为不等式约束，将精细化调速系统达到最优一次调频性能作为最终目标，采用多目标多约束非线性优化方法寻得使一次调频性能最优的调速系统参数，建立优化后的发电机调速系统模型。
【IPC分类】H02P9/00, H02J3/00, H02P101/20
【公开号】CN105490266
【申请号】CN201510988585
【发明人】智勇, 秦睿, 梁福波, 拜润卿, 郑伟, 陈仕彬, 肖仕武, 刘军
【申请人】国网甘肃省电力公司电力科学研究院, 国网甘肃省电力公司, 国家电网公司, 华北电力大学
【公开日】2016年4月13日
【申请日】2015年12月24日
【技术领域】
[0001] 本发明属于电力系统分析与控制技术领域，具体设及一种基于多变量拟合的发电 机调速系统参数优化建模方法。
【背景技术】
[0002] -次调频性能的研究重点在于控制系统的瞬态性能指标，其应达到超调量尽可能 小，上升时间和调节时间尽可能快，控制系统稳、准、快的要求。一次调频性能指标受调节系 统参数的影响，但目前并没有一次调频性能指标与控制参数之间的定量关系。传统调节系 统参数的选择通常根据经验，缺乏具体的理论计算依据，往往得不到最优控制参数。
【发明内容】
[0003] 本发明的目的在于克服已有技术存在的不足，提供一种基于多变量拟合的发电机 调速系统参数优化建模方法。
[0004] 本发明的目的通过W下技术方案实现。
[0005] 基于多变量拟合的发电机调速系统参数优化建模方法，包括W下步骤：
[0006] (1)系统超调量由四个线性控制参数，单位调节功率、油动机时间常数、汽轮机蒸 汽时间常数和发电机惯性时间常数共同决定，上升时间只与发电机惯性时间常数有关，调 节时间与汽轮机蒸汽时间常数和发电机惯性时间常数有关，故将系统超调量01%、系统第 二个周波超调量〇2 %及调节时间ts进行多变量拟合，对上升时间tr进行单变量拟合。
[0007] 多变量拟合的表达式形式需根据单变量拟合的结果来确定，故首先假设：目标因 变量T因单变量X变化时的多项式拟合形式为式(1)如下：
[000引
[0009]式中，T为超调量、上升时间或调节时间；X为单位调节功率K、油动机时间常数Ts、 汽轮机蒸汽时间常数Tch或发电机惯性时间常数化。
[0010]根据单变量拟合的结果，假设T与K、Ts、Tch、Th拟合出的多项式的次数分别为i，j， k, 1，指数变化范围分别为m = 0, i ;n = 0, j ;p = 0，k;q = 0,1则多变量拟合因变量Τ如（2)所 示：
[0011]
[0012]将上述方法运用到调速系统性能指标的多变量拟合中去，于是对于系统超调量 01%的拟合过程就有:首先进行单变量拟合，得其随1(、了3、1'。4、1'4均成线性关系变化，根据此 结果，可知式(2)中i = j = k=l = l，于是可得如式(3)所示的〇1%的多变量拟合式；
[0016]同理可得系统第二个周波超调量02%、调节时间ts、上升时间tr的拟合表达式如式 (4) -(6)，其中，02%随1(、了3、1'。4、化均成线性关系变化，*3随了3成线性关系变化，随化成平方 关系变化，tr随调速系统线性参数化变化。
[0022] 其中，K为单位调节功率，Ts为油动机时间常数，Tc功汽轮机蒸汽时间常数，Th为发 电机惯性时间常数，aia = l，2,...，n)为每一项的系数，通过仿真获得大量离散数据，将运 些数据带入各拟合表达式中，将所得到的方程联立，就可W求得ai的值；；
[0023] (2)调速系统稳定运行时，给定调速系统控制参数取值范围及其他约束条件，建立 调节系统参数的不等式约束方程：
[0024]
[002引其中，8<01/02。0为对调速系统中衰减率的限制；50| ΔΡ(1+0ι%化I <0.2为对 系统中频率最大波动值50 ± 0.2Ηζ的限制条件；
[0026] (3)在步骤（1)和(2)的基础上，建立W系统一次调频的性能指标最为优化的目标 函数，即式(8)，W下两式表示要使超调量、上升时间和调节时间取到最小的值，使系统的频 率响应更为优化，
[0027]
[0028] (4)基于W上条件，应用目标多约束非线性优化方法进行寻优求解过程：
[0029] 1)取初始惩罚因子'^>〇，允许误差6>〇;
[0030] 2)选取迭代初值，取值如式(9)所示；
[0031] x=(〇i,tr,ts,K,Ts,Tch,Th) = (0.25,20,23,0.11,0.2,7,0.05) (9)
[0032] 3)构造惩罚函数
如式（10)所示，
[0033]
[0034] 其中，rW是正数序列，逐个递减，称惩罚因子：
産罚因子的取值为: =1.0,0.1，0.01，0.001，1(^(°)>'(1)>'口)>1>卢)>卢")>1>0，邑。口)为式(7)中的不等 式约束条件，各惩罚函数为并列关系；
[0035] 4)从点出发用无约束优化方法求解惩罚函数的的极值点
[0036] 5)设定相对误差默认值为10-6作为检查迭代是否终止的准则，若误差在默认值范 围内，则停止迭代计算，并为原目标函数f(X)的约束最优解，否则迭代次数加1， 转入步骤3)，继续求解。
[0037] 进一步，上述步骤建立了一次调频性能指标与控制参数的数量关系，在此基础上 W调频性能指标最优为优化目标，并满足相应的不等式约束条件，通过多目标多约束非线 性优化方法，将一次调频指标与控制参数之间关系量化，将调节系统参数优化问题转化为 求解式(8)最小值问题，即可得到最优的调速系统参数。
[0038] 上述数学模型是一个非线性规划模型，通过多目标多约束非线性优化方法可求得 其最优解，既满足有约束的条件，又能使目标函数取得最小值。
[0039] 与现有技术相比，本发明具有W下优点和有益效果：
[0040] 本发明所述方法W超调量、上升时间、调节时间作为一次调频的性能评价指标，考 虑调速系统四个线性参数(单位调节功率、油动机时间常数、汽轮机蒸汽时间常数或发电机 惯性时间常数)对指标的影响，采用多变量拟合的方法，利用仿真数据拟合控制参数和性能 评价指标之间的定量关系，得出它们之间的表达式。在此基础上，通过插值拟合建立一次调 频指标与控制参数之间的定量关系，将调节系统参数的优化问题转化为求解非线性规划的 最优解问题，改善了传统靠经验确定控制参数不理想的问题。通过本申请方法可得到一次 调频性能最优的一组调速系统控制参数，提高调速器一次调频性能，使超调量尽可能小，上 升时间和调节时间尽可能快，达到控制系统稳、准、快的要求。
【附图说明】
[0041] 图1是某地区大规模电力外送系统结构图。
[0042] 图2是超调量指标的单变量拟合图。
[0043] 图3是系统第二个周波超调量02%的单变量拟合图。
[0044] 图4是上升时间指标的单变量拟合图。
[0045] 图5是调节时间的单变量拟合图。
[0046] 图6是系统发生单相接地后的发电机转速偏差曲线图；
[0047] 图7是系统发生两相短路后的发电机转速偏差曲线图；
[0048] 图8是系统切除国能一台机后的其他发电机转速偏差曲线图。
【具体实施方式】
[0049] 下面结合附图对本发明进一步详细说明。本发明提出的基于多变量拟合的发电机 调速系统参数优化建模方法说明如下：
[0050] 某地区大规模电力外送系统结构如图1所示。±800kV特高压直流配套火电电源共 6600MW，其中神华国能电厂4X660MW、国投电厂、国电电厂和瑞虹哈密电厂均为2X660MW。 10台火电机组通过500kV线路接入新哈换500kV交流变电站，新哈换500kV交流母线连接± 800kV特高压直流整流侧换流变压器，连接750/500kV联络变压器，通过750/500kV联络变压 器将500kV送端电网接入750kV地区主网。系统的运行方式为直流孤岛运行方式，即切除该 地区送端交流网络，仅投直流配套电源。
[0051] (1)依据电厂试验报告所提供的参数和多变量拟合方法得到拟合出的系统超调 量、系统第二个周波超调量、调节时间和上升时间的多变量拟合式。
[0052] 系统超调量σι%随调速 系统线性参数1(^3、1'。4、化变化，对系统超调量进行单变量 拟合，拟合结果如图2所示。
[0053] 由单变量拟合的结果可见:系统超调量01%随1(、了3、1'。4、化均成线性关系变化，根据 此结果，写出01%的多变量拟合式如式(3)。
[0054] 通过仿真获得大量离散数据，将运些数据带入各拟合表达式中，将所得到的方程 联立，就可W求得各系数日1的值。
[0055] σι %的多变量拟合表达式的各系数如表1所示。
[0056] 表1超调量的多变量拟合表达式系数 [0化7]
[0059] 系统第二个周波超调量02%的单变量拟合图像如图3所示，由拟合结果可见，系统 第二个周波超调量02%也随1(、了3、1'。4、扣均成线性关系变化，根据系统超调量〇1%的拟合式 书写过程，可W得到系统第二个周波超调量02%的拟合式如式(4)。同理，我们可W通过仿 真得到大量离散数据求得拟合式中各系数的值，系数值如表2所示。
[0060] 表2第二周期超调量的多变量拟合表达式系数
[0061]
[0062] 上升时间tr随调速系统线性参数化变化。对上升时间进行单变量拟合，拟合结果如 图4所示。可见上升时间tr随调速系统线性参数化呈线性关系变化。
[0063] 则上升时间tr的拟合表达式如式(5)所示。将仿真所得离散数据代入式(5)，联立 方程求解各系数，得拟合出的上升时间拟合式系数如表3所示。
[0064] 表3上升时间的多变量拟合表达式系数
[00 化]
[0066] 调节时间ts随调速系统线性参数Ts、Th变化。利用MTLAB对调节时间进行多变量拟 合。图5为ts随Ts成线性关系变化，随化成平方关系变化。
[0067] 根据ts的单变量拟合结果，得到调节时间ts的拟合表达式为式(6)所示。将仿真所 得离散数据代入式(6)，联立方程求解各系数，得拟合出的上升时间拟合式系数如表4所示。
[0068] 表4调节时间的多变量拟合表达式系数
[0069]
[0070] (2)得到一次调频调速系统参数寻优目标函数规划，其中等式约束为步骤(1)中得 到的一次调频的多变量拟合式，不等式约束为控制参数取值范围及其他条件；
[0071] (3)步骤（1)和(2)的基础上，建立W系统一次调频的性能指标最为优化的目标函 数如下式：
[0072]
[0073] (4)基于步骤(2)的约束及(3)的目标函数，应用多目标多约束非线性优化方法进 行寻优求解过程：
[0074] 1)取初始惩罚因子rW>〇,允许误差ε>〇;
[0075] 2)选取迭代初值，取值如式(9)所示；
[0076] x=(〇i,tr,ts,K,Ts,Tch,Th) = (0.25,20,23,0.11,0.2,7,0.05) (9)
[0077] 3)构造惩罚函数口 (义，卢1)如式（10)所示；
[007引
[0079] 其中，是正数序列，逐个递减，称惩罚因子，
惩罚因子的取值为: = 1.0,0.1，0.01，0.001，K，rW>r(l)>r(リ>L>r(k)>r(k")>L>0，gu(X)为式(7)中的不等 式约束条件，各惩罚函数为并列关系；
[0080] 4)从xfk-il点出发用无约束优化方法求解惩罚函数的极值点Χ*0-?);
[0081] 5)设定相对误差默认值为10-6作为检查迭代是否终止的准则，若误差在默认值范 围内，则停止迭代计算，并为原目标函数f(X)的约束最优解，否则迭代次数加1， 转入步骤3)，继续求解。求解满足要求后，得到的控制参数取值结果如表5所示：
[0082] 表5内点法寻优的参数取值结果
[0083]
[0084] (5)将得到的最优控制参数及系统内电厂采用的控制参数分别用于该系统的调速 系统模型进行仿真，得到系统在直流孤岛运行情况下的仿真结果对比如图6-图8所示。
[0085] 由仿真结果来看，在直流孤岛运行方式下，系统中发生单相接地短路、两相相间短 路及切除国能一台机时，寻优所得参数的优势体现为频响超调量小、稳定时间快，一次调频 的性能优势体现明显。
【主权项】
1.基于多变量拟合的发电机调速系统参数优化建模方法，其特征在于，包括以下步骤： (1) 将系统超调量、系统第二个周波超调量σ2%及调节时间ts进行多变量拟合，对 上升时间tr进行单变量拟合，建立一次调频性能指标与控制参数之间的数量关系式（1)~ (4)：其中，K为单位调节功率，Ts为油动机时间常数，Td1为汽轮机蒸汽时间常数，Th为发电机 惯性时间常数，ai(i = l，2,. . .，n)为每一项的系数，通过大量离散实验数插值拟合求得； (2) 调速系统稳定运行时，给定调速系统控制参数取值范围及其他约束条件，建立调节 系统参数的不等式约束方程：其中，^ 10为对调速系统中衰减率的限制；50 I ΔΡ(1+σι%)Κ I <0.2为对系统 中频率最大波动值50 ± 0.2Hz的限制条件； (3) 在步骤（1)和(2)的基础上，建立以系统一次调频的性能指标最为优化的目标函数， 即式(6)，(4) 应用目标多约束非线性优化方法迭代求解函数最优解，具体算法如下： 1) 取初始惩罚因子r(())>0,允许误差ε>0; 2) 选取迭代初值，取值如式(7)所示： X = (〇i,tr,ts,K, Ts,Tch, Th) = (0.25,20,23,0.11,0.2,7,0.05) (7) 3) ?? ? I8)^,其中，r(k)是正数序列，逐个递减，称惩罚因子，惩罚因子的取值为:r(k) = 1.0,0.1,0.01,0.001，K，r(0)>r(1)>r(2)>L>r( k)>r(k+1)>L>0，gu(X)为式(5)中的不等式 约束条件，各惩罚函数为并列关系； 4) 从点出发用无约束优化方法求解惩罚函数...................|的极值点浐(r(k));5) 设定相对误差默认值为HT6作为检查迭代是否终止的准则，若误差在默认值范围内， 则停止迭代计算，并以浐(r (k))为原目标函数f (X)的约束最优解，否则迭代次数加1，转入步 骤3)，继续求解。2.根据权利1所述的基于多变量拟合的发电机调速系统参数优化建模方法，其特征在 于，建立一次调频性能指标与控制参数的数量关系，在此基础上以调频性能指标最优为优 化目标，并满足相应的不等式约束条件，通过多目标多约束非线性优化方法，将一次调频指 标与控制参数之间关系量化，将调节系统参数优化问题转化为求解式(6)最小值问题，即可 得到最优的调速系统参数。
【专利摘要】本发明公开了一种基于多变量拟合的发电机调速系统参数优化建模方法。该方法首先采用多变量拟合方法，拟合发电机一次调频三个性能指标超调量、上升时间、调节时间与调速系统四个控制参数单位调节功率、油动机时间常数、汽轮机蒸汽时间常数和发电机惯性时间常数之间的关系，然后将得到的一次调频的多变量拟合式作为等式约束，将控制参数取值范围及其他条件作为不等式约束，将精细化调速系统达到最优一次调频性能作为最终目标，采用多目标多约束非线性优化方法寻得使一次调频性能最优的调速系统参数，建立优化后的发电机调速系统模型。
【IPC分类】H02P9/00, H02J3/00, H02P101/20
【公开号】CN105490266
【申请号】CN201510988585
【发明人】智勇, 秦睿, 梁福波, 拜润卿, 郑伟, 陈仕彬, 肖仕武, 刘军
【申请人】国网甘肃省电力公司电力科学研究院, 国网甘肃省电力公司, 国家电网公司, 华北电力大学
【公开日】2016年4月13日
【申请日】2015年12月24日
最新回复(0)