各向异性介质弹性波传播的三维不规则网格有限差分方法

合集下载

基于交错网格有限差分弹性波正演模拟及波场特征分析

基于交错网格有限差分弹性波正演模拟及波场特征分析

基于交错网格有限差分弹性波正演模拟及波场特征分析【摘要】为研究和认识多种储层中弹性波的波场特征,以利于多波地震资料解释,高精度数值模拟是有效的方法之一。

本文在弹性波方程基础上,采用高阶交错网格有限差分技术模拟地震波在各向同性介质和各向异性介质中的传播,可得到不同类型介质的弹性波场。

同时,文中也分析了各向异性系数对多波波场特征的影响。

通过对高精度数值模拟得到的波场快照对比研究表明,该方法可有效获得高精度弹性波正演结果,为研究各种复杂介质中弹性波的波场特征和传播规律奠定了基础。

【关键词】多波多分量波场特征各向异性弹性波正演1 引言随着油气田勘探技术的不断发展[1][2],人们对地震资料的认识也不断加深,纵波地震资料在含油气的显示上存在一些不确定性,单一纵波资料解释的多解性问题尤为突出。

在地震勘探领域中,过去一直把各向同性弹性体理论作为研究地下介质的前提,但是在实际地层中普遍存在各向异性,地下介质的各向异性(如周期薄互层引起的各向异性、以及裂隙引起的各向异性)产生的弹性波场与各向同性介质产生的弹性波场存在着不可忽略的差异。

由此,多波地震勘探作为油储地球物理的主要方法之一应运而生。

在多波资料解释过程中,要求搞清楚储层的岩性与多波的波场特征之间的关系,因此,多波波场数值模拟技术显得非常重要。

高精度数值模拟技术是联系地震、地质、测井以及油藏工程的纽带,其作用主要体现在提高人们对各种复杂介质中地震波传播规律的认知,并可为新技术、新方法提供试验数据,以满足方法技术研究的需要,同时也可以检验解释结果的正确性。

弹性波波动方程高精度数值模拟可以得到全波场信息,包含了地震波的动力学和运动学特点,为准确描述地震波场特征和波的传播规律奠定基础,本文在弹性波方程基础上,采用高阶交错网格有限差分技术模拟地震波在各向同性介质和各向异性介质中的传播,比较地震波在各向同性介质和各向异性介质中的波场响应异同,并分析了各向异性系数对多波波场特征的影响,这对研究各种复杂介质中弹性波的波场特征和传播规律有着重要的意义。

基于旋转交错网格的双相各向异性介质二维三分量波场模拟

基于旋转交错网格的双相各向异性介质二维三分量波场模拟

基于旋转交错网格的双相各向异性介质二维三分量波场模拟林朋;彭苏萍;卢勇旭;王攀【摘要】基于Biot双相介质模型,推导了双相TTI介质二维三分量一阶速度-应力弹性波方程,采用旋转交错网格(RSG)技术建立了各向异性孔隙介质波动方程的二维三分量高精度有限差分格式(FDTD),并引入不分裂卷积完全匹配层(CPML)作为吸收边界条件.为了验证算法可行性,对均匀双相TTI介质中的弹性波场进行了模拟.结果表明:使用旋转交错网格有限差分技术能够模拟出双相TTI介质中存在的快横波、慢横波、快纵波和慢纵波;双相各向异性介质中存在明显的横波分裂、波前面尖角和三分叉现象;不分裂卷积完全匹配层对边界反射的吸收效果较好.通过对比传统交错网格(SSG)和旋转交错网格有限差分技术,证明了旋转交错网格有限差分算法稳定性较强,精度较高,是一种实用的地震波场数值模拟方法.【期刊名称】《煤炭学报》【年(卷),期】2016(041)005【总页数】9页(P1203-1211)【关键词】双相各向异性介质;旋转交错网格;传统交错网格;不分裂卷积完全匹配层【作者】林朋;彭苏萍;卢勇旭;王攀【作者单位】中国矿业大学(北京)煤炭资源与安全开采国家重点实验室,北京100083;中国矿业大学(北京)煤炭资源与安全开采国家重点实验室,北京100083;中国矿业大学(北京)煤炭资源与安全开采国家重点实验室,北京100083;中国矿业大学(北京)煤炭资源与安全开采国家重点实验室,北京100083【正文语种】中文【中图分类】P631林朋,彭苏萍,卢勇旭,等.基于旋转交错网格的双相各向异性介质二维三分量波场模拟[J].煤炭学报,2016,41(5):1203-1211. doi:10. 13225/ j. cnki. jccs. 2015. 1561Lin Peng,Peng Suping,Lu Yongxu,et al. Study on 2D/3C wave propagation in two-phase anisotropic media using the rotated staggered-grid method[J]. Journal of China Coal Society,2016,41(5):1203-1211. doi:10. 13225/ j. cnki. jccs. 2015. 1561地球内部分布不均的岩石裂缝和孔隙在地震尺度上往往表现为各向异性性质,而传统的地震波传播理论仅仅适用于单相介质。

各向异性介质弹性波传播的三维不规则网格有限差分方法

各向异性介质弹性波传播的三维不规则网格有限差分方法

r1 2
+
r2
3
+
54 < 5 x4
1 4!
r1 2
+
r2
4
+ …,
(11)
由上述方程组
,55
<表示为
x
<
的线性组合
5 5
<
x
= η1 <m +1
-
η2 <m
+ η3 <m+2
-
η4 <m- 1
,
(12)
η1
=
1 4
r21
+ 4 r1 r2 + 4 s1 r2 r1 r2 ( r1 + s1 )
,
(13)
1 2!
r1 2
+ s1
2
-
53 < 5 x3
1 3!
r1 2
+
s1
3
+
54 < 5 x4
1 4!
r1 2
+
s1
4
+ …,
(8)
<m
= <i
-
5 < r1 5x 2
+
52 < 5 x2
1 2!
r1 2
2
-
53 < 1 5 x3 3 !
r1 2
3
+
54 < 5 x4
1 4!
r1 2
4
+ …,
2 理论公式
Hale Waihona Puke 方晶系介质弹性系数矩阵 C (9 个独立系数) 定义如

弹性波交错网格高阶有限差分法波场分离数值模拟

弹性波交错网格高阶有限差分法波场分离数值模拟

© 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved.
本页已使用福昕阅读器进行编辑。
福昕软件(C)2005-2007,版权所有,
5 12
石油地球物理勘探
Pi-
1 2
,j
k
Pi+
1 2
, j-
1
k
Pi-
1 2
, j- 1
+ [ T - T - 2 T + k
Si+
1 2
,
j +1
k
Si-
1 2
, j +1
k
Si+
1 2
,j
+ 2 T + T - T ]} k
Si-
1 2
,j
k
Si+
1 2
,
j-
1
k
Si-
1 2
, j- 1
(5)
其他方程的精度为 O (Δt4 +Δx2 N ) 的差分方程同理 可得 ,这里从略 。
3 山东省东营市中国石油大学 (华东) 地球资源与信息学院 ,257061 本文于 2007 年 1 月 30 日收到 ,修改稿于同年 5 月 9 日收到 。 本研究得到国家自然科学基金 (40474041) 、国家 863 专题 (2006AA06Z206) 、CN PC 中青年创新基金 (04 E7040) 、中原油田博士后科研工作站 和 CN PC 物探重点实验室中国石油大学 (华东) 研究室资助 。
+
5vz 5z
5σPzz 5t

三维方位各向异性介质数值模拟及波场分析

三维方位各向异性介质数值模拟及波场分析

三维方位各向异性介质数值模拟及波场分析王宁;吕希华;闫磊;苗长盛【摘要】采用交错网格技术将速度-应力方程中的速度对时间的导数转化为应力对空间的导数,将弹性波动方程表示为与二阶双曲方程等价的一阶应力-速度公式,以实现三维三分量地震波场模拟.对获得的波场快照、VSP记录图像和地面记录图像进行波场分析,发现在波场图像中存在明显的拟P波、拟快横波和慢横波,还出现了横波分裂、横波分裂盲区及波面三分叉等地震波在方位各向异性介质传播时产生的特殊现象.研究结果表明,用交错网格三维模拟方法研究方位各向异性介质对方位各向异性介质中地震波传播反演能起到很强的辅助作用.【期刊名称】《世界地质》【年(卷),期】2010(029)001【总页数】9页(P130-138)【关键词】三维各向异性介质;交错网格;数值模拟;TIH波场分析【作者】王宁;吕希华;闫磊;苗长盛【作者单位】吉林大学,地球探测与信息技术学院,长春,130061;内蒙古地勘局,呼和浩特,010020;吉林大学,地球科学学院,长春,130061;吉林大学,地球科学学院,长春,130061【正文语种】中文【中图分类】P631.440 引言各向异性介质中地震波的传播是勘探地震学中的前沿课题之一,研究地震波在各向异性介质中的传播规律,能为矿产资源勘探和工程勘察等提供方法指导和技术支撑。

19世纪 20年代,Navier,Cauchy和 Poisson等人建立了弹性力学的第一个通用方程组[1-4];60年代,Kraut研究了应用于各向异性介质中的格林函数并得出了定性的波场解[5]; 20世纪,各向异性的研究取得了巨大的进展,1935年,Bruggeman第一次指出层状固体均具有轴对称性,确定此类固体的性质需要 5个弹性参数[6]; Hess对在太平洋上获得的几条折射剖面的 Pn速度进行的分析研究,第一次证实了在整个上地幔中都存在各向异性[7];80年代后,各向异性介质中正演模拟方法的研究取得了很大进展,Crampin et al.对VTI介质用差分方法进行模拟,提出横向各向同性介质中 SH波运动与 P波和 SV波运动解耦合的认识[8];Dablain和Virieux最早将高阶交错网格有限差分法应用于各向同性介质[9];Crampin et al.采用谱法模拟了三维各向异性介质中的弹性波传播[10];Igel et al.使用交错网格有限差分法对各向异性介质进行了波场模拟[11];Ramos和 Ortega使用交错网格有限差分及一阶速度--应力方程模拟了两组斜交垂直裂隙介质中的横波分裂现象[12];傅旦丹和何樵登进行了方位各向异性介质地震弹性波场的伪谱法正演模拟[13];徐天吉等运用Hartley变换模拟了地震波在方位各向异性介质中的传播[14]。

三维非均匀介质波动方程有限差分python开源代码

三维非均匀介质波动方程有限差分python开源代码

文章标题:探索三维非均匀介质波动方程有限差分python开源代码1. 简介在地质勘探、医学成像和地震监测等领域,对三维非均匀介质波动方程的研究与应用日益重要。

而有限差分方法在数值求解波动方程中具有广泛的应用。

在本文中,我们将探讨如何利用Python编程语言实现三维非均匀介质波动方程的有限差分方法,并开源共享相应的代码,以便更多人能够深入理解和应用这一重要领域。

2. 三维非均匀介质波动方程简介三维非均匀介质波动方程描述了波在非均匀介质中的传播规律,是地震勘探、医学成像等领域中常见的数学模型之一。

该方程的数值求解通常采用有限差分方法,通过离散网格化空间和时间来逼近连续的微分方程,从而得到数值解。

3. 有限差分方法有限差分方法是数值求解微分方程的一种常见方法,其基本思想是将微分方程中的导数用差分近似代替,从而将连续的问题转化为离散的问题。

在三维非均匀介质波动方程中,有限差分方法可以有效地模拟波的传播过程,并得到波场的数值解。

4. Python编程实现利用Python编程语言实现三维非均匀介质波动方程的有限差分方法具有许多优势,如简洁易读的代码、丰富的科学计算库等。

在实现过程中,我们可以利用NumPy库进行数组操作,使用Matplotlib库进行波场可视化,并通过SciPy库进行数值求解等。

5. 开源代码共享在本文中,我们将共享我们编写的三维非均匀介质波动方程有限差分Python开源代码,包括空间离散化、时间离散化、边界条件处理、波场更新等关键部分。

我们也会附上详细的注释和使用说明,以便感兴趣的读者能够下载并运行我们的代码,深入理解和学习有限差分方法在波动方程中的应用。

6. 个人观点和理解通过编写三维非均匀介质波动方程的有限差分Python开源代码,我深刻体会到数值模拟在地质勘探、医学成像等领域中的重要作用。

Python作为一种强大的科学计算语言,为我们提供了丰富的工具和库,使得数值模拟变得更加高效和灵活。

三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟

三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟

三维双相各向异性介质弹性波方程交错网格高阶有限差分法模

裴正林
【期刊名称】《中国石油大学学报(自然科学版)》
【年(卷),期】2006(030)002
【摘要】基于Biot理论,给出了三维双相各向异性介质应力-速度弹性波方程交错网格任意偶阶精度有限差分解法,并对三维双相各向异性(TIH)介质中弹性波场进行了模拟.结果表明,弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横波qS2和慢纵波qP2,并清楚地观测到了横波分裂、横波分裂盲点、波面三分叉等特殊现象.快纵波在固相和流相中的相位相同;而慢纵波在固相和流相中的相位相反,即慢纵波在流相中振幅大,而在固相中的振幅较小.另外,三维双相各向异性介质中弹性波场的快纵波和快横波的耦合关系、波的类型、能量分布和相位等都是在三维空间中变化的.
【总页数】5页(P16-20)
【作者】裴正林
【作者单位】中国石油大学,中国石油天然气集团公司物探重点实验室,北
京,100083
【正文语种】中文
【中图分类】P631.443
【相关文献】
1.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟 [J], 裴正林
2.双相各向异性介质弹性波传播交错网格高阶有限差分法模拟 [J], 裴正林
3.任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟 [J], 裴正林
4.三维各向同性介质弹性波方程交错网格高阶有限差分法模拟 [J], 裴正林
5.任意倾斜各向异性介质中弹性波波场交错网格高阶有限差分法模拟 [J], 裴正林;王尚旭
因版权原因,仅展示原文概要,查看原文内容请购买。

三维波动方程有限差分正演方法

三维波动方程有限差分正演方法

三维波动方程有限差分正演方法三维波动方程有限差分正演方法是地球物理勘探领域中常用的数值计算方法之一,其主要应用于地震波传播与反演等领域。

一、三维波动方程有限差分正演方法原理三维波动方程的一般形式可以表示为:\[ \frac{\partial^2 p}{\partial t^2} = \nabla^2 p +f(x,y,z,t) \]其中$p$表示波场,$f(x,y,z,t)$表示源项函数,$\nabla^2$表示拉普拉斯算子。

对于三维波动方程的有限差分正演方法,其基本的数值离散形式如下:\[ \frac{p_{i,j,k}^{n+1} - 2p_{i,j,k}^n + p_{i,j,k}^{n-1}}{\Delta t^2} = c_x^2 \frac{p_{i+1,j,k}^n - 2p_{i,j,k}^n + p_{i-1,j,k}^n}{\Delta x^2} + c_y^2 \frac{p_{i,j+1,k}^n -2p_{i,j,k}^n + p_{i,j-1,k}^n}{\Delta y^2} + c_z^2\frac{p_{i,j,k+1}^n - 2p_{i,j,k}^n + p_{i,j,k-1}^n}{\Delta z^2} + f_{i,j,k}^n \]其中$p_{i,j,k}^n$表示波场在离散网格点$(i,j,k)$处的值,$\Delta t,\Delta x,\Delta y,\Delta z$分别表示时间和空间的离散步长,$c_x,c_y,c_z$分别表示波速在$x,y,z$方向上的离散形式,$f_{i,j,k}^n$表示源项在离散网格点$(i,j,k)$处的值。

该有限差分正演方法可以通过迭代求解,即根据当前时刻$t^n$的波场值$p_{i,j,k}^n$,计算当前时刻$t^{n+1}$的波场值$p_{i,j,k}^{n+1}$。

在迭代过程中,需要进行边界条件处理和源项的更新等操作,以确保该方法的数值计算精度和稳定性。

起伏地表弹性波传播有限差分法数值模拟

起伏地表弹性波传播有限差分法数值模拟

起伏地表弹性波传播有限差分法数值模拟董良国;郭晓玲;吴晓丰;马在田【期刊名称】《天然气工业》【年(卷),期】2007(027)010【摘要】有限差分是进行地震波传播数值模拟的最常用方法,但该方法处理起伏的自由边界比较困难.为此,通过对不同地形起伏情况下自由边界的具体分析,将整个二维空间离散点划分为24类,对每一类自由边界处的网格点选择了合理的表现方式,实现了起伏地表自由边界条件的数值化.该方法可以模拟出地表起伏情况下弹性波复杂的传播现象,为进行起伏地表地震波传播规律研究、山地地震勘探野外观测系统设计、山地地震勘探干扰波分析和识别、以及静校正研究提供了正演模拟工具.模拟实例表明,地形起伏引起面波、体波等地震波型之间的相互转化,产生了大量的散射P波、散射S波和散射面波,尤其是由沿地表传播的强能量面波,在地表起伏及近地表物性突变处产生了大量的强能量散射面波,同时也产生了相对较弱的散射P 波和散射S波,这是造成山地地震资料信噪比低的主要原因.【总页数】4页(P38-41)【作者】董良国;郭晓玲;吴晓丰;马在田【作者单位】同济大学海洋地质国家重点实验室;同济大学海洋地质国家重点实验室;同济大学海洋地质国家重点实验室;同济大学海洋地质国家重点实验室【正文语种】中文【中图分类】P61【相关文献】1.双相各向异性介质弹性波传播交错网格高阶有限差分法模拟 [J], 裴正林2.任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟 [J], 裴正林3.起伏地表条件下2.5维声波方程有限差分法数值模拟 [J], 齐鹏;孙建国4.起伏地表弹性波传播的间断Galerkin有限元数值模拟方法 [J], 薛昭;董良国;李晓波;刘玉柱5.有限积分法与有限差分法在弹性波数值模拟中的对比分析 [J], 李明智;熊章强;张大洲因版权原因,仅展示原文概要,查看原文内容请购买。

三维非均匀介质波动方程有限差分python开源代码

三维非均匀介质波动方程有限差分python开源代码

三维非均匀介质波动方程有限差分python开源代码(实用版)目录一、引言1.1 背景介绍1.2 三维非均匀介质波动方程有限差分算法的意义二、相关理论知识2.1 有限差分算法原理2.2 三维非均匀介质中波动方程的求解三、Python 开源代码实现3.1 代码结构概述3.2 关键代码模块解析四、结论与展望4.1 对比现有方法的优势4.2 潜在改进方向正文一、引言1.1 背景介绍在地球物理勘探、海洋工程等领域,三维非均匀介质波动方程的研究具有重要意义。

为了更好地理解和模拟地壳构造、地震波传播等现象,学者们致力于寻求高效的数值计算方法来求解这类偏微分方程。

有限差分方法是一种广泛应用于求解偏微分方程的数值方法,它将连续的空间和时间离散化,从而降低问题的复杂度。

1.2 三维非均匀介质波动方程有限差分算法的意义对于三维非均匀介质中的波动方程,有限差分算法可以提供一种有效的求解手段。

通过将连续的空间和时间离散化,可以降低问题的复杂度,使得求解过程更加高效。

此外,有限差分算法具有较好的稳定性和收敛性,可以较准确地模拟波动现象。

因此,研究三维非均匀介质波动方程有限差分算法具有重要的理论和实际意义。

二、相关理论知识2.1 有限差分算法原理有限差分算法是一种常用的数值计算方法,它通过将连续的空间和时间离散化,将偏微分方程转化为离散形式的代数方程组。

这种离散化的过程需要满足稳定性条件,以保证数值解的准确性和稳定性。

在有限差分算法中,通常采用中心差分格式来对偏微分方程进行离散化。

2.2 三维非均匀介质中波动方程的求解在三维非均匀介质中,波动方程的求解可以采用有限差分算法。

通过将连续的空间和时间离散化,可以将波动方程转化为离散形式的代数方程组。

然后,可以通过迭代法等数值方法求解这个代数方程组,得到数值解。

这种方法可以有效地模拟三维非均匀介质中的波动现象。

三、Python 开源代码实现3.1 代码结构概述本文提供的 Python 开源代码实现了三维非均匀介质波动方程有限差分算法。

文献检索报告(石油大学)

文献检索报告(石油大学)

2014—2015学年第1学期《计算机信息检索》试卷专业班级勘查11级(卓越班)姓名王嘉骏学号 ********开课系室地球物理系考试日期 2014年11月20日题号总分得分阅卷人《计算机信息检索》报告要求一、报告内容选择石油地球物理勘探中重要的方法、技术为主题开展计算机信息检索,并针对检索结果开展全面深入的分析,同时在文献分析的基础上基于文献管理专业软件撰写出一篇格式规范的科技论文。

要求采用《地球物理学报》的最新投稿格式撰写论文,要求正文内容不少于10页、能正确的参照指定格式完成论文撰写、参考文献格式正确、文献种类多样化,参考文献著录规则参照《文后参考文献著录规则》(GB-T_7714-2005)。

二、检索主题建议选择选择如下主题开展计算机信息检索与分析,也可选择与石油地球物理专业相关的其它主题。

地震资料采集、地震资料处理、地震资料解释、动校正、静校正、速度分析、组合、叠加、深度偏移、时间偏移、逆时偏移RTM、AVO、弹性阻抗反演、波阻抗反演、全波形反演FWI、地震反演、地震属性、井间地震、VSP、逆VSP、随钻地震、各向异性、粘弹性介质、双相介质、地震子波、时深转换、多波多分量、转换波、多次波、面波、去噪、衰减、正演模拟、广角反射、地震照明、反褶积、观测系统、分辨率、震源、电法、磁法、微地震、微测井、低速带、层析成像三、检索要求1、检索数据库范围:中文CNKI、英文SCI、EI或SEG数据库;2、通过CNKI检索查询后列出该主题在国内的主要研究单位和研究人员分布情况,并用图表形式展示近20年内该领域的文章发表数量,结合图表分析其发展趋势;3、通过SCI、EI或SEG数据库检索查询后列出国际英文期刊上发表的该主题文章的主要研究单位和作者,并用图表形式展示近20年该领域文章发表的数量,并分析其趋势;4、分别阅读中、英文数据库中引用次数最多的两篇文献、最新的两篇文献,并尽量分别阅读一篇该领域内的中文和英文文献综述,根据阅读情况总结该方向的研究现状,对研究现状的分析需要加上相应的参考文献;5、根据计算机信息检索结果深入分析该研究方向今后的发展趋势,并结合自己的兴趣特点来分析如果让你从事相关研究工作的话应该选择什么样的研究课题,根据分析结果提出你本科毕业设计希望从事的具体研究方向,并根据该方向拟定出一个初步的研究方案;6、根据上述任务完成一个报告,报告内容至少包括但不限于报告内容格式中规定的范围(见附页);7、报告格式参照《地球物理学报》论文投稿格式执行,篇幅10页以上。

三维波动方程有限差分正演方法

三维波动方程有限差分正演方法

三维波动方程有限差分正演方法摘要:本文主要介绍了三维波动方程有限差分正演方法的基本原理和实现过程,并对其优缺点进行了分析。

通过数值实验验证了该方法的可行性和准确性,为地震勘探、地下水文学等领域的研究提供了参考。

关键词:三维波动方程、有限差分、正演方法、数值实验一、引言三维波动方程是地震勘探、地下水文学等领域的基础理论,其求解方法对于相关领域的研究具有重要意义。

有限差分正演方法是求解三维波动方程常用的数值方法之一,其具有计算速度快、精度高等优点,在实际应用中得到了广泛的应用。

本文将介绍三维波动方程有限差分正演方法的基本原理和实现过程,并对其优缺点进行分析,同时通过数值实验验证该方法的可行性和准确性。

二、三维波动方程有限差分正演方法的基本原理三维波动方程可以表示为:![image.png](attachment:image.png)其中,u(x,y,z,t)为波动场,c(x,y,z)为介质速度,ρ(x,y,z)为介质密度,t为时间。

有限差分正演方法通过将空间和时间离散化,将三维波动方程转化为差分方程,进而求解波动场在不同时刻的数值解。

具体而言,有限差分正演方法将空间和时间分别离散化,将空间网格点和时间网格点相结合,构成一个四维网格空间,其中每个网格点对应一个波动场的数值解。

有限差分正演方法的求解过程可以分为以下几个步骤:1. 空间离散化对于三维空间,可以将其分为x、y、z三个方向,分别进行离散化。

假设在x方向上,空间步长为Δx,在y方向上,空间步长为Δy,在z方向上,空间步长为Δz。

则可以将空间网格点表示为(xi,yj,zk),其中i=1,2,...,nx,j=1,2,...,ny,k=1,2,...,nz,nx、ny、nz分别为x、y、z方向的网格数。

2. 时间离散化假设时间步长为Δt,则在t时刻波动场的数值解可以表示为ui,j,k^n,其中n表示时间步数,i、j、k表示空间网格点的索引。

3. 有限差分近似将三维波动方程中的导数项用有限差分近似表示,例如:![image-2.png](attachment:image-2.png)其中,ui,j,k^n表示在t时刻,在(xi,yj,zk)处的波动场数值解,c(xi,yj,zk)表示在(xi,yj,zk)处的介质速度,ρ(xi,yj,zk)表示在(xi,yj,zk)处的介质密度,Δx、Δy、Δz、Δt分别为空间和时间步长。

几类各向异性界面问题的有限元-有限差分方法

几类各向异性界面问题的有限元-有限差分方法

几类各向异性界面问题的有限元-有限差分方法几类各向异性界面问题的有限元-有限差分方法引言:在实际的工程及科学计算中,很多问题涉及到界面及其附近区域的模拟与分析。

针对各向异性介质中的界面,选择合适的数值模拟方法是非常重要的。

本文将探讨几类各向异性界面问题,并介绍有限元-有限差分方法在这些问题中的应用。

一、各向异性界面问题的数学建模各向异性材料的界面问题在许多科学领域中都有重要的应用,如材料科学、固体力学、地球物理学等。

各向异性界面问题的建模通常基于弹性力学理论,其中关键的数学方程是拉普拉斯方程和散度方程。

拉普拉斯方程描述了界面附近的势场分布,而散度方程描述了场的散度。

通过适当的边界条件,我们可以对各向异性材料中的界面问题进行建模。

二、有限元-有限差分方法的基本原理有限元-有限差分方法是一种常用的数值计算方法,广泛应用于工程学和科学计算中。

它通过将问题的定义域划分为多个小区域,并在每个小区域内求解偏微分方程,最终得到整个系统的数值解。

该方法结合了有限元方法和有限差分方法的优点,能够有效地处理各向异性问题。

三、几类典型各向异性界面问题及其求解方法1. 各向异性界面的热传导问题在热传导问题中,考虑材料的各向异性性质对界面附近的温度分布和热流动态的影响是非常重要的。

有限元-有限差分方法可以在微观尺度上模拟热传导过程,并考虑界面处的各向异性性质。

通过适当的边界条件和材料参数,可以获得准确的界面温度分布和热流动态。

2. 各向异性界面的电磁场问题在电磁场问题中,各向异性界面通常对电场和磁场的分布产生重要影响。

有限元-有限差分方法可以将问题离散化,并在界面处进行适当的数值处理。

通过求解麦克斯韦方程组,可以得到各向异性材料中界面附近的电场和磁场分布情况。

3. 各向异性界面的流体力学问题在流体力学问题中,各向异性界面可以对流体流动产生重要影响,如湍流的发生和分离现象。

有限元-有限差分方法能够将流体力学方程离散化,并在界面处考虑各向异性的影响。

三维波动方程双变网格有限差分并行模拟方法研究的开题报告

三维波动方程双变网格有限差分并行模拟方法研究的开题报告

三维波动方程双变网格有限差分并行模拟方法研究的开题报告一、研究背景与意义地震勘探是石油勘探中不可或缺的一种手段。

在地震勘探中,波动方程数值模拟是一种快速、经济、精确的方法。

在三维区域中,波动方程数值模拟需要处理大量数据,因此需要高效的并行算法来降低计算时间。

然而,双变网格有限差分方法因其高精度和效率等特点,成为三维波动方程数值模拟中的一种主流方法。

因此,研究双变网格有限差分并行算法具有重要意义。

二、研究内容本文将研究三维波动方程双变网格有限差分方法的并行化实现。

具体来说,我们将研究以下内容:1. 双变网格有限差分方法基础理论的研究和分析,包括双变网格有限差分方法的基本原理、误差分析、收敛性和稳定性分析等。

2. 并行算法的研究,包括并行算法的设计和实现,以及并行算法的效率分析等。

3. 数值实验的设计和实现,验证所提出的算法的正确性和稳定性,并进行性能测试,比较不同算法之间的差异。

三、研究方法1. 理论分析方法:通过理论推导和分析,研究双变网格有限差分方法的基本原理、误差分析、收敛性和稳定性等。

2. 并行算法设计和实现:通过编程实现,设计出适合双变网格有限差分方法的并行算法,并进行大规模数据的性能测试。

3. 数值实验:通过使用已知的三维波动方程数值模拟数据,验证所提出的算法的正确性和稳定性,并进行性能测试,比较不同算法之间的差异。

四、研究现状双变网格有限差分方法是一种高效的数值模拟方法,近年来受到越来越多的关注。

在并行算法的研究方面,已有一些研究成果,如OpenMP 并行算法、MPI并行算法等。

但这些算法存在一些局限性,在处理大规模数据上效果不是很好。

因此,还需要进一步研究和发展更加高效的并行算法。

五、研究计划时间节点和工作安排如下:1. 第1-2周:阅读文献,熟悉已有研究成果和方法,初步确定研究内容和算法设计。

2. 第3-4周:学习双变网格有限差分方法的基础理论和并行算法设计。

3. 第5-6周:完成双变网格有限差分方法的理论推导和分析,并结合实际案例进行验证。

三维波动方程有限差分正演方法

三维波动方程有限差分正演方法

三维波动方程有限差分正演方法
首先,将三维波动方程进行离散化处理,将连续的空间和时间域离散
为网格点的集合。

在空间域中,将地下介质划分为多个均匀网格,每个网
格点上对应一个地震波场的值。

在时间轴上,将时间域离散为多个时间步长,每个时间步长表示一个时刻。

然后,利用有限差分算子将三维波动方程转化为离散形式。

常用的有
限差分格式有常系数差分格式和变系数差分格式。

常系数差分格式适用于
各向同性介质,而变系数差分格式适用于各向异性介质。

在三维波动方程有限差分正演方法中,需要考虑各向异性介质的模型。

各向异性介质与各向同性介质不同,其波速和波阻抗在不同方向上具有不
同的值。

因此,在离散化过程中,需要引入特殊的差分算子来考虑各向异
性的影响。

最后,利用迭代求解的方法,按时间步长依次求解离散化的波动方程。

利用差分算子和初始条件,根据时空的变化逐步更新波场的数值,并计算
出每个网格点的波动值。

通过不断迭代求解,最终可以得到地震波在地下
的传播和反射结果。

三维波动方程有限差分正演方法在地震勘探中具有重要的应用价值。

它可以模拟地震波在地下介质中的传播和反射过程,帮助研究人员了解地
下结构、油气储层等。

同时,它也可以用来研究地震波在各向异性介质中
的传播规律,提高地震勘探的准确性和效率。

总之,三维波动方程有限差分正演方法是一种常用的数值模拟方法,
可以用来模拟地震波在地下介质中的传播和反射。

它具有广泛的应用价值,对地震勘探、岩石物理、地质测井等领域具有重要的意义。

一种三维正交方位各向异性介质岩石物理建模及弹性波正演模拟方法

一种三维正交方位各向异性介质岩石物理建模及弹性波正演模拟方法

一种三维正交方位各向异性介质岩石物理建模及弹性波正演模拟方法李雨生;吴国忱【摘要】In this paper,by rock physical equivalent theory and linear slip theo-ry,two sets of orthorhombic fractures are equivalent to a new three-dimensional azimuthally anisotropic medium of orthorhombic symmetry for rock physical modeling.And then elastic wave equations are solved to simulate the propaga-tion process in this medium by high-order staggered-grid finite-difference method.During the process of modeling,with the variation of physical parame-ters,we analyze CSP (common shot points)sets and the wave field characteris-tics on the conditions of different fracture densities,as well as azimuthal charac-teristics of orthorhombic anisotropy.The research results show that the aniso-tropic strengths enhance with fracture density increasing,which is reflected in the CSP sets and forwarding wavefield.%通过线性滑动理论和岩石物理等效理论,将两组正交直立裂隙介质等效为一种正交方位各向异性介质进行三维岩石物理建模,通过高阶交错网格有限差分求解弹性波动方程模拟地震波在该种介质中的传播过程。

各向异性介质弹性波高阶交错网格有限差分模拟

各向异性介质弹性波高阶交错网格有限差分模拟
An i s o t r o p i c me d i a e l a s t i c wa v e h i g h - o r d e r s t a g g e r e d - g r i d f i n i t e d i fe r e n c e s i mu l a t i o n
Ab s t r a c t : Th e e l a s t i c wa v e e q u a t i o n i s s i mu l a t e d b y t h e h i g h — o r d e r s t a g g e r e d — g r i d i f n i t e d i f f e r e n c e , t h e s t a b i l i t y a n d a s t r i n g e n c y a r e a n a l y z e d , i n wh i c h t h e a b s o r b i n g b o u n d a r y c o n d i t i o n s a n d a t t e n u a t i o n r e g i o n s a r e a d d e d i n . T he r e - s ui t s o f t h e ni a s o t r o p i c a n d i s o t r o p i c me d i a s h o w t h a t t h e h i g h - o r d e r d i f f e r e n c e e l a s t i c wa v e s i mu l a t i o n c a n ma k e a g o o d e f f e c t b e c a u s e o f l e s s g r i d d i s p e r s i o n a n d h i g h e r a c c u r a c y .

三维各向同性介质弹性波方程交错网格高阶有限差分法模拟

三维各向同性介质弹性波方程交错网格高阶有限差分法模拟

三维各向同性介质弹性波方程交错网格高阶有限差分法模拟裴正林
【期刊名称】《石油物探》
【年(卷),期】2005(044)004
【摘要】给出了三维各向同性介质中一阶应力速度弹性波方程交错网格任意偶数阶精度有限差分格式,推导出了三维各向同性弹性介质完全匹配层吸收边界条件公式和相应的交错网格高阶有限差分格式.对三维French模型进行了弹性波模拟,结果表明,该方法模拟精度高,边界吸收效果好.在三维French模型的xoz平面和yoz 平面的弹性波场快照中可以见到断面反射波、侧面反射波和散射波等波场特征.多次侧面反射波和多次散射波说明不同三维构造所形成的波场是相互影响的.
【总页数】8页(P308-315)
【作者】裴正林
【作者单位】石油大学CNPC物探重点实验室,北京,100083
【正文语种】中文
【中图分类】P631.4
【相关文献】
1.三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟 [J], 裴正林
2.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟 [J], 裴正林
3.任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟 [J], 裴正林
4.基于BISQ模型的各向同性孔隙介质弹性波三维交错网格高阶有限差分数值模拟
[J], 李红星;刘财;陶春辉
5.双相TI介质中弹性波交错网格高阶有限差分法数值模拟 [J], 尹学爱;邱光辉因版权原因,仅展示原文概要,查看原文内容请购买。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

η2
=
1 4
r21
+ 4 r1 s1 + r1 s1 ( r1 +
4 r2 s1 r2 )
,
(14)
η3 = -
1 4
r2 ( r1
+
r21 r2 ) ( r1
+
r2
+
s1 )
,
(15)
η4 = -
1 4
s1 ( r1
+
r21 s1 ) ( r1
+
r2
+
s1 )
,
(16)
同理 ,容易得到55 <y和55 <z . 如果是规则网格 ( s1 = r1 = r2 =Δx) ,则
图 1 交错网格有限差分示意图 Fig. 1 Staggered grids finite difference
笛卡儿坐标系斜方晶系各向异性模型弹性波动
方程速度 - 应力公式为[18]
ρv = D ·T ,
(1)
T = C ·DT ·v ,
(2)
其中 ρ( x) 为 介 质 密 度 , 弹 性 介 质 振 动 速 度 矢 量
第 47 卷 第 2 期
地球物理学报
Vol. 47 , No. 2
2004 年 3 月
CHINESE JOURNAL OF GEOPHYSICS
Mar. , 2 0 0 4
Sun W T , Yang H Z. A 32D finite difference method using irregular grids for elastic wave propagation in anisotropic media. Chinese J . Geo2 phys. (in Chinese) , 2004 , 47 (2) :332~337
r1 2
+
r2
3
+
54 < 5 x4
1 4!
r1 2
+
r2
4
+ …,
(11)

由上述方程组
,55
<表示为
x
<
的线性组合
5 5
<
x
= η1 <m +1
-
η2 <m
+ η3 <m+2
-
η4 <m- 1
,
(12)
η1
=
1 4
r21
+ 4 r1 r2 + 4 s1 r2 r1 r2 ( r1 + s1 )
,
(13)
1 2!
r1 2
+ s1
2
-
53 < 5 x3
1 3!
r1 2
+
s1
3
+
54 < 5 x4
1 4!
r1 2
+
s1
4
+ …,
(8)
<m
= <i
-
5 < r1 5x 2
+
52 < 5 x2
1 2!
r1 2
2
-
53 < 1 5 x3 3 !
r1 2
3
+
54 < 5 x4
1 4!
r1 2
4
+ …,
s1 = xm - xm- 1 ,
r1 = xm +1 - xm ,
r2 = xm +2 - xm+1 .
(7)
节点 m 和 m + 1 的中点命名为节点 i . 节点
m - 1 、m 、m + 1 、m + 2 处的波场值 < 写成级数形式
<m - 1 = <i
-
5< 5x
r1 2
+
s1
+
52 < 5 x2
(4)
5 5x
0
0
0
55 5z 5y
D=
0
5 5y
0
5 5z
0
5 5x
,
(5)
0
0
555 5z 5y 5x
0
c11 c12 c13 0 0 0
c12 c22 c23 0 0 0
c13 c23 c33 0 0 0
C=
. (6)
0 0 0 c44 0 0
0 0 0 0 c55 0
0 0 0 0 0 c66 波动方程中包含速度和应力的一阶导数 ,下面 推导具有四阶空间精度和二阶时间精度的不规则网 格差分算子. 在交错网格 (见图 1) 上离散一阶偏微 分算子 ,笛卡儿坐标轴上四节点示意图见图 2 ,它们 之间的间距分别是 s1 、r1 、和 r2 .
基金项目 中国石油天然气集团公司基金资助 (2002CXKF24) . 作者简介 孙卫涛 ,男 ,1975 年生 ,1998 年毕业于大连理工大学工程力学系 ,2003 年在清华大学工程力学系获硕士 、博士学位 ,现在清华大学计
算机系做博士后研究 ,主要从事弹性波动力学 、高性能计算的理论和方法研究. E2mail : sunwt @mail . tsinghua. edu. cn
η1
= η2
=
9 8Δx
,
(17)
η3
= η4
=-
1 24Δx
,
(18)
5 5
<
x
=
9 8
( <m +1 Δx
<m )
-
1 24
( <m +2 - <m- 1 ) Δx
,
(19)
这是 Levander[5] 给出的规则网格四阶差分算子. 本
文提出具有更一般形式的不规则网格公式 ,模拟弹
(9)
<m +1
= <i
+
5 5
<
x
r1 2
+
52 < 5 x2
1 2!
r1 2
2
+
53 < 5 x3
1 3!
r1 2
3
+
54 < 5 x4
1 4!
r1 2
4
+ …, (10)
<m +2
= <i
+
5 5
<
x
r1 2
+
r2
+
52 < 5 x2
1 2!
r1 2
+
r2
2
+
53 < 5 x3
1 3!
2 期 孙卫涛等 :各向异性介质弹性波传播的三维不规则网格有限差分方法
333
1 引 言
20 世纪 70 年代以来 ,在地震波传播和地震地 表运动的数值模拟研究中 ,有限差分方法成为求解 波动方程的最有力工具之一. 地震波正演模拟中有 限差分方法的早期研究见 Alterman[1] , Kelly 等[2] 和 Virieux[3 ,4] 等学者的著作 ,Levander[5] 在 P2SV 波模拟 中引入了四阶空间差分算子 , Graves[6] 给出了等效物 质参数的三维四阶速度 - 应力有限差分方法 , Dablain[7] 提出了高阶差分算子的方法. 这些工作全 部基于笛卡儿坐标系中的规则网格 ,用普通网格模 拟曲线界面时出现“阶梯状”边界 ,在地形构造模型 中产生虚假衍射波. 另外 ,局部物理参数的变化也 会要求加密整个模型网格 ,导致计算量的增加. Shortley 等[8] 首先研究了 Laplace 方程中的不规则网 格有限差分方法 , Jastram 等[9] 和 Falk 等[10] 给出了 交错网格上的可变网格差分方法 , Tessmer 等[11] 和 Hestholm 等[12] 用变形的矩形网格模拟曲线边界 , Ivo Oprsal 等[13] 研究了非均匀介质中的二阶波动方程的 不规则网格差分方法 , Pitarka[14] 提出了各向同性介 质中不规则网格的有限差分方法 ,董良国等[15 ,16] 研 究了交错网格高阶差分方法 ,张剑锋[17] 通过坐标映 射得到非规则网格差分法. 但是 ,以往变形网格有 限差分方法需要波动方程在不同坐标系之间的正反 变换 ,正演模拟的计算量相当大 ;对于具有非平面界 面的三维非均匀各向异性介质模型 ,不规则网格有 限差分方法目前还鲜有相应的研究 ,普通有限差分 方法仍然面临着计算规模和计算精度的限制. 高精 度 、高计算速度正演成为三维复杂介质弹性波传播 问题的研究热点. 本文给出一种具有高阶精度 、非 均匀各向异性波动方程不规则网格有限差分方法. 这种方法无需坐标变换和网格间的插值 ,简单易行 而且占用内存少 、计算量小 ,适合于三维非均匀各向 异性介质模型的弹性波正演问题.
Abstract This paper presents a new 3D finite2difference ( FD) method using spatially irregular grids to sim2 ulate elastic wave propagation in heterogeneous anisotropic media with topographic structures. The method ap2 proximates the first2order elastic wave equations by the finite difference operators on irregular grids with sec2 ond2order time precise and fourth2order spatial precise . Unlike the multi2grid scheme , this method has no in2 terpolation between the fine and coarse grids. All grids are computed at the same spatial iteration. Complicat2 ed geometrical structures like rough submarine interfaces , faults and nonplanar interfaces are treated with fine irregular grids. Theoretical analysis and numerical simulations show that this method saves considerable memo2 ry and computing time , at the same time , has satisfactory stability and accuracy. The proposed scheme is more efficient than conventional methods in simulating seismic wave propagation in complex topographic struc2 tures. Key words Seismic wave , Irregular grid , Finite difference ,Anisotropy media.
相关文档
最新文档