富水砂层冻结壁形成过程声场响应特征研究

张基伟; 刘书杰; 张松; 李方政; 韩玉福; 王磊

doi:10.11779/CJGE202012009

富水砂层冻结壁形成过程声场响应特征研究 English Version

张基伟^{1, 2, 3,},
刘书杰^{1, 2, 3, 4, ,},
张松^{1, 2, 3, 5},
李方政^{1, 2, 3},
韩玉福^{1, 2, 3},
王磊^{1, 2, 3}

1.
煤炭科学研究总院建井研究分院,北京　100013
2.
北京中煤矿山工程有限公司,北京　100013
3.
天地科技股份有限公司,北京　100013
4.
北京科技大学土木与资源工程学院,北京　100083
5.
石家庄铁道大学土木工程学院,河北　石家庄　050043

基金项目:

国家自然科学基金青年基金项目 51804157

国家自然科学基金项目 51774183

十三五国家重点研发计划 2016YFC0600904

详细信息

作者简介:
张基伟(1988—),男,副研究员,博士,主要从事地下工程特殊施工方法的科研工作。E-mail：357102668@qq.com

通讯作者:
刘书杰, E-mail：125701782@qq.com

中图分类号: TU43
计量
- 文章访问数: 231
- HTML全文浏览量: 19
- PDF下载量: 121
出版历程
- 收稿日期: 2020-02-19
- 网络出版日期: 2022-12-05
- 刊出日期: 2020-11-30

Response characteristics of sound fields of stratum frozen wall of water-rich sand during developing process

ZHANG Ji-wei^{1, 2, 3,},
LIU Shu-jie^{1, 2, 3, 4, ,},
ZHANG Song^{1, 2, 3, 5},
LI Fang-zheng^{1, 2, 3},
HAN Yu-fu^{1, 2, 3},
WANG Lei^{1, 2, 3}

1.
Shaft Branch of China Coal Research Institute, Beijing 100013, China
2.
Beijing China Coal Mining Engineering Co., Ltd., Beijing 100013, China
3.
Tian Di Science & Technology Co., Ltd., Beijing 100013, China
4.
School of Resources & Civil Engineering, University of Science and Technology Beijing, Beijing 100083, China
5.
School of Civil Engineering, Shijiazhuang Tiedao University, Shijiazhang 050043, China

摘要

摘要: 为了掌握富水砂层冻结壁形成过程超声波响应特征,利用NM-4A型非金属声波检测仪对广州地铁三号线北延伸冻结段-16 m层位进行超声波检测。基于此,根据孔隙介质热传导与压力声学理论,提出了人工冻结壁热-声耦合数值模拟方法。获得了富水砂层冻结过程中温度场-声场的演化特征与声学特征值变化规律。掌握了不同冻结锋面位置、孔间未交圈距离与波速之间的定量关系。总结了热-声耦合模型在冻结施工中应用思路与流程。①应用热-声耦合数值模拟方法求解出J₁—J₂检测孔间的波速准确率为89.7%～96.97%；②超声波由冻土向未冻土的传播过程中反射与折射现象非常明显,根据声波理论计算出J₁—J₂检测孔间声波平均反射系数为0.278、折射系数为1.278；③冻结早期在未冻土-冻土交界面附近出现明显的声压增高区域,声压场分布不均匀,且声线发生偏移。冻结后期无声压集中现象且声压场分布较均匀；④冻结壁发育过程中,波速随着冻结时间的增加而增加,而峰值声压、声强均随着冻结时间增加呈先减小后增加的趋势；⑤冻结锋面位置、冻结孔间未交圈距离均与波速之间呈非线性单调递增趋势,满足二次函数关系,相关性很好；⑥总结了热-声耦合模型应用思路与流程,可实现冻结壁发育状况的“早期预报、过程监测、效果评价”。
- 热-声耦合 /
- 热传导 /
- 峰值声压 /
- 声强 /
- 冻结锋面
Abstract: The ultrasonic tests on artificial frozen wall of water-rich sand stratum at depth of 16 m in Guangzhou Metro Line 3 are conducted by using NM-4A nonmetal ultrasonic test meter to study the response characteristics of sound fields during developing process of frozen wall. A thermoacoustic coupling numerical simulation method is proposed based on the heat conduction and pressure acoustic theory. The variation rule of temperature-sound field and acoustic characteristic value are obtained. The relationship among freezing front, unclosed distance and wave velocity is grasped. The application process of the thermoacoustic coupling model is summarized. The results are as follows: (1) The accuracy range of wave velocity is 89.7% to 96.97% betweenJ₁and J₂ holes by the numerical model. (2) The reflection and refraction are obvious when ultrasonic waves transmit from unfrozen soil to frozen soil. The average reflection and refraction coefficients are 0.278 and 1.278 betweenJ₁ and J₂holes. (3) The concentration area of acoustic pressure appears near the interface between unfrozen and frozen soils, and the sound field is inhomogeneous during early freezing period. Conversely, the sound field is homogeneous, and there is no pressure concentration area. (4) The wave velocity increases with the increasing freezing time, while both the acoustic pressure and the sound intensity increase first, then decrease as the freezing time increases. (5) Both the freezing front and the unclosed distance increase with the increasing wave velocity, and they meet quadratic function relations, indicating that the correlation is very good. (6) The application process of the thermoacoustic coupling model is summarized, which can achieve good effect on developing evaluation of frozen wall.
- thermoacoustic coupling /
- heat conduction /
- absolute acoustic pressure /
- sound intensity /
- freezing front

HTML全文

0. 引言

人工冻结法是利用人工制冷技术,将待建结构周围地层冻结为连续封闭冻土体的特殊地层加固方法。该工法具有隔绝地下水、增加土体强度、提高稳定性、无污染等优势,已经被广泛应用于矿井建设、地铁建设、基坑及综合管廊建设中^[1-3]。但是文献[4, 5]均指出人工冻结法在流速较大的富水砂层中应用时,往往会出现冻结壁不交圈、冻结壁发展过慢等现象,容易造成重大安全事故。因此掌握冻结壁发育状况对人工冻结法安全性施工具有重要意义。

目前,广泛应用的冻结壁发育状况检测方法包括图乘法^[6]、经验公式法^[6]、数值模拟法^[7]等。然而,以上3种方法均基于测温孔内实时温度监测结果,在距测温孔较远位置冻结壁发育状况判断中往往会忽略局部冻结壁开窗、强度不足等关键信息,在一定程度上增加了该工法在高渗透性富水砂层中应用风险。近年来,为了整体掌握冻结壁发育状况,宋雷等^[8]提出了地质雷达探测冻结壁发育状况方法。但该方法探测深度较小（3～30 m）,结果受到冻结管（材质一般为无缝钢管）的干扰非常大。

超声波检测技术是通过测定超声透射岩土后的声学信号变化来间接地反映岩土体的物理力学参数、应力状态特征。目前已被广泛应用于工程岩土分类分级、岩土体物理力学性质测定等领域中,具有灵敏度较高、检测周期短、受金属影响较低等优点^[9]。自20世纪70年代,国内外学者长期致力于建立波速与冻土物理力学性质相关关系,将超声波应用于冻土物理力学性质室内测定中。Kurfurst等^[10]提出了超声波声学参数可以反映冻土物理力学性质。Nakano等^[11]通过试验发现冻结砂土和粉土的波速与温度之间符合指数关系。Thimus等^[12]提出利用声波法测定超固结粘土中未冻水含量的方法。Martin等^[13]建立了冻土声学参数与物理力学性质之间的关系,并利用超声波技术评价未冻水含量和冻土的弹性常数。国内学者^[14-17]利用室内试验掌握了不同负温、含水量、密度的声学参数与冻土抗压、抗拉强度、弹模与纵、横波速、振幅衰减及动弹模之间的关系。为超声波检测冻土物理力学性质奠定了坚实的理论基础。

但是,以上室内试验成果只能对比不同负温等级养护条件下的冻土声学参数差异,不能反应超声波在冻土温度梯度场中的传播规律与演化特征。当在人工冻结壁形成过程中进行超声波检测时,声波相当于在温度梯度场中传播。由“费马原理”可知^[18],温度的不均匀分布会导致声波发生“折射效应”,进而导致声学参数随着温度梯度的变化而发生时空演变。可见,目前缺乏对超声波在冻土温度梯度场中的传播规律与演化特征的相关研究工作,制约了超声波法应用于人工冻结壁发育状况判定中的应用。

本文基于对广州地铁三号线北延伸段富水砂层冻结壁超声波原位检测结果,建立了热-声间接耦合数值计算方法。利用该方法对实际工程冻结过程中温度场、声场的演化特征与声学特征值变化规律进行分析。基于该数值模型讨论不同冻结锋面位置、冻结孔间未交圈距离的声学参数响应特征。总结出利用声-热耦合模型判定施工现场冻结壁发育状况的思路及应用流程。揭示了超声波在冻土温度梯度场中的传播规律与演化特征,为超声波法评价富水砂层冻结壁发育状况研究工作提供参考。

1. 工程概况与超声波原位检测

1.1 工程概况

广州市轨道交通三号线北延段机场南站—机场北站区间约58.80 m,局部曲线段采用冻结法加固施工。设计冻结段隧道拱顶覆土厚度7.63～8.08 m,地下水位埋深约1.0 m,暗挖隧道宽12.6 m,高8.82 m,断面面积达90.273 m²,底板埋深约为16.0 m。为了尽量控制冻胀融沉,减小冻结体量,采用垂直+分段保温冻结方法（图1）。其中埋深14.0～20.1 m的细砂层厚度大,水量丰富,为富水砂层。A排—B排设计孔、排间距1800 mm,C排—D排设计孔、排间距2250 mm,设计积极冻结时间60 d。

图 1 地层柱状及冻结孔布置立面图

Figure 1. Elevation of stratigraphic columns and arrangement of freezing holes

下载: 全尺寸图片幻灯片

1.2 超声波检测孔布置

超声波声测孔布置在A排冻结孔附近,故本文选取A排—C排冻结孔分析。如图2所示,由于现场施工条件限制,A排垂直冻结孔间距为1.8～1.96 m。Y₆—Y₁₁为防止冷量不足而在A,B排之间补加的冻结孔。在原测温孔T₈附近设两个超声波检测孔J₁,J₂。其中J₁—T₈测管间距为1360 mm,平行A排冻结孔轴线方向布置,旨在测量不同负温等级下冻土超声波波速。J₁—J₂测管间距为1470 mm,垂直轴线布置在A₂₅—A₂₆冻结孔之间。旨在测量冻土温度梯度场条件下不同冻结锋面位置的超声波波速。利用NM-4A非金属超声波探测装置分别在不同冻结时期进行现场超声波检测（其中J₁—J₂孔间超声波仅在积极冻结期60 d内测量）。

图 2 冻结孔、检测孔平面布置图（局部）

Figure 2. Plane layout of freezing and detection holes

下载: 全尺寸图片幻灯片

1.3 超声波检测参数与结果

在超声波测试前需要利用十字交叉法标定超声波检测仪从发射至接收仪器系统产生的系统延迟时间τ₀。发射电压选用1000 V,采样周期选择1.6 μs。将超声波发射与接收探头分别放入两检测孔J₁,J₂或J₁,T₈中,下放至目标层位-16 m,保证发射与接收探头在同一水平上,盐水作为耦合剂。

现场超声波检测结果如图3所示,垂直于A排冻结孔轴线方向布置的J₁—J₂孔波速-冻结天数变化趋势与平行于轴线方向布置的J₁—T₈孔波速-温度趋势基本一致,均为缓慢增加→急剧增加→缓慢增加,但是波速值差异较大。

图 3 J₁—J₂孔间与J₁—T₈孔间波速变化曲线

Figure 3. Variation of wave velocities of holesJ₁-J₂andJ₁-T₈

下载: 全尺寸图片幻灯片

其中J₁—J₂孔间富水砂土相当于单向冻结模式（冻结壁由冻土侧（J₂孔）向未冻土侧（J₁孔）不断发展）,温度梯度非常明显。在积极冻结前期（0～30 d）,J₁—J₂孔间呈冻土区域小于未冻土区域的特点,故当声波在未冻水分降低中传播时反射、折射现象非常显著,波速上升幅度较小。积极冻结中期（30～42 d）冻结壁逐渐靠近J₁孔,此时孔间冻土区域大于未冻土区。故J₁—J₂孔间平均温度下降幅度较大,大量冰颗粒的增加使冻土弹性性能增强,波速增加梯度较高。积极冻结后期（42～60 d）J₁—J₂孔间温度梯度减小,基本可视为冻土,平均温度低于-5℃,故两孔之间弹性模量、抗压、抗折强度上升梯度下降,波速值上升趋势趋于缓慢。

J₁—T₈孔间相当于恒温养护模式（距离较近的J₁—T₈孔间土体温度基本相同）,孔间的温度梯度可忽略不计。根据Fermat原理^[18],超声波在等温场中一定是按照检测孔间传播时间最短的路径传播,基本不会发生折射效应。在孔间土体温度为0℃以上时,未冻土波速值上升幅度非常小,基本保持不变。孔间冻土温度为0～-5℃,土体波速随温度降低而急剧上升,其原因为达到结冰温度后,土颗粒中部分水变为冰晶,充填砂粒之间的孔隙使颗粒之间胶结在一起,形成连续的整体。同时该时期内未冻土至冻土过程中弹性模量、抗压、抗折强度均急剧升高。孔间冻土温度为-5℃～-20℃,弹性模量、抗压、抗折强度上升梯度下降,波速值上升趋势趋于缓慢,与国内外学者室内试验趋势基本一致^[11-17]。

由以上分析可知,由于受到“折射效应”的影响,声波在冻土温度梯度场中传播与恒温养护冻土中传播规律存在较大差异。证明了恒温养护冻土相关试验成果不能准确反应超声波在冻土温度梯度场中的传播规律与演化特征的结论。

2. 人工冻结壁热-声耦合数值模拟方法

2.1 基本假设

含水人工冻土是由土骨架,冰晶体和未冻水共同组成的多孔介质土体,具有多个声阻抗区域。与此同时,由于人工冻结过程中温度场实时变化导致温度梯度场中声阻抗区位置与面积不断改变。由此可见,人工冻土温度场-声场耦合问题是一个非常复杂的多场耦合计算问题,本研究涉及的物理量进行如下假设：

（1）土体是连续、均匀、各向同性的饱和孔隙介质,孔隙间相互连通,满足Whitaker连续型假设^[19]。声波不受土体粒径不同、介质突变而发生折射现象。

（2）相比于相变潜热忽略水分迁移引起的对流传热、导热通量变化,且不考虑渗流作用。

（3）水冰相变仅在[T_m,0]范围内发生相变,T_m为土体结冰点的温度。

（4）设计冻结壁边界以外视为无限大区域且模型边界对超声波检测无反射影响。

（5）热量与声波均在同一标高延径向传递,忽略热量轴向传递,即可将检测水平简化为二维模型。

2.2 控制方程

根据孔隙介质的热传导原理,均质各向同性体遵循导热方程可由下式表示：

$\begin{matrix} d_{z} ρ C_{p} (\frac{\partial T}{\partial t} + u \cdot \nabla T) + \nabla \cdot q + \nabla \cdot q_{r} = d_{z} Q, \\ q = - d_{z} k \nabla T 。 \end{matrix}}$

(1)

式中ρ为介质密度（kg/m³）； $C_{p}$ 为等效体积热容（kJ/(kg·K)）；T为表面温度（℃）； $u$ 为节点平移运动的速度矢量（m/s）； $q$ 为热传导的热通量（W/m^-2）； $q_{r}$ 为热辐射的热通量（W/m^-2）；Q为热源（W/m^-3）；k为等效导热系数（W/(m·K)）。冻结过程中不考虑辐射换热,故 $q_{r}$ 项取为0；将超声波-16 m检测层位的温度场问题简化为平面二维模型,d_z取1 mm。

式（1）中由于利用Comsol软件自带的相变材料模拟水→冰的相变,将模型等效体积分为水→冰的相变体积分数与土体骨架体积分数两部分,则C_p与相变材料的密度 $ρ_{x}$ 、热容C_x、导热系数k_x分别为

$C_{p} = θ_{g} ρ_{g} C_{g} + θ_{x} ρ_{x} C_{x}$ ,

(2)

$ρ_{x} = θ_{w} ρ_{w} + θ_{i} ρ_{i}$ ,

(3)

$C_{x} = \frac{1}{ρ_{x}} (θ_{w} ρ_{w} C_{w} + θ_{i} ρ_{i} C_{i}) + \frac{\partial α_{m}}{\partial T}$ ,

(4)

其中,

$\begin{array}{l} α_{m} = \frac{1}{2} \frac{θ_{i} ρ_{i} - θ_{w} ρ_{w}}{θ_{i} ρ_{i} + θ_{w} ρ_{w}}, \\ k_{x} = θ_{w} k_{w} + θ_{i} k_{i} 。 \end{array}}$

(5)

式中θ_g,θ_w,θ_i分别为土体骨架、水、冰的体积分数； $ρ_{g}$ , $ρ_{w}$ , $ρ_{i}$ , $ρ_{x}$ 分别为土体骨架、水、冰、相变材料的密度（kg·m^-3）；C_g,C_x,C_w,C_i分别为土体骨架、相变材料、水、冰的热容（kJ/(kg·K)）；k_w,k_i,k_x分别为水、冰、相变材料的导热系数（W/(m·K)）。其中相变温度为0℃,水→冰的转变间隔为10 K,相变潜热为333 kJ/kg。

根据压力声学理论^[20],声波在无损介质中的波动方程可由下式表示：

$\frac{1}{ρ c^{2}} \frac{\partial^{2} p}{\partial t^{2}} + \nabla \cdot (- \frac{1}{ρ} (\nabla ρ - q_{d})) = Q_{m}$ 。

(6)

式中c为声波在介质中的传播速度（km/s）；p为自变量声压（Pa）；t为时间（μs）； $q_{d}$ 为偶极子声源（N/m³）；Q_m为单极子声源（1/s²）。

2.3 热-声耦合模拟方法

首先,利用式（1）进行人工冻结热力学分析,获得不同冻结阶段的温度场。其次,将本文中J₁—T₈孔间测得不同冻结时期无温度梯度冻土的波速与温度非线性关系c (T,t)代入下式获得不同冻结时刻声阻抗场：

$\frac{\partial Z}{\partial t} = \frac{\partial c (T, t)}{\partial t} \cdot ρ$ ,

(7)

式中,Z为介质声阻抗（N·s/m³）,c (T,t)为恒温冻结模式冻土波速与温度之间非线性关系式。

利用PDE自定义模块导入不同冻结时期声阻抗场作为声场分析的初始条件。设置声学边界条件,设置偶极子声源 $q_{d}$ ,利用波动方程（6）进行声压场的时域瞬态计算。最终,可获得不同冻结时期声压场分布。

2.4 热-声参数选取与边界条件

（1）热物理参数选取与温度场边界条件

数值模型根据对-16 m冻结细砂土层位进行二维简化处理,考虑到热力学、声学边界条件求解域长×宽=17710 mm×10740 mm。冻结孔、检测孔位布置如图2所示。模拟冻结时间为100 d。土体骨架、水、冰的热物理参数如表1所示。

表 1 热物理参数

Table 1. Thermophysical parameters

名称		密度/(kg·m^-3)	导热系数/(W/m·K)	热容kJ/(kg·K)
土体骨架温度	-30℃	1870	1.785	1.386
	-1℃	1870	1.785	1.386
	0℃	1870	1.56	2.102
	18℃	1870	0.86	2.102
水		1000	1.56	4.2
冰		910	2.24	2.1

下载: 导出CSV

| 显示表格

根据现场实测,设置模型土体初始温度为23.8℃。所有冻结管外壁温度荷载近似为实测盐水温度（表2）。模型4个边界不受边界以外土体传热影响,等于土体初始温度。

表 2 盐水温度

Table 2. Brine temperatures

时间/d	温度/℃	时间/d	温度/℃
0	8	20	-25
2	-5	30	-27
5	-12	42	-28
10	-18	60	-28
15	-20	70	-22

下载: 导出CSV

| 显示表格

（2）波速-温度关系表达式

根据J₁—T₈检测孔测得的现场土体波速-温度的非线性关系（图3）可以利用线性函数分段拟合：

$c = {\begin{cases} - 0.0031 T + 1.81425 (24 ℃ \leq T \leq 1 ℃) \\ - 0.5403 T + 2.38372 (1 ℃ < T \leq - 4.25 ℃) \\ - 0.02635 T + 4.54634 (- 4.25 ℃ < T < - 21 ℃) \end{cases}$ 。

(8)

相关系数分别为1,1和0.992,可准确应用于声阻抗场计算。

（3）声学相关参数选取与边界条件

由于现场应用KON-HNQ-100纵向换能器进行超声波检测,故声源激励源 $q_{d}$ 选用柱面波辐射,为了符合实际检测所用的超声脉冲波且便于分析,选取单频雷克子波作为声源函数,其函数表达式为

$S (t) = 2 e^{- 2 \frac{t - T_{0}}{T_{0} / 2}} \sin (2 π f_{0} t)$ ,

(9)

式中,t为声时（μs）, $f_{0}$ 为声源中心频率（kHz）,T₀为周期（μs）。根据现场超声波检测参数选取声源中心频率 $f_{0}$ 为50 kHz,则T₀为20 μs。

由于外边界土体基本对声场无反射作用,故选取辐射边界模拟超声波在无限大区域内的传播,消除边界效应。声场分析中单元网格大小由声源决定,在有意义的波长内含有不少于6～8个的空间步长可以满足计算精度,本文选取网格最大尺寸小于4 mm。计算时间步长与声源周期有关,本文选取T₀/10,计算时间范围100*T₀。

3. 温度场-声场数值模拟结果分析

3.1 温度场演化规律

提取测温孔J₁,J₂,T₈所在位置不同时间的温度数据,将模拟结果与现场实测结果进行比较。如图4所示,模拟温度数据与实测温度随时间的变化曲线虽在个别区段有较小偏差,但发展趋势一致。故认为本研究所选岩土热物理性质参数合理,温度计算结果可有效反应工程实际情况。

图 4 检测孔模拟结果与实测结果对比

Figure 4. Comparison between simulated and measured results

下载: 全尺寸图片幻灯片

如图5所示,根据室内试验结果,以结冰温度-0.2℃为界,判定冻结锋面位置。随着冻结时间的推移,冻结锋面不断扩展,冻结20～30 d内J₁—J₂检测孔之间土体由未冻土转变为局部冻土,冻结42 d冻结锋面基本到达J₁附近。冻结60～100 dJ₁—J₂测管之间冻土温度梯度逐渐减小。证明J₁—J₂检测孔之间可视为单向冻结养护模式冻土。J₁—T₈检测孔之间温度基本一致,无温度梯度,可视为恒温冻结模式冻土。

图 5 不同冻结时刻温度场分布云图

Figure 5. Temperature fields at different freezing time

下载: 全尺寸图片幻灯片

3.2 不同冻结时期波速响应特征

由以上声压场分析可见,冻结区内声压扩展面积远大于未冻区,其原因为冻结区内声波传播速度远大于未冻结区。本文利用互相关函数法求解平均波速,其前提为发射、接收超声波换能器同轴线、同检测平面^[21]。发射信号x (t)、接收信号y (t)峰值位置所对应的时间位移τ即为信号在该系统中的传递时间,相关运算如下式所示,该方法具有一定的去噪能力：

$R_{x y} (τ) = \frac{1}{N} \int_{0}^{T} x (t) y (t + τ) d t$ ,

(10)

式中,N为激励脉冲波形中正弦波的数量,T为接收超声波声时（μs）,τ为时间位移。

将数值模拟计算出不同冻结时间（12,24,28,32,36,40,42,52,60 d）J₁—J₂检测孔间平均波速结果与现场检测结果进行比较（图6）。

图 6 J₁—J₂检测孔间波速模拟结果与实测结果对比

Figure 6. Comparison between simulated and measured results of wave velocity of holeJ₁-J₂

下载: 全尺寸图片幻灯片

如图6所示,现场实测的J₁—J₂检测孔波速与数值模拟均随着冻结时间增加呈缓慢增加→急剧增加→缓慢增加的趋势。总体而言,冻结后期模拟精度高于冻结前期。J₁—J₂检测孔波速计算准确率在89.7%～96.97%范围内,基本满足工程精度。证明本文所选声学参数合理,热-声耦合模拟结果可有效反应工程实际超声波检测情况。

3.3 声阻抗场演化规律

如图7所示,不同冻结时期声阻抗在3.3×10⁶～9.4×10⁶ N·s/m³范围内。其中T>0℃的未冻土区声阻抗值范围为3.3×10⁶～4.6×10⁶ N·s/m³,T<0℃（局部冻土和冻土）区域内声阻抗值范围为4.6×10⁶～9.4×10⁶ N·s/m³。温度为0℃～-4℃范围内声阻抗值变化梯度较大,为4.6×10⁶～8×10⁶ N·s/m³。原因为该区域内含未冻土→局部冻土→冻土三类声阻抗差异较大的土体。

图 7 不同冻结时刻声阻抗场分布云图

Figure 7. Acoustic impedance fields at different freezing time

下载: 全尺寸图片幻灯片

由于0℃～-4℃范围内的声阻抗值变化梯度较大,声波在局部冻土内部传播过程中会在冻结锋面附近发生反射和折射现象。根据声波理论可知^[25],声波垂直入射的反射系数r_φ与折射系数t_φ可以通过两种介质的声阻抗进行计算：

$| r_{φ} | = | \frac{Z_{2} - Z_{1}}{Z_{2} + Z_{1}} | = | \frac{ρ_{2} c_{2} - ρ_{1} c_{1}}{ρ_{2} c_{2} + ρ_{1} c_{1}} |$ ,

(11)

$| t_{φ} | = (\frac{ρ_{1}}{ρ_{2}}) | \frac{2 ρ_{2} c_{2}}{ρ_{2} c_{2} + ρ_{1} c_{1}} |$ ,

(12)

式中,Z₁,Z₂分别为未冻土区与冻土区平均声阻抗,（N·s/m³）,ρ₁,ρ₂分别为未冻土区与冻土区平均密度,（kg/m³）,c₁,c₂分别为未冻土区与冻土区平均波速（km/s）。

将未冻区（T>0℃）内平均声阻抗值3.95×10⁶ N·s/m³与冻土区（T<-0.2℃）内平均声阻抗值7×10⁶ N·s/m³代入公式内,求得J₁—J₂检测孔之间平均 $r_{φ}$ , $t_{φ}$ 分别为0.278,1.278。可见,超声波由冻土向未冻土的传播过程中反射与折射现象较明显,阻碍了声波在检测孔间的传递。

3.4 声压场演化规律

以不同冻结时间24,60 d为例,设置声波激励源,得到各个冻结时间不同声时（400,500,600 μs）的声压场分布云图。

如图8所示,冻结24,60 d在相同声时的声压场分布有以下区别：

图 8 不同冻结时刻声场分布云图

Figure 8. Sound fields at different freezing time

下载: 全尺寸图片幻灯片

（1）冻结24 d相比于冻结60 d在J₁,J₂检测孔之间的未冻土-冻土交界面附近出现明显的声压增高区域。其原因为冻结24 d,J₁,J₂检测孔之间大部分为未冻土,未冻土相比于冻土在声学性质上较“软”,使得入射波质点在冻结锋面附近会产生非弹性碰撞^[25],此类碰撞会产生明显的折射、反射现象,声压等值线向温度较低的方向弯曲。在此过程中,反射声波与入射声波相互叠加造成声压在冻结锋面附近的声压增高,声压值基本等于同一声时60 d的两倍。冻结60 d,J₁,J₂检测孔之间已经完全形成冻土,且温度梯度较小,故基本不出现声波折射现象。

（2）冻结24 d相比于冻结60 d声压场分布不均匀,且声线（声压等值线法线方向）发生偏移。该现象符合Fermat声波在温度场中的传播路径选择遵从捷线理论。声波在温度梯度场中传播时,温度的不均匀会导致声速改变,声速的不均匀分布又会改变声线的方向。冻结24 d时,声线沿温度梯度最大的方向传播即冻结孔主面方向传播远快于界面方向传播,故在600 μs时冻结孔主面方向的声波传播至到达J₁孔,而界面方向声波仍然停滞在声波反射区内,最终在冻结前期形成不均匀的声压场分布。随着冻结壁不断发展至60 d时,冻结孔界面的冻土温度不断降低,温度梯度降低,声波折射效应不明显。声线沿冻结孔主面方向偏移程度逐渐减小,故在该期间声压场分布较均匀。

3.5 不同冻结时期声场特征值响应特征

声场特征值主要包括峰值声压、声强,其中峰值声压是指某一点于某一瞬时在超声场中的压强与同一点不存在超声波作用时的静态压强差。声强是指超声波在介质中传播时,其传播方向上单位截面中单位时间通过的声能^[25]：

$I = \frac{P_{t}^{2}}{ρ c}$ ,

(13)

式中,I为声强值（W/m²）,P_t为不同冻结时间峰值声压（接收波声压曲线中最大值）（Pa）,ρ为介质密度（kg/m³）,c为不同冻结时间到达接收孔的声速（km/s）。

如图9所示,峰值声压、声强两个声学参数均随着冻结时间增加呈先减小后增加的趋势。积极冻结前期（0～30 d）时,J₁,J₂检测孔间冻结壁逐渐发育导致声阻抗梯度增加,未冻土-冻土界面产生散乱反射引起超声波衰减的现象。而该时期内波速增加幅度较小,导致其增加峰值声压作用远小于折射产生的峰值声压减小作用。积极冻结中后期（30～60 d）,声速上升梯度增加导致折射的减小峰值声压、声强作用逐渐减小,表现为峰值声压、声强逐渐增加。

图 9 不同冻结时间声压、声强变化曲线

Figure 9. Variation of acoustic pressure and sound intensity at different freezing time

下载: 全尺寸图片幻灯片

4. 热-声耦合模型应用思路

由以上分析可知,本文提出的热-声耦合数值模拟方法可有效反应工程实际温度及超声波检测结果。超声波场声学特征值变化规律为：波速随着冻结时间的增加而增加,而峰值声压、声强均随着冻结时间增加呈先减小后增加的趋势。可见,只有波速可作为判断冻结壁发育状况的敏感声学参数。

4.1 冻结锋面位置与波速关系

在数值模型中测量J₁—J₂检测孔之间冻结锋面（-0.2℃等温线位置）与J₂距离。如图10所示,随着冻结锋面的不断发展,波速与冻结锋面位置呈非线性单调递增趋势,可利用二次函数拟合。相关系数值达到了0.983,可见冻结锋面位置与波速的相关性很强。证明利用超声波波速可以推算出两检测孔之间冻结锋面位置。

图 10 不同冻结锋面位置波速变化曲线

Figure 10. Variation of wave velocity at different freezing fronts

下载: 全尺寸图片幻灯片

4.2 冻结孔间未交圈距离与波速关系

定量计算冻结孔间未交圈距离是冻结壁发育状况评价理论长期未解决的难题之一。基于本文数值模型,将声波激励源设置在图2中A₂₆冻结孔作为发射孔,A₂₅冻结孔作为接收孔,模拟A₂₅,A₂₆两冻结孔交圈过程中的未交圈距离与波速之间的关系。

如图11所示,未交圈距离与波速之间呈非线性单调递增趋势,未交圈距离L>0.961的斜率远小于L<0.961的斜率。其原因为L>0.961时,两冻结孔间未冻土区域占A₂₅,A₂₆孔间距一半以上,声波在未冻土内传播速度远小于冻土（图3）。未交圈距离与波速关系可利用二次函数拟合。相关系数值达到了0.983,可见二者相关性很强。证明利用超声波波速可以推算出两冻结孔间未交圈距离。

图 11 冻结孔间未交圈距离波速变化曲线

Figure 11. Variation of wave velocity at different unclosed distances between two freezing holes

下载: 全尺寸图片幻灯片

4.3 热-声耦合模型应用思路

富水砂层两孔冻结模型如图12所示,A,B为冻结孔,C,D为温度和超声检测孔。总结热-声耦合模型应用思路如下：

图 12 富水砂层两孔间冻结壁发育模型

Figure 12. Development model for frozen wall between two holes of water-rich sand stratum

下载: 全尺寸图片幻灯片

（1）冻结施工准备阶段,进行两孔冻结模型实验（图12）。利用C—D检测孔获得如式（8）中所列出的土体波速与温度回归公式。室内测试冻土热物理参数。

（2）冻结孔钻孔施工后,立即根据冻结孔实际偏斜情况建立热-声耦合数值模型。

（3）在冻结施工前期,利用实测测温孔数据变化对温度场数值模型进行优化,获得准确的热-声耦合预测模型。利用热-声耦合模型可获得正常冻结壁发育状况下冻结区域任意位置与时间的冻结锋面、未交圈距离与计算波速之间定量关系。

（4）在冻结施工过程中,现场定期对任意两孔间（冻结孔间、冻结孔-检测孔）进行超声波测试。将现场检测波速代入冻结锋面位置、未交圈距离与计算波速之间定量关系表达式中。推算该时期孔间冻结锋面位置h₁或未交圈距离h₂。据此分析冻结壁厚度和冻结孔间未交圈距离,结合测温孔数据综合判断出冻结壁发育状况。当出现冻结锋面位置、未交圈距离异常时,及时分析孔间出现的问题,合理地选取补冻结孔、降低盐水温度、注浆堵水等工程对策。

（5）在冻结施工验收阶段,现场可对距离测温孔较远区域的冻结孔间（如A—B）进行逐孔排查检测,评价任意两孔间是否仍然存在由于地层渗流等因素造成冻结壁局部“开窗”现象。结合现场测温数据、泄压孔压力数据、钻孔取芯情况,综合评价冻结效果。

综上,热-声耦合数值模拟方法可作为评价富水砂层冻结壁发育状况的手段之一,可实现冻结壁发育状况的“早期预报、过程监测、效果评价”。

5. 结论

本文基于对广州地铁三号线北延伸段第二冻结段富水砂层冻结壁超声波原位检测结果,提出了热-声间接耦合数值计算方法。利用该方法对实际工程冻结过程中温度场、超声场的演化特征与声学特征值变化规律进行了研究。

（1）现场超声波检测结果可知,声波在温度梯度较大的单向冻结土体中传播时反射、折射现象非常显著。无法应用恒温养护冻土相关试验成果反应声波在冻土温度梯度场中的传播规律与演化特征。

（2）根据孔隙介质的热传导原理与压力声学理论,基于J₁—T₈检测孔测得的波速-温度关系表达式,提出了热-声耦合数值模拟方法,应用现场原位测试结果证明该方法计算结果较准确,J₁—J₂检测孔波速计算准确率在89.7%～96.97%。

（3）超声波由冻土向未冻土的传播过程中反射与折射现象非常明显,根据声波理论计算出J₁—J₂检测孔间声波平均反射系数为0.278、折射系数为1.278。

（4）冻结早期在未冻土-冻土交界面附近出现明显的声压增高区域。声压场分布不均匀,且声线发生偏移。冻结后期声压场分布较均匀,无声压集中现象。

（5）冻结壁发育过程中,波速随着冻结时间的增加而增加,而峰值声压、声强均随着冻结时间增加呈先减小后增加的趋势。

（6）冻结锋面位置、孔间未交圈距离均与波速之间呈非线性单调递增趋势,满足二次函数关系,相关性很好。据此提出了热-声耦合模型应用思路,可实现冻结壁发育状况的“早期预报、过程监测、效果评价”。

图 1 地层柱状及冻结孔布置立面图

Figure 1. Elevation of stratigraphic columns and arrangement of freezing holes

下载: 全尺寸图片幻灯片

图 2 冻结孔、检测孔平面布置图（局部）

Figure 2. Plane layout of freezing and detection holes

下载: 全尺寸图片幻灯片

图 3 J₁—J₂孔间与J₁—T₈孔间波速变化曲线

Figure 3. Variation of wave velocities of holesJ₁-J₂andJ₁-T₈

下载: 全尺寸图片幻灯片

图 4 检测孔模拟结果与实测结果对比

Figure 4. Comparison between simulated and measured results

下载: 全尺寸图片幻灯片

图 5 不同冻结时刻温度场分布云图

Figure 5. Temperature fields at different freezing time

下载: 全尺寸图片幻灯片

图 6 J₁—J₂检测孔间波速模拟结果与实测结果对比

Figure 6. Comparison between simulated and measured results of wave velocity of holeJ₁-J₂

下载: 全尺寸图片幻灯片

图 7 不同冻结时刻声阻抗场分布云图

Figure 7. Acoustic impedance fields at different freezing time

下载: 全尺寸图片幻灯片

图 8 不同冻结时刻声场分布云图

Figure 8. Sound fields at different freezing time

下载: 全尺寸图片幻灯片

图 9 不同冻结时间声压、声强变化曲线

Figure 9. Variation of acoustic pressure and sound intensity at different freezing time

下载: 全尺寸图片幻灯片

图 10 不同冻结锋面位置波速变化曲线

Figure 10. Variation of wave velocity at different freezing fronts

下载: 全尺寸图片幻灯片

图 11 冻结孔间未交圈距离波速变化曲线

Figure 11. Variation of wave velocity at different unclosed distances between two freezing holes

下载: 全尺寸图片幻灯片

图 12 富水砂层两孔间冻结壁发育模型

Figure 12. Development model for frozen wall between two holes of water-rich sand stratum

下载: 全尺寸图片幻灯片

表 1 热物理参数

Table 1 Thermophysical parameters

名称密度/(kg·m^-3) 导热系数/(W/m·K) 热容kJ/(kg·K)

土体骨架温度 -30℃ 1870 1.785 1.386
-1℃ 1870 1.785 1.386
0℃ 1870 1.56 2.102
18℃ 1870 0.86 2.102
水 1000 1.56 4.2
冰 910 2.24 2.1

下载: 导出CSV

表 2 盐水温度

Table 2 Brine temperatures

时间/d 温度/℃ 时间/d 温度/℃

0 8 20 -25
2 -5 30 -27
5 -12 42 -28
10 -18 60 -28
15 -20 70 -22

下载: 导出CSV

参考文献(20)

[1]	ARMAGHANI D J, AMIN M F M, YAGIZ S, et al. Prediction of the uniaxial compressive strength of sandstone using various modeling techniques[J]. Int J Rock Mech Mining Sci 2016, 85: 174-186. doi: 10.1016/j.ijrmms.2016.03.018
[2]	VITEL M, ROUABHI A, TIJANI M, et al. Thermo-hydraulic modeling of artificial ground freezing: application to an underground mine in fractured sandstone[J]. Computers and Geotechnics, 2016, 75: 80-92. doi: 10.1016/j.compgeo.2016.01.024
[3]	杨平, 陈瑾, 张尚贵, 等. 软弱地层联络通道冻结法施工温度及位移场全程实测研究[J]. 岩土工程学报, 2017, 39(12): 2226-2234. doi: 10.11779/CJGE201712011 YANG Ping, CHEN Jin, ZHANG Shang-gui, et al. Whole range monitoring for temperature and displacement fields of cross passage in soft soils by AGF[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(12): 2226-2234. (in Chinese) doi: 10.11779/CJGE201712011
[4]	张基伟, 刘志强, 单仁亮, 等. 复杂地层井筒冻结壁异常状况监测技术研究现状与展望[J]. 煤炭科学技术, 2019, 47(1): 103-109. https://www.cnki.com.cn/Article/CJFDTOTAL-MTKJ201901050.htm ZHANG Ji-wei, LIU Zhi-qiang, SHAN Ren-liang, et al. Review and prospect of abnormal condition of shaft frozen wall monitoring technique in complex formation condition[J]. Coal science and technology, 2019, 47(1): 103-109. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTKJ201901050.htm
[5]	ALZOUBI M A, MADISEH A, HASSANI F P, et al. Heat transfer analysis in artificial ground freezing under high seepage: validation and heatlines visualization[J]. International Journal of Thermal Sciences, 2019, 139: 232-245. doi: 10.1016/j.ijthermalsci.2019.02.005
[6]	胡向东. 直线形单排管冻土帷幕平均温度计算方法[J]. 冰川冻土, 2010, 32(4): 142-149. https://www.cnki.com.cn/Article/CJFDTOTAL-BCDT201004019.htm HU Xiang-dong. Average temperature calculation for the straight single-row-pipe frozen soil wall[J]. Journal of Glaciology and Geocryology, 2010, 32(4): 142-149. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BCDT201004019.htm
[7]	蔡海兵, 程桦, 姚直书, 等. 基于冻土正交各向异性冻胀变形的隧道冻结期地层位移数值分析[J]. 岩石力学与工程学报, 2015, 34(8): 1667-1676. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201508019.htm CAI Hai-bing, CHENG Hua, YAO Zhi-shu, et al. Numerical analysis of ground displacement due to orthotropic frost heave of frozen soil in freezing period of tunnel[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(8): 1667-1676. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201508019.htm
[8]	宋雷, 张小俊, 杨维好, 等. 人工冻结工程地质雷达模型试验研究[J]. 岩土工程学报, 2012(1): 115-122. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201201011.htm SONG Lei, ZHANG Xiao-jun, YANG Wei-hao, et al. Experimental study on GPR model for artificial freezing projects[J]. Journal of Geotechnical Engineering, 2012(1): 115-122. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201201011.htm
[9]	赵建铧, 周崇良. 声波透射法检测对桩基质量的分析与判别[J]. 岩土工程技术, 2012, 26(4): 200-203. doi: 10.3969/j.issn.1007-2993.2012.04.010 ZHAO Jian-hua, ZHOU Chong-liang. Waves penetrating tests of the quality analysis and discretion in pile foundation[J]. Chinese Journal of Geotechnical Engineering, 2012, 26(4): 200-203. (in Chinese) doi: 10.3969/j.issn.1007-2993.2012.04.010
[10]	KURFURST P J. Ultrasonic wave measurements on frozen soils at permafrost temperatures[J]. Canadian Journal of Earth Sciences, 1976, 13(11): 1571-1576.
[11]	NAKANO Y, MARTIN A J, SMITH M. Ultrasonic velocities of the dilatational and shear waves in frozen soils[J]. Water Resources Research, 1972, 8(4): 1024-1030.
[12]	THIMUS J Fr, Aguirre-Puente J, Cohen-Tenoudji FR. Determination of unfrozen water content of an over consolidated clay down to −160℃ by sonic approaches—Comparison with classical methods[C]//International Symposium on Ground Freezing, 1991, Rotterdam.
[13]	MARTIN Christ, PARK Jun-Boum. Ultrasonic technique as tool for determining physical and mechanical properties of frozen soils[J]. Cold Regions Science and Technology, 2009, 3(58): 136-142.
[14]	LI Dong-qing, XING H, FENG M, et al. The impact of unfrozen water content on ultrasonic wave velocity in frozen soils[J]. Procedia Engineering, 2016, 143(1): 1210-1217.
[15]	WAND D Y, ZHU Y L, MA W, et al. Application of ultrasonic technology for physical-mechanical properties of frozen soils[J]. Cold Regions Science and Technology, 2006, 44(1): 12-19.
[16]	杨平, 李强, 郁楚侯. 人工冻土声波参数实验研究[J]. 冰川冻土, 1997, 19(2): 149-153. YANG Ping, LI Qiang, YU Chu-hou. An experimental study on the acoustic wave parameters of artifical frozen soil[J]. Journal of Glaciology and Geocryology, 1997, 19(2): 149-153. (in Chinese)
[17]	黄星, 李东庆, 明锋, 等. 冻结粉质黏土声学特性与物理力学性质试验研究[J]. 岩石力学与工程学报, 2015, 34(7): 1489-1496. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201507022.htm HUANG Xing, LI Dong-qing, MING Feng, et al. Experimental study on acoustic characteristics and physical-mechnical properties off frozen slity clay[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(7): 1489-1496. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201507022.htm
[18]	李言钦, 姬会东, 周俊杰, 等. 声波法测量炉内温度场有限元模拟研究[J]. 动力工程学报, 2014, 34(11): 837-842, 866. https://www.cnki.com.cn/Article/CJFDTOTAL-DONG201411001.htm LI Yan-qin, JI Hui-dong, ZHOU Jun-jie, et al. Finite element simulation on acoustic measurement of the in-furnace temperature field[J]. Journal of Chinese Society of Power Engineering, 2014, 34(11): 837-842, 866. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DONG201411001.htm
[19]	邹大鹏, 吴百海, 卢博, 等. 温度对水饱和孔隙介质压缩波速度的影响[J]. 声学技术, 2008(4): 492-496. https://www.cnki.com.cn/Article/CJFDTOTAL-SXJS200804003.htm ZOU WU Da-peng, HU Bai-hai, LU Bo, et al. Thermal effect on compressional speed of water-saturated porous medium[J]. Technical Acoustics, 2008(4): 492-496. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SXJS200804003.htm
[20]	杜功焕, 朱哲民, 龚秀芬. 声学基础[M]. 南京: 南京大学出版社, 2001. DU Gong-huan, ZHU Zhe-min, GONG Xiu-fen. Basic of Acoustics[M]. Nanjing: Nanjing University Press, 2001. (in Chinese)

施引文献(14)

期刊类型引用(6)

1.	李方政，喻新皓，张基伟. 深部冻结井壁早期温度应力计算方法研究. 金属矿山. 2023(03): 86-93 . 百度学术
2.	王涛，马骏，周国庆，许大晴，季雨坤. 冻土地层三维空间变异性表征及冻结帷幕温度特征值演化过程研究. 岩石力学与工程学报. 2022(10): 2094-2108 . 百度学术
3.	黄宝龙，韩玉福，杨志刚，张松，温汉宏，魏可东，杨宁. 大断面长距离通道冻结效果分析研究. 建井技术. 2022(06): 32-40 . 百度学术
4.	张基伟，李方政，喻新皓，丁航，孔令辉. 深部冻结井筒内壁早期温度-应力场演化特征研究. 煤炭科学技术. 2021(02): 69-76 . 百度学术
5.	汪恩良，任志凤，韩红卫，田雨，胡胜博，刘兴超. 超低温冻结黏土单轴抗压力学性质试验研究. 岩土工程学报. 2021(10): 1851-1860 . 本站查看
6.	鲁先龙，陈湘生，陈曦. 人工地层冻结法风险预控. 岩土工程学报. 2021(12): 2308-2314 . 本站查看

其他类型引用(8)

资源附件(0)

图(12) / 表(2)

计量

文章访问数: 231
HTML全文浏览量: 19
PDF下载量: 121
被引次数: 14

0. 引言
1. 工程概况与超声波原位检测
1.1 工程概况
1.2 超声波检测孔布置
1.3 超声波检测参数与结果
2. 人工冻结壁热-声耦合数值模拟方法
2.1 基本假设
2.2 控制方程
2.3 热-声耦合模拟方法
2.4 热-声参数选取与边界条件
3. 温度场-声场数值模拟结果分析
3.1 温度场演化规律
3.2 不同冻结时期波速响应特征
3.3 声阻抗场演化规律
3.4 声压场演化规律
3.5 不同冻结时期声场特征值响应特征
4. 热-声耦合模型应用思路
4.1 冻结锋面位置与波速关系
4.2 冻结孔间未交圈距离与波速关系
4.3 热-声耦合模型应用思路
5. 结论

0. 引言
1. 工程概况与超声波原位检测
1.1 工程概况
1.2 超声波检测孔布置
1.3 超声波检测参数与结果
2. 人工冻结壁热-声耦合数值模拟方法
2.1 基本假设
2.2 控制方程
2.3 热-声耦合模拟方法
2.4 热-声参数选取与边界条件
3. 温度场-声场数值模拟结果分析
3.1 温度场演化规律
3.2 不同冻结时期波速响应特征
3.3 声阻抗场演化规律
3.4 声压场演化规律
3.5 不同冻结时期声场特征值响应特征
4. 热-声耦合模型应用思路
4.1 冻结锋面位置与波速关系
4.2 冻结孔间未交圈距离与波速关系
4.3 热-声耦合模型应用思路
5. 结论

参考文献(20)

施引文献

资源附件(0)

富水砂层冻结壁形成过程声场响应特征研究 English Version

作者简介: 张基伟(1988—),男,副研究员,博士,主要从事地下工程特殊施工方法的科研工作。E-mail：357102668@qq.com

通讯作者: 刘书杰, E-mail：125701782@qq.com

计量

出版历程

Response characteristics of sound fields of stratum frozen wall of water-rich sand during developing process

0. 引言

1. 工程概况与超声波原位检测

1.1 工程概况

1.2 超声波检测孔布置

1.3 超声波检测参数与结果

2. 人工冻结壁热-声耦合数值模拟方法

2.1 基本假设

2.2 控制方程

2.3 热-声耦合模拟方法

2.4 热-声参数选取与边界条件

3. 温度场-声场数值模拟结果分析

3.1 温度场演化规律

3.2 不同冻结时期波速响应特征

3.3 声阻抗场演化规律

3.4 声压场演化规律

3.5 不同冻结时期声场特征值响应特征

4. 热-声耦合模型应用思路

4.1 冻结锋面位置与波速关系

4.2 冻结孔间未交圈距离与波速关系

4.3 热-声耦合模型应用思路

5. 结论

期刊类型引用(6)

其他类型引用(8)

计量

出版历程

目录

0. 引言

1. 工程概况与超声波原位检测

1.1 工程概况

1.2 超声波检测孔布置

1.3 超声波检测参数与结果

2. 人工冻结壁热-声耦合数值模拟方法

2.1 基本假设

2.2 控制方程

2.3 热-声耦合模拟方法

2.4 热-声参数选取与边界条件

3. 温度场-声场数值模拟结果分析

3.1 温度场演化规律

3.2 不同冻结时期波速响应特征

3.3 声阻抗场演化规律

3.4 声压场演化规律

3.5 不同冻结时期声场特征值响应特征

4. 热-声耦合模型应用思路

4.1 冻结锋面位置与波速关系

4.2 冻结孔间未交圈距离与波速关系

4.3 热-声耦合模型应用思路

5. 结论

作者简介:
张基伟(1988—),男,副研究员,博士,主要从事地下工程特殊施工方法的科研工作。E-mail：357102668@qq.com

通讯作者:
刘书杰, E-mail：125701782@qq.com