CN106950595A - 一种角变换成像方法与装置 - Google Patents

一种角变换成像方法与装置 Download PDF

Info

Publication number
CN106950595A
CN106950595A CN201710071542.1A CN201710071542A CN106950595A CN 106950595 A CN106950595 A CN 106950595A CN 201710071542 A CN201710071542 A CN 201710071542A CN 106950595 A CN106950595 A CN 106950595A
Authority
CN
China
Prior art keywords
angular transformation
angular
ray
angle
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
Application number
CN201710071542.1A
Other languages
English (en)
Other versions
CN106950595B (zh
Inventor
周青春
谢桂生
康德拉什科夫
姜瑞友
陈相君
田迎春
张云波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
New JCM Group Co.,Ltd.
Original Assignee
Landocean Group Ltd By Share Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Landocean Group Ltd By Share Ltd filed Critical Landocean Group Ltd By Share Ltd
Priority to CN201710071542.1A priority Critical patent/CN106950595B/zh
Publication of CN106950595A publication Critical patent/CN106950595A/zh
Application granted granted Critical
Publication of CN106950595B publication Critical patent/CN106950595B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/34Displaying seismic recordings or visualisation of seismic data or attributes
    • G01V1/345Visualisation of seismic data or attributes, e.g. in 3D cubes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/70Other details related to processing
    • G01V2210/74Visualisation of seismic data

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)

Abstract

本发明提出了一种角变换成像方法与装置,具体为通过扫描下行波与上行波的夹角、地层倾角两个参数,实现地震叠加成像,消除实际地震波射线速度与椭圆展开方法估算的法向成像速度不同导致的聚焦成像不理想问题,改善成像效果并提高成像精度,并为后续处理提供更好、更多的基础性资料。

Description

一种角变换成像方法与装置
技术领域
本发明涉及地震勘探中反射地震资料处理技术领域,特别涉及一种角变换成像方法与装置。
背景技术
随着地震勘探的不断发展,实际地震资料处理生产中,主要采用常规共中心点(CMP)道集动校正(NMO)后叠加的方法。它假定地下为水平层状均匀介质,采用CMP抽道集、双曲线动校正、多次覆盖同相叠加的方式实现叠加成像。这在地下构造较为简单,地层倾角较小的情况下,可以得到比较好的叠加剖面。但当地层倾斜或弯曲时,反射点发散,该方法难以实现反射波共反射点同相叠加,使得叠加剖面难以取得较真实的结果。而且所求得的速度为反射界面曲率的函数,没有实际地质意义。为了弥补常规CMP叠加方法的缺陷,国内外的地球物理学者做了大量的科学研究和实践工作,形成了一些较好的成像技术,如倾角校正(DMO)叠加成像技术、共反射面元叠加(CRS)技术、椭圆展开成像技术等,取得了一定成果。DMO叠加成像技术是对NMO叠加公式引入地层倾角校正项,将非零炮检距地震记录转化为自激自收地震记录,保证在任何倾斜地层的情况下都能实现共反射点叠加的一种重要技术。但是,DMO技术有诸多假设条件限制,如常速度、倾角或炮检距的限制等。
20世纪90年代,Hubral等人提出了一种不依赖于宏观速度模型的反射成像方法——共反射面元叠加(CRS)技术。该方法可大幅度地提高地震资料信噪比和分辨率,被视为今后深层和复杂地区地震资料处理方法的重要发展方向。Hocht与Bergler推导了适合三维情形的CRS叠加表达式,该表达式需要八个参数确定走时关系,计算量大,效率低,不实用。此外,该方法存在大偏移距成像精度较差的问题和数据驱动策略导致CRS叠加剖面上出现“倾角歧视现象”的问题,导致动力学特征失真并且影响叠后偏移的质量。椭圆展开成像方法的理论基础是射线理论。均匀介质条件下,利用该方法可以对任意形状的反射界面求取唯一正确的共反射点叠加速度,得到真正的零偏移距时间剖面。从而比较好地解决了DMO和CRS技术等无法彻底解决的准确求取叠加速度和反射点发散等问题。非均匀介质条件下,上行波速度和下行波速度存在差异,基于均匀介质导出的椭圆展开成像方法不再适应,需做改进。
发明内容
本发明的目的旨在至少解决所述的技术缺陷之一。
为此,本发明的一个目的在于提出一种角变换成像方法与装置,通过扫描下行波与上行波的夹角、地层倾角两个参数,消除实际地震波射线速度与椭圆展开方法估算的法向成像速度不同导致的聚焦成像不理想问题,改善成像效果并提高成像精度,并为后续处理提供更好、更多的基础性资料。
为了实现上述目的,本发明一方面的实施例提供一种角变换成像方法,其特征在于,包括以下步骤:
步骤S1,利用炮点、检波点和反射点画圆,将所述圆反射点轨迹拟合曲线;
步骤S2,在所述圆内添加辅助线,所述辅助线包括每个反射点处的角平分线,经几何理论推导可知,该角平分线是倾角为的反射界面的法线。所有的法线相交于一点,该交点也是圆和纵向垂直的直径的交点,在弧上的其他任意一点圆都高于炮检线;相等的角α/2将圆弧分成相等的两段弧;α为入射线和反射线之间的夹角;
步骤S3,根据椭圆展开理论、椭圆展开变换方程和界面倾角与射线参数的关系建立角变换方程
其中,-下行射线与反射线之间夹角的一半,-地下反射界面的倾角,或者中间射线与垂向线的夹角;
步骤S4,根据角变换方程,调整参数,求取角变换叠加剖面。
优选的,步骤S2中所述的椭圆展开变换方程为
其中,l为炮点到检波点的距离,l0为反射点法线在地面出射点与炮点的距离;t0为沿法向的双程旅行时间。
优选的,步骤S2中所述的界面倾角与射线参数的关系式为:
其中,l为炮点到检波点的距离,l0为反射点法线在地面出射点与炮点的距离;t0为沿法向的双程旅行时间。
优选的,步骤S4中,获取角变换叠加剖面时,通过各个展开等时线相切干涉叠加的方法,获得零偏移距剖面。
本发明的另一个目的在于提供一种角变换成像装置,包括轨迹拟合模块、分析辅助模块、角变换模块、角变换叠加剖面模块;
所述轨迹拟合模块,利用炮点、检波点和反射点画圆,用于将所述圆反射点轨迹拟合曲线;
所述分析辅助模块,用于在所述圆内添加辅助线,所述辅助线包括每个反射点处的角平分线,经几何理论推导可知,该角平分线是倾角为的反射界面的法线。所有的法线相交于一点,该交点也是圆和纵向垂直的直径的交点,在弧上的其他任意一点圆都高于炮检线;相等的角α/2将圆弧分成相等的两段弧;α为入射线和反射线之间的夹角;
所述角变换模块,用于根据椭圆展开理论、椭圆展开变换方程和界面倾角与射线参数的关系建立角变换方程
其中-下行射线与反射线之间夹角的一半,-地下反射界面的倾角,或者中间射线与垂向线的夹角;
所述角变换叠加剖面模块,用于根据角变换方程,调整参数,求取角变换叠加剖面。
优选的,所述角变换模块所述的PRO变换方程为
其中,l为炮点到检波点的距离,l0为反射点法线在地面出射点与炮点的距离;t0为沿法向的双程旅行时间。
优选的,所述角变换模块中,所述的界面倾角与射线参数的关系式为:
其中,l为炮点到检波点的距离,l0为反射点法线在地面出射点与炮点的距离;t0为沿法向的双程旅行时间。
优选的,角变换叠加剖面模块中,获取角变换叠加剖面时,通过各个展开等时线相切干涉叠加的方法,获得零偏移距剖面。
根据本发明实施例的提供的角变换成像方法与装置,通过扫描下行波与上行波的夹角、地层倾角两个参数,实现地震叠加成像,消除实际地震波射线速度与椭圆展开方法估算的法向成像速度不同导致的聚焦成像不理想问题,改善成像效果并提高成像精度,并为后续处理提供更好、更多的基础性资料。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为本发明一种角变换成像方法的流程图;
图2为本发明一种角变换成像装置的连接结构框图;
图3为地震波射线传播示意图;
图4为本发明实施例中均匀介质时角变换射线图;
图5为本发明实施例中建立的参数模型示意图;
图6为本发明实施例中参数模型中地层倾角扫描结果示意图;
图7为本发明实施例参数模型中入射角α/2扫描结果示意图;
图8本发明实施例参数模型中地层倾角为45°时角变换叠加剖面。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
如图1所示,本发明实施例的提出一种角变换成像方法,包括以下步骤:
步骤S1,利用炮点、检波点和反射点画圆,将所述圆反射点轨迹拟合曲线。
步骤S2,在所述圆内添加辅助线,所述辅助线包括每个反射点处的角平分线,经几何理论推导可知,该角平分线是倾角为的反射界面的法线。所有的法线相交于一点,该交点也是圆和纵向垂直的直径的交点,在弧上的其他任意一点圆都高于炮检线;相等的角α/2将圆弧分成相等的两段弧;α为入射线和反射线之间的夹角。
步骤S3,根据椭圆展开理论、椭圆展开变换方程和界面倾角与射线参数的关系建立角变换方程
从椭圆展开理论已知,椭圆展开变换方程为
界面倾角与射线参数的关系式为:
其中,l为炮点到检波点的距离,l0为反射点法线在地面出射点与炮点的距离;t0为沿法向的双程旅行时间。
进而变换得出角变换方程
其中,-下行射线与反射线之间夹角的一半,-地下反射界面的倾角,或者中间射线与垂向线的夹角;
角变换方程,对于均匀介质来说,它绝对准确,通过该方程,我们可以实现角变换地震数据成像。
步骤S4,根据角变换方程,调整参数,求取角变换叠加剖面。在步骤S4中,获取角变换叠加剖面时,通过各个展开等时线相切干涉叠加的方法,获得零偏移距剖面。
如图5所示,模型示意图有六组炮检对,地下有一个反射地层,地层倾角45°,这六组炮检对在该反射地层上有同一个反射点。由于反射点相同,因此图中直线l0为六组炮检的公共法线;反射地层右下方的曲线为共反射点时距曲线。
图6-图8为上述模型的处理结果。图6为地层倾角扫描结果,可以看到各时距曲线相交时的地层倾角为45°与模型一致。需要说明的是地层倾角为本模型中选取的值;根据本发明提供的方法建立的本模型,可以根据实际情况设定任意地层倾角。图7为入射角α/2扫描结果,其中α/2=25°。图8为地层倾角为45°时角变换叠加剖面,各个展开等时线相切干涉叠加,成功获得零偏移距剖面。
本发明的另一个目的在于提供一种角变换成像装置,包括轨迹拟合模块1、分析辅助模块2、角变换模块3、角变换叠加剖面模块4;
所述轨迹拟合模块1,利用炮点、检波点和反射点画圆,用于将所述圆反射点轨迹拟合曲线;
所述分析辅助模块2,用于在所述圆内添加辅助线,所述辅助线包括每个反射点处的角平分线,经几何理论推导可知,该角平分线是倾角为的反射界面的法线。所有的法线相交于一点,该交点也是圆和纵向垂直的直径的交点,在弧上的其他任意一点圆都高于炮检线;相等的角α/2将圆弧分成相等的两段弧;α为入射线和反射线之间的夹角;
所述角变换模块3,用于根据椭圆展开理论、椭圆展开变换方程和界面倾角与射线参数的关系建立角变换方程;
从椭圆展开理论已知,椭圆展开变换方程为参照上述公式1;界面倾角与射线参数的关系式为参照上述公式2;角变换方程为参照上述公式3。
角变换叠加剖面模块4,根据角变换方程,调整参数,通过各个展开等时线相切干涉叠加的方法,求取角变换叠加剖面,即获得零偏移距剖面。
结合附图3-8可知,参考上述实施例中的模型示意图有六组炮检对,地下有一个反射地层,地层倾角45°,这六组炮检对在该反射地层上有同一个反射点。图6-图8为模型的处理结果。图6为地层倾角扫描结果,可以看到各时距曲线相交时的地层倾角为45°,与模型一致。图7为入射角α/2扫描结果,其中α/2=25°。图8为地层倾角为45°时角变换叠加剖面,各个展开等时线相切干涉叠加,成功获得零偏移距剖面。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。

Claims (8)

1.一种角变换成像方法,其特征在于,包括以下步骤
步骤S1,利用炮点、检波点和反射点画圆,将所述圆反射点轨迹拟合曲线;
步骤S2,在所述圆内添加辅助线,所述辅助线包括每个反射点处的角平分线,经几何理论推导可知,该角平分线是倾角为的反射界面的法线;所有的法线相交于一点,该交点也是圆和纵向垂直的直径的交点,在弧上的其他任意一点圆都高于炮检线;相等的角α/2将圆弧分成相等的两段弧;α为入射线和反射线之间的夹角;
步骤S3,根据椭圆展开理论、椭圆展开变换方程和界面倾角与射线参数的关系建立角变换方程
其中-下行射线与反射线之间夹角的一半,-地下反射界面的倾角,或者中间射线与垂向线的夹角;
步骤S4,根据角变换方程,调整参数,求取角变换叠加剖面。
2.根据权利要求1所述的角变换成像方法。其特征在于,步骤S2中所述的椭圆展开变换方程为
其中,l为炮点到检波点的距离,l0为反射点法线在地面出射点与炮点的距离;t0为沿法向的双程旅行时间。
3.根据权利要求1所述的角变换成像方法,其特征在于,步骤S2中所述的界面倾角与射线参数的关系式为:
其中,l为炮点到检波点的距离,l0为反射点法线在地面出射点与炮点的距离;t0为沿法向的双程旅行时间。
4.根据权利要求1所述的角变换成像方法,其特征在于,步骤S4中,获取角变换叠加剖面时,通过各个展开等时线相切干涉叠加的方法,获得零偏移距剖面。
5.一种角变换成像装置,其特征在于,包括轨迹拟合模块、分析辅助模块、角变换模块、角变换叠加剖面模块;
所述轨迹拟合模块,利用炮点、检波点和反射点画圆,用于将所述圆反射点轨迹拟合曲线;
所述分析辅助模块,用于在所述圆内添加辅助线,所述辅助线包括每个反射点处的角平分线,经几何理论推导可知,该角平分线是倾角为的反射界面的法线,所有的法线相交于一点,该交点也是圆和纵向垂直的直径的交点,在弧上的其他任意一点圆都高于炮检线;相等的角α/2将圆弧分成相等的两段弧;α为入射线和反射线之间的夹角;
所述角变换模块,用于根据椭圆展开理论、椭圆展开变换方程和界面倾角与射线参数的关系建立角变换方程
其中-下行射线与反射线之间夹角的一半,-地下反射界面的倾角,或者中间射线与垂向线的夹角;
所述角变换叠加剖面模块,用于根据角变换方程,调整参数,求取角变换叠加剖面。
6.根据权利要求5所述的角变换成像装置,其特征在于,所述角变换模块中所述的椭圆展开变换方程为
其中,l为炮点到检波点的距离,l0为反射点法线在地面出射点与炮点的距离;t0为沿法向的双程旅行时间。
7.根据权利要求5所述的角变换成像装置,其特征在于,所述角变换模块中,所述的界面倾角与射线参数的关系式为:
其中,l为炮点到检波点的距离,l0为反射点法线在地面出射点与炮点的距离;t0为沿法向的双程旅行时间。
8.根据权利要求5所述的角变换成像装置,其特征在于,角变换叠加剖面模块中,获取角变换叠加剖面时,通过各个展开等时线相切干涉叠加的方法,获得零偏移距剖面。
CN201710071542.1A 2017-02-09 2017-02-09 一种角变换成像方法与装置 Active CN106950595B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710071542.1A CN106950595B (zh) 2017-02-09 2017-02-09 一种角变换成像方法与装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710071542.1A CN106950595B (zh) 2017-02-09 2017-02-09 一种角变换成像方法与装置

Publications (2)

Publication Number Publication Date
CN106950595A true CN106950595A (zh) 2017-07-14
CN106950595B CN106950595B (zh) 2019-05-21

Family

ID=59466179

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710071542.1A Active CN106950595B (zh) 2017-02-09 2017-02-09 一种角变换成像方法与装置

Country Status (1)

Country Link
CN (1) CN106950595B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112346122A (zh) * 2020-11-06 2021-02-09 中国地震灾害防御中心 一种地震数据处理vda双参数成像方法
CN112346124A (zh) * 2020-11-06 2021-02-09 中国地震灾害防御中心 一种地震数据处理via双参数成像方法
CN112346125A (zh) * 2020-11-06 2021-02-09 中国地震灾害防御中心 一种地震数据处理vda双参数分析方法
CN112346123A (zh) * 2020-11-06 2021-02-09 中国地震灾害防御中心 一种地震数据处理via双参数分析方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040220743A1 (en) * 2003-04-30 2004-11-04 Conocophillips Company Method for determining shear-wave velocity model for depth migration of mode-converted data
CN101900832A (zh) * 2010-06-25 2010-12-01 恒泰艾普石油天然气技术服务股份有限公司 地震资料处理真地表条件下的椭圆展开成像方法和装置
CN101900831A (zh) * 2010-06-25 2010-12-01 恒泰艾普石油天然气技术服务股份有限公司 地震资料处理真地表条件下的椭圆展开成像方法和装置
CN104570069A (zh) * 2013-10-12 2015-04-29 中国石油集团东方地球物理勘探有限责任公司 一种从起伏地表直接成像的地震成像方法和装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040220743A1 (en) * 2003-04-30 2004-11-04 Conocophillips Company Method for determining shear-wave velocity model for depth migration of mode-converted data
CN101900832A (zh) * 2010-06-25 2010-12-01 恒泰艾普石油天然气技术服务股份有限公司 地震资料处理真地表条件下的椭圆展开成像方法和装置
CN101900831A (zh) * 2010-06-25 2010-12-01 恒泰艾普石油天然气技术服务股份有限公司 地震资料处理真地表条件下的椭圆展开成像方法和装置
CN104570069A (zh) * 2013-10-12 2015-04-29 中国石油集团东方地球物理勘探有限责任公司 一种从起伏地表直接成像的地震成像方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周青春 等: "椭圆展开共反射点叠加方法的应用研究", 《地球物理学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112346122A (zh) * 2020-11-06 2021-02-09 中国地震灾害防御中心 一种地震数据处理vda双参数成像方法
CN112346124A (zh) * 2020-11-06 2021-02-09 中国地震灾害防御中心 一种地震数据处理via双参数成像方法
CN112346125A (zh) * 2020-11-06 2021-02-09 中国地震灾害防御中心 一种地震数据处理vda双参数分析方法
CN112346123A (zh) * 2020-11-06 2021-02-09 中国地震灾害防御中心 一种地震数据处理via双参数分析方法
CN112346125B (zh) * 2020-11-06 2023-03-14 中国地震灾害防御中心 一种地震数据处理vda双参数分析方法
CN112346124B (zh) * 2020-11-06 2023-03-14 中国地震灾害防御中心 一种地震数据处理via双参数成像方法
CN112346122B (zh) * 2020-11-06 2023-03-14 中国地震灾害防御中心 一种地震数据处理vda双参数成像方法

Also Published As

Publication number Publication date
CN106950595B (zh) 2019-05-21

Similar Documents

Publication Publication Date Title
WO2015196779A1 (zh) 利用叠加速度求取高精度地震波速度的方法
US6014342A (en) Method of evaluating a subsurface region using gather sensitive data discrimination
CN106950595A (zh) 一种角变换成像方法与装置
CN102879819B (zh) 保持地震波场运动学特征的地震数据处理方法及装置
CN107367761B (zh) 一种大沙漠区深层碳酸盐岩储层广角反射波成像方法
CN105093292A (zh) 一种地震成像的数据处理方法和装置
CN104297784A (zh) 一种基于地震纵波方位各向异性的裂缝预测方法
CN106199704B (zh) 一种三维三分量海底电缆地震资料速度建模方法
CN108196305A (zh) 一种山地静校正方法
CN106970417B (zh) 椭圆展开转换波速度分析方法与系统
CN109839660A (zh) 一种利用叠前道集数据建立速度深度模型的方法
CN104533396A (zh) 一种远探测声波的处理方法
CN107656308B (zh) 一种基于时间深度扫描的共散射点叠前时间偏移成像方法
CN107390266A (zh) 基于角道集的速度更新方法及叠前深度偏移速度建模方法
CN106125139B (zh) 一种三维地震数据处理方法及系统
CN102053260B (zh) 获得地震纵波的方位速度的方法及处理地震数据的方法
CN104422959A (zh) 一种检测储层边界的曲率属性方法
CN106990434B (zh) 椭圆展开转换波成像方法及系统
CN106526664A (zh) 一种时空域几何变换水平叠加地震成像方法
CN105510958A (zh) 一种适用复杂介质的三维vsp观测系统设计方法
CN105425300B (zh) 一种剩余静校正方法
CN101963674A (zh) 共转换点抽道集方法
CN104516016B (zh) 一种确定三维转换波地震数据方位速度的方法及装置
CN107179548A (zh) 一种基于真地表的叠前地震成像方法
CN109307887A (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
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 100094 Room 401, building 4, courtyard 3, fengxiu Middle Road, Haidian District, Beijing

Patentee after: New JCM Group Co.,Ltd.

Address before: 100094 Room 401, Building 4, Courtyard 3, Fengxiu Middle Road, Haidian District, Beijing

Patentee before: LANDOCEAN ENERGY SERVICES CO.,LTD.