Study of elastic wave propagation in two-phase anisotropic media by numerical modeling of pseudospectral method

When there exists anisotropy in underground media, elastic parameters of the observed coordinate possibly do not coincide with that of the natural coordinate. According to the theory that the density of potential energy, dissipating energy is independent of the coordinate, the relationship of elastic parameters between two coordinates is derived for two-phase anisotropic media. Then, pseudospectral method to solve wave equations of two-phase anisotropic media is derived. At last, we use this method to simulate wave propagation in two-phase anisotropic media, four types of waves are observed in the snapshots, i.e., fast P wave and slow P wave, fast S wave and slow S wave. Shear wave splitting, SV wave cusps and elastic wave reflection and transmission are also observed.     

含黏滞流体各向异性孔隙介质中弹性波的频散和衰减

基于Biot理论,考虑液相的黏弹性变形和固液相接触面上的相对扭转,提出了含黏滞流体VTI孔隙介质模型.从理论上推导出,在该模型中除存在快P波、慢P波、SV波、SH波以外,还将存在两种新横波－慢SV波和慢SH波.数值模拟分析了6种弹性波的相速度、衰减、液固相振幅比随孔隙度、频率的变化规律以及快P波、快SV波的衰减随流体性质、渗透率、入射角的变化规律.结果表明慢SV波和慢SH波主要在液相中传播,高频高孔隙度时,速度较高;大角度入射时,快P波衰减表现出明显的各向异性,而快SV波的衰减则基本不变;储层纵向和横向渗透率存在差异时,快SV波衰减大的方向渗透率高.     

双相各向异性介质中弹性波方程的有限元解法及波场模拟

基于双相各向异性介质模型,首先推导了双相各向异性介质中弹性波传播的动力学方程及其Galerkin变分方程和有限元运动方程,然后给出了孔隙弹性波方程的有限元数值解法以及二维双相PTL介质中波场模拟的人为吸收边界条件. 最后,利用本文给出的有限元方法对双相PTL介质和双相各向同性介质中的弹性波传播进行了数值模拟. 结果表明:有限元方法和吸收边界条件有效、可行,在理想相界条件下,不论是从固体位移,还是从流体位移的波场快照都能看到明显的慢速拟P波;在黏滞相界情况下,能否观察到慢速拟P波,与含流体地层介质的耗散性质有关.对实际含流体介质,从流体位移分量的波场快照比从固体位移波场快照更容易观察到慢速拟P波.     

部分饱和岩石的黏弹性Chapman-Kelvin模型的正演模拟方法研究

在讨论地震波传播理论时, 大部分情况下是把地震波看作弹性波.事实上, 地下介质是非完全弹性介质, 这使得传统的均匀、完全弹性介质理论受到了严重挑战.地震波在裂隙岩石中传播受裂隙系统和流体含量的影响较大, 以往许多关于频变各向异性的理论大多局限于单相流体假设, 但是几乎所有储层中通常被一种以上的流体部分饱和或完全饱和.本文在Chapman理论基础上, 提出了基于部分饱和的黏弹性Chapman-Kelvin模型的正演模拟新方法, 以提高对部分饱和岩石的黏弹性地震波频变各向异性的认识.该方法基于Chapman和Kelvin理论模型, 计算了双相不混溶流体饱和裂隙岩石中的黏弹性波频变各向异性弹性系数的表达式, 提出包含喷射流和斑块效应的统一地震波传播的黏弹性Chapman-Kelvin新模型.通过对新模型进行数值试验, 讨论了无裂缝和存在裂缝两种情况下, 含裂隙储层部分饱和岩石中耦合的喷射流和斑块效应对黏弹性介质频变各向异性的影响.试验结果反映出黏弹性介质下地震波频变各向异性的变化规律, 验证了本文提出的新方法和新模型正确.本文将黏弹性介质各向异性与裂缝中的流体流动参数相联系, 有利于提高对含裂缝储层部分饱和岩石的黏弹性波频变地震各向异性的认识, 以及地震学与油藏工程的结合程度.

     

On the applicability of the Biot solution to the problems of elastic propagation in saturated porous media

In this paper, the solution of the system of homogeneous Biot equations, which was derived by Biot for the displacement vectors of plane monochrome elastic waves propagating in a homogeneous infinite two-phase medium, is expanded to the case where the propagation area of the elastic waves is limited and the wavefront is a piecewise smooth curved surface. It is shown that the arbitrary system of homogeneous Biot equations for the displacement vectors of the solid and liquid phases can be reduced to three different equations pertaining to the class of Helmholtz equations. From this, irrespective of the geometry of the seismic wavefront and the boundaries of the studied two-phase medium, there is the following. (1) Each displacement vector (of the solid and liquid phase) splits into three independent vectors satisfying three different Helmholtz equations. Two of these vectors correspond to the two types of compressional waves, namely, fast waves (waves of the first kind) and slow waves (waves of the second kind). The third vector describes shear waves. (2) The similar (related to the same wave type) components of the displacement vector in the solid and liquid phases satisfy the same Helmholtz equation and are linked with each other through a corresponding scalar factor that is expressed in terms of the coefficients of the Biot equations. Taking into account the established properties of the displacement vectors in the solid and liquid phases seems to be helpful in the problems dealing with calculation of elastic fields of arbitrary sources in piecewise-homogeneous two-phase media.     

基于GSLS模型TI介质衰减拟声波方程

随着计算机硬件技术的发展以及高分辨率勘探需求的增加,我们希望能够更准确地模拟地下介质,得到更丰富的地层信息.然而,传统的声学假设并不能描述实际地层所存在各向异性和黏滞性,使得成像分辨率较低.为了实现深部储层的高精度成像,本文同时考虑了介质的各向异性和黏滞性,从TI介质弹性波的基本理论出发,结合各向异性GSLS理论,并通过声学近似方法导出基于GSLS模型的各向异性衰减拟声波方程.数值模拟表明该方程既能准确地描述各向异性介质下的准P波运动学规律,又能体现地层的吸收衰减效应;模型逆时偏移结果表明,在实现成像过程中考虑各向异性和黏滞性的影响,能对高陡构造清晰成像,且剖面振幅相对均衡,分辨率较高.     

双相各向异性研究、问题与应用前景

地球内部流体的存在和岩石各向异性是地下介质的两大表征,考虑地下流体和介质各向异性问题的双相各向异性理论是当今地震学和地球物理学理论与应用研究的前沿和难题之一。孔隙流体的存在、固体和流体之间的相互作用会弱化或硬化岩石的力学属性,上起声波或弹性声波速度的频散的振幅的衰减,并产生第二类压缩波。裂缝或裂隙的定向分布、岩层的旋性沉积、应力场的定向排列,都会引起传播速度的各向异性、横波发生分裂等重要现象,这些     

各向异性介质中地震波前面的偏微分方程

从含２１个弹性参数的各向异性介质中关于位移分量ｕ_ｘ、ｕ_ｙ与ｕ_ｚ的偏微分波动方程组出发，通过假定平面波位移函数解，导出准Ｐ波、准ＳＶ波与准ＳＨ波的波前面偏微分控制方程，进而对各类特殊各向异性介质（横向各向同性介质、椭圆及立方体各向异性介质）中地震波前面偏微分方程进行了讨论．以上结果为研究各向异性介质中地震波传播规律以及进行正、反演研究奠定了理论基础．     

各向异性介质qP波传播描述I：伪纯模式波动方程

地球介质相对于地震波波长尺度的定向非均匀性会导致波速的各向异性,进而影响地震波场的运动学与动力学特征．各向异性弹性波动方程是描述该类介质波场传播的基本工具,在正演模拟、偏移成像与参数反演中起着关键作用．为了面向实际应用构建灵活、简便的各向异性波场传播算子,人们一直在寻求简化的各向异性波动方程．本文借鉴各向异性弹性波波型分离思想,通过对平面波形式的弹性波方程(即Christoffel方程)实施一种代表向波矢量方向投影的相似变换,推导出了一种适应任意各向异性介质、运动学上与原始弹性波方程完全等价,在动力学上突出qP波的新方程,即qP波伪纯模式波动方程．文中以横向各向同性(TI)介质为例,给出了相应的qP波伪纯模式波动方程及其声学与各向同性近似,并在此基础上开展了正演模拟和逆时偏移试验,展示了这种描述各向异性波场传播的新方程的特点与优势．     

双相横向各向同性介质分界面上弹性波反射与透射问题研究

从双相横向各向同性介质弹性波波动方程出发 ,利用边界上的 4个连续性条件 ,计算双相横向各向同性介质分界面上弹性波反射和透射系数 .计算表明 ,快纵波在双相横向各向同性介质分界面上 ,要产生反射快纵波、反射转换 qSV波、反射转换慢纵波和透射快纵波、透射转换qSV波、透射转换慢纵波 .反射转换慢纵波振幅和透射转换慢纵波振幅均较小 .频率、耗散和各向异性大小影响着快纵波反射系数的大小 .     

Effects of petrophysical parameters on attenuation and dispersion of seismic waves in the simplified poroelastic theory: Comparison between the Biot's theory and viscosity‐extended theory

The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly.     

垂向裂隙各向异性煤层地震波响应

煤层中存在的裂隙会导致介质表现为各向异性,本文以HTI型煤层为例,结合各向异性介质弹性矩阵和各向异性裂隙理论,推导出不同充填物的垂直裂隙中各向异性参数表达式,将其应用于地震波响应分析;通过改进的交错网格差分法和各向异性Christoffel方程波场分解法,得到地震波合成记录和分解后的P波和SV波记录;将Thomsen群速度与相速度公式,经过坐标轴旋转变换,得到HTI型煤层中不同各向异性参数的地震波速度响应表达式;建立不同类型煤层地质模型,分析了裂隙密度、裂隙充填物以及煤层厚度等参数变化时的地震波响应特征.研究结果为分析垂向裂隙各向异性薄煤层地震波传播规律提供工具,为选用相应地震数据进行地震波各向异性参数反演提供依据.

     

双相介质中地震波的频率-空间域数值模拟

本文利用优化的25点频率-空间域有限差分算法对基于BISQ模型双相各向同性介质中的地震波进行了数值模拟.通过与经典的Biot模型理论模拟结果进行对比,分析了Biot流动（宏观流体流动）和Squirt流动（微观流体流动）耦合作用对地震波在孔隙介质中传播特性的影响.数值模拟在地震频段进行,结果显示：在理想相界和黏滞相界情况下,Squirt流动机制都比Biot流动机制产生了更大的速度频散和能量衰减.其中,在Biot流动和Squirt流动耦合作用下的快P波的速度和振幅小于仅考虑Biot流动影响下快P波速度和振幅,而且慢P波的衰减也更加强烈.本文还研究了地震波在双层双相各向同性介质分界面处的反射和透射特征,双相介质中波的反射与透射现象类似于单相介质的情况.模拟结果表明,利用优化25点频率-空间域有限差分法模拟双相孔隙介质中的地震波场是可行的,这为开展双相孔隙介质全波形反演问题的研究提供了可能.     

地震各向异性流体检测技术在胜利油田垦71井区的应用(英文)

根据Chapman理论模型,在各向异性介质(如HTI介质)中,当入射角在0-45。范围内,慢横波会发生较大的衰减和频散,且对流体粘度敏感,而P波和快横波则比较小。对于沿裂隙法向传播的慢横波,其振幅受流体影响很大。因此,在P波响应对流体不敏感的情况下,可利用慢横波来获得裂隙型油气藏的流体信息。本文分析了胜利油田垦71地区三维三分量地震数据,检测出的慢横波振幅和旅行时异常与该区的测井资料十分吻合。分析结果还发现,与含油区相比,含水区会产生更高的横波分裂。在含水区,慢横波振幅会产生明显变化,而在含油区则几乎没有变化。     

Seismoelectric numerical simulation in 2D vertical transverse isotropic poroelastic medium

Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ.     

求解3D弹性波方程的并行WNAD方法及其TI介质中的波场模拟

准确模拟TTI介质中弹性波的传播是研究地震各向异性、AVO反演的基础. 在二维加权近似解析离散化(WNAD)算法的基础上, 本文发展的并行WNAD算法是一种研究三维横向各向同性(TI)介质中弹性波传播的、快速高效的数值模拟方法. 我们首先介绍三维WNAD方法的构造过程, 然后与经典的差分格式——交错网格(SG)算法进行了比较. 理论分析和数值算例表明, WNAD算法比交错网格算法更适合在高性能计算机上进行大规模弹性波场模拟. 同时, 本文利用并行的WNAD方法研究了弹性波在TTI介质中的传播规律, 观测了TI介质中弹性波传播的重要特征:横波分离、体波耦合和速度各向异性等. 在TTI介质分界面处, 弹性波产生更加复杂的折射、反射和波型转化, 使得波场非常复杂, 研究和辨别不同类型的波能够加深我们对由裂隙诱导的各向异性介质的认识.     

Analytical representation of Green’s functions of the Biot equations for elastic waves in a homogeneous two-phase medium

The system of Biot vector equations in the frequency space includes two elliptic-type vector partial differential equations with unknown displacement vectors in the solid and liquid phases. Considering the Biot equations, alongside with Pride’s equations, the key approaches to the theoretical study of the elastic waves in the two-phase fluid-saturated media, the author suggests an analytical solution for the inhomogeneous Biot equations in the frequency space, which is reduced to finding its fundamental solution (Green’s function). The solution of this problem consists of solutions for two systems of Biot equations. In the first system, only the first equation is inhomogeneous, while in the second system, only the second equation is inhomogeneous and, as it is shown, its right-hand side is exclusively a potential function. The fundamental solution of the full system of inhomogeneous Biot equations (in which both equations are inhomogeneous) is represented in the form of Green’s matrix-tensor, for the scalar elements of which the analytical relations are presented. The obtained formulas describing the elastic displacements of both the solid and liquid phases reflect three wave types, namely, compressional waves of the first and the second kind (the fast and the slow waves, respectively) and shear waves. Similar terms (those describing the same type of the elastic waves in the solid and liquid phases) in the expressions for Green’s functions are linked with each other through the coefficient that links the components of the displacement vectors of the solid and liquid phases corresponding to the given wave type.

     

Analytical representation of Green’s functions of the Biot equations for elastic waves in a homogeneous two-phase medium

The system of Biot vector equations in the frequency space includes two elliptic-type vector partial differential equations with unknown displacement vectors in the solid and liquid phases. Considering the Biot equations, alongside with Pride??s equations, the key approaches to the theoretical study of the elastic waves in the two-phase fluid-saturated media, the author suggests an analytical solution for the inhomogeneous Biot equations in the frequency space, which is reduced to finding its fundamental solution (Green??s function). The solution of this problem consists of solutions for two systems of Biot equations. In the first system, only the first equation is inhomogeneous, while in the second system, only the second equation is inhomogeneous and, as it is shown, its right-hand side is exclusively a potential function. The fundamental solution of the full system of inhomogeneous Biot equations (in which both equations are inhomogeneous) is represented in the form of Green??s matrix-tensor, for the scalar elements of which the analytical relations are presented. The obtained formulas describing the elastic displacements of both the solid and liquid phases reflect three wave types, namely, compressional waves of the first and the second kind (the fast and the slow waves, respectively) and shear waves. Similar terms (those describing the same type of the elastic waves in the solid and liquid phases) in the expressions for Green??s functions are linked with each other through the coefficient that links the components of the displacement vectors of the solid and liquid phases corresponding to the given wave type.     

Finite element equations and numerical simulation of elastic wave propagation in two—phase anisotropic media

Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotopic media are derived in this paper.Numerical solution of finite element equations is given.Finally,properties of elastic wave propagation are observed and analyzed through FEM modeling.     

双相各向异性介质中弹性波传播有限元方程及数值模拟

基于Biot双相各向异性介质理论和动态问题的哈密顿原理，推导出任意双相各向异性介质中弹性波传播的有限元方程，并给出双相各向异性介质中弹性波有限元方程的数值解法．最后进行有限元法的数值模拟，对双相各向异性介质中弹性波传播特征进行了模拟与分析． 