克拉玛依陆相油砂注水过程中的剪胀诱导渗透率模型

高彦芳 王晓阳 任战利 陈勉 姜海龙

高彦芳,王晓阳,任战利,陈勉,姜海龙. 克拉玛依陆相油砂注水过程中的剪胀诱导渗透率模型[J]. 石油钻采工艺,2021,43(2):244-249 doi:  10.13639/j.odpt.2021.02.017
引用本文: 高彦芳,王晓阳,任战利,陈勉,姜海龙. 克拉玛依陆相油砂注水过程中的剪胀诱导渗透率模型[J]. 石油钻采工艺,2021,43(2):244-249 doi:  10.13639/j.odpt.2021.02.017
GAO Yanfang, WANG Xiaoyang, REN Zhanli, CHEN Mian, JIANG Hailong. A model of shear dilation induced permeability evolution in the process of water injection in Karamay continental oil sand[J]. Oil Drilling & Production Technology, 2021, 43(2): 244-249 doi:  10.13639/j.odpt.2021.02.017
Citation: GAO Yanfang, WANG Xiaoyang, REN Zhanli, CHEN Mian, JIANG Hailong. A model of shear dilation induced permeability evolution in the process of water injection in Karamay continental oil sand[J]. Oil Drilling & Production Technology, 2021, 43(2): 244-249 doi:  10.13639/j.odpt.2021.02.017

克拉玛依陆相油砂注水过程中的剪胀诱导渗透率模型

doi: 10.13639/j.odpt.2021.02.017
基金项目: 国家自然科学基金项目“超深井井筒安全构建工程基础理论与方法”(编号:U1762215);陕西省自然科学基础研究计划项目“超压气藏垂直裂缝井产能来源分析和渗流机理研究”(编号:2021JM-407);陕西省教育厅专项科研计划项目“基于流固耦合的超深高压高产气井渗流特征和出砂机理研究”(编号:20JK0843)
详细信息
    作者简介:

    高彦芳(1991-),2020年毕业于中国石油大学(北京)油气井工程专业,获博士学位,现从事非常规油气储层地质力学研究工作,讲师。通讯地址:(710069)陕西省西安市太白北路西北大学地质学系。E-mail:gaoyanfang@nwu.edu.cn

  • 中图分类号: TE345

A model of shear dilation induced permeability evolution in the process of water injection in Karamay continental oil sand

  • 摘要: 特超稠油油藏普遍采用蒸汽辅助重力泄油(SAGD)方式开发,在前期注水改造阶段储层剪胀扩容导致剪切带产生剪切裂缝,孔隙度和含水/油饱和度发生变化。为了研究注水过程中水的有效渗透率的动态演化规律,考虑孔隙度和含水饱和度随扩容体积应变的变化,通过Kozeny-Poiseuille方程和Touhidi-Baghini方程关联绝对渗透率和孔隙度的关系,通过油、水两相相渗曲线方程关联水的有效渗透率和含水饱和度的关系,建立了3个剪胀诱导渗透率演化数学模型,并根据实验数据对各个模型的预测效果进行了评价分析。研究结果表明,传统的Kozeny-Carman方程只适用于孔隙性颗粒介质,不适用于克拉玛依陆相油砂这种由沥青基底式胶结形成的复合材料骨架,也不适用描述剪切扩容带内剪切裂缝引起的渗透率演化。基于Touhidi-Baghini方程建立的剪胀诱导渗透率模型,可以准确预测注水过程中水的有效渗透率演化规律,为渗透率的相关研究提供参考。
  • 图  1  油、水两相相对渗透率随含水饱和度的变化

    Figure  1.  Variation of oil and water relative permeability with water saturation

    图  2  剪胀或剪缩过程中水的有效渗透率改善程度随体积应变的变化数据

    Figure  2.  Change data of the evolution of the improvement degree of water effective permeability with volumetric strain in the process of shear dilation or shear contraction

    图  3  不同放大倍数下克拉玛依油砂微观结构

    Figure  3.  Microstructure of Karamay oil sand under different magnifications

    表  1  克拉玛依油砂和Alberta油砂的基本物理性质

    Table  1.   Basic physical properties of Karamay oil sand and Alberta oil sand

    物性参数克拉玛依油砂Alberta油砂
    密度/(103 kg · m−3) 1.96 1.75~2.40
    孔隙度/% 33 16~46
    水有效渗透率/(10−3 μm2) 2.886 0.01~1
    沥青体积分数/% 15.51 18.6
    水体积分数/% 17.49 7.9
    矿物及黏土体积分数/% 67 73.5
    20 °C原油黏度/mPa · s 5.06×106 6.42×106
    100 °C原油黏度/mPa · s 846 130
    150 °C原油黏度/mPa · s 92 32
    200 °C原油黏度/mPa · s 22 12
    下载: 导出CSV

    表  2  油砂三轴剪切前、后的水的有效渗透率改善程度

    Table  2.   Improvement degree of water effective permeability before and after the triaxial shear of oil sand

    温度/
    °C
    有效围压/
    MPa
    最大体积
    扩容量/%
    改造后水的有效
    渗透率/(10−3 μm2)
    渗透率
    改善程度
    20 0.5 −7.1 5.996 2.08
    1.0 −4.9 4.866 1.69
    2.0 −3.5 3.144 1.09
    5.0 2.7 1.996 0.69
    70 0.5 −5.9 5.274 1.83
    5.0 3.01 1.663 0.58
    下载: 导出CSV

    表  3  克拉玛依油砂储层油、水两相相渗关系拟合参数

    Table  3.   Oil and water relative permeability relation fitting parameters of Karamay oil sand reservoir

    模型参数取值
    束缚水饱和度Swr/% 29.6
    残余油饱和度Sor/% 30.6
    实验参数krw 0.079 5
    实验参数nw 1.398 7
    实验参数kro 1
    实验参数no 2.368 8
    下载: 导出CSV
  • [1] 高彦芳, 陈勉, 林伯韬, 等. 温度对油砂力学性质的影响规律研究[J]. 岩石力学与工程学报, 2018, 37(11):2520-2535. doi:  10.13722/j.cnki.jrme.2018.0635

    GAO Yanfang, CHEN Mian, LIN Botao, et al. Thermal influences on mechanical properties of oil sands[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(11): 2520-2535. doi:  10.13722/j.cnki.jrme.2018.0635
    [2] PANG H W, JIN Y, GAO Y F. Evaluation of elastic property changes in Karamay oil sand reservoir during thermal stimulation[J]. Energy Science and Engineering, 2019, 7(4): 1233-1253. doi:  10.1002/ese3.342
    [3] 高彦芳, 陈勉, 林伯韬, 等. 稠油油藏SAGD微压裂阶段储层压缩系数研究−以新疆风城陆相储层重1区齐古组为例[J]. 石油科学通报, 2017, 2(2):240-250. doi:  10.3969/j.issn.2096-1693.2017.02.022

    GAO Yanfang, CHEN Mian, LIN Botao, et al. Study on compressibility during micro-fracturing in continental ultra-heavy oil sand reservoirs-Taking the Qigu Formation of Xinjiang Fengcheng Oilfield Z1 Block for instance[J]. Petroleum Science Bulletin, 2017, 2(2): 240-250. doi:  10.3969/j.issn.2096-1693.2017.02.022
    [4] AZAD A, CHALATURNYK R J. A mathematical improvement to SAGD using geomechanical modeling[J]. Journal of Canadian Petroleum Technology, 2010, 49(10): 53-64. doi:  10.2118/141303-PA
    [5] TOUHIDI-BAGHINI A. Absolute permeability of McMurray Formation oil sands at low confining stresses[D]. Edmonton: University of Alberta, 1998.
    [6] OLDAKOWSKI K. Stress induced permeability changes of Athabasca oil sands[D]. Edmonton: University of Alberta, 1994.
    [7] YALE D P, MAYER T, WANG J. Geomechanics of oil sands under injection[R]. ARMA 10-257, 2010.
    [8] 陈森, 林伯韬, 金衍, 等. SAGD井微压裂储层渗透率变化规律研究[J]. 西南石油大学学报(自然科学版), 2018, 40(1):141-148. doi:  10.11885/j.issn.1674-5086.2016.06.21.05

    CHEN Sen, LIN Botao, JIN Yan, et al. Study on patterns of change in oil reserve permeability during microfracturing of SAGD wells[J]. Journal of Southwest Petroleum University(Science & Technology Edition), 2018, 40(1): 141-148. doi:  10.11885/j.issn.1674-5086.2016.06.21.05
    [9] WONG R C K, LI Y. A deformation-dependent model for permeability changes in oil sand due to shear dilation[J]. Journal of Canadian Petroleum Technology, 2001, 40(8): 37-44. doi:  10.2118/01-08-03
    [10] GAO Y F, PANG H W, JIN Y, et al. Evaluation of shear dilation capability/potential and permeability changes in Karamay oil sands under water injection[J]. Geofluids, 2019, 2019: 7245081. doi:  10.1155/2019/7245081
    [11] GAO Y F, CHEN M, LIN B T, et al. An analytical model of hydraulic dilation area for Karamay oil sand reservoir under water injection in SAGD wells[J]. Journal of Petroleum Science and Engineering, 2019, 179: 1090-1101. doi:  10.1016/j.petrol.2019.05.008
    [12] 刘江涛, 廖东良, 葛新民. 基于Kozeny-Carman方程的水相相对渗透率计算方法[J]. 科学技术与工程, 2012, 12(29):7500-7504. doi:  10.3969/j.issn.1671-1815.2012.29.003

    LIU Jiangtao, LIAO Dongliang, GE Xinmin. Water phase relative permeability calculation based on Kozeny-Carman equation[J]. Science Technology and Engineering, 2012, 12(29): 7500-7504. doi:  10.3969/j.issn.1671-1815.2012.29.003
    [13] 高彦芳, 陈勉, 林伯韬, 等. SAGD井挤液扩容水力波及范围模型[J]. 新疆石油地质, 2018, 39(2):202-208. doi:  10.7657/XJPG20180212

    GAO Yanfang, CHEN Mian, LIN Botao, et al. An analytical model for water swept area in SAGD wells[J]. Xinjiang Petroleum Geology, 2018, 39(2): 202-208. doi:  10.7657/XJPG20180212
    [14] GAO Y F, CHEN M, LIN B T, et al. Modeling of reservoir temperature upon preheating in SAGD wells considering phase change of bitumen[J]. International Journal of Heat and Mass Transfer, 2019, 144: 118650. doi:  10.1016/j.ijheatmasstransfer.2019.118650
    [15] 李兆敏, 王壮壮, 李松岩, 等. 温度对油砂沥青相渗规律影响研究[J]. 特种油气藏, 2015, 22(1):92-94. doi:  10.3969/j.issn.1006-6535.2015.01.020

    LI Zhaomin, WANG Zhuangzhuang, LI Songyan, et al. Research on influence of temperature on rules of relative permeability of asphalt and oil sands[J]. Special Oil & Gas Reservoirs, 2015, 22(1): 92-94. doi:  10.3969/j.issn.1006-6535.2015.01.020
  • [1] 孙涛, 孟祥娟, 王静, 黎真, 吴红军, 阿克巴尔·卡得拜, 赵静.  注水水质对裂缝性油藏储层的影响 . 石油钻采工艺, 2022, 44(2): 199-203, 210. doi: 10.13639/j.odpt.2022.02.010
    [2] 邢德钢, 程红晓, 赵长喜, 李洲, 毛为成.  春光油田低温稀油油藏泡沫+树脂固砂技术 . 石油钻采工艺, 2021, 43(2): 259-264. doi: 10.13639/j.odpt.2021.02.020
    [3] 李洪建, 张鹏, 常菁铉, 赵利庆, 兰晓龙, 王苗苗, 马文娟.  酒东油田注入水对岩心渗透率的损害率实验 . 石油钻采工艺, 2019, 41(6): 756-762. doi: 10.13639/j.odpt.2019.06.013
    [4] 李伟峰, 刘云, 于小龙, 魏浩光.  致密油储层岩石孔喉比与渗透率、孔隙度的关系 . 石油钻采工艺, 2017, 39(2): 125-129. doi: 10.13639/j.odpt.2017.02.001
    [5] 刘浩, 蔡记华, 肖长波, 王济君, 陈宇.  热处理提高煤岩渗透率的机理 . 石油钻采工艺, 2012, 34(4): 96-99.
    [6] 孙铭新.  东营北带深层砂砾岩体储层伤害机理研究 . 石油钻采工艺, 2011, 33(6): 45-48.
    [7] 贺丰果, 岳湘安, 李良川, 孙玉龙, 魏浩光, 刘怀珠.  底水油藏非均质性对水平井开采特性的影响 . 石油钻采工艺, 2009, 31(4): 63-66.
    [8] 刘鹏飞, 姜汉桥, 徐晖, 胡庆贺, 涂兴万.  缝洞型油藏开发室内模拟研究 . 石油钻采工艺, 2009, 31(5): 72-76.
    [9] 付静, 周晓君, 孙宝江, 张洪茂.  低渗裂缝性油藏渗透率的影响因素研究 . 石油钻采工艺, 2008, 30(2): 69-71.
    [10] 范学平, 吴宏.  近井地应力场及其对渗透率的伤害 . 石油钻采工艺, 2001, 23(4): 42-45,84,85. doi: 10.3969/j.issn.1000-7393.2001.04.015
    [11] 王增林, 辛林涛, 崔玉海, 何云, 唐高峰.  埕岛浅海油田注水管柱及配套工艺技术 . 石油钻采工艺, 2001, 23(3): 64-67,86. doi: 10.3969/j.issn.1000-7393.2001.03.021
    [12] 冯瑞林, 李仰民, 李凤群, 李运娥, 王进宝, 翟光华.  普通稠油油藏周期注水提高水驱效率技术研究与应用 . 石油钻采工艺, 2001, 23(1): 63-66,85. doi: 10.3969/j.issn.1000-7393.2001.01.019
    [13] 毛琼, 袁建国, 于俊吉, 田平, 张洪茂.  任丘雾迷山组油藏恢复注水数值模拟研究 . 石油钻采工艺, 2001, 23(2): 50-53,85. doi: 10.3969/j.issn.1000-7393.2001.02.017
    [14] 杜殿发, 陈月明, 汪庐山, 张代森.  底水潜山油藏出水规律研究 . 石油钻采工艺, 2000, 22(6): 42-45. doi: 10.3969/j.issn.1000-7393.2000.06.013
    [15] 黄春, 汤志强, 蒋官澄.  润湿性对恳东29块稠油油藏注水采收率的影响 . 石油钻采工艺, 1999, 21(3): 92-94. doi: 10.3969/j.issn.1000-7393.1999.03.017
    [16] 李洪建, 李海涛, 赵敏, 宫立武, 张玉涛.  矿化度梯度注水工艺新技术研究 . 石油钻采工艺, 1998, 20(6): 82-85. doi: 10.3969/j.issn.1000-7393.1998.06.024
    [17] 朱文, 朱华银.  支撑剂裂缝渗透率差异及其优质问题对压裂后经济净现值的影响 . 石油钻采工艺, 1996, 18(4): 85-89. doi: 10.3969/j.issn.1000-7393.1996.04.019
    [18] 余维初, 孙天锡, 李淑廉.  渗透率梯度测试仪的研制 . 石油钻采工艺, 1995, 17(5): 82-86,89. doi: 10.3969/j.issn.1000-7393.1995.05.019
    [19] 石步乾, 周蔚云, 王元第, 司光礼.  确定油井污染带渗透率及各种表皮因子的方法 . 石油钻采工艺, 1993, 15(5): 63-65,75. doi: 10.3969/j.issn.1000-7393.1993.05.013
    [20] 陈广.  渤南油田低渗透区块的压裂引效 . 石油钻采工艺, 1990, 12(2): 69-76. doi: 10.3969/j.issn.1000-7393.1990.02.009
  • 加载中
图(3) / 表 (3)
计量
  • 文章访问数:  6
  • HTML全文浏览量:  108
  • PDF下载量:  29
  • 被引次数: 0
出版历程
  • 修回日期:  2021-01-11
  • 网络出版日期:  2021-06-21
  • 刊出日期:  2021-06-21

克拉玛依陆相油砂注水过程中的剪胀诱导渗透率模型

doi: 10.13639/j.odpt.2021.02.017
    基金项目:  国家自然科学基金项目“超深井井筒安全构建工程基础理论与方法”(编号:U1762215);陕西省自然科学基础研究计划项目“超压气藏垂直裂缝井产能来源分析和渗流机理研究”(编号:2021JM-407);陕西省教育厅专项科研计划项目“基于流固耦合的超深高压高产气井渗流特征和出砂机理研究”(编号:20JK0843)
    作者简介:

    高彦芳(1991-),2020年毕业于中国石油大学(北京)油气井工程专业,获博士学位,现从事非常规油气储层地质力学研究工作,讲师。通讯地址:(710069)陕西省西安市太白北路西北大学地质学系。E-mail:gaoyanfang@nwu.edu.cn

  • 中图分类号: TE345

摘要: 特超稠油油藏普遍采用蒸汽辅助重力泄油(SAGD)方式开发,在前期注水改造阶段储层剪胀扩容导致剪切带产生剪切裂缝,孔隙度和含水/油饱和度发生变化。为了研究注水过程中水的有效渗透率的动态演化规律,考虑孔隙度和含水饱和度随扩容体积应变的变化,通过Kozeny-Poiseuille方程和Touhidi-Baghini方程关联绝对渗透率和孔隙度的关系,通过油、水两相相渗曲线方程关联水的有效渗透率和含水饱和度的关系,建立了3个剪胀诱导渗透率演化数学模型,并根据实验数据对各个模型的预测效果进行了评价分析。研究结果表明,传统的Kozeny-Carman方程只适用于孔隙性颗粒介质,不适用于克拉玛依陆相油砂这种由沥青基底式胶结形成的复合材料骨架,也不适用描述剪切扩容带内剪切裂缝引起的渗透率演化。基于Touhidi-Baghini方程建立的剪胀诱导渗透率模型,可以准确预测注水过程中水的有效渗透率演化规律,为渗透率的相关研究提供参考。

English Abstract

高彦芳,王晓阳,任战利,陈勉,姜海龙. 克拉玛依陆相油砂注水过程中的剪胀诱导渗透率模型[J]. 石油钻采工艺,2021,43(2):244-249 doi:  10.13639/j.odpt.2021.02.017
引用本文: 高彦芳,王晓阳,任战利,陈勉,姜海龙. 克拉玛依陆相油砂注水过程中的剪胀诱导渗透率模型[J]. 石油钻采工艺,2021,43(2):244-249 doi:  10.13639/j.odpt.2021.02.017
GAO Yanfang, WANG Xiaoyang, REN Zhanli, CHEN Mian, JIANG Hailong. A model of shear dilation induced permeability evolution in the process of water injection in Karamay continental oil sand[J]. Oil Drilling & Production Technology, 2021, 43(2): 244-249 doi:  10.13639/j.odpt.2021.02.017
Citation: GAO Yanfang, WANG Xiaoyang, REN Zhanli, CHEN Mian, JIANG Hailong. A model of shear dilation induced permeability evolution in the process of water injection in Karamay continental oil sand[J]. Oil Drilling & Production Technology, 2021, 43(2): 244-249 doi:  10.13639/j.odpt.2021.02.017
  • 与常规油气不同,特超稠油(油砂)油藏具有原油黏度极高,沥青基底式胶结,骨架疏松,具有相嵌互锁结构等特点[1-2]。对于注水或注蒸汽过程中的储层渗透率演化规律,国内外学者进行了大量研究。为了缩短循环预热周期,在SAGD技术正式实施之前,通过双水平井对向储层挤液,能够诱导井周产生剪切裂缝,增加储层渗透率,提高热对流效应[3]。在蒸汽腔发育过程中,腔外泄油区在热膨胀和孔隙压力作用下,同样能够产生剪切裂缝,诱导渗透率动态演化[4]。针对油砂剪胀过程中的体积扩容-渗透率演化问题,Touhidi-Baghini[5]测试了Athabasca油砂在低围压下进行三轴剪切过程中的体应变-绝对渗透率关系。Oldakowski[6]测试了Athabasca油砂的渗透率应力敏感性,研究了含水饱和度和水的有效渗透率之间的关系。Yale等[7]基于Alberta油砂三轴压缩实验和渗透率实验,提出了剪胀过程中的渗透率分析模型。陈森等[8]通过实验发现,克拉玛依油砂的绝对渗透率和水的有效渗透率随孔隙度或体应变的变化基本遵循Kozeny-Poiseuille方程。Wong等[9]通过颗粒堆积理论,建立了骨架颗粒在弹性变形和剪胀过程中的渗透率和应力/应变之间的数学关系。Gao等 [10]提出剪胀扩容和渗透率改善潜力的概念,描述注水过程中的孔隙度和渗透率演化,并基于颗粒堆积理论及克拉玛依油砂微观结构建立数学预测模型;同时也给出了克拉玛依油砂三轴剪切实验过程中,水的有效渗透率随有效围压和体应变的变化规律。但目前剪胀诱导渗透率实验普遍采用海相Alberta油砂,渗透率模型普遍假设油砂为颗粒堆积组成的孔隙性介质,忽略了体积扩容过程中的含水/油饱和度变化。而克拉玛依油砂为未经冰川压实作用的陆相油砂,与海相油砂孔隙结构明显不同[11]。在实际注水过程中,沥青基底式胶结油砂剪胀导致剪切带产生剪切裂缝,孔隙性岩石的Kozeny-Carman方程不再适用[12]。另外,注水扩容过程中的含水/油饱和度变化将影响两相有效渗透率的大小[7]。因此,有必要建立考虑剪胀扩容的裂缝性岩石有效渗透率预测模型,分析注水扩容过程中体积应变和水的有效渗透率之间的耦合关系。

    研究考虑剪胀扩容过程中孔隙度和含水饱和度随体积应变的变化,通过Kozeny-Poiseuille方程和Touhidi-Baghini方程关联绝对渗透率和孔隙度的关系,通过油、水两相相渗曲线方程关联水的有效渗透率和含水饱和度的关系,建立了3个剪胀诱导渗透率演化数学模型,并根据实验数据对各个模型的预测效果进行了评价。

    • 假设油砂剪胀前的代表性体积单元(RVE)体积为V0,初始孔隙体积为Vp0;注水剪胀扩容后,RVE体积变为V,孔隙体积为Vp

      假设基质颗粒不可压缩,则RVE体积增加量等于孔隙体积增加量。RVE膨胀体应变εv写为

      $${\varepsilon _{\rm{v}}} = \frac{{V - {V_0}}}{{{V_0}}} = \frac{{{V_{\rm{p}}} - {V_{{\rm{p0}}}}}}{{{V_0}}}$$ (1)

      岩石剪胀后的孔隙度$\phi $写为

      $$\phi = \frac{{{V_{\rm{p}}}}}{V} = \frac{{{V_{{\rm{p0}}}} + {\varepsilon _{\rm{v}}}{V_0}}}{{{V_{\rm{0}}} + {\varepsilon _{\rm{v}}}{V_0}}} = \frac{{{\phi _0}{V_{\rm{0}}} + {\varepsilon _{\rm{v}}}{V_0}}}{{{V_{\rm{0}}} + {\varepsilon _{\rm{v}}}{V_0}}} = \frac{{{\phi _0} + {\varepsilon _{\rm{v}}}}}{{1 + {\varepsilon _{\rm{v}}}}}$$ (2)

      由Kozeny-Carman方程可知[8],孔隙性岩石的绝对渗透率ka与孔隙度$\phi $的关系式为

      $${k_{\rm{a}}} = \frac{{C{\phi ^3}}}{{{{\left( {1 - \phi } \right)}^2}}} = \frac{{C{{\left( {{\phi _0} + {\varepsilon _{\rm{v}}}} \right)}^3}}}{{\left( {1 + {\varepsilon _{\rm{v}}}} \right){{\left( {1 - {\phi _0}} \right)}^2}}}$$ (3)

      假设C为不随孔隙度变化而变化的常数,则初始绝对渗透率ka0

      $${k_{{\rm{a0}}}} = \frac{{C\phi _0}^3}{{{{\left( {1 - {\phi _0}} \right)}^2}}}$$ (4)

      式(3)除以式(4),可得绝对渗透率随体积应变的演化关系为

      $$\frac{{{k_{\rm{a}}}}}{{{k_{{\rm{a0}}}}}} = \frac{{{{\left( {{\phi _0} + {\varepsilon _{\rm{v}}}} \right)}^3}}}{{{\phi _0}^3\left( {1 + {\varepsilon _{\rm{v}}}} \right)}}$$ (5)

      油砂注水剪胀的微观力学机理为,砂粒的翻滚带动粒间黏土和沥青的挤压和错动,在剪切扩容带产生剪切裂缝[3]。对于裂缝性介质,采用Kozeny-Carman方程描述其渗透率-孔隙度关系时可能存在一定误差。Touhidi-Baghini[5]通过大量实验研究,认为油砂剪胀过程中的绝对渗透率与体应变的关系可以表示为

      $$\frac{{{k_{\rm{a}}}}}{{{k_{{\rm{a}}0}}}} = {{\rm{e}}^{n{\varepsilon _{\rm{v}}}}}$$ (6)

      式中,${\phi _0}$为初始孔隙度(沥青视为孔隙流体);C为实验参数;n为实验参数。

      式(5)和式(6)为注水过程中考虑剪胀体积扩容的绝对渗透率演化方程的2种形式。

      假设含气饱和度忽略不计,则注水过程中的RVE孔隙体积增加量等于水体积增加量。含水饱和度Sw写为

      $${S_{\rm{w}}} = \frac{{{V_{\rm{w}}}}}{{{V_{\rm{p}}}}} = \frac{{{V_{{\rm{p0}}}}{S_{{\rm{w0}}}} + {V_0}{\varepsilon _{\rm{v}}}}}{{{V_{{\rm{p0}}}} + {V_0}{\varepsilon _{\rm{v}}}}} = \frac{{{V_{\rm{0}}}{\phi _0}{S_{{\rm{w0}}}} + {V_0}{\varepsilon _{\rm{v}}}}}{{{V_{\rm{0}}}{\phi _0} + {V_0}{\varepsilon _{\rm{v}}}}} = \frac{{{\phi _0}{S_{{\rm{w0}}}} + {\varepsilon _{\rm{v}}}}}{{{\phi _0} + {\varepsilon _{\rm{v}}}}}$$ (7)

      剪胀扩容过程中,水的相对渗透率krw和水的有效渗透率kw随含水饱和度Sw的变化写为[7]

      $${k_{{\rm{rw}}}} = \frac{{{k_{\rm{w}}}}}{{{k_{\rm{a}}}}} = {k_{{\rm{rw}}}}^\prime {\left( {\frac{{{S_{\rm{w}}} - {S_{{\rm{wr}}}}}}{{1 - {S_{{\rm{wr}}}} - {S_{{\rm{or}}}}}}} \right)^{{n_{\rm{w}}}}}$$ (8)

      同理,油的相对渗透率kro和油的有效渗透率ko随含水饱和度Sw的变化写为

      $${k_{{\rm{ro}}}} = \frac{{{k_{\rm{o}}}}}{{{k_{\rm{a}}}}} = {k_{{\rm{ro}}}}^\prime{\left( {\frac{{1 - {S_{\rm{w}}} - {S_{{\rm{or}}}}}}{{1 - {S_{{\rm{wr}}}} - {S_{{\rm{or}}}}}}} \right)^{{n_{\rm{o}}}}}$$ (9)

      式中,Vw为剪胀后的孔隙水体积,m3Sw0为剪胀前的初始含水饱和度;${ k_{{\rm{rw}}}}^\prime$nw为实验参数;Swr为束缚水饱和度;Sor为残余油饱和度;${k_{{\rm{ro}}}}^\prime$no为实验参数。

      联立式(5)和式(8)或联立式(6)和式(8),可得水的有效渗透率与膨胀体应变的关系为

      $${k_{\rm{w}}} = {k_{{\rm{a0}}}}\frac{{{{\left( {{\phi _0} + {\varepsilon _{\rm{v}}}} \right)}^3}}}{{{\phi _0}^3\left( {1 + {\varepsilon _{\rm{v}}}} \right)}}{k_{{\rm{rw}}}}^\prime{\left( {\dfrac{{\dfrac{{{\phi _0}{S_{{\rm{w0}}}} + {\varepsilon _{\rm{v}}}}}{{{\phi _0} + {\varepsilon _{\rm{v}}}}} - {S_{{\rm{wr}}}}}}{{1 - {S_{{\rm{wr}}}} - {S_{{\rm{or}}}}}}} \right)^{{n_{\rm{w}}}}}$$ (10)

      $${k_{\rm{w}}} = {k_{{\rm{a0}}}}{{\rm{e}}^{n{\varepsilon _{\rm{v}}}}}{k_{{\rm{rw}}}}^\prime{\left( {\dfrac{{\dfrac{{{\phi _0}{S_{{\rm{w0}}}} + {\varepsilon _{\rm{v}}}}}{{{\phi _0} + {\varepsilon _{\rm{v}}}}} - {S_{{\rm{wr}}}}}}{{1 - {S_{{\rm{wr}}}} - {S_{{\rm{or}}}}}}} \right)^{{n_{\rm{w}}}}}$$ (11)

      εv=0,式(10)和式(11)写为

      $${k_{{\rm{w0}}}} = {k_{{\rm{a0}}}}{k_{{\rm{rw}}}}^\prime{\left( {\frac{{{S_{{\rm{w0}}}} - {S_{{\rm{wr}}}}}}{{1 - {S_{{\rm{wr}}}} - {S_{{\rm{or}}}}}}} \right)^{{n_{\rm{w}}}}}$$ (12)

      式(10)除以式(12),可得

      $$\frac{{{k_{\rm{w}}}}}{{{k_{{\rm{w0}}}}}} = \frac{{{{\left( {{\phi _0} + {\varepsilon _{\rm{v}}}} \right)}^3}}}{{{\phi _0}^3\left( {1 + {\varepsilon _{\rm{v}}}} \right)}}{\left( {\dfrac{{\dfrac{{{\phi _0}{S_{{\rm{w0}}}} + {\varepsilon _{\rm{v}}}}}{{{\phi _0} + {\varepsilon _{\rm{v}}}}} - {S_{{\rm{wr}}}}}}{{{S_{{\rm{w0}}}} - {S_{{\rm{wr}}}}}}} \right)^{{n_{\rm{w}}}}}$$ (13)

      式(11)除以式(12),可得

      $$\frac{{{k_{\rm{w}}}}}{{{k_{{\rm{w0}}}}}} = {{\rm{e}}^{n{\varepsilon _{\rm{v}}}}}{\left( {\dfrac{{\dfrac{{{\phi _0}{S_{{\rm{w0}}}} + {\varepsilon _{\rm{v}}}}}{{{\phi _0} + {\varepsilon _{\rm{v}}}}} - {S_{{\rm{wr}}}}}}{{{S_{{\rm{w0}}}} - {S_{{\rm{wr}}}}}}} \right)^{{n_{\rm{w}}}}}$$ (14)

      式(13)和式(14)为注水过程中考虑剪胀体积扩容的水的有效渗透率演化方程的2种形式。

      在注水阶段,沥青处于不流动状态,起到胶结物的作用[13]。若将沥青视为骨架,则孔隙中只含有水相。根据式(5)可以估算水的有效渗透率为

      $$\frac{{{k_{\rm{w}}}}}{{{k_{{\rm{w0}}}}}} = \frac{{{{\left( {{\phi _1} + {\varepsilon _{\rm{v}}}} \right)}^3}}}{{{\phi _1}^3\left( {1 + {\varepsilon _{\rm{v}}}} \right)}}$$ (15)

      式中,kw0为油砂剪胀扩容前的初始水的有效渗透率,10−3 μm2${\phi _1}$为不包括沥青在内的孔隙度。

      式(2)和式(15)是不矛盾的。式(2)中的孔隙度用来计算绝对渗透率,而式(15)中的不包括沥青在内的孔隙度用来计算水的有效渗透率。在将油砂视为孔隙性介质且把沥青视为骨架的条件下,根据式(5)可以近似估计水的有效渗透率。虽然这种方法比较简单,但是结果比较粗略。陈森等[8]即采用了这种方法粗略估计渗透率的演化规律。

    • 采用GCTS RTR-1500型高温高压岩石三轴仪和瞬时脉冲渗透率仪,测试油砂在剪切过程中的体积应变和水的有效渗透率。取风城油田重1区某井下储层油砂样品,实验温度为20~70 ℃,孔隙压力为5 MPa,围压为5.5~10 MPa,最大加载轴向应变为9%。表1为该区块储层与Alberta油砂储层基本物理及力学性质对比[14]。由表可知,与Alberta油砂相比,克拉玛依油砂沥青黏度、孔隙度和水的有效渗透率较高,而储层密度和含油饱和度较低。克拉玛依油砂包含沥青在内的初始孔隙度为33%,不包括沥青的初始孔隙度为16%,初始含油饱和度47%。

      表 1  克拉玛依油砂和Alberta油砂的基本物理性质

      Table 1.  Basic physical properties of Karamay oil sand and Alberta oil sand

      物性参数克拉玛依油砂Alberta油砂
      密度/(103 kg · m−3) 1.96 1.75~2.40
      孔隙度/% 33 16~46
      水有效渗透率/(10−3 μm2) 2.886 0.01~1
      沥青体积分数/% 15.51 18.6
      水体积分数/% 17.49 7.9
      矿物及黏土体积分数/% 67 73.5
      20 °C原油黏度/mPa · s 5.06×106 6.42×106
      100 °C原油黏度/mPa · s 846 130
      150 °C原油黏度/mPa · s 92 32
      200 °C原油黏度/mPa · s 22 12

      实验结果表明,当温度为20~70 ℃时,在5 MPa有效围压下只发生了剪缩,剪缩量约为2%~3%;在0.5~2 MPa有效围压下先发生了剪缩,后发生了明显的剪胀,剪胀量约为3%~7%。随着有效围压逐渐降低,剪胀体积应变增加,水的有效渗透率增加。温度对剪胀体积的影响不大。表2列出了不同温度和有效围压下三轴剪切实验结束后,最终体积扩容(或压缩)量和有效渗透率改善(或降低)程度的关系。渗透率改善程度为剪胀(或剪缩)后的渗透率与初始渗透率之比[10]。由表可知,在20 ℃和0.5 MPa有效围压下,最大体积扩容量为7.1% (负值表示体积膨胀),最大渗透率改善程度为2.08倍;在70 ℃和5 MPa有效围压下,最大体积压缩量为3.01% (正值表示体积压缩),最大渗透率降低倍数为0.58倍。

      表 2  油砂三轴剪切前、后的水的有效渗透率改善程度

      Table 2.  Improvement degree of water effective permeability before and after the triaxial shear of oil sand

      温度/
      °C
      有效围压/
      MPa
      最大体积
      扩容量/%
      改造后水的有效
      渗透率/(10−3 μm2)
      渗透率
      改善程度
      20 0.5 −7.1 5.996 2.08
      1.0 −4.9 4.866 1.69
      2.0 −3.5 3.144 1.09
      5.0 2.7 1.996 0.69
      70 0.5 −5.9 5.274 1.83
      5.0 3.01 1.663 0.58
    • 将式(14)、式(13)和式(15)描述的模型分别称为模型一、模型二和模型三。利用3种模型预测油砂注水剪胀扩容过程中水的有效渗透率随体积应变的动态演化规律,并结合实验室实测数据,对3种模型的预测效果进行了评价分析。

      图1所示散点为室内实验测得的研究区块油藏在200 ℃热水驱条件下的油、水两相相渗关系。绝对渗透率是岩石本身的一种属性,不随通过其中的流体性质而变化,因此温度对绝对渗透率的影响可以忽略。不同温度下,水的相对渗透率随含水饱和度的变化基本相同[15]。因此,采用200 ℃热水驱条件下的油、水两相相渗结果可以近似描述20~70 ℃条件下的水相有效渗透率的动态变化。

      图  1  油、水两相相对渗透率随含水饱和度的变化

      Figure 1.  Variation of oil and water relative permeability with water saturation

      根据式(8)和式(9)拟合图1中的实验数据,计算得到的拟合参数如表3所示。将表中的模型参数分别代入到式(14)、式(13)和式(15)进行计算,根据3种模型预测油砂剪缩或剪胀过程中水的有效渗透率的动态演化规律,如图2所示。由图可知,体积应变从−0.04 (体积压缩)到0.08 (体积膨胀)变化时,3种模型预测的水的有效渗透率改善程度均逐渐增加。与实测数据对比,模型一的预测效果最好,模型二和模型三的预测结果均偏大。因此,建议采用模型一预测陆相克拉玛依油砂剪胀扩容过程中的水的有效渗透率演化过程。

      图3(a)所示为克拉玛依油砂原始微观结构,可知克拉玛依油砂疏松程度高,颗粒与颗粒之间的接触点/面稀少,粒间充填大量的沥青和黏土混合物。如图3(b)所示为克拉玛依油砂在20 ℃和0.5 MPa有效围压下发生剪胀后,不同放大倍数下的剪切扩容带的微观结构。由图可知,剪切带发育明显,角砾状砂粒显著翻转,形成粒间大孔隙,增加孔隙度并形成优势渗流通道。油砂注水剪胀的微观力学机理为,砂粒的翻滚带动粒间黏土和沥青的挤压和错动,在剪切扩容带产生了剪切裂缝[10]。研究中建立的是渗流参数与体积变形之间的关系,沥青的变形主要为剪应变,因此沥青的变形是不需要考虑的。

      表 3  克拉玛依油砂储层油、水两相相渗关系拟合参数

      Table 3.  Oil and water relative permeability relation fitting parameters of Karamay oil sand reservoir

      模型参数取值
      束缚水饱和度Swr/% 29.6
      残余油饱和度Sor/% 30.6
      实验参数krw 0.079 5
      实验参数nw 1.398 7
      实验参数kro 1
      实验参数no 2.368 8

      图  2  剪胀或剪缩过程中水的有效渗透率改善程度随体积应变的变化数据

      Figure 2.  Change data of the evolution of the improvement degree of water effective permeability with volumetric strain in the process of shear dilation or shear contraction

      图  3  不同放大倍数下克拉玛依油砂微观结构

      Figure 3.  Microstructure of Karamay oil sand under different magnifications

      模型二预测结果偏大的原因:Kozeny-Carman方程适用于孔隙性颗粒介质[12],不适用于这种由沥青基底式胶结形成的复合材料骨架,也不适用描述剪切扩容带的剪切裂缝引起的渗透率演化。

      模型三预测结果较大的原因:除模型二提到的原因外,模型三将油、水两相流体简化为单相流体(水),没有考虑油和水之间的相互作用。由于水的有效渗透率和油的有效渗透率之和一定小于绝对渗透率,因此这种简化产生了很大的误差。

    • (1)克拉玛依陆相油砂注水剪胀的微观力学机理为:在低有效围压下砂粒的翻滚带动粒间黏土和沥青的挤压和错动,在剪切扩容带产生剪切裂缝,增加孔隙度并形成优势渗流通道。

      (2)传统的Kozeny-Carman方程只适用于孔隙性颗粒介质,不适用于克拉玛依陆相油砂这种由沥青基底式胶结形成的复合材料骨架,也不适用描述剪切裂缝引起的渗透率演化。

      (3)考虑孔隙度和含水饱和度随体积应变的变化,通过Touhidi-Baghini方程关联绝对渗透率和孔隙度的关系,通过油、水两相相渗曲线方程关联水的有效渗透率和含水饱和度的关系,建立克拉玛依油砂剪胀诱导渗透率模型,可以准确预测注水过程中水的有效渗透率演化规律。

参考文献 (15)

目录

    /

    返回文章
    返回