四重旋转对称含裂隙双相各向异性介质中的地震波

基于各向异性介质和含裂隙双相介质中的地震波传播理论,导出了四重旋转对称含裂隙双相各向异性介质的本构关系与地震波传播的运动方程,并以平面波为例,进行了初步的分析研究．      

含裂隙双相各向异性介质中的地震波传播

基于各向异性介质和含裂隙双相介质中的地震波传播理论,导出了具有四阶对称轴的含裂隙双相各向异性介质的本构关系与波传播的运动方程．指出在含裂隙双相各向异性介质中可能传播６种类型的准纵横波,即快纵波ＱＰ１,慢纵波ＱＰ２,两个分裂的快横波ＱＳＶ１,ＱＳＨ１和两个分裂的慢横波ＱＳＶ２,ＱＳＨ２．并以平面波传播为例作了进一步分析     

弹性介质各向异性研究沿革、现状与问题

各向异性介质中地震波的传播研究是当今地震学研究领域中的前沿课题之一,同时也是地震学中难题之一.由于地下岩石的各向异性主要表现在:地震波速度随传播方向发生变化;不同类型体波间相互耦合;横渡发生分裂:面波速度频散依赖于传播方向等.薄互层与裂隙定向分布等产生视各向异性,放在石油地震勘探、地震预测和岩石层物理与动力学研究中有极大潜力和应用前景,并受到广泛重视.为此,本文较为详尽地讨论了弹性介质中地震波各向异性研究的沿革,简述了国内外现今取得的主要研究成果以及目前尚存在的和有待解决的一些主要问题.最后对我国在各向异性介质中地震波动、检测和应用等研究提出了一些看法.     

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

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

Studies on seismic waves

The development of seismic wave study in China in the past four years is reviewed. The discussion is divided into several aspects, including seismic wave propagation in laterally homogeneous media, laterally heterogeneous me-dia, anisotropie and porous media, surface wave and seismic wave inversion, and seismic wave study in prospect-ing and logging problems. Important projects in the current studies on seismic wave is suggested as the develop-ment of high efficient numerical methods, and applying them to the studies of excitation and propagation of seis-mic waves in complex media and strong ground motion, which will form a foundation for refined earthquake haz-ard analysis and prediction.     

地震波研究://car

回顾了近4年来中国地球物理学家在地震波研究方面取得的主要进展. 分别简要地讨论了层状均匀介质、横向不均匀介质和多孔介质中的地震波传播，地震面波及波形反演，以及地震勘探、测井问题中的地震波研究方面取得的新成果. 指出当前地震波研究的重要课题是发展高效的数值方法，并利用该方法开展复杂介质中的地震波激发和传播，以及强地面运动的研究，为精细的地震危险分析与预测奠定基础.      

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

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

各向异性介质qP波传播描述II:分离纯模式标量波

在各向异性地震波场中,qP波与qS波常常是耦合在一起的.多分量地震数据处理中一个关键环节就是波型分离(即模式解耦),以纵波成分为主的常规单分量地震数据的成像则需要合理描述标量qP波的传播算子.本文作者曾构建了在运动学上同弹性波动方程等价,动力学上突出标量qP波的伪纯模式波动方程.为了彻底消除qS波残余,本文根据波矢量与qP波偏振矢量之间的偏差,提出从伪纯模式波场提取纯模式标量qP波的方法.数值分析展示了投影偏差算子在波数域和空间域的特征.基于不同复杂程度理论模型的试验结果表明,联合"伪纯模式传播算子"与"投影偏差校正"可为各向异性介质分离模式波场传播过程提供一种简便的描述工具.     

基于分数阶时间导数的黏弹性衰减VTI介质中平面波传播

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推导出准压缩波qP、准剪切波qSV和纯剪切波SH的复速度、相速度、能量速度以及品质因子的解析表达式.对模型的正确性进行了数值验证,并分析了qP,qSV和SH波在介质中的传播特性.数值试验结果表明:本模型能够实现理想的恒定Q行为,表现了品质因子和速度的各向异性特征,显示出黏弹性增强将导致能量速度和相速度的频散曲线变化剧烈;速度和衰减各向异性参数与传播角度之间的耦合效应对qP,qSV和SH波的速度和能量影响明显;qP,qSV和SH波的频散曲线和波前面随着衰减各向异性强度的改变发生显著变化,其中耦合在一起的qP和qSV波变化趋势相同,而SH波与它们呈现相反的变化规律.本研究为从常Q模型角度分析地震波在衰减各向异性黏弹性介质中的传播特征奠定了理论基础.     

裂缝诱导双相HTI介质地震波场错格伪谱法模拟与波场特征分析

裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.     

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 δ.     

双相介质中地震波衰减的物理机制

High-frequency seismic attenuation is conventionally attributed to anelastic absorption. In this paper, I present three studies on high-frequency seismic attenuation and propose that the physical mechanism results from the interference of elastic microscopic multiple scattering waves. First, I propose a new theory on wave propagation in a two-phase medium which is based on the concept that the basic unit for wave propagation is a nano- mass point. As a result of the elasticity variations of pore fluid and rock framework, micro multiple scattering waves would emerge at the wavelength of the seismic waves passing through the two-phase medium and their interference and overlap would generate high- frequency seismic attenuation. Second, I present a study of the frequency response of seismic transmitted waves by modeling thin-layers with thicknesses no larger than pore diameters. Results indicate that high-frequency seismic waves attenuate slightly in a near-surface water zone but decay significantly in a near-surface gas zone. Third, I analyze the seismic attenuation characteristics in near-surface water and gas zones using dual-well shots in the Songliao Basin, and demonstrate that the high-frequency seismic waves attenuate slightly in water zones but in gas zones the 160-1600 Hz propagating waves decay significantly. The seismic attenuation characteristics from field observations coincide with the modeling results. Conclusions drawn from these studies theoretically support seismic attenuation recovery.     

Propagation of seismic waves in a viscoelastic medium and the effect of topography

In this paper we have simulated the propagation of seismic waves in viscoelastic medium and calculated the effect of irregular surface topography by using finite element method. Several types of elastic and viscoelastic medium models, such as block structure, vidges and/or valley with and without a soft layer underneath have been studied. The distribution of maximum amplitude of displacement and acceleration on surface has been calculated in the case of vertically incident SH and P waves from the basement. The possible application in earthquake engineering is also discussed. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 442–449, 1991.     

利用介质的不均匀分布构建各向异性模型

在复杂区域应力场的作用下, 大范围扩容各向异性(extensive-dilatancy anisotropy, EDA)造成的裂隙定向排列形成了地壳内部介质的不均匀性。 因此, 地壳介质各向异性的描述不仅限于横向各向同性(transversely isotropic, TI)的理论。 我们首先利用自相关函数随机扰动得到水平向小尺度速度不均匀, 然后利用这种不均匀性描述地壳介质中EDA裂隙定向排列所造成的介质各向异性, 利用数值模拟弹性波在该介质中的传播过程验证了利用介质的不均匀性构建介质的各向异性是一个有效、 可靠的方法。 可以将此方法应用到复杂地壳介质结构中, 了解地震波在复杂地壳各向异性介质中的传播特征。     

Q value of anelastic S-wave attenuation in Yunnan region

The study of seismic attenuation property is a major subject in seismology. Seismic waves recorded by seismic stations （seismographs） contain source effect, seismic wave propagation effect, site response of seismic stations and instrumental response. The path effect of seismic wave propagation, site response of seismic stations and instrumental response must be taken out in the study of source property with seismic data. The path effect of seismic wave propagation （seismic attenuation） involves an important influential factor, the anelastic attenuation of medium, which is measured with quality factor Q, apart from geometric attenuation with the distance. As a basic physical parameter of the Earth medium, Q value is essential for quantitative study of earthquakes and source property （e.g. determination of source parameters）, which is widely used in earthquake source physics and engineering seismology.     

Anisotropic elastic full‐waveform inversion of walkaway vertical seismic profiling data from the Arabian Gulf

Borehole seismic addresses the need for high‐resolution images and elastic parameters of the subsurface. Full‐waveform inversion of vertical seismic profile data is a promising technology with the potential to recover quantitative information about elastic properties of the medium. Full‐waveform inversion has the capability to process the entire wavefield and to address the wave propagation effects contained in the borehole data—multi‐component measurements; anisotropic effects; compressional and shear waves; and transmitted, converted, and reflected waves and multiples. Full‐waveform inversion, therefore, has the potential to provide a more accurate result compared with conventional processing methods. We present a feasibility study with results of the application of high‐frequency (up to 60 Hz) anisotropic elastic full‐waveform inversion to a walkaway vertical seismic profile data from the Arabian Gulf. Full‐waveform inversion has reproduced the majority of the wave events and recovered a geologically plausible layered model with physically meaningful values of the medium.     

正交各向异性介质中qP波入射的二阶近似反射系数与透射系数

地震波在各向异性介质中以一个准P波(qP)和两个准S波(qS1和qS2)的形式传播.研究三种波的相速度、群速度以及偏振方向等传播性质能够为各向异性介质中的正反演问题提供有效支撑.具有比横向各向同性(TI)介质更一般对称性的正交各向异性介质通常需要9个独立参数对其进行描述,这使得对传播特征的计算更为复杂.当两个准S波速度相近时具有耦合性,从而令慢度的计算产生奇异性.因此,奇异点(慢度面的鞍点和交叉点)附近的反射与透射(R/T)系数的求解不稳定,会导致波场振幅不准确.本文首次通过结合耦合S波射线理论和基于迭代的各向异性相速度与偏振矢量的高阶近似解,得到了适用于正交各向异性介质以qP波入射所产生的二阶R/T系数的计算方法.与基于一阶近似的结果相比,基于二阶近似的方法提高了qP波R/T系数的精度,能得到一阶耦合近似无法表达的准确的qP-qS转换波的R/T系数解,且方法适用于较强的各向异性介质.     

热弹性介质中波传播特征

深层-超深层油气地震勘探涉及高温介质地震波传播问题,热弹介质参数对地震波传播有重要影响.含弛豫时间修正项的Lord-Shulman双曲型耦合热弹波动方程从理论上预测了热弹性介质中存在快纵波、慢纵波(一种准静态慢纵波,简称热波)和横波的传播,两个纵波为热耗散衰减波而横波不受介质热特性的影响.本文结合平面波频散分析和格林函数法数值模拟,详细研究两个热耗散衰减波的频散和衰减特征,着重分析热导率、热膨胀系数及比热的变化对波速和衰减的影响.研究表明热导率作为主要参数决定了波速与衰减的临界变化,热膨胀系数对波速和衰减的幅度有明显影响,比热则兼顾了前两个热弹系数的影响特征.最后,利用热弹性动力学频率域的二阶格林函数进行波场快照数值模拟,展示热弹性介质中纵波、横波和热波的传播行为.     

Propagation characteristics of converted refracted wave and its application in static correction of converted wave

Three-component seismic exploration through P-wave source and three-component geophone is an effective technique used in complicated reservoir exploration. In three-component seismic exploration data processing,one of the difficulties is static correction of converted wave. This paper analyzes propagation characteristics of non-converted and converted refracted waves,and discovers a favor-able condition for the formation of converted refracted wave,i.e. the velocity of overlaying medium S wave is much lower than that of underlying medium S wave. In addition,the paper proposes the static correction method of converted wave based on PPS converted refracted wave,and processes the real three-component seismic data with better results of static correction of converted wave.     

双相各向异性介质中偶数阶精度有限差分数值模拟

To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.