CN102749644B - 搜寻反射界面法线角获得最佳成像效果叠前偏移成像方法 - Google Patents
搜寻反射界面法线角获得最佳成像效果叠前偏移成像方法 Download PDFInfo
- Publication number
- CN102749644B CN102749644B CN201210194803.6A CN201210194803A CN102749644B CN 102749644 B CN102749644 B CN 102749644B CN 201210194803 A CN201210194803 A CN 201210194803A CN 102749644 B CN102749644 B CN 102749644B
- Authority
- CN
- China
- Prior art keywords
- angle
- migration
- point
- stack
- cmp
- 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.)
- Expired - Fee Related
Links
- 238000013508 migration Methods 0.000 title claims abstract description 56
- 230000005012 migration Effects 0.000 title claims abstract description 56
- 238000003384 imaging method Methods 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000012545 processing Methods 0.000 description 6
- 101100464782 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) CMP2 gene Proteins 0.000 description 5
- 101100464779 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) CNA1 gene Proteins 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 101000957815 Culex pipiens Alpha-glucosidase Proteins 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于搜寻反射界面法线角获得最佳成像效果的叠前偏移成像方法,该方法首先从O点向一系列CMP点做射线追踪,在O点得到一系列射线与垂直线的交角,即倾角;再将该系列CMP道集上的地震数据在O点按不同倾角范围进行偏移叠加,获得不同角度的叠加能量,比较系列角度对应的偏移叠加能量大小,找到代表最大偏移能量的角度,确认该角度为界面面元法线角F;通过界面面元的法线角F确定最佳偏移叠加图像的数据之共中心点CMP;应用确定的CMP道集和其周围一定范围内道集的地震数据在O点进行叠前偏移成像。该方法相对于传统的叠前偏移成像方法,能够获得更为清晰的图像。
Description
技术领域
本发明属于反射地震勘探领域,特别是涉及利用地震叠前偏移道集数据确定反射界面法线对应的地震数据共中心点位置,从而获得更为清晰的地震偏移图像的方法。
背景技术
对地质构造进行成像的技术,比较常见的的是利用反射地震叠前偏移道集数据进行成像,如果地下介质是水平层状和各向同性的,地下界面法线应该垂直于水平地面,形成最佳成像图的地震数据之共中心点位于法线在地面的投影点。但是介质若为非水平层状,即界面如果是倾斜的,按常规,我们不知道界面某点的法线角和形成该点最佳图像的地震数据之共中心点,所以将所有可能在该点形成反射的地震数据全部累加,并假设噪音是随机的,得到成像结果。但该方法不能确切了解地下界面法线角和反射点对应的共地震数据中心点位置,所以在叠加时将所有可能涉及的数据全部累加,且要求叠加次数很大以便于降低噪音,但实际数据的叠加次数是有限的,很多情况下不能满足需要,应用速度的不准确和非该点反射数据在该点的叠加反而降低了成像的分辨率。
发明内容
在反射地震勘探领域,应用叠前偏移数据确定地下反射界面法线,进而确定对应的地震数据共中心点位置,到目前为止没有有效的方法。而如果在偏移成像过程中,了解地下成像点的界面法线角,进而知道在该点形成最佳叠加图像的地震数据之共中心点,就可以专门对该共中心点附近的地震数据进行偏移叠加成像,从而得到较高分辨率的地震剖面。
本发明根据地下反射点的界面法线角对应的CMP道集数据具有最大成像振幅的原理,对数据进行选择性叠加,仅仅应用少量地震数据成像,从而避免了由于大量非该点反射数据在该点的叠加和速度偏差带来叠加噪音,从而分辨率降低的弊病。
如图1所示,地下倾斜界面上一点O,由两个地震数据组CMP1和CMP2进行成像作业。CMP1道集数据由炮检点对S1-R1,S2-R2组成,其中S1、S2代表放炮点,RI、R2代表检波点。显然,CMP1道集数据不能在O点形成良好的像,原因是S1-R1和S2-R2点对在O点均没有有效的高振幅反射,CMP1与O点的连线不是界面在O点的法线。相反,由于CMP2与O点的连线是界面在O点的法线,CMP2的炮-检点对围绕法线形成完整的反射系列,所以CMP2道集数据在O点可以形成良好的像。如果在成像过程中,我们丢弃CPM1数据而将CMP2数据挑选出来单独成像,成像精度将大大提高,从而避免了大量非O点反射数据叠加带来的分辨率降低的问题。如果能够把地下划分为有限个界元,对每个界元按照CMP2一样,找到其法线,得到围绕法线的完整反射系列,并依照该道集数据分别成像,就能得到高精度的整个地下介质的图像。
在具体检测中,采取的模式是:
假设地下介质为各向同性、速度为V的缓慢变化的介质,炮点S(xs)和检波点R(xr)位于地面,地下反射点O(x,z)所处界面面元的法线方向为F(图2)。其中入射角Ф是指入射波与法线的夹角;t0为O点到地面的双程旅行时。V和t0为实验可得数据。
如图2所示,据DSR方程,α角和β角可以由下式得到[1],即
和
理论上,入射角度Φ应该是两个角度平均值,即
Φ=(-α-β)/2, ……(3)
如图所示,法线角度F为
F=α+Ф ……(4)
在均匀介质情况下,(1)和(2)式的α角和β角可以由下式计算,
和
在实际地震数据处理中,真正的F角和Φ角是难于得到的,只有在确切知道某时刻的信号来源就是从该面元反射的情况下,才能唯一性的得到,而这种情况几乎没有可能。
但是,我们知道,只有在应用(4)式得到的F角与地下界面面元的真实法线一致的情况下,该点的反射振幅最大,偏移叠加得到的能量最强。由此,一个可行的方案是,从一系列CMP点向O点做射线追踪,在O点得到一系列角度;再将一系列CMP道集上的地震数据在该点各自进行偏移叠加,获得偏移后的一系列像能量;比较系列角度对应的偏移叠加能量大小,找到代表最大偏移能量的角度,确认该角度为界面面元法线角F。同时,我们也可以确知,造成最大偏移叠加能量的地震数据的CMP点即为形成最佳偏移叠加成像所对应的地震数据共中心点。
具体步骤为:
(1)首先从O点向一系列CMP点做射线追踪,在O点得到一系列射线与垂直线的交角,即倾角;
(2)再将该系列CMP道集上的地震数据在O点按不同倾角范围进行偏移叠加,获得不同角度的叠加能量;
(3)比较系列角度对应的偏移叠加能量大小,找到代表最大偏移能量的角度,确认该角度为界面面元法线角F;
(4)通过界面面元的法线角F确定最佳偏移叠加图像的数据之共中心点CMP;
(5)应用依据步骤(4)确定的CMP道集和其周围一定范围内道集的地震数据在O点进行叠前偏移成像。
常规Kirchhoff叠前时间偏移不应用射线追踪方法获得走时信息,因此,利用射线追踪和扫描估计得到的F角就近于奢侈。一般情况,若地下介质近于光滑和各向同性,作为近似,F角的扫描估计可以直接利用CMP到O点的连线来近似实现,亦即
图3是利用二维地震偏移数据对法线角进行扫描的实际范例。图左侧是偏移叠加剖面,在2700-2800毫秒处,地层倾斜角度大约为35度,通过将地面上以CMP为单位的地震数据在地下成像点进行偏移叠加,形成法线角道集(图3中的右侧图,横坐标为角度,纵坐标为时间)后,可以见到,右侧的偏移叠加能量集中在法线角为25-45度间,正好对应地层的倾斜度。对图中右侧数据进行能量扫描,最大能量对应35度左右。此说明地层界面的法线角完全可以通过对角度叠加道集数据能量的扫描过程获取。
如果找到界面面元的法线角和确定了可形成最佳偏移叠加图像的数据之共中心点CMP,选择性的挑选数据进行叠前偏移成像,成像剖面的分辨率将大大提高,从而实现聚焦叠加成像。在CMP道集上,炮点和检波点对CMP点是平面几何对称的。在各向同性均匀介质假设前提下,地下成像点的反射波只有在界面法向所指的CMP位置上才有记录,其它CMP点由于接受不到该点的反射,没有反射记录。然而,按照Fresnel波动理论,该法向角所指的数据共中心点CMP位置周围还应该具备接受反射波的条件,与CMP点不同,它们接受的是能量较弱的可称为衍射的地震波。因此,应用Fresnell第一反射带所包围的范围(图4),其中的阴影区是Fresnell带范围,W为叠加权。将地震数据进行加权偏移叠加,应用高斯权系数,可以获得地下界面面元最佳偏移—聚焦成像成果。
偏移叠加的范围可以应用Fresnell带界定,也可以应用其它方式界定。如果地下介质反射面是水平的,该范围的界定等效于偏移孔径大小的界定。
为了处理问题的方便,实际偏移过程中,当扫描得到法线F角后,可以直接应用F角的前后容忍度△F来界定叠加范围,即
Fmin=F-ΔF ……(8)
与
Fmax=F-ΔF ……(9)
其中:
Fmin是最小参与叠加的法线角;
Fmax是最大可参与叠加的法线角。
对地震数据进行聚焦叠加,意味着:
其中:
Image(0)为地下O点的像;SF为对应某CMP数据在该点的叠加像函数,该CMP与O点的法向角为F;WF为叠加权系数。
将上述地震数据进行加权偏移叠加,可以获得地下界面面元最佳偏移-聚焦成像成果。
地震数据若为三维数据,F角的选择还要在不同方位上进行,亦即要在不同方位角范围内选择F角,再进行上述偏移叠加成像。
附图说明
图1是地震CMP道集数据与地下反射面关系图。
图2是界面面元法线和激发点与接收点几何关系图。
图3是叠前偏移剖面与法线角道集数据关系图示。
图4是反射数据在成像点的偏移叠加示意图。
图5是应用常规技术(左)和应用本发明的技术(右)对地震数据进行成像结果对比图。
图6是本发明技术应用于叠前时间偏移获得的成像剖面与常规叠前时间偏移获得的成像剖面之对比图。
图7是国内某研究院应用常规的叠前时间偏移成像方法得到的地震反射成像剖面图。
图8应用本技术开发的软件对图7地震数据进行的叠前时间偏移成像处理对比图。
实施例1
应用本发明方法,我们对某地区地震数据进行了成像处理,获得了良好的成像效果。图5左侧显示了应用常规叠前时间偏移技术得到的最终偏移成像剖面成果,右侧显示了应用本发明获得的叠前时间偏移成像成果剖面。可以看出,应用本发明获得的地震图像具有较高的分辨率,偏移噪音明显较少,频率明显提高,特别是大套地层间的小层清晰可辨。
实施例2
下图6地震数据来源于某火山岩地区,信噪比中等,但是偏移叠加后同相轴连续性不理想。分析原因为部分干扰波能量较强,在偏移过程中由于覆盖次数少,所以噪音的叠加衰减不理想;为此,应用本发明技术,获得成像点界面面元法线角后,对法线角对应的地震数据共中心点CMP和周围若干点地震数据进行选择性加权叠加成像,获得了良好的更高分辨率的叠前偏移成像剖面。图6左侧图像应用常规叠加偏移方法,右侧图像应用本发明的技术进行的叠前时间偏移方法。比较两侧图像,新的技术确实提高了偏移剖面的分辨率,使得图像更加清晰,同相轴连续,图像整体噪音水平也降低了许多。
实施例3
某地震数据处理项目,要求对地震数据进行精细成像处理,断层和地层分辨清晰。此前,国内某研究院曾经对该数据进行过精细处理,没有达到甲方的要求。我们在认真分析了数据特点以后,认为新开发的成像技术有能力达到甲方的要求,所以应用本专利技术开发的成像软件对该数据进行了成像处理。
图7是国内某研究单位对该地震数据成像结果。图中可见主要断层和地层都比较清晰,但是较为细节的结构和主构造周围的地层、特别是浅部构造有些模糊,解释中似乎可以认为存在一条平卧的主断层并已经直接将地层错断。图8是应用我们开发的技术软件对该数据进行的成像处理结果。可见主要的结构与图7基本相同,但细结构有了很大的变化,控制本区上部地层发育的断层并没有成为将地层错断的大型断层从而促使上部发育为一个大型滑脱体,它实际上仍旧可以解释为一个大型平卧褶皱的核部。此处理成果得到甲方的充分肯定,为我们争得了荣誉,也为本技术的应用提供了强有力的市场支持。
Claims (2)
1.基于搜寻反射界面法线角获得最佳成像效果的叠前偏移成像方法,其特征在于包含以下步骤:
(1)首先从O点向一系列CMP点做射线追踪,在O点得到一系列射线与垂直线的交角,即倾角;
(2)再将该系列CMP道集上的地震数据在O点按不同倾角范围进行偏移叠加,获得不同角度的叠加能量;
(3)比较系列角度对应的偏移叠加能量大小,找到代表最大偏移能量的角度,确认该角度为界面面元法线角F,在地下介质近于光滑和各向同性假设情况下,法线角F直接利用CMP到O点的连线来近似实现,亦即
(4)通过界面面元的法线角F确定最佳偏移叠加图像的数据之共中心点CMP;
(5)应用依据步骤(4)确定的CMP道集和其周围一定范围内道集的地震数据在O点进行叠前偏移成像;
将在步骤(4)和(5)中Fresnell第一反射带范围内所得数据应用高斯系数进行加权偏移叠加成像;将在步骤(4)和(5)中法线角F前后一定范围△F内所得数据进行加权偏移叠加成像。
2.根据权利要求1所述基于搜寻反射界面法线角获得最佳成像效果的叠前偏移成像方法,其特征在于步骤(4)和(5)中的地震数据若为三维数据,F角的选择还要在不同方位上进行,亦即要在不同方位角范围内选择F角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210194803.6A CN102749644B (zh) | 2012-06-13 | 2012-06-13 | 搜寻反射界面法线角获得最佳成像效果叠前偏移成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210194803.6A CN102749644B (zh) | 2012-06-13 | 2012-06-13 | 搜寻反射界面法线角获得最佳成像效果叠前偏移成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102749644A CN102749644A (zh) | 2012-10-24 |
CN102749644B true CN102749644B (zh) | 2015-01-14 |
Family
ID=47029987
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210194803.6A Expired - Fee Related CN102749644B (zh) | 2012-06-13 | 2012-06-13 | 搜寻反射界面法线角获得最佳成像效果叠前偏移成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102749644B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103823240B (zh) * | 2014-03-11 | 2017-01-04 | 西南石油大学 | 一种基于crp的弯线采集方法 |
CN108693554A (zh) * | 2017-04-11 | 2018-10-23 | 中国石油化工股份有限公司 | 倾角域共成像点道集求取方法及计算机可读存储介质 |
CN110646840B (zh) * | 2018-06-27 | 2021-07-20 | 中国石油化工股份有限公司 | 角道集提取方法及*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004090575A1 (en) * | 2003-04-09 | 2004-10-21 | Norsar | Method for simulating local prestack depth migrated seismic images |
WO2005026776A1 (en) * | 2003-09-16 | 2005-03-24 | Geosystem S.R.L. | Wide-offset-range pre-stack depth migration method for seismic exploration |
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
CN102103216A (zh) * | 2009-12-16 | 2011-06-22 | 中国石油天然气集团公司 | 二维高斯射线束叠前偏移方法 |
CN102176053A (zh) * | 2011-01-27 | 2011-09-07 | 中国科学院地质与地球物理研究所 | 提升波动方程叠前深度偏移成像效果的方法 |
CN102183786A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 双复杂条件下保真振幅高斯束叠前深度偏移方法 |
CN102313903A (zh) * | 2011-07-01 | 2012-01-11 | 中国海洋石油总公司 | 基于波动方程外推算子的vti介质中叠前时间偏移方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6904368B2 (en) * | 2002-11-12 | 2005-06-07 | Landmark Graphics Corporation | Seismic analysis using post-imaging seismic anisotropy corrections |
-
2012
- 2012-06-13 CN CN201210194803.6A patent/CN102749644B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004090575A1 (en) * | 2003-04-09 | 2004-10-21 | Norsar | Method for simulating local prestack depth migrated seismic images |
WO2005026776A1 (en) * | 2003-09-16 | 2005-03-24 | Geosystem S.R.L. | Wide-offset-range pre-stack depth migration method for seismic exploration |
CN101630014A (zh) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | 一种利用垂直地震剖面数据对各向异性介质成像的方法 |
CN102103216A (zh) * | 2009-12-16 | 2011-06-22 | 中国石油天然气集团公司 | 二维高斯射线束叠前偏移方法 |
CN102176053A (zh) * | 2011-01-27 | 2011-09-07 | 中国科学院地质与地球物理研究所 | 提升波动方程叠前深度偏移成像效果的方法 |
CN102183786A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 双复杂条件下保真振幅高斯束叠前深度偏移方法 |
CN102313903A (zh) * | 2011-07-01 | 2012-01-11 | 中国海洋石油总公司 | 基于波动方程外推算子的vti介质中叠前时间偏移方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102749644A (zh) | 2012-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xie et al. | Wave-equation-based seismic illumination analysis | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
Nishizawa et al. | Seismic structure of rifting in the Okinawa Trough, an active backarc basin of the Ryukyu (Nansei-Shoto) island arc–trench system | |
US4611311A (en) | Method of seismic processing involving converted P- or S-wave data | |
Nowack et al. | Application of Gaussian-beam migration to multiscale imaging of the lithosphere beneath the Hi-CLIMB array in Tibet | |
CN103091710A (zh) | 一种逆时偏移成像方法及装置 | |
CN105093292A (zh) | 一种地震成像的数据处理方法和装置 | |
CN102854526B (zh) | 一种多分量地震资料的处理方法 | |
CN102749644B (zh) | 搜寻反射界面法线角获得最佳成像效果叠前偏移成像方法 | |
WO2004017097A1 (en) | Seismic survey design technique | |
Xu et al. | Seismic evidence for lateral asthenospheric flow beneath the northeastern Tibetan Plateau derived from S receiver functions | |
CN106950595B (zh) | 一种角变换成像方法与装置 | |
Babuška et al. | Mantle fabric of western Bohemian Massif (central Europe) constrained by 3D seismic P and S anisotropy | |
Saito | Anisotropic traveltime tomography at the Buckhorn Test Facility in Illinois | |
Tsutsui et al. | Seismic reflectors beneath the central cones of Aso Volcano, Kyushu, Japan | |
Stroujkova et al. | A magma mass beneath Casa Diablo? Further evidence from reflected seismic waves | |
CN107728196A (zh) | 获取角度域共成像道集的方法及*** | |
CN109655879B (zh) | 目的层信噪比一致性能量照明的观测***优化方法及装置 | |
AU721765B2 (en) | Method of generating a fold distribution and of evaluating seismic survey | |
CN104076395A (zh) | 基于滤波组合的镜面能量提取与成像方法 | |
Shen et al. | Research progress on layered seismic anisotropy-A review | |
Yang et al. | Separation of split shear waves based on a hodogram analysis of HTI media | |
Nakata et al. | Seismic imaging of a carbonate rock block with fracture-induced anisotropy | |
Tsingas et al. | Fracture detection by diffraction imaging | |
Pahlavanloo et al. | Improving seismic image in complex structures by new solving strategies in the CO-CRS and the CO-CDS methods |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150114 |