DE-GWO algorithm to optimize SVM inversion mechanical parameters of soft soil
-
摘要: 岩土物理力学参数是岩土工程设计及施工的主要依据。为得到更符合实际的岩土力学参数,研究提出了一种更加快捷准确的参数确定方法,即基于支持向量机(SVM)的位移反分析法。基于发展于自然仿生学的群智能算法-灰狼算法(GWO),结合差分进化算法(DE)的交叉变异操作对灰狼算法进行进化,得到一种差分进化灰狼混合优化算法(DEGWO),用其对支持向量机参数进行优化后反演出岩土力学参数。结合昆明某地铁隧道项目,采用Midas/GTS有限元软件构建模拟隧道施工的三维仿真模型。以某段隧道顶地表实测沉降数据为依据,用差分进化灰狼混合算法优化后的支持向量机模型反演岩体力学参数,并据反演参数及地勘报告建议参数计算得出的、Peck法计算的与实际地表监测的沉降值进行对比分析,结果表明：据反分析得出的参数计算的地表沉降值与实测值基本吻合,另外两种方法的计算值与实测值相差较大。通过以上分析表明本方法合理可行,反演出参数合理,用其进行计算的结果准确更符合实际,本方法能为岩土力学参数的获取、校核提供新的途径。Abstract: The geotechnical mechanical parameters are the main basis for geotechnical engineering design and construction.In order to obtain the geotechnical parameter information that is more in line with the actual situation, a displacement back analysis method based on the support vector machines, faster and more accurate parameter determination method, is proposed.Based on the grey wolf algorithm, a swarm intelligence algorithm developed in nature bionics, and combined with the crossover mutation operation of the differential evolutionary algorithm, the grey wolf algorithm is evolved, and the hybrid grey wolf optimization algorithm is obtained, which is used to reverse the geotechnical mechanical parameters after optimization of the support vector machine parameters.Based on the subway tunnel project in Kunming, the Midas/GTS finite element is used to establish a three-dimensional simulation model to simulate the tunnel construction process.Based on the measured surface deformation of the top section of a tunnel, the support vector machine model optimized by the hybrid differential evolution grey wolf algorithm is used to invert the mechanical parameters of rock mass.The settlements calculated by the back analysis method and the traditional Peck formula method are compared with the meas surface ones.The results show that the calculated surface settlements based on the back analysis are basically consistent with the measured ones, and the calculated values by the other two methods differ greatly from the measured ones.The above analysis shows that the proposed method is reasonable and feasible, the reverse performance parameters are reasonable, and the calculated results are more accurate and in line with the reality.The proposed method can provide a new way for the acquisition and verification of the parameters of rock and soil mechanics.
-
0. 引言
微生物诱导碳酸钙沉积（microbially induced calcite precipitation，简称MICP）是本世纪初提出的一项岩土工程实用技术，主要利用以碳酸钙为产物的微生物矿化过程，在岩土材料中形成类似水泥的胶凝材料，进而对岩土材料起到加固作用[1-2]。由于基于尿素水解的MICP技术具有反应过程容易控制、能够在短时间内产生大量CO32-、生成碳酸钙结晶的反应过程和胶凝强度容易控制以及所生成的结晶物在一定条件下可替代传统胶凝材料等特点[3-7]，国内外研究人员通过开展天然岩石节理现场试验和模拟岩石节理室内试验对MICP技术封堵岩石节理的适用性进行了广泛研究[8-16]，试验结果表明：采用MICP方法对岩石节理进行多轮次灌注后，由于MICP反应产生的碳酸钙沉积在岩石节理中导致其隙宽显著下降，MICP封堵后岩石节理渗透系数最大降幅达4个数量级，封堵效果良好。
尽管如此，Lambert等[17]针对加拿大南安大略地区地下16 m处水平灰岩裂隙开展原位激发MICP封堵试验，结果发现注入营养液52 d后灰岩裂隙导水系数最大降幅达62%，封堵效果较为理想，但是230 d后复测发现灰岩裂隙导水系数出现一定程度恢复，这表明岩体裂隙MICP封堵效果在复杂地下环境中可能出现退化。Minto等[18]针对固碳储层岩体CO2泄露问题，分析了Berea砂岩岩芯MICP封堵效果及所生成的碳酸钙弱酸性侵蚀情况，结果表明尽管Berea砂岩岩芯中优势渗流通道被MICP封堵，阻止了酸性地下水在岩芯中的渗流，抑制了碳酸钙侵蚀，但在岩芯边缘与酸性地下水接触的碳酸钙产生了明显侵蚀现象。上述研究成果表明耐久性是岩体MICP封堵需要关注的关键问题，同时由于MICP产生的碳酸钙封堵层对酸性环境较为敏感，酸性侵蚀是影响MICP封堵耐久性的最不利因素，因此抗酸性侵蚀性能是影响MICP封堵耐久性的关键因素。尽管近年来已有学者关注碳酸钙酸性溶蚀行为对MICP胶结砂土长期性能的影响[19-20]，目前关于MICP封堵处理岩石节理抗酸性侵蚀方面的研究成果还鲜见报导。
鉴于此，本文以单一岩石节理为研究对象，开展MICP封堵岩石节理试件在不同PH值侵蚀溶液中侵蚀不同时间的酸性侵蚀试验，研究酸性侵蚀后MICP封堵岩石节理碳酸钙封堵层溶蚀情况及其在酸性侵蚀前后的渗透性演化规律，在此基础上分析岩石节理MICP封堵抗酸性侵蚀性能，为MICP技术在岩体封堵工程中的应用提供科学依据。
1. MICP封堵岩石节理试件
1.1 3D打印类岩石节理试件制备
由于本次MICP封堵岩石节理酸性侵蚀试验需要大量具有相同表面形貌特征的岩石节理试件，在此以人工劈裂花岗岩节理试件（平面尺寸长×宽=100 mm×50 mm）为原型，先采用三维激光扫描仪获取其表面三维形貌数据，构建人工劈裂花岗岩节理数字模型，再采用透明光敏树脂（型号Crysta-8QEF1）为打印材料，利用3D打印技术打印透明树脂类岩石节理试件，制备流程如图 1所示。固化后的透明光敏树脂表面无孔隙，几乎不吸水，物理特性与原型花岗岩近似。本次试验共打印了9组透明树脂类岩石节理试件（分别编号为1-1~1-9），用于后续MICP封堵试验和酸性侵蚀试验。
![]()
图 1 3D打印透明树脂类岩石节理试件制备过程Figure 1. Preparation of transparent resin rock joint specimens完整花岗岩基本物理力学参数、人工劈裂花岗岩节理JRC值以及透明树脂基本物理参数汇总如表 1所示。其中，完整花岗岩基本物理力学参数通过3组标准圆柱体花岗岩试件室内单轴压缩试验结果得到；人工劈裂花岗岩节理JRC值计算先利用表面三维形貌数据获得沿节理试件长度方向的节理剖面线二维坐标数据，再按照Tse等[21]提出的JRC计算公式计算：
JRC=32.2+32.47lgZ2。 (1) 式中：Z2为岩石节理剖面线坡度均方根一阶导数。
表 1 花岗岩和固化透明光敏树脂基本物理力学参数Table 1. Physical and mechanical parameters of granite and resin透明光敏树脂材料 花岗岩和花岗岩节理 邵氏硬度 弯曲模量/
GPa拉伸模量/
GPa拉伸强度/
MPa热变形温度/
℃密度/
(g·cm-3)单轴抗压强度/MPa 弹性模量/
GPa节理JRC值 76~82 2.57~2.86 2.49~2.6 44~60 40~51 2.586 155.0 9.0 11.56 1.2 3D打印类岩石节理试件MICP封堵处理
当3D打印类岩石节理试件制备好后，为了避免MICP封堵试验中试件入口因为碳酸钙沉积而造成堵塞，首先在节理表面4个角点处放置不同厚度小垫片形成具有不同隙宽值的岩石节理试件开展MICP封堵预试验，最终确定隙宽为2 mm时试件入口处基本不会发生堵塞。因此，正式MICP封堵试验中在节理表面4个角点处放置2 mm厚小垫片以形成隙宽2 mm的节理空腔，沿节理长度方向用纳米胶密封试件两侧，并在试件宽度一侧用纳米胶固定一个“Y”型接头作为MICP反应灌注入口，另一侧为灌注出口自然出流并贴上滤纸防止MICP反应生成的碳酸钙流出，密封好的3D打印类岩石节理试件如图 2所示。为了测试3D打印类岩石节理试件初始渗透系数K0，在进行MICP封堵处理前先采用如图 3所示装置对其进行沿节理长度方向的单向渗流试验，在此基础上采用立方定理估算其渗透系数，最终测得3D打印类岩石节理试件初始渗透系数K0为9.814×10-2 m/s。
3D打印类岩石节理试件MICP封堵试验采用单向非循环灌注方式，试验装置如图 4所示，包括蠕动泵、试件固定支架、连接软管、菌液和反应液储罐以及废液收集桶等。其中，试件固定支架具有角度调节功能，角度调节范围为0°~90°。本次试验菌种采用购自上海保藏生物技术中心的巴氏芽孢杆菌（编号为SHBCC D10472），菌种经活化培养后测试菌液OD600值为0.824。反应液采用1.0 moL/L等摩尔尿素-氯化钙混合溶液，其中尿素为细菌生长提供氮源和能量，氯化钙为MICP反应提供钙源。
![]()
图 4 岩石节理MICP封堵试验装置Figure 4. Test devices for sealing rock joint specimens by MICP method岩石节理MICP封堵试验过程中，首先将固定支架调整为水平，将密封好的岩石节理试件置于支架上并接入蠕动泵，通过蠕动泵向岩石节理空腔注满菌液并静置2 h，保证细菌充分附着在岩石节理表面。然后开始菌液和反应液单向非循环灌注，灌注前把固定支架调整为向下倾斜一定角度，先以0.5 mL/min速率连续泵送菌液和反应液60 min，静置2 h；再以1 mL/min速率连续泵送菌液和反应液30 min，再次静置2 h；最后再以0.5 mL/min速率泵送菌液和反应液直至灌注入口处无法泵送或出现漏液情况，停止灌注。为了避免岩石节理空腔内部碳酸钙填充不密实，在岩石节理单向非循环灌注结束后将固定支架调整为向下倾斜90°，通过蠕动泵以0.3 mL/min速率泵入约20 mL菌液和20 mL胶结液并静置3 h，重复5次以保证岩石节理空腔能够被碳酸钙完全填充密实。
按照上述步骤得到的MICP封堵处理3D打印类岩石节理试件（以下简称“MICP封堵岩石节理试件”）如图 5所示，将其置于40℃烘箱中连续烘干72 h后取出冷却称重，再继续烘干6 h后重新取出冷却称重，若两次称重差值不超过3 mg，则认为试件中碳酸钙层已经完全干燥固化，可以用于后续酸性侵蚀试验，最终得到酸性侵蚀试验前9组MICP封堵岩石节理试件中碳酸钙质量（记为M1）数值汇总如表 2所示。此外，MICP封堵岩石节理试件在进行酸性侵蚀试验之前先接入图 3所示自制渗流试验装置，测试试件在MICP封堵后的渗透系数K1，具体数值如表 2所示。
![]()
图 5 MICP封堵3D打印类岩石节理试件Figure 5. MICP-treated 3D printed artificial rock joint specimens表2 MICP封堵前后以及酸性侵蚀后3D打印类岩石节理试件碳酸钙质量和渗透系数
Table 2. Calcium carbonates and permeabilities of artificial rock joint before and after MICP treatment and acid erosion tests, respectively试件编号 试件中碳酸钙质量/g 试件渗透系数/(10-3m·s-1) M1 M2 ΔM K0 K1 K2 1-1 7.406 7.203 0.203 98.14 2.64 2.98 1-2 8.140 7.957 0.183 3.43 3.72 1-3 8.822 8.709 0.113 2.93 3.00 1-4 7.341 7.021 0.320 3.21 4.14 1-5 6.483 6.227 0.256 3.69 4.52 1-6 7.943 7.823 0.120 2.78 2.89 1-7 9.141 8.271 0.870 2.34 3.57 1-8 8.75 8.147 0.603 2.44 3.48 1-9 7.568 7.452 0.116 2.98 3.11 2. MICP封堵岩石节理酸性侵蚀试验
根据《工程地质手册》（第五版）[22]关于酸性和弱酸性地下水的分类，本次试验考虑地下水为酸性和弱酸性的情况，通过对浓度0.01 mol/L的盐酸溶液进行稀释并以0.01 mol/L的NaOH溶液作为调节液配制得到pH值为4和5的稀盐酸溶液，以此模拟酸性地下水，同时配制pH值为7的中性溶液作为对照组。此外，由于酸性地下水对碳酸钙封堵层的侵蚀程度与侵蚀时间密切相关，因此在本次MICP封堵岩石节理酸性侵蚀试验中综合考虑试验效率和侵蚀效果，设置7，14，28 d三种侵蚀时间。由此设置MICP封堵岩石节理试件酸性侵蚀试验的9种工况及其对应的试件编号如表 3所示。
表 3 MICP封堵岩石节理试件酸性侵蚀试验工况Table 3. Test conditions of MICP-treated artificial rock joint specimens during acid erosion experiments试件
编号1-1 1-2 1-3 1-4 1-5 1-6 1-7 1-8 1-9 pH值 4 5 7 4 5 7 4 5 7 侵蚀时间/d 7 7 7 14 14 14 28 28 28 试验中，将按3种预设pH值配制好的酸性侵蚀溶液分别倒入3个聚丙烯透明浸泡箱中，然后将图 5所示9组MICP封堵岩石节理试件（拆除入口纳米胶和“Y”型接头）放入浸泡箱并保证试件完全浸没在酸性溶液液面以下。浸泡过程中每隔8 h采用手持式pH酸度计测定酸性溶液的pH值，若溶液pH值增大则及时补充稀盐酸溶液以保证其pH值恒定。当浸泡时间达到7 d时取出3组试件（试件编号为1-1~1-3），按照1.2节所述烘干称重方法测定酸性侵蚀后MICP封堵岩石节理试件中碳酸钙质量（记为M2），在此基础上结合酸性侵蚀前试件中碳酸钙质量M1即可求得酸性侵蚀试验中碳酸钙溶蚀量ΔM。对于浸泡时间分别为14，28 d的MICP封堵岩石试件（对应试件编号分别为1-4~1-6和1-7~1-9）同样按照上述方法计算碳酸钙溶蚀量，汇总结果见表 2。
此外，为了分析酸性侵蚀后MICP封堵岩石节理渗透性变化规律，将经历不同酸性侵蚀时间的MICP封堵岩石节理试件接入图 3所示渗流试验装置，测试对应的渗透系数K2，结果见表 2。
3. 试验结果
3.1 碳酸钙溶蚀率
为了定量表征MICP封堵岩石节理试件碳酸钙封堵层在酸性侵蚀试验中的溶蚀情况，在此提出碳酸钙溶蚀率（记为δm）概念，其定义为MICP封堵岩石节理试件中碳酸钙溶蚀量ΔM与酸性侵蚀前该试件中碳酸钙质量M1之比，由此得到MICP封堵岩石节理试件在3种不同pH值溶液中侵蚀不同时间的碳酸钙溶蚀率如图 6所示。由图 6可知：①当酸性侵蚀溶液pH值一定时，MICP封堵岩石节理试件碳酸钙溶蚀率随侵蚀时间增长呈显著增加趋势，以酸性侵蚀溶液pH=4为例，当侵蚀时间由7 d增加到28 d时，试件碳酸钙溶蚀率从2.74%增长至9.52%，增幅达247%，但在中性溶液中试件碳酸钙溶蚀率非常小且随侵蚀时间增长并无明显变化；②当侵蚀时间相同时，MICP封堵岩石节理试件碳酸钙溶蚀率随侵蚀溶液酸性增强而增大，以试件酸性侵蚀28 d为例，当侵蚀溶液由弱酸性（pH=5）变为酸性（pH=4）时，试件碳酸钙溶蚀率从6.89%增长至9.52%，增幅为38.2%，而侵蚀溶液由中性（pH=7）变为弱酸性（pH=5）时，试件碳酸钙溶蚀率从1.53%增长至6.89%，增幅更是高达350%。
![]()
图 6 MICP封堵岩石节理试件碳酸钙溶蚀率变化情况Figure 6. Dissolution rates of MICP-treated rock joint specimens after acid erosion tests3.2 酸性侵蚀后MICP封堵岩石节理渗透系数变化
为了分析MICP封堵岩石节理试件在酸性侵蚀前后的渗透系数变化规律，根据表 2数据先计算MICP封堵前后岩石节理试件渗透系数变化情况，得到9组3D打印类岩石节理试件在MICP封堵后其渗透系数平均降幅为97%；再计算9组MICP封堵岩石节理试件在酸性侵蚀前后对应的渗透系数变化情况，得到不同pH值侵蚀溶液作用下MICP封堵岩石节理试件渗透系数变化量随侵蚀时间的变化曲线如图 7所示。
由图 7可知：①经MICP封堵后3D打印类岩石节理试件渗透系数平均降幅为97%，说明MICP方法对岩石节理具有显著的封堵效果；②MICP封堵岩石节理试件经酸性侵蚀后，其渗透系数比未侵蚀前均有增大趋势，增幅在2.39%~52.56%，且试件渗透系数增幅与侵蚀溶液pH值和侵蚀时间密切相关，即：当侵蚀溶液pH值一定时，酸性侵蚀后MICP封堵岩石节理试件渗透系数增幅随侵蚀时间增长而显著增加，而当MICP封堵岩石节理试件侵蚀时间相同时，对应的渗透系数增幅随侵蚀溶液酸性增强而增大。
4. 分析与讨论
4.1 碳酸钙溶蚀分析
MICP封堵岩石节理酸性侵蚀本质上是MICP生成的碳酸钙在酸性溶液侵蚀作用下发生溶蚀，进而影响MICP封堵岩石节理的渗透性。为了研究酸性侵蚀后MICP封堵岩石节理试件渗透系数变化与碳酸钙溶蚀量之间的定量关系，根据表 2所示酸性侵蚀后岩石节理试件碳酸钙溶蚀量和酸性侵蚀前后岩石节理试件渗透系数变化量，建立酸性侵蚀后MICP封堵岩石节理试件渗透系数增幅与碳酸钙溶蚀率δm之间的经验关系如式（2）所示，拟合曲线见如图 8。由图 8可知：酸性侵蚀后MICP封堵岩石节理试件渗透系数增幅随碳酸钙溶蚀量变化基本呈线性增长趋势。
η=5.635δm(R2=0.9804)， (2a) η=K2−K1K1×100%。 (2b) 式中：η为MICP封堵岩石节理试件渗透系数增幅；δm为碳酸钙溶蚀率，定义同3.1节；K1和K2含义同表 2。
![]()
图 8 酸性侵蚀后MICP封堵岩石节理试件渗透系数增幅与碳酸钙溶蚀率拟合曲线Figure 8. Fitting curve between permeability increment and CaCO3 dissolution rate of MICP-treated rock joint specimen after acid erosion需要说明的是，由于式（2）时针对表 3所示9种酸性侵蚀试验工况结果回归分析得到的，其不适用MICP封堵岩石节理试件中碳酸钙100%酸性溶蚀的情况。
同时，为了查明酸性侵蚀试验中MICP封堵岩石节理试件碳酸钙溶蚀部位，在利用图 3所示渗流试验装置测试酸性侵蚀后MICP封堵岩石节理试件渗透系数时采用红墨水进行渗流试验，由此获得酸性侵蚀后MICP封堵岩石节理试件在渗流试验过程中的渗流区域。以酸性侵蚀28 d的MICP封堵岩石节理试件为例，3组试件稳定渗流过程中最终形成的渗流区域如图 9所示，图中红色区域即为水流渗流区域。
![]()
图 9 酸性侵蚀28 d后MICP封堵岩石节理试件渗流区域Figure 9. Preferential flow domain in MICP-treated rock joint specimens after acid erosion for 28 days由图 9可知：①对于1-7试件，其经历pH=4酸性溶液侵蚀28 d，渗流试验中水流从试件渗流入口处开始渗流，经由试件上下两侧边缘区域形成向渗流出口的渗流流径，渗流区域面积较大且出口处渗流区域几乎连通，这说明在酸性溶液侵蚀作用下，1-7试件渗流入口附近区域、上下侧边缘附近区域以及渗流出口处的碳酸钙封堵层均有溶蚀，从而为水流渗流提供通道。②对于1-8试件，其经历pH=5弱酸性溶液侵蚀作用，渗流试验中水流从渗流入口处至试件长度二分之一范围内均有渗流，但随后主要沿试件上侧边缘附近区域形成向渗流出口的渗流流径，渗流出口处渗流区域也主要集中在试件上侧边缘附近区域，说明1-8试件中碳酸钙溶蚀位置与1-7试件有一定差异；同时与1-7试件相比，1-8试件中水流渗流区域总面积也有显著减小，这说明1-8试件中碳酸钙溶蚀量也有减少，这与图 6所示碳酸钙溶蚀试验结果是一致的。③对于1-9试件，其在pH=7中性溶液中侵蚀，渗流试验中水流从渗流入口处靠近上侧区域开始，主要沿试件渗流上侧边缘附近区域形成向渗流出口的渗流流径，但是在渗流出口附近区域水流突然从试件上侧附近区域朝试件右下角区域渗流，试件中渗流区域面积进一步缩小。
综上所述，本次酸性侵蚀试验中MICP封堵岩石节理试件碳酸钙侵蚀部位主要集中在试件左侧渗流入口附近以及试件上下侧附近等易与酸性侵蚀溶液接触区域，而试件中部区域溶蚀程度较低。
4.2 MICP封堵岩石节理抗酸性侵蚀性能分析
文献[17]关于MICP封堵Berea砂岩岩芯稀盐酸侵蚀试验结果表明，采用MICP方法封堵Berea砂岩岩芯后其渗透系数相比封堵前下降了95.53%，由于稀盐酸溶液溶解了砂岩岩芯中封堵孔隙的碳酸钙，导致稀盐酸溶液侵蚀后MICP封堵Berea砂岩岩芯渗透系数略有回升，其渗透系数下降幅度收窄至94.63%。但是由于文献[17]针对MICP封堵Berea砂岩岩芯的稀盐酸溶液侵蚀时间相对较短，酸性侵蚀前后MICP封堵Berea砂岩岩芯渗透系数变化不明显。
鉴于此，本文MICP封堵岩石节理试件酸性侵蚀试验设置了7，14，28 d三种侵蚀时间。根据图 7所示MICP封堵岩石节理试件酸性侵蚀前后的渗透系数变化规律可知：MICP封堵岩石节理试件经酸性侵蚀后，其渗透系数相比未侵蚀前均有增大趋势，增幅在2.39%~52.56%，且MICP封堵岩石节理试件渗透系数增幅与侵蚀时间密切相关，即：当侵蚀溶液pH值一定时，酸性侵蚀后MICP封堵岩石节理试件渗透系数增幅随侵蚀时间增长而显著增加，如侵蚀溶液为酸性（pH=4）时，当侵蚀时间由7 d变为28 d，对应的试件渗透系数增幅则由12.88%增长至52.56%。因此，为了分析MICP封堵岩石节理试件抗酸性侵蚀性能，将本次酸性侵蚀试验中在不同pH值酸性溶液中经历侵蚀时间最长（28 d）的1-7~1-9试件作为分析对象，将3组试件对应的初始渗透系数k0、MICP封堵后渗透系数k1以及酸性侵蚀28 d后渗透系数k2试验结果进行比较，得到岩石节理试件在经历MICP封堵和酸性侵蚀作用后的渗透系数变化情况如图 10所示。由图 10可知：1-7~1-9等3组岩石节理试件经MICP封堵后的渗透系数分别为2.34×10-3，2.44×10-3，2.98×10-3 m/s，相比MICP封堵前3组试件渗透系数降幅分别为97.62%，97.51%，96.96%，而后3组试件分别在3种不同pH值酸性溶液中侵蚀28 d后对应的渗透系数分别为3.57×10-3，3.48×10-3，3.11×10-3 m/s，相比MICP封堵前3组试件在酸性侵蚀28 d后的渗透系数降幅分别为96.36%，96.45%，96.83%，其原因在于酸性溶液溶解了MICP封堵岩石节理试件中的部分碳酸钙，导致其渗透系数略有回升，这与文献[17]研究成果是一致的。尽管如此，经历不同pH值酸性溶液侵蚀28 d后3组MICP封堵岩石节理试件渗透系数降幅依然高达95%以上，说明在历经最长28 d酸性侵蚀后MICP封堵岩石节理试件碳酸钙封堵层总体抗渗效果良好。
![]()
图 10 MICP封堵岩石节理试件经历酸性侵蚀28 d后的渗透系数变化情况Figure 10. Permeability change of MICP-treated rock joint specimens after acid erosion for 28 days在此基础上再结合图 9所示酸性侵蚀后MICP封堵岩石节理试件渗流流径可知：碳酸钙酸性溶蚀是引起MICP封堵岩石节理试件渗透系数回升的主要原因，而岩石节理中碳酸钙溶蚀首先发生在碳酸钙封堵层与酸性溶液接触位置，且酸性溶液会沿着岩石节理试件碳酸钙溶蚀区发生渗流，这又导致不断有新的碳酸钙被溶蚀，形成新的碳酸钙溶蚀区，为酸性溶液渗流提供新的通道，即：碳酸钙酸性溶蚀和酸性溶液渗流之间存在耦合作用，这正是影响MICP封堵岩石节理试件抗酸性侵蚀性能的关键。因此，如何防止碳酸钙封堵层发生初始酸性溶蚀是保证MICP封堵岩石节理抗酸性侵蚀耐久性需要重点解决的技术问题。工程实践中当采用MICP方法封堵处理岩石节理后，应当密切监测处理区域地下水pH值变化情况，一旦地下水呈现弱酸性就要采取措施提升其pH值，保证碳酸钙封堵层附近地下水呈中性或碱性。
5. 结论
通过开展MICP封堵岩石节理试件在不同pH值和侵蚀时间下的酸性侵蚀试验，研究MICP封堵岩石节理中碳酸钙封堵层溶蚀情况以及其在酸性侵蚀前后的渗透性演化规律，分析其抗酸性侵蚀性能，得到以下3点结论。
（1）采用MICP方法封堵岩石节理后，与封堵前相比其渗透系数显著下降，平均下降幅度达97%，MICP封堵效果明显。
（2）在酸性侵蚀作用下，MICP封堵岩石节理试件中碳酸钙封堵层有轻微溶蚀现象，碳酸钙溶蚀量随侵蚀溶液酸性增强和侵蚀时间增长而增大，且试件中碳酸钙侵蚀部位主要集中在试件四周边缘附近区域；
（3）MICP封堵岩石节理试件在经历酸性侵蚀后其渗透系数会略微回升，且酸性侵蚀后试件渗透系数增幅与碳酸钙溶蚀率呈线性增长关系，但经历不同pH值酸性溶液侵蚀28 d的MICP封堵岩石节理试件渗透系数降幅均大于95%，MICP封堵效果依然良好。
此外，本次MICP封堵岩石节理试件酸性侵蚀试验设置的最长侵蚀时间仅为28 d，从分析岩石节理碳酸钙封堵层酸性侵蚀长期耐久性角度来看是不够的，需要进一步开展更长侵蚀时间的室内侵蚀试验甚至是现场试验，同时还应关注MICP封堵岩石节理试件CaCO3封堵层酸性侵蚀后对试件剪切力学特性的弱化影响，这将是本文后续研究拟重点开展的研究内容。
-
表 1 模型土层参数取值表
Table 1 Values of parameters of model soil
土类 c/kPa φ/(°) Es/MPa γ/(kN·m-3) μ e 素填土 10 5.0 5.0 17.0 0.42 1.00 黏质粉土 16 7.0 8.5 18.9 0.30 0.75 泥炭质土 21 5.3 2.5 13.2 0.39 3.10 黏土 30 11.0 4.5 17.8 0.36 1.10 加固土 46 12.5 15.0 18.0 0.35 0.60 
下载: 导出CSV
表 2 力学参数反演范围表
Table 2 Inversion ranges of mechanical parameters
指标 Es/MPa μ e ρ/(g·cm-3) 泥炭质土 1.5~6 0.35~0.45 1.8~5.78 1.09~1.51
下载: 导出CSV
表 3 正态分布学习样本表
Table 3 Orthogonal experiment learning samples
样本号 Es/MPa μ e ρ/(g·cm-3) BP31/mm BP32/mm BP33/mm BP34/mm … BP42/mm BP43/mm BP44/mm 1 2.4815 0.4498 3.7968 1.3384 -26.4095 -40.9652 -24.6060 -7.1573 … -2.7355 -4.8089 -3.9057 2 4.7481 0.3831 4.2207 1.3524 -18.4219 -45.5279 -18.7600 -10.0022 … -2.4166 -4.4952 -4.2558 3 4.7432 0.4176 3.6040 1.2932 -17.7189 -32.7310 -26.6711 -10.6636 … -4.2959 -4.7174 -4.2234 ︙ ︙ ︙ ︙ ︙ ︙ ︙ ︙ ︙ … ︙ ︙ ︙ 48 4.0956 0.4062 4.0299 1.3321 -26.5025 -32.2280 -20.8879 -6.0910 … -2.5950 -4.7490 -4.2098 49 3.6584 0.4165 3.3979 1.2334 -18.5138 -34.1621 -19.6161 -6.0939 … -2.0894 -3.5371 -2.9590 50 2.5926 0.4074 3.0930 1.1614 -23.8329 -29.1038 -24.8184 -11.4369 … -2.0067 -5.3758 -4.9246
下载: 导出CSV
表 4 均匀试验测试样本表
Table 4 Uniform test samples
样本 Es/MPa μ e ρ/(g·cm-3) BP31/mm BP32/mm BP33/mm BP34/mm … BP42/mm BP43/mm BP44/mm 1 1.5 0.3944 3.1266 1.5103 -22.5433 -41.1159 -22.1167 -8.9535 … -4.8597 -3.3338 -3.3443 2 2.0 0.4499 2.6844 1.4636 -25.9139 -39.9132 -16.5514 -9.1734 … -2.7393 -3.2539 -3.0694 3 2.5 0.4166 1.8000 1.4169 -23.4582 -38.0973 -16.9451 -12.6672 … -4.3348 -3.1331 -4.0029 ︙ ︙ ︙ ︙ ︙ ︙ ︙ ︙ ︙ … ︙ 8 5.0 0.3611 4.0110 1.1834 -21.4858 -29.2345 -25.7267 -10.2678 … -4.1049 -5.0648 -3.4669 9 5.5 0.3833 5.7798 1.1367 -26.2777 -39.2444 -21.3229 -10.2249 … -3.6693 -4.2848 -4.0764 10 6.0 0.3500 4.8954 1.0900 -25.7881 -45.0308 -24.9058 -6.7625 … -4.9053 -5.4120 -3.2425
下载: 导出CSV
表 6 实测地表沉降数据表
Table 6 Data of surface settlement
监测点位 沉降值/mm 监测点位 沉降值/mm BP31 -22.43 BP38 -22.38 BP32 -38.17 BP39 -22.60 BP33 -22.20 BP40 -13.51 BP34 -8.97 BP41 -3.35 BP35 -3.95 BP42 -3.67 BP36 -2.27 BP43 -3.94 BP37 -2.84 BP44 -3.49
下载: 导出CSV
表 7 地表沉降值对比差异表
Table 7 Comparison of differences in settlements
地表监测点位 BP31 BP32 BP33 BP34 BP35 BP36 BP37 BP38 BP39 BP40 BP41 BP42 BP43 BP44 监测值/mm -22.40 -38.10 -22.2 -8.97 -3.95 -2.27 -2.84 -22.30 -22.60 -13.50 -3.35 -3.67 -3.94 -3.49 反演值/mm -21.10 -35.30 -21.6 -9.51 -4.18 -2.56 -2.69 -21.00 -21.40 -14.00 -3.18 -3.26 -3.59 -3.12 相对误差/% 5.75 7.34 2.3 -6.02 -5.82 -12.70 5.28 6.08 5.13 -3.77 5.07 11.17 8.88 10.60 建议值计算/mm -18.20 -29.40 -17.6 -8.01 -2.85 -1.56 -1.71 -18.40 -18.80 -11.10 -2.78 -2.92 -3.01 -2.73 相对误差/% 18.77 22.80 20.3 10.70 27.85 31.28 39.70 17.40 16.70 17.60 17.01 20.44 23.60 21.78 Peck计算值/mm -31.20 -49.70 -30.9 -15.30 -9.61 -5.32 -5.61 -31.10 -31.10 -20.40 -7.78 -7.96 -8.32 -7.82 相对误差/% -39.10 -30.30 -39.5 -70.90 -143.00-134.00 -97.50 -39.20 -38.00 -51.60-132.00 -116 -111.00 -124.00
下载: 导出CSV
-
[1] 王文杰, 叶世伟. 人工智能原理与应用[M]. 北京: 人民邮电出版社, 2004. WANG Wen-jie, YE Shi-wei. Principles and applications of artificial intelligence[M]. Beijing: Posts&Telecom Press, 2004. (in Chinese)
[2] CERNYCH M, SATAS A, BRAZAITIS M. Post-sauna recovery enhances brain neural network relaxation and improves cognitive economy in oddball tasks[J]. International Journal of Hyperthermia, 2018, 35(1): 375-382. doi: 10.1080/02656736.2018.1504992
[3] 李宁, 王李管, 贾明涛, 等. 改进遗传算法和支持向量机的岩体结构面聚类分析[J]. 岩土力学, 2014, 35(增刊2): 405-411. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2014S2058.htm LI Ning, WANG Li-guan, JIA Ming-tao, et al. Application of improved genetic algorithm and support vector machine to clustering analysis of rock mass structural plane[J]. Rock and Soil Mechanics, 2014, 35(S2): 405-411. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2014S2058.htm
[4] WANG Ding-wei, WANG Jun-wei, WANG Hong-feng. Intelligent OptimizationMethods[M]. Beijing: Higher Education Press, 2007.
[5] 漆祖芳, 姜清辉, 周创兵, 等. 基于v-SVR和MVPSO算法的边坡位移反分析方法及其应用[J]. 岩石力学与工程学报, 2013, 32(6): 1185-1196. doi: 10.3969/j.issn.1000-6915.2013.06.012 QI Zu-fang, JIANG Qinghui, ZHOU Chuang-bing, et al. A new slope displacement back analysis method based ON v-SVR and mvpso algorithm and its application[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(6): 1185-1196. (in Chinese) doi: 10.3969/j.issn.1000-6915.2013.06.012
[6] Khatibinia Mohsen, Feizbakhsh Abdosattar, Ehsan Mohseni and Malek Mohammad Ranjbar. Modelingmechanical strength of self-compacting mortar containing nanoparticles using wavelet-based support vector machine[J]. Computers and Concrete, 2016, 18(6): 1065-1082.
[7] 张宝珍, 张尧, 林凌雪. 基于改进差分进化算法的估计等值法[J]. 华南理工大学学报(自然科学版), 2014, 42(4): 7-12, 18. doi: 10.3969/j.issn.1000-565X.2014.04.002 ZHANG Bao-zhen, ZHANG Yao, LIN Ling-xue. Estimation equivalence method based on modified differential evolution algorithm[J]. Journal of South China University of Technology(Natural Science Edition), 2014, 42(4): 7-12, 18. (in Chinese) doi: 10.3969/j.issn.1000-565X.2014.04.002
[8] YANG X S, DEB S. Cuckoo search via lévy flights[C]//2009 World Congress on Nature&Biologically Inspired Computing(NaBIC). Coimbatore, 2009: 210-214
[9] 郭振洲, 刘然, 拱长青, 等. 基于灰狼算法的改进研究[J]. 计算机应用研究, 2017, 34(12): 3603-3606, 3610. doi: 10.3969/j.issn.1001-3695.2017.12.019 GUO Zhen-zhou, LIU Ran, GONG Chang-qing, et al. Studyonimprovementofgraywolfalgorithm[J]. Application Research of Computers, 2017, 34(12): 3603-3606, 3610. (in Chinese) doi: 10.3969/j.issn.1001-3695.2017.12.019
[10] HEIDARI A A, PAHLAVANI P. An efficient modified grey wolf optimizer with Lévy flight for optimization tasks[J]. Applied Soft Computing, 2017, 60: 115-134. doi: 10.1016/j.asoc.2017.06.044
[11] LU C, GAO L, YI J. Grey wolf optimizer with cellular topological structure[J]. Expert Systems With Applications, 2018, 107: 89-114. doi: 10.1016/j.eswa.2018.04.012
[12] 方开泰, 马长兴. 正交与均匀试验设计[M]. 北京: 科学出版社, 2001. FANG Kai-tai, MA Chang-xing. Orthogonal and uniform experiment design[M]. Beijing: Science Press, 2001. (in Chinese)
[13] 梅松华, 盛谦, 冯夏庭. 均匀设计在岩土工程中的应用[J]. 岩石力学与工程学报, 2004, 23(16): 2694-2697. doi: 10.3321/j.issn:1000-6915.2004.16.007 MEI Song-hua, SHENG Qian, FENG Xia-ting. Geotechnical engineering[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(16): 2694-2697. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.16.007
[14] 韩煊, 李宁, Standing J. R.. Peck公式在中国隧道施工地面变形预测中的适用性分析[J]. 岩土力学, 2007, 28(1): 23-28, 35. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200701005.htm HAN Xuan, LI Ning, STANDING J R. An adaptability study of Gaussian equation applied to predicting ground settlements induced by tunneling in China[J]. Rock and Soil Mechanics, 2007, 28(1): 23-28, 35. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200701005.htm
-
期刊类型引用(14)
1. 肖时辉，王峰，何柏延，徐国元，黄伟真，李伟鹏. 基于LHS-WOA-ELM的隧道围岩参数反演分析. 矿业研究与开发. 2025(03): 91-100 . 
百度学术
2. 刘月华，朱庆闯，毕乃晨. 基于量子粒子群优化的支持向量机模型反演土力学参数. 建筑施工. 2025(03): 361-365 . 百度学术
3. 李书兆，孙国栋，申辰，李文逵，罗进华，王教龙. 基于深度学习和地震数据的海上风电场CPT预测研究. 工程地球物理学报. 2025(02): 216-226 . 百度学术
4. 阮永芬，李鹏辉，施虹，吴龙，李飞鹏，肖潇. 基于机器学习的软土压缩模量预测及沉降分析. 成都理工大学学报(自然科学版). 2024(02): 258-268+280 . 百度学术
5. 许晓逸，张贵科，陈世壮，邓韶辉，王如宾，张海龙. 基于优化支持向量回归的磨古倾倒变形体参数反分析研究. 三峡大学学报(自然科学版). 2024(03): 15-22 . 百度学术
6. 李向海，杨玲，魏静. 基于优化多维支持向量机回归模型的土体参数反演. 公路交通科技. 2024(05): 27-34+78 . 百度学术
7. 张坤勇，苏政凯，聂美军，郭楷文，沈小锐. 基于ISSA-ELM的土体参数反演与变形预测方法. 应用基础与工程科学学报. 2024(05): 1434-1448 . 百度学术
8. 李谦，周治刚，邓光宏，刘绪勇，丁晔. 基于钻进参数实时预测土体力学性质的Stacking集成模型. 钻探工程. 2024(S1): 61-69 . 百度学术
9. 闫浩，张吉雄，周楠，时培涛. 基于DA-DE-SVM智能模型的煤岩体SC-CO_2压裂效果预测. 岩土工程学报. 2023(02): 362-368+443 . 本站查看
10. 何晓莹，林潮宁，杜效鹄，袁力，周兴波，刘晓青. 考虑渗流场作用的重力坝坝基物理力学参数有限元智能反演分析. 三峡大学学报(自然科学版). 2023(03): 8-13 . 百度学术
11. 王述红，董福瑞. 基于变形预测和参数反演的山岭隧道围岩稳定性分析. 岩土工程学报. 2023(05): 1024-1035 . 本站查看
12. 高新博，谢渊，姜育科，高玮. 基于微扰动的软土区综合管廊结构与地层参数智能反演研究. 北京交通大学学报. 2023(03): 35-43 . 百度学术
13. 杨鹏，顾春生，董教社，胡洋，陈洋. 基于多元回归模型与信息熵的岩土参数预测. 地下空间与工程学报. 2023(06): 1782-1790+1808 . 百度学术
14. 郑晓亮，谢晓贤，王强. 基于灰狼优化的埋地管道泄漏双波谱定位方法. 仪器仪表学报. 2022(08): 204-214 . 百度学术
其他类型引用(14)
计量
- 文章访问数: 219
- HTML全文浏览量: 32
- PDF下载量: 61
- 被引次数: 28
