渗透系数的测定 *** 有哪些？
GB/T 1038-2000 塑料薄膜和薄片气体透过性试验 *** 压差法
ASTM D1434 测定塑料薄膜和薄片透气性能的测试 ***
ISO 2556 塑料.在大气压力下薄膜和薄板的气体透过率的测定.测压法
ISO 15105-1 塑料.薄膜和薄板.气体传输率的测定.第1部分压差法
JIS K7126 塑料薄膜及薄板的透气性试验 ***
YBB00082003气体透过量测定法
CNAS包装安全检测实验室兰光检测多种气体多种材质的渗透系数
渗透系数的测定 ***
渗透系数的测定 *** 主要分“实验室测定”和“野外现场测定“两大类。
1.实验室测定法
目前在实验室中测定渗透系数 k 的仪器种类和试验 *** 很多，但从试验原理上大体可分为”常水头法“和变水头法两种。
常水头试验法就是在整个试验过程中保持水头为一常数，从而水头差也为常数。 如图：
试验时，在透明塑料筒中装填截面为A，长度为L的饱和试样，打开水阀，使水自上而下流经试样，并自出水口处排出。待水头差△h和渗出流量Q稳定后，量测经过一定时间 t 内流经试样的水量V，则
V = Q*t = ν*A*t
根据达西定律，v = k*i，则
V = k*（△h/L）*A*t
从而得出
k = q*L / A*△h=Q*L /( A*△h)
常水头试验适用于测定透水性大的沙性土的渗透参数。粘性土由于渗透系数很小，渗透水量很少，用这种试验不易准确测定，须改用变水头试验。
变水头试验法就是试验过程中水头差一直随时间而变化，其装置如图：水从一根直立的带有刻度的玻璃管和U形管自下而上流经土样。试验时，将玻璃管充水至需要高度后，开动秒表，测记起始水头差△h1，经时间 t 后，再测记终了水头差△h2，通过建立瞬时达西定律，即可推出渗透系数 k 的表达式。
设试验过程中任意时刻 t 作用于两段的水头差为△h，经过时间dt后，管中水位下降dh，则dt时间内流入试样的水量为
dVe = －a dh
式中 a 为玻璃管断面积；右端的负号表示水量随△h的减少而增加。
根据达西定律，dt时间内流出试样的渗流量为：
dVo = k*i*A*dt = k*（△h/L）*A*dt
式中，A——试样断面积；L——试样长度。
根据水流连续原理， 应有dVe = dVo，即得到
k = （a*L/A*t）㏑（△h1/△h2）
或用常用对数表示，则上式可写为
k = 2.3*（a*L/A*t）lg（△h1/△h2）
2. 野外现场测定法
渗水试验(infiltration test)一般采用试坑渗水试验，是野外测定包气带松散层和岩层渗透系数的简易 *** 。试坑渗水试验常采用的是试坑法、单环法、和双环法。 是试坑底嵌入两个铁环，增加一个内环，形成同心环，外环直径可取0.5米, 内环直径可取0.25米。试验时往铁环内注水，用马利奥特瓶控制外环和内环的水柱都保持在同一高度上，（例如10厘米）。根据内环取的的资料按上述 *** 确定松散层、岩层的渗透系数值。由于内环中的水只产生垂直方向的渗入，排除了侧向渗流带的误差，因此，比试坑法和单环法精确度高。内外环之间渗入的水，主要是侧向散流及毛细管吸收，内环则是松散层和岩层在垂直方向的实际渗透。
当渗水试验进行到渗入水量趋于稳定时，可按下式精确计算渗透系数（考虑了毛细压力的附加影响）：K(渗透系数)= QL/ F(H+Z+L)。
式中：
Q-----稳定的渗入水量(立方厘米/分)；
F------试坑内环的渗水面积(平方厘米)；
Z-----试坑内环中的水厚度(厘米)；
H-----毛细管压力（一般等于岩土毛细上升高度的一半）(厘米)；
L-----试验结束时水的渗入深度（试验后开挖确定）(厘米)。
与饱和土相比，非饱和土的渗透系数有何特点？如何测试非饱和土的渗透系数
土壤水分渗透速率与入渗时间的关系，在入渗的初始阶段速率较快，到一定时间后速率趋于稳定，此时的渗透速率称为稳定入渗率，可用以表征土壤的渗透特性。
温度变化，水中有机物、无机物的成分和含量多少，均对渗透系数有影响。
在均质 含水层中,不同地点具有相同的渗透系数;在非均质含水层中，渗透系数与水流方向无关，而在各向异性含水层中，同一地点当水流方向不同时，具有不同的渗透系数值。一般说来，对于同一性质的 地下水饱和带中一定地点的渗透系数是常数；而非饱和带的渗透系数随岩土含水量而变，含水量减少时渗透系数急剧增大。
土石混合体渗透性能的试验研究
周中1 傅鹤林1 刘宝琛1 谭捍华2 龙万学2 罗强2
（1.中南大学土木建筑学院 湖南 长沙 410075
2.贵州省交通规划勘察设计研究院 贵州 贵阳 550001）
摘要 土石混合体属于典型的多孔介质，其渗透特性与砾石的百分含量关系密切。通过自制的常水头渗透仪，测定了不同含砾量时土石混合体渗透系数值，研究发现含砾量与土石混合体渗透系数之间存在指数关系；基于幂平均法，提出了土石混合体复合渗透系数的计算公式，并通过试验结果验证了该式的正确性，为土石混合体渗透系数的理论计算提供了一个简明有用的计算工具。
关键词 土石混合体 多孔介质 渗透性能 复合渗透系数 经验公式
土石混合体一般是由作为骨料的砾石或块石与作为充填料的粘土或砂组成，它是介于土体与岩体之间的一种特殊的地质体，是土和石块的介质耦合体［1］。因为土石混合体具有物质组成的复杂性、结构分布的不规则性以及试样的难以采集性等特殊的性质，从而给研究带来极大的困难，目前人们对于它的研究仍处于探索之中［2］。渗透与强度和变形特性都是土力学中所要研究的主要力学性质，其在土木工程的各个领域都有重要的作用［3］。土石混合体属于典型的非均质多孔介质［4］，其渗透系数是由高渗透性的砾石和低渗透性的土体复合而成的。土的渗透系数可以通过室内试验由达西定理计算得出，然而土石混合体的渗透系数却难以确定，主要原因是：取样困难；难以进行常规的渗透试验；大尺度的渗透试验不仅造价高、准确性差，而且试验结果离散度大，难以掌握其规律性。因此能够求出土石混合体复合渗透系数的计算公式具有重要的理论意义和工程应用价值。
土石混合体中土与砾石粒径的界限值为5mm，即将粒径小于5mm的颗粒称为土、大于5mm的颗粒称为石，砾石含量用P5表示［1］。利用自制的常水头渗透仪，研究砾石体积百分含量P5从0%逐步过渡到100%（间隔10%）时土石混合体的渗透系数，每种配比作平行试验3次，共33次渗透试验。
1 土石混合体渗透性能试验
1.1 试样的基本物理力学性质
试验所取土样为正在修建的上瑞高速公路贵州段晴隆隧道出口处典型性土石混合体，其天然状态土的物理指标及颗粒级配曲线见表1和图1。由图1可知现场取回土样的不均匀系数Cu为12.31，说明土样中包含的粒径级数较多，粗细粒径之间差别较大，颗粒级配曲线的曲率系数Cc为1.59，级配优良。
表1 天然状态土的基本物理指标
图1 天然状态土的颗粒级配曲线
1.2 大型渗透仪的研制
《土工试验规程》（SL237—1999）规定粗粒土的室内渗透系数需由常水头渗透仪测试，国内常用的常水头渗透仪是70型渗透仪。70型渗透仪的筒身内径为9.44cm，试验材料的更大粒径为2cm，规范［5］要求筒身内径应为更大粒径的8～10倍，因此70型渗透仪的筒身内径过小，有必要研制大尺寸的渗透仪。自制渗透仪的内径和试样高度至少应为更大颗粒粒径的8倍，即至少应为16cm，另外，考虑到边界效应，试样的上下两头分别增加2cm，因此，自制渗透仪的内径和试样高分别取为16cm和20cm。考虑到土石混合体的渗透性较强，选取进排水管的口径为2cm。自制的大型常水头渗透仪如图2和图3所示。
图2 自行研制的渗透仪
图3 常水头渗透仪示意图
数据单位为cm
1.3 试验步骤
首先，将由现场取回的土样烘干、过筛，并根据粒径的大小分为0～5 mm的土和5～20mm的砾石两部分。然后，按照试验要求的砾石体积百分含量P5，以10%的初始含水量配制试样，静置24 h。试验时，将配制好的试样分层装入圆桶中，每层装料厚度30mm左右，分层压实，记录每层的击实数。按上述步骤逐层装样，至试样顶部高出测压孔约3cm为止。测出装样高度，准确至0.1cm。在试样顶部铺一层2cm厚的细砾石作缓冲层。之后，由进水管注入蒸馏水，直至出水孔有水流出，静置24 h使试样充分饱和。用量筒从渗透水出口测定渗透量，同时用温度计测量水温，用秒表测记经一定时间的渗水量，共测读6次，取其平均值，6次结果相差不得超过7%，否则需重新测定。
1.4 试验数据
按照试验设计的各种砾石体积百分含量P5共需作11组试验，每组试验作平行试验三次，取3次测量的平均值，并乘以温度校正系数
，即可求出每组试验20℃时的渗透系数，渗透系数的测量结果见表2。
表2 渗透系数测定结果
2 试验结果分析
2.1 渗透系数与砾石含量的关系
不同含砾量的颗粒级配曲线如图4所示，由图4可以求出各曲线的粒径特征系数及不均匀系数Cu和曲率系数Cc。
图4 试样的颗粒级配曲线
图5为土石混合体砾石含量P5与20℃时渗透系数的关系曲线。从图5可以看出，随着含砾量的增加，渗透系数急剧增加，可见，在设计中可以通过调节砾石的含量来控制土石混合体的宏观渗透性能。
图5 粗粒含量与渗透系数的关系
从图5还可以发现，土石混合体中砾石的含量P5与渗透系数k之间存在指数关系，与文献［6］的研究成果相似，即
土石混合体
式中：k0为P5＝0时土的初始渗透系数；n为与土石混合体本身性质相关的常数。对于文中试验值，k0与n分别为0.0006cm/s和8.82。在工程中可以通过少量试验来确定k0，n值，以此来预测不同级配土石混合体的渗透性。
2.2 土石混合体的复合渗透系数
近几十年来，许多学者在揭示影响和决定土的渗透系数内在因素及其相互关系方面进行了大量工作，并取得了有益的成果［7～12］，被认为依然有效且目前常用的确定渗透系数的半经验、半理论公式有：
（1）水利水电科学研究院公式［7］：
土石混合体
式中：k10，k20分别为温度为10℃和20℃时的渗透系数（cm/s）；η10/η20为温度为10℃和20℃的粘滞系数比；n为孔隙率；d20为等效粒径（mm）。
（2）泰勒（Taylor）［9］用毛管流的哈根-伯努力（Hange-Poiseuille）方程导出渗透系数的表达式：
土石混合体
式中：ds为当量圆球直径，可以用等效粒径d20代替；γw为液体容重；μ为液体粘滞度；e为孔隙比；C为形状系数，通常取C＝0.2。
式（2）和式（3）均是针对土体的渗透特性提出的半经验、半理论公式，然而对于非均质性更强、粒径差别更大的土石混合体来说，其适用性不是很强。土石混合体中砾石形成骨架，细颗粒充填孔隙，其渗透系数是由低渗透介质土体的渗透系数kS和高渗透性介质砾石的渗透系数kG复合而成。土石混合体复合渗透系数不是按体积百分含量的简单复合，而是高低渗透性介质的耦合。在参考相关文献［10～12］的基础上，基于幂平均法，本文提出的土石混合体复合渗透系数k复合的表达式为
土石混合体
式中：P5为砾石的体积百分含量，%；kG为砾石的渗透系数，cm/s；kS为土的渗透系数，cm/s；f为系数。
砾石的体积百分含量P5可以由筛分法求出；土的渗透系数kS和砾石的渗透系数kG可以由室内试验直接求出或参考相关资料确定；系数f可以通过少量试验回归分析确定，因此可以说（4）式是一个简明实用的土石混合体复合渗透系数计算公式。
图6 不同计算 *** 结果比较
为进一步验证（4）式，我们将试验测得的k值与用（2），（3），（4）式计算得到的k值进行对比分析。结果见图6，具体数值见表3。由图6和表3可知据水利水电科学研究院公式和泰勒公式计算结果均高于实测值，尤其是当P5≤30%时，（2）式计算结果和（3）式计算结果比实测值大2～3个数量级，与实测值相差较大。而用本文 *** 得到的土石混合体的渗透系数最接近实测值，平均相对误差仅为0.6%，能够作为土石混合体渗透系数定量预测的有效工具。在工程设计中，可以根据工程对土石混合体渗透性的要求，依据本文提供的经验公式，调整土石混合体中砾石的含量，达到控制土石混合体渗透能力的目的。
表3 土石混合体渗透系数及相关参数
3 结论
（1）利用自制的常水头渗透仪，测定了不同含砾量时土石混合体的渗透系数值，并指出含砾量与土石混合体渗透系数之间存在指数关系。在工程设计中可以通过合理调整土石混合体中砾石的含量，达到控制其渗透性能的目的。
（2）指出土石混合体的渗透系数是一种由高渗透性的砾石和低渗透性的土体复合而成的，给出了土石混合体复合渗透系数的计算公式，并通过试验结果验证了计算公式的正确性，为土石混合体渗透系数的定量预测提供了一个简明有用的计算工具。
参考文献
［1］油新华.土石混合体随机结构模型及其应用研究.北方交通大学博士论文.2001：1～18
［2］油新华，汤劲松.土石混合体野外水平推剪试验研究.岩石力学与工程学报.2002，21（10）：1537～1540，60～129
［3］黄文熙.土的工程性质.北京：水利电力出版社.1984：60～129
［4］薛定谔 A E.多孔介质中的渗流物理.北京：石油工业出版社.1984：141～173
［5］中华人民共和国水利部.土工试验规程（SL237—1999）.北京：中国水利水电出版社，1999：114～120
［6］邱贤德，阎宗岭，刘立等.堆石体粒径特征对其渗透性的影响.岩土力学，2004，25（6）：950～954
［7］刘杰.土的渗透稳定与渗流控制.北京：水利电力出版社，1992：1～20
［8］ Wen X H，Gomez-Hernandez J J.Upscaling hydraulic conductivities in heterogeneous media：An overview.Journal of Hydrology，1996，183：ix～xxxii
［9］ Taylor D W.Fundamentals of soil mechanics.John Wiley ＆ SONS，Inc.，1948
［10］ Brown W F.Solid mixture permitivities.Journal of Chemical Physic，1955，23（8）：1514～1517
［11］ Dagan G.Analysis of flow through heterogeneous random aquifers by the method of embedding matrix—1：Steady flow.Water Resources Research，1981，17（1）：107～122
［12］ Noetinger B.The effective permeability of a heterogeneous porous medium.Transport in Porous Media，1994，15：99～127
数值模拟流程
不同的软件进行数值模拟时所需的参数、计算 *** 、剖分格式等不尽相同，数值模拟的过程也不同，但大致相同，本文以TOUGHREACT为例介绍CO2地质储存数值模拟的流程。
（一）研究范围的确定
一般情况下，独立的天然水文地质系统是计算区更好的选择，它具有自然边界，便于较准确地利用其真实的边界条件，避免人为边界在资料提供上的困难和误差。但是在实际工作中，常常不能完全利用自然边界，这就需要充分利用勘察和长期观测资料等建立人为边界。在确定计算区域时，除了保证范围足够大以外，还应使假定的边界条件尽可能接近真实状态。
计算范围的划定应充分考虑研究目的、区域地质构造、储层岩性、储层岩石矿物组成及地下水化学成分等多方面因素。数值模拟时间根据研究目的不同具有不同的时间尺度。就CO2地质储存数值模拟而言，如果不考虑地球化学作用，封存系统在1000年数量级的模拟时间内基本上已达到平衡或稳定。在划定边界时还应考虑CO2在储层中的扩散距离，与研究区地质模型的孔隙度，渗透率等参数关系密切。为了保证所选模型范围边界在模拟期内不影响模拟结果，尽量通过具有相同地质条件的天然CO2气田（藏）进行类比，确定大体的计算范围的边界。如果考虑地球化学反应，由于CO2注入引发的水－岩－气反应对围岩岩性改变较显著，制约着CO2注人的速度和径向运移的距离等。
（二）明确研究目的
在进行数值模拟以前首先要明确利用数值模拟技术要解决什么样的问题。对于CO2地质储存工程而言，进行数值模拟的目的主要是在CO2地质储存工程实施前，通过数值模拟技术对工程的选址、方案设计进行优化，工程实施期技术指导、运行期监测及后期CO2泄漏的风险评估等进行预测，以指导项目科学、合理地实施，将CO2泄漏风险降至更低。
研究目的决定着前期资料的收集类型、地质建模的侧重点、地质模型离散的精密度以及初始、边界条件的处理方式等过程。
（三）资料的收集整理
1）通过遥感、综合地质调查、物探、钻探和各类样品测试分析等手段获取场地深部地层岩性、地质构造、水文地质、水文地球化学、岩石矿物资料和数据；
2）搜集和分析CO2地质储存场地地质岩性、区域构造格架、活动断层与地震活动情况等；
3）采用钻井岩心、测井和地震反射 *** ，调查CO2地质储存场地目标储层和盖层的空间分布形态，埋深、厚度和规模等；
4）使用X射线衍射、扫描电镜等 *** 研究分析封存场地岩石矿物组成、孔隙结构特征及其物理化学性质；
5）通过采取浅部、深部含水层水样进行水质全分析，获得储盖层地层水及浅部含水层初始水化学成分。
不同的数值模拟软件其数学模型的数值解法不同，空间离散方式也不尽相同，所需的模型参数也有一定的差异，表9-1即为TOUGHREACT数值模拟所需要的主要参数。
表9-1 CO2地质储存模拟过程中需要的主要参数（以TOUGHREACT为例）
（四）概念模型的建立
利用所掌握研究区的地形地貌、气象水文、地层岩性、地质构造、水文地质、水文地球化学等资料，根据研究目的初步建立地质模型雏形。为了解决实际问题，往往要对所建模型及初始条件、边界条件等进行必须要的概化，概化内容包括：
1）含水层系统结构的概化：应依据含水层的类型、岩性、厚度、渗透系数等，将内部结构概化为均质、非均质各向同性或者异性的含水层；
2）侧向边界和顶，底部边界的概化；根据研究区边界的划定，将侧向边界和顶、底部边界概化为一、二和三类边界条件；
3）根据所要解决问题的难易程度，将地质模型概化为一维、二维或者三维模型，并进行合理精密级的网格剖分；
4）源汇项的处理：储存库一般位于地表800m以下，上覆有致密的盖层岩石，CO2难以通过越流方式得到补充和排泄，在理想的条件下整个地质储库系统的CO2通过径向对流或抽取井的方式得到补充和排泄。
（五）模拟方案的设计
针对不同的研究问题，可以根据设计者技术水平、场地经验等设计出不同的方案。利用数值模拟技术对多种方案进行模拟分析，评估方案实施的可行性，并对方案进行优化，最终得出经济、合理的方案。如进行CO2灌注场地的遴选时，要对几个目标场区进行灌注能力，储存潜力、CO2在储层中的扩散速度及距离以及后期CO2泄漏风险等环节进行模拟，然后对模拟输出的数据进行后期处理。通过分析、统计这些数据确定更佳的灌注场地，同时合理划分储盖层的层位、厚度等。对场区灌注能力和储存潜力评估是数值模拟技术所能解决的CO2地质储层领域的又一难题。不同的灌注方式导致CO2进入储层的速度、流量等不同。根据设计者所设计的所有可行的注入方式进行模拟，分析比较不同方案下注入能力和注入量的大小，确定更优的灌注方式。
模拟方案的设计要根据研究问题不同视情况而定，而且方案设计的合理与否依赖于设计者自身理论和实际场地经验，不同的研究着对于同一个问题可能设计出不同的方案。我们可以针对这些方案建立不同的模型，通过模拟技术判断其合理性与可行性，最后确定更佳的方案。
（六）数值模型及模拟软件的选择
数值模拟的关键是地质模型概化、计算精度和计算速度。由于计算的精度取决于离散的程度，而离散的程度又决定了计算的速度，这是一对矛盾，要根据解决问题的需要来选择离散化的程度和计算速度。
CO2在储层中的运移、溶解以及与围岩的化学反应形成了一个多相、多组分的反应体系，涉及的主要数学方程有超临界CO2－水的两相流体运动控制方程、溶质运移控制方程和化学反应方程等。在建立数值模型时，通常采用的 *** 有有限差分法、有限元法和积分有限差法等。
由于实际应用时多采用已有的数值模拟软件对CO2地质储存的全过程进行模拟.不涉及软件的开发及程序代码的编写，只需根据研究的需要选择合适的软件进行模拟预测，而软件一旦选定，数学模型和数值模型基本上已经确定。以TOUGHREACT为例，基于上述所建的概念模型建立数学模型。气相和液相统一的偏微分方程为（9-1），咸水的偏微分方程为（9-2），超临界方程为（9-3）。方程中所涉及的相关字符和角标见表9-2。
二氧化碳地质储存技术 *** 概论
二氧化碳地质储存技术 *** 概论
表9-2 数学模型中所涉及的符号的意义
（七）数值模型的建立
1.网格剖分
建立了地质模型以后要对研究区离散化，即进行网格的剖分。首先确定离散点，即把所研究的区域按照某种几何形状（如矩形、任意多边形等）分割成网格系统。研究区的边界可以用最接近它的格线近似表示。当网格划分得足够小时，曲折的格线也能够很好地刻画出边界的形状。此过程又称研究区域的离散化（剖分）。在离散化时要遵循以下两条基本原则。
1）几何相似。要求物理模拟模型从几何形状方面接近真实被模拟体。
2）物理相似。要求离散单元的特性从物理性质方面（含水层结构、水流状态）近似于真实结构在这个区域的物理性质。
网格类型大体上可分为规则和不规则两种。规则网格剖分包括矩形、三角形和其他规则图形的剖分（图9-3），而不规则网格则包括不规则多边形等。网格的形状主要根据研究区形状而定。
图9-3 网格剖分
网格剖分对计算的精度及计算的效率有很重要的影响。精度越高对模拟结果刻画的越精细，但是数据的计算量越大，对计算机的要求也越高。建议在进行地质模型剖分时先采用较粗的网格剖分，如果这种剖分方式下模拟结果合理然后再进行精细化剖分，用于对模拟结果更加详细的刻画。
2.参数和初始条件
初始条件是指在初始时刻（t=0）时研究区内求解数学模型主要状态变量的初始值。选择的应用软件不同所需的状态变量数量、种类不同。如TOUGHREACT所需的初始主要状态变量包括压力、温度和组分浓度的空间分布。地质参数包括孔隙度、渗透率、密度、压力、温度、毛细压力等参数值。这些数值一部分采用室内实验测得，另一部分采用参考文献的经验值；地层水的化学成分的初始值采用实际地层水的化学分析，主要是8大离子的浓度、盐度和pH 等。如果研究区深部地层中的水样难以获得，如盖层，则采用静态平衡的 *** ，利用具有与储层相同盐度的咸水与含有原生矿物的地层岩石在原地层环境下进行化学反应，获取平衡状态下的地层水化学成分的初始值；通过岩矿分析、电子扫描、Ⅹ衍射等手段，获得组成CO2地质储层盖层的原生矿物成分体积含量初始值，并根据原生矿物的组成合理判断次生矿物。
从原则上讲，初始时刻是可以任意取定的，只要该时刻所需的参数和状态变量值已知即可。因此我们不应该把初始条件理解为研究系统的初始状态。具体如何取，应该视问题的需要、资料来源、计算方便与否等因素而定。
3.边界条件
边界条件是某一实际问题数学模型具有定解的必要条件之一。地下水流问题和溶质运移问题边界条件的定义不尽相同，但一般概化为以下三种。
（1）一类边界条件（Dirichlet条件）
解决水流问题时，此类边界条件为在边界上所有点的水头是给定的；对于溶质运移问题，一类边界条件是指研究区边界上的溶质浓度分布已知。解决CO2—水两相流动问题时，此类边界条件为在边界上所有点的压力是给定的。
（2）二类边界条件（Neumann条件）
当已知某一边界的单位面积流入或流出的流量时，可视作解决流动问题的二类边界；相对溶质运移来讲，此类边界又称给定弥散通量边界，即边界上的弥散通量随时间变化规律已知。
（3）三类边界（Cauchy条件）
当研究区一部分满足一类Dirichlet条件，而另一部分满足二类Neumann条件时，这类问题称为混合边界问题，称为三类边界。对溶质运移而言，此类边界为边界上溶质通量随时间变化规律已知。
在CO2地质储存数值模拟过程中，由于储层地层多在800m以下，地质模型的顶部和底部根据实际需要可以处理为不透边界；为了避免边界对模拟结果的影响，研究区的范围一般比实际CO2所能运移到的范围大得多，因此，在处理四周边界时一般设置为无穷一类边界或不透边界。在确定边界条件时，应根据水文地质条件以及现有的资料来综合考虑。
4.源汇项处理
在多孔介质中流动和溶质运移的问题中，对流、水动力弥散和溶质源或／和汇，是决定含水层中任一内点上溶质质量时变率的两大因素。源汇项问题在水质与水量计算中以及正确处理对流－弥散方程和渗流基本微分方程中占有重要地位。作为源汇项的方式很多，如越流补给、含水层弹性释放补给以及抽（注）井的补给等。
对于深部咸水层CO2地质储存系统而言，系统的顶部一般为具有低渗、低孔的泥岩、页岩等致密性岩层，越流补给较难发生。整个CO2地质储库系统的源汇项主要指对流（如侧向边界）和抽（注）井。
（八）模型的校正与验证
模型识别是建立地下流体数值模型最重要的环节之一，正确理解和进行拟合对于提高数值模型的仿真性是至关重要的。在有实测结果的情况下如示范工程，可将模拟结果与实测结果进行比较，对相关参数进行适当合理的调整，使模拟结果在给定的误差范围内与实测结果吻合。若误差较大，应该重新检验概念模型的可靠性，甚至重新建立概念模型。在识别校正以后，应采用校正好的模型继续计算，并与未用来识别校正的实际数据比较，验证模型的准确性和可靠性。若存在较大误差，需重复前面的过程。在没有实测结果的情况下，数值模型的可靠性可通过类比相关资料或根据个人经验和理论判断。
（九）模拟预测
模型预测是实施数值模拟技术的主要目的。对于CO2地质储存工程而言，由于CO2地质储存技术的提出为时尚短，针对CO2在深部咸水层中的运移、扩散、与地层水和围岩产生的化学反应，以及由于CO2灌注引起的储盖层物理、化学性质变化研究均处于研究和发展阶段。因此，在工程实施过程中急需具有技术指导性的工具产生，避免造成投资的浪费和CO2泄漏等风险的出现。
利用经过识别校正与验证过的数值模型对CO2地质储存过程进行模拟预测，有针对性地对模拟数据进行后期处理，如统计分析、比较等手段对结果进行解译，以此达到场地的优选，目标储层灌注能力、储存潜力的评估，CO：扩散运移途径和速度、不同捕集方式封存量及它们之间的时空转化等过程的详细刻画与模拟仿真等目的。同时可以预测CO2在已有、重新激活或新生成的裂隙中逃逸的可能性及时间、CO2泄漏风险评估以及评价CO2泄漏对浅层地下水的水质、水量及对地表环境的影响等。
上述结果的分析只是数值模拟技术所能解决问题的冰山一角。对于数值模拟结果的处理要根据所研究的目的进行有针对性的提取和解译。通过对处理后的数据进行总结分析，发现问题从而解决问题，并掌握内在规律，为CO2地质储存工程的前期设计、工程实施、中期监测管理提供理论支持和科学的技术指导，并可以提前开展风险预测，尽早制定预案防范CO2地质储存工程实施及运行过程中可能出现的隐患。
0条大神的评论