CN105785439B - 小尺度非均匀地质体空间分布位置的预测方法和装置 - Google Patents
小尺度非均匀地质体空间分布位置的预测方法和装置 Download PDFInfo
- Publication number
- CN105785439B CN105785439B CN201610069160.0A CN201610069160A CN105785439B CN 105785439 B CN105785439 B CN 105785439B CN 201610069160 A CN201610069160 A CN 201610069160A CN 105785439 B CN105785439 B CN 105785439B
- Authority
- CN
- China
- Prior art keywords
- msub
- mover
- diffracted wave
- small scale
- mrow
- 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.)
- Active
Links
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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
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)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种小尺度非均匀地质体空间分布位置的预测方法和装置。该方法和装置利用小尺度非均匀地质体的绕射地震波的特点,从叠前共偏移距道集中提取绕射波信息,并对绕射波信息进行偏移成像,接着对绕射波成像体有利于小尺度非均匀地质体空间分布预测的属性分析,最终提取到小尺度非均匀地质体的空间分布参数,根据提取到的空间分布参数从而实现对小尺度非均匀地质体空间分布位置的预测。因此,该方法解决了小尺度非均匀地质体空间分布的预测,是复杂构造地区的地质体预测的有效工具。
Description
技术领域
本发明涉及地质勘探领域,尤其涉及一种小尺度非均匀地质体预测方法和装置。
背景技术
随着地质勘探开发技术的发展,一些小尺度非均匀地质体,如灭尘点、小断块、裂缝等的探测越来越受到工业界重视。这些异常体的位置信息对识别油气储层有着重要意义。
然而,在地震勘探开发领域,传统的地震数据处理通常只利用反射波来获取地下地质构造信息,这种地震数据处理方式对大尺度连续地质体例如层状沉积地层成像较好,而针对尺度较小的地质体,其反射波成像中往往表现为杂乱反射,与噪音难以区分开来。
因此,利用传统的地震数据处理方式很难探测小尺度非均匀地质体的空间分布情况。
发明内容
有鉴于此,本发明的发明人研究发现,在小尺度非均匀地质体区域,绕射地震波发育明显,具有能量较弱,但成像分辨率高、识别能力强的特点,因此,通过对绕射波运动学和动力学特征分析,本发明提供了一种基于绕射波信息的小尺度非均匀地质体空间分布位置的预测方法和装置,以解决利用传统的地震数据处理方法很难探测小尺度非均匀地质体的空间分布情况。
为了达到上述发明目的,本发明采用了如下技术方案:
一种小尺度非均匀地质体空间分布位置的预测方法,包括:
从叠前共偏移距道集中提取绕射波信息;
对所述绕射波信息进行偏移成像,形成绕射波成像体;
对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析;
提取小尺度非均匀地质体的空间分布参数;
根据所述小尺度非均匀地质体的空间分布参数预测所述小尺度非均匀地质体的空间分布位置。
可选地,所述对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析,具体包括:
通过蚂蚁体技术和绕射波能量检测技术对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析。
可选地,所述从叠前共偏移距道集中提取绕射波信息,包括:
获取地震数据;
将所述地震数据分选为共中心点道集数据;
对共中心点道集数据进行动校正;
将动校正后的共中心点道集数据分选为共偏移距道集;
对所述共偏移距道集进行规则化处理;
在规则化处理后的共偏移距道集上,利用反射波同相轴局部倾角和相邻地震道数据,预测反射波同相轴;
从规则化处理后的共偏移距道集上对反射波进行自适应衰减,残余的数据即为动校正处理后的绕射波数据;
对所述动校正处理后的绕射波数据进行反动校正处理,获取最终分离的绕射波信息。
可选地,所述在规则化处理后的共偏移距道集上,利用反射波同相轴局部倾角和相邻地震道数据,预测反射波同相轴之前,还包括:获取反射波同相轴局部倾角。
可选地,所述获取反射波同相轴局部倾角,具体包括:
对地震数据中不同偏移距的地震道进行正常时差校正;
对正常时差校正后的地震数据进行叠加;
对叠加后的地震数据求取反射波同相轴局部倾角。
可选地,所述对所述绕射波信息进行偏移成像,形成所述绕射波成像体,具体包括:
基于克希霍夫积分偏移理论对绕射波信息进行偏移成像,形成所述绕射波成像体;
所述克希霍夫积分偏移理论的公式如下:
其中,为地表观测得到的波场值;
分别为炮点和检波点坐标;
为反射点处在t时刻的波场;
rs和rr分别代表炮点到反射点、检波点到反射点的距离;
vd和vu分别代表下行波和上行波沿射线路径的均方根速度;
系数A为保幅处理因子,作为振幅比例因子,实现保幅处理。
一种小尺度非均匀地质体空间分布位置的预测装置,包括:
第一提取单元,用于从叠前共偏移距道集中提取绕射波信息;
偏移成像单元,用于对所述绕射波信息进行偏移成像,形成绕射波成像体;
属性分析单元,用于对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析;
第二提取单元,用于提取小尺度非均匀地质体的空间分布参数;
预测单元,用于根据所述小尺度非均匀地质体的空间分布参数预测所述小尺度非均匀地质体的空间分布位置。
可选地,所述第一提取单元具体包括:
第一获取子单元,用于获取地震数据;
第一分选子单元,用于将所述地震数据分选为共中心点道集数据;
动校正子单元,用于对共中心点道集数据进行动校正;
第二分选子单元,用于将动校正后的共中心点道集数据分选为共偏移距道集;
规则化处理子单元,用于对所述共偏移距道集进行规则化处理;
预测子单元,用于在规则化处理后的共偏移距道集上,利用反射波同相轴局部倾角和相邻地震道数据,预测反射波同相轴;
自适应衰减子单元,用于从规则化处理后的共偏移距道集上对反射波进行自适应衰减,残余的数据即为动校正处理后的绕射波数据;
第二获取子单元,用于对所述动校正处理后的绕射波数据进行反动校正处理,获取最终分离的绕射波信息。
可选地,所述第一提取单元还包括:
第三获取子单元,用于获取反射波同相轴局部倾角,并触发所述预测子单元执行预测反射波同相轴的操作。
可选地,所述第三获取子单元具体包括:
正常时差校正子单元,用于对地震数据中不同偏移距的地震道进行正常时差校正;
叠加子单元,用于对正常时差校正后的地震数据进行叠加;
计算子单元,用于对叠加后的地震数据求取反射波同相轴局部倾角。
可选地,所述偏移成像单元具体包括基于克希霍夫积分偏移理论对绕射波信息进行偏移成像,形成所述绕射波成像体的子单元;
其中,所述克希霍夫积分偏移理论的公式如下:
其中,为地表观测得到的波场值;
分别为炮点和检波点坐标;
为反射点处在t时刻的波场;
rs和rr分别代表炮点到反射点、检波点到反射点的距离;
vd和vu分别代表下行波和上行波沿射线路径的均方根速度;
系数A为保幅处理因子,作为振幅比例因子,实现保幅处理。
相较于现有技术,本发明具有以下有益效果:
本发明提供的小尺度非均匀地质体空间分布位置的预测方法,利用小尺度非均匀地质体的绕射地震波的特点,从叠前共偏移距道集中提取绕射波信息,并对绕射波信息进行偏移成像,接着对绕射波成像体有利于小尺度非均匀地质体空间分布预测的属性分析,最终提取到小尺度非均匀地质体的空间分布参数,根据提取到的空间分布参数从而实现对小尺度非均匀地质体空间分布位置的预测。因此,该方法解决了小尺度非均匀地质体空间分布的预测,是复杂构造地区的地质体预测的有效方法。
附图说明
为了清楚地理解本发明的具体实施方式,下面将描述本发明的具体实施方式时用到的附图做一简要说明。显而易见地,这些附图仅是本发明的部分实施例,本领域技术人员在未付出创造性劳动的前提下,还可以获得其它附图。
图1是本发明实施例提供的小尺度非均匀地质体空间分布位置的预测方法流程示意图;
图2是本发明实施例提供的步骤S11的一个具体实施方式流程示意图;
图3是本发明实施例提供的小尺度非均匀地质体空间分布位置的预测装置结构示意图;
图4是本发明实施例提供的第一提取单元的具体结构示意图。
具体实施方式
为使本发明的发明目的、技术方案和技术效果更加清楚、完整,下面结合附图对本发明的具体实施方式进行详细描述。
首先需要说明的是,本发明实施例所述的小尺度非均匀地质体包括但不限于尖灭点、小断块和裂缝。
图1是本发明实施例提供的小尺度非均匀地质体空间分布位置的预测方法流程示意图。如图1所示,该预测方法包括以下步骤:
S11、从叠前共偏移距道集中提取绕射波信息:
在叠前共偏移距道集数据上,反射波和绕射波特征表现为:反射波同相轴局部光滑连续,能量腔,子波波形稳定,可预测;绕射波能量弱,波形不稳定,可预测性差。
根据上述特征差异,可以设计一滤波器,例如平面波分解滤波器,估算反射波同相轴局部倾角,然后利用滤波算子从全波场中取出强反射波,将连续性较差的绕射波信息保留下来。作为本发明的一个具体实施例,步骤S11的具体实现方式如图2所示,其包括以下步骤:
S111、获取地震数据:
需要说明的是,在本发明实施例中,地震数据一般为炮集数据。炮集数据为同一炮点激发,不同检波点接收的所有道上的数据。
S112、将所述地震数据分选为共中心点道集数据:
基于共中心点将地震数据分选为共中心点道集数据。共中心点道集数据是地震资料处理中最常用的一种道集数据形式。其道集中的所有道来自于同一个中心点。
S113、对共中心点道集数据进行动校正:
为了提高空间分布位置预测的精准度,本步骤对共中心点道集数据进行动校正处理。
S114、将动校正后的共中心点道集数据分选为共偏移距道集:
基于共偏移距将动校正后的共中心点道集数据分选为共偏移距道集。此时,每个共偏移距道集数据体都可以看作零偏移距数据体。
S115、对所述共偏移距道集进行规则化处理:
由于偏移距局部可能分布不均匀,存在缺道、坏道等问题,因此,需要对分选后的共偏移距道集数据进行规则化处理。规则化处理后的共偏移距道集数据才能应用局部倾角滤波。
S116、获取发射波同相轴局部倾角:
作为一个具体实施方式,获取反射波同相轴局部倾角的过程可以包括以下步骤:
A1、对地震数据中不同偏移距的地震道进行正常时差校正:
需要说明的是,正常时差校正即为动校正。
在本发明实施例中,为了做到各道同相叠加,需要对不同偏移距的地震道进行正常时差校正。其中,偏移距是指炮点到检波点之间的距离。
A2、对正常时差校正后的地震数据进行叠加:
需要说明的是,为了减少求取反射波同轴向局部倾角的计算量,本发明实施例优选对正常时差校正后的地震数据进行叠加。通过叠加处理,把来自同一个反射点上的反射波同相地叠加在一起,如此可以提高地震记录信噪比,有利于更精确地提取反射波同相轴倾角信息。
A3、对叠加后的地震数据求取反射波同相轴局部倾角:
需要说明的是,作为示例,本发明可以利用平面解构滤波算法从叠加后的地震数据中求取反射波同相轴局部倾角。此外,作为本发明实施例的扩展,本发明也可以采用其它滤波算法从叠加后的地震数据中求取反射波同轴相局部倾角。
S117、在规则化处理后的共偏移距道集上,利用反射波同相轴局部倾角和相邻地震数据,预测反射波同相轴:
步骤S117可以具体为:在规则化处理后的共偏移距道集上,利用反射波同相轴局部倾角和相邻地震数据,通过局部倾角滤波预测反射波同相轴。
S118、从规则化处理后的共偏移距道集上对反射波进行自适应衰减,残余的数据即为动校正处理后的绕射波数据。
S119、对所述动校正处理后的绕射波数据进行反动校正处理,获取最终分离的绕射波信息。
S12、对所述绕射波信息进行偏移成像,形成绕射波成像体:
反射波剖面断裂构造(如断层、裂缝系统等)不能很好成像,聚焦性差,同相轴凌乱;绕射波剖面中,绕射波在断裂构造区域性发育的位置能量分布集中,通过绕射波成像在断裂构造位置有明确显示,可应用其特性寻找并确定断裂系统的位置。因此,为了预测小尺度非均匀地质体的空间分布位置,需要形成绕射波成像体。
作为本发明的一个具体实施例可以基于克希霍夫积分偏移理论对上述提取到的绕射波信息进行偏移成像,形成所述绕射波成像体;
所述克希霍夫积分偏移理论的公式如下:
其中,为地表观测得到的波场值;
分别为炮点和检波点坐标;
为反射点处在t时刻的波场;
rs和rr分别代表炮点到反射点、检波点到反射点的距离;
vd和vu分别代表下行波和上行波沿射线路径的均方根速度;
系数A为保幅处理因子,作为振幅比例因子,实现保幅处理。
S13、对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析:
需要说明的是,由于绕射波成像体的信噪比较低,通过直接分析绕射成像体很难将小尺度非均匀地质体如小断层、裂缝系统等准确地分辨出来,这就需要将绕射波成像体进行有利于小尺度非君度地质体空间分布预测的属性分析。
在本发明实施例中,可以通过蚂蚁体技术和绕射波能量检测技术对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析。
其中,蚂蚁体技术又称断裂系统自动追踪技术,对小尺度非均匀地质体的识别非常有效。由于常规反射波成像方法只对大尺度层状模型成像较好,例如反射波成像层对层状沉积地层能够很好成像。当地质体的尺度较小时,在反射波成像中难以进行识别。反射波成像方法对小尺度不规则、非均匀地质体的裂缝带则成像难度较大,这些小尺度不规则、非均匀地质体的反射波成像体往往表现为杂乱反射,与噪音难以区分开来,识别难度大。针对裂缝,孔洞等小尺度非均匀地质体,绕射波具有比反射波强的识别能力。因此我们利用绕射数据进行蚂蚁体裂缝追踪,最终获得一个低噪音、具有清晰断裂痕迹的数据体。
经过研究发现,通过反射波和绕射波成像剖面的分析,在地层沉积稳定区域,构造破坏影响较小,绕射波能量相对较弱;而在构造发育区,受多期构造作用,地层破坏严重,反射波聚焦差,同相轴相对凌乱,绕射波对应于这些断裂系统,波场较发育。通过绕射信息分离及成像后,这些区域绕射能量相对较强,指示该地区裂缝发育程度较高。因此,绕射波能量检测技术是有利于小尺度非均匀地质体空间分布位置预测的属性分析手段。
S14、提取小尺度非均匀地质体的空间分布参数:
当小尺度非均匀地质体为裂缝时,所述小尺度非均匀地质体的空间分布参数包括裂缝密度和方向。
S15、根据所述小尺度非均匀地质体的空间分布参数预测所述小尺度非均匀地质体的空间分布位置。
以上为本发明实施例提供的小尺度非均匀地质体空间分布位置的预测方法的具体实施方式。在该具体实施方式中,利用小尺度非均匀地质体的绕射地震波的特点,从叠前共偏移距道集中提取绕射波信息,并对绕射波信息进行偏移成像,接着对绕射波成像体有利于小尺度非均匀地质体空间分布预测的属性分析,最终提取到小尺度非均匀地质体的空间分布参数,根据提取到的空间分布参数从而实现对小尺度非均匀地质体空间分布位置的预测。因此,该方法解决了小尺度非均匀地质体空间分布的预测,是复杂构造地区的地质体预测的有效方法。
基于上述实施例提供的小尺度非均匀地质体空间分布位置的预测方法,本发明实施例还提供了一种小尺度非均匀地质体空间分布位置的预测装置。具体参见以下实施例。
图3是本发明实施例提供的小尺度非均匀地质体空间分布位置的预测装置结构示意图。如图3所示,该预测装置包括以下单元:
第一提取单元31,用于从叠前共偏移距道集中提取绕射波信息;
偏移成像单元32,用于对所述绕射波信息进行偏移成像,形成绕射波成像体;
属性分析单元33,用于对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析;
第二提取单元34,用于提取小尺度非均匀地质体的空间分布参数;
预测单元35,用于根据所述小尺度非均匀地质体的空间分布参数预测所述小尺度非均匀地质体的空间分布位置。
作为本发明的一个具体实施例,如图4所示,所述第一提取单元31可以具体包括:
第一获取子单元311,用于获取地震数据;
第一分选子单元312,用于将所述地震数据分选为共中心点道集数据;
动校正子单元313,用于对共中心点道集数据进行动校正;
第二分选子单元314,用于将动校正后的共中心点道集数据分选为共偏移距道集;
规则化处理子单元315,用于对所述共偏移距道集进行规则化处理;
预测子单元316,用于在规则化处理后的共偏移距道集上,利用反射波同相轴局部倾角和相邻地震道数据,预测反射波同相轴;
自适应衰减子单元317,用于从规则化处理后的共偏移距道集上对反射波进行自适应衰减,残余的数据即为动校正处理后的绕射波数据;
第二获取子单元318,用于对所述动校正处理后的绕射波数据进行反动校正处理,获取最终分离的绕射波信息。
作为本发明的更具体实施例,上述所述的第一提取单元31还可以包括:
第三获取子单元319,用于获取反射波同相轴局部倾角,并触发所述预测子单元执行预测反射波同相轴的操作。
作为本发明的更进一步具体实施例,所述第三获取子单元319可以具体包括:
正常时差校正子单元,用于对地震数据中不同偏移距的地震道进行正常时差校正;
叠加子单元,用于对正常时差校正后的地震数据进行叠加;
计算子单元,用于对叠加后的地震数据求取反射波同相轴局部倾角。
作为本发明的另一具体实施例,所述偏移成像单元32具体包括基于克希霍夫积分偏移理论对绕射波信息进行偏移成像,形成所述绕射波成像体的子单元;
其中,所述克希霍夫积分偏移理论的公式如下:
其中,为地表观测得到的波场值;
分别为炮点和检波点坐标;
为反射点处在t时刻的波场;
rs和rr分别代表炮点到反射点、检波点到反射点的距离;
vd和vu分别代表下行波和上行波沿射线路径的均方根速度;
系数A为保幅处理因子,作为振幅比例因子,实现保幅处理。
以上为本发明实施例提供的小尺度非均匀地质体空间分布位置的预测装置的具体实施方式。在该具体实施方式中,利用小尺度非均匀地质体的绕射地震波的特点,从叠前共偏移距道集中提取绕射波信息,并对绕射波信息进行偏移成像,接着对绕射波成像体有利于小尺度非均匀地质体空间分布预测的属性分析,最终提取到小尺度非均匀地质体的空间分布参数,根据提取到的空间分布参数从而实现对小尺度非均匀地质体空间分布位置的预测。因此,该装置解决了小尺度非均匀地质体空间分布的预测,是复杂构造地区的地质体预测的有效工具。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (9)
1.一种小尺度非均匀地质体空间分布位置的预测方法,其特征在于,包括:
从叠前共偏移距道集中提取绕射波信息;对所述绕射波信息进行偏移成像,形成绕射波成像体;
对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析;
提取小尺度非均匀地质体的空间分布参数;
根据所述小尺度非均匀地质体的空间分布参数预测所述小尺度非均匀地质体的空间分布位置;
所述从叠前共偏移距道集中提取绕射波信息,包括:
获取地震数据;
将所述地震数据分选为共中心点道集数据;
对共中心点道集数据进行动校正;
将动校正后的共中心点道集数据分选为共偏移距道集;
对所述共偏移距道集进行规则化处理;
在规则化处理后的共偏移距道集上,利用反射波同相轴局部倾角和相邻地震道数据,预测反射波同相轴;
从规则化处理后的共偏移距道集上对反射波进行自适应衰减,残余的数据即为动校正处理后的绕射波数据;
对所述动校正处理后的绕射波数据进行反动校正处理,获取最终分离的绕射波信息。
2.根据权利要求1所述的方法,其特征在于,所述对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析,具体包括:
通过蚂蚁体技术和绕射波能量检测技术对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析。
3.根据权利要求2所述的方法,其特征在于,所述在规则化处理后的共偏移距道集上,利用反射波同相轴局部倾角和相邻地震道数据,预测反射波同相轴之前,还包括:获取反射波同相轴局部倾角。
4.根据权利要求3所述的方法,其特征在于,所述获取反射波同相轴局部倾角,具体包括:
对地震数据中不同偏移距的地震道进行正常时差校正;
对正常时差校正后的地震数据进行叠加;
对叠加后的地震数据求取反射波同相轴局部倾角。
5.根据权利要求1-4任一项所述的方法,其特征在于,所述对所述绕射波信息进行偏移成像,形成所述绕射波成像体,具体包括:
基于克希霍夫积分偏移理论对绕射波信息进行偏移成像,形成所述绕射波成像体;
所述克希霍夫积分偏移理论的公式如下:
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&Integral;</mo>
<mi>A</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mi>t</mi>
<mo>+</mo>
<mfrac>
<msub>
<mover>
<mi>r</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msub>
<msub>
<mi>v</mi>
<mi>d</mi>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mover>
<mi>r</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<msub>
<mi>v</mi>
<mi>u</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>d</mi>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msub>
<mi>d</mi>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>;</mo>
</mrow>
其中,为地表观测得到的波场值;
分别为炮点和检波点坐标;
为反射点处在t时刻的波场;
rs和rr分别代表炮点到反射点、检波点到反射点的距离;
vd和vu分别代表下行波和上行波沿射线路径的均方根速度;
系数A为保幅处理因子,作为振幅比例因子,实现保幅处理。
6.一种小尺度非均匀地质体空间分布位置的预测装置,其特征在于,包括:
第一提取单元,用于从叠前共偏移距道集中提取绕射波信息;
偏移成像单元,用于对所述绕射波信息进行偏移成像,形成绕射波成像体;
属性分析单元,用于对所述绕射波成像体进行有利于小尺度非均匀地质体空间分布预测的属性分析;
第二提取单元,用于提取小尺度非均匀地质体的空间分布参数;
预测单元,用于根据所述小尺度非均匀地质体的空间分布参数预测所述小尺度非均匀地质体的空间分布位置;
所述第一提取单元具体包括:
第一获取子单元,用于获取地震数据;
第一分选子单元,用于将所述地震数据分选为共中心点道集数据;
动校正子单元,用于对共中心点道集数据进行动校正;
第二分选子单元,用于将动校正后的共中心点道集数据分选为共偏移距道集;
规则化处理子单元,用于对所述共偏移距道集进行规则化处理;
预测子单元,用于在规则化处理后的共偏移距道集上,利用反射波同相轴局部倾角和相邻地震道数据,预测反射波同相轴;
自适应衰减子单元,用于从规则化处理后的共偏移距道集上对反射波进行自适应衰减,残余的数据即为动校正处理后的绕射波数据;
第二获取子单元,用于对所述动校正处理后的绕射波数据进行反动校正处理,获取最终分离的绕射波信息。
7.根据权利要求6所述的装置,其特征在于,所述第一提取单元还包括:
第三获取子单元,用于获取反射波同相轴局部倾角,并触发所述预测子单元执行预测反射波同相轴的操作。
8.根据权利要求7所述的装置,其特征在于,所述第三获取子单元具体包括:
正常时差校正子单元,用于对地震数据中不同偏移距的地震道进行正常时差校正;
叠加子单元,用于对正常时差校正后的地震数据进行叠加;
计算子单元,用于对叠加后的地震数据求取反射波同相轴局部倾角。
9.根据权利要求6-8任一项所述的装置,其特征在于,所述偏移成像单元具体包括基于克希霍夫积分偏移理论对绕射波信息进行偏移成像,形成所述绕射波成像体的子单元;
其中,所述克希霍夫积分偏移理论的公式如下:
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&Integral;</mo>
<mi>A</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>z</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mi>t</mi>
<mo>+</mo>
<mfrac>
<msub>
<mover>
<mi>r</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msub>
<msub>
<mi>v</mi>
<mi>d</mi>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mover>
<mi>r</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
<msub>
<mi>v</mi>
<mi>u</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>d</mi>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>s</mi>
</msub>
<mi>d</mi>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>r</mi>
</msub>
</mrow>
其中,为地表观测得到的波场值;
分别为炮点和检波点坐标;
为反射点处在t时刻的波场;
rs和rr分别代表炮点到反射点、检波点到反射点的距离;
vd和vu分别代表下行波和上行波沿射线路径的均方根速度;
系数A为保幅处理因子,作为振幅比例因子,实现保幅处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610069160.0A CN105785439B (zh) | 2016-02-01 | 2016-02-01 | 小尺度非均匀地质体空间分布位置的预测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610069160.0A CN105785439B (zh) | 2016-02-01 | 2016-02-01 | 小尺度非均匀地质体空间分布位置的预测方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105785439A CN105785439A (zh) | 2016-07-20 |
CN105785439B true CN105785439B (zh) | 2018-05-18 |
Family
ID=56402356
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610069160.0A Active CN105785439B (zh) | 2016-02-01 | 2016-02-01 | 小尺度非均匀地质体空间分布位置的预测方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105785439B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443785B (zh) * | 2016-11-03 | 2018-05-29 | 中国矿业大学(北京) | 绕射波场获取方法和装置 |
CN106772592B (zh) * | 2016-11-10 | 2018-08-07 | 中国矿业大学(北京) | 绕射波聚焦能量的分析方法及装置 |
CN108226999A (zh) * | 2018-01-19 | 2018-06-29 | 中国石油化工股份有限公司 | 碳酸盐岩小尺度缝洞体信息的处理方法 |
CN109116423B (zh) * | 2018-07-11 | 2020-03-17 | 北京奥能恒业能源技术有限公司 | 一种绕射多次波压制方法及装置 |
CN112083490B (zh) * | 2019-06-12 | 2024-03-01 | 中国石油天然气集团有限公司 | 地震数据噪声衰减方法及装置 |
CN111025383B (zh) * | 2019-11-21 | 2021-09-24 | 徐州工程学院 | 一种基于绕射横波定性判断隧道前方溶洞充水情况的方法 |
CN111239832B (zh) * | 2020-03-10 | 2020-12-22 | 中国地质大学(北京) | 3d3c-vsp成像处理方法、装置及设备 |
CN113671573A (zh) * | 2020-05-14 | 2021-11-19 | 中国石油化工股份有限公司 | 地震数据叠前绕射波分离方法、装置、电子设备及介质 |
CN112415591A (zh) * | 2020-10-30 | 2021-02-26 | 中国石油天然气集团有限公司 | 绕射波的成像方法、装置、电子设备及存储介质 |
CN115343763B (zh) * | 2021-05-14 | 2023-08-01 | 中国石油化工股份有限公司 | 基于地震同相轴零相位判别的蚂蚁体属性倾角校正方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102455439A (zh) * | 2010-11-02 | 2012-05-16 | 中国石油大学(北京) | 基于克希霍夫积分法的绕射波场分离方法 |
CN103675896A (zh) * | 2012-08-30 | 2014-03-26 | 中国石油化工股份有限公司 | 一种绕射波与反射波分离成像方法 |
CN104730571A (zh) * | 2015-03-11 | 2015-06-24 | 中国科学院地质与地球物理研究所 | 一种利用绕射再聚焦识别小尺度地质体的方法与装置 |
CN104730572A (zh) * | 2015-03-11 | 2015-06-24 | 中国科学院地质与地球物理研究所 | 一种基于l0半范数的绕射波成像方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8209125B2 (en) * | 2007-03-12 | 2012-06-26 | Geomage (2003) Ltd. | Method for identifying and analyzing faults/fractures using reflected and diffracted waves |
-
2016
- 2016-02-01 CN CN201610069160.0A patent/CN105785439B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102455439A (zh) * | 2010-11-02 | 2012-05-16 | 中国石油大学(北京) | 基于克希霍夫积分法的绕射波场分离方法 |
CN103675896A (zh) * | 2012-08-30 | 2014-03-26 | 中国石油化工股份有限公司 | 一种绕射波与反射波分离成像方法 |
CN104730571A (zh) * | 2015-03-11 | 2015-06-24 | 中国科学院地质与地球物理研究所 | 一种利用绕射再聚焦识别小尺度地质体的方法与装置 |
CN104730572A (zh) * | 2015-03-11 | 2015-06-24 | 中国科学院地质与地球物理研究所 | 一种基于l0半范数的绕射波成像方法及装置 |
Non-Patent Citations (1)
Title |
---|
"绕射波叠前时间偏移速度分析及成像";刘玉金 等;《地球物理学进展》;20131231;第28卷(第6期);第3023-3028页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105785439A (zh) | 2016-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105785439B (zh) | 小尺度非均匀地质体空间分布位置的预测方法和装置 | |
RU2694621C1 (ru) | Способ и устройство для обработки сейсмических данных | |
US9766357B2 (en) | Seismic image dip decomposition estimation and recomposition | |
CN105372705B (zh) | 一种基于多波资料的地层切片方法 | |
Palm et al. | A reflection seismic image of Caledonian structure in central Sweden | |
CN102590862B (zh) | 补偿吸收衰减的叠前时间偏移方法 | |
CN102053261A (zh) | 一种地震数据处理方法 | |
CN111880219B (zh) | 基于方位角-倾角道集的绕射波成像方法及装置 | |
EP1141750A1 (en) | Hydrocarbon edge detection using seismic amplitude | |
CN104459794A (zh) | 共反射点道集时变时间差值的校正方法及装置 | |
CN105425299B (zh) | 确定地层裂缝分布的方法和装置 | |
CN106772593B (zh) | 绕射波的成像方法及装置 | |
CN106526678A (zh) | 一种反射声波测井的波场分离方法及装置 | |
CN102053260B (zh) | 获得地震纵波的方位速度的方法及处理地震数据的方法 | |
CN106443791B (zh) | 求取倾斜地层或各向异性地层横波剩余静校正量的方法 | |
CN102636813B (zh) | 一种用于处理物探地震数据的动校正拉伸切除方法 | |
CN107942405A (zh) | 预测砂泥岩薄互层砂体累积厚度的方法 | |
CN104422955B (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
Craft et al. | Azimuthal anisotropy analysis from P-wave seismic traveltime data | |
CN102236101A (zh) | 多层裂缝预测方法和装置 | |
CN105093318A (zh) | 一种自适应波动方程波场延拓静校正方法 | |
CN108121010B (zh) | 基于孔巷联合的井下独头巷道槽波超前探测方法及系统 | |
CN103941289B (zh) | 一种绕射点成像的方法和设备 | |
Rodriguez-Tablante et al. | Cross-profile acquisition and cross-dip analysis for extracting 3D information from 2D surveys, a case study from the western Skellefte District, northern Sweden | |
CN104570115A (zh) | 一种面波衰减方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |