本发明属于地震资料解释技术领域，更具体地讲，涉及一种基于波形指示反演的滩体预测方法。
背景技术：
滩体是有利的油气储集体，通过地震数据可以对滩体的规模、空间展布特征进行描述，其准确划分对于油气勘探意义重大。
传统的滩体预测大多采用定性的人工识别方法，随着科学技术水平的发展和地震资料采集技术的不断提高，地震数据中包含的信息更加丰富，地震属性分析方法被不断地应用于滩体预测中。滩体在地震剖面上多表现为地层增厚、滩体底波峰反射等特征，现阶段对滩体的识别主要为定性识别，主要通过对地震剖面波形特征、地震波峰属性等地震信息进行分析，达到滩体定性识别的效果。但是，目前的预测方法只能通过波形特征和波峰振幅属性等地震属性切片，定性地划分出滩体横向上的边界，对于多个滩体纵向叠置的情况，无法有效地识别滩体边界，因此解释结果存在精度不高、多解性强等问题。
技术实现要素：
针对现有技术中存在的问题，本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如，本发明的目的在于解决传统的定性划分滩体横向边界时，对于多个滩体纵向叠置的情况，无法有效地识别滩体边界，导致解释结果精度不高、多解性强的问题。
为了实现上述目的，本发明提供了一种基于波形指示反演的滩体预测方法，所述预测方法包括：获取地震偏移数据、地震层位数据、钻井数据和测井数据；获得准确的时深关系，将深度域测井曲线转换为时间域测井曲线；获取敏感时间域测井曲线，通过所述敏感时间域测井曲线获得初始框架模型；确认适合滩体预测的反演参数，通过所述地震偏移数据和初始框架模型获得反演成果数据体；通过所述敏感时间域测井曲线获得滩体对应的响应值域，根据所述反演成果数据体和滩体对应的响应值域，提取出所述反演成果数据体中代表滩体的数据，获得滩体预测成果，其中，所述获得初始框架模型包括：对所述敏感时间域测井曲线进行标准化处理；以目的层段顶、底地震层位为界面，根据地震反射特征和/或沉积规律确定地层接触关系。
在本发明的一个示例性实施例中，所述时深关系可以通过井震标定获得。
在本发明的一个示例性实施例中，所述获得初始框架模型步骤中还可以包括在标准化处理前对敏感时间域测井曲线中的异常值进行剔除。
在本发明的一个示例性实施例中，所述标准化处理可以采用直方图分析法。
在本发明的一个示例性实施例中，所述敏感时间域测井曲线的确定可以通过直方图分析法和/或交汇图分析法。
在本发明的一个示例性实施例中，所述反演参数可以包括波阻抗曲线、有效样本数、平滑系数和截止频率。
在本发明的一个示例性实施例中，所述有效样本数可以为3～5。
在本发明的一个示例性实施例中，所述平滑系数可以为1～3。
在本发明的一个示例性实施例中，所述截止频率中低频部分的高通频率可以为5～8hz，高截频率可以为8～10hz，高频部分的低截频率可以为40～60hz，低通频率可以为60～80hz，高通频率可以为150～300hz，高截频率可以为250～400hz。
与现有技术相比，本发明的有益效果包括：基于波形指示反演的滩体预测方法能够识别多个滩体纵向叠置情况下的滩体边界；能够提高解释结果的精度；能够降低多解性；能够更加准确的预测滩体；能够对于储层和油气预测起到重要的指导意义。
附图说明
图1示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的一个流程图。
图2示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的另一个流程图。
图3示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的预测剖面图。
图4示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的预测平面图。
具体实施方式
在下文中，将结合示例性实施例和附图来详细说明本发明的一种基于波形指示反演的滩体预测方法。
图1示出了根据本发明的示例性实施例的基于波形指示反演的滩体预测方法的流程图。如图1所示，在一个示例性实施例中，本发明的基于波形指示反演的滩体预测方法可通过如下步骤实现：
步骤(1)，获取地震偏移数据、地震层位数据、测井数据和钻井数据。例如，获取处理好的四川盆地某三维工区地震偏移数据，作为波形指示反演的基础数据；获取目的层段(例如，龙王庙组)的顶、底地震层位，作为波形指示反演的约束时窗；获取该三维工区内的钻井数据和测井数据，作为井-震标定的基础数据。该三维工区的地震偏移数据可以通过野外数据采集和叠前时间偏移处理得到。
步骤(2)，通过井-震标定方法获得准确的时深关系，将深度域测井曲线转换为时间域测井曲线。步骤(1)中获取的测井数据为深度域，需要将深度域测井数据转换为时间域测井数据用于波形指示反演。
步骤(3)，分析得到敏感时间域测井曲线，通过敏感时间域测井曲线获得初始框架模型。对该三维工区中的时间域测井曲线进行分析，确定出存在异常值的时间域测井曲线，剔除时间域测井曲线中的异常值，然后对正常的时间域测井曲线进行标准化处理，例如，可以采用直方图分析法进行标准化处理；例如，可以根据时间域测井曲线中的声波速度、密度、伽马、孔隙度、波阻抗等，选取目的层段样点，根据实际钻井情况，将样点分为滩体和非滩体两类，可以通过直方图分析方法和/或交汇图分析方法，对目的层段时间域测井曲线进行分析，确定出能够区分出滩体和非滩体的时间域测井曲线，得到敏感时间域测井曲线，例如，针对该三维工区的敏感时间域测井曲线可以为波阻抗曲线；可以目的层段顶、底地震层位为界面，可以根据地震反射特征和/或沉积规律确定地层接触关系，可以井点敏感时间域测井曲线为基础，进行全区线性内插，计算得到初始框架模型。
步骤(4)，分析得到适合滩体预测的反演参数，通过地震偏移数据和初始框架模型获得反演成果数据体。通过试验得出适合该三维工区滩体预测的反演参数，例如，反演参数可以包括输入曲线、有效样本数、平滑系数和截止频率，其中，输入曲线可以为波阻抗曲线，有效样本数可以为3～5，优选地，有效样本数为5；平滑系数可以为1～3，优选地，平滑系数为1；截止频率中低频部分的高通频率可以为5～8hz，高截频率可以为8～10hz，高频部分的低截频率可以为40～60hz，低通频率可以为60～80hz，高通频率可以为150～300hz，高截频率可以为250～400hz，优选地，最佳截止频率中低频部分的高通频率可以为5hz，高截频率可以为10hz，高频部分的低截频率可以为60hz，低通频率可以为80hz，高通频率可以为200hz，高截频率可以为300hz。通过波形指示反演获得反演成果数据体。
步骤(5)，通过敏感时间域测井曲线获得滩体对应的响应值域，根据步骤(4)获得的反演成果数据体和滩体对应的响应值域，提取出反演成果数据体中代表滩体的数据，获得滩体预测成果，用于储层含油气性的预测，如图3、图4所示。如图3所示为龙王庙组滩体预测剖面图，图中的a、b、c区域分别表示预测出的三个滩体，虚线表示滩体的边界。如图4所示为采用本方法对四川盆地某三维工区龙王庙组开展的滩体预测效果图，图中磨溪以及宝龙加数字代表井名，钻井钻遇滩体的井均位于滩体发育区，滩体预测符合率100％，滩体预测成果对于储层与含油气性解释有很好的定性指导作用。
如图2所示，在本发明的另一个示例性实施例中，基于波形指示反演的滩体预测方法可由以下步骤实现：
步骤(1)，获取地震偏移数据、地震层位数据、钻井数据和测井数据。首先，以本领域公知的方法进行野外勘探采集，得到原始采集的数据，经过本领域公知的叠前时间偏移处理方法进行地震资料处理，得到地震偏移数据，作为波形指示反演的基础数据；然后，获取目的层段的顶、底地震层位，作为波形指示反演的约束时窗；最后，获得工区内钻、测井数据，作为井-震标定的基础数据。例如，应用本方法可对四川盆地某三维工区龙王庙组开展滩体预测。
步骤(2)，井-震标定。通过本领域公知的井-震标定方法，对目的层段滩体进行标定，从而确定滩体的地震响应特征，并为下一步波形指示反演提供准确的时深关系。
步骤(3)，波形指示反演。地震有效频带一般为10～60hz，无法直接获取高频成分，但地震波形的横向变化反映了沉积环境的变化，而相似的沉积环境具有可类比的波形特征，因此可以充分利用地震波形的横向变化开展高频成分估计。波形指示反演和传统统计学反演最大的区别就在于统计样本的筛选，传统方法是基于空间域变差函数的，只能粗略表达空间变异程度，无法体现相变特征，而且对样本分布均匀程度要求比较苛刻，限制该方法的应用领域。而波形指示反演利用沉积学基本原理，充分利用地震波形的横向变化来反应储层空间的相变特征，进而分析储层垂向岩性组合高频结构特征，更好地体现了相控的思想，是一种真正的井震结合高频模拟方法，使反演结果从完全随机到逐步确定，同时对井位分布的均匀性没有严格要求。反演步骤如下：
①数据预处理
a、对工区内井数据进行分析，找到存在异常值的测井曲线，剔除曲线中的异常值，并存储为新曲线名备用。
b、可以采用本领域公知的直方图分析法对测井曲线进行标准化处理。
②敏感参数分析
可以采用测井曲线中的声波速度、密度、伽马、孔隙度、波阻抗等曲线，选取目的层段样点，根据钻井实际情况，将样点分为滩体和非滩体两类，可以通过直方图、交汇图分析方法，对目的层段测井曲线进行分析，找出能区分出滩体和非滩体的测井曲线，例如，在本工区中波阻抗曲线能够区分出滩体和非滩体。
③框架模型建立
以目的层段顶、底地震层位为界面，可以依据地震反射特征和沉积规律设置地层接触关系，可以井点波阻抗曲线值为基础，进行全区线性内插，计算出初始框架模型。
④地震反演
通过试验，可以确定适合本工区滩体预测的反演参数，完成波形指示反演，获得反演成果数据体。反演参数可以包括波阻抗曲线、有效样本数、平滑系数和截止频率，其中，输入曲线可以为波阻抗曲线，有效样本数可以为3～5，优选地，有效样本数为5；平滑系数可以为1～3，优选地，平滑系数为1；截止频率中低频部分的高通频率可以为5～8hz，高截频率可以为8～10hz，高频部分的低截频率可以为40～60hz，低通频率可以为60～80hz，高通频率可以为150～300hz，高截频率可以为250～400hz，优选地，最佳截止频率中低频部分的高通频率可以为5hz，高截频率可以为10hz，高频部分的低截频率可以为60hz，低通频率可以为80hz，高通频率可以为200hz，高截频率可以为300hz。通过波形指示反演获得反演成果数据体。
步骤(4)，滩体预测。根据测井曲线敏感参数分析结果，以滩体对应响应值域为提取范围，提取出反演成果数据体中代表滩体的数据，在反演剖面上将滩体顶、底界对比出来，如图3所示为龙王庙组滩体预测剖面图，图中a、b、c区域分别表示预测出的三个滩体，虚线表示滩体的边界。
如图4所示为采用本方法对四川盆地某三维工区龙王庙组开展的滩体预测效果图，图中磨溪以及宝龙加数字代表井名，钻井钻遇滩体的井均位于滩体发育区，滩体预测符合率100％，滩体预测成果对于储层与含油气性解释有很好的定性指导作用。
综上所述，本发明提出了一种基于波形指示反演的滩体预测方法，以地震偏移数据为输入，采用波形指示反演技术，获得目的层段波阻抗反演数据体。根据实际钻井资料和井-震标定，获得钻井已钻遇滩体的波阻抗值域范围，从根据波形指示反演得到的波阻抗反演数据体中提取出滩体值域范围值的数据，可以有效地实现对滩体的预测，实践表明该方法对于该工区的储层和油气预测具有重要的指导意义。
尽管上面已经结合附图和示例性实施例描述了本发明，但是本领域普通技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述实施例进行各种修改。