一种基于群体智能优化算法的高聚物浆液参数识别方法与流程

本发明属于化学注浆领域，尤其涉及一种基于群体智能优化算法的高聚物浆液参数识别方法。

背景技术：

注浆作为工程基础设施加固修复的重要手段，在交通、水利、市政等领域得到了广泛应用。发展快捷高效的注浆加固修复新材料、新技术和新理论始终是人们关注的热点，尤其是近年来，随着现代高分子工业的进步，具有自膨胀特性的双组分聚氨酯类高聚物注浆材料及其高压注射技术在国际上发展十分迅速，应用日趋广泛，成为化学灌浆领域较为活跃的发展方向之一，已广泛应用于矿井、隧道等地下工程以及水害防治、地基加固、道路维修等方面，成为岩土工程领域具有明显特色的发展方向。

研究高聚物浆液在被注介质中的扩散机理对于注浆设计和施工具有重要指导意义。借助仿真手段模拟高聚物在被注介质中的膨胀扩散过程是研究高聚物浆液扩散机理的重要途径。高聚物浆液的内在化学反应机制是其膨胀扩散的原动力，实现对高聚物浆液化学反应过程的准确求解是借助仿真手段研究浆液扩散机理的前提，其关键在于准确确定高聚物浆液化学反应动力学参数(指前因子和活化能)。

确定化学反应动力学参数的传统方法是通过试验测定，但该方法只能针对单一化学反应实施，由于高聚物浆液的膨胀扩散过程包含胶凝反应和发泡反应两个并行的化学反应，受两个同步发生的化学反应的耦合作用以及外加剂等因素的影响，与单一反应相比，浆液化学反应动力学参数会发生变化，进而改变反应速率。cn110470559a公开了一种固体高温高压化学反应动力学参数获取方法与装置，其获取方法为：电炉升至指定温度并通入高压反应气体，将固体样品送入电炉炉膛内发生反应并实时测量样品质量，通过对样品质量曲线的微分获取固体高温高压化学反应动力学参数。该方法也是通过试验并采集试验数据以及结合样品质量曲线推算出化学反应动力学参数，并且此方法也并不适合存在多个同步反应条件下的高聚物浆液化学反应动力学参数的获取。如何实现对存在多个同步反应条件下的高聚物浆液化学反应动力学参数的准确辨识，目前尚缺乏有效的手段，成为制约高聚物注浆理论研究和技术发展的瓶颈。

技术实现要素：

针对现有技术的不足，本发明的目的是提供一种基于群体智能优化算法的高聚物浆液参数识别方法，以解决现有技术中对存在多个同步反应条件下的高聚物浆液化学反应动力学参数的识别缓慢且不准确的问题。

为了解决上述技术问题，本发明采用了如下的技术方案：

本发明提供了一种基于群体智能优化算法的高聚物浆液参数识别方法，包括以下步骤：

s10、预处理：确定所述智能优化算法的种群规模n、个体维数d、最大迭代次数t、收敛精度ε；

s20、初始化：按照所述智能优化算法规则生成种群，为种群中每个个体赋初始值，每个个体代表一组高聚物浆液化学反应动力学参数，计算每个个体对应的初始目标函数值；赋予每个个体的个体历史最优值和相应的个体历史最优目标函数值；赋予整个群体的群体历史最优值和相应的群体历史最优目标函数值；为当前迭代次数赋初值；

s30、迭代求解：按照所述智能优化算法规则更新种群每个个体值，计算每个个体相应的优化目标函数值，根据个体当前所述优化目标函数值与所述个体历史最优目标函数值及所述群体历史最优目标函数值的关系，确定是否更新所述个体历史最优值和相应的所述个体历史最优目标函数值以及所述群体历史最优值和相应的所述群体历史最优目标函数值；

s40、判断是否执行了预定次数的运算或所述群体历史最优目标函数值是否达到了精度要求，若是，则执行步骤s50；若否，则执行步骤s30；

s50、输出参数识别结果：输出所述群体历史最优目标函数值对应的化学反应动力学参数，即为识别得到的化学反应动力学参数。

进一步的，所述智能优化算法为粒子群算法、遗传算法、蚁群算法、鱼塘算法。

进一步的，所述步骤s20具体包括：

s21、以化学反应动力学参数取值可行域作为求解空间，在求解空间内随机产生种群的每个个体的初始值形式为xi＝(anco,enco,aoh,eoh)，i＝1:n，所述anco代表每个个体在x轴上离原点的距离，所述enco代表每个个体在y轴上离原点的距离，所述aoh代表每个个体在z轴上离原点的距离，所述eoh代表每个个体在时间轴上离原点的距离；

s22、用表达式(1)计算每个个体对应的所述初始目标函数值fi；

所述表达式(1)为：

所述表达式(1)中，为实测温度值，为计算温度值，ρk为实测密度值，ρ′k为计算密度值，k为记录点序列号，m为记录点总数，α、β为权重系数；

s23、赋予每个个体的所述个体历史最优值和相应的所述个体历史最优目标函数值：将每个个体的所述初始值xi和相应的所述初始目标函数值fi，分别作为各所述个体历史最优值pi和相应的所述个体历史最优目标函数值fpi的初值；

s24、赋予整个群体的所述群体历史最优值和相应的所述群体历史最优目标函数值：在各所述个体历史最优目标函数值fpi的初值中选择其中的最小值作为所述群体历史最优目标函数值fg的初值，相应的所述个体历史最优值pi初值中的最小值作为所述群体历史最优值pg初值；

s25、设置所述当前迭代次数t＝1。

进一步的，所述步骤s22中，通过联立求解高聚物浆液胶凝反应速率方程(2)、发泡反应速率方程(3)和能量守恒方程(4)得到，ρ′k通过表达式(5)计算得到，和ρk分别通过开展高聚物浆液化学反应试验过程中对不同时刻浆液的温度和体积的测量得到；

所述高聚物浆液胶凝反应速率方程(2)为：

所述发泡反应速率方程(3)为：

所述能量守恒方程(4)为：

所述表达式(5)为：

所述高聚物浆液胶凝反应速率方程(2)、所述发泡反应速率方程(3)、所述能量守恒方程(4)、所述表达式(5)中，aoh为凝胶反应的指前因子，eoh为凝胶反应的活化能，aw为发泡反应的指前因子，ew为发泡反应的活化能，rg为理想气体常数，t为浆液绝对温度，xoh为胶凝反应的转化率，xw为发泡反应的转化率，cnco,0、coh,0、cw,0分别为异氰酸酯、多元醇以及水组分在反应混合物中的当前摩尔浓度，ρp、ρbl、ρw分别为异氰酸酯与多元醇混合液、液态物理发泡剂、水的密度，c为比热，λ为蒸发热，(△h)oh为凝胶反应的反应热，(△h)w为发泡反应的反应热，ρ为浆液密度，rw,0是水组分的初始质量，rbl,0是物理发泡剂溶解度的初始值，rw、rbg、rbl分别为各组分的当前质量，t是反应时间。

进一步的，所述步骤s30具体包括：

s31：按照所述智能优化算法规则更新种群每个个体当前值xi；

s32、利用所述表达式(1)计算每个个体相应的优化目标函数值fi；

s33、比较个体当前目标函数值fi^t和所述个体历史最优目标函数值fpi，如果fi^t<fpi，则更新fpi为fpi＝fi^t，所述个体历史最优值pi用个体当前值替换；

如果fpi更新，进一步比较其与所述群体历史最优目标函数值fg的大小，若fpi<fg，则更新fg＝fpi，所述群体历史最优值pg也相应替换为pi。

进一步的，所述步骤s40具体包括：

当t＜t且fg>ε时，执行步骤s30，同时t＝t+1；

如果t＞t或fg≤ε，则计执行步骤s50；

其中所述ε为设定的收敛精度。

本发明提供的一种基于群体智能优化算法的高聚物浆液参数识别方法与现有技术相比，至少具有如下有益效果：

通过利用群体智能优化算法从而实现了对存在多个同步化学反应条件下的高聚物浆液化学反应动力学参数的快速准确识别。

附图说明

为了更清楚地说明本申请的方案，下面将对实施例描述中所需要使用的附图作一个简单的介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明提供的一种基于群体智能优化算法的高聚物浆液参数识别方法的总流程示意图；

图2为本发明提供的一种基于群体智能优化算法的高聚物浆液参数识别方法的实施例一实测温度与基于识别参数得到的计算温度的对比图；

图3为本发明提供的一种基于群体智能优化算法的高聚物浆液参数识别方法的实施例一实测密度与基于识别参数得到的计算密度的对比图。

具体实施方式

为了便于理解本申请，下面将参照相关附图对本申请进行更全面的描述。附图中给出了本申请的较佳实施例。但是，本申请可以以许多不同的形式来实现，并不限于本文所描述的实施例。相反地，提供这些实施例的目的是使对本申请的公开内容的理解更加透彻全面。除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中在本申请的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本申请。

实施例一

本实施例中，采用基于群体智能优化算法的自膨胀高聚物浆液化学反应动力学参数识别方法，结合自膨胀高聚物浆液化学反应过程中的温度、密度测试数据，对浆液化学反应动力学参数实施辨识。首先设定自膨胀高聚物浆液的组分参数及相关物理参数：浆液组分浓度coh,0为4079.6mol/m³，cnco,0为4557.7mol/m³，cw,0为294.7mol/m³，密度ρp为1100kg/m³，ρw为1000kg/m³，反应热(-△h)w为8.6×10⁴j·mol^-1，(-△h)oh为7.075×10⁴j·mol^-1，rbl,0为1228kg/m³，为4.4×10⁴，cp为1800j/(kg·k)，为836.6j/(kg·k)，cbl为1159j/(kg·k)，cbg为1000j/(kg·k)，λ为206.8kj/kg。

本实施例中的群体智能优化算法具体采用粒子群算法。

s10、预处理：设置种群规模n＝60，粒子维数d＝4，最大迭代次数t＝50，收敛精度ε＝0.00001；

s20、初始化：按照所述智能优化算法规则生成种群，为种群中每个个体赋初始值，每个个体代表一组高聚物浆液化学反应动力学参数，计算每个个体对应的初始目标函数值；赋予每个个体的个体历史最优值和相应的个体历史最优目标函数值；赋予整个群体的群体历史最优值和相应的群体历史最优目标函数值；为当前迭代次数赋初值；

s21、以化学反应动力学参数取值可行域作为求解空间，在求解空间内随机产生种群的每个个体的初始值形式为xi＝(anco,enco,aoh,eoh)，i＝1:n，所述anco代表每个个体在x轴上离原点的距离，所述enco代表每个个体在y轴上离原点的距离，所述aoh代表每个个体在z轴上离原点的距离，所述eoh代表每个个体在时间轴上离原点的距离；

s22、用表达式(1)计算每个个体对应的所述初始目标函数值fi；

所述表达式(1)为：

所述表达式(1)中，为实测温度值，为计算温度值，ρk为实测密度值，ρ′k为计算密度值，k为记录点序列号，m为记录点总数，α、β为权重系数；

所述步骤s22中，通过联立求解高聚物浆液胶凝反应速率方程(2)、发泡反应速率方程(3)和能量守恒方程(4)得到，ρ′k通过表达式(5)计算得到，和ρk分别通过开展高聚物浆液化学反应试验过程中对不同时刻浆液的温度和体积的量测得到。

所述表达式(2)为：

所述表达式(3)为：

所述表达式(4)为：

所述表达式(5)为：

所述高聚物浆液胶凝反应速率方程(2)、所述发泡反应速率方程(3)、所述能量守恒方程(4)、所述表达式(5)中，aoh为凝胶反应的指前因子，eoh为凝胶反应的活化能，aw为发泡反应的指前因子，ew为发泡反应的活化能，rg为理想气体常数，t为绝对温度，xoh为胶凝反应的转化率，xw为发泡反应的转化率，cnco,0、coh,0、cw,0分别为异氰酸酯、多元醇以及水组分在反应混合物中的当前摩尔浓度，ρp、ρbl、ρw分别为混合浆液、液态物理发泡剂、水的密度，c为比热，λ为蒸发热，(△h)oh为凝胶反应的反应热，(△h)w为发泡反应的反应热，ρ为密度，rw,0是水组分的初始质量，rbl,0是溶解度的初始值，rw、rbg、rbl分别为各组分的当前质量，t是反应时间。

s23、赋予每个个体的所述个体历史最优值和相应的所述个体历史最优目标函数值：将每个个体的所述初始值xi和相应的所述初始目标函数值fi，分别作为各所述个体历史最优值pi和相应的所述个体历史最优目标函数值fpi的初值；

s24、赋予整个群体的所述群体历史最优值和相应的所述群体历史最优目标函数值：在各所述个体历史最优目标函数值fpi的初值中选择其中的最小值作为所述群体历史最优目标函数值fg的初值，相应的所述个体历史最优值pi初值中的最小值作为所述群体历史最优值pg初值；

s25、设置所述当前迭代次数t＝1。

s30、迭代求解：按照所述智能优化算法规则更新种群每个个体值，计算每个个体相应的优化目标函数值，根据个体当前所述优化目标函数值与所述个体历史最优目标函数值及所述群体历史最优目标函数值的关系，确定是否更新所述个体历史最优值和相应的所述个体历史最优目标函数值以及所述群体历史最优值和相应的所述群体历史最优目标函数值；

s30具体步骤如下：

s31、按表达式(6)更新每个粒子的位置；

所述表达式(6)为：

s32、利用表达式(1)计算每个粒子的目标函数值fi；

s33、对每个粒子进行比较和替换操作：

比较粒子当前目标函数值fi^t和个体历史最优目标函数值fpi，如果fi^t<fpi，则更新fpi为fpi＝fi^t，个体历史最优值pi用当前位置替换；如果fpi更新，进一步比较其与群体历史最优目标函数值fg的大小，若fpi<fg，则更新fg＝fpi，对应的群体历史最优值pg也相应替换为个体历史最优值pi；

s34、按表达式(7)更新每个粒子的速度。

所述表达式(7)为：

所述表达式(7)中，t表示当前迭代次数，t-1表示当前迭代次数的上一次迭代次数；r1表示第一随机因子、r2表示第二随机因子，取值范围在[0,1]；w为惯性权重，用来控制速度，平衡算法的整体和局部搜索能力；c1为算法的第一学习因子、c2为算法的第二学习因子，通常取两个相同的非负常数，取值范围在[0,4]，本实施例中，第一学习因子c1＝0.5，第二学习因子c2＝0.5；pi为粒子个体历史最优值，pg为群体历史最优值。

所述表达式(7)中，惯性权重w可用表达式(8)计算；

所述表达式(8)为：

w＝(t-t)×1.0/t(8)

s40、判断是否执行了预定次数的运算或所述群体历史最优目标函数值是否达到了精度要求，若是，则执行步骤s50；若否，则执行步骤s30；

s50、输出参数识别结果：输出所述群体历史最优目标函数值对应的化学反应动力学参数，即为识别得到的化学反应动力学参数。

表1为本发明提供的一种基于群体智能优化算法的高聚物浆液参数识别方法的实施例一最终识别得到的化学反应动力学参数。

表1

如表1所示可知，本发明提供的一种基于群体智能优化算法的高聚物浆液参数识别方法的实施例一最终得到的化学反应动力学参数在参数取值范围内，从图2、图3可以看出，本发明提供的一种基于群体智能优化算法的高聚物浆液参数识别方法中实测温度与计算温度几乎相同、实测密度与计算密度几乎相同，说明本发明提供的一种基于群体智能优化算法的高聚物浆液参数识别方法所得到的的化学反应动力学参数是较准确的。

显然，以上所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例，附图中给出了本申请的较佳实施例，但并不限制本申请的专利范围。本申请可以以许多不同的形式来实现，相反地，提供这些实施例的目的是使对本申请的公开内容的理解更加透彻全面。尽管参照前述实施例对本申请进行了详细的说明，对于本领域的技术人员来而言，其依然可以对前述各具体实施方式所记载的技术方案进行修改，或者对其中部分技术特征进行等效替换。凡是利用本申请说明书及附图内容所做的等效结构，直接或间接运用在其他相关的技术领域，均同理在本申请专利保护范围之内。