一种基于多波资料的地层切片方法
【技术领域】
[0001]本发明涉及石油地球物理勘探领域，尤指一种基于多波资料的地层切片方法。
【背景技术】
[0002]目前，地震沉积学是继地震地层学和层序地层学之后出现的一个新的学科。近年来国内外学者开展了广泛的研究。概括地讲，地震沉积学是一门主要利用地震信息和技术研究有关沉积岩及其形成过程的学科。具体来讲，它侧重于研究地震岩石学、地震地貌学、沉积结构、沉积史、沉积体系和沉积相平面展布(曾洪流，朱筱敏，朱如凯等，陆相坳陷型盆地地震沉积学研究规范，石油勘探与开发，2012，39 (3) =275-284 ;董春梅，张宪国，林承焰，有关地震沉积学若干问题的探讨，石油地球物理勘探，2006，41 (4) =405-409 ;刘化清，卫平生，李相博等，对地震沉积学理论创新的认识，岩性油气藏，2012，24 (1):7-11)。其主要的研究方法为:90度相位转换、分频解释和地层切片技术。90度相位转换的目的是相位转换后同相轴具有了岩性意义，从而使地震道与岩性测井更加具有一致性。基于地震数据的频率控制了同相轴倾角和内部反射结构这一理念，分频解释可以适当调整地震数据的频率确保解释的等时性。利用地层切片技术进行沉积时间界面的地震地貌学研究。
[0003]现行的所有的地层切片方法都是基于90度相移数据及叠后反演数据。之所以选择这两种数据，就是因为这两种数据可能具备区分岩性的功能，但实际上在复杂岩性区、致密油气区等这两种数据并不能有效区分岩性，那么基于这两种数据的地层切片进行地震沉积学研究就失去了意义。因此，现有的地层切片方法还存在一定的不足。
【发明内容】
[0004]为解决基于90度相移数据及叠后反演数据的地层切片不适用于复杂岩性区、致密油气区等地区的问题，本发明提供一种基于多波资料的地层切片方法，首次将多波资料应用到地震沉积学研究中，可以有效克服现有技术的不足，提高地层切片在不同地区的适用性。
[0005]本发明提出的基于多波资料的地层切片方法包括:步骤1，对包含PP波及PS波的地震数据进行保幅处理，获取叠前道集的PP波分偏移距数据和PS波分偏移距数据；步骤2，联合全波列测井资料，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体；步骤3，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层；步骤4，将所述弹性参数体进行转换获得地质年代体；步骤5，根据所述地质年代体，在所述等时标志层之间以等比例内插的方式对地层进行切片处理，获得具有近似年代地层意义的地层切片，通过多张所述地层切片的平面相解释，获得反映砂组尺度的沉积特征，恢复沉积演化过程。
[0006]进一步的，在步骤1中，至少获取3个以上的叠前道集的PP波分偏移距数据和PS波分偏移距数据。
[0007]进一步的，步骤2，联合全波列测井资料，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体，包括:利用所述全波列测井资料分别对所述PP波分偏移距数据和所述PS波分偏移距数据进行标定，将所述PS波分偏移距数据对齐到所述PP波分偏移距数据的时间域，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体。
[0008]进一步的，步骤3，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层，包括:对所述PP波分偏移距数据进行追踪解释，获取大套标志层；对所述大套标志层进行小波分频，对小波分频后的数据进行地震倾角估算，获得高、低频的倾角；当高、低频的倾角差值低于阈值时，判定地震反射轴是等时的，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层。
[0009]进一步的，步骤4，将所述弹性参数体进行转换获得地质年代体，包括:将所述弹性参数体进行Wheeler域变换，获得地质年代体。
[0010]通过本发明提出的基于多波资料的地层切片方法，可以在地震PP+PS分偏移距资料的基础上，进行叠前PP+PS联合同时反演方法求取高分辨率高精度的纵波阻抗、横波阻抗、泊松比等弹性参数体，创新性地在弹性参数体上进行地层切片研究，可以解决复杂岩性区、致密油气区等地区利用90度相移数据或者叠后反演数据的地层切片不能区分岩性等问题，并将多波资料获得的纵波阻抗、横波阻抗、泊松比等弹性参数体应用到地震沉积学研究中，可以提高地震沉积学技术在复杂岩性区、致密油气区等不同地区的适用性。
【附图说明】
[0011]此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中:
[0012]图1为本发明一实施例的基于多波资料的地层切片方法流程图。
[0013]图2为某三维工区Line753多波叠前反演获得的弹性参数剖面。
[0014]图3为横波阻抗弹性参数体典型地层切片。
[0015]图4为泊松比弹性参数体典型地层切片。
【具体实施方式】
[0016]以下配合图示及本发明的较佳实施例，进一步阐述本发明为达成预定发明目的所采取的技术手段。
[0017]图1为本发明一实施例的基于多波资料的地层切片方法流程图。如图1所示，该方法包括:
[0018]步骤1，对包含PP波及PS波的地震数据进行保幅处理，获取叠前道集的PP (纵波)波分偏移距数据和PS(转换波)波分偏移距数据；
[0019]步骤2，联合全波列测井资料，根据所述PP波分偏移距数据和所述PS波分偏移距数据联合同时反演方法获得至少包括纵波阻抗、横波阻抗、泊松比在内的弹性参数体；
[0020]步骤3，在等时性分析的控制下，对所述PP波分偏移距数据追踪获取等时标志层；
[0021]步骤4，将所述弹性参数体进行转换获得地质年代体；
[0022]步骤5，根据所述地质年代体，在所述等时标志层之间以等比例内插的方式对地层进行切片处理，获得具有近似年代地层意义的地层切片，通过多张所述地层切片的平面相解释，获得反映砂组尺度的沉积特征，恢复沉积演化过程。
[0023]具体而言，在步骤1中，获取保幅处理的η (η彡3)个叠前道集的ΡΡ波分偏移距数据和η (η多3)个PS波分偏移距数据。
[0024]在步骤2中，利用所述全波列测井资料分别对所述ΡΡ波分偏移距数据和所