Image‐based modelling and analysis of microstructures for two‐scale problems in geomechanics

This paper focuses on the geometry modelling and numerical analysis of microstructures of geomaterials employing the concept of image‐based engineering. The novel modelling and analysis techniques with digital images are incorporated with the mathematical homogenization method to study the interaction between individual phases, each of whose shape and spatial distribution are irregular. Owing to the distinctive features of these computational techniques, the evaluation of homogenized properties for geomaterials provides the reliable information about the micro‐ or macroscopic mechanical behaviours for engineering practice. It is, naturally, inevitable that engineers' demands on safety and efficient design place emphasis on quantitative estimates for these values. Thus, calibration accompanied with actual measurements comes within the scope of this study so that these properties would be realistic and practical from the engineering viewpoints. Copyright © 2002 John Wiley & Sons, Ltd.     

Multiscale poro‐creep model for cement‐based materials

Experimental observations clearly show that the relative humidity (h_r) conditions influence significantly the creep behavior of cement‐based materials, indicating that the water present within these materials plays a crucial role. This work presents a creep model for hardened cement pastes (HCP), based on a multiscale homogenization approach. It takes into account both free and adsorbed water contained in the porosity and investigates their effects on the HCP macroscopic creep behavior. The calcium silicate hydrate phase is assumed to be linear viscoelastic, and the Mori–Tanaka scheme is applied in the Laplace–Carson space to the composite formed of porosity, calcium silicate hydrate, and the other main hydrated compounds (which behavior is linearly elastic) by making use of the correspondence principle. With this model, estimations of the evolution of the macroscopic creep behavior of HCP submitted to constant external loading are examined under different h_r and compared with available experimental data. Finally, a method for implementing the model in a finite element code is proposed, and simulations of standard creep tests are performed to assess its validity. Copyright © 2011 John Wiley & Sons, Ltd.     

Compressibility and permeability of sand–kaolin mixtures. Experiments versus non‐linear homogenization schemes

This study deals with the behaviour of mixtures of sand and saturated kaolin paste considered as composite materials made of permeable and deformable (with non‐linear behaviour) matrix (the kaolin paste) with rigid and impervious inclusions (the sand grains). Oedometric and permeability tests highlight the key role of the state of the clay paste, and show the existence of a threshold of sand grain concentration above which a structuring effect influences both compressibility and permeability. At the light of these experiments two homogenization schemes (with simplifying assumptions to make the problem manageable) are considered to model these two parameters. Qualitative and quantitative comparisons with experimental data point out their respective domain of interest and limitations: a tangent homogenization scheme is shown to be sufficient to describe the macroscopic properties for dilute sand concentration; above the concentration threshold, the structuring effect is captured by the new homogenization scheme developed in this paper. Copyright © 2010 John Wiley & Sons, Ltd.     

利用激发极化驰豫时间谱确定渗透率的实验研究

激发极化驰豫时间是孔隙大小的函数,采用奇异值分解方法对激发极化衰减曲线进行多指数反演,得到光滑的驰豫时间谱,该谱呈现多峰特性.实验结果表明,溶质对该谱的影响主要是由于溶质离子的扩散系数不同所致.对于同一块岩心,衰减谱随扩散系数的增大而向小驰豫时间方向移动.结合岩心地层因素F,研究了激发极化驰豫时间几何平均值Tg与岩石的渗透率K的关系,结果表明,logK和log(T1/2g/F)、log(Tg/F)呈现良好的线性关系.     

复杂应力路径下大理岩三轴渗透试验研究

在隧道施工过程中,围岩的应力条件非常复杂,研究复杂应力路径下围岩渗透性能的变化规律对在高水压地区修建隧道具有重要意义。以锦屏水电枢纽二级电站交通辅助洞大理岩为研究对象,进行了常规三轴渗透试验与控制轴向应变（简称应变 ）、围压先升后降的三轴渗透试验,探讨了轴压与围压之差的绝对值 （简称应力差）与渗透率的关系。试验结果表明: (1)常规渗透试验大理岩渗透率的变化过程有3个重要的特征点：渗透率最低点、渗透率峰值点以及渗透率稳定点;(2)大理岩的渗透率随着应力差的增大而减小,两者呈负指数关系;(3)对某一固定的应力差,升围压阶段测得的渗透率大于降围压阶段测得的渗透率;(4)对任何试验岩样总存在一阈值,当应力差小于该值时,应力差的改变对渗透率有显著影响;(5)应力差减小过程中岩样渗透率的“恢复能力”随着岩样轴向应变的增大而逐渐减弱。     

On computation of strain‐dependent permeability of rocks and rock‐like porous media

Determination of transport properties of geomaterials is an important issue in many fields of engineering analysis and design. For example, in petroleum engineering, in situ permeability of an oil reservoir may be crucial in establishing its viability for exploitation, whilst prevention of leakage from underground storage facilities for oil and gas, nuclear waste as well as viability of CO₂ sequestration projects crucially depends on its long‐term values. Permeability is indirectly related to the porosity, pore‐size distribution and pore architecture of the porous media. These parameters evolve when a strain field is imposed. Physical measurement of permeability under a strain field in laboratory conditions is difficult, expensive and prone to a number of uncertainties. In the past, pore network models have been used to compute permeability of materials under stress/strain‐free conditions. In this paper, we propose an enhanced pore network model to compute permeability of rocks and rock‐like porous media under a stress/strain field. Data of pore‐size distribution obtained from mercury intrusion porosimetry are used to compute permeability of rock samples from various unspecified oilfields in the world. It is shown that the two permeabilities can be predicted from the model with sufficient accuracy. A hypothesis for change in porosity, pore‐size distribution and pore architecture as a result of imposed mechanical strains is then proposed. Based on this, permeability is computed again for one of the rock samples under uniaxial and triaxial compressive and tensile strain fields. It is shown that depending on the state of strain field imposed, permeability evolves in an anisotropic manner. Permeability under tensile strain field increases dramatically compared with the reduction that takes place under compressive strain field of the same magnitude. Copyright © 2014 John Wiley & Sons, Ltd.     

不饱和膨润土中气体渗透研究

作者对不饱和膨润土中气体渗透情况进行了研究，对实验所使用的膨润土和砂进行了密度、容重、粒重、孔隙度及含水率等特性的表征。实验所建立的气体渗透研究装置可用于测定气体渗透系数，并将其作为材料气体渗透性能的表征，实验结果表明：膨润土的孔隙度和含水率对其渗透系数有十分显著的影响。     

Numerical modelling on water retention and permeability of compacted GMZ bentonite under free-swelling conditions

In this study, a solid–liquid–gas coupled equation was established to simulate water retention characteristics of highly compacted GMZ bentonite. Then, modelling results were compared with laboratory test results. Results indicate that GMZ bentonite has a strong moisture expansion (or a limit drying shrinkage) characteristic. The control equation can simulate the water absorption and deformation characteristics very well at high relative humidity (or low suction). Environmental scanning electron microscope (ESEM) observation reveals the course grain soil texture of the surface under low relative humidity (RH), while the surface of GMZ bentonite becomes smooth (more fine-grained soil texture) as RH increases. Differences were found between the porosities calculated by macroexperiment results and microscopic observations with ESEM method. This is because only the interaggregate pores can be observed by ESEM photographs. Additionally, we find that the simulated effective porosities are close to the results calculated by microscopic tests, while the effective porosity is considered as the main flow channel of flow. Further, the intrinsic permeability, the effective water and gas permeability are calculated based on the proposed model. The modelling results coincide well with the laboratory experimental results and support the reliability of the proposed model.     

3‐D mesoscale simulation of crack‐permeability coupling in the Brazilian splitting test

In this paper, mesoscale hydromechanical simulations are performed to study (1) fracture features and (2) crack‐gas permeability coupling evolution in the context of the tensile splitting test. The mesostructure is based on a 2‐phase 3‐D representation of heterogeneous materials, such as concrete, where stiff aggregates are embedded into a mortar matrix. To take into account these heterogeneities without any mesh adaptation, a weak discontinuity is introduced into the strain field. In addition, a strong discontinuity is also added to take into account microcracking. This mechanical model is cast into the framework of the enhanced finite element method. Concerning the coupling with gas permeability, a double‐porosity method is used to simulate the flow through the cracks and the porosity. The apparent gas permeability is afterwards evaluated by a homogenization method. On the basis of finite element simulations, influence of aggregate size on ultimate crack opening, macroscopic ultimate tensile stress, total dissipated energy, and gas permeability evolution is numerically investigated. Furthermore, gas permeability evolution is also compared with experimental results from the literature. In addition, in the spirit of a sequential multiscale approach, macroscale gas permeability equations are identified from the hydromechanical results coming from the mesoscale computations. These equations lead to a relation between macroscale gas permeability evolution and crack opening. Besides, we show how the aggregate size influences the percolation threshold and that after this threshold, a cubic relation between macroscale gas permeability and crack opening is obtained.     

PPTS-1压力传递型泥页岩渗透性测试系统

井壁失稳一直是困扰油气勘探开发的难题,而泥页岩渗透性是评价泥页岩井壁稳定性的重要指标之一。介绍了一套基于LabVIEW软件平台的PPTS-1压力传递实验系统平台,采用岩心分析法和压力传递测试技术精确测试泥页岩的渗透率。实验结果表明,该实验平台操作简单,性能稳定,测试精度高,为井壁稳定性深入研究提供了支撑。     

基于电缆地层测试的储层径向有效渗透率计算新方法

试井方法是求取地层径向有效渗透率最直接有效的方法,但其测试费用较高。提出了基于电缆地层测试结合核磁共振测井获得地层径向有效渗透率的新方法。应用电缆地层测试仪器求取地层纵向有限测试点的有效渗透率水平分量,进而应用各测试点有效渗透率水平分量对核磁共振渗透率进行刻度,计算刻度后的核磁测井渗透率平均值,可获得既能反映测井表现的地层纵向非均质,又能表现电缆地层测试渗透率的有效性。该方法应用于现场测试,与DST (钻柱地层测试)测试结果对比具有较高的计算精度,同时可以节省大量的测试费用。     

周期振荡法在低渗透测量中的应用研究

周期振荡法是一种新兴的多孔介质渗透率测量方法,具有测量周期短、稳定性好、精度高等优点。然而,由于缺乏对周期振荡法的系统研究,该方法未能在低渗透测量领域得到推广应用。采用数值模拟方法,分别研究了不同孔隙率、渗透率对周期振荡法气压测量结果的影响,以及不同形式周期波在周期振荡法渗透率测量中的适用性。结果表明:无论是孔隙率,还是渗透率都会影响气体压力传递过程;周期振荡法渗透率测量时可以不需要初始气压平衡过程,一定时间后测量结果几乎不受初始气压状态影响;与正弦波相比,方波对下游气体压力响应的影响更明显,而且方波、三角波和锯齿波具有容易加载的特点,建议采用方波取代正弦波进行周期振荡法试验。研究成果对指导周期振荡法试验和解决低渗透测量难题有一定帮助。     

One‐dimensional analytical solution for mesoscopic flow induced damping in a double‐porosity dual‐permeability material

Heterogeneities, such as fractures and cracks, are ubiquitous in porous rocks. Mesoscopic heterogeneities, that is, heterogeneities on length scales much larger than typical pore size but much smaller than the wavelength, are increasingly believed to be responsible for significant wave energy loss in the seismic frequency band. When a compressional wave stresses a material containing mesoscopic heterogeneities, the more compliant parts of the material (e.g., fractures and cracks) respond with a greater fluid pressure than the stiffer portions (e.g., matrix pores). The induced fluid flow, resulting from the pressure gradients developed on such scale, is called mesoscopic flow. In the present study, the double‐porosity dual‐permeability model is adopted to incorporate mesoscopic heterogeneities into rock models to account for the attenuation of wave energy. Based on the model, the damping effect due to mesoscopic flow in a one‐dimensional porous structure is investigated. Analytical solutions for several boundary‐value problems are obtained in the frequency domain. The dynamic responses of infinite and finite porous layer are examined. Numerical calculations show that the damping effect of mesoscopic flow is significant on the pore pressure response and the resulting effective stress. For the displacement, the effect is seen only at the very low frequency range or near the resonance frequencies. Copyright © 2017 John Wiley & Sons, Ltd.     

基于孔压静力触探确定土体的渗透系数

近年来许多研究者致力于通过孔压静力触探(CPTU)确定土体的渗透系数,但是存在偏差大,而且不能从机制上清楚地解释孔压静力触探确定土体渗透系数的原因。为了应用孔压静力触探现场快速地确定土体的渗透系数,在总结前人对孔压静力触探确定土体渗透系数方法的研究基础上,提出孔压静力触探时锥头处呈球冠形孔隙水流动的计算模型,假定锥头处初始超孔隙水压力为负指数型分布,推导出土体渗透系数的表达式,并且结合现场实测数据,将所提出的方法与已有的方法进行对比分析。研究发现：所提出的方法扩大了现有方法的应用范围;渗透系数计算结果受锥头角度的影响较大,随着锥头角度的增大而增大;针对国际上通用的60°角的标准锥头,所提出的方法计算的渗透系数大于前人的方法,而且更接近室内外试验结果。     

超临界CO2温变对低渗透煤层孔渗变化的实验研究

为解决我国高瓦斯煤层渗透性差导致瓦斯抽采率低的难题,利用超临界二氧化碳强扩散和溶解增透等独特优点,采用自制三轴渗透实验装置,开展不同温度下超临界二氧化碳作用后煤的宏观增透实验,在宏观增透实验基础上进行煤微观扫描成像实验。结果表明：恒定体积应力和孔隙压力条件下,不同温度超临界二氧化碳作用后,煤的渗透率较增透前提高一个数量级,但在二氧化碳的超临界温度范围内,煤的渗透率随温度增加呈负指数变化规律。超临界二氧化碳作用后,煤微观孔隙率较增透前提高两个数量级,随着温度增加,煤微观孔裂隙的演化速率减慢,孔隙率随温度增加呈负指数变化规律。宏微观实验数据同时表明,煤宏观渗透率随微观孔隙率增加而增大。超临界二氧化碳增透过程中,孔隙压力对低渗透煤层的增透效果起主控作用。     

煤体瓦斯渗透率的PSO-LSSVM预测模型

结合有关煤体渗透率的众多研究成果,总结出影响煤体渗透率的3个主要因素为有效应力、温度和瓦斯压力。采用粒子群优化算法(PSO)对最小二乘支持向量机(LSSVM)的参数进行优化选择,并以上述3个因素和抗压强度为输入值,以渗透率为目标输出值,建立煤体瓦斯渗透率的PSO-LSSVM预测模型。利用25组数据进行PSO-LSSVM模型与BP神经网络及支持向量机的比较实验,PSO-LSSVM预测结果与实际值拟合程度优于其他两个模型,且具有更小的误差。实验结果表明,采用PSO-LSSVM模型可由有效应力、温度和瓦斯压力对渗透率进行较高精度的预测。     

非均匀多孔介质等效渗透率的普适表达式

多孔介质渗流是普遍的物理过程,涉及地下工程、地热开采、环境工程等各行各业,尤其是工程建设,常面临防渗问题。由于地质条件的复杂性,工程区域地层受到成岩、压实、风化、生物作用等各种影响,故渗流性质复杂,常需要对建设区域的渗流状况进行数值模拟,从而为工程的设计施工提供决策依据。数值仿真结果依赖于对地层介质关键参数的选取,但目前工程多将其视为均匀介质处理,对于介质的非均匀特性考虑较少。文章旨在研究非均质多孔介质渗透率空间分布与等效渗透率的关系。基于连续介质假定、达西定律以及非均匀多孔介质渗透率空间分布函数,建立一维到三维的达西渗流问题模型,通过求解偏微分方程和理论推导,得到基于渗透率空间分布函数的等效渗透率理论表达式,并与有限元计算的数值解进行对比分析,结果表明理论值和数值解误差很小,证明等效渗透率的表达式的合理性。利用该成果可通过多点局部渗透率的测定构建渗透率空间分布函数,从而对整体渗流区域的渗透性质进行快速计算和评估,从而简化异常复杂的工程地质模型以减少计算量需求,对于工程仿真的快速计算和结果评估有重要意义。