CN104459780A - 一种利用波动方程获取地震波路径的方法 - Google Patents
一种利用波动方程获取地震波路径的方法 Download PDFInfo
- Publication number
- CN104459780A CN104459780A CN201410747999.6A CN201410747999A CN104459780A CN 104459780 A CN104459780 A CN 104459780A CN 201410747999 A CN201410747999 A CN 201410747999A CN 104459780 A CN104459780 A CN 104459780A
- Authority
- CN
- China
- Prior art keywords
- frequency
- seismic
- seismic wave
- point
- path
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种利用波动方程获取地震波路径的方法,其包括以下步骤:1)给定初始参数;2)将待勘探的复杂非均匀地下介质模型划分为网状格点;3)将震源置于模型的震源点S处,逐点计算整个空间内所求频段内某个频率下的地震波场;4)将震源置于模型的接收点R处,逐点计算整个空间内所求频段内某个频率下的地震波场;5)在频率域内,将同一频率下在震源点S和接收点R处计算的地震波场相乘并与前一频率下得到的相乘结果进行累加;6)在所求频段内进行频率循环,若频率循环完成,则存储地震波路径。本发明避免了传统方法的高频近似和多路径问题,可以广泛应用于石油勘探复杂构造的地震模拟、地震照明和地震偏移技术中。
Description
技术领域
本发明涉及一种石油勘探领域中的地震波路径的获取方法,特别是关于一种利用波动方程获取地震波路径的方法。
背景技术
目前,地震波的路径通常由Kirchhoff(克希霍夫)理论计算得到,其进行高频近似,即认为地震波能量集中在地下介质两点间的某条线上(称为“射线”)。然而,实际上地震波是有限频带的信号。根据Fermat(费马特)原理,地震波沿着旅行时最小的路径传播,即地震波从震源点S传播到接收点R所用时间最短的路径。射线法借鉴最短路径算法求出震源和接收器两点间最短的旅行时。实际中,由于地下介质十分复杂,且计算地震波路径时需要对地下介质进行离散,即将地下介质划分成一系列的小单元并形成网格,因此,计算射线路径的精度与网格的大小或稀疏程度有密切的关系。另外,当地下介质较为复杂时,射线路径的计算会出现多路径的问题,即从震源点S到接收点R地震波旅行时最小的路径有两条或以上。精确计算地震波的传播路径对于理解地震波在复杂介质中的传播规律以及地震偏移成像有着十分重要的意义。
如前所述,计算地震波传播路径的方法大多基于Kirchhoff理论,其主要存在高频近似(实际中地震波为有限频带信号)和多路径的问题。随着我国油气勘探的不断深入和勘探程度的不断提高,石油勘探的目标已逐步转向复杂隐蔽储层。这些复杂沉积环境下的储层具有厚度薄、岩性高度非均质、有效储层分布分散、岩石物理关系复杂和储层岩性对比弱等特征。传统的地震波路径计算方法无法满足复杂构造勘探的要求。
发明内容
针对上述问题,本发明的目的是提供一种精确度更高、适用不同频段地震波的可广泛应用于复杂构造勘探的利用波动方程计算地震波路径的实现方法。
为实现上述目的,本发明采取以下技术方案:一种利用波动方程获取地震波路径的方法,包括以下步骤:1)给定初始参数,包括待勘探的复杂非均匀地下介质模型、地震波场的频率以及震源点S和接收点R的空间位置;2)将待勘探的复杂非均匀地下介质模型,以地表为起始点向着待勘探的方向按照不同的深度间隔和水平间隔划分为网状格点;3)将震源置于地下介质模型的震源点S处,利用波动方程逐点计算整个空间内所求频段内某个频率下震源点S处的地震波场uS(x,ω);4)将震源置于地下介质模型的接收点R处,利用波动方程逐点计算整个空间内所求频段内某个频率下接收点R处的地震波场uR(x,ω);5)在频率域内,将同一频率下震源点S处的地震波场uS(x,ω) 和接收点R处的地震波场uR(x,ω)相乘并与前一频率下得到的相乘结果进行累加;6)在所求频段内进行频率循环,即重复步骤3)-5):若频率循环未完成,即未累加上所求频段的最后一个频率值,则返回步骤3)继续进行计算;若频率循环完成,即已累加上所求频段的最后一个频率值,则存储地震波路径。
所述步骤3)和4)中,所使用的波动方程为以下二者之一:
其中,x=(x,y,z)为待勘探的复杂非均匀地下介质模型某一点处的空间位置,ω为地震波场的频率,u(x,ω)为x点处的频率域地震波场,为拉普拉斯算子,k为波数,kz为垂直波数,i为虚数单位。
所述步骤6)中,地震波路径为
本发明由于采取以上技术方案,其具有以下优点:1、本发明由于通过频率域波动方程计算不同位置的地震波路径,与常规的射线方法相比,能够更加准确的获得地震波路径;2、本发明方法得出的地震波路径由于是基于波动方程,计算精度高,克服了常规方法多路径的固有缺点;3、本发明方法由于在频率域中进行数值实施,可计算不同频率/频段的地震波路径,克服了常规方法高频近似的固有缺点。4、本发明可以广泛应用于石油勘探复杂构造的地震模拟、地震照明和地震偏移技术中。
附图说明
图1是本发明的技术实现流程示意图
图2是本发明计算的0-20Hz频段的地震波路径
图3是本发明计算的20-30Hz频段的地震波路径
图4是本发明计算的40-50Hz频段的地震波路径
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明提供的利用波动方程获取地震波路径的方法,其包括以下步骤:
1)给定初始参数,包括待勘探的复杂非均匀地下介质模型、地震波场的频率以及震源点S和接收点R的空间位置。
2)将待勘探的复杂非均匀地下介质模型,以地表为起始点向着待勘探的方向按照不同的深度间隔和水平间隔划分为网状格点。
3)将震源置于地下介质模型的震源点S处,并利用下述波动方程(1)或(2)逐点计 算整个空间内所求频段内某个频率下震源点S处的地震波场uS(x,ω)。
其中,x=(x,y,z)为待勘探的复杂非均匀地下介质模型某一点处的空间位置,ω为地震波场的频率,u(x,ω)为x点处的频率域地震波场,为Laplace(拉普拉斯)算子,k为波数,kz为垂直波数,i为虚数单位。
4)将震源置于地下介质模型的接收点R处,并利用上述波动方程(1)或(2)逐点计算整个空间内所求频段内某个频率下接收点R处的地震波场uR(x,ω)。
5)在频率域内,将同一频率下震源点S处的地震波场uS(x,ω)和接收点R处的地震波场uR(x,ω)相乘并与前一频率下得到的相乘结果进行累加。
6)在所求频段内(如0-20Hz)进行频率循环,即重复步骤3)-5)。若频率循环未完成,即未累加上所求频段的最后一个频率值,则返回步骤3)继续进行计算;若频率循环完成,即已累加上所求频段的最后一个频率值,则存储地震波路径 图2显示了0-20Hz频段的地震波路径(其中震源点S与接收点R之间的直线为射线路径)。
图3和图4分别显示了利用本发明方法获取的20-30Hz和40-50Hz频段的地震波路径,这进一步说明本发明方法可以在不同频段内精确计算地震波路径。从图中可以看出,采用本发明计算的地震波路径避免了传统方法的高频近似和多路径问题。
上述各实施例仅用于说明本发明,其中的一些方法和步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (3)
1.一种利用波动方程获取地震波路径的方法,包括以下步骤:
1)给定初始参数,包括待勘探的复杂非均匀地下介质模型、地震波场的频率以及震源点S和接收点R的空间位置;
2)将待勘探的复杂非均匀地下介质模型,以地表为起始点向着待勘探的方向按照不同的深度间隔和水平间隔划分为网状格点;
3)将震源置于地下介质模型的震源点S处,利用波动方程逐点计算整个空间内所求频段内某个频率下震源点S处的地震波场uS(x,ω);
4)将震源置于地下介质模型的接收点R处,利用波动方程逐点计算整个空间内所求频段内某个频率下接收点R处的地震波场uR(x,ω);
5)在频率域内,将同一频率下震源点S处的地震波场uS(x,ω)和接收点R处的地震波场uR(x,ω)相乘并与前一频率下得到的相乘结果进行累加;
6)在所求频段内进行频率循环,即重复步骤3)-5):若频率循环未完成,即未累加上所求频段的最后一个频率值,则返回步骤3)继续进行计算;若频率循环完成,即已累加上所求频段的最后一个频率值,则存储地震波路径。
2.如权利要求1所述的一种利用波动方程获取地震波路径的方法,其特征在于:所述步骤3)和4)中,所使用的波动方程为以下二者之一:
其中,x=(x,y,z)为待勘探的复杂非均匀地下介质模型某一点处的空间位置,ω为地震波场的频率,u(x,ω)为x点处的频率域地震波场,为拉普拉斯算子,k为波数,kz为垂直波数,i为虚数单位。
3.如权利要求1或2所述的一种利用波动方程获取地震波路径的方法,其特征在于:所述步骤6)中,地震波路径为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410747999.6A CN104459780B (zh) | 2014-12-09 | 2014-12-09 | 一种利用波动方程获取地震波路径的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410747999.6A CN104459780B (zh) | 2014-12-09 | 2014-12-09 | 一种利用波动方程获取地震波路径的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104459780A true CN104459780A (zh) | 2015-03-25 |
CN104459780B CN104459780B (zh) | 2017-02-22 |
Family
ID=52906125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410747999.6A Active CN104459780B (zh) | 2014-12-09 | 2014-12-09 | 一种利用波动方程获取地震波路径的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104459780B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106569265A (zh) * | 2015-10-12 | 2017-04-19 | 中国石油化工股份有限公司 | 基于波动方程目的层逆向照明的vsp观测方法及系统 |
CN106896402A (zh) * | 2017-03-21 | 2017-06-27 | 中国矿业大学(北京) | 基于地质单元体的地震正演方法及装置 |
CN109490947A (zh) * | 2018-10-16 | 2019-03-19 | 中国科学院地质与地球物理研究所 | 一种高温介质地震波传播模拟方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
US20120300585A1 (en) * | 2011-05-27 | 2012-11-29 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
CN103675894A (zh) * | 2013-12-24 | 2014-03-26 | 中国海洋石油总公司 | 一种基于三维高斯束射线追踪和频率域合成地震记录方法 |
WO2014151519A1 (en) * | 2013-03-15 | 2014-09-25 | Westerngeco Llc | Frequency-sparse seismic data acquisition and processing |
-
2014
- 2014-12-09 CN CN201410747999.6A patent/CN104459780B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
US20120300585A1 (en) * | 2011-05-27 | 2012-11-29 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
WO2014151519A1 (en) * | 2013-03-15 | 2014-09-25 | Westerngeco Llc | Frequency-sparse seismic data acquisition and processing |
CN103675894A (zh) * | 2013-12-24 | 2014-03-26 | 中国海洋石油总公司 | 一种基于三维高斯束射线追踪和频率域合成地震记录方法 |
Non-Patent Citations (4)
Title |
---|
任浩然 等: ""沿射线路径的波动方程延拓吸收与衰减补偿方法"", 《石油物探》 * |
管西竹 等: ""复杂地表边界元-体积元波动方程数值模拟"", 《地球物理学报》 * |
蒋韬 等: ""基于波动方程的目标导向观测系统设计方法研究"", 《地球物理学进展》 * |
陈生昌 等: ""地震波传播路径追踪的波动方程方法"", 《计算物理》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106569265A (zh) * | 2015-10-12 | 2017-04-19 | 中国石油化工股份有限公司 | 基于波动方程目的层逆向照明的vsp观测方法及系统 |
CN106896402A (zh) * | 2017-03-21 | 2017-06-27 | 中国矿业大学(北京) | 基于地质单元体的地震正演方法及装置 |
CN109490947A (zh) * | 2018-10-16 | 2019-03-19 | 中国科学院地质与地球物理研究所 | 一种高温介质地震波传播模拟方法 |
CN109490947B (zh) * | 2018-10-16 | 2020-01-10 | 中国科学院地质与地球物理研究所 | 一种高温介质地震波传播模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104459780B (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107203002A (zh) | 反演速度模型及其建立方法和地下结构的像的获得方法 | |
CN103984012B (zh) | 基于叠前高斯束深度偏移的绕射波场分离方法 | |
CN102636809B (zh) | 一种传播角度域共成像点道集的生成方法 | |
US11221429B2 (en) | Coherent transmit and receiver bi-static electromagnetic geophysical tomography | |
CN103573250A (zh) | 一种计算水平井井眼到地层上下界面距离的方法 | |
CN103412333B (zh) | 一种静校正基准面确定方法 | |
CN103728660A (zh) | 基于地震数据的多道匹配追踪方法 | |
CN106291725A (zh) | 一种快速反演地下地质体空间位置的方法 | |
CN104459780B (zh) | 一种利用波动方程获取地震波路径的方法 | |
CN105388523A (zh) | 一种高精度的品质因子提取方法 | |
CN104181593B (zh) | 一种三维无射线追踪回折波层析成像方法及装置 | |
CN104459768A (zh) | 一种基于可视化的三维空间目标地质体追踪方法 | |
CN105549077A (zh) | 基于多级多尺度网格相似性系数计算的微震震源定位方法 | |
CN104597488A (zh) | 非等边长网格波动方程有限差分模板优化设计方法 | |
CN105277976A (zh) | 基于岩石露头雷达探测的地震正演模拟方法 | |
Shen | Subsurface focusing measurement of diving waves and its application to reflection tomography | |
CN102385066B (zh) | 一种叠前地震定量成像方法 | |
CN106199694A (zh) | 基于深变子波的合成记录制作方法 | |
CN107942405A (zh) | 预测砂泥岩薄互层砂体累积厚度的方法 | |
CN103777241A (zh) | 基于时间域广义Hilbert变换的三维地震资料快速边缘检测方法 | |
CN114764151B (zh) | 一种大地电磁分频层析反演方法 | |
CN103543467B (zh) | 地震弱信号增强处理方法 | |
CN106908838A (zh) | 构建目标区域地层非弹性衰减品质因子立体模型的方法 | |
CN105572732A (zh) | 一种属性变化率逐次晋级的裂缝发育带检测方法 | |
CN107085236B (zh) | 最大炮检距的确定方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |