共查询到19条相似文献，搜索用时 140 毫秒
1.
裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征. 相似文献
2.
3.
伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值算法,由于它的微分算子的全局性,一般认为该方法不适于并行计算. 本文介绍了并行计算非均匀介质中地震波传播的重叠区域分解算法,给出了一种基于傅里叶伪谱法的并行算法. 文中给出的算法将介质划分为相互重叠的若干区域,在各个子域上单独求解,利用重叠部分的解的传递,将各个子域连接起来,实现了伪谱法在分布式并行处理机上的计算. 文中给出了一个将二维区域分解的算例,比较了并行算法和整体算法的结果,分析了并行算法的计算精度. 结果表明,并行算法会有效降低计算时间,并且保证计算精度. 该方法在大规模三维非均匀介质的地震波场模拟方面有应用价值. 相似文献
4.
油气勘探开发实践证明,裂缝常常是油气藏存储的空间或运移的通道,因此,裂缝各向异性介质中地震波场的研究越来越倍受关注,国内外很多岩石物理学者、地球物理专家等对裂缝信息的描述提出了很多理论认识与方法技术.本文根据Eshelby-Cheng各向异性裂缝介质模型理论,求取各向异性裂缝介质的弹性参数,并建立Eshelby-Cheng各向异性裂缝介质的波动方程,利用时间错格伪谱法对含流体裂缝介质进行数值模拟,模拟结果表明,采用时间错格伪谱法能有效解决各向异性介质的波场传播,利用时间错格有限差分算子替代普通的差分算子来求解时间导数,利用快速傅氏变换求解空间导数,大大提高了正演模拟的计算精度与计算效率.并且与各向同性介质相比,地震波在含流体裂缝各向异性介质中的传播要复杂得多,各向同性介质层中的波是纯的,其横波不会发生分裂,而在各向异性介质层中,横波将发生分裂. 相似文献
5.
6.
7.
横向各向同性介质是地球内部广泛存在的一种各向异性介质，因此为了能够更好地认识地震波在这种介质中的传播特征，用数值方法进行地震波模拟显得十分必要.本文采用谱元法对横向同性介质中的地震波进行模拟，该方法基于弹性力学方程弱形式基础之上，具有有限元适应任意复杂介质模型的韧性和伪谱法的精度.文中阐述了基于Legendre多项式的谱元法的理论和推导过程，该方法可以形成全局对角质量矩阵，在时间域使用显式的差分算法，提高运算效率，最后通过横向各向同性介质的数值计算,模拟结果表明该方法是一种有效的数值模拟方法. 相似文献
8.
伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值方法,但是由于它的微分算子的全局性,使得该方法不适用于分散内存的并行计算.本文将有限差分算子的局部性和伪谱法算子的高效、高精度相结合,发展基于两种方法的伪谱/有限差分混合方法.该方法在一个空间坐标方向上利用交错网格高阶有限差分算子,在另外的空间坐标方向上利用交错网格伪谱法算子,既保留了后者的高效、高精度优势,又便于在PC集群上实现并行计算.对二维模型的计算显示,混合方法能有效处理介质不连续面,在保证伪谱法计算精度的情况下,提供了一种并行计算的可能途径. 相似文献
9.
本文导出了一种由单程波方程利用反射／透射算子的可分表示方法模拟复杂介质中一次反射地震波的数值算法. 文中利用算子可分表示理论将反射／透射算子分解成适合于双域（空间域和波数域）运算的表达形式，使得本文得到的地震波数值模拟算法可适应于一定程度横向非均匀介质和界面起伏情况，在入射角小于45°时能够准确模拟振幅随入射角(AVA)的关系. 就模拟一次反射地震波而言，与前人研究的双程波动方程伪谱法地震波模拟相比，本文算法具有足够高的模拟精度，且计算效率成倍地提高. 相似文献
10.
弹性波在离散介质中传播的贝塞尔函数解Ⅰ 总被引：3，自引：0，他引：3
对弹性波在离散介质中的传播问题进行了新的探索研究,本文用贝塞尔函数的线性组合构造了弹性波在非均匀离散介质中传播解。借助于贝塞尔函数的良好性质,其研究结果可用于地震波传播问题的研究。 相似文献
11.
12.
Orthorhombic models are often used in the seismic industry nowadays to describe azimuthal and polar anisotropy and reasonably realistic in capturing the features of the earth interior. It is challenging to handle so many model parameters in the seismic data processing. In order to reduce the number of the parameters for P wave, the acoustic orthorhombic medium is proposed by setting all on-axis S wave velocities to zero. However, due to the coupled behaviour for P and S waves in the orthorhombic model, the ‘S wave artefacts’ are still remained in the acoustic orthorhombic model, which kinematics needs to be defined and analysed. In this paper, we analyse the behaviour of S wave in acoustic orthorhombic media. By analysis of the slowness surface in acoustic orthorhombic media, we define the S waves (or S wave artefacts) that are more complicated in shape comparing to the one propagating in an acoustic transversely isotropic medium with a vertical symmetry axis. The kinematic properties of these waves are defined and analysed in both phase and group domain. The caustics, amplitude and the multi-layered case for S wave in acoustic orthorhombic model are also discussed. It is shown that there are two waves propagating in this acoustic orthorhombic medium. One of these waves is similar to the one propagating in acoustic vertical symmetry axis media, whereas another one has a very complicated shape consisting of two crossing surfaces. 相似文献
13.
In this work, we develop the indirect boundary element method (IBEM) to simulate the seismic site response in a realistic, large-scale 3-D sedimentary basin. Most previous applications of boundary element method have used full-space Green’s functions for wave propagation between element points. We use half-space Green’s functions, which include the seismic wavefield interactions at the free surface and require only the boundary elements of the basin interface. In this way, the size of the matrix equation for solution in the IBEM can be reduced to approximately a quarter of that using full-space Green’s functions. The site response modeling of the Granada basin in southern Spain using the IBEM shows that the basin-induced scattering waves were identified as propagating back and forth inside the basin. The scattered waves also generate surface waves that are weakly propagated outside of the basin. The wave propagation inside and outside of the basin shows different patterns. We observe that the scattered wave is locally amplified, and its propagation direction deviates from that of the incident waves propagation direction. Therefore, the computed seismic response in the basin could provide us with good estimates of the seismic motion. 相似文献
14.
15.
《Soil Dynamics and Earthquake Engineering》1998,17(7-8):509-517
A new method of analyzing seismic response of soil-building systems is introduced. The method is based on the discrete-time formulation of wave propagation in layered media for vertically propagating plane shear waves. Buildings are modeled as an extension of the layered soil media by assuming that each story in the building is another layer. The seismic response is expressed in terms of wave travel times between the layers, and the wave reflection and transmission coefficients at layer interfaces. The calculation of the response is reduced to a pair of simple finite-difference equations for each layer, which are solved recursively starting from the bedrock. Compared with commonly used vibration formulation, the wave propagation formulation provides several advantages, including the ability to incorporate soil layers, simplicity of the calculations, improved accuracy in modeling the mass and damping, and better tools for system identification and damage detection. 相似文献
16.
Ray path of head waves with irregular interfaces 总被引：1，自引：0，他引：1
Head waves are usually considered to be the refracted waves propagating along flat interfaces with an underlying higher velocity.
However, the path that the rays travel along in media with irregular interfaces is not clear. Here we study the problem by
simulation using a new approach of the spectral-element method with some overlapped elements (SEMO) that can accurately evaluate
waves traveling along an irregular interface. Consequently, the head waves are separated from interface waves by a time window.
Thus, their energy and arrival time changes can be analyzed independently. These analyses demonstrate that, contrary to the
case for head waves propagating along a flat interface, there are two mechanisms for head waves traveling along an irregular
interface: a refraction mechanism and transmission mechanism. That is, the head waves may be refracted waves propagating along
the interface or transmitted waves induced by the waves propagating in the higher-velocity media. Such knowledge will be helpful
in constructing a more accurate inversion method, such as head wave travel-time tomography, and in obtaining a more accurate
model of subsurface structure which is very important for understanding the formation mechanism of some special areas, such
as the Tibetan Plateau. 相似文献
17.
Modified image algorithm to simulate seismic channel waves in 3D tunnel model with rugged free surfaces 总被引：1，自引：0，他引：1 
下载免费PDF全文
下载免费PDF全文 The existence of rugged free‐surface three‐dimensional tunnel conditions in the coal seams, caused either by geological or mining processes, will inevitably influence wave propagation characteristics when the seismic waves go through the coal mines. Thus, a modified image algorithm has been developed to account for seismic channel waves propagating through this complicated topography with irregular free surfaces. Moreover, the seismic channel waves commonly exhibit damped and dispersive signatures, which is not only because of their own unique sandwich geometry of rock–coal–rock but also because of the viscoelastic behavior of coal. Considering the complexity of programming in three‐dimensional tunnel models with rugged free surfaces, an optimized vacuum grid search algorithm, enabling to model highly irregular topography and to compute efficiently, is also proposed when using high‐order staggered finite‐difference scheme to simulate seismic channel wave propagations in viscoelastic media. The numerical simulations are implemented to investigate the accuracy and stability of the method and the impact of coal's viscoelastic behavior on seismic channel wave propagation characteristics. The results indicate that the automatic vacuum grid search algorithm can be easily merged into high‐order staggered finite‐difference scheme, which can efficiently be applied to calculate three‐dimensional tunnel models with rugged free surfaces in the viscoelastic media. The simulation also suggests that the occurrence of a three‐dimensional tunnel with free surfaces has a remarkable influence on the seismic channel wave propagation characteristics and elastic energy distribution. 相似文献
18.
本文从地球介质的复杂性与非均匀性的普遍性以及地震波在其中的传播过程的复杂性出发,对地震波散射理论及方法研究的历史进行了简要的回顾,并对该领域从20世纪90年代以来的某些进展作了简要的介绍与评述.本文分两大部分进行阐述.(1)．历史回顾：简介并回顾了地震波散射理论几十年的发展历程、存在的问题及对解决这些问题的潜在的基本途径.(2)．20世纪90年代以来的若干进展：简要介绍了弱散射理论体系的形成及强散射问题研究的部分成果,提出了散射理论未来发展的可能方向之一,即在无近似的拟微局部分析理论的基础上,建立并发展可处理高度复杂的地球固体介质中地震波传播问题的地震波散射及全波理论. 相似文献
19.
We have modeled the effect of a direct current (DC) electric field on the propagation of seismic waves by the pseudospectral
time domain (PSTD) method, based on a set of governing equations for the poroelastic media. This study belongs to the more
general term of the seismoelectric coupling effect. The set of physical equations consists of the poroelastodynamic equations
for the seismic waves and the Maxwell's equations for the electromagnetic waves; the magnitude of the seismoelectric coupling
effect is characterized by the charge density, the electric conductivity, the Onsager coefficient, a function of the dielectric
permittivity, the fluid viscosity, and the zeta potential. The poroelastodynamic vibration of a solid matrix generates an
electric oscillation with the form of streaming current via the fluctuation of pore pressure. Meanwhile, fluctuating pore
pressure also causes oscillatory variation of the electric resistivity of the solid matrix. The simulated poroelastic wave
propagation and electric field variation with an existing background DC electric field are compared with the results of a
physical experiment carried out in an oilfield. The results show that the DC electric field can significantly affect the propagating
elastic energy through the seismoelectric coupling in a wide range of the seismic frequency band. 相似文献