CN113110065A - 基于双rbf神经网络的反渗透膜组压力优化控制方法 - Google Patents

基于双rbf神经网络的反渗透膜组压力优化控制方法 Download PDF

Info

Publication number
CN113110065A
CN113110065A CN202110519454.XA CN202110519454A CN113110065A CN 113110065 A CN113110065 A CN 113110065A CN 202110519454 A CN202110519454 A CN 202110519454A CN 113110065 A CN113110065 A CN 113110065A
Authority
CN
China
Prior art keywords
pressure
reverse osmosis
membrane
osmosis membrane
neural network
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110519454.XA
Other languages
English (en)
Other versions
CN113110065B (zh
Inventor
褚晓广
王恬
王铭涛
孔英
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qufu Normal University
Original Assignee
Qufu Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qufu Normal University filed Critical Qufu Normal University
Priority to CN202110519454.XA priority Critical patent/CN113110065B/zh
Publication of CN113110065A publication Critical patent/CN113110065A/zh
Application granted granted Critical
Publication of CN113110065B publication Critical patent/CN113110065B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/124—Water desalination
    • Y02A20/131—Reverse-osmosis

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Separation Using Semi-Permeable Membranes (AREA)

Abstract

本发明公开了一种基于双RBF神经网络的反渗透膜组压力优化控制方法，构建反渗透膜组首段膜系统压力动态模型及其它各段膜稳态模型，建立考虑膜组各段膜压力综合最优的优化目标函数，根据实际滤饼污染层厚度在线调整压力优化值，采用RBF神经网络以实际膜组反渗透效率逼近最优反渗透效率为目标对压力优化值在线调整，在状态反馈控制器的基础上以RBF神经网络自适应补偿控制抑制系统的未知扰动，实现膜组各段压力和反渗透效率双重意义下的压力最优、变工况下优化压力的快速跟踪及反渗透膜组的稳定运行。

Description

基于双RBF神经网络的反渗透膜组压力优化控制方法
技术领域
本发明提出了一种基于双RBF神经网络的反渗透膜组压力优化控制方法，是一种应用于多膜反渗透膜组变工况下压力实时优化及跟踪控制的方法，解决了系统参数变化及未知干扰导致的优化压力不准确和系统运行不稳定的问题，属于反渗透海水淡化系统控制领域。
背景技术
多膜反渗透膜组海水淡化系统是目前主流的海水淡化系统，存在运行参数多、参数时变、系统耦合性强、未知干扰多、工况变化频繁等特点。因此，在线优化系统压力，降低压力波动是反渗透膜组优化控制研究的关键。为此，曲阜师范大学新能源研究所提出基于双RBF神经网络的反渗透膜组压力优化控制方法，提高了优化压力的准确性和系统运行的稳定性。由于反渗透膜组中膜的数量多，参数实时变化，系统的优化压力难以实时达到最优；由于系统受到膜组弹性干扰及电磁阀绕组温升干扰的影响，系统在运行过程汇总压力不稳定，造成反渗透膜损伤。传统的状态反馈控制方法虽可实现优化压力的跟踪，但压力跟踪速度慢，且由于状态反馈控制不具备抗干扰能力，在变工况时压力波动大，稳态运行时的稳定性较差，系统能耗较高，膜承受严重的压力冲击。
发明内容
本发明提供了一种基于双RBF神经网络的反渗透膜组压力优化控制方法，包括变工况下优化压力的获取、基于RBF神经网络的优化压力在线调整、采用RBF神经网络自适应补偿的优化压力控制。所述优化压力的获取以变海水盐度下额定优化压力为单膜优化压力初始值，以各段膜压力综合最优为目标构建优化目标函数，采用拉格朗日乘数法获取反渗透膜组首段膜系统压力优化值；所述优化压力在线调整以膜组反渗透效率为性能评价指标，采用RBF神经网络令实际反渗透效率逼近最优反渗透效率，在线调整首段膜系统压力优化值；所述RBF神经网络自适应补偿控制采用基于状态反馈的主压力跟踪控制器和基于RBF神经网络自适应的干扰补偿控制器。所述主压力跟踪控制器采用状态反馈控制器，以优化压力跟踪误差及其导数作为状态反馈控制输入；所述RBF神经网络补偿控制器以优化压力跟踪误差及其导数为控制输入，自适应地调整神经网络权值，以输出的干扰估计值逼近系统中的未知扰动。
步骤1构建用于优化控制的首段膜系统压力动态模型及其它各段膜压力稳态模型、用于优化性能评价的各段膜产水流量稳态模型。
Figure BDA0003063362810000011
其中，ρ为浓水区液体密度，Pi(i＝1，2，……，n)为各段膜的系统压力，ηs为粘度系数，Vvr浓水阀液体体积，evr为浓水阀阻，vr为浓水流速，Qf为进水流量，λ为沿程阻力系数，l为膜元件轴向长度，d为浓水区管道当量直径，Am为膜面积，Ap为水管截面积，
Figure BDA0003063362810000021
为因膜组弹性造成的浓水区平均流量变化量，CH为膜组液容，Qpi(i＝1,2,3,…，n)为各段膜产水流量，Rvi(i＝1,2,3,4)为各段膜渗透液阻，Δπi(i＝1,2,3,4)为各段膜渗透液压差，Ki(i＝1,2,…,7)、C1为膜间关联系数。
步骤2反渗透膜组首段膜系统压力在线优化
第一步，以不同海水盐度下单膜额定压力PNref为输入优化目标函数的单膜优化压力初始值。
PNref＝0.02117S5-3.281S4+204.6S3-6637S2+241500S+653700 (2)
其中，S为海水盐度。
第二步，以膜组各段膜压力综合最优为目标构建优化目标函数。
minG＝γ1(P1-PNref)2+γ2(P2-PNref)2+…+γn(Pn-PNref)2 (3)
其中Pi(i＝1，2，……，n)为第i段膜的系统压力，γi(i＝1，2，……，n)为第i段膜的优化权值。
第三步，设定反渗透膜组压力优化的约束条件：
Figure BDA0003063362810000022
其中，Pf、Qf为膜组进水压力及进水流量，Pdi(i＝1,2,3,…，n)表示n个反渗透膜的管压降，Qi(i＝1,2,3,…，n)为各段膜的轴向流量，Pr为浓水出口处的压力，Qr为浓水出口的流量，Psysmin、Psysmax为反渗透膜所允许的系统压力最大值和最小值，Smin、Smax为送入反渗透膜组海水的最小盐度和最大盐度。
第四步，采用拉格朗日乘数法获得反渗透膜组首段膜额定系统压力优化值。
Figure BDA0003063362810000023
第五步，根据实时采集的反渗透膜组产水流量计算反渗透膜滤饼污染层厚度。
Figure BDA0003063362810000024
其中，Qpi为第i段膜的产水流量，α为滤饼层累积系数，Am为反渗透膜面积，
Figure BDA0003063362810000025
为海水平均污染浓度。
第六步，根据实时获得的滤饼层厚度，更新反渗透膜渗透液阻，获得变海水盐度下的反渗透膜组反渗透效率和最大反渗透效率ηrM对应的首段膜系统压力
Figure BDA0003063362810000026
第七步，根据实时监测的反渗透膜组各段膜产水流量计算膜组的实际反渗透效率。
Figure BDA0003063362810000027
其中，Qpi为第i个膜的产水流量，Δπi为渗透压差，Pf为膜组的进水压力，Qf为膜组的进水流量。
第八步，采用RBF神经网络实现反渗透膜组压力在线优化，计算最大反渗透效率逼近误差eη＝ηrM-ηr，以eηr和
Figure BDA0003063362810000028
为神经网络输入,即
Figure BDA0003063362810000029
第九步，计算RBF神经网络的高斯基函数：
Figure BDA0003063362810000031
其中，cj为高斯基函数中心，bj为高斯基函数宽度。
第十步，得到RBF神经网络输出的优化压力补偿值：
Figure BDA0003063362810000032
其中，ωj为输出层权值，hj为高斯基函数。
第十一步，将额定压力优化值与优化压力补偿值相加，获得首段膜系统压力优化值：
P1ref＝P1Nref+ΔP1ref (10)
第十二步，设定RBF神经网络权值在线优化最小目标函数：
Figure BDA0003063362810000033
第十三步，根据梯度下降法，在线调整RBF神经网络权值。
Figure BDA0003063362810000034
wj(N)＝wj(N-1)+Δwj(N)+k(wj(N-1)-wj(N-2)) (13)
其中，Δwj(N)为RBF神经网络权值的修正值，η为学习率，满足η∈[0,1]。
步骤3构建包含首段膜系统压力动态模型与电磁阀线圈电流模型的四膜反渗透膜组模型
Figure BDA0003063362810000035
其中，evrmax为最大阀阻，imax为电磁阀最大励磁电流，
Figure BDA0003063362810000036
为膜组弹性形变扰动，ΔRi为电磁阀绕组温升扰动。上式中，第四段膜的系统压力P4、膜组浓水流速vr以及函数l(P1)形如
Figure BDA0003063362810000037
上述模型中的Ki(i＝1,2,3,…，23)、C1、C2、C3均为膜间关联系数。
步骤4反渗透膜组海水淡化系统模型转化
第一步，将如式(14)所示的反渗透膜组模型简化如下：
Figure BDA0003063362810000038
其中，
Figure BDA0003063362810000041
Figure BDA0003063362810000042
第二步，求取反渗透膜组海水淡化系统平衡点
忽略电流模型中的扰动，令电流的变化率为零，即求解式(19)，得到的电流即为平衡点处的电流i0。
Figure BDA0003063362810000043
其中，uin0为电磁阀外电路输入电压uin的稳态值。
将求得的系统平衡点电流i0带入式(20)计算平衡点的压力，得到平衡点处的压力P1,0。
Figure BDA0003063362810000044
由上式求得的压力即为平衡点的压力P1,0，该系统的平衡点即为(P1,0,i0)。
第三步，在反渗透膜组系统平衡点处得到线性化的反渗透膜组模型为：
Figure BDA0003063362810000045
其中，干扰项归结如下：
Figure BDA0003063362810000046
步骤5主压力跟踪控制器设计
第一步，以P1为状态变量x1，以
Figure BDA0003063362810000047
为状态变量x2，将模型中的扰动项统一归结为干扰项fxp，写出反渗透膜组淡化系统的状态空间方程。
Figure BDA0003063362810000048
其中，扰动项为
Figure BDA0003063362810000049
第二步，计算压力跟踪误差e＝P1ref-P1＝P1ref-x1，
Figure BDA00030633628100000410
定义虚拟变量
Figure BDA00030633628100000411
写出系统压力跟踪误差增广模型：
Figure BDA0003063362810000051
其中，虚拟控制输入为
Figure BDA0003063362810000052
第三步，对式(25)进行极点配置，采用状态反馈法设计主压力跟踪控制器，写出反渗透膜组淡化系统闭环特征多项式：
Figure BDA0003063362810000053
其中，KE＝[-kp -kd]为状态反馈增益矩阵。
第四步，通过选取合适的闭环极点，得到反馈增益矩阵KE的参数值，则此时的控制输入为：
Figure BDA0003063362810000054
步骤6基于RBF神经网络自适应的干扰补偿控制器设计
第一步，在状态反馈的基础上引入RBF神经网络自适应补偿，设计控制律为：
Figure BDA0003063362810000055
第二步，将(29)代入式(25)，得到状态空间方程：
Figure BDA0003063362810000056
第三步，定义最优权值为
Figure BDA0003063362810000057
式中，Ω表示具有希望边界的子集，即Ω＝{W|||f||≤M}，M为设计参数。
第四步，以RBF神经网络输出的自适应补偿
Figure BDA0003063362810000058
逼近干扰项fxp。
Figure BDA0003063362810000059
第五步，定义模型逼近误差为
Figure BDA00030633628100000510
第六步，将式(33)代入式(30)，将状态空间方程表示为：
Figure BDA00030633628100000511
其中，
Figure BDA0003063362810000061
第七步，将式(32)代入式(34)，将闭环系统状态空间方程改写为：
Figure BDA0003063362810000062
第八步，设计Lyapunov函数为
Figure BDA0003063362810000063
其中，γ是一个正实数，P为对称正定矩阵，且满足Lyapunov方程ΛTP+PΛ＝-Q，Q为对称正定矩阵。
第九步，取
Figure BDA0003063362810000064
对V1求导，可得
Figure BDA0003063362810000065
其中
Figure BDA0003063362810000066
则有
Figure BDA0003063362810000067
第十步，对V2求导，可得
Figure BDA0003063362810000068
第十一步，计算V的导数：
Figure BDA0003063362810000069
第十二步，为确保
Figure BDA00030633628100000610
设计权值自适应律：
Figure BDA00030633628100000611
其中，γ是一个正实数，P为对称正定矩阵，且满足ΛTP+PΛ＝-Q，Q为对称正定矩阵。
本发明的有益效果是：
1)提出了基于RBF神经网络的优化压力在线调整，依靠神经网络强大的学习能力，在线调整目标函数的结构，实现各段膜系统压力与膜组反渗透效率的综合最优。
2)构建污染层厚度与优化压力的关系，实现了考虑污染程度变化的反渗透膜组压力最优，扩大了优化方法的适用范围。
3)提出了基于RBF神经网络自适应补偿的优化压力控制，使系统具备抑制干扰的能力，提升了系统压力跟踪的响应速度，大幅度降低了变工况及稳态运行时的压力波动，降低系统能耗，减少因压力冲击造成的膜损伤，延长了反渗透膜的使用寿命。
附图说明
图1为本发明基于双RBF神经网络的反渗透膜组压力优化控制方法的额定工况下单膜最大反渗透效率曲线。
图2为本发明基于双RBF神经网络的反渗透膜组压力优化控制方法的反渗透膜组最大反渗透效率曲线。
图3为本发明基于双RBF神经网络的反渗透膜组压力优化控制方法的反渗透膜组控制结构图。
图4为本发明基于双RBF神经网络的反渗透膜组压力优化控制方法的控制与PID控制下的变工况压力跟踪实验图。
图5为本发明基于双RBF神经网络的反渗透膜组压力优化控制方法的控制与PID控制下的变工况膜组反渗透效率实验图。
图中：1-额定工况下单膜优化压力获取，2-反渗透膜组首段膜系统压力离线优化，3-膜组最优反渗透效率及最优压力实时计算，4-膜组实际反渗透效率计算，5-基于RBF神经网络的优化模型在线调整，6-状态反馈控制，7-RBF神经网络自适应补偿控制，8-反渗透膜组海水淡化系统模型，9-反渗透膜滤饼层厚度实时计算。
具体实施方式
本发明提供了一种基于双RBF神经网络的反渗透膜组压力优化控制方法，包括变工况下优化压力的获取、基于RBF神经网络的优化压力在线调整、采用RBF神经网络自适应补偿的优化压力控制。所述优化压力的获取以变海水盐度下额定优化压力为单膜优化压力初始值，以各段膜压力综合最优为目标构建优化目标函数，采用拉格朗日乘数法获取反渗透膜组首段膜系统压力优化值；所述优化压力在线调整以膜组反渗透效率为性能评价指标，采用RBF神经网络令实际反渗透效率逼近最优反渗透效率，在线调整首段膜系统压力优化值；所述RBF神经网络自适应补偿控制采用基于状态反馈的主压力跟踪控制器和基于RBF神经网络自适应的干扰补偿控制器。所述主压力跟踪控制器采用状态反馈控制器，以优化压力跟踪误差及其导数作为状态反馈控制输入；所述RBF神经网络补偿控制器以优化压力跟踪误差及其导数为控制输入，自适应地调整神经网络权值，以输出的干扰估计值逼近系统中的未知扰动。
步骤1构建用于优化控制的首段膜系统压力动态模型及其它各段膜压力稳态模型、用于优化性能评价的各段膜产水流量稳态模型。
Figure BDA0003063362810000081
其中，ρ为浓水区液体密度，Pi(i＝1，2，……，n)为各段膜的系统压力，ηs为粘度系数，Vvr浓水阀液体体积，evr为浓水阀阻，vr为浓水流速，Qf为进水流量，λ为沿程阻力系数，l为膜元件轴向长度，d为浓水区管道当量直径，Am为膜面积，Ap为水管截面积，
Figure BDA0003063362810000082
为因膜组弹性造成的浓水区平均流量变化量，CH为膜组液容，Qpi(i＝1,2,3,…，n)为各段膜产水流量，Rvi(i＝1,2,3,4)为各段膜渗透液阻，Δπi(i＝1,2,3,4)为各段膜渗透液压差，Ki(i＝1,2,…,7)、C1为膜间关联系数。
步骤2反渗透膜组首段膜系统压力在线优化
第一步，以不同海水盐度下单膜额定压力PNref为输入优化目标函数的单膜优化压力初始值。
PNref＝0.02117S5-3.281S4+204.6S3-6637S2+241500S+653700 (2)
其中，S为海水盐度。
第二步，以膜组各段膜压力综合最优为目标构建优化目标函数。
minG＝γ1(P1-PNref)2+γ2(P2-PNref)2+…+γn(Pn-PNref)2 (3)
其中Pi(i＝1，2，……，n)为第i段膜的系统压力，γi(i＝1，2，……，n)为第i段膜的优化权值。
第三步，设定反渗透膜组压力优化的约束条件：
Figure BDA0003063362810000083
其中，Pf为膜组进水压力，Pdi(i＝1,2,3,…，n)表示n个反渗透膜的管压降，Qi(i＝1,2,3,…，n)为各段膜的轴向流量，Pr为浓水出口处的压力，Qr为浓水出口的流量，Psysmin、Psysmax为反渗透膜所允许的系统压力最大值和最小值，Smin、Smax为送入反渗透膜组海水的最小盐度和最大盐度。
第四步，采用拉格朗日乘数法获得反渗透膜组首段膜额定系统压力优化值。
Figure BDA0003063362810000091
第五步，根据实时采集的反渗透膜组产水流量计算反渗透膜滤饼污染层厚度。
Figure BDA0003063362810000092
其中，Qpi为第i段膜的产水流量，α为滤饼层累积系数，Am为反渗透膜面积，
Figure BDA0003063362810000093
为海水平均污染浓度。
第六步，根据实时获得的滤饼层厚度，更新反渗透膜渗透液阻，获得变海水盐度下的反渗透膜组反渗透效率和最大反渗透效率ηrM对应的首段膜系统压力
Figure BDA0003063362810000094
第七步，根据实时监测的反渗透膜组各段膜产水流量计算膜组的实际反渗透效率。
Figure BDA0003063362810000095
其中，Qpi为第i个膜的产水流量，Δπi为渗透压差，Pf为膜组的进水压力，Qf为膜组的进水流量。
第八步，采用RBF神经网络实现反渗透膜组压力在线优化，计算最大反渗透效率逼近误差eη＝ηrM-ηr，以eηr和
Figure BDA0003063362810000096
为神经网络输入,即
Figure BDA0003063362810000097
第九步，计算RBF神经网络的高斯基函数：
Figure BDA0003063362810000098
其中，cj为高斯基函数中心，bj为高斯基函数宽度。
第十步，得到RBF神经网络输出的优化压力补偿值：
Figure BDA0003063362810000099
其中，ωj为输出层权值，hj为高斯基函数。
第十一步，将额定压力优化值与优化压力补偿值相加，获得首段膜系统压力优化值：
P1ref＝P1Nref+ΔP1ref (10)
第十二步，设定RBF神经网络权值在线优化最小目标函数：
Figure BDA00030633628100000910
第十三步，根据梯度下降法，在线调整RBF神经网络权值。
Figure BDA00030633628100000911
wj(N)＝wj(N-1)+Δwj(N)+k(wj(N-1)-wj(N-2)) (13)
其中，Δwj(N)为RBF神经网络权值的修正值，η为学习率，满足η∈[0,1]。
步骤3构建包含首段膜系统压力动态模型与电磁阀线圈电流模型的四膜反渗透膜组模型
Figure BDA0003063362810000101
其中，evrmax为最大阀阻，imax为电磁阀最大励磁电流，
Figure BDA0003063362810000102
为膜组弹性形变扰动，ΔRi为电磁阀绕组温升扰动。上式中，第四段膜的系统压力P4、膜组浓水流速vr以及函数l(P1)形如
Figure BDA0003063362810000103
上述模型中的Ki(i＝1,2,3,…，23)、C1、C2、C3均为膜间关联系数。
步骤4反渗透膜组海水淡化系统模型转化
第一步，将如式(14)所示的反渗透膜组模型简化如下：
Figure BDA0003063362810000104
其中
Figure BDA0003063362810000105
Figure BDA0003063362810000106
第二步，求取反渗透膜组海水淡化系统平衡点
忽略电流模型中的扰动，令电流的变化率为零，即求解式(19)，得到的电流即为平衡点处的电流i0。
Figure BDA0003063362810000107
其中，uin0为电磁阀外电路输入电压uin的稳态值。
将求得的系统平衡点电流i0带入式(20)计算平衡点的压力，得到平衡点处的压力P1,0。
Figure BDA0003063362810000108
由上式求得的压力即为平衡点的压力P1,0，该系统的平衡点即为(P1,0,i0)。
第三步，在反渗透膜组系统平衡点处得到线性化的反渗透膜组模型为：
Figure BDA0003063362810000109
其中，干扰项归结如下：
Figure BDA0003063362810000111
步骤5主压力跟踪控制器设计
第一步，以P1为状态变量x1，以
Figure BDA0003063362810000112
为状态变量x2，将模型中的扰动项统一归结为干扰项fxp，写出反渗透膜组淡化系统的状态空间方程。
Figure BDA0003063362810000113
其中，扰动项为
Figure BDA0003063362810000114
第二步，计算压力跟踪误差e＝P1ref-P1＝P1ref-x1，
Figure BDA0003063362810000115
定义虚拟变量
Figure BDA0003063362810000116
写出系统压力跟踪误差增广模型：
Figure BDA0003063362810000117
其中，虚拟控制输入为
Figure BDA0003063362810000118
第三步，对式(25)进行极点配置，采用状态反馈法设计主压力跟踪控制器，写出反渗透膜组淡化系统闭环特征多项式：
Figure BDA0003063362810000119
其中，KE＝[-kp -kd]为状态反馈增益矩阵。此时系统闭环特征多项式为
第四步，通过选取合适的闭环极点，得到反馈增益矩阵KE的参数值，则此时的控制输入为：
Figure BDA00030633628100001110
步骤6基于RBF神经网络自适应的干扰补偿控制器设计
第一步，在状态反馈的基础上引入RBF神经网络自适应补偿，设计控制律为：
Figure BDA00030633628100001111
第二步，将(29)代入式(25)，得到状态空间方程：
Figure BDA0003063362810000121
第三步，定义最优权值为
Figure BDA0003063362810000122
式中，Ω表示具有希望边界的子集，即Ω＝{W|||f||≤M}，M为设计参数。
第四步，以RBF神经网络输出的自适应补偿
Figure BDA0003063362810000123
逼近干扰项fxp。
Figure BDA0003063362810000124
第五步，定义模型逼近误差为
Figure BDA0003063362810000125
第六步，将式(33)代入式(30)，将状态空间方程表示为：
Figure BDA0003063362810000126
其中，
Figure BDA0003063362810000127
第七步，将式(32)代入式(34)，将闭环系统状态空间方程改写为：
Figure BDA0003063362810000128
第八步，设计Lyapunov函数为
Figure BDA0003063362810000129
其中，γ是一个正实数，P为对称正定矩阵，且满足Lyapunov方程ΛTP+PΛ＝-Q，Q为对称正定矩阵。
第九步，取
Figure BDA00030633628100001210
对V1求导，可得
Figure BDA00030633628100001211
其中
Figure BDA00030633628100001212
则有
Figure BDA0003063362810000131
第十步，对V2求导，可得
Figure BDA0003063362810000132
第十一步，计算V的导数：
Figure BDA0003063362810000133
第十二步，为确保
Figure BDA0003063362810000134
设计权值自适应律：
Figure BDA0003063362810000135
其中，γ是一个正实数，P为对称正定矩阵，且满足ΛTP+PΛ＝-Q，Q为对称正定矩阵。
下面结合附图以及实例，对本发明作进一步详细说明。反渗透膜组海水淡化系统参数如表1所示，反渗透膜面积为15.6m2，膜透水系数为2.6×1012m/(Pa·s)，水管截面积为0.000127m2，膜元件容积为0.04m3，膜轴向长度为1m，浓水区当量直径为0.01272m，为说明本发明的效果，进行压力跟踪实验。
表1反渗透膜性能参数
Figure BDA0003063362810000136
变海水盐度优化压力跟踪实验对比图如图4所示，首段膜系统压力初始值设定为6MPa，初始工况海水盐度为28，在60s时海水盐度变为33，120s时海水盐度变回28。分析首段膜系统压力跟踪的调节时间，超调量及稳态压力波动值，三个工况下，RBF神经网络自适应补偿控制的调节时间比传统的PID控制快19.5s以上，其稳态压力波动值降低至PID控制的0.047倍以下。在两次变工况时，PID控制的超调量分别为27.6kPa、26.8kPa，而RBF神经网络自适应补偿控制实现压力无超调。显然，RBF神经网络自适应控制提高了动态响应速度，大幅度降低了变工况时的压力波动和稳态压力波动，提高了反渗透效率，减轻膜损伤。变海水盐度下的优化压力跟踪性能对比如表2所示。
表2不同控制策略下压力跟踪仿真实验结果
Figure BDA0003063362810000137
图5为膜组反渗透效率的对比图，在第一次变工况，态反馈控制控制和RBF神经网络自适应控制的系统分别在20.6s和1.1s达到最优反渗透效率，最优反渗透效率为0.2155；在第二次变工况，态反馈控制控制和RBF神经网络自适应控制的系统分别在22.2s和1.7s达到最优反渗透效率，最优反渗透效率为0.2136。可见RBF神经网络自适应控制幅度提高了反渗透效率达到最优的速度，且系统运行稳定性强。

Claims (2)

1.基于双RBF神经网络的反渗透膜组压力优化控制方法，其特征在于：包括变工况下优化压力的获取、基于RBF神经网络的优化压力在线调整、采用RBF神经网络自适应补偿的优化压力控制；所述优化压力的获取以变海水盐度下额定优化压力为单膜优化压力初始值，以各段膜压力综合最优为目标构建优化目标函数，采用拉格朗日乘数法获取反渗透膜组首段膜系统压力优化值；所述优化压力在线调整以膜组反渗透效率为性能评价指标，采用RBF神经网络令实际反渗透效率逼近最优反渗透效率，在线调整首段膜系统压力优化值；所述RBF神经网络自适应补偿控制采用基于状态反馈的主压力跟踪控制器和基于RBF神经网络自适应的干扰补偿控制器；所述主压力跟踪控制器采用状态反馈控制器，以优化压力跟踪误差及其导数作为状态反馈控制输入；所述干扰补偿控制器采用RBF神经网络自适应控制，以优化压力跟踪误差及其导数为控制输入，自适应地调整神经网络权值，以输出的干扰估计值逼近系统中的未知扰动。
2.根据权利要求1所述的基于双RBF神经网络的反渗透膜组压力优化控制方法，包括以下步骤：
步骤1构建用于优化控制的首段膜系统压力动态模型及其它各段膜压力稳态模型、用于优化性能评价的各段膜产水流量稳态模型
Figure FDA0003063362800000011
其中，ρ为浓水区液体密度，Pi(i＝1，2，……，n)为各段膜的系统压力，ηs为粘度系数，Vvr浓水阀液体体积，evr为浓水阀阻，vr为浓水流速，Qf为进水流量，λ为沿程阻力系数，l为膜元件轴向长度，d为浓水区管道当量直径，Am为膜面积，Ap为水管截面积，
Figure FDA0003063362800000012
为因膜组弹性造成的浓水区平均流量变化量，CH为膜组液容，Qpi(i＝1,2,3,…，n)为各段膜产水流量，Rvi(i＝1,2,3,4)为各段膜渗透液阻，Δπi(i＝1,2,3,4)为各段膜渗透液压差，Ki(i＝1,2,…,7)、C1为膜间关联系数；
步骤2反渗透膜组首段膜系统压力在线优化
第一步，以不同海水盐度下单膜额定压力PNref为输入优化目标函数的单膜优化压力初始值；
PNref＝0.02117S5-3.281S4+204.6S3-6637S2+241500S+653700 (2)
其中，S为海水盐度；
第二步，以膜组各段膜压力综合最优为目标构建优化目标函数；
min G＝γ1(P1-PNref)2+γ2(P2-PNref)2+…+γn(Pn-PNref)2 (3)
其中，γi(i＝1，2，……，n)为第i段膜的优化权值；
第三步，设定反渗透膜组压力优化的约束条件：
Figure FDA0003063362800000021
其中，Pf为膜组进水压力，Pdi(i＝1,2,3,…，n)表示n个反渗透膜的管压降，Qi(i＝1,2,3,…，n)为各段膜的轴向流量，Pr为浓水出口处的压力，Qr为浓水出口的流量，Psysmin、Psysmax为反渗透膜所允许的系统压力最大值和最小值，Smin、Smax为送入反渗透膜组海水的最小盐度和最大盐度；
第四步，采用拉格朗日乘数法获得反渗透膜组首段膜额定系统压力优化值：
Figure FDA0003063362800000022
第五步，根据实时采集的反渗透膜组产水流量计算反渗透膜滤饼污染层厚度：
Figure FDA0003063362800000023
其中，α为滤饼层累积系数，Am为反渗透膜面积，
Figure FDA0003063362800000024
为海水平均污染浓度；
第六步，根据实时获得的滤饼层厚度，更新反渗透膜渗透液阻，获得变海水盐度下的反渗透膜组反渗透效率和最大反渗透效率ηrM对应的首段膜系统压力P1ref new；
第七步，根据实时监测的反渗透膜组各段膜产水流量计算膜组的实际反渗透效率：
Figure FDA0003063362800000025
第八步，采用RBF神经网络实现反渗透膜组压力在线优化，计算最大反渗透效率逼近误差eη＝ηrM-ηr，以eηr和
Figure FDA0003063362800000029
为神经网络输入，即
Figure FDA0003063362800000026
第九步，计算RBF神经网络的高斯基函数：
Figure FDA0003063362800000027
其中，cj为高斯基函数中心，bj为高斯基函数宽度；
第十步，得到RBF神经网络输出的优化压力补偿值：
Figure FDA0003063362800000028
其中，ωj为输出层权值，hj为高斯基函数；
第十一步，将额定压力优化值与优化压力补偿值相加，获得首段膜系统压力优化值：
P1ref＝P1Nref+ΔP1ref (10)
第十二步，设定RBF神经网络权值在线优化最小目标函数：
Figure FDA0003063362800000031
第十三步，根据梯度下降法，在线调整RBF神经网络权值：
Figure FDA0003063362800000032
wj(N)＝wj(N-1)+Δwj(N)+k(wj(N-1)-wj(N-2)) (13)
其中，Δwj(N)为RBF神经网络权值的修正值，η为学习率，满足η∈[0,1]；
步骤3构建包含首段膜系统压力动态模型与电磁阀线圈电流模型的四膜反渗透膜组模型
Figure FDA0003063362800000033
其中，evrmax为最大阀阻，imax为电磁阀最大励磁电流，
Figure FDA0003063362800000034
为膜组弹性形变扰动，ΔRi为电磁阀绕组温升扰动，上式中，第四段膜的系统压力P4、膜组浓水流速vr以及函数l(P1)形如
Figure FDA0003063362800000035
上述模型中的Ki(i＝1,2,3,…，23)、C1、C2、C3均为膜间关联系数；
步骤4反渗透膜组海水淡化系统模型转化
第一步，将如式(14)所示的反渗透膜组模型简化如下：
Figure FDA0003063362800000036
其中
Figure FDA0003063362800000037
Figure FDA0003063362800000038
第二步，求取反渗透膜组海水淡化系统平衡点
忽略电流模型中的扰动，令电流的变化率为零，即求解式(19)，得到的电流即为平衡点处的电流i0：
Figure FDA0003063362800000039
其中，uin0为电磁阀外电路输入电压uin的稳态值；
将求得的系统平衡点电流i0带入式(20)计算平衡点的压力，得到平衡点处的压力P1,0，
Figure FDA0003063362800000041
由上式求得的压力即为平衡点的压力P1,0，该系统的平衡点即为(P1,0,i0)；
第三步，在反渗透膜组系统平衡点处得到线性化的反渗透膜组模型为：
Figure FDA0003063362800000042
其中，干扰项归结如下：
Figure FDA0003063362800000043
步骤5主压力跟踪控制器设计
第一步，以P1为状态变量x1，以
Figure FDA00030633628000000411
为状态变量x2，将模型中的扰动项统一归结为干扰项fxp，写出反渗透膜组淡化系统的状态空间方程：
Figure FDA0003063362800000044
其中，扰动项为
Figure FDA0003063362800000045
第二步，计算压力跟踪误差e＝P1ref-P1＝P1ref-x1，
Figure FDA0003063362800000046
定义虚拟变量
Figure FDA0003063362800000047
写出系统压力跟踪误差增广模型：
Figure FDA0003063362800000048
其中，虚拟控制输入为
Figure FDA0003063362800000049
第三步，对式(25)进行极点配置，采用状态反馈法设计主压力跟踪控制器，写出反渗透膜组淡化系统闭环特征多项式：
Figure FDA00030633628000000410
其中，KE＝[-kp -kd]为状态反馈增益矩阵；
第四步，通过选取合适的闭环极点，得到反馈增益矩阵KE的参数值，则此时的控制输入为：
Figure FDA0003063362800000051
步骤6基于RBF神经网络自适应的干扰补偿控制器设计
第一步，在状态反馈的基础上引入RBF神经网络自适应补偿，设计控制律为：
Figure FDA0003063362800000052
第二步，将(29)代入式(25)，得到状态空间方程：
Figure FDA0003063362800000053
第三步，定义最优权值为
Figure FDA0003063362800000054
其中，Ω表示具有希望边界的子集，即Ω＝{W|||f||≤M}，M为设计参数；
第四步，以RBF神经网络输出的自适应补偿
Figure FDA0003063362800000055
逼近干扰项fxp
Figure FDA0003063362800000056
第五步，定义模型逼近误差为
Figure FDA0003063362800000057
第六步，将式(33)代入式(30)，将状态空间方程表示为：
Figure FDA0003063362800000058
其中，
Figure FDA0003063362800000059
第七步，将式(32)代入式(34)，将闭环系统状态空间方程改写为：
Figure FDA00030633628000000510
第八步，设计Lyapunov函数为
Figure FDA00030633628000000511
其中，γ是一个正实数，P为对称正定矩阵，且满足Lyapunov方程ΛTP+PΛ＝-Q，Q为对称正定矩阵；
第九步，取
Figure FDA00030633628000000512
对V1求导，可得
Figure FDA0003063362800000061
其中，
Figure FDA0003063362800000062
则有
Figure FDA0003063362800000063
第十步，对V2求导，可得
Figure FDA0003063362800000064
第十一步，计算V的导数：
Figure FDA0003063362800000065
第十二步，为确保
Figure FDA0003063362800000066
设计权值自适应律：
Figure FDA0003063362800000067
其中，γ是一个正实数，P为对称正定矩阵，且满足ΛTP+PΛ＝-Q，Q为对称正定矩阵。
CN202110519454.XA 2021-05-13 2021-05-13 基于双rbf神经网络的反渗透膜组压力优化控制方法 Active CN113110065B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110519454.XA CN113110065B (zh) 2021-05-13 2021-05-13 基于双rbf神经网络的反渗透膜组压力优化控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110519454.XA CN113110065B (zh) 2021-05-13 2021-05-13 基于双rbf神经网络的反渗透膜组压力优化控制方法

Publications (2)

Publication Number Publication Date
CN113110065A true CN113110065A (zh) 2021-07-13
CN113110065B CN113110065B (zh) 2023-08-01

Family

ID=76722026

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110519454.XA Active CN113110065B (zh) 2021-05-13 2021-05-13 基于双rbf神经网络的反渗透膜组压力优化控制方法

Country Status (1)

Country Link
CN (1) CN113110065B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001062255A (ja) * 1999-08-27 2001-03-13 Toray Ind Inc 逆浸透膜プラントおよびその製造、運転方法ならびにそのための記憶媒体
KR20130022132A (ko) * 2011-08-25 2013-03-06 고려대학교 산학협력단 분산분석 공정을 이용한 해수담수화용 역삼투막의 성능분석 방법 및 이를 수행하는 분석장치
CN108549229A (zh) * 2018-04-20 2018-09-18 曲阜师范大学 一种桥式吊车神经网络自适应控制器及其设计方法
CN109081399A (zh) * 2018-08-01 2018-12-25 江苏凯纳水处理技术有限公司 反渗透系统的预测及诊断方法及其装置
KR20190113421A (ko) * 2018-03-28 2019-10-08 광주과학기술원 역삼투막 오염 예측 모델의 표본 입력 데이터 양 결정 방법 및 이를 이용한 장치
CN111173678A (zh) * 2020-01-09 2020-05-19 曲阜师范大学 一种rbf神经网络俯仰干扰补偿的风力机舱悬浮控制方法
CN111709132A (zh) * 2020-06-10 2020-09-25 杭州电子科技大学 一种全流程反渗透海水淡化装置及其优化方法
CN112464471A (zh) * 2020-11-25 2021-03-09 国网辽宁省电力有限公司 一种反渗透海水淡化系统的建模方法
CN112591852A (zh) * 2020-12-04 2021-04-02 河北工业大学 一种分流部分二级卷式膜元件反渗透海水淡化系统的操作优化方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001062255A (ja) * 1999-08-27 2001-03-13 Toray Ind Inc 逆浸透膜プラントおよびその製造、運転方法ならびにそのための記憶媒体
KR20130022132A (ko) * 2011-08-25 2013-03-06 고려대학교 산학협력단 분산분석 공정을 이용한 해수담수화용 역삼투막의 성능분석 방법 및 이를 수행하는 분석장치
KR20190113421A (ko) * 2018-03-28 2019-10-08 광주과학기술원 역삼투막 오염 예측 모델의 표본 입력 데이터 양 결정 방법 및 이를 이용한 장치
CN108549229A (zh) * 2018-04-20 2018-09-18 曲阜师范大学 一种桥式吊车神经网络自适应控制器及其设计方法
CN109081399A (zh) * 2018-08-01 2018-12-25 江苏凯纳水处理技术有限公司 反渗透系统的预测及诊断方法及其装置
CN111173678A (zh) * 2020-01-09 2020-05-19 曲阜师范大学 一种rbf神经网络俯仰干扰补偿的风力机舱悬浮控制方法
CN111709132A (zh) * 2020-06-10 2020-09-25 杭州电子科技大学 一种全流程反渗透海水淡化装置及其优化方法
CN112464471A (zh) * 2020-11-25 2021-03-09 国网辽宁省电力有限公司 一种反渗透海水淡化系统的建模方法
CN112591852A (zh) * 2020-12-04 2021-04-02 河北工业大学 一种分流部分二级卷式膜元件反渗透海水淡化系统的操作优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张彪;邢健峰;纪志成;: "基于优化SVM的反渗透脱盐水故障诊断", 系统仿真学报, no. 05, pages 141 - 147 *

Also Published As

Publication number Publication date
CN113110065B (zh) 2023-08-01

Similar Documents

Publication Publication Date Title
CN111474965B (zh) 基于模糊神经网络的串联输水渠道水位预测与控制方法
CN111812968B (zh) 基于模糊神经网络pid控制器的阀位串级控制方法
CN105425612A (zh) 一种水轮机调节系统控制参数的优选方法
CN101598109B (zh) 一种风力发电机偏航系统的智能控制方法
CN104808494B (zh) 基于自适应蚁群遗传混合算法pid参数整定方法
CN204595644U (zh) 基于神经网络的铝棒加热炉燃烧温度自动控制装置
CN106054596B (zh) 一种基于改进性能指标的pid控制器参数优化整定方法
CN110008575B (zh) 循环冷却水系统工艺介质多温度目标设定值切换多参数预测控制算法
CN108181802A (zh) 一种性能可控pid控制器参数优化整定方法
CN103744443A (zh) 一种基于蓄量动态调节的大型明渠水位自动控制方法
CN102411308A (zh) 基于递归神经网络模型的溶解氧的自适应控制方法
CN105388764A (zh) 基于动态矩阵前馈预测的电液伺服pid控制方法及系统
CN109473182B (zh) 一种基于深度信念网络的mbr膜透水率智能检测方法
CN106485064A (zh) 一种抽水蓄能机组水轮机工况智能开机方法
KR101187416B1 (ko) 예측 모델을 이용한 해수담수화 방법 및 이를 수행하는 장치
CN113267314A (zh) 一种暂冲式风洞的超声速流场总压控制系统
CN113110065A (zh) 基于双rbf神经网络的反渗透膜组压力优化控制方法
CN113213585B (zh) 基于电液相似原理的多膜组建模及性能分析方法
CN108549207A (zh) 一种水轮机调节系统控制参数的方法
CN104123348A (zh) 一种基于arx的供水管网在线建模方法
CN112856014B (zh) 一种智能阀门电动执行机构的控制系统构建方法
CN113219833B (zh) 基于优化压力自适应补偿的反渗透膜组优化控制方法
Mahmud Neural network based PID control analysis
CN102139769B (zh) 基于自组织cmac的挠性卫星快速稳定控制方法
CN112766565B (zh) 梯级水库群防洪优化调度分数阶逐步优化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant