一种寒冷区域内多孔介质中LNAPL迁移模拟方法与流程

一种寒冷区域内多孔介质中lnapl迁移模拟方法
技术领域
1.本发明涉及污染水文地质领域，尤其涉及一种寒冷区域内多孔介质中lnapl迁移模拟方法。


背景技术：

2.随着石油化工业的发展，非相有机污染物(napl)在其炼制、贮存及运输过程中泄漏至环境，对人类健康造成严重威胁，是土壤及地下水污染治理的重要对象。我国季节性冻土面积广阔，北方相当长的时间处于冬季冻结状态。建立冻结条件下饱和-非饱和带的lnapl运移模型有利于预测冻土区lnapl泄漏后在土壤-地下水中的分布迁移，有利于寒区的lnapl污染修复治理。
3.数值模拟作为一种预测napl在土壤—地下水中迁移分布的重要手段，在napl污染防治及油气开采中都得到了广泛应用，目前关于三相流在饱和-非饱和带中的模拟已经相对成熟。有关在饱和带中不同降雨速率、介质非均质性及napl泄漏速率、地下水流速下的napl两相流模型已比较完善[施小清,吴吉春,刘德朋,等.饱和介质中重非水相液体运移的数值模拟及敏感性分析[j].南京大学学报(自然科学版),2011,47(3):299-307.yang z b,zandin h,niemi a,et al.the role of geological heterogeneity and variability in water infiltration on non-aqueous phase liquid migration[j].environmental earth sciences,2013,68(7):2085-2097.郑菲,高燕维,施小清,等.地下水流速及介质非均质性对重非水相流体运移的影响[j].水利学报,2015,46(8):925-933.]。而在非饱和带，一些二维实验研究了napl的运移，并对运移过程进行了数值模拟。napl在包气带的迁移涉及到水相-气相-napl相三相共存，通常的模拟方法是使用三相的达西定律并将两两相间的压强关系通过毛细压力联系起来。求解三相流方程所需的相对渗透率、有效饱和度与毛细压力的关系模型一般采用vg经验关系式以及brooks-corey关系式。对于三相间的物质质量转换则采用局部平衡假设关系式。lenhard等通过建立自由相、残余相和截留相napl与流体水头的函数关系式建立了napl在非饱和带-饱和带的迁移模型，预测了其分布及渗透率,但未进行该模型实验校正[lenhard r j,rayner j l,davis g b.a practical tool for estimating subsurface lnapl distributions and transmissivity using current and historical fluid levels in ground water wells:effects of entrapped and residual lnapl[j].journal of contaminant hydrology,2017,205:1
–
11.]。之后进行的补充实验验证中考虑了滞留作用，comsol multiphysics软件被用于三相流的数值模拟。
[0004]
迄今，寒区冻结条件下lnapl的迁移研究仍相对较少。少数学者开展了实验研究。研究表明lnapl在冻结条件下泄漏后在土壤包气带及地下水中的迁移因为冰晶影响会与非冻结条件下的迁移图像有所不同，冰晶有效阻止石油垂直迁移，加强石油沿优先流路径下渗及在地表浅层水平流动，但即使完全冰饱和的土壤也表现出一定的冰渗透性[biggar k w.site investigations of fuel spill migration into permafrost[j].journal of cold regions engineering,1998,12(2):84
–
104.chuvilin e m,naletova n s,
miklyaeva e c,et al.factors affecting spreadability and transportation of oil in regions of frozen ground[j].polar record,2001,37(202):229-238.barnes d l,wolfe s m.influence of ice on the infiltration of petroleum into frozen coarse-grained soil[j].petroleum science and technology,2008,26(7-8):856-867.]。且在冻结条件下，土壤包气带伴随冻结滞水的形成，饱和水汽压作用减小，而毛管作用得以加强[徐树林,那平山,武俊英.关于冻结滞水的探讨[j].内蒙古林学院报,1994,(2):46-51.],地下水中石油污染物及lnapl ganglia在冻结势和基质势作用下会向冻结峰积聚，且有部分石油汇聚于表层[singh k,niven r k.non-aqueous phase liquid spills in freezing and thawing soils:critical analysis of pore-scale processes[j].critical reviews in environmental science and technology,2013,43(6):551-597.李兴柏,李国玉.温度梯度对多年冻土区石油迁移影响的研究[j].甘肃科学学报,2013,25(01):73-76.]。国内也有少数学者通过室内冻结箱实验开展lnapl于冻结条件下的迁移规律研究[多孔介质内油、水介质迁移对比低温实验[j].当代化工,2016,45(06):1108-1111.]。然而由于冻结固相(冰)的存在，需将温度场与lnapl迁移所涉及的渗流场、化学场耦合，同时多相多场耦合问题求解往往不易收敛，有关冻结条件下气相、水相、napl相三相运移模拟尚未见相关报道。
[0005]
常见的刻化常温非冻结条件下napl于饱和带、非饱和带中迁移分布的模拟方法。由于冻结固相(冰)的存在，需将温度场与lnapl迁移所涉及的渗流场、化学场耦合，同时多相多场耦合问题求解往往不易收敛，现有的模拟方法往往不考虑冻结作用对lnapl迁移的影响，仅适用于我国南方非冻结地区的napl污染的修复指导，无法为北方冻土区冬季napl泄漏的修复提供指导。


技术实现要素：

[0006]
发明目的：本发明旨在提供一种固相-气相-水相-napl相四相共存，温度场、流场、化学场三场耦合的寒冷区域内多孔介质中lnapl迁移模拟方法。
[0007]
技术方案：本发明的寒冷区域内多孔介质中lnapl迁移模拟方法，基于comsol multiphysics模拟多孔介质水-气-napl的三相流动模型，包括以下步骤：
[0008]
(1)通过建立lnapl迁移和lnapl的存在状态及相互转化过程的控制方程构建流场；所述控制方程包括基于质量守恒定律的四相控制方程，自由相lnapl在水中的溶解的过程表达方程，自由相lnapl挥发于空气的过程表达方程，溶于水的lnapl挥发于空气的过程表达方程和溶于水中的lnapl吸附于土壤骨架的过程表达方程；基于van genuchten方程和parker模型的毛细压力、相对渗透率、饱和度关系方程，
[0009]
(2)通过建立多孔介质传热方程构建温度场，在所述三相流的源项增设冰的融化潜热项并基于热量守恒建立多孔介质传热方程；
[0010]
(3)通过未冻水和冻结冰的经验关系式对温度场和流场的耦合；所述经验关系式为：
[0011][0012]
式中：为比例系数，为起始冻结温度(绝对温度)，系数1.1为水与冰的密度之比；
[0013]
(4)通过建立溶解相-挥发相lnapl质量平衡方程构建化学场并通过质量转换项实现与流场的耦合；
[0014]
(5)采用comsolmultiphysics对步骤(1)～(4)中所述的表达方程计算求解，得到lnapl迁移规律。
[0015]
进一步地，步骤(1)中，所述四相控制方程为：
[0016][0017][0018][0019]
式中：ε、ε1为冻结时和未冻时的孔隙度，sw、sn、sg分别为水相、lnapl相、气相的饱和度，si为冰相相对饱和度，ρi、ρw、ρn、ρg分别为冰相、水相、lnapl相、气相的密度，vw、vn、vg分别为水相、lnapl相、气相的达西流速，分别lnapl溶于水、lnapl挥发至气相、溶于水的lnapl挥发于气相、溶于水的lnapl吸附于土壤骨架的速率。
[0020]
进一步地，步骤(1)中，所述毛细压力、相对渗透率、饱和度关系方程中k-s关系式为：
[0021][0022]
式中，k
rw
、k
rg
、k
rn
分别为水、气、lnapl的相对渗透率，srw、srg、srn分别为水、气、lnapl的相对饱和度，srtw为总的液相饱和度。
[0023]
进一步地，步骤(1)中，所述毛细压力、相对渗透率、饱和度关系方程中s-p关系式为：
[0024][0025]
[0026][0027][0028]
式中，s
wr
、s
nr
为残余水饱和度及残余lnapl饱和度，γ
gn
、γ
nw
为缩放因子，δ
gn
、δ
nw
分别为气-lnapl表面张力及lnapl-水表面张力，p
cnw
为napl—水毛管压力，p
cgn
为napl—气毛管压力。α、n、m为vg方程参数。
[0029]
进一步地，步骤(1)中，自由相lnapl在水中的溶解的过程表达方程为：
[0030][0031]
所述自由相lnapl挥发于空气的过程表达方程为：
[0032][0033]
所述溶于水的lnapl挥发于空气的过程表达方程为：
[0034][0035]
所述溶于水中的lnapl吸附于土壤骨架的过程表达方程为：
[0036][0037]
式中，ε为孔隙度，sw为气相的相对饱和度，sg为气相的相对饱和度，λ
nw
、λ
ng
、为传质系数，为lnapl在水中的溶解度，c1为lnapl溶于水的浓度；为lnapl在空气中的溶解度，c2为lnapl在空气中的浓度，h为亨利常数，ρs为土体密度，m为lnapl的质量分数，kd为分布系数。
[0038]
更进一步地，所述分布系数kd的经验公式为：
[0039]
kd＝f
ockoc
[0040]
式中，f
oc
为土体中有机碳的质量分数，k
oc
为有机碳分配系数。
[0041]
进一步地，步骤(2)中，所述多孔介质传热方程为：
[0042][0043]
式中，ρ、ρ
p
分别为流体和固体骨架密度，c
p
分别比定压热容，u为达西流速，t为温度，k
eff
为热导率，lf为融化潜热。
[0044]
更进一步地，步骤(4)中，所述质量平衡方程包括lnapl溶解相的溶质运移方程和lnapl挥发相的溶质运移方程，所述lnapl溶解相的溶质运移方程为：
[0045][0046][0047]
式中，m为lnapl的摩尔质量分数，dw、dg为水动力弥散系数。
[0048]
本发明通过在以往三相流运移方程的水相运移质量守恒方程中加入冰项，通过未冻水和冻结冰的经验关系式实现温度场和流场的耦合，并将冰视为土壤固体骨架来实现水-气-napl-固四相向水-气-napl三相的转换，将传热方程中的比热容和导热系数设为水-气-napl-固四相体积占比的加权函数。从而建立了一套固相-气相-水相-napl相四相共存，温度场、流场、化学场三场耦合模拟方法。
[0049]
有益效果：与现有技术相比，本发明具有如下显著优点：本发明可用于模拟及定量刻化lnapl泄漏过程中发生冻结后lnapl的迁移分布，为寒区lnapl泄漏后的修复提供参考。
附图说明
[0050]
图1为本发明的模型设置与初始条件、边界条件；
[0051]
图2为本发明的温度场随时间的演变过程；
[0052]
图3为本发明的含冰率随时间的演变过程；
[0053]
图4为本发明的冻结条件下lnapl饱和度变化图；
[0054]
图5为本发明的冻结条件下水饱和度变化图；
[0055]
图6为本发明的冻结条件下空气饱和度变化图；
[0056]
图7为本发明的lnapl运移量分布图；
[0057]
图8为本发明的常温条件下lnapl饱和度变化图；
[0058]
图9为本发明的冻结条件下lnapl溶解相浓度分布图
[0059]
图10为本发明的冻结条件下lnapl挥发相浓度分布图；
[0060]
图11为本发明的冻结条件不同状态lnapl含量变化图；
[0061]
图12为本发明的常温条件下lnapl溶解相浓度分布图；
[0062]
图13为本发明的常温条件下lnapl挥发相浓度分布图；
[0063]
图14为本发明的常温条件不同状态lnapl含量变化图。
具体实施方式
[0064]
下面结合附图对本发明的技术方案作进一步说明。
[0065]
本发明首先构建冻结条件下多孔介质中lnapl迁移的数学模型，包括描述流动场的水-气-napl-固(冰)四相控制方程、刻化温度场的多孔介质传热方程以及刻化化学场的溶解相-挥发相lnapl质量平衡方程。在数学模型的构建中，因冻结条件下多孔介质lnapl迁移因有冻结的固相即冰存在，需要在以往三相流运移方程的水相运移质量守恒方程中加入冰项，并通过未冻水和冻结冰的经验关系式实现温度场和流场的耦合，然后将其在comsol软件中进行设置，并根据模拟情景设置相关参数、初始条件和边界条件进行求解。
[0066]
本发明的寒冷区域内多孔介质中lnapl迁移模拟方法，即构建冻结条件下多孔介质中lnapl迁移的数学模型，具体步骤如下：
[0067]
(1)构建流场
[0068]
冻结条件下lnapl的迁移包括两个过程，一个是lnapl自由相以连续相在水-气-lnapl-冰四相共存时的毛细压力作用下迁移，另一个是溶于水的lnapl及其挥发在空气中的lnapl的溶质迁移过程，以对流-弥散作用为主。该过程涉及水-冰相变、水-气-napl三相流迁移、多孔介质温度传导及溶质运移模型耦合，其控制方程如下：
[0069]
1)水-气-napl-固(冰)四相控制方程
[0070]
对于水-气-napl-固(冰)相控制流动方程据各相的质量守恒定律得到：
[0071][0072][0073][0074]
式中：ε、ε1为冻结时和未冻时的孔隙度，sw、sn、sg分别为水相、lnapl相、气相的饱和度，si为冰相相对饱和度，ρi、ρw、ρn、ρg分别为冰相、水相、lnapl相、气相的密度，vw、vn、vg分别为水相、lnapl相、气相的达西流速，分别lnapl溶于水、lnapl挥发至气相、溶于水的lnapl挥发于气相、溶于水的lnapl吸附于土壤骨架的速率。lnapl在非饱和带土壤的迁移涉及到水相、气相、lnapl相、土壤固相的质量转换过程，考虑包括lnapl自由相溶于水，lnapl自由相挥发于空气，溶于水中的lnapl挥发于空气及溶于水中的lnapl吸附于土壤骨架的四个质量转换过程。
[0075]
自由相lnapl在水中的溶解过程可用下式表达：
[0076][0077]
式中：ε为孔隙度，sw为气相的相对饱和度，λ
nw
为传质系数，为lnapl在水中的溶解度，c1为lnapl溶于水的浓度。
[0078]
自由相lnapl挥发于空气可用下式表达：
[0079][0080]
式中：ε为孔隙度，sg为气相的相对饱和度，λ
ng
为传质系数，为lnapl在空气中的溶解度，c2为lnapl在空气中的浓度。
[0081]
溶于水的lnapl挥发于空气可用下式表达：
[0082][0083]
式中：ε为孔隙度，sw为水相的相对饱和度，为传质系数，c1为lnapl溶于水的浓度，c2为lnapl在空气中的浓度，h为亨利常数。
[0084]
溶于水中的lnapl吸附于土壤骨架可用下式表示：
[0085][0086]
式中，ρs为土体密度，sw为水相的相对饱和度，c1为lnapl溶于水的浓度，m为lnapl的质量分数，kd为分布系数，kd的经验公式为：
[0087]
kd＝f
ockoc
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0088]
式中，f
oc
为土体中有机碳的质量分数，k
oc
为有机碳分配系数。
[0089]
而对于毛管压力p
cnw
、p
cgn
以及相对渗透率k
rw
、k
rn
、k
rg
，需构建k-s-p关系式。本模型k-s关系将van genuchten方程拓展到三相流，得到：
[0090][0091]
式中，srw、srg、srn分别为水、气、lnapl的相对饱和度，srtw为总的液相饱和度。其与毛管压力的s-p关系式采用parker模型，表达式如下：
[0092][0093][0094][0095][0096]
式中，s
wr
、s
nr
为残余水饱和度及残余lnapl饱和度，γ
gn
、γ
nw
为缩放因子，δ
gn
、δ
nw
分别为气-lnapl表面张力及lnapl-水表面张力，p
cnw
为napl—水毛管压力，p
cgn
为napl—气毛管压力。。α、n、m为vg方程参数。
[0097]
而对于饱和渗透率k，它随冻结过程孔隙度的变化而变化(将冰视为土壤固体骨架来实现水-气-napl-固四相向水-气-napl三相的转换)，采用kozeny-carman模型，表达式如下：
[0098][0099]
式中，d为颗粒直径。
[0100]
(2)构建温度场
[0101]
多孔介质传热方程即热量守恒方程如下：
[0102][0103]
式中，ρ、ρ
p
分别为流体和固体骨架密度，c
p
分别比定压热容，u为达西流速，t为温度，k
eff
为热导率，lf为融化潜热(zhang等，2016)。
[0104]
(3)温度场和流场的耦合
[0105]
为了实现温度场与三相流场的耦合，一方面需要引入冻结条件下未冻水和冻结冰的经验比例关系式，另一方面，通过将冰视为土壤固体骨架来实现水-气-napl-固四相向
水-气-napl三相的转换，即冻结过程中固体骨架的孔隙度将随冰的形成而不断变化，且传热方程中的比热容和导热系数是水-气-napl-固四相体积占比的加权函数。
[0106]
冻结场和三相流场的耦合关键即未冻水与冻结冰的经验比例式如下：
[0107][0108]
式中，b(t)为比例系数，tf为起始冻结温度(绝对温度)，系数1.1为水与冰的密度之比。
[0109]
(4)构建化学场
[0110]
建立溶解相-挥发相lnapl质量平衡方程，所述溶解相-挥发相lnapl质量平衡方程，包括lnapl溶解相的溶质运移方程和lnapl挥发相的溶质运移方程。
[0111]
lnapl溶解相的溶质运移方程为：
[0112][0113]
lnapl挥发相的溶质运移方程为：
[0114][0115]
式中，m为lnapl的摩尔质量分数，dw、dg为水动力弥散系数。
[0116]
非冻结条件下lnapl迁移的控制方程除水相运移方程(方程1)不包含含冰项及无需进行多孔介质温度传导模块设置外，其余控制方程与冻结条件下lnapl迁移控制方程一致。
[0117]
(5)模型求解
[0118]
采用comsol multiphysics的多孔介质多相流模块、多孔介质温度传导模块、多孔介质物质传递模块及达西定律模块输入并求解冻结条件下lnapl迁移控制方程。comsol multiphysics是一款强大的多物理场耦合软件，常用于多场多相问题的求解。而冻结条件下lnapl的迁移，涉及到水、气、napl三相在多孔介质中运移，是温度场、流场及化学场高度耦合的问题。comsol中水相方程并未包含含冰项，故在模拟冻结时需要在水相质量守恒方程中加入含冰项随时间变化(方程1)。而气相、lnapl相则采用其内置方程。comsol内置的多孔介质相传递模块可用来模拟三相流运移，但其未考虑四相质量转换过程，需要在模块中增设(方程4-8)。毛管压力以及相对渗透率也需要在comsol中自行构建其与相对饱和度间的本构关系(方程9-13)。多孔介质温度传导模块可用来模拟介质传热，但其热传导方程缺少冰的融化潜热项，需要在源项增设(方程15)。而且其缺少冻结方程设置模块，需将冻结条件下未冻结水和冻结冰的经验比例关系式(方程16)在comsol中进行设置。
[0119]
对于非冻条件下lnapl迁移模型，可采用comsol的多孔介质多相流模块、多孔介质稀物质传递模块及达西定律模块输入并求解lnapl迁移控制方程。同样将四相质量转换过
程在多孔介质相传递模块中进行增设。而三相流运移控制方程均采用comsol内置方程。
[0120]
comsol通过有限单元法离散上述模型，以求解冻结、非冻结条件下lnapl迁移的非稳定问题。冻结、非冻结条件下lnapl迁移模拟输出结果中自由相lnapl的含量用饱和度表示，溶解相及挥发相lnapl含量用质量浓度表示。
[0121]
通过以下算例来验证本发明的有效性：
[0122]
(1)算例设置
[0123]
1)算例概述及初始边界条件
[0124]
冻结条件下lnapl在非饱和带中二维迁移数值算例见图1。研究区为二维xy剖面，长1.5m，厚1.5m，下部0.5m为地下水饱和区，取地下水面为基准面。区域网格剖分采用comsol的物理场细化网格剖分。
[0125]
至于初始条件，对于非饱和带与饱和带分别进行如下设置：
[0126][0127]
模型的边界条件如下设置：
[0128]
①
流场
[0129][0130]
②
化学场
[0131][0132]
③
温度场
[0133][0134]
所有的初始边界条件如图1。而对于非冻结条件下，无需设置温度场边界及初始条件，流场和溶质场边界条件及初始条件与冻结条件下保持一致。
[0135]
2)基本参数
[0136]
模型中采用的参数值如表1，包括土壤、水、lnapl、冰的特征参数以及三相质量转换的相关参数。lnapl在此以甲苯为模拟对象。
[0137]
表1模型中的参数值
[0138]
[0139][0140]
3)情景设置
[0141]
设置冻结条件与非冻结条件lnapl下渗两种情景进行模拟，并通过模拟结果对比
来探究冻结条件对lnapl迁移影响。非冻结条件下除温度场外，其它条件与冻结条件下保持一致。
[0142]
(2)结果与分析
[0143]
1)温度场
[0144]
冻结过程中温度随冻结时间增长在地表浅层由0℃逐渐降低至-14℃(图2)，冰饱和度在地表浅层逐渐增加，冻结12h时，地表浅层冰饱和度达到0.8(图3)。因为温度只在上部薄层发生变化，下部恒为5℃，故取土壤表层0.1m厚的温度场模拟结果进行分析。可知冻结场在很短时间内达到稳定，温度在土壤表层骤降，在土壤表层会形成薄冰层。这再现了低温条件下多孔介质lnapl、水迁移实验所得实验结果：低温环境中lnapl相迁移时高温区域面积较大，等温线在土壤表层较为密集，温度梯度变化较快。说明lnapl泄漏会对土壤的对流传热产生影响，对土壤起保温作用。lnapl流经薄冰层时，顺着薄冰层做层流流动。而当流体做层流流动时，在垂直于流动方向的热量传递，主要以热传导方式进行。lnapl流体的导热系数较小，传热热阻主要集中在该层中，温差也集中在该层。且达西数较小时，也会明显削弱对流传热。
[0145]
2)流场
[0146]
①
冻结条件下lnapl迁移特征
[0147]
冻结过程中lnapl饱和度变化如图4所示，可知，冻结前5h lnapl在重力作用下向下迁移，在含冰层下端有突变界面，lnapl含量在泄漏点的冻层处聚集，即冻层对lnapl的下渗有明显阻挡作用。而随着冻结作用增强，lnapl不再继续向下迁移。冻结5h后，其在冻结势的作用下克服重力势作用向上迁移，迁移程度较小，随后达到稳定。而对于包气带中的水分，其在冻结前1h内在随lnapl下渗产生的压力作用下发生向上运移，随后其饱和度随lnapl污染晕的移动而变化(图5)。含冰层处水分转化为冰，水饱和度相对较低。包气带中的气体则随lnapl的下渗逐渐向边缘扩散，直至整个区域被lnapl和水分填充。且其在含冰层下部出现模拟高值，发生聚集(图6)。
[0148]
冻结及非冻条件下lnapl运移5h和12h相比的量分布差异图以及冻结过程向下泄漏量分布剖面变化如图7所示，由冻结5h和冻结12h的量分布剖面图可知，lnapl向下泄漏，其量由泄漏表面至深部逐渐减少。而且lnapl从冻结5h至冻结12h量分布曲线上移，剖面上lnapl量均有不同程度减小，在冻结势作用下向地表方向迁移，暴露在表层。以冻结层浓度减小最明显，最高可达5mg/l。而在冻结层下部量减少较小，这是因为下部lnapl向上迁移，会受到冻层阻滞，形成累积，这与李兴柏等学者的冻结石油迁移实验结果一致(李兴柏和李国玉，2013)。故下部迁移的lnapl累积量抵消了部分lnapl的向上迁移量。由模拟结果计算得冻结条件下前5h lnapl向下迁移量为21.24ml，而冻结5h至冻结12h lnapl向上迁移至土壤表层的量为1.04ml。
[0149]
②
非冻结条件下lnapl迁移特征
[0150]
非冻结条件lnapl下渗过程中饱和度变化如图8所示。lnapl随时间推移逐渐向下渗透，前5h变化程度大，随后达到区域饱和稳定状态。非冻结条件下lnapl运移5h和12h相比的量分布差异图以及向下泄漏量分布剖面变化如图7所示，可知lnapl在非冻结条件下在重力势作用下向下渗，其量同样由泄漏表面至深部逐渐减少。在第5h至第12h，上部土壤水带和中间带lnapl浓度略有减少，lnapl向下迁移，在下部地下水饱和带以及毛细水带累积，而
且泄露源对上部土壤水带和中间带的补给量小于其向地下水饱和带及毛细水带的下渗量。由模拟结果计算得非冻结条件下前5h lnapl向下迁移量为38.59ml，至第12h达到41.30ml。
[0151]
非冻结条件与冻结条件下lnapl下渗模拟结果相比，因其不存在冻结层的阻挡作用，下渗速度快，下渗量远大于冻结时下渗量。12h后地下水面处体积分数可达0.040(图8)，而冻结条件下地下水面土壤处lnapl体积分数仅达0.017(图4)。非冻结条件lnapl向地下水面迁移，不存在向上迁移现象，而冻结条件下lnapl在冻结5h后在冻结势作用下发生向上迁移(图7)。冻结条件下lnapl相的下渗量相比非冻结条件下的量产生-51.1％相对误差。
[0152]
3)化学场
[0153]
①
冻结条件lnapl溶解挥发相迁移特征
[0154]
冻结条件下lnapl溶解相浓度分布如图9所示，随lnapl下渗，其在水中的溶解度逐渐增加，随后趋于稳定状态。冻结1h时，lnapl溶解相分布图像呈椭圆型，在地下水中浓度达116.5mg/l。迁移1h时，lnapl挥发相在包气带中浓度可达118.6mg/l。冻结稳定后，在lnapl以0.138cm/min的泄漏速度下，其在包气带中浓度可达186.7mg/l(lnapl在水中的溶解度为515mg/l)。而冻结条件下lnapl挥发相随气相逐渐向两边扩散(图10)，在土壤表层浓度可达101.5mg/l(lnapl在空气中的溶解度为139mg/l)。
[0155]
模拟可得冻结条件下不同状态lnapl含量关系图，如图11所示。可知冻结条件下冻结1h时总溶解量为0.304ml，挥发量为0.195ml。冻结稳定时lnapl总下渗量为23.74ml，溶于水的量为0.48ml，气相lnapl含量为0.165ml，而lnapl吸附量为2.71ml(图11a)。冻结条件下12h时可自由移动lnapl量(溶解相lnapl与自由相lnapl的和)为20.68ml，其含量冻结前5h逐渐增加，从冻结5h后的21.72ml逐渐减少，减少量为冻结势作用下向土壤表层迁移量。而吸附量随冻结时间增长增加到2.71ml后保持稳定(图11b)。
[0156]
②
非冻结条件lnapl溶解挥发相迁移特征
[0157]
非冻结条件下lnapl溶解相浓度分布如图12所示，水中的溶解相随lnapl下渗逐渐增加，随后稳定。迁移1h时，lnapl溶解相分布图像呈楔形，其在垂直方向扩散速度大于水平方向，在地下水饱和区浓度可达117.2mg/l。迁移1h时，lnapl挥发相在包气带中浓度可达143.2mg/l。下渗稳定后，在lnapl以0.138cm/min的泄漏速度下，lnapl溶解相在包气带中浓度可达185.8mg/l(lnapl在水中的溶解度为515mg/l)。而非冻结条件下lnapl挥发相随气相逐渐向两边扩散(图13)，在土壤表层浓度可达101.9mg/l(lnapl在空气中的溶解度为139mg/l)。
[0158]
非冻结下不同状态lnapl含量关系图如图14所示。可知非冻结条件下下渗1h时总溶解量为0.196ml，挥发量为0.251ml。稳定后lnapl总下渗量为41.30ml，溶于水的量为0.48ml，气相lnapl含量为0.168ml，lnapl吸附量为2.90ml(图14a)。吸附量占泄漏总量7％，吸附量大于溶解量，这与kim和corapcioglu的lnapl泄漏实验结果吻合(kim和corapcioglu,2003)。非冻结条件lnapl从泄漏点开始不断下渗，至第7h达到38.43ml稳定状态，吸附量随冻结时间增长增加到2.90ml后保持稳定(图14b)。
[0159]
对比lnapl于冻结、非冻结条件下1h时溶解相、挥发相浓度分布可知，lnapl溶解相在冻结条件下在水平方向上扩散速度大于非冻结条件，地下水饱和区溶解相浓度差为0.7mg/l，土壤非饱和带中lnapl挥发相最大浓度之差为24.6mg/l。在下渗过程中lnapl溶解量与挥发量冻结条件相比非冻结条件会产生2.4％、-24.9％的最大相对误差。稳定时lnapl
在水中及空气中的溶解扩散量浓度差别在1mg/l内，其受冰冻层影响小。又由不同状态lnapl含量变化图可知两种情景下lnapl溶于水的量相等，而气相lnapl含量差仅为0.003ml，相对误差为-1.8％。即含冰层对lnapl的溶解量和挥发量造成的影响很小，这是因为冻结条件下溶解量和挥发量主要通过分子扩散作用通过冰晶孔隙，含冰层的存在对分子扩散作用影响小。而对比吸附量，非冻结条件lnapl吸附量高于冻结条件lnapl吸附量，相对误差为-6.6％。这是因为一方面lnapl在非冻结条件下下渗量大于冻结条件，另一方面土壤表面被冰包裹会减少土壤的比表面积从而降低其吸附性能。