一种砂泥岩薄互层层位的精细标定方法
【技术领域】
[0001]本发明有关一种石油勘探和开发地质领域中的滩坝砂岩的分析标定方法，特别是指一种以岩相划分为基础，结合测井数据与地震数据的对应关系，在对地质分层及地震层位特征进行分析的基础上，通过不同类型子波的应用和不同标定方法的研宄，对砂泥岩薄互层层位进行精细标定的方法。
【背景技术】
[0002]目前，岩相、测井和地震资料的综合解释和分析是油气勘探和开发过程中最基本、最重要的一种综合分析手段，以岩相划分为基础，结合测井数据与地震数据的对应关系，对砂泥岩薄互层进行精确标定，是精确进行构造解释的关键因素，对薄储层的识别和描述也是一种挑战。
[0003]目前，针对砂泥岩薄互层精细标定的方法很多，针对反射系数序列改进的方法主要有:拟声波曲线重构方法、测井曲线由深度域转换到时间域重采样的方法等；针对子波改进的方法主要有:子波参数调整方法、AGC (自动增益控制)窗口调整方法和多时窗子波提取方法等。以上两类方法从合成记录形成的两个参数出发，有一定的合理性，却忽略了根据子波形成的特点，从标定过程进行研宄，尤其对于砂泥岩薄互层层位来说，标定准确度不高，影响油气勘探的进展。
【发明内容】
[0004]有鉴于此，本发明的主要目的在于提供一种以岩相划分为基础，结合测井数据与地震数据的对应关系，对砂泥岩薄互层层位进行精细标定的方法。
[0005]为达到上述目的，本发明提供一种砂泥岩薄互层层位的精细标定方法，该方法包括:
步骤一、测井分层和测井相特征的识别和描述；
步骤二、地震、测井数据的加载:建立数据库，分别向计算机、工作站软件系统中加载三维高精度地震数据和整理好的测井数据；
步骤三、判断地震剖面极性；
步骤四、识别和描述标志层的地震响应特征及研宄区地震相特征；
步骤五、利用标志层的地震剖面特征与子波提取对某一目的层及全井段完成精细标定;
步骤六、根据振幅与反射系数序列包络的关系，对标定结果进行验证，达到砂泥岩薄互层层位精细标定的目的。
[0006]所述步骤一具体包括:在岩相及测井相模式的指导下，根据岩性及测井曲线特征，确定标志层及砂泥岩薄互层的测井分层深度、测井相特征。
[0007]所述步骤三具体包括:根据研宄区地质特点，在无特殊岩性体的工区可根据地震子波极性判断方法，在有特殊岩性体(如基岩顶面、火成岩顶面等)的工区可根据特殊岩性体地震响应与地震子波极性判断相结合的方法，来进行剖面极性判断。
[0008]在所述步骤四中具体包括，在岩相及测井相模式的指导下，在地震剖面上，识别并分析标志层的地震剖面特征及目的层沉积体系的地震相特征。
[0009]所述步骤五包括:利用时不变子波，运用“二次标定法则”对目的层完成精细标定，及对“长井段、多目的层”的井位完成全井段的精细标定。
[0010]所述利用时不变子波，运用“二次标定法则”对目的层完成精细标定具体包括:
(1)提取统计性子波:根据统计性子波频率分析图，由统计性子波确定井旁道地震数据的主频fm;
(2)提取主频为fm的雷克子波:根据统计性子波确定的主频，提取该频率的雷克子波，并在地震相模式的指导下，选用该主频的雷克子波进行标志层的“初次标定”，达到合成记录与地震剖面主要标志层的吻合；
(3)选用确定性子波进行层位的“二次标定”:运用“二次标定法则”进行薄互层层位精细标定，在此过程中仅需在初次标定成果的基础上进行微调。
[0011]所述对“长井段、多目的层”的井位完成全井段的精细标定包括:运用时变子波，在二次标定的基础之上进行分段调整，从而达到全井段的精细标定。
[0012]所述步骤六具体包括:
典型地震相特征的指导作用:根据研宄区沉积相类型的地震相特征，对合成记录的标定过程进行指导，与标志层的地震响应特征结合，确定人工合成记录的漂移方向和大小，从而达到层位的正确标定；
根据反射系数序列与振幅的对应关系，对层位精细标定结果进行验证:根据反射系数序列与子波的极性，对合成记录的标定结果进行检验。
[0013]本发明是一种新的对砂泥岩薄互层层位进行精细标定的方法，该方法以岩相划分为基础，结合测井信息与地震剖面的对应关系，以地震相为指导，采用多种子波交替综合应用来解决砂泥岩薄互层层位精细标定的问题，为薄互层层位精细构造解释和储层预测提供技术支持。
【附图说明】
[0014]图1为本发明砂泥岩薄互层层位的精细标定方法的流程图。
【具体实施方式】
[0015]为便于对本发明的方法及达到的效果有进一步的了解，现结合附图并举较佳实施例详细说明如下。
[0016]本文利用子波的不同特点，对砂泥岩薄互层层位精细标定过程进行研宄，来提高砂泥岩薄互层层位精细标定的精度，其关键技术包括:①地震剖面极性判断技术；②不同类型子波的提取与应用技术“二次标定技术”;④地震相特征的指导作用；⑤反射系数序列分析技术。以上述技术手段为基础，在实现砂泥岩薄互层层位精细标定的过程中，关键是在剖面极性确定的前提下，在标志层的框架内，如何利用地震相模式的指导作用，实现人工合成记录和地震剖面的精确匹配，从而达到对地震响应特征较弱的砂泥岩薄互层进行精确识别的目的。
[0017]R.E.Sheriff (1980)指出:“极性是表示一个波阻抗界面增加或减少的反射系数(或正或负)的符号。”在地震资料分辨率的范围之内，强波阻抗是否能产生强振幅值与剖面极性相关。因此，在对砂泥岩薄互层层位精细标定方法进行探讨之前，首先需要对剖面进行极性判断。用于地震剖面极性判断的方法很多，根据实例区块的特点及方法的适用性，本发明拟采用根据特殊岩性体的地震响应和根据地震子波的极性两种方法，进行剖面极性的判断。
[0018]根据特殊岩性体进行剖面极性判断的方法如下:首先根据区域地质分析或钻、测井地质资料获取地层界面信息，然后根据反射系数界面的正、负与地震剖面响应的对应关系，来确定地震剖面的极性。如果特殊岩性体(如基岩顶面、火成岩顶面等)的反射系数界面为正，对应波峰，在地震剖面上表现为单轨强峰反射，则表明为正极性剖面，反之，则为负极性剖面；如大套油页岩顶面，则表现为双轨强峰反射，则为正极性剖面，反之，则为负极性剖面。这样，就可以根据这些特殊岩性体的地震响应，来确定地震剖面为正极性剖面还是负极性剖面。
[0019]根据确定性子波正负进行剖面极性判断是另外一种较为常用的地震剖面极性判别方法。利用反射系数序列和井旁道信息提取的确定性子波的极性可确定地震剖面的极性。如子波极性为正，则为正极性剖面；反之，则为负极性剖面。
[0020]子波提取原理不同，其应用也有所不同。根据子波时窗的变化与否，分为时不变子波和时变子波，其中时不变子波除标准雷克子波之外，还包括统计性子波和确定性子波。统计性子波在假定子波的自相关系数与提取时窗段的地震道自相关相同的前提下，给定输入子波的相位及提取时窗和地震道范围提出，该类型子波形状标准，能反映井旁地震道的频率信息；雷克子波为标准零相位子波，根据统计性子波确定的主频，用该频率的雷克子波可进行层位初次标定；确定性子波既用到了测井曲线得到的反射系数序列，又考虑到井旁地震道信息，与测井数据和地震资料均具有较强的相关性，可用于层位的精细标定(即二次标定)。“ 二次标定法则”主要用于单个目的层标定，如果井段较长且存在多目的层，就可采用时变子波，该子波主要通过对地震数据加窗实现。
[0021 ] 地震相指的是沉积相的地震响应，沉积类型不同，其地震响应特征也有所不同。本发明实施例的研宄区为滩坝砂岩沉积，滩坝砂是典型的薄互层砂体，地震剖面特征为“短轴状、中强振幅、不连续”，且该研宄目的层在油页岩集中段的下部，因此，对砂泥岩薄互层有明显的指示作用。坝砂厚度较大，根据坝砂的地震相特征，可以确定漂移的方向，通过对坝砂地震相特征进行准确标定，达到对薄互层精细标定的目的。
[0022]反射系数序列对层位标定的结果起到验证作用。反射系数序列由测井资料获得，具体过程如下，先由声波测井资料可知波在岩石中的传播速度，再根据Gardnar公式计算出密度值，最终得到反射系数序列。反射系数序列具有较高的采样率，如果子波为正极性子波，则反射系数序列的包络与地震剖面的波峰相对应，由此达到对层位标定结果进行验证的目的。
[0023]本发明综合运用以上方法技术，根据砂泥岩薄互层的特点，采用特定技术流程，实现砂泥岩薄互层层位的精细标定，具体步骤流程如下:
步骤一、测井分层和测井相特征的识别和描述:在岩相及测井相模式的指导下，根据岩性及测井曲线特征，确定标志层及砂泥岩薄互层的测井分层深度、测井相特征； 步骤二、地震、测井数据的加载:建立数据库，分别向计算机、工作站相应软件系统中加载三维高精度地震数据和整理好的测井数据；
步骤三、判断地震剖面极性:根据研宄区地质特点，特殊岩性体(如基岩、火成岩顶面、油页岩集中段等)的工区，可根据特殊岩性体地震响应与地震子波极性判断相结合的方法，来进行剖面极性判断，对于其他无特殊岩性体的工区，可根据地震子波极性判断方法，；步骤四、识别和描述标志层的地震响应特征及研宄区(砂泥岩薄互层)地震相特征:在岩相及测井相模式的指导下，在地震剖面上，识别并分析标志层的地震剖面特征及目的层沉积体系的地震相特征；
步骤五、利用标志层的地震剖面特征与子波提取对某一目的层及全井段完成精细标定;
步骤六、根据振幅与反射系数序列包络的关系，对标定结果进行验证，达到砂泥岩薄互层层位精细标定的目的。