一种识别薄储层隐蔽岩性油气藏的方法
【技术领域】
[0001] 本发明涉及一种识别薄储层隐蔽岩性油气藏的方法，属于石油勘探开发技术领 域。
【背景技术】
[0002] 油气勘探中，对隐蔽岩性油气藏的勘探是发展进程的必然结果，特别是薄储层岩 性油气藏，很难进行储层预测。地震资料品质直接制约着储层预测，必须采取一系列的处理 手段得到高保真、高分辨率的地震资料，然后采用合适的解释方法与流程，提高隐蔽岩性油 气藏的识别效果。
[0003] 专利《基于生烃期古构造的油藏有利区的识别方法》（【申请号】201310064704. 0) 提出通过分析生烃期古凸起（斜坡）构造或古油气的运聚动力来寻找油藏有利区，该专利 主要是单一的通过地质手段进行综合研究，未与地震资料结合进行分析。专利《薄层岩性 储层识别与水平钻井跟踪技术》（【申请号】201310618364. 1)公开了利用高产井反射特征 追踪有利砂体，刻画薄层岩性砂体并实现水平井钻井跟踪技术，该专利只是从高产井的地 震反射特征出发，未考虑到其它区域的适用性，具有一定局限性。专利《河道砂岩性油藏的 识别方法》（【申请号】201410844281.9)提出了河道砂体识别、描述技术，该专利没有对地 震资料的处理过程，未考虑到地震资料对工区的适用性。专利《一种基于演化进程的曲流 河砂体储层建筑结构分析方法》（【申请号】201410211934X)公开一种基于演化进程的曲流 河砂体储层建筑结构分析方法，该专利单一的通过地质手段对储层进行研究，未与地震资 料结合进行分析。专利《利用测井约束波阻抗反演预测砂体厚度的方法和装置》（【申请号】 2013103468502)提出通过测井约束波阻抗反演可提高地质条件复杂地区砂体厚度预测的 可靠性，该专利主要是建立模型预测井之间储层厚度的变化，不能进一步在平面上反映储 层展布。
【发明内容】
[0004] 本发明的目的是提供一种识别薄储层隐蔽岩性油气藏的方法，以解决目前薄储层 隐蔽岩性油气藏识别存在的上述问题。
[0005] 本发明为解决上述技术问题提供了一种识别薄储层隐蔽岩性油气藏的方法，该方 法包括以下步骤：
[0006] 1)计算地层的品质因子，根据该品质因子利用反Q滤波对叠后地震资料进行二次 高频补偿处理，结合测井资料选取地震资料主频，以得到理想地震资料；
[0007] 2)以得到的理想地震资料对目的层进行地震响应特征研究，明确储层波形特征， 依据储层顶底面波形反射特征制作地震资料骨架剖面；
[0008] 3)在地震资料骨架剖面基础上精细刻画砂体的平面分布，利用地震属性分析确定 砂体边界，以此确定薄储层隐蔽性油藏。
[0009] 所述步骤1)中地层品质因子是基于VSP下行波记录采用质心偏移算法计算得到。
[0010] 所述步骤1)中地震主频的选取的原则为：所选取的主频数据体与合成记录匹配； 与已钻井的钻井资料情况吻合。
[0011] 所述步骤2)是通过地震正演模拟确定储层地震响应特征，根据所确定的地震响 应特征对地震数据体进行运算，抽取代表储层展布的地震波形骨架信息。
[0012] 所述步骤3)中砂体的精细刻画采用地震线描法，包括以下过程：
[0013] 以砂层组为单元，根据波形分析识别砂体，用砂体顶面解释线条表示砂体发育分 布；
[0014] 进行逐道解释，刻画砂体的平面分布。
[0015] 所述步骤3)中地震属性分析是针对实际钻探油气藏情况，根据地震的瞬时振幅 属性、瞬时相位属性、瞬时频率属性，优选振幅属性平面图，并对振幅分布平面图所指示的 意义进行综合分析，识别出研究区内已存在及可能存在油气藏的区域。
[0016] 本发明的有益效果是：本发明通过计算地层的品质因子，根据该品质因子利用反 Q滤波对叠后地震资料进行二次高频补偿处理，以得到理想地震资料；以得到的理想地震 资料对目的层进行地震响应特征研究，明确储层波形特征，依据储层顶底面波形反射特征 制作地震资料骨架剖面；在地震资料骨架剖面基础上精细刻画砂体的平面分布，利用地震 属性分析确定砂体边界，以此确定薄储层隐蔽性油藏。并通过具体实例进行了验证，本发明 能够可靠的识别隐蔽岩性油气藏，提高薄储层隐蔽岩性油气藏升温识别率。
【附图说明】
[0017] 图1零偏VSP直达波记录；
[0018] 图2-a是直达波的主频随深度变化图；
[0019] 图2-b是品质因子的理论值与反演结果示意图；
[0020] 图3是本发明实施例中春27井地层Q值曲线图；
[0021] 图4是本发明实施例中春33井地层Q值曲线图；
[0022] 图5是本发明实施例中反Q滤波前后（南部）剖面及目的层频谱图；
[0023] 图6是本发明实施例中反Q滤波前后（中部）剖面及目的层频谱图；
[0024] 图7是本发明实施例中反Q滤波前后（北部）剖面及目的层频谱图；
[0025] 图8是本发明实施例中不同主频地震资料时间剖面图；
[0026] 图9是本发明实施例中春29井合成记录示意图；
[0027] 图10是本发明实施例中主频50HZ时过春33-春27地震剖面图；
[0028] 图11是本发明实施例中主频80HZ时过春33-春27地震剖面图；
[0029] 图12是本发明实施例中地震剖面线描解释示意图；
[0030] 图13是本发明实施例中春光油田沙湾组2砂组3小层振幅属性平面图；
[0031] 图14是本发明实施例中春光油田沙湾组2砂组3小层瞬时频率属性平面图；
[0032] 图15是本发明实施例中春光油田沙湾组2砂组3小层瞬时相位属性平面图。
【具体实施方式】
[0033] 下面结合附图对本发明的【具体实施方式】做进一步的说明。
[0034] 本发明借助已有的地震、测井资料，运用井控处理提高分辨率技术、地震线描法和 属性优选法识别油气藏的综合研究方法。该方法首先采用井控处理技术，该技术能最大程 度的利用已有的测井或VSP等资料，将井点数据和地面地震数据进行匹配，从而生成一系 列的匹配属性（包括基于频率的可信度、可预测度、测井资料的固有子波和地震子波的传 递函数、零相位算子等等），并用这些属性来检验地震资料与井资料的匹配程度以达到最佳 匹配，最终便可得到高保真、高分辨率的地震资料；然后以地震响应特征为基础制作地震资 料骨架剖面，在明确储层地震波形特征的基础上，对地震数据体进行运算，抽取代表储层展 布的地震波形骨架信息，地震资料骨架剖面叠合在原始地震资料上能较好的指导对储层顶 底界面反射同相轴的追踪解释，有助于地震线描的准确开展；再对地震资料进行精细解释， 刻画砂体展布范围；最后利用地震属性优选技术识别油气藏。
[0035] 下面以某工区的三维地震资料为例进行说明，该工区三维地震资料的面积 1023km2,以往地震资料主频较低，为50HZ，地震资料品质不好，同相轴相对较少，单砂体在 地震剖面上无法分辨，很难满足储层预测的需求。本发明基于叠后资料进行反Q滤波实验， 对地震资料进行提高分辨率处理，结合测井资料等，将地震资料主频提至80HZ，然后对地震 资料采取地震线描法进行精细解释，优选地震属性，提取地震属性平面图，以识别油气藏可 能发育的位置。该方法的具体过程如下：
[0036] 1.计算地层的品质因子Q，置后反Q滤波处理
[0037] 通过叠前地表一致性反褶积处理后，地震资料中各种频率成分的能量得到了一定 的均衡补偿，但由于高低频反射成份的信噪比不同，多次覆盖叠加的结果势必导致优势频 率向低频方向滑动，因此需要对地震资料进行叠后反Q滤波处理，该步骤首先需要计算地 层的品质因子，然后再通过反Q滤波技术对叠后资料进行二次高频补偿处理。
[0038] 目前常用的计算衰减属性的方法按时间域和频率域分为两类，一类是时间域方 法，一类是频率域方法，其中时间域方法有振幅衰减法、上升时间法、子波模拟法和解析信 号法，频率域方法有谱比法、匹配技术和谱模拟，本实施例中采用比较稳健的质心频率偏移 算法计算品质因子。
[0039] 为了验证本实施例中所采用的质心频率偏移法的可行性，将其应用于一维层状零 偏VSP模型，该模型中震源位于井口，检波器均匀布设在井口，则每个检波器接收到的信号 路径是一条直线，可将上一个检波器接收到的信号作为下一个检波器的输入信号，接收点 的质心频率为fR
[0040]
[0041 ]
[0042] ai为相邻两个检波器之间的地层平均吸收因子，Δf1=f「fi+1为相邻两道质心 频率之差，AΖι为道间距，<为第i道频谱的方差，σ〗是震源频谱的方差。
[0043] 一维零偏VSP模型有六层，深h= 660米，共布置了 45个检波器，道间距为15米， 其中第2层和第5层为薄层，分别厚30米和20米，各层的厚度和品种因子理论值和反演结 果如表1所示。首先，利用一维水平层状粘弹介质VSP正演模拟方法，模拟零偏移距的直达 波记录，如图1所示，所使用的子波为峰值频率40Hz的雷克子波，然后在频率域中求取每一 道的主频和方差。图1是每道的主