基于三维地震数据的砂泥岩孔隙类型反演方法

文档序号:8471565阅读:238来源:国知局
基于三维地震数据的砂泥岩孔隙类型反演方法
【技术领域】
[0001] 本发明专利涉及测量技术,特别涉及基于三维地震数据的砂泥岩孔隙类型反演方 法。
【背景技术】
[0002] 在现有的技术中,只有基于测井资料的孔隙扁率反演方法,主要是针对单井的岩 石等效孔隙纵横比进行反演,该方法只能获得井点处的孔隙结构信息,不能有效地获得反 映研宄区孔隙结构信息的三维数据。

【发明内容】

[0003] 本发明的目的在于提供基于三维地震数据的砂泥岩孔隙类型反演方法,以有效地 获得反映研宄区孔隙结构信息的三维数据。
[0004] 本发明提供一种基于三维地震数据的砂泥岩孔隙类型反演方法,包括:设定当前 采样点的线号、道号和时间采样点号,提取当前采样点的孔隙度Φ'、基质矿物的等效体积 模量K' ^、基质矿物的等效剪切模量μ' ^、砂岩的百分比含量v's、泥岩的百分比含量或泥 质含量V'。、密度P '和纵波速度t s;设定当前采样点的砂岩孔隙扁率a ' s的初始值公 式、搜索范围和搜索间隔;并设定每个采样点的泥岩孔隙的扁率α '。的初始值、搜索范围 和搜索间隔;设定完当前采样点的砂岩孔隙扁率〇、和泥岩孔隙扁率α'。的初始值后,计 算孔隙度为Φ'时的干岩石剪切模量μ' ty;计算当前采样点的横波速度V' s >进行当 前采样点的计算的两种横波速度V' s (l和V' s的误差分析,并保存当前砂岩孔隙的扁率 a ' s和泥岩孔隙的扁率α '。的值所对应的误差值σ。
[0005] 优选的,还包括:修改当前采样点的砂岩孔隙扁率α ' 3和泥岩孔隙扁率α '。的 值,直到搜索范围内所有的扁率值计算完毕;寻找砂岩孔隙的扁率α ' 3和泥岩孔隙的扁 率α'。的值所对应的最小误差值〇,并输出对应的砂岩孔隙的扁率a' s、泥岩孔隙的扁 率α '。和φ s'、Φ。'作为当前采样点的计算结果。
[0006] 优选的,还包括:修改当前采样点的线号、道号和时间采样点号,直到所有的采样 点计算完毕,并输出砂岩孔隙的扁率a s及其百分比含量Φ s、泥岩孔隙的扁率α。及其百分 比含量Φ。的三维数据体;对a s、φ3、α。、Φ。进行三维中值滤波,去除异常值的影响,并作 为最终的结果输出。
[0007] 优选的,还包括:输入砂泥岩研宄区的测井参数:纵波速度\、横波速度Vs、孔隙度 Φ和密度P ;统计得到\与P的关系式,拟合得到系数α和β ;再利用双波法统计φ与 \和V 3的关系式,拟合得到系数A、B、C。
[0008] 本发明的有益效果在于:基于砂泥岩研宄区的三维地震资料和测井资料,计算反 映砂泥岩研宄区的孔隙类型的三维数据体,获得地下每个采样点的孔隙类型信息(包括组 成孔隙的各种孔隙类型的扁率及其百分比含量),从而有效地描述砂泥岩储层的孔隙类型, 为砂泥岩储层预测和钻探提供重要的参考数据。
【附图说明】
[0009] 图1是基于三维地震数据的砂泥岩孔隙类型反演方法的流程图。
【具体实施方式】
[0010] 为使本发明的目的、技术方案和优点更加清楚明白,下文中将结合实施例进行详 细说明。
[0011] 如图1所示,本发明基于三维地震数据的砂泥岩孔隙类型反演方法具体包括如下 步骤:
[0012] 步骤(1),输入砂泥岩研宄区的测井参数:纵波速度(Vp)、横波速度(v s)、孔隙度 (Φ)和密度(P)。利用Gardner公式分别统计得到\与P的关系式(见公式4-1),拟合 得到系数α和β ;再利用双波法统计φ与VjPV 3的关系式(见公式4-2),拟合得到系数 A、B、C ;
【主权项】
1. 一种基于三维地震数据的砂泥岩孔隙类型反演方法,其特征在于,包括: 设定当前采样点的线号、道号和时间采样点号,提取当前采样点的孔隙度巾'、基质矿 物的等效体积模量K'^、基质矿物的等效剪切模量y'^、砂岩的百分比含量v's、泥岩的百 分比含量或泥质含量v'。、密度P'和纵波速度V' s; 设定当前采样点的砂岩孔隙扁率a's的初始值公式、搜索范围和搜索间隔;并设定 每个采样点的泥岩孔隙的扁率a'。的初始值、搜索范围和搜索间隔; 设定完当前采样点的砂岩孔隙扁率a、和泥岩孔隙扁率a'。的初始值后,计算孔隙度 为小'时的干岩石剪切模量U'd,y; 计算当前采样点的横波速度Vs' > 进行当前采样点的计算的两种横波速度V's(l和's的误差分析,并保存当前砂岩孔隙的 扁率a's和泥岩孔隙的扁率a'。的值及其所对应的误差值〇。
2. 如权利要求1所述的基于三维地震数据的砂泥岩孔隙类型反演方法,其特征在于, 还包括: 修改当前采样点的砂岩孔隙扁率a' 3和泥岩孔隙扁率a'。的值,直到搜索范围内 所有的扁率值计算完毕; 寻找砂岩孔隙的扁率a's和泥岩孔隙的扁率a'。的值所对应的最小误差值〇,并 输出对应的砂岩孔隙的扁率a's、泥岩孔隙的扁率a'。和(})'s、巾'。作为当前采样点的计 算结果。
3. 如权利要求2所述的基于三维地震数据的砂泥岩孔隙类型反演方法,其特征在于, 还包括: 修改当前采样点的线号、道号和时间采样点号,直到所有的采样点计算完毕,并输出砂 岩孔隙的扁率a3及其百分比含量巾s、泥岩孔隙的扁率a。及其百分比含量巾。的三维数 据体; 对as、(})s、a。、(}>。进行三维中值滤波,去除异常值的影响,并作为最终的结果输出。
4. 如权利要求3所述的基于三维地震数据的砂泥岩孔隙类型反演方法,其特征在于, 还包括:输入砂泥岩研宄区的测井参数:纵波速度Vp、横波速度Vs、孔隙度巾和密度P;统 计得到\与P的关系式,拟合得到系数a和0 ;再利用双波法统计巾与VjPVs的关系 式,拟合得到系数A、B、C。
5. 如权利要求4所述的基于三维地震数据的砂泥岩孔隙类型反演方法,其特征在于, 还包括:输入砂泥岩研宄区的基于叠后反演得到的伽玛GR三维地震数据,并采计算泥质含 量即泥岩的百分比含量V。的三维地震数据。
6. 如权利要求5所述的基于三维地震数据的砂泥岩孔隙类型反演方法,其特征在于, 还包括:输入基于叠前弹性反演得到纵波阻抗IP、横波阻抗Is的三维地震数据;计算Vp、Vs 和P的三维地震数据;计算孔隙度巾的三维地震数据。
7. 如权利要求6所述的基于三维地震数据的砂泥岩孔隙类型反演方法,其特征在于, 还包括:输入研宄区的砂岩和泥岩的剪切模量,采计算基质矿物的等效体积模量Ktl和等效 剪切模量Utl的三维数据体。
8. 如权利要求7所述的基于三维地震数据的砂泥岩孔隙类型反演方法,其特征在于, 还包括:计算砂岩孔隙在总孔隙中的百分比含量巾 3和泥岩孔隙在总孔隙中的百分比含量 的三维数据体。
【专利摘要】本发明提供一种基于三维地震数据的砂泥岩孔隙类型反演方法,包括:设定当前采样点的线号、道号和时间采样点号,提取当前采样点的孔隙度φ'、基质矿物的等效体积模量K0'、基质矿物的等效剪切模量μ'0、砂岩的百分比含量v's、泥岩的百分比含量即泥质含量v'c、密度ρ'和纵波速度V's;设定当前采样点的砂岩孔隙扁率α's的初始值公式、搜索范围和搜索间隔;并设定每个采样点的泥岩孔隙的扁率α'c的初始值、搜索范围和搜索间隔;设定完当前采样点的砂岩孔隙扁率α's和泥岩孔隙扁率α'c的初始值后,计算孔隙度为φ'时的干岩石剪切模量μ'dry;计算当前采样点的横波速度V's_0;进行当前采样点的计算的两种横波速度V's_0和V's的误差分析,并保存当前砂岩孔隙的扁率α's和泥岩孔隙的扁率α'c的值及其所对应的误差值σ。
【IPC分类】G01V1-30, G01N15-08
【公开号】CN104792684
【申请号】CN201510189404
【发明人】熊晓军, 谭峰, 侯秋平, 李良泉, 黄智勇
【申请人】成都理工大学
【公开日】2015年7月22日
【申请日】2015年4月18日
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1