致密储层人工裂缝与天然裂缝夹角对水驱油效果的影响

王晓明 陈军斌 孙晨 柳文欣

王晓明,陈军斌,孙晨,柳文欣. 致密储层人工裂缝与天然裂缝夹角对水驱油效果的影响[J]. 石油钻采工艺,2022,44(1):70-76 doi:  10.13639/j.odpt.2022.01.013
引用本文: 王晓明,陈军斌,孙晨,柳文欣. 致密储层人工裂缝与天然裂缝夹角对水驱油效果的影响[J]. 石油钻采工艺,2022,44(1):70-76 doi:  10.13639/j.odpt.2022.01.013
WANG Xiaoming, CHEN Junbin, SUN Chen, LIU Wenxin. Effect of the angle between hydraulic fracture and natural fracture on water flooding in tight reservoirs[J]. Oil Drilling & Production Technology, 2022, 44(1): 70-76 doi:  10.13639/j.odpt.2022.01.013
Citation: WANG Xiaoming, CHEN Junbin, SUN Chen, LIU Wenxin. Effect of the angle between hydraulic fracture and natural fracture on water flooding in tight reservoirs[J]. Oil Drilling & Production Technology, 2022, 44(1): 70-76 doi:  10.13639/j.odpt.2022.01.013

致密储层人工裂缝与天然裂缝夹角对水驱油效果的影响

doi: 10.13639/j.odpt.2022.01.013
基金项目: 国家自然科学基金“陆相页岩井周天然裂隙力学活动性评价方法基础研究”(编号:51874239)
详细信息
    作者简介:

    王晓明(1990-),博士研究生,从事油气渗流理论与提高采收率技术。通讯地址:(710054)陕西省西安市雁塔路58号西安科技大学地质与环境学院。E-mail: wxm18392177016@sina.com

  • 中图分类号: TE357.2

Effect of the angle between hydraulic fracture and natural fracture on water flooding in tight reservoirs

  • 摘要: 为研究人工裂缝与天然裂缝夹角对水驱油效果的影响,考虑重力和裂缝对水驱油效果的影响,建立了非活塞式水驱油油水两相单向渗流模型,利用有限元分析方法,通过产量和含油率定量分析人工裂缝与天然裂缝夹角处于0°~90°时对水驱油效果的影响。研究表明:人工裂缝与天然裂缝夹角越大,油井产量越高,当人工裂缝与天然裂缝夹角为90°时,油井累积产量最高,且油井产量递减速度最小;人工裂缝与天然裂缝夹角增大,扩大了注入水波及基质系统的范围,驱油效果显著,当人工裂缝与天然裂缝夹角为90°时,裂缝系统和基质系统的驱油效果显著优于其他角度的裂缝分布,注入水波及范围最广,剩余油最少,但水驱油效率整体较低。该研究结果有助于认识裂缝分布对于水驱开发效果的影响,并预防水窜的发生。
  • 图  1  双重孔隙介质物理模型及非活塞式水驱油单向流模型

    Figure  1.  Physical model of dual porous media and one-way flow model of non piston water flooding

    图  2  网格剖分结果

    Figure  2.  Mesh generation results

    图  3  不同裂缝夹角下储层累积产量随时间变化

    Figure  3.  Variation of reservoir cumulative production with time under different fracture angles

    图  4  第2 000 d储层含油率分布

    Figure  4.  Distribution of oil-bearing rate of the 2 000 d reservoir

  • [1] 赵阳, 曲志浩, 刘震. 裂缝水驱油机理的真实砂岩微观模型实验研究[J]. 石油勘探与开发, 2002, 29(1):116-119. doi:  10.3321/j.issn:1000-0747.2002.01.031

    ZHAO Yang, QU Zhihao, LIU Zhen. Experimental study on water/oil displacement mechanisms in fractured reservoir by real sand stone micro-models[J]. Petroleum Exploration and Development, 2002, 29(1): 116-119. doi:  10.3321/j.issn:1000-0747.2002.01.031
    [2] 时佃海, 陈光梅, 邢正岩, 等. 牛20低渗断块裂缝的分布及对开发效果的影响[J]. 大庆石油学院学报, 2004, 28(3):10-13,117.

    SHI Dianhai, CHEN Guangmei, XING Zhengyan, et al. Fracture distribution and effect in Niu 20 low permeability faulted block on exploitation efficiency[J]. Journal of Northeast Petroleum University, 2004, 28(3): 10-13,117.
    [3] 刘建军, 刘先贵. 裂缝分布对水驱油效果的影响[J]. 新疆石油地质, 2002, 23(2):146-147,86. doi:  10.3969/j.issn.1001-3873.2002.02.019

    LIU Jianjun, LIU Xiangui. Study of numerical simulation of rational well patterns for development of fractured oilfield[J]. Xinjiang Petroleum Geology, 2002, 23(2): 146-147,86. doi:  10.3969/j.issn.1001-3873.2002.02.019
    [4] 周娟, 薛惠, 郑德温, 等. 裂缝油藏水驱油渗流机理[J]. 重庆大学学报(自然科学版), 2000, 23(增刊1):65-67.

    ZHOU Juan, XUE Hui, ZHENG Dewen, et al. The mechanisms of water displacing oil in fractured reservoir[J]. Journal of Chongqing University, 2000, 23(S1): 65-67.
    [5] ZHANG R H, ZHANG L H, LUO J X, et al. Numerical simulation of water flooding in natural fractured reservoirs based on control volume finite element method[J]. Journal of Petroleum Science and Engineering, 2016(146): 1211-1225.
    [6] 吕文雅, 曾联波, 刘静, 等. 致密低渗透储层裂缝研究进展[J]. 地质科技情报, 2016, 35(4):74-83.

    LV Wenya, ZENG Lianbo, LIU Jing, et al. Fractured research progress in low permeability tight reservoir[J]. Bulletin of Geological Science and Technology, 2016, 35(4): 74-83.
    [7] JON K, RICHARD A M. Oil production from fractured reservoirs by water displacement[R]. SPE 5084, 1974.
    [8] GU S H, LIU Y T, CHEN X, et al. A method for evaluation of water flooding performance in fractured reservoirs[J]. Journal of Petroleum Science and Engineering, 2014(120): 130-140.
    [9] YUAN B, XU C, WANG K, et al. Enhance horizontal well performance by optimising multistage hydraulic fracture and water flooding[J]. Int. J. of Oil, Gas and Coal Technology, 2017, 15(1): 25-46. doi:  10.1504/IJOGCT.2017.083857
    [10] 陈守雨, 李文权, 杜林麟, 等. 水力裂缝与天然裂缝相交作用准则研究[J]. 中外能源, 2014, 19(10):37-43.

    CHEN Shouyu, LI Wenquan, DU Linlin, et al. Research on the acting criteria when hydraulic fracture intersects with natural fracture[J]. Sino-Global Energy, 2014, 19(10): 37-43.
    [11] BLANTON T L. Propagation of hydraulically and dynamically induced fractures in naturally fractured reservoirs[R]. SPE 15261, 1986.
    [12] 李璗, 陈军斌. 油气渗流力学[M]. 1版. 北京: 石油工业出版社, 2011: 126-133.

    LI Dang, CHEN Junbin. Mechanics of oil and gas flow in porous media[M]. 1st Edition. Beijing: Petroleum Industry Press, 2011: 126-133.
    [13] 亓倩, 朱维耀. 复杂压裂缝网页岩气储层压力传播动边界研究[J]. 工程科学学报, 2019, 41(11):1387-1396.

    QI Qian, ZHU Weiyao. Moving boundary analysis of fractured shale gas reservoir[J]. Chinese Journal of Engineering, 2019, 41(11): 1387-1396.
    [14] 康永梅, 李明, 王联国, 等. 鄂尔多斯盆地演武油田S1区块延8油层组微观水驱渗流特征[J]. 油气地质与采收率, 2019, 26(5):48-57.

    KANG Yongmei, LI Ming, WANG Lianguo, et al. Microscopic water flooding percolation characteristics of Yan 8 Formation in Block S1, Yanwu Oilfield, Ordos Basin[J]. Petroleum Geology and Recovery Efficiency, 2019, 26(5): 48-57.
    [15] 叶仲并. 提高采收率原理[M]. 2版. 北京: 石油工业出版社, 2013: 34-36.

    YE Zhongbing. Principles of EOR[M]. 2nd Edition. Beijing: Petroleum Industry Press, 2013: 34-36.
  • [1] 刘学伟, 田福春, 李东平, 杨立永, 阴启武, 余嵚鑫, 闫伟.  沧东凹陷孔二段页岩压裂套变原因分析及预防对策 . 石油钻采工艺, 2022, 44(1): 77-82. doi: 10.13639/j.odpt.2022.01.011
    [2] 孙业恒, 殷代印, 刘凯.  致密油藏不同沉积岩人工裂缝导流能力实验 . 石油钻采工艺, 2021, 43(1): 104-109. doi: 10.13639/j.odpt.2021.01.016
    [3] 许建国, 刘光玉, 王艳玲.  致密储层缝内暂堵转向压裂工艺技术 . 石油钻采工艺, 2021, 43(3): 374-378. doi: 10.13639/j.odpt.2021.03.015
    [4] 肖阳, 张丽英, 邓虎成, 邓勇, 何文, 王盾.  水力裂缝及天然裂缝相交延伸规律及真三轴物理模拟实验 . 石油钻采工艺, 2021, 43(1): 83-88. doi: 10.13639/j.odpt.2021.01.013
    [5] 郭建春, 任冀川, 苟波, 王世彬.  厚层强非均质砂岩储层酸液立体放置优化设计 . 石油钻采工艺, 2019, 41(5): 649-656. doi: 10.13639/j.odpt.2019.05.016
    [6] 叶亮, 邹雨时, 赵倩云, 李四海, 丁勇, 马新星.  致密砂岩储层CO2压裂裂缝扩展实验研究 . 石油钻采工艺, 2018, 40(3): 361-368. doi: 10.13639/j.odpt.2018.03.015
    [7] 朱金智, 游利军, 张震, 康毅力, 徐三峰, 林冲.  聚磺混油钻井液对深层裂缝性致密储层的保护能力评价 . 石油钻采工艺, 2018, 40(3): 311-317. doi: 10.13639/j.odpt.2018.03.007
    [8] 张健, 张国祥, 邹雨时, 赵文韬, 王金意.  致密储层CO2压裂裂缝扩展规律数值模拟 . 石油钻采工艺, 2018, 40(3): 354-360. doi: 10.13639/j.odpt.2018.03.014
    [9] 乔领良.  “NEW-DRILL”提速新工具在元坝气田的应用 . 石油钻采工艺, 2016, 38(1): 14-17. doi: 10.13639/j.odpt.2016.01.004
    [10] 刘立峰, 冉启全, 王欣, 李冉.  致密储层水平井体积压裂段间距优化方法 . 石油钻采工艺, 2015, 37(3): 84-87. doi: 10.13639/j.odpt.2015.03.019
    [11] 崔明月, 刘玉章, 修乃领, 丁云宏, 鄢雪梅, 窦晶晶, 付海峰.  形成复杂缝网体积(ESRV)的影响因素分析 . 石油钻采工艺, 2014, 36(2): 82-87. doi: 10.13639/j.odpt.2014.02.021
    [12] 王晓东, 赵振峰, 李向平, 马兵, 杨晓刚, 闫永萍.  鄂尔多斯盆地致密油层混合水压裂试验   . 石油钻采工艺, 2012, 34(5): 80-84.
    [13] 吴奇, 胥云, 刘玉章, 丁云宏, 王晓泉, 王腾飞.  美国页岩气体积改造技术现状及对我国的启示 . 石油钻采工艺, 2011, 33(2): 1-7.
    [14] 吴小平, 陶良军, 张雪莲, 李伟, 张仁俊, 余芳.  双重介质低渗透储层压裂技术研究与应用 . 石油钻采工艺, 2010, 32(4): 128-130.
    [15] 谢润成, 韩振林, 周文, 张生军, 葛善良, 王辛.  深层致密砂岩气藏高异常破裂压力影响因素分析 . 石油钻采工艺, 2009, 31(5): 60-64.
    [16] 姜晶, 李春兰, 杨敏.  低渗透油藏压裂水平井裂缝优化研究 . 石油钻采工艺, 2008, 30(4): 50-52.
    [17] 李中锋, 何顺利, 门成全, 杨文新.  非均质三维模型水驱剩余油试验研究 . 石油钻采工艺, 2005, 27(4): 41-44. doi: 10.3969/j.issn.1000-7393.2005.04.013
    [18] 何顺利, 李中锋, 杨文新, 门成全.  非均质线性模型水驱油试验研究 . 石油钻采工艺, 2005, 27(5): 49-52. doi: 10.3969/j.issn.1000-7393.2005.05.014
    [19] 李勇明, 赵金洲, 郭建春, 崔慧芝, 王万迅.  裂缝性低渗透储层压裂液滤失计算新模型 . 石油钻采工艺, 2004, 26(5): 44-46,49. doi: 10.3969/j.issn.1000-7393.2004.05.014
    [20] 宋洪才.  油层渗透率及厚度对聚合物驱效果的影响 . 石油钻采工艺, 1999, 21(1): 100-101. doi: 10.3969/j.issn.1000-7393.1999.01.025
  • 加载中
图(4)
计量
  • 文章访问数:  0
  • HTML全文浏览量:  57
  • PDF下载量:  32
  • 被引次数: 0
出版历程
  • 修回日期:  2021-12-05
  • 网络出版日期:  2022-07-06
  • 刊出日期:  2022-07-06

致密储层人工裂缝与天然裂缝夹角对水驱油效果的影响

doi: 10.13639/j.odpt.2022.01.013
    基金项目:  国家自然科学基金“陆相页岩井周天然裂隙力学活动性评价方法基础研究”(编号:51874239)
    作者简介:

    王晓明(1990-),博士研究生,从事油气渗流理论与提高采收率技术。通讯地址:(710054)陕西省西安市雁塔路58号西安科技大学地质与环境学院。E-mail: wxm18392177016@sina.com

  • 中图分类号: TE357.2

摘要: 为研究人工裂缝与天然裂缝夹角对水驱油效果的影响,考虑重力和裂缝对水驱油效果的影响,建立了非活塞式水驱油油水两相单向渗流模型,利用有限元分析方法,通过产量和含油率定量分析人工裂缝与天然裂缝夹角处于0°~90°时对水驱油效果的影响。研究表明:人工裂缝与天然裂缝夹角越大,油井产量越高,当人工裂缝与天然裂缝夹角为90°时,油井累积产量最高,且油井产量递减速度最小;人工裂缝与天然裂缝夹角增大,扩大了注入水波及基质系统的范围,驱油效果显著,当人工裂缝与天然裂缝夹角为90°时,裂缝系统和基质系统的驱油效果显著优于其他角度的裂缝分布,注入水波及范围最广,剩余油最少,但水驱油效率整体较低。该研究结果有助于认识裂缝分布对于水驱开发效果的影响,并预防水窜的发生。

English Abstract

王晓明,陈军斌,孙晨,柳文欣. 致密储层人工裂缝与天然裂缝夹角对水驱油效果的影响[J]. 石油钻采工艺,2022,44(1):70-76 doi:  10.13639/j.odpt.2022.01.013
引用本文: 王晓明,陈军斌,孙晨,柳文欣. 致密储层人工裂缝与天然裂缝夹角对水驱油效果的影响[J]. 石油钻采工艺,2022,44(1):70-76 doi:  10.13639/j.odpt.2022.01.013
WANG Xiaoming, CHEN Junbin, SUN Chen, LIU Wenxin. Effect of the angle between hydraulic fracture and natural fracture on water flooding in tight reservoirs[J]. Oil Drilling & Production Technology, 2022, 44(1): 70-76 doi:  10.13639/j.odpt.2022.01.013
Citation: WANG Xiaoming, CHEN Junbin, SUN Chen, LIU Wenxin. Effect of the angle between hydraulic fracture and natural fracture on water flooding in tight reservoirs[J]. Oil Drilling & Production Technology, 2022, 44(1): 70-76 doi:  10.13639/j.odpt.2022.01.013
    • 致密低渗透储层天然裂缝发育,天然裂缝的存在不仅会改变储层流体的渗流路径,还决定着储层注水开发效果。对于前者,赵阳等[1]采用真实砂岩模型进行水驱替实验,并发现在裂缝型储层水驱油过程中,注入水将以多条线路向出口端非均匀推进,即发生非活塞式水驱油。对于后者,时佃海等[2]通过矿场实验研究了天然裂缝对水驱开发效果的影响,并得出在油水井间的裂缝会造成部分油井含水上升快,甚至水淹,从而影响注水效果。为进一步研究裂缝如何影响水驱效果,刘建军等[3]采用薄板模型物理模拟裂缝方向与流动方向平行和垂直的水驱油效果,并得出与渗流方向平行的裂缝容易导致注入水突破,见水时间早,采收率低;而与渗流方向垂直的裂缝能够使无法运移的油活动起来,采收率高。为细化研究裂缝方向与流动方向关系对水驱油效果的影响,周娟等[4]采用二维玻璃微模型研究了裂缝与流动方向夹角为0°~90°时的水驱油效果,并得出剩余油会随裂缝与主流线角度的增加而减小,水驱油效果则逐渐变好。Zhang R H等[5]采用数值模拟方法研究了裂缝方向与流动方向夹角为0°和±45°时的水驱效果,得出裂缝方向与流动方向夹角为−45°时采收率最高,储层含水率最低,驱油效果最好。因此,天然裂缝的存在及其方向是影响水驱开发效果的重要因素。

      尽管致密低渗透储层发育了一些天然微裂缝,并且单条微裂缝的渗透率大于基质渗透率,但由于其仅占储层总体积的极小部分,原始状态下对储层整体渗透率的贡献几乎等同于基质部分,在开发过程中表现为油井不压裂就无自然产能或自然产能低[6-7]。因此,对致密低渗透储层仍需要进行压裂开采。同时,全球采油经验表明:裂缝发育型油藏水驱开发效果要优于常规储层水驱开发效果[8]。结合这2种提高采收率方式,Yuan B等[9]针对天然裂缝发育致密储层,将压裂和水驱结合,通过数值模拟开展多级压裂水驱实验研究,并得出这种方法能够提高产量和最终采收率。然而,在该研究过程中忽略了天然裂缝的分布对水驱油的影响,天然裂缝的分布也会在压裂后影响注入水的波及范围,认识不清则同样会导致水窜,浪费注入水。

      考虑储层天然裂缝发育、储层和流体的可压缩性以及流体重力对渗流的影响,基于非活塞式水驱油理论,建立了油水两相渗流数学模型,通过有限元方法定量分析天然裂缝分布对多级压裂储层水驱油开发效果的影响,揭示油水两相流体在多尺度孔隙结构中的渗流规律和剩余油分布规律,为后期寻找剩余油、堵水调剖及制定采油开发方案奠定基础。

    • 储层单井模拟区域为480 m×400 m×70 m,根据其对称性,选用水平井筒以上储层区域进行研究,该井段采用3级分段压裂技术,水力裂缝在高应力和高逼近角条件下于交点处直接穿过天然裂缝,继续沿原方向延伸[10-11],水力裂缝与天然裂缝夹角为a,在近井筒附近形成典型的基质-裂缝系统,物理模型如图1(a)所示,其中$ {\Gamma _{\text{1}}} $$ {\Gamma _{\text{2}}} $为储层边界;$ {\Gamma _{\text{3}}} $为人工裂缝;$ {\Gamma _{\text{4}}} $为天然裂缝。同时,压裂后采用人工注水保持压力,假设油水界面非均匀单向向前推进,即非活塞式水驱油,非活塞式水驱油单向流模型如图1(b)所示。

      图  1  双重孔隙介质物理模型及非活塞式水驱油单向流模型

      Figure 1.  Physical model of dual porous media and one-way flow model of non piston water flooding

    • 基于双重介质模型提出假设条件:(1)储层压裂后渗流通道主要为基质和裂缝,裂缝不连续;(2)储层均质各向同性、微可压缩且压缩系数不随时间变化;(3)考虑重力对致密油流动的影响;(4)油水两相之间相互独立,二者之间无相间传质现象发生;(5)忽略基岩与裂缝间的窜流;(6)忽略油井的不完善性;(7)储层中油水两相流动均符合达西渗流规律,且渗流过程等温;(8)储层流体仅包含油相和水相,流体微可压缩但压缩系数不随时间变化;(9)油水两相流体黏度与压力无关;(10)油水两相相对渗透率仅与饱和度有关。

    • 考虑油水两相流的基质系统连续性方程[12]

      $$ \dfrac{\partial }{\partial t}({\phi }_{\text{m}}{\rho }_{{\text{s}}_{{i}}}{S}_{{i}})+\nabla ·({\rho }_{{\text{s}}_{{i}}}{v}_{\text{m}{i}})=0 $$ (1)

      基质系统油水两相运动方程

      $$\left\{ \begin{array}{l} {\nu _{{\rm{mo}}}} = - \dfrac{{{k_{{{\rm{s}}_{\rm{o}}}}}}}{{{\mu _{\rm{o}}}}}\nabla ({p_{{{\rm{s}}_{\rm{o}}}}} - {\rho _{{{\rm{s}}_{\rm{o}}}}}g)\\ {\nu _{{\rm{mw}}}} = - \dfrac{{{k_{{{\rm{s}}_{\rm{w}}}}}}}{{{\mu _{\rm{w}}}}}\nabla( {p_{{{\rm{s}}_{\rm{w}}}}} - {\rho _{{{\rm{s}}_{\rm{w}}}}}g) \end{array} \right. $$ (2)

      基质系统固体状态方程

      $$ {\phi _{\rm{m}}} = {\phi _{{\rm{mo}}}} + {C_{\rm{m}}}({p_{\rm{m}}} - {p_{\rm{i}}})$$ (3)

      基质系统油水两相液体状态方程

      $$ \left\{ \begin{array}{l} {\rho _{{{\rm{s}}_{\rm{o}}}}} = {\rho _{\rm{o}}}\left[ {1 + {C_{\rm{o}}}\left( {{p_{\rm{m}}} - {p_{{{\rm{s}}_{\rm{o}}}}}} \right)} \right]\\ {\rho _{{{\rm{s}}_{\rm{w}}}}} = {\rho _{\rm{w}}}\left[ {1{\rm{ + }}{C_{\rm{w}}}\left( {{p_{\rm{m}}} - {p_{{{\rm{s}}_{\rm{w}}}}}} \right)} \right] \end{array} \right.$$ (4)

      基质系统油水两相相对渗透率

      $$ \left\{ \begin{array}{l} {k_{{\rm{r}}{{\rm{s}}_{\rm{o}}}}} = \dfrac{{{k_{{{\rm{s}}_{\rm{o}}}}}}}{{{k_{\rm{m}}}}} = {k_{{\rm{r}}{{\rm{s}}_{\rm{o}}}}}\left( {{S_{ \rm{o}}}} \right)\\ {k_{{\rm{r}}{{\rm{s}}_{\rm{w}}}}} = \dfrac{{{k_{{{\rm{s}}_{\rm{w}}}}}}}{{{k_{\rm{m}}}}} = {k_{\rm{r}}}_{{{\rm{s}}_{\rm{w}}}}\left( {{S_{ \rm{w}}}} \right) \end{array} \right.$$ (5)

      将式(2)~式(5)代入式(1)化简得出流出量和流入量公式。

      (1)流出量公式。

      $$ \begin{split} \dfrac{\partial }{{\partial t}}({\phi _{\rm{m}}}{\rho _{{{\rm{s}}_{{i}}}}}{S_{{i}}}){\rm{ = }}&\dfrac{\partial }{{\partial t}}\left( {{\phi _{\rm{m}}}{\rho _{{{\rm{s}}_{\rm{o}}}}}{S_{ \rm{o}}}} \right) + \dfrac{\partial }{{\partial t}}\left( {{\phi _{\rm{m}}}{\rho _{{{\rm{s}}_{\rm{w}}}}}{S_{ \rm{w}}}} \right){\rm{ = }}\\ & \dfrac{\partial }{{\partial t}}\left[ {{\rho _{\rm{o}}}{\phi _{{\rm{mo}}}}\left( {{S_{ \rm{o}}}{\rm{ + }}{C_{\rm{o}}}{S_{ \rm{o}}}} \right){\rm{ + }}{\rho _{\rm{o}}}{C_{\rm{m}}}{S_{ \rm{o}}}\left( {{p_{\rm{m}}} - {p_{\rm{i}}}} \right){\rm{ + }}} \right.\\ & \left. {{\rho _{\rm{o}}}{C_{\rm{m}}}{C_{\rm{o}}}{S_{ \rm{o}}}\left( {{p_{\rm{m}}} - {p_{\rm{i}}}} \right)\left( {{p_{\rm{m}}} - {p_{{{\rm{s}}_{\rm{o}}}}}} \right)} \right] + \\ & \dfrac{\partial }{{\partial t}}\left[ {{\rho _{\rm{w}}}{\phi _{{\rm{mo}}}}\left( {{S_{ \rm{w}}}{\rm{ + }}{C_{\rm{w}}}{S_{ \rm{w}}}} \right){\rm{ + }}{\rho _{\rm{w}}}{C_{\rm{m}}}{S_{ \rm{w}}}\left( {{p_{\rm{m}}} - {p_{\rm{i}}}} \right){\rm{ + }}} \right.\\ & \left. {{\rho _{\rm{w}}}{C_{\rm{m}}}{C_{\rm{w}}}{S_{ \rm{w}}}\left( {{p_{\rm{m}}} - {p_{\rm{i}}}} \right)\left( {{p_{\rm{m}}} - {p_{{{\rm{s}}_{\rm{w}}}}}} \right)} \right] \end{split}$$ (6)

      由于$ {C_{\rm{m}}}{C_{\rm{o}}}$${C_{\rm{m}}}{C_{\rm{w}}} $的乘积极小忽略不计,式(6)化简为

      $$ \begin{split}& \dfrac{\partial }{{\partial t}}({\phi _{\rm{m}}}{\rho _{{{\rm{s}}_{{i}}}}}{S_{{i}}}){\rm{ = }}\dfrac{\partial }{{\partial t}}[ {{\rho _{\rm{o}}}{\phi _{{\rm{mo}}}}\left( {{S_{ \rm{o}}}{\rm{ + }}{C_{\rm{o}}}{S_{ \rm{o}}}} \right){\rm{ + }}{\rho _{\rm{o}}}{C_{\rm{m}}}{S_{ \rm{o}}}({p_{\rm{m}}} - {p_{\rm{i}}}}) ]{\rm{ + }}\\&\qquad \dfrac{\partial }{{\partial t}}\left[ {{\rho _{\rm{w}}}{\phi _{{\rm{mo}}}}\left( {{S_{ \rm{w}}}{\rm{ + }}{C_{\rm{w}}}{S_{ \rm{w}}}} \right){\rm{ + }}{\rho _{\rm{w}}}{C_{\rm{m}}}{S_{ \rm{w}}}({p_{\rm{m}}} - {p_{\rm{i}}}}) \right] \end{split}$$ (7)

      储层中任意一点油水两相饱和度关系为${S_{ \rm{w}}} + {S_{ \rm{o}}} = 1 $,因此式(7)简化为

      $$ \begin{split}& \dfrac{\partial }{{\partial t}}({\phi _{\rm{m}}}{\rho _{{{\rm{s}}_{{i}}}}}{S_{{i}}}){\rm{ = }}\dfrac{\partial }{{\partial t}}\left[ {{\rho _{\rm{o}}}{\phi _{{\rm{mo}}}}\left( {{S_{ \rm{o}}}{\rm{ + }}{C_{\rm{o}}}{S_{ \rm{o}}}} \right){\rm{ + }}{\rho _{\rm{o}}}{C_{\rm{m}}}{S_{ \rm{o}}}\left( {{p_{\rm{m}}} - {p_{\rm{i}}}} \right)} \right]{\rm{ + }}\\ &\qquad \dfrac{\partial }{{\partial t}}\left[ { - {\rho _{\rm{w}}}{\phi _{{\rm{mo}}}}{S_{ \rm{o}}}\left( {{\rm{1 + }}{C_{\rm{w}}}} \right){\rm{ + }}{\rho _{\rm{w}}}{C_{\rm{m}}}\left( {1 - {S_{ \rm{o}}}} \right)\left( {{p_{\rm{m}}} - {p_{\rm{i}}}} \right)} \right]\\[-15pt] \end{split}$$ (8)

      (2)流入量公式。

      $$ \begin{split}& \nabla · ({\rho _{{{\rm{s}}_{{i}}}}}{\nu _{{{{\rm{m}}}{i}}}}) = \dfrac{{ - {k_{{{\rm{s}}_{\rm{o}}}}}}}{{{\mu _{\rm{o}}}}}\nabla ·\left( {{p_{{{\rm{s}}_{\rm{o}}}}} - {\rho _{{{\rm{s}}_{\rm{o}}}}}g} \right) + \dfrac{{ - {k_{{{\rm{s}}_{\rm{w}}}}}}}{{{\mu _{\rm{w}}}}}\nabla· \left( {{p_{{{\rm{s}}_{\rm{w}}}}} - {\rho _{{{\rm{s}}_{\rm{w}}}}}g} \right)\\ &\qquad {\rm{ = }}\dfrac{{ - {k_{{\rm{r}}{{\rm{s}}_{\rm{o}}}}}}}{{{\mu _{\rm{o}}}}}{k_{\rm{m}}}\nabla \left\{ {{p_{{{\rm{s}}_{\rm{o}}}}} - {\rho _{\rm{o}}}\left[ {1 + {C_{\rm{o}}}\left( {{p_{\rm{m}}} - {p_{{{\rm{s}}_{\rm{o}}}}}} \right)} \right]g} \right\}+\\ &\qquad \dfrac{{ - {k_{{\rm{r}}{{\rm{s}}_{\rm{w}}}}}}}{{{\mu _{\rm{w}}}}}{k_{\rm{m}}}\nabla \left\{ {{p_{{{\rm{s}}_{\rm{w}}}}} - {\rho _{\rm{w}}}\left[ {1 + {C_{\rm{w}}}\left( {{p_{\rm{m}}} - {p_{{{\rm{s}}_{\rm{w}}}}}} \right)} \right]g} \right\}\\[-17pt] \end{split}$$ (9)

      流入量与流出量相等,式(1)最终简化为

      $$ \begin{split}& \dfrac{\partial }{{\partial t}}\left[ {{\rho _{\rm{o}}}{\phi _{{\rm{mo}}}}{S_{ \rm{o}}}\left( {{\rm{1 + }}{C_{\rm{o}}}} \right){\rm{ + }}{\rho _{\rm{o}}}{C_{\rm{m}}}{S_{ \rm{o}}}({p_{\rm{m}}} - {p_{\rm{i}})}} \right]{\rm{ + }}\\ &\qquad \dfrac{\partial }{{\partial t}}\left[ { - {\rho _{\rm{w}}}{\phi _{{\rm{mo}}}}{S_{ \rm{o}}}\left( {{\rm{1 + }}{C_{\rm{w}}}} \right){\rm{ + }}{\rho _{\rm{w}}}{C_{\rm{m}}}\left( {1 - {S_{ \rm{o}}}} \right)({p_{\rm{m}}} - {p_{\rm{i}})}} \right] = \\ &\qquad \dfrac{{{k_{{\rm{r}}{{\rm{s}}_{\rm{o}}}}}}}{{{\mu _{\rm{o}}}}}{k_{\rm{m}}}\nabla \left[ {{p_{\rm{o}}} - {C_{\rm{o}}}{\rho _{\rm{o}}}g\left( {{p_{\rm{m}}} - {p_{\rm{o}}}} \right)} \right] + \\ &\qquad \dfrac{{{k_{{\rm{r}}{{\rm{s}}_{\rm{w}}}}}}}{{{\mu _{\rm{w}}}}}{k_{\rm{m}}}\nabla \left[ {{p_{\rm{w}}} - {C_{\rm{w}}}{\rho _{\rm{w}}}g\left( {{p_{\rm{m}}} - {p_{\rm{w}}}} \right)} \right]\\[-16pt] \end{split}$$ (10)
    • 考虑油水两相流裂缝系统连续性方程[15]

      $$ {d_{\rm{f}}}\dfrac{\partial }{{\partial t}}({\phi _{\rm{f}}}{\rho _{{{\rm{s}}_{{i}}}}}{S_{{i}}}) + {d_{\rm{f}}}\nabla · ({\rho _{{{\rm{s}}_{{i}}}}}{\nu _{{{{\rm{f}}i}}}}) = 0$$ (11)

      裂缝系统油水两相运动方程

      $$ \left\{ \begin{array}{l} {\nu _{{\rm{fo}}}} = - \dfrac{{{k_{{{\rm{s}}_{\rm{o}}}}}}}{{{\mu _{\rm{o}}}}}\nabla ({p_{{{\rm{s}}_{\rm{o}}}}} - {\rho _{{{\rm{s}}_{\rm{o}}}}}g)\\ {\nu _{{\rm{fw}}}} = - \dfrac{{{k_{{{\rm{s}}_{\rm{w}}}}}}}{{{\mu _{\rm{w}}}}}\nabla ({p_{{{\rm{s}}_{\rm{w}}}}} - {\rho _{{{\rm{s}}_{\rm{w}}}}}g) \end{array} \right.$$ (12)

      裂缝系统固体状态方程

      $$ {\phi _{\rm{f}}} = {\phi _{{\rm{fo}}}} + {C_{\rm{f}}}({p_{\rm{f}}} - {p_{\rm{i}}})$$ (13)

      裂缝系统油水两相液体状态方程

      $$ \left\{ \begin{array}{l} {\rho _{{{\rm{s}}_{\rm{o}}}}} = {\rho _{\rm{o}}}\left[ {1 + {C_{\rm{o}}}\left( {{p_{\rm{f}}} - {p_{{{\rm{s}}_{\rm{o}}}}}} \right)} \right]\\ {\rho _{{{\rm{s}}_{\rm{w}}}}} = {\rho _{\rm{w}}}\left[ {1{\rm{ + }}{C_{\rm{w}}}\left( {{p_{\rm{f}}} - {p_{{{\rm{s}}_{\rm{w}}}}}} \right)} \right] \end{array} \right.$$ (14)

      裂缝系统油水两相相对渗透率

      $$ \left\{ \begin{array}{l} {k_{{\rm{r}}{{\rm{s}}_{\rm{o}}}}} = \dfrac{{{k_{{{\rm{s}}_{\rm{o}}}}}}}{{{k_{\rm{f}}}}} = {k_{{\rm{r}}{{\rm{s}}_{\rm{o}}}}}\left( {{S_{ \rm{o}}}} \right)\\ {k_{{\rm{r}}{{\rm{s}}_{\rm{w}}}}} = \dfrac{{{k_{{{\rm{s}}_{\rm{w}}}}}}}{{{k_{\rm{f}}}}}{\rm{ = }}{k_{{\rm{r}}{{\rm{s}}_{\rm{w}}}}}\left( {{S_{ \rm{w}}}} \right) \end{array} \right.$$ (15)

      将式(12)~式(15)代入式(11)化简得

      $$ \begin{split}& {d_{\rm{f}}}\dfrac{\partial }{{\partial t}}\left[ {{\rho _{\rm{o}}}{\phi _{{\rm{fo}}}}\left( {{S_{ \rm{o}}}{\rm{ + }}{C_{\rm{o}}}{S_{ \rm{o}}}} \right){\rm{ + }}{\rho _{\rm{o}}}{C_{\rm{m}}}{S_{ \rm{o}}}({p_{\rm{f}}} - {p_{\rm{i}}}}) \right]{\rm{ + }}\\&\qquad {d_{\rm{f}}}\dfrac{\partial }{{\partial t}}\left[ { - {\rho _{\rm{w}}}{\phi _{{\rm{fo}}}}{S_{ \rm{o}}}\left( {{\rm{1 + }}{C_{\rm{w}}}} \right){\rm{ + }}{\rho _{\rm{w}}}{C_{\rm{m}}}\left( {1 - {S_{ \rm{o}}}} \right){p_{\rm{f}}} - {p_{\rm{i}}}} \right]{\rm{ = }}\\&\qquad {\rm{ }}{d_{\rm{f}}}\dfrac{{{k_{{\rm{r}}{{\rm{s}}_{\rm{o}}}}}}}{{{\mu _{\rm{o}}}}}{k_{\rm{f}}}\nabla \left[ {{p_{\rm{o}}} - {C_{\rm{o}}}{\rho _{\rm{o}}}g\left( {{p_{\rm{f}}} - {p_{\rm{o}}}} \right)} \right] + {\rm{ }}\\&\qquad {d_{\rm{f}}}\dfrac{{{k_{{\rm{r}}{{\rm{s}}_{\rm{w}}}}}}}{{{\mu _{\rm{w}}}}}{k_{\rm{f}}}\nabla \left[ {{p_{\rm{w}}} - {C_{\rm{w}}}{\rho _{\rm{w}}}g\left( {{p_{\rm{f}}} - {p_{\rm{w}}}} \right)} \right]\\[-17pt] \end{split}$$ (16)
    • 任意时刻通过两相区内任一截面水、油流量为

      $$ {q_{\rm{w}}} = - \dfrac{{k{k_{{\rm{rs_{{\rm{w}}}}}}}}}{{{\mu _{\rm{w}}}}}A\nabla \left( {{p_{\rm{w}}} - {\rho _{\rm{w}}}g} \right)$$ (17)
      $$ {q_{\rm{o}}} = - \dfrac{{k{k_{{\rm{rs_{{\rm{o}}}}}}}}}{{{\mu _{\rm{o}}}}}A\nabla \left( {{p_{\rm{o}}} - {\rho _{\rm{o}}}g} \right)$$ (18)

      任意时刻该断面上的含水率、含油率为

      $$ {f_{\rm{w}}} = \dfrac{{{q_{\rm{w}}}}}{{{q_{\rm{w}}} + {q_{\rm{o}}}}}$$ (19)
      $$ {f_{\rm{o}}} = 1 - {f_{\rm{w}}}$$ (20)
    • (1)外边界条件。

      $$ {\left. p \right|_{{\Gamma _1}}} = {p_{\rm{e}}}$$ (21)

      (2)内边界条件。储层内边界条件为

      $$ {\left. {\dfrac{{\partial p\left( {x,y,t} \right)}}{{\partial n}}} \right|_{{\Gamma _2}}} = 0, \left( {x,y} \right) \in {\Gamma _2}$$ (22)
      $$ {\left. {p\left( {x,y,t} \right)} \right|_{{\Gamma _3}}} = {p_{\rm{w}}},\left( {x,y} \right) \in {\Gamma _3}$$ (23)

      裂缝内边界条件为

      $$ {p_{\rm{m}}}\left( {x,y,t} \right) = {\left. {{p_{\rm{f}}}\left( {x,y,t} \right)} \right|_{{\Gamma _4}}},\left( {x,y} \right) \in {\Gamma _4}$$ (24)
      $$ {\left. {{k_{\rm{m}}}\dfrac{{\partial {p_{\rm{m}}}}}{{\partial n}}} \right|_{{\Gamma _4}}} = {\left. {{k_{\rm{f}}}\dfrac{{\partial {p_{\rm{f}}}}}{{\partial n}}} \right|_{{\Gamma _4}}},\left( {x,y} \right) \in {\Gamma _4}$$ (25)

      油水边界条件为

      $$ {\left. {{p_{\rm{o}}}} \right|_{\rm{f}}} = {\left. {{p_{\rm{w}}}} \right|_{\rm{f}}}$$ (26)
      $$ {\left. {{u_{\rm{o}}}} \right|_{\overrightarrow n }} = {\left. {{u_{\rm{w}}}} \right|_{\overrightarrow n }}$$ (27)

      (3)初始条件。

      $$ {p_{\rm{m}}}\left( {x,y,t{\rm{ = }}0} \right){\rm{ = }}{p_{\rm{f}}}\left( {x,y,t{\rm{ = }}0} \right) = {p_{\rm{o}}}\left( {x,y,t{\rm{ = }}0} \right)$$ (28)
    • 井深4 000 m,井眼半径0.01 m,储层初始压力42.32 MPa,井底压力12 MPa,油的饱和压力11.34 MPa,油藏长度400 m,油藏宽度240 m,油藏厚度70 m,水相黏度1 mPa·s,油相黏度2.25 mPa·s,水相密度1 000 kg/m3,油相密度775 kg/m3,基质初始孔隙度0.0164,人工裂缝初始孔隙度0.45,天然裂缝初始孔隙度0.45,基质渗透率0.004 μm2,人工裂缝渗透率1.5 μm2,天然裂缝渗透率0.04 μm2,人工裂缝宽度5 mm,天然裂缝宽度2 mm,基质压缩系数3×10−4 MPa−1,人工裂缝压缩系数4×10−4 MPa−1,天然裂缝压缩系数2.5×10−4 MPa−1,油相压缩系数7×10−3 MPa−1,水相压缩系数3.7×10−3 MPa−1,注水速度30 m3/d,油相相对渗透率为${k_{{{{\rm{rs}}_{{\rm{o}}}}}}} = 35.408 - 240.688{S_{\text{w}}}+ {\text{ 618}}{\text{.365}}S_{\text{w}}^2 - 707.125S_{\text{w}}^3 + 302.464S_{\text{w}}^4$,水相相对渗透率为${k_{{{{\rm{rs}}_{{\rm{w}}}}}}} = 0.457 - 2.737{S_{\text{w}}} + 5.001S_{\text{w}}^2 - 2.508S_{\text{w}}^3$

    • 选取瞬态油相流量和储层含油率为基本变量,采用有限元数值解法对水驱油数学方程进行离散求解。首先对求解域进行网格划分,选用三角形网格,并在裂缝区域进行网格加密,通过自由剖分方式离散化求解域,网格剖分结果如图2。将求解参数赋值于基质系统和裂缝系统,基于水驱油数学方程,采用向后差分方法,最终获得水驱开发条件下油田的瞬态产量和储层含油率。

      图  2  网格剖分结果

      Figure 2.  Mesh generation results

    • (1)累积产量。图3为不同的人工裂缝与天然裂缝夹角下储层累积产量随时间变化,可以看出,采用水驱开发压裂后的储层能大幅提高单井产量。

      图  3  不同裂缝夹角下储层累积产量随时间变化

      Figure 3.  Variation of reservoir cumulative production with time under different fracture angles

      水驱初始阶段,基质系统中首先进行水驱油,从供给边界到原始油水界面只有水在运动,水驱前缘在基质系统中匀速向前推进,而在近井筒附近人工裂缝中的原油依托高导流裂缝进入水平井筒,水驱前缘到水平井筒间为纯油流,进而油井累积产量迅速增加,累积产油量从小到大依次为人工裂缝与天然裂缝夹角为45°、75°、60°、15°、90°和30°时的产量,且油井产量递减速度较小。随着水驱开采时间延长,水驱前缘进入裂缝区域,由于裂缝渗流阻力较小且水的黏度较小,注入水沿裂缝突进,同时,基质系统中的油在重力和供给边界压力及水驱动力的作用下向裂缝系统中运移,此时水驱前缘到水平井筒之间仍为纯油流,但原始油水界面到水驱前缘为油水两相流。然而,由于基质的渗透率极低,造成压力无法向储层基质含油区域较远边界传播[13],同时,油水两相的存在也增大了油相的渗流阻力,使得基质中的油向裂缝运移的速度变慢,无法快速运移至水平井筒,油井产量递减速度较大,致使油井累积产量增加速度减慢,油井产量递减速度从小到大依次为人工裂缝与天然裂缝夹角为90°、60°、45°、30°、75°和15°时的产量,所以人工裂缝与天然裂缝夹角为90°时油井产量最高。水驱油开发后期,油井仍未见水,受水驱前缘动力和储层含油量影响,油井以较低产量进行生产,日产量在5m3/d波动,累积产油量从小到大依次为人工裂缝与天然裂缝夹角为15°、30°、75°、60°、45°和90°时的产量,人工裂缝与天然裂缝夹角相对越大,油井累积产量相对越高。

      (2)含油率。图4为压裂储层水驱开发第2 000 d时不同裂缝分布下储层剩余油分布,总体而言,水驱油在基质和裂缝中同时发生,水驱油效率较低且水驱前缘向前推进不均匀。基质系统渗透率相对裂缝系统较低,注入水难以大范围波及基质系统,注入水则发生绕流并沿着阻力较小的裂缝系统向前突进,水窜较快[14-15]。因此,基质中的水驱油效率极低,裂缝中的水驱油效率相对基质中的水驱油效率较高。与此同时,不同的人工裂缝和天然裂缝夹角显著影响水驱油开发效果。当人工裂缝与天然裂缝夹角小于等于30°时,基质系统和裂缝系统的驱油效果相当;当人工裂缝与天然裂缝夹角大于30°时,裂缝系统的驱油效果明显优于基质系统;当人工裂缝与天然裂缝夹角为90°时,裂缝系统和基质系统的驱油效果显著优于其他角度的裂缝分布。当人工裂缝与天然裂缝夹角为90°时,人工裂缝与主流线平行,增大了储层渗透率,天然裂缝与主流线方向垂直,阻挡了注入水沿着主流线方向的推进,注入水发生部分绕流,增大了注入水在基质系统的波及范围,使得基质中无法运移油活动起来。同时,裂缝间干扰对水驱油效果影响显著,不同的人工裂缝与天然裂缝夹角条件下裂缝间的干扰程度不同,人工裂缝与天然裂缝夹角越大,裂缝间的干扰越明显,注入水更能够波及裂缝间难以驱替的基质区域,水驱开发后剩余油较少。因此,当人工裂缝与天然裂缝夹角为90°时,注入水波及范围最大,剩余油最少。

      图  4  第2 000 d储层含油率分布

      Figure 4.  Distribution of oil-bearing rate of the 2 000 d reservoir

    • (1)基于非活塞式水驱油机理,考虑流体的重力和裂缝对水驱油效果的影响,建立了压裂后水驱开发致密油藏的油水两相渗流模型,通过累积产量和含油率定量分析人工裂缝与天然裂缝夹角对压裂储层水驱开发效果。

      (2)人工裂缝与天然裂缝夹角显著影响压裂储层单井产量,人工裂缝与天然裂缝夹角越大,油井产量相对越高。水驱初始阶段,供给边界到原始油水界面为纯水流,油井累积产量迅速增加,累积产油量从小到大依次为人工裂缝与天然裂缝夹角为45°、75°、60°、15°、90°和30°时的产量,且油井产量递减速度较小;随着水驱开采时间的延长,水驱前缘进入裂缝区域,注入水沿着裂缝突进,原始油水界面到水驱前缘为油水两相流,油井累积产量增加速度减小,而油井产量递减速度增大,从小到大依次为人工裂缝与天然裂缝夹角为90°、60°、45°、30°、75°和15°时的产量;水驱油开发后期,从油水前缘到生产井为纯油流,油井以较低产量进行生产,累积产油量从小到大依次为人工裂缝与天然裂缝夹角为15°、30°、75°、60°、45°和90°时的产量。

      (3)人工裂缝与天然裂缝夹角显著影响压裂储层剩余油分布,夹角增大,扩大了注入水波及基质系统的范围,驱油效果显著,但水驱油效率整体较低。当人工裂缝与天然裂缝夹角小于等于30°时,基质系统和裂缝系统的驱油效果相当;当人工裂缝与天然裂缝夹角大于30°时,裂缝系统的驱油效果明显优于基质系统;当人工裂缝与天然裂缝夹角为90°时,裂缝系统和基质系统的驱油效果显著优于其他角度的裂缝分布,注入水波及范围最大,剩余油最少。

      符号说明:

      A为液体通过的截面面积,m2$ {C_{\text{f}}} $$ {C_{\text{m}}} $分别为裂缝系统、基质系统孔隙压缩系数,MPa−1$ {C_{\text{o}}} $$ {C_{\text{w}}} $分别为油相、水相压缩系数,MPa−1${d_{\text{f}}}$为裂缝宽度,mm;${f_{\rm{o}}}$${f_{\rm{w}}}$分别为含油率和含水率,无量纲;$ {k_{\text{f}}} $$ {k_{\text{m}}} $分别为裂缝系统和基质系统渗透率,μm2$ {k_{{{\text{s}}_{\text{o}}}}} $$ {k_{{{\text{s}}_{\text{w}}}}} $分别为油相和水相渗透率,μm2$ {p_{\text{f}}} $为裂缝系统压力,MPa;$ {p_{\text{i}}} $为原始储层压力,MPa;$ {p_{\text{m}}} $为基质系统压力,MPa; po为储层油相压力,MPa; pw为储层水相压力,MPa; $ {p_{{{\text{s}}_{\text{o}}}}} $为油相饱和度为So时储层油相压力,MPa; $ {p_{{{\text{s}}_{\text{w}}}}} $为水相饱和度为Sw时储层水相压力,MPa;$ {q_{\text{o}}} $$ {q_{\text{w}}} $分别为油藏条件下油相、水相流量,m3/d;${S_{{i}}}$为某一液相(油相或水相)饱和度,%;${S_{\text{o}}}$为油相饱和度,%;${S_{\text{w}}}$为水相饱和度,%;$ {v_{\text{m}}} $为基质系统流体渗流速度,cm/s;$ {v_{{\text{mo}}}} $$ {v_{{\text{fo}}}} $分别为基质系统、裂缝系统油相渗流速度,cm/s;$ {v_{{\text{mw}}}} $$ {v_{{\text{fw}}}} $分别为基质系统、裂缝系统水相渗流速度,cm/s;${v_{{{{\rm{f}}i}}}}$为裂缝系统某一液相(油相或水相)渗流速度,cm/s;$ {\phi _{\text{f}}} $为裂缝孔隙度,%;$ {\phi _{{\text{fo}}}} $为裂缝系统初始孔隙度,%;$ {\phi _{\text{m}}} $为基质孔隙度,%;$ {\phi _{{\text{mo}}}} $为基质系统初始孔隙度,%;$ {\rho _{{{\text{s}}_{\text{o}}}}} $为储层条件下油相密度,kg/m3$ {\rho _{{{\text{s}}_{\text{w}}}}} $为储层条件下水相密度,kg/m3$ {\rho _{\text{o}}} $$ {\rho _{\text{w}}} $分别为原始储层油相、水相密度,kg/m3${\rho _{{{\text{s}}_{{i}}}}}$为储层条件下某一液相(油相或水相)密度,kg/m3$ {\mu _{\text{o}}} $$ {\mu _{\text{w}}} $分别为储层条件下油相、水相黏度,Pa·s;$ \nabla $为梯度算符,无量纲。

参考文献 (15)

目录

    /

    返回文章
    返回