CN102879816A - 一种地震多次波偏移方法 - Google Patents
一种地震多次波偏移方法 Download PDFInfo
- Publication number
- CN102879816A CN102879816A CN2012102476511A CN201210247651A CN102879816A CN 102879816 A CN102879816 A CN 102879816A CN 2012102476511 A CN2012102476511 A CN 2012102476511A CN 201210247651 A CN201210247651 A CN 201210247651A CN 102879816 A CN102879816 A CN 102879816A
- Authority
- CN
- China
- Prior art keywords
- imaging
- wave
- multiple reflection
- record
- point
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种地震多次波偏移方法,本发明是在逆时偏移方法的基础上进行了改进,将多次波融合到逆时偏移方程中,使得多次波归位到地下正确的位置上。本发明的多次波偏移方法主要分为两步。第一步,利用一些已知的预测多次波的方法得到多次波,第二步,利用逆时偏移方法进行盐下深度域成像。本发明利用多次波信息,并将多次波融合到逆时偏移方程中,利用一次波和多次波的地震记录波场做成像,与传统技术相比,本发明方可以实现多次波盐下成像,达到更大范围、更加均衡的地下照明度,极大地提高了盐下成像的分辨率。本发明可以广泛用于油气勘探领域中,特别是对盐下结构成像的效果更加明显。
Description
技术领域
本发明涉及一种地下结构成像方法,特别是关于一种用于油气勘探领域的地震多次波偏移方法。
背景技术
几十年来,地质学家们认为盐体可以作为一些油气藏的天然盖层,其下方可能存在超大规模的油气藏。然而盐体是高速度地质体,这些地质体像声音透镜一样,能使地震波信号发散或者聚集,造成盐下成像面临很大的挑战。近期国际上提出了很多针对高速盐体问题地震成像技术,其中很多方法可以提高成像质量。然而上述方法都是基于反射波成像,在复杂盐下区域,反射波照明范围小,照明度低,很难获得满意的盐下成像结果。
有别于反射波,多次波最有可能实现盐下成像,多次波是经一次以上的地下界面反射或海水面反射,最终在地面被接收的地震波。一般来说,与一次反射波相比,多次波在地下传播的射线路径更长,覆盖的区域更广。多次波可以传播到地球内部,可照明到一次反射波无法到达的阴影区。另外,多次波的反射角一般比一次反射波小,能提供很好的地球内部结构信息。传统的成像方法只对一次波偏移,把多次波视为噪声,而多次波偏移则利用多次波信息,提供地下界面构造成像。从多次波中直接提取有用信息的观点可以追溯到Claerbout(克莱尔博特)在地震干涉测量方面的研究,该研究说明了如何通过被动震源激发产生的自相关道,获得地球表面格林函数(例如一个脉冲点源)。被动震源的空间位置和激发时间未知,所以可以认为是一种被动地震方法。基于这种虚震源方式,可以运用运动学方法将多次波转换为由地表虚震源生成的一次波。多次波偏移,最直观的方法就是将多次波转换为反射波,然后采用传统偏移方法。如运用地震干涉测量方法,将三维垂直地震剖面的地震多次波数据转换成反射波数据。另一种可以用到的多次波偏移方法是最小平方偏移,其是利用海底电缆数据对一次反射波和鬼波同时成像。还有一种是利用反馈模型和波场延拓反演,将多次波转换为一次反射波实现多次波成像。
然而,很多研究证明,将多次波直接转换成反射波,只是对垂直地震数据成立,对传统的表面采集的地震数据并不成功。
发明内容
针对上述问题,本发明的目的是提供一种通过改进传统逆时偏移方法,将多次波归位到地下正确位置的地震多次波偏移方法。
为实现上述目的,本发明采取以下技术方案:一种地震多次波偏移方法,包括以下步骤:1)对所有炮的地震记录,先用相关多次波预测方法,预测表面的多次波和层间的多次波,用最小二乘法匹配预测的多次波;2)将每一单炮中任一点P的检波器记录到的地震记录作正向传播,将另一点Mg的检波器记录到的一阶多次波记录作反向传播;3)对从P点正向传播的地震记录,通过双程波波动方程,计算出在地下任何一个位置的地震记录正向传播的格林函数P(g);4)对从Mg点反向传播的一阶多次波记录,通过双程波波动方程计算出在地下任何一个位置的一阶多次波反向传播的格林函数Mg(g);5)在地下任何一个网格点,对P(g)与Mg(g)作互相关,形成可能的成像点;6)如果在可能的成像点确实有地下反射界面,步骤5)中的两个波场的格林函数互相关,就会在该地下反射界面成像;7)对Mg点的一阶多次波记录正向传播,计算出在地下任何一个位置的一阶多次波记录正向传播的格林函数Mg(g)’;8)对再一点Mr的检波器记录的二阶多次波记录反向传播,通过双程波波动方程计算出在地下任何一个位置的二阶多次波逆时传播的格林函数Mr(g);9)在地下任何一个网格点,对Mg(g)’与Mr(g)作互相关,形成可能的成像点;10)如果在可能的成像点确实有地下反射界面,步骤9)中的两个波场的格林函数的互相关就会在该地下反射界面成像;11)对完成步骤6)和步骤10)的单炮成像结果进行拉普拉斯滤波,以消除成像过程中的低频成像噪音;12)对所有单炮记录重复步骤2)至步骤11),将所有单炮成像结果进行叠加,得到最终成像结果。
本发明由于采取以上技术方案,其具有以下优点:1、本发明方法用包含一次反射波,表层多次波和层间多次波的地震记录代替震源子波,用预测多次波代替一次反射波作为输入数据,将地表多次波的记录逆时延拓到各个深度层,将包含一次波和多次波的地震记录沿时间延拓到相同深度层,通过互相关成像条件,使得多次波归位到地下正确的位置上。2、本发明方法采用的多次反射波的传播路径与传统的一次反射波传播路径不同,本发明方法可以提高盐下油气藏的照明度。与传统的一次反射波偏移比较,多次波偏移可提供更好的盐体成像。3、本发明利用多次波信息,并将多次波融合到逆时偏移方程中,利用一次波和多次波的地震记录波场做成像,与传统技术相比,本发明方可以实现多次波盐下成像,达到更大范围、更加均衡的地下照明度,极大地提高了盐下成像的分辨率。本发明可以广泛用于油气勘探领域中,特别是对盐下结构成像的效果更加明显。
附图说明
图1是本发明原理图
图2是一个三层速度模型
图3是多次波数据200炮中的一炮
图4是多次波偏移结果
图5是著名的验证各种偏移方法的Sigsbee2B盐丘模型
图6是采用传统反射波偏移方法得到的结果
图7是采用本发明多次波偏移方法得到的结果
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明提出的地震多次波偏移方法,不需要将多次波变换为一次反射波,而是用标准偏移算子处理表面数据。偏移成像方法一般为Kirchhoff(克希霍夫)成像、单程波波动方程和逆时偏移等。其中,Kirchhoff成像方法最为常用,其是假定波动方程的解是高频近似解,对速度模型不敏感,但是对横向速度变化巨大的复杂地质体成像不够准确,存在不容易处理多走时问题和横向速度变化等缺陷。单程波波动方程方法,计算效率快,可以对横向变化剧烈的构造成像,但是对倾角有限制,不能处理大倾角断层。基于双程波波动方程的逆时偏移方法(简称逆时偏移方法)理论上效果最好,既能对横向速度变化剧烈的构造成像,也能对大倾角地层成像,但是计算量大。本发明是在逆时偏移方法的基础上进行了改进,将多次波融合到逆时偏移方程中,使得多次波归位到地下正确的位置上。本发明的多次波偏移方法主要分为两步。第一步,利用一些已知的预测多次波的方法得到多次波,如Radon(拉东变换)变换、逆散射层间多次波和表层相关多次波预测方法等。第二步,利用逆时偏移方法进行盐下深度域成像。用包含一次反射波,表层多次波和层间多次波的地震记录代替震源子波,用多次波代替一次反射波。在逆时偏移处理流程中,将地表多次波的记录逆时延拓到各个深度层,将包含一次波和多次波的地震记录沿时间延拓到相同深度层。成像条件是在各个深度层上自相关两个波场。
具体说,本发明提供的地震多次波偏移方法,包括以下步骤:
1)对所有炮的地震记录,先用相关多次波预测方法,预测表面的多次波和层间的多次波,用最小二乘法匹配预测的多次波;
2)如图1所示,图中P点、Mg点和Mr点均为单炮记录中某一检波点,将在P点的检波器记录到的地震记录作正向传播,将在Mg点的检波器记录到的一阶多次波记录作反向传播;
3)对从P点正向传播的地震记录,通过双程波波动方程,计算出在地下任何一个位置的地震记录正向传播的格林函数P(g);
4)对从Mg点反向传播的一阶多次波记录,通过双程波波动方程计算出在地下任何一个位置的一阶多次波反向传播的格林函数Mg(g);
5)在地下任何一个网格点,对P(g)与Mg(g)作互相关,形成可能的成像点;
6)如果在地下X2,X3和X4点确实有地下反射界面,步骤5)中的两个波场的格林函数互相关就会在该地下反射界面成像;
7)对Mg点的一阶多次波记录正向传播,计算出在地下任何一个位置的一阶多次波记录正向传播的格林函数Mg(g)’;
8)对Mr点的二阶多次波记录反向传播,通过双程波波动方程计算出在地下任何一个位置的二阶多次波逆时传播的格林函数Mr(g);
9)在地下任何一个网格点,对Mg(g)’与Mr(g)作互相关,形成可能的成像点;
10)如果在地下X5点确实有地下反射界面,步骤9)中的两个波场的格林函数的互相关就会在该地下反射界面成像;
11)对步骤6)和步骤10)得到的单炮成像结果进行Laplace(拉普拉斯)滤波,以消除成像过程中的低频成像噪音;
12)对所有单炮记录重复步骤2)至步骤11),并将所有单炮成像结果进行叠加,得到最终成像结果。
为验证本发明方法的可行性和有效性,下面举两个实例
实例1:
如图2所示,对于一地下3层速度模型,采用有限差分声波正演模拟,可以得到单炮带有多次波信号的地震记录(如图3所示),如果采用本发明的多次波偏移方法(如图4所示),可以看到能量最强的两个同像轴C和D,与3层速度模型的速度界面位置和形态相吻合,证明了采用本发明方法能够得到正确的成像,同时也证明了本发明的可行性。
实例2:
如图5所示,是著名的Sigsbee2B盐丘模型,该模型含有高速盐丘体和单个速度异常体,模型复杂,是验证各种偏移方法成像效果的国际标准模型之一。如图6所示,采用传统逆时偏移方法只利用一次反射波成像,会将多次波当作噪音消去,同时由于检波点范围的限制,一次反射波在盐下照明度不强,会导致盐下成像效果不好,如图中椭圆虚线标出的区域,反射波无法形成。如图7所示,由于本发明方法利用了多次波相对一次反射波走时长,路径丰富的特征,因此可以看到图中椭圆虚线标示的盐下区域成像清楚,效果明显区别并优于传统技术。
上述描述和各实施例仅用于说明本发明,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (1)
1.一种地震多次波偏移方法,包括以下步骤:
1)对所有炮的地震记录,先用相关多次波预测方法,预测表面的多次波和层间的多次波,用最小二乘法匹配预测的多次波;
2)将每一单炮中任一点P的检波器记录到的地震记录作正向传播,将另一点Mg的检波器记录到的一阶多次波记录作反向传播;
3)对从P点正向传播的地震记录,通过双程波波动方程,计算出在地下任何一个位置的地震记录正向传播的格林函数P(g);
4)对从Mg点反向传播的一阶多次波记录,通过双程波波动方程计算出在地下任何一个位置的一阶多次波反向传播的格林函数Mg(g);
5)在地下任何一个网格点,对P(g)与Mg(g)作互相关,形成可能的成像点;
6)如果在可能的成像点确实有地下反射界面,步骤5)中的两个波场的格林函数互相关,就会在该地下反射界面成像;
7)对Mg点的一阶多次波记录正向传播,计算出在地下任何一个位置的一阶多次波记录正向传播的格林函数Mg(g)’;
8)对再一点Mr的检波器记录的二阶多次波记录反向传播,通过双程波波动方程计算出在地下任何一个位置的二阶多次波逆时传播的格林函数Mr(g);
9)在地下任何一个网格点,对Mg(g)’与Mr(g)作互相关,形成可能的成像点;
10)如果在可能的成像点确实有地下反射界面,步骤9)中的两个波场的格林函数的互相关就会在该地下反射界面成像;
11)对完成步骤6)和步骤10)的单炮成像结果进行拉普拉斯滤波,以消除成像过程中的低频成像噪音;
12)对所有单炮记录重复步骤2)至步骤11),将所有单炮成像结果进行叠加,得到最终成像结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210247651.1A CN102879816B (zh) | 2012-07-17 | 2012-07-17 | 一种地震多次波偏移方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210247651.1A CN102879816B (zh) | 2012-07-17 | 2012-07-17 | 一种地震多次波偏移方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102879816A true CN102879816A (zh) | 2013-01-16 |
CN102879816B CN102879816B (zh) | 2015-08-05 |
Family
ID=47481209
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210247651.1A Active CN102879816B (zh) | 2012-07-17 | 2012-07-17 | 一种地震多次波偏移方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102879816B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926623A (zh) * | 2014-05-06 | 2014-07-16 | 王维红 | 一种压制逆时偏移低频噪音的方法 |
CN103984023A (zh) * | 2014-05-15 | 2014-08-13 | 中国科学院地质与地球物理研究所 | 一种地震噪声去除方法及装置 |
CN104280768A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | 一种适用于逆时偏移的吸收边界条件方法 |
CN104570082A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于格林函数表征的全波形反演梯度算子的提取方法 |
CN105223607A (zh) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | 基于反稳相滤波的绕射波成像方法 |
CN105974470A (zh) * | 2016-07-04 | 2016-09-28 | 中国石油大学(华东) | 一种多分量地震资料最小二乘逆时偏移成像方法及系统 |
CN107402405A (zh) * | 2016-05-18 | 2017-11-28 | 中国石油化工股份有限公司 | 静相位虚源道集构建方法 |
CN108415073A (zh) * | 2018-03-06 | 2018-08-17 | 中国科学院测量与地球物理研究所 | 角度域逆散射偏移成像方法及装置 |
CN110858000A (zh) * | 2018-08-24 | 2020-03-03 | 中国石油天然气股份有限公司 | 地震数据重构方法、装置、计算机设备及存储介质 |
CN111480097A (zh) * | 2017-12-15 | 2020-07-31 | 沙特阿拉伯石油公司 | 用于解释员的盐下成像工具 |
CN112748463A (zh) * | 2020-12-21 | 2021-05-04 | 中国科学技术大学 | 一种基于深度学习照明分析的局部偏移成像方法 |
CN115993650A (zh) * | 2023-03-22 | 2023-04-21 | 中国石油大学(华东) | 一种基于棱柱波的地震干涉成像方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101017204A (zh) * | 2006-02-09 | 2007-08-15 | Pgs地球物理公司 | 三维双向声波动方程叠前成像系统和方法 |
-
2012
- 2012-07-17 CN CN201210247651.1A patent/CN102879816B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101017204A (zh) * | 2006-02-09 | 2007-08-15 | Pgs地球物理公司 | 三维双向声波动方程叠前成像系统和方法 |
Non-Patent Citations (4)
Title |
---|
OONG K. YOUN ET AL.: "Depth imaging with multiples", 《GEOPHYSICS》 * |
李鹏,等: "多次波问题的研究进展", 《地球物理学进展》 * |
郭书娟,等: "基于广义的炮偏移方法实现地表多次波和一次波联合成像", 《地球物理学报》 * |
陈聪: "叠前逆时偏移及成像", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104280768A (zh) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | 一种适用于逆时偏移的吸收边界条件方法 |
CN104280768B (zh) * | 2013-07-12 | 2017-03-15 | 中国石油天然气集团公司 | 一种适用于逆时偏移的吸收边界条件方法 |
CN104570082A (zh) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于格林函数表征的全波形反演梯度算子的提取方法 |
CN104570082B (zh) * | 2013-10-29 | 2017-04-19 | 中国石油化工股份有限公司 | 一种基于格林函数表征的全波形反演梯度算子的提取方法 |
CN103926623A (zh) * | 2014-05-06 | 2014-07-16 | 王维红 | 一种压制逆时偏移低频噪音的方法 |
CN103926623B (zh) * | 2014-05-06 | 2015-02-18 | 东北石油大学 | 一种压制逆时偏移低频噪音的方法 |
CN103984023B (zh) * | 2014-05-15 | 2017-01-25 | 中国科学院地质与地球物理研究所 | 一种地震噪声去除方法及装置 |
CN103984023A (zh) * | 2014-05-15 | 2014-08-13 | 中国科学院地质与地球物理研究所 | 一种地震噪声去除方法及装置 |
CN105223607A (zh) * | 2014-07-04 | 2016-01-06 | 中国石油化工股份有限公司 | 基于反稳相滤波的绕射波成像方法 |
CN107402405A (zh) * | 2016-05-18 | 2017-11-28 | 中国石油化工股份有限公司 | 静相位虚源道集构建方法 |
CN105974470B (zh) * | 2016-07-04 | 2017-06-16 | 中国石油大学(华东) | 一种多分量地震资料最小二乘逆时偏移成像方法及系统 |
CN105974470A (zh) * | 2016-07-04 | 2016-09-28 | 中国石油大学(华东) | 一种多分量地震资料最小二乘逆时偏移成像方法及系统 |
CN111480097A (zh) * | 2017-12-15 | 2020-07-31 | 沙特阿拉伯石油公司 | 用于解释员的盐下成像工具 |
CN108415073A (zh) * | 2018-03-06 | 2018-08-17 | 中国科学院测量与地球物理研究所 | 角度域逆散射偏移成像方法及装置 |
CN108415073B (zh) * | 2018-03-06 | 2019-09-13 | 中国科学院测量与地球物理研究所 | 角度域逆散射偏移成像方法及装置 |
CN110858000A (zh) * | 2018-08-24 | 2020-03-03 | 中国石油天然气股份有限公司 | 地震数据重构方法、装置、计算机设备及存储介质 |
CN110858000B (zh) * | 2018-08-24 | 2021-07-02 | 中国石油天然气股份有限公司 | 地震数据重构方法、装置、计算机设备及存储介质 |
CN112748463A (zh) * | 2020-12-21 | 2021-05-04 | 中国科学技术大学 | 一种基于深度学习照明分析的局部偏移成像方法 |
CN115993650A (zh) * | 2023-03-22 | 2023-04-21 | 中国石油大学(华东) | 一种基于棱柱波的地震干涉成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102879816B (zh) | 2015-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102879816B (zh) | 一种地震多次波偏移方法 | |
CN106597533B (zh) | 一种用于山前带地震资料处理的深度域速度建模方法 | |
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
CN104570125A (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN103308943B (zh) | 一种海洋地震资料处理中层间多次波衰减的方法及装置 | |
CN105093292A (zh) | 一种地震成像的数据处理方法和装置 | |
US20120116680A1 (en) | Microseismic source location estimation method with high resolution using green's functions | |
CN105093301B (zh) | 共成像点反射角角道集的生成方法及装置 | |
CN104533396A (zh) | 一种远探测声波的处理方法 | |
CN104360381A (zh) | 一种地震资料的偏移成像处理方法 | |
CN106125139B (zh) | 一种三维地震数据处理方法及系统 | |
CN105301636A (zh) | 速度模型的建立方法和装置 | |
CN103913768A (zh) | 基于地震波资料对地表中浅层进行建模的方法及装置 | |
CN105093296A (zh) | 一种优化观测系统的方法及装置 | |
CN101738636A (zh) | 一种三维vsp高斯束法多波联合偏移成像方法 | |
CN102520443B (zh) | 一种绕射波叠前偏移的方法 | |
CN105572735B (zh) | 一种提高叠前深度偏移成像精度的方法及装置 | |
Song et al. | Experiment of 3D seismic reflection technique for forward probing on TBM tunnel Face | |
CA2875984C (en) | Seismic imaging | |
CN115236730B (zh) | 一种层间多次波傅里叶有限差分的地震波场偏移成像方法 | |
CN104267434A (zh) | 一种三维多分量地震观测系统炮检距分布获取方法及装置 | |
Wu et al. | Offset VSP P wave and converted S wave salt proximity migration study | |
CN109765617B (zh) | 基于运动学的垂向双震源压制多次反射折射波的方法 | |
Buur et al. | Salt interpretation enabled by reverse-time migration | |
CN112147686B (zh) | 多期发育火成岩叠前深度偏移成像速度建模方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |