流体饱和致密砂岩的实验测量与岩石物理建模应用

与孔隙几何相对简单的高孔隙度、高渗透率砂岩不同,致密砂岩的岩石骨架和孔隙结构通常表现出不同尺度上的非均质性。为了研究这种非均质性和饱和流体对致密砂岩性质的影响,测量了具有9块不同孔隙度的致密砂岩在不同围压、干燥和水饱和条件下的超声波速度。利用速度随有效压力的变化曲线估计软孔隙的纵横比分布,进而利用喷射流模型进行流体替换效应解释。基于超声测试数据和局部流机制,建立了考虑软孔隙随纵横比分布的喷射流模型,实现了致密砂岩跨频段岩石物理建模,利用此模型分析了超声频段和地震频段内饱和致密砂岩的流体敏感性。基于超声频段测量的流体敏感度被高估,实际上地震频段的流体敏感度可能与超声分析结果存在较大差异,因此实现实验室低频直接测量对于准确的岩石物理建模和应用具有关键性的意义。     

云冈石窟砂岩循环冻融试验研究

风化是云冈石窟目前所面临的严重的地质病害之一,温度和水分的变化是造成石窟岩体风化的重要原因,尤其是在循环冻融条件下岩体更易风化,因此,利用室内试验研究循环冻融条件下云冈石窟砂岩的物理力学性质,对于石窟岩体的稳定性评价和保护具有重要的意义。将取自云冈石窟的砂岩岩样分为饱水组、干燥组和对比组3组,通过对饱水组和干燥组岩样进行35次循环冻融试验,模拟云冈石窟砂岩的风化过程。在冻融循环开始前以及每5次冻融循环结束后,量测岩样的质量、体积,并利用超声检测分析仪对各岩样进行超声纵波测试;利用INSTRON－1346岩石伺服试验机对上述3组砂岩岩样进行单轴压缩试验,并对试验后的岩样进行SEM微观结构分析。通过试验研究,得到不同含水状态下云冈石窟砂岩岩样的冻融破坏特征以及不同循环冻融次数后岩样体积、质量、超声波纵波波速、砂岩的单轴应力-应变全过程曲线、抗压强度、抗冻系数以及微观结构的变化,分析归纳出循环冻融条件下云冈石窟砂岩的主要物理力学特性。     

受酸腐蚀砂岩的物理力学性质探讨

考虑化学作用的岩土力学性质研究及其稳定性问题已成为国内外学者关注的焦点。采用酸性环境下的加速试验,研究了钙质胶结砂岩受酸腐蚀后在不同时段单轴压缩应力作用下表现出来力学特性,揭示了砂岩腐蚀过程和受酸腐蚀机制,其目的是为酸雨环境下,研究岩体的物理力学性质提供可以借鉴的思路和方法。     

高压孔隙CO2/水作用下完整四川三叠系砂岩剪切特性的试验研究

越来越多的地下岩土工程面临着高压孔隙流体的作用,研究高压孔隙CO2/水作用下岩石的剪切特性对工程设计及安全施工、运营具有重要意义。研究了完整四川三叠系砂岩在干燥、饱和水、高压孔隙水和高压孔隙CO24种不同条件下的剪切特性。研究发现:4种不同条件下砂岩的剪切强度和残余剪切强度都会随着预剪切面上有效正应力的增加而增加;在干燥条件下随正应力的增加砂岩的剪切刚度也会增大。孔隙压力对砂岩剪切强度的影响遵守Terzaghi有效应力原理。水和CO2对砂岩的剪切强度具有显著的弱化作用,水会明显降低内摩擦角φ,而CO2对φ几乎不影响;而水对黏结力c的影响小于CO2的影响。其弱化作用可以归结为水/CO2-岩石之间的相互作用,由于砂岩中含有黏土成分,水和CO2都会对黏土产生软化作用。     

泌阳凹陷地质流体对砂岩储集层中黏土矿物形成和分布的控制作用

运用扫描电镜观察和X射线衍射分析等技术,在对泌阳凹陷核桃园组砂岩储集层中黏土矿物分布特征研究的基础上,从黏土矿物的形成条件和形成过程入手,详细分析了地质流体的性质、组成、流动方式、迁移速度等对砂岩储集层中黏土矿物形成和分布的影响,指出砂岩储集层中参与成岩反应的碎屑活性组分和地质流体对砂岩储集层中黏土矿物的形成和分布具有控制作用,并讨论了它们的形成机理.研究表明,黏土矿物的形成和分布能够不同程度地影响砂岩储集层的孔渗性和含油气性;高岭石发育层段与次生孔隙的发育和渗透率的增加以及油气的富集层段具有良好的对应关系.     

深海相中的砂质碎屑流沉积--以西藏特提斯喜马拉雅侏罗-白垩系为例

阐述了砂质碎屑流沉积特征、研究现状及其与其他深海碎屑沉积的区别 ,认为西藏特提斯喜马拉雅上侏罗统 -下白垩统深海沉积背景下的块状砂岩具有砂质碎屑流沉积性质 ,指出深海相中的块状砂岩可以预测。     

含砂岩石力学特性及其致灾机制研究

含砂岩石是发生突水溃砂灾害前在高位关键层形成的特殊岩石,其强度与力学性质均与普通岩石不同,决定着高位关键层的稳定性。研究发现：不同裂隙角的裂隙岩石与含砂岩石具有不同的特征应力,且随着裂隙角的增加,裂隙岩石与含砂岩石的起裂应力、损伤应力和峰值应力均增加,双峰应力先增加后减小。相同裂隙角下的含砂岩石各特征应力均小于裂隙岩石,说明砂体对岩石特征应力具有弱化效应。从破坏形态来看,裂隙岩石易呈现翼形拉伸裂隙,含砂岩石在低裂隙角（30°）条件下形成拉伸裂隙,高裂隙角（60°）条件下易形成剪切裂隙,表明砂体进入岩石裂隙后对岩石具有剪切效应。同时建立了充砂力学模型,指出了含砂岩石强度小于裂隙岩石的原因是砂体降低了岩石的摩擦系数。根据声发射累计振铃计数定义了岩石损伤量并分析了含砂岩石致灾机制,现场溃砂灾害可分为4个阶段：弹性变形阶段、裂隙扩展阶段、蓄砂储能阶段、溃砂释能阶段。最后利用PFC2D验证了裂隙岩石与含砂岩石的差异性,分析了不同类型岩石的能量演化规律。研究结果可作为煤矿顶板突水溃砂现象的前兆信息识别,有助于指导突水溃砂工作面的安全生产。     

煤岩应力敏感性的有限元数值模拟

利用有限元软件对比了煤岩与砂岩的裂缝界面应变分布规律。有限元法模拟表明:煤岩与砂岩的应力分布有较大差异,煤岩裂缝的最大应变约为砂岩裂缝最大应变的25倍,且煤岩裂缝更容易闭合;结合煤岩裂缝在不同有效应力条件下的应变规律,给出了渗透率随有效应力变化的拟合公式。参照砂岩应力敏感实验程序,对不同围压下煤岩渗透率变化规律的实验研究发现,渗透率拟合公式计算值与实验实测值变化趋势一致,部分煤样渗透率实验值与预测值平均误差小于10%,拟合精度高。研究结果对预测实际煤岩的渗透率和不同有效应力下煤岩渗透率的变化规律具有较高的实用价值。      

升温速率对高温作用后砂岩的宏细观性质影响

砂岩受到高温影响后,宏细观性质会发生不同程度的变化,400 ℃ 和1 000 ℃ 是砂岩宏细观性质变化的两个重要节点。为进一步研究升温速率对高温砂岩物理性质、力学特性以及内部微观破裂的影响,以350 ℃和 950 ℃为目标温度,开展了不同升温速率下砂岩的单轴压缩试验,单轴压缩全程采用声发射监测,采用扫描电镜分析了砂岩破坏后的细观形貌。研究结果表明：经不同升温速率处理后,350 ℃砂岩的质量、体积、密度以及纵波波速变化较小;950 ℃砂岩的质量、密度及纵波波速显著降低,体积明显增大,且升温速率越高变化率越小。350 ℃砂岩的应力−应变曲线、应力−体积应变曲线、抗压强度及弹性模量受升温速率的影响较小;950 ℃ 下砂岩试样应力−应变曲线随升温速率的增加向上偏移,抗压强度和弹性模量则先增加后达到恒定。350 ℃ 及 950 ℃砂岩的声发射振铃累计数都有随升温速率增大而减小的趋势,950 ℃ 下升温速率为1 ℃/min 时砂岩的振铃累计数最大,突发性声发射活动区域最多。对不同升温速率砂岩进行扫描电镜分析发现,升温速率对350 ℃ 砂岩的微观形貌影响较小,950 ℃ 砂岩试件则会随着升温速率的下降出现微裂纹和微孔隙数量增多、体积扩大的情况。     

单轴压缩下高温后砂岩的声发射特征

通过在单轴压缩下实施的声发射测试,研究焦作砂岩受20～1 200 ℃温度作用后的声发射演变过程;结合不同温度下砂岩的力学性质,通过声发射参数分析研究砂岩在不同受力阶段的声发射特点。研究表明：400 ℃以内温度对砂岩的声发射影响不太明显,在100 ℃后和600 ℃后声发射振铃累计数均发生急剧变化,100 ℃是砂岩裂纹扩展发育的门槛值,600 ℃后砂岩内部结构成分发生了变化,声发射现象较为明显。600～1 200 ℃时,砂岩呈现出明显的脆塑性转变现象,高温导致声发射信号的时间有所推迟,声发射信号增长率不断上升。1 200 ℃后,砂岩释放密集的声发射信号,呈现出塑性破坏特征。     

酸性环境下类砂岩材料物理性质的试验研究

采用酸性环境下的加速试验方法,研究了不同时段类砂岩材料-砂浆的物理性质,测试了试样的波速,监测了溶液中H+、Ca2+浓度变化,发现在不同腐蚀阶段砂浆表现出来的物理性质相差较大,其反应过程具有阶段性和Ca2+离子溶出具有滞后性的特点。同时,对不同尺寸试样受酸腐蚀的影响程度进行了对比分析,为酸雨地区岩石和混凝土性质的研究提供了可借鉴的思路。     

古海洋活性磷埋藏记录及其在氧气地球化学循环研究中的运用

在综合现代大洋活性磷循环特征的基础上,提出大洋中活性磷主要来源于有CO2 参与的地表岩石的化学风化,输入通量在很大程度上受到气候与构造等因素的控制;进入大洋之后的活性磷经光合作用而进入有机物,并通过一系列的转化环节而最终埋藏下来;从而可以根据古海洋活性磷埋藏记录,推断地质历史中不同时间尺度的全球变化情况。对比160 Ma以来的古海洋总磷和100 Ma活性磷埋藏记录可以发现:(1)总磷和活性磷的埋藏速率与长期的海平面变化明显关联,在温室气候条件下,两者基本呈正相关关系;在冰室气候条件下则为反相关关系;(2)总磷和活性磷埋藏速率的变化,能够在一定程度上指示古海洋营养条件的变化;(3)白垩纪中期大洋缺氧事件对应活性磷的埋藏峰值,表明大洋活性磷循环与大洋缺氧事件之间存在着非常密切的关系。此外还例举了古海洋活性磷埋藏记录在氧气地球化学循环研究中的运用,指出活性磷循环对大气氧含量的稳定起到重要的作用。文章最后对古海洋活性磷埋藏记录研究中存在的一些问题进行了讨论。     

Reactive transport benchmark of MoMaS

We present here the definition of the reactive transport benchmark of Groupement Mathematical Modeling and Numerical Simulation for Nuclear Waste Management Problems. The aim of this benchmark is to propose a challenging test for numerical methods used for reactive transport modeling in porous media. In order to focus on numerical methods, the problem presented here is of quite a small size, both from a hydrodynamical and from a geochemical point of view. Though the chemical coefficients used in this benchmark are not taken from a real chemical system, they are realistic, and the test case is quite challenging.     

Diffusive Reactive Transport of Multicomponent Chemicals Under Coupled Thermal, Hydraulic, Chemical and Mechanical Conditions

This paper presents an investigation of the reactive transport of multicomponent chemicals in clays under coupled thermal, hydraulic, chemical and mechanical framework, considering the diffusion processes in detail. More specifically, combined effects due to the electrochemical and the thermal diffusion potentials are investigated. A theoretical framework for coupling thermal diffusion, i.e. the Soret effect, with electrochemical diffusion in a multi-ionic system is provided. An explicit form of a definition for the thermal diffusion coefficient in a multicomponent chemical transport model is developed. Chemical transport is linked to an advanced geochemical model, PHREEQC (version 2), in order to include chemical reactions. The effects of the combined diffusion potentials on the reactive transport of multicomponent chemicals are investigated by a series of numerical simulations of coupled thermal, hydraulic and chemical behaviour.     

A global method for coupling transport with chemistry in heterogeneous porous media

Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection–diffusion PDEs coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper, a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton–Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that one be able to solve chemical equilibrium problems (and compute derivatives) without having to know the solution method. An additional advantage of the Newton–Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.     

Adsorption of reactive dyes from aqueous solutions by tannery sludge developed activated carbon: Kinetic and equilibrium studies

Adsorption kinetic and equilibrium studies of two reactive dyes, namely, Reactive Red 31 and Reactive Red 2 were conducted. The equilibrium studies were conducted for various operational parameters such as initial dye concentration, pH, agitation speed, adsorbent dosage and temperature. The initial dye concentration was varied from 10 - 60 mg/L, pH from 2–11, agitation speed from 100–140 rpm, adsorbent dosage from 0.5 g to 2.5 g and temperature from 30 °C -50 °C respectively. The activated carbon of particle size 600 μm was developed from preliminary tannery sludge. The dye removal capacity of the two reactive red dyes decreased with increasing pH. The zero point charge for the sludge carbon was 9.0 and 7.0 for the two dyes, respectively. Batch kinetic data investigations on the removal of reactive dyes using tannery sludge activated carbon have been well described by the lagergren plots. It was suggested that the Pseudo second order adsorption mechanism was predominant for the sorption of the reactive dyes onto the tannery sludge based carbon. Thus, the adsorption phenomenon was suggested as a chemical process. The adsorption data fitted well with Langmuir model than the Freundlich model. The maximum adsorption capacity(q⁰) from Langmuir isotherm were found to have increased in the range of 23.15–39.37 mg/g and 47.62–55.87 mg/g for reactive dyes reactive red 31 and reactive red 2, respectively.     

First results of operating and monitoring an innovative design of a permeable reactive barrier for the remediation of chromate contaminated groundwater

An innovative setup of a permeable reactive barrier (PRB) was installed in Willisau, Switzerland to remediate chromate contaminated groundwater. Instead of a conventional continuous barrier, this PRB consists of cylinders installed in rows: a single row for lower expected Cr^VI-concentrations and an offset double row for higher expected Cr^VI-concentrations. The cylinders are filled with reactive grey cast-Fe shavings mixed with gravel to prevent extensive precipitation of secondary phases in the pore space. The treatment of the contaminants takes place both within the cylinders and in the dissolved Fe^II plume generated downstream of the barrier. Monitoring of the contamination situation over a period of 3 a provided evidence of the mobilization, transport and behavior of the contaminants in the aquifer. Groundwater and reactive material were sampled upstream, within and downstream of the barrier by a Multi-Port Sampling System (MPSS) that revealed the geochemical processes as a function of time and space. Comprehensive chemical analyses included sensitive parameters such as Cr^VI, Fe^II/Fe^III, redox potential, dissolved O₂ and pH. Several campaigns using multiple optical tracers revealed a rather complex hydrological regime at different scales, thereby complicating the barrier performance.     

Estimation of reactive mineral surface area during water-rock interaction using fluid chemical data

Mineral dissolution and precipitation reactions actively participate to control fluid chemistry during water-rock interaction. However, it is difficult to estimate and normalize bulk reaction rates if the mineral surface area effectively participating in the reactions is unknown. In this study, we evaluated the changing of the reactive mineral surface area during the interaction between CO₂-rich fluids and albitite rock reacting under flow-through conditions. Our methodology, adopting an inverse modelling approach, is based on the measured chemical fluid composition as raw data. We estimated the rates of dissolution and the reactive surface areas of the different minerals by reconstructing the chemical evolution of the interacting fluids. This was done by a reaction process schema that was defined by a fractional degree of advance of the irreversible mass-transfer process and by attaining the continuum limit during the water-rock interaction. Calculations were carried out for albite, microcline, biotite and calcite assuming that the ion activity of dissolved silica and aluminium ions was limited by the equilibrium with quartz and kaolinite.We found that the absolute dissolution rate of albite, microcline, biotite and calcite remains essentially constant as a function of time, and the calcite dissolution rate is orders of magnitude higher than silicate minerals. On the contrary, the reactive surface area of the parent minerals varied by more than two orders of magnitude during the observed reaction time, especially for albite. We propose that the reactive surface area depends mainly on the stability of the secondary mineral coating that may passivate the effective reactive surface area of the parent minerals.     

Three-Dimensional Numerical Method of Moments for Linear Equilibrium-Adsorbing Solute Transport in Physically and Chemically Nonstationary Formations

A Lagrangian perturbation method is applied to develop a method of moments for reactive solute flux through a three-dimensional, nonstationary flow field. The flow nonstationarity may stem from medium nonstationarity, finite domain boundaries, and/or fluid pumping and injecting. The reactive solute flux is described as a space–time process where time refers to the solute flux breakthrough in a control plane at some distance downstream of the solute source and space refers to the transverse displacement distribution at the control plane. The analytically derived moments equations for solute transport in a nonstationary flow field are too complicated to solve analytically; therefore, a numerical finite difference method is implemented to obtain the solutions. This approach combines the stochastic model with the flexibility of the numerical method to boundary and initial conditions. The approach provides a tool to apply stochastic theory to reactive solute transport in complex subsurface environments. Several case studies have been conducted to investigate the influence of the physical and chemical heterogeneity of a medium on the reactive solute flux prediction in nonstationary flow field. It is found that both physical and chemical heterogeneity significantly affect solute transport behavior in a nonstationary flow field. The developed method is also applied to an environmental project for predicting solute flux in the saturated zone below the Yucca Mountain Project area, demonstrating the applicability of the method in practical environmental projects.     

Effects of medium and pore‐fluid compressibility on chemical‐dissolution front instability in fluid‐saturated porous media

Theoretical analysis and computational simulations have been carried out to investigate how medium and pore‐fluid compressibility affects the chemical‐dissolution front propagation, which is associated with a fully‐coupled nonlinear problem between porosity, pore‐fluid pressure, pore‐fluid density and reactive chemical‐species transport within a deformable and fluid‐saturated porous medium. When the fully‐coupled nonlinear system is in a subcritical state, some analytical solutions have been derived for a special case, in which the ratio of the equilibrium concentration to the solid molar density of the chemical species is approaching zero. To investigate the effect of either medium compressibility or pore‐fluid compressibility on the evolutions of chemical dissolution fronts in supercritical chemical dissolution systems, numerical algorithms and procedures have been also proposed. The related theoretical and numerical results have demonstrated that: (i) not only can pore‐fluid compressibility affect the propagating speeds of chemical dissolution fronts in both subcritical and supercritical systems, but also it can affect the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; (ii) medium compressibility may have a little influence on the propagating speeds of chemical dissolution fronts, but it can have significant effects on the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; and (iii) both medium and pore‐fluid compressibility may stabilize irregular chemical‐dissolution‐fronts in supercritical chemical dissolution systems. Copyright © 2011 John Wiley & Sons, Ltd.