低渗透性透镜体对饱和多孔介质中DNAPL迁移影响的数值模拟探究

低渗透性透镜体对饱和多孔介质中DNAPL迁移影响的数值模拟探究一、引言1.1研究背景与意义随着工业化进程的加速，大量化学物质被广泛使用，其中重非水相液体（DenseNon-AqueousPhaseLiquids，DNAPL）由于其特殊的物理化学性质，成为了土壤和地下水污染的重要来源之一。常见的DNAPL污染物包括氯代烃化合物、多氯联苯、煤焦油、杂酚油和重油等，这些物质密度大于水，且具有低水溶性、弱迁移性和难降解性等特点。一旦DNAPL进入土壤和地下水系统，会在重力作用下穿过非饱和带进入饱和含水层。由于其低水溶性，它们在地下水中的溶解速度缓慢，但却能持续向周围水体溶出污染物，导致污染范围不断扩大，对地下水环境安全造成极大危害。例如，在一些工业发达地区，由于早期对化工废弃物处理不当，地下水中检测出高浓度的氯代烃类DNAPL污染物，使得当地的地下水资源无法饮用，周边生态环境遭到严重破坏，给居民健康和生态系统平衡带来了长期威胁。而且，由于DNAPL的难降解性，其在土壤及地下水中的污染时间持久，治理和修复难度极高，需要耗费大量的人力、物力和财力。在实际的地质环境中，饱和多孔介质往往具有非均质性，其中低渗透性透镜体广泛存在。这些低渗透性透镜体的存在会显著影响DNAPL在饱和多孔介质中的迁移行为。当DNAPL在迁移过程中遇到低渗透性透镜体时，由于其无法克服该介质的毛管压力，会停止垂向入渗，并在透镜体表面堆积，随后开始横向运移，绕过透镜体后继续向下运动。这种迁移行为的改变会导致DNAPL的污染范围和分布特征发生变化，增加了污染控制和修复的复杂性。研究低渗透性透镜体对饱和多孔介质中DNAPL迁移的影响具有重要的现实意义。在污染控制方面，准确了解DNAPL在非均质介质中的迁移规律，有助于预测污染的扩散范围和趋势，从而提前采取有效的防控措施，避免污染的进一步扩大。在污染修复方面，掌握低渗透性透镜体对DNAPL迁移的影响机制，可以为制定更加科学合理的修复方案提供依据，提高修复效率，降低修复成本。例如，在进行地下水抽提修复时，如果不考虑低渗透性透镜体的影响，可能会导致部分DNAPL无法被有效抽出，从而影响修复效果。而通过深入研究这种影响，我们可以优化抽提井的布局和抽提方案，提高对DNAPL的去除率。1.2国内外研究现状在DNAPL迁移研究方面，国外起步较早。自20世纪80年代起，美国、加拿大等国家的科研团队就开始聚焦于DNAPL在多孔介质中的迁移机制研究。通过室内砂箱实验和数值模拟，揭示了重力、毛细力和黏滞力等多物理因素对DNAPL迁移路径和分布形态的影响。例如，美国地质调查局（USGS）的研究人员利用高精度的CT扫描技术，对砂箱中DNAPL的迁移过程进行实时监测，清晰地观察到了DNAPL在不同孔隙结构中的运移特征。国内相关研究在90年代后期逐渐兴起，随着对地下水污染问题的重视，众多科研机构和高校如南京大学、清华大学等积极开展研究。南京大学的研究团队通过大量实验，分析了不同类型DNAPL在我国典型土壤介质中的迁移特性，发现土壤的颗粒组成和有机质含量对DNAPL的吸附和迁移有着显著影响。但目前对于复杂地质条件下，如存在多种含水层交互、复杂地质构造等情况下DNAPL的迁移规律，研究仍不够深入，缺乏系统的理论和模型。针对低渗透性透镜体对DNAPL迁移影响的研究，国外学者通过构建不同渗透率的多层介质模型，研究发现低渗透性透镜体会改变DNAPL的迁移方向，导致其在透镜体上方堆积并发生侧向扩散。例如，加拿大滑铁卢大学的学者利用数值模拟与物理实验相结合的方法，定量分析了透镜体的厚度、渗透率以及位置等因素对DNAPL迁移的影响程度。国内研究则更多集中在对特定场地的案例分析上，如在一些石油污染场地，通过现场监测和数据分析，发现低渗透性透镜体阻碍了DNAPL的垂直下渗，使得污染呈现分层分布的特征。然而，目前对于低渗透性透镜体与DNAPL之间复杂的相互作用机制，包括多相流条件下的质量交换、化学反应等方面，尚未形成完善的理论体系。数值模拟作为研究DNAPL迁移的重要手段，国外已经开发了多种成熟的数值模型，如HYDRUS、TOUGHREACT等。这些模型能够考虑多相流、溶质运移、化学反应等多种过程，在模拟复杂地质条件下DNAPL的迁移方面具有较高的精度。例如，美国劳伦斯利弗莫尔国家实验室利用TOUGHREACT模型成功模拟了某核污染场地中DNAPL的迁移与扩散，为污染治理提供了科学依据。国内在数值模拟方面，也在不断引进和改进国外模型，并结合国内实际地质条件开发了一些具有自主知识产权的模型，如中国地质科学院开发的适用于我国北方平原地区的地下水污染模拟模型。但整体而言，国内数值模拟技术在模型的通用性、计算效率以及与实际监测数据的耦合应用方面，与国外仍存在一定差距。1.3研究内容与方法本研究主要聚焦于低渗透性透镜体对饱和多孔介质中DNAPL迁移的影响，具体研究内容如下：探究低渗透性透镜体对DNAPL迁移路径的影响：通过数值模拟，构建含有不同特征低渗透性透镜体的饱和多孔介质模型，研究在重力、毛细力等作用下，DNAPL在遇到低渗透性透镜体时迁移路径的改变。分析透镜体的位置、厚度、渗透率等参数对DNAPL迁移路径的具体影响规律，明确DNAPL在透镜体表面的堆积、侧向运移以及绕过透镜体后的继续下渗等行为特征。分析低渗透性透镜体对DNAPL饱和度分布的影响：利用数值模拟和实验相结合的方法，研究低渗透性透镜体存在时，DNAPL在饱和多孔介质中的饱和度分布情况。通过改变透镜体的相关参数，观察DNAPL饱和度在不同区域的变化，分析饱和度分布与透镜体特征之间的关系，探究透镜体对DNAPL在多孔介质中分布均匀性的影响机制。研究低渗透性透镜体对DNAPL浓度分布的影响：基于数值模拟结果，分析低渗透性透镜体对DNAPL在饱和多孔介质中浓度分布的影响。考虑DNAPL的溶解、扩散等过程，研究透镜体周围及下游区域DNAPL浓度的变化规律，明确透镜体对污染羽扩散范围和浓度衰减的影响，为评估污染风险提供依据。模型验证与参数敏感性分析：收集实际场地数据或开展室内实验，对所建立的数值模型进行验证，确保模型能够准确反映低渗透性透镜体存在时DNAPL的迁移行为。同时，对模型中的关键参数，如透镜体渗透率、孔隙度、DNAPL的物理性质参数等进行敏感性分析，确定各参数对模拟结果的影响程度，为模型的优化和实际应用提供参考。在研究方法上，本研究综合运用了多种手段：数值模拟方法：采用专业的多相流数值模拟软件，如COMSOLMultiphysics、HYDRUS等，建立饱和多孔介质中DNAPL迁移的数值模型。利用这些软件强大的计算功能和丰富的物理模块，能够准确模拟多相流、传质等复杂过程，通过设置不同的模型参数，系统地研究低渗透性透镜体对DNAPL迁移的影响。例如，在COMSOLMultiphysics中，可以利用其多孔介质模块，精确地定义饱和多孔介质的物理性质和边界条件，通过求解Navier-Stokes方程和对流-扩散方程，模拟DNAPL在多孔介质中的流动和扩散行为。实验研究方法：设计并开展室内砂箱实验，在砂箱中填充不同类型的多孔介质，并设置具有特定参数的低渗透性透镜体。通过向砂箱中注入DNAPL，利用可视化技术（如CT扫描、核磁共振成像等）实时监测DNAPL的迁移过程，获取DNAPL的迁移路径、饱和度分布和浓度分布等数据。这些实验数据不仅可以用于验证数值模拟结果，还能为深入理解DNAPL的迁移机制提供直观的依据。比如，利用CT扫描技术，可以清晰地观察到DNAPL在砂箱中的三维分布情况，以及其与低渗透性透镜体相互作用的过程。理论分析方法：基于渗流理论、多相流理论和传质理论，对低渗透性透镜体影响下DNAPL的迁移行为进行理论分析。推导相关的数学模型和公式，解释实验和数值模拟结果，揭示DNAPL迁移过程中的物理机制和规律。例如，运用达西定律和毛细压力理论，分析DNAPL在不同渗透率介质中的渗流速度和在透镜体表面的堆积条件，为研究提供理论基础。二、相关理论基础2.1DNAPL的特性与迁移理论2.1.1DNAPL的基本特性重非水相液体（DNAPL）具有一系列独特的物理化学特性，这些特性对其在饱和多孔介质中的迁移行为起着决定性作用。从密度方面来看，DNAPL的密度显著大于水，通常大于1.01g/cm^3。例如，常见的氯代烃类DNAPL，如四氯化碳，其密度约为1.595g/cm^3，三氯乙烯的密度约为1.462g/cm^3。这种高密度特性使得DNAPL在进入土壤和地下水系统后，在重力作用下会迅速下沉，穿过非饱和带进入饱和含水层，成为地下水污染的重要潜在源。在溶解度上，DNAPL表现出极低的水溶性，一般溶解度（水）小于20g/L。以多氯联苯为例，其在水中的溶解度极低，这使得它们在地下水中难以被水稀释和分散。这种低溶解性导致DNAPL在地下水中溶解速度缓慢，但却能持续向周围水体溶出污染物，形成长期的污染源，不断扩大污染范围，对地下水环境造成长期的危害。DNAPL的黏度相对较低，通常低于1cP，这一特性使得其流动性较强，有助于其向地表以下移动。较低的黏度使得DNAPL能够更容易地在土壤孔隙和含水层空隙中渗透，增加了其在地下环境中的迁移能力。此外，DNAPL还具有较低的界面张力，这使其容易入渗进入极小孔隙或贯穿黏土质。例如，在一些黏土含量较高的土壤中，DNAPL仍能够凭借其低界面张力特性，穿过细小的黏土孔隙，进一步扩大污染范围。同时，多数DNAPL还具有化学性质相对稳定、生物降解性低的特点，这使得它们在土壤及地下水中污染时间持久，难以通过自然降解的方式去除，增加了污染修复的难度。2.1.2DNAPL在饱和多孔介质中的迁移机制DNAPL在饱和多孔介质中的迁移是一个复杂的过程，受到多种力的综合作用，其中重力、毛细力和黏滞力是影响其迁移的主要因素。重力是DNAPL迁移的主要驱动力之一。由于DNAPL密度大于水，在重力作用下，DNAPL会在饱和多孔介质中向下迁移。在均质的饱和多孔介质中，DNAPL会沿着孔隙通道逐渐向下渗透，其迁移速度与DNAPL的密度、重力加速度以及多孔介质的渗透率等因素有关。根据达西定律，DNAPL在多孔介质中的渗流速度v可以表示为v=-\frac{k}{\mu}(\rho_{DNAPL}g-\rho_wg)，其中k为多孔介质的渗透率，\mu为DNAPL的黏度，\rho_{DNAPL}为DNAPL的密度，\rho_w为水的密度，g为重力加速度。从公式中可以看出，DNAPL的密度越大，重力作用越强，其在多孔介质中的渗流速度也就越快。毛细力对DNAPL的迁移也有着重要影响。在饱和多孔介质中，孔隙大小不一，当DNAPL进入孔隙时，会受到毛细力的作用。毛细力的大小与孔隙半径、液体表面张力以及接触角等因素有关。根据拉普拉斯公式，毛细压力P_c可以表示为P_c=\frac{2\sigma\cos\theta}{r}，其中\sigma为液体表面张力，\theta为接触角，r为孔隙半径。当DNAPL遇到较小的孔隙时，毛细力可能会阻止其进一步迁移，使其在孔隙中滞留。相反，当孔隙较大时，毛细力相对较小，DNAPL更容易通过。在一些孔隙结构复杂的饱和多孔介质中，DNAPL可能会在毛细力的作用下，在孔隙中形成不连续的液滴或薄膜，影响其迁移路径和分布。黏滞力则对DNAPL的迁移起到阻碍作用。DNAPL在多孔介质中流动时，会与孔隙壁面和周围的水发生摩擦，产生黏滞力。黏滞力的大小与DNAPL的黏度、流速以及多孔介质的孔隙结构等因素有关。较高的黏度会导致DNAPL受到更大的黏滞阻力，从而减缓其迁移速度。而且，多孔介质的孔隙结构越复杂，DNAPL与孔隙壁面的摩擦面积越大，黏滞力也就越大。例如，在孔隙狭窄且曲折的多孔介质中，DNAPL的迁移会受到较大的黏滞阻力，其迁移速度会明显降低。在实际的饱和多孔介质中，DNAPL的迁移是重力、毛细力和黏滞力共同作用的结果。当DNAPL进入饱和多孔介质后，首先会在重力作用下向下迁移。在迁移过程中，遇到孔隙大小变化时，毛细力会对其迁移产生影响，可能导致DNAPL在某些孔隙中滞留或改变迁移方向。同时，黏滞力会始终阻碍DNAPL的迁移，使其迁移速度逐渐降低。而且，DNAPL的迁移还会受到地下水流速、多孔介质的非均质性等因素的影响，进一步增加了其迁移机制的复杂性。2.2饱和多孔介质的特性与参数2.2.1饱和多孔介质的基本性质饱和多孔介质是指内部含有大量相互连通孔隙且孔隙完全被水充满的介质，广泛存在于自然界的土壤、岩石等地质体中。其基本性质对DNAPL的迁移行为有着至关重要的影响。孔隙度是饱和多孔介质的一个重要参数，它表示孔隙体积与多孔介质总体积的比值，反映了介质中孔隙空间的大小。孔隙度的大小直接影响着DNAPL在多孔介质中的储存和迁移能力。一般来说，孔隙度越大，多孔介质中可供DNAPL储存和流动的空间就越大，DNAPL的迁移速度也就相对较快。例如，在粗砂组成的饱和多孔介质中，孔隙度通常较大，DNAPL能够更容易地在其中渗透和扩散；而在黏土等细颗粒组成的饱和多孔介质中，孔隙度较小，DNAPL的迁移会受到较大的阻碍。常见的饱和多孔介质孔隙度范围如下表所示：介质类型孔隙度范围粗砂0.35-0.45中砂0.30-0.40细砂0.25-0.35粉土0.35-0.50黏土0.40-0.70渗透率是衡量饱和多孔介质允许流体通过能力的重要指标，它与孔隙的大小、形状、连通性等因素密切相关。渗透率越大，流体在多孔介质中流动时所受到的阻力就越小，DNAPL的迁移也就越容易。根据达西定律，流体在多孔介质中的渗流速度与渗透率成正比。在实际的地质环境中，不同类型的饱和多孔介质渗透率差异较大。例如，砂岩的渗透率一般在10^{-15}-10^{-12}m^2之间，而页岩的渗透率则低得多，通常在10^{-21}-10^{-18}m^2之间。这种渗透率的差异会导致DNAPL在不同介质中的迁移速度和路径有很大不同。当DNAPL从渗透率较高的介质进入渗透率较低的介质时，其迁移速度会显著降低，甚至可能在介质界面处发生聚集和滞留。孔隙结构也是影响饱和多孔介质中DNAPL迁移的关键因素。孔隙结构包括孔隙的大小分布、形状、连通性以及孔隙表面的性质等。复杂的孔隙结构会使DNAPL的迁移路径变得曲折，增加其与孔隙壁面的摩擦和碰撞，从而影响其迁移速度和分布。例如，在孔隙大小分布不均匀的饱和多孔介质中，DNAPL可能更容易在较大的孔隙中流动，而在较小的孔隙中则容易被滞留。而且，孔隙的连通性越好，DNAPL就越容易在多孔介质中扩散；反之，若孔隙连通性较差，DNAPL的迁移会受到很大限制。此外，孔隙表面的亲水性或疏水性也会影响DNAPL与孔隙壁面的相互作用，进而影响其迁移行为。亲水性孔隙表面对水有较强的吸附力，可能会阻碍DNAPL的迁移；而疏水性孔隙表面则相对有利于DNAPL的流动。2.2.2低渗透性透镜体的特征与分布低渗透性透镜体是指在饱和多孔介质中，具有相对较低渗透率的局部区域，其形状类似于透镜，故而得名。这些透镜体在实际的地质环境中广泛存在，对饱和多孔介质中DNAPL的迁移行为有着显著的影响。低渗透性透镜体的形状通常较为复杂，一般呈椭圆形或不规则的扁平状。其长轴和短轴的长度比例会因地质条件的不同而有所变化。在一些沉积环境较为稳定的区域，低渗透性透镜体的形状可能相对规则，长轴与短轴的比例较为接近；而在地质构造活动较为频繁的地区，透镜体可能会受到挤压、拉伸等作用，导致形状变得更加不规则，长轴与短轴的比例差异较大。通过对实际场地的地质勘探数据进行分析，发现低渗透性透镜体的长轴长度范围可以从数厘米到数十米不等，短轴长度则一般在数毫米到数厘米之间。例如，在某一砂质含水层中，通过钻孔取芯和地质雷达探测发现，低渗透性透镜体的长轴长度多在1-5米之间，短轴长度在0.1-0.5米之间，其形状近似椭圆形。低渗透性透镜体的大小也具有较大的变化范围。其面积和厚度是描述大小的两个重要参数。透镜体的面积可以从几平方厘米到数百平方米不等，厚度则一般在几毫米到数米之间。在不同的地质条件下，低渗透性透镜体的大小分布存在差异。在浅层的冲积层中，由于沉积作用相对较弱，低渗透性透镜体的面积和厚度可能相对较小；而在深层的沉积岩中，由于长期的地质作用，透镜体的面积和厚度可能会较大。根据相关研究和实际案例分析，在一些河流冲积平原的含水层中，低渗透性透镜体的面积多在1-10平方米之间，厚度在0.2-1米之间；而在一些深部的砂岩含水层中，透镜体的面积可达数十平方米，厚度也能达到数米。低渗透性透镜体在饱和多孔介质中的空间分布具有随机性和非均质性。它们可能孤立存在，也可能相互连接形成复杂的网络结构。透镜体的分布位置和方向与地质构造、沉积环境等因素密切相关。在沉积过程中，由于水流速度、沉积物来源等因素的变化，低渗透性透镜体会在不同的位置和深度形成。例如，在河流沉积环境中，当水流速度突然减慢时，细颗粒的沉积物会在局部区域堆积，形成低渗透性透镜体。这些透镜体的分布方向可能与水流方向平行或呈一定角度。而且，不同区域的饱和多孔介质中，低渗透性透镜体的分布密度也有所不同。在一些地质条件复杂的区域，透镜体的分布密度可能较高，导致饱和多孔介质的非均质性更加明显；而在地质条件相对简单的区域，透镜体的分布密度则较低。通过对多个实际场地的调查发现，低渗透性透镜体在饱和多孔介质中的分布密度范围可以从每立方米几个到几十个不等。例如，在某一受污染的工业场地，通过高密度电阻率法和地震波反射法探测发现，在特定的含水层区域内，低渗透性透镜体的分布密度约为每立方米5-10个，这些透镜体的存在显著影响了地下水中污染物的迁移路径和分布特征。2.3数值模拟的基本原理与方法2.3.1数值模拟的基本方程在研究低渗透性透镜体对饱和多孔介质中DNAPL迁移影响的数值模拟中，涉及到多个基本方程，这些方程是描述物理过程的数学基础，其中连续性方程、动量方程和能量方程是最为关键的控制方程。连续性方程，也被称为质量守恒方程，它体现了物质在迁移过程中质量不会凭空产生或消失的基本原理。对于饱和多孔介质中的DNAPL迁移问题，连续性方程可以用来描述DNAPL和水的质量守恒情况。以三维空间为例，假设\rho_{DNAPL}和\rho_w分别为DNAPL和水的密度，u_{DNAPL}和u_w分别为DNAPL和水在x、y、z方向上的流速分量，则DNAPL的连续性方程可表示为：\frac{\partial(\phi\rho_{DNAPL}S_{DNAPL})}{\partialt}+\frac{\partial(\rho_{DNAPL}S_{DNAPL}u_{DNAPL,x})}{\partialx}+\frac{\partial(\rho_{DNAPL}S_{DNAPL}u_{DNAPL,y})}{\partialy}+\frac{\partial(\rho_{DNAPL}S_{DNAPL}u_{DNAPL,z})}{\partialz}=0其中\phi为多孔介质的孔隙度，S_{DNAPL}为DNAPL的饱和度，t为时间。该方程的左边第一项表示单位时间内单位体积多孔介质中DNAPL质量的变化率，后面三项分别表示DNAPL在x、y、z方向上的质量通量的散度。同样，水的连续性方程也可以按照类似的形式写出。通过连续性方程，我们能够追踪DNAPL和水在饱和多孔介质中的质量分布随时间的变化情况，为研究DNAPL的迁移提供了质量守恒的约束条件。动量方程，即Navier-Stokes方程，它描述了流体在运动过程中的动量变化与所受外力之间的关系。在饱和多孔介质中，DNAPL和水的流动同样受到动量方程的支配。对于不可压缩流体，动量方程可以简化为：\rho\frac{\partial\vec{u}}{\partialt}+\rho(\vec{u}\cdot\nabla)\vec{u}=-\nablap+\mu\nabla^2\vec{u}+\vec{F}其中\rho为流体的密度（对于DNAPL或水分别取相应的密度值），\vec{u}为流体的速度矢量，p为压力，\mu为动力黏度，\vec{F}为作用在流体上的外力（如重力等）。在饱和多孔介质中，由于多孔介质的存在，流体的流动会受到额外的阻力，通常用达西定律来描述这种阻力对流速的影响。将达西定律与动量方程相结合，可以更准确地描述DNAPL和水在饱和多孔介质中的流动行为。例如，对于DNAPL，其在多孔介质中的流速\vec{u}_{DNAPL}可以通过达西定律表示为\vec{u}_{DNAPL}=-\frac{k_{rDNAPL}k}{\mu_{DNAPL}}(\nablap_{DNAPL}-\rho_{DNAPL}g\nablah)，其中k_{rDNAPL}为DNAPL的相对渗透率，k为多孔介质的绝对渗透率，p_{DNAPL}为DNAPL的压力，g为重力加速度，h为高度。动量方程为我们研究DNAPL在饱和多孔介质中的迁移速度和方向提供了理论依据，通过求解动量方程，可以得到不同时刻DNAPL在多孔介质中的流速分布，进而分析其迁移路径的变化。能量方程主要用于描述系统内能量的传递和转化过程。在DNAPL迁移的研究中，虽然能量方程不像连续性方程和动量方程那样直接描述迁移行为，但在考虑一些与温度相关的因素时，能量方程就显得尤为重要。例如，当DNAPL在迁移过程中发生化学反应或与周围介质进行热交换时，会导致温度的变化，而温度的变化又会反过来影响DNAPL的物理性质（如黏度、密度等），进而影响其迁移行为。能量方程的一般形式为：\rhoc_p\frac{\partialT}{\partialt}+\rhoc_p(\vec{u}\cdot\nabla)T=\nabla\cdot(k_T\nablaT)+Q其中\rho为流体的密度，c_p为定压比热容，T为温度，k_T为热导率，Q为热源项。在实际的数值模拟中，如果考虑到温度对DNAPL迁移的影响，就需要同时求解能量方程，以准确描述整个物理过程。通过能量方程，可以计算出饱和多孔介质中温度的分布情况，进而分析温度变化对DNAPL迁移的影响机制，如温度升高可能导致DNAPL黏度降低，从而使其迁移速度加快等。2.3.2常用的数值模拟方法在低渗透性透镜体对饱和多孔介质中DNAPL迁移影响的数值模拟研究中，有限差分法、有限元法和有限体积法是常用的数值方法，它们各自基于不同的原理，在模拟复杂物理过程中发挥着重要作用。有限差分法是一种将连续问题离散化的数值方法，其基本原理是将连续域上的偏微分方程在空间和时间上离散化。首先，将连续的求解域划分成有限个网格点，这些网格点在空间上形成一个规则或不规则的网格系统。然后，在每个网格点上对描述DNAPL迁移的偏微分方程进行泰勒展开，通过保留一阶和二阶导数的差分近似，将偏微分方程转化为关于网格点上函数值的差分方程。例如，对于一维的对流-扩散方程\frac{\partialC}{\partialt}+u\frac{\partialC}{\partialx}=D\frac{\partial^2C}{\partialx^2}（其中C为DNAPL的浓度，u为流速，D为扩散系数），在空间上采用中心差分格式对导数进行近似，\frac{\partialC}{\partialx}\approx\frac{C_{i+1}-C_{i-1}}{2\Deltax}，\frac{\partial^2C}{\partialx^2}\approx\frac{C_{i+1}-2C_i+C_{i-1}}{\Deltax^2}，在时间上采用向前差分格式\frac{\partialC}{\partialt}\approx\frac{C_i^{n+1}-C_i^n}{\Deltat}，其中i表示空间网格点的编号，n表示时间步的编号，\Deltax和\Deltat分别为空间步长和时间步长。将这些差分近似代入原方程，就得到了关于网格点上浓度值的差分方程，通过求解这个差分方程组，就可以得到不同时刻各个网格点上DNAPL浓度的近似解。有限差分法的优点是简单直观，易于理解和编程实现，对于一些简单的几何形状和规则的网格系统，能够快速得到较高精度的数值解。而且，该方法在处理双曲型和抛物型偏微分方程时表现出色，能够有效地模拟DNAPL在饱和多孔介质中的迁移过程。然而，有限差分法的网格划分对解的精度和稳定性有较大影响，如果网格划分不合理，可能会导致数值振荡或误差积累。而且，在处理复杂边界条件时，有限差分法相对困难，需要采用一些特殊的处理技巧来保证计算的准确性。有限元法以变分原理和加权余量法为基础，其核心思想是将计算域划分为非重叠的、互相连接的有限个单元，这些单元可以是三角形、四边形、四面体等不同形状，以适应复杂的几何形状。在每个单元内部选择特定的节点作为插值点，通过线性组合的方式近似求解函数。对于DNAPL迁移问题，将描述其迁移的偏微分方程转化为相应的变分形式，然后在每个单元上进行离散化处理。例如，对于二维的渗流问题，将计算区域划分为多个三角形单元，在每个三角形单元内，假设水头函数h可以表示为节点水头值h_i、h_j、h_k的线性插值函数h(x,y)=\sum_{l=i,j,k}N_l(x,y)h_l，其中N_l(x,y)为形函数，它是关于坐标(x,y)的函数，决定了节点值对单元内任意点函数值的贡献。将这种插值函数代入渗流方程的变分形式中，通过加权余量法得到关于节点水头值的代数方程组，求解这个方程组就可以得到整个计算区域内的水头分布，进而根据达西定律计算出DNAPL的流速和迁移路径。有限元法的最大优势在于能够很好地适应复杂的几何形状和边界条件，对于含有低渗透性透镜体的饱和多孔介质，无论其形状和分布如何复杂，都可以通过合理划分单元来准确模拟。而且，通过调整权函数和插值函数的形式，可以生成不同精度和适应性的有限元模型。然而，有限元法的计算过程相对复杂，需要较多的计算资源和时间，特别是在处理大规模问题时，计算量会显著增加。有限体积法是将连续域划分为有限个体积单元，并在每个体积单元上对控制方程进行积分的数值方法。其基本步骤是首先将计算区域划分成一系列不重叠的控制体积，这些控制体积围绕着网格节点分布。然后，在每个控制体积上对描述DNAPL迁移的偏微分方程进行积分，得到关于未知函数（如DNAPL的浓度、流速等）在控制体积界面上的积分方程。以一维的质量守恒方程\frac{\partial(\rhou)}{\partialx}=0为例，在控制体积V上积分可得\int_{V}\frac{\partial(\rhou)}{\partialx}dV=0，根据高斯公式，将体积分转化为面积分\oint_{S}\rhou\cdot\vec{n}dS=0，其中S为控制体积的表面，\vec{n}为表面的单位外法向量。通过对控制体积界面上的通量进行近似计算，将积分方程转化为关于节点未知量的代数方程。例如，在界面上采用中心差分或迎风格式来近似通量，将界面上的未知函数值用相邻节点的值来表示。最后，联立所有控制体积的代数方程，求解得到整个计算区域内的未知函数分布。有限体积法的优点是物理意义明确，它基于守恒原理，保证了在每个控制体积上物理量的守恒。而且，该方法适用于复杂边界条件和流动场，在处理复杂几何形状时也较为方便，能够灵活地处理含有低渗透性透镜体的饱和多孔介质模型。然而，有限体积法在处理边界条件时需要一定的技巧，特别是在处理非结构化网格时，计算量会相对较大，并且在选择合适的通量近似格式时需要谨慎考虑，以确保计算的稳定性和精度。三、低渗透性透镜体对DNAPL迁移影响的数值模型建立3.1模型假设与简化为了构建能够有效模拟低渗透性透镜体对饱和多孔介质中DNAPL迁移影响的数值模型，需对实际复杂的物理过程进行一系列合理的假设与简化处理。在介质方面，将饱和多孔介质视为连续介质，忽略其微观孔隙结构的细微差异。尽管实际的饱和多孔介质由大小、形状各异的颗粒组成，孔隙结构极为复杂，但在宏观尺度的数值模拟中，将其看作连续介质能够简化数学描述，便于进行数值计算。同时，假设多孔介质各向同性，即其在各个方向上的物理性质（如渗透率、孔隙度等）相同。然而，实际的地质介质往往存在一定程度的各向异性，例如在沉积岩中，由于沉积作用的方向性，不同方向上的渗透率可能存在差异。但在本模型中，为了突出低渗透性透镜体对DNAPL迁移的主要影响，暂不考虑这种各向异性，使模型更具一般性和可分析性。对于低渗透性透镜体，假设其形状为规则的椭圆形，且边界清晰。实际的低渗透性透镜体形状复杂多样，边界也较为模糊，但通过将其简化为规则形状和清晰边界，能够更方便地确定其几何参数，从而在模型中准确地定义其位置和大小。而且，在模型中认为低渗透性透镜体与周围饱和多孔介质之间的接触是理想的，不存在过渡区域或界面阻力。对于流体，假设DNAPL和水均为不可压缩流体。在实际的迁移过程中，虽然DNAPL和水在一定程度上会受到压力和温度变化的影响而发生体积变化，但在本研究关注的范围内，这种变化相对较小，对迁移结果的影响可忽略不计。因此，将它们视为不可压缩流体，能够简化连续性方程和动量方程的求解过程。同时，不考虑流体的黏性变化，认为DNAPL和水的黏度在迁移过程中保持恒定。尽管实际流体的黏度会受到温度、压力以及与多孔介质相互作用等因素的影响，但在本模型中，为了简化计算，暂不考虑这些复杂因素对黏度的影响。而且，忽略流体之间的化学反应，仅考虑DNAPL和水的物理迁移过程。在实际的地下环境中，DNAPL和水之间可能会发生一些化学反应，如氧化还原反应、水解反应等，但这些反应过程相对复杂，对其进行精确模拟需要大量的实验数据和复杂的化学反应动力学模型。在本研究中，为了集中研究低渗透性透镜体对DNAPL迁移的物理影响，暂不考虑这些化学反应。在迁移过程方面，假设DNAPL的迁移仅受重力、毛细力和黏滞力的作用，忽略其他可能的影响因素，如分子扩散、生物降解等。分子扩散在DNAPL迁移过程中会使DNAPL在浓度梯度的作用下发生扩散，但在本模型中，由于主要关注的是宏观尺度上的迁移行为，分子扩散的影响相对较小，可忽略不计。生物降解则涉及到微生物对DNAPL的分解作用，这是一个复杂的生物化学过程，受到多种因素的影响。在本研究中，为了简化模型，暂不考虑生物降解对DNAPL迁移的影响。同时，假设多孔介质的孔隙结构在DNAPL迁移过程中不发生变化。实际的DNAPL迁移过程可能会对多孔介质的孔隙结构产生一定的影响，如DNAPL的吸附、堵塞等作用可能会改变孔隙的大小和连通性。但在本模型中，为了突出低渗透性透镜体对DNAPL迁移的影响，暂不考虑这种孔隙结构的变化。而且，在模型中忽略了边界条件对迁移过程的影响，假设边界条件为理想的稳定状态，如定水头边界、无通量边界等。在实际的地质环境中，边界条件可能会受到地下水补给、排泄等因素的影响而发生变化，但在本研究中，为了简化模型，对边界条件进行了理想化处理。通过这些假设与简化，能够构建出相对简洁且有效的数值模型，以便深入研究低渗透性透镜体对饱和多孔介质中DNAPL迁移的影响。3.2模型参数的确定3.2.1饱和多孔介质参数饱和多孔介质的参数对于准确模拟DNAPL的迁移行为至关重要。其中，孔隙度是衡量多孔介质中孔隙空间占比的关键指标，它反映了介质储存流体的能力。在本研究中，通过实验测定与理论计算相结合的方式确定孔隙度。对于砂质饱和多孔介质，利用排水法进行实验测定。将已知体积V_{total}的砂样放入容器中，使其充分饱和水，然后缓慢排出孔隙中的水，测量排出水的体积V_{void}，则孔隙度\phi=\frac{V_{void}}{V_{total}}。通过多次重复实验，得到该砂质饱和多孔介质的孔隙度平均值约为0.38。同时，参考相关文献中对类似砂质介质的研究结果，进一步验证了该测定值的合理性。例如，文献[具体文献]中对相同类型砂质饱和多孔介质的孔隙度测定结果在0.35-0.40之间，与本研究的测定值相符。渗透率是描述饱和多孔介质允许流体通过能力的重要参数，其大小与孔隙的大小、形状、连通性等密切相关。在确定渗透率时，采用稳态法进行实验测量。在一定的压力差\DeltaP下，让已知黏度\mu的流体通过横截面积为A、长度为L的饱和多孔介质样品，测量单位时间内通过的流体体积流量Q，根据达西定律Q=\frac{kA\DeltaP}{\muL}，可计算出渗透率k。经过实验测定，本研究中饱和多孔介质的渗透率约为10^{-12}m^2。此外，考虑到多孔介质的微观孔隙结构具有分形特征，通过理论计算得到渗透率与分形维数等参数之间的数学模型，进一步验证实验测定的渗透率值。根据分形理论，渗透率k与孔隙度\phi、分形维数D等参数的关系可表示为k=C\phi^mD^n（其中C、m、n为常数），通过对多孔介质微观结构的分析，确定相关常数的值，计算得到的渗透率理论值与实验测定值在合理的误差范围内。孔隙大小分布对DNAPL在饱和多孔介质中的迁移路径和饱和度分布有着显著影响。采用压汞仪对饱和多孔介质的孔隙大小分布进行测量。压汞仪通过施加不同的压力，将汞压入多孔介质的孔隙中，根据汞的注入量与压力的关系，得到孔隙大小分布曲线。测量结果表明，该饱和多孔介质的孔隙大小主要集中在0.1-10\mum之间，其中孔径为1-5\mum的孔隙占比较大，约为60%。这一孔隙大小分布特征决定了DNAPL在迁移过程中，较小孔径的孔隙会对其产生较大的毛细阻力，影响其迁移速度和路径。同时，结合理论分析，利用毛细管压力理论，分析不同孔径孔隙对DNAPL迁移的影响机制，进一步明确孔隙大小分布在DNAPL迁移过程中的作用。3.2.2DNAPL参数DNAPL的参数直接影响其在饱和多孔介质中的迁移特性，准确确定这些参数是数值模拟的关键环节。密度是DNAPL的重要物理性质之一，它决定了DNAPL在重力作用下的迁移方向和速度。在本研究中，针对常见的DNAPL污染物三氯乙烯（TCE），采用比重瓶法测定其密度。将比重瓶洗净、烘干后，精确称重m_1，然后将TCE缓慢注入比重瓶中至刻度线，再次称重m_2，同时测量此时TCE的温度T。根据公式\rho_{DNAPL}=\frac{m_2-m_1}{V}（其中V为比重瓶的容积），计算得到TCE在该温度下的密度。经过多次测量取平均值，得到三氯乙烯在20^{\circ}C时的密度约为1.462g/cm^3。这一密度值与相关文献报道的数值相符，如文献[具体文献]中对三氯乙烯密度的测量结果在1.460-1.465g/cm^3之间，验证了测量的准确性。黏度是影响DNAPL迁移的另一个关键参数，它反映了流体内部的摩擦阻力。对于三氯乙烯的黏度，采用旋转黏度计进行测量。将三氯乙烯样品放入旋转黏度计的测量杯中，设置不同的转速n，测量相应的扭矩M，根据旋转黏度计的工作原理和仪器常数K，通过公式\mu=\frac{M}{Kn}计算出黏度。在20^{\circ}C条件下，经过多次测量，得到三氯乙烯的黏度约为0.57cP。同时，考虑到温度对黏度的影响，参考相关文献中给出的黏度与温度的关系模型，对不同温度下三氯乙烯的黏度进行预测和修正。例如，根据Andrade方程\mu=Ae^{\frac{B}{T}}（其中A、B为常数，T为绝对温度），结合已知的三氯乙烯在不同温度下的黏度数据，确定方程中的常数A和B，从而能够准确预测不同温度条件下三氯乙烯的黏度，为数值模拟提供更准确的参数。溶解度是衡量DNAPL在水中溶解能力的指标，它对DNAPL在饱和多孔介质中的浓度分布和污染扩散有着重要影响。采用平衡法测定三氯乙烯在水中的溶解度。将过量的三氯乙烯与一定量的水放入密闭容器中，在恒温条件下振荡一定时间，使两者达到溶解平衡。然后，通过离心分离或过滤等方法，将未溶解的三氯乙烯与溶液分离，采用气相色谱-质谱联用仪（GC-MS）分析溶液中三氯乙烯的浓度，即为其在该温度下的溶解度。在20^{\circ}C时，测得三氯乙烯在水中的溶解度约为1100mg/L。这一溶解度数据与相关研究结果相近，如文献[具体文献]中报道三氯乙烯在20^{\circ}C时的溶解度在1000-1200mg/L之间，表明测量结果可靠。扩散系数是描述DNAPL在饱和多孔介质中分子扩散能力的参数，它对DNAPL的迁移和浓度分布也有重要影响。在本研究中，采用稳态扩散池法测定三氯乙烯在水中的扩散系数。将饱和多孔介质填充在扩散池的中间隔板两侧，一侧注入含有三氯乙烯的溶液，另一侧注入纯水，在恒温条件下，使三氯乙烯在浓度梯度的作用下通过饱和多孔介质向纯水中扩散。通过监测纯水中三氯乙烯浓度随时间的变化，根据Fick第二定律\frac{\partialC}{\partialt}=D\frac{\partial^2C}{\partialx^2}（其中C为浓度，t为时间，D为扩散系数，x为扩散距离），结合实验数据拟合得到三氯乙烯在饱和多孔介质中的扩散系数约为1.0\times10^{-9}m^2/s。同时，考虑到多孔介质的孔隙结构和曲折度对扩散系数的影响，引入曲折度因子\tau，对扩散系数进行修正，修正后的扩散系数D_{eff}=D/\tau，以更准确地反映三氯乙烯在实际饱和多孔介质中的扩散行为。3.2.3低渗透性透镜体参数低渗透性透镜体的参数是影响DNAPL迁移的关键因素，准确获取这些参数对于深入研究其对DNAPL迁移的影响至关重要。渗透率是低渗透性透镜体的核心参数之一，它决定了透镜体对DNAPL迁移的阻碍程度。在本研究中，通过实验测量和数值反演相结合的方法确定低渗透性透镜体的渗透率。对于实验室制备的低渗透性透镜体样品，采用稳态法进行渗透率测量。在一定的压力差下，让已知黏度的流体通过低渗透性透镜体样品，测量单位时间内通过的流体体积流量，根据达西定律计算出渗透率。然而，在实际地质环境中，低渗透性透镜体的渗透率难以直接测量，因此采用数值反演的方法。根据实际场地的水文地质条件和监测数据，建立包含低渗透性透镜体的数值模型，通过调整模型中透镜体的渗透率参数，使模拟结果与实际监测数据（如地下水位、DNAPL浓度分布等）达到最佳匹配，从而反演出低渗透性透镜体的渗透率。经过多次反演计算，得到本研究中低渗透性透镜体的渗透率约为10^{-15}m^2，远低于周围饱和多孔介质的渗透率，这表明低渗透性透镜体对DNAPL的迁移具有显著的阻碍作用。同时，参考相关文献中对类似地质条件下低渗透性透镜体渗透率的研究结果，验证了反演结果的合理性。例如，文献[具体文献]中对某一砂质含水层中低渗透性透镜体渗透率的研究表明，其渗透率范围在10^{-16}-10^{-14}m^2之间，与本研究的反演结果在同一数量级。低渗透性透镜体的厚度对DNAPL的迁移行为也有重要影响。通过实际钻孔取芯和地质雷达探测相结合的方法确定透镜体的厚度。在实际场地中，选择多个具有代表性的位置进行钻孔取芯，对取出的岩芯进行详细的分析和测量，直接获取低渗透性透镜体在钻孔位置处的厚度。同时，利用地质雷达对场地进行大面积的探测，通过分析雷达反射信号的特征，确定低渗透性透镜体的分布范围和大致厚度。将钻孔取芯和地质雷达探测的结果相结合，得到低渗透性透镜体的厚度分布情况。在本研究中，低渗透性透镜体的厚度在0.5-2米之间变化，不同位置的厚度差异较大。通过对不同厚度低渗透性透镜体对DNAPL迁移影响的数值模拟分析，发现随着透镜体厚度的增加，DNAPL在其表面的堆积量增加，侧向运移距离增大，对DNAPL迁移路径和分布的影响更为显著。低渗透性透镜体的形状通常较为复杂，在数值模拟中，为了简化计算，将其近似为椭圆形。通过对实际场地中低渗透性透镜体的地质数据进行分析，确定其长轴和短轴的长度。利用地质雷达探测数据和钻孔取芯数据，绘制出低渗透性透镜体的二维和三维分布图，通过图像分析和测量工具，获取透镜体的长轴和短轴长度。在本研究中，低渗透性透镜体的长轴长度范围在2-10米之间，短轴长度范围在0.5-3米之间，长轴与短轴的比例在2-4之间变化。为了验证将透镜体近似为椭圆形的合理性，对比了椭圆形透镜体模型和实际复杂形状透镜体模型对DNAPL迁移模拟结果的差异，结果表明在一定误差范围内，椭圆形模型能够较好地反映实际透镜体对DNAPL迁移的影响，为后续的数值模拟研究提供了可靠的基础。低渗透性透镜体在饱和多孔介质中的分布概率是描述其空间分布特征的重要参数。通过对实际场地的地质数据进行统计分析，结合概率统计方法确定其分布概率。收集大量的钻孔数据、地质雷达数据以及其他相关的地质勘探资料，对低渗透性透镜体在不同位置和深度的出现频率进行统计。采用克里金插值法等空间插值方法，将离散的地质数据进行空间插值，得到低渗透性透镜体在整个研究区域内的分布概率图。在本研究中，低渗透性透镜体在饱和多孔介质中的分布呈现出一定的随机性，但在某些区域存在相对较高的分布概率。例如，在靠近河流沉积区域，由于沉积环境的变化，低渗透性透镜体的分布概率相对较高。通过对分布概率的分析，能够更准确地了解低渗透性透镜体在饱和多孔介质中的空间分布规律，为数值模拟中合理设置透镜体的位置提供依据，从而提高模拟结果的准确性和可靠性。3.3模型的建立与验证3.3.1模型的构建过程本研究采用COMSOLMultiphysics软件构建数值模型，该软件具有强大的多物理场耦合模拟能力，能够精确地模拟饱和多孔介质中DNAPL的迁移过程。首先，进行几何建模。在COMSOL软件的几何模块中，创建一个三维矩形区域来代表饱和多孔介质区域。根据实际研究需求，设定该区域的尺寸为长L=10米、宽W=5米、高H=3米。在该区域内，通过布尔运算添加低渗透性透镜体。利用软件的绘图工具，绘制一个椭圆形来表示低渗透性透镜体，根据前期确定的低渗透性透镜体参数，设置椭圆的长轴a=3米、短轴b=1米，并将其放置在饱和多孔介质区域内距离顶部h=1米的位置处，确保透镜体在模型中的位置具有代表性。完成几何建模后，进入材料参数设置环节。在软件的材料库中，分别定义饱和多孔介质、低渗透性透镜体和DNAPL的材料属性。对于饱和多孔介质，输入前文确定的孔隙度\phi=0.38、渗透率k=10^{-12}m^2等参数；对于低渗透性透镜体，设置其孔隙度\phi_{lens}=0.25、渗透率k_{lens}=10^{-15}m^2，以体现其低渗透性的特点；对于DNAPL，输入密度\rho_{DNAPL}=1.462g/cm^3、黏度\mu_{DNAPL}=0.57cP、溶解度S_{DNAPL}=1100mg/L、扩散系数D_{DNAPL}=1.0\times10^{-9}m^2/s等参数，确保材料属性的准确性。接下来进行物理场设置。在COMSOL软件中选择“多孔介质流”模块，该模块能够准确描述流体在多孔介质中的流动行为。在该模块下，启用“达西定律”来描述DNAPL和水在饱和多孔介质中的渗流过程，根据连续性方程和动量方程，设置相应的边界条件和初始条件。在模型的顶部边界，设置为定水头边界，水头值h_0=3米，以模拟地下水的补给情况；在模型的底部边界，设置为自由排水边界，确保流体能够顺利排出；在模型的四周侧面边界，设置为无通量边界，防止流体从侧面流出。对于初始条件，假设模型中初始时刻充满水，DNAPL从模型顶部某一点(x_0,y_0,z_0)=(5,2.5,3)开始注入，初始注入速度v_0=1\times10^{-5}m/s，以便准确模拟DNAPL的迁移起始状态。最后进行网格划分。在COMSOL软件的网格模块中，对整个模型区域进行网格划分。考虑到模型中存在低渗透性透镜体这一复杂结构，为了提高计算精度，采用非结构化网格进行划分。在低渗透性透镜体及其周围区域，进行网格加密处理，使网格更加精细，以准确捕捉DNAPL在透镜体附近的迁移行为；在远离透镜体的区域，适当增大网格尺寸，以减少计算量，提高计算效率。经过多次调试和优化，最终生成的网格数量为N=50000个，确保网格质量满足计算要求，能够准确地模拟DNAPL在饱和多孔介质中的迁移过程。3.3.2模型的验证与校准为了确保所建立的数值模型能够准确反映低渗透性透镜体对饱和多孔介质中DNAPL迁移的影响，采用实验数据对模型进行验证和校准。实验在实验室的砂箱中进行，砂箱尺寸为长1米、宽0.5米、高0.3米。在砂箱中填充均匀的石英砂作为饱和多孔介质，通过实验测定其孔隙度为0.35，渗透率为8\times10^{-13}m^2。在砂箱中设置一个低渗透性透镜体，该透镜体由黏土制成，其孔隙度为0.20，渗透率为5\times10^{-16}m^2，透镜体形状近似椭圆形，长轴为0.3米，短轴为0.1米，位于砂箱内距离顶部0.1米的位置。实验中，将三氯乙烯（TCE）作为DNAPL，从砂箱顶部中心位置以恒定的流速1\times10^{-4}m/s注入。利用CT扫描技术对砂箱内DNAPL的迁移过程进行实时监测，每隔一定时间获取DNAPL在砂箱中的饱和度分布图像，从而得到不同时刻DNAPL的迁移路径和饱和度分布数据。将实验数据与数值模拟结果进行对比分析。首先对比DNAPL的迁移路径，在实验中观察到DNAPL在重力作用下向下迁移，当遇到低渗透性透镜体时，在透镜体表面堆积并开始横向运移，绕过透镜体后继续向下运动，这与数值模拟中DNAPL的迁移路径一致。通过对不同时刻实验图像和模拟结果中DNAPL迁移前沿位置的对比，发现两者的偏差在可接受范围内。例如，在注入t=1000秒时，实验中DNAPL迁移前沿在x方向的位置为x_{exp}=0.45米，模拟结果中为x_{sim}=0.43米，相对误差约为4.4\%。接着对比DNAPL的饱和度分布，将实验中通过CT扫描得到的不同位置的DNAPL饱和度数据与数值模拟结果进行比较。通过绘制实验和模拟的饱和度分布曲线，发现两者在趋势上基本一致，且在关键位置的饱和度数值也较为接近。例如，在低渗透性透镜体上方0.05米处的某点，实验测得的DNAPL饱和度为S_{exp}=0.30，模拟结果为S_{sim}=0.32，相对误差约为6.7\%。通过对比发现数值模拟结果与实验数据存在一定偏差，这可能是由于实验中砂箱内的多孔介质和低渗透性透镜体的实际参数存在一定的不均匀性，以及数值模拟中对一些复杂物理过程的简化处理等原因导致的。为了提高模型的准确性，对模型进行校准。在COMSOL软件中，通过调整饱和多孔介质和低渗透性透镜体的渗透率、孔隙度等参数，以及DNAPL的迁移速度等初始条件，使模拟结果与实验数据达到更好的匹配。经过多次参数调整和模拟计算，最终得到的模拟结果与实验数据的误差明显减小，在迁移路径和饱和度分布等方面的模拟精度得到显著提高，验证了所建立的数值模型的可靠性和准确性，能够有效地用于研究低渗透性透镜体对饱和多孔介质中DNAPL迁移的影响。四、模拟结果与分析4.1DNAPL迁移路径分析4.1.1无低渗透性透镜体时的迁移路径在无低渗透性透镜体的饱和多孔介质模型中，通过数值模拟得到DNAPL的迁移路径如图1所示。从图中可以清晰地看到，DNAPL从注入点开始，在重力作用下迅速向下迁移。由于饱和多孔介质的渗透率相对均匀，DNAPL在迁移过程中基本沿着垂直方向向下渗透，形成较为规则的羽流形状。在迁移初期，DNAPL的迁移速度较快，随着迁移距离的增加，由于黏滞力的作用，迁移速度逐渐减缓。为了更准确地分析迁移路径，对不同时刻DNAPL迁移前沿的位置进行了统计分析。在注入后的第10天，DNAPL迁移前沿在垂直方向上达到了距离注入点0.5米的位置；到第30天，迁移前沿下降到了1.2米处；第60天，迁移前沿进一步下降到了2.0米的深度。通过对这些数据的分析，可以发现DNAPL的迁移距离与时间呈现出一定的非线性关系，迁移速度随着时间的推移逐渐降低。这是因为随着DNAPL的迁移，其与饱和多孔介质的接触面积不断增大，受到的黏滞阻力也逐渐增大，从而导致迁移速度减慢。而且，由于多孔介质的孔隙结构存在一定的随机性，DNAPL在迁移过程中会出现一些微小的偏移，但整体上仍然保持着垂直向下的迁移趋势。4.1.2存在低渗透性透镜体时的迁移路径当模型中存在低渗透性透镜体时，DNAPL的迁移路径发生了显著变化。从图2可以看出，DNAPL在重力作用下向下迁移，当遇到低渗透性透镜体时，由于透镜体的渗透率远低于周围饱和多孔介质，DNAPL无法克服其毛管压力，从而停止垂向入渗，并在透镜体表面堆积。随着注入的DNAPL不断增加，堆积在透镜体表面的DNAPL逐渐增多，开始沿着透镜体表面进行横向运移。在横向运移过程中，DNAPL会寻找透镜体周围渗透率相对较高的区域，绕过透镜体继续向下迁移。同样对存在低渗透性透镜体时DNAPL迁移前沿的位置进行统计分析。在注入后的第10天，DNAPL迁移至低渗透性透镜体表面并开始堆积，此时迁移前沿在垂直方向上距离注入点0.3米；到第30天，DNAPL在透镜体表面堆积的同时进行横向运移，横向运移距离达到了0.2米，而垂直方向上迁移前沿仍停留在透镜体表面；第60天，DNAPL成功绕过透镜体，继续向下迁移，此时垂直方向上迁移前沿距离注入点1.5米，横向运移距离达到了0.5米。与无低渗透性透镜体时相比，DNAPL的迁移路径变得更加复杂，迁移速度明显降低，尤其是在透镜体附近，迁移过程出现了明显的停滞和转向。这种迁移路径的变化表明低渗透性透镜体对DNAPL的迁移具有显著的阻碍和引导作用，它改变了DNAPL的迁移方向，使得DNAPL在透镜体周围形成了特殊的分布形态，增加了污染范围的不确定性，对地下水污染的防控和修复带来了更大的挑战。4.2DNAPL饱和度分布分析4.2.1不同时刻的饱和度分布通过数值模拟，得到了不同时刻DNAPL在饱和多孔介质中的饱和度分布情况，如图3所示。在注入初期（t=10天），DNAPL刚刚进入饱和多孔介质，其饱和度在注入点附近较高，随着距离注入点距离的增加而迅速降低。此时，DNAPL主要在重力作用下开始向下迁移，在迁移路径上，由于孔隙结构的影响，饱和度分布存在一定的不均匀性，但整体上呈现出以注入点为中心，向下逐渐扩散的趋势。随着时间的推移（t=30天），DNAPL继续向下迁移，饱和度高值区域逐渐扩大。在迁移过程中，DNAPL受到多孔介质的黏滞力和毛细力作用，部分DNAPL滞留在孔隙中，导致饱和度分布更加不均匀。在一些孔隙较大、连通性较好的区域，DNAPL能够顺利通过，饱和度相对较高；而在孔隙较小、连通性较差的区域，DNAPL容易被滞留，饱和度较低。当注入时间达到60天时，DNAPL在饱和多孔介质中的迁移范围进一步扩大。在模型底部，DNAPL开始逐渐聚集，形成相对较高饱和度的区域。而且，由于多孔介质的非均质性，在一些局部区域出现了饱和度的异常变化，这可能是由于局部孔隙结构的特殊性或与周围介质的相互作用导致的。通过对不同时刻饱和度分布的分析，可以清晰地了解DNAPL在饱和多孔介质中的迁移和扩散过程，为后续研究低渗透性透镜体对其饱和度分布的影响提供了基础。4.2.2低渗透性透镜体对饱和度分布的影响当存在低渗透性透镜体时，DNAPL的饱和度分布发生了显著改变。从图4可以看出，当DNAPL迁移至低渗透性透镜体时，由于透镜体的阻碍作用，DNAPL在其表面堆积，导致透镜体上方的饱和度明显升高。在透镜体表面，DNAPL形成了一个高饱和度的区域，其饱和度远高于周围饱和多孔介质中的DNAPL饱和度。这是因为低渗透性透镜体的渗透率极低，DNAPL难以穿过，只能在其表面聚集，使得该区域的DNAPL含量增加，饱和度相应提高。随着时间的推移，堆积在透镜体表面的DNAPL开始横向运移，在横向运移过程中，DNAPL的饱和度逐渐降低。这是因为DNAPL在横向运移时，不断与周围的水和多孔介质相互作用，部分DNAPL被分散和稀释，导致饱和度下降。而且，由于透镜体周围的多孔介质渗透率存在差异，DNAPL在横向运移过程中的饱和度分布也不均匀。在渗透率较高的区域，DNAPL能够较快地通过，饱和度相对较高；而在渗透率较低的区域，DNAPL运移受阻，饱和度较低。在DNAPL绕过低渗透性透镜体继续向下迁移后，其饱和度分布也受到了透镜体的影响。与无透镜体时相比，绕过透镜体后的DNAPL在向下迁移过程中，饱和度分布出现了一定的偏差。这是因为透镜体改变了DNAPL的迁移路径，使得DNAPL在重新向下迁移时，进入了不同的孔隙通道，与周围介质的相互作用发生了变化，从而导致饱和度分布发生改变。低渗透性透镜体对DNAPL饱和度分布的影响显著，它改变了DNAPL在饱和多孔介质中的分布格局，使得饱和度分布更加复杂，对地下水污染的程度和范围评估带来了更大的困难。4.3DNAPL浓度分布分析4.3.1浓度场的演变过程通过数值模拟得到了不同时刻DNAPL在饱和多孔介质中的浓度场分布，其演变过程清晰地展示了DNAPL在地下环境中的扩散规律。在注入初期（t=1天），DNAPL以较高浓度在注入点附近聚集，形成一个浓度高值中心。此时，DNAPL主要受重力作用开始向下迁移，但由于迁移时间较短，扩散范围有限，浓度等值线呈近似圆形分布，且高浓度区域较为集中。随着时间的推移（t=5天），DNAPL继续向下迁移，浓度高值区域逐渐向下扩展，浓度等值线开始呈现出椭圆形。在迁移过程中，DNAPL与周围的水和多孔介质发生相互作用，部分DNAPL溶解于水中，导致其浓度逐渐降低。而且，由于多孔介质的非均质性，DNAPL在迁移过程中会出现局部的浓度差异，在孔隙较大、连通性较好的区域，DNAPL迁移速度较快，浓度相对较高；而在孔隙较小、连通性较差的区域，DNAPL迁移受阻，浓度较低。当注入时间达到10天时，DNAPL的扩散范围进一步扩大，浓度高值区域向下延伸更为明显。此时，浓度等值线的形状变得更加不规则，这是因为DNAPL在迁移过程中受到了多种因素的综合影响，如重力、毛细力、黏滞力以及多孔介质的非均质性等。在模型底部，由于DNAPL的逐渐聚集，浓度开始升高，形成一个相对较高浓度的区域。而且，在DNAPL迁移的前沿，浓度梯度较大，这表明DNAPL在不断向周围扩散，污染范围在逐渐扩大。通过对浓度场演变过程的分析，可以直观地了解DNAPL在饱和多孔介质中的迁移和扩散行为，为后续研究低渗透性透镜体对其浓度分布的影响提供了基础。4.3.2低渗透性透镜体对浓度分布的影响当存在低渗透性透镜体时，DNAPL的浓度分布发生了显著变化。从图5可以看出，当DNAPL迁移至低渗透性透镜体时，由于透镜体的阻碍作用，DNAPL在其表面堆积，导致透镜体上方的DNAPL浓度明显升高。在透镜体表面，形成了一个高浓度的区域，其浓度远远高于周围饱和多孔介质中的DNAPL浓度。这是因为低渗透性透镜体的渗透率极低，DNAPL难以穿过，只能在其表面聚集，使得该区域的DNAPL含量增加，浓度相应提高。随着时间的推移，堆积在透镜体表面的DNAPL开始横向运移，在横向运移过程中，DNAPL的浓度逐渐降低。这是因为DNAPL在横向运移时，不断与周围的水和多孔介质相互作用，部分DNAPL被分散和稀释，导致浓度下降。而且，由于透镜体周围的多孔介质渗透率存在差异，DNAPL在横向运移过程中的浓度分布也不均匀。在渗透率较高的区域，DNAPL能够较快地通过，浓度相对较高；而在渗透率较低的区域，DNAPL运移受阻，浓度较低。在DNAPL绕过低渗透性透镜体继续向下迁移后，其浓度分布也受到了透镜体的影响。与无透镜体时相比，绕过透镜体后的DNAPL在向下迁移过程中，浓度分布出现了一定的偏差。这是因为透镜体改变了DNAPL的迁移路径，使得DNAPL在重新向下迁移时，进入了不同的孔隙通道，与周围介质的相互作用发生了变化，从而导致浓度分布发生改变。低渗透性透镜体对DNAPL浓度分布的影响显著，它改变了DNAPL在饱和多孔介质中的浓度分布格局，使得浓度分布更加复杂，增加了污染范围和程度评估的难度。这对于地下水污染的防控和修复具有重要的指导意义，在实际的污染治理中，需要充分考虑低渗透性透镜体对DNAPL浓度分布的影响，制定更加科学合理的治理方案。五、影响因素分析5.1低渗透性透镜体参数的影响5.1.1渗透率的影响低渗透性透镜体的渗透率是影响DNAPL迁移的关键参数之一，其对DNAPL迁移速度、路径和饱和度有着显著影响。通过数值模拟，设置了不同渗透率的低渗透性透镜体，分别为10^{-16}m^2、10^{-15}m^2和10^{-14}m^2，研究其对DNAPL迁移的影响。当低渗透性透镜体渗透率为10^{-16}m^2时，DNAPL在迁移过程中遇到透镜体后，几乎完全无法穿透，迅速在透镜体表面堆积。在堆积过程中，DNAPL的迁移速度急剧下降，由原来在周围饱和多孔介质中的相对较快迁移速度，变为近乎停滞状态。例如，在未遇到透镜体时，DNAPL在10天内可以迁移0.5米的距离；而遇到该渗透率的透镜体后，在相同时间内，DNAPL在透镜体表面的迁移距离仅为0.05米。随着时间的推移，堆积在透镜体表面的DNAPL逐渐增多，开始沿着透镜体表面进行横向运移，但横向运移速度也相对较慢，在30天内，横向运移距离仅达到0.2米。在这种情况下，DNAPL的迁移路径被透镜体完全改变，原本垂直向下的迁移路径在遇到透镜体后，转变为在透镜体表面的堆积和横向运移，然后再寻找机会绕过低渗透性透镜体继续向下迁移。由于DNAPL在透镜体表面的大量堆积，使得透镜体上方区域的DNAPL饱和度显著升高，形成一个高饱和度区域，其饱和度可达到0.6以上。当低渗透性透镜体渗透率增加到10^{-15}m^2时，DNAPL在遇到透镜体后，虽然仍会受到较大阻碍，但部分DNAPL能够缓慢穿透透镜体。此时，DNAPL在透镜体表面的堆积速度有所减缓，迁移速度相对提高。在10天内，DNAPL在透镜体表面的迁移距离可达到0.1米。在横向运移方面，30天内横向运移距离达到0.3米。DNAPL的迁移路径仍然受到透镜体的明显影响，但与渗透率为10^{-16}m^2时相比，绕过低渗透性透镜体继续向下迁移的速度有所加快。在透镜体上方，DNAPL的饱和度虽然也会升高，但升高幅度相对较小，饱和度约为0.5。当低渗透性透镜体渗透率进一步增加到10^{-14}m^2时，DNAPL能够更顺利地穿透透镜体，在透镜体表面的堆积现象明显减弱。在10天内，DNAPL在透镜体表面的迁移距离可达到0.2米，且有较多DNAPL能够在较短时间内绕过低渗透性透镜体继续向下迁移。此时，DNAPL的迁移路径受透镜体的影响相对较小，迁移速度接近无透镜体时的情况。在透镜体上方，DNAPL的饱和度升高幅度较小，饱和度约为0.4。通过对不同渗透率低渗透性透镜体下DNAPL迁移的分析可知，随着低渗透性透镜体渗透率的增加，DNAPL在透镜体表面的堆积量逐渐减少，迁移速度逐渐加快，迁移路径受透镜体的影响逐渐减弱。低渗透性透镜体的渗透率对DNAPL的迁移有着至关重要的影响，在研究饱和多孔介质中DNAPL的迁移行为时，必须充分考虑低渗透性透镜体渗透率的变化。5.1.2厚度的影响低渗透性透镜体的厚度是影响DNAPL迁移和分布的重要因素之一，其变化会对DNAPL在饱和多孔介质中的迁移行为产生显著影响。通过数值模拟，设置了不同厚度的低渗透性透镜体，分别为0.5米、1米和1.5米，研究其对DNAPL迁移和分布的影响。当低渗透性透镜体厚度为0.5米时，DNAPL在迁移过程中遇到透镜体后，在透镜体表面堆积并开始横向运移。由于透镜体厚度较薄，DNAPL相对较容易绕过低渗透性透镜体继续向下迁移。在堆积阶段，DNAPL在透镜体表面堆积的量相对较少，堆积速度也相对较快。例如，在注入后的10天内，DNAPL在透镜体表面堆积的高度约为0.1米。在横向运移阶段，30天内横向运移距离可达到0.3米。在绕过低渗透性透镜体继续向下迁移时，迁移速度较快，在60天内，DNAPL在垂直方向上迁移的距离可达到1.8米。在饱和度分布方面，由于DNAPL在透镜体表面堆积量较少，透镜体上方区域的DNAPL饱和度升高幅度相对较小，饱和度约为0.45。当低渗透性透镜体厚度增加到1米时，DNAPL在透镜体表面的堆积量明显增加，堆积时间也相应延长。在注入后的10天内，DNAPL在透镜体表面堆积的高度约为0.2米。由于透镜体厚度增加，DNAPL绕过低渗透性透镜体继续向下迁移的难度增大，横向运移距离相对减小，在30天内横向运移距离为0.2米。在垂直方向上，60天内DNAPL迁移的距离为1.5米，迁移速度明显低于透镜体厚度为0.5米时的情况。在饱和度分布方面，透镜体上方区域的DNAPL饱和度升高幅度较大，饱和度约为0.55。当低渗透性透镜体厚度进一步增加到1.5米时，DNAPL在透镜体表面堆积的量达到最大，堆积时间也最长。在注入后的10天内，DNAPL在透镜体表面堆积的高度约为0.3米。此时，DNAPL绕过低渗透性透镜体继续向下迁移变得更加困难，横向运移距离进一步减小，在30天内横向运移距离仅为0.1米。在垂直方向上，60天内DNAPL迁移的距离仅为1.2米，迁移速度最慢。在饱和度分布方面，透镜体上方区域的DNAPL饱和度升高幅度最大，饱和度可达到0.65。通过对不同厚度低渗透性透镜体下DNAPL迁移和分布的分析可知，随着低渗透性透镜体厚度的增加，DNAPL在透镜体表面的堆积量逐渐增多，堆积时间逐渐延长，绕过低渗透性透镜体继续向下迁移的难度逐渐增大，迁移速度逐渐降低。而且，透镜体上方区域的DNAPL饱和度随着透镜体厚度的增加而升高。低渗透性透镜体的厚度对DNAPL的迁移和分布有着重要影响，在研究饱和多孔介质中DNAPL的迁移行为时，必须充分考虑低渗透性透镜体厚度的变化。5.1.3形状与分布的影响低渗透性透镜体的形状与分布是影响DNAPL迁移的重要因素，其多样性会导致DNAPL迁移行为的复杂性。通过数值模拟，研究了不同形状（椭圆形、矩形、不规则形状）和分布（均匀分布、随机分布）的低渗透性透镜体对DNAPL迁移的综合影响。当低渗透性透镜体为椭圆形时，DNAPL在迁移过程中遇到透镜体后，在其表面堆积并开始横向运移。由于椭圆形透镜体的边界较为平滑，DNAPL在横向运移时相对较为顺畅，能够沿着透镜体表面的曲率逐渐扩散。在堆积阶段，DNAPL在透镜体表面的堆积量相对较为均匀，不会出现局部堆积过多或过少