被动式含水层管理补给脱盐系统的设计与表征：多层材料优化与技术经济建模

《Results in Engineering》：Design and Characterization of Passive Managed Aquifer Recharge Desalination: Multilayer Material Optimization and Techno-Economic Modelling

　　

淡水短缺，特别是在超干旱沉积盆地中，日益威胁区域水与粮食安全。不断上升的区域温度及竞争性行业需求正迅速耗尽可获取的储量，而后者仅占全球资源的不足1%。浅层非承压含水层储量的快速下降，促使研究人员对替代水源策略产生兴趣，尤其是利用非常规水源进行含水层管理补给（MAR）。全球范围内，不可持续的地下水开采支撑着超过20亿人口及近40%的农业生产。这种水文地质压力在北非西北撒哈拉含水层系统（NWSAS）和乍得盆地等超干旱、跨边界环境中尤为严峻，无管控的水资源开采导致地下水位下降、土壤盐渍化及不可逆的水质退化。为缓解这些地下水供给赤字，与枯竭油气藏相关的深层含水层代表了巨大但未被充分利用的水体体积。然而，利用提取的次储层盐水支持MAR需要高能效的脱盐途径以降低高浓度溶解离子。尽管海水反渗透（SWRO）等常规地表技术提供了可靠的离子截留能力，但其在偏远内陆井口的大规模部署仍受限于巨大的热力学功需求、持续的化学药剂投加以及严重的环境盐水处理瓶颈。通过土壤-含水层处理的被动式地下水抽取、入渗及脱盐提供了一种由静水压力驱动的低能耗、重力驱动型替代方案。然而，此类系统的工程设计受限于对复杂地质构造中溶质运移动力学机制的有限理解。研究人员提出假设：地下脱盐效率主要由增加中间岩性单元内流体停留时间的垂向渗透率差异控制，而工程改良剂通过引入局部化学非平衡汇来放大这些传输控制过程。本研究旨在量化岩性非均质性和工程改良剂如何控制MAR期间的被动式地下脱盐，采用定制的面向对象Python脚本求解器。具体而言，研究人员开发了一个耦合的块中心数值框架，用于模拟分层土壤-含水层系统中的水提取、压力降（ΔP）、平流-地球化学溶质运移及质量平衡衰减阵列。此外，该框架评估了工程反应性改良剂通过耦合二元阳离子交换和微孔扩散截留增强盐度去除的矢量化效能。最后，应用基于方差的全局敏感性分析和蒙特卡洛不确定性量化来确定主导脱盐参数（Cfinal, tsat）的统计变异，以检验层方差贡献和改良剂驱动汇，最终通过净现值（NPV）评估模拟出水对直接回用水分配的适用性。该集成模型旨在将这一地下资产的多种资源回收置于具体背景下，提供一个快速、可扩展的框架来评估其经济可行性及向真实世界环境的转化运行。

论文解读

一、研究背景与现存问题

淡水短缺，尤其在超干旱沉积盆地中，正严重威胁区域水与粮食安全。随着气温上升及各行业用水需求竞争加剧，可获取的淡水资源日益枯竭。在此背景下，深层含水层，特别是与枯竭油气藏相关的深层地层，因其巨大的水体体积而被视为潜在的替代水源。然而，从这些地层中提取的盐水通常具有极高的总溶解固体（TDS）浓度，必须经过高效脱盐处理才能用于含水层管理补给（MAR）。

传统的海水反渗透（SWRO）等表面脱盐技术在偏远内陆井口的应用中面临诸多挑战，包括高昂的热力学能耗、持续的化学药剂投加以及严重的环境盐水处置难题。相比之下，被动式地下脱盐——即通过土壤-含水层处理实现的重力驱动、低能耗脱盐方式——展现出显著优势。然而，当前对该技术的认知存在关键空白：天然地质结构（如不同岩性层的渗透率差异）与工程反应性改良剂如何共同控制高盐流体的水动力学行为、非平衡溶质突破动力学以及空间盐度衰减规律尚不清晰。为此，研究人员提出了核心假设：地下脱盐效率主要由垂向渗透率差异决定，这些差异能增加流体在中低渗透性岩性单元中的停留时间；同时，嵌入的工程改良剂可通过引入局部化学非平衡汇来进一步强化这些传输控制过程。

二、研究方法与技术路线

为验证上述假设，研究人员开发了一套定制的面向对象Python脚本求解器，构建了耦合的块中心数值框架。该框架主要包括以下几项关键技术方法：

1. 1.
   空间域离散与岩石物性参数化：将深层枯竭油气藏简化为五区深度分辨结构连续体，并将这些参数化存储阵列映射至一个分六层的浅层土壤-入渗剖面。各层厚度、孔隙度、流体饱和度及渗透率等关键属性均基于公开文献数据进行赋值。该系统采用显式代数传递循环而非隐式偏微分方程求解器，并通过直接体积平衡约束保证整个模拟期内（20年）的质量守恒。
2. 2.
   平流-弥散-地球化学耦合运移模块：该模块模拟了盐水在分层多孔介质中的物理迁移（包括平流、机械弥散和分子扩散）以及与基质间的化学相互作用，包括二元阳离子交换和可逆吸附等非平衡反应动力学。针对氯离子（Cl?）因阴离子排斥而导致的早期穿透问题，引入了经季铵基团修饰的颗粒活性炭（GAC）作为永久性化学吸附汇，其去除效率通过经验缩放因子η_AC进行量化。
3. 3.
   基于方差的全局敏感性分析与蒙特卡洛不确定性量化：通过对关键输出变量（最终出水浓度C_final和系统饱和寿命t_sat）进行统计分析，该方法能够定量评估不同岩性层参数及改良剂效应对脱盐性能的贡献度，从而识别出最敏感的控制因素。
4. 4.
   技术经济生命周期评估：将模拟得到的出水水质和系统寿命数据输入确定性净现值（NPV）计算模型，以评估该被动脱盐方案在经济上的可行性和实际推广潜力。

三、主要研究结果

3.1 岩性非均质性对水动力和停留时间的影响

研究发现，分层含水层系统中相邻岩性单元间巨大的水力传导率（K）和有效孔隙度差异是控制局部地下水流速矢量的关键因素。低渗透性的粉砂或粘土透镜体会迫使流体在其内部经历更长的停留时间，从而最大化可动大孔隙流体与不动基质块体之间的扩散传质窗口。相反，在高渗透性的砾石或粗砂岩单元中，结构旁通流会显著缩短停留时间，限制地球化学交换，导致快速的溶质突破和对单价离子的弱衰减能力。

3.2 工程反应性改良剂的脱盐效能

通过将斜发沸石、热解生物炭和改性颗粒活性炭（GAC）战略性嵌入中等渗透性地层中，可以显著增强离子去除动力学。沸石凭借其开放的三维铝硅酸盐骨架结构和高负电荷密度，通过竞争性二元阳离子交换机制高效捕获钠离子（Na?）。生物炭则通过其去质子化的含氧表面官能团（如羧基-COOH、羟基-OH）进一步增加净负表面电荷密度和阳离子交换容量（CEC）。对于氯离子（Cl?），传统材料表面的静电排斥效应会导致其过早突破，而经季铵基团改性的GAC能将这种物理滞后转变为高容量的不可逆化学吸附，有效延长了氯离子的截留时间。

3.3 全局敏感性分析揭示主导控制因素

基于方差的全局敏感性分析明确显示，最终出水浓度（C_final）和系统饱和寿命（t_sat）的输出方差主要归因于粘土-粉砂层的参数变化，证实了低渗透性夹层在调控整体脱盐性能中的核心作用。同时，工程改良剂的引入被证明是一个强有力的“扰动”因素，能够显著改变原有的参数敏感性排序，将控制权从纯物理传输转向化学非平衡过程。

四、讨论与结论

研究表明，被动式地下脱盐的效率并非单一因素决定，而是天然地质结构与人工干预措施协同作用的结果。岩性非均质性通过调节流体停留时间为地球化学反应创造了必要条件，而工程改良剂则在此基础上通过引入强大的化学非平衡汇来放大这一效果。这种耦合机制使得即便在处理极高盐度的油田采出水时，也能获得满足灌溉水质标准的出水。

研究结论明确指出：分层地质介质的垂向渗透率对比是控制被动脱盐效率的首要因素，它决定了流体与活性材料的接触时间；而合理设计的工程反应性改良剂层能够有效弥补天然地层在特定离子（尤其是Cl?）去除上的不足。通过集成的技术经济模型评估，该被动脱盐方案在20年的运营周期内展现出了正的净现值（NPV），证明了其在干旱地区进行多资源回收的经济可行性。这一工作为将枯竭油气藏转化为可持续的淡水供应源提供了一个快速、可扩展的计算评估框架，具有重要的理论指导意义和实际应用价值。