集成ERT、岩石物理与钻孔测井技术在埃及拉斯埃尔赫克马海岸带岩土分层与可持续海岸开发中的应用

《Scientific Reports》：Integrated ERT, petrophysics and borehole logging for geotechnical zonation and sustainable coastal development in Ras El-Hekma, Egypt

　　

埃及西北部地中海沿岸的快速城市扩张需要具备可靠浅地表表征能力的研究，以同时支撑岩土安全性与可持续地下水管理。研究人员整合87条电气电阻率层析成像（Electrical Resistivity Tomography, ERT）剖面与包括自然伽马（Gamma Ray, GR）、自然电位（Spontaneous Potential, SP）及粒度分析在内的钻孔测井数据，为埃及拉斯埃尔赫克马（Ras El-Hekma）海岸带建立高分辨率岩土与水文地层框架。集成数据集分辨出四个具有不同电阻率特征、岩性组成及工程特性的主要浅地表单元，并通过共址钻孔记录与实验室粒度分析进行了标定。统计分析显示电阻率与孔隙度（r ≈ ?0.52；厚度校正后 r ≈ ?0.85）及黏土含量（r ≈ ?0.78；95% CI (?0.85, ?0.69)）之间存在强负相关关系，而GR值与黏土百分比呈强正相关（r ≈ 0.95；p < 0.001）。第三电性层以富含黏土的上新世页岩为主，是最具问题的岩土单元，其电阻率值低于12 Ω·m，黏土含量超过34%，厚度约在1.5至3.5 m之间。空间分析表明存在与埋藏鲕粒脊、萨布哈（sabkha）环境及局部盐水入侵相关的显著NE–SW向非均质性。基于集成岩石物理特征，研究区被划分为四个承载力分区（A–D），每个分区均关联特定的地基与地基改良建议。所提出的工作流程通过结合地球物理成像、钻孔标定与统计解译，为干旱海岸区域提供了一个可迁移且具有成本效益的框架，以支持可持续城市规划、基础设施开发与浅地表风险评估。

研究背景与意义

埃及西北部地中海沿岸目前正处于加速城市化进程中，拉斯埃尔赫克马地区作为重要的战略发展走廊，正在开展大规模的住宅、旅游及投资项目。然而，该地区浅地表具有极强的横向非均质性，由未固结碳酸盐砂、鲕粒灰岩脊、萨布哈沉积、裂隙碳酸盐岩及富黏土上新世单元构成。这种变异性引发了差异沉降、土壤膨胀性、岩溶化、盐水入侵及含水层脆弱性等重大岩土与水文地质问题。传统仅依赖孤立钻孔与室内试验的岩土调查难以捕捉海岸沉积环境中典型的强烈横向非均质性，可能遗漏不连续黏土层、溶蚀带及盐羽流等关键特征；而缺乏地质标定的地球物理调查则可能导致非唯一性解译，限制工程适用性。既往沿埃及西北海岸的研究多聚焦水文地质解译，未建立直接适用于城市规划与地基工程的集成岩土分类系统，特别是拉斯埃尔赫克马区域内膨胀性黏土电性层的横向连续性、厚度变化及其工程意义尚未明确约束，导致保守的地基设计与建设成本增加。因此，研究人员通过开展集成高密度地球物理成像与标定钻孔及粒度数据集的研究，旨在生成实用且可复现的海岸岩土评估工作流程，为类似环境的可持续开发提供科学基础。该论文发表于《Scientific Reports》。

主要关键技术方法

研究人员在拉斯埃尔赫克马地区布设87个ERT测点与共址钻孔，获取浅至约40 m深度的数据集。采用温纳-施卢贝格尔（Wenner-Schlumberger）阵列进行二维ERT测量，使用稳健平滑约束最小二乘算法进行电阻率反演。钻孔数据涵盖SP、GR、地层厚度及岩性描述，采集分层土样进行孔隙度与粒度分布（砂%、粉砂%、黏土%）实验室分析，计算伽马射线指数（IGR）以归一化天然放射性。基于阿尔奇方程（Archie’s law）建模电阻率-孔隙度关系，并采用线性回归模型根据GR估算黏土体积（V_clay）。所有数据集集成至地理空间数据库，运用Pearson相关矩阵、层状分布统计及双变量可视化进行分析，利用Python科学库进行数据处理与空间分布制图，通过Moran’s I指数评估空间自相关性。

研究结果

空间分布模式

研究人员通过地理空间分布图揭示显著的横向非均质性。第四电性层高电阻率异常（>1000 Ω·m）集中于研究区东北部，对应于灰岩基底高程较高处；第三黏土电性层厚度在1.5 m至3.5 m间变化且无简单东西向梯度，暗示受下伏灰岩古地形控制；第一表层最大厚度（达15 m）见于NW-SE向古干谷通道。空间自相关分析（Moran’s I = 0.72，p < 0.01）表明500 m滞后距下采样密度充足，支持空间插值可靠性，指示横向相变及可能的构造控制影响地基稳定性与含水层分隔。

一般统计关系

基于87组钻孔-ERT配对数据的统计分析显示明确岩石物理模式与强参数间相关性。厚度与电阻率呈中度正相关（r ≈ 0.62），厚度与孔隙度呈强负相关（r ≈ ?0.85），电阻率与孔隙度呈显著负相关（原始值 r ≈ ?0.52，厚度校正后 r ≈ ?0.85，95% CI [?0.89, ?0.80]），符合Archie定律。SP、GR与IGR相互呈极强正相关（r > 0.84）；GR及IGR与黏土含量呈极强正相关（r ≈ 0.95），与砂含量呈强负相关（r ≈ ?0.95），验证GR测井可作为本地沉积序列中黏土体积极佳代理指标。经中心对数比（centred log-ratio, CLR）变换后粉砂与黏土呈强负相关（r = ?0.73，p < 0.001），修正了原始组分数据的虚假相关。

按电性地层的地层区分

箱线图显示四层间统计显著差异。第一表层（全新统未固结碳酸盐砂）电阻率41.5–1408 Ω·m，孔隙度32.0%–37.4%，砂含量96.1%–98.0%，黏土1.2%–2.4%，GR 24.0–30.8 API，厚度1.5–15.0 m（均值≈6.5 m）。第二电性层（浅海鲕粒灰岩脊）电阻率202–557 Ω·m，孔隙度37.2%–39.1%，黏土0.44%–0.70%，GR最低18.5–20.4 API，厚度0.75–7.0 m（均值≈3.5 m）。第三电性层（富黏土上新世页岩）电阻率极低（2–12 Ω·m，均值≈9 Ω·m），孔隙度18.3%–19.3%，砂含量降至39.7%–42.9%，黏土升至34.3%–36.2%，GR极高（99.5–101.5 API），IGR 0.80–0.82，SP≈0 mV，厚度1.5–3.5 m（均值≈2.5 m）。第四电性层（裂隙中新世灰岩）电阻率最高（816–1287 Ω·m，多>1000 Ω·m），孔隙度最低（9.4%–11.4%），砂含量87.7%–89.6%，黏土6.2%–7.4%，GR 48.0–52.3 API，厚度13.5–33.8 m（均值≈24 m）。

交叉图分析

电阻率-孔隙度散点图确认各层间强反比关系：第四电性层聚于高电阻率（>800 Ω·m）低孔隙度（~10%），第三电性层位于极低电阻率（2–12 Ω·m）中等孔隙度（18–19%），第一表层居中（40–1400 Ω·m，32–37%）。砂与黏土散点图清晰分离三组分场：第三层富黏土区（砂<43%，黏土>34%），第二层洁净碳酸盐区（砂>98.8%，黏土<0.7%），第一与第四层混合砂主导区（砂87–98%，黏土1–7.5%）。GR与黏土含量呈极佳线性相关（r ≈ 0.95，p < 0.001），尤其第三层GR>99 API对应黏土>34%，该关系独立于实验室粒度分析，验证了天然伽马测井用于非侵入式黏土量化的可靠性。SP-电阻率图显示层间电化学特征：第一与第四层SP负值为?20至?22 mV及?9至?11 mV，反映渗透性与流体流动性；第三层SP≈0 mV，符合低渗透富黏土层受限电化学活动。

讨论部分总结

研究人员在讨论中指出，岩石物理关系与Archie定律一致，但强调Archie定律严格适用于洁净无黏土地层，含黏土沉积需考虑表面电导，故对洁净单元（层1、2、4）采用经验标定胶结指数（m），层3电阻率视为黏土定性代理而非定量孔隙度估计。层1同孔隙度下电阻率系统性低于层4归因于孔隙水盐度差异及细粒中导电黏土矿物表面传导，层2高孔隙度（37–39%）中等电阻率（202–557 Ω·m）与层4低孔隙度高电阻率差异体现m与迂曲度（a）具岩性依赖性。GR-黏土近完美线性（r ≈ 0.95）源于研究区放射性强几乎专属于伊利石与蒙脱石，无钾长石等干扰，使V_clay= 0.88 × IGR + 0.05模型具备高实用性，可从测井快速量化黏土以评估膨胀土灾害。基于集成数据提出四级承载力分区：A区（ρ > 1000 Ω·m，黏土<2.5%）适宜浅基础但需ERT排查岩溶；B区（ρ = 200–1000 Ω·m，黏土2.5–6.5%）推荐筏板基础并配排水；C区（ρ = 50–200 Ω·m，黏土+粉砂10–20%）需桩基或深层地基改良及抗硫酸盐混凝土；D区（ρ < 50 Ω·m，黏土>34%）为膨胀黏土危害区，需锚入A区深桩、格构板或石灰-粉煤灰化学改良。空间上高阻呈NE–SW不连续集群反映鲕粒脊与萨布哈沉积，层1与层2低阻异常（<50 Ω·m）初步指示局部盐水入侵。研究虽限于浅至中深域且受海岸湿度盐度季节变化影响，但通过87点标定显著提升可靠性，分区为区域框架，重大工程前仍需详勘。

结论部分翻译

本研究为拉斯埃尔赫克马海岸带这一埃及极具战略意义且快速城市化的地中海走廊提供了高分辨率多源浅地表表征。通过在87条电气电阻率层析成像（ERT）剖面共址集成钻孔测井与粒度分析，研究超越常规仅钻孔评估，在地质异质区域内交付了空间连续且定量约束的浅地表状况图景。自上而下识别出四个主要电性单元/层：全新统碳酸盐砂（42–1408 Ω·m）从松散沙丘砂至胶结海滩岩；下伏更新统鲕粒灰岩脊（202–557 Ω·m）构成最有利岩土基底；富黏土上新世页岩（2–12 Ω·m；黏土34–36%；厚1.5–3.5 m）构成主要危害，为具差异沉降潜力的膨胀压缩单元；最深裂隙中新世灰岩（816–1287 Ω·m）提供稳固基底并容纳区域含水层。统计关系稳健：电阻率与孔隙度呈反比（r ≈ ?0.85）验证Archie定律；伽马射线与黏土含量近完美相关（r ≈ 0.95）得出场地方程 V_clay= 0.88 × IGR + 0.05，可从测井快速估黏土。空间上浅地表呈NE-SW向展布，中新世基底成椭圆簇，黏土层厚度独立变化暗示古地形控制，故需密集地球物理覆盖，稀疏钻孔不足以可靠表征。引入四级承载力分类（A–D）：A区（ρ > 1000 Ω·m；黏土<2.5%）允浅扩展基础但需ERT筛查岩溶，1000 Ω·m阈值符合埃及建筑规范ECC 203–2020胜任灰岩标准与ASTM D2487非塑性边界，2.5%黏土阈值符合Atterberg非塑性准则；B区需带排水筏基；C区（ρ < 100 Ω·m；黏土+粉砂10–20%）需深桩或地基改良；D区（膨胀上新世黏土）需锚入A区深桩、格构板或石灰-粉煤灰改良（塑性指数降达50%）。水文地质上低阻异常标示盐水入侵，页岩近零自然电位确认为区域隔水层，分隔第四系与中新世系统，对可持续管理与腐蚀风险评估至关重要。该工作流程（高密度ERT、钻孔标定、岩石物理建模、地统计学）可复用于其他干旱海岸大项目（如新阿拉曼、加拉拉城）及类似地中海/红海区域。未来工作含时序ERT探水分迁移、激发极化辨黏土束缚水、机器学习辅助三维反演建预测模型。综上，ERT集成钻孔地球物理与沉积学是复杂海岸负责任开发的实用必要手段；拉斯埃尔赫克马浅地表是风险与机遇并存的空间结构化景观，所述承载力分区、岩石物理关系与空间图为塑造埃及最关键发展前沿的规划者、工程师与水文地质学家提供了科学基础。

要不要我帮你把论文中四个岩土承载力分区（A–D）的具体工程建议与阈值整理成一段简明扼要的要点总结？