CN112415591A - 绕射波的成像方法、装置、电子设备及存储介质 - Google Patents
绕射波的成像方法、装置、电子设备及存储介质 Download PDFInfo
- Publication number
- CN112415591A CN112415591A CN202011193675.4A CN202011193675A CN112415591A CN 112415591 A CN112415591 A CN 112415591A CN 202011193675 A CN202011193675 A CN 202011193675A CN 112415591 A CN112415591 A CN 112415591A
- Authority
- CN
- China
- Prior art keywords
- diffracted
- gather
- wave
- ovt
- waves
- 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 69
- 238000003860 storage Methods 0.000 title claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 76
- 238000013508 migration Methods 0.000 claims abstract description 47
- 230000005012 migration Effects 0.000 claims abstract description 47
- 238000012937 correction Methods 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000001914 filtration Methods 0.000 claims abstract description 30
- 238000004458 analytical method Methods 0.000 claims abstract description 14
- 238000007781 pre-processing Methods 0.000 claims abstract description 14
- 238000004590 computer program Methods 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 7
- 230000006870 function Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 13
- 238000004891 communication Methods 0.000 description 7
- 239000011324 bead Substances 0.000 description 6
- 238000000926 separation method Methods 0.000 description 5
- 239000000872 buffer Substances 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000001028 reflection method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation 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
本发明公开了一种绕射波的成像方法、装置、电子设备及存储介质,其中,方法包括:获取地震资料,并对地震资料进行预处理,地震资料包括:反射波和绕射波;根据预处理后的地震资料生成共中心点(CMP)道集;对预处理后的地震资料进行反射波速度分析操作,以得到反射波叠加速度和叠前时间偏移速度场;通过反射波叠加速度对CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成炮检距矢量片(OVT)道集;对OVT道集进行频率波数域滤波处理,以得到绕射波,并对绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。通过本发明,可以较好地分离绕射波,从而提高绕射波成像的分辨率,得到较高分辨率的绕射波成像剖面。
Description
技术领域
本发明涉及油气田勘探技术领域,具体涉及一种绕射波的成像方法、装置、电子设备及存储介质。
背景技术
现有反射法在解决石油、天然气、煤田等资源勘探与开发,以及近地表工程等问题中做出了重大贡献,为社会带来了巨大的经济财富。但对于常规反射地震勘探来说,随着地质问题的复杂化程度的提高,其弊端和问题也越来越多,反射勘探技术主要利用反射波等来研究地下地质情况,而实际地震波的类型是十分丰富的,反射波仅仅携带了一部分地质信息,只利用反射波反演地质属性特征得到的地质信息必将是不完整的。
绕射波是非水平弹性界面或非均匀体的存在而产生的,它们与地下某些地质现象存在必然联系,是检验断层、盐丘边界及碳酸盐岩缝洞储集体等小尺度地质体的有效波,能够用于高分辨率探测,它在识别裂缝、研究非均匀性问题等方面比反射波更为有效。绕射波一方面是确定地下异常体的重要依据,另一方面它也使得地震记录复杂化,一定程度上干扰反射波的成像。因此,十分有必要单独利用绕射波获取小尺度异常地质体信息。
在常规地震资料处理中,绕射波常常会因采用了不当的滤波器而被压制掉,使得地下小异常体无法精确地判断出来,降低了成像的分辨率。因此,要想有效地利用绕射波信息解决地质问题,首要任务是将绕射波完整地分离出来。目前,国内外对绕射波的分离方法主要依据反射波和绕射波的走时差异,但构造异常点位置的反射波与绕射波常常相切,仅根据走时很难保证绕射波的能量;并且叠前数据中绕射波能量分布范围大,波场干涉严重,分离的难度很大。当介质速度变化强烈时,绕射波的形态更为复杂,无法很好地将它们分离。
发明内容
有鉴于此,本发明提供一种绕射波的成像方法、装置、电子设备及存储介质,以解决上述提及的至少一个问题。
根据本发明的第一方面,提供一种绕射波的成像方法,所述方法包括:
获取地震资料,并对所述地震资料进行预处理,所述地震资料包括:反射波和绕射波;
根据所述预处理后的地震资料生成共中心点CMP道集;
对预处理后的地震资料进行反射波速度分析操作,以得到反射波叠加速度和叠前时间偏移速度场;
通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成炮检距矢量片OVT道集;
对所述OVT道集进行频率波数域滤波处理,以得到绕射波,并对所述绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
根据本发明的第二方面,提供一种绕射波的成像装置,所述装置包括:
数据获取单元,用于获取地震资料,所述地震资料包括:反射波和绕射波;
预处理单元,用于对所述地震资料进行预处理;
CMP道集生成单元,用于根据所述预处理后的地震资料生成共中心点CMP道集;
反射波速度分析单元,用于对预处理后的地震资料进行反射波速度分析操作,以得到反射波叠加速度和叠前时间偏移速度场;
OVT道集生成单元,用于通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成OVT道集;
绕射波获得单元,用于对所述OVT道集进行频率波数域滤波处理,以得到绕射波;
绕射波成像单元,用于对所述绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
根据本发明的第三方面,提供一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现上述方法的步骤。
根据本发明的第四方面,提供一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现上述方法的步骤。
由上述技术方案可知,通过对获取的地震数据进行预处理后生成CMP道集,并对预处理后的地震资料进行反射波速度分析操作,得到反射波叠加速度和叠前时间偏移速度场,之后通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成OVT道集,之后对OVT道集进行频率波数域滤波处理,得到绕射波,并对绕射波进行叠前深度偏移处理,以生成绕射波成像剖面,相比于现有技术,本技术方案可以较好地分离绕射波,从而提高绕射波成像的分辨率,得到较高分辨率的绕射波成像剖面。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是根据本发明实施例的绕射波成像方法的流程图;
图2(a)是OVT域频率波数域滤波前的截图;
图2(b)是根据本发明实施例的OVT域频率波数域滤波后的截图;
图3(a)是任意线原始叠前深度偏移后的截图;
图3(b)是根据本发明实施例的绕射波叠前深度偏移后的截图;
图4(a)是原始储层预测平面图;
图4(b)是根据本发明实施例的绕射波储层预测平面图;
图5是根据本发明实施例的绕射波成像装置的结构框图;
图6为本发明实施例的电子设备600的系统构成的示意框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
针对现有技术中,无法较好地将绕射波从地震数据中分离出来,降低了绕射波成像的分辨率,本发明实施例提供一种绕射波的成像方案,该方案可以将绕射波从地震数据中较好地分离出来,从而提高绕射波成像的分辨率。以下结合附图来详细描述本发明实施例。
图1是根据本发明实施例的绕射波成像方法的流程图,如图1所示,该流程包括:
步骤101,获取地震资料,并对所述地震资料进行预处理,所述地震资料包括:反射波和绕射波。
这里的预处理可以是抑制噪声等操作。
步骤102,根据所述预处理后的地震资料生成共中心点(Common Middle Point,CMP)道集。
步骤103,对预处理后的地震资料进行反射波速度分析操作,以得到反射波叠加速度和叠前时间偏移速度场。
步骤104,通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成炮检距矢量片(Offset Vector Tile,OVT)道集。
具体地,可以先通过所述反射波叠加速度对所述CMP道集进行动校正处理,之后对动校正处理后的CMP道集通过分选方位角和炮检距方式生成多个炮检距矢量片OVT,以生成所述OVT道集。
步骤105,对所述OVT道集进行频率波数域滤波处理,以得到绕射波,并对所述绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
具体地,对所述OVT道集进行频率波数域滤波处理,以对所述反射波进行压制处理,保留所述绕射波。之后,通过所述反射波叠加速度对所述绕射波进行反动校正处理,并根据所述叠前时间偏移速度场对反动校正处理后的绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
上述叠前深度偏移是实现地质构造空间归位的一项处理技术,具体实现可以参见现有技术,本发明对此不作限制。
通过对获取的地震数据进行预处理后生成CMP道集,并对预处理后的地震资料进行反射波速度分析操作,得到反射波叠加速度和叠前时间偏移速度场,之后通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成OVT道集,之后对OVT道集进行频率波数域滤波处理,得到绕射波,并对绕射波进行叠前深度偏移处理,以生成绕射波成像剖面,相比于现有技术,本发明实施例可以较好地分离绕射波,从而提高绕射波成像的分辨率,得到较高分辨率的绕射波成像剖面。
为了更好地理解本发明,以下结合实例来详细描述本发明实施例的绕射波成像流程,该流程包括如下步骤1-3:
步骤1、建立反射波叠加速度模型和叠前时间偏移速度场。
具体地,在地震资料常规处理的基础上,进行反射波速度分析。检查输出的CMP道集和CRP(Common Reflection Point,共反射点道集)道集中的有效同相轴是否被拉平判断速度是否准确,形成反射波叠加速度模型和叠前时间偏移速度场。
步骤2、通过频率波数域滤波分离绕射波。
步骤1):CMP道集动校正。利用反射波叠加速度对CMP道集进行动校正,通过动校正后可以将反射波校平,绕射波由于校正量不足,仍然是曲线。
步骤2):OVT道集选排。将动校后的CMP道集通过分选方位角和炮检距的方式将相邻检波点和相邻炮点之间两两组合,形成闭合的一个矩形区域,该小矩形区域即是一个炮检距矢量片,所有的炮检距矢量片构成OVT道集。
步骤3):OVT道集频率波数域滤波。滤波窗内绕射波信号与反射波信号视速度不同,能量集中方向不同,因而可以通过压制强能量反射波,使得绕射波信号得以保留。
这里的频率波数域滤波,是指对于地震数据中的一个共炮点道集或一个CMP道集,可以对它进行二维傅里叶变换,即把反射时间t和道位置x表示的函数f(x,t)变换为以频率f和空间波数k表示的函数,在频率-波数域平面上为一条通过原点的直线,其斜率值f/k是一个常数,表示在时间-空间域上相应同相轴的视速度,这里视速度可表示为v*=f/k。倾角越大,在频率-波数域中此直线越靠近波数轴。倾角分量则在频率轴上,倾角相当于零波数,示波的视速度为无穷。有效信号和干扰信号在频率和视速度上是可以分开的,为利用频率-域滤波压制各种干扰提供了条件。对于地震剖面上不同倾角的同相轴,频率-波数域上表现为不同斜率的直线。
3、绕射波叠前深度偏移成像。
步骤1):将分离的绕射波利用叠加速度进行反动校正。
步骤2):在OVT域进行叠前深度偏移,即在深度域对每个偏移距矢量片先用已建立的叠前时间偏移速度场进行叠前深度偏移,之后从每个偏移后的偏移距矢量片抽取共反射点道集。OVT域叠前深度偏移后得到包含所有方位的CIP(Common Image Point,共成像点)道集,对这一CIP道集统一进行RMO(剩余时差)拾取、统一射线追踪建立方程组,最后联立求解方程组,进入速度更新循环。
步骤3):通过Poisson算法对OVT偏移的CIP道集进行统一剩余时差拾取,然后将这些深度剩余时差输入层析反演,反演时被转换为旅行时差得到速度更新量。进一步优化速度模型,调整叠前成像参数,直到获取高分辨率绕射波成像剖面。
由以上面描述可知,通过在CMP道集动校正后分选到OVT道集,在OVT道集分离绕射波并进行成像,提高了地震资料横向分辨率。
在本发明实施例中,绕射波成像的核心是准确分离绕射波,在OVT域动校正后直接利用反射波与绕射波的相干差异,对全波场数据进行滤波处理,频率波数域滤波的滤波窗内绕射波信号与反射波信号视速度不同,能量集中方向不同,通过压制强能量反射波,使得绕射波信号得以保留,得到了较为完整的绕射波波场。采用OVT域叠前深度偏移方法对绕射波进行成像处理,可以保留更多叠前地震信息,使得绕射波准确成像,获取更加丰富的构造内幕信息。
以下以果勒地区奥陶系鹰2段为例,来描述本发明实施例的绕射波分离成像流程,具体流程及结果如下:
1)OVT域道集选排,用叠前偏移速度对OVT域道集进行动校正后,进行频率波数域滤波分离绕射波,具体参见图2(a)和图2(b),其中,图2(a)是OVT域频率波数域滤波前的截图,图2(b)是OVT域频率波数域滤波后的截图。
2)OVT域叠前深度偏移,得到包含所有方位的CIP道集,对这一CIP道集统一进行RMO拾取、统一射线追踪建立方程组,最后联立求解方程组,进入速度更新循环。通过Poisson算法对OVT偏移的CIP道集进行统一剩余时差拾取,然后将这些深度剩余时差输入层析反演,反演时被转换为旅行时差得到速度更新量。进一步优化速度模型,调整叠前成像参数,直到获取高分辨率绕射波成像剖面。具体参见图3(a)和图3(b),其中,图3(a)是任意线原始叠前深度偏移后的截图,图3(b)是绕射波叠前深度偏移后的截图。
3)奥陶系鹰2段储层预测平面图对比分析,绕射波偏移资料串珠形态清晰,数量增加,绕射波分离成像技术效果明显,具体参见图4(a)和图4(b),其中,图4(a)是原始储层预测平面图,图4(b)是绕射波储层预测平面图。
塔里木盆地塔北、塔中地区奥陶系碳酸盐岩受多期岩溶叠加作用及断裂作用影响,缝洞型储层特别发育,是塔里木油田勘探开发的重要领域。缝洞体规模达到地震可识别的范围时,在地震叠后数据体上常表现为绕射波,偏移成像后呈“串珠”状,它是缝洞体中裂缝和溶洞的综合响应。“串珠”的位置基本对应了储层发育位置。因此,“串珠”的准确归位是实现碳酸盐岩储层准确预测的重要基础。
通过本发明实施例,奥陶系鹰2段绕射波偏移资料串珠形态清晰,数量增加,对断裂的刻画更加清楚,为储层预测和油藏开发提供了有力支持。
基于相似的发明构思,本发明实施例还提供一种绕射波成像装置,该装置优选地用于实现上述方法实施例中的流程。
图5是该绕射波成像装置的结构框图,如图5所示,该装置包括:数据获取单元51、预处理单元52、CMP道集生成单元53、反射波速度分析单元54、OVT道集生成单元55、绕射波获得单元56和绕射波成像单元57,其中:
数据获取单元51,用于获取地震资料,所述地震资料包括:反射波和绕射波;
预处理单元52,用于对所述地震资料进行预处理;
CMP道集生成单元53,用于根据所述预处理后的地震资料生成CMP道集;
反射波速度分析单元54,用于对预处理后的地震资料进行反射波速度分析操作,以得到反射波叠加速度和叠前时间偏移速度场;
OVT道集生成单元55,用于通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成OVT道集;
绕射波获得单元56,用于对所述OVT道集进行频率波数域滤波处理,以得到绕射波;
绕射波成像单元57,用于对所述绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
通过CMP道集生成单元53对获取的地震数据进行预处理后生成CMP道集,反射波速度分析单元54对预处理后的地震资料进行反射波速度分析操作,得到反射波叠加速度和叠前时间偏移速度场,之后OVT道集生成单元55通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成OVT道集,绕射波获得单元56对OVT道集进行频率波数域滤波处理,得到绕射波,绕射波成像单元57对绕射波进行叠前深度偏移处理,以生成绕射波成像剖面,相比于现有技术,本发明实施例可以较好地分离绕射波,从而提高绕射波成像的分辨率,得到较高分辨率的绕射波成像剖面。
具体地,上述OVT道集生成单元55包括:动校正处理模块和OVT道集生成模块,其中:动校正处理模块,用于通过所述反射波叠加速度对所述CMP道集进行动校正处理;OVT道集生成模块,用于对动校正处理后的CMP道集通过分选方位角和炮检距方式生成多个炮检距矢量片OVT,以生成所述OVT道集。
上述绕射波获得单元56具体用于:对所述OVT道集进行频率波数域滤波处理,以对所述反射波进行压制处理,保留所述绕射波。
上述绕射波成像单元57包括:反动校正处理模块和绕射波成像模块,其中:反动校正处理模块,用于通过所述反射波叠加速度对所述绕射波进行反动校正处理;绕射波成像模块,用于根据所述叠前时间偏移速度场对反动校正处理后的绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
上述各单元、各模块的具体执行过程,可以参见上述方法实施例中的描述,此处不再赘述。
在实际操作中,上述各单元、各模块可以组合设置、也可以单一设置,本发明不限于此。
本实施例还提供一种电子设备,该电子设备可以是台式计算机、平板电脑及移动终端等,本实施例不限于此。在本实施例中,该电子设备可以参照上述方法实施例进行实施及绕射波成像装置的实施例进行实施,其内容被合并于此,重复之处不再赘述。
图6为本发明实施例的电子设备600的系统构成的示意框图。如图6所示,该电子设备600可以包括中央处理器100和存储器140;存储器140耦合到中央处理器100。值得注意的是,该图是示例性的;还可以使用其他类型的结构,来补充或代替该结构,以实现电信功能或其他功能。
一实施例中,绕射波成像功能可以被集成到中央处理器100中。其中,中央处理器100可以被配置为进行如下控制:
获取地震资料,并对所述地震资料进行预处理,所述地震资料包括:反射波和绕射波;
根据所述预处理后的地震资料生成共中心点CMP道集;
对预处理后的地震资料进行反射波速度分析操作,以得到反射波叠加速度和叠前时间偏移速度场;
通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成炮检距矢量片OVT道集;
对所述OVT道集进行频率波数域滤波处理,以得到绕射波,并对所述绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
从上述描述可知,本申请实施例提供的电子设备,通过对获取的地震数据进行预处理后生成CMP道集,并对预处理后的地震资料进行反射波速度分析操作,得到反射波叠加速度和叠前时间偏移速度场,之后通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成OVT道集,之后对OVT道集进行频率波数域滤波处理,得到绕射波,并对绕射波进行叠前深度偏移处理,以生成绕射波成像剖面,相比于现有技术,本发明实施例可以提高绕射波成像的分辨率,得到较高分辨率的绕射波成像剖面。
在另一个实施方式中,绕射波成像装置可以与中央处理器100分开配置,例如可以将绕射波成像装置配置为与中央处理器100连接的芯片,通过中央处理器的控制来实现绕射波成像功能。
如图6所示,该电子设备600还可以包括:通信模块110、输入单元120、音频处理单元130、显示器160、电源170。值得注意的是,电子设备600也并不是必须要包括图6中所示的所有部件;此外,电子设备600还可以包括图6中没有示出的部件,可以参考现有技术。
如图6所示,中央处理器100有时也称为控制器或操作控件,可以包括微处理器或其他处理器装置和/或逻辑装置,该中央处理器100接收输入并控制电子设备600的各个部件的操作。
其中,存储器140,例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息,此外还可存储执行有关信息的程序。并且中央处理器100可执行该存储器140存储的该程序,以实现信息存储或处理等。
输入单元120向中央处理器100提供输入。该输入单元120例如为按键或触摸输入装置。电源170用于向电子设备600提供电力。显示器160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器,但并不限于此。
该存储器140可以是固态存储器,例如,只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器,其即使在断电时也保存信息,可被选择性地擦除且设有更多数据,该存储器的示例有时被称为EPROM等。存储器140还可以是某种其它类型的装置。存储器140包括缓冲存储器141(有时被称为缓冲器)。存储器140可以包括应用/功能存储部142,该应用/功能存储部142用于存储应用程序和功能程序或用于通过中央处理器100执行电子设备600的操作的流程。
存储器140还可以包括数据存储部143,该数据存储部143用于存储数据,例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器140的驱动程序存储部144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块110即为经由天线111发送和接收信号的发送机/接收机110。通信模块(发送机/接收机)110耦合到中央处理器100,以提供输入信号和接收输出信号,这可以和常规移动通信终端的情况相同。
基于不同的通信技术,在同一电子设备中,可以设置有多个通信模块110,如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)110还经由音频处理器130耦合到扬声器131和麦克风132,以经由扬声器131提供音频输出,并接收来自麦克风132的音频输入,从而实现通常的电信功能。音频处理器130可以包括任何合适的缓冲器、解码器、放大器等。另外,音频处理器130还耦合到中央处理器100,从而使得可以通过麦克风132能够在本机上录音,且使得可以通过扬声器131来播放本机上存储的声音。
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时以实现上述绕射波成像方法的步骤。
综上所述,本发明实施例充分利用地震绕射波真实的传播规律反演地质属性,针对绕射波自身的规律和特征,提供了一种针对绕射波的分离及成像方案,该方案在地震资料处理中提取绕射波并单独成像,获得比传统反射波法更加丰富的异常地质信息,利用绕射波携带的地质信息,可以提高地质异常体成像精度,进而提高储层钻遇率,为油气田高效开发提供有利技术支撑。
以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的,因此权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外,由于本领域的技术人员容易想到很多修改和改变,因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作,而是可以涵盖落入其范围内的所有合适修改和等同物。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种绕射波的成像方法,其特征在于,所述方法包括:
获取地震资料,并对所述地震资料进行预处理,所述地震资料包括:反射波和绕射波;
根据所述预处理后的地震资料生成共中心点CMP道集;
对预处理后的地震资料进行反射波速度分析操作,以得到反射波叠加速度和叠前时间偏移速度场;
通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成炮检距矢量片OVT道集;
对所述OVT道集进行频率波数域滤波处理,以得到绕射波,并对所述绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
2.根据权利要求1所述的方法,其特征在于,通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成炮检距矢量片OVT道集包括:
通过所述反射波叠加速度对所述CMP道集进行动校正处理;
对动校正处理后的CMP道集通过分选方位角和炮检距方式生成多个炮检距矢量片OVT,以生成所述OVT道集。
3.根据权利要求1所述的方法,其特征在于,对所述OVT道集进行频率波数域滤波处理,以得到绕射波包括:
对所述OVT道集进行频率波数域滤波处理,以对所述反射波进行压制处理,保留所述绕射波。
4.根据权利要求1所述的方法,其特征在于,对所述绕射波进行叠前深度偏移处理,以生成绕射波成像剖面包括:
通过所述反射波叠加速度对所述绕射波进行反动校正处理;
根据所述叠前时间偏移速度场对反动校正处理后的绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
5.一种绕射波的成像装置,其特征在于,所述装置包括:
数据获取单元,用于获取地震资料,所述地震资料包括:反射波和绕射波;
预处理单元,用于对所述地震资料进行预处理;
CMP道集生成单元,用于根据所述预处理后的地震资料生成共中心点CMP道集;
反射波速度分析单元,用于对预处理后的地震资料进行反射波速度分析操作,以得到反射波叠加速度和叠前时间偏移速度场;
OVT道集生成单元,用于通过所述反射波叠加速度对所述CMP道集进行动校正处理、并通过分选方位角和炮检距方式生成OVT道集;
绕射波获得单元,用于对所述OVT道集进行频率波数域滤波处理,以得到绕射波;
绕射波成像单元,用于对所述绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
6.根据权利要求5所述的装置,其特征在于,所述OVT道集生成单元包括:
动校正处理模块,用于通过所述反射波叠加速度对所述CMP道集进行动校正处理;
OVT道集生成模块,用于对动校正处理后的CMP道集通过分选方位角和炮检距方式生成多个炮检距矢量片OVT,以生成所述OVT道集。
7.根据权利要求5所述的装置,其特征在于,所述绕射波获得单元具体用于:
对所述OVT道集进行频率波数域滤波处理,以对所述反射波进行压制处理,保留所述绕射波。
8.根据权利要求5所述的装置,其特征在于,所述绕射波成像单元包括:
反动校正处理模块,用于通过所述反射波叠加速度对所述绕射波进行反动校正处理;
绕射波成像模块,用于根据所述叠前时间偏移速度场对反动校正处理后的绕射波进行叠前深度偏移处理,以生成绕射波成像剖面。
9.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现权利要求1至4中任一项所述方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至4中任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011193675.4A CN112415591B (zh) | 2020-10-30 | 2020-10-30 | 绕射波的成像方法、装置、电子设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011193675.4A CN112415591B (zh) | 2020-10-30 | 2020-10-30 | 绕射波的成像方法、装置、电子设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112415591A true CN112415591A (zh) | 2021-02-26 |
CN112415591B CN112415591B (zh) | 2024-08-27 |
Family
ID=74828245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011193675.4A Active CN112415591B (zh) | 2020-10-30 | 2020-10-30 | 绕射波的成像方法、装置、电子设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112415591B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113640872A (zh) * | 2021-08-12 | 2021-11-12 | 中国矿业大学(北京) | 绕射波分离方法、装置和电子设备 |
CN114325821A (zh) * | 2021-12-23 | 2022-04-12 | 西安交通大学 | 基于3d-snacnn网络的叠前地震资料中强散射噪声压制方法及系统 |
CN115616660A (zh) * | 2022-12-20 | 2023-01-17 | 广州海洋地质调查局三亚南海地质研究所 | 利用绕射波监测海域碳封存项目co2泄露情况的方法及装置 |
CN117233844A (zh) * | 2023-09-19 | 2023-12-15 | 中国地质科学院地球物理地球化学勘查研究所 | 一种断层发育区地震高精度成像的数据处理方法及系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100131205A1 (en) * | 2007-03-12 | 2010-05-27 | Geomage (2003) Ltd | Method for identifying and analyzing faults/fractures using reflected and diffracted waves |
CN102778693A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油化工股份有限公司 | 一种基于反射波层拉平提取并消除的绕射波分离处理方法 |
CN103675897A (zh) * | 2012-08-30 | 2014-03-26 | 中国石油化工股份有限公司 | 一种地震绕射波分离成像方法 |
CN105785439A (zh) * | 2016-02-01 | 2016-07-20 | 北京中科联衡科技有限公司 | 小尺度非均匀地质体空间分布位置的预测方法和装置 |
CN107918147A (zh) * | 2017-11-20 | 2018-04-17 | 中国矿业大学(北京) | 绕射波成像方法和装置 |
CN108226999A (zh) * | 2018-01-19 | 2018-06-29 | 中国石油化工股份有限公司 | 碳酸盐岩小尺度缝洞体信息的处理方法 |
US20180246241A1 (en) * | 2016-02-02 | 2018-08-30 | China University Of Mining & Technology, Beijing | Method and apparatus for identifying collapsed coal column |
US20180292553A1 (en) * | 2017-01-10 | 2018-10-11 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method and Apparatus for Separating Seismic Diffracted Wave |
CN109507732A (zh) * | 2018-11-07 | 2019-03-22 | 中石化石油工程技术服务有限公司 | 一种基于成像道集的绕射波分离成像方法 |
CN111474580A (zh) * | 2020-05-13 | 2020-07-31 | 中国海洋石油集团有限公司 | 一种基于炮检距矢量片的方位角道集提取方法和系统 |
-
2020
- 2020-10-30 CN CN202011193675.4A patent/CN112415591B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100131205A1 (en) * | 2007-03-12 | 2010-05-27 | Geomage (2003) Ltd | Method for identifying and analyzing faults/fractures using reflected and diffracted waves |
CN102778693A (zh) * | 2011-05-13 | 2012-11-14 | 中国石油化工股份有限公司 | 一种基于反射波层拉平提取并消除的绕射波分离处理方法 |
CN103675897A (zh) * | 2012-08-30 | 2014-03-26 | 中国石油化工股份有限公司 | 一种地震绕射波分离成像方法 |
CN105785439A (zh) * | 2016-02-01 | 2016-07-20 | 北京中科联衡科技有限公司 | 小尺度非均匀地质体空间分布位置的预测方法和装置 |
US20180246241A1 (en) * | 2016-02-02 | 2018-08-30 | China University Of Mining & Technology, Beijing | Method and apparatus for identifying collapsed coal column |
US20180292553A1 (en) * | 2017-01-10 | 2018-10-11 | Institute Of Geology And Geophysics, Chinese Academy Of Sciences | Method and Apparatus for Separating Seismic Diffracted Wave |
CN107918147A (zh) * | 2017-11-20 | 2018-04-17 | 中国矿业大学(北京) | 绕射波成像方法和装置 |
CN108226999A (zh) * | 2018-01-19 | 2018-06-29 | 中国石油化工股份有限公司 | 碳酸盐岩小尺度缝洞体信息的处理方法 |
CN109507732A (zh) * | 2018-11-07 | 2019-03-22 | 中石化石油工程技术服务有限公司 | 一种基于成像道集的绕射波分离成像方法 |
CN111474580A (zh) * | 2020-05-13 | 2020-07-31 | 中国海洋石油集团有限公司 | 一种基于炮检距矢量片的方位角道集提取方法和系统 |
Non-Patent Citations (6)
Title |
---|
KONG XUE 等: "Diffraction separation by plane-wave prediction filtering", APPLIED GEOPHYSICS, vol. 14, no. 03, 30 September 2017 (2017-09-30), pages 399 - 405 * |
仲旭艳 等: "中值滤波及均值滤波在地震异常波提取中的应用", 地下水, vol. 39, no. 02, 31 March 2017 (2017-03-31), pages 222 - 224 * |
朱生旺 等: "局部倾角滤波和预测反演联合分离绕射波", 地球物理学报, vol. 56, no. 01, 31 January 2013 (2013-01-31), pages 280 - 288 * |
王志辉 等: "地震绕射波分离方法研究进展", 地球物理学进展, vol. 34, no. 01, 31 December 2019 (2019-12-31), pages 221 - 228 * |
薛明喜 等: "碳酸盐岩缝洞型储层三种不同地震成像效果对比", 科学技术与工程, vol. 18, no. 18, 30 June 2018 (2018-06-30), pages 7 - 15 * |
赵娟娟 等: "F-K滤波方法分离地震绕射波和反射波", 能源技术与管理, no. 03, 31 December 2010 (2010-12-31), pages 16 - 17 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113640872A (zh) * | 2021-08-12 | 2021-11-12 | 中国矿业大学(北京) | 绕射波分离方法、装置和电子设备 |
CN113640872B (zh) * | 2021-08-12 | 2022-03-08 | 中国矿业大学(北京) | 绕射波分离方法、装置和电子设备 |
CN114325821A (zh) * | 2021-12-23 | 2022-04-12 | 西安交通大学 | 基于3d-snacnn网络的叠前地震资料中强散射噪声压制方法及系统 |
CN115616660A (zh) * | 2022-12-20 | 2023-01-17 | 广州海洋地质调查局三亚南海地质研究所 | 利用绕射波监测海域碳封存项目co2泄露情况的方法及装置 |
CN115616660B (zh) * | 2022-12-20 | 2023-05-19 | 广州海洋地质调查局三亚南海地质研究所 | 利用绕射波监测海域碳封存项目二氧化碳泄露情况的方法及装置 |
CN117233844A (zh) * | 2023-09-19 | 2023-12-15 | 中国地质科学院地球物理地球化学勘查研究所 | 一种断层发育区地震高精度成像的数据处理方法及系统 |
CN117233844B (zh) * | 2023-09-19 | 2024-06-04 | 中国地质科学院地球物理地球化学勘查研究所 | 一种断层发育区地震高精度成像的数据处理方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112415591B (zh) | 2024-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112415591B (zh) | 绕射波的成像方法、装置、电子设备及存储介质 | |
CN103454678B (zh) | 一种地震切片等时性的确定方法及系统 | |
CN105116448B (zh) | 一种转换波方位各向异性校正方法及装置 | |
CN103635831A (zh) | 在3d各向异性介质中的基于敏感度内核的偏移速度分析 | |
WO2017035490A1 (en) | Optimal survey design | |
CN105301636A (zh) | 速度模型的建立方法和装置 | |
CN104199103A (zh) | 转换波静校正方法及装置 | |
CN109696705A (zh) | 深海obs上下行波场成像数据校正方法及装置 | |
CN105425298A (zh) | 一种消除有限差分正演过程中数值频散的方法和装置 | |
CN106054252B (zh) | 一种叠前时间偏移的方法及装置 | |
CN113433587B (zh) | 岩石物理反演方法及系统 | |
CN112684498A (zh) | 一种基于宽方位地震数据的储层裂缝预测方法及系统 | |
CN109459790B (zh) | 针对煤系地层地震速度场建立方法及系统 | |
CN112817042B (zh) | 沉积储层中油气识别方法及装置 | |
CN112711064B (zh) | 初始层速度建模方法及装置 | |
CN113933901B (zh) | 深度域地震数据拼接方法及装置 | |
CN112305585B (zh) | 三维数据体获取方法及装置 | |
CN109884701B (zh) | 地质体散射角度导向深度成像方法 | |
CN114114417B (zh) | 地震成像方法及装置 | |
CN112394408A (zh) | 各向异性介质裂缝预测方法及装置 | |
CN113075747B (zh) | 储层裂缝发育区域的预测方法及装置 | |
CN104166160B (zh) | 一种地震数据组合方法、装置及系统 | |
CN103852786B (zh) | 一种应用于陆上地震数据的逆时偏移成像的方法及系统 | |
CN112505748B (zh) | 炮域反射波拾取方法及装置 | |
CN110632650B (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 |