• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 12
    • 下载费用:30 金币  

    重庆时时彩龙虎算法: 一种裂缝性砂岩储层含油饱和度的计算方法.pdf

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

    CN201611090890.5

    申请日:

    2016.12.01

    公开号:

    CN106772674A

    公开日:

    2017.05.31

    当前法律状态:

    授权

    有效性:

    有权

    法律详情: 授权|||实质审查的生效IPC(主分类):G01V 11/00申请日:20161201|||公开
    IPC分类号: G01V11/00 主分类号: G01V11/00
    申请人: 中国石油天然气股份有限公司
    发明人: 朱光有; 赵斌
    地址: 100007 北京市东城区东直门北大街9号
    优先权:
    专利代理机构: 北京三友知识产权代理有限公司 11127 代理人: 张德斌;姚亮
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201611090890.5

    授权公告号:

    ||||||

    法律状态公告日:

    2018.11.16|||2017.06.23|||2017.05.31

    法律状态类型:

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

    摘要

    本发明提供了种裂缝性砂岩储层含油饱和度的计算方法。方法包括:(1)测定砂岩的弹性模量、泊松比、裂缝法向刚度、裂缝间距以及储层初始有效孔隙度;(2)确定储层主地应力方向,确定裂缝倾向与最大主地应力方向夹角与裂缝面倾角,计算储层垂向应力、最大水平主地应力与最小水平主地应力;(3)测定储层孔隙压力;(4)计算储层三个主地应力方向上的有效孔隙;(5)测定地层水的电阻率,储层在三个主地应力方向上的电阻率,计算含100%水时储层的电阻率:(6)计算储层的含油饱和度,测定不同含水饱和度的砂岩样品的电阻率,使用最小二乘法原理拟合出曲线,求解不同有效孔隙度砂岩的含油饱和度。

    权利要求书

    1.一种裂缝性砂岩储层含油饱和度的计算方法,其中,所述方法包括如下步骤:
    (1)测定裂缝性砂岩储层的弹性模量、泊松比、裂缝的法向刚度、裂缝间距以及砂岩的
    初始有效孔隙度;并测定不包含裂缝的砂岩在三个主地应力方向上的弹性模量;
    (2)确定裂缝性砂岩储层主地应力方向,确定裂缝倾向与最大主地应力方向的夹角α、
    裂缝面的倾角β,计算储层垂向应力、最大水平主地应力与最小水平主地应力;
    (3)测定裂缝性砂岩储层的孔隙压力;
    (4)按照如下公式计算油气开采过程中裂缝性砂岩储层三个主地应力方向上的有效有
    效孔隙度
    <mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>&phi;</mi> <mi>p</mi> <mn>1</mn> </msubsup> <mo>=</mo> <msub> <mi>&phi;</mi> <mn>0</mn> </msub> <mo>+</mo> <mfrac> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>1</mn> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>1</mn> </msubsup> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <msub> <mi>&mu;</mi> <mn>12</mn> </msub> <mo>+</mo> <msub> <mi>&mu;</mi> <mn>13</mn> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&sigma;</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>sin</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>2</mn> </msubsup> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>sin</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>&mu;</mi> <mn>23</mn> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&sigma;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>3</mn> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&beta;</mi> </mrow> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>3</mn> </msubsup> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&beta;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>&mu;</mi> <mn>32</mn> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&sigma;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
    <mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>&phi;</mi> <mi>p</mi> <mn>2</mn> </msubsup> <mo>=</mo> <msub> <mi>&phi;</mi> <mn>0</mn> </msub> <mo>-</mo> <mfrac> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>1</mn> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>1</mn> </msubsup> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>&mu;</mi> <mn>13</mn> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&sigma;</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>sin</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>2</mn> </msubsup> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>sin</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <msub> <mi>&mu;</mi> <mn>21</mn> </msub> <mo>+</mo> <msub> <mi>&mu;</mi> <mn>23</mn> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&sigma;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>3</mn> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&beta;</mi> </mrow> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>3</mn> </msubsup> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&beta;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>&mu;</mi> <mn>31</mn> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&sigma;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
    <mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>&phi;</mi> <mi>p</mi> <mn>3</mn> </msubsup> <mo>=</mo> <msub> <mi>&phi;</mi> <mn>0</mn> </msub> <mo>-</mo> <mfrac> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>1</mn> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>1</mn> </msubsup> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>&mu;</mi> <mn>13</mn> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&sigma;</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>sin</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>2</mn> </msubsup> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>sin</mi> <mi>&alpha;</mi> <mi>sin</mi> <mi>&beta;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>&mu;</mi> <mn>21</mn> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&sigma;</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mfrac> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>3</mn> </msubsup> <mo>+</mo> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&beta;</mi> </mrow> <mrow> <msubsup> <mi>E</mi> <mi>s</mi> <mn>3</mn> </msubsup> <msub> <mi>K</mi> <mi>f</mi> </msub> <msub> <mi>s</mi> <mi>f</mi> </msub> <mi>cos</mi> <mi>&beta;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <mrow> <msub> <mi>&mu;</mi> <mn>31</mn> </msub> <mo>+</mo> <msub> <mi>&mu;</mi> <mn>32</mn> </msub> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&sigma;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>p</mi> <mi>p</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
    式中,为最大主地应力方向上的有效孔隙度,无量纲;为中间主地应力方向上的有
    效孔隙度,无量纲;为最小主地应力方向上的有效孔隙度,无量纲;φ0为裂缝性砂岩储层
    的初始有效孔隙度,无量纲;为完整砂岩在最大主地应力方向上的弹性模量,MPa;为完
    整砂岩在中间主地应力方向上的弹性模量,MPa;为完整砂岩在最小主地应力方向上的弹
    性模量,MPa;α为裂缝倾向与最大主地应力方向的夹角,°;β为裂缝面的倾角,°;Kf为裂缝的
    法向刚度,MPa·m-1;sf为裂缝间距,m;μij为泊松比,表征由i方向的单位拉应变在j方向上产
    生的压缩应变,i,j=1,2,3,无量纲;pp为裂缝性砂岩储层孔隙压力,无量纲;σ1为最大主地
    应力,MPa;σ2为中间主地应力,MPa;σ3为最小主地应力,MPa;储层垂向应力σv、最大水平主地
    应力σH、最小水平主地应力σh依大小排序,最大值赋给σ1,中间值赋给σ2,最小值赋给σ3;
    (5)测定地层水的电阻率Rw,裂缝性砂岩储层在三个主地应力方向上的电阻率,依据如
    下公式计算100%含水时裂缝性砂岩储层的电阻率:
    <mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>R</mi> <mn>0</mn> <mi>i</mi> </msubsup> <mo>=</mo> <msub> <mi>a</mi> <mi>i</mi> </msub> <msub> <mi>R</mi> <mi>w</mi> </msub> <mo>/</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>&phi;</mi> <mi>p</mi> <mi>i</mi> </msubsup> <mo>)</mo> </mrow> <msub> <mi>m</mi> <mi>i</mi> </msub> </msup> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> </mrow> </mtd> </mtr> </mtable> </mfenced>
    式中,为100%含水时裂缝性砂岩在三个主地应力方的电阻率,Ω·m;ai为三个主地
    应力方向的岩性系数,无量纲;mi为三个主地应力方向的孔隙指数,i=1,2,3,无量纲;为
    裂缝性砂岩在三个主地应力方向上的有效孔隙度,无量纲;Rw为地层水的电阻率,Ω·m;
    (6)计算裂缝性砂岩储层的含油饱和度,测定不同含水饱和度的砂岩样品的电阻率
    使用最小二乘法原理拟合出曲线,依据曲线方程求解不同有效孔隙度砂岩
    的含油饱和度,其中i=1,2,3。
    2.根据权利要求1所述的计算方法,其中,步骤(1)是按照《岩石力学试验建议方法》测
    量砂岩储层的弹性模量、泊松比、裂缝的法向刚度、裂缝间距以及砂岩的初始有效孔隙度,
    和不包含裂缝的砂岩储层在三个主地应力方向上的弹性模量。
    3.根据权利要求1所述的计算方法,其中,步骤(1)是通过钻井取心,使用单轴压缩试验
    测定不含裂缝的完整砂岩块在三个主地应力方向的弹性模量并测定包含裂缝
    在内的砂岩储层在三个主地应力方向相关的泊松比μ;选择含裂缝的砂岩块,测定裂缝的法
    向刚度Kf和有效孔隙度φ0,利用裂缝成像技术和岩心观测方法测定裂缝性砂岩储层裂缝间
    距sf。
    4.根据权利要求1所述的计算方法,其中,步骤(2)是利用井壁裂缝成像分析技术确定
    裂缝性砂岩储层的主地应力方向;使用密度积分方法计算垂向应力的大小,使用组合弹簧
    模型计算最大水平主地应力和最小水平主地应力的大小。
    5.根据权利要求1所述的计算方法,其中,步骤(2)中计算垂向应力是使用密度积分法,
    公式如下
    <mrow> <msub> <mi>&sigma;</mi> <mi>v</mi> </msub> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>&rho;</mi> <mi>i</mi> </msub> <msub> <mi>gh</mi> <mi>i</mi> </msub> </mrow>
    式中,σv为垂向地应力,MPa;ρi为地层岩石密度,g/cm3;g为重力加速度,取9.8m/s2;hi为
    地层岩石厚度,m。
    6.根据权利要求1所述的计算方法,其中,步骤(2)中使用组合弹簧模型计算最大水平
    主地应力与最小水平主地应力,公式如下:
    <mrow> <msub> <mi>&sigma;</mi> <mi>H</mi> </msub> <mo>=</mo> <mfrac> <mi>E</mi> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>&mu;</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msub> <mi>&epsiv;</mi> <mi>H</mi> </msub> <mo>+</mo> <mfrac> <mrow> <mi>&mu;</mi> <mi>E</mi> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>&mu;</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msub> <mi>&epsiv;</mi> <mi>h</mi> </msub> <mo>+</mo> <mfrac> <mi>&mu;</mi> <mrow> <mn>1</mn> <mo>-</mo> <mi>&mu;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>&sigma;</mi> <mi>v</mi> </msub> <mo>-</mo> <msub> <mi>&alpha;p</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&alpha;p</mi> <mi>p</mi> </msub> </mrow>
    <mrow> <msub> <mi>&sigma;</mi> <mi>h</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&mu;</mi> <mi>E</mi> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>&mu;</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msub> <mi>&epsiv;</mi> <mi>H</mi> </msub> <mo>+</mo> <mfrac> <mi>E</mi> <mrow> <mn>1</mn> <mo>-</mo> <msup> <mi>&mu;</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msub> <mi>&epsiv;</mi> <mi>h</mi> </msub> <mo>+</mo> <mfrac> <mi>&mu;</mi> <mrow> <mn>1</mn> <mo>-</mo> <mi>&mu;</mi> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>&sigma;</mi> <mi>v</mi> </msub> <mo>-</mo> <msub> <mi>&alpha;p</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&alpha;p</mi> <mi>p</mi> </msub> </mrow>
    式中,σH为水平最大主地应力,MPa;σh为水平最小主地应力,MPa;pp为储层孔隙压力,E
    为裂缝性砂岩储层的弹性模量,MPa;μ为岩石泊松比;σv为垂向应力,MPa;α为Biot系数;εH最
    大水平主地应力方向的应变,无量纲;εh为最小水平主地应力方向的应变,无量纲。
    7.根据权利要求1所述的计算方法,其中,步骤(3)中采用重复地层测试(RFT)测定裂缝
    性砂岩储层的孔隙压力。
    8.根据权利要求1所述的计算方法,其中,步骤(5)的三个主地应力方向上的岩性系数
    ai,和三个主地应力方向上的孔隙指数mi,其中i=1,2,3,通过如下方法测定:在三个主地应
    力方向钻取岩心,测定地层水的电阻率Rw,测定不同有效孔隙度时的电阻率,作
    曲线,使用最小二乘法拟合曲线,由图线求取裂缝性砂岩储层三个主地应力方
    向上的岩性系数ai,以及三个主地应力方向上的孔隙指数mi。
    9.根据权利要求1所述的计算方法,其中,步骤(5)的岩性系数的测定方法为:采集多组
    岩样,分别测定三个主地应力方向的R0,以及三个主地应力方向的φp,在双对数坐标系中拟
    合出直线纵坐标截距为logai,ai即为岩性系数。
    10.根据权利要求1所述的计算方法,其中,步骤(6)按照如下公式计算裂缝性砂岩储层
    的含油饱和度:
    <mrow> <msubsup> <mi>S</mi> <mrow> <mi>o</mi> <mi>g</mi> </mrow> <mi>i</mi> </msubsup> <mo>=</mo> <mn>1</mn> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>b</mi> <mi>i</mi> </msub> <msubsup> <mi>R</mi> <mn>0</mn> <mi>i</mi> </msubsup> </mrow> <msubsup> <mi>R</mi> <mi>t</mi> <mi>i</mi> </msubsup> </mfrac> <mo>)</mo> </mrow> <mfrac> <mn>1</mn> <msup> <mi>n</mi> <mi>i</mi> </msup> </mfrac> </msup> </mrow>
    式中,为含100%水时裂缝性砂岩在三个主地应力方向的电阻率,Ω·m;为含油气
    裂缝性砂岩在三个主地应力方向上的电阻率,Ω·m;bi为实验拟合参数,无量纲;ni为实验
    拟合参数,无量纲,其中i=1,2,3。

    说明书

    一种裂缝性砂岩储层含油饱和度的计算方法

    技术领域

    本发明涉及油气勘探开发技术领域,具体地说,是涉及一种裂缝性砂岩储层含油
    饱和度的计算方法。

    背景技术

    含油饱和度的准确测定是判断储层品质和调整开发方案的基础资料,裂缝性砂岩
    储层的岩石力学性质与电性具有显著的各向异性特征,使用电阻率测井评价裂缝性储层含
    油性的生产实践中发现,在构造挤压区和张拉区,电阻率值出现较为明显的偏差,造成储层
    含油性评价的失误。裂缝作为一种不连续面,其存在会导致储层在各个方向上的弹性模量、
    泊松比以及孔渗性能产生较大差异,最终造成含油性的各向异性。到目前,还没有一种能够
    反映这种含油性各向异性差异的定量计算方法。

    发明内容

    本发明的一个目的在于提供一种裂缝性砂岩储层含油饱和度的计算方法。本发明
    通过综合考虑垂向应力、构造侧向挤压作用,以及采油时间等因素对裂缝性砂岩储层有效
    孔隙度变化的影响,发明了计算裂缝性砂岩储层含油饱和度的计算方法,可用于生产过程
    中裂缝性砂岩储层含油饱和度的定量计算,预测裂缝性砂岩储层的含油性,解决了地应力
    环境对电阻率测井精度的影响问题,降低了裂缝性砂岩储层含油饱和度测定的成本,填补
    了地应力对测井电阻率数据影响研究的国际空白,为砂岩油气勘探与开发提供技术支持。

    为达上述目的,一方面,本发明提供了一种裂缝性砂岩储层含油饱和度的计算方
    法,其中,所述方法包括如下步骤:

    (1)测定裂缝性砂岩储层的弹性模量、泊松比、裂缝的法向刚度、裂缝间距以及砂
    岩储层的初始有效孔隙度;并测定不包含裂缝的砂岩在三个主地应力方向上的弹性模量;

    (2)确定裂缝性砂岩储层主地应力方向,确定裂缝倾向与最大主地应力方向的夹
    角α、裂缝面的倾角β,计算储层垂向应力、最大水平主地应力与最小水平主地应力;

    (3)测定裂缝性砂岩储层的孔隙压力;

    (4)按照如下公式计算油气开采过程中裂缝性砂岩储层三个主地应力方向上的有
    效有效孔隙度




    式中,为最大主地应力方向上的有效孔隙度,无量纲;为中间主地应力方向上
    的有效孔隙度,无量纲;为最小主地应力方向上的有效孔隙度,无量纲;φ0为裂缝性砂岩
    储层的初始有效孔隙度,无量纲;为完整砂岩在最大主地应力方向上的弹性模量,MPa;
    为完整砂岩在中间主地应力方向上的弹性模量,MPa;为完整砂岩在最小主地应力方向上
    的弹性模量,MPa;α为裂缝倾向与最大主地应力方向的夹角,°;β为裂缝面的倾角,°;Kf为裂
    缝的法向刚度,MPa·m-1;sf为裂缝间距,m;μij为泊松比,表征由i方向的单位拉应变在j方向
    上产生的压缩应变,i,j=1,2,3,无量纲;pp为裂缝性砂岩储层孔隙压力,无量纲;σ1为最大
    主地应力,MPa;σ2为中间主地应力,MPa;σ3为最小主地应力,MPa;储层垂向应力σv、最大水平
    主地应力σH、最小水平主地应力σh依大小排序,最大值赋给σ1,中间值赋给σ2,最小值赋给σ3;

    (5)测定地层水的电阻率Rw,裂缝性砂岩储层在三个主地应力方向上的电阻率,依
    据如下公式计算含100%水时裂缝性砂岩储层的电阻率:


    式中,为含100%水时裂缝性砂岩在三个主地应力方的电阻率,Ω·m;ai为三个
    主地应力方向的岩性系数,无量纲;mi为三个主地应力方向的孔隙指数,i=1,2,3,无量纲;
    为裂缝性砂岩在三个主地应力方向上的有效孔隙度,无量纲;Rw为地层水的电阻率,Ω·
    m;

    (6)计算裂缝性砂岩储层的含油饱和度,测定不同含水饱和度的砂岩样品的电
    阻率使用最小二乘法原理拟合出曲线,依据曲线方程求解不同有效孔隙度
    砂岩的含油饱和度,其中i=1,2,3。

    根据本发明一些具体实施方案,其中,步骤(1)是按照《岩石力学试验建议方法》测
    定砂岩的弹性模量、泊松比、裂缝的法向刚度、裂缝间距以及砂岩的初始有效孔隙度,和不
    包含裂缝的砂岩在三个主地应力方向上的弹性模量。

    根据本发明一些具体实施方案,其中,步骤(1)是通过钻井取心,使用单轴压缩试
    验测定不含裂缝的完整砂岩块在三个主地应力方向的弹性模量并测定包含
    裂缝的砂岩在三个主地应力方向相关的泊松比μ;选择含裂缝的砂岩块,测定裂缝的法向刚
    度Kf和有效孔隙度φ0,利用裂缝成像技术和岩心观测方法测定裂缝性砂岩储层裂缝间距
    sf。

    根据本发明一些具体实施方案,其中,步骤(2)是利用井壁裂缝成像分析技术确定
    裂缝性砂岩储层的主地应力方向;使用密度积分方法计算垂向应力的大小,使用组合弹簧
    模型计算最大水平主地应力和最小水平主地应力的大小。

    根据本发明一些具体实施方案,其中,步骤(2)中计算垂向应力是使用密度积分
    法,公式如下


    式中,σv为垂向地应力,MPa;ρi为地层岩石密度,g/cm3;g为重力加速度,取9.8m/
    s2;hi为地层岩石厚度,m。

    根据本发明一些具体实施方案,其中,步骤(2)中使用组合弹簧模型计算最大水平
    主地应力与最小水平主地应力,公式如下:



    式中,σH为水平最大主地应力,MPa;σh为水平最小主地应力,MPa;pp为储层孔隙压
    力,E为岩石弹性模量,MPa;μ为岩石泊松比;σv为垂向应力,MPa;α为Biot系数;εH最大水平主
    地应力方向的应变,无量纲;εh为最小水平主地应力方向的应变,无量纲。

    根据本发明一些具体实施方案,其中,步骤(3)中采用重复地层测试(RFT)测定裂
    缝性砂岩储层的孔隙压力。

    根据本发明一些具体实施方案,其中,步骤(5)的三个主地应力方向上的岩性系数
    ai,和三个主地应力方向上的孔隙指数mi,其中i=1,2,3,通过如下方法测定:在三个主地应
    力方向钻取岩心,测定地层水的电阻率Rw,测定不同有效有效孔隙度时的电阻率,作
    曲线,使用最小二乘法拟合曲线,由图线求取裂缝性砂岩储层三个主地应力方
    向上的岩性系数ai,以及三个主地应力方向上的孔隙指数mi。

    根据本发明一些具体实施方案,其中,步骤(5)所述的岩性系数的测定方法为:采
    集多组岩样,分别测定三个主地应力方向的R0,以及三个主地应力方向的φp,在双对数坐标
    系中拟合出直线纵坐标截距为logai,ai即为岩性系数。

    根据本发明一些具体实施方案,其中,步骤(6)按照如下公式计算裂缝性砂岩储层
    的含油饱和度:


    式中,为含100%水时裂缝性砂岩在三个主地应力方向的电阻率,Ω·m;为含
    油气裂缝性砂岩在三个主地应力方向上的电阻率,Ω·m;bi为实验拟合参数,无量纲;ni为
    实验拟合参数,无量纲,其中i=1,2,3。

    综上所述,本发明提供了一种裂缝性砂岩储层含油饱和度的计算方法。本发明的
    方法具有如下优点:

    本发明提供了一种考虑裂缝对储层含油性影响的储层含油饱和度计算方法,该方
    法在主地应力空间中评价裂缝性砂岩储层的各向异性力学性质与电性特征,将裂缝性砂岩
    储层三个主地应力方向的有效孔隙度作为关于主地应力、储层砂岩和裂缝力学参数,以及
    孔隙压力的函数,可用于计算不同地应力环境下储层的有效孔隙度,也可用于计算油气井
    不同生产时期储层的有效孔隙度。省去了用声波速度测井、密度测井或岩性-密度测井等方
    法测定储层有效孔隙度的工序,降低了含油饱和度测定的成本。而且,解决了由于构造挤压
    和张拉作用导致的储层含油性评价失误问题,提高了储层含油性测定的精度。为勘探开发
    储层品质的判别提供了方法技术。

    附图说明

    图1为实施例1的裂缝性砂岩储层取心示意图;

    说明:以较大尺寸从储层中钻取岩心(图1左侧),以较小尺寸在三个主地应力方向
    上钻取岩心(图1右侧),加工成标准试件,进行岩体力学参数测试。测定裂缝性砂岩储层中
    砂岩(不含裂缝)在三个主地应力方向上的弹性模量测定裂缝性砂岩储层中
    砂岩(含裂缝)的泊松比μij(i,j=1,2,3)。

    图2为实施例1的主地应力空间中裂缝面的倾向、倾角与裂缝面间距示意图;

    说明:1轴、2轴、3轴分别表示最大主地应力方向、中间主地应力方向、最小主地应
    力方向,1、2、3坐标系表示主地应力坐标系,其表征的空间为主地应力空间。使用井壁裂缝
    成像技术与岩心裂缝观测等方法,确定主地应力方向,裂缝面倾向与最大主地应力方向的
    夹角α,裂缝的倾角β,并测定裂缝性砂岩储层裂缝间距sf。裂缝间距为相邻裂缝面之间的距
    离,储层中构造裂缝的产状一般具有较强的规律性,可以测定一组或多组构造裂缝的间距。

    图3为实施例1的裂缝性砂岩储层裂缝法向刚度测定示意图;

    说明:测定裂缝性砂岩储层裂缝的法向刚度Kf,垂直裂缝面方向施加应力σn,监测
    裂缝面闭合位移δn,以定义式Kf=σn/δn计算裂缝的法向刚度。

    图4为实施例1的储层深度-有效孔隙度曲线;

    说明:计算含一组构造裂缝情况下裂缝性砂岩储层在三个主地应力方向的有效孔
    隙度,储层深度6600-6800m,由曲线可知,裂缝性砂岩储层有效孔隙度具有显著的各向异性
    特征。

    图5为实施例1的储层深度-含油饱和度曲线;

    说明:计算含一组构造裂缝情况下裂缝性砂岩储层在三个主地应力方向的含油饱
    和度,储层深度6600-6800m,由曲线可知,裂缝性砂岩储层含油性具有显著的各向异性特
    征。

    具体实施方式

    以下通过具体实施例详细说明本发明的实施过程和产生的有益效果,旨在帮助阅
    读者更好地理解本发明的实质和特点,不作为对本案可实施范围的限定。

    实施例1

    (1)对某开发井钻井至裂缝性砂岩储层6600-6800m深度时采集岩心,统计裂缝间
    距、倾向与倾角,使用声发射或古地磁方法确定最大主地应力方向,统计裂缝倾角与最大主
    地应力方向的夹角,裂缝面与主地应力方向的关系如图2所示。在岩心柱不同深度位置沿三
    个主地应力方向钻取直径25mm,长度50mm的岩样,如图1;

    裂缝产状与力学性能表

    裂缝倾角/°
    夹角/°
    裂缝法向刚度/MPa·m-1
    裂缝间距/m
    60
    60
    20000
    0.1

    裂缝性砂岩储层在三个主地应力方向上的岩性系数、孔隙指数、拟合系数值

    方向
    a
    b
    m
    n
    最大水平主地应力方向
    1.00
    1.50
    1.80
    4.00
    最小水平主地应力方向
    1.00
    0.45
    2.00
    3.20
    垂向地应力方向
    1.00
    0.65
    0.65
    2.80

    (2)参考国际岩石力学学会《岩石力学试验建议方法》测定岩样的弹性模量与泊松
    比,并使用含裂缝岩样测定裂缝的法向刚度,如图3;

    (3)使用压汞法测定裂缝性砂岩的有效孔隙度;

    (4)采用多层位测试器(FMT)测试,测定储层不同深度的孔隙压力;

    (5)采用密度积分法计算储层不同深度的垂向地应力,使用组合弹簧模型计算最
    大水平主地应力与最小水平主地应力;

    (6)使用权利要求1步骤四中的公式,计算三个主地应力方向上的孔隙度,计算结
    果如图4所示;

    (7)使用权利要求13中的公式,计算三个主地应力方向上的含油饱和度,计算结果
    如图5所示。

    关于本文
    本文标题:一种裂缝性砂岩储层含油饱和度的计算方法.pdf
    链接地址://www.4mum.com.cn/p-6003102.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

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


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 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
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • 彩票大乐透周一走势图 体彩双色球开奖直播现场直播 任选9场胜负查询 快乐8登录首页 吉林时时彩走势图开奖 时时彩平台哪个好 山东手机彩票投注 幸运赛车视频直播观看 秒速时时彩正规吗 体彩排列三走势图带连线 内蒙古时时彩11选5开奖结果 甘肃十一选五任五遗漏 体彩顶呱刮中奖图片 双色球奇偶分布图100期 山东省彩票开奖时间表 三分彩全天计划手机版