两相饱和介质半无限空间位错位移场

弹性位错理论是研究地震震源正演和反演的重要工具之一. 单相介质位错理论求得的地面变形场是稳定的、不变的，而实际地震后的地面变形场是变化的. 本文利用集中力作用下无限空间中两相饱和介质Green函数的柱函数表达式，通过弹性运动方程积分解的Helmholtz方程的Hansen矢量变换，再叠加自由表面影响场的方法求得两相饱和介质半无限空间位错位移场. 本文结果结合Nur的DD模式，在进水排水的假说下能解释地震后地面变形场的变化.     

Solution of a Green’s function for a saturated porous medium in a half-space subjected to a torsional force

Based on the solutions of the Green’s function for a saturated porous medium obtained by the authors,and using transformation of axisymmetric coordinates,Sommerfeld integrals and superposition of the influence field on a free surface,the authors have obtained displacement solutions of a saturated porous medium subjected to a torsional force in a half-space.The relationship curves of the displacement solutions and various parameters(permeability,frequency,etc.)under action of a unit of torque are also given in this paper.The results are consistent with previous Reissner’s solutions,where a two-phase medium decays to a single-phase medium.The solution is useful in solving relevant dynamic problems of a twophase saturated medium in engineering.     

关于两相饱和介质中流相Green函数的解析解

自然界中许多动力学问题需要用两相饱和介质Green函数求得解答,本文在两相饱和介质位移场Green函数解答的基础上,求得的有关饱和流体孔隙压力,进、排流体流量等三个Green函数.从而使两相饱和介质涉及流体问题的Lamb方程能完整地使用.论文在推导过程中详细地阐述了流体三个Green函数的物理意义,并将所得结果与Chen（1994）的纯数学的结果作了对比.最后文章讨论了解析解的稳定性.     

两相饱和介质中的集中力点源位移场解与应用

通过两相饱和介质Ｂｉｏｔ方程变换,利用势场分解较好地解决了饱和两相介质中Ｐ１和Ｐ２波的耦合问题．再由Ｐｏｉｓｓｏｎ方程和Ｈｅｌｍｈｏｔｚ方程的特性,求解得到两相饱和介质在集中力点源作用下的位移场Ｇｒｅｅｎ函数,进而通过阶跃函数和Ｆｏｕｒｉｅｒ逆变换,求得两相介质在集中力．作用下的波场．由上述结果,根据扩容进水模型,合理地解释了孕震前各主要阶段在地震记录中垂直向振幅与振幅比统计值的变化．最后结合工作实例,利用两相介质波场理论求得由一般弹性介质理论不易求得的骨架固体力学参数．     

Two-phase pore-fluid distribution in fractured media: acoustic wave velocity vs saturation

The acoustic wave velocity varies with fluid saturation and pore-fluid distribution. We use a P-wave source and the staggered grid finite-difference method, with second-order accuracy in time and eighth-order accuracy in space, to simulate the acoustic wave field in a fractured medium that is saturated with a two-phase pore fluid (gas & water). Further, we analyze the variation of acoustic wave velocity with saturation for different pore-fluid distribution modes. The numerical simulation method is simple and yields accurate results.     

处于海洋环境的桩柱在瑞利波作用下的响应

本文以一海洋深水环境桩柱为例,考虑到海洋地基为两相饱和土介质,基于Biot两相饱和多孔介质理论,通过数值算例分析了饱和介质中瑞利波加速度、位移分布规律并与单相弹性介质地基条件下进行了对比,结果表明,瑞利波主要作用于地表附近且振动幅度随深度衰减很快,饱和地基下的瑞利波影响深度要大于单相弹性地基,在瑞利波地震作用下,考虑两相饱和介质下的桩柱位移响应大于单相弹性介质地基。并且由于桩柱处海洋深水环境,本文对环境水体对桩柱振动响应的影响进行了研究,结果表明,流固耦合效应可明显降低桩柱运动响应,在海洋深水环境地震响应分析时,有必要考虑饱和地基以及流固耦合效应对桩柱运动响应的影响。     

无限饱水孔隙弹性空间中由运动源所产生的位移场

基于Ｂｉｏｔ的饱和孔隙弹性介质的运动方程和利用复变函数方法，本文研究了无限饱水孔隙弹性空间中由常速运动源所产生的位移场。考虑了两种类型源：ａ．沿无限空间水平轴运动的斜向集中力源；ｂ．运动双力偶源。关于力源的运动速度，考虑了４种情形：ａ．力源的运动速度Ｕ小于饱水孔隙弹性介质的三种体波速度一亚音速情形；ｂ．速度Ｕ小于介质的第一纵波过度和横波速度，但大于第二纵波速度─—弱跨音速情形；ｃ．速度Ｕ小于第一纵波速度，但大于横波速度和第二级波速度─—强跨音速情形；ｄ．速度Ｕ大于介质的所有三种体波速度─—超音速情形。结果表明，在跨音达和超音速情形里，解呈现出与力源相联系的平面冲击波特征，位移出现了相应的跳跃。     

两相饱和介质中的集中力点源Green函数

通过对两相饱和介质 Ｂｉｏｔ 方程的变换,利用 Ｐｏｉｓｓｏｎ 方程和 Ｈｅｌｍ ｈｏｌｔｚ 方程性质,求解得到两相饱和介质在集中力点源作用下的位移场 Ｇｒｅｅｎ 函数     

基于Comsol的饱和多孔介质动力方程的数值模拟及应用

两相饱和多孔介质的动力响应问题在地震工程领域具有重要的研究意义,由于涉及到固相和液相的动力耦合,使得该问题的求解尤为复杂。本文利用Comsol在求解多场耦合问题上的优点,针对Biot饱和多孔介质u-U耦合形式下的波动方程特征,经过一系列微分算子运算和矩阵变换得到导数形式下的波动方程,基于Comsol Multiphysics提供的广义偏微分方程模式对变形后的波动方程进行求解,并把改进后的无限元边界应用到无限域动力问题的模拟中。通过与饱和多孔介质动力响应的解析解进行对比,验证模型求解技术的可行性和正确性,并在此基础上讨论饱和土地基中空沟隔振效果与饱和土体参数孔隙率、泊松比的关系。通过研究分析,可以为饱和土地基中空沟隔振设计提供一些有价值的参考。     

双相各向同性介质弹性波高精度波场分离数值模拟方法

为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法.采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出同相和流相分量相应的散度场(纯纵波场)和旋度场...     

Biot流动和喷射流动耦合作用下波传播的FCT紧致差分模拟

本文从含流体多孔隙介质中同时包含Biot流动和喷射流动两种力学机制的BISQ(Biot-Squirt)方程出发,利用FCT(Flux-Corrected Transport)紧致差分方法数值模拟了在Biot流和喷射流共同作用下的波在含流体多孔隙各向同性介质中的传播.通过与仅受Biot流动作用下的波场结果对比,我们研究了...     

饱水介质波动分析的谐波输入初始间断的处理方法

本文以饱水两相介质的土力学模型为研究对象,在假定两相介质为弹性介质条件下,采用了显式有限元法和透射边界进行了饱和弹性半空间动力响应问题的研究。为避免谐波输入初始间断的影响,文中提出了一个处理函数,并以弹性半空间为算例,对饱水介质和单相介质分别进行了在底边界P波垂直入射时的动力响应分析,验证了该处理函数的有效性和实用性。     

Inverting reservoir parameters in a two-phase fractured medium with a niche genetic algorithm

In this paper we calculate a synthetic medium surface displacement response that is consistent with real measurement data by applying the least-square principle and a niche genetic algorithm to the parameters inversion problem of the wave equation in a two-phase medium. We propose a niche genetic multi-parameter (including porosity, solid phase density and fluid phase density) joint inversion algorithm based on a two-phase fractured medium in the BISQ model. We take the two-phase fractured medium of the BISQ model in a two-dimensional half space as an example, and carry out the numerical reservoir parameters inversion. Results show that this method is very convenient for solving the parameters inversion problem for the wave equation in a two-phase medium, and has the advantage of strong noise rejection. Relative to conventional genetic algorithms, the niche genetic algorithm based on a sharing function can not only significantly speed up the convergence, but also improve the inversion precision.     

Transmitting boundary for water-saturated transversely isotropic strata based on the u–U formulation

A new transmitting boundary in a cylindrical coordinate system has been developed for modeling the elastic waves radiating out to an infinite boundary in water-saturated transversely isotropic soil strata over a rigid bedrock. The saturated soil strata are assumed to consist of a porous material and modeled as a transversely isotropic two-phase medium, based on the u−U formulation. The newly developed transmitting boundary is combined with the finite-elements model of the near-field region, using the same u−U formulation, and applied to the study of the dynamics of a rigid circular foundation in porous isotropic or transversely isotropic layered strata, either fully or partly saturated with water. The verification and application examples give valuable insights into new and interesting aspects of the dynamic behavior of rigid circular foundations in fully or partly saturated two-phase ground in terms of permeability, transverse anisotropy, and ground-water table level.     

A RESPONSE OF A RANDOMLY LAYERED EARTH TO AN ELECTRIC POINT OR DIPOLE SOURCE*

A point source or a dipole source of electric current, placed on a randomly layered semi-infinite medium, produces an inhomogeneous random potential field on the surface. The variance of the random potential, normalized with reference to the normal field (that is, a field due to a point source or a dipole source on a homogeneous medium), falls off inversely as the distance from the source. The conductivity of the layers is assumed to vary randomly about a mean value (unity) such that the variations can be represented by a zero mean homogeneous random function. Further the variations are assumed to be small compared to the mean so that the first order perturbation is adequate. The analysis shows that the dipole field is more sensitive than the single pole field to the conductivity variations.     

Exact solutions to radially symmetric two-phase flow for an arbitrary diffusivity

Over the past decade there have been a variety of exact solutions developed for one-dimensional two-phase flow, however when higher dimensions are considered there is a distinct scarcity of solutions. In this paper we consider the problem of radially symmetric two-phase flow, into an infinite medium of uniform initial saturation, subject to a constant flux V from a line source at the origin. We show that in the absence of gravity and when the two-phase diffusivity D is related to the fraction flow function f by βD = V df/dθ, where θ is the water content and β is a constant of proportionality, a new class of exact solutions can be found. In particular, when β = 2, we show that the solution is given by a simple quadrature for arbitrary D, and is fully integrable for specific functional forms of D. It has been shown by Weeks et al. [Weeks SW, Sander GC, Parlange J-Y. n-Dimensional first integral and similarity solutions for two-phase flow. ANZIAM J 2003;44:365–80] that when D obeys the above relation, a saturated zone does not grow around the line of injection, consequently we find for β = 1, the flow equation maps to one-dimensional single-phase flow under a saturated boundary condition. Consequently solutions developed for one-dimensional single-phase flow (exact or approximate) apply to radially symmetric two-phase flow. Solutions for β = 1 or 2 can be derived for either a wetting fluid displacing a non-wetting fluid, or a non-wetting fluid displacing a wetting fluid, however for arbitrary β numerical methods are required.     

One-dimensional dispersion in unsteady flow in an adsorbing porous medium: An analytical solution

Analytical solutions describing the concentration distribution along one-dimensional unsteady seepage flow through adsorbing saturated finite porous medium have been obtained. An exponential function concentration is enforced at the source of the dispersion, while the change in the concentration is zero at the other boundary. A new time variable has been introduced to solve the unsteady flow problem and the solution is illustrated.     

横观各向同性层状半空间中的弹性位错

本文在柱向量函数系下,利用传播矩阵法求解了层状横观各向同性半空间由内部点源位锚引起的变形;对六个基本点源位错,以等价体力法推出了横观各向同性情形下的点源函数,并且给出了内部任意剪切位错源引起的地表位移的积分表达式。为研究地球的层状结构,特别是其上部的横观各向同性对地表的地震位移、应变以及倾斜场的影响提供了计算公式。     

Numerical simulation of spatial-temporal evolution characteristics of subsurface fluid based on strong body seismogenic model

Using finite element technique of the plane-strain problem in solid-liquid two-phase medium, we studied the characteristics of "field precursors" and "focus precursors" of subsurface fluid and their spatial-temporal evolution in case of dip-slip earthquake. The results show that: ① the change of ground fluid is slow and the anomaly is not prominent in the early period which is of elastic accumulation and non-linear; ②dilatancy emerges and anomaly focus mainly in the source region in the moderate period which is hardening and of local dilatancy. In the period the focus precursors emerge earlier than the field precursors; ③anomalies spread continuously in the source area and new regions with big anomaly emerge out of the source region in the middle-short period which is of large-scale dilatancy.     

基于反演基质矿物模量和多约束条件的双相介质AVO正演方法（英文）

基于双相介质理论的AVO正演技术是储层性质描述和流体预测的有效技术手段之一,但是输入参数中基质矿物模量的准确性和双相介质模型的的合理性极大地影响双相介质AVO正演效果。因此,本文采用基于流体因子的基质矿物模量反演方法,自适应反演基质矿物体积模量。引入具有岩石物理意义的多约束条件,使得流体替换技术制作的双相介质模型具有岩石物理意义。保证获得的双相介质AVO特征反映实际地层响应,真实可靠。通过不同岩性岩样的对比分析,说明反演方法的优越性和准确性。同时LH地区实际资料应用,获得孔隙度和流体饱和度等重要岩性参数变化时双相介质AVO特征,特别是不同储层孔隙度在同一入射角对应快纵波和横波反射系数幅值的大小差异和突变角差异是分辨储层孔隙度大小的依据。