本发明属于石油地质工程领域，特别涉及一种重力流水道砂体中优质储层的预测方法。
背景技术：
在中国东部陆相盆地中，大面积分布着重力流水道砂体。由于重力流是沉积物-流体混合物在重力作用下整体顺斜坡向下运动、以悬移方式搬运为主的高密度流体流动，其形成属于事件性沉积作用，在洪水、地震、海啸、火山喷发等阵发性因素直接或间接诱发下导致沉积物的滑动形成重力流沉积。重力流水道砂体中碎屑颗粒大小混杂，分选差，储层物性差异大，如何应用地震预测优质储层难度较大。
目前国内外石油地质专家一般采用地震资料中的地震属性预测重力流水道砂体储层的分布，对重力流水道砂体储层进行定性描述和分析。
在实现本发明的过程中，本发明人发现现有技术中至少存在以下问题：
现有的技术只能区分出储层，不能定量的区别优质储层和差储层，同时，目前还没有一种用来预测重力流水道砂体中优质储层的有效方法，钻井作业比较盲目，钻井成功率低。
技术实现要素：
鉴于此，本发明提供一种能够定量预测重力流水道砂体中优质储层的预测方法。
具体而言，包括以下的技术方案：
一种重力流水道砂体中优质储层的预测方法，所述方法包括：
获取多个单井的重力流水道砂体的测井资料，得到每个所述单井的重力流水道砂体对应的孔隙度、岩石密度和纵波速度；
将每个所述单井的所述岩石密度和所述纵波速度相乘，得到每个所述单井的重力流水道砂体对应的波阻抗；
根据多个所述单井的重力流水道砂体对应的所述孔隙度和所述波阻抗，建立孔隙度-波阻抗交互图版，并在该图版上确定重力流水道砂体中的优质储层；
对目标重力流水道砂体储层进行波阻抗反演，并结合所述孔隙度-波阻抗交互图版所确定的所述重力流水道砂体的优质储层信息，得到所述目标重力流水道砂体储层的定量分布。
进一步地，所述根据多个所述单井的重力流水道砂体对应的所述孔隙度和所述波阻抗，建立孔隙度-波阻抗交互图版，并在该图版上确定重力流水道砂体中的优质储层包括：根据所述孔隙度-波阻抗交互图版，分别选用主频对所述单井进行正演，来确定所述重力流水道砂体中的优质储层。
优选地，所选用的所述主频为20Hz、30Hz、40Hz。
进一步地，所述获取单井的重力流水道砂体的测井资料，得到每个所述单井的重力流水道砂体对应的孔隙度、岩石密度和纵波速度中，所述测井资料包括声波时差。
优选地，所述孔隙度的计算方式为：Φ＝(Δt-Δtma)/[(Δtf-Δtma)*Cp]，其中：Φ为孔隙度；Δt为声波时差；Δtma为岩石骨架声波时差值；Δtf为岩石孔隙流体声波时差值；Cp为压实系数。
优选地，所述孔隙度的计算方式为：Φ＝(ρma-ρb)/(ρma-ρf)，其中：Φ为孔隙度；ρma是岩石密度骨架值；ρb是岩石密度测量值；ρf是岩石孔隙流体密度值。
进一步地，所述纵波速度的计算方式为：V＝1×106/△t，其中：V为纵波速度。
进一步地，所述孔隙度-波阻抗交互图版包括：当孔隙度大于12％时，波阻抗值小于9900g.cm-3.m.s-1；当孔隙度小于8％时，波阻抗值大于11600g.cm-3.m.s-1。
优选地，所述对目标重力流水道砂体储层进行波阻抗反演，并结合所述孔隙度-波阻抗交互图版所确定的重力流水道砂体的优质储层信息，得到所述目标重力流水道砂体储层的定量分布包括：
在所述目标重力流水道砂体储层的地震数据体中解释所述目标重力流水道砂体储层的顶底界面；
提取所述顶底界面内的地震属性，定性预测所述目标重力流水道砂体储层中优质储层的分布情况；
根据所述分布情况，进行波阻抗约束反演，并结合所述孔隙度-波阻抗交互图版，定量预测所述目标重力流水道砂体储层中优质储层的分布。
本发明实施例提供的技术方案的有益技术效果如下：
本发明实施例提供了一种能够定量预测重力流水道砂体中优质储层的预测方法。具体地，首先通过运用单井的重力流水道砂体储层的测井资料，得到单井的重力流水道砂体对应的波阻抗，建立孔隙度-波阻抗交互图版，并确定重力流水道砂体中的优质储层，再进行井震结合，通过目标重力流水道砂体储层的地震属性和测井约束反演进行优质储层的定性与定量预测。本发明实施例提供的预测方法可定量预测重力流水道砂体优质储层的分布，克服了钻井的盲目性，提高钻井成功率。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种重力流水道砂体中优质储层的预测方法的方法流程图。
具体实施方式
为使本发明的技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。除非另有定义，本发明实施例所用的所有技术术语均具有与本领域技术人员通常理解的相同的含义。
本发明提供一种重力流水道砂体优质储层的预测方法，参见图1，该方法包括如下步骤：
步骤101：获取多个单井的重力流水道砂体的测井资料，得到每个单井的重力流水道砂体对应的孔隙度Φ、岩石密度和纵波速度V；
根据已钻单井的测井资料(包括声波时差△t等)，求取各单井中重力流水道砂体对应的孔隙度Φ和纵波速度V，优选地，可以用两种方式计算孔隙度Φ：
一种是采用声波时差孔隙度的计算方式：
Φ＝(Δt-Δtma)/[(Δtf-Δtma)*Cp]，其中：Δt为声波时差；Δtma为岩石骨架声波时差值；Δtf为岩石孔隙流体声波时差值；Cp为压实系数。
另一种是采用密度测井孔隙度计算方式：
Φ＝(ρma-ρb)/(ρma-ρf)，其中：ρma是岩石密度骨架值；ρb是岩石密度测量值；ρf是岩石孔隙流体密度值。
纵波速度V＝1×106/Δt，声波时差Δt单位为μs/m。
由此，根据多个单井的重力流水道砂体储层的孔隙度Φ和纵波速度V的值，可建立反映孔隙度与纵波速度的对应关系的地质模型，如下：
孔隙度Φ大于12％的重力流水道砂体对应低速层，纵波速度V＝4000-4350m.s-1；孔隙度Φ小于8％的重力流水道砂体对应高速层，纵波速度V＝5400-5500m.s-1；即优质储层对应低速层，差储层对应高速层。
步骤102：将每个单井的岩石密度和纵波速度相乘，得到每个单井的重力流水道砂体对应的波阻抗；
用岩石密度的测量值ρb与纵波速度V相乘，求取各单井重力流水道砂体对应的波阻抗(impedanee)。
步骤103：根据多个单井的重力流水道砂体对应的孔隙度Φ和波阻抗，建立孔隙度-波阻抗交互图版，并在该图版上确定重力流水道砂体中的优质储层；
根据多个单井的重力流水道砂体储层的孔隙度Φ和波阻抗的值，可建立反映孔隙度与波阻抗的线性关系交互图版，由关系图版可得出，孔隙度Φ大于12％的中高孔隙度区，波阻抗值小于9900g.cm-3.m.s-1；孔隙度Φ小于8％的致密区，波阻抗值大于11600g.cm-3.m.s-1，从而可知，中低波阻抗对应的砂体储层物性较好，为优质储层，高阻抗对应的砂体储层物性相对较差，为致密砂岩发育区。
进一步地，在孔隙度-波阻抗交互图版的基础上，分别选用主频对各单井进行正演，以进一步确定重力流水道砂体中的优质储层。结果表明，在厚度相近的情况下，中波阻抗对应的砂体储层物性较好，高阻抗对应的砂体储层物性相对较差。
且由于在地震采集时，主要目的层的主频一般在20-30Hz之间，而提高主频可以更精确的推测砂体积尖灭点的位置，但是太高了又容易造成假象，因此优选地，用20Hz、30Hz和40Hz进行正演。
步骤104：对目标重力流水道砂体储层进行波阻抗反演，并结合孔隙度-波阻抗交互图版所确定的重力流水道砂体的优质储层信息，得到目标重力流水道砂体储层的定量分布。
具体地，在目标重力流水道砂体储层的地震数据体中解释目标重力流水道砂体储层的顶底界面；
提取顶底界面内的地震属性，定性预测目标重力流水道砂体储层中优质储层的分布情况；
根据分布情况，进行波阻抗约束反演，并结合孔隙度-波阻抗交互图版，定量预测所述目标重力流水道砂体储层中优质储层的分布。
本发明实施例提供的预测方法可定量预测重力流水道砂体优质储层的分布，经5口井钻探实施后，优质储层钻遇率由原来的50％提高到100％，单井日产油由原来的0.01-5.2吨提高到10.5--40.8吨，克服了钻井的盲目性，提高钻井成功率。
以上所述仅是为了便于本领域的技术人员理解本发明的技术方案，并不用以限制本发明。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。