CN104656139A - 一种转换波的静校正方法及设备 - Google Patents
一种转换波的静校正方法及设备 Download PDFInfo
- Publication number
- CN104656139A CN104656139A CN201310594371.2A CN201310594371A CN104656139A CN 104656139 A CN104656139 A CN 104656139A CN 201310594371 A CN201310594371 A CN 201310594371A CN 104656139 A CN104656139 A CN 104656139A
- Authority
- CN
- China
- Prior art keywords
- data
- shot
- correction
- wave
- converted 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种转换波的静校正方法及设备,包括:采集三分量地震数据;对所述的三分量地震数据进行预处理,得到预处理后的转换波数据;获取纵波野外炮点静校正量;根据所述预处理后的转换波数据以及所述的纵波野外炮点静校正量得到动校正后的转换波共检波点道集;根据动校正后的转换波共检波点道集得到炮点剩余静校正后的转换波数据;根据所述的炮点剩余静校正后的转换波数据得到炮点和检波点静校正后的转换波数据。解决了现有技术中由于转换波近地表比较复杂而存在的适用条件严格以及计算精度差的问题,实现了在一定程度上适应转换波检波点上静校正影响大于1/2个周期的情况,有效提高了转换波静校正的精度和效率。
Description
技术领域
本发明关于地球物理勘探技术领域,特别是关于多波地震勘探技术,具体的讲是一种转换波的静校正方法及设备。
背景技术
多波地震勘探是综合利用纵波、横波、转换波等多种信息对油气藏进行精细勘察的一种勘探技术。这种技术可以识别薄互层、小断层、低幅度构造,准确预测储层边界和储层物性、监测流体,甚至还可以直接用于油气预测等。
目前,随着地震勘探技术的不断提高和成熟,转换波勘探日益增多。转换波地震资料的信噪比在一般情况下都比纵波低,再加上静校正量比纵波大,炮点校正量的精度和效果严重依赖于纵波静校正量的精度和效果,对于陆上复杂山地地震勘探,纵波静校正问题仍然是地震勘探的一个难点,对于复杂山地多波地震勘探,转换波静校正问题是多波地震资料处理的一个瓶颈。通常,多波地震勘探中的转换波静校正效果是衡量转换波地震资料处理技术水平高低的重要标志,转换波静校正量是由炮点的纵波静校正量和接收点转换波静校正量组成。由于在同一介质中转换波的传播速度比纵波速度低,且转换波不受孔隙流体影响,这就造成同一位置的转换波静校正量大于纵波静校正量,一般情况转换波校正量可能是纵波静校正量的2到10倍。
转换波的初至波静校正,由于转换波传播路径不对称,使得常规静校正算法难以在转换波静校正中直接应用。再加上转换波记录信噪比通常比较低,初至附近纵波折射、转换波折射等信息混杂在一起,使得转换波初至难以分辨和拾取,这在很大程度上加大了转换波静校正的难度。尤其是在许多复杂勘探区域,纵波静校正问题目前还没有完全解决,而转换波低速带又更加复杂,有时可能还要受到方位各向异性的影响,依靠初至波静校正方法去解决转换波静校证问题就很难得到实现。
转换波的反射波静校正,当前地震资料处理转换波静校正解决方案还主要是依靠共检波点叠加相关的静校正方法,这种静校正方法有两个假设条件1)假设转换波炮点校正量精度已经达到了转换波资料处理的要求;2)转换波检波点上静校正影响不大于1/2个周期。由于转换波近地表比较复杂,尽管这种算法在实际应用中也取得了一定的效果,但是存在着诸如适用条件严格和计算精度差等问题。
发明内容
为了克服现有技术存在的上述问题,本发明提供了一种转换波的静校正方法及设备,解决了现有技术中由于转换波近地表比较复杂而存在的适用条件严格以及计算精度差的问题,实现了在一定程度上适应转换波检波点上静校正影响大于1/2个周期的情况,有效提高了转换波静校正的精度和效率。
本发明的目的之一是,提供一种转换波的静校正方法,包括:采集三分量地震数据;对所述的三分量地震数据进行预处理,得到预处理后的转换波数据;获取纵波野外炮点静校正量;根据所述预处理后的转换波数据以及所述的纵波野外炮点静校正量得到动校正后的转换波共检波点道集;根据动校正后的转换波共检波点道集得到炮点剩余静校正后的转换波数据;根据所述的炮点剩余静校正后的转换波数据得到炮点和检波点静校正后的转换波数据。
本发明的目的之一是,提供了一种转换波的静校正设备,包括:三分量地震数据采集装置,用于采集三分量地震数据;预处理装置,用于对所述的三分量地震数据进行预处理,得到预处理后的转换波数据;野外炮点静校正量获取装置,用于获取纵波野外炮点静校正量;转换波共检波点道集生成装置,用于根据所述预处理后的转换波数据以及所述的纵波野外炮点静校正量得到动校正后的转换波共检波点道集;第一转换波数据生成装置,用于根据动校正后的转换波共检波点道集得到炮点剩余静校正后的转换波数据;第二转换波数据生成装置,用于根据所述的炮点剩余静校正后的转换波数据得到炮点和检波点静校正后的转换波数据。
本发明的有益效果在于,提供了一种转换波的静校正方法及设备,解决了现有技术中由于转换波近地表比较复杂而存在的适用条件严格以及计算精度差的问题,实现了后续常规转换波静校正处理不受转换波静校正影响大于1/2个周期情况的限制,转换波反射波静校正可以在本发明获得数据的基础上继续进行,在一定程度上提高了转换波静校正的精度和效果。
为让本发明的上述和其他目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附图式,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种转换波的静校正方法的流程图;
图2为图1中的步骤S102的具体流程图;
图3为图1中的步骤S104的具体流程图;
图4为图1中的步骤S105的具体流程图;
图5为图1中的步骤S106的具体流程图;
图6为本发明实施例提供的一种转换波的静校正设备的结构框图;
图7为本发明实施例提供的一种转换波的静校正设备中预处理装置200的具体结构框图;
图8为本发明实施例提供的一种转换波的静校正设备中转换波共检波点道集生成装置400的结构框图;
图9为本发明实施例提供的一种转换波的静校正设备中第一转换波数据生成装置500的结构框图;
图10为本发明实施例提供的一种转换波的静校正设备中第二转换波数据生成装置600的结构框图;
图11为本发明提供的实施例一中动校正后的转换波单炮模型记录示意图;
图12为本发明提供的实施例一中对图11的数据处理后得到的转换波单炮模型记录示意图;
图13为本发明提供的实施例二中动校正后的转换波共检波点道集实际地震记录示意图;
图14为发明提供的实施例二中对图13的数据处理后得到的共检波点道集地震记录示意图;
图15为根据图14形成的共检波点道集分选成共炮点道集录示意图;
图16为对图15的数据处理后的共炮点道集地震记录示意图;
图17为对图15的数据处理时计算出的检波点校正量示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种转换波的静校正方法及设备,解决了现有技术中由于转换波近地表比较复杂而存在的适用条件严格以及计算精度差的问题,实现了在一定程度上适应转换波检波点上静校正影响大于1/2个周期的情况,有效提高了转换波静校正的精度和效率。
图1为本发明提供的一种转换波的静校正方法的流程图,由图1可知,该方法具体包括:
S101:采集三分量地震数据。
S102:对所述的三分量地震数据进行预处理,得到预处理后的转换波数据。图2为步骤S102的具体流程图,由图2可知,步骤S102具体包括:
S201:对采集的三分量地震数据进行观测系统加载。观测系统加载为现有技术中已有的模块,这个模块主要是对地震数据的一些特征进行描述和记录,包括炮点坐标,检波点坐标地震道在道集中的位置,这个地震道属于哪个面元等。
S202:对观测系统加载后的三分量地震数据进行面波衰减;
S203:对面波衰减后的三分量地震数据进行规则噪声衰减;
S204:对规则噪声衰减后的三分量地震数据进行强能量压制;
S205:对强能量压制后的三分量地震数据进行能量补偿;
S206:对能量补偿后的三分量地震数据进行振幅恢复,得到预处理后的转换波数据。
由图1可知,该方法还包括:
S103:获取纵波野外炮点静校正量。该步骤中纵波野外炮点静校正量的获取方式有两个途径,可以从野外提供的近地表信息中直接获得炮点校正量,也可以在室内通过在纵波上拾取初至波后重新计算炮点校正量。
如果野外提供的近地表信息中炮点校正量精度比较高,可以在转换波数据上直接应用炮点校正量,无需室内再计算野外校正量。如果野外提供的校正量不可靠,就必须在室内重新计算炮点校正量,计算校正量时初至的拾取必须是在加载观测系统后的纵波数据上拾取。静校正方法可以选择常规纵波静校正方法,如层析静校正和折射波静校正两种方法。地表比较复杂的一般选择层析静校正;地表高程起伏不大,近地表横向速度变化缓慢的可以选择折射波静校正的方法获得校正量。
S104:根据所述预处理后的转换波数据以及所述的纵波野外炮点静校正量得到动校正后的转换波共检波点道集。图3为步骤S104的具体流程图,由图3可知,步骤S104具体包括:
S301:将预处理后的转换波数据应用于所述的纵波野外炮点静校正量,得到炮点校正后的转换波数据。
S302:对所述炮点校正后的转换波数据进行共中心点道集分选,得到转换波速度。在具体的实施方式中,对步骤S301形成的炮点校正后的转换波数据进行共中心点道集(CCP)分选,然后可做转换波速度分析,进而获得转换波速度。
S303:对所述炮点校正后的转换波数据进行共检波点道集分选,得到分选后的转换波共检波点道集数据。在具体的实施方式中,对步骤S301形成的炮点校正后的转换波数据进行共检波点道集分选,得到分选后的转换波共检波点道集数据。
S304:根据所述的转换波速度对分选后的转换波共检波点道集数据进行动校正,得到动校正后的转换波共检波点道集。
由图1可知,该方法还包括:
S105:根据动校正后的转换波共检波点道集得到炮点剩余静校正后的转换波数据。图4为步骤S105的具体流程图,由图4可知,步骤S105具体包括:
S401:对动校正后的转换波共检波点道集进行互相关处理,得到综合相关时差。
互相关处理在具体的实施方式中可通过如下方式实现:对于一个转换波共检波点道集的一个地震道,先选定一个时间窗诸如200毫秒,200毫秒内需要有一组较强的反射信号同相轴,窗的大小根据资料实际情况定,校正量越大,相应的时间窗选择就大。再给定一个地震道个数诸如5,也可以是3,4,6,7,8等,用这个地震道与其左边的4个地震道依次做相关,获得5个最大相关值时的时差,时差分别为t1、t2、t3、t4。再用这个地震道与其右边的4个地震道依次做相关,获得最大相关值时的时差分别为t6、t7、t8、t9,然后可以求得一个最大相关值时的综合相关时差:
t5=(1*t1+2*t2+3*t3+4*t4+4*t6+3*t7+2*t8+1*t9)/8
t5就是这个选定地震道的炮点校正量,也是步骤S401的综合相关时差。用步骤S401对每个转换波共检波点道集的每个地震道都进行相同处理,每个地震道都获得一个综合相关时差。
S402:对所述的综合相关时差进行地表一致性处理,得到转换波炮点剩余静校正量。
地表一致性处理在具体的实施方式中可通过如下方式实现:把每个地震道获得的综合相关时差,按炮点进行统计,把地震道属于同一炮的相关时差相加后求平均值,本发明把这个平均时差当成这个炮点剩余静校正量。这里的剩余静校正量主要是针对步骤S301已经用过一次纵波炮点校正量,由于纵波炮点校正量精度不够,还剩余了一部分炮点校正量。如果步骤S301使用的炮点的校正量精度比较高,炮点上的剩余静校正量就可以忽略,本发明的求炮点剩余静校正量的步骤S401和步骤S402就可以不被采用。
S403:将所述的转换波炮点剩余静校正量应用于所述的转换波共检波点道集上,得到炮点剩余静校正后的转换波数据。
由图1可知,该方法还包括:
S106:根据所述的炮点剩余静校正后的转换波数据得到炮点和检波点静校正后的转换波数据。图5为图1中的步骤S106的具体流程图,由图5可知,步骤S106具体包括:
S501:将所述的炮点剩余静校正后的转换波数据分选成共炮点道集;
S502:对所述的共炮点道集进行互相关处理,得到所述共炮点道集对应的综合相关时差。该步骤的地表一致性是把每个地震道获得的综合相关时差,按检波点进行统计,把地震道属于同一检波点的综合相关时差相加后求平均值,本发明把这个平均时差当成这个检波点的校正量。
S503:对所述共炮点道集对应的综合相关时差进行地表一致性处理,得到转换波检点校正量;
S504:将所述的转换波检点校正量应用于所述炮点剩余静校正后的转换波数据,得到炮点和检波点静校正后的转换波数据。
如上所述,步骤S105对这个转换波共检波点道集数据做互相关处理,获得一个综合相关时差,对综合时差进行地表一致性分解,可以获得炮点剩余静校正量,在应用炮点剩余静校正量数据的基础上可以用类似的方法获得检波点校正量。应用剩余炮点校正量后再应用检波点校正量,就可以获得本发明转换波校正后的数据,然后就可以进行后续常规的转换波静校正处理,这样可以在一定程度上避免校正量大于1/2个周期时对常规转换波静校正造成的影响,达到本发明提高转换波校正静校正精度和效率的目的。
图6为本发明实施例提供的一种转换波的静校正设备的结构框图,由图6可知,该设备具体包括:
三分量地震数据采集装置100,用于采集三分量地震数据。
预处理装置200,用于对所述的三分量地震数据进行预处理,得到预处理后的转换波数据。图7为转换波的静校正设备中预处理装置200的具体结构框图,由图7可知,预处理装置200具体包括:
加载模块201,用于对采集的三分量地震数据进行观测系统加载。
面波衰减模块202,用于对观测系统加载后的三分量地震数据进行面波衰减;
噪声衰减模块203,用于对面波衰减后的三分量地震数据进行规则噪声衰减;
强能量压制模块204,用于对规则噪声衰减后的三分量地震数据进行强能量压制;
能量补偿模块205,用于对强能量压制后的三分量地震数据进行能量补偿;
振幅恢复模块206,用于对能量补偿后的三分量地震数据进行振幅恢复,得到预处理后的转换波数据。
由图6可知,该设备还包括:
野外炮点静校正量获取装置300,用于获取纵波野外炮点静校正量。纵波野外炮点静校正量的获取方式有两个途径,可以从野外提供的近地表信息中直接获得炮点校正量,也可以在室内通过在纵波上拾取初至波后重新计算炮点校正量。
如果野外提供的近地表信息中炮点校正量精度比较高,可以在转换波数据上直接应用炮点校正量,无需室内再计算野外校正量。如果野外提供的校正量不可靠,就必须在室内重新计算炮点校正量,计算校正量时初至的拾取必须是在加载观测系统后的纵波数据上拾取。静校正方法可以选择常规纵波静校正方法,如层析静校正和折射波静校正两种方法。地表比较复杂的一般选择层析静校正;地表高程起伏不大,近地表横向速度变化缓慢的可以选择折射波静校正的方法获得校正量。
转换波共检波点道集生成装置400,用于根据所述预处理后的转换波数据以及所述的纵波野外炮点静校正量得到动校正后的转换波共检波点道集。图8为转换波的静校正设备中转换波共检波点道集生成装置400的结构框图,由图8可知,转换波共检波点道集生成装置具体包括:
第一应用模块401,用于将预处理后的转换波数据应用于所述的纵波野外炮点静校正量,得到炮点校正后的转换波数据。
共中心点道集分选模块402,用于对所述炮点校正后的转换波数据进行共中心点道集分选,得到转换波速度。在具体的实施方式中,对步骤S301形成的炮点校正后的转换波数据进行共中心点道集(CCP)分选,然后可做转换波速度分析,进而获得转换波速度。
共检波点道集分选模块403,用于对所述炮点校正后的转换波数据进行共检波点道集分选,得到分选后的转换波共检波点道集数据。在具体的实施方式中,对步骤S301形成的炮点校正后的转换波数据进行共检波点道集分选,得到分选后的转换波共检波点道集数据。
动校正模块404,用于根据所述的转换波速度对分选后的转换波数据进行动校正,得到动校正后的转换波共检波点道集。
由图6可知,该设备还包括:
第一转换波数据生成装置500,用于根据动校正后的转换波共检波点道集得到炮点剩余静校正后的转换波数据。图9为转换波的静校正设备中第一转换波数据生成装置500的结构框图,由图9可知,第一转换波数据生成装置具体包括:
第一互相关处理模块501,用于对动校正后的转换波共检波点道集进行互相关处理,得到综合相关时差。
互相关处理在具体的实施方式中可通过如下方式实现:对于一个转换波共检波点道集的一个地震道,先选定一个时间窗诸如200毫秒,200毫秒内需要有一组较强的反射信号同相轴,窗的大小根据资料实际情况定,校正量越大,相应的时间窗选择就大。再给定一个地震道个数诸如5,也可以是3,4,6,7,8等,用这个地震道与其左边的4个地震道依次做相关,获得5个最大相关值时的时差,时差分别为t1、t2、t3、t4。再用这个地震道与其右边的4个地震道依次做相关,获得最大相关值时的时差分别为t6、t7、t8、t9,然后可以求得一个最大相关值时的综合相关时差:
t5=(1*t1+2*t2+3*t3+4*t4+4*t6+3*t7+2*t8+1*t9)/8
t5就是这个选定地震道的炮点校正量,也是步骤S401的综合相关时差。用步骤S401对每个转换波共检波点道集的每个地震道都进行相同处理,每个地震道都获得一个综合相关时差。
第一一致性处理模块502,用于对所述的综合相关时差进行地表一致性处理,得到转换波炮点剩余静校正量。
地表一致性处理在具体的实施方式中可通过如下方式实现:把每个地震道获得的综合相关时差,按炮点进行统计,把地震道属于同一炮的相关时差相加后求平均值,本发明把这个平均时差当成这个炮点剩余静校正量。这里的剩余静校正量主要是针对步骤S301已经用过一次纵波炮点校正量,由于纵波炮点校正量精度不够,还剩余了一部分炮点校正量。如果步骤S301使用的炮点的校正量精度比较高,炮点上的剩余静校正量就可以忽略,本发明的求炮点剩余静校正量的步骤S401和步骤S402就可以不被采用。
第二应用模块503,用于将所述的转换波炮点剩余静校正量应用于所述的转换波共检波点道集上,得到炮点剩余静校正后的转换波数据。
由图6可知,该设备还包括:
第二转换波数据生成装置600,用于根据所述的炮点剩余静校正后的转换波数据得到炮点和检波点静校正后的转换波数据。图10为转换波的静校正设备中第二转换波数据生成装置600的结构框图,由图10可知,第二转换波数据生成装置具体包括:
分选模块601,用于将所述的炮点剩余静校正后的转换波数据分选成共炮点道集;
第二相关处理模块602,用于对所述的共炮点道集进行互相关处理,得到所述共炮点道集对应的综合相关时差。该步骤的地表一致性是把每个地震道获得的综合相关时差,按检波点进行统计,把地震道属于同一检波点的综合相关时差相加后求平均值,本发明把这个平均时差当成这个检波点的校正量。
第二一致性处理模块603,用于对所述共炮点道集对应的综合相关时差进行地表一致性处理,得到转换波检点校正量;
第三应用模块604,用于将所述的转换波检点校正量应用于所述炮点剩余静校正后的转换波数据,得到炮点和检波点静校正后的转换波数据。
如上所述,第一转换波数据生成装置对这个转换波共检波点道集数据做互相关处理,获得一个综合相关时差,对综合时差进行地表一致性分解,可以获得炮点剩余静校正量,在应用炮点剩余静校正量数据的基础上可以用类似的方法获得检波点校正量。应用剩余炮点校正量后再应用检波点校正量,就可以获得本发明转换波校正后的数据,然后就可以进行后续常规的转换波静校正处理,这样可以在一定程度上避免校正量大于1/2个周期时对常规转换波静校正造成的影响,达到本发明提高转换波校正静校正精度和效率的目的。
下面结合具体的实施例,详细介绍本发明的技术方案。
实施例一
本发明提供一种提高转换波静校正精度和效率的方法,运用于转换波数据静校正,对常规转换波静校正进行优化处理。图11为本发明提供的实施例一中动校正后的转换波单炮模型记录示意图,炮点上已经没有校正量,只剩下检波点上的校正量,所以可以不进行炮点剩余静校正量的求取,只在动校正后的炮点上计算检波点的校正量。从图11可以看到转换波同相轴抖动剧烈,说明检波点校正量很大,静校正的影响超过了半个排列(横向坐标是道号,纵向坐标是时间)。本发明首先在图11的单炮模型记录上选择了0.4秒到1.6秒的相关时间窗,把单炮模型记录上的三个同相轴都包括在内,求得综合相关时差,进行地表一致性分解,求得转换波检波点校正量,再把这个校正量应用于图11,获得本发明处理后的转换波数据,如图12。从图12可以发现检波点上的校正量基本消失,转换波检波点上静校正的影响也远远小于了1/2个周期,达到了本发明优化转换波静校正的目的。
实施例二
以本发明对实际地震记录的处理结果为例。图13是一个动校正后的转换波单炮实际地震记录,从图13上可以看出炮点上还存在一些剩余的校正量,这些校正量不解决会影响本方发明求取检波点道集转换波静校正量的精度。所以在求检波点静校正前,必须进行炮点剩余静校正量的求取和应用。图14就是本发明在检波点上求炮点剩余静校正量并应用的结果。图13上存在明显的同相轴抖动,本发明选择的时间窗是1700毫秒到2200毫秒,这个时间窗把2000毫秒左右的同相轴包含在内。经过本发明处理,求得综合相关时差,进行地表一致性分解,求得转换波炮点上剩余静校正量,再把这个校正量应用于图13的数据,本发明处理后炮点上剩余静校正量基本上被解决,同相轴变得很光滑,如图14。把炮点剩余静校正量应用后共检波点数据分选成共炮点道集数据,如图15。图16为对图15的数据处理后的共炮点道集地震记录示意图。在炮点上用同样的方法计算检波点的校正量,相关时间窗选择1.5秒到2.5秒,确保了剧烈抖动的同相轴都被包括在时间窗内。图15中同相轴抖动剧烈说明检波点校正量很大,在炮记录显示上可以明显看出静校正影响超过了一个排列,本发明计算出的检波点校正量显示图也证明了这一点,如图17,图17纵坐标为综合相关时差,横坐标为地震道顺序号。本方法求得综合相关时差,进行地表一致性分解,求得转换波检波点上的校正量,如图17,再把这个校正量应用于图15数据,本发明处理后检波点上校正量基本上得到解决,同相轴变得光滑,转换波检波点上静校正的影响也已经远远小于了1/2个周期,已经达到了本发明优化转换波静校正的目的。本发明处理后,检波点上可能还有一部分剩余静校正量,这种剩余静校正量完全可以在后续转换波静校正中用常规的转换波静校正方法进行校正处理解决。
综上所述,本发明的有益成果是:提供了一种转换波的静校正方法及设备,解决了现有技术中由于转换波近地表比较复杂而存在的适用条件严格以及计算精度差的问题,实现了后续常规转换波静校正处理不受转换波静校正影响大于1/2个周期情况的限制,转换波反射波静校正可以在本发明获得数据的基础上继续进行,在一定程度上提高了转换波静校正的精度和效果。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一般计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random AccessMemory,RAM)等。
本领域技术人员还可以了解到本发明实施例列出的各种功能是通过硬件还是软件来实现取决于特定的应用和整个系统的设计要求。本领域技术人员可以对于每种特定的应用,可以使用各种方法实现所述的功能,但这种实现不应被理解为超出本发明实施例保护的范围。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种转换波的静校正方法,其特征是,所述的方法具体包括:
采集三分量地震数据;
对所述的三分量地震数据进行预处理,得到预处理后的转换波数据;
获取纵波野外炮点静校正量;
根据所述预处理后的转换波数据以及所述的纵波野外炮点静校正量得到动校正后的转换波共检波点道集;
根据动校正后的转换波共检波点道集得到炮点剩余静校正后的转换波数据;
根据所述的炮点剩余静校正后的转换波数据得到炮点和检波点静校正后的转换波数据。
2.根据权利要求1所述的方法,其特征是,对所述的三分量地震数据进行预处理,得到预处理后的转换波数据具体包括:
对采集的三分量地震数据进行观测系统加载;
对观测系统加载后的三分量地震数据进行面波衰减;
对面波衰减后的三分量地震数据进行规则噪声衰减;
对规则噪声衰减后的三分量地震数据进行强能量压制;
对强能量压制后的三分量地震数据进行能量补偿;
对能量补偿后的三分量地震数据进行振幅恢复,得到预处理后的转换波数据。
3.根据权利要求1所述的方法,其特征是,根据所述预处理后的转换波数据以及所述的纵波野外炮点静校正量得到动校正后的转换波共检波点道集具体包括:
将预处理后的转换波数据应用于所述的纵波野外炮点静校正量,得到炮点校正后的转换波数据;
对所述炮点校正后的转换波数据进行共中心点道集分选,得到转换波速度;
对所述炮点校正后的转换波数据进行共检波点道集分选,得到分选后的转换波共检波点道集数据;
根据所述的转换波速度对分选后的转换波共检波点道集数据进行动校正,得到动校正后的转换波共检波点道集。
4.根据权利要求1或3所述的方法,其特征是,根据动校正后的转换波共检波点道集得到炮点剩余静校正后的转换波数据具体包括:
对动校正后的转换波共检波点道集进行互相关处理,得到综合相关时差;
对所述的综合相关时差进行地表一致性处理,得到转换波炮点剩余静校正量;
将所述的转换波炮点剩余静校正量应用于所述的转换波共检波点道集上,得到炮点剩余静校正后的转换波数据。
5.根据权利要求4所述的方法,其特征是,根据所述的炮点剩余静校正后的转换波数据得到炮点和检波点静校正后的转换波数据具体包括:
将所述的炮点剩余静校正后的转换波数据分选成共炮点道集;
对所述的共炮点道集进行互相关处理,得到所述共炮点道集对应的综合相关时差;
对所述共炮点道集对应的综合相关时差进行地表一致性处理,得到转换波检点校正量;
将所述的转换波检点校正量应用于所述炮点剩余静校正后的转换波数据,得到炮点和检波点静校正后的转换波数据。
6.一种转换波的静校正设备,其特征是,所述的设备具体包括:
三分量地震数据采集装置,用于采集三分量地震数据;
预处理装置,用于对所述的三分量地震数据进行预处理,得到预处理后的转换波数据;
野外炮点静校正量获取装置,用于获取纵波野外炮点静校正量;
转换波共检波点道集生成装置,用于根据所述预处理后的转换波数据以及所述的纵波野外炮点静校正量得到动校正后的转换波共检波点道集;
第一转换波数据生成装置,用于根据动校正后的转换波共检波点道集得到炮点剩余静校正后的转换波数据;
第二转换波数据生成装置,用于根据所述的炮点剩余静校正后的转换波数据得到炮点和检波点静校正后的转换波数据。
7.根据权利要求6所述的设备,其特征是,所述的预处理装置具体包括:
加载模块,用于对采集的三分量地震数据进行观测系统加载;
面波衰减模块,用于对观测系统加载后的三分量地震数据进行面波衰减;
噪声衰减模块,用于对面波衰减后的三分量地震数据进行规则噪声衰减;
强能量压制模块,用于对规则噪声衰减后的三分量地震数据进行强能量压制;
能量补偿模块,用于对强能量压制后的三分量地震数据进行能量补偿;
振幅恢复模块,用于对能量补偿后的三分量地震数据进行振幅恢复,得到预处理后的转换波数据。
8.根据权利要求6所述的设备,其特征是,所述的转换波共检波点道集生成装置具体包括:
第一应用模块,用于将预处理后的转换波数据应用于所述的纵波野外炮点静校正量,得到炮点校正后的转换波数据;
共中心点道集分选模块,用于对所述炮点校正后的转换波数据进行共中心点道集分选,得到转换波速度;
共检波点道集分选模块,用于对所述炮点校正后的转换波数据进行共检波点道集分选,得到分选后的转换波共检波点道集数据;
动校正模块,用于根据所述的转换波速度对分选后的转换波共检波点道集数据进行动校正,得到动校正后的转换波共检波点道集。
9.根据权利要求6或8所述的设备,其特征是,所述的第一转换波数据生成装置具体包括:
第一互相关处理模块,用于对动校正后的转换波共检波点道集进行互相关处理,得到综合相关时差;
第一一致性处理模块,用于对所述的综合相关时差进行地表一致性处理,得到转换波炮点剩余静校正量;
第二应用模块,用于将所述的转换波炮点剩余静校正量应用于所述的转换波共检波点道集上,得到炮点剩余静校正后的转换波数据。
10.根据权利要求9所述的设备,其特征是,所述的第二转换波数据生成装置具体包括:
分选模块,用于将所述的炮点剩余静校正后的转换波数据分选成共炮点道集;
第二相关处理模块,用于对所述的共炮点道集进行互相关处理,得到所述共炮点道集对应的综合相关时差;
第二一致性处理模块,用于对所述共炮点道集对应的综合相关时差进行地表一致性处理,得到转换波检点校正量;
第三应用模块,用于将所述的转换波检点校正量应用于所述炮点剩余静校正后的转换波数据,得到炮点和检波点静校正后的转换波数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310594371.2A CN104656139B (zh) | 2013-11-21 | 2013-11-21 | 一种转换波的静校正方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310594371.2A CN104656139B (zh) | 2013-11-21 | 2013-11-21 | 一种转换波的静校正方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104656139A true CN104656139A (zh) | 2015-05-27 |
CN104656139B CN104656139B (zh) | 2018-04-06 |
Family
ID=53247486
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310594371.2A Active CN104656139B (zh) | 2013-11-21 | 2013-11-21 | 一种转换波的静校正方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104656139B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105866839A (zh) * | 2016-06-28 | 2016-08-17 | 中国矿业大学(北京) | 基于共姿态道集的静校正方法及装置 |
CN106855639A (zh) * | 2016-11-22 | 2017-06-16 | 中国石油天然气集团公司 | 静校正方法和装置 |
CN107179548A (zh) * | 2017-07-10 | 2017-09-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于真地表的叠前地震成像方法 |
CN107229069A (zh) * | 2016-03-24 | 2017-10-03 | 中国石油化工股份有限公司 | 一种对共检波点数据道集进行速度分析的方法 |
CN109581503A (zh) * | 2017-09-28 | 2019-04-05 | 中国石油化工股份有限公司 | 一种解决沙漠探区层析静校正边界问题的方法及系统 |
CN110879413A (zh) * | 2018-09-05 | 2020-03-13 | 中国石油化工股份有限公司 | 一种射线参数域转换波静校正方法及系统 |
CN111352157A (zh) * | 2018-12-20 | 2020-06-30 | 中国石油天然气集团有限公司 | 一种横波静校正方法及系统 |
CN112824939A (zh) * | 2019-11-20 | 2021-05-21 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据的剩余静校正方法、装置、设备及可读存储介质 |
CN112946745A (zh) * | 2019-12-11 | 2021-06-11 | 中国石油天然气集团有限公司 | 一种基于引力拟合的转换波剩余静校正方法及系统 |
CN113093279A (zh) * | 2020-01-08 | 2021-07-09 | 中国石油天然气集团有限公司 | 转换波静校正方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4611311A (en) * | 1983-04-20 | 1986-09-09 | Chevron Research Company | Method of seismic processing involving converted P- or S-wave data |
CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
-
2013
- 2013-11-21 CN CN201310594371.2A patent/CN104656139B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4611311A (en) * | 1983-04-20 | 1986-09-09 | Chevron Research Company | Method of seismic processing involving converted P- or S-wave data |
CN1773310A (zh) * | 2004-11-09 | 2006-05-17 | 中国石油天然气集团公司 | 转换波延迟时静校正方法 |
Non-Patent Citations (5)
Title |
---|
张旭东 等: "几种转换波静校正方法讨论", 《岩性油气藏》 * |
张旭东: "转换波静校正方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
李彦鹏 等: "利用陆上三分量数据改善气云区构造成像", 《石油地球物理勘探》 * |
李彦鹏 等: "厚风化层覆盖区转换波静校正方法", 《地球物理学报》 * |
潘树林 等: "共炮(检)点剩余静校正方法", 《石油地球物理勘探》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107229069B (zh) * | 2016-03-24 | 2019-04-02 | 中国石油化工股份有限公司 | 一种对共检波点数据道集进行速度分析的方法 |
CN107229069A (zh) * | 2016-03-24 | 2017-10-03 | 中国石油化工股份有限公司 | 一种对共检波点数据道集进行速度分析的方法 |
CN105866839A (zh) * | 2016-06-28 | 2016-08-17 | 中国矿业大学(北京) | 基于共姿态道集的静校正方法及装置 |
CN106855639B (zh) * | 2016-11-22 | 2018-09-04 | 中国石油天然气集团公司 | 静校正方法和装置 |
CN106855639A (zh) * | 2016-11-22 | 2017-06-16 | 中国石油天然气集团公司 | 静校正方法和装置 |
CN107179548A (zh) * | 2017-07-10 | 2017-09-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于真地表的叠前地震成像方法 |
CN107179548B (zh) * | 2017-07-10 | 2019-03-29 | 中国石油集团东方地球物理勘探有限责任公司 | 一种基于真地表的叠前地震成像方法 |
CN109581503A (zh) * | 2017-09-28 | 2019-04-05 | 中国石油化工股份有限公司 | 一种解决沙漠探区层析静校正边界问题的方法及系统 |
CN110879413A (zh) * | 2018-09-05 | 2020-03-13 | 中国石油化工股份有限公司 | 一种射线参数域转换波静校正方法及系统 |
CN111352157A (zh) * | 2018-12-20 | 2020-06-30 | 中国石油天然气集团有限公司 | 一种横波静校正方法及系统 |
CN112824939A (zh) * | 2019-11-20 | 2021-05-21 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据的剩余静校正方法、装置、设备及可读存储介质 |
CN112824939B (zh) * | 2019-11-20 | 2022-08-26 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据的剩余静校正方法、装置、设备及可读存储介质 |
CN112946745A (zh) * | 2019-12-11 | 2021-06-11 | 中国石油天然气集团有限公司 | 一种基于引力拟合的转换波剩余静校正方法及系统 |
CN113093279A (zh) * | 2020-01-08 | 2021-07-09 | 中国石油天然气集团有限公司 | 转换波静校正方法及装置 |
CN113093279B (zh) * | 2020-01-08 | 2024-05-28 | 中国石油天然气集团有限公司 | 转换波静校正方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104656139B (zh) | 2018-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104656139A (zh) | 一种转换波的静校正方法及设备 | |
Grigoli et al. | Automated seismic event location by travel‐time stacking: An application to mining induced seismicity | |
CN104375188B (zh) | 一种地震波透射衰减补偿方法及装置 | |
Rentsch et al. | Fast location of seismicity: A migration-type approach with application to hydraulic-fracturing data | |
CN107918156B (zh) | 检测海底节点采集地震数据极性的方法及装置 | |
CN104199103B (zh) | 转换波静校正方法及装置 | |
CN109001813B (zh) | 一种压制多次波的方法、装置及系统 | |
CN107193045B (zh) | 一种地震数据处理方法及装置 | |
CN108828668B (zh) | 一种叠前时间偏移数据处理方法及装置 | |
CN106054251A (zh) | 一种初至波拾取方法及装置 | |
CN110879412A (zh) | 地下横波速度反演方法、装置、计算设备及存储介质 | |
CN105954797A (zh) | 地震数据的断裂识别方法和装置 | |
CN103210323A (zh) | 处理地震数据的系统和方法 | |
CN106772593A (zh) | 绕射波的成像方法及装置 | |
Kuge | Source modeling using strong-motion waveforms: toward automated determination of earthquake fault planes and moment-release distributions | |
CN108562936B (zh) | 一种裂缝预测方法、系统、存储介质及终端 | |
CN105093300A (zh) | 一种地质体边界识别方法及装置 | |
CN104216016A (zh) | 一种纵波约束扫描的转换波动校正方法以及系统 | |
CN111812716A (zh) | 页岩气储层总有机碳含量叠前定量预测方法、装置及设备 | |
Matsuzawa et al. | Love‐and Rayleigh‐wave microseisms excited by migrating ocean swells in the North Atlantic detected in Japan and Germany | |
CN107422375B (zh) | 海底节点的共反射点道集的确定方法和装置 | |
CN112630832B (zh) | 基于衰减因子随入射角变化的含气性预测方法和装置 | |
CN112255679A (zh) | 地震资料绕射深度偏移处理方法及装置 | |
CN103698812B (zh) | 利用叠前地震道集计算地层品质因数的方法及装置 | |
CN113296157B (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 |