CN104199103A - 转换波静校正方法及装置 - Google Patents
转换波静校正方法及装置 Download PDFInfo
- Publication number
- CN104199103A CN104199103A CN201410403429.5A CN201410403429A CN104199103A CN 104199103 A CN104199103 A CN 104199103A CN 201410403429 A CN201410403429 A CN 201410403429A CN 104199103 A CN104199103 A CN 104199103A
- Authority
- CN
- China
- Prior art keywords
- data
- correction
- road collection
- wave
- evened
- 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
一种转换波静校正方法及装置,该方法包括:对多波地震数据进行预处理;获取野外炮点校正量及检波点校正量;对地震数据进行校正;对校正后的地震数据进行分选及速度分析,生成地震波传播速度;对道集进行动校正及切除操作,生成动校正道集数据;将PP波CMP道集数据转换到PS域,将转换到PS域的PP波CMP道集数据进行叠加,生成叠加地震数据;在叠加地震数据上选择一个目标层,对目标层进行层位拾取,获得层位时间;预设基准面,计算层位时间与预设水平基准面时间的时差;将动校正道集数据校正到预设基准面上,生成构造拉平后的PP波动校正道集数据;将构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的共检波点道集叠加数据。
Description
技术领域
本发明涉及地球物理多波地震勘探技术领域,尤其涉及一种转换波静校正方法及装置。
背景技术
多波地震勘探是利用纵波、横波、转换波等多种信息对油气藏进行精细勘察的一种勘探技术。这种技术可以识别薄互层、小断层、低幅度构造,准确预测储层边界和储层物性、监测流体,甚至还可以直接用于油气预测等。由于转换波地震数据的信噪比较纵波低,其静校正量比纵波大,特别是山区转换波静校正已经成为影响转换波地震勘探效果重要的关键技术之一,是公认的陆上转换波地震勘探的难点,转换波静校正效果已成为衡量转换波地震数据处理技术水平高低的重要标准。
转换波的静校正量是由震源的PP波静校正量和接收点静校正量组成。由于在同一介质中转换波的传播速度比PP波速度低,并且转换波传播不受孔隙流体影响,这就造成了同一位置的转换波静校正量大于PP波静校正量。转换波信噪比通常较低,初至附近纵波折射、转换波折射等信息混杂在一起,使得转换波初至不易分辨,难以拾取,这就造成了折射静校正算法难以在转换波静校正中发挥作用。
目前反射波共检波点叠加转换波静校正方法在实际生产中已经得到了广泛的应用,呈现了很好的实际生产效果。目前使用较多、效果较好的转换波静校正方法是应用PP波和转换波共检波点叠加进行转换波静校正量计算,该方法中,炮点的静校正量可以通过常规的纵波静校正计算获得,应用炮点静校正量后共检波点道集中每个地震道的静校正剩余量相同,均为对应检波点的静校正量。把经过炮点静校正后的共检波点道集数据进行动校正后叠加,每个共检波点道集叠加后形成的叠加道,包涵了对应检波点的转换波静校正量。如果反射界面水平,可以直接把PP波共检波点叠加转换道PS域(不需要拉平,因为反射界面已经水平),通过拾取共检波点叠加剖面上的同一个层位,获得PP波对应的层位时间,同样在PS波共检波点叠加剖面上拾取层位(PP波拾取层位必须与PS波拾取层位相对应),获得PS波共检波点叠加剖面上对应层位的拾取时间,然后把PP波与PS波共检波点层位时间的差作为PS波检波点上的静校正量。这种常规交互转换波静校正方法的前提条件是:1)假设PS炮点上静校正已经完全解决,剩余校正量都是检波点上的静校正量,而且共检波点上剩余校正量相同;2)拾取作为计算PS波检波点静校正的层位是水平的。第2个前提条件对于有较大构造起伏的转换波勘探区域而言,这种方法的静校正效果不理想,主要是CRP道集在较大构造起伏的地质情况下层位不水平,不能进行精确的同相叠加,CRP叠加后有效信号的信噪比较低,没法拾取层位或者拾取层位精度不高。
发明内容
本发明提供一种转换波静校正方法及装置,可以有效提高地下构造起伏大的多波地震数据转换波静校正的精度和效率。
为了实现上述目的,本发明提供一种转换波静校正方法,所述的方法包括:
采集多波地震数据,并对所述多波地震数据进行预处理,生成预处理后的地震数据,所述预处理后的地震数据包括:预处理后的PP波地震数据及预处理后的PS波地震数据;
获取野外炮点校正量及检波点校正量;
应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据;
对所述校正后的地震数据进行分选,并进行速度分析,生成PP和PS地震波传播速度;
利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,所述的动校正道集数据包括:PP波CMP道集数据及PS波CCP道集数据;
将所述PP波CMP道集数据转换到PS域,并将转换到PS域的所述PP波CMP道集数据进行叠加,生成叠加地震数据;
在所述叠加地震数据上选择一个目标层,对所述目标层进行层位拾取,获得层位时间;
预设基准面,并计算所述层位时间与所述预设水平基准面的时间的时差;
根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正道集数据;
将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的共检波点道集叠加数据。
在一实施例中,所述获取野外炮点校正量及检波点校正量,包括:获取PP波的野外炮点校正量和检波点校正量及PS波的野外炮点校正量及检波点校正量。
在一实施例中,应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据,包括:
应用PP波的炮点校正量和检波点校正量,对所述预处理后的PP波数据进行静校正,生成静校正后的PP波数据。
在一实施例中,应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据,包括:
应用PS波的炮点校正量,对所述预处理后的PS波数据进行静校正,生成静校正后的PS波数据。
在一实施例中,对所述校正后的地震数据进行分选,并进行速度分析,生成PP波地震波传播速度,包括:
对所述静校正后的PP波数据进行CMP道集分选,然后进行PP波速度分析,生成PP速度。
在一实施例中,对所述校正后的地震数据进行分选,并进行速度分析,生成PS波地震波传播速度,包括:
对所述静校正后的PS波数据进行CCP道集分选,然后进行PS波速度分析,生成PS速度。
在一实施例中,利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,包括:
利用所述PP速度对所述CMP道集进行动校正及拉伸切除,生成动校正后的PP波CMP道集数据。
在一实施例中,利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,包括:
利用所述PS速度对所述CCP道集进行动校正及拉伸切除,生成动校正后的PS波CCP道集数据。
在一实施例中,根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的波动校正道集数据,包括:
根据所述时差将所述PP波CMP道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正CMP道集数据。
在一实施例中,所述的方法还包括:利用所述时差将所述PS波动校正CCP道集数据进行构造拉平操作,生成构造拉平后的PS波动校正CCP道集数据。
在一实施例中,将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的道集叠加数据,包括:
对PP波动校正CMP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PP波共检波点道集叠加数据。
在一实施例中,将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的道集叠加数据,包括:
对PS波动校正CCP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PS波共检波点道集叠加数据。
为了实现上述目的,本发明提供一种转换波静校正装置,所述的装置包括:
数据预处理单元,用于采集多波地震数据,并对所述多波地震数据进行预处理,生成预处理后的地震数据,所述预处理后的地震数据包括:预处理后的PP波地震数据及预处理后的PS波地震数据;
获取单元,用于获取野外炮点校正量及检波点校正量;
数据校正单元,用于应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据;
速度分析单元,用于对所述校正后的地震数据进行分选,并进行速度分析,生成获得PP和PS地震波传播速度;
数据处理单元,用于利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,所述的动校正道集数据包括:PP波CMP道集数据及PS波CCP道集数据;
数据叠加单元,用于将所述PP波CMP道集数据转换到PS域,并将转换到PS域的所述PP波CMP道集数据进行叠加,生成叠加地震数据;
层位时间生成单元,用于在所述叠加地震数据上选择一个信噪比高于预设值的目标层,对所述目标层进行层位拾取,生成层位时间;
时差计算单元,用于预设基准面,并计算所述层位时间与所述预设水平基准面的时间的时差;
构造拉平单元,用于根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正道集数据,将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的共检波点道集叠加数据。
在一实施例中,所述获取单元具体用于:获取PP波的野外炮点校正量和检波点校正量及PS波的野外炮点校正量及检波点校正量。
在一实施例中,所述的数据校正单元包括:
PP波数据校正模块,用于应用PP波的炮点校正量和检波点校正量,对所述预处理后的PP波数据进行静校正,生成静校正后的PP波数据。
在一实施例中,所述的数据校正单元包括:
PS波数据校正模块,用于应用PS波的炮点校正量,对所述预处理后的PS波数据进行静校正,生成静校正后的PS波数据。
在一实施例中,所述的速度分析单元包括:
PP速度分析模块,用于对所述静校正后的PP波数据进行CMP道集分选,然后进行PP波速度分析,生成PP速度。
在一实施例中,所述的速度分析单元包括:
PS速度分析模块,用于对所述静校正后的PS波数据进行CCP道集分选,然后进行PS波速度分析,生成PS速度。
在一实施例中,所述的数据处理单元包括:
PP速度数据处理模块,用于利用所述PP速度对所述CMP道集进行动校正及拉伸切除,生成动校正后的PP波CMP道集数据。
在一实施例中,所述的数据处理单元包括:
PS速度数据处理模块,用于利用所述PS速度对所述CCP道集进行动校正及拉伸切除,生成动校正后的PS波CCP道集数据。
在一实施例中,所述的构造拉平单元包括:
PP波构造拉平模块,用于根据所述时差将所述PP波CMP道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正CMP道集数据。
在一实施例中,所述的构造拉平单元包括:
PS波构造拉平模块,用于利用所述时差将所述PS波动校正CCP道集数据进行构造拉平操作,生成构造拉平后的PS波动校正CCP道集数据。
在一实施例中,所述的PP波构造拉平模块还用于对PP波动校正CMP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PP波共检波点道集叠加数据。
在一实施例中,所述的PS波构造拉平模块还用于对PS波动校正CCP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PS波共检波点道集叠加数据。
本发明实施例的有益效果在于,本发明能适应地下构造起伏大的转换波静校正情况,可有效提高地下构造起伏大的多波地震资料转换波静校正的精度和效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的转换波静校正方法示意图;
图2为本发明实施例的转换波静校正装置的结构框图;
图3为本发明实施例的数据校正单元203的结构框图;
图4为本发明实施例的速度分析单元204的结构框图;
图5为本发明实施例的数据处理单元205的结构框图;
图6为本发明实施例的构造拉平单元209的结构框图;
图7为本发明实施例转换波静校正构造拉平流程图;
图8为本发明实施例预处理后的PP波CMP道集(未做动校正);
图9为本发明实施例将图2数据动校正后的PP波CMP道集;
图10为本发明实施例将图3数据动转换到PS域的CMP道集;
图11为本发明实施例PP波CMP道集叠加剖面;
图12为本发明实施例对图5叠加数据进行构造拉平;
图13为本发明实施例构造拉平前的PP波CRP道集叠加剖面;
图14为本发明实施例构造拉平后的PP波CRP道集叠加剖面;
图15为本发明实施例构造拉平前的PS波CRP道集叠加剖面;
图16为本发明实施例构造拉平后的PS波CRP道集叠加剖面;
图17为本发明实施例构造拉平后的PP波CRP道集叠加剖面层位拾取示意图;
图18为本发明实施例构造拉平后的PS波CRP道集叠加剖面层位拾取示意图;
图19为本发明实施例转换波静校正后PS波CMP道集叠加结果;
图20现有技术中交互转换波静校正后PS波CMP道集叠加结果。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供一种转换波静校正方法,该转换波静校正方法包括:
步骤101:采集多波地震数据,并对所述多波地震数据进行预处理,生成预处理后的地震数据,所述预处理后的地震数据包括:预处理后的PP波地震数据及预处理后的PS波地震数据;
步骤102:获取野外炮点校正量及检波点校正量;
步骤103:应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据;
步骤104:对所述校正后的地震数据进行分选,并进行速度分析,生成PP和PS地震波传播速度;
步骤105:利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,所述的动校正道集数据包括:PP波CMP道集数据及PS波CCP道集数据;
步骤106:将所述PP波CMP道集数据转换到PS域,并将转换到PS域的所述PP波CMP道集数据进行叠加,生成叠加地震数据;
步骤107:在所述叠加地震数据上选择一个目标层,对所述目标层进行层位拾取,获得层位时间;
步骤108:预设基准面,并计算所述层位时间与所述预设水平基准面的时间的时差;
步骤109:根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正道集数据;
步骤110:将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的共检波点道集叠加数据。
本发明首先野外采集地震数据,并对采集的地震数据进行预处理,应用所述炮点校正量和/或检波点校正量对预处理后的地震数据进行校正、分选、速度分析、动校正及切除操作,生成包括PP波CMP道集数据及PS波CCP道集数据的动校正道集数据;将PP波CMP道集数据转换到PS域,将转换到PS域的所述PP波CMP道集数据进行叠加,生成叠加地震数据;在叠加地震数据上选择一个目标层,对该目标层进行层位拾取获得层位时间;计算层位时间与预设水平基准面的时间的时差,根据该时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正道集数据;将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的共检波点道集叠加数据。通过上述过程,本发明能适应地下构造起伏大的转换波静校正情况,可有效提高地下构造起伏大的多波地震数据转换波静校正的精度和效率。
具体实施时,步骤101中对地震数据进行预处理,包括加载观测系统,地震资料叠前去噪以及振幅恢复等。
预处理后的地震数据包括:预处理后的PP波地震数据及预处理后的PS波地震数据,故步骤102具体实施时,需要获取PP波的野外炮点校正量和检波点校正量及PS波的野外炮点校正量及检波点校正量。
步骤103具体实施时,应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据,包括:应用PP波的炮点校正量和检波点校正量,对所述预处理后的PP波数据进行静校正,生成静校正后的PP波数据;还包括:应用PS波的炮点校正量,对所述预处理后的PS波数据进行静校正,生成静校正后的PS波数据。
步骤104具体实施时,对所述校正后的地震数据进行分选,并进行速度分析,生成PP波地震波传播速度,包括:对所述静校正后的PP波数据进行CMP道集分选,然后进行PP波速度分析,生成PP速度;还包括:对所述静校正后的PS波数据进行CCP道集分选,然后进行PS波速度分析,生成PS速度。
步骤105具体实施时,利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,包括:利用所述PP速度对所述CMP道集进行动校正及拉伸切除,生成动校正后的PP波CMP道集数据;还包括:利用所述PS速度对所述CCP道集进行动校正及拉伸切除,生成动校正后的PS波CCP道集数据。
步骤109具体实施时,根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的波动校正道集数据,包括:根据所述时差将所述PP波CMP道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正CMP道集数据;还包括:利用所述时差将所述PS波动校正CCP道集数据进行构造拉平操作,生成构造拉平后的PS波动校正CCP道集数据。
步骤110具体实施时,将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的道集叠加数据,包括:对PP波动校正CMP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PP波共检波点道集叠加数据;还包括:对PS波动校正CCP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PS波共检波点道集叠加数据。
本发明通过野外采集多波地震数据。对采集的多波数据进行观测系统加载;预处理并计算野外炮点校正量和检波点校正量。对预处理后的PP波数据,应用炮点静校正量和检波点校正量,获得应用炮点校正量后的PP波数据,对静校正后的PP波数据进行共中心点道集(CMP)分选,然后做PP波速度分析,获得PP速度,用PP波速度对PP波CMP道集做动校正,获得动校正后的PP波CMP道集数据,把动校正后的PP波CMP道集数据转换到PS域,然后进行叠加,获得PP波转换到PS域叠加后的地震数据,在这个叠加数据上选择一个信噪比较高的目标层,对目标层进行层位拾取,获得层位时间,对获得的层位时间,选定一个水平基准面时间,计算层位时间到基准面时间的差,利用这个时差把获得的PP波转换到PS域的CMP道集数据校正到这个水平基准面上,这样可以得到构造拉平后的PP波动校正CMP道集在PS域的数据。对PS数据采取与PP波数据处理类似的处理方法,其中PS波只应用炮点校正量,PS波本来就在PS域,不必转换,也不用做CCP叠加和拾取层位(PP波是CMP),可以直接应用PP波到水平基准面的时间把动校正后的PS波构造拉平,得到校正到平到水平基准面上的PS波CCP数据。然后把校正到水平基准面上的数据PP波数据和PS波数据都分选成CRP道集,然后做CRP道集叠加,叠加后的数据就可以做常规的转换波静校正处理了。
常规交互转换波静校正方法是直接在PP波(PS域)和PS波CRP叠加剖面上进行层位拾取(拾取同一层位),然后计算PP波层位时间与PS波层位时间的差,把这个时间差作为PS波对应检波点上的校正量,对PS波动校正前的数据进行校正,就完成了常规交互转换波静校正的计算和应用,达到了消除PS波检波点上静校正的目的。
如图2所示,本发明提供一种转换波静校正装置,该转换波静校正装置包括:数据预处理单元201,获取单元102,数据校正单元203,速度分析单元204,数据处理单元205,数据叠加单元206,层位时间生成单元207,时差计算单元208及构造拉平单元209。
数据预处理单元201用于采集多波地震数据,并对所述多波地震数据进行预处理,生成预处理后的地震数据,所述预处理后的地震数据包括:预处理后的PP波地震数据及预处理后的PS波地震数据。
获取单元202用于获取野外炮点校正量及检波点校正量。
数据校正单元203用于应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据。
速度分析单元204用于对所述校正后的地震数据进行分选,并进行速度分析,生成获得PP和PS地震波传播速度。
数据处理单元205用于利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,所述的动校正道集数据包括:PP波CMP道集数据及PS波CCP道集数据。
数据叠加单元206用于将所述PP波CMP道集数据转换到PS域,并将转换到PS域的所述PP波CMP道集数据进行叠加,生成叠加地震数据。
层位时间生成单元207用于在所述叠加地震数据上选择一个信噪比高于预设值的目标层,对所述目标层进行层位拾取,生成层位时间。
时差计算单元208用于预设基准面,并计算所述层位时间与所述预设水平基准面的时间的时差。
构造拉平单元209用于根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正道集数据,将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的共检波点道集叠加数据。
具体实施时,获取单元102获取PP波的野外炮点校正量和检波点校正量及PS波的野外炮点校正量及检波点校正量。
在一实施例中,如图3所示,数据校正单元203包括:PP波数据校正模块301及PS波数据校正模块302。
PP波数据校正模块301用于应用PP波的炮点校正量和检波点校正量,对所述预处理后的PP波数据进行静校正,生成静校正后的PP波数据。PS波数据校正模块302用于应用PS波的炮点校正量,对所述预处理后的PS波数据进行静校正,生成静校正后的PS波数据。
在一实施例中,如图4所示,速度分析单元204包括:PP速度分析模块401及PS速度分析模块402。PP速度分析模块401可以用于对所述静校正后的PP波数据进行CMP道集分选,然后进行PP波速度分析,生成PP速度。PS速度分析模块402可以用于对所述静校正后的PS波数据进行CCP道集分选,然后进行PS波速度分析,生成PS速度。
在一实施例中,如图5所示,数据处理单元205包括:PP速度数据处理模块501及PS速度数据处理模块502。
PP速度数据处理模块501用于利用所述PP速度对所述CMP道集进行动校正及拉伸切除,生成动校正后的PP波CMP道集数据。PS速度数据处理模块502用于利用所述PS速度对所述CCP道集进行动校正及拉伸切除,生成动校正后的PS波CCP道集数据。
在一实施例中,如图6所示,构造拉平单元209包括:PP波构造拉平模块601及PS波构造拉平模块602。
PP波构造拉平模块601用于根据所述时差将所述PP波CMP道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正CMP道集数据。PS波构造拉平模块602用于利用所述时差将所述PS波动校正CCP道集数据进行构造拉平操作,生成构造拉平后的PS波动校正CCP道集数据。
在另一实施例中,PP波构造拉平模块601还用于对PP波动校正CMP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PP波共检波点道集叠加数据。
在另一实施例中,PS波构造拉平模块602还用于对PS波动校正CCP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PS波共检波点道集叠加数据。本发明实施例的有益效果在于,本发明能适应地下构造起伏大的转换波静校正情况,可有效提高地下构造起伏大的多波地震资料转换波静校正的精度和效率。
本发明实施例的有益效果在于,本发明能适应地下构造起伏大的转换波静校正情况,可有效提高地下构造起伏大的多波地震资料转换波静校正的精度和效率。
为了详细说明本发明,下面结合具体的例子,对本发明的转换波静校正方法及取得的效果进行说明。
对已经进行PP波炮点和检波点校正后的数据,以及对PS波进行炮点校正后的数据,可以用流程图来表示本发明的具体实施方式的过程,如图7所示。
本发明运用于转换波数据静校正,适用于构造起伏的转换波静校正计算。图8是动校正前的PP波CMP道集模型记录示意图,此时已经实施了步骤101至步骤103,完成了炮点和检波点校正量。
对图8中的数据实施步骤104及步骤105,完成PP波速度分析,并进行动校正处理,获得动校正、切除后的数据,如图9。
将图9数据转换到PS域,转换到PS域数据,如图10。
将图10数据做叠加处理,得到叠加后的数据,如图11。
从图11中可以看出这个模型数据构造起伏大,同相轴不水平。本发明在图11上拾取了一个层位,获得层位时间,如图11中虚线所示,选定1000毫秒作为水平基准面。
计算拾取的层位时间到水平基准面的时间差,用这个时差对图11数据进行校正,这样可以把拾取的层位构造拉平到1000毫秒处,拉平构造后的数据如图12所示。构造拉平的主要目的是保证PP波CRP道集的同相轴能同相叠加,提高PP波CRP道集叠加数据信噪比,便于CRP叠加剖面上层位准确拾取。
图13是构造拉平前的CRP道集叠加剖面示意图,箭头部分成像质量很差,层位拾取时不能准确拾取。
图14是构造拉平后的CRP道集叠加剖面示意图,选择的水平基准面为1000毫秒,拉平构造后,CRP道集数据在1000毫秒处的同相轴实现了同相叠加,1000毫秒处同相轴连续性比造拉平前的CRP道集叠加数据好,对比图13,这一点在图14的红色箭头标识的区域可以明显识别出来,这样就达到了本发明的目的。本发明把PP波计算出来的拉平构造的时差应用到PS波炮点校正后的CCP动校正道集上,把PS波与PP波相同层位的构造拉平,然后进行CRP道集分选和叠加。
图15为构造拉平前的PS波CRP道集叠加剖面示意图,图16为造拉平后的PS波CRP道集叠加剖面示意图。对比图16及图15,可以得出与PP波CRP道集构造拉平前后对比相同的结论,1000毫秒处同相轴连续性比造拉平前的CRP道集叠加数据好,信噪比也相对较高,便于常规交付转换波层位正确拾取。
本发明的PP波(PS域)动校正后的CMP道集和PS波动校正后的CCP道集构造拉平后就可以进行后续的常规交付转换波静校正处理,把构造拉平后的PP波和PS波数据分选成CRP道集进行叠加,在CRP道集叠加后的数据上拾取PP波与PS波的同一层位,获得PP波和PS波的两个层位时间。图17为PP波CRP道集叠加剖面和层位拾取示意图,见1000毫秒层位处的水平拾取线,图18为PS波CRP道集叠加剖面和层位拾取示意图,见1000毫秒层位处的弯曲拾取线。用PP波层位时间减去PS波层位时间,就可以得到一个层位时间差,这个层位时间差就是PS波检波点上的校正量,把计算出的PS波校正量应用到PS波进行炮点校正后的数据上,就完成了本发明的一种改善地下构造起伏大的转换波静校正效果的方法。
本发明转换波静校正结果,如图19,常规交付转换波静校正结果如图20。比较图19和图20,图19同相轴的连续性比图20好,信噪比、分辨率也比图20高。达到了本发明的改善地下构造起伏大的转换波静校正的效果。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (24)
1.一种转换波静校正方法,其特征在于,所述的方法包括:
采集多波地震数据,并对所述多波地震数据进行预处理,生成预处理后的地震数据,所述预处理后的地震数据包括:预处理后的PP波地震数据及预处理后的PS波地震数据;
获取野外炮点校正量及检波点校正量;
应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据;
对所述校正后的地震数据进行分选,并进行速度分析,生成PP和PS地震波传播速度;
利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,所述的动校正道集数据包括:PP波CMP道集数据及PS波CCP道集数据;
将所述PP波CMP道集数据转换到PS域,并将转换到PS域的所述PP波CMP道集数据进行叠加,生成叠加地震数据;
在所述叠加地震数据上选择一个目标层,对所述目标层进行层位拾取,获得层位时间;
预设基准面,并计算所述层位时间与所述预设水平基准面的时间的时差;
根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正道集数据;
将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的共检波点道集叠加数据。
2.根据权利要求1所述的转换波静校正方法,其特征在于,所述获取野外炮点校正量及检波点校正量,包括:获取PP波的野外炮点校正量和检波点校正量及PS波的野外炮点校正量及检波点校正量。
3.根据权利要求2所述的转换波静校正方法,其特征在于,应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据,包括:
应用PP波的炮点校正量和检波点校正量,对所述预处理后的PP波数据进行静校正,生成静校正后的PP波数据。
4.根据权利要求2所述的转换波静校正方法,其特征在于,应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据,包括:
应用PS波的炮点校正量,对所述预处理后的PS波数据进行静校正,生成静校正后的PS波数据。
5.根据权利要求3所述的转换波静校正方法,其特征在于,对所述校正后的地震数据进行分选,并进行速度分析,生成PP波地震波传播速度,包括:
对所述静校正后的PP波数据进行CMP道集分选,然后进行PP波速度分析,生成PP速度。
6.根据权利要求4所述的转换波静校正方法,其特征在于,对所述校正后的地震数据进行分选,并进行速度分析,生成PS波地震波传播速度,包括:
对所述静校正后的PS波数据进行CCP道集分选,然后进行PS波速度分析,生成PS速度。
7.根据权利要求5所述的转换波静校正方法,其特征在于,利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,包括:
利用所述PP速度对所述CMP道集进行动校正及拉伸切除,生成动校正后的PP波CMP道集数据。
8.根据权利要求6所述的转换波静校正方法,其特征在于,利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,包括:
利用所述PS速度对所述CCP道集进行动校正及拉伸切除,生成动校正后的PS波CCP道集数据。
9.根据权利要求7所述的转换波静校正方法,其特征在于,根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的波动校正道集数据,包括:
根据所述时差将所述PP波CMP道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正CMP道集数据。
10.根据权利要求8所述的转换波静校正方法,其特征在于,根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的波动校正道集数据,包括:
利用所述时差将所述PS波动校正CCP道集数据进行构造拉平操作,生成构造拉平后的PS波动校正CCP道集数据。
11.根据权利要求9所述的转换波静校正方法,其特征在于,将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的道集叠加数据,包括:
对PP波动校正CMP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PP波共检波点道集叠加数据。
12.根据权利要求10所述的转换波静校正方法,其特征在于,将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的道集叠加数据,包括:
对PS波动校正CCP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PS波共检波点道集叠加数据。
13.一种转换波静校正装置,其特征在于,所述的装置包括:
数据预处理单元,用于采集多波地震数据,并对所述多波地震数据进行预处理,生成预处理后的地震数据,所述预处理后的地震数据包括:预处理后的PP波地震数据及预处理后的PS波地震数据;
获取单元,用于获取野外炮点校正量及检波点校正量;
数据校正单元,用于应用所述炮点校正量和/或检波点校正量对所述预处理后的地震数据进行校正,生成校正后的地震数据;
速度分析单元,用于对所述校正后的地震数据进行分选,并进行速度分析,生成获得PP和PS地震波传播速度;
数据处理单元,用于利用所述地震波传播速度对道集进行动校正及切除操作,生成动校正道集数据,所述的动校正道集数据包括:PP波CMP道集数据及PS波CCP道集数据;
数据叠加单元,用于将所述PP波CMP道集数据转换到PS域,并将转换到PS域的所述PP波CMP道集数据进行叠加,生成叠加地震数据;
层位时间生成单元,用于在所述叠加地震数据上选择一个目标层,对所述目标层进行层位拾取,生成层位时间;
时差计算单元,用于预设基准面,并计算所述层位时间与所述预设水平基准面的时间的时差;
构造拉平单元,用于根据所述时差将所述动校正道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正道集数据,将所述构造拉平后的动校正道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的共检波点道集叠加数据。
14.根据权利要求13所述的转换波静校正装置,其特征在于,所述获取单元具体用于:获取PP波的野外炮点校正量和检波点校正量及PS波的野外炮点校正量及检波点校正量。
15.根据权利要求14所述的转换波静校正装置,其特征在于,所述的数据校正单元包括:
PP波数据校正模块,用于应用PP波的炮点校正量和检波点校正量,对所述预处理后的PP波数据进行静校正,生成静校正后的PP波数据。
16.根据权利要求14所述的转换波静校正装置,其特征在于,所述的数据校正单元包括:
PS波数据校正模块,用于应用PS波的炮点校正量,对所述预处理后的PS波数据进行静校正,生成静校正后的PS波数据。
17.根据权利要求15所述的转换波静校正装置,其特征在于,所述的速度分析单元包括:
PP速度分析模块,用于对所述静校正后的PP波数据进行CMP道集分选,然后进行PP波速度分析,生成PP速度。
18.根据权利要求16所述的转换波静校正装置,其特征在于,所述的速度分析单元包括:
PS速度分析模块,用于对所述静校正后的PS波数据进行CCP道集分选,然后进行PS波速度分析,生成PS速度。
19.根据权利要求17所述的转换波静校正装置,其特征在于,所述的数据处理单元包括:
PP速度数据处理模块,用于利用所述PP速度对所述CMP道集进行动校正及拉伸切除,生成动校正后的PP波CMP道集数据。
20.根据权利要求18所述的转换波静校正装置,其特征在于,所述的数据处理单元包括:
PS速度数据处理模块,用于利用所述PS速度对所述CCP道集进行动校正及拉伸切除,生成动校正后的PS波CCP道集数据。
21.根据权利要求19所述的转换波静校正装置,其特征在于,所述的构造拉平单元包括:
PP波构造拉平模块,用于根据所述时差将所述PP波CMP道集数据校正到所述预设基准面上,并进行构造拉平操作,生成构造拉平后的PP波动校正CMP道集数据。
22.根据权利要求20所述的转换波静校正装置,其特征在于,所述的构造拉平单元包括:
PS波构造拉平模块,用于利用所述时差将所述PS波动校正CCP道集数据进行构造拉平操作,生成构造拉平后的PS波动校正CCP道集数据。
23.根据权利要求21所述的转换波静校正装置,其特征在于,所述的PP波构造拉平模块还用于对PP波动校正CMP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PP波共检波点道集叠加数据。
24.根据权利要求22所述的转换波静校正装置,其特征在于,所述的PS波构造拉平模块还用于对PS波动校正CCP道集数据进行共检波点道集分选及道集叠加,生成构造拉平后的PS波共检波点道集叠加数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410403429.5A CN104199103B (zh) | 2014-08-15 | 2014-08-15 | 转换波静校正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410403429.5A CN104199103B (zh) | 2014-08-15 | 2014-08-15 | 转换波静校正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104199103A true CN104199103A (zh) | 2014-12-10 |
CN104199103B CN104199103B (zh) | 2017-01-25 |
Family
ID=52084415
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410403429.5A Active CN104199103B (zh) | 2014-08-15 | 2014-08-15 | 转换波静校正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104199103B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104635270A (zh) * | 2015-02-16 | 2015-05-20 | 中国地质大学(北京) | 基于合成记录约束的转换波静校正方法 |
CN106324677A (zh) * | 2015-07-08 | 2017-01-11 | 中国石油化工股份有限公司 | 一种位移场剩余动校正方法及装置 |
CN106772601A (zh) * | 2016-12-21 | 2017-05-31 | 成都理工大学 | 一种基于叠前地震数据的裂缝储层快速预测方法 |
CN106855639A (zh) * | 2016-11-22 | 2017-06-16 | 中国石油天然气集团公司 | 静校正方法和装置 |
CN107255833A (zh) * | 2017-07-20 | 2017-10-17 | 中国石油天然气集团公司 | 地震叠加剖面的确定方法和装置 |
CN112649873A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 地震数据分频动校正处理方法及系统 |
CN113138410A (zh) * | 2020-01-16 | 2021-07-20 | 中国石油天然气集团有限公司 | 纵横波速度获取方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102998704A (zh) * | 2012-12-13 | 2013-03-27 | 石颖 | 一种地球物理勘探地震资料处理方法 |
CN103116184A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种表层调查中计算横波静校正量的方法 |
CN103487834A (zh) * | 2013-09-09 | 2014-01-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 转换波共检波点叠加静校正方法 |
CN103630934A (zh) * | 2012-08-23 | 2014-03-12 | 中国石油天然气集团公司 | 一种确定转换波检波点大的横波静校正量的方法 |
-
2014
- 2014-08-15 CN CN201410403429.5A patent/CN104199103B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103116184A (zh) * | 2011-11-16 | 2013-05-22 | 中国石油集团东方地球物理勘探有限责任公司 | 一种表层调查中计算横波静校正量的方法 |
CN103630934A (zh) * | 2012-08-23 | 2014-03-12 | 中国石油天然气集团公司 | 一种确定转换波检波点大的横波静校正量的方法 |
CN102998704A (zh) * | 2012-12-13 | 2013-03-27 | 石颖 | 一种地球物理勘探地震资料处理方法 |
CN103487834A (zh) * | 2013-09-09 | 2014-01-01 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 转换波共检波点叠加静校正方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104635270A (zh) * | 2015-02-16 | 2015-05-20 | 中国地质大学(北京) | 基于合成记录约束的转换波静校正方法 |
CN104635270B (zh) * | 2015-02-16 | 2017-03-08 | 中国地质大学(北京) | 基于合成记录约束的转换波静校正方法 |
CN106324677A (zh) * | 2015-07-08 | 2017-01-11 | 中国石油化工股份有限公司 | 一种位移场剩余动校正方法及装置 |
CN106855639A (zh) * | 2016-11-22 | 2017-06-16 | 中国石油天然气集团公司 | 静校正方法和装置 |
CN106855639B (zh) * | 2016-11-22 | 2018-09-04 | 中国石油天然气集团公司 | 静校正方法和装置 |
CN106772601A (zh) * | 2016-12-21 | 2017-05-31 | 成都理工大学 | 一种基于叠前地震数据的裂缝储层快速预测方法 |
CN106772601B (zh) * | 2016-12-21 | 2018-10-23 | 成都理工大学 | 一种基于叠前地震数据的裂缝储层快速预测方法 |
CN107255833A (zh) * | 2017-07-20 | 2017-10-17 | 中国石油天然气集团公司 | 地震叠加剖面的确定方法和装置 |
CN107255833B (zh) * | 2017-07-20 | 2019-06-28 | 中国石油天然气集团公司 | 地震叠加剖面的确定方法和装置 |
CN112649873A (zh) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | 地震数据分频动校正处理方法及系统 |
CN112649873B (zh) * | 2019-10-10 | 2024-04-09 | 中国石油化工股份有限公司 | 地震数据分频动校正处理方法及系统 |
CN113138410A (zh) * | 2020-01-16 | 2021-07-20 | 中国石油天然气集团有限公司 | 纵横波速度获取方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104199103B (zh) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104199103A (zh) | 转换波静校正方法及装置 | |
CN104181599B (zh) | 一种基于近地表层的折射波静校正处理方法以及系统 | |
CN104536043A (zh) | 一种深度域整体速度模型融合方法及装置 | |
CN105182408A (zh) | 一种合成地震记录的制作方法和装置 | |
CN107765298A (zh) | 一种确定近地表速度模型的方法及装置 | |
CN104570079A (zh) | 一种纵波、转换横波地震资料的时间匹配方法 | |
CN105093319A (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN104483705A (zh) | 一种三维剩余静校正方法 | |
CN107817525B (zh) | 一种近地表常速填充方法及装置 | |
CN102879819A (zh) | 保持地震波场运动学特征的地震数据处理方法及装置 | |
CN103713318A (zh) | 对三维叠后地震数据进行拼接的方法 | |
CN103499836A (zh) | 空变-多时窗融合高精度剩余静校正方法 | |
CN112415591A (zh) | 绕射波的成像方法、装置、电子设备及存储介质 | |
CN109696705A (zh) | 深海obs上下行波场成像数据校正方法及装置 | |
CN103210323A (zh) | 处理地震数据的系统和方法 | |
CN104749623B (zh) | 一种地震资料成像处理方法 | |
CN106054252B (zh) | 一种叠前时间偏移的方法及装置 | |
CN104614762B (zh) | 疏松砂岩气藏边界确定方法及装置 | |
CN106707345B (zh) | 一种角度弹性参数岩性识别方法和装置 | |
CN112764102B (zh) | 地震资料采集方法及装置 | |
CN110376644A (zh) | 叠前地震数据拼接方法及装置 | |
CN107422375B (zh) | 海底节点的共反射点道集的确定方法和装置 | |
CN109581494A (zh) | 叠前偏移方法及装置 | |
CN113933901B (zh) | 深度域地震数据拼接方法及装置 | |
CN105425290A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |