CN107918147A - 绕射波成像方法和装置 - Google Patents
绕射波成像方法和装置 Download PDFInfo
- Publication number
- CN107918147A CN107918147A CN201711162018.1A CN201711162018A CN107918147A CN 107918147 A CN107918147 A CN 107918147A CN 201711162018 A CN201711162018 A CN 201711162018A CN 107918147 A CN107918147 A CN 107918147A
- Authority
- CN
- China
- Prior art keywords
- data
- imaging
- target data
- diffraction
- seismic
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 174
- 238000000034 method Methods 0.000 claims abstract description 15
- 230000002889 sympathetic effect Effects 0.000 claims description 15
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 2
- 230000008901 benefit Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 210000002706 plastid Anatomy 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000006854 communication Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000003325 tomography Methods 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
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
技术领域
本发明涉及绕射波成像技术领域,尤其是涉及一种绕射波成像方法和装置。
背景技术
小尺度的不连续地质体,如断层、陷落柱、裂缝等,往往和矿产资源分布具有紧密的关系,因而,对非均匀的不连续地质体进行精确定位,可有效提高对矿产资源的勘探成功率,降低成本,避免可能发生的地质灾害,规避风险。
绕射波是小尺度地质体的地震响应,绕射波携带了小尺度地质体的构造信息,因而,绕射波可以被用来精确定位非均匀的不连续地质体,提供更强的地下空间照明度。但相对反射波而言,绕射波在传播过程中衰减较快,能量较弱,覆盖次数较低,容易被反射波掩盖,因此,绕射波对小尺度地质体定位的精确性较低。
发明内容
有鉴于此,本发明的目的在于提供一种绕射波成像方法和装置,以缓解了绕射波对小尺度地质体定位的精确性较低的技术问题。
第一方面,本发明实施例提供了一种绕射波成像方法,包括:
获取待处理区域的叠前地震单炮数据;
对所述待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,所述网格离散剖分后地表面的一个网格对应一个所述绕射成像点;
从所述叠前地震单炮数据中选取目标数据,其中,一个所述目标数据为一个所述绕射成像点的共成像点道集数据;
对所述目标数据进行时差校正,得到绕射波成像数据;
通过所述绕射波成像数据进行绕射波成像。
结合第一方面,本发明实施例提供了第一方面的第一种可能的实施方式,其中,从所述叠前地震单炮数据中选取目标数据,包括:
基于多个目标地震道组建所述绕射成像点的共成像点道集,其中,所述目标地震道中所包含的检波点与所述绕射成像点的间隔距离小于预设长度;
从所述叠前地震单炮数据中抽取所述共成像点道集的叠前地震道集,并将所述共成像点道集的叠前地震道集确定为所述目标数据。
结合第一方面的第一种可能的实施方式,本发明实施例提供了第一方面的第二种可能的实施方式,其中,所述预设长度为
其中,r表示所述预设长度,v表示所述待处理区域内地震波的传播速度,t0表示所述绕射成像点的垂直双程走时,f表示地震波震源子波主频。
结合第一方面,本发明实施例提供了第一方面的第三种可能的实施方式,其中,对所述目标数据进行时差校正,得到绕射波成像数据,包括:
计算第一水平距离、第二水平距离,其中,所述第一水平距离为检波器与所述目标数据相应绕射成像点之间的水平距离,所述第二水平距离为炮点与所述目标数据相应绕射成像点之间的水平距离;
根据所述第一水平距离对所述目标数据对应的共成像点道集进行排列,并按照所述排列顺序通过目标关系对所述目标数据相应共成像点道集进行时差校正,得到校正数据,其中,所述目标关系为所述目标数据相应地震道的地震旅行时与所述第一水平距离、所述第二水平距离的关系;
对所述校正数据进行求和,得到所述绕射波成像数据。
结合第一方面的第三种可能的实施方式,本发明实施例提供了第一方面的第四种可能的实施方式,其中,所述目标关系通过以下公式进行表示:
其中,t表示所述地震旅行时,t0表示所述目标数据相应地震道成像点的垂直双程走时,v表示所述待处理区域内地震波的传播速度,rd表示所述第一水平距离;sd表示所述第二水平距离。
第二方面,本发明实施例还提供一种绕射波成像装置,包括:
获取模块,用于获取待处理区域的叠前地震单炮数据;
确定模块,用于对所述待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,所述网格离散剖分后地表面的一个网格对应一个所述绕射成像点;
选取模块,用于从所述叠前地震单炮数据中选取目标数据,其中,一个所述目标数据为一个所述绕射成像点的共成像点道集数据;
校正模块,用于对所述目标数据进行时差校正,得到绕射波成像数据;
成像模块,用于通过所述绕射波成像数据进行绕射波成像。
本发明实施例带来了以下有益效果:
该绕射波成像方法包括:获取待处理区域的叠前地震单炮数据;对待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,一个网格对应一个绕射成像点;从叠前地震单炮数据中选取目标数据,其中,一个目标数据为一个绕射成像点的共成像点道集数据;对目标数据进行时差校正,得到绕射波成像数据;通过绕射波成像数据进行绕射波成像。
本发明实施例由绕射成像点的共成像点道集进行叠加成像,由于叠加成像所利用的道集数据更加丰富,因而提高了覆盖次数,增强了数据信噪比,从而实现了对地下空间构造的高精度成像,缓解了绕射波对小尺度地质体定位的精确性较低的技术问题。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一提供的一种绕射波成像方法的流程图;
图2为地震波产生-接收示意图;
图3为本发明实施例一提供的一种对目标数据进行时差校正,得到绕射波成像数据的方法流程图;
图4为本发明实施例二提供的一种绕射波成像装置的结构框图;
图5为本发明实施例二提供的一种校正模块的结构框图。
图标:100-获取模块;200-确定模块;300-选取模块;400-校正模块;401-计算单元;402-校正单元;403-求和单元;500-成像模块。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
绕射波可以用来定位非均匀不连续地质体,提供更强的地下空间照明度。但相对反射波而言,绕射波在传播过程中衰减较快,能量较弱,覆盖次数较低,容易被反射波掩盖,使得绕射波对小尺度地质体定位的精确性较低。基于此,本发明实施例提供的一种绕射点成像方法和装置,可以提高覆盖次数,缓解绕射波对小尺度地质体定位的精确性较低的技术问题。
实施例一
本发明实施例提供的一种绕射波成像方法,如图1所示,包括:
步骤S102,获取待处理区域的叠前地震单炮数据。
具体地,地震数据,一般以地震道为单位进行组织,该叠前地震单炮数据为包含时间域地震数据的地震单炮数据,可以是包含有标准的SEGY道头信息,其中,SEGY道头中存储炮点坐标、检波点坐标、成像点坐标。SEGY格式是由勘探地球物理学家协会(Society ofExploration Geophysicists,简称SEG)提出的标准磁带数据格式之一,它是地震数据的最为普遍的格式之一。
步骤S104,对待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,网格离散剖分后地表面的一个网格对应一个绕射成像点。
具体地,在绕射点成像过程中,需要在地下区域选取一个成像点假定为反射点,然后以这个成像点进行成像。在通过地震单炮数据进行成像时,需要使成像点遍历所需勘探的待处理区域。由于待处理区域是一个三维空间,每个成像点的位置由地表面位置和地下深度位置结合表征。本发明实施例中,把待处理区域的地表面划分成大量很小的网格,一个网格对应一个绕射成像点,一个绕射成像点是具有相同地表面位置、不同地下深度位置的大量成像点的代表。参见图2,点O表示本发明实施例中的绕射成像点,点D表示上述成像点,点S表示炮点,点R1和点R2表示两个不同的检波点。
步骤S106,从叠前地震单炮数据中选取目标数据,其中,一个目标数据为一个绕射成像点的共成像点道集数据。
具体地,一个绕射成像点具有一组共成像点道集,叠前地震单炮数据中的一组共成像点道集的地震数据即为一个目标数据,也就是说,一个绕射成像点关联一个目标数据。叠前地震单炮数据中包括多个目标数据,从叠前地震单炮数据中选取出各个目标数据。
步骤S108,对目标数据进行时差校正,得到绕射波成像数据。
步骤S110,通过绕射波成像数据进行绕射波成像。
在本发明实施例中,由绕射成像点的共成像点道集进行叠加成像,由于叠加成像所利用的道集数据更加丰富,因而提高了覆盖次数,从而很好地缓解了绕射波在传播过程中衰减较快、能量较弱、覆盖次数较低的问题,实现了增强数据信噪比和对地下空间构造的高精度成像的目的,缓解了绕射波对小尺度地质体定位的精确性较低的技术问题。
本发明实施例的一个可选实施方式中,步骤S106,从叠前地震单炮数据中选取目标数据,包括:
基于多个目标地震道组建绕射成像点的共成像点道集,其中,目标地震道中所包含的检波点与绕射成像点的间隔距离小于预设长度;
从叠前地震单炮数据中抽取共成像点道集的叠前地震道集,并将共成像点道集的叠前地震道集确定为目标数据。
具体地,本发明实施例的另一个可选实施方式中,预设长度为
其中,r表示预设长度,v表示待处理区域内地震波的传播速度,t0表示目标数据相应地震道成像点的垂直双程走时,f表示地震波震源子波主频。
需要说明的是,目标数据包含多个地震道的地震数据,每个地震到都有一个成像点,对目标数据进行时差校正时,需要对各个地震道的地震数据进行时差校正。若将当前时差校正涉及的地震道确定为目标地震道,则当前时差校正中目标数据相应地震道即指目标地震道。参见图1,如果当前时差校正中目标数据相应地震道的成像点为点D,则t0表示地震波在点O和点D之间传播的双程走时。
此外,对于任一绕射成像点而言,预设长度由地震勘探的横向分辨率决定,在预设长度选为的情况下,地震波旅行时差不大于二分之一个子波周期,在该预设长度以内的地震道被认为来源于地下同一成像点。
本发明实施例的另一个可选实施方式中,如图3所示,步骤S108,对目标数据进行时差校正,得到绕射波成像数据,包括:
步骤S301,计算第一水平距离、第二水平距离,其中,第一水平距离为检波器与目标数据相应绕射成像点之间的水平距离,第二水平距离为炮点与目标数据相应绕射成像点之间的水平距离。
具体地,参见图1,第一水平距离即指点O到点R1的距离,或者,点O到点R2的距离;第二水平距离即指点O到点S的距离。
可以通过第一公式计算第一水平距离,通过第二公式计算第二水平距离,其中,
第一公式为:
第二公式为:
cdpx表示成像点处的横向坐标,cdpy表示成像点处的纵向坐标,sx表示炮点位置的横向坐标,sy表示炮点位置的纵向坐标,rx表示检波点位置处的横向坐标,ry表示检波点位置处的纵向坐标。需要说明的是,地震波产生-接收过程中,炮和检波器一般放在地表面,若将地表面(即,点S、点R1、点R2确定的平面)确定为三维坐标系的XY平面,则上述横向坐标即XY平面中X方向的坐标,纵向坐标即XY平面中Y方向的坐标。
步骤S302,根据第一水平距离对目标数据对应的共成像点道集进行排列,并按照排列顺序通过目标关系对目标数据相应共成像点道集进行时差校正,得到校正数据,其中,目标关系为目标数据相应地震道的地震旅行时与第一水平距离、第二水平距离的关系。
具体地,按照排列顺序通过目标关系对目标数据相应共成像点道集进行时差校正,从而拉平道集内的同相轴。
步骤S303,对校正数据进行求和,得到绕射波成像数据。
具体地,对校正数据进行求和,得到绕射波成像数据可用公式进行表示,其中,w表示所述绕射波成像数据,u表示所述校正数据。需要说明的是,本发明实施例中的叠前地震单炮数据是包含时间域地震数据和炮点坐标、检波点坐标、成像点坐标的数据,即叠前地震单炮数据包括sd、rd和原始的时间数据,校正数据包括sd、rd和t,其中,sd、rd和t的表示意义如本发明实施例中前文所述。
本发明实施例的另一个可选实施方式中,目标关系通过以下公式进行表示:
其中,t表示地震旅行时,t0表示目标数据相应地震道成像点的垂直双程走时,v表示待处理区域内地震波的传播速度,rd表示第一水平距离;sd表示第二水平距离。
实施例二
本发明实施例提供的一种绕射波成像装置,如图4所示,包括:
获取模块100,用于获取待处理区域的叠前地震单炮数据;
确定模块200,用于对待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,网格离散剖分后地表面的一个网格对应一个绕射成像点;
选取模块300,用于从叠前地震单炮数据中选取目标数据,其中,一个目标数据为一个绕射成像点的共成像点道集数据;
校正模块400,用于对目标数据进行时差校正,得到绕射波成像数据;
成像模块500,用于通过绕射波成像数据进行绕射波成像。
在本发明实施例中,获取模块100获取待处理区域的叠前地震单炮数据;确定模块200对待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,一个网格对应一个绕射成像点;选取模块300从叠前地震单炮数据中选取目标数据,其中,一个目标数据为一个绕射成像点的共成像点道集数据;校正模块400对目标数据进行时差校正,得到绕射波成像数据;成像模块500通过绕射波成像数据进行绕射波成像。
本发明实施例由绕射成像点的共成像点道集进行叠加成像,由于叠加成像所利用的道集数据更加丰富,因而提高了覆盖次数,增强了数据信噪比,从而实现了对地下空间构造的高精度成像,缓解了绕射波对小尺度地质体定位的精确性较低的技术问题。
本发明实施例的一个可选实施方式中,选取模块用于:
基于多个目标地震道组建绕射成像点的共成像点道集,其中,目标地震道中所包含的检波点与绕射成像点的间隔距离小于预设长度;
从叠前地震单炮数据中抽取共成像点道集的叠前地震道集,并将共成像点道集的叠前地震道集确定为目标数据。
本发明实施例的另一个可选实施方式中,预设长度为其中,r表示预设长度,v表示待处理区域内地震波的传播速度,t0表示目标数据相应地震道成像点的垂直双程走时,f表示地震波震源子波主频。
本发明实施例的另一个可选实施方式中,校正模块,如图5所示,包括:
计算单元401,用于计算第一水平距离、第二水平距离,其中,第一水平距离为检波器与目标数据相应绕射成像点之间的水平距离,第二水平距离为炮点与目标数据相应绕射成像点之间的水平距离;
校正单元402,用于根据第一水平距离对目标数据对应的共成像点道集进行排列,并按照排列顺序通过目标关系对目标数据相应共成像点道集进行时差校正,得到校正数据,其中,目标关系为目标数据相应地震道的地震旅行时与第一水平距离、第二水平距离的关系;
求和单元403,用于对校正数据进行求和,得到绕射波成像数据。
本发明实施例的另一个可选实施方式中,目标关系通过以下公式进行表示:
其中,t表示地震旅行时,t0表示目标数据相应地震道成像点的垂直双程走时,v表示待处理区域内地震波的传播速度,rd表示第一水平距离;sd表示第二水平距离。
本发明实施例所提供的绕射波成像方法及装置的计算机程序产品,包括存储了程序代码的计算机可读存储介质,程序代码包括的指令可用于执行前面方法实施例中的方法,具体实现可参见方法实施例,在此不再赘述。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的***和装置的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
另外,在本发明实施例的描述中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种绕射波成像方法,其特征在于,包括:
获取待处理区域的叠前地震单炮数据;
对所述待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,所述网格离散剖分后地表面的一个网格对应一个所述绕射成像点;
从所述叠前地震单炮数据中选取目标数据,其中,一个所述目标数据为一个所述绕射成像点的共成像点道集数据;
对所述目标数据进行时差校正,得到绕射波成像数据;
通过所述绕射波成像数据进行绕射波成像。
2.根据权利要求1所述的方法,其特征在于,从所述叠前地震单炮数据中选取目标数据,包括:
基于多个目标地震道组建所述绕射成像点的共成像点道集,其中,所述目标地震道中所包含的检波点与所述绕射成像点的间隔距离小于预设长度;
从所述叠前地震单炮数据中抽取所述共成像点道集的叠前地震道集,并将所述共成像点道集的叠前地震道集确定为所述目标数据。
3.根据权利要求2所述的方法,其特征在于,所述预设长度为
<mrow>
<mi>r</mi>
<mo>=</mo>
<mi>v</mi>
<msqrt>
<mfrac>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mi>f</mi>
</mfrac>
</msqrt>
<mo>,</mo>
</mrow>
其中,r表示所述预设长度,v表示所述待处理区域内地震波的传播速度,t0表示所述目标数据相应地震道成像点的垂直双程走时,f表示地震波震源子波主频。
4.根据权利要求1所述的方法,其特征在于,对所述目标数据进行时差校正,得到绕射波成像数据,包括:
计算第一水平距离、第二水平距离,其中,所述第一水平距离为检波器与所述目标数据相应绕射成像点之间的水平距离,所述第二水平距离为炮点与所述目标数据相应绕射成像点之间的水平距离;
根据所述第一水平距离对所述目标数据对应的共成像点道集进行排列,并按照所述排列的顺序通过目标关系对所述目标数据相应共成像点道集进行时差校正,得到校正数据,其中,所述目标关系为所述目标数据相应地震道的地震旅行时与所述第一水平距离、所述第二水平距离的关系;
对所述校正数据进行求和,得到所述绕射波成像数据。
5.根据权利要求4所述的方法,其特征在于,所述目标关系通过以下公式进行表示:
<mrow>
<mi>t</mi>
<mo>=</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>rd</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>sd</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,t表示所述地震旅行时,t0表示所述目标数据相应地震道成像点的垂直双程走时,v表示所述待处理区域内地震波的传播速度,rd表示所述第一水平距离;sd表示所述第二水平距离。
6.一种绕射波成像装置,其特征在于,包括:
获取模块,用于获取待处理区域的叠前地震单炮数据;
确定模块,用于对所述待处理区域内的地下空间进行网格离散剖分,以确定绕射成像点,其中,所述网格离散剖分后地表面的一个网格对应一个所述绕射成像点;
选取模块,用于从所述叠前地震单炮数据中选取目标数据,其中,一个所述目标数据为一个所述绕射成像点的共成像点道集数据;
校正模块,用于对所述目标数据进行时差校正,得到绕射波成像数据;
成像模块,用于通过所述绕射波成像数据进行绕射波成像。
7.根据权利要求6所述的装置,其特征在于,所述选取模块用于:
基于多个目标地震道组建所述绕射成像点的共成像点道集,其中,所述目标地震道中所包含的检波点与所述绕射成像点的间隔距离小于预设长度;
从所述叠前地震单炮数据中抽取所述共成像点道集的叠前地震道集,并将所述共成像点道集的叠前地震道集确定为所述目标数据。
8.根据权利要求7所述的装置,其特征在于,所述预设长度为
<mrow>
<mi>r</mi>
<mo>=</mo>
<mi>v</mi>
<msqrt>
<mfrac>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mi>f</mi>
</mfrac>
</msqrt>
<mo>,</mo>
</mrow>
其中,r表示所述预设长度,v表示所述待处理区域内地震波的传播速度,t0表示所述目标数据相应地震道成像点的垂直双程走时,f表示地震波震源子波主频。
9.根据权利要求6所述的装置,其特征在于,所述校正模块,包括:
计算单元,用于计算第一水平距离、第二水平距离,其中,所述第一水平距离为检波器与所述目标数据相应绕射成像点之间的水平距离,所述第二水平距离为炮点与所述目标数据相应绕射成像点之间的水平距离;
校正单元,根据所述第一水平距离对所述目标数据对应的共成像点道集进行排列,并按照所述排列的顺序通过目标关系对所述目标数据相应共成像点道集进行时差校正,得到校正数据,其中,所述目标关系为所述目标数据相应地震道的地震旅行时与所述第一水平距离、所述第二水平距离的关系;
求和单元,用于对所述校正数据进行求和,得到所述绕射波成像数据。
10.根据权利要求9所述的装置,其特征在于,所述目标关系通过以下公式进行表示:
<mrow>
<mi>t</mi>
<mo>=</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mi>rd</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>sd</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,t表示所述地震旅行时,t0表示所述目标数据相应地震道成像点的垂直双程走时,v表示所述待处理区域内地震波的传播速度,rd表示所述第一水平距离;sd表示所述第二水平距离。
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 true CN107918147A (zh) | 2018-04-17 |
CN107918147B CN107918147B (zh) | 2018-12-14 |
Family
ID=61896660
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711162018.1A Active CN107918147B (zh) | 2017-11-20 | 2017-11-20 | 绕射波成像方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107918147B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110850469A (zh) * | 2019-11-20 | 2020-02-28 | 李志勇 | 一种基于克希霍夫积分解的地震槽波深度偏移的成像方法 |
CN111929730A (zh) * | 2020-08-20 | 2020-11-13 | 中国矿业大学(北京) | 小尺度地质异常体检测方法及装置 |
CN112415591A (zh) * | 2020-10-30 | 2021-02-26 | 中国石油天然气集团有限公司 | 绕射波的成像方法、装置、电子设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070203673A1 (en) * | 2005-11-04 | 2007-08-30 | Sherrill Francis G | 3d pre-stack full waveform inversion |
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 | 中国矿业大学(北京) | 绕射波的成像方法及装置 |
-
2017
- 2017-11-20 CN CN201711162018.1A patent/CN107918147B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070203673A1 (en) * | 2005-11-04 | 2007-08-30 | Sherrill Francis G | 3d pre-stack full waveform inversion |
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)
Title |
---|
李正伟 等: "基于倾角-偏移距域道集的绕射波成像", 《地球物理学报》 * |
黄建平 等: "基于PWD的绕射波波场分离成像方法综述", 《地球物理学进展》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110850469A (zh) * | 2019-11-20 | 2020-02-28 | 李志勇 | 一种基于克希霍夫积分解的地震槽波深度偏移的成像方法 |
CN111929730A (zh) * | 2020-08-20 | 2020-11-13 | 中国矿业大学(北京) | 小尺度地质异常体检测方法及装置 |
CN112415591A (zh) * | 2020-10-30 | 2021-02-26 | 中国石油天然气集团有限公司 | 绕射波的成像方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107918147B (zh) | 2018-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Martínez-Moreno et al. | Integrated geophysical methods for studying the karst system of Gruta de las Maravillas (Aracena, Southwest Spain) | |
CN107817523B (zh) | 绕射波偏移速度的分析方法及装置 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN106291725B (zh) | 一种快速反演地下地质体空间位置的方法 | |
EA008733B1 (ru) | Способ сейсморазведочных построений в сложных геологических формациях | |
CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN107918147B (zh) | 绕射波成像方法和装置 | |
CN104570068A (zh) | 一种基于三维gis的观测***crp属性分析方法 | |
CN101984366A (zh) | 基于共孔径面的起伏地表叠前时间偏移方法及设备 | |
CN106959467A (zh) | 地震波速度反演方法和装置 | |
CN109188520A (zh) | 薄储层厚度预测方法及装置 | |
CN104391319A (zh) | 一种地震资料采集***的确定方法及装置 | |
CN106646602A (zh) | 一种基于多种震源模型的震后快速生成震动图的方法 | |
CN101551464B (zh) | 一种复杂地表三维勘探初至波时间确定方法 | |
Lü et al. | Crustal-scale structure and deformation of Lu-Zong ore district: Joint interpretation from integrated geophysical data | |
CN106772593A (zh) | 绕射波的成像方法及装置 | |
CN104570070B (zh) | 一种建立二维近地表地质模型的方法及设备 | |
CN106054252A (zh) | 一种叠前时间偏移的方法及装置 | |
CN105093318A (zh) | 一种自适应波动方程波场延拓静校正方法 | |
CN104155690B (zh) | 基于椭球展开的三维地震数据叠加速度求取方法 | |
CN108051858B (zh) | 优化的地震勘探方法及装置 | |
Prado et al. | Geophysical investigation of the Colônia structure, Brazil | |
CN104950330B (zh) | 气顶油藏深度域成像的速度建模方法及*** | |
CN105093324B (zh) | 一种剩余静校正方法 | |
CN105259577B (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 |