CN109387868A - 一种基于地震波同相轴斜率信息的三维层析成像方法 - Google Patents
一种基于地震波同相轴斜率信息的三维层析成像方法 Download PDFInfo
- Publication number
- CN109387868A CN109387868A CN201811136635.9A CN201811136635A CN109387868A CN 109387868 A CN109387868 A CN 109387868A CN 201811136635 A CN201811136635 A CN 201811136635A CN 109387868 A CN109387868 A CN 109387868A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- slope
- shot
- obtains
- geophone station
- 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.)
- Pending
Links
- 238000004587 chromatography analysis Methods 0.000 title claims abstract description 19
- 238000003384 imaging method Methods 0.000 title claims abstract description 17
- 239000012504 chromatography matrix Substances 0.000 claims abstract description 14
- 230000005012 migration Effects 0.000 claims abstract description 5
- 238000013508 migration Methods 0.000 claims abstract description 5
- 238000001228 spectrum Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000000034 method Methods 0.000 claims description 10
- 238000007689 inspection Methods 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims description 2
- 238000003325 tomography Methods 0.000 description 8
- 238000011208 chromatographic data Methods 0.000 description 6
- 238000000605 extraction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 241001282153 Scopelogadus mizolepis Species 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
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/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/514—Post-stack
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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)
Abstract
本发明涉及一种基于地震波同相轴斜率信息的三维层析成像方法,其特征在于包括以下步骤:1)对三维共炮点道集和共检波点道集进行有效数据拾取,得到炮点和检波点处的射线参数水平分量信息;2)根据得到的炮点和检波点处射线参数水平分量信息建立三维斜率层析矩阵;3)在建立的三维斜率层析矩阵线性方程组中加入正则化项,得到大规模稀疏线性方程组;4)对大规模稀疏线性方程组进行求解,获得斜率层析反演的宏观速度模型,用于三维偏移成像或作为三维全波形反演的初始模型。本发明可以广泛应用于勘探地震学领域。
Description
技术领域
本发明涉及勘探地震学领域,特别是关于一种基于地震波同相轴斜率信息的三维层析成像方法。
背景技术
宏观速度模型估计是勘探地震学中最为重要的参数提取工作,面向全带宽速度建模的全波形反演方法已是近年来国际研究热点。然而,在勘探地震学领域,任何先进的算法都必须经受实际数据的严格检验,才能够最终判定其价值。由于实际地震数据中低频信息的缺失、噪音干扰、子波估计以及各种现有理论尚难以描述的复杂波现象的存在,导致波形反演的目标函数凸性较差,在推向实用中依然遭遇各种棘手问题。反观基于地震数据运动学信息的层析成像算法,由于运动学信息相对易于提取以及程函方程本身的稳健可靠,数十年来依然是工业界在低频背景速度估算和偏移速度建模中最为倚重的建模工具。
斜率层析成像是一种极具特色的层析成像方法。从数据空间的利用方面来说,斜率层析成像巧妙利用了地震波场中通常为人所忽视的斜率信息。除了旅行时信息,局部相干同相轴在炮道集与检波点道集内的射线参数水平分量、炮点坐标和检波点坐标都被纳入到数据空间;除了速度,它还将地下反射点坐标,炮点一侧与检波点一侧的散射角也纳入反演的模型空间。这个特色使得斜率层析成像成为射线类层析成像方法中唯一一种可以同时反演速度结构、反射点位置与反射层局部形态的方法。
然而,目前仅仅提出了二维斜率层析成像算法用于实际地震数据的宏观速度建模,未曾将该算法推广到三维,而三维斜率层析成像需要解决更多问题,比如三维地震数据中斜率信息的高效提取和质量检查问题、三维斜率层析矩阵的建立问题、以及更强有力的正则化项求取问题。
发明内容
针对上述问题,本发明的目的是提供一种基于地震波同相轴斜率信息的三维层析成像方法,该方法基于倾斜叠加能量谱进行斜率信息的拾取,然后利用拾取的斜率信息建立三维斜率层析矩阵,并引入正则化项对三为斜率层析矩阵进行求解。
为实现上述目的,本发明采取以下技术方案:一种基于地震波同相轴斜率信息的三维层析成像方法,其包括以下步骤:1)对三维共炮点道集和共检波点道集进行有效数据拾取,得到炮点和检波点处的射线参数水平分量信息;2)根据得到的炮点和检波点处射线参数水平分量信息建立三维斜率层析矩阵;3)在建立的三维斜率层析矩阵线性方程组中加入正则化项,得到大规模稀疏线性方程组;4)对大规模稀疏线性方程组进行求解,获得斜率层析反演的宏观速度模型,用于三维偏移成像或作为三维全波形反演的初始模型。
进一步地,所述步骤1)中,对三维共炮点道集和共检波点道集进行有效数据拾取,得到炮点和检波点处的射线参数水平分量信息的方法,包括以下步骤:1.1)对于三维共炮点道集,计算其倾斜叠加谱,通过人工或者半自动方式拾取倾斜叠加谱中能量大于炮点能量阈值的炮点位置,并获得相应的炮点处的射线参数水平分量prx,pry;1.2)对于三维共检波点道集,计算其倾斜叠加谱,通过人工或者半自动方式拾取倾斜叠加谱中能量大于检波点能量阈值的检波点位置,并获得相应的检波点处的射线参数水平分量psx,psy。
进一步地,所述步骤3)中,建立的三维斜率层析矩阵为:
式中,psx,psy为炮点坐标;prx,pry为检波点坐标,x0,y0,z0为地下反射点坐标,θs,φs代表从炮点一侧出射的地下张角与方位角,θr,φr代表从检波点一侧出射的地下张角与方位角,v代表地下介质速度信息。
进一步地,所述步骤4)中,所述大规模稀疏线性方程组为:
式中,m0为先验模型,α0为常规阻尼正则化因子,H1,H2,H3分别为X,Y,Z方向上的二阶光滑因子,△d为数据空间的残差,△m为模型空间的残差。
本发明由于采取以上技术方案,其具有以下优点:1、本发明基于倾斜叠加能量谱的拾取引入了人工干预,可以有效避开相干性的干扰波和多次波,将斜率信息的拾取聚焦于一次反射波,拾取的数据空间可靠性更高;2、本发明利用拾取的斜率信息建立三维斜率层析矩阵,相比于二维斜率层析矩阵,不仅可以反演宏观速度模型,同时也可以反演反射点位置与构造倾角,这是传统的旅行时反演做不到的;3、本发明在求解层析线性方程组的过程中引入了更强有力的正则化项进行约束,改善了现有方程组的条件数,保证了线性方程组的稳定求解,进一步提高了求解的准确性。因此,本发明可以广泛应用于勘探地震学领域中。
附图说明
图1是二维斜率层析数据空间各分量与模型空间各分量示意图,其中,图1(a)是斜率层析的模型与数据分量,图1(b)是同相轴斜率与射线理论的慢度矢量。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1(a)和图1(b)所示,本发明首先对二维斜率层析成像的实现过程进行描述。图1(a)显示了二维斜率层析成像的模型与数据分量示意图,根据实际地震观测数据分别建立层析数据空间和地下模型空间。其中图1(a)中观测数据是从实际地震数据中提取出来的,正演数据则是在反演过程中基于射线正演获得的数据空间,层析反演的本质就是利用数据残差求取模型残差,当模型残差作用到初始模型上,初始模型就得到了更新,当前轮反演工作完成;在更新后的模型上实施射线正演就可以获得新的数据空间dcal,于是就进入下一轮迭代。注意斜率层析中所需要的层析数据空间不包含旅行时和射线地表出射坐标,为d=[psx,prx],其中psx和prx分别为炮点和检波点处的射线参数水平分量;建立的地下的模型空间为m=(x0,z0,θs,θr,v),其中x0,z0为地下反射点X的坐标;θs和θr分别表示从炮点一侧与检波点一侧出射的地下张角,v代表地下介质速度信息。
注意数据空间中的Psx可以在共检波点道集内利用倾斜叠加在倾斜叠加能量谱内拾取得到,Prx可以在共炮点道集内利用倾斜叠加在倾斜叠加能量谱内拾取得到,原理如下:图1(b)显示了一个共炮道集,如果考察某一局部相干同相轴,由几何关系以及射线理论中慢度矢量的定义得:
式(1)中,pslope为局部相干同相轴的斜率,prx为检波点处射线参数水平分量,△t为叠前炮道集在时间方向的采样率,△x为叠前炮道集在空间方向的采样率。式(1)表明共炮点道集内同相轴的斜率对应于检波点处射线参数水平分量prx,因此,由炮检互易原理易知,共检波点道集内同相轴的斜率必对应于炮点处射线参数水平分量psx。需要指出,正是由于共炮道集和共检波点道集数据内同相轴的斜率(即prx,psx)被引入层析数据空间,使得在层析反演中反射波在炮检点处的局部传播方向得到强有力的约束。具体的优势体现在斜率层析反演中无需再去拾取连续的反射同相轴,只要选择连续性好、信噪比高的局部波包即可。同时,将反射问题化为透射问题,也回避了传统反射层析中速度与反射层深度耦合的问题。
但是,由于在层析数据空间中引入了prx,psx,模型空间也随之扩大,使得斜率层析在反演速度的同时也可以同时反演反射点的深度位置以及反射层的倾角。公式(2)展示了二维斜率层析反演的层析矩阵(或FRECHET偏导数矩阵),其中σ为针对层析矩阵每一行物理量的数量级不同设置的均衡系数。可以看到常规层析反演由于其模型空间一般仅由速度V构成,其数据空间一般仅考虑旅行时T,可认为其层析矩阵仅仅相当于公式(2)中的某一个子矩阵的规模。因此斜率层析矩阵无论从规模还是其稀疏性方面都超过了常规层析反演矩阵,因此高质量的数据空间提取对降低其稀疏性、降低其对规则化的要求、进而提高求解的精度是非常有意义的。
选择斜率层析数据空间的标准就是优先考虑连续性好、信噪比高的一次反射局部波包。从数据空间的特点我们不难看出,只要prx和psx参数选择好了,其他参数如炮检点横坐标和旅行时都可以轻易得到。而从公式(1)可以知道,prx和psx参数估计等价于同相轴局部斜率。因此如何准确有效地估计共炮点道集与共检波点道集内同相轴的斜率就是斜率层析反演中最为重要的数据准备工作。
基于上述分析,本发明提供的一种基于地震波同相轴斜率信息的三维层析成像方法,包括以下步骤:
1)对三维共炮点道集和共检波点道集进行有效数据拾取,得到炮点和检波点处的射线参数水平分量信息;具体的包括以下步骤:
1.1)对于三维共炮点道集,计算其倾斜叠加谱,通过人工或者半自动方式拾取倾斜叠加谱中能量大于炮点能量阈值的炮点位置,并获得相应的炮点处的射线参数水平分量prx,pry;其中,炮点能量阈值根据实际情况设定。
1.2)将数据分选为三维共检波点道集,计算其倾斜叠加谱,通过人工或者半自动方式拾取倾斜叠加谱中能量大于检波点能量阈值的检波点位置,并获得相应的检波点处的射线参数水平分量psx,psy;其中,检波点能量阈值根据实际情况设定。
2)根据得到的炮点和检波点处射线参数水平分量信息建立三维斜率层析矩阵。
三维情形下,射线参数具有两个水平分量,因此炮点与检波点处的射线参数水平分量都变成了两个,于是数据空间变为d=(psx,psy,prx,pry),模型空间变为m=(x0,y0,z0,θs,φs,θr,φr,v),其中,x0,y0,z0为地下反射点坐标,θs,φs代表从炮点一侧出射的地下张角与方位角,θr,φr代表从检波点一侧出射的地下张角与方位角,v代表地下介质速度信息。本发明利用降阶汉米尔顿算子在直角坐标系下推导层析矩阵,将公式(2)所示的二维斜率层析矩阵转换为如公式(3)所示的三维斜率层析矩阵:
3)在建立的三维斜率层析矩阵线性方程组中加入更为强有力的正则化项,得到大规模稀疏线性方程组;
三维斜率层析矩阵的基本构建思路与二维并无区别,都是将从炮点S→反射点X→检波点R的反射过程分解为从X→S以及X→R这样两个透射过程。其中数据空间的残差△d与模型空间残差△m之间的关系由射线扰动理论建立(FARRA,1987),这里不妨简写为:
△d=F△m (4)
在三维情形下,由于三维斜率层析矩阵F更为庞大稀疏,仅仅依靠阻尼正则化因子已经无法获得稳定的解,必须附加更为强有力的光滑正则化项,才能够获得稳定的解,此时得到的大规模稀疏线性方程组为:
式中,m0为先验模型,α0为常规阻尼正则化因子,H1,H2,H3分别为X,Y,Z方向上的二阶光滑因子。
4)对大规模稀疏线性方程组进行求解,获得斜率层析反演的宏观速度模型,用于三维偏移成像或作为三维全波形反演的初始模型。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和制作工艺等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (4)
1.一种基于地震波同相轴斜率信息的三维层析成像方法,其特征在于包括以下步骤:
1)对三维共炮点道集和共检波点道集进行有效数据拾取,得到炮点和检波点处的射线参数水平分量信息;
2)根据得到的炮点和检波点处射线参数水平分量信息建立三维斜率层析矩阵;
3)在建立的三维斜率层析矩阵线性方程组中加入正则化项,得到大规模稀疏线性方程组;
4)对大规模稀疏线性方程组进行求解,获得斜率层析反演的宏观速度模型,用于三维偏移成像或作为三维全波形反演的初始模型。
2.如权利要求1所述的一种基于地震波同相轴斜率信息的三维层析成像方法,其特征在于:所述步骤1)中,对三维共炮点道集和共检波点道集进行有效数据拾取,得到炮点和检波点处的射线参数水平分量信息的方法,包括以下步骤:
1.1)对于三维共炮点道集,计算其倾斜叠加谱,通过人工或者半自动方式拾取倾斜叠加谱中能量大于炮点能量阈值的炮点位置,并获得相应的炮点处的射线参数水平分量prx,pry;
1.2)对于三维共检波点道集,计算其倾斜叠加谱,通过人工或者半自动方式拾取倾斜叠加谱中能量大于检波点能量阈值的检波点位置,并获得相应的检波点处的射线参数水平分量psx,psy。
3.如权利要求1所述的一种基于地震波同相轴斜率信息的三维层析成像方法,其特征在于:所述步骤3)中,建立的三维斜率层析矩阵为:
式中,psx,psy为炮点坐标;prx,pry为检波点坐标,x0,y0,z0为地下反射点坐标,θs,φs代表从炮点一侧出射的地下张角与方位角,θr,φr代表从检波点一侧出射的地下张角与方位角,v代表地下介质速度信息。
4.如权利要求1所述的一种基于地震波同相轴斜率信息的三维层析成像方法,其特征在于:所述步骤4)中,所述大规模稀疏线性方程组为:
式中,m0为先验模型,α0为常规阻尼正则化因子,H1,H2,H3分别为X,Y,Z方向上的二阶光滑因子,△d为数据空间的残差,△m为模型空间的残差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811136635.9A CN109387868A (zh) | 2018-09-28 | 2018-09-28 | 一种基于地震波同相轴斜率信息的三维层析成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811136635.9A CN109387868A (zh) | 2018-09-28 | 2018-09-28 | 一种基于地震波同相轴斜率信息的三维层析成像方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109387868A true CN109387868A (zh) | 2019-02-26 |
Family
ID=65418170
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811136635.9A Pending CN109387868A (zh) | 2018-09-28 | 2018-09-28 | 一种基于地震波同相轴斜率信息的三维层析成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109387868A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110780348A (zh) * | 2019-11-01 | 2020-02-11 | 中国石油大学(华东) | 基于立体成像条件的一次波和多次波联合成像方法及*** |
CN111781639A (zh) * | 2020-06-04 | 2020-10-16 | 同济大学 | 针对obs多分量数据的炮检互易弹性波全波形反演方法 |
CN112782758A (zh) * | 2019-11-11 | 2021-05-11 | 中国石油天然气集团有限公司 | 稀疏采样观测***的状态确定方法和装置 |
CN114839675A (zh) * | 2021-01-31 | 2022-08-02 | 中国石油化工股份有限公司 | 一种建立三维速度模型的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020183980A1 (en) * | 2001-03-05 | 2002-12-05 | Patrice Guillaume | Methods of tomographically inverting picked events on migrated seismic data |
CN104714251A (zh) * | 2015-03-24 | 2015-06-17 | 西安交通大学 | 用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法 |
CN105319589A (zh) * | 2014-07-25 | 2016-02-10 | 中国石油化工股份有限公司 | 一种利用局部同相轴斜率的全自动立体层析反演方法 |
CN107505651A (zh) * | 2017-06-26 | 2017-12-22 | 中国海洋大学 | 地震初至波和反射波联合斜率层析成像方法 |
CN107515420A (zh) * | 2017-04-28 | 2017-12-26 | 西安石油大学 | 一种用于局部相关同相轴的走时与梯度精确拾取方法 |
-
2018
- 2018-09-28 CN CN201811136635.9A patent/CN109387868A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020183980A1 (en) * | 2001-03-05 | 2002-12-05 | Patrice Guillaume | Methods of tomographically inverting picked events on migrated seismic data |
CN105319589A (zh) * | 2014-07-25 | 2016-02-10 | 中国石油化工股份有限公司 | 一种利用局部同相轴斜率的全自动立体层析反演方法 |
CN104714251A (zh) * | 2015-03-24 | 2015-06-17 | 西安交通大学 | 用于同相轴自动拾取的倾斜叠加峰值振幅处边缘检测法 |
CN107515420A (zh) * | 2017-04-28 | 2017-12-26 | 西安石油大学 | 一种用于局部相关同相轴的走时与梯度精确拾取方法 |
CN107505651A (zh) * | 2017-06-26 | 2017-12-22 | 中国海洋大学 | 地震初至波和反射波联合斜率层析成像方法 |
Non-Patent Citations (2)
Title |
---|
倪瑶 等: "立体层析反演方法理论分析与应用测试", 《石油物探》 * |
杨锴 等: "基于非降阶汉密尔顿算子的三维立体层析反演", 《地球物理学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110780348A (zh) * | 2019-11-01 | 2020-02-11 | 中国石油大学(华东) | 基于立体成像条件的一次波和多次波联合成像方法及*** |
CN110780348B (zh) * | 2019-11-01 | 2021-07-09 | 中国石油大学(华东) | 基于立体成像条件的一次波和多次波联合成像方法及*** |
CN112782758A (zh) * | 2019-11-11 | 2021-05-11 | 中国石油天然气集团有限公司 | 稀疏采样观测***的状态确定方法和装置 |
CN111781639A (zh) * | 2020-06-04 | 2020-10-16 | 同济大学 | 针对obs多分量数据的炮检互易弹性波全波形反演方法 |
CN114839675A (zh) * | 2021-01-31 | 2022-08-02 | 中国石油化工股份有限公司 | 一种建立三维速度模型的方法 |
CN114839675B (zh) * | 2021-01-31 | 2023-09-05 | 中国石油化工股份有限公司 | 一种建立三维速度模型的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108037526B (zh) | 基于全波波场vsp/rvsp地震资料的逆时偏移方法 | |
CN109387868A (zh) | 一种基于地震波同相轴斜率信息的三维层析成像方法 | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN107505651B (zh) | 地震初至波和反射波联合斜率层析成像方法 | |
US6388947B1 (en) | Multi-crosswell profile 3D imaging and method | |
US10317548B2 (en) | Reflection seismic data Q tomography | |
CN106483559B (zh) | 一种地下速度模型的构建方法 | |
US9869783B2 (en) | Structure tensor constrained tomographic velocity analysis | |
Yang et al. | Illumination compensation for image-domain wavefield tomography | |
CN105093281B (zh) | 一种反演框架下的地震多波建模方法 | |
CN107219554B (zh) | 陆地地震资料的剩余静校正量的自动获取方法 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN105319589B (zh) | 一种利用局部同相轴斜率的全自动立体层析反演方法 | |
CN109116413A (zh) | 成像域立体层析速度反演方法 | |
CN107656308B (zh) | 一种基于时间深度扫描的共散射点叠前时间偏移成像方法 | |
CN113156500B (zh) | 数据驱动的快速构造约束叠前地震多道反演方法 | |
Zhu et al. | Recent applications of turning-ray tomography | |
Guo et al. | Becoming effective velocity-model builders and depth imagers, Part 2—The basics of velocity-model building, examples and discussions | |
CN106338760B (zh) | 误差补偿的起伏地表偏移方法 | |
Alerini et al. | Prestack depth imaging of ocean-bottom node data | |
Wu et al. | Offset VSP P wave and converted S wave salt proximity migration study | |
Dong et al. | Automatic migration velocity estimation for prestack time migration | |
Xing et al. | Application of 3D stereotomography to the deep-sea data acquired in the South China Sea: a tomography inversion case | |
Ni* et al. | Preliminary practice of stereotomography | |
CN111538081A (zh) | 一种地震数据初至波的外推方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190226 |
|
RJ01 | Rejection of invention patent application after publication |