一种碎屑岩等效孔隙结构反演方法及系统技术方案

本发明专利技术提供了一种碎屑岩等效孔隙结构反演方法及系统，属于岩石物理学及地球物理测井领域。该碎屑岩等效孔隙结构反演方法利用目的层段的多个岩芯样品的等效孔隙结构参数和测井数据中的纵波速度、横波速度、孔隙度、泥质含量、含水饱和度反演出该目的层段的等效孔隙结构参数。本发明专利技术所需要的参数少，只需要纵、横波速度、孔隙度、泥质含量和含水饱和度，就可以预测出岩石等效孔隙结构，避免了传统的经验公式中因统计关系而存在的不确定性，应用范围广，具有计算速度快、稳定性好的优点，可以为工区内测井横波速度预测提供重要参数。内测井横波速度预测提供重要参数。内测井横波速度预测提供重要参数。

【技术实现步骤摘要】
一种碎屑岩等效孔隙结构反演方法及系统


[0001]本专利技术属于岩石物理学及地球物理测井领域，具体涉及一种碎屑岩等效孔隙结构反演方法及系统。

技术介绍

[0002]流体是直接赋存在岩石孔隙中的，所以基于地震资料的孔隙结构研究在油气/流体识别方面有时比弹性三参数(纵、横波速度和岩石密度)反演更直接、有效。但目前对于从测井资料或者从地震资料预测岩石孔隙结构的研究较少。
[0003]叠前AVO分析技术和叠前反演技术在储层预测中发挥了重要作用，这些技术方法均需要准确的纵波及横波速度测井曲线，然而在实际测井资料中往往缺乏横波速度信息，这给储层预测工作带来了一定的困难。为此，许多地球物理工作者对测井横波速度预测进行了大量的研究工作，发展了经验公式法和理论模型法两大类测井横波速度预测方法。其中经验公式法基于统计规律，其应用效果依赖于统计样品，受地域性影响，不具普遍性。理论模型法基于岩石物理方程，由于具有较高的预测精度，取得了广泛的应用。其中Xu-White速度模型法(或称为计算流程)预测横波速度精度较高，但是由于该模型需要输入的参数较多，有些参数对预本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种碎屑岩等效孔隙结构反演方法，其特征在于：所述方法利用目的层段的多个岩芯样品的等效孔隙结构参数和测井数据中的纵波速度、横波速度、孔隙度、泥质含量、含水饱和度反演出该目的层段的等效孔隙结构参数。2.根据权利要求1所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述方法包括：(1)根据岩芯样品测试数据统计等效孔隙结构参数的数值分布范围及分布规律；(2)利用所述数值分布范围及分布规律获得随机模拟等效孔隙结构参数，并计算得到纵波速度和横波速度；(3)利用计算得到的纵波速度和横波速度以及测井数据中的纵波速度和横波速度建立误差函数f(α
s0
)；(4)判断f(α
s0
)≤ε是否成立，如果是，则转入步骤(5)，如果否，则返回步骤(2)；(5)将所述随机模拟等效孔隙结构参数作为等效孔隙结构参数的反演结果进行输出。3.根据权利要求2所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述步骤(1)的操作包括：输入测井数据，所述测井数据包括：纵波速度、横波速度、孔隙度、泥质含量、含水饱和度；统计岩芯样品的等效孔隙结构参数的数值分布范围及分布规律，将所述数值分布范围表示为[P
l
,P
u
]，其中，P
l
为岩芯样品的等效孔隙结构参数中的最小值，P
u
为岩芯样品的等效孔隙结构参数中的最大值。4.根据权利要求3所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述步骤(2)的操作包括：按照所述分布规律在数值分布范围[P
l
,P
u
]内随机生成一个随机模拟等效孔隙结构参数α
s0
；利用所述随机模拟等效孔隙结构参数α
s0
和孔隙度、泥质含量、含水饱和度，由岩石物理模型计算得到纵波速度V
p
(α
s0
)、横波速度V
s
(α
s0
)。5.根据权利要求4所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述岩石物理模型采用的是适用于碎屑岩的岩石物理模型。6.根据权利要求2所述的碎屑岩等效孔隙结构反演方法，其特征在于：所述步骤(3)的操作包括：利用公式(1)建立误差函数f(α
s0
)：f(α
s0
)＝||V
p
(α
s0
)-V
pm
||+||V
s
(α
s0
)-V
sm
||+|V
p
(α
s0
)/V
s
(α
s0
)-V
pm
/V
sm
|
ꢀꢀꢀ
(1)公式(1)中，V
pm
为测井数据中的纵波速度，V
sm
为测井数据中的横波速度。7.一种实现权利要求1-6任一项所述的碎屑岩等效孔隙结构反演方法的系统，其特征在于：所述系统包括：测井数据输入单元：用于输入测井数据，测井数据包括：纵波速度、横波速度、孔隙度、泥质含量、含水饱和度；统...

【专利技术属性】
技术研发人员：白俊雨，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：