CN109490956B - 一种基于交错网格的声波波动方程正演模拟方法及装置 - Google Patents
一种基于交错网格的声波波动方程正演模拟方法及装置 Download PDFInfo
- Publication number
- CN109490956B CN109490956B CN201811354411.5A CN201811354411A CN109490956B CN 109490956 B CN109490956 B CN 109490956B CN 201811354411 A CN201811354411 A CN 201811354411A CN 109490956 B CN109490956 B CN 109490956B
- Authority
- CN
- China
- Prior art keywords
- wave
- wave equation
- time
- sound wave
- difference
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 117
- 239000006185 dispersion Substances 0.000 claims abstract description 69
- 238000004088 simulation Methods 0.000 claims abstract description 36
- 238000010521 absorption reaction Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 14
- 239000000126 substance Substances 0.000 claims description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000002123 temporal effect Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000013508 migration Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013211 curve analysis Methods 0.000 description 2
- 238000013213 extrapolation Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及地震波场数值模拟技术领域,公开了一种基于交错网格的声波波动方程正演模拟方法及装置,该方法包括:获取地震参数;建立基于交错网格的声波波动方程;采用时空域有限差分法计算所述声波波动方程的频散关系式;根据所述频散关系式获取波场模拟所满足的稳定条件;采用吸收边界条件对声波波动方程进行波场延拓,获取波场及地震记录;通过交错网格建立新的差分结构,求取的差分系数能够使声波在更大的波数范围内压制数值频散,进一步提高了声波波动方程的模拟精度。
Description
技术领域
本发明涉及地震波场数值模拟技术领域,尤其涉及一种基于交错网格的声波波动方程正演模拟方法及装置。
背景技术
地震波成像和反演需要一种高效率和高精度的算法来模拟波的正向和反向传播,例如当前热门的逆时偏移和全波形反演技术。因此研究高精度和高效率的数值模拟算法很有必要。
有限差分法是一种灵活且简便的数值算法,已被广泛应用在波动方程的数值求解之中。对于有限差分法而言,同时压制波动方程的时间和空间频散是最大的挑战之一。1986年,Dablain通过运用LAX-WENDROFF方法来提高声波波动方程的模拟精度,即将波动方程中高阶时间偏导使用空间偏导替代,但这种方法计算量将显著增加;传统的有限差分法的空间差分系数是基于空间域频散关系推导得到的,2007年,Finkelstein等提出在时间空间域确定空间差分系数,该方法使得时空域频散方程在指定的若干个频率上严格成立,从而获得若干个方程,进而求解出空间差分系数,虽然获得的差分系数能够减小频散关系式,但因为不同的频率需要不同的差分系数因此难以用于实际。为了在不显著增加计算机内存的条件下提高模拟精度,Liu等基于时间空间域频散关系式获得了声波波动方程的有限差分系数,相比于基于空间域频散关系获得的差分系数,在相同的离散条件下,一维模拟无条件稳定且时间空间均为阶精度,二维模拟在波场传播的8个方向时间达到阶精度,三维情况48个方向达到阶精度,而在其他的传播方向上时间精度仍为二阶;为了使波传播的各个方向空间和时间达到阶精度,2013年,Liu等发展了一种菱形差分格式,但计算效率大大地降低。为了同时提高模拟精度和效率,2014年,Tan等基于新差分结构发展了时间四阶、空间任意偶数阶精度的交错网格有限差分法,其差分系数是经泰勒级数展开获得的,这种新的差分格式保留了传统交错网格有限差分控制空间频散关系式的优势,同时增强了其压制时间频散关系式的能力;为了在保持时间四阶精度的同时进一步的提高空间精度,同年,Tan等提出采用非线性优化方法来求得优化的差分系数,然而,优化方法需要重复的迭代计算比较耗时,而Chen等则通过在最小二乘意义下使得交错网格有限差分算子与一阶波数-空间算子的波数响应误差最小化来优化差分系数。随后,在2016年,张保庆等发展了交错网格剖分的空间任意阶,时间四阶精度的有限差分法。但模拟精度和稳定性都不够高。
发明内容
本发明的主要目的在于提出一种基于交错网格的声波波动方程正演模拟方法及装置,通过交错网格建立新的差分结构,求取的差分系数能够使声波在更大的波数范围内压制数值频散,进一步提高了声波波动方程的模拟精度。
为实现上述目的,本发明提供的一种基于交错网格的声波波动方程正演模拟方法,包括:
获取地震参数;
建立基于交错网格的声波波动方程;
采用时空域有限差分法计算所述声波波动方程的频散关系式;
根据所述频散关系式获取波场模拟所满足的稳定条件;
采用吸收边界条件对声波波动方程进行波场延拓,获取波场及地震记录。
可选地,所述地震参数包括:正演模拟所需的速度场文件、差分算子阶数、震源函数及其主频、正演所采用的时间与空间步长及地震记录时长,海绵吸收边界条件的参数。
可选地,所述建立基于交错网格的声波波动方程包括:
时间偏导数采用二阶差分离散;其公式为:
其中,ρ代表介质密度,u和w分别代表水平和垂直方向的速度波场,x,z为直角坐标系的两个轴,t为时间,v为地震波的传播速度,p代表声波波场;
空间偏导数的公式为:
可选地,所述采用时空域有限差分法计算所述声波波动方程的频散关系式包括:
计算声波波动方程的离散形式;
获取时空域有限差分法与LAX-WENDROFF方案的关系式;
采用褶积微分算子法求解所述差分系数;
利用平面波理论,得到时空域频散关系式。
可选地,所述声波波动方程的离散形式计算公式为:
作为本发明的另一方面,提供的一种基于交错网格的声波波动方程正演模拟装置,包括:
获取模块,用于获取地震参数;
建模模块,用于建立基于交错网格的声波波动方程;
差分模块,用于采用时空域有限差分法计算所述声波波动方程的频散关系式;
模拟模块,用于根据所述频散关系式获取波场模拟所满足的稳定条件;
延拓模块,用于采用吸收边界条件对声波波动方程进行波场延拓,获取波场及地震记录。
可选地,所述地震参数包括:正演模拟所需的速度场文件、差分算子阶数、震源函数及其主频、正演所采用的时间与空间步长及地震记录时长,海绵吸收边界条件的参数。
可选地,所述建立基于交错网格的声波波动方程包括:
时间偏导数采用二阶差分离散;其公式为:
其中,ρ代表介质密度,u和w分别代表水平和垂直方向的速度波场,x,z为直角坐标系的两个轴,t为时间,v为地震波的传播速度,p代表声波波场;
空间偏导数的公式为:
可选地,所述采用时空域有限差分法计算所述声波波动方程的频散关系式包括:
计算声波波动方程的离散形式;
获取时空域有限差分法与LAX-WENDROFF方案的关系式;
采用褶积微分算子法求解所述差分系数;
利用平面波理论,得到时空域频散关系式。
可选地,所述声波波动方程的离散形式计算公式为:
本发明提出的一种基于交错网格的声波波动方程正演模拟方法及装置,该方法包括:获取地震参数;建立基于交错网格的声波波动方程;采用时空域有限差分法计算所述声波波动方程的频散关系式;根据所述频散关系式获取波场模拟所满足的稳定条件;采用吸收边界条件对声波波动方程进行波场延拓,获取波场及地震记录;通过交错网格建立新的差分结构,求取的差分系数能够使声波在更大的波数范围内压制数值频散,进一步提高了声波波动方程的模拟精度。
附图说明
图1为本发明实施例一提供的一种基于交错网格的声波波动方程正演模拟方法的流程图;
图2为图1中步骤S30的方法流程图;
图3为本发明实施例一提供的频散误差对比示意图;
图4为本发明实施例一提供的稳定性对比图;
图5为本发明实施例一提供的均匀速度模型中0.6s时刻的波场快照图;
图6为本发明实施例一提供的复杂的Marmousi速度模型;
图7为本发明实施例一提供的Marmousi速度模型中的4.0s时刻的波场快照图;
图8为本发明实施例二提供的另一种基于交错网格的声波波动方程正演模拟装置的示范性结构框图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
在后续的描述中,使用用于表示元件的诸如“模块”、“部件”或“单元”的后缀仅为了有利于本发明的说明,其本身并没有特定的意义。因此,"模块"与"部件"可以混合地使用。
实施例
如图1所示,在本实施例中,一种基于交错网格的声波波动方程正演模拟方法,包括:
S10、获取地震参数;
S20、建立基于交错网格的声波波动方程;
S30、采用时空域有限差分法计算所述声波波动方程的频散关系式;
S40、根据所述频散关系式获取波场模拟所满足的稳定条件;
S50、采用吸收边界条件对声波波动方程进行波场延拓,获取波场及地震记录。
在本实施例中,通过交错网格建立新的差分结构,求取的差分系数能够使声波在更大的波数范围内压制数值频散,进一步提高了声波波动方程的模拟精度。
在本实施例中,所述地震参数包括:正演模拟所需的速度场文件、差分算子阶数、震源函数及其主频、正演所采用的时间与空间步长及地震记录时长,海绵吸收边界条件的参数。
在本实施例中,所述建立基于交错网格的声波波动方程包括:
时间偏导数采用二阶差分离散;其公式为:
其中,ρ代表介质密度,u和w分别代表水平和垂直方向的速度波场,x,z为直角坐标系的两个轴,t为时间,v为地震波的传播速度,p代表声波波场;
运用二阶差分对时间偏导数进行离散:
在本实施例中,空间偏导数的公式为:
如图2所示,在本实施例中,所述步骤S30包括:
S31、计算声波波动方程的离散形式;
S32、获取时空域有限差分法与LAX-WENDROFF方案的关系式;
S33、采用褶积微分算子法求解所述差分系数;
S34、利用平面波理论,得到时空域频散关系式。
在本实施例中,所述声波波动方程的离散形式计算公式为:
在本实施例中,新差分结构的时空域有限差分法与LAX-WENDROFF方案的关系式为:
其中,px、pxxx、pxxxxx为声波波场p关于x的1阶、3阶、5阶偏导数;pxzz为声波波场p关于x和z的混合偏导数;r=vΔt/h。
在本实施例中,褶积微分算子的基本原理是:将偏导数算子的傅里叶变换进行傅里叶反变换得到的;本实施例推导得到的差分系数表达式为:
对声波波动方程的离散形式计算公式的两边进行傅里叶变换,得到声波波动方程基于新差分结构的时空域有限差分法的频散关系表达式:
由差分系数表达式得到新差分结构的时空域有限差分法的频散误差表达式为:
如图3所示,为频散误差对比示意图,其中,图3(a)和图3(b)为传统有限差分法,图3(c)和图3(d)为本发明方法;图中纵坐标为数值相速度与真实速度比值,如果该比值接近于1,则频散误差小,反之频散误差大,图中横坐标表示归一化波数,算子长度2M=16,如图3(a)所示,当r=0.2,传统有限差分法仍然存在明显的频散误差;当r=0.4(实际上是时间步长Δt增大一倍)传统有限差分法的频散显著增加;分别对比图3(a)和图3(c)、图3(b)和图3(d)可以看出本发明方法可以在更大的波数范围内保持频散误差小,特别是在高波数情况下,表明本发明方法的模拟精度优于传统有限差分法。
由频散关系得到本实施例的有限差分法满足:
r2G≤1
取那奎思波特波数(kx,kz)=(π/h,π/h)代入频散误差表达式得到稳定条件为:
采用海绵吸收边界条件对边界反射进行吸收,利用求取的差分系数进行声波波场外推,得到任意时刻的波场及整个地震记录。记录波场快照切片,输出地震记录。
图4为本发明方法与传统有限差分法稳定性对比图,图中纵坐标为稳定性因子,该值越大表明方法更稳定,横坐标为差分算子长度2M,对比发现,当差分算子长度相同时,本方明方法的稳定性条件均大于传统有限差分法,这意味着本发明方法稳定性更好,可以采用更大的时间步长进行波场延拓。
图5为均匀速度模型中1s时刻的波场快照图。介质的地震波传播速度为2000m/s,模型网格区域为401×401点,网格的纵横空间点距为20m,震源位置在模型正中央,震源主频15Hz,差分算子的长度2M=16。对比图5(a)和图5(b)可知,当r由0.2增至0.4时(实际上是时间步长Δt由2ms增大至4ms),传统有限差分法的数值频散显著增强;再对比图5中(a)和(c)、(b)和(d)可知,本发明方法的模拟精度优于传统有限差分法,这与图3中的频散误差曲线分析一致,证实了本发明方法模拟精度更高。
图6为复杂的Marmousi速度模型,模型网格区域为2721×701点,网格的纵横空间点距为10m,震源主频为15Hz,位于模型地表的正中央,差分算子长度2M=16。
图7为Marmousi速度模型中的2.5s时刻的波场快照图,(a)为采用传统有限差分法,(b)为采用本发明方法图,从图7中黑色矩形框的放大图可以明显的看出,本发明方法提高了波场模拟精度,说明了本文发展的方法的有效性与稳健性。
在本实施例中,基于交错网格的差分结构除了包含轴向2M个离散点外,还包含4个非轴向离散点。与传统的差分结构相比,引入了额外的计算量。然而,新差分结构改善了传统时间2阶方法的精度,且通过稳定性分析可知,基于新差分结构的时空域有限差分法可以采用更大的时间步长。从推导出来的差分系数表达式可以看出,差分系数随着介质速度的变化而改变,本实施例的实现策略是预先计算并存储给定速度模型范围内差分系数。例如,对于给定的模型速度范围为1000m/s至5000m/s,本实施例从1000m/s开始以1m/s的增量计算并存储这些差分系数,时间4阶空间16阶差分格式只需要大约0.1MB存储这些差分系数。因此,本实施例的有限差分法和传统时间二阶方案计算内存几乎相等,本方法能够为计算量巨大的逆时偏移成像和全波形反演提供更为准确的波场信息,具有重要的实际用途。
实施例二
如图8所示,在本实施例中,一种基于交错网格的声波波动方程正演模拟装置,包括:
获取模块10,用于获取地震参数;
建模模块20,用于建立基于交错网格的声波波动方程;
差分模块30,用于采用时空域有限差分法计算所述声波波动方程的频散关系式;
模拟模块40,用于根据所述频散关系式获取波场模拟所满足的稳定条件;
延拓模块50,用于采用吸收边界条件对声波波动方程进行波场延拓,获取波场及地震记录。
在本实施例中,通过交错网格建立新的差分结构,求取的差分系数能够使声波在更大的波数范围内压制数值频散,进一步提高了声波波动方程的模拟精度。
在本实施例中,所述地震参数包括:正演模拟所需的速度场文件、差分算子阶数、震源函数及其主频、正演所采用的时间与空间步长及地震记录时长,海绵吸收边界条件的参数。
在本实施例中,所述建立基于交错网格的声波波动方程包括:
时间偏导数采用二阶差分离散;其公式为:
其中,ρ代表介质密度,u和w分别代表水平和垂直方向的速度波场,x,z为直角坐标系的两个轴,t为时间,v为地震波的传播速度,p代表声波波场;
运用二阶差分对时间偏导数进行离散:
空间偏导数的公式为:
在本实施例中,所述采用时空域有限差分法计算所述声波波动方程的频散关系式包括:
计算声波波动方程的离散形式;
获取时空域有限差分法与LAX-WENDROFF方案的关系式;
采用褶积微分算子法求解所述差分系数;
利用平面波理论,得到时空域频散关系式。
在本实施例中,所述声波波动方程的离散形式计算公式为:
在本实施例中,新差分结构的时空域有限差分法与LAX-WENDROFF方案的关系式为:
在本实施例中,褶积微分算子的基本原理是:将偏导数算子的傅里叶变换进行傅里叶反变换得到的;本实施例推导得到的差分系数表达式为:
对声波波动方程的离散形式计算公式的两边进行傅里叶变换,得到声波波动方程基于新差分结构的时空域有限差分法的频散关系表达式:
由差分系数表达式得到新差分结构的时空域有限差分法的频散误差表达式为:
如图3所示,为频散误差对比示意图,其中,图3(a)和图3(b)为传统有限差分法,图3(c)和图3(d)为本发明方法;图中纵坐标为数值相速度与真实速度比值,如果该比值接近于1,则频散误差小,反之频散误差大,图中横坐标表示归一化波数,算子长度2M=16,如图3(a)所示,当r=0.2,传统有限差分法仍然存在明显的频散误差;当r=0.4(实际上是时间步长Δt增大一倍)传统有限差分法的频散显著增加;分别对比图3(a)和图3(c)、图3(b)和图3(d)可以看出本发明方法可以在更大的波数范围内保持频散误差小,特别是在高波数情况下,表明本发明方法的模拟精度优于传统有限差分法。
由频散关系得到本实施例的有限差分法满足:
r2G≤1
取那奎思波特波数(kx,kz)=(π/h,π/h)代入频散误差表达式得到稳定条件为:
采用海绵吸收边界条件对边界反射进行吸收,利用求取的差分系数进行声波波场外推,得到任意时刻的波场及整个地震记录。记录波场快照切片,输出地震记录。
图4为本发明方法与传统有限差分法稳定性对比图,图中纵坐标为稳定性因子,该值越大表明方法更稳定,横坐标为差分算子长度2M,对比发现,当差分算子长度相同时,本方明方法的稳定性条件均大于传统有限差分法,这意味着本发明方法稳定性更好,可以采用更大的时间步长进行波场延拓。
图5为均匀速度模型中1s时刻的波场快照图。介质的地震波传播速度为2000m/s,模型网格区域为401×401点,网格的纵横空间点距为20m,震源位置在模型正中央,震源主频15Hz,差分算子的长度2M=16。对比图5(a)和图5(b)可知,当r由0.2增至0.4时(实际上是时间步长Δt由2ms增大至4ms),传统有限差分法的数值频散显著增强;再对比图5中(a)和(c)、(b)和(d)可知,本发明方法的模拟精度优于传统有限差分法,这与图3中的频散误差曲线分析一致,证实了本发明方法模拟精度更高。
图6为复杂的Marmousi速度模型,模型网格区域为2721×701点,网格的纵横空间点距为10m,震源主频为15Hz,位于模型地表的正中央,差分算子长度2M=16。
图7为Marmousi速度模型中的2.5s时刻的波场快照图,(a)为采用传统有限差分法,(b)为采用本发明方法图,从图7中黑色矩形框的放大图可以明显的看出,本发明方法提高了波场模拟精度,说明了本文发展的方法的有效性与稳健性。
在本实施例中,基于交错网格的差分结构除了包含轴向2M个离散点外,还包含4个非轴向离散点。与传统的差分结构相比,引入了额外的计算量。然而,新差分结构改善了传统时间2阶方法的精度,且通过稳定性分析可知,基于新差分结构的时空域有限差分法可以采用更大的时间步长。从推导出来的差分系数表达式可以看出,差分系数随着介质速度的变化而改变,本实施例的实现策略是预先计算并存储给定速度模型范围内差分系数。例如,对于给定的模型速度范围为1000m/s至5000m/s,本实施例从1000m/s开始以1m/s的增量计算并存储这些差分系数,时间4阶空间16阶差分格式只需要大约0.1MB存储这些差分系数。因此,本实施例的有限差分法和传统时间二阶方案计算内存几乎相等,本方法能够为计算量巨大的逆时偏移成像和全波形反演提供更为准确的波场信息,具有重要的实际用途。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中,包括若干指令用以使得一台终端设备(可以是手机,计算机,服务器,空调器,或者网络设备等)执行本发明各个实施例所述的方法。
以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (6)
1.一种基于交错网格的声波波动方程正演模拟方法,其特征在于,包括:
获取地震参数;
建立基于交错网格的声波波动方程;
采用时空域有限差分法计算所述声波波动方程的频散关系式;
根据所述频散关系式获取波场模拟所满足的稳定条件;
采用吸收边界条件对声波波动方程进行波场延拓,获取波场及地震记录;
其中,所述建立基于交错网格的声波波动方程包括:
时间偏导数采用二阶差分离散;其公式为:
其中,ρ代表介质密度,u和w分别代表水平和垂直方向的速度波场,x,z为直角坐标系的两个轴,t为时间,v为地震波的传播速度,p代表声波波场;
空间偏导数的公式为:
其中,所述采用时空域有限差分法计算所述声波波动方程的频散关系式包括:
计算声波波动方程的离散形式;
获取时空域有限差分法与LAX-WENDROFF方案的关系式;
采用褶积微分算子法求解所述差分系数;
利用平面波理论,得到时空域频散关系式。
2.根据权利要求1所述的一种基于交错网格的声波波动方程正演模拟方法,其特征在于,所述地震参数包括:正演模拟所需的速度场文件、差分算子阶数、震源函数及其主频、正演所采用的时间与空间步长及地震记录时长,海绵吸收边界条件的参数。
4.一种基于交错网格的声波波动方程正演模拟装置,其特征在于,包括:
获取模块,用于获取地震参数;
建模模块,用于建立基于交错网格的声波波动方程;
差分模块,用于采用时空域有限差分法计算所述声波波动方程的频散关系式;
模拟模块,用于根据所述频散关系式获取波场模拟所满足的稳定条件;
延拓模块,用于采用吸收边界条件对声波波动方程进行波场延拓,获取波场及地震记录;
其中,所述建立基于交错网格的声波波动方程包括:
时间偏导数采用二阶差分离散;其公式为:
其中,ρ代表介质密度,u和w分别代表水平和垂直方向的速度波场,x,z为直角坐标系的两个轴,t为时间,v为地震波的传播速度,p代表声波波场;
空间偏导数的公式为:
其中,所述采用时空域有限差分法计算所述声波波动方程的频散关系式包括:
计算声波波动方程的离散形式;
获取时空域有限差分法与LAX-WENDROFF方案的关系式;
采用褶积微分算子法求解所述差分系数;
利用平面波理论,得到时空域频散关系式。
5.根据权利要求4所述的一种基于交错网格的声波波动方程正演模拟装置,其特征在于,所述地震参数包括:正演模拟所需的速度场文件、差分算子阶数、震源函数及其主频、正演所采用的时间与空间步长及地震记录时长,海绵吸收边界条件的参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811354411.5A CN109490956B (zh) | 2018-11-14 | 2018-11-14 | 一种基于交错网格的声波波动方程正演模拟方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811354411.5A CN109490956B (zh) | 2018-11-14 | 2018-11-14 | 一种基于交错网格的声波波动方程正演模拟方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109490956A CN109490956A (zh) | 2019-03-19 |
CN109490956B true CN109490956B (zh) | 2020-12-08 |
Family
ID=65695936
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811354411.5A Active CN109490956B (zh) | 2018-11-14 | 2018-11-14 | 一种基于交错网格的声波波动方程正演模拟方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109490956B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111487677A (zh) * | 2020-03-31 | 2020-08-04 | 深圳市勘察研究院有限公司 | 一种声波波动方程叠前逆时偏移成像方法及装置 |
CN113589362B (zh) * | 2020-04-30 | 2024-03-19 | 中国石油化工股份有限公司 | 三维陆上耦合波正演模拟方法 |
CN114089419B (zh) * | 2020-08-24 | 2024-04-30 | 中国石油化工集团有限公司 | 优化的变网格地震正演模拟方法 |
CN112255675B (zh) * | 2020-10-07 | 2023-02-21 | 长安大学 | 一种地震资料震源波场重构方法、系统、设备、介质及应用 |
CN112526605B (zh) * | 2020-12-24 | 2022-09-02 | 广州海洋地质调查局 | 一种采用地震数值模拟勘探天然气水合物的方法 |
CN113281808B (zh) * | 2021-04-22 | 2023-10-20 | 南方海洋科学与工程广东省实验室(湛江) | 一种抗频散地震波正演方法、系统、装置及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459773A (zh) * | 2014-08-08 | 2015-03-25 | 中国石油天然气集团公司 | 基于交错网格Lowrank分解的无条件稳定地震波场延拓方法 |
WO2016108896A1 (en) * | 2014-12-31 | 2016-07-07 | Landmark Graphics Corporation | Seismic elastic wave simulation for tilted transversely isotropic media using adaptive lebedev staggered grid |
CN105974471A (zh) * | 2016-07-19 | 2016-09-28 | 中国地质大学(北京) | 一种基于异步流的地震数据多gpu快速正演计算方法 |
CN106443793A (zh) * | 2016-11-10 | 2017-02-22 | 中国石油化工股份有限公司 | 一种时空双变正演模拟方法 |
-
2018
- 2018-11-14 CN CN201811354411.5A patent/CN109490956B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104459773A (zh) * | 2014-08-08 | 2015-03-25 | 中国石油天然气集团公司 | 基于交错网格Lowrank分解的无条件稳定地震波场延拓方法 |
WO2016108896A1 (en) * | 2014-12-31 | 2016-07-07 | Landmark Graphics Corporation | Seismic elastic wave simulation for tilted transversely isotropic media using adaptive lebedev staggered grid |
CN105974471A (zh) * | 2016-07-19 | 2016-09-28 | 中国地质大学(北京) | 一种基于异步流的地震数据多gpu快速正演计算方法 |
CN106443793A (zh) * | 2016-11-10 | 2017-02-22 | 中国石油化工股份有限公司 | 一种时空双变正演模拟方法 |
Non-Patent Citations (1)
Title |
---|
声波和弹性波波动方程有限差分正反演方法研究;任志明;《中国博士学位论文全文数据库》;20180228(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109490956A (zh) | 2019-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109490956B (zh) | 一种基于交错网格的声波波动方程正演模拟方法及装置 | |
CN109490955B (zh) | 一种基于规则网格的声波波动方程正演模拟方法及装置 | |
CN104122585B (zh) | 基于弹性波场矢量分解与低秩分解的地震正演模拟方法 | |
Monteiller et al. | Three-dimensional full waveform inversion of short-period teleseismic wavefields based upon the SEM–DSM hybrid method | |
Lombard et al. | Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves | |
CN103139907B (zh) | 一种利用指纹法的室内无线定位方法 | |
Luo et al. | Fast Huygens sweeping methods for Helmholtz equations in inhomogeneous media in the high frequency regime | |
US20120316844A1 (en) | System and method for data inversion with phase unwrapping | |
CN103149585B (zh) | 一种弹性偏移地震波场构建方法及装置 | |
KR20140021584A (ko) | 스펙트럼 형상을 이용하는 전 파동장 역산의 수렴 레이트 | |
CN112596107B (zh) | 弹性参数反演方法及装置 | |
WO2012170201A2 (en) | System and method for seismic data inversion by non-linear model update | |
Long et al. | A temporal fourth-order scheme for the first-order acoustic wave equations | |
Habashy et al. | Source-receiver compression scheme for full-waveform seismic inversion | |
JP7126612B1 (ja) | 物理的プロセスのシミュレーションのための空き空間ドメイン・デコンポジション | |
Ren et al. | Multiscale viscoacoustic waveform inversion with the second generation wavelet transform and adaptive time–space domain finite-difference method | |
CN107102359B (zh) | 地震数据保幅重建方法和系统 | |
Chemingui et al. | Seismic data reconstruction by inversion to common offset | |
Ha et al. | 3D Laplace-domain waveform inversion using a low-frequency time-domain modeling algorithm | |
CN108828659B (zh) | 基于傅里叶有限差分低秩分解的地震波场延拓方法及装置 | |
Bai et al. | Gaussian beam reconstruction of seismic data | |
CN113221395A (zh) | 基于分层介质的地震走时参数化网格模型构建方法及应用 | |
Takahashi et al. | Convolutional equivalent layer for gravity data processing | |
CN110888166B (zh) | 基于l-bfgs算法的最小二乘偏移成像方法及装置 | |
Amodei et al. | Computation of uniform wave forms using complex rays |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |