时序三维瞬变电磁成像揭示滑坡断裂带水力压裂损伤演化

《Water》：Time-Lapse Three-Dimensional TEM Imaging and Damage Evolution of Hydraulic Fracturing in Landslide Fracture Zones

　　

摘要：滑坡断裂带在边坡失稳中起着关键作用，然而其水力压裂过程难以使用常规地球物理方法进行监测。本研究提出一种时序三维瞬变电磁法（Transient Electromagnetic Method, TEM）成像方法，用于表征中国昭通市小龙洞滑坡断裂带的损伤演化特征。研究人员将重复TEM观测与全球导航卫星系统（GNSS）位移监测和环境测量相结合，通过三维反演重建动态电阻率场。典型的电阻率前兆异常表现为3–8%的升高后快速下降。在水力压裂过程中，电阻率从30–50 Ω·m下降至8–15 Ω·m，反映了渐进式的裂缝连通性。雨季渗透系数达到10-4–10-3 m/s，比旱季高出一到两个数量级。观测到的位移中有相当一部分发生在季节性温度波动和降雨相关水文变化主导的时期内。然而，现有观测数据尚无法定量区分冻融过程、降雨入渗、重力蠕变及其他驱动机制的贡献。时序TEM观测揭示了滑坡断裂带电阻率结构的显著时间变化。低阻异常的发育及其空间连通性与季节性水文变化相关联，表明重复TEM调查在表征地下水分再分布和断裂带演化方面具有潜力。然而，现有TEM调查的时间分辨率不足以识别分钟尺度的前兆信号。所提出的框架为断裂演化提供了定量可视化手段，并为滑坡稳定性评估与早期预警提供了有效的地球物理约束条件。

研究背景与科学问题

滑坡断裂带是边坡失稳的关键控制因素，但其内部水力压裂过程的时空演化难以通过传统点状监测手段有效捕捉。尽管已有研究表明热致开裂、孔隙水压力升高及地质结构对滑坡失稳具有重要作用，但地下裂缝网络的空间异质性和时间演化特征仍缺乏有效的非侵入式表征手段。瞬变电磁法（Transient Electromagnetic Method, TEM）虽已应用于复杂地质调查，但在滑坡断裂带长期动态监测中的应用尚未充分建立。针对这一不足，研究人员在中国云南省昭通市小龙洞滑坡区开展了一项系统性研究，旨在通过时序三维TEM成像，结合全球导航卫星系统（GNSS）位移监测与环境测量，揭示断裂带水力压裂过程中的电阻率响应规律及其与变形演化的耦合关系。

研究方法与技术路线

研究人员采用的主要关键技术包括：基于双线圈600D集成系统的时序TEM数据采集，该系统采用零偏移距共点配置，有效消除了浅层探测盲区，提升了0–30 m深度范围内低阻含水结构的探测分辨率；基于非结构化四面体网格的三维时域有限元正演模拟与反演算法，融合光滑约束与小波稀疏约束以平衡反演结果的平滑性与局部异常增强需求；地形约束下的网格化测线布设，沿坡向布置10条主测线和3条横向连接测线，共187个物理测点，并在密集裂隙区加密点距至0.5 m；此外，还集成了GNSS位移监测站、自动雨量计及温湿度传感器等环境监测设备，建立了多参数协同观测体系。

研究结果

电阻率结构的季节性演变特征：旱季初期，断裂带视电阻率范围为180–320 Ω·m，低阻异常呈孤立斑块状分布。雨季中期，受降雨入渗影响，电阻率急剧下降55–70%，降至45–95 Ω·m，低阻异常沿NE方向逐渐连通形成宽约15 m的连续带。冻融期，裂缝边界3–8 m范围内电阻率上升15–25%，反映冻胀引起的结构性变化。次年旱季电阻率恢复至120–180 Ω·m，但仍较同期低20–40%，表明存在累积性不可逆损伤。

电阻率变化与裂缝孔径的定量关系：基于平行导电模型推导出的理论关系，研究人员利用非线性最小二乘法拟合现场数据，得到岩石基质与裂隙水电阻率之比为185、代表性单元长度为3.2 m，决定系数R2达0.92。电阻率变化率Δρ随裂缝孔径d增大呈非线性增长：d小于0.5 mm时Δρ缓慢增加至0.12–0.20；d在0.5–2.5 mm区间内Δρ加速增长，每增加0.5 mm对应Δρ增加0.15–0.25；d超过2.5 mm后Δρ趋于饱和，最大值约为0.75–0.85。

各向异性电阻率模型与渗透张量演化：研究人员构建了考虑裂缝发育方向的各向异性电阻率模型，将各向同性反演结果与GNSS位移的空间相关性从0.61提升至0.78。裂缝区内各向异性系数λ介于2.5–6.8之间，最大电阻率方向垂直于主裂缝面，最小电阻率方向平行于裂缝走向。基于时序TEM反演结果获得的渗透张量显示：旱季最大渗透系数为2.5×10^-6–8.6×10^-6 m/s；雨季中期增至1.2×10^-4–6.8×10^-4 m/s，提高约1.5个数量级；雨季末期回落至3.5×10^-5–1.2×10^-4 m/s，但仍高于旱季水平，表明水力压裂导致不可逆渗透性增强。

损伤变量与归一化电阻率的关系：研究人员建立了损伤变量D与归一化视电阻率指数R_n之间的线性关系，表达式为D = 0.38R_n + 0.09，可作为早期裂缝损伤的定量指示指标。

讨论与结论

讨论部分指出，裂缝介质电阻率响应的非唯一性是本研究的核心局限之一。除裂缝连通性外，体积含水量、孔隙水电导率、温度、黏土矿物含量及水冰相变等因素均可能影响实测电阻率。因此，研究中识别的低阻异常应首先被解释为地下介质整体电学性质的变化，其与充水裂缝及渐进损伤的关联需结合地质背景与变形观测加以佐证，不能视为水力裂缝扩展的唯一证据。未来研究应整合孔隙水压力、土壤含水量、温度和水电导率等多参数观测，以降低解释不确定性。

研究结论如下：本研究采用时序三维瞬变电磁法成像技术，探究了小龙洞滑坡断裂带内水力压裂与损伤过程的季节性演化特征。结果显示，滑坡断裂带电性结构表现出显著的季节变异性。旱季初期断裂带视电阻率较高，低阻异常呈孤立斑块状分布；雨季入渗后电阻率下降约55–70%，低阻异常沿NE方向逐步连通，最显著的电性异常集中于5–15 m深度，与识别出的滑动面高度吻合。冻融期电阻率在裂缝边界3–8 m范围内上升15–25%。GNSS观测表明，超过电阻率阈值的区域累计位移约为周边稳定区的2.3倍，显示出电性演化与滑坡变形的强耦合关系。定量分析进一步证实，电阻率变化率Δρ与裂缝孔径发展高度相关，决定系数达0.93，可可靠反演0.1–2.5 mm范围内的裂缝孔径。电阻率各向异性比介于2.3–7.1之间，最大电阻率方向始终垂直于主裂缝面。损伤变量D与归一化视电阻率R_n之间存在显著线性关系。渗透张量反演进一步表明，雨季最大渗透系数可达10^-4–10^-3 m/s，较旱季提高一至两个数量级，主渗透方向与主裂缝走向一致，渗透响应滞后降雨事件约2–3天。时序TEM方法在约20 m深度范围内对浅层裂缝网络具有良好的成像能力，水平定位误差小于1 m，垂直误差小于2 m。裂缝传播前持续观测到电阻率先升后降3–8%的典型前兆信号。卸荷阶段部分裂缝自愈合导致残余电阻率降低20–40%。总体而言，所提出的时序三维TEM框架为量化裂缝演化、评估滑坡稳定性及在地质活跃山区支持早期预警系统提供了有效工具。