基于热-声耦合算法的冻结壁发育状况判断预警方法.pdf

《基于热-声耦合算法的冻结壁发育状况判断预警方法.pdf》由会员分享，可在线阅读，更多相关《基于热-声耦合算法的冻结壁发育状况判断预警方法.pdf（21页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 202010439568.9 (22)申请日 2020.05.22 (71)申请人 北京中煤矿山工程有限公司 地址 100013 北京市朝阳区和平里青年沟 路5号64号楼 (72)发明人 张基伟李方政韩玉福张松 高伟许舒荣崔兵兵王磊 孔令辉郑新赟孙佳 (74)专利代理机构 北京冠榆知识产权代理事务 所(特殊普通合伙) 11666 代理人 王道川 (51)Int.Cl. G06F 30/20(2020.01) G06F 30/13(2020.01) (54)发明名称 一种基于热-声。

2、耦合算法的冻结壁发育状况 判断预警方法 (57)摘要 本发明公开一种基于热-声耦合算法的冻结 壁发育状况判断预警方法， 冻结施工准备阶段， 建立富水砂层两孔冻结的室内物理试验模型， 并 绘制波速-温度关系曲线； 利用波速-温度关系曲 线建立波速-温度的非线性关系表达式； 得到考 虑相变的人工冻土热-声耦合算法； 利用热声耦 合算法求解设计冻结壁发育的平均波速； 计算不 同冻结锋面、 未交圈距离与计算的平均波速之间 定量关系； 通过对比实际测量与计算的冻结锋 面、 未交圈距离判断及预警冻结壁发育状况。 利 用该标准判断及预警任意两孔间冻结壁发育异 常情况， 并及时处理。 可做到对冻结壁交圈状况 。

3、的 “早期预报、 过程监测、 效果评价” 。 权利要求书5页 说明书10页 附图5页 CN 111680403 A 2020.09.18 CN 111680403 A 1.一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其特征在于， 包括如下步 骤： 步骤1、 冻结施工准备阶段， 建立富水砂层两孔冻结的室内物理试验模型， 并绘制波速- 温度关系曲线； 步骤2、 利用步骤1中波速-温度关系曲线建立波速-温度的非线性关系表达式； 步骤3、 得到考虑相变的人工冻土热-声耦合算法； 步骤4、 利用热声耦合算法求解设计冻结壁发育的平均波速； 步骤5、 计算不同冻结锋面、 未交圈距离与计算的平均波速。

4、之间定量关系； 步骤6、 通过对比实际测量与计算的冻结锋面、 未交圈距离判断及预警冻结壁发育状 况。 2.根据权利要求1所述的一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其 特征在于， 在步骤1中 (1-1)选取与现场检测土性热力学参数相同的土样， 在室内进行两孔冻结物理模型实 验， 实验箱周围用聚乙烯材料进行保温； (1-2)其中两冻结孔分别为冻结孔A和冻结孔B， 冻结孔A和冻结孔B的盐水流量、 盐水温 度下降梯度、 制冷系统、 清水循环系统、 土体冻结系统条件与现场情况保持一致； 冻结锋面距离冻结孔A或冻结孔B的厚度为h1， 冻结孔A和冻结孔B之间的未交圈距离为 h2； (1-3。

5、)在平行于冻结孔A和冻结孔B的位置分别布置的检测孔C和检测孔D， 用于超声波 检测； (1-4)在检测孔C和检测孔D、 冻结孔A和冻结孔B内注入盐水， 开机试循环， 当冻结管A和 冻结管B内盐水温度为-25-28后， 打开检测孔C和检测孔D的封头帽， 将探头至于检测孔 C和检测孔D， 下放至检测位置， 同步提升发射与接收探头， 保证发射与接收探头在同一水平 上， 检测不同冻结时期的超声波波速， 每隔1-2小时测一次， 然后绘制波速-温度关系曲线。 3.根据权利要求2所述的一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其 特征在于， 检测孔C和检测孔D的间距小于1.5m， 上、 下相邻两。

6、测点的间距为250600mm。 4.根据权利要求3所述的一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其 特征在于， 在步骤(2)中： (2-1)利用线性函数分段拟合波速-温度变化曲线， 获得波速-温度的非线性关系表达 式， 如式(1)所示： 式(1)中， c为波速， km/s； t为时间， h； a1、 a2、 b1、 b2、 d1、 d2为待定系数； T1、 T2、 T3、 T4为温度 曲线中温度拐点， ； 权利要求书 1/5 页 2 CN 111680403 A 2 (2-2)将波速c(T,t)代入式(1)中获得不同冻结时刻声阻抗场， 如式(2)所示： 式中：为不同冻结时刻声阻抗。

7、场。 5.根据权利要求4所述的一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其 特征在于， 在步骤(3)中： 将含水人工冻土视为由土骨架， 冰晶体和未冻水共同组成的多孔介质土体， 具有多个 声阻抗区域； 均质各向同性体含水人工冻土遵循导热方程可由下式表示： 式中： 为密度； Cp是等效体积热容， kJ/(kgK)； T为表面温度， ； u是节点平移运动的 速度矢量， m/s； q是热传导的热通量， W/m-2； qr是热辐射的热通量， W/m-2； Q是热源， W/m-3； k是 等效导热系数， W/(mK)； 冻结过程中不考虑辐射换热， 故qr项取为0； 将超声波检测层位的 温度场问。

8、题简化为平面二维模型， dz取1mm； 模型等效体积分为水至冰的相变体积分数与土体骨架体积分数两部分， 则式(3)中Cp与 相变材料的密度 x、 热容Cx、 导热系数kx分别为式(4)、 式(5)、 式(6)和式(7)： Cp ggCg+ xxCx (4)； x ww+ ii (5)； 其中： kx wkw+ iki； (7) 式中： g、 w、 i分别为土体骨架、 水、 冰的体积分数； g、 w、 i、 x分别为土体骨架、 水、 冰、 相变材料的密度， kgm-3； Cg、 Cx、 Cw、 Ci分别为土体骨架、 相变材料、 水、 冰的热容， kJ/(kg K)； kw、 ki、 kx分别为水。

9、、 冰、 相变材料的导热系数， W/(mK)； 其中相变温度为0， 水到冰的 转变间隔为10K， 相变潜热为333kJ/kg； 则考虑相变因素的人工冻土热力学分析公式为式(8)所示： 利用式(8)进行人工冻结热力学分析， 获得冻结壁正常发育状况下不同冻结阶段的温 权利要求书 2/5 页 3 CN 111680403 A 3 度场， 即式(1)中 将式(8)求解得出的温度场导入式(1)， 求解c(T,t)， 再代入式(2)， 求解作为声场 分析的初始声阻抗条件； 设置声学边界条件， 声源激励源选用柱面波辐射， 选取单频雷克子波作为声源函数， 其 函数表达式为： 式中， t为声时， s； f0为声。

10、源中心频率， kHz； T0为周期， s； 将以上所述声阻抗与声源边界条件波动方程可由下式表示： 式中： p为自变量声压,Pa； t为时间, s； Qm为单极子声源1/s2； 利用式(10)可以求解出不同冻结时期声压场分布， 即 6.根据权利要求5所述的一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其 特征在于， 在步骤4中： 提取冻结前与冻结后某时刻的声压信号传播曲线， 对比两者发射信号x(t)、 接收信号y (t)峰值位置所对应的时间位移 ， 代入到式(11)中， 求解出信号在冻结壁中的传递时间Rxy ( )； 之后， 将两检测孔距离S,传递时间Rxy( )代入到式(12)求解出正。

11、常冻结壁发育状况下 不同冻结时期的平均波速Vp； 7.根据权利要求6所述的一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其 特征在于， 在步骤5中： 利用热-声耦合算法获得正常冻结壁发育状况下冻结区域任意位置的冻结锋面h1、 未交 圈距离h2与平均波速Vp之间定量关系； 回归冻结壁正常发育条件下不同冻结锋面h1、 未交圈 距离h2与计算平均波速Vp之间定量关系符合二次函数关系式(13)、 (14)： 式中， h1为冻结锋面距离冻结孔厚度， mm； h2为两冻结孔间未交圈距离， mm； f1、 f2、 f3、 g1、 权利要求书 3/5 页 4 CN 111680403 A 4 g2、 。

12、g3为待定系数。 8.根据权利要求7所述的一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其 特征在于， 在步骤6中： 在现场冻结项目中， 冻结壁发育状况超声波检测分为两个检测阶段： 冻结过程阶段和冻结验收阶段； 冻结前： 在超声波检测之前应当对现场冻结孔、 检测孔进行验收： 1)检查其管内是否有 异物堵塞， 同时检查管身是否有裂纹、 弯曲或压扁情况； 2)管子之间应基本上保持平行， 不 平行度应控制在小于1； 3)管道安装完毕后应将上口封闭， 以免落入异物， 致使孔道堵塞； 4)管道高于地表300mm500mm， 保证外露高度宜相同； 冻结过程阶段： 在冻结现场， 设置冻结孔A26和冻。

13、结孔A27、 检测孔J1和测温孔T8； 在冻结现场超声波检测冻结孔A26与检测孔J1波速vp1， 代入到公式(13)计算现场实际 冻结锋面距离H1， 7天/次， 直至冻结壁达到冻结边界； 在冻结现场超声波检测任意两冻结孔如A26和冻结孔A27之间波速vp2， 检测时期为设计 冻结壁交圈日， 设计冻结壁交圈日为205天， 代入到公式(14)计算未交圈距离H2； 冻结验收阶段： 同步提升发射与接收探头， 保证发射与接收探头在同一水平上； 记录冻 结孔A26与检测孔J1波速vp1。 9.根据权利要求8所述的一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其 特征在于， 在步骤6中的判断预警方法。

14、为： 将算法中冻结壁正常发育状况下平均波速Vp1和Vp2代入到式(13)、 (14)， 得到冻结壁正 常发育状况下的两检测孔之间冻结锋面位置h1， 与冻结壁正常发育状况下的冻结孔之间未 交圈距离h2； 将现场检测波速vp1和vp2代入到式(13)、 (14)， 得到现场冻结孔A26与检测孔J1 之间冻结锋面位置H1， 与冻结孔A26和冻结孔A27之间未交圈距离H2； 针对出现现场实测判断的同一冻结时间冻结锋面位置H1比热-声耦合算法得到的冻结 壁正常发育状况下冻结锋面位置h1小30以上时， 判断该检测区域为冻结壁发育高风险区 域； 针对设计冻结壁交圈时间内， 某处现场实测波速vp2计算出两冻结。

15、孔未交圈距离H20， 而热-声耦合算法获得的冻结壁正常发育状况下该时刻冻结孔之间未交圈距离h20， 同时 现场检测其它位置两冻结孔未交圈距离H20， 则判断该检测区域为冻结壁 “开窗” 高风险 区域。 10.根据权利要求9所述的一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 其特征在于， 发育高风险区域和 “开窗” 高风险区域的异常情况为： (a)渗流速度过大或者其它地质因素； (b)附近冻结孔冷量不足、 堵塞、 流量小； (c)附近冻结孔偏斜过大。 逐项排查以上所述可能情况， 并及时进行补救， 补救措施： 如发现情况(a)， 应及时对高风险区进行注浆堵水、 减小渗流措施； 如发现情况(。

16、b)， 应及时增加冷量、 疏导冻结孔内异物； 权利要求书 4/5 页 5 CN 111680403 A 5 如发现情况(c)， 应及时进行补冻结孔或者加大冷量措施。 权利要求书 5/5 页 6 CN 111680403 A 6 一种基于热-声耦合算法的冻结壁发育状况判断预警方法 技术领域 0001 本发明涉及冻结壁发育状况检测技术领域。 具体地说是一种基于热-声耦合算法 的冻结壁发育状况判断预警方法。 背景技术 0002 人工冻结法是利用人工制冷技术， 将待建结构周围地层冻结为连续封闭冻土体的 特殊地层加固方法。 该工法具有隔绝地下水、 增加土体强度、 提高稳定性、 无污染等优势， 已 经被广。

17、泛应用于矿井建设、 地铁建设、 基坑及综合管廊建设中。 但是人工冻结法在流速较大 的富水砂层中应用时， 往往会出现冻结壁不交圈、 冻结壁发展过慢等现象， 容易造成重大安 全事故。 因此掌握冻结壁发育状况对人工冻结法安全性施工具有重要意义。 0003 目前， 国内外广泛应用的冻结壁发育状况检测包括： 图乘法、 经验公式法、 数值模 拟法等， 均能够准确的获得局部冻结壁发育状况。 然而， 以上三种方法均基于测温孔内实时 温度监测结果， 在距测温孔较远位置冻结壁发育状况判断中往往会忽略局部冻结壁开窗、 强度不足等关键信息， 在一定程度上增加了该工法在在地下水流速大于5m/d时富水砂层中 应用风险。 。

18、不能实现对冻结壁交圈状况的 “早期预报、 过程监测、 效果评价” 。 发明内容 0004 为此， 本发明所要解决的技术问题在于提供一种可以在地下水流速大于 5m/d时， 能准确判断和预警的基于热-声耦合算法的冻结壁发育状况判断方法及预警方法。 0005 为解决上述技术问题， 本发明提供如下技术方案： 0006 一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 包括如下步骤： 0007 步骤1、 冻结施工准备阶段， 建立富水砂层两孔冻结的室内物理试验模型， 并绘制 波速-温度关系曲线； 0008 步骤2、 利用步骤1中波速-温度关系曲线建立波速-温度的非线性关系表达式； 0009 步骤3、 。

19、得到考虑相变的人工冻土热-声耦合算法； 0010 步骤4、 利用热声耦合算法求解设计冻结壁发育的平均波速； 0011 步骤5、 计算不同冻结锋面、 未交圈距离与计算的平均波速之间定量关系； 0012 步骤6、 通过对比实际测量与计算的冻结锋面、 未交圈距离判断及预警冻结壁发育 状况。 0013 上述一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 在步骤1 中 0014 (1-1)选取与现场检测土性热力学参数相同的土样， 在室内进行两孔冻物理结模 型实验， 实验箱周围用聚乙烯材料进行保温； 0015 (1-2)其中两冻结孔分别为冻结孔A和冻结孔B， 冻结孔A和冻结孔B的盐水流量、 盐 水温。

20、度下降梯度、 制冷系统、 清水循环系统、 土体冻结系统条件与现场情况保持一致； 0016 冻结锋面距离冻结孔A或冻结孔B的厚度为h1， 冻结孔A和冻结孔 B之间的未交圈 距离为h2； 说明书 1/10 页 7 CN 111680403 A 7 0017 (1-3)在平行于冻结孔A和冻结孔B的位置分别布置的检测孔C 和检测孔D， 用于超 声波检测； 0018 (1-4)在检测孔C和检测孔D、 冻结孔A和冻结孔B内注入盐水， 开机试循环， 当冻结 管A和冻结管B内盐水温度为-25-28后， 打开检测孔C和检测孔D的封头帽， 将探头至于 检测孔C和检测孔D， 下放至检测位置， 同步提升发射与接收探头。

21、， 保证发射与接收探头在同 一水平上， 检测不同冻结时期的超声波波速， 每隔1-2小时测一次， 然后绘制波速-温度关系 曲线。 0019 上述一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 检测孔C 和检测 孔D的间距小于1.5m， 上、 下相邻两测点的间距为250600mm。 0020 上述一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 在步骤(2) 中： 0021 (2-1)利用线性函数分段拟合波速-温度变化曲线， 获得波速-温度的非线性关系 表达式， 如式(1)所示： 0022 0023 式(1)中， c为波速， km/s； t为时间， h； a1、 a2、 b1、 b2、 。

22、d1、 d2为待定系数； T1、 T2、 T3、 T4为 温度曲线中温度拐点， ； 0024 (2-2)将波速c(T,t)代入式(1)中获得不同冻结时刻声阻抗场， 如式(2)所示： 0025 0026式中：为不同冻结时刻声阻抗场。 0027 上述一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 在步骤(3) 中： 0028 将含水人工冻土视为由土骨架， 冰晶体和未冻水共同组成的多孔介质土体， 具有 多个声阻抗区域； 0029 均质各向同性体含水人工冻土遵循导热方程可由下式表示： 0030 0031 式中： 为密度； Cp是等效体积热容， kJ/(kgK)； T为表面温度， ； u是节点平移。

23、运 动的速度矢量， m/s； q是热传导的热通量， W/m-2； qr是热辐射的热通量， W/m-2； Q是热源， W/m -3； k是等效导热系数， W/(mK)； 冻结过程中不考虑辐射换热， 故qr项取为0； 将超声波检测 层位的温度场问题简化为平面二维模型， dz取1mm； 0032 模型等效体积分为水至冰的相变体积分数与土体骨架体积分数两部分， 则式 (3) 说明书 2/10 页 8 CN 111680403 A 8 中Cp与相变材料的密度 x、 热容Cx、 导热系数kx分别为式(4)、 式(5)、 式(6)和式(7)： 0033 Cp ggCg+ xxCx (4)； 0034 x w。

24、w+ ii (5)； 0035 0036其中： 0037 kx wkw+ iki； (7) 0038 式中： g、 w、 i分别为土体骨架、 水、 冰的体积分数； g、 w、 i、 x分别为土体骨架、 水、 冰、 相变材料的密度， kgm-3； Cg、 Cx、 Cw、 Ci分别为土体骨架、 相变材料、 水、 冰的热容， kJ/ (kgK)； kw、 ki、 kx分别为水、 冰、 相变材料的导热系数， W/(mK)； 其中相变温度为0， 水到 冰的转变间隔为10K， 相变潜热为333kJ/kg； 0039 则考虑相变因素的人工冻土热力学分析公式为式(8)所示： 0040 0041 利用式(8)进。

25、行人工冻结热力学分析， 获得冻结壁正常发育状况下不同冻结阶段 的温度场， 即式(1)中 0042将式(8)求解得出的温度场导入式(1)， 求解c(T,t)， 再代入式(2)， 求解作为 声场分析的初始声阻抗条件； 0043 设置声学边界条件， 声源激励源选用柱面波辐射， 选取单频雷克子波作为声源函 数， 其函数表达式为： 0044 0045 式中， t为声时， s； f0为声源中心频率， kHz； T0为周期， s； 0046 将以上所述声阻抗与声源边界条件波动方程可由下式表示： 0047 0048 式中： p为自变量声压,Pa； t为时间, s； Qm为单极子声源1/s2； 0049利用式(。

26、10)可以求解出不同冻结时期声压场分布， 即 0050 上述一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 在步骤4 中： 说明书 3/10 页 9 CN 111680403 A 9 0051 提取冻结前与冻结后某时刻的声压信号传播曲线， 对比两者发射信号x(t)、 接收 信号y(t)峰值位置所对应的时间位移 ， 代入到式(11)中， 求解出信号在冻结壁中的传递时 间Rxy( )； 0052 0053 之后， 将两检测孔距离S,传递时间Rxy( )代入到式(12)求解出正常冻结壁发育状 况下不同冻结时期的平均波速Vp； 0054 0055 上述一种基于热-声耦合算法的冻结壁发育状况判断预。

27、警方法， 在步骤5 中： 0056 利用热-声耦合算法获得正常冻结壁发育状况下冻结区域任意位置与时间的冻结 锋面h1、 未交圈距离h2与平均波速Vp之间定量关系； 冻结壁正常发育条件下不同冻结锋面 h1、 未交圈距离h2与计算平均波速Vp之间定量关系符合二次函数关系式(13)、 (14)： 0057 0058 0059 式中， h1为冻结锋面距离冻结孔厚度， mm； h2为两冻结孔间未交圈距离， mm； f1、 f2、 f3、 g1、 g2、 g3为待定系数。 0060 上述一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 在步骤6 中： 在现 场冻结项目中， 冻结壁发育状况超声波检测分为。

28、两个检测阶段： 冻结过程阶段和冻结验收 阶段； 0061 冻结前： 在超声波检测之前应当对现场冻结孔和检测孔进行验收： 1)检查其管内 是否有异物堵塞， 同时检查管身是否有裂纹、 弯曲或压扁情况； 2) 管子之间应基本上保持 平行， 不平行度应控制在小于1； 3)管道安装完毕后应将上口封闭， 以免落入异物， 致使孔 道堵塞； 4)管道高于地表300mm 500mm， 保证外露高度宜相同； 0062 冻结过程阶段： 0063 在冻结现场， 设置冻结孔A26和冻结孔A27、 检测孔J1和测温孔T8； 0064 在冻结现场超声波检测冻结孔A26与检测孔J1波速vp1， 代入到公式(13) 计算现场 。

29、实际冻结锋面距离H1， 7天/次， 直至冻结壁达到冻结边界； 0065 在冻结现场超声波检测任意两冻结孔如A26和冻结孔A27之间波速vp2， 检测时期为 设计冻结壁交圈日， 设计冻结壁交圈日为205天， 代入到公式(14) 计算未交圈距离H2 0066 冻结验收阶段： 同步提升发射与接收探头， 保证发射与接收探头在同一水平上； 记 录冻结孔A26与检测孔J1波速vp1。 0067 上述一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 在步骤 6中的判 断预警方法为： 0068 将算法中冻结壁正常发育状况下计算平均波速Vp1和Vp2代入到式(13)、 (14)， 得 到冻结壁正常发育状况下。

30、的两检测孔之间冻结锋面位置h1， 与冻结壁正常发育状况下的冻 结孔之间未交圈距离h2； 将现场检测波速vp1和vp2代入到式(13)、 (14)， 得到现场冻结孔A26 说明书 4/10 页 10 CN 111680403 A 10 与检测孔J1之间冻结锋面位置H1， 与冻结孔A26和冻结孔A27之间未交圈距离H2； 0069 针对出现现场实测判断的同一冻结时间冻结锋面位置H1比热-声耦合算法得到的 冻结壁正常发育状况下冻结锋面位置h1小30以上时， 判断该检测区域为冻结壁发育高风 险区域； 0070 针对设计冻结壁交圈时间内， 某处现场实测波速vp2计算出两冻结孔未交圈距离H2 0， 而热-。

31、声耦合算法获得的冻结壁正常发育状况下该时刻冻结孔之间未交圈距离h20， 同时现场检测其他位置两冻结孔未交圈距离H20， 则判断该检测区域为冻结壁 “开窗” 高 风险区域。 0071 上述一种基于热-声耦合算法的冻结壁发育状况判断预警方法， 0072 发育高风险区域和 “开窗” 高风险区域的异常情况为： 0073 (a)渗流速度过大或者其它地质因素； 0074 (b)附近冻结孔冷量不足、 堵塞、 流量小； 0075 (c)附近冻结孔偏斜过大。 0076 逐项排查以上所述可能情况， 并及时进行补救， 补救措施： 0077 如发现情况(a)， 应及时对高风险区进行注浆堵水、 减小渗流措施； 0078。

32、 如发现情况(b)， 应及时增加冷量、 疏导冻结孔内异物； 0079 如发现情况(c)， 应及时进行补冻结孔或者加大冷量措施。 0080 本发明的技术方案取得了如下有益的技术效果： 0081 本发明的一种基于热-声耦合算法的冻结壁发育状况判断及预警方法， 能够通过 对平行于冻结孔布置的两检测孔超声波检测结果， 获得波速-温度关系表达式。 进而获得 热-声耦合公式。 增加激励源之后， 求解正常发育状况下不同冻结时期的声学参数曲线作为 冻结壁发育状况判断标准。 现场检测孔波速计算准确率在89.796.97范围内， 表明 热-声耦合算法计算波速结果可有效反应工程实际超声波检测情况。 利用该标准判断及。

33、预 警任意两孔间冻结壁发育异常情况， 并及时处理， 该方法能够准确判断冻结壁发育状况， 特 别是在地下水流速大于5m/d时的富水砂层中， 可以对局部冻结壁发育状况， 如高风险区域 的冻结壁开窗、 强度不足等问题进行及时地判断和预警， 预报准确率可达到80以上， 减少 了在地下水流速大于5m/d时富水砂层中应用风险， 可做到对冻结壁交圈状况的 “早期预报、 过程监测、 效果评价” ， 为冻结工程施工保驾护航。 附图说明 0082 图1本发明富水砂层两孔间冻结壁发育模型结构示意图； 0083 图2恒温养护模式冻土波速-温度变化曲线； 0084 图3冻结孔、 检测孔平面布置图(局部)； 0085 图。

34、4检测孔模拟结果与实测结果对比 0086 图5不同冻结时刻温度场分布云图： (a)为24d， (b)为42d， (c)为 60d； 0087 图6不同冻结时刻声阻抗场分布云图： (a)为24d， (b)为42d， (c) 为60d； 0088 图7不同冻结时刻声场分布云图： (a)24d(400 s)； (b)24d(500 s)； (c)24d(600 s)； (d)60d(400 s)； (e)60d(500 s)； (f)60d(600 s)； 0089 图8检测孔间波速模拟结果与实测结果对比； 说明书 5/10 页 11 CN 111680403 A 11 0090 图9不同冻结锋面位。

35、置波速变化曲线。 具体实施方式 0091 本实施例中的基于热-声耦合算法的冻结壁发育状况判断预警方法应用到广州市 轨道交通三号线北延段机场南站机场北站区间约58.80m， 局部曲线段采用冻结法加固施 工。 0092 1、 背景： 广州市轨道交通三号线北延段机场南站机场北站区间约 58.80m， 局部 曲线段采用冻结法加固施工。 设计冻结段隧道拱顶覆土厚度7.638.08m， 地下水位埋深约 1.0m， 暗挖隧道宽12.6m， 高8.82m， 断面面积达90.273m2， 底板埋深约为16.0m。 为了尽量控 制冻胀融沉， 减小冻结体量， 采用垂直+分段保温冻结方法。 其中埋深14.020.1m。

36、的细砂层 厚度大， 水量丰富， 为富水砂层。 A排B排设计孔、 排间距1800mm， C 排D排设计孔、 排间距 2250mm， 设计积极冻结时间60d。 0093 超声波声测孔布置在A排冻结孔附近， 故本实施例中选取A排-C排冻结孔分析。 如 图3所示， 由于现场施工条件限制， A排垂直冻结孔间距为1.81.96m之间。 Y6-Y11为防止冷 量不足而在A、 B排之间补加的冻结孔。 在原测温孔T8附近设两个超声波检测孔J1、 J2。 J1-J2 测管间距1470mm， 垂直轴线布置在A25A26冻结孔之间。 旨在测量冻土温度梯度场条件下 不同冻结锋面位置的超声波波速。 利用NM-4A非金属超。

37、声波探测装置分别在不同冻结时期 进行现场超声波检测(其中J1-J2孔间超声波仅在积极冻结期60d内测量)。 0094 2、 建立室内试验模型， 包括如下步骤： 0095 步骤1、 冻结施工准备阶段， 建立富水砂层两孔冻结的室内物理试验模型， 并绘制 波速-温度关系曲线； 0096 (1-1)选取与现场检测土性热力学参数相同的土样， 在室内进行两孔冻物理结模 型实验， 实验箱周围用聚乙烯材料进行保温； 如图1所示； 0097 (1-2)其中两孔分别为冻结孔A和冻结孔B， 冻结孔A和冻结孔B 的盐水流量、 盐水 温度下降梯度、 制冷系统、 清水循环系统、 土体冻结系统条件与现场情况保持一致， 即与。

38、广 州市轨道交通三号线北延段机场南站机场北站区间约58.80m处的情况保持一致； 0098 冻结锋面距离冻结孔A或冻结孔B的厚度为h1， 冻结孔A和冻结孔 B之间的未交圈 距离为h2； 0099 (1-3)在平行于冻结孔A和冻结孔B的位置分别布置的检测孔C 和检测孔D， 用于超 声波检测； 0100 (1-4)在检测孔C和检测孔D、 冻结孔A和冻结孔B内注入盐水， 开机试循环， 当冻结 管A和冻结管B内盐水温度为-25-28后， 打开检测孔C和检测孔D的封头帽， 将探头至于 检测孔C和检测孔D， 下放至检测位置， 同步提升发射与接收探头， 保证发射与接收探头在同 一水平上， 检测不同冻结时期的。

39、超声波波速， 每隔1-2小时测一次， 然后绘制波速- 温度关 系曲线， 如图2所示。 检测孔C和检测孔D的间距小于1.5m， 上、 下相邻两测点的间距为250 600mm； 0101 步骤2、 利用步骤1中波速-温度关系曲线建立波速-温度的非线性关系表达式； 0102 利用线性函数分段拟合图2波速-温度变化曲线， 获得波速-温度的非线性关系表 达式。 说明书 6/10 页 12 CN 111680403 A 12 0103 0104 式(1)中， c为波速， km/s； t为时间， h； a1、 a2、 b1、 b2、 d1、 d2为待定系数； T1、 T2、 T3、 T4为 温度曲线中温度拐。

40、点， ； 0105 则利用图2中的恒温养护模式冻土波速-温度变化曲线拟合的波速- 温度的非线 性关系表达式如下式所示： 0106 0107(2-2)不同冻结阶段的温度场计算， 即式(1)中通过图4可以看出模拟温度数 据与实测温度随时间的变化曲线虽在个别区段有较小偏差， 但发展趋势一致。 考虑相变的 多孔介质传热算法计算出温度结果可有效反应工程实际情况， 得到如图5所示的不同冻结 时刻温度场分布云图。 0108不同冻结阶段的声阻抗场计算， 即求解如图6所示， 0109 将波速c(T,t)代入式(1)中获得不同冻结时刻声阻抗场， 如式(2) 所示： 0110 0111式中：为不同冻结时刻声阻抗场。。

41、 0112 步骤3、 得到考虑相变的人工冻土热-声耦合算法： 不同冻结时期声压场分布(以冻 结24d和60d为例)， 即如图7所示的不同冻结时刻声场分布云图； 0113 将含水人工冻土视为由土骨架， 冰晶体和未冻水共同组成的多孔介质土体， 具有 多个声阻抗区域； 0114 均质各向同性体含水人工冻土遵循导热方程可由下式表示： 0115 0116 式中： 为密度； Cp是等效体积热容， kJ/(kgK)； T为表面温度， ； u是 说明书 7/10 页 13 CN 111680403 A 13 0117 故qr项取为0； 将超声波检测层位的温度场问题简化为平面二维模型， dz取1mm； 0118。

42、 模型等效体积分为水至冰的相变体积分数与土体骨架体积分数两部分， 则式(3) 中Cp与相变材料的密度 x、 热容Cx、 导热系数kx分别为式(4)、 式(5)、 式(6)和式(7)： 0119 Cp ggCg+ xxCx (4)； 0120 x ww+ ii (5)； 0121 0122其中： 0123 kx wkw+ iki； (7) 0124 式中： g、 w、 i分别为土体骨架、 水、 冰的体积分数； g、 w、 i、 x分别为土体骨架、 水、 冰、 相变材料的密度， kgm-3； Cg、 Cx、 Cw、 Ci分别为土体骨架、 相变材料、 水、 冰的热容， kJ/ (kgK)； kw、 。

43、ki、 kx分别为水、 冰、 相变材料的导热系数， W/(mK)； 其中相变温度为0， 水到 冰的转变间隔为10K， 相变潜热为333kJ/kg； 0125 则考虑相变因素的人工冻土热力学分析公式为式(8)所示： 0126 0127 利用式(8)进行人工冻结热力学分析， 获得冻结壁正常发育状况下不同冻结阶段 的温度场， 即式(1)中 0128将式(8)求解得出的温度场导入式(2)， (1)， 求解c(T,t)， 再代入式(2)， 求解 作为声场分析的初始声阻抗条件； 0129 设置声学边界条件， 声源激励源选用柱面波辐射， 选取单频雷克子波作为声源函 数， 其函数表达式为： 0130 0131。

44、 式中， t为声时， s； f0为声源中心频率， kHz； T0为周期， s； 0132 将以上所述声阻抗与声源边界条件波动方程可由下式表示： 0133 0134 式中： p为自变量声压,Pa； t为时间, s； Qm为单极子声源1/s2； 0135利用式(10)可以求解出不同冻结时期声压场分布， 即如图7 所示的不同冻结 说明书 8/10 页 14 CN 111680403 A 14 时刻声场分布云图。 0136 步骤4、 利用热声耦合算法求解设计冻结壁发育的平均波速； 0137 提取冻结前与冻结后某时刻的声压信号传播曲线， 对比两者发射信号 x(t)、 接收 信号y(t)峰值位置所对应的时。

45、间位移 ， 代入到式(11)中， 求解出信号在冻结壁中的传递时 间Rxy( )； 0138 0139 之后， 将两检测孔距离S,传递时间Rxy( )代入到式(12)求解出正常冻结壁发育状 况下不同冻结时期的平均波速Vp； 0140 0141 如图8所示， J1-J2检测孔波速计算准确率在89.796.97范围内， 基本满足工 程精度。 证明热-声耦合算法计算波速结果可有效反应工程实际超声波检测情况。 0142 步骤5、 计算不同冻结锋面、 未交圈距离与计算的平均波速之间定量关系； 0143 利用热-声耦合算法获得正常冻结壁发育状况下冻结区域任意位置与时间的冻结 锋面h1、 未交圈距离h2与平均。

46、波速Vp之间定量关系； 回归冻结壁正常发育条件下不同冻结锋 面h1、 未交圈距离h2与计算平均波速Vp之间定量关系符合二次函数关系式(13)、 (14)： 0144 0145 0146 式中， h1为冻结锋面距离冻结孔厚度， mm； h2为两冻结孔间未交圈距离， mm； f1、 f2、 f3、 g1、 g2、 g3为待定系数。 0147 利用热-声耦合算法后的式(13)、 式(14)的计算方法计算出不同冻结锋面、 未交圈 距离与计算波速之间定量关系如(15)， (16)所示。 0148 0149 0150 步骤6、 冻结壁发育状况的判断及预警。 0151 冻结过程阶段和冻结验收阶段； 0152。

47、 冻结前： 在超声波检测之前应当对冻结孔现场冻结孔、 检测孔进行验收： 1)检查其 管内是否有异物堵塞， 同时检查管身是否有裂纹、 弯曲或压扁等情况； 2)管子之间应基本上 保持平行， 不平行度应控制在小于1； 3)管道安装完毕后应将上口封闭， 以免落入异物， 致 使孔道堵塞； 4)管道高于地表300mm500mm， 保证外露高度宜相同； 0153 冻结过程阶段： 如图3所示， 针对冻结孔A26和冻结孔A27、 检测孔J1和测温孔T8； 0154 在冻结现场超声波检测冻结孔A26与检测孔J1波速vp1， 代入到公式 (13)计算现场 实际冻结锋面距离H1， 7天/次， 直至冻结壁达到冻结边界；。

48、 0155 在冻结现场超声波检测任意两冻结孔如A26和冻结孔A27之间波速 vp2， 检测时期 为设计冻结壁交圈日， 设计冻结壁交圈日为205天， 代入到公式(14)计算未交圈距离H2； 说明书 9/10 页 15 CN 111680403 A 15 0156 冻结验收阶段： 同步提升发射与接收探头， 保证发射与接收探头在同一水平上； 记 录冻结孔A26与检测孔J1波速vp1。 0157 对于波速异常、 通行时间与其它同水平检测结果不一致情况， 需要及时分析地质 原因， 排查空洞、 地下水流速等原因， 达到早期预报的目的。 待检测完毕， 应当及时拧紧冻结 孔封头帽， 并检查封头帽垫圈封水情况，。

49、 防止盐水泄露。 声测孔应当也设置封头帽， 防止异 物落入孔中。 0158 在步骤6中的判断预警方法为： 0159 将算法中冻结壁正常发育状况下平均波速Vp1和Vp2代入到式(13)、 (14)， 得到冻 结壁正常发育状况下的两检测孔之间冻结锋面位置h1， 与冻结壁正常发育状况下的冻结孔 之间未交圈距离h2； 将现场检测波速vp1和vp2代入到式(13)、 (14)， 得到现场冻结孔A26与检 测孔J1之间冻结锋面位置H1， 与冻结孔A26和冻结孔A27之间未交圈距离H2； 0160 针对出现现场实测判断的同一冻结时间的冻结锋面位置H1比设计冻结锋面位置h1 小30以上时， 判断该检测区域为冻。

50、结壁发育高风险区域； 0161 针对设计冻结壁交圈时间内， 某处现场实测波速vp2计算出两冻结孔未交圈距离H2 0， 而热-声耦合算法获得的冻结壁正常发育状况下该时刻冻结孔之间未交圈距离h20， 同时现场检查其它位置两冻结孔未交圈距离H20， 则判断该检测区域为冻结壁 “开窗” 高 风险区域。 0162 发育高风险区域和 “开窗” 高风险区域的异常情况为： 0163 (a)渗流速度过大或者其它地质因素； 0164 (b)附近冻结孔冷量不足、 堵塞、 流量小等； 0165 (c)附近冻结孔偏斜过大。 0166 逐项排查以上所述可能情况， 并及时进行补救， 补救措施： 0167 如发现情况(a)，。