[0040] (3-4)根据当前采样周期中每一步预测得到的控制律增量、机组频率预测值与给定值的频率偏差,获得状态量预测序列。
[0041] 优选地,上述的抽水蓄能机组调速系统非线性预测模型,建立抽水蓄能机组调速系统的非线性预测模型的方法,包括如下子步骤:
[0042] (1-1)建立执行机构模型;包括描述接力器工作速率限制的饱和环节1,描述主接力器积分环节和描述接力器运动行程限制的饱和环节2的反馈闭环结构;
[0043] 该执行机构模型考虑了各饱和限幅环节和主接力器时间常数,能精确地描述抽蓄机组控制系统执行机构的非线性运动特性;
[0044] (1-2)采用考虑引水管道弹性水击效应的弹性水击模型,建立引水系统模型,其函数关系为
[0045] 其中,Tr为水击相长,hw为管路特性系数,q(s)是指流量偏差相对值,h(s)指水头偏差相对值,s是指拉普拉斯算子;
[0046] (1-3)根据电站实测的水泵水轮机全特性曲线数据建立水泵水轮机的非线性插值模型;具体为转矩-转速特性曲线Mt11=f1(α,n11)和流量-转速特性曲线Q11=f2(α,n11);
[0047] 其中,Mt11为单位转矩,Q11为单位流量,α为导叶开度,n11为单位转速;
[0048] 为了克服水泵水轮机的“驼峰”特性与“S”区域对插值的影响,对全特性曲线采用改进Suter变换处理以消除其插值的多值性问题,可充分描述水泵水轮机复杂的工况变化特性;
[0049] (1-4)建立发电机模型,其函数关系为
[0050] 其中,Ta为抽蓄机组的惯性常数;en为水轮发电机组自调节系数;ω为转速偏差相对值;mt为水泵水轮机机械转矩偏差相对值;mg0为负载转矩变动偏差相对值;
[0051] 发电机模型为采用反应转子运动特性与机组自调节能力的一阶发电机模型,在保证机组模型转速求解精度前提下,简化发电机电磁特性,提高模型计算效率;
[0052] (1-5)建立抽水蓄能机组调速系统的非线性预测模型,包括执行机构模型、引水系统模型、水泵水轮机模型和发电机模型;;
[0053] 其中,模型输入为机组即时控制律u,输出为机组实时转速w;模型状态变量包括导叶开度y,流量q,水压h,机械转矩mt和转速w。
[0054] 总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
[0055] (1)本发明提供的抽水蓄能机组调速系统的快速非线性模糊预测控制方法,充分考虑调速系统各组成部分的复杂非线性特性,构建了抽水蓄能机组非线性预测模型;
[0056] 由于模型中直接运用抽水蓄能机组实测全特性曲线并对其进行数学变换处理作为水泵水轮机非线性模型,较传统采用状态微分方程、六参数模型等局部线性化模型的方法更能全面反映其复杂非线性特性,且更能适应控制过程中机组工况的大范围变化;而液压执行机构模型重点考虑主接力器动作速率、行程上的饱和限制,而忽略掉配压阀,辅助接力器,油系统等对调速控制过程影响微弱的部件,简化模型阶次;对与水泵水轮机同轴的发电/电动机模型,忽略比调速暂态过程动作快很多的电磁暂态,将其简化为一根刚性转动的轴,既保留了机组转速变化的主要特性又大大减小模型计算量;由此,本发明使用的调速系统非线性模型为抽蓄机组调速系统预测控制的实施建立了模型基础;
[0057] (2)本发明提供的抽水蓄能机组调速系统的快速非线性模糊预测控制方法,借助模糊控制进行在线滚动预测,该方法在获取控制律增量预测序列的过程中,由于控制律的求解过程借助模糊PID控制规律,在预测时域内用滚动预测代替了滚动优化,避免了传统预测控制算法中复杂的非线性规划问题的在线求解,因此较传统预测控制的滚动优化方法极大程度地减少了在线计算量,缩短了在线计算时间,保证了预测控制的实时性;
[0058] (3)本发明提供的抽水蓄能机组调速系统的快速非线性模糊预测控制方法,区别于传统预测控制中只保留预测控制律序列的第一项(即时控制律)而抛弃其余时刻的预测控制律的做法,本发明通过对模糊PID控制下预测时域内各预测采样时刻的机组频率偏差预测序列以非线性递减形式进行加权处理得到累积预测频率偏差,并将其作为PID控制器的输入来求取预测控制器的即时控制律,从而使机组即时控制律的设置同时包含机组的实时状态偏差信息和对未来控制过程的预测信息;由于控制器实时控制律中已经考虑了未来系统可能出现的状态量偏差及其后续控制规律,故控制偏差较传统控制方法能够得到进一步减小;
[0059] (4)本发明提供的抽水蓄能机组调速系统的快速非线性模糊预测控制方法,由于在每个采样周期,模糊推理机都会根据实时工况信息中的转速偏差及其偏差变化率,分别通过三个控制参数的模糊推理原则,并在预测时域内考虑到PID参数的变化速率和系统状态偏差变化趋势间的关系,不断对实时PID参数进行调整,达到控制器参数工况自适应的效果,进而运用PID控制根据滚动预测原理预测机组未来预测采样时刻状态量信息,从而提高控制品质。
附图说明
[0060] 图1抽蓄机组调速控制系统非线性模糊预测控制流程;
[0061] 图2是实施例中抽水蓄能机组调速系统的结构框图;
[0062] 图3是实施例中抽水蓄能机组调速系统的执行机构非线性模型;
[0063] 图4是实施例中抽水蓄能机组调速系统的水泵水轮机全特性曲线;其中,(a)是转矩特性曲线,(b)是流量特性曲线;
[0064] 图5是实施例中抽水蓄能机组调速系统的改进Suter变换后的水泵水轮机全特性曲线;其中,(a)是WM曲线,(b)是WH曲线;
[0065] 图6是实施例中抽水蓄能机组调速系统非线性模糊预测控制框图;
[0066] 图7是实施例中模糊PID参数自适应原理图;
[0067] 图8是实施例中快速滚动预测流程图;
[0068] 图9是实施例中本发明控制方法下的抽蓄技能机组发电方向开机过程仿真结果示意图。
具体实施方式
[0069] 为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
[0070] 本发明实施例以某抽水蓄能电站机组的开机过程控制为例,对本发明提供的抽水蓄能机组调速系统的快速非线性模糊预测控制方法进行阐述。
[0071] 实施例中的抽水蓄能电站机组的设计参数如表1所列;
[0072] 表1抽水蓄能机组设计参数
[0073]转轮直径(m) 额定转速(r/min) 额定水头(m) 额定流量(m3/s) 额定机械转矩Nm
3.850 500 540 62.09 5.8446×106
[0074] 结合抽水蓄能机组实际运行经验,设定实施例的抽水蓄能电站机组的开机过程仿真整体控制策略:开机过程初期,在机组转速低于90%额定转速时,水轮机导叶以一定速度直线开启,至空载限制开度后保持不变;待机组转速到达90%额定转速时,机组转入实施例提供的预测控制;其中,预测控制器及仿真重要参数如表2所列;
[0075] 表2控制器及仿真参数
[0076]Kp0 Ki0 Kd0 N 投入频率 仿真时间(s) 开机时刻(s)
1.2 0.5 0.1 5 90%额定 60 5
[0077] 在开机后机组频率达到90%额定频率(即45Hz)后,机组转为由预测控制控制水轮机导叶开度的变化;抽水蓄能机组调速系统的非线性模糊预测控制的流程如图1所示,具体包括如下步骤:
[0078] 步骤1:建立抽水蓄能机组调速系统非线性预测模型;实施例中,抽水蓄能机组调速系统结构如图2所示,除控制器外,系统由执行机构、引水系统、水泵水轮机和发电/电动机四个部分组成;建立的抽水蓄能机组调速系统非线性预测模型,包括执行机构模型、弹性水击模型、水泵水轮机全特性曲线插值模型和一阶发电模型;其中,执行机构模型考虑各饱和限幅环节和主接力器时间常数,其结构如图3所示,饱和环节1描述接力器动作速率的限制,饱和环节2描述接力器运动行程的限制,Ty表示主接力器时间常数。
[0079] 引水系统模型采用考虑引水管道弹性水击效应的弹性水击模型,其函数式如下[0080]
[0081] 其中,Tr为水击相长,hw为管路特性系数;
[0082] 为了充分描述水泵水轮机复杂的工况变化特性,实施例中,抽水蓄能机组采用电站实测的水泵水轮机全特性曲线数据作为其非线性插值模型,如图4所示,图4(a)为不同典型开度下机组的单位转速-单位转矩关系曲线;图4(b)为不同典型开度下机组的单位转速-单位流量关系曲线。
[0083] 机组全特性曲线包含了机组的所有可能工况,在任意工况下通过确定的导叶开度和机组转速可以通过插值从图中分别得到转矩和流量。然而由于水泵水轮机存在“驼峰”效应和反“S”特性,故曲线在某些区域存在交叉和插值的多样性问题。
[0084] 为了克服水泵水轮机的“驼峰”特性与“S”区域对插值的影响,对全特性曲线采用改进Suter变换处理以消除其插值的多值性问题,如图5所示;其中,图5(a)为使用改进Suter变换处理后的WM曲线,图5(b)为使用改进Suter变换处理后的WH曲线。从图中可以看出,数学变换后的机组全特性曲线中的曲线交叠被完全解开,且不同开度下的曲线分布均匀,利于插值处理。
[0085] 发电机采用反应转子运动特性与机组自调节能力的一阶发电机模型,其函数式为[0086] 其中,Ta为抽蓄机组的惯性常数;en为水轮发电机组自调节系数;ω为转速偏差相对值;mt为水泵水轮机机械转矩偏差相对值;mg0为负载转矩变动偏差相对值;
[0087] 实施例中的抽水蓄能机组调速系统非线性模糊预测控制结构如图6所示,整个调速预测控制系统由预测控制器,执行机构,水泵水轮机及其引水系统和发电电动机四个部分组成。其中,模型输入为机组即时控制律u,输出为机组实时转速w;模型状态变量包括导叶开度y,流量q,水压h,机械转矩mt和转速w。
[0088] 实施例中,对采用函数表示的数学模型进行离散化,便于计算机实现。
[0089] 步骤2:根据初始状态对调速系统状态变量设置初值,包括水头偏差相对值h0、流量偏差相对值q0、导叶开度偏差相对值y0、机械转矩偏差相对值mt0和机组频率偏差相对值f0;
[0090] 根据初始运行工况设置PID参数初值Kp0,Ki0,Kd0,并设定预测控制器采样步长Ts与预测步数Np;实施例中,控制步数Nc与预测步数Np取值相等;
[0091] 步骤3:根据经验值与现场状态给定ΔKp,ΔKi,ΔKd三个控制参数增量的模糊隶属度函数和模糊规则矩阵,并根据机组频率偏差及其增量对其进行模糊PID参数自适应,如图7所示;图中,频率偏差和频率偏差增量信号通过模糊推理机,经过模糊化,模糊规则推理和去模糊三个步骤后级联自适应速率变化因子环节输出PID控制器的三个控制参数增量ΔKp,ΔKi和ΔKd。
[0092] 为了控制器实现的方便性与运算的快速性,上述三个控制参数的模糊隶属度函数均取为三角形函数,频率偏差e和偏差增量ec的模糊子集均为{NB,NM,NS,ZO,PS,PM,PB}七个状态;
[0093] 实施例中,ΔKp,ΔKi,ΔKd三个参数的模糊规则分别如表3至表5所列;
[0094] 表3ΔKp模糊规则表
[0095]
[0096] 表4ΔKi模糊规则表
[0097]
[0098] 表5ΔKd模糊规则表
[0099]
[0100] 步骤4:设置离散时间的机组频率给定值序列fref。
[0101] 步骤5:对控制系统状态进行在线滚动预测,其具体实施流程如图8所示;具体地,根据实时机组频率偏差由模糊PID控制器求得调速器控制律,通过快速滚动预测中的预测模型和机组实时状态,采用离散的预测算法可以求得下一时刻的状态量预测值;
[0102] 用某一时刻机组频率给定fref(k)与机组频率f(k)相减得到当前时刻频率偏差e(k),并计算频率偏差增量ec(k)=e(k)-e(k-1);
[0103] 以e(k)和ec(k)作为模糊推理输入,通过控制器模糊规则矩阵,分别得到当前时刻PID参数增量ΔKp(k),ΔKi(k),ΔKd(k);
[0104] 为避免PID参数过快变化造成的控制过程中出现不期望的波形抖动,根据系统存储的频率偏差变化趋势,设置PID参数速率变化因子λ(k);
[0105] 其中,
[0106] 将其与PID参数初值相加得到当前时刻自适应的PID参数
[0107] Kp(k),Ki(k),Kd(k);
[0108] 其中,
[0109] 步骤6:获取预测时域内各预测采样时刻预测控制律增量,具体包含以下子步骤:
[0110] (6-1)令
[0111] 获取对应控制律预测值u(k/k)=u(k-1)+Δu(k/k);
[0112] 其中,即时PID控制律增量输出
[0113] Δu(k/k)=[KP(k)+Ki(k)+Kd(k)]·e(k)+[-KP(k)-2Kd(k)]·e(k-1)+Kd(k)·e(k-2);
[0114] (6-2)根据抽蓄机组调速系统离散状态方程组、即时PID控制律u(k/k)与当前系统状态X(k),获得系统在k时刻对(k+1)时刻的预测状态X(k+1/k);
[0115] 其中,抽蓄机组调速系统离散状态方程组为X(k+1/k)=f(X(k),u(k/k))[0116] 其中,X(k)=[x1(k),x2(k),...,xn(k)]表示k时刻系统状态向量,[0117] Y(k)=[y1(k),y2(k),...,yn(k)]表示k时刻系统输出向量;
[0118] (6-3)将(k+1)时刻的状态预测量X(k+1/k)中机组频率预测值f(k+1/k)与其给定值fref(k+1)进行比较,获得k时刻对(k+1)时刻的频率偏差e(k+1/k)及频率偏差增量ec(k+1/k);
[0119] 其中,e(k+1/k)=fref(k+1)-f(k+1/k);ec(k+1/k)=e(k+1/k)-e(k/k);
[0120] 采用e(k+1/k)和ec(k+1/k)重复步骤5,获得k时刻对(k+1)时刻的预测PID自适应参数Kp(k+1/k),Ki(k+1/k),Kd(k+1/k);
[0121] 获取k时刻对(k+1)时刻的预测PID控制律增量Δu(k+1/k),
[0122]
[0123] 其中,对应的控制律u(k+1/k)=u(k/k)+Δu(k+1/k);
[0124] 利用X(k+1/k)和u(k+1/k),获得k时刻对(k+2)时刻的状态量预测值X(k+2/k)=f(X(k+1),u(k+1/k));
[0125] (6-4)自Δu(k+2/k)开始,求解k时刻的预测时域内的PID预测控制律增量Δu(k+i/k)不再与本时刻(k时刻)之前的机组频率偏差量相关;
[0126] 故预测时域内剩余预测采样时刻控制律增量预测值求解公式统一表示为:
[0127]
[0128] 其中,i=2,3,…,Np-1;(k+i)时刻的PID自适应参数预测值采用步骤(5)的方法根据e(k+i/k)以及ec(k+i/k)获取;
[0129] 对应的控制律u(k+i/k)=u(k+i-1/k)+Δu(k+i/k);
[0130] 根据X(k+i/k)和u(k+i/k)获取k时刻对(k+i+1)时刻的状态量预测值X(k+i+1/k)=f(X(k+i),u(k+i/k)),并从中计算对应时刻的机组频率与频率给定的偏差e(k+i+1/k)=f(k+i+1/k)-f(k+i/k);其中,i=2,3,…,Np-1;
[0131] (6-5)对k时刻控制器预测时域内预测采样时刻进行滚动预测完毕后,获取预测时域Np内每个预测采样时刻的机组频率的预测偏差值。
[0132] 步骤7:对远离当前时刻的控制器预测时域内未来时刻机组频率偏差以非线性递减权重进行加权得到预测时域内累计误差
[0133] 利用步骤5中PID参数,由
[0134] Δu(k)=[KP(k)+KI(k)+Kd(k)]·e∑(k)+[-KP(k)-2Kd(k)]·e∑(k-1)+Kd(k)·e∑(k-2)计算预测控制下的控制律增量。
[0135] 步骤8:获取当前时刻预测控制器的实际控制律取值
[0136] u(k)=u(k-1)+Δu(k)。
[0137] 实施例中,用本发明提供的预测控制方法对抽水蓄能机组调速系统从停机状态进行开机控制直到机组进入空载稳态的过程仿真,抽蓄机组各状态变量变化过程如图9所示,如图中所示的,机组由停机状态在5s开机启动,在机组转速达到0.9倍额定转速前,导叶开度以恒定速率逐渐开启;在转速达到0.9倍额定转速时,机组控制转入模糊预测控制,根据机组实时运行状态和控制算法预测得到的机组未来状态信息对导叶开度变化规律进行控制,从图中可以看出,导叶转速相对值从0逐步上升至1的过程平滑,基本无超调量,开机规律符合机组导叶开度给定曲线趋势且系统进入稳态运行时间短;流量、开度、水压等状态变量的暂态过程总体变化平稳,过程曲线符合正常开机规律。
[0138] 本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。