CN112678904B - 一种多级闪蒸海水淡化***结垢清洗装置及方法 - Google Patents

一种多级闪蒸海水淡化***结垢清洗装置及方法 Download PDF

Info

Publication number
CN112678904B
CN112678904B CN202011216530.1A CN202011216530A CN112678904B CN 112678904 B CN112678904 B CN 112678904B CN 202011216530 A CN202011216530 A CN 202011216530A CN 112678904 B CN112678904 B CN 112678904B
Authority
CN
China
Prior art keywords
brine
flash
module
flash chamber
stage
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.)
Active
Application number
CN202011216530.1A
Other languages
English (en)
Other versions
CN112678904A (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.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi 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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN202011216530.1A priority Critical patent/CN112678904B/zh
Publication of CN112678904A publication Critical patent/CN112678904A/zh
Application granted granted Critical
Publication of CN112678904B publication Critical patent/CN112678904B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

Landscapes

  • Heat Treatment Of Water, Waste Water Or Sewage (AREA)
  • Vaporization, Distillation, Condensation, Sublimation, And Cold Traps (AREA)

Abstract

本发明涉及一种多级闪蒸海水淡化***结垢清洗装置及方法，通过人机交互模块给定多级闪蒸海水淡化***的结构参数、必要系数和数据采集周期；利用数据采集模块采集、通过A/D和D/A转换模块将上述参数发送给中央处理模块；中央处理模块根据内部存储的问题模型，调用模拟计算模块，将上述计算后参数发送到参数决策模块。参数决策模块估算得出当前***的结垢程度以及***下次清洗时间，并且模拟计算模块和参数决策模块将其计算数据发送至中央处理模块，中央处理模块将数据进行整理发送至显示模块。进行实时数据采集，并由模型进行计算，得出***结垢系数从而判断***的结垢情况，指导工厂能够及时进行相应的结垢处理，保持***的优良性能。

Description

一种多级闪蒸海水淡化***结垢清洗装置及方法
技术领域
本发明属于化工生产过程控制技术领域，尤其涉及一种多级闪蒸海水淡化***结垢清洗装置及方法。
背景技术
多级闪蒸是海水淡化工业中最成熟，运行安全性最高，弹性大，适合于大型和超大型淡化装置的一项技术。在未来的海水淡化领域中仍将继续发挥重要的作用。
MSF海水淡化装置有两种经典结构，分别为直通式多级闪蒸(One ThroughMultistage Flash,OT-MSF)海水淡化装置和盐水循环式多级闪蒸(Brine RecirculationMultistage Flash,BR-MSF)海水淡化装置。其中OT-MSF装置结构较简单，而BR-MSF因为其具有更好的综合效果而得到更加普遍的应用。BR-MSF装置由盐水加热器、热回收段、热排放段以及混合分离模块四大部分组成。多级闪蒸海水淡化***经典结构流程如图1所示。原海水首先进入热排放段，自身被预热的同时冷凝闪蒸室中产生的蒸汽。进料海水经过热排放段的预热后，被分流为两部分，一部分返回大海，一部分与循环盐水混合在一起，然后被泵入热回收段的末端，当它从右到左流过一系列的热交换器时，自身被逐渐加热，同时闪蒸室中闪蒸出来的蒸汽得到冷凝。最后，海水从热回收段的第一级闪蒸室内的换热器内出来，流入盐水加热器被进一步加热，并以热盐水的形式流入第一级闪蒸室。由于该闪蒸室中的压力控制在低于热盐水温度所对应的饱和蒸气压，故热盐水进入闪蒸室后即成为过热水而急速地部分气化，所产生的蒸汽遇到冷凝管，被冷凝后滴入淡水托盘。如此重复，直到最后一级，浓盐水被排放，淡水被抽出。
多级闪蒸过程中的循环盐水经过预热器加热，由于海水中的盐度、硬度、总固溶物及其他杂质的含量均较高,某些组分的溶解度可能达到过饱和，导致在多级闪蒸海水淡化装置的换热面上易结垢，这已成为热法海水淡化面临的主要问题之一。污垢的存在不仅增大流体阻力，增加能量消耗，降低传热效率，严重时还会堵塞管路，使设备瘫痪。但多级闪蒸海水淡化***性能在运行过程中不仅受到污垢参数影响，还受到进料海水流量、盐度、循环浓盐水量、盐水最高温度、盐水加热器参数等影响，导致产水性能多变。***在运行过程中，其结垢特性无法直接测量，其清洗保养一般都是基于经验而进行操作，缺乏科学依据。
本发明可根据***实时运行的输入输出运行数据，通过模型计算出当前***的结垢系数的大小，从而得出***结垢程度，并给出相应的***维护和结垢处理方法。
发明内容
本发明的目的是通过在线估测结垢特征的结垢系数，给出清洗操作指导。
传统的多级闪蒸海水淡化装置的除垢依据，是产水量和造水比远低于标准值或***性能出现严重问题，其对装置寿命周期可能会产生不可逆的危害，对设备再次投入生产也会产生考验。本发明对热法海水淡化***结垢特性参数的估计和维保决策，涉及在数理结合建模下通过对传热系数和产水性能的计算，对结垢参数进行软测量并指导***清洗操作。进行实时数据采集，并由模型进行计算，得出***结垢系数从而判断***的结垢情况，指导工厂能够及时进行相应的结垢处理，保持***的优良性能。
一种多级闪蒸海水淡化***结垢清洗装置包括传感器模块、数据采集模块、人机交互模块、模拟计算模块、A/D和D/A转换模块、中央处理模块、参数决策模块、显示模块；传感器模块包括五路流量传感器、四路温度传感器、两路海水含盐量传感器；数据采集模块用于采集传感器模块获得的进料海水温度和含盐量、淡水流量、循环盐水和热排放海水流量，以及加热蒸汽流量和温度；人机交互模块用于设定多级闪蒸海水淡化***的各级闪蒸室结构参数、盐水加热器参数、设定数据采集周期；A/D和D/A转换模块用于将接收到的模拟量转换为相应的数字量或将接收到的数字量转换为相应的模拟量；中央处理模块用于存储多级闪蒸海水淡化***结垢系数处理问题模型和所需的物性参数，通过A/D和D/A转换模块接收并存储数据采集模块采集的数据以及接收并存储模拟计算模块和参数决策模块处理得到的数据，然后将数据传输到显示模块；模拟计算模块根据当前采集到的数据，调用***程序进行计算，得出当前盐水加热器的结垢系数fBH、闪蒸室的结垢系数fj((j＝1,2,…,n)；n为***闪蒸室级数)、造水比GOR、日运行成本TOC，并将上述参数传至显示模块；参数决策模块根据模拟计算模块的之前所存储的数据进行分析，估计出***结垢程度，和下次***清洗时间，并将上述参数传至显示模块；显示模块显示模拟计算模块和参数决策模块的计算结果。
一种多级闪蒸海水淡化***结垢清洗方法包括如下步骤：首先由工程师通过人机交互模块设定多级闪蒸海水淡化***的结构参数，给定运行周期。此后该装置以给定的周期采集所需数据，并通过A/D和D/A转换模块发给中央处理模块，中央处理模块调用模拟计算模块。模拟计算模块通过计算中央处理模块中的模型问题，得出当前***盐水加热器的结垢系数fBH和闪蒸室的结垢系数fj、造水比GOR、日运行成本TOC。参数决策模块根据模拟计算模块所得参数，进行比较分析，得出当前***的结垢程度和下次***清洗时间并将相关结论发送至***显示模块。该装置在下一个运行周期，重复数据采集、模拟计算、比较分析过程，不断地实时估测***结垢程度。
具体包括如下步骤：
步骤A1：操作员或者工程师通过人机交互模块给定多级闪蒸海水淡化***的结构参数、物性系数和数据采集周期Tc；
步骤A2：利用数据采集模块采集当前时间点的进入热排放段的冷却盐水流量WF、循环海水流量WRe、加热蒸汽流量Wsteam、闪蒸盐水流量WB、闪蒸淡水流量WD、进料海水温度Tsea、加热蒸汽温度Tsteam、闪蒸盐水温度TB、闪蒸淡水温度TD、进料盐水浓度CF、闪蒸盐水浓度CB。记录当前时间，令变量T1等于当前时间，然后通过A/D和D/A转换模块将上述参数发送给中央处理模块；
步骤A3：中央处理模块根据内部存储的多级闪蒸海水淡化装置稳态模拟和日运行成本模型，调用模拟计算模块，从而计算出在当前T1时刻盐水加热器的结垢系数fBH、闪蒸室的结垢系数fj、造水比GOR和日运行成本TOC，然后将上述参数发送到参数决策模块。
步骤A4：参数决策模块将此周期得到的参数集L1(fBH,fj,GOR,TOC)与之前各周期得到的参数集Li(fBH,fj,GOR,TOC)((i＝1,2,…,m)；m为当前周期数)进行比较分析，估算得出当前***的结垢程度以及***下次清洗时间，并且模拟计算模块和参数决策模块将其计算数据发送至中央处理模块。
步骤A5：中央处理模块将数据进行整理发送至显示模块，显示模块显示出当前时间T1、盐水加热器的结垢系数fBH、闪蒸室的结垢系数fj、造水比GOR、日运行成本TOC、***结垢程度、***待清洗时间Tim。
步骤A6：记录当前时间为T2，如果T2-T1<Tc，则继续等待；否则转步骤A2，重新进行数据采集。
所述的中央处理模块中的多级闪蒸海水淡化装置稳态模拟和日运行成本的模型由式(1)～(38)组成；
多级闪蒸海水淡化装置的稳态模型由闪蒸室方程、盐水加热器方程、混合分离方程、物性参数方程组成。对于第j级闪蒸室，闪蒸室模块模型由下列式(1)～(8)组成：
WBj-1+WDj-1＝WBj+WDj (1)；
WBj-1CBj-1＝WBjCBj (2)；
WBj-1hBj-1＝WBjhBj+VBjhVj (3)；
WBj-1-WBj＝VBj (4)；
Figure BDA0002760582980000041
Figure BDA0002760582980000042
其中，WBj表示第j级的闪蒸盐水质量流量，WBj-1表示第j-1级的闪蒸盐水质量流量，WDj表示第j级产出的淡水质量流量，WDj-1表示第j-1级的闪蒸盐水质量流量。CBj表示第j级的闪蒸盐水浓度，CBj-1表示第j-1级的闪蒸盐水浓度。hBj表示第j级闪蒸盐水的比焓，hBj-1表示第j-1级闪蒸盐水的比焓，hVj表示第j级闪蒸蒸汽的比焓，VBj表示第j级闪蒸室中盐水蒸发量。WF表示进入热排放段的冷却盐水的质量流量，CPRj表示离开第j级闪蒸室的冷却盐水热容，CPDj表示离开第j级闪蒸室的淡水比热容，CPDj-1表示离开第j-1级闪蒸室的淡水比热容，CPBj表示离开第j级闪蒸室的闪蒸盐水比热容，CPBj-1表示离开第j-1级闪蒸室的闪蒸盐水比热容，TFj表示离开第j级闪蒸室的冷却盐水温度，TFj+1表示进入第j级闪蒸室的冷却盐水温度，TBj表示离开第j级闪蒸室的盐水温度，TBj-1表示离开第j-1级闪蒸室的盐水温度，TDj表示离开第j级闪蒸室的淡水温度，TDj-1表示离开第j-1级闪蒸室的淡水温度，T*表示理想状态下闪蒸参考温度。Aj表示第j级闪蒸室的传热面积，
Figure BDA0002760582980000043
表示第j级闪蒸室的传热系数，其中W指的是第j级闪蒸室的进料流量，对于热排放段来说W＝WF，表示进料海水流量，而对于热回收段来说W＝WR，表示热回收段的进料流量，
Figure BDA0002760582980000044
表示闪各级蒸室外直径，
Figure BDA0002760582980000045
表示各级闪蒸室冷凝管内直径，Di表示内直径，其在物性方程部分有具体公式展开。
TBj＝TDj+ΔBPEj+ΔNETDj+ΔTLj (7)；
TVj＝TDj+ΔTLj (8)；
其中，ΔBPEj表示第j级盐水沸点升高，ΔNETDj表示非平衡余量，ΔTLj表示经过除雾器和冷凝器的温度损失。TVj表示第j级闪蒸室的闪蒸蒸汽温度。
盐水加热器模块模型由式(9)～(12)组成：
WB0＝WR (9)；
CB0＝CR (10)；
WRCPRH(TB0-TF1)＝Wsteamλs (11)；
其中，
Figure BDA0002760582980000051
Figure BDA0002760582980000052
其中，WB0表示离开盐水加热器的闪蒸盐水质量流量，WR表示进入热回收段的冷却盐水质量流量。CB0表示进入盐水加热器的闪蒸盐水浓度，CR表示进入热回收段的冷却盐水浓度。CPRH表示进入盐水加热器的冷却盐水热容，TB0表示离开盐水加热器的闪蒸盐水温度，Wsteam表示为加热蒸汽质量流量，λs表示蒸汽潜热，Tsteam表示为盐水加热器蒸汽温度。AH表示盐水加热器的传热面积，
Figure BDA0002760582980000053
表示盐水加热器的传热系数，其中W指的是盐水加热器的进料流量WR，
Figure BDA0002760582980000054
表示盐水加热器外直径，
Figure BDA0002760582980000055
表示盐水加热器冷凝管内直径，其在物性方程部分有具体公式展开。
混合分离模块模型由式(13)～(20)组成：
WBD＝WBN-WRe (13)；
Wm＝WF-Wr (14)；
S＝Wr-Cw (15)；
WR＝WRe+Wm (16)；
WR·CR＝WRe·CRe+Wm·Cm (17)；
WR·hR＝WRe·hRe+Wm·hm (18)；
WF＝S+WS (19)；
Figure BDA0002760582980000056
其中WBD表示废弃海水质量流量，WBN表示离开最后一级闪蒸室的盐水质量流量，WRe表示循环盐水质量流量。Wm表示补充盐水质量流量，Wr表示热排放海水质量流量。S表示热排放段回流流量，Cw表示热排放段排出盐水的质量流量。CRe表示循环盐水浓度，Cm表示补充盐水浓度。hR表示进入热回收段的盐水比焓，hRe表示循环盐水比焓，hm表示补充盐水比焓。WS表示进料海水质量流量。
Figure BDA0002760582980000057
表示热排放段入口盐水比焓，hS表示热排放段回流盐水比焓，
Figure BDA0002760582980000058
表示进料海水比焓。
物性参数方程模型由式(21)～(31)组成：
Figure BDA0002760582980000059
CPB＝[1-CB(0.011311-1.146×10-5TB)]×CPD (22)；
hV＝596.912+0.46694TS-0.000460256TS 2 (23)；
hB＝CPB·TB (24)；
hD＝CPD·TD (25)；
Figure BDA0002760582980000061
其中，CC＝(19.819CB)/(1-CB)，CC表示浓度转换，CB表示闪蒸盐水浓度。
Figure BDA0002760582980000062
其中，ωj＝WF/wj，ωj表示第j级进入热排放段的冷却盐水单位长度上的质量流量，wj表示第j级闪蒸室冷凝管的宽度。
TL＝exp(1.885-0.02063TD)/1.8 (28)；
U＝4.8857/(y+z+4.8857f) (29)；
y＝[0.0013(v×Di)0.2]/[(0.2018+0.0031×T)v] (30)；
Figure BDA0002760582980000063
GOR＝WDN/Wsteam (32)；
其中CPD表示闪蒸室的淡水比热容。CPB表示闪蒸室的盐水比热容，CB表示闪蒸盐水浓度，TB表示闪蒸盐水温度。hv表示蒸汽焓值，TS为蒸汽温度。hB表示盐水焓值。hD表示淡水焓值，TD表示淡水温度。BPE表示盐水沸点升高，其中，T表示处于闪蒸室或者盐水加热器的盐水温度。NETD表示非平衡温差，Hj表示闪蒸盐水液位，wj表示第j级闪蒸室冷凝管的宽度，ΔTB表示两级间的盐水温差。TL表示经过除雾器和冷凝管的温度损失。U表示盐水加热器或者各级闪蒸室的传热系数。y为有关管内流速、盐水温度和冷凝管内直径的中间表达式，z为有关闪蒸室淡水温度的中间表达式，y和z仅为简化计算，f为盐水加热器或者闪蒸室的结垢系数。v表示管内流速，Di表示冷凝管内直径，T表示在换热器出口处的盐水温度。WDN表示闪蒸室总淡水的质量流量，GOR表示造水比(Gained Output Ratio)是反应***性能的一个关键指标。
日运行成本模块模型由式(33)～(38)组成：
C1＝22×Wsteam×[(Tsteam-40)/85]×0.00415 (33)；
C2＝22×(DN/ρw)×0.109 (34)；
C3＝22×(DN/ρw)×0.082 (35)；
C4＝22×(DN/ρb)×0.024 (36)；
C5＝22×(DN/ρw)×0.1 (37)；
TOC＝C1+C2+C3+C4+C5 (38)；
其中，C1表示加热蒸汽费用，C2表示装置消耗的电能费用，C3表示维护与空闲费用，C4表示为预处理费用，C5表示为人工费用，TOC表示***日运行费用。DN表示总淡水产量，ρw表示纯水密度，ρb表示浓盐水密度。
模拟计算模块将计算式(1)～(38)构成的稳态模拟问题运用拟牛顿法进行求解，确定当前***盐水加热器结垢系数fBH、热回收段和热排放段闪蒸室结垢系数fj、造水比GOR和日运行成本TOC，其具体步骤如下所示：
步骤B1：将上式(1)～(38)中有待求未知量的方程整理为非线性方程组，其具体表现形式如式(39)所示：
Figure BDA0002760582980000071
向量的表示形式为：
G(x)＝0 (40)；
这里，x＝(x1,x2,…,xp)T，G＝(g1,g2,…,gq)T，gi(i＝1,2,…,q):Rn→R。
其中，x1,x2,…,xp表示非线性方程组的p个未知数，x表示方程组的未知数构成的p×1维的矩阵。g1,g2,…,gq表示构成非线性方程组的q个等式，G表示非线性方程组的等式构成的q×p维的矩阵。
步骤B2：令初始迭代次数k＝0，设置初始值x0∈Rn，初始拟牛顿矩阵Bk为p×p维单位矩阵，函数值计算精度为ε，最小步长精度为δ。
步骤B3：计算迭代次数k下搜索方向向量dk，以得到下一个点的值：
xk+1＝xk+dk (41)；
dk＝-BkΔG(xk) (42)；
步骤B4：计算G(xk)和G(xk+1)，若||G(xk+1)||＜ε或者||xk+1-xk||＜δ，则终止迭代，得出非线性方程组的解及盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR、日运行成本TOC；否则，转至步骤B5。
步骤B5：通过修正Bk得到Bk+1：
Figure BDA0002760582980000081
其中，sk＝xk+1-xk，yk＝ΔG(xk+1)-ΔG(xk)。xk表示第k次迭代所求值，xk+1表示第k+1次迭代所求值，sk表示两次所求值之差。ΔG(xk+1)表示G(xk+1)的梯度向量，ΔG(xk)表示G(xk)的梯度向量，yk表达第k+1次迭代和k次迭代的梯度向量差。
步骤B6：令k:＝k+1，转至步骤B3进行计算。
参数决策模块根据模拟计算模块所得出的所有时间点的盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR和日运行费用TOC，采用最小二乘法进行拟合并进行分析，对***结垢程度进行合理的判断，具体步骤如下所示：
步骤C1：继承模拟计算模块从开始进行数据采集时的所有数据，并且针对于盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR和日运行费用TOC，对盐水加热器结垢系数fBH设定上限H1，闪蒸室结垢系数fj设定上限H2、造水比GOR设定下限X1、日运行费用TOC设定上限X2。
步骤C2：将上述数据分成四部分分别进行拟合，首先对盐水加热器结垢系数fBH采用最小二乘法进行曲线拟合。假设给定采集的数据点(xxi,yyi)(i＝0,1,…,m)，xxi表示数据采集时间点Ti，yyi表示每个时间点采集的盐水加热器结垢系数fBH(i)。拟合多项式为
Figure BDA0002760582980000082
为了保证模块计算速度令多项式阶数nn＝3，使得
Figure BDA0002760582980000083
其中I为a0,a1,a2,a3的四元函数，因此上述问题即为求I＝I(a0,a1,a2,a3)的极值问题。
步骤C3：由多元函数求极值的必要条件，得
Figure BDA0002760582980000091
即
Figure BDA0002760582980000092
步骤C4：将式(46)改写成正规方程组形式
Figure BDA0002760582980000093
可以证明正规方程组(47)的系数矩阵是一个对称正定矩阵，故存在唯一解。从式(47)中解得akk(kk＝0,1,…,3)，从而可得多项式
Figure BDA0002760582980000094
步骤C5：重复步骤C2至C4分别得出盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR和日运行费用TOC关于***运行时间的拟合曲线，为P1(t)、P2(t)、P3(t)和P4(t)，其中t代表时间。
步骤C6：根据参数所设上下限，和所求得的四条拟合曲线，判断当前***是否需要清洗，并根据拟合方程推出上述4个指标下的***待清洗时间点Tim1、Tim2、Tim3、Tim4。
步骤C7：将Tim1、Tim2、Tim3、Tim4进行比较得出最近的时间节点记为Tim，该时间点为***下次待清洗时间点。
步骤C8：对盐水加热器的结垢系数fBH、闪蒸室结垢系数fj的大小进行划分，分别对应不同等级的结垢程度，依据当前数据采集点的数值大小可得出***结垢程度。
附图说明
图1为本发明针对的多级闪蒸海水淡化***结构图；
图2本发明的整体结构示意图；
图3为本发明的运行流程图；
图4为模拟计算模块流程图；
图5为参数决策模块流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的分析。
如图2所示，多级闪蒸海水淡化***结垢处理装置，包括传感器模块、数据采集模块、人机交互模块、模拟计算模块、A/D和D/A转换模块、中央处理模块、参数决策模块、显示模块；传感器模块包括五路流量传感器、四路温度传感器、两路海水含盐量传感器；数据采集模块用于采集传感器模块获得的进料海水温度和含盐量、淡水流量、循环盐水和热排放海水流量，以及加热蒸汽流量和温度；人机交互模块用于设定多级闪蒸海水淡化***的各级闪蒸室结构参数、盐水加热器参数、设定数据采集周期；A/D和D/A转换模块用于将接收到的模拟量转换为相应的数字量或将接收到的数字量转换为相应的模拟量；中央处理模块用于存储多级闪蒸海水淡化***结垢系数处理问题模型和所需的物性参数，通过A/D和D/A转换模块接收并存储数据采集模块采集的数据以及接收并存储模拟计算模块和参数决策模块处理得到的数据，然后将数据传输到显示模块；模拟计算模块根据当前采集到的数据，调用***程序进行计算，得出当前盐水加热器的结垢系数fBH、闪蒸室的结垢系数fj((j＝1,2,…,n)；n为***闪蒸室级数)、造水比GOR、日运行成本TOC，并将上述参数传至显示模块；参数决策模块根据模拟计算模块的之前所存储的数据进行分析，估计出***结垢程度，和下次***清洗时间，并将上述参数传至显示模块；显示模块显示模拟计算模块和参数决策模块的计算结果。
如图3所示，对于采用本发明的某多级闪蒸海水淡化***，为了在线估测***的结垢程度并做出相应处理，需要执行以下步骤：
步骤A1:用户或者工程师通过人机交互模块给定多级闪蒸海水淡化***的结构参数，包括热回收段级数NR＝16，热排放段级数NJ＝3；盐水加热器、闪蒸室尺寸等重要系数；设定数据采集周期Tc＝2小时；
步骤A2：利用数据采集模块采集当前时间点的进入热排放段的冷却盐水流量WF、循环海水流量WRe、加热蒸汽流量Wsteam、闪蒸盐水流量WB、闪蒸淡水流量WD、进料海水温度Tsea、加热蒸汽温度Tsteam、闪蒸盐水温度TB、闪蒸淡水温度TD、进料盐水浓度CF、闪蒸盐水浓度CB。记录当前时间，令变量T1等于当前时间，然后通过A/D和D/A转换模块将上述参数发送给中央处理模块；
步骤A3：中央处理模块根据内部存储的多级闪蒸海水淡化装置稳态模拟和日运行成本模型，调用模拟计算模块，从而计算出在当前T1时刻盐水加热器的结垢系数fBH、闪蒸室的结垢系数fj、造水比GOR和日运行成本TOC，然后将上述参数发送到参数决策模块。
步骤A4：参数决策模块将此周期得到的参数集L1(fBH,fj,GOR,TOC)与之前各周期得到的参数集Li(fBH,fj,GOR,TOC)((i＝1,2,…,m)；m为当前周期数)进行比较分析，估算得出当前***的结垢程度以及***下次清洗时间，并且模拟计算模块和参数决策模块将其计算数据发送至中央处理模块。
步骤A5：中央处理模块将数据进行整理发送至显示模块，显示模块显示出当前时间T1、盐水加热器的结垢系数fBH、闪蒸室的结垢系数fj、造水比GOR、日运行成本TOC、***结垢程度、***待清洗时间Tim。
步骤A6：记录当前时间为T2，如果T2-T1<Tc，则继续等待；否则转步骤A2，重新进行数据采集。
所述的中央处理模块中的多级闪蒸海水淡化装置稳态模拟和日运行成本的模型由式(1)～(38)组成：
多级闪蒸海水淡化装置的稳态模型由闪蒸室方程、盐水加热器方程、混合分离方程、物性参数方程组成。对于第j级闪蒸室，闪蒸室模块模型由下列式(1)～(8)组成
WBj-1+WDj-1＝WBj+WDj (1)；
WBj-1CBj-1＝WBjCBj (2)；
WBj-1hBj-1＝WBjhBj+VBjhVj (3)；
WBj-1-WBj＝VBj (4)；
Figure BDA0002760582980000111
Figure BDA0002760582980000112
其中，WBj表示第j级的闪蒸盐水质量流量，WBj-1表示第j-1级的闪蒸盐水质量流量，WDj表示第j级产出的淡水质量流量，WDj-1表示第j-1级的闪蒸盐水质量流量。CBj表示第j级的闪蒸盐水浓度，CBj-1表示第j-1级的闪蒸盐水浓度。hBj表示第j级闪蒸盐水的比焓，hBj-1表示第j-1级闪蒸盐水的比焓，hVj表示第j级闪蒸蒸汽的比焓，VBj表示第j级闪蒸室中盐水蒸发量。WF表示进入热排放段的冷却盐水的质量流量，CPRj表示离开第j级闪蒸室的冷却盐水热容，CPDj表示离开第j级闪蒸室的淡水比热容，CPDj-1表示离开第j-1级闪蒸室的淡水比热容，CPBj表示离开第j级闪蒸室的闪蒸盐水比热容，CPBj-1表示离开第j-1级闪蒸室的闪蒸盐水比热容，TFj表示离开第j级闪蒸室的冷却盐水温度，TFj+1表示进入第j级闪蒸室的冷却盐水温度，TBj表示离开第j级闪蒸室的盐水温度，TBj-1表示离开第j-1级闪蒸室的盐水温度，TDj表示离开第j级闪蒸室的淡水温度，TDj-1表示离开第j-1级闪蒸室的淡水温度，T*表示理想状态下闪蒸参考温度。Aj表示第j级闪蒸室的传热面积，
Figure BDA0002760582980000121
表示第j级闪蒸室的传热系数，其中W指的是第j级闪蒸室的进料流量，对于热排放段来说W＝WF，表示进料海水流量，而对于热回收段来说W＝WR，表示热回收段的进料流量，
Figure BDA0002760582980000122
表示闪各级蒸室外直径，
Figure BDA0002760582980000123
表示各级闪蒸室冷凝管内直径，Di表示内直径，其在物性方程部分有具体公式展开。
TBj＝TDj+ΔBPEj+ΔNETDj+ΔTLj (7)；
TVj＝TDj+ΔTLj (8)；
其中，ΔBPEj表示第j级盐水沸点升高，ΔNETDj表示非平衡余量，ΔTLj表示经过除雾器和冷凝器的温度损失。TVj表示第j级闪蒸室的闪蒸蒸汽温度。
盐水加热器模块模型由式(9)～(12)组成：
WB0＝WR (9)；
CB0＝CR (10)；
WRCPRH(TB0-TF1)＝Wsteamλs (11)；
其中，
Figure BDA0002760582980000124
Figure BDA0002760582980000125
其中WB0表示离开盐水加热器的闪蒸盐水质量流量，WR表示进入热回收段的冷却盐水质量流量。CB0表示进入盐水加热器的闪蒸盐水浓度，CR表示进入热回收段的冷却盐水浓度。CPRH表示进入盐水加热器的冷却盐水热容，TB0表示离开盐水加热器的闪蒸盐水温度，Wsteam表示为加热蒸汽质量流量，λs表示蒸汽潜热，Tsteam表示为盐水加热器蒸汽温度。AH表示盐水加热器的传热面积，
Figure BDA0002760582980000126
表示盐水加热器的传热系数，其中W指的是盐水加热器的进料流量WR，
Figure BDA0002760582980000127
表示盐水加热器外直径，
Figure BDA0002760582980000128
表示盐水加热器冷凝管内直径，其在物性方程部分有具体公式展开。
混合分离模块模型由式(13)～(20)组成：
WBD＝WBN-WRe (13)；
Wm＝WF-Wr (14)；
S＝Wr-Cw (15)；
WR＝WRe+Wm (16)；
WR·CR＝WRe·CRe+Wm·Cm (17)；
WR·hR＝WRe·hRe+Wm·hm (18)；
WF＝S+WS (19)；
Figure BDA0002760582980000131
其中WBD表示废弃海水质量流量，WBN表示离开最后一级闪蒸室的盐水质量流量，WRe表示循环盐水质量流量。Wm表示补充盐水质量流量，Wr表示热排放海水质量流量。S表示热排放段回流流量，Cw表示热排放段排出盐水的质量流量。CRe表示循环盐水浓度，Cm表示补充盐水浓度。hR表示进入热回收段的盐水比焓，hRe表示循环盐水比焓，hm表示补充盐水比焓。WS表示进料海水质量流量。
Figure BDA0002760582980000132
表示热排放段入口盐水比焓，hS表示热排放段回流盐水比焓，
Figure BDA0002760582980000133
表示进料海水比焓。
物性参数方程模型由式(21)～(31)组成：
Figure BDA0002760582980000134
CPB＝[1-CB(0.011311-1.146×10-5TB)]×CPD (22)；
hV＝596.912+0.46694TS-0.000460256TS 2 (23)；
hB＝CPB·TB (24)；
hD＝CPD·TD (25)；
Figure BDA0002760582980000135
其中，CC＝(19.819CB)/(1-CB)，CC表示浓度转换，CB表示闪蒸盐水浓度，。
Figure BDA0002760582980000136
其中，ωj＝WF/wj，ωj表示第j级进入热排放段的冷却盐水单位长度上的质量流量，wj表示第j级闪蒸室冷凝管的宽度。
TL＝exp(1.885-0.02063TD)/1.8 (28)；
U＝4.8857/(y+z+4.8857f) (29)；
y＝[0.0013(v×Di)0.2]/[(0.2018+0.0031×T)v] (30)；
Figure BDA0002760582980000141
GOR＝WDN/Wsteam (32)；
其中CPD表示闪蒸室的淡水比热容。CPB表示闪蒸室的盐水比热容，CB表示闪蒸盐水浓度，TB表示闪蒸盐水温度。hv表示蒸汽焓值，TS为蒸汽温度。hB表示盐水焓值。hD表示淡水焓值，TD表示淡水温度。BPE表示盐水沸点升高，其中，T表示处于闪蒸室或者盐水加热器的盐水温度。NETD表示非平衡温差，Hj表示闪蒸盐水液位，wj表示第j级闪蒸室冷凝管的宽度，ΔTB表示两级间的盐水温差。TL表示经过除雾器和冷凝管的温度损失。U表示盐水加热器或者各级闪蒸室的传热系数。y为有关管内流速、盐水温度和冷凝管内直径的中间表达式，z为有关闪蒸室淡水温度的中间表达式，y和z仅为简化计算，f为盐水加热器或者闪蒸室的结垢系数。v表示管内流速，Di表示冷凝管内直径，T表示在换热器出口处的盐水温度。WDN表示闪蒸室总淡水的质量流量，GOR表示造水比(Gained Output Ratio)是反应***性能的一个关键指标。
日运行成本模块模型由式(33)～(38)组成：
C1＝22×Wsteam×[(Tsteam-40)/85]×0.00415 (33)；
C2＝22×(DN/ρw)×0.109 (34)；
C3＝22×(DN/ρw)×0.082 (35)；
C4＝22×(DN/ρb)×0.024 (36)；
C5＝22×(DN/ρw)×0.1 (37)；
TOC＝C1+C2+C3+C4+C5 (38)；
其中，C1表示加热蒸汽费用，C2表示装置消耗的电能费用，C3表示维护与空闲费用，C4表示为预处理费用，C5表示为人工费用，TOC表示***日运行费用。DN表示总淡水产量，ρw表示纯水密度，ρb表示浓盐水密度。
模拟计算模块将计算式(1)～(38)构成的稳态模拟问题运用拟牛顿法进行求解，确定当前***盐水加热器结垢系数fBH、热回收段和热排放段闪蒸室结垢系数fj、造水比GOR和日运行成本TOC，其具体步骤如下所示：
步骤B1：将上式(1)～(38)中有待求未知量的方程整理为非线性方程组，其具体表现形式如式(39)所示：
Figure BDA0002760582980000151
向量的表示形式为：
G(x)＝0 (40)；
这里，x＝(x1,x2,…,xp)T，G＝(g1,g2,…,gq)T，gi(i＝1,2,…,q):Rn→R。
其中，x1,x2,…,xp表示非线性方程组的p个未知数，x表示方程组的未知数构成的p×1维的矩阵。g1,g2,…,gq表示构成非线性方程组的q个等式，G表示非线性方程组的等式构成的q×p维的矩阵。
步骤B2：令初始迭代次数k＝0，设置初始值x0∈Rn，初始拟牛顿矩阵Bk为p×p维单位矩阵，函数值计算精度为ε，最小步长精度为δ。
步骤B3：计算迭代次数k下搜索方向向量dk，以得到下一个点的值：
xk+1＝xk+dk (41)；
dk＝-BkΔG(xk) (42)；
步骤B4：计算G(xk)和G(xk+1)，若||G(xk+1)||＜ε或者||xk+1-xk||＜δ，则终止迭代，得出非线性方程组的解及盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR、日运行成本TOC；否则，转至步骤B5。
步骤B5：通过修正Bk得到Bk+1：
Figure BDA0002760582980000152
其中，sk＝xk+1-xk，yk＝ΔG(xk+1)-ΔG(xk)。xk表示第k次迭代所求值，xk+1表示第k+1次迭代所求值，sk表示两次所求值之差。ΔG(xk+1)表示G(xk+1)的梯度向量，ΔG(xk)表示G(xk)的梯度向量，yk表达第k+1次迭代和k次迭代的梯度向量差。
步骤B6：令k:＝k+1，转至步骤B3进行计算。
参数决策模块根据模拟计算模块所得出的所有时间点的盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR和日运行费用TOC，采用最小二乘法进行拟合并进行分析，对***结垢程度进行合理的判断，具体步骤如下所示：
步骤C1：继承模拟计算模块从开始进行数据采集时的所有数据，并且针对于盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR和日运行费用TOC，对盐水加热器结垢系数fBH设定上限H1，闪蒸室结垢系数fj设定上限H2、造水比GOR设定下限X1、日运行费用TOC设定上限X2。
步骤C2：将上述数据分成四部分分别进行拟合，首先对盐水加热器结垢系数fBH采用最小二乘法进行曲线拟合。假设给定采集的数据点(xxi,yyi)(i＝0,1,…,m)，xxi表示数据采集时间点Ti，yyi表示每个时间点采集的盐水加热器结垢系数fBH(i)。拟合多项式为
Figure BDA0002760582980000161
为了保证模块计算速度令多项式阶数nn＝3，使得
Figure BDA0002760582980000162
其中I为a0,a1,a2,a3的四元函数，因此上述问题即为求I＝I(a0,a1,a2,a3)的极值问题。
步骤C3：由多元函数求极值的必要条件，得
Figure BDA0002760582980000163
即
Figure BDA0002760582980000164
步骤C4：将式(46)改写成正规方程组形式
Figure BDA0002760582980000171
可以证明正规方程组(47)的系数矩阵是一个对称正定矩阵，故存在唯一解。从式(47)中解得akk(kk＝0,1,…,3)，从而可得多项式
Figure BDA0002760582980000172
步骤C5：重复步骤C2至C4分别得出盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR和日运行费用TOC关于***运行时间的拟合曲线，为P1(t)、P2(t)、P3(t)和P4(t)，其中t代表时间。
步骤C6：根据参数所设上下限，和所求得的四条拟合曲线，判断当前***是否需要清洗，并根据拟合方程推出上述4个指标下的***待清洗时间点Tim1、Tim2、Tim3、Tim4。
步骤C7：将Tim1、Tim2、Tim3、Tim4进行比较得出最近的时间节点记为Tim，该时间点为***下次待清洗时间点。
步骤C8：对盐水加热器的结垢系数fBH、闪蒸室结垢系数fj的大小进行划分，分别对应不同等级的结垢程度，依据当前数据采集点的数值大小可得出***结垢程度。
以上内容是结合具体的实例给出的实施方式对本发明的进一步说明，不能认定本发明的具体实施只限于这些说明。对于本发明所属技术领域的技术人员来说，在不脱离发明构思的前提下还可以做出一定程度的简单推演或者替换，都应当视为属于本发明的保护范围。

Claims (2)

1.一种多级闪蒸海水淡化***结垢清洗方法，其特征在于：
包括如下步骤：首先由工程师通过人机交互模块设定多级闪蒸海水淡化***的结构参数，给定运行周期；此后多级闪蒸海水淡化***结垢清洗装置以给定的周期采集所需数据，并通过A/D和D/A转换模块发给中央处理模块，中央处理模块调用模拟计算模块；模拟计算模块通过计算中央处理模块中的模型问题，得出当前***盐水加热器的结垢系数fBH和闪蒸室的结垢系数fj、造水比GOR、日运行成本TOC；参数决策模块根据模拟计算模块所得参数，进行比较分析，得出当前***的结垢程度和下次***清洗时间并将相关结论发送至***显示模块；多级闪蒸海水淡化***结垢清洗装置在下一个运行周期，重复数据采集、模拟计算、比较分析过程，不断地实时估测***结垢程度；
具体包括如下步骤：
步骤A1：操作员或者工程师通过人机交互模块给定多级闪蒸海水淡化***的结构参数、物性系数和数据采集周期Tc；
步骤A2：利用数据采集模块采集当前时间点的进入热排放段的冷却盐水流量WF、循环海水流量WRe、加热蒸汽流量Wsteam、闪蒸盐水流量WB、闪蒸淡水流量WD、进料海水温度Tsea、加热蒸汽温度Tsteam、闪蒸盐水温度TB、闪蒸淡水温度TD、进料盐水浓度CF、闪蒸盐水浓度CB；记录当前时间，令变量T1等于当前时间，然后通过A/D和D/A转换模块将上述参数发送给中央处理模块；
步骤A3：中央处理模块根据内部存储的多级闪蒸海水淡化装置稳态模拟和日运行成本模型，调用模拟计算模块，从而计算出在当前T1时刻盐水加热器的结垢系数fBH、闪蒸室的结垢系数fj、造水比GOR和日运行成本TOC，然后将上述参数发送到参数决策模块；
所述的中央处理模块中的多级闪蒸海水淡化装置稳态模拟和日运行成本的模型由式(1)～(38)组成；
多级闪蒸海水淡化装置的稳态模型由闪蒸室方程、盐水加热器方程、混合分离方程、物性参数方程组成；对于第j级闪蒸室，闪蒸室模块模型由下列式(1)～(8)组成：
WBj-1+WDj-1＝WBj+WDj (1)；
WBj-1CBj-1＝WBjCBj (2)；
WBj-1hBj-1＝WBjhBj+VBjhVj (3)；
WBj-1-WBj＝VBj (4)；
Figure FDA0004051947610000021
Figure FDA0004051947610000022
其中，WBj表示第j级的闪蒸盐水质量流量，WBj-1表示第j-1级的闪蒸盐水质量流量，WDj表示第j级产出的淡水质量流量，WDj-1表示第j-1级的闪蒸盐水质量流量；CBj表示第j级的闪蒸盐水浓度，CBj-1表示第j-1级的闪蒸盐水浓度；hBj表示第j级闪蒸盐水的比焓，hBj-1表示第j-1级闪蒸盐水的比焓，hVj表示第j级闪蒸蒸汽的比焓，VBj表示第j级闪蒸室中盐水蒸发量；WF表示进入热排放段的冷却盐水的质量流量，CPRj表示离开第j级闪蒸室的冷却盐水热容，CPDj表示离开第j级闪蒸室的淡水比热容，CPDj-1表示离开第j-1级闪蒸室的淡水比热容，CPBj表示离开第j级闪蒸室的闪蒸盐水比热容，CPBj-1表示离开第j-1级闪蒸室的闪蒸盐水比热容，TFj表示离开第j级闪蒸室的冷却盐水温度，TFj+1表示进入第j级闪蒸室的冷却盐水温度，TBj表示离开第j级闪蒸室的盐水温度，TBj-1表示离开第j-1级闪蒸室的盐水温度，TDj表示离开第j级闪蒸室的淡水温度，TDj-1表示离开第j-1级闪蒸室的淡水温度，T*表示理想状态下闪蒸参考温度；Aj表示第j级闪蒸室的传热面积，
Figure FDA0004051947610000023
表示第j级闪蒸室的传热系数，其中W指的是第j级闪蒸室的进料流量，对于热排放段来说W＝WF，表示进料海水流量，而对于热回收段来说W＝WR，表示热回收段的进料流量，
Figure FDA0004051947610000024
表示闪各级蒸室外直径，
Figure FDA0004051947610000025
表示各级闪蒸室冷凝管内直径，Di表示内直径，其在物性方程部分有具体公式展开；
TBj＝TDj+ΔBPEj+ΔNETDj+ΔTLj (7)；
TVj＝TDj+ΔTLj (8)；
其中，ΔBPEj表示第j级盐水沸点升高，ΔNETDj表示非平衡余量，ΔTLj表示经过除雾器和冷凝器的温度损失；TVj表示第j级闪蒸室的闪蒸蒸汽温度；
盐水加热器模块模型由式(9)～(12)组成：
WB0＝WR (9)；
CB0＝CR (10)；
WRCPRH(TB0-TF1)＝Wsteamλs (11)；
其中，
Figure FDA0004051947610000026
Figure FDA0004051947610000027
其中，WB0表示离开盐水加热器的闪蒸盐水质量流量，WR表示进入热回收段的冷却盐水质量流量；CB0表示进入盐水加热器的闪蒸盐水浓度，CR表示进入热回收段的冷却盐水浓度；CPRH表示进入盐水加热器的冷却盐水热容，TB0表示离开盐水加热器的闪蒸盐水温度，Wsteam表示为加热蒸汽质量流量，λs表示蒸汽潜热，Tsteam表示为盐水加热器蒸汽温度；AH表示盐水加热器的传热面积，
Figure FDA0004051947610000031
表示盐水加热器的传热系数，其中W指的是盐水加热器的进料流量WR，
Figure FDA0004051947610000032
表示盐水加热器外直径，
Figure FDA0004051947610000033
表示盐水加热器冷凝管内直径，其在物性方程部分有具体公式展开；
混合分离模块模型由式(13)～(20)组成：
WBD＝WBN-WRe (13)；
Wm＝WF-Wr (14)；
S＝Wr-Cw (15)；
WR＝WRe+Wm (16)；
WR·CR＝WRe·CRe+Wm·Cm (17)；
WR·hR＝WRe·hRe+Wm·hm (18)；
WF＝S+WS (19)；
WF·hWF＝S·hS+WS·hWS (20)；
其中WBD表示废弃海水质量流量，WBN表示离开最后一级闪蒸室的盐水质量流量，WRe表示循环盐水质量流量；Wm表示补充盐水质量流量，Wr表示热排放海水质量流量；S表示热排放段回流流量，Cw表示热排放段排出盐水的质量流量；CRe表示循环盐水浓度，Cm表示补充盐水浓度；hR表示进入热回收段的盐水比焓，hRe表示循环盐水比焓，hm表示补充盐水比焓；WS表示进料海水质量流量；
Figure FDA0004051947610000034
表示热排放段入口盐水比焓，hS表示热排放段回流盐水比焓，
Figure FDA0004051947610000035
表示进料海水比焓；
物性参数方程模型由式(21)～(31)组成：
Figure FDA0004051947610000036
CPB＝[1-CB(0.011311-1.146×10-5TB)]×CPD (22)；
hV＝596.912+0.46694TS-0.000460256TS 2 (23)；
hB＝CPB·TB (24)；
hD＝CPD·TD (25)；
Figure FDA0004051947610000037
其中，CC＝(19.819CB)/(1-CB)，CC表示浓度转换，CB表示闪蒸盐水浓度；
Figure FDA0004051947610000041
其中，ωj＝WF/wj，ωj表示第j级进入热排放段的冷却盐水单位长度上的质量流量，wj表示第j级闪蒸室冷凝管的宽度；
TL＝exp(1.885-0.02063TD)/1.8 (28)；
U＝4.8857/(y+z+4.8857f) (29)；
y＝[0.0013(v×Di)0.2]/[(0.2018+0.0031×T)v] (30)；
Figure FDA0004051947610000042
GOR＝WDN/Wsteam (32)；
其中CPD表示闪蒸室的淡水比热容；CPB表示闪蒸室的盐水比热容，CB表示闪蒸盐水浓度，TB表示闪蒸盐水温度；hv表示蒸汽焓值，TS为蒸汽温度；hB表示盐水焓值；hD表示淡水焓值，TD表示淡水温度；BPE表示盐水沸点升高，其中，T表示处于闪蒸室或者盐水加热器的盐水温度；NETD表示非平衡温差，Hj表示闪蒸盐水液位，wj表示第j级闪蒸室冷凝管的宽度，ΔTB表示两级间的盐水温差；TL表示经过除雾器和冷凝管的温度损失；U表示盐水加热器或者各级闪蒸室的传热系数；y为有关管内流速、盐水温度和冷凝管内直径的中间表达式，z为有关闪蒸室淡水温度的中间表达式，y和z仅为简化计算，f为盐水加热器或者闪蒸室的结垢系数；v表示管内流速，Di表示冷凝管内直径，T表示在换热器出口处的盐水温度；WDN表示闪蒸室总淡水的质量流量，GOR表示造水比是反应***性能的一个关键指标；
日运行成本模块模型由式(33)～(38)组成：
C1＝22×Wsteam×[(Tsteam-40)/85]×0.00415 (33)；
C2＝22×(DN/ρw)×0.109 (34)；
C3＝22×(DN/ρw)×0.082 (35)；
C4＝22×(DN/ρb)×0.024 (36)；
C5＝22×(DN/ρw)×0.1 (37)；
TOC＝C1+C2+C3+C4+C5 (38)；
其中，C1表示加热蒸汽费用，C2表示装置消耗的电能费用，C3表示维护与空闲费用，C4表示为预处理费用，C5表示为人工费用，TOC表示***日运行费用；DN表示总淡水产量，ρw表示纯水密度，ρb表示浓盐水密度；
模拟计算模块将计算式(1)～(38)构成的稳态模拟问题运用拟牛顿法进行求解，确定当前***盐水加热器结垢系数fBH、热回收段和热排放段闪蒸室结垢系数fj、造水比GOR和日运行成本TOC，其具体步骤如下所示：
步骤B1：将上式(1)～(38)中有待求未知量的方程整理为非线性方程组，其具体表现形式如式(39)所示：
Figure FDA0004051947610000051
向量的表示形式为：
G(x)＝0 (40)；
这里，x＝(x1,x2,…,xp)T，G＝(g1,g2,…,gq)T，gi(i＝1,2,…,q):Rn→R；
其中，x1,x2,…,xp表示非线性方程组的p个未知数，x表示方程组的未知数构成的p×1维的矩阵；g1,g2,…,gq表示构成非线性方程组的q个等式，G表示非线性方程组的等式构成的q×p维的矩阵；
步骤B2：令初始迭代次数k＝0，设置初始值x0∈Rn，初始拟牛顿矩阵Bk为p×p维单位矩阵，函数值计算精度为ε，最小步长精度为δ；
步骤B3：计算迭代次数k下搜索方向向量dk，以得到下一个点的值：
xk+1＝xk+dk (41)；
dk＝-BkΔG(xk) (42)；
步骤B4：计算G(xk)和G(xk+1)，若||G(xk+1)||＜ε或者||xk+1-xk||＜δ，则终止迭代，得出非线性方程组的解及盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR、日运行成本TOC；否则，转至步骤B5；
步骤B5：通过修正Bk得到Bk+1：
Figure FDA0004051947610000052
其中，sk＝xk+1-xk，yk＝ΔG(xk+1)-ΔG(xk)；xk表示第k次迭代所求值，xk+1表示第k+1次迭代所求值，sk表示两次所求值之差；ΔG(xk+1)表示G(xk+1)的梯度向量，ΔG(xk)表示G(xk)的梯度向量，yk表达第k+1次迭代和k次迭代的梯度向量差；
步骤B6：令k:＝k+1，转至步骤B3进行计算；
步骤A4：参数决策模块将此周期得到的参数集L1(fBH,fj,GOR,TOC)与之前各周期得到的参数集Li(fBH,fj,GOR,TOC)，i＝1,2,…,m，m为当前周期数，进行比较分析，估算得出当前***的结垢程度以及***下次清洗时间，并且模拟计算模块和参数决策模块将其计算数据发送至中央处理模块；
步骤A5：中央处理模块将数据进行整理发送至显示模块，显示模块显示出当前时间T1、盐水加热器的结垢系数fBH、闪蒸室的结垢系数fj、造水比GOR、日运行成本TOC、***结垢程度、***待清洗时间Tim；
步骤A6：记录当前时间为T2，如果T2-T1<Tc，则继续等待；否则转步骤A2，重新进行数据采集。
2.如权利要求1所述的一种多级闪蒸海水淡化***结垢清洗方法，其特征在于：参数决策模块根据模拟计算模块所得出的所有时间点的盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR和日运行费用TOC，采用最小二乘法进行拟合并进行分析，对***结垢程度进行合理的判断，具体步骤如下所示：
步骤C1：继承模拟计算模块从开始进行数据采集时的所有数据，并且针对于盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR和日运行费用TOC，对盐水加热器结垢系数fBH设定上限H1，闪蒸室结垢系数fj设定上限H2、造水比GOR设定下限X1、日运行费用TOC设定上限X2；
步骤C2：将上述数据分成四部分分别进行拟合，首先对盐水加热器结垢系数fBH采用最小二乘法进行曲线拟合；假设给定采集的数据点(xxi,yyi)，i＝0,1,…,m，xxi表示数据采集时间点Ti，yyi表示每个时间点采集的盐水加热器结垢系数fBH(i)；拟合多项式为
Figure FDA0004051947610000061
为了保证模块计算速度令多项式阶数nn＝3，使得：
Figure FDA0004051947610000062
其中I为a0,a1,a2,a3的四元函数，因此即为求I＝I(a0,a1,a2,a3)的极值问题；
步骤C3：由多元函数求极值的必要条件，得：
Figure FDA0004051947610000063
即
Figure FDA0004051947610000064
步骤C4：将式(46)改写成正规方程组形式：
Figure FDA0004051947610000071
能够证明正规方程组(47)的系数矩阵是一个对称正定矩阵，故存在唯一解；从式(47)中解得akk(kk＝0,1,…,3)，从而可得多项式：
Figure FDA0004051947610000072
步骤C5：重复步骤C2至C4分别得出盐水加热器结垢系数fBH、闪蒸室结垢系数fj、造水比GOR和日运行费用TOC关于***运行时间的拟合曲线，为P1(t)、P2(t)、P3(t)和P4(t)，其中t代表时间；
步骤C6：根据参数所设上下限，和所求得的四条拟合曲线，判断当前***是否需要清洗，并根据拟合方程推出上述4个参数下的***待清洗时间点Tim1、Tim2、Tim3、Tim4；
步骤C7：将Tim1、Tim2、Tim3、Tim4进行比较得出最近的时间节点记为Tim，该时间点为***下次待清洗时间点；
步骤C8：对盐水加热器的结垢系数fBH、闪蒸室结垢系数fj的大小进行划分，分别对应不同等级的结垢程度，依据当前数据采集点的数值大小可得出***结垢程度。
CN202011216530.1A 2020-11-04 2020-11-04 一种多级闪蒸海水淡化***结垢清洗装置及方法 Active CN112678904B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011216530.1A CN112678904B (zh) 2020-11-04 2020-11-04 一种多级闪蒸海水淡化***结垢清洗装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011216530.1A CN112678904B (zh) 2020-11-04 2020-11-04 一种多级闪蒸海水淡化***结垢清洗装置及方法

Publications (2)

Publication Number Publication Date
CN112678904A CN112678904A (zh) 2021-04-20
CN112678904B true CN112678904B (zh) 2023-03-17

Family

ID=75445842

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011216530.1A Active CN112678904B (zh) 2020-11-04 2020-11-04 一种多级闪蒸海水淡化***结垢清洗装置及方法

Country Status (1)

Country Link
CN (1) CN112678904B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113158466A (zh) * 2021-04-23 2021-07-23 浙江捷创智能技术有限公司 一种磷化工稀磷酸供酸***管堵检修预测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1866140A (zh) * 2005-02-10 2006-11-22 通用电气公司 优化联合循环/联合生产设备的方法和装置
CN102706387A (zh) * 2012-05-31 2012-10-03 苏州市金翔钛设备有限公司 冷凝器污垢系数检测装置
CN205426201U (zh) * 2016-02-15 2016-08-03 北京洁禹通环保科技有限公司 一种冷凝器污垢积累程度检测装置
CN106872514A (zh) * 2017-02-21 2017-06-20 沈阳艾柏瑞环境科技有限公司 稳态换热过程传热系数及污垢热阻值在线监测***和方法
CN110143635A (zh) * 2019-04-03 2019-08-20 杭州电子科技大学 一种多级闪蒸海水淡化***的优化运行方法及装置
CN111059957A (zh) * 2019-12-28 2020-04-24 上海九瀚机电设备有限公司 一种实时监测板式换热器污垢系数的***及方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10508043B2 (en) * 2013-12-20 2019-12-17 Massachusetts Institute Of Technology Thermal desalination for increased distillate production

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1866140A (zh) * 2005-02-10 2006-11-22 通用电气公司 优化联合循环/联合生产设备的方法和装置
CN102706387A (zh) * 2012-05-31 2012-10-03 苏州市金翔钛设备有限公司 冷凝器污垢系数检测装置
CN205426201U (zh) * 2016-02-15 2016-08-03 北京洁禹通环保科技有限公司 一种冷凝器污垢积累程度检测装置
CN106872514A (zh) * 2017-02-21 2017-06-20 沈阳艾柏瑞环境科技有限公司 稳态换热过程传热系数及污垢热阻值在线监测***和方法
CN110143635A (zh) * 2019-04-03 2019-08-20 杭州电子科技大学 一种多级闪蒸海水淡化***的优化运行方法及装置
CN111059957A (zh) * 2019-12-28 2020-04-24 上海九瀚机电设备有限公司 一种实时监测板式换热器污垢系数的***及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"多级闪蒸海水淡化***的动态模拟与优化";高寒寒等;《计算机与应用化学》;20190628;第36卷(第3期);第212-222页 *

Also Published As

Publication number Publication date
CN112678904A (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
Alatiqi et al. Process control in water desalination industry: an overview
Ettouney Design of single-effect mechanical vapor compression
CN109211439B (zh) 一种汽轮机低压缸排汽焓值在线监测***及方法
CN112678904B (zh) 一种多级闪蒸海水淡化***结垢清洗装置及方法
Bayrak et al. Energy and exergy analyses of sugar production stages
Shakib et al. A new approach for process optimization of a METVC desalination system
Abdul‐Wahab et al. Development of a steady‐state mathematical model for multistage flash (MSF) desalination plant
Shen et al. Study of steam parameters on the performance of a TVC-MED desalination plant
Li et al. On-line fouling monitoring model of condenser in coal-fired power plants
CN106096310A (zh) 一种抽汽供热机组供热工况热力性能计算方法
CN111781904A (zh) 基于bp神经网络的植物提取的生产过程控制***
Jradi et al. Application of an artificial neural networks method for the prediction of the tube-side fouling resistance in a shell-and-tube heat exchanger
CN101825869A (zh) 基于数据驱动的过热器模型参数辨识方法
CN110143635B (zh) 一种多级闪蒸海水淡化***的优化运行方法及装置
El-Nashar Predicting part load performance of small MED evaporators-a simple simulation program and its experimental verification
Palenzuela et al. Modeling of the heat transfer of a solar multi-effect distillation plant at the Plataforma Solar de Almería
CN113685797B (zh) 一种余热锅炉省煤器变工况热力计算方法
Hassanean et al. Simulation of multi-stage flash with brine circulating desalination plant
CN106894853B (zh) 凝汽式汽轮机冷端节能诊断方法
Sadri et al. Multi objective optimization of the MED-TVC system with exergetic and heat transfer analysis
CN109189029B (zh) 火力发电厂低温省煤器节能量在线监测***及方法
Marmouch et al. Effect of a cooling tower on a solar desalination system
CN114925529A (zh) 凝汽器数字孪生模型的状态校正方法、***、设备及介质
Alatiqi et al. Measurements of dynamic behavior of a multistage flash water desalination system
Ibrahimi et al. Embodiment design of a mechanical vapor compression desalination unit

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