多孔介质中非水相流体运移的数值模拟

针对多孔介质中水、气和非水相流体(NAPLs)的多相流动特点,建立了非水相流体(NAPLs)污染物迁移模型,分析了非水相流体在土壤非饱和区和地下水系统中的运移规律。通过有限元数值解对轻非水相流体和重非水相流体在土壤系统中的迁移过程进行模拟,得到了污染物的时空分布特征和污染范围。计算结果表明,数值模拟方法能够合理地描述非水相流体的运移过程和污染特征。土体渗透性和污染物残余饱和度是其重要影响因素。     

非均质裂隙介质中重非水相流体运移

由于裂隙介质具有强烈非均质性,使得重非水相流体(DNAPLs)在裂隙介质中的运移行为较孔隙介质更为复杂.基于指示模拟算法模拟裂隙介质的二元变量(粗糙面接触和裂隙开口),综合模拟退火算法构建裂隙渗透率随机场.采用T2VOC模拟DNAPLs在裂隙介质的运移,探讨粗糙面接触的相关长度、各向异性比、倾角以及非均质性程度对DNAPLs运移分布的影响.数值分析结果表明,粗糙面接触的相关性越差,DNAPLs污染范围越大;粗糙面接触的各向异性比和倾角增大,将导致更多DNAPLs残留;而相关长度和各向异性比增大以及倾角减小,都会引起DNAPLs锋面运移速率增大;随着渗透率非均质性增强,会导致局部粗糙面接触上蓄积的DNAPLs饱和度增大,运移速率减小.     

土壤水分运动空间变异性尺度效应的染色示踪入渗试验研究

探讨了土壤水非均匀流动特性和描述方法,通过染色示踪剂调查了三种试验尺度条件下非均匀流动模式,并采用随机层叠模型对不同实验尺度条件下非均匀流动模式进行了模拟。随机层叠模型中具有对数正态分布性质的随机层叠发生器被用来描述水流入渗过程,不同的方法被用于模型参数求解。试验观测和模拟计算结果均表明,尺度特性是非均匀流动的重要影响因素之一,准确的描述不同研究尺度下的非均匀流动特征,须同时考虑流动在水平和垂直方向的变异性。随着研究尺度的增加,流动的非均匀性变异程度更加明显。     

石羊河流域平原区土壤入渗特性空间变异的研究

在野外原位入渗试验的基础上,建立了简化Philip入渗模型及其土壤转换函数,借助ARCVIEW软件生成了石羊河流域平原区土壤入渗特性参数的空间分布图,并分析了入渗特性的空间变异特征。研究表明,简化Philip入渗模型是点面转化中比较理想的模型,模型参数α的土壤转换函数可以简便地估算土壤入渗特性,两者的结合是研究大尺度土壤入渗空间变异的一种有效方法。     

非均质介质中重非水相污染物运移受泄漏速率影响数值分析

重非水相污染物(DNAPL)在地下介质中运移和分布受多种因素控制,包括DNAPL本身的物理化学性质,土的性质,泄漏条件等等。由于介质的非均质性,使得多相流运移行为更为复杂。基于地下水随机理论构建渗透率随机场,采用蒙特卡罗方法探讨泄漏速率对非均质饱和介质中DNAPL运移的影响。数值结果表明,在泄漏总量一定的情况下,泄漏速率越低,介质非均质性对DNAPL运移的影响程度越高。反之,DNAPL的渗漏速率越高,小尺度地层的非均质性影响越低。由于DNAPL运移过程中在垂直方向受重力的影响,污染羽在空间上的质心位置(一阶矩)以及展布范围(二阶矩)在垂直方向上的变异程度要高于水平方向。     

注入丁醇调节重非水液相密度的微空隙试验模拟

密度大于水的重非水液相(DNAPLs)有机污染物在重力作用下向地下介质深部迁移从而增加污染范围.前人通过一维砂柱和二维砂箱试验发现利用密度调节技术可降低DNAPLs向下迁移的风险,但目前缺乏微观尺度上密度调节影响DNAPLs迁移的定量观测.本研究试验模拟丁醇注入微空隙调节四氯乙烯(PCE)的密度,通过建立非水相中染色P...     

水—气二相流及污染物运移数值模拟研究

本文通过对近些年来水－气二相流数值模拟进展情况的分析,总结了在水－气二相流模拟过程中,各个主要参数或物理过程的数学概化方法,以及近年来在求解水气二相流方程和污染运移方程的方法上的改进情况。并在此基础上提出了需进一步研究解决的问题。     

塔里木盆地构造变形迁移与三大地质事件的关系及其地质意义

塔里木盆地经历多期构造运动,形成多期断裂系统和隆凹构造,在地质历史的演化过程中,均发生了横向迁移。通过大量二维地震资料解释和断裂分析,结合地层残余厚度及分布状况,参考重点层系构造图,经研究发现,不同地质时期的断裂系统具有按照先后顺序有规律的横向迁移特征。同时,不同时期塔里木盆地隆凹格局也发生了巨大变化。随着塔里木盆地独有的早期深埋、中期二叠纪火山活动及晚期快速沉降三大地质事件的发生,地层沉积沉降和构造隆升的不均衡性,使得烃源岩在平面上的演化程度呈现不均一性。按现今构造单元可分为古隆起生烃终止区、古隆起阶段生烃区、凹陷生烃衰竭区、凹陷阶段生烃区和凹陷周期持续生烃区等5类烃源岩演化区。现今台盆区仅有阿满过渡带、塔北隆起中段、草湖次凹、麦盖提斜坡及塔西南山前等5个地区发育有效烃源岩,持续生烃或阶段生烃,为塔里木盆地台盆区晚期成藏提供了不竭的烃源。

     

土壤入渗特性的空间变异性及土壤转换函数

以在杨陵一级阶地上进行的入渗试验为基础,利用多重分形和联合多重分形方法,研究分析了土壤入渗特性在多尺度上的空间变异性及与土壤物理特性的相关性,并在此基础上建立了土壤入渗特性的土壤转换函数。研究表明:土壤入渗特性具有多重分形特征,观测尺度上,稳定入渗率、前30 min累积入渗量的空间变异都主要由粗粉粒和粘粒含量的空间变异造成。在多尺度上,稳定入渗率的空间变异受土壤容重、粗粉粒和粘粒含量的空间分布影响显著,前30 min累积入渗量的空间变异受粗粉粒和粘粒含量的空间分布影响显著,与观测尺度上的影响因素相同。基于联合多重分形分析建立的稳定入渗率的土壤转换函数的计算误差较小,前30 min累积入渗量的土壤转换函数的计算误差较大。     

非饱和冻土水汽迁移与相变过程的光滑粒子法模拟

为了研究多年冻土表层的水热分布情况,在非饱和冻土的能量守恒方程和水分迁移的质量控制方程的基础上考虑冰水相变和水汽相变过程,并考虑水汽运移传热及温度势对水汽迁移的影响,建立了非饱和冻土的水-热-汽耦合模型。采用光滑粒子流体动力学(smoothed particle hydrodynamics,简称SPH)方法可方便地计算它们的演化过程。为此,在计算中先求解能量守恒方程的含冰量及气态水含量,再对未冻水含量和温度场进行求解,从而实现了温度场与水汽场的耦合。在此基础上,模拟计算了第1类热边界条件下半无限空间介质内非稳态温度场、体积含水率及水汽通量的分布情况,并将计算结果与未考虑耦合的解析解进行比较,结果显示水汽耦合的作用不容忽略。最后,针对处于季节性周期温度边界下路基的水热场的分布情况进行计算。研究表明,相比于水-热耦合模型,所建立的水-热-汽耦合模型得到的计算结果更为接近实际监测结果,可很好地揭示非饱和冻土中的水热汽迁移特征及其相变过程。     

The role of water in opal

The occurrence of water in natural opaline silicas has been studied by differential thermal analysis, thermogravimetric analysis, infra‐red analysis and nuclear magnetic resonance. The results show that in the “crystalline” opals studied, some 90 per cent or more of the total water is physically adsorbed whereas in “amorphous” opals, at least 20 per cent but perhaps much more of the total water is held as hydroxyl groups chemically bonded to the silica surface. The rate of water loss on heating is also different, being chiefly controlled by the pore structure in “crystalline” opals but to a significant extent by the surface structure in “amorphous”.     

Three-phase numerical model of water migration in partially frozen geological media: model formulation,validation, and applications

Water in the subsurface of the Earth’s cold regions—and possibly the subsurface of Mars—resides in the liquid, vapor, and ice phases. However, relatively few simulations addressing full three-phase, nonisothermal water dynamics in below-freezing porous media have been undertaken. This paper presents a nonisothermal, three-phase approach to modeling water migration in partially frozen porous media. Conservation equations for water (as ice, liquid, and vapor) and a single gas species (in the gas phase and dissolved in water) are coupled to a heat transport equation and solved by a finite-volume method with fully implicit time stepping. Particular attention is given to the method of spatial differencing when the pore space is partially filled with ice. The numerical model is able to reproduce freezing-induced water redistribution observed in laboratory experiments. Simulations of Earth permafrost dynamics and of the formation and evolution of a planetary-scale cryosphere on Mars demonstrate the new capabilities.     

Pore-scale removal mechanisms of residual light non-aqueous phase liquids in porous media

Glass-etched micromodel and visualization technology were used to research the removal mechanisms of residual light non-aqueous phase liquids (LNAPLs) droplets at pore scale. In this research, n-hexadecane was selected as the LNAPL model substance. During the injection of de-ionized water into the pore channel, residual n-hexadecane droplet was removed quickly at first, then gradually became slow and finally remained nearly invariable; the droplet size after a definite time depends on its initial value. The relationship between total volume of residual n-hexadecane and time can be described by a logarithmic equation: V = V ₀ − 0.0152lnt + 0.0360 (0 < t < 9,420, R ² = 0.9621). Based on percolation theory, the relationship between fluid percolation velocity in pores and time was determined through the analysis of residual n-hexadecane droplet size and the flowing characteristics of injected fluid. Gaussian model was used to fit the cumulated dissolved n-hexadecane, the correlation coefficient R ² = 0.97269. The removal process of n-hexadecane in network model experiment is similar to that of the sand column flushing experiment, which indicates that the mass transfer micro-mechanism of network model experiment could be used to interpret the results of column flushing experiment.     

土壤入渗特性的空间变异规律及其变异源

以封丘地区典型的4种土地利用类型(传统耕作地、免耕地、金银花地和杨树林地)为例,系统研究了不同土地利用方式下土壤入渗特性的空间变异强度、空间变异结构及其主要变异源。结果表明,频繁的耕作会降低土壤异质性,弱化土壤特性空间变异结构,扰动较小的免耕地和金银花地的土壤入渗特性则具有相对明显的空间变异结构,扰动最小的杨树林地土壤的入渗特性具有最为明显的空间变异结构和最小的变程。容重、有机质含量和初始含水量是耕作地入渗的最主要变异源,初始含水量则是金银花地入渗的唯一显著性变异源;杨树林地入渗变异的主要变异源是有机质和初始含水量;影响免耕地入渗特性变异的变异源并不显著。在有些情况下,水温也可能成为入渗特性空间变异研究的重要干扰因素。     

地表水入渗对垃圾填埋场水质水量影响的数值模拟分析

将填埋垃圾视作一种特殊土体,建立了综合描述垃圾填埋场内渗滤液水分的饱和-非饱和渗流、垃圾土体沉降变形、有机污染物的释放、迁移、转化,以及微生物的生长、衰减等物理化学过程的多场耦合仿真分析模型,将任意的Lagrangian- Eulerian（ALE）方法引入到模型求解中,基于自行开发的仿真分析程序,开展了地表水入渗对垃圾填埋场水质、水量影响的数值仿真分析研究。模拟结果表明：300 d的非饱和水分入渗及内源水产生可使填埋单元内渗滤液饱和液面最高达6 m;地表水入渗对水相可溶性可降解有机物浓度和厌氧微生物浓度具有显著的影响,并间接地影响固相可溶性可降解有机物和微生物的浓度分布。     

The role of eclogite in the mantle heterogeneity at Cape Verde

The Cape Verde hotspot, like many other Ocean Island Basalt provinces, demonstrates isotopic heterogeneity on a 100–200 km scale. The heterogeneity is represented by the appearance of an EM1-like component at several of the southern islands and with a HIMU-like component present throughout the archipelago. Where the EM1-like component is absent, a local DMM-like component replaces the EM1-like component. Various source lithologies, including peridotite, pyroxenite and eclogite have been suggested to contribute to generation of these heterogeneities; however, attempts to quantify such contributions have been limited. We apply the minor elements in olivine approach (Sobolev et al. in Nature 434:590–597, 2005; Science, doi: 10.1126/science.1138113, 2007), to determine and quantify the contributions of peridotite, pyroxenite and eclogite melts to the mantle heterogeneity observed at Cape Verde. Cores of olivine phenocrysts of the Cape Verde volcanics have low Mn/FeO and low Ni*FeO/MgO that deviate from the negative trend of the global array. The global array is defined by mixing between peridotite and pyroxenite, whereas the Cape Verde volcanics indicate contribution of an additional eclogite source. Eclogite melts escape reaction with peridotite either by efficient extraction in an area of poor mantle flow or by reaction of eclogite melts with peridotite, whereby an abundance of eclogite can seal off the melt from further reaction. Temporal trends of decreasing Mn/FeO indicate that the supply of eclogite melts is increasing. Modelling suggests the local DMM-like end-member is formed from a relatively peridotite-rich melt, while the EM1-like end-member has a closer affinity to a mixed peridotite–pyroxenite–eclogite melt. Notably the HIMU-like component ranges from pyroxenite–peridotite-rich melt to one with up to 77 % eclogite melt as a function of time, implying that sealing of melt pathways is becoming more effective.     

The role of geochemical self-organization in the migration and trapping of hydrocarbons

Many sedimentary basins of the world contain abnormally pressured “compartments” surrounded by impermeable seals, which may serve as hydrocarbon traps. Observations of the many types of seals that may exist point to an origin due to self-organization. Self-organization is the development of repetitive or other low symmetry patterns from an initially unpatterned system in the absence of an externally imposed template, and involves feedback generated by the interaction of transport, chemical reactions, and grain growth and dissolution. Herein we focus on the self-organizing mechanisms of reaction driven advection as they pertain to the formation of diagenetic traps.