• 浪潮孙丕恕从信息时代到智能时代 人工智能价值将爆发式释放 2019-12-21
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
    • / 18
    • 下载费用:30 金币  

    重庆时时彩源代码: 一种含气饱和度确定方法和装置.pdf

    关 键 词:
    一种 饱和度 确定 方法 装置
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201510452356.3

    申请日:

    2015.07.28

    公开号:

    CN106405638A

    公开日:

    2017.02.15

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01V 1/28申请日:20150728|||公开
    IPC分类号: G01V1/28; G01V1/30 主分类号: G01V1/28
    申请人: 中国石油化工股份有限公司; 中国石油化工股份有限公司石油勘探开发研究院
    发明人: 王瑞; 董宁; 韩磊
    地址: 100728 北京市朝阳区朝阳门北大街22号
    优先权:
    专利代理机构: 北京思创毕升专利事务所 11218 代理人: 刘明华
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201510452356.3

    授权公告号:

    ||||||

    法律状态公告日:

    2018.09.18|||2017.03.15|||2017.02.15

    法律状态类型:

    授权|||实质审查的生效|||公开

    摘要

    本公开提供了一种含气饱和度确定方法和装置,包括:对分角度叠加的数据体进行射线弹性阻抗反演得到多个角度的射线弹性阻抗反演体;根据岩性和孔隙度值进行岩相划分;根据提取的井旁道数据及测井岩相进行计算得到AVO阻抗体;对AVO阻抗体的阻抗值域的划分,计算形成地震岩相体;再次进行射线弹性阻抗反演得到储层弹性参数体;在地震岩相体的控制下,对储层弹性参数体进行散点交汇分析,获得对含气饱和度敏感的弹性参数体;提取井旁道数据,进行压实分析形成干岩石骨架体;根据利用干岩石骨架体计算得到的流体体积??楹土魈迕芏?,定量计算含气饱和度。本公开实现了排除现有叠前反演含气性预测方法的多解性,准确的预测有效储层中的含气饱和度效果。

    权利要求书

    1.一种含气饱和度确定方法,其特征在于,包括:
    对分角度叠加的数据体进行射线弹性阻抗反演得到多个角度的射线弹性阻
    抗反演体;
    根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分,得
    到测井岩相;
    根据提取的井旁道数据及测井岩相进行计算得到AVO阻抗体;
    对AVO阻抗体的阻抗值域的划分,计算形成地震岩相体;
    将多个射线弹性阻抗反演体进行射线弹性阻抗反演得到储层弹性参数体;
    在地震岩相体的控制下,对储层弹性参数体进行散点交汇分析,获得对含
    气饱和度敏感的弹性参数体;
    对敏感弹性参数体提取井旁道数据,进行压实分析形成干岩石骨架体;
    根据利用干岩石骨架体计算得到的流体体积??楹土魈迕芏?,定量计算含
    气饱和度。
    2.根据权利要求1所述的含气饱和度确定方法,其特征在于,所述根据岩
    性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分,得到测井岩相
    包括:
    从多个角度的射线弹性阻抗反演体中提取小角度射线弹性阻抗和大角度射
    线弹性阻抗;
    提取小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性
    信息;
    根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度信息及大角度射线
    弹性阻抗中的岩性信息进行岩相划分,得到测井岩相。
    3.根据权利要求2所述的含气饱和度确定方法,其特征在于,所述根据提
    取的井旁道数据及划分出的测井岩相进行计算得到AVO阻抗体包括:
    根据抽取的井旁道数据及划分出的测井岩相,进行散点交汇分析;
    利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进
    行坐标旋转,得到坐标转化公式;
    根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度
    射线弹性阻抗进行计算得到AVO阻抗体。
    4.根据权利要求1所述的含气饱和度确定方法,其特征在于,所述对敏感
    弹性参数体提取井旁道数据,进行压实分析形成干岩石骨架体包括:
    对敏感弹性参数体提取井旁道数据;
    对井旁道数据进行压实分析;
    根据压实分析的结果获取干岩石骨架的属性变化趋势;
    根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体。
    5.根据权利要求1所述的含气饱和度确定方法,其特征在于,所述对含气
    饱和度敏感的弹性参数体包括密度和体积模量。
    6.一种含气饱和度确定装置,其特征在于,包括:
    第一反演???,用于对分角度叠加的数据体进行射线弹性阻抗反演,得到
    多个角度的射线弹性阻抗反演体;
    岩相划分???,用于根据岩性和孔隙度值对多个角度的射线弹性阻抗反演
    体进行岩相划分,得到测井岩相;
    阻抗体计算???,用于根据提取的井旁道数据及划分出的测井岩相进行计
    算得到AVO阻抗体;
    岩相体计算???,用于对AVO阻抗体的阻抗值域的划分,计算形成地震岩
    相体;
    第二反演???,用于将多个射线弹性阻抗反演体进行射线弹性阻抗反演,
    得到储层弹性参数体;
    交汇分析???,用于在地震岩相体的控制下,对储层弹性参数体进行散点
    交汇分析,获得对含气饱和度敏感的弹性参数体;
    压实分析???,用于对敏感弹性参数体提取井旁道数据,进行压实分析形
    成干岩石骨架体;
    饱和度计算???,用于根据利用干岩石骨架体计算得到的流体体积??楹?br />流体密度,定量计算含气饱和度。
    7.根据权利要求6所述的含气饱和度确定装置,其特征在于,所述岩相划
    分??榘ǎ?br />
    第一提取子???,用于从多个角度的射线弹性阻抗反演体中提取小角度射
    线弹性阻抗和大角度射线弹性阻抗;
    第二提取子???,用于提取小角度射线弹性阻抗中的孔隙度信息及大角度
    射线弹性阻抗中的岩性信息;
    划分子???,用于根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度
    信息及大角度射线弹性阻抗中的岩性信息进行岩相划分,得到测井岩相。
    8.根据权利要求7所述的含气饱和度确定装置,其特征在于,所述阻抗体
    计算??榘ǎ?br />
    第一分析子???,用于根据抽取的井旁道数据及划分出的测井岩相,进行
    散点交汇分析;
    坐标旋转子???,用于利用属性投影的方法对小角度射线弹性阻抗体和大
    角度射线弹性阻抗体进行坐标旋转,得到坐标转化公式;
    计算子???,用于根据散点交汇分析的结果及坐标转化公式对小角度射线
    弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。
    9.根据权利要求6所述的含气饱和度确定装置,其特征在于,所述压实分
    析??榘ǎ?br />
    数据提取子???,用于对敏感弹性参数体提取井旁道数据;
    压实分析子???,用于对井旁道数据进行压实分析;
    获取子???,用于根据压实分析的结果获取干岩石骨架的属性变化趋势;
    插值子???,用于根据干岩石骨架的属性变化趋势进行多属性井插值形成
    干岩石骨架体。
    10.根据权利要求6所述的含气饱和度确定装置,其特征在于,所述对含
    气饱和度敏感的弹性参数体包括密度和体积模量。

    说明书

    一种含气饱和度确定方法和装置

    技术领域

    本发明涉及地球物理勘探地震解释领域领域,更具体地,涉及一种基于岩
    相体控制的含气饱和度确定方法及系统。

    背景技术

    二十世纪90年代以来,我国加到了天然气勘探力度,探明储量每年快速稳
    步增长,但是其探明程度依然很低,仅为6.7%。随着勘探和开发程度的增加,
    储层评价和含气性预测成为工作中的核心部分。如何从地震数据中获得饱和度、
    孔隙度等储层物性参数是进行储层孔隙流体预测的关键。目前,利用地震数据进
    行油气预测基本还是针对储层所含流体的定性预测。

    然而地下储层岩性的复杂性决定了不能用一种属性或是一个岩石物理方程
    来描述地下储层分布,通过相关预测方法得到的流体因子或是反演得到的储层
    物性参数都会受到非有效储层的干扰,预测精度低。

    发明内容

    本公开提出了一种排除非有效储层的干扰、预测精度高的含气饱和度确定
    方法及系统。

    根据本公开的一方面,提出了一种含气饱和度确定方法,包括:

    对分角度叠加的数据体进行射线弹性阻抗反演得到多个角度的射线弹性阻
    抗反演体;

    根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分,得
    到测井岩相;

    根据提取的井旁道数据及测井岩相进行计算得到AVO阻抗体;

    对AVO阻抗体的阻抗值域的划分,计算形成地震岩相体;

    将多个射线弹性阻抗反演体进行射线弹性阻抗反演得到储层弹性参数体;

    在地震岩相体的控制下,对储层弹性参数体进行散点交汇分析,获得对含
    气饱和度敏感的弹性参数体;

    对敏感弹性参数体提取井旁道数据,进行压实分析形成干岩石骨架体;

    根据利用干岩石骨架体计算得到的流体体积??楹土魈迕芏?,定量计算含
    气饱和度。

    根据本公开的第一方面,在本公开的第一种可能的实现方式中,所述根据
    岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相划分,得到测井岩
    相包括:

    从多个角度的射线弹性阻抗反演体中提取小角度射线弹性阻抗和大角度射
    线弹性阻抗;

    提取小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗中的岩性
    信息;

    根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度信息及大角度射线
    弹性阻抗中的岩性信息进行岩相划分,得到测井岩相。

    根据本公开第一方面第一种可能的实现方式,在本公开的第二种可能的实
    现方式中,所述根据提取的井旁道数据及划分出的测井岩相进行计算得到AVO
    阻抗体包括:

    根据抽取的井旁道数据及划分出的测井岩相,进行散点交汇分析;

    利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体进
    行坐标旋转,得到坐标转化公式;

    根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻抗和大角度
    射线弹性阻抗进行计算得到AVO阻抗体。

    根据本公开的第一方面,在本公开的第三种可能的实现方式中,所述对敏
    感弹性参数体提取井旁道数据,进行压实分析形成干岩石骨架体包括:

    对敏感弹性参数体提取井旁道数据;

    对井旁道数据进行压实分析;

    根据压实分析的结果获取干岩石骨架的属性变化趋势;

    根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石骨架体。

    根据本公开的第一方面,所述对含气饱和度敏感的弹性参数体包括密度和
    体积模量。

    根据本公开的第二方面,一种含气饱和度确定装置,包括:

    第一反演???,用于对分角度叠加的数据体进行射线弹性阻抗反演,得到
    多个角度的射线弹性阻抗反演体;

    岩相划分???,用于根据岩性和孔隙度值对多个角度的射线弹性阻抗反演
    体进行岩相划分,得到测井岩相;

    阻抗体计算???,用于根据提取的井旁道数据及划分出的测井岩相进行计
    算得到AVO阻抗体;

    岩相体计算???,用于对AVO阻抗体的阻抗值域的划分,计算形成地震岩
    相体;

    第二反演???,用于将多个射线弹性阻抗反演体进行射线弹性阻抗反演,
    得到储层弹性参数体;

    交汇分析???,用于在地震岩相体的控制下,对储层弹性参数体进行散点
    交汇分析,获得对含气饱和度敏感的弹性参数体;

    压实分析???,用于对敏感弹性参数体提取井旁道数据,进行压实分析形
    成干岩石骨架体;

    饱和度计算???,用于根据利用干岩石骨架体计算得到的流体体积??楹?br />流体密度,定量计算含气饱和度。

    根据本公开的第二方面,在第二方面第一种可能的实现方式中,所述岩相
    划分??榘ǎ?br />

    第一提取子???,用于从多个角度的射线弹性阻抗反演体中提取小角度射
    线弹性阻抗和大角度射线弹性阻抗;

    第二提取子???,用于提取小角度射线弹性阻抗中的孔隙度信息及大角度
    射线弹性阻抗中的岩性信息;

    划分子???,用于根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度
    信息及大角度射线弹性阻抗中的岩性信息进行岩相划分,得到测井岩相。

    根据本公开的第二方面第一种可能的实现方式,在第二方面第二种可能的
    实现方式中,所述阻抗体计算??榘ǎ?br />

    第一分析子???,用于根据抽取的井旁道数据及划分出的测井岩相,进行
    散点交汇分析;

    坐标旋转子???,用于利用属性投影的方法对小角度射线弹性阻抗体和大
    角度射线弹性阻抗体进行坐标旋转,得到坐标转化公式;

    计算子???,用于根据散点交汇分析的结果及坐标转化公式对小角度射线
    弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。

    根据本公开的第二方面,在第二方面第三种可能的实现方式中,所述压实
    分析??榘ǎ?br />

    数据提取子???,用于对敏感弹性参数体提取井旁道数据;

    压实分析子???,用于对井旁道数据进行压实分析;

    获取子???,用于根据压实分析的结果获取干岩石骨架的属性变化趋势;

    插值子???,用于根据干岩石骨架的属性变化趋势进行多属性井插值形成
    干岩石骨架体。

    根据本公开的第二方面,在第二方面第四种可能的实现方式中,所述对含
    气饱和度敏感的弹性参数体包括密度和体积模量。

    本公开的各方面通过利用射线弹性阻抗交汇形成新的AVO阻抗反演地震岩
    相体,预测储层有利分布,在岩相体控制下排除非储层影响,进行定量含气饱
    和度转化,转化过程中,利用Gassmann方程的简化式,将孔隙岩石理解为干燥
    骨架与孔隙流体两部分,在去除干燥部分,得到的流体密度与含气饱和度之间
    存在线性关系,实现了排除现有叠前反演含气性预测方法的多解性,准确的预
    测有效储层中的含气饱和度效果。

    附图说明

    通过结合附图对本公开示例性实施方式进行更详细的描述,本公开的上述
    以及其它目的、特征和优势将变得更加明显,其中,在本公开示例性实施方式
    中,相同的参考标号通常代表相同部件。

    图1示出了本公开一示例性实施例提供的一种含气饱和度确定方法的流程
    示意图。

    图2示出了本公开一示例性实施例示出的步骤102的一种流程示意图。

    图3示出了本公开一示例性实施例示出的步骤103的一种流程示意图。

    图4示出了本公开一示例性实施例示出的步骤107的一种流程示意图。

    图5示出了本公开一示例性实施例提供的一种含气饱和度确定装置的结构
    示意图。

    图6示出了本公开一示例性实施例示出的岩相划分???02的一种结构示
    意图。

    图7示出了本公开一示例性实施例示出的阻抗体计算???03的一种结构
    示意图。

    图8示出了本公开一示例性实施例示出的压实分析???07的一种结构示
    意图。

    图9示出了本公开涉及的岩石物理模型示意图。

    具体实施方式

    下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了
    本公开的优选实施方式,然而应该理解,可以以各种形式实现本公开而不应被
    这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本公开更加透
    彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。

    图1示出了根据本发明的一个实施例的含气饱和度确定方法的流程图,如
    图1所示,该方法包括以下步骤。

    步骤101,对分角度叠加的数据体进行射线弹性阻抗反演得到多个角度的射
    线弹性阻抗反演体。

    所述步骤101在本公开实施例中可具体实现为,由于多个角度的射线弹性
    阻抗反演体中包括小角度射线弹性阻抗和大角度射线弹性阻抗,而小角度射线
    弹性阻抗含有孔隙度信息,大角度射线弹性阻抗含有岩性信息,所以可以通过
    因此可以通过大角度射线弹性阻抗(REI25)和小角度射线弹性阻抗(REI5)区分岩
    相。角度的数值范围一般看实际资料的角度范围而定,例如一般叠前道集入射
    角范围为0~30,平均分为三份,小角度为5,中角度为15,大角度为25。

    步骤102,根据岩性和孔隙度值对多个角度的射线弹性阻抗反演体进行岩相
    划分,得到测井岩相。

    步骤103,根据提取的井旁道数据及测井岩相进行计算得到AVO阻抗体。

    所述步骤103在本公开实施例中可具体实现为,交汇图中散点形状,可以
    通过测井岩相解释获得,如可令三角形表示泥岩,正方形为干砂岩储层,圆形
    为砂岩储层等。

    岩相体的数值范围应该是岩相分类值域(1,2,3…n)。在本实施例中,两个角
    度的射线弹性阻抗通过坐标转换,形成新的阻抗体(AVOImp),通过AVOImp与
    密度交汇图显示可以看出,AVO阻抗在值域上区分储层、干层及泥岩类型,通
    过阈值划分的方式,可以将AVO阻抗转化成岩性体(定义1—泥岩,2—干层,3
    —储层)。

    步骤104,对AVO阻抗体的阻抗值域的划分,计算形成地震岩相体。

    所述步骤104在本公开实施例中可具体实现为,岩相体LITH与AVO阻抗
    体应满足如下公式:


    value1和value2为区分泥岩的AVOImp的最小值、最大值;value1和value2
    为区分干层的AVOImp的最小值、最大值;value1和value2为区分储层的AVOImp
    的最小值、最大值。

    其中,AVO阻抗与输入的大小角度射线弹性阻抗体满足如下公式:

    AVOImp=A*X+B*Y+C

    式中,X是小角度射线弹性阻抗体,Y是大角度射线弹性阻抗体,A、B和
    C为拟合参数。

    步骤105,将多个射线弹性阻抗反演体进行射线弹性阻抗反演得到储层弹性
    参数体。

    所述步骤105在本公开实施例中可具体实现为,在岩相体的控制下,储层
    弹性参数可以用Gassmann方程进行描述,如图9所示,当储层为饱和岩石时:


    μs=μd



    ρs=(1-φ)ρm+φρf

    式中:Ks是饱和岩石的体积模量;Kd是干燥岩石骨架的体积模量;Km是岩
    石固体基质的体积模量;Kf是孔隙流体的体积模量;φ是孔隙度;μs是饱和岩
    石的剪切模量;μd是干燥岩石骨架的剪切模量;ρs是饱和岩石的密度。

    这里Gassmann方程孔隙流体只涉及到单相流体。

    步骤106,在地震岩相体的控制下,对储层弹性参数体进行散点交汇分析,
    获得对含气饱和度敏感的弹性参数体。

    所述步骤106在本公开实施例中可具体实现为,所述对含气饱和度敏感的
    弹性参数体包括密度和体积模量,如果假设流体包含气和水两相,流体密度可
    写成两相流体的加权和,流体的体积模量可用Wood方程来表示:


    或采用Brie经验流体混合方程,e为经验值,一般取2.5,值越小逼近Wood
    方程,越大逼近斑块饱和关系流体混合方程:

    Kf=(Kw-Kg)Sge+Kg

    流体密度ρf可写为:

    ρf=(1-Sg)ρw+Sgρg

    通过以上三个公式可知,当孔隙中含两相流体时,其中一相流体的饱和度
    可由流体体积模量用非线性公式来表达,流体密度用线性公式来表达。

    步骤107,对敏感弹性参数体提取井旁道数据,进行压实分析形成干岩石骨
    架体。

    所述步骤107在本公开实施例中可具体实现为,对敏感弹性参数体提取井
    旁道数据,由于Kd、ρd与孔隙度及成岩作用有关,随深度、压力趋势变化,可以
    通过测井曲线压实分析,形成干岩石骨架体。

    步骤108,根据利用干岩石骨架体计算得到的流体体积??楹土魈迕芏?,定
    量计算含气饱和度。

    通过对Gassmann方程重新整理,可以将流体密度和流体体积模量写成与饱
    和岩石、骨架、基质和孔隙度有关的量:




    从上式可以看出,要求得流体密度和流体体积模量,需要知道的物理量有:
    干岩石的体积模量及密度Kd、ρd,固体基质的体积模量及密度Km、ρm,以及孔隙
    度φ;由于孔隙度φ随深度变化,在特定物相分类中,可定义为常数;Km、ρm排
    除与泥质含量的关系,与孔隙度无关,不随深度变化,岩样数据测试及查表获
    得。

    本实施例通过利用射线弹性阻抗交汇形成新的AVO阻抗反演地震岩相体,
    预测储层有利分布,在岩相体控制下排除非储层影响,进行定量含气饱和度转
    化,转化过程中,利用Gassmann方程的简化式,将孔隙岩石理解为干燥骨架与
    孔隙流体两部分,在去除干燥部分,得到的流体密度与含气饱和度之间存在线
    性关系,实现了排除现有叠前反演含气性预测方法的多解性,准确的预测有效
    储层中的含气饱和度效果。

    由于小角度射线弹性阻抗含有孔隙度信息,大角度射线弹性阻抗含有岩性
    信息,因此可以通过大、小角度射线弹性阻抗区分岩相,如图2所示,所述步
    骤102包括:

    步骤1021,从多个角度的射线弹性阻抗反演体中提取小角度射线弹性阻抗
    和大角度射线弹性阻抗;

    步骤1022,提取小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻
    抗中的岩性信息;

    步骤1023,根据岩性和孔隙度值、小角度射线弹性阻抗中的孔隙度信息及
    大角度射线弹性阻抗中的岩性信息进行岩相划分,得到测井岩相。

    本实施例通过小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗
    中的岩性信息进行岩相划分,可以准确的得到测井岩相。

    为了能够在AVO阻抗上可以排除泥岩和干层的干扰,直接找到孔隙度与
    AVO阻抗的关系,如图3所示,所述步骤103包括:

    步骤1031,根据抽取的井旁道数据及划分出的测井岩相,进行散点交汇分
    析;

    步骤1032,利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹
    性阻抗体进行坐标旋转,得到坐标转化公式;

    步骤1033,根据散点交汇分析的结果及坐标转化公式对小角度射线弹性阻
    抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。

    本实施例通过根据抽取的井旁道数据及划分出的测井岩相,进行散点交汇
    分析,利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体
    进行坐标旋转,得到坐标转化公式,根据散点交汇分析的结果及坐标转化公式
    对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体,可以
    实现利用计算得到的AVO阻抗排除泥岩和干层的干扰,直接找到孔隙度与AVO
    阻抗的关系。

    为了能够去除压实作用的岩石骨架,获得反映流体的含气饱和度体,达到
    定量预测储层含气性目的,如图4所示,所述步骤107包括:

    步骤1071,对敏感弹性参数体提取井旁道数据;

    步骤1072,对井旁道数据进行压实分析;

    步骤1073,根据压实分析的结果获取干岩石骨架的属性变化趋势;

    步骤1074,根据干岩石骨架的属性变化趋势进行多属性井插值形成干岩石
    骨架体。

    本实施例通过对敏感弹性参数体提取井旁道数据,对井旁道数据进行压实
    分析,根据压实分析的结果获取干岩石骨架的属性变化趋势,根据干岩石骨架
    的属性变化趋势进行多属性井插值形成干岩石骨架体,实现了去除压实作用的
    岩石骨架,获得反映流体的含气饱和度体。

    另外,作为对上述各实施例的实现,本公开实施例还提供了一种含气饱和
    度确定装置,如图5所示,包括:第一反演???01、岩相划分???02、阻抗
    体计算???03、岩相体计算???04、第二反演???05、交汇分析???06、
    压实分析???07和饱和度计算???08。

    第一反演???01,用于对分角度叠加的数据体进行射线弹性阻抗反演,得
    到多个角度的射线弹性阻抗反演体。

    岩相划分???02,用于根据岩性和孔隙度值对多个角度的射线弹性阻抗反
    演体进行岩相划分,得到测井岩相。

    阻抗体计算???03,用于根据提取的井旁道数据及划分出的测井岩相进行
    计算得到AVO阻抗体。

    岩相体计算???04,用于对AVO阻抗体的阻抗值域的划分,计算形成地
    震岩相体。

    第二反演???05,用于将多个射线弹性阻抗反演体进行射线弹性阻抗反
    演,得到储层弹性参数体。

    交汇分析???06,用于在地震岩相体的控制下,对储层弹性参数体进行散
    点交汇分析,获得对含气饱和度敏感的弹性参数体,所述对含气饱和度敏感的
    弹性参数体包括密度和体积模量。

    压实分析???07,用于对敏感弹性参数体提取井旁道数据,进行压实分析
    形成干岩石骨架体;

    饱和度计算???08,用于根据利用干岩石骨架体计算得到的流体体积???br />和流体密度,定量计算含气饱和度。

    本实施例通过利用射线弹性阻抗交汇形成新的AVO阻抗反演地震岩相体,
    预测储层有利分布,在岩相体控制下排除非储层影响,进行定量含气饱和度转
    化,转化过程中,利用Gassmann方程的简化式,将孔隙岩石理解为干燥骨架与
    孔隙流体两部分,在去除干燥部分,得到的流体密度与含气饱和度之间存在线
    性关系,实现了排除现有叠前反演含气性预测方法的多解性,准确的预测有效
    储层中的含气饱和度效果。

    由于小角度射线弹性阻抗含有孔隙度信息,大角度射线弹性阻抗含有岩性
    信息,因此可以通过大、小角度射线弹性阻抗区分岩相,如图6所示,所述岩
    相划分???02包括:第一提取子???021、第二提取子???022和划分子模
    块2023。

    第一提取子???021,用于从多个角度的射线弹性阻抗反演体中提取小角
    度射线弹性阻抗和大角度射线弹性阻抗。

    第二提取子???022,用于提取小角度射线弹性阻抗中的孔隙度信息及大
    角度射线弹性阻抗中的岩性信息。

    划分子???023,用于根据岩性和孔隙度值、小角度射线弹性阻抗中的孔
    隙度信息及大角度射线弹性阻抗中的岩性信息进行岩相划分,得到测井岩相。

    本实施例通过小角度射线弹性阻抗中的孔隙度信息及大角度射线弹性阻抗
    中的岩性信息进行岩相划分,可以准确的得到测井岩相。

    为了能够在AVO阻抗上可以排除泥岩和干层的干扰,直接找到孔隙度与
    AVO阻抗的关系,如图7所示,所述阻抗体计算???03包括:第一分析子模
    块2031、坐标旋转子???032和计算子???033。

    第一分析子???031,用于根据抽取的井旁道数据及划分出的测井岩相,
    进行散点交汇分析。

    坐标旋转子???032,用于利用属性投影的方法对小角度射线弹性阻抗体
    和大角度射线弹性阻抗体进行坐标旋转,得到坐标转化公式。

    计算子???033,用于根据散点交汇分析的结果及坐标转化公式对小角度
    射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体。

    本实施例通过根据抽取的井旁道数据及划分出的测井岩相,进行散点交汇
    分析,利用属性投影的方法对小角度射线弹性阻抗体和大角度射线弹性阻抗体
    进行坐标旋转,得到坐标转化公式,根据散点交汇分析的结果及坐标转化公式
    对小角度射线弹性阻抗和大角度射线弹性阻抗进行计算得到AVO阻抗体,可以
    实现利用计算得到的AVO阻抗排除泥岩和干层的干扰,直接找到孔隙度与AVO
    阻抗的关系。

    为了能够去除压实作用的岩石骨架,获得反映流体的含气饱和度体,达到
    定量预测储层含气性目的,如图8所示,所述压实分析???07包括:数据提
    取子???071、压实分析子???072、获取子???073和插值子???074。

    数据提取子???071,用于对敏感弹性参数体提取井旁道数据。

    压实分析子???072,用于对井旁道数据进行压实分析。

    获取子???073,用于根据压实分析的结果获取干岩石骨架的属性变化趋
    势。

    插值子???074,用于根据干岩石骨架的属性变化趋势进行多属性井插值
    形成干岩石骨架体。

    本实施例通过对敏感弹性参数体提取井旁道数据,对井旁道数据进行压实
    分析,根据压实分析的结果获取干岩石骨架的属性变化趋势,根据干岩石骨架
    的属性变化趋势进行多属性井插值形成干岩石骨架体,实现了去除压实作用的
    岩石骨架,获得反映流体的含气饱和度体。

    应用示例

    为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。
    本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非
    意在以任何方式限制本发明。

    川西某地区蓬莱镇组发育多套气层,展示了很好的勘探潜力。但在开发过
    程中,发现这种低孔低渗储层气藏气水关系复杂,不同含气丰度气层的地震响
    应特征差异不明显,导致常规油气检测方法难以解决本区气水预测问题,针对
    本区地质特点及地震资料现状,利用岩相控制下含气饱和度预测解决了当前的
    开发难点。

    研究区定义地震岩相综合考虑了沉积作用和成岩作用、孔隙度,并结合了
    测井储层评价的标准,对研究区岩相划分出三类岩相(表1),通过反演获得地
    震岩相体。通过岩相体划分,将I、II类物相孔隙度定义为

    表1‐研究区岩相划分标准

    岩相
    储层级别
    划分标准
    I类
    1类储层、2类储层
    7<Por<15
    II类
    干层
    Por<7
    III类
    泥岩
    Por<7

    在岩相计算基础上,采用四维交汇岩石物理分析,即密度(横坐标)、纵横
    波速度比(纵坐标)、孔隙度(颜色)以及地震岩相(符号)四种参数进行交汇
    分析,通过交会分析可以知道,在有利储层区,在一定孔隙度控制下密度和含
    气饱和度呈线性关系。

    利用密度参数作为我们本次预测含气丰度的主要参数。在此基础上将流体
    密度通过Gassman方程变形式转化为含气饱和度。图中红、黄色代表气层,绿
    色代表气水区,蓝色为水层,白色为储层不发育段,预测结果与各测试井结果
    基本吻合,有效验证了方法的适用性。

    本公开可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括
    计算机可读存储介质,其上载有用于使处理器实现本公开的各个方面的计算机
    可读程序指令。

    计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令
    的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、
    磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合
    适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:便携
    式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可
    编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压
    缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编
    码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合
    适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身,诸如
    无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波
    (例如,通过光纤电缆的光脉冲)、或者通过电线传输的电信号。

    这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个
    计算/处理设备,或者通过网络、例如因特网、局域网、广域网和/或无线网下载
    到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传
    输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设
    备中的网络适配卡或者网络接口从网络接收计算机可读程序指令,并转发该计
    算机可读程序指令,以供存储在各个计算/处理设备中的计算机可读存储介质中。

    用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构
    (ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、
    或者以一种或多种编程语言的任意组合编写的源代码或目标代码,所述编程语
    言包括面向对象的编程语言—诸如Smalltalk、C++等,以及常规的过程式编程语
    言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计
    算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分
    在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上
    执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络—包
    括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计
    算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中,通过
    利用计算机可读程序指令的状态信息来个性化定制电子电路,例如可编程逻辑
    电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA),该电子电路可
    以执行计算机可读程序指令,从而实现本公开的各个方面。

    这里参照根据本公开实施例的方法、装置(系统)和计算机程序产品的流
    程图和/或框图描述了本公开的各个方面。应当理解,流程图和/或框图的每个方
    框以及流程图和/或框图中各方框的组合,都可以由计算机可读程序指令实现。

    这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编
    程数据处理装置的处理器,从而生产出一种机器,使得这些指令在通过计算机
    或其它可编程数据处理装置的处理器执行时,产生了实现流程图和/或框图中的
    一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令
    存储在计算机可读存储介质中,这些指令使得计算机、可编程数据处理装置和/
    或其他设备以特定方式工作,从而,存储有指令的计算机可读介质则包括一个
    制造品,其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的
    各个方面的指令。

    也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、
    或其它设备上,使得在计算机、其它可编程数据处理装置或其它设备上执行一
    系列操作步骤,以产生计算机实现的过程,从而使得在计算机、其它可编程数
    据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方
    框中规定的功能/动作。

    附图中的流程图和框图显示了根据本公开的多个实施例的系统、方法和计
    算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图
    中的每个方框可以代表一个???、程序段或指令的一部分,所述???、程序段
    或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有
    些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺
    序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以
    按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图
    中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能
    或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的
    组合来实现。

    以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,
    并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的
    情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
    本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场
    中的技术的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露
    的各实施例。

       内容来自专利网重庆时时彩单双窍门 www.4mum.com.cn转载请标明出处

    关于本文
    本文标题:一种含气饱和度确定方法和装置.pdf
    链接地址://www.4mum.com.cn/p-6014397.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 浪潮孙丕恕从信息时代到智能时代 人工智能价值将爆发式释放 2019-12-21
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 球探即时指数 彩票站分布图 安徽快三 顶级信誉棋牌平台 百人牛牛攻略 江苏微乐棋牌麻将 三肖两码中特网 江苏十一选五任三遗漏一定牛 彩票开奖设计 华夏彩票首页 广东26选5最后一期开奖结果 秒速飞艇开奖走势 球探体育比分4.8 甘肃快三遗漏数据统计 江苏快3下载安装 黄金城棋牌游戏官方网站