饱和多孔介质化学-热-弹性模型及应用

本文基于热局部非平衡(LTNE)条件和加权平均温度概念,并假设孔隙流体由溶质和溶剂两组元组成,对页岩（饱和多孔介质）,推导给出了一种LTNE条件下的化学-热-弹性模型,同时讨论了耦合方程组的解耦求解问题．作为模型的应用,考虑无限大平面含一圆形孔的情况,研究了冷/热对流以及溶质摩尔分数突变边界条件下圆孔附近的孔隙压力和化-热应力问题,用Laplace变换得到了平面轴对称情况下有关力学变量的表达式．数值分析了圆孔边界上冷/热对流的Biot数和溶质摩尔分数改变量对圆孔附近孔隙压力和化-热应力的影响．结果表明：在Biot数为中等值(1~5)范围内,LTNE效应是非常明显的;化学作用对孔隙压力和固相应力的影响不可忽视．     

一维流体饱和粘弹性多孔介质层的动力响应

本文研究了不可压流体饱和粘弹性多孔介质层的一维动力响应问题。基于粘弹性理论和多孔介质理论,在流相和固相微观不可压、固相骨架服从粘弹性积分型本构关系和小变形的假定下,建立了不可压流体饱和粘弹性多孔介质层一维动力响应的数学模型,利用Laplace变换,求得了原初边值问题在变换空间中的解析解,并利用Laplace逆变换的Crump数值反演方法,得到原动力响应问题的数值解。数值研究了饱和标准线性粘弹性多孔介质层的动力响应,分析了固相位移、渗流速度、孔隙压力及固相有效应力等的响应特征。结果表明,与不可压流体饱和弹性多孔介质相同,不可压流体饱和粘弹性多孔介质中亦只存在一个纵波,并且固相骨架的粘性对动力行为有显著的影响。     

A semi-analytic solution of a wellbore in a non-isothermal low-permeability porous medium under non-hydrostatic stresses

The transient stress, displacement, pore pressure and temperature fields around a wellbore in a thermo-poro-elastic (THM) medium subject to non-hydrostatic remote stresses are analyzed under non-isothermal plane-strain conditions. The linear THM model proposed by Coussy (1989) is adopted in the analysis with a focus on thermal effects in low-permeability saturated rocks, characterized by a latent heat associated with local changes of fluid mass content. Non-dimensionalized parameters are identified by reformulating the fully-coupled governing equations and boundary conditions. The wellbore problem is simplified by decomposing it into axisymmetric and deviatoric loading cases. The corresponding analytical solutions are obtained in Laplace space. The inverse Laplace transforms are performed numerically to find the time-dependent distributions of field variables in the rock mass around the wellbore. These numerical results show that although the pore pressure diffusion has little influence on temperature and stress, temperature changes can strongly affect the pore pressure and stress around the wellbore. The temperature change can lead to changes in near-well stresses and the resulting significant change in wellbore breakdown pressure illustrates the importance of considering the THM coupling.     

Stress analysis of borehole subjected to fluid injection in transversely isotropic poroelastic medium

This paper proposes an analytical solution, using the combined Laplace–Fourier integral transform technique, for a borehole drilled in transversely isotropic porous medium and subjected to a fluid discharge over a finite length of its surface. Especially, the coupled boundary condition between the total radial stress and injection-induced pore pressure at the borehole surface is addressed in a rigorous fashion, which leads essentially to a set of dual integral equations that can be solved through standard numerical procedure. The study focuses on the calculation of stress fields around the borehole with particular attention given to the time-dependent effective tangential stress and pore pressure distributions. Numerical solutions are presented for verification with those recently derived for the limiting case of an isotropic medium and, more importantly, to investigate the influences of material anisotropy on the stress responses of the borehole and porous medium.     

Convective cooling/heating induced thermal stresses in a fluid saturated porous medium undergoing local thermal non-equilibrium

Local thermal non-equilibrium (LTNE) may have profound effects on the pore pressure and thermal stresses in fluid saturated porous media under transient thermal loads. This work investigates the temperature, pore pressure, and thermal stress distributions in a porous medium subjected to convective cooling/heating on its boundary. The LTNE thermo-poroelasticity equations are solved by means of Laplace transform for two fundamental problems in petroleum engineering and nuclear waste storage applications, i.e., an infinite porous medium containing a cylindrical hole or a spherical cavity subjected to symmetrical thermo-mechanical loads on the cavity boundary. Numerical examples are presented to examine the effects of LTNE under convective cooling/heating conditions on the temperature, pore pressure and thermal stresses around the cavities. The results show that the LTNE effects become more pronounced when the convective heat transfer boundary conditions are employed. For the cylindrical hole problem of a sandstone formation, the thermally induced pore pressure and the magnitude of thermal stresses are significantly higher than the corresponding values in the classical poroelasticity, which is particularly true under convective cooling with moderate Biot numbers. For the spherical cavity problem of a clay medium, the LTNE effect may become significant depending on the boundary conditions employed in the classical theory.     

Thermo-Osmosis Effect in Saturated Porous Medium

In this paper, the thermo-poroelasticity theory is used to investigate the quasi-static response of temperatures, pore pressure, stress, displacement, and fluid flux around a cylindrical borehole subjected to impact thermal and mechanical loadings in an infinite saturated poroelastic medium. It has been reported in literatures that coupled flow known as thermo-osmosis by which flux is driven by temperature gradient, can significantly change the fluid flux in clay, argillaceous and many other porous materials whose permeability coefficients are very small. This study presents a mathematical model to investigate the coupled effect of thermo-osmosis in saturated porous medium. The energy balance equations presented here fulfill local thermal non-equilibrium condition (LTNE) which is different from the local thermal equilibrium transfer theory, accounting for that temperatures of solid and fluid phases are not the same and governed by different heat transfer equations. Analytical solutions of temperatures, pore pressure, stress, displacement, and fluid flux are obtained in Laplace transform space. Numerical results for a typical clay are used to investigate the effect of thermo-osmosis. The effects of LTNE on temperatures, pore pressure, and stress are also studied in this paper.     

Analytical Determination of Flow Patterns in a Double Porosity Medium Containing Ellipsoidal Occluded Macro-Voids

This article presents the analytical study of fluid flow in a porous medium presenting pores of two different length scales: at the smallest or microscopic scale, the presence of connected voids confers a porous medium structure to the material investigated, while at the upper or mesoscopic scale, occluded macro-pores are present. This microstructure is employed to represent the progressive opening of inter-aggregate pore spaces observed in natural compacted montmorillonites polluted by heavy metal ions. Three-dimensional analytical expressions are rigorously derived for the pore fluid velocity and excess pore fluid pressure within the porous matrix, around an occluded ellipsoidal inter-aggregate void. The eccentricity ratio is employed to characterize the geometrical shape of the ellipsoidal void, while its size is characterized by the macro-porosity. Confrontations are made with numerical solutions in order to investigate the applicability of the analytical pressure and velocity solutions to microstructures of finite size.     

粘弹性饱和土体中半封闭圆形隧洞的动力响应分析

基于Biot波动方程,研究分析了粘弹性饱和土体中半封闭圆形隧洞的动力响应问题．假定衬砌材料为多孔介质,引入了更符合工程实际的半封闭边界条件．通过引入势函数,在Laplace变换域中得到隧洞边界上作用轴对称荷载和流体压力条件下应力、位移和超孔隙水压力的解答．利用Laplace数值逆变换得到时域中的解,分析了隧洞边界的半透水特性对隧洞动力响应问题的影响,结果表明：隧洞边界的半透水特性对应力、位移场的变化和超孔隙水压力的消散有很大的影响,透水和不透水下条件的解仅是本文的两个特例。     

Relaxation Effects of a Saturated Porous Media Using the Two-Dimensional Generalized Thermoelastic Theory

A model of the equations of a generalized thermoelasticity (GT) with relaxation times for a saturated porous medium is given in this article. The formulation can be applied to the GT theories: Lord–Shulman theory, Green–Lindsay theory, and Coupled theory for the porous medium. A two-dimensional thermoelastic problem that is subjected to a time-dependent thermal/mechanical source is investigated with the model of the generalized porous thermoelasticity. By using the Laplace transform and the Fourier transform technique, solutions for the displacement, temperature, pore pressure, and stresses are obtained with a semi-analytical approach in the transform domain. Numerical results are also performed for portraying the nature of variations of the field variables. In addition, comparisons are presented with the corresponding four theories.     

Analysis of hydrodynamic and thermal dispersion in porous media by means of a local approach

A pore scale analysis is implemented in this numerical study to investigate the behavior of microscopic inertia and thermal dispersion in a porous medium with a periodic structure. The macroscopic characteristics of the transport phenomena are evaluated with an averaging technique of the controlling variables at a pore scale level in an elementary cell of the porous structure. The Darcy–Forchheimer model describes the fluid motion through the porous medium while the continuity and Navier–Stokes equations are applied within the unit cell. An average energy equation is employed for the thermal part of the porous medium. The macroscopic pressure loss is computed in order to evaluate the dominant microscopic inertial effects. Local fluctuations of velocity and temperature at the pore scale are instrumental in the quantification of the thermal dispersion through the total effective thermal diffusivity. The numerical results demonstrate that microscopic inertia contributes significantly to the magnitude of the macroscopic pressure loss, in some instances with as much as 70%. Depending on the nature of the porous medium, the thermal dispersion may have a marked bearing on the heat transfer, particularly in the streamwise direction for a highly conducting fluid and certain values of the Peclet number.     

The Flow Analysis of Viscoelastic Fluid with Fractional Order Derivative in Horizontal Well

The fractional calculus approach is introduced into the seepage mechanics. A three-dimensional relaxation model of viscoelastic fluid is built. The models based on four boundary conditions of exact solution in Laplace space for some unsteady flows in an infinite reservoir is obtained by using the Laplace transform and Fourier sine and cosine integral transform. The pressure transient behavior of non-Newtonian viscoelastic fluid is studied by using Stehfest method of the numerical Laplace transform inversion and Gauss–Laguerre numerical integral formulae. The viscoelastic fluid is very sensitive to the order of the fractional derivative. The change rules of pressure are discussed when the parameters of the models change. The plots of type pressure curves are given, and the results can be provided to theoretical basis and well-test method for oil field.     

Filtration in a Porous Granular Medium: 1. Simulation of Pore-Scale Particle Deposition and Clogging

This paper presents a numerical model for simulating the pore-scale transport and infiltration of dilute suspensions of particles in a granular porous medium under the action of hydrodynamic and gravitational forces. The formulation solves the Stokes’ flow equations for an incompressible fluid using a fixed grid, multigrid finite difference method and an embedded boundary technique for modeling particle–fluid coupling. The analyses simulate a constant flux of the fluid suspension through a cylindrical model pore. Randomly generated particles are collected within the model pore, initially through contact and attachment at the grain surface (pore wall) and later through mounding close to the pore inlet. Simple correlations have been derived from extensive numerical simulations in order to estimate the volume of filtered particles that accumulate in the pore and the differential pressure needed to maintain a constant flux through the pore. The results show that particle collection efficiency is correlated with the Stokes’ settling velocity and indirectly through the attachment probability with the particle–grain surface roughness. The differential pressure is correlated directly with the maximum mound height and indirectly with particle size and settling velocity that affect mound packing density. Simple modification factors are introduced to account for pore length and dip angle. These parameters are used to characterize pore-scale infiltration processes within larger scale network models of particle transport in granular porous media in a companion paper. Articlenote: Currently at GZA GeoEnvironmental Inc., 1 Edgewater Drive, Norwood, MA 02062, U.S.A.     

Plane strain consolidation of a multi-layered soil with anisotrpic permeability and compressible constituents

对多层地基的平面应变固结问题进行了研究,并同时考虑了土体的渗透各向异性和孔隙 流体的可压缩性. 从平面应变Biot固结的控制方程出发,对时间t, 坐标z和x进行 Laplace和Fourier变换,建立了地基表面(z=0)和任意深度z处的基本量 在Laplace-Fourier变换域内的传递矩阵关系. 利用传递矩阵 法,结合土层连续条件和边界条件,并应用Laplace-Fourier逆变换技术,推导出渗透各向 异性可压缩多层地基平面应变固结的理论解. 基于该解,编制了计算程序,并进行了 数值计算. 讨论了土体的渗透各向异性、孔隙流体的可压缩性以及地基的分层特性对地基固 结的影响,分析结果表明：土体的渗透各向异性、孔隙流体的可压缩性,以及地基的分层特 性对地基的固结行为有着重要的影响.     

饱和土中球空腔的瞬态动力响应

采用工程上通用的饱和土力学模型,考虑流体与固全之间的耦合作用,利用Ｌａｐｌａｃｅ变换求解了饱和土中球空腔的瞬态动力响应问题,得到了变换域内的解析解,借助数值Ｌａｐｌａｃｅ换求解了饱和土中球空腔瞬态动力响应的位移、应力及孔压的变化规律,为分析地下结构瞬态动力响应提供了一种有效的方法,模型符合工程实际,有一定的工程应用价值。     

State space approach to unsteady magnetohydrodynamics natural convection heat and mass transfer through a porous medium saturated with a viscoelastic fluid

A technique of the state space approach and the inversion of the Laplace transform method are applied to dimensionless equations of an unsteady one-dimensional boundary-layer flow due to heat and mass transfer through a porous medium saturated with a viscoelastic fluid bounded by an infinite vertical plate in the presence of a uniform magnetic field is described. Complete analytical solutions for the temperature, concentration, velocity, and induced magnetic and electric fields are presented. The inversion of the Laplace transforms is carried out by using a numerical approach. The proposed method is used to solve two problems: boundary-layer flow in a viscoelastic fluid near a vertical wall subjected to the initial conditions of a stepwise temperature and concentration and viscoelastic fluid flow between two vertical walls. The solutions are found to be dependent on the governing parameters including the Prandtl number, the Schmidt number, the Grashof number, reaction rate coefficient, viscoelastic parameter, and permeability of the porous medium. Effects of these major parameters on the transport behavior are investigated methodically, and typical results are illustrated to reveal the tendency of the solutions. Representative results are presented for the velocity, temperature, concentration, and induced magnetic and electric field distributions, as well as the local skin-friction coefficient and the local Nusselt and Sherwood numbers.     

Dynamic response of an axially loaded rigid sphere embedded in a saturated poroelastic medium

In this paper, an analytical solution for the response of a rigid sphere embedded in a full space poroelastic medium subjected to a dynamic lateral load is derived. The solution is obtained using Biots theory for acoustic waves. In this solution, the displacements of the solid skeleton and the pore pressure are expressed in terms of three scalar potentials. These potentials correspond to the wave velocities of the slow and fast compressional waves and to the shear wave. The governing equation for the dynamic motion is expressed in the frequency domain using Fourier transformation. Different boundary and load conditions were investigated. Curves showing variation in the fluid pressure and solid displacements with the loads frequency were plotted in non-dimensional forms.     

Asymptotic analysis of the equilibrium equation of a fluid-saturated porous medium by the homogenization method

The homogenization of static elasticity equations describing the stress strain state of fluid-saturated porous medium is considered. In this paper, the homogenization method is used to determine the pore pressure transfer tensor, which (a coefficient in the isotropic case) is an important parameter influencing the stress-strain state of fluid-saturated rocks. It shows what a part of the pressure in the fluid is “active” in the formation of macroscopic strains.The pore pressure transfer tensor is calculated for model and real geological specimens. The dependence of this tensor on the porosity, pore shape, and Poisson ratio is investigated. The use of the computational technique for determining the effective properties of rocks shows that it is practically important in the engineering geology.     

黏弹性饱和多孔介质中圆柱孔洞的频域响应

研究了无限黏弹性饱和多孔介质中圆柱孔洞(有衬砌)表面受轴对称简谐荷载和流体压力 作用下的频域响应问题. 引入Carcione提出的本构模型来描述介质的流变和松弛性质. 考 虑衬砌和介质的相对渗透性，孔洞处于半封闭状态，边界半透水. 引入两个势函数，在频域 中得到了应力、位移和超孔隙水压力响应解答. 并进行了算例分析，讨论了反映介质黏弹 性性质的最小质量因子，反映孔洞边界半透水性质的渗透性参数及衬砌和介质的相对刚度对 问题的影响. 分析结果表明：以上参数对圆柱孔洞的频域响应有很大影响.     

New exact solutions of Stokes' second problem for an MHD second grade fluid in a porous space

We investigate a problem describing the oscillating flow of an incompressible magnetohydrodynamic (MHD) second grade fluid in a porous half space. Exact solutions for sine and cosine oscillations are developed by applying the Laplace transform method. The total obtained solution is a sum of steady and transient solutions. Particular attention is given to the effects of magnetic and porous medium parameters on the velocity. It is shown that previous results for a non-porous medium and hydrodynamic fluid are the limiting cases of the present problem. The results for velocity are plotted and discussed carefully.     

Impact on the boundary of a compressible two-layer fluid

Within the framework of the acoustic approximation a solution of the plane nonstationary problem of impact on a fluid boundary is found. The fluid occupies the lower half-plane and consists of two layers with given speeds of sound and densities. The upper layer has a constant depth and is bounded above by a plate with a given normal velocity. The solution is constructed using the Fourier and Laplace integral transforms. Numerical calculations are performed for piston impact across a rigid screen and the impact of a jet with an aerated head on a rigid wall. It is shown that the presence of an interlayer with reduced speed of sound and/or density considerably changes the evolution of the hydrodynamic pressure distribution over the impacting surface: the absolute pressure maximum decreases but pressures of significant amplitude are maintained for a longer time than for a homogeneous fluid.