集成直流电阻率与地震折射成像用于埃及新米尼亚市基岩不连续面填图

《Scientific Reports》：Integrated DC resistivity and seismic refraction imaging for bedrock discontinuity mapping in New Minia city, Egypt

　　

本研究通过整合电阻率层析成像（ERT）、地震折射和直流电（DC）电阻率测深与钻孔数据，解决了埃及上埃及新米尼亚市裂隙石灰岩环境中的地质技术挑战。研究目标是描绘浅层地下非均质性和结构不连续性，并建立地球物理参数与地质技术参数之间的相关性。研究人员应用遗传算法（GA）优化对直流电阻率反演，减少了非唯一性并改善了模型收敛性。17个钻孔（深度10米）为校准和验证提供了约束。结果揭示了三个主要地下序列：（1）含黏土夹层的砂质—粉质土壤，（2）非均质中间沉积物，（3）裂隙石灰岩基岩。主要构造特征，包括南北走向断层，被清晰识别。该集成方法增强了地下表征，减少了解释不确定性，并为复杂碳酸盐岩地形中的灾害评估和城市规划提供了稳健框架。

论文解读文章

**研究背景与问题**
在喀斯特化石灰岩地区，基岩不连续面（如断层、裂隙、洞穴）严重影响工程稳定性，传统钻孔方法难以探测亚垂直结构。埃及新米尼亚市（New Minia city）位于上埃及尼罗河东岸，区域石灰岩基岩高度裂隙化且含大量洞穴，NW-SE和NE-SW走向断层控制排水和地形。随着城市扩张，亟需详细场地调查以评估地基适宜性。然而，单一地球物理方法（如直流电阻率或地震折射）存在分辨率局限：1D电阻率测深受等值性和压制效应影响，地震折射难以识别高速层下的低速薄层。因此，研究人员提出集成直流电阻率测深、2D电阻率层析成像（ERT）和地震折射成像，结合钻孔数据，以降低不确定性并全面描绘基岩不连续面。

**研究内容与结论**
研究人员在11 km2区域内采集了34个施伦贝谢（Schlumberger）垂直电阻率测深点、5条2D-ERT剖面（Wenner Beta阵列，最大电极距5 m，数据层数n=15）和25条地震折射测线（24道，道间距5 m），同时利用17个10 m深钻孔的地质技术数据（包括岩石质量指标RQD、无侧限抗压强度UCS、含水量WC%、标准贯入试验SPT-N值）进行校准。通过MATLAB代码预处理电阻率数据，先后应用常规线性滤波反演和遗传算法（GA）反演（种群大小50，世代数50），GA反演降低了非唯一性并提高了与钻孔的吻合度（如VES15的RMS误差低至3.46）。2D-ERT数据采用Res2dinv和DC2DInvRes软件反演，地震折射数据采用SeisImager/2D和ZondST2D软件进行走时层析成像（RMS误差2.5–3.4）。最终通过Rockworks软件构建3D模型。

**研究结果**
**Bedrock discontinuities from resistivity sounding**（基岩不连续面来自电阻率测深）：通过拼接1D反演电阻率剖面（如A-A’剖面），识别出五个地电层：表层（电阻率251–61,412 Ω·m，厚0.5–1.5 m）、冲积层（50–1,051 Ω·m，厚2–30 m）、裂隙化石石灰岩（1,000–5,995 Ω·m，厚5–30 m）、泥质石灰岩（41–704 Ω·m，厚5–20 m）和致密白色石灰岩（5,034–11,428 Ω·m）。电阻率急剧变化区域对应断层，与钻孔和露头对比一致。3D栅栏模型展示横向和垂向非均质性，受构造控制。

**Imaging of rock discontinuities**（岩石不连续面成像）：2D-ERT（如ERT1-5）和地震折射（如S1-5）联合揭示：高电阻率（>1,500 Ω·m）与高波速（约2.9 km/s）对应致密石灰岩，低电阻率（100–600 Ω·m）与低波速（约1.2 km/s）对应冲积层或泥质石灰岩。在ERT3和S3剖面中，西侧浅部高阻垂向边界与NW-SE断层吻合，并识别出低阻洞穴（如25 m和50 m处，深2 m）。地震折射在高速层下无法探测低速泥质石灰岩（隐藏层问题），而ERT弥补了该缺陷。

**Geotechnical characterization of bedrock and soil**（基岩与土壤的地质技术表征）：相关分析表明：土壤（砂质黏土）电阻率<150 Ω·m，含水量约29%，SPT-N值10–30，波速0.7–1.3 km/s；裂隙石灰岩电阻率1,070–3,946 Ω·m，波速2.034–4.357 km/s，UCS值22–523 kg/cm2，RQD值10%–25%（极差至差）。黏土和细粒沉积物（高塑性指数PI=8–19）集中在北部和西部浅层（1–4 m），存在中等膨胀性风险。

**总结讨论**
研究人员指出，各方法在精度上优势互补：GA反演的DC电阻率测深与钻孔吻合度高（RMS低至3.46）；2D-ERT提供高分辨率横向成像，成功识别断层和洞穴；地震折射约束基岩深度和速度对比，但分辨率在高度非均质区下降。钻孔校准锚定地球物理异常，但10 m深度限制下方主要依赖反演。集成工作流通过迭代解决数据冲突（如电阻率示黏土而地震示石灰岩时以钻孔优先），降低了不确定性，为地基设计、灾害评估和城市规划提供了稳健框架。

**研究结论翻译**
该研究表明，集成直流电阻率测深、2D-ERT和地震折射数据，并利用钻孔信息进行校准，为裂隙石灰岩中地下不连续面填图提供了可靠框架。该组合方法成功描绘了岩性变异和构造特征，包括NW-SE走向断层和可能的落水洞带。电阻率和地震波速对比的证据显示，中部和西部区段下方为相对致密的石灰岩，而东部和北部区段含有细粒土壤（高塑性），构成地质技术风险。基于GA的反演应用通过减少非唯一性并增强与钻孔数据的相关性，提高了电阻率模型的可靠性。总体而言，该集成方法降低了解释不确定性，并为新米尼亚市等快速发展区域的基础设计、灾害评估和城市规划提供了稳健基础。