地震波波动方程数值模拟方法
地震波波动方程数值模拟方法
地震波波动方程数值模拟方法主要包括克希霍夫积分法、傅里叶变换法、有限元法和有限差分法等。
克希霍夫积分法引入射线追踪过程，本质上是波动方程积分解的一个数值计算，在某种程度上相当于绕射叠加。该方法计算速度较快，但由于射线追踪中存在着诸如焦散、多重路径等问题，故其一般只能适合于较简单的模型，难以模拟复杂地层的波场信息。
傅里叶变换法是利用空间的全部信息对波场函数进行三角函数插值，能更加精确地模拟地震波的传播规律，同时，利用快速傅里叶变换(FFT)进行计算，还可以提高运算效率，其主要优点是精度高，占用内存小，但缺点是计算速度较慢，对模型的适用性差，尤其是不适应于速度横向变化剧烈的模型.
波动方程有限元法的做法是：将变分法用于单元分析，得到单元矩阵，然后将单元矩阵总体求和得到总体矩阵，最后求解总体矩阵得到波动方程的数值解；其主要优点是理论上可适宜于任意地质体形态的模型，保证复杂地层形态模拟的逼真性，达到很高的计算精度，但有限元法的主要问题是占用内存和运算量均较大，不适用于大规模模拟，因此该方法在地震波勘探中尚未得到广泛地应用。。
相对于上述几种方法，有限差分法是一种更为快速有效的方法。虽然其精度比不上有限元法，但因其具有计算速度快，占用内存较小的优点，在地震学界受到广泛的重视与应用。
声波方程的有限差分法数值模拟
对于二维速度-深度模型，地下介质中地震波的传播规律可以近似地用声波方程描述：
)()(2222222tSzuxuvtu （4-1）
(,)vxz是介质在点（x , z）处的纵波速度，u为描述速度位或者压力的波场，)(ts为震源函数。
为求式（4-1）的数值解，必须将此式离散化，即用有限差分来逼近导数，用差商代替微商。为此，先把空间模型网格化（如图4-1所示）。
设x、z方向的网格间隔长度为h，t为时间采样步长，则有： z,ij1,ij2,ij1,ijhix （i为正整数）
hjz (j为正整数)tnt (n为正整数)
kjiu, 表示在(i,j)点，k时刻的波场值。
将1,kjiu在(i,j)点k时刻用Taylor展式展开：
)(*21*22*22*,1,tottuttuuutkttktkjikji (4-2)
将1,kjiu在(i,j)点k时刻用Taylor展式展开：
)(*21*22*22*,1,tottuttuuutkttktkjikji (4-3)
1）将上两式相加，略去高阶小量，整理得(i,j)点k时刻的二阶时间微商为：
21,,1,222tuuutukjikjikji (4-4)
同理可得(i,j)点k时刻的二阶空间微商分别为：
(4-5)
(4-6)
这就实现了用网个点波场值的差商代替了偏微分方程的微商，将上三个式子代入(4-1)式中得：
(4-7)
2）采用四阶精度差分格式，（以X方向为例）即将1,2kjiu、1,1kjiu、1,1kjiu、1,1kjiu分别在(i,j)点k时刻展开到四阶小量，消除四阶小量并解出二阶微分得：
}25][34][121{1,,1,1,2,2222kjikjikjikjikjiuuuuuxxu (4-8)
}25][34][121{1,1,1,2,2,222kjikjikjikjikjiuuuuuzzu (4-9)
这就实现了用网个点波场值的差商代替了偏微分方程的微商，代入(4-1)式中得可得时间二阶、空间四阶精度的声波方程差分格式为：：
}25][34][121{2,,1,1,2,22221,,1,kjikjikjikjikjikjikjikjiuuuuuhtvuuu 2,1,,1222huuuxukjikjikji21,,1,222huuuzukjikjikji)(*)(*)()42()(*)/*(001,,22221,1,,1,121,jjiiksuuhtvhuuuuhtvukjikjikjikjikjikjikji kjikjikjikjikjiuuuuuhtv,,1,1,2,222225][34][121{}
)(**)(*)(00jjiits （4-10）
式中),(jiv为介质速度的空间离散值，h是空间离散步长，t为时间离散步长，)(ks为震源函数，关于)(ks一般使用一个理论的雷克型子波代替，即：
fttftse2cos)/2()(22 （4-11）
上式中，t为时间， f为中心频率，一般取为20-40HZ，为控制频带宽度的参数，一般取3-5。在实际计算过程中，需把此震源函数离散，参与波场计算。)(*)(00jjii确定震源位置。
稳定性条件：对于特定的偏微分方程只有特定的几种有限差分格式是无条件或有条件稳定的，（4-7）、（4-10）式即是已被证明的有条件稳定格式，其稳定性条件分别为：
83/*maxhtv (4-12)
这里maxv表示的是地下介质的最大波速；若地下介质网格间隔、最小速度、及时间采样间隔不符合(4-12)式时，递推求解（4-10）式，波场值会出现误差（高阶小量）累积，出现不稳定现象。
频散关系式： )/(minNGfvh (4-13)
同时，在差分计算过程中，如果空间和时间采样间隔不当，就会导致波形畸变，甚至派生出多个同相轴，这种现象称为频散现象。
偏微分方程本身没有频散，网格频散是由于差分方程近似替代微分方程引起的。当波场按照波动方程所表示的微分方程传播时，波场的传播速度就是波动方程中的速度，但当波场按照波动方程离散化后的差分方程传播时，波场的传播速度就不再是波动方程中的速度了，而是与波的频率和波数有关的函数，具有不同频率和波数的波有不同的传播速度，因而在传播过程中会出现频散，发生畸变，且随走时的增加而增加。
式中minv为最小速度，Nf为Nyquist频率。一般取震源子波中的主频f的2倍值参与计算，G为每个波长所占的网格点数，对于空间二阶差分、时间二阶差分G取8，而对于空间为四阶差分的情况则G取4方能有效减少频散。
边界条件：在地震波场正演模拟中，必须引入人工边界来界定计算区域。人工边界若不做特殊处理，就会随着波场的递推计算在边界上产生虚假反射波从而扰乱波场，人工边界的处理是地震波场正演数值模拟的一个重要课题。 2/2/*maxhtv
2.5维非均匀介质中的地震波数值模拟
第２ ３卷 第４期
２００８年８月（页码：１０９２～１０９８） 地球物理学进展
ＰＲＯＧＲＥＳＳ ＩＮ ＧＥＯＰＨ ＹＳＩＣＳ Ｖｏ【．２３ Ｎｏ．４
Ａｕｇ． ２００８
２．５维非均匀介质中的地震波数值模拟
武 晔ｈ ， 李小凡 ， 顾观文。， 龙桂华 ， 赵晓燕 ， 武巴特尔。
（１．中科院地质与地球物理研究所，北京１０００２９； ２．防灾科技学院，燕郊１０１６ｏｌ； ３．中国地质科学院物化探研究所，廊坊０６５０００）
摘要 目前，完全的３Ｄ弹性波数值模拟计算仍然需要庞大的计算资源，而２Ｄ弹性波数值模拟，又不能准确地近似
３Ｄ数值模拟．因此，本文在较高数值精度的一阶应力～速度弹性波动方程的基础上，采用时间错格差分拟谱方法在
２Ｄ模型介质中计算３Ｄ弹性波场，即２．５Ｄ数值模拟．并将３Ｄ计算结果与２．５Ｄ和２Ｄ结果进行比较，发现２．５Ｄ模拟
结果在波场形态和振幅方面比２Ｄ模拟结果更接近３Ｄ模拟结果，从而证实了２．５Ｄ数值模拟比２Ｄ数值模拟可以更
好地近似３Ｄ波场数值模拟，同时证实了时间错格差分拟谱方法是一种高精度、高效率的正演模拟方法．
关键词 弹性波，傅立叶变换，数值模拟，拟谱，２．５维，非均匀介质
中图分类号Ｐ３１５ 文献标识码 Ａ 文章编号 １００４—２９０３（２００８）０４—１０９２—０７
２．５Ｄ Ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｓｅｉｓｍｉｃ ｗａｖｅ ｉｎ ｉｎｈｏｍｏｇｅｎｅｏｕｓ ｍｅｄｉａ
ＷＵ Ｙｅ ～， ＬＩ Ｘｉａｏ—ｆａｎ ， ＧＵ Ｇｕａｎ—Ｗｅｎ。， ＬＯＮＧ Ｇｕｉ—Ｈｕａ ，
ＺＨＡＯ Ｘｉａｏ—ｙａｎ ， ＷＵ Ｂａ—ｔｅ—ｅｒ
（１．Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｇｅｏｌｏｇｙ ａｎｄ Ｇｅｏｐｈｙｓｉｃｓ Ｃｈｉｎｅｓｅ Ａｃａｄｅｍｙ ｏｆ Ｓｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ １０００２９，Ｃｈｉｎａ； ２．Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｄｉｓａｓｔｅｒ Ｐｒｅｖｅｎｔｉｏｎ Ｓｃｉｅｎｃｅ Ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ，Ｙａｎｊｉａｏ １０１６０ｌ，Ｃｈｉｎａ； ３．Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｇｅｏｐｈｙｓｉｃａｌ＆Ｇｅｏｃｈｅｍｉｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ，Ｌａｎｇｆａｎｇ １０２８００，Ｃｈｉｎａ）
可变网格多尺度有限差分模拟方法
可变网格多尺度有限差分模拟方法
摘要院针对典型各项同性地质模型，地震波波动方程和有限差分法对交错网格地震波场进行模拟仿真研究。为了得到一阶速度应力弹性波方程，首先对地震波波动方程原始形式进行推演，然后用时空二元高阶交错差分网格的方法对其进行离散化，最后对地质模型模拟仿真实验，并获取了对应的波场快照。最后由仿真结果表明，算法计算精度准确，计算效率高，而且可变网格的变换对地质模型的研究起到了很重要的作用。
Abstract: Aiming at the typical gay geological model, the seismic wave equation
and finite difference method are used to do simulationresearch to the staggered grid
seismic wave field. In order to get a first order rate of seismic wave equation, the
original form of the seismicwave equation is deduced firstly. Then space-time dual
higher order difference staggered grid method is used to do discretization. Finallythe
geological model simulation experiment is done and the corresponding wave field
snapshots are obtained. The simulation results showthat the algorithm is accurate with
high computational efficiency, and the varision of variable grid is of important
基于波动方程的地震波形反演与成像方法研究
基于波动方程的地震波形反演与成像方法研究
地震波形反演与成像是地球物理学中重要的领域，它通过分析地震波形数据来研究地下的地质结构和介质参数。这项研究对于地质勘探、地震灾害预测、地震工程以及地球内部结构的理解具有重要意义。基于波动方程的地震波形反演与成像方法可以提供更准确的地下模型和地震源的参数。
1. 地震波形反演方法
地震波形反演是通过分析地震波形数据，推测地下的地质结构和介质参数。波形反演方法有许多种，其中最常用的是基于波动方程的全波形反演方法。全波形反演方法通过求解正问题和反问题来估计地下介质模型。正问题是根据已知的地下介质模型和地震源参数，计算出模拟地震波数据。反问题是根据观测到的地震波数据，反推估计出地下介质模型参数。 全波形反演方法是一种迭代方法，通过多次迭代求解正问题和反问题来逐步优化地下模型的估计值。在正问题的求解中，需要使用波动方程模拟地震波传播过程。波动方程是描述地震波传播的基本方程，它是一个偏微分方程，可以通过数值方法（如有限差分法、有限元法等）进行求解。在反问题的求解中，需要使用优化算法进行参数估计，最常用的方法包括共轭梯度法、拟牛顿法等。
2. 地震波形成像方法
地震波形成像是通过分析地震波形数据，进行地下介质的成像。它与波形反演方法类似，但是波形成像方法更注重于地下结构的成像，而不太关注参数估计。波形成像方法有许多种，常用的方法包括偏移成像、反射成像和散射成像。
偏移成像是一种常用的波形成像方法，它利用地震波的走时信息来定位地下结构。在偏移成像中，首先需要进行资料处理，包括去噪、去除仪器响应等。然后根据速度模型对地震波数据进行偏移处理，得到反射面在地下的位置。偏移成像的优点是处理速度快，适用于大规模数据。但是它对速度模型的准确性要求较高。
反射成像是一种基于地震波反射的成像方法。它通过分析地震波在地下发生反射的位置和特征，来推测地下的反射面。反射成像常用的方法有叠前偏移和叠后偏移等。叠前偏移是对地震道数据进行偏移处理，得到反射面的位置和形状。叠后偏移是在反射波坐标之后，对地震道数据进行偏移处理。反射成像的优点是可以提供较高分辨率的成像结果，但是对噪音和多路径效应比较敏感。
起伏地表复杂介质波动方程有限元数值模拟方法
第２ ２卷 第２期
２００７年４月（页码：５２２－－５２９） 地球物理学进展
ＰＲ０ＧＲＥＳＳ １Ｎ ＧＥ０ＰＨＹＳ１ＣＳ ＶｏＩ．２２ Ｎｏ．２
Ａｐｒ． ２００７
起伏地表复杂介质波动方程有限元数值模拟方法
薛东川 ， 王尚旭 ， 焦淑静。
（１．中国石油大学（北京）ＣＮＰＣ物探重点实验室，北京１０２２４９；２．中国石油大学（北京）机电工程学院，北京１０２２４９）
摘要 波动方程数值模拟是深入研究地震波传播规律的有效方法．有限差分法因其方法简单、精度高而得到了广
泛的应用．但其缺点是不能准确模拟具有复杂几何形态的物性界面．因而当遇到起伏地表或复杂构造时，求解精度 低．为了准确模拟起伏地形、复杂构造和复杂介质条件下的地震波场，本文采用有限元法模拟二维声波方程．用三角
形单元模拟地形和速度界面；把单元内的场和波速均看作单元上的线性函数，以适应复杂介质压制边角散射；采用吸
收边界条件去除来自截断边界上的反射；采用集中质量矩阵和集中阻尼矩阵使得显式时间递推无需对矩阵求逆，提 高了计算效率．对模型的计算表明该方法正确有效．
关键词起伏地表，复杂介质，有限元，数值模拟，吸收边界条件，集中质量矩阵
中图分类号Ｐ６３１ 文献标识码 Ａ 文章编号 １００４—２９０３（２００７）０２—０５２２—０８
Ｗａｖｅ ｅｑｕａｔｉｏｎ ｆｉｎｉｔｅ－ｅｌｅｍｅｎｔ ｍｏｄｅｌｉｎｇ ｉｎｃｌｕｄｉｎｇ ｒｕｇｇｅｄ
ｔｏｐｏｇｒａｐｈｙ ａｎｄ ｃｏｍｐｌｉｃａｔｅｄ ｍｅｄｉｕｍ
ＸＵＥ Ｄｏｎｇ—ｃｈｕａｎ ，ＷＡＮＧ Ｓｈａｎｇ—ＸＵ 。ＪＩＡＯ Ｓｈｕ—ｊｉｎｇ
（１．ＣＮＰＣＧｅｏｐｈｙｓｉｃｓ Ｋｅｙ Ｌａｂ．，Ｃｈｉｎａ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｐｅｔｒｏｌｅｕｍ，Ｂｅｉｊｉｎｇ １０２２４９，Ｃｈｉｎａ： ２．Ｆａｃｕｌｔｙ ｏｆ Ｍｅｃｈａｎｉｃａｌ ａｎｄ Ｅｌｅｃｔｒｏｎｉｃ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｃｈｉｎａ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｐｅｔｒｏｌｅｕｍ，Ｂｅｉｊｉｎｇ １０２２４９，Ｃｈｉｎａ）
波动方程pml matlab
波动方程PML（Perfectly Matched Layer）是一种用于模拟波传播的有效方法，特别适用于处理波在介质边界处的反射和折射现象。而Matlab作为一种强大的数学建模和仿真软件，被广泛应用于波动方程的数值求解和仿真研究中。在本文中，我将深入探讨波动方程PML在Matlab中的应用，旨在帮助您更全面地理解这一主题。
1. 波动方程的基本概念
波动方程是描述波动现象的数学模型，广泛应用于声波、电磁波和地震波等领域。其一般形式为：
\[ \nabla^2 p - \dfrac{1}{c^2} \dfrac{\partial^2 p}{\partial t^2} =
S \]
其中，\(p\)表示波场的物理量，\(c\)为介质中的波速，\(S\)表示波源的项。波动方程的数值求解需要考虑介质边界的影响，而PML方法正是为了处理波在边界处的吸收和反射而被提出。
2. PML方法的原理与应用
PML方法是通过在模拟区域周围引入特殊的吸收层，使得波在边界处能够被有效吸收，从而减小边界效应对波场的影响。PML方法在波动方程的数值模拟中起着至关重要的作用，能够更真实地模拟波在复杂介质中的传播过程。
在Matlab中，可以通过编写PML吸收边界条件的自定义函数，结合有限差分或有限元等数值方法，实现对波动方程PML的数值求解。Matlab提供了丰富的数学函数和绘图工具，能够方便地实现PML方法的仿真研究。
3. Matlab中波动方程PML的实现
在Matlab中，首先需要定义模拟区域的网格和介质参数，设置好波源的位置和波形类型。通过调用自定义的PML吸收边界条件函数，将吸收层的参数和边界条件添加到波动方程的数值求解过程中。利用Matlab的绘图功能，可以直观地展示波场在介质中的传播和反射情况。
4. 我对波动方程PML在Matlab中的理解
在实际工程和科研中，波动方程PML在Matlab中的应用具有重要意义。通过对复杂介质中波的传播行为进行数值模拟，可以更深入地理解波与介质之间的相互作用，为声波成像、地震勘探和无损检测等领域的研究提供有力支持。我也认为在学习和研究过程中，深入理解和掌握PML方法的原理和在Matlab中的实现技巧，能够为我们提升数值模拟的能力和水平，从而更好地应用于实际问题的求解和分析中。
地下空洞地震波场的旋转交错网格有限差分数值模拟-精选文档
地下空洞地震波场的旋转交错网格有限差分数值模拟
1 引言
天然形成或人工挖凿形成的地下空洞，对城市基础建设具有潜在巨大隐患。地震勘探在探测地下空洞方面取得了显著的进展，为了更好的理解地震波在地下存在空洞时传播的特征，本文采用旋转交错网格进行了高阶旋转交错网格空洞波场数值模拟，并分析了震源子波频率、地下空洞埋深及大小和形状对地震记录的影响。在旋转交错网格一阶速度-应力波动方程波场数值模拟中，相同物理量的不同分量都被定义在单元网格中的一个相同位置，应力（质点速度）定义在单元网格的顶点，质点速度（应力）在相应的对角节点上。质点密度和弹性常数定义在与应力相同的节点上，当模型存在物性分界面时，只需要对密度进行平均即可，因此在模拟裂隙介质和地形空洞等物性差异较大的非均匀介质时更具优势。另外，为了减少人工边界产生的反射波的影响，在模型的人工边界外侧采用了完全匹配层。经过分析不同模型参数下的地震记录及波场快照，加深了对地下空洞存在情况下地震波传播特征的理解。
2 旋转交错网格差分格式
3 一阶应力速度弹性波方程
二维各向同性介质一阶应力―速度弹性波方程可以表示为： （5）
其中， 为质点速度分量， 为密度的倒数， 为应力分量， 为拉梅系数。
4 模型算例
本文对埋深为2m和5m，边长为2m和5m的方形空洞模型分别在30Hz和80Hz雷克子波震源下进行了数值模拟。空洞周围介质的纵横波速度及密度分别为866m/s、500m/s、2000Kg/m3，模型的左边界距震源40m，网格间距为0.1m，时间步长为0.00005s。
为了进一步研究震源频率、空洞大小及埋深对地震记录的影响，本文计算了含空洞和不含空洞地震记录之差，如图5所示。模型的左边界距震源40m，最小偏移距为5m，即空洞的左边界位于检波器排列的第35道处。
5 结论
通过将CZ真空法和高阶旋转交错网格相结合，较稳定而精确的模拟了地下存在空洞时的地震波传播。通过模拟不同震源子波主频、不同深度和不同大小空洞情况下的地震记录及波场快照，分析了震源子波主频、深度及大小对地震波传播的影响。这有助于探测地下空洞及管道时更合理地进行观测系统及采集参数优化设计，即可根据空洞的大小及埋深适当选择震源频率，以达到分辨率及左右边界同相轴可分辨的效果，同时也为识别空腔地震波提供了一种模拟工具。
地震勘探中波动方程的推导过程
地震勘探中波动方程的推导过程
地震勘探是通过记录和分析地震波在地下传播的情况，来获取地下结构和物性信息的一种方法。地震波动方程是描述地震波在地下传播过程中的数学模型，它是地震勘探研究中的重要理论基础。本文将通过推导地震波动方程的过程，介绍地震波在地下传播的基本原理，让读者对地震勘探有一个更深入的理解。
1.地震波动方程的基本概念
地震波动方程是描述地震波在地下传播的数学模型，它是通过物理规律和方程推导出来的，用来描述地震波在地下传播的基本规律。地震波动方程通常是一个偏微分方程，描述了地震波在地下传播的速度、能量耗散和波阻尼等物理过程。通过地震波动方程，可以推导出地震波在地下的传播速度、能量耗散和反射、折射等现象，从而获取地下结构和物性信息。
2.地震波动方程的推导过程
地震波动方程的推导过程是通过物理规律和方程推导出来的，主要涉及弹性力学、波动方程和偏微分方程等知识。地震波动方程的基本形式是弹性波动方程，描述了地震波在地下传播的速度和能量耗散等物理过程。下面我们将通过推导地震波动方程的过程，介绍地震波在地下传播的基本原理。
首先，我们需要了解弹性力学的基本概念。在地震波动方程的推导过程中，我们需要用到弹性力学的基本原理和方程。弹性力学是研究物体在外力作用下的形变和应力的力学学科，它描述了物体受力后的形变和应力分布。在地震波动方程的推导过程中，我们将利用弹性力学的基本原理和方程，推导出地震波在地下传播的数学模型。
其次，我们需要了解波动方程的基本概念。波动方程是描述波动过程的数学模型，它描述了波动在空间和时间上的传播规律。在地震波动方程的推导过程中，我们将利用波动方程的基本原理和方程，推导出地震波在地下的传播速度和能量耗散等物理过程。
最后，我们需要了解偏微分方程的基本概念。偏微分方程是描述多维空间中变量的变化规律的数学模型，它描述了变量在空间和时间上的变化规律。在地震波动方程的推导过程中，我们将利用偏微分方程的基本原理和方程，推导出地震波在地下传播的速度和能量耗散等物理过程。 通过以上的基本概念的介绍，我们可以开始推导地震波动方程了。首先，我们需要假设地震波是一种弹性波，即在地下传播时，地震波具有弹性形变和能量耗散等物理过程。其次，我们需要利用弹性力学的基本原理和方程，推导出地震波在地下传播的速度和能量耗散等物理过程。最后，我们需要利用波动方程的基本原理和方程，推导出地震波在地下传播的数学模型。
地震波横波分裂现象的数值模拟特征
地震波横波分裂现象的数值模拟特征
全红娟;朱光明;李桂花
【摘 要】开展了横波分裂的机理、横波分裂的影响因素及横波分裂的可识别条件等3个方面的研究.采用交错网格高阶有限差分法,对垂直裂缝介质中地震波的传播特征进行了数值模拟.研究结果表明:由于入射横波在裂缝介质中存在两个正交极化分量导致横波分裂现象的发生.横波分裂的因素与震源、传播介质、分界面、观测系统等有关.横波分裂可识别的条件是:采用剪切力源或集中力源及VSP或井间地震观测系统,更容易识别横波分裂现象;对于非均匀各向异性介质,如果观测方位合适,会观测到横波分裂和再分裂现象,该项工作对实际中裂隙检测和地震波多波多分量勘探有一定的理论指导意义.
【期刊名称】《西北大学学报（自然科学版）》
【年(卷),期】2014(044)005
【总页数】6页(P745-750)
【关键词】裂缝介质;横波分裂;井间地震
【作 者】全红娟;朱光明;李桂花
【作者单位】长安大学地质工程与测绘学院,陕西西安710054;长安大学理学院,陕西西安710064;长安大学地质工程与测绘学院,陕西西安710054;山东科技大学地质科学与工程学院,山东青岛266590
【正文语种】中 文
【中图分类】P315.3
地球各向异性普遍存在，横波较之纵波对各向异性介质更为敏感，尤其对于裂缝性储集层，故近些年横波勘探方法被认为是行之有效的油气藏勘探方法。横波分裂现象可以携带很多裂缝信息，所以人们对裂缝介质中的横波分裂做了很多物理模拟和数值模拟工作。董敏煜等［1］用有限差分法模拟了弹性波在EDA介质中的传播，获得了反射波、透射波及绕射波等波型的VSP剖面，揭示了横波双折射现象;甘文权等［2］用高阶有限差分法对含裂隙介质中横波分裂现象进行了数值模拟，得出横波分裂现象与裂缝走向方位角有直接关系，快横波比慢横波衰减慢的结论;吴松翰等［3］通过EDA物理模型获得了快慢横波记录，计算出快慢横波的走时时差;裴正林［4］用交错网格高阶有限差分法对层状各向异性介质中横波分裂和再分裂进行了数值模拟，揭示了快横波的多值性和波面三分叉现象;郭桂红等［5］用伪谱法模拟了EDA介质中横波分裂的地面三分量记录，研究了由于裂缝方位和密度引起的横波分裂等。对于实际勘探工作者来讲，最想了解的是横波分裂产生的机理是什么;究竟在何种条件下可以观测到横波分裂现象;横波分裂究竟与哪些因素有关系等;这些问题，首先应考虑横波分裂与震源激励、传播介质、分界面和观测系统等之间的关系，本文针对此，用交错网格高阶有限差分法对垂直裂缝介质中横波分裂现象进行了数值模拟，并对模拟结果展开了详细的分析和讨论。
波动方程_精品文档
博学笃行 自强不息
1
波动方程
波动方程是描述波动现象的数学模型。它是最基本的物理方程之一，广泛应用于各个领域，包括物理学、工程学、地球科学等。波动方程描述了波动传播的机制和特性，是许多领域中研究和分析波动现象的重要工具。
波动方程的一般形式可以表示为：
∇²u = (1/c²) * ∂²u/∂t²
其中，u是波动的物理量，∇²代表拉普拉斯算子，c是波速，∂²u/∂t²是波动量的二阶时间导数。
波动方程的解决了初值问题：给定初始条件下，求解在给定时间和空间范围内波动的传播和变化情况。对于简单的一维情况，波动方程可以简化为：
∂²u/∂x² = (1/c²) * ∂²u/∂t²
这是常用的一维波动方程，描述了波沿着x轴的传播行为。根据边界条件和初值条件，可以求解出特定系统下的波动解。 博学笃行 自强不息
2
波动方程描述了各种类型的波动现象，包括机械波、电磁波、声波等。在物理学中，波动方程常被用于研究弹性体的传播行为，如声波在空气中的传播、地震波在地壳中的传播等。在工程学中，波动方程可以用于分析结构中的振动问题，如桥梁、建筑物等的振动特性。在地球科学中，波动方程被广泛应用于地震勘探和地震波传播等研究。
波动方程的研究可以帮助我们理解和预测波动现象的行为。通过求解波动方程，我们可以得到波的传播速度、波的形状、波的幅度等信息。这些信息对于研究和应用波动现象都非常重要。
除了一维波动方程外，波动方程还可以推广到二维和三维情况。在二维情况下，波动方程可以表示为：
∇²u = (1/c²) * ∂²u/∂t²
这是二维波动方程，描述了波沿着平面的传播行为。在三维情况下，波动方程可以表示为：
∇²u = (1/c²) * ∂²u/∂t²
博学笃行 自强不息
3
这是三维波动方程，描述了波沿着空间的传播行为。对于二维和三维情况，波动方程的求解相对复杂，但同样具有重要的应用价值。
波动方程在现代科学中扮演着重要的角色。它不仅有助于理解波的传播机制和特性，还能用于模拟和预测各种波动现象。通过数值方法和计算机模拟，可以求解波动方程并得到更准确的结果。这为各个领域的研究和应用提供了有力的支持。
地震波场数值模拟中的信号恢复问题
第４７卷第４期 ２００８年７月 石油物探 ＧＥ０ＰＨＹＳＩＣＡＬ ＰＲＯＳＰＥＣＴＩＮＧ ＦＯＲ ＰＥＴＲＯＩ ＥＵＭ Ｖｏ１．４７，Ｎｏ．４ Ｊｕ１．，２００８ 文章编号：１０００—１４４１（２００８）０４—０３５７—０４ 地震波场数值模拟中的信号恢复问题 何宝庆，王秋成，蔡锡伟 （中国石油集团东方地球物理公司采集技术支持部，河北涿州０７２７５１） 摘要：基于波动方程的数值模型正演模拟方法已经被广泛应用于地震数据采集、处理和解释的各个方面。为了 选择最佳的网格尺寸以及以最小的计算量完整地恢复出信号，针对地震波场正演模拟中经常使用的Ｒｉｃｋｅｒ子 波，根据空间采样定理和地震波在不同速度介质中的传播规律，分析了有限差分算法中的假频和信号恢复问题， 推导出模拟网格尺寸、模型地层的最小速度和能够恢复的子波主频之间的近似定量关系，并用数值算例验证了 这一关系。 关键词：正演模拟；Ｒｉｃｋｅｒ子波；有限差分；信号恢复 中图分类号：Ｐ６３１．４ 文献标识码：Ａ 随着油气勘探难度的增加和地球物理勘探水 平的提高，基于射线的地震勘探理论越来越难以解 决比较复杂的地质构造和各向异性问题，地震波正 演模拟技术在地震资料采集、处理和解释过程中正 发挥越来越重要的作用。但在正演模拟过程中，子 波主频与模拟网格大小和模型速度之间存在着微 妙的关系，如果不考虑这种关系，模拟结果就可能 出现假频，从而不能完整地恢复信号。为了进一步 指导波动方程正演模拟技术在实际生产过程中的 应用，本文对信号恢复问题进行了讨论。 １ Ｒｉｃｋｅｒ子波概述 Ｒｉｃｋｅｒ子波是地震勘探理论中常用的子波， 现在已经广泛被人们所接受。俞寿朋提出了宽带 Ｒｉｃｋｅｒ子波 ，其后云美厚 在时间、频率域内对 Ｒｉｃｋｅｒ子波作出了对比和总结，提出了独到的见 解。在地震波场正演模拟中，Ｒｉｃｋｅｒ子波往往作 为外力项（也就是我们通常所说的震源子波）被引 入，所以有必要对Ｒｉｃｋｅｒ子波的特性进行分析。 Ｒｉｃｋｅｒ子波在时间域的表达式为 。 ｒ（￡）一（１—２７ｃ ￡ ）ｅ ｚ （１） 而它在频率域内的表达式为 ９以 ￡ Ｒ（厂）一｛ ｅ （２） 式中：ｔ为时间，厂为频率， 为峰值频率。 图１和图２分别为Ｒｉｃｋｅｒ子波的时间域和频 率域图像。时间域图像从左到右的主频依次是 ７０，５０和３０ Ｈｚ；频率域图像从左到右的主频依次 是３０，５０和７０ Ｈｚ。从时间域和频率域图像可以 看出，随着主频的增加，时间域的子波延迟时间减 小；频率域的频带宽度增加，振幅的最大值减小。