发明内容
本发明为了克服以上技术的不足,提供了一种基于专家内模控制的串级回路PID控制器参数整定方法。
本发明的方法,无论是开环的数据还是闭环的数据,只要回路数据具有一定程度的激励,即可进行串级回路的PID控制器参数整定。且本发明适用于大部分串级回路,具有普适性。本发明采用基于专家内模控制的串级回路PID控制器参数整定方法,具体是在副回路采用基于专家经验的内模控制进行PI控制器参数的整定,在主回路采用基于专家经验的二自由度内模控制进行PID控制器参数的整定。通过该方法使得在不改变串级回路采用PID控制器的前提下,实现了和专家内模控制等效的控制效果,可以同时满足串级回路控制的鲁棒性和跟踪性要求。
术语解释:
1、PID控制器:PID(比例-积分-微分)控制器是工业控制应用中常见的一类控制器,其主要参数包括比例K、积分时间Ti和微分时间Td。
2、内模控制(Internal Model Control,简称IMC):是一种基于过程数学模型进行控制器设计的新型控制策略。
本发明克服其技术问题所采用的技术方案是:
一种基于专家内模控制的串级回路PID控制器参数整定方法,串级回路包括主回路和副回路,至少包括:
利用串级回路的数据,采用网格化搜索寻优,得到主副回路的辨识模型;
采用基于专家经验的内模控制方法,进行副回路控制器参数的整定;
通过副回路闭环阶跃仿真特征参数与主回路的辨识模型,推出等效串级回路开环模型;
采用基于专家经验的二自由度内模控制方法,进行主回路控制器参数的整定。
进一步地,利用串级回路的数据,采用网格化搜索寻优,得到主副回路的辨识模型,包括如下过程:
建立主副回路的过程模型;
对于主副回路的过程模型中待辨识的模型参数,进行动态分割并构成网络;
通过网格遍历,寻找最优的辨识模型参数,得到主副回路的辨识模型。
进一步地,用基于专家经验的内模控制方法,进行副回路控制器参数的整定,包括如下过程:
设计基于专家经验的副回路内模控制器:对副回路的内部过程模型进行分解;对副回路控制器进行确定,得到副回路内模控制器;
进行副回路控制器参数的整定:根据副回路内模控制器和副回路控制器的等价关系,得到副回路控制器的参数。
进一步地,通过副回路闭环阶跃仿真特征参数与主回路的辨识模型,推出等效串级回路开环模型,包括如下过程:
根据副回路控制器和副回路辨识模型,进行副回路的闭环阶跃响应仿真,得到副回路的闭环近似模型特征参数,包括静态增益和上升时间;
结合主回路的辨识模型,进行串级回路近似模型转换。
进一步地,采用基于专家经验的二自由度内模控制方法,进行主回路控制器参数的整定,包括如下过程:
设计基于专家经验二自由度的主回路内模控制器:对主回路的内部过程模型进行分解;对主回路控制器进行确定,得到主回路内模控制器;
进行主回路控制器参数的整定:根据主回路内模控制器和主回路控制器的等价关系,得到主回路控制器的参数。
进一步地,利用串级回路的数据,采用网格化搜索寻优,得到主副回路的辨识模型,具体包括:
通过二阶加时滞模型,将过程模型P表示为:
上式中,s为拉普拉斯算子;K为模型增益;T1为时间常数1;T2为时间常数2,当T2=0时,此过程模型为一阶模型;τ为时滞时间;
需要说明的是:本发明中,P代表用于辨识的过程模型,M代表内部过程模型,两者实际上都是过程模型,其对应的模型是一样的,是为了在不同计算阶段做区分使用;P在辨识前称为过程模型,辨识后成为辨识模型,M称为内部过程模型。
对于待辨识的模型参数,设定上下限范围,在设定的上下限范围内进行动态分割并构成网格,再进行网格遍历以寻找最优的辨识模型参数,从而实现模型辨识的过程;
在第i(i=1,2,3,...,n)次网格计算中,给定时滞时间τi、时间常数T1i和T2i,求出最优的模型增益Ki;在每个网格对应的参数下,计算辨识模型输出和真实输出的误差,寻找误差最小的网格参数和模型增益作为最优模型参数;
对于给定的时滞时间τ、时间常数T1和T2,其优化函数表示为:
上式中,yP为过程真实输出,yM为辨识模型输出,为根据已知条件得到的无增益模型输出;
令:
由此得到最优增益Kopt:
据此得到主副回路的辨识模型;
副回路的辨识模型为:
上式中,K2为副回路的模型增益;T21为副回路的时间常数1;T22为副回路的时间常数2,τ2为副回路的时滞时间;
主回路的辨识模型为:
上式中,K1为主回路的模型增益;T11为主回路的时间常数1;T12为主回路的时间常数2,τ1为主回路的时滞时间。
进一步地,采用基于专家经验的内模控制方法,进行副回路控制器参数的整定,具体包括:
基于副回路的跟随性要求,设副回路采用PI控制器,则副回路控制器的结构为:
上式中,Kp2为副回路控制器的比例系数,Ti2为副回路控制器的积分时间;
当副回路控制器C2和副回路内模控制器C2 *满足以下关系时,副回路内模控制器C2 *下的副回路控制效果与副回路控制器C2下的副回路控制效果是等价的;
首先,设计基于专家经验的副回路内模控制器:
1)对副回路的内部过程模型M2进行分解:
M2=M2+M2-
上式中,M2+包含所有纯滞后环节和s平面右半边的零点(所述s平面,是指通过拉普拉斯变换,将时域函数变换到复频域中,并分别以实轴、虚轴为横纵坐标轴组成的平面,通过分析零、极点在s平面的分布情况可判断系统稳定性,s平面右半边的零点可能会导致系统不稳定),并设其稳态增益为1;
则此处,M2+和M2-分别为:
2)对副回路控制器进行确定:
在副回路控制器中,增加一个低通滤波器f2:
上式中,λ2为低通滤波器f2的时间常数,m为正整数;
则副回路控制器C2 *为:
通过专家经验,结合副回路的内部过程模型M2中的时间常数和时滞时间与时间常数的比值,给定λ2的数值;其主要调节原理为:当模型时间常数较大时(即慢过程),调小λ2;当模型时间常数较小时(快过程),调大λ2;当模型时滞时间与时间常数的比值较大时(即大时滞过程),适当调小λ2;当模型时滞时间与时间常数的比值较小时(即小时滞或无时滞过程),适当调大λ2。
然后,进行副回路控制器参数的整定:根据副回路内模控制器C2 *和副回路控制器C2的等价关系,得出副回路控制器的参数为:
当副回路的时滞时间τ2≥0.5(T21+T22)时,采用一阶派德近似:
副回路控制器的参数为:
当时滞时间τ2<0.5(T21+T22)时,采用一阶泰勒近似:
副回路控制器的参数为:
进一步地,通过副回路闭环阶跃仿真特征参数与主回路的辨识模型,推出等效串级回路开环模型,具体包括:
根据副回路控制器和副回路辨识模型,进行副回路的闭环阶跃响应仿真,提取闭环响应的特性参数,得到副回路的闭环近似模型:
t2=tr/2.2
上式中,A2为副回路闭环响应静态增益,t2为副回路闭环时间常数,tr为副回路闭环响应上升时间;
一阶系统的单位阶跃响应为:
c(t)=1-e-t/T
上式中,T为一阶系统的时间常数;
副回路闭环响应上升时间tr定义为从阶跃响应终值的10%上升到终值90%所需的时间;
据此计算得到:
tr=2.20T
结合主回路的辨识模型,进行串级回路近似模型转换,如下:
将串级回路模型等效为二阶滞后模型:
当主回路为比例模型时:串级回路等效增益K1等于主回路模型P1的增益K1和副回路闭环响应静态增益A2的乘积:
K=K1·A2
串级回路等效时间常数计算如下:
T1=t2
T2=0
串级回路等效时滞时间等于主副回路时滞时间之和:
τ=τ1+τ2
当主回路为一阶模型时:串级回路等效增益K等于主回路的模型增益K1和副回路闭环响应静态增益A2的乘积:
K=K1·A2
串级回路等效时间常数计算如下:
T1=max{T11,t2}
T2=min{T11,t2}
串级回路等效时滞时间等于主回路的时滞时间与副回路的时滞时间之和:
τ=τ1+τ2
当主回路为二阶模型时:串级回路等效增益K等于主回路的模型增益K1和副回路闭环响应静态增益A2的乘积:
K=K1·A2
串级回路等效时间常数T1和T2为主回路的两个时间常数T11、T12和副回路的闭环时间常数t2之中,三者取较大的两个;
串级回路等效时滞时间等于主回路的时滞时间与副回路的时滞时间之和:
τ=τ1+τ2。
进一步地,采用基于专家经验的二自由度内模控制方法,进行主回路控制器参数的整定,具体包括:
设主回路采用PID控制器,则主回路控制器的结构为:
上式中,Kp1为主回路控制器的比例系数,Ti1为主回路控制器的积分时间,Td1为主回路控制器的微分时间;
当主回路控制器C1和主回路内模控制器C1 *满足以下关系时,主回路内模控制器C1 *下的串级回路控制效果与主回路控制器C1下的串级回路控制效果是等价的;
首先,设计基于专家经验二自由度的主回路内模控制器:
1)对主回路的内部过程模型M1进行分解:
M1=M1+M1-
上式中,M1+包含所有纯滞后环节和s平面右半边的零点,并设其稳态增益为1;
则此处,M1+和M1-分别为:
M1+=e-τs
2)对主回路控制器进行确定:
在主回路控制器中,增加一个低通滤波器f1:
上式中,λ1为低通滤波器f1的时间常数;n为正整数;
则主回路控制器C1 *为:
为了使串级回路的鲁棒性和跟踪性能够独立调整,设计滤波控制器F1为:
其中,f滤波器为:
上式中,λ为f滤波器的时间常数;r为正整数;
通过专家经验,结合主回路的内部过程模型M1中的时间常数和时滞时间与时间常数的比值,给定λ1和λ的数值;其主要调节原理为:当模型时间常数较大时(即慢过程),调小λ1和λ;当模型时间常数较小时(快过程),调大λ1和λ;当模型时滞时间与时间常数的比值较大时(即大时滞过程),适当调小λ1和λ;当模型时滞时间与时间常数的比值较小时(即小时滞或无时滞过程),适当调大λ1和λ。
然后,进行主回路控制器参数的整定:根据主回路内模控制器C1 *和主回路控制器C1的等价关系,得出主回路控制器的参数:
当时滞时间τ≥0.5(T1+T2)时,采用一阶派德近似:
主回路控制器的参数为:
当时滞时间τ<0.5(T1+T2)时,采用一阶泰勒近似:
e-τs≈1-τs
主回路控制器的参数为:
本发明的有益效果是:
1、本发明采用开环或闭环任意状态下的具有激励的串级回路数据进行整定,应用场景更广泛。
2、本发明提出的串级回路整定方法具有普遍适用性,通过副回路闭环阶跃仿真的特征参数结合主回路模型得到的等效串级回路开环模型在保证等效模型的稳定性的同时,降低了模型的复杂度,方便运算。
3、采用基于专家经验的内模整定和二自由度内模整定方法通过模型参数结合专家经验调整滤波器时间常数,避免了用户手动调整该项参数,提高了用户易用度。
4、本发明提出的基于专家经验二自由度内模控制的串级回路PID控制器参数整定方法适用于串级回路的参数整定,无需改变当前工况,即可同时得到主副回路的PID控制器整定参数,整定效率高。
具体实施方式
为了便于本领域人员更好的理解本发明,下面结合附图和具体实施例对本发明做进一步详细说明,下述仅是示例性的不限定本发明的保护范围。
下面通过一个具体实例脱碳九塔中的一个串级回路为例展开对本发明所述的基于专家内模控制的串级回路PID控制器参数整定方法进行详细说明:
主回路对塔下部温度进行调节,副回路对再沸器入口处过热蒸汽的进口流量进行调节。如图2所示,本实施例所述的串级回路中,C1为主回路的PID控制器,P1为主回路的辨识模型,C2为副回路的PID控制器,P2为副回路的辨识模型,SV1为主回路的塔下部温度设定值,SV2为副回路的进口流量设定值,OP1为主回路控制器的输出值,OP2为副回路控制器的输出值,PV1为主回路塔下部温度的测量值,PV2为副回路进口流量的测量值。副回路的设定值SV2等于主回路控制器的输出值OP1。PV1的量程为[0,200℃],OP1的量程为[0,10000t/h],PV2的量程为[0,10000t/h],OP2的量程为[0,100%]。
该串级回路的趋势线如图5所示,其中,图5(a)为整定前主回路温度设定值SV1和主回路温度测量值PV1的趋势图,图5(b)为整定前主回路控制器输出OP1的趋势图,图5(c)为整定前副回路流量设定值SV2和副回路流量测量值PV2的趋势图,图5(d)为整定前副回路控制器输出OP1的趋势图。从该组趋势图中可知,整定前该串级回路的副回路阀门的平稳性较差,存在频繁的振荡,容易降低阀门使用寿命;过快的响应速度,使副回路流量存在微小的超调;且主回路的温度控制存在一定的超调,不符合工艺上需要温度控制平稳的要求。基于这样的背景,下面展开对该串级回路的PID参数整定。
对主副回路均采用一阶模型进行辨识,其模型形式为:
主回路辨识模型:
副回路辨识模型:
通过前述的网格遍历寻优,得到的主副回路的模型参数分别为:
主回路模型辨识参数:
K1=0.017,T11=5.96min,τ1=1min
副回路模型辨识参数:
K2=97.13,T21=0.54min
首先根据专家经验内模控制整定副回路的PID控制器C2的参数。对于副回路,通过内部过程模型M2和控制器输出OP2,可以得到模型响应PV2 *。真实的响应PV2与模型响应PV2*的差PV2﹣PV2*作为内模控制器C2*的反馈输入信号。通常,由于存在模型误差M2≠G2和未被考虑在模型中的未知扰动D2,因此PV2≠PV2 *。在这个过程中,通过专家经验,依据内部过程模型M2对C2 *中的参数进行调节。
基于副回路的跟随性要求,副回路考虑采用PI控制器:
可以发现,当C2和C2 *满足以下关系时,图3和图4的副回路的框图是等价的。
下面设计基于专家经验的副回路内模控制器:
首先对副回路内部过程模型进行分解:
M2=M2+M2-
M2+=1
然后对副回路控制器进行确定:
为了增加控制器的鲁棒性,加入一个低通滤波器f2:
其中,时间常数初始值取λ2=3(T21+τ2)=97.2。
据此设计副回路控制器为:
滤波器的时间常数λ2的取值通过专家经验,结合内部过程模型M2中的时间常数和时滞时间与时间常数的比值的不同给定合适的值。本实施例中,对λ2的具体调节如表1和表2所示。
表1流量回路的过程响应速度与λ2调节量之间的关系
流量回路 |
慢 |
较慢 |
中 |
较快 |
快 |
λ2调节量 |
-0.8(T21+τ2) |
-0.4(T21+τ2) |
0 |
0.4(T21+τ2) |
0.8(T21+τ2) |
表2流量回路的时滞程度与λ2调节量之间的关系
流量回路 |
无时滞 |
小时滞 |
较大时滞 |
大时滞 |
λ2调节量 |
0.2(T1+τ) |
0.1(T1+τ) |
-0.1(T1+τ) |
-0.2(T1+τ) |
本实施例所述的该流量回路作为副回路,属于较慢过程,无时滞,根据上表调节得到最终的λ2=90.68。
根据等价关系,采用一阶泰勒近似导出PI控制器参数为:
根据副回路控制器和副回路辨识模型,进行副回路的闭环阶跃响应仿真,如图6所示。
提取闭环响应的特性参数,得到副回路闭环响应静态增益A2=1和副回路闭环响应上升时间tr=100s,从而得到副回路的闭环近似模型:
t2=tr/2.2=45.45
将串级回路模型等效为二阶滞后模型,则等效串级回路开环模型为:
最后展开基于专家经验二自由度内模的串级回路PID控制器参数整定。在内模控制器C1 *的基础上,增加F1用于调节系统的跟踪性,其作用等效于对SV1进行了滤波处理。M1为串级回路内部过程模型,通过OP1可以得到串级回路的模型响应PV1 *。真实的响应PV1与模型响应PV1 *的差PV1﹣PV1 *作为内模控制器C1 *的反馈输入信号。通常,由于存在模型误差M1≠G’和未被考虑在模型中的未知扰动D1,因此PV1≠PV1 *。在这个过程中,通过专家经验,依据内部过程模型M1对C1 *和F1中的参数进行调节。
对主回路采用PI控制器,其控制器结构为:
当C1和C1 *满足以下关系时,图3和图4的串级回路整体的框图是等价的。
下面设计基于专家经验二自由度的主回路内模控制器:
首先分解主回路的内部过程模型:
M1=M1+M1-
M1+=e-τs=e-60s
然后确定主回路的控制器:
为了增加控制器的鲁棒性,加入一个低通滤波器f1:
其中,时间常数初始值取λ1=2(T1+τ)=835.2。
据此设计主回路控制器为:
设计滤波器f为:
其中,时间常数初始值取λ=(T1+τ)=417.6。
滤波控制器F1为:
在上述主回路控制器设计中,两个滤波器的时间常数λ1和λ的取值通过专家经验,结合模型M1中的时间常数和时滞时间与时间常数的比值的不同给定合适的值。本实施例中,λ1和λ的具体调节如表3和表4所示。
表3温度回路的过程响应速度与λ1、λ调节量之间的关系
温度回路 |
慢 |
较慢 |
中 |
较快 |
快 |
λ1调节量 |
-(T1+τ) |
-0.5(T1+τ) |
0 |
0.5(T1+τ) |
(T1+τ) |
λ调节量 |
-0.4(T1+τ) |
-0.2(T1+τ) |
0 |
0.2(T1+τ) |
0.4(T1+τ) |
表4温度回路的时滞程度与λ1、λ调节量之间的关系
温度回路 |
无时滞 |
小时滞 |
较大时滞 |
大时滞 |
λ1调节量 |
0.2(T1+τ) |
0.1(T1+τ) |
-0.1(T1+τ) |
-0.2(T1+τ) |
λ调节量 |
0.08(T1+τ) |
0.04(T1+τ) |
-0.04(T1+τ) |
-0.08(T1+τ) |
本实施例所述的该温度回路作为主回路,属于慢过程,小时滞,根据上表调节参数得到最终的λ1=459.4,λ=267.3。
根据等价关系,采用一阶泰勒近似导出主回路PI控制器参数为:
经过量程转换和比例度换算,最终得到整定后的主副回路控制器参数分别为:
最终的整定效果如图7所示,其中,图7(a)为整定后主回路温度设定值SV1和主回路温度测量值PV1的趋势图,图7(b)为整定后主回路控制器输出OP1的趋势图,图7(c)为整定后副回路流量设定值SV2和副回路流量测量值PV2的趋势图,图7(d)为整定后副回路控制器输出OP1的趋势图。整定后,消除了整定前副回路的阀门高频率的振荡,同时流量控制保持了一定的跟随性,且在相同外界扰动的情况下,使温度控制较稳定且具有一定的跟随性。
以上仅描述了本发明的基本原理和优选实施方式,本领域人员可以根据上述描述做出许多变化和改进,这些变化和改进应该属于本发明的保护范围。