买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

龙图腾网获悉河北工业大学申请的专利一种考虑不确定性的脱碳海水进料反渗透脱硼海水淡化系统优化方法获国家发明授权专利权，本发明授权专利权由国家知识产权局授予，授权公告号为：CN115331744B 。

龙图腾网通过国家知识产权局官网在2025-10-10发布的发明授权授权公告中获悉：该发明授权的专利申请号/专利号为：202210979398.2，技术领域涉及：G16C20/10；该发明授权一种考虑不确定性的脱碳海水进料反渗透脱硼海水淡化系统优化方法是由杜亚威;魏敬欣;李亚新;邓会宁;张少峰设计研发完成，并于2022-08-16向国家知识产权局提交的专利申请。

本一种考虑不确定性的脱碳海水进料反渗透脱硼海水淡化系统优化方法在说明书摘要公布了：本发明公开了一种考虑不确定性的脱碳海水进料反渗透脱硼海水淡化系统优化方法。该方法采用严格机理的方式建立反渗透脱硼传递机理的数值模型，将海水中的碳酸脱除后结垢倾向大幅度降低，各级反渗透均提高进水pH值，硼截留率大幅度提高，建立了反渗透系统的超结构模型，考虑了沿压力容器轴向盐水压力、浓度、流量和pH值的变化，物料分配器和混合器考虑了物料、盐度、硼浓度和pH值平衡，加入操作条件约束保证系统安全运行。采用二级随机优化模型考虑运行过程中的不确定性因素，使得产水满足饮用水中硼含量的标准，根据不同情景调节流程和操作条件进一步降低能耗，更具有鲁棒性，为其工业化应用提供理论基础和技术参考，具有非常好的应用前景。

本发明授权一种考虑不确定性的脱碳海水进料反渗透脱硼海水淡化系统优化方法在权利要求书中公布了：1.一种考虑不确定性的脱碳海水进料反渗透脱硼海水淡化系统优化方法，其特征在于该方法包括以下步骤： 步骤1：建立卷式膜元件反渗透脱硼海水淡化过程模型； 根据反渗透过程机理，采用微分方程描述压力容器内盐度、硼浓度、压力、流量沿轴向变化，将有限差分方法微分方程离散化，得到如下方程： Cboric,l＝α0,lCTB,l8 Kboric，l＝2Kl13 其中A、B、BTB分别为纯水、盐和硼透过系数，P表示压力，C表示盐度，π表示渗透压，ρp和Vw分别表示淡水的密度和流速，Jw和Js分别为纯水通量和盐透过通量，Vw为渗透流速，T为温度，K为传质系数，de为进料流道的当量直径，Sl为膜元件一个微分单元的面积，Sl＝Sm·nmL，Sm为单支膜元件的面积，nm为压力容器膜元件的个数，L为总的微分单元节点数，Re为雷诺准数，Re＝ρVdeμ，其中μ为动力粘度，Sc为施密特准数，Sc＝μρDs，Ds为盐的扩散系数，Q为流量，V为进料流速，V＝Q3600Sfcsεsp，Sfcs为进料流道横截面积，εsp为进料流道隔网孔隙率，σ为反射系数，α0和α1分别为硼酸和硼酸盐离子的分率，pKa为硼酸一级电离常数，Kλ为摩擦系数，FFd为污染系数，e为膜的活化能T≤298K时，e＝25,000Jmol-1，T298K时，e＝22,000Jmol-1，R是气体常数，Bin为每年盐透过增加率，Nmlp为反渗透膜的平均寿命，NRO表示总的反渗透级，kλ表示膜元件内摩擦系数，下标ch为膜元件的进料或产水流道，b为浓盐水，boric为硼酸分子，borate为硼酸盐离子，f为进料海水，p为产水，mw为膜表面，ref为T0在298K时没有污染时膜的参数，l为微分单元节点，j表示第j个反渗透级； 有限差分法的边界条件： z＝0，V＝Vin，Q＝Qin，Cboric＝Cboric,in，Cborate＝Cborate,in，C＝Cin，P＝Pin； 其中in表示压力容器进口； 产水分流设计将压力容器前端的产水直接输送至最终产水，后端的产水进入下一级反渗透再脱盐，压力容器两端加入流量调节阀调节两端出水比例，其模型表述为： Qch,b,l+1＝Qch,b,l-3600Vw，lSl15 Qch,b，l+1Cch,b,l+1-Qch，b，lCch，b，l＝-3600Vw,lSlCch,p，l16 Qch,b,l+1CTB,ch,b,l+1-Qch,b,lCTB,ch,b,l＝-3600Vw,lSlCTB,ch,p,l17 Qf,n＝Qb，n+Qp,n,lc+Qp,n,hc24 Qf，nCf,n＝Qb，nCb,n+Qp,n,lcCp,n,lc+Qp,n,hcCp,n,hc25 Qf,nCTB,f,n＝Qb，nCTB,b,n+Qp,n,lcCTB,p,n,lc+Qp,n,hcCTB，p，n，hc26 Yl-Yl+1≥027 其中二元变量Y表示压力容器内微分单元的产水流向，下标lc表示压力容器前端产水，hc表示压力容器后端产水，n表示第n个压力容器；公式27表示在压力容器内前端或后端每个微分单元产水有一致的流向； 盐水渗透压π、动力粘度μ和盐的扩散系数Ds由下列拟合公式计算： π＝4.54047103CMsρ0.98728 μ＝1.4757×10-3+2.4817×10-6C+9.3287×10-9C2exp-0.02008T29 Ds＝6.725×10-6exp0.1546×10-3C-2513T+273.1530 其中Ms为溶质的摩尔质量； 步骤2.建立反渗透超结构模型； 反渗透系统包括海水取水和前处理、产水后处理、反渗透膜组、泵、功交换器pressureexchanger，PX、物流混合器及分离器，由于海水中存在碳酸根离子，在反渗透操作过程中易形成沉淀，将海水加强酸后调节pH值到4.0，通过风机将海水中游离的CO2脱除，脱碳海水可有效降低结垢倾向，加入强碱将pH值提高至碱性有效提高硼截留率，反渗透超结构中包含NPS个增压级和NRO个反渗透级构成，总共有NPS+2个物流节点，2是指最终离开反渗透系统的盐水和淡水，NPS个物流节点中的每一个节点表示有一股物流经过高压泵增压后或不经高压泵不增压直接进入1个反渗透单元，每一个反渗透级由多个并联压力容器构成，每个压力容器由2～8个膜元件串联且在相同条件下运行，离开反渗透级的每一股盐水和淡水都可以进入NPS+2个物流节点，每一个物流表示为流量、盐度、硼浓度和压力的函数，每一个进料MIN经过物流分配器后分为MOUT个物流，一些物流经过物流混合器后汇成一股物流，物流分配器和混合器表示为： Cin,out＝Cinout＝1，...MOUT32 CTB,in,out＝CTB,inout＝1，...MOUT33 Pin,out＝Pinout＝1，...MOUT34 pHin,out＝pHinout＝1，...MOUT35 0＝Pin-PoutQin,outin＝1，...MIN40 公式31-35表示物流分配器，公式36-39表示物流混合器，公式40表示等压混合约束，允许反渗透级的产水与系统最终产水、反渗透级的高压浓盐水泄压后与系统进料混合，总的硼浓度CTB为硼酸分子Cboric和硼酸盐离子Cborate浓度的和，Qin、Cin、CTB,in和Pin分别表示物流分配器进口的流量、盐度、硼浓度和压力，Qout、Cout，CTB,out和Pout分别表示物流混合器出口的流量、盐度、硼浓度和压力，Qin,out表示物流分配器出口流量，Cin,out、CTB,in,out和Pin,out表示物流分配器出口盐度、硼浓度和压力，QacidandQbase分别为酸和碱的添加量，CacidandCbase分别为酸和碱的浓度，分别为酸和碱液中氢离子浓度，下标in，out分别表示进口和出口； 高压泵与功交换器的物料平衡方程为： Qps，1＝Qhpp+Qpxlin41 Qps,1Cps，1＝QhppChpp+QpxlinCpxlin42 Qps，1CTB，ps，1＝QhppCTB,hpp+QpxlinCTB,pxlin43 QRO，1＝Qhpp+Qpxhout44 QRO，1CRO，1＝QhppChpp+QpxhoutCpxhout45 QRO，1CTB，RO，1＝QhppCTB,hpp+QpxhoutCTB,pxhout46 Qpxhout＝Qpxlin47 Qpxhin＝Qpxlout48 Lpxpxhin100＝Qpxhin-Qpxhout49 px[％]＝0.3924+0.01238Ppxhin50 pxhout＝MixCpxhin-Cpxlin+Cpxlin51 TB,pxhout＝MixCTB,pxhin-CTB,pxlin+CTB，pxlin52 Mix＝6.0057-0.3559OF+0.0084OF53 OF[％]＝100×Qpxhin,-QpxhoutQpxhin54 pxloutpxlout＝Qpxlinpxlin+Qpxhinpxhin-Qpxhoutpxhout55 TB,pxloutpxlout＝QpxlinTB，pxlin+QpxhinTB,pxhin-QpxhoutTB,pxhout56 其中LPX为泄漏率，Mix为体积混合率，OF-10％≤OF≤15％为润滑流量，下标hpp、pxhin、pxlin、pxhout和pxhin分别表示高压泵、进入功交换器的低压进料海水和高压浓盐水、离开功交换器的加压海水和泄压浓盐水，下标ps表示增压级，RO表示反渗透级； 离开第i个增压级的物流直接进入第j个反渗透级，同一级反渗透压力容器内的采用相同种类的膜元件，其特性如纯水透过系数、溶质透过系数、膜面积和进料隔网厚度保持不变，由如下公式确定第j级反渗透内压力容器采用的膜元件型号k： -Pk，max≤U1-yj,kj＝1，2，...，NRO,k＝1，2，...，Kt58 引入二元变量yj，k，当它取1时表示在第j级反渗透选取第k种类型的膜元件，否则取0；公式58限定了膜元件允许的最大进水压力，U是一个足够大的数，Kt是反渗透膜元件的种类集合； 引入公式60和61确定反渗透级的浓盐水和产水进入下一个反渗透级，如果二元变量βb，i，j等于1，则第j级的浓盐水进入下一个反渗透级，当βp，i,j等于1，则第j级的产水进入下一个反渗透级，公式62可避免浓盐水和产水同时进入下一个反渗透级； 随着水透过膜，其浓盐水和产水的pH值也会发生变化，采用下列公式对其进行描述； 其中表示第j级反渗透的回收率，rl,j表示第j级反渗透压力容器内每个微分单元的回收率，可以确定第j个反渗透级的进料海水是上一个反渗透级或反渗透段； 海水的pH值为8.2，为有效去除海水中的CO2，使用Qacid,0+Qacid，1的强酸将海水调节到pH值为4.0，然后通过吹扫塔将海水中的CO2脱除，之后加入Qbase,0+Qbase,1的强碱将海水pH值提高以增加硼脱除率， Qf+Qacid,010-7+Qacid,1·Cacid＝Qf+Qacid,0+Qacid,1]10-469 Qf+Qacid,0+Qacid,110-4＝Qbase，0Cbase70 其中公式71计算提高反渗透级进料pH值所需的强碱量； 整个反渗透网络满足如下物料平衡关系，以及产品水需求约束： Qf＝Qb+Qp72 QfCf＝QbCb+QpCp73 QfCTB，f＝QbCTB，b+QpCTB,p74 Qp≥Qp，lo81 Cp≤Cp,up82 CTB，p≤CTB,p,up83 式中Qb、Cb和CTB，b分别表示离开反渗透网络的盐水流量、盐度和硼浓度，Qp、Cp和CTB,p分别表示产品水的流量、盐度和硼浓度，下标lo和up分别表示最低需求值和最大允许值； 步骤3.系统流量和操作条件约束； 为保障反渗透系统安全运行，在模型中设定以下约束：浓差极化因子是膜表面盐浓度Cch.mw.l与主体溶液盐度Cch,b.l的比值，脱碳海水一级反渗透浓差极化因子极限值为1.22，传统海水进料一级反渗透浓差极化因子极限值为1.2，二级反渗透由于其进水含盐量已经显著降低，浓差极化因子最大为1.4；单支压力容器最大压力损失为0.35MPa，第一级和第二级平均产水通量最大值分别为20Lm2·h和40Lm2·h，第一级和第二级首支膜元件最大产水通量分别为35Lm2·h和48Lm2·h，第一级和第二级压力容器内最小浓盐水流量分别为3.6m3h和2.4m3h，浓盐水浓度小于90kgm3，反渗透进料pH值最高为9.5，反渗透级进料pH值最高为11.0； 步骤4.建立反渗透系统优化设计模型 反渗透系统的优化设计问题表达为一个混合整数非线性规划，以公式79为目标函数、满足过程热力学、单元操作、设计要求约束： 目标函数Ew考虑了系统能耗、总取水量、膜组件规模和酸碱试剂添加量，ΔPSWIP、ΔPhpp、ΔPbp和ΔPbppx分别表示海水取水泵、高压泵、级间泵和增压泵的压差，Qf、Qhpp、Qbp和Qbppx分别表示进料海水、高压泵、级间泵和增压泵的流量，η为设备的效率，fc为反渗透装置的载荷系数，nm，j和npv，j表示第j级反渗透压力容器内膜元件个数和第j级反渗透压力容器个数，下标moter为泵的电机； 步骤5.对形成的反渗透系统优化问题进行求解； 优化目标为Ew公式84， 约束方程为：反渗透过程模型公式1-30， 反渗透超结构模型公式31-83， 系统流量和操作条件约束； 采用二级随机模型优化不确定性条件下反渗透系统的流程和运行条件，随机优化模型中的不确定性参数表达为概率分布Ω＝{ω1，ω2}，则两级随机模型表示为： 约束： 约束： 前两个方程表示为第一级模型，后两个方程为第二级模型，和为决策变量，包括反渗透模型中的各级流量Q、压力P、pH、物流分配器出口变量Qin,out、表示第j级反渗透压力容器内膜元件个数nm,j、第j级反渗透压力容器个数npv,j、二元变量Yl和二元变量yj,k；cT是目标函数中的参数，是第二级模型最优解的数学期望，公式85表示反渗透系统优化问题中的约束方程，其中是约束方程的系数，和是约束方程的右边项，TωandWω分别为转置和补偿矩阵；在第一级模型中，nm，j和npv,j保持不变，其他决策变量在第二级模型中根据不确定性因素进行调节； 采用数学规划软件求解以上混合整数非线性规划问题，通过给变量赋不同的初值，从多个初始点出发进行迭代，可获得系统优化的流程和操作条件。

如需购买、转让、实施、许可或投资类似专利技术，可联系本专利的申请人或专利权人河北工业大学，其通讯地址为：300130 天津市北辰区西平道5340号；或者联系龙图腾网官方客服，联系龙图腾网可拨打电话0551-65771310或微信搜索“龙图腾网”。

以上内容由龙图腾AI智能生成。