CN107918147B - 绕射波成像方法和装置 - Google Patents

绕射波成像方法和装置 Download PDF

Info

Publication number
CN107918147B
CN107918147B CN201711162018.1A CN201711162018A CN107918147B CN 107918147 B CN107918147 B CN 107918147B CN 201711162018 A CN201711162018 A CN 201711162018A CN 107918147 B CN107918147 B CN 107918147B
Authority
CN
China
Prior art keywords
imaging
data
target data
diffraction
point
Prior art date
Application number
CN201711162018.1A
Other languages
English (en)
Other versions
CN107918147A (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 CN201711162018.1A priority Critical patent/CN107918147B/zh
Publication of CN107918147A publication Critical patent/CN107918147A/zh
Application granted granted Critical
Publication of CN107918147B publication Critical patent/CN107918147B/zh

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction

Abstract

本发明提供了一种绕射波成像方法和装置,涉及绕射波成像的技术领域,该方法包括:获取待处理区域的叠前地震单炮数据;对待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,网格离散剖分后地表面的一个网格对应一个绕射成像点;从叠前地震单炮数据中选取目标数据,其中,一个目标数据为一个绕射成像点的共成像点道集数据;对目标数据进行时差校正,得到绕射波成像数据;通过绕射波成像数据进行绕射波成像。本发明缓解了绕射波对小尺度地质体定位的精确性较低的技术问题。

Description

绕射波成像方法和装置

技术领域

[0001] 本发明涉及绕射波成像技术领域,尤其是涉及一种绕射波成像方法和装置。

背景技术

[0002] 小尺度的不连续地质体,如断层、陷落柱、裂缝等,往往和矿产资源分布具有紧密 的关系,因而,对非均勾的不连续地质体进彳丁精确定位,可有效提尚对矿广资源的勘探成功 率,降低成本,避免可能发生的地质灾害,规避风险。

[0003] 绕射波是小尺度地质体的地震响应,绕射波携带了小尺度地质体的构造信息,因 而,绕射波可以被用来精确定位非均匀的不连续地质体,提供更强的地下空间照明度。但相 对反射波而言,绕射波在传播过程中衰减较快,能量较弱,覆盖次数较低,容易被反射波掩 盖,因此,绕射波对小尺度地质体定位的精确性较低。

发明内容

[0004] 有鉴于此,本发明的目的在于提供一种绕射波成像方法和装置,以缓解了绕射波 对小尺度地质体定位的精确性较低的技术问题。

[0005] 第一方面,本发明实施例提供了 一种绕射波成像方法,包括:

[0006] 获取待处理区域的叠前地震单炮数据;

[0007] 对所述待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,所 述网格离散剖分后地表面的一个网格对应一个所述绕射成像点;

[0008] 从所述叠前地震单炮数据中选取目标数据,其中,一个所述目标数据为一个所述 绕射成像点的共成像点道集数据;

[0009] 对所述目标数据进行时差校正,得到绕射波成像数据;

[0010] 通过所述绕射波成像数据进行绕射波成像。

[0011] 结合第一方面,本发明实施例提供了第一方面的第一种可能的实施方式,其中,从 所述叠前地震单炮数据中选取目标数据,包括:

[0012] 基于多个目标地震道组建所述绕射成像点的共成像点道集,其中,所述目标地震 道中所包含的检波点与所述绕射成像点的间隔距离小于预设长度;

[0013] 从所述叠前地震单炮数据中抽取所述共成像点道集的叠前地震道集,并将所述共 成像点道集的叠前地震道集确定为所述目标数据。

[0014] 结合第一方面的第一种可能的实施方式,本发明实施例提供了第一方面的第二种 可能的实施方式,其中,所述预设长度为

Figure CN107918147BD00041

[0015] 其中,r表示所述预设长度,V表示所述待处理区域内地震波的传播速度,to表示所 述绕射成像点的垂直双程走时,f表示地震波震源子波主频。

[0016] 结合第一方面,本发明实施例提供了第一方面的第三种可能的实施方式,其中,对 所述目标数据进行时差校正,得到绕射波成像数据,包括:

[0017] 计算第一水平距离、第二水平距离,其中,所述第一水平距离为检波器与所述目标 数据相应绕射成像点之间的水平距离,所述第二水平距离为炮点与所述目标数据相应绕射 成像点之间的水平距离;

[0018] 根据所述第一水平距离对所述目标数据对应的共成像点道集进行排列,并按照所 述排列顺序通过目标关系对所述目标数据相应共成像点道集进行时差校正,得到校正数 据,其中,所述目标关系为所述目标数据相应地震道的地震旅行时与所述第一水平距离、所 述第二水平距离的关系;

[0019] 对所述校正数据进行求和,得到所述绕射波成像数据。

[0020] 结合第一方面的第三种可能的实施方式,本发明实施例提供了第一方面的第四种 可能的实施方式,其中,所述目标关系通过以下公式进行表示:

[0021]

Figure CN107918147BD00051

[0022] 其中,t表示所述地震旅行时,to表示所述目标数据相应地震道成像点的垂直双程 走时,V表示所述待处理区域内地震波的传播速度,rd表示所述第一水平距离;sd表示所述 第二水平距离。

[0023] 第二方面,本发明实施例还提供一种绕射波成像装置,包括:

[0024] 获取模块,用于获取待处理区域的叠前地震单炮数据;

[0025] 确定模块,用于对所述待处理区域内的地下空间进行网格离散剖分,以确定绕射 成像点,其中,所述网格离散剖分后地表面的一个网格对应一个所述绕射成像点;

[0026] 选取模块,用于从所述叠前地震单炮数据中选取目标数据,其中,一个所述目标数 据为一个所述绕射成像点的共成像点道集数据;

[0027] 校正模块,用于对所述目标数据进行时差校正,得到绕射波成像数据;

[0028] 成像模块,用于通过所述绕射波成像数据进行绕射波成像。

[0029] 本发明实施例带来了以下有益效果:

[0030] 该绕射波成像方法包括:获取待处理区域的叠前地震单炮数据;对待处理区域内 的地下空间进行网格离散剖分,以确定绕射成像点,其中,一个网格对应一个绕射成像点; 从叠前地震单炮数据中选取目标数据,其中,一个目标数据为一个绕射成像点的共成像点 道集数据;对目标数据进行时差校正,得到绕射波成像数据;通过绕射波成像数据进行绕射 波成像。

[0031] 本发明实施例由绕射成像点的共成像点道集进行叠加成像,由于叠加成像所利用 的道集数据更加丰富,因而提高了覆盖次数,增强了数据信噪比,从而实现了对地下空间构 造的高精度成像,缓解了绕射波对小尺度地质体定位的精确性较低的技术问题。

[0032] 本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变 得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书 以及附图中所特别指出的结构来实现和获得。

[0033] 为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合 所附附图,作详细说明如下。

附图说明

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

[0035] 图1为本发明实施例一提供的一种绕射波成像方法的流程图;

[0036] 图2为地震波产生-接收示意图;

[0037] 图3为本发明实施例一提供的一种对目标数据进行时差校正,得到绕射波成像数 据的方法流程图;

[0038] 图4为本发明实施例二提供的一种绕射波成像装置的结构框图;

[0039] 图5为本发明实施例二提供的一种校正模块的结构框图。

[0040] 图标:100-获取模块;200-确定模块;300-选取模块;400-校正模块;401-计算单 元;402-校正单元;403-求和单元;500-成像模块。

具体实施方式

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

[0042] 绕射波可以用来定位非均匀不连续地质体,提供更强的地下空间照明度。但相对 反射波而言,绕射波在传播过程中衰减较快,能量较弱,覆盖次数较低,容易被反射波掩盖, 使得绕射波对小尺度地质体定位的精确性较低。基于此,本发明实施例提供的一种绕射点 成像方法和装置,可以提高覆盖次数,缓解绕射波对小尺度地质体定位的精确性较低的技 术问题。

[0043] 实施例一

[0044] 本发明实施例提供的一种绕射波成像方法,如图1所示,包括:

[0045] 步骤S102,获取待处理区域的叠前地震单炮数据。

[0046] 具体地,地震数据,一般以地震道为单位进行组织,该叠前地震单炮数据为包含时 间域地震数据的地震单炮数据,可以是包含有标准的SEGY道头信息,其中,SEGY道头中存储 炮点坐标、检波点坐标、成像点坐标。S E G Y格式是由勘探地球物理学家协会(S 〇 c i e t y 〇 f Exploration Geophysicists,简称SEG)提出的标准磁带数据格式之一,它是地震数据的最 为普遍的格式之一。

[0047] 步骤S104,对待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其 中,网格离散剖分后地表面的一个网格对应一个绕射成像点。

[0048] 具体地,在绕射点成像过程中,需要在地下区域选取一个成像点假定为反射点,然 后以这个成像点进行成像。在通过地震单炮数据进行成像时,需要使成像点遍历所需勘探 的待处理区域。由于待处理区域是一个三维空间,每个成像点的位置由地表面位置和地下 深度位置结合表征。本发明实施例中,把待处理区域的地表面划分成大量很小的网格,一个 网格对应一个绕射成像点,一个绕射成像点是具有相同地表面位置、不同地下深度位置的 大量成像点的代表。参见图2,点0表示本发明实施例中的绕射成像点,点D表示上述成像点, 点S表示炮点,点Ri和点R2表示两个不同的检波点。

[0049] 步骤S106,从叠前地震单炮数据中选取目标数据,其中,一个目标数据为一个绕射 成像点的共成像点道集数据。

[0050] 具体地,一个绕射成像点具有一组共成像点道集,叠前地震单炮数据中的一组共 成像点道集的地震数据即为一个目标数据,也就是说,一个绕射成像点关联一个目标数据。 叠前地震单炮数据中包括多个目标数据,从叠前地震单炮数据中选取出各个目标数据。

[0051] 步骤S108,对目标数据进行时差校正,得到绕射波成像数据。

[0052] 步骤SI 10,通过绕射波成像数据进行绕射波成像。

[0053] 在本发明实施例中,由绕射成像点的共成像点道集进行叠加成像,由于叠加成像 所利用的道集数据更加丰富,因而提高了覆盖次数,从而很好地缓解了绕射波在传播过程 中衰减较快、能量较弱、覆盖次数较低的问题,实现了增强数据信噪比和对地下空间构造的 高精度成像的目的,缓解了绕射波对小尺度地质体定位的精确性较低的技术问题。

[0054] 本发明实施例的一个可选实施方式中,步骤S106,从叠前地震单炮数据中选取目 标数据,包括:

[0055] 基于多个目标地震道组建绕射成像点的共成像点道集,其中,目标地震道中所包 含的检波点与绕射成像点的间隔距离小于预设长度;

[0056] 从叠前地震单炮数据中抽取共成像点道集的叠前地震道集,并将共成像点道集的 叠前地震道集确定为目标数据。

[0057] 具体地,本发明实施例的另一个可选实施方式中,预设长度为

[0058]

Figure CN107918147BD00071

[0059] 其中,r表示预设长度,V表示待处理区域内地震波的传播速度,to表示目标数据相 应地震道成像点的垂直双程走时,f表示地震波震源子波主频。

[0060] 需要说明的是,目标数据包含多个地震道的地震数据,每个地震到都有一个成像 点,对目标数据进行时差校正时,需要对各个地震道的地震数据进行时差校正。若将当前时 差校正涉及的地震道确定为目标地震道,则当前时差校正中目标数据相应地震道即指目标 地震道。参见图1,如果当前时差校正中目标数据相应地震道的成像点为点D,则to表示地震 波在点〇和点D之间传播的双程走时。

[0061] 此外,对于任一绕射成像点而言,预设长度由地震勘探的横向分辨率决定,在预设 长度选为

Figure CN107918147BD00072

的情况下,地震波旅行时差不大于二分之一个子波周期,在该预设长度 以内的地震道被认为来源于地下同一成像点。

[0062] 本发明实施例的另一个可选实施方式中,如图3所示,步骤S108,对目标数据进行 时差校正,得到绕射波成像数据,包括:

[0063] 步骤S301,计算第一水平距离、第二水平距离,其中,第一水平距离为检波器与目 标数据相应绕射成像点之间的水平距离,第二水平距离为炮点与目标数据相应绕射成像点 之间的水平距离。

[0064] 具体地,参见图1,第一水平距离即指点0到点R1的距离,或者,点0到点此的距离;第 二水平距离即指点0到点S的距离。

[0065] 可以通过第一公式计算第一水平距离,通过第二公式计算第二水平距离,其中,

[0066] 第一公式为:

Figure CN107918147BD00081

[0067] 第二公式为:

Figure CN107918147BD00082

[0068] cdpx表示成像点处的横向坐标,cdpy表示成像点处的纵向坐标,sx表示炮点位置 的横向坐标,sy表示炮点位置的纵向坐标,rx表示检波点位置处的横向坐标,ry表示检波点 位置处的纵向坐标。需要说明的是,地震波产生-接收过程中,炮和检波器一般放在地表面, 若将地表面(即,点S^gR1、点R2确定的平面)确定为三维坐标系的XY平面,则上述横向坐标 即XY平面中X方向的坐标,纵向坐标即XY平面中Y方向的坐标。

[0069] 步骤S302,根据第一水平距离对目标数据对应的共成像点道集进行排列,并按照 排列顺序通过目标关系对目标数据相应共成像点道集进行时差校正,得到校正数据,其中, 目标关系为目标数据相应地震道的地震旅行时与第一水平距离、第二水平距离的关系。

[0070] 具体地,按照排列顺序通过目标关系对目标数据相应共成像点道集进行时差校 正,从而拉平道集内的同相轴。

[0071] 步骤S303,对校正数据进行求和,得到绕射波成像数据。

[0072]具体地,对校正数据进行求和,得到绕射波成像数据可用公式

Figure CN107918147BD00083

进行表示,其中,W表示所述绕射波成像数据,u表示所述校正数据。 需要说明的是,本发明实施例中的叠前地震单炮数据是包含时间域地震数据和炮点坐标、 检波点坐标、成像点坐标的数据,即叠前地震单炮数据包括sd、rd和原始的时间数据,校正 数据包括sd、rd和t,其中,sd、rd和t的表示意义如本发明实施例中前文所述。

[0073] 本发明实施例的另一个可选实施方式中,目标关系通过以下公式进行表示:

[0074]

Figure CN107918147BD00084

[0075] 其中,t表示地震旅行时,to表示目标数据相应地震道成像点的垂直双程走时,V表 示待处理区域内地震波的传播速度,rd表示第一水平距离;sd表示第二水平距离。

[0076] 实施例二

[0077] 本发明实施例提供的一种绕射波成像装置,如图4所示,包括:

[0078] 获取模块100,用于获取待处理区域的叠前地震单炮数据;

[0079] 确定模块200,用于对待处理区域内的地下空间进行网格离散剖分,以确定绕射成 像点,其中,网格离散剖分后地表面的一个网格对应一个绕射成像点;

[0080] 选取模块300,用于从叠前地震单炮数据中选取目标数据,其中,一个目标数据为 一个绕射成像点的共成像点道集数据;

[0081] 校正模块400,用于对目标数据进行时差校正,得到绕射波成像数据;

[0082] 成像模块500,用于通过绕射波成像数据进行绕射波成像。

[0083] 在本发明实施例中,获取模块100获取待处理区域的叠前地震单炮数据;确定模块 200对待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,一个网格对 应一个绕射成像点;选取模块300从叠前地震单炮数据中选取目标数据,其中,一个目标数 据为一个绕射成像点的共成像点道集数据;校正模块400对目标数据进行时差校正,得到绕 射波成像数据;成像模块500通过绕射波成像数据进行绕射波成像。

[0084] 本发明实施例由绕射成像点的共成像点道集进行叠加成像,由于叠加成像所利用 的道集数据更加丰富,因而提高了覆盖次数,增强了数据信噪比,从而实现了对地下空间构 造的高精度成像,缓解了绕射波对小尺度地质体定位的精确性较低的技术问题。

[0085] 本发明实施例的一个可选实施方式中,选取模块用于:

[0086] 基于多个目标地震道组建绕射成像点的共成像点道集,其中,目标地震道中所包 含的检波点与绕射成像点的间隔距离小于预设长度;

[0087] 从叠前地震单炮数据中抽取共成像点道集的叠前地震道集,并将共成像点道集的 叠前地震道集确定为目标数据。

[0088] 本发明实施例的另一个可选实施方式中,预设长度为

Figure CN107918147BD00091

,其中,r表示预设 长度,V表示待处理区域内地震波的传播速度,to表示目标数据相应地震道成像点的垂直双 程走时,f表示地震波震源子波主频。

[0089] 本发明实施例的另一个可选实施方式中,校正模块,如图5所示,包括:

[0090] 计算单元401,用于计算第一水平距离、第二水平距离,其中,第一水平距离为检波 器与目标数据相应绕射成像点之间的水平距离,第二水平距离为炮点与目标数据相应绕射 成像点之间的水平距离;

[0091] 校正单元402,用于根据第一水平距离对目标数据对应的共成像点道集进行排列, 并按照排列顺序通过目标关系对目标数据相应共成像点道集进行时差校正,得到校正数 据,其中,目标关系为目标数据相应地震道的地震旅行时与第一水平距离、第二水平距离的 关系;

[0092] 求和单元403,用于对校正数据进行求和,得到绕射波成像数据。

[0093] 本发明实施例的另一个可选实施方式中,目标关系通过以下公式进行表示:

[0094]

Figure CN107918147BD00092

[0095] 其中,t表示地震旅行时,to表示目标数据相应地震道成像点的垂直双程走时,V表 示待处理区域内地震波的传播速度,rd表示第一水平距离;sd表示第二水平距离。

[0096] 本发明实施例所提供的绕射波成像方法及装置的计算机程序产品,包括存储了程 序代码的计算机可读存储介质,程序代码包括的指令可用于执行前面方法实施例中的方 法,具体实现可参见方法实施例,在此不再赘述。

[0097] 所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统 和装置的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。

[0098] 另外,在本发明实施例的描述中,除非另有明确的规定和限定,术语“安装”、“相 连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可 以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是 两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本 发明中的具体含义。

[0099] 所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以 存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说 对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计 算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个 人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。 而前述的存储介质包括:U盘、移动硬盘、只读存储器®0M,Read-0nly Memory)、随机存取存 储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

[0100] 在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、 “水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了 便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、 以特定的方位构造和操作,因此不能理解为对本发明的限制。

[0101] 此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相 对重要性。

[0102] 最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明 的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发 明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员 在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻 易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使 相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护 范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。

Claims (8)

1. 一种绕射波成像方法,其特征在于,包括: 获取待处理区域的叠前地震单炮数据; 对所述待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,所述网 格离散剖分后地表面的一个网格对应一个所述绕射成像点; 从所述叠前地震单炮数据中选取目标数据,其中,一个所述目标数据为一个所述绕射 成像点的共成像点道集数据; 对所述目标数据进行时差校正,得到绕射波成像数据; 通过所述绕射波成像数据进行绕射波成像; 其中,基于多个目标地震道组建所述绕射成像点的共成像点道集,其中,所述目标地震 道中所包含的检波点与所述绕射成像点的间隔距离小于预设长度; 从所述叠前地震单炮数据中抽取所述共成像点道集的叠前地震道集,并将所述共成像 点道集的叠前地震道集确定为所述目标数据。
2. 根据权利要求1所述的方法,其特征在于,所述预设长度为
Figure CN107918147BC00021
其中,r表示所述预设长度,V表示所述待处理区域内地震波的传播速度,to表示所述目 标数据相应地震道成像点的垂直双程走时,f表示地震波震源子波主频。
3. 根据权利要求1所述的方法,其特征在于,对所述目标数据进行时差校正,得到绕射 波成像数据,包括: 计算第一水平距离、第二水平距离,其中,所述第一水平距离为检波器与所述目标数据 相应绕射成像点之间的水平距离,所述第二水平距离为炮点与所述目标数据相应绕射成像 点之间的水平距离; 根据所述第一水平距离对所述目标数据对应的共成像点道集进行排列,并按照所述排 列的顺序通过目标关系对所述目标数据相应共成像点道集进行时差校正,得到校正数据, 其中,所述目标关系为所述目标数据相应地震道的地震旅行时与所述第一水平距离、所述 第二水平距离的关系; 对所述校正数据进行求和,得到所述绕射波成像数据。
4. 根据权利要求3所述的方法,其特征在于,所述目标关系通过以下公式进行表示:
Figure CN107918147BC00022
其中,t表示所述地震旅行时,to表示所述目标数据相应地震道成像点的垂直双程走时, V表示所述待处理区域内地震波的传播速度,rd表示所述第一水平距离;sd表示所述第二水 平距离。
5. —种绕射波成像装置,其特征在于,包括: 获取模块,用于获取待处理区域的叠前地震单炮数据; 确定模块,用于对所述待处理区域内的地下空间进行网格离散剖分,以确定绕射成像 点,其中,所述网格离散剖分后地表面的一个网格对应一个所述绕射成像点; 选取模块,用于从所述叠前地震单炮数据中选取目标数据,其中,一个所述目标数据为 一个所述绕射成像点的共成像点道集数据; 校正模块,用于对所述目标数据进行时差校正,得到绕射波成像数据; 成像模块,用于通过所述绕射波成像数据进行绕射波成像; 其中,所述选取模块用于: 基于多个目标地震道组建所述绕射成像点的共成像点道集,其中,所述目标地震道中 所包含的检波点与所述绕射成像点的间隔距离小于预设长度; 从所述叠前地震单炮数据中抽取所述共成像点道集的叠前地震道集,并将所述共成像 点道集的叠前地震道集确定为所述目标数据。
6. 根据权利要求5所述的装置,其特征在于,所述预设长度为
Figure CN107918147BC00031
其中,r表示所述预设长度,V表示所述待处理区域内地震波的传播速度,to表示所述目 标数据相应地震道成像点的垂直双程走时,f表示地震波震源子波主频。
7. 根据权利要求5所述的装置,其特征在于,所述校正模块,包括: 计算单元,用于计算第一水平距离、第二水平距离,其中,所述第一水平距离为检波器 与所述目标数据相应绕射成像点之间的水平距离,所述第二水平距离为炮点与所述目标数 据相应绕射成像点之间的水平距离; 校正单元,根据所述第一水平距离对所述目标数据对应的共成像点道集进行排列,并 按照所述排列的顺序通过目标关系对所述目标数据相应共成像点道集进行时差校正,得到 校正数据,其中,所述目标关系为所述目标数据相应地震道的地震旅行时与所述第一水平 距离、所述第二水平距离的关系; 求和单元,用于对所述校正数据进行求和,得到所述绕射波成像数据。
8. 根据权利要求7所述的装置,其特征在于,所述目标关系通过以下公式进行表示:
Figure CN107918147BC00032
其中,t表示所述地震旅行时,to表示所述目标数据相应地震道成像点的垂直双程走时, V表示所述待处理区域内地震波的传播速度,rd表示所述第一水平距离;sd表示所述第二水 平距离。
CN201711162018.1A 2017-11-20 2017-11-20 绕射波成像方法和装置 CN107918147B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711162018.1A CN107918147B (zh) 2017-11-20 2017-11-20 绕射波成像方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711162018.1A CN107918147B (zh) 2017-11-20 2017-11-20 绕射波成像方法和装置

Publications (2)

Publication Number Publication Date
CN107918147A CN107918147A (zh) 2018-04-17
CN107918147B true CN107918147B (zh) 2018-12-14

Family

ID=61896660

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711162018.1A CN107918147B (zh) 2017-11-20 2017-11-20 绕射波成像方法和装置

Country Status (1)

Country Link
CN (1) CN107918147B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102866421A (zh) * 2012-09-04 2013-01-09 中国科学院地质与地球物理研究所 识别小断距断点的散射波叠前成像方法
CN103941289A (zh) * 2014-05-09 2014-07-23 中国科学院地质与地球物理研究所 一种绕射点成像的方法和设备
CN104237940A (zh) * 2014-09-29 2014-12-24 中国石油天然气股份有限公司 一种基于动力学特征的绕射波成像方法及装置
CN105607121A (zh) * 2016-02-02 2016-05-25 中国矿业大学(北京) 一种煤炭陷落柱识别方法和装置
CN106772593A (zh) * 2016-11-10 2017-05-31 中国矿业大学(北京) 绕射波的成像方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7373252B2 (en) * 2005-11-04 2008-05-13 Western Geco L.L.C. 3D pre-stack full waveform inversion

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102866421A (zh) * 2012-09-04 2013-01-09 中国科学院地质与地球物理研究所 识别小断距断点的散射波叠前成像方法
CN103941289A (zh) * 2014-05-09 2014-07-23 中国科学院地质与地球物理研究所 一种绕射点成像的方法和设备
CN104237940A (zh) * 2014-09-29 2014-12-24 中国石油天然气股份有限公司 一种基于动力学特征的绕射波成像方法及装置
CN105607121A (zh) * 2016-02-02 2016-05-25 中国矿业大学(北京) 一种煤炭陷落柱识别方法和装置
CN106772593A (zh) * 2016-11-10 2017-05-31 中国矿业大学(北京) 绕射波的成像方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于PWD的绕射波波场分离成像方法综述;黄建平 等;《地球物理学进展》;20121215;第27卷(第6期);第2499-2510页 *
基于倾角-偏移距域道集的绕射波成像;李正伟 等;《地球物理学报》;20160415;第61卷(第4期);第1447-1459页 *

Also Published As

Publication number Publication date
CN107918147A (zh) 2018-04-17

Similar Documents

Publication Publication Date Title
EP2560030B1 (fr) Méthode pour construire un modèle géologique comportant une mise en position de depot d'unites stratigraphiques
Reyners et al. Imaging subduction from the trench to 300 km depth beneath the central North Island, New Zealand, with Vp and Vp/Vs
Dreger et al. Kinematic finite‐source model for the 24 August 2014 South Napa, California, earthquake from joint inversion of seismic, GPS, and InSAR data
Kent et al. Distribution of magma beneath the East Pacific Rise between the Clipperton transform and the 9 17′ N Deval from forward modeling of common depth point data
CN104237940B (zh) 一种基于动力学特征的绕射波成像方法及装置
Tronicke et al. Multivariate analysis of cross‐hole georadar velocity and attenuation tomograms for aquifer zonation
Beres et al. Integration of ground-penetrating radar and microgravimetric methods to map shallow caves
CN1950721B (zh) 关于碳氢化合物储藏的电磁勘探
US20180143334A1 (en) Seismic data processing
Zelt Modelling strategies and model assessment for wide-angle seismic traveltime data
US6643590B2 (en) Method for computing finite-frequency seismic migration traveltimes from monochromatic wavefields
EP2250520B1 (en) Method for indexing a subsurface volume for the purpose of inferring geologic information
CN102939546B (zh) 用于在地震处理中的局部属性匹配的系统和方法
CN100487488C (zh) 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法
US20130218539A1 (en) Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces
CN105492723A (zh) 勘探和开采环境下的电震勘测
CA2671434C (en) Systems and methods for measuring sea-bed resistivity
Méric et al. Seismic noise-based methods for soft-rock landslide characterization
CN102967882B (zh) 地层的层速度模型的建模方法
AU2009215000B2 (en) Method for determining geological information related to a subsurface volume of interest
CN102282481B (zh) 基于地震能见度分析的数据采集和叠前偏移
US20010032051A1 (en) Automated feature identification in data displays
CN102914791B (zh) 一种起伏地表地震数据处理的克希霍夫叠前时间偏移方法
EP2920615B1 (en) Process for creating image gathers
US20150066460A1 (en) Stratigraphic function

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