CN106772592A - 绕射波聚焦能量的分析方法及装置 - Google Patents
绕射波聚焦能量的分析方法及装置 Download PDFInfo
- Publication number
- CN106772592A CN106772592A CN201610988823.9A CN201610988823A CN106772592A CN 106772592 A CN106772592 A CN 106772592A CN 201610988823 A CN201610988823 A CN 201610988823A CN 106772592 A CN106772592 A CN 106772592A
- Authority
- CN
- China
- Prior art keywords
- diffracted wave
- data
- wave
- common offset
- diffracted
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 32
- 239000013598 vector Substances 0.000 claims abstract description 91
- 238000003384 imaging method Methods 0.000 claims abstract description 41
- 238000000034 method Methods 0.000 claims abstract description 37
- 238000013508 migration Methods 0.000 claims abstract description 33
- 230000005012 migration Effects 0.000 claims abstract description 33
- 238000000605 extraction Methods 0.000 claims description 21
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 239000000284 extract Substances 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 241001269238 Data Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 239000002360 explosive Substances 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 210000002706 plastid Anatomy 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000001914 filtration Methods 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
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000017105 transposition Effects 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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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
技术领域
本发明涉及高分辨率地震勘探技术领域,尤其是涉及一种绕射波聚焦能量的分析方法及装置。
背景技术
在地震勘探中,对于勘探人员来说,对定位断层、不整合、小断块、裂缝和陷落柱等非均匀不连续地质体的勘探是一项具有挑战的任务。常规的处理技术,是在成像的剖面上利用相干性分析、反模糊滤波及构造导向滤波等方法提高地震分辨率。然而,由于该类方法属于图像驱动型处理领域,易产生假构造,并且其理论建立在地震带限子波与无限薄反射系数褶积(多维)模型,因此通过该方法无法突破瑞丽准则限制。本质上讲,常规的处理方法都不是基于地震波物理规律实现的,相反,绕射波是这些小尺度地质异常体地震响应的可靠载体,具有高分辨或者超分辨探测能力。当地震波在传播过程中遇到反射界面的突变点(如断层的断棱,地层尖灭点,不整合面突起点、裂缝等)时,该突变点会成为次震源,发出球面波。这种地震波在地震勘探中称为绕射波,在研究小尺度地质体上具有重要意义,根据绕射波处理方法也随之被大家广泛使用。
目前,大多数的绕射波技术仅仅关注于绕射波的分离与成像,而对于绕射波聚焦分析的研究很少,工业界中常用的速度分析方法是基于反射波理论建立的,不利于绕射波能量聚焦,绕射波聚焦质量差。现有技术中还没有提出一种能够准确的求取绕射波聚焦能量的方法,绕射波聚焦问题仍然存在。
发明内容
有鉴于此,本发明的目的在于提供一种绕射波聚焦能量的分析方法及装置,以缓解现有技术中在求取绕射波聚焦能量时,精确度较低的技术问题。
第一方面,本发明实施例提供了一种绕射波聚焦能量的分析方法,包括:
在采集到的地震共偏移距数据中提取共偏移距绕射波数据,其中,所述共偏移距绕射波数据为偏移距相同的绕射波数据;
根据所述地震共偏移距数据的地震波偏移速度文件和所述共偏移距绕射波数据,获取角道集地震数据;
从所述角道集地震数据中提取绕射波目标振幅值,并将所述提取得到的所述目标振幅值作为绕射波能量向量组中的绕射波能量向量,其中,所述目标振幅值的数量为多个;
计算所述绕射波能量向量组中的每个所述绕射波能量向量的方差值,得到多个方差值;
在所述多个方差值中查找最大方差值,并将所述最大方差值作为目标聚焦能量,其中,在通过所述目标聚焦能量进行绕射波聚焦成像时,成像结果的清晰度高于通过除所述目标聚焦能量之外的任意一个聚焦能量进行绕射波聚焦成像时,成像结果的清晰度。
结合第一方面,本发明实施例提供了第一方面的第一种可能的实施方式,其中,所述从所述角道集地震数据中提取绕射波目标振幅值,并将所述提取得到的所述目标振幅值作为绕射波能量向量组中的绕射波能量向量包括:
获取预先设置的多个绕射波能量聚焦参数,其中,所述多个绕射波能量聚焦参数中任意两个参数的参数值不相同;
根据所述多个绕射波能量聚焦参数中的每个绕射波能量聚焦参数确定绕射波角道集曲线方程,得到多个所述绕射波角道集曲线方程;
从所述角道集地震数据中提取与每个所述绕射波角道集曲线方程的空间和时间位置相对应的所述目标振幅值,并将所述目标振幅值作为所述绕射波能量向量组中的所述绕射波能量向量。
结合第一方面,本发明实施例提供了第一方面的第二种可能的实施方式,其中,所述绕射波角道集曲线方程为:其中,τ为对应于角度α的绕射波时间值,τ0为对应于角度α=0时的绕射波初始时间值,γ为所述绕射波能量聚焦参数。
结合第一方面,本发明实施例提供了第一方面的第三种可能的实施方式,其中,在所述多个方差值中查找最大方差值之后,所述方法还包括:
确定与所述最大方差值相对应的绕射波能量向量;
将用于计算所述相对应的绕射波能量向量的绕射波能量聚焦参数作为绕射波聚焦目标参数。
结合第一方面,本发明实施例提供了第一方面的第四种可能的实施方式,其中,所述在采集到的地震共偏移距数据中提取共偏移距绕射波数据包括:
获取所述地震共偏移距数据,其中,所述地震共偏移距数据包括:共偏移距反射波数据、共偏移距绕射波数据;
对所述地震共偏移距数据进行平面波分解,得到所述共偏移距反射波数据;
根据所述地震共偏移距数据和所述共偏移距反射波数据,计算得到所述共偏移距绕射波数据。
结合第一方面,本发明实施例提供了第一方面的第五种可能的实施方式,其中,所述根据所述地震共偏移距数据的地震波偏移速度文件和所述共偏移距绕射波数据,获取角道集地震数据包括:
获取所述地震波偏移速度文件;
根据所述地震波偏移速度文件和所述共偏移距绕射波数据,通过Kirchhoff角道集叠前时间偏移算法计算得到所述角道集地震数据。
第二方面,本发明实施例还提供一种绕射波聚焦能量的分析装置,所述装置包括:
第一提取模块,用于在采集到的地震共偏移距数据中提取共偏移距绕射波数据,其中,所述共偏移距绕射波数据为偏移距相同的绕射波数据;
获取模块,用于根据所述地震共偏移距数据的地震波偏移速度文件和所述共偏移距绕射波数据,获取角道集地震数据;
第二提取模块,用于从所述角道集地震数据中提取绕射波目标振幅值,并将所述提取得到的所述目标振幅值作为绕射波能量向量组中的绕射波能量向量,其中,所述目标振幅值的数量为多个;
计算模块,用于计算所述绕射波能量向量组中的每个所述绕射波能量向量的方差值,得到多个方差值;
查找模块,用于在所述多个方差值中查找最大方差值,并将所述最大方差值作为目标聚焦能量,其中,在通过所述目标聚焦能量进行绕射波聚焦成像时,成像结果的清晰度高于通过除所述目标聚焦能量之外的任意一个聚焦能量进行绕射波聚焦成像时,成像结果的清晰度。
结合第二方面,本发明实施例提供了第二方面的第一种可能的实施方式,其中,所述第二提取模块包括:
第一获取子模块,用于获取预先设置的多个绕射波能量聚焦参数,其中,所述多个绕射波能量聚焦参数中任意两个参数的参数值不相同;
确定子模块,用于根据所述多个绕射波能量聚焦参数中的每个绕射波能量聚焦参数确定绕射波角道集曲线方程,得到多个所述绕射波角道集曲线方程;
提取子模块,用于从所述角道集地震数据中提取与每个所述绕射波角道集曲线方程的空间和时间位置相对应的所述目标振幅值,并将所述目标振幅值作为所述绕射波能量向量组中的所述绕射波能量向量。
结合第二方面,本发明实施例提供了第二方面的第二种可能的实施方式,其中,所述绕射波角道集曲线方程为:其中,τ为对应于角度α的绕射波时间值,τ0为对应于角度α=0时的绕射波初始时间值,γ为所述绕射波能量聚焦参数。
结合第二方面,本发明实施例提供了第二方面的第三种可能的实施方式,其中,所述装置还包括:
确定模块,用于确定与所述最大方差值相对应的绕射波能量向量;
设定模块,用于将用于计算所述相对应的绕射波能量向量的绕射波能量聚焦参数作为绕射波聚焦目标参数。
本发明实施例带来了以下有益效果:本发明实施例提供的一种绕射波聚焦能量的分析方法及装置,该方法包括在采集到的地震共偏移距数据中提取共偏移距绕射波数据,其中,共偏移距绕射波数据为偏移距相同的绕射波数据;根据地震共偏移距数据的地震波偏移速度文件和共偏移距绕射波数据,获取角道集地震数据;从角道集地震数据中提取绕射波目标振幅值,并将提取得到的目标振幅值作为绕射波能量向量组中的绕射波能量向量,其中,目标振幅值的数量为多个;计算绕射波能量向量组中的每个绕射波能量向量的方差值,得到多个方差值;在多个方差值中查找最大方差值,并将最大方差值作为目标聚焦能量,其中,在通过目标聚焦能量进行绕射波聚焦成像时,成像结果的清晰度高于通过除目标聚焦能量之外的任意一个聚焦能量进行绕射波聚焦成像时,成像结果的清晰度;
与现有技术中采用的速度分析方法进行绕射波聚焦能量分析相比,其通过获取角道集地震数据并从获取得到的角道集地震数据中提取绕射波目标振幅值,同时将提取得到的目标振幅值作为绕射波能量向量组中的绕射波能量向量,最终,计算绕射波能量向量组中的每个绕射波能量向量的方差值,得到的最大方差值作为目标聚焦能量。现有技术中采用的速度分析方法是基于反射波理论建立的,用于绕射波聚焦能量分析中不够准确,分析过程中绕射波能量损耗严重,不利于绕射波能量聚焦。本发明根据角道集地震数据实现了角度域绕射波能量聚焦,从角道集地震数据中提取绕射波目标振幅值,并将目标振幅值作为绕射波能量向量组中的绕射波能量向量,然后求取每个绕射波能量向量的方差值,以得到最大方差值并将其作为目标聚焦能量,能够准确的求取绕射波聚焦能量,实现绕射波聚焦成像,达到精确定位地下空间小尺度非均匀地质体的目的,本发明在资源勘探和地质灾害防治中具有重要应用价值。缓解了现有技术中在求取绕射波聚焦能量时,精确度较差的技术问题。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种绕射波聚焦能量的分析方法的流程图;
图2为本发明实施例提供的提取共偏移距绕射波数据的流程图;
图3为本发明实施例提供的获取角道集地震数据的流程图;
图4为本发明实施例提供的从角道集地震数据中提取绕射波目标振幅值,并将提取得到的目标振幅值作为绕射波能量向量组中的绕射波能量向量的流程图;
图5为本发明实施例提供的确定绕射波聚焦目标参数的流程图;
图6为本发明实施例提供的一种绕射波聚焦能量的分析装置的结构示意图。
图标:
1-绕射波聚焦能量的分析装置;11-第一提取模块;12-获取模块;13-第二提取模块;14-计算模块;15-查找模块;16-确定模块;17-设定模块;131-第一获取子模块;132-确定子模块;133-提取子模块;111-第二获取子模块;112-平面波分解子模块;113-第一计算子模块;121-第三获取子模块;122-第二计算子模块。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前,现有技术中采用的速度分析方法是基于反射波理论建立的,绕射波能量损耗严重,不利于绕射波能量聚焦。基于此,本发明实施例提供的一种绕射波聚焦能量的分析方法及装置,可根据角道集地震数据实现角度域绕射波能量聚焦,从角道集地震数据中提取绕射波目标振幅值,并将目标振幅值作为绕射波能量向量组中的绕射波能量向量,然后求取每个绕射波能量向量的方差值,以得到最大方差值并将其作为目标聚焦能量,能够准确的求取绕射波聚焦能量,实现绕射波聚焦成像。缓解了现有技术中在求取绕射波聚焦能量时,精确度较差的技术问题。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种绕射波聚焦能量的分析方法进行详细介绍。
一种绕射波聚焦能量的分析方法,参考图1,该方法包括:
S101、在采集到的地震共偏移距数据中提取共偏移距绕射波数据,其中,共偏移距绕射波数据为偏移距相同的绕射波数据;
在本发明实施例中,地震共偏移距数据为偏移距相同的地震波数据,而地震波数据是通过地震观测系统采集得到的,然后地震波数据经过地震软件的处理得到了用于分析的地震共偏移距数据。在物理学中,把波绕过障碍物而传播的现象称为绕射,绕射时,波的路径发生了改变或弯曲,符合惠更斯原理。而在地震勘探工作中,把地震波在传播中遇到地层剧烈变化的地方,例如断层的断点、断棱,地层尖灭点,不均匀体、侵入体和地下礁的边缘等所引起的波为绕射波,所以绕射波是地质信息的载体。为了便于分析,常常在采集到的地震共偏移距数据中提取共偏移距绕射波数据,然后通过对共偏移距绕射波数据的准确分析能够精确的定位地下空间小尺度的非均匀地质体,减少灾害的发生。
S102、根据地震共偏移距数据的地震波偏移速度文件和共偏移距绕射波数据,获取角道集地震数据;
具体的,地震波偏移速度文件用于指示地震波在地质体内的传播速度;最终获取得到的角道集地震数据是不仅包含绕射波的振幅值,还包括角度、传播时间等参数,其中角度可以认为是空间位置,这样角道集地震数据就是包含绕射波的振幅值、空间和时间位置参数。
S103、从角道集地震数据中提取绕射波目标振幅值,并将提取得到的目标振幅值作为绕射波能量向量组中的绕射波能量向量,其中,目标振幅值的数量为多个;
具体的,在所有得到的角道集地震数据中,只有一部分的数据和绕射波角道集曲线方程的时间和空间(即角度)位置相对应,提取得到的绕射波目标振幅值即为和绕射波角道集曲线方程的时间和空间(即角度)位置相对应的角道集地震数据的绕射波振幅值,将得到的目标振幅值作为绕射波能量向量组中的绕射波能量向量,以供后续的分析计算。
S104、计算绕射波能量向量组中的每个绕射波能量向量的方差值,得到多个方差值;
具体的,方差值的计算公式为:其中,方差值φ用于判断绕射波聚焦能量,si为绕射波能量向量,N为绕射波振幅个数。
S105、在多个方差值中查找最大方差值,并将最大方差值作为目标聚焦能量,其中,在通过目标聚焦能量进行绕射波聚焦成像时,成像结果的清晰度高于通过除目标聚焦能量之外的任意一个聚焦能量进行绕射波聚焦成像时,成像结果的清晰度。
在本发明实施例中,在计算得到多个方差值后,在多个方差值中查找最大方差值。具体的,可以通过最大值查找算法进行查找,也可以通过其它方式进行查找,本发明实施例对其不做具体限制。将查找得到的最大方差值作为目标聚焦能量,其中,目标聚焦能量指的是最优聚焦能量。在通过目标聚焦能量进行绕射波聚焦成像时,成像结果的清晰度高于通过除目标聚焦能量之外的任意一个聚焦能量进行绕射波聚焦成像时,成像结果的清晰度。
本发明实施例提供的一种绕射波聚焦能量的分析方法,与现有技术中采用的速度分析方法进行绕射波聚焦能量分析相比,其通过获取角道集地震数据并从获取得到的角道集地震数据中提取绕射波目标振幅值,同时将提取得到的目标振幅值作为绕射波能量向量组中的绕射波能量向量,最终,计算绕射波能量向量组中的每个绕射波能量向量的方差值,得到的最大方差值作为目标聚焦能量。现有技术中采用的速度分析方法是基于反射波理论建立的,用于绕射波聚焦能量分析中不够准确,分析过程中绕射波能量损耗严重,不利于绕射波能量聚焦。本发明根据角道集地震数据实现了角度域绕射波能量聚焦,从角道集地震数据中提取绕射波目标振幅值,并将目标振幅值作为绕射波能量向量组中的绕射波能量向量,然后求取每个绕射波能量向量的方差值,以得到最大方差值并将其作为目标聚焦能量,能够准确的求取绕射波聚焦能量,实现绕射波聚焦成像,达到精确定位地下空间小尺度非均匀地质体的目的,缓解了现有技术中在求取绕射波聚焦能量时,精确度较差的技术问题,本发明在资源勘探和地质灾害防治中具有重要应用价值。
上述内容对绕射波聚焦能量的分析方法进行了整体简要的描述,下面对其中涉及的内容进行具体描述。参考图2,在采集到的地震共偏移距数据中提取共偏移距绕射波数据包括:
S201、获取地震共偏移距数据,其中,地震共偏移距数据包括:共偏移距反射波数据、共偏移距绕射波数据;
具体的,首先需要采集地震波数据。在进行采集时,需要在预设地面部署地震观测系统,包括:设置炮点位置以及设置检波点位置。在预设地面设置好炮点位置后,在该位置埋藏炸药;同时,在检波点位置放置检波器以收集地震波数据,其中,检波点位置为多个,同时对应的检波器也为多个。测试开始后,引爆炸药,产生的地震波经地面进入地质体内部,在到达不连续非均匀地质体后,地震波一部分经过反射以反射波的形式到达地面,另一部分经过绕射以绕射波的形式到达地面。置于地表面的检波器就可以接收通过地质体的地震波数据,这些地震波数据携带有地质信息。
然后,将上述检波器中接收到的地震波数据导入地震软件中,通过地震软件完成加载数据,以使炮点位置、检波点位置和采集到的地震波数据在坐标系中对应起来,也就是完成建立坐标的过程,通过地震软件对采集到的地震波数据去除噪声,静校正以及分选,最终可得到地震共偏移距数据,可用于后续的分析处理。
需要说明的是:上述炸药也可以为其它能够对地面产生地震波的物质,本发明实施例对其不做具体限制。
S202、对地震共偏移距数据进行平面波分解,得到共偏移距反射波数据;
具体的,平面波分解的过程为:d=C(σ)s,其中,d=[d1,d2,...dN]T为共偏移距反射波数据的矩阵表示形式,T为矩阵的转置运算,s为地震共偏移距数据矩阵表示形式,C为平面波分解算子,形式如下:
其中,Ι单位矩阵,P为低阶平面波预测算子,其标量形式如下:
其中,Zt,Zx分别为变量t,x的Z变换;σi表示地震倾角斜率。
最终得到共偏移距反射波数据d=[d1,d2,…dN]T。
S203、根据地震共偏移距数据和共偏移距反射波数据,计算得到共偏移距绕射波数据。
进一步的,由地震共偏移距数据(包括共偏移距反射波数据和共偏移距绕射波数据)减去步骤S202中得到的共偏移距反射波数据,计算得到共偏移距绕射波数据。
参考图3,根据地震共偏移距数据的地震波偏移速度文件和共偏移距绕射波数据,获取角道集地震数据包括:
S301、获取地震波偏移速度文件;
具体的,地震波偏移速度文件用于指示地震波在地质体中的传播速度。偏移速度文件是由地震波数据经过偏移速度分析处理得到的;其中,上述偏移速度文件,对应于地层速度参数,用于计算地震波传播旅行时(即时间),其是由对地震波数据进行偏移速度分析获得;而上述偏移速度分析是一种速度建模技术,该速度建模技术通过分析地震采集数据中的地震波聚焦性确定速度参数,即获得地震波的偏移速度文件。
S302、根据地震波偏移速度文件和共偏移距绕射波数据,通过Kirchhoff角道集叠前时间偏移算法计算得到角道集地震数据。
在本发明实施例中,Kirchhoff角道集叠前时间偏移为:
其中,射线由炮点位置s(x,y,z)激发,经地质体成像点到达检波点位置r(x,y,z)。分别为由炮点到成像点和由检波点到成像点的射线参数,方向分别由表示。θ,分别为出射角和方位角。θm,为假设反射界面法向量,式中,nx=(1,0,0),W(τ0,r)为振幅补偿因子。
具体的,由地震波偏移速度文件可以确定上述公式中的由共偏移距绕射波数据可以确定上述公式中的W(τ0,r)。这样得到的角道集地震数据为时间和角度(即空间)的关系方程,另外角道集地震数据还包括绕射波振幅值。
参考图4,从角道集地震数据中提取绕射波目标振幅值,并将提取得到的目标振幅值作为绕射波能量向量组中的绕射波能量向量包括:
S401、获取预先设置的多个绕射波能量聚焦参数,其中,多个绕射波能量聚焦参数中任意两个参数的参数值不相同;
具体的,绕射波能量聚焦参数可用γ表示,可以根据需要对γ的值自行设定,不同的绕射波能量聚焦参数可以构成不同的绕射波角道集曲线方程。进而,通过角道集地震数据和绕射波角道集曲线方程可以确定出多个绕射波目标振幅值。
S402、根据多个绕射波能量聚焦参数中的每个绕射波能量聚焦参数确定绕射波角道集曲线方程,得到多个绕射波角道集曲线方程;
具体的,绕射波角道集曲线方程为:其中,τ为对应于角度α的绕射波时间值,τ0为对应于角度α=0时的绕射波初始时间值,γ为绕射波能量聚焦参数;
这样,根据多个绕射波能量聚焦参数,就能确定出对应的多个绕射波角道集曲线方程。
S403、从角道集地震数据中提取与每个绕射波角道集曲线方程的空间和时间位置相对应的目标振幅值,并将目标振幅值作为绕射波能量向量组中的绕射波能量向量。
具体的,绕射波角道集曲线方程为角度(即空间位置)和时间的关系曲线方程,步骤S402中确定出多个绕射波能量聚焦参数对应的绕射波角道集曲线方程后,从角道集地震数据中提取与每个绕射波角道集曲线方程的空间和时间位置相对应的目标振幅值,即提取满足角道集曲线方程中角度(即空间位置)和时间关系的角道集地震数据的绕射波目标振幅值,并将目标振幅值作为绕射波能量向量组中的绕射波能量向量,这样绕射波能量向量组中的每个绕射波能量向量都有对应的绕射波能量聚焦参数。
参考图5,在多个方差值中查找最大方差值之后,方法还包括:
S501、确定与最大方差值相对应的绕射波能量向量;
具体的,在得到最大方差值后,可以同时确定与最大方差值相对应的绕射波能量向量组中的某个绕射波能量向量。
S502、将用于计算相对应的绕射波能量向量的绕射波能量聚焦参数作为绕射波聚焦目标参数。
具体的,在得到与最大方差值相对应的绕射波能量向量后,可以进一步得到与该绕射波能量向量所对应的绕射波能量聚焦参数,并将该绕射波能量聚焦参数作为绕射波聚焦目标参数。上述绕射波聚焦目标参数对应于绕射波目标聚焦能量。
本发明提供的一种绕射波聚焦能量的分析方法,提高了绕射波聚焦质量,从而可准确探测与油气储集有关的缝洞体和与地质灾害有关的断裂和陷落柱等不连续非均匀地质体,在地质勘探中,意义重大。
本发明实施例还提供了一种绕射波聚焦能量的分析装置1,参考图6,该装置包括:
第一提取模块11,用于在采集到的地震共偏移距数据中提取共偏移距绕射波数据,其中,共偏移距绕射波数据为偏移距相同的绕射波数据;
获取模块12,用于根据地震共偏移距数据的地震波偏移速度文件和共偏移距绕射波数据,获取角道集地震数据;
第二提取模块13,用于从角道集地震数据中提取绕射波目标振幅值,并将提取得到的目标振幅值作为绕射波能量向量组中的绕射波能量向量,其中,目标振幅值的数量为多个;
计算模块14,用于计算绕射波能量向量组中的每个绕射波能量向量的方差值,得到多个方差值;
查找模块15,用于在多个方差值中查找最大方差值,并将最大方差值作为目标聚焦能量,其中,在通过目标聚焦能量进行绕射波聚焦成像时,成像结果的清晰度高于通过除目标聚焦能量之外的任意一个聚焦能量进行绕射波聚焦成像时,成像结果的清晰度。
其中,第二提取模块13包括:
第一获取子模块131,用于获取预先设置的多个绕射波能量聚焦参数,其中,多个绕射波能量聚焦参数中任意两个参数的参数值不相同;
确定子模块132,用于根据多个绕射波能量聚焦参数中的每个绕射波能量聚焦参数确定绕射波角道集曲线方程,得到多个绕射波角道集曲线方程;
提取子模块133,用于从角道集地震数据中提取与每个绕射波角道集曲线方程的空间和时间位置相对应的目标振幅值,并将目标振幅值作为绕射波能量向量组中的绕射波能量向量。
绕射波角道集曲线方程为:其中,τ为对应于角度α的绕射波时间值,τ0为对应于角度α=0时的绕射波初始时间值,γ为绕射波能量聚焦参数。
绕射波聚焦能量的分析装置1还包括:
确定模块16,用于确定与最大方差值相对应的绕射波能量向量;
设定模块17,用于将用于计算相对应的绕射波能量向量的绕射波能量聚焦参数作为绕射波聚焦目标参数。
其中,第一提取模块11包括:
第二获取子模块111,用于获取地震共偏移距数据,其中,地震共偏移距数据包括:共偏移距反射波数据、共偏移距绕射波数据;
平面波分解子模块112,用于对地震共偏移距数据进行平面波分解,得到共偏移距反射波数据;
第一计算子模块113,用于根据地震共偏移距数据和共偏移距反射波数据,计算得到共偏移距绕射波数据。
进一步的,获取模块12包括:
第三获取子模块121,用于获取地震波偏移速度文件;
第二计算子模块122,用于根据地震波偏移速度文件和共偏移距绕射波数据,通过Kirchhoff角道集叠前时间偏移算法计算得到角道集地震数据。
本发明实施例所提供的绕射波聚焦能量的分析装置,包括存储了程序代码的计算机可读存储介质,所述程序代码包括的指令可用于执行前面方法实施例中所述的方法,具体实现可参见方法实施例,在此不再赘述。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的装置的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
另外,在本发明实施例的描述中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种绕射波聚焦能量的分析方法,其特征在于,包括:
在采集到的地震共偏移距数据中提取共偏移距绕射波数据,其中,所述共偏移距绕射波数据为偏移距相同的绕射波数据;
根据所述地震共偏移距数据的地震波偏移速度文件和所述共偏移距绕射波数据,获取角道集地震数据;
从所述角道集地震数据中提取绕射波目标振幅值,并将所述提取得到的所述目标振幅值作为绕射波能量向量组中的绕射波能量向量,其中,所述目标振幅值的数量为多个;
计算所述绕射波能量向量组中的每个所述绕射波能量向量的方差值,得到多个方差值;
在所述多个方差值中查找最大方差值,并将所述最大方差值作为目标聚焦能量,其中,在通过所述目标聚焦能量进行绕射波聚焦成像时,成像结果的清晰度高于通过除所述目标聚焦能量之外的任意一个聚焦能量进行绕射波聚焦成像时,成像结果的清晰度。
2.根据权利要求1所述的方法,其特征在于,所述从所述角道集地震数据中提取绕射波目标振幅值,并将所述提取得到的所述目标振幅值作为绕射波能量向量组中的绕射波能量向量包括:
获取预先设置的多个绕射波能量聚焦参数,其中,所述多个绕射波能量聚焦参数中任意两个参数的参数值不相同;
根据所述多个绕射波能量聚焦参数中的每个绕射波能量聚焦参数确定绕射波角道集曲线方程,得到多个所述绕射波角道集曲线方程;
从所述角道集地震数据中提取与每个所述绕射波角道集曲线方程的空间和时间位置相对应的所述目标振幅值,并将所述目标振幅值作为所述绕射波能量向量组中的所述绕射波能量向量。
3.根据权利要求2所述的方法,其特征在于,所述绕射波角道集曲线方程为:其中,τ为对应于角度α的绕射波时间值,τ0为对应于角度α=0时的绕射波初始时间值,γ为所述绕射波能量聚焦参数。
4.根据权利要求2所述的方法,其特征在于,在所述多个方差值中查找最大方差值之后,所述方法还包括:
确定与所述最大方差值相对应的绕射波能量向量;
将用于计算所述相对应的绕射波能量向量的绕射波能量聚焦参数作为绕射波聚焦目标参数。
5.根据权利要求1所述的方法,其特征在于,所述在采集到的地震共偏移距数据中提取共偏移距绕射波数据包括:
获取所述地震共偏移距数据,其中,所述地震共偏移距数据包括:共偏移距反射波数据、共偏移距绕射波数据;
对所述地震共偏移距数据进行平面波分解,得到所述共偏移距反射波数据;
根据所述地震共偏移距数据和所述共偏移距反射波数据,计算得到所述共偏移距绕射波数据。
6.根据权利要求1所述的方法,其特征在于,所述根据所述地震共偏移距数据的地震波偏移速度文件和所述共偏移距绕射波数据,获取角道集地震数据包括:
获取所述地震波偏移速度文件;
根据所述地震波偏移速度文件和所述共偏移距绕射波数据,通过Kirchhoff角道集叠前时间偏移算法计算得到所述角道集地震数据。
7.一种绕射波聚焦能量的分析装置,其特征在于,所述装置包括:
第一提取模块,用于在采集到的地震共偏移距数据中提取共偏移距绕射波数据,其中,所述共偏移距绕射波数据为偏移距相同的绕射波数据;
获取模块,用于根据所述地震共偏移距数据的地震波偏移速度文件和所述共偏移距绕射波数据,获取角道集地震数据;
第二提取模块,用于从所述角道集地震数据中提取绕射波目标振幅值,并将所述提取得到的所述目标振幅值作为绕射波能量向量组中的绕射波能量向量,其中,所述目标振幅值的数量为多个;
计算模块,用于计算所述绕射波能量向量组中的每个所述绕射波能量向量的方差值,得到多个方差值;
查找模块,用于在所述多个方差值中查找最大方差值,并将所述最大方差值作为目标聚焦能量,其中,在通过所述目标聚焦能量进行绕射波聚焦成像时,成像结果的清晰度高于通过除所述目标聚焦能量之外的任意一个聚焦能量进行绕射波聚焦成像时,成像结果的清晰度。
8.根据权利要求7所述的装置,其特征在于,所述第二提取模块包括:
第一获取子模块,用于获取预先设置的多个绕射波能量聚焦参数,其中,所述多个绕射波能量聚焦参数中任意两个参数的参数值不相同;
确定子模块,用于根据所述多个绕射波能量聚焦参数中的每个绕射波能量聚焦参数确定绕射波角道集曲线方程,得到多个所述绕射波角道集曲线方程;
提取子模块,用于从所述角道集地震数据中提取与每个所述绕射波角道集曲线方程的空间和时间位置相对应的所述目标振幅值,并将所述目标振幅值作为所述绕射波能量向量组中的所述绕射波能量向量。
9.根据权利要求8所述的装置,其特征在于,
所述绕射波角道集曲线方程为:其中,τ为对应于角度α的绕射波时间值,τ0为对应于角度α=0时的绕射波初始时间值,γ为所述绕射波能量聚焦参数。
10.根据权利要求8所述的装置,其特征在于,所述装置还包括:
确定模块,用于确定与所述最大方差值相对应的绕射波能量向量;
设定模块,用于将用于计算所述相对应的绕射波能量向量的绕射波能量聚焦参数作为绕射波聚焦目标参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610988823.9A CN106772592B (zh) | 2016-11-10 | 2016-11-10 | 绕射波聚焦能量的分析方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610988823.9A CN106772592B (zh) | 2016-11-10 | 2016-11-10 | 绕射波聚焦能量的分析方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772592A true CN106772592A (zh) | 2017-05-31 |
CN106772592B CN106772592B (zh) | 2018-08-07 |
Family
ID=58972807
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610988823.9A Active CN106772592B (zh) | 2016-11-10 | 2016-11-10 | 绕射波聚焦能量的分析方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772592B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107817523A (zh) * | 2017-10-30 | 2018-03-20 | 中国矿业大学(北京) | 绕射波偏移速度的分析方法及装置 |
CN108226999A (zh) * | 2018-01-19 | 2018-06-29 | 中国石油化工股份有限公司 | 碳酸盐岩小尺度缝洞体信息的处理方法 |
CN108614292A (zh) * | 2018-08-08 | 2018-10-02 | 中国海洋石油集团有限公司 | 一种基于绕射波数据的岩浆岩信息提取方法 |
CN108873067A (zh) * | 2018-09-26 | 2018-11-23 | 中国矿业大学(北京) | 绕射系数求解方法和装置 |
CN108919350A (zh) * | 2018-09-26 | 2018-11-30 | 中国矿业大学(北京) | 绕射波分离方法及装置 |
CN108919351A (zh) * | 2018-05-15 | 2018-11-30 | 中国海洋石油集团有限公司 | 基于逆时聚焦原理进行观测系统双向聚焦性的评价方法 |
CN111025383A (zh) * | 2019-11-21 | 2020-04-17 | 徐州工程学院 | 一种基于绕射横波定性判断隧道前方溶洞充水情况的方法 |
CN112255679A (zh) * | 2020-10-26 | 2021-01-22 | 中国石油天然气集团有限公司 | 地震资料绕射深度偏移处理方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101105537A (zh) * | 2006-07-12 | 2008-01-16 | 中国石油集团东方地球物理勘探有限责任公司 | 一种高精度的深度域叠前地震数据反演方法 |
US20110313667A1 (en) * | 2007-02-13 | 2011-12-22 | Levin Stewart A | Analyzing Borehole Seismic Data to Locate Subsurface Diffractors |
CN102520444A (zh) * | 2011-12-13 | 2012-06-27 | 中国科学院地质与地球物理研究所 | 一种叠后地震波中绕射波信息提取方法 |
CN104076395A (zh) * | 2014-04-16 | 2014-10-01 | 孙赞东 | 基于滤波组合的镜面能量提取与成像方法 |
CN104237940A (zh) * | 2014-09-29 | 2014-12-24 | 中国石油天然气股份有限公司 | 一种基于动力学特征的绕射波成像方法及装置 |
CN104730572A (zh) * | 2015-03-11 | 2015-06-24 | 中国科学院地质与地球物理研究所 | 一种基于l0半范数的绕射波成像方法及装置 |
WO2016032353A1 (ru) * | 2014-08-26 | 2016-03-03 | Закрытое Акционерное Общество "Технологии Обратных Задач" | Способ поиска залежей углеводородов, приуроченных к трещинно-кавернозным коллекторам |
CN105785439A (zh) * | 2016-02-01 | 2016-07-20 | 北京中科联衡科技有限公司 | 小尺度非均匀地质体空间分布位置的预测方法和装置 |
-
2016
- 2016-11-10 CN CN201610988823.9A patent/CN106772592B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101105537A (zh) * | 2006-07-12 | 2008-01-16 | 中国石油集团东方地球物理勘探有限责任公司 | 一种高精度的深度域叠前地震数据反演方法 |
US20110313667A1 (en) * | 2007-02-13 | 2011-12-22 | Levin Stewart A | Analyzing Borehole Seismic Data to Locate Subsurface Diffractors |
CN102520444A (zh) * | 2011-12-13 | 2012-06-27 | 中国科学院地质与地球物理研究所 | 一种叠后地震波中绕射波信息提取方法 |
CN104076395A (zh) * | 2014-04-16 | 2014-10-01 | 孙赞东 | 基于滤波组合的镜面能量提取与成像方法 |
WO2016032353A1 (ru) * | 2014-08-26 | 2016-03-03 | Закрытое Акционерное Общество "Технологии Обратных Задач" | Способ поиска залежей углеводородов, приуроченных к трещинно-кавернозным коллекторам |
CN104237940A (zh) * | 2014-09-29 | 2014-12-24 | 中国石油天然气股份有限公司 | 一种基于动力学特征的绕射波成像方法及装置 |
CN104730572A (zh) * | 2015-03-11 | 2015-06-24 | 中国科学院地质与地球物理研究所 | 一种基于l0半范数的绕射波成像方法及装置 |
CN105785439A (zh) * | 2016-02-01 | 2016-07-20 | 北京中科联衡科技有限公司 | 小尺度非均匀地质体空间分布位置的预测方法和装置 |
Non-Patent Citations (5)
Title |
---|
刘玉金 等: ""绕射波叠前时间偏移速度分析及成像"", 《地球物理学进展》 * |
朱生旺 等: ""局部倾角滤波和预测反演联合分离绕射波"", 《地球物理学报》 * |
袁茂林 等: ""构造倾角约束的有效能量叠加成像方法及应用"", 《石油地球物理勘探》 * |
赵惊涛 等: ""基于绕射波的储层预测方法及其应用"", 《石油学报》 * |
黄建平 等: ""基于PWD的绕射波波场分离成像方法综述"", 《地球物理学进展》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107817523A (zh) * | 2017-10-30 | 2018-03-20 | 中国矿业大学(北京) | 绕射波偏移速度的分析方法及装置 |
CN108226999A (zh) * | 2018-01-19 | 2018-06-29 | 中国石油化工股份有限公司 | 碳酸盐岩小尺度缝洞体信息的处理方法 |
CN108919351A (zh) * | 2018-05-15 | 2018-11-30 | 中国海洋石油集团有限公司 | 基于逆时聚焦原理进行观测系统双向聚焦性的评价方法 |
CN108614292A (zh) * | 2018-08-08 | 2018-10-02 | 中国海洋石油集团有限公司 | 一种基于绕射波数据的岩浆岩信息提取方法 |
CN108873067A (zh) * | 2018-09-26 | 2018-11-23 | 中国矿业大学(北京) | 绕射系数求解方法和装置 |
CN108919350A (zh) * | 2018-09-26 | 2018-11-30 | 中国矿业大学(北京) | 绕射波分离方法及装置 |
CN111025383A (zh) * | 2019-11-21 | 2020-04-17 | 徐州工程学院 | 一种基于绕射横波定性判断隧道前方溶洞充水情况的方法 |
CN111025383B (zh) * | 2019-11-21 | 2021-09-24 | 徐州工程学院 | 一种基于绕射横波定性判断隧道前方溶洞充水情况的方法 |
CN112255679A (zh) * | 2020-10-26 | 2021-01-22 | 中国石油天然气集团有限公司 | 地震资料绕射深度偏移处理方法及装置 |
CN112255679B (zh) * | 2020-10-26 | 2023-09-26 | 中国石油天然气集团有限公司 | 地震资料绕射深度偏移处理方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106772592B (zh) | 2018-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772592A (zh) | 绕射波聚焦能量的分析方法及装置 | |
CN104237940B (zh) | 一种基于动力学特征的绕射波成像方法及装置 | |
CN107817523B (zh) | 绕射波偏移速度的分析方法及装置 | |
CN105607121A (zh) | 一种煤炭陷落柱识别方法和装置 | |
CN104730571A (zh) | 一种利用绕射再聚焦识别小尺度地质体的方法与装置 | |
CN101900833B (zh) | 一种地震散射p-p波成像速度分析方法 | |
CN102937721A (zh) | 利用初至波走时的有限频层析成像方法 | |
CN105425292A (zh) | 一种油气预测方法及装置 | |
CN106772593B (zh) | 绕射波的成像方法及装置 | |
CN101923175B (zh) | 波动方程偏移直接产生角道集方法 | |
US20190361138A1 (en) | Methods and systems for automated sonic imaging | |
CN112305591B (zh) | 隧道超前地质预报方法、计算机可读存储介质 | |
CN102914799A (zh) | 非等效体波场正演模拟方法及装置 | |
US11125900B2 (en) | Estimating an earth response | |
CN106990436A (zh) | 陷落柱的识别方法和装置 | |
CN103513279B (zh) | 一种基于地震波波动方程的照明分析计算方法及计算装置 | |
CN113703044B (zh) | 古河道宽度的校正方法、装置、电子设备及存储介质 | |
CN105093300A (zh) | 一种地质体边界识别方法及装置 | |
CN112444859A (zh) | 协同异变蚂蚁体的页岩储层断裂识别方法及系统 | |
CN107422375B (zh) | 海底节点的共反射点道集的确定方法和装置 | |
CN105277981B (zh) | 基于波场延拓补偿的非一致性时移地震面元匹配方法 | |
CN108051858A (zh) | 优化的地震勘探方法及装置 | |
CN103020999A (zh) | 一种形变区域识别的方法及系统 | |
CN108897046B (zh) | 远探测声波逆时偏移成像方法及装置 | |
CN107085236B (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 |