页岩气泄漏易受影响的地下水中类采出水化学特征的几率更高

《Nature Sustainability》：Higher odds of produced-water-like chemistry in groundwater vulnerable to shale gas spills

　　

非常规油气开发（Unconventional Oil and Gas Development，UOGD）是指通过水平钻井和水力压裂从页岩及其他低渗透性地层中开采碳氢化合物，该做法提升了美国的国内能源生产和经济水平。然而，这些效益也伴随着公众对饮用水资源、环境和人类健康影响的担忧。本研究探讨了阿巴拉契亚盆地北部观测到的地下水化学特征与UOGD已记录地表泄漏之间的关联。研究人员利用已发表的超过16,000份地下水样本数据，以及高分辨率的、基于物理过程的地下水流和溶质运移模型，发现即使考虑了其他可能解释此类化学印迹的因素，地表泄漏水文脆弱性较高的区域，在家用井水中检测到类似于采出水或深部地层卤水的化学特征的优势比（odds）仍增加三倍。在整个研究区域内，研究人员估计超过120,000名依赖家用井水的人口居住在地下水易受UOGD活动污染的地区。

论文解读：基于物理模型的页岩气泄漏地下水脆弱性评估及其与水质关联研究

非常规油气开发（Unconventional Oil and Gas Development, UOGD）通过水平钻井和水力压裂技术，使美国阿巴拉契亚盆地成为重要的天然气产区，极大推动了国内能源自给与经济发展。然而，伴生的环境风险，尤其是对饮用水资源的潜在威胁，始终是社会关注的焦点。早期研究多集中于甲烷泄漏对地下水的影响，但UOGD过程中涉及的大量淡水、压裂液及返排液、采出水等废水的管理，使更多具有健康风险的化学物质（如疑似致癌物、内分泌干扰物）可能通过地表泄漏释放到环境中。此类泄漏常发生在井场、地面蓄水池及运输管道等环节，被认为是民用井水质恶化的潜在来源之一。然而，污染物来源的归因十分困难，因为常规油气开发、天然卤水上升等也能产生相似的化学印迹。传统化学指纹识别方法虽能区分来源，却只能事后追溯，无法预测未来风险。因此，亟需一种基于物理机制的区域性评估框架，以识别哪些地区在发生泄漏后更易导致地下水污染，从而为预防性管理提供依据。

针对上述问题，研究人员在覆盖宾夕法尼亚、俄亥俄和西弗吉尼亚州、面积达305,000 km2的研究区域内，构建了基于物理过程的地下水脆弱性评估框架，并将其与大规模地下水化学观测数据相结合。研究人员利用超过16,000份民用井水样本和公开的UOGD井场泄漏记录，评估了地下水对地表泄漏的脆弱性（即污染源释放后对饮用水受体的潜在影响程度），并进一步量化了脆弱性与类采出水化学特征检出之间的统计关联。该研究发表于《Nature Sustainability》，其重要意义在于：首次在大区域尺度上，通过水文连通性而非简单距离阈值来识别高风险人群，为制定饮用水保护政策提供了科学依据。

在技术方法上，研究人员主要采用了以下关键手段。首先，基于MODFLOW-6构建了研究区非承压浅层地下水流的稳态数值模型，模型网格为250 m×250 m，采用Zell和Sanford开发的美国本土模型参数化方案，并选取覆盖研究区的四个水文单元域。其次，通过PESTPP-IES迭代集合平滑器进行参数校准，生成包含1,000个水力传导系数场实现的集合，以表征模型不确定性。第三，利用MODPATH-7对每个实现进行前向粒子追踪，模拟从UOGD井场泄漏的保守性溶质在地下水中的稳态运移路径线，并据此计算每个网格单元的脆弱性值（V），其定义为路径线经过该网格的实现比例。第四，结合Shale Network和已发表研究收集的俄亥俄、宾夕法尼亚、西弗吉尼亚州民用井水化学数据，采用Ca/Sr与Ca/Mg离子比值分类准则识别类采出水印迹。最后，建立逻辑回归模型，控制地形位置、距地质构造距离以及常规油气井邻近性等协变量，分析脆弱性与水质指标之间的关联。

研究结果分为三个部分。第一部分“对地下水和污染物路径的地形控制”表明，模拟的地下水位随地表高程变化，地下水总体从高地流向低地并排泄至河流，符合该地区浅层水文地质概念模型。集合校准后的平均绝对误差为0.1～1.0 m，均方根误差为0.9～1.6 m，证明模拟性能良好。第二部分“通过地下水脆弱性绘制的风险热点”中，粒子追踪结果显示，路径线长度中位数为532 m，均值为671 m，95%分位数为1,790 m。仅考虑已记录泄漏的井场时，V≥0.1的面积达590 km2；而将全部UOGD井场作为潜在源时，该面积扩大至3,538 km2。人口风险评估显示，约22,000人生活在对已记录泄漏高脆弱（V≥0.5）的区域；若考虑所有井场，这一数字增至123,000人。热点区域集中在宾夕法尼亚东北部及三州交界地区。第三部分“脆弱性与UOGD类印迹之间的关联”通过逻辑回归发现：脆弱性每增加0.10，检出类采出水化学特征的优势比增加14%；调整协变量后仍为13%。高脆弱样本（V≥0.5）检出类采出水印迹的优势比是低脆弱样本的3倍以上（优势比3.24，P<0.001）。敏感性分析通过随机排列数据验证了结果的稳健性。此外，钡（Ba）浓度超过一级最大污染物水平（MCL，2 mg l^?1）的样本占总体2%，高脆弱样本中为5%；氯离子（Cl）超标也与高脆弱性显著相关。

在讨论部分，研究人员指出，模拟的污染物运移距离远超现行监管缓冲区。宾夕法尼亚、俄亥俄和西弗吉尼亚州的设定距离分别为152 m、15 m和76 m。建议中的762 m和1,000 m缓冲距离分别被34%和21%的模拟路径线超过。与基于固定半径的风险评估相比，该框架排除了缺乏水文连通性的区域，因而估计的受影响人口低于传统方法（1 km与2 km半径分别对应16万和38.9万人），但更为精确。研究同时承认模型假设了稳态流场和保守性溶质运移，忽略了弥散、吸附和降解，且未预测溶质运移时间，因此结果应视为提示性而非确定性。本研究的结论是：基于物理过程的脆弱性评估能够有效识别UOGD地表泄漏对地下水的高风险区域，且脆弱性与民用井水中类采出水化学特征的存在显著相关，支持泄漏影响地下水水质的假说。该框架可作为饮用水保护决策的补充工具，并强调了改进泄漏记录和加强脆弱地区水井监测的必要性。