CN103941289A - 一种绕射点成像的方法和设备 - Google Patents

一种绕射点成像的方法和设备 Download PDF

Info

Publication number
CN103941289A
CN103941289A CN201410195993.2A CN201410195993A CN103941289A CN 103941289 A CN103941289 A CN 103941289A CN 201410195993 A CN201410195993 A CN 201410195993A CN 103941289 A CN103941289 A CN 103941289A
Authority
CN
China
Prior art keywords
imaging
diffraction
point
data
described current
Prior art date
Application number
CN201410195993.2A
Other languages
English (en)
Other versions
CN103941289B (zh
Inventor
王真理
Original Assignee
中国科学院地质与地球物理研究所
北京中科联衡科技有限公司
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 中国科学院地质与地球物理研究所, 北京中科联衡科技有限公司 filed Critical 中国科学院地质与地球物理研究所
Priority to CN201410195993.2A priority Critical patent/CN103941289B/zh
Publication of CN103941289A publication Critical patent/CN103941289A/zh
Application granted granted Critical
Publication of CN103941289B publication Critical patent/CN103941289B/zh

Links

Abstract

本发明公开了一种绕射点成像的方法和设备。该方法包括:获取当前炮集数据,并计算地震波传播走时;选取一个未选用过的成像点作为当前绕射成像点;依据地震波传播走时,针对当前绕射成像点的地下深度,对在各个炮检距处的当前炮集数据进行时差校正、高通滤波及叠加,得到当前炮集数据在当前绕射成像点处的绕射点成像数据;返回执行所述选取一个未选用过的成像点作为当前绕射成像点,直至得到当前炮集数据在各个成像点处的绕射点成像数据;对各个成像点处的绕射点成像数据进行绕射点成像。通过本发明实施例的技术方案,可以使绕射点处绕射波信号突出,从而利用突出绕射点信号的数据成像之后可以更容易、更准确地识别出绕射点的位置。

Description

一种绕射点成像的方法和设备

技术领域

[0001] 本发明涉及地震勘探中的信号处理技术,特别是涉及一种绕射点成像的方法和设备。

背景技术

[0002] 在物理学中,把波绕过障碍物而传播的现象称为绕射,或衍射。在地震勘探工作中,地震波在传播过程中,遇到一些规模较小、物性横向变化较快、与突变相关的特殊地质体等地层剧烈变化的地方时,会产生绕射现象而产生绕射波。这些地层剧烈变化的地方称之为绕射点,例如断层的断点、地层尖灭点、溶洞、裂缝等都属于能够使地震波产生绕射波的绕射点。绕射点可视为一个新的震源,绕射波可视为绕射点位置的振动而产生的。由于绕射点所在的位置通常位于地下小尺度的地质体,而油气藏往往也与这些地质体相关联,因此,寻找地下绕射点的位置在地震勘探工作中是十分重要的。

[0003] 在地震勘探工作中,通常地面某一激发点(也称炮点)振动而在地下激发出地震波,再通过地面其他点接收到地震信号(即该炮点的炮集数据)来分析地下地质体的。其中,接收到的地震信号不仅包括绕射波,还包括有反射波等其他形式的地震波。由于地震信号中反射波的信号较强,目前最常用于分析地震信号的方法是反射波成像法。但是,正是由于绕射波的信号较弱且难以成像的很好,所以在对接收到的地震信号进行反射波成像的成像结果中,绕射点难以被发现。

[0004] 因此,为了识别绕射点位置,现有技术中主要采用的是先对接收到的地震信号区分反射波与绕射波,再对绕射波进行成像处理,最终实现从绕射波成像数据中分析绕射点位置。但是,由于绕射波与反射波在绕射点处相切,沿绕射波轨迹的一段空间内绕射波会被反射波干扰,这使得现有技术一直无法准确分离绕射波与反射波,分离出的绕射波信号一直不够准确,从而导致无法准确识别绕射点位置。

发明内容

[0005] 本发明所要解决的技术问题是,提供一种绕射点成像的方法和设备,通过先对炮集数据进行信号处理而得到突出绕射点位置的信号数据,再对处理后的信号数据进行成像处理,从而无需再对地震信号分离反射波与绕射波,以避免按照现有技术中反射波与绕射波分离不准确而导致的无法准确识别绕射点的位置。

[0006] 为解决上述技术问题,本发明提供了一种绕射点成像的方法和设备,该方法包括:

[0007] 获取一个未处理过的炮集数据作为当前炮集数据,并计算地震波传播走时;

[0008] 选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点;

[0009] 依据所述地震波传播走时,针对所述当前绕射成像点的地下深度,对在各个炮检距处的所述当前炮集数据进行时差校正、高通滤波及叠加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据;[0010] 返回执行所述选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点,直至得到所述当前炮集数据在各个成像点处的绕射点成像数据;

[0011] 对各个成像点处的绕射点成像数据进行绕射点成像。

[0012] 可选的,所述依据各个炮检距的地震波传播走时,针对所述当前绕射成像点的地下深度,对在各个炮检距的所述当前炮集数据进行动校正、高通滤波及叠加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据,包括:

[0013] 依据所述地震波传播走时,将在各个炮检距处的所述当前炮集数据时差校正到所述当前绕射成像点的双程垂直旅行时,得到所述当前炮集数据在各个炮检距处的校正数据;

[0014] 对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理,得到所述当前炮集数据在各个炮检距处的绕射点成像数据;

[0015] 将所述当前炮集数据在各个炮检距处的绕射点成像数据相加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据。

[0016] 可选的,还包括:

[0017] 响应于所述当前炮集数据具有多条接收线且每条接收线具有多个炮检距的接受传感器,选用二维滤波器执行所述对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理。

[0018] 可选的,还包括:

[0019] 响应于得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据,将所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据与所述当前绕射成像点处的总绕射点成像数据相加,重新得到所述当前绕射成像点处的总绕射点成像数据;所述当前绕射成像点处的总绕射点成像数据表示已处理过的所有炮集数据在所述当前绕射成像点处的绕射点成像数据之和;

[0020] 响应于得到所述当前炮集数据在各个成像点处的绕射点成像数据,返回执行所述获取一个未处理过的炮集数据作为当前炮集数据,直至得到所有炮集数据处理完毕,再执行所述对各个成像点处的绕射点成像数据进行绕射点成像。

[0021] 此外,本发明还提供了一种绕射点成像的设备,该设备包括:

[0022] 炮集选取模块,用于获取一个未处理过的炮集数据作为当前炮集数据;

[0023] 走时计算模块,用于计算地震波传播走时;

[0024] 成像点选取模块,用于选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点;

[0025] 信号处理模块,用于依据所述地震波传播走时,针对所述当前绕射成像点的地下深度,对在各个炮检距处的所述当前炮集数据进行时差校正、高通滤波及叠加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据;

[0026] 成像点选取触发模块,用于返回触发所述成像点选取模块再选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点,直至得到所述当前炮集数据在各个成像点处的绕射点成像数据;

[0027] 成像模块,用于对各个成像点处的绕射点成像数据进行绕射点成像。

[0028] 可选的,所述信号处理模块包括:[0029] 时差校正子模块,用于依据所述地震波传播走时,将在各个炮检距处的所述当前炮集数据时差校正到所述当前绕射成像点的双程垂直旅行时,得到所述当前炮集数据在各个炮检距处的校正数据;

[0030] 高通滤波子模块,用于对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理,得到所述当前炮集数据在各个炮检距处的绕射点成像数据。

[0031] 叠加子模块,用于将所述当前炮集数据在各个炮检距处的绕射点成像数据相加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据。

[0032] 可选的,还包括:

[0033] 滤波器确定模块,用于响应于所述当前炮集数据具有多条接收线且每条接收线具有多个炮检距的接受传感器,选用二维滤波器执行所述对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理。

[0034] 可选的,还包括:

[0035] 炮集叠加模块,用于响应于得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据,将所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据与所述当前绕射成像点处的总绕射点成像数据相加,重新得到所述当前绕射成像点处的总绕射点成像数据;所述当前绕射成像点处的总绕射点成像数据表示已处理过的所有炮集数据在所述当前绕射成像点处的绕射点成像数据之和;

[0036] 炮集选取触发模块,用于响应于得到所述当前炮集数据在各个成像点处的绕射点成像数据,返回执行所述获取一个未处理过的炮集数据作为当前炮集数据,直至得到所有炮集数据处理完毕,再执行所述对各个成像点处的绕射点成像数据进行绕射点成像。

[0037] 与现有技术相比,本发明具有以下优点:

[0038] 根据本发明实施例的技术方案,先对炮集数据进行时差校正、高通滤波及叠加的信号处理,然后再对处理后的数据进行成像。一方面,在绕射点处绕射波信号的相位会发生180度反转、振幅迅速增大,呈现出高频变化的特点,而在非绕射点处绕射波信号的振幅随空间变化则非常缓慢,呈现出低频变化的特点,另一方面,反射波在空间的连续性很强、振幅随空间的变化也呈现出低频变化的特点,因此,通过对炮集数据进行高通滤波,可以压制作为低频信号的反射波信号以及绕射点两边的绕射波信号,使作为高频信号的绕射点处绕射波信号得以突出,从而使得利用突出绕射点的数据进行成像之后,可以更容易、更准确地识别出绕射点的位置。

附图说明

[0039] 为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。

[0040] 图1为本发明中绕射点成像的方法实施例1的流程图;

[0041] 图2为绕射波随空间的振幅变化示意图;

[0042] 图3为本发明实施例中一种高通滤波后的绕射波振幅变化示例示意图;

[0043] 图4为本发明实施例中对炮集数据的信号处理过程一实施方式的流程图;[0044] 图5为本发明中绕射点成像的方法实施例2的流程图;

[0045] 图6为本发明中绕射点成像的设备实施例1的结构图;

[0046] 图7为本发明实施例中信号处理模块604 —实施方式的结构图;

[0047] 图8为本发明中绕射点成像的设备实施例2的结构图。

具体实施方式

[0048] 为了使本技术领域的人员更好地理解本申请方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。

[0049] 发明人经过研究发现,现有技术中之所以先分离反射波和绕射波再从绕射波中寻找绕射点位置,原因在于地震信号中的绕射波信号较弱而使得绕射点难以被发现。但是,由于绕射波与反射波在绕射点处是相切的,绕射轨迹的一段空间内绕射波会被反射波干扰,可见,反射波与绕射波实际上是难以被准确地分开的。与此同时,发明人还发现,在绕射点处绕射波信号的相位会发生180度反转、振幅迅速增大,呈现出高频变化的特点,而在非绕射点处绕射波信号的振幅随空间变化则非常缓慢,呈现出低频变化的特点,并且,反射波在空间的连续性很强、振幅随空间的变化也呈现出低频变化的特点。因此,地震信号中,在绕射点处的绕射波信号与绕射点两边的绕射波信号、整个空间的反射波信号呈现出不同的频谱特性,据此可以分离出绕射点处的绕射波信号。

[0050] 基于上述分析,本发明的主要思想是:先对炮集数据进行时差校正、高通滤波及叠加的信号处理,然后再对处理后的数据进行成像。这样,通过对炮集数据进行高通滤波,使作为低频信号的反射波信号以及绕射点两边的绕射波信号得以压制,并使作为高频信号的绕射点处绕射波信号得以突出,从而使得利用突出绕射点的数据进行成像之后,可以更容易、更准确地识别出绕射点的位置。

[0051] 在介绍了本发明的主要思想之后,下面结合附图,通过实施例来详细说明本发明中绕射点成像的方法和设备的实现方式。

[0052] 参见图1,示出了本发明中绕射点成像的方法实施例1的流程图。在本实施例中,例如具体可以包括如下步骤:

[0053] S101、获取一个未处理过的炮集数据作为当前炮集数据,并计算地震波传播走时。

[0054] 其中,本实施例中所指的一个炮集数据表示的是由同一个炮点激发地震波而在各个炮检距的接收传感器上接收的所有数据。一般地,地震勘探时会在多个不同的炮点多次激发地震波,每个炮点激发地震波以后可以在多个不同炮检距的接收传感器接收传播来的地震信号,因此,每一个炮点就会对应一组地震信号,即该炮点的炮集数据。其中,任意一个炮点的炮集数据都可以全面地反映地下绕射点的情况,即可利用任一炮集数据进行绕射点成像。可以理解的是,本实施例中,可以选用未处理过的炮集数据来作为当前炮集数据对绕射点进行成像,此时,在通过叠加多次绕射点成像而放大绕射点成像结果时,就可以利用不同炮集数据对同一地质体进行多次绕射点成像,从而减小激发位置对绕射点成像的影响。

[0055] 在获取当前炮集数据之后,由于各个炮检距之间的距离是不同的,各个炮检距的接收传感器接收到的地震信号数据存在时差。为了后续对各个炮检距处的数据进行信号处理时消除数据之间的时差,此时需要计算地震波传播走时。例如具体地,可以计算一个地震波传播走时表,即用表格形式表现的地震波传播距离与所需时间的关系。

[0056] S102、选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点。

[0057] 在绕射点成像过程中,需要在地下区域选取一个成像点假定为绕射点,然后对这个假定的绕射点进行成像,本实施例中所选取的当前绕射成像点即是这一假定的绕射点。

[0058] 需要说明的是,在对一个炮集数据进行成像时,需要使成像点遍历所需勘探的整个地下区域,为此,本实施例中需要针对整个地下区域的每个成像点而对数据进行信号处理。可以理解的是,由于需要遍历的地下区域是一个三维空间,该三维空间中的成像点位置可以是通过地面位置与地下深度位置两者来描述。因此,在本实施例的一些实施方式中,可以把地面所在的平面划分成大量很小的网格点,而在选取当前绕射成像点,可以先选择一个未信号处理过的网格点,然后遍历该网格点所有的地下深度位置来分别进行信号处理,再选择一个未信号处理过网格点,直至遍历所有的网格点为止,以此来完成遍历地下区域所有的成像点。其中,地下深度位置可以通过成像点处的地震波传播到地面所需的时间来表示,具体地,可以是以成像点处的地震波垂直旅行时。

[0059] S103、依据所述地震波传播走时,针对所述当前绕射成像点的地下深度,对在各个炮检距处的所述当前炮集数据进行时差校正、高通滤波及叠加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据。

[0060] 本实施例中,对数据的信号处理至少包括时差校正、高通滤波和叠加。其中,对不同炮检距的数据进行时差校正是为了校正各个炮检距的炮集数据之间存在的时差,从而使得各个炮检距的炮集数据可以进行叠加。而对不同炮检距的数据进行叠加则是为了提高信噪比而消除绕射点两边的噪声信号。

[0061] 需要说明的是,通常的成像过程中,主要是通过对炮集数据进行时差校正、叠加和偏移来完成的。在沿绕射轨迹求和的叠加过程中,由于绕射点两边的绕射波信号振幅符号相反(如图2所示)而互相抵消、互相压制,使得绕射点处的绕射波信号得以保留。但是,由于通常炮集数据的采集是关于绕射点不对称的,因此,绕射点以外的绕射波信号在成像过程中实际上并不能完全抵消,这样会使绕射点成像信噪比降低,影响成像质量,这也是反射成像数据中绕射信息弱的原因之一。而本实施例中,为了突出绕射点处的绕射波信号,可以先对时差校正过的炮集数据进行高通滤波,然后再对滤波后的数据进行叠加。由于绕射波信号的振幅仅在绕射点处处于高频变化而在绕射点两边处于低频变化,同时反射波的振幅随空间变化一直处于低频,因此,经过高通滤波的数据中反射波信号以及绕射点两边的绕射波都会被压制,而绕射点处的绕射波信号因被保留而得以突出。例如,图3示出了一种高通滤波后的绕射波振幅变化示例示意图,从图3中可以看出,滤波后绕射波信号仅在绕射点处具有振幅变化。其中,本实施例中所涉及的高通滤波,表示的是具有使低频受到压制而使高频相对放大或通过的信号滤波处理方式。

[0062] 可以理解的是,在本实施例的一些实施方式中,可以先对炮集数据进行时差校正,然后对校正过的数据进行高通滤波,再对滤波后的数据进行叠加,从而得到成像数据。其中,时差校正具体可以是将各个炮检距处的炮集数据都校正到当前绕射成像点的双程垂直旅行时。具体地,如图4所示,对炮集数据的信号处理过程,例如可以包括:[0063] S401、依据所述地震波传播走时,将在各个炮检距处的所述当前炮集数据时差校正到所述当前绕射成像点的双程垂直旅行时,得到所述当前炮集数据在各个炮检距处的校正数据。

[0064] 由于此时的信号处理是要为了针对当前绕射成像点进行成像,因此,可以将各个炮检距处的炮集数据都校正到该当前绕射成像点的双程垂直旅行时,从而实现将曲线拉平为直线。可以理解的是,由于各个炮检距处所接收到的数据需要校正的时差是不同的,因此,对炮集数据的时差校正也称为动校正。

[0065] S402、对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理,得到所述当前炮集数据在各个炮检距处的绕射点成像数据。

[0066] 其中,对校正数据的高通滤波处理,可以是通过预先设定为高通滤波器的数字滤波器来实现,该高通滤波器具体是指具有使低频受到压制而使高频相对放大或通过的滤波器。可以理解的是,对于不同的地质,地震信号的振幅空间变化可能有很大差别,并且,炮集数据采集的方式不同,也可能造成不同的地震信号的处理方式不同,因此,在对校正数据进行高通滤波处理时,可以依据炮集数据的情况来选择不同设置的数字滤波器。例如,一些实施方式中,还可以响应于所述当前炮集数据具有多条接收线且每条接收线具有多个炮检距的接受传感器,选用二维滤波器执行所述对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理。

[0067] S403、将所述当前炮集数据在各个炮检距处的绕射点成像数据相加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据。

[0068] 在对炮集数据进行高通滤波之后,可以将各个炮检距处的滤波后数据相加而得到绕射点成像数据,以提高成像数据的信噪比,使得绕射点处的信号通过叠加得以进一步突出,而使绕射点以外的噪声信号通过叠加得以消除。

[0069] 接着返回图1。在S103执彳丁完成以后,接着进入S104的执打。

[0070] S104、返回执行所述选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点,直至得到所述当前炮集数据在各个成像点处的绕射点成像数据。

[0071] 在绕射点成像过程中,需要遍历所有的成像点来获得各个成像点的绕射点成像数据。因此,在得到所述当前炮集数据在当前绕射成像点处的绕射点成像数据之后,还需要返回S102来重新选取一个未针对该炮集数据选用过的成像点,从而得到所述当前炮集数据在重新选取的成像点处的绕射点成像数据,直至得到所述当前炮集数据在各个成像点处的绕射点成像数据为止。

[0072] 需要说明的是,在前述以网格点及深度位置表示成像点位置的实施方式基础上,可以先针对同一个网格点选择不同的深度位置来获得绕射点成像数据,在遍历该网格点的所有深度位置之后,再选择另一个网格点来获得绕射点成像数据,直至遍历所有网格点,从而得到各个深度位置的各个网格点处的绕射点成像数据。

[0073] S105、对各个成像点处的绕射点成像数据进行绕射点成像。

[0074] 可以理解的是,由于绕射点成像数据经过了高通滤波,因此,绕射点两边的信号可以得到压制,而绕射点处的信号得以突出,从而使得成像结果中的绕射点位置得以突出。

[0075] 在本实施的一些实施方式中,由于单一炮点的炮集数据中绕射点处的绕射波信号往往较弱,即使经过了前述高通滤波处理,绕射点处的绕射波信号可能依然难以与绕射点以外的噪声信号区分。为了进一步提高绕射点处信号的信噪比,可以将多个不同炮点的炮集数据经过高通滤波之后而在同一成像点处的成像数据进行叠加,从而使得绕射点处的信号进一步突出,绕射点以外的噪声信号进一步消除。具体地,例如在S103执行完成之后,还可以响应于得到所述当前炮集数据在所述当前绕射成像点处的总绕射点成像数据,将所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据与所述当前绕射成像点处的绕射点成像数据相加,重新得到所述当前绕射成像点处的总绕射点成像数据;而在S104执行完成之后,还可以响应于得到所述当前炮集数据在各个成像点处的绕射点成像数据,返回执行SlOl获取一个未处理过的炮集数据作为当前炮集数据,直至得到所有炮集数据处理完毕,再执行S105对各个成像点处的绕射点成像数据进行绕射点成像。其中,所述当前绕射成像点处的总绕射点成像数据表示已处理过的所有炮集数据在所述当前绕射成像点处的绕射点成像数据之和。可以理解的是,此时S105用于绕射点成像的数据为各个成像点的总绕射点成像数据。

[0076] 本实施例的技术方案,通过对炮集数据进行高通滤波,使作为低频信号的反射波信号以及绕射点两边的绕射波信号得以压制,并使作为高频信号的绕射点处绕射波信号得以突出,从而使得利用突出绕射点的数据进行成像之后,可以更容易、更准确地识别出绕射点的位置。

[0077] 为了使本领域技术人员更加清楚地理解本发明的应用方式,下面以一个应用场景来对本发明的实施方式进行示例性地介绍。在该应用场景中,绕射点的成像是利用多个不同炮点的炮集数据叠加而实现的,而成像点是由预先离散化地质体而得到的网格点和深度位置来表示的。

[0078] 参见图5,示出了本发明中绕射点成像的方法实施例2的流程图。在本实施例中,例如具体可以包括如下步骤:

[0079] S501、获取用于绕射点成像的所有炮集数据,计算地震波传播走时表,并选择高通滤波器。

[0080] S502、获取一个未处理过的炮集数据作为当前炮集数据。

[0081 ] S503、选取一个未针对当前炮集数据而选用过的网格点作为当前网格点。

[0082] S504、选取一个未针对当前炮集数据及当前网格点而选用过的深度作为当前网格点的当前深度。

[0083] S505、依据地震波传播走时,将在各个炮检距处的当前炮集数据时差校正到在当前深度下当前网格点的双程垂直旅行时,得到当前炮集数据在各个炮检距处的校正数据。

[0084] S506、采用对当前炮集数据在各个炮检距处的校正数据进行高通滤波处理,得到当前炮集数据在各个炮检距处的绕射点成像数据。

[0085] S507、将当前炮集数据在各个炮检距处的绕射点成像数据相加,得到当前炮集数据在当前深度下当前网格点的绕射点成像数据并存储。

[0086] S508、将当前炮集数据在当前深度下当前网格点的绕射点成像数据与当前深度下当前网格点的总绕射点成像数据相加,重新得到当前深度下当前网格点的总绕射点成像数据。

[0087] S509、迭代S504〜S508,直至对当前网格点遍历所有深度。

[0088] S510、迭代S503〜S509,直至对当前炮集数据遍历所有当前网格点。[0089] S511、迭代S502〜S510,直至遍历所有炮集数据。

[0090] S512、对各深度下各个网格点的总绕射点成像数据进行绕射点成像。

[0091] 本实施例的技术方案,一方面,通过对炮集数据进行高通滤波,使作为低频信号的反射波信号以及绕射点两边的绕射波信号得以压制,并使作为高频信号的绕射点处绕射波信号得以突出,另一方面,通过将多个不同炮点的炮集数据在同一成像点的成像数据进行了叠加,提高了绕射点处信号的信噪比,进一步突出了绕射点处的信号,从而使得利用突出绕射点的数据进行成像之后,可以更容易、更准确地识别出绕射点的位置。

[0092] 对应于方法实施例,本发明还提供了一种绕射点成像设备。

[0093] 参见图6,示出了本发明中绕射点成像的设备实施例1的结构图。在本实施例中,所述设备例如具体可以包括:

[0094] 炮集选取模块601,用于获取一个未处理过的炮集数据作为当前炮集数据;

[0095] 走时计算模块602,用于计算地震波传播走时;

[0096] 成像点选取模块603,用于选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点;

[0097] 信号处理模块604,用于依据所述地震波传播走时,针对所述当前绕射成像点的地下深度,对在各个炮检距处的所述当前炮集数据进行时差校正、高通滤波及叠加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据;

[0098] 成像点选取触发模块605,用于返回触发所述成像点选取模块再选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点,直至得到所述当前炮集数据在各个成像点处的绕射点成像数据;

[0099] 成像模块606,用于对各个成像点处的绕射点成像数据进行绕射点成像。

[0100] 参见图7,示出了本发明实施例中信号处理模块604—实施方式的结构图。在本实施方式中,所述信号处理模块604例如可以包括:

[0101] 时差校正子模块701,用于依据所述地震波传播走时,将在各个炮检距处的所述当前炮集数据时差校正到所述当前绕射成像点的双程垂直旅行时,得到所述当前炮集数据在各个炮检距处的校正数据;

[0102] 高通滤波子模块702,用于对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理,得到所述当前炮集数据在各个炮检距处的绕射点成像数据;

[0103] 叠加子模块703,用于将所述当前炮集数据在各个炮检距处的绕射点成像数据相力口,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据。

[0104] 进一步而言,可选的,在图7所示的实施方式基础上,所述设备例如还可以包括:

[0105] 滤波器确定模块,用于响应于所述当前炮集数据具有多条接收线且每条接收线具有多个炮检距的接受传感器,选用二维滤波器执行所述对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理。

[0106] 参见图8,示出了本发明中绕射点成像的设备实施例2的结构图。在本实施例中,除了图6所示的所有结构以外,所述设备例如还可以包括:

[0107] 炮集叠加模块801,用于响应于得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据,将所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据与所述当前绕射成像点处的总绕射点成像数据相加,重新得到所述当前绕射成像点处的总绕射点成像数据;所述当前绕射成像点处的总绕射点成像数据表示已处理过的所有炮集数据在所述当前绕射成像点处的绕射点成像数据之和;

[0108] 炮集选取触发模块802,用于响应于得到所述当前炮集数据在各个成像点处的绕射点成像数据,返回执行所述获取一个未处理过的炮集数据作为当前炮集数据,直至得到所有炮集数据处理完毕,再执行所述对各个成像点处的绕射点成像数据进行绕射点成像。

[0109] 本发明设备实施例的技术方案,通过对炮集数据进行高通滤波,使作为低频信号的反射波信号以及绕射点两边的绕射波信号得以压制,并使作为高频信号的绕射点处绕射波信号得以突出,从而使得利用突出绕射点的数据进行成像之后,可以更容易、更准确地识别出绕射点的位置。

[0110] 需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

[0111] 对于设备实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的设备实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。

[0112] 以上所述仅是本申请的具体实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本申请原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本申请的保护范围。

Claims (8)

1.一种绕射点成像的方法,其特征在于,包括: 获取一个未处理过的炮集数据作为当前炮集数据,并计算地震波传播走时; 选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点; 依据所述地震波传播走时,针对所述当前绕射成像点的地下深度,对在各个炮检距处的所述当前炮集数据进行时差校正、高通滤波及叠加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据; 返回执行所述选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点,直至得到所述当前炮集数据在各个成像点处的绕射点成像数据; 对各个成像点处的绕射点成像数据进行绕射点成像。
2.根据权利要求1所述的方法,其特征在于,所述依据各个炮检距的地震波传播走时,针对所述当前绕射成像点的地下深度,对在各个炮检距的所述当前炮集数据进行动校正、高通滤波及叠加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据,包括: 依据所述地震波传播走时,将在各个炮检距处的所述当前炮集数据时差校正到所述当前绕射成像点的双程垂直旅行时,得到所述当前炮集数据在各个炮检距处的校正数据;对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理,得到所述当前炮集数据在各个炮检距处的绕射点成像数据; 将所述当前炮集数据在各个炮检距处的绕射点成像数据相加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据。
3.根据权利要求2所述的方法,其特征在于,还包括: 响应于所述当前炮集数据具有多条接收线且每条接收线具有多个炮检距的接受传感器,选用二维滤波器执行所述对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理。
4.根据权利要求1所述的方法,其特征在于,还包括: 响应于得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据,将所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据与所述当前绕射成像点处的总绕射点成像数据相加,重新得到所述当前绕射成像点处的总绕射点成像数据;所述当前绕射成像点处的总绕射点成像数据表示已处理过的所有炮集数据在所述当前绕射成像点处的绕射点成像数据之和; 响应于得到所述当前炮集数据在各个成像点处的绕射点成像数据,返回执行所述获取一个未处理过的炮集数据作为当前炮集数据,直至得到所有炮集数据处理完毕,再执行所述对各个成像点处的绕射点成像数据进行绕射点成像。
5.一种绕射点成像的设备,其特征在于,包括: 炮集选取模块,用于获取一个未处理过的炮集数据作为当前炮集数据; 走时计算模块,用于计算地震波传播走时; 成像点选取模块,用于选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点; 信号处理模块,用于依据所述地震波传播走时,针对所述当前绕射成像点的地下深度,对在各个炮检距处的所述当前炮集数据进行时差校正、高通滤波及叠加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据; 成像点选取触发模块,用于返回触发所述成像点选取模块再选取一个未针对所述当前炮集数据而选用过的成像点作为当前绕射成像点,直至得到所述当前炮集数据在各个成像点处的绕射点成像数据; 成像模块,用于对各个成像点处的绕射点成像数据进行绕射点成像。
6.根据权利要求5所述的设备,其特征在于,所述信号处理模块包括: 时差校正子模块,用于依据所述地震波传播走时,将在各个炮检距处的所述当前炮集数据时差校正到所 述当前绕射成像点的双程垂直旅行时,得到所述当前炮集数据在各个炮检距处的校正数据; 高通滤波子模块,用于对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理,得到所述当前炮集数据在各个炮检距处的绕射点成像数据。 叠加子模块,用于将所述当前炮集数据在各个炮检距处的绕射点成像数据相加,得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据。
7.根据权利要求6所述的设备,其特征在于,还包括: 滤波器确定模块,用于响应于所述当前炮集数据具有多条接收线且每条接收线具有多个炮检距的接受传感器,选用二维滤波器执行所述对所述当前炮集数据在各个炮检距处的校正数据进行高通滤波处理。
8.根据权利要求5所述的设备,其特征在于,还包括: 炮集叠加模块,用于响应于得到所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据,将所述当前炮集数据在所述当前绕射成像点处的绕射点成像数据与所述当前绕射成像点处的总绕射点成像数据相加,重新得到所述当前绕射成像点处的总绕射点成像数据;所述当前绕射成像点处的总绕射点成像数据表示已处理过的所有炮集数据在所述当前绕射成像点处的绕射点成像数据之和; 炮集选取触发模块,用于响应于得到所述当前炮集数据在各个成像点处的绕射点成像数据,返回执行所述获取一个未处理过的炮集数据作为当前炮集数据,直至得到所有炮集数据处理完毕,再执行所述对各个成像点处的绕射点成像数据进行绕射点成像。
CN201410195993.2A 2014-05-09 2014-05-09 一种绕射点成像的方法和设备 CN103941289B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410195993.2A CN103941289B (zh) 2014-05-09 2014-05-09 一种绕射点成像的方法和设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410195993.2A CN103941289B (zh) 2014-05-09 2014-05-09 一种绕射点成像的方法和设备

Publications (2)

Publication Number Publication Date
CN103941289A true CN103941289A (zh) 2014-07-23
CN103941289B CN103941289B (zh) 2017-01-25

Family

ID=51189034

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410195993.2A CN103941289B (zh) 2014-05-09 2014-05-09 一种绕射点成像的方法和设备

Country Status (1)

Country Link
CN (1) CN103941289B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237940A (zh) * 2014-09-29 2014-12-24 中国石油天然气股份有限公司 一种基于动力学特征的绕射波成像方法及装置
CN105806949A (zh) * 2016-03-07 2016-07-27 杭州华新检测技术股份有限公司 一种超声ct检测中的射线走时校正方法
CN107918147A (zh) * 2017-11-20 2018-04-17 中国矿业大学(北京) 绕射波成像方法和装置
CN109116423A (zh) * 2018-07-11 2019-01-01 北京奥能恒业能源技术有限公司 一种绕射多次波压制方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102455439A (zh) * 2010-11-02 2012-05-16 中国石油大学(北京) 基于克希霍夫积分法的绕射波场分离方法
CN102520444A (zh) * 2011-12-13 2012-06-27 中国科学院地质与地球物理研究所 一种叠后地震波中绕射波信息提取方法
CA2841775A1 (en) * 2011-07-19 2013-01-24 Halliburton Energy Services, Inc. System and method for moment tensor migration imaging
CN103675897A (zh) * 2012-08-30 2014-03-26 中国石油化工股份有限公司 一种地震绕射波分离成像方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102455439A (zh) * 2010-11-02 2012-05-16 中国石油大学(北京) 基于克希霍夫积分法的绕射波场分离方法
CA2841775A1 (en) * 2011-07-19 2013-01-24 Halliburton Energy Services, Inc. System and method for moment tensor migration imaging
CN102520444A (zh) * 2011-12-13 2012-06-27 中国科学院地质与地球物理研究所 一种叠后地震波中绕射波信息提取方法
CN103675897A (zh) * 2012-08-30 2014-03-26 中国石油化工股份有限公司 一种地震绕射波分离成像方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
KLKOV A, ET AL.: ""Separation and imaging of seismic diffractions in the dip-angle domain"", 《GEOPHYSICS》, vol. 77, no. 6, 17 September 2012 (2012-09-17), pages 131 - 143 *
刘玉金 等: ""绕射波叠前时间偏移速度分析成像"", 《地球物理学进展》, vol. 28, no. 6, 31 December 2013 (2013-12-31), pages 3023 - 3029 *
孔雪 等: ""基于平面波记录的绕射目标成像方法研究"", 《石油地球物理勘探》, vol. 47, no. 4, 15 August 2012 (2012-08-15), pages 674 - 681 *
朱生旺 等: ""通过压制共散射点道集映射噪声改善绕射波成像分辨率"", 《石油物探》, vol. 49, no. 2, 15 March 2010 (2010-03-15), pages 107 - 114 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237940A (zh) * 2014-09-29 2014-12-24 中国石油天然气股份有限公司 一种基于动力学特征的绕射波成像方法及装置
CN104237940B (zh) * 2014-09-29 2017-03-08 中国石油天然气股份有限公司 一种基于动力学特征的绕射波成像方法及装置
CN105806949A (zh) * 2016-03-07 2016-07-27 杭州华新检测技术股份有限公司 一种超声ct检测中的射线走时校正方法
CN105806949B (zh) * 2016-03-07 2019-04-02 杭州华新检测技术股份有限公司 一种超声ct检测中的射线走时校正方法
CN107918147A (zh) * 2017-11-20 2018-04-17 中国矿业大学(北京) 绕射波成像方法和装置
CN107918147B (zh) * 2017-11-20 2018-12-14 中国矿业大学(北京) 绕射波成像方法和装置
CN109116423A (zh) * 2018-07-11 2019-01-01 北京奥能恒业能源技术有限公司 一种绕射多次波压制方法及装置

Also Published As

Publication number Publication date
CN103941289B (zh) 2017-01-25

Similar Documents

Publication Publication Date Title
Baeten et al. The use of low frequencies in a full‐waveform inversion and impedance inversion land seismic case study
Foti et al. Guidelines for the good practice of surface wave analysis: A product of the InterPACIFIC project
Duputel et al. W phase source inversion for moderate to large earthquakes (1990–2010)
Sabbione et al. Automatic first-breaks picking: New strategies and algorithms
AU2011206597B2 (en) Method to process marine seismic data
Yuan et al. Multiscale adjoint waveform tomography for surface and body waves
Ross et al. Automatic picking of direct P, S seismic phases and fault zone head waves
Taner et al. Velocity spectra—Digital computer derivation applications of velocity functions
Poggi et al. Estimating Rayleigh wave particle motion from three-component array analysis of ambient vibrations
Forbriger Inversion of shallow-seismic wavefields: I. Wavefield transformation
Yao et al. Surface-wave array tomography in SE Tibet from ambient seismic noise and two-station analysis—I. Phase velocity maps
Mayne Common reflection point horizontal data stacking techniques
Chambers et al. Testing the ability of surface arrays to monitor microseismic activity
Arai et al. S-wave velocity profiling by inversion of microtremor H/V spectrum
EP0889331B1 (fr) Méthode pour modéliser en 3D l'impédance d'un milieu hétérogène
Gajewski et al. Reverse modelling for seismic event characterization
Grigoli et al. Automated seismic event location by waveform coherence analysis
EP1096270B1 (fr) Méthode pour modeliser en 2d ou 3d un milieu hétérogène tel que le sous-sol décrit par un ou plusieurs paramètres physiques
Cassidy et al. Ground penetrating radar data processing, modelling and analysis
Draganov et al. Seismic exploration‐scale velocities and structure from ambient seismic noise (> 1 Hz)
Slob et al. Interferometric electromagnetic Green's functions representations using propagation invariants
Picozzi et al. Characterization of shallow geology by high-frequency seismic noise tomography
CN101086535B (zh) 表征油藏随时间演变的方法和程序
Xia et al. Advantages of using multichannel analysis of Love waves (MALW) to estimate near-surface shear-wave velocity
Xia et al. Feasibility of detecting near-surface feature with Rayleigh-wave diffraction

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20150522

Address after: 100029 Beijing city Chaoyang District Beitucheng West Road No. 19

Applicant after: Institute of Geology and Geophysics, Chinese Academy of Sciences

Address before: 100029 Beijing city Chaoyang District Beitucheng West Road No. 19

Applicant before: Institute of Geology and Geophysics, Chinese Academy of Sciences

Applicant before: Sino-Lionheart Technology Co., Ltd.

ASS Succession or assignment of patent right

Free format text: FORMER OWNER: BEIJING ZHONGKE LIANHENG TECHNOLOGY CO., LTD.

Effective date: 20150522

C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20181029

Address after: 100083 room 707, block B, jincode building, 38, Xue Qing Road, Haidian District, Beijing.

Patentee after: Guo Kan Petrotechnik Ltd.

Address before: No. 19, North Tu Cheng West Road, Chaoyang District, Beijing

Patentee before: Institute of Geology and Geophysics, Chinese Academy of Sciences