• 全国中文核心期刊
  • 中国科技核心期刊
  • 中国核心学术期刊(RCCSE)
  • PA收录期刊
  • CA收录期刊

基于分形理论的页岩气分支水平井产能数学模型

李忠厚 吴小斌 DUZhongwei 任茜莹

李忠厚, 吴小斌, DUZhongwei, 任茜莹. 基于分形理论的页岩气分支水平井产能数学模型[J]. 石油钻采工艺, 2017, 39(1): 7-13. doi: 10.13639/j.odpt.2017.01.002
引用本文: 李忠厚, 吴小斌, DUZhongwei, 任茜莹. 基于分形理论的页岩气分支水平井产能数学模型[J]. 石油钻采工艺, 2017, 39(1): 7-13. doi: 10.13639/j.odpt.2017.01.002
LI Zhonghou, WU Xiaobin, Zhongwei DU, REN Qianying. A mathematical model for productivity calculation of shale-gas multilateral horizontal wells based on fractal theory[J]. Oil Drilling & Production Technology, 2017, 39(1): 7-13. doi: 10.13639/j.odpt.2017.01.002
Citation: LI Zhonghou, WU Xiaobin, Zhongwei DU, REN Qianying. A mathematical model for productivity calculation of shale-gas multilateral horizontal wells based on fractal theory[J]. Oil Drilling & Production Technology, 2017, 39(1): 7-13. doi: 10.13639/j.odpt.2017.01.002

基于分形理论的页岩气分支水平井产能数学模型

doi: 10.13639/j.odpt.2017.01.002
基金项目: 国家科技重大专项“井震结合油藏精细结构表征技术”(编号:2011ZX05010-001);陕西省高水平大学建设专项资金资助项目(编号:2013SXTS03);延安大学专项科研基金项目(编号:YDK2012-1)
详细信息
    作者简介:

    李忠厚(1982-),2008年毕业于中国石油大学(华东)油气田开发工程专业,现为延安大学讲师,从事物理法采油,非常规油气田开发方面的研究工作。通讯地址:(716000)中国陕西延安市圣地路580号延安大学石油工程与环境工程学院。电话:18329909987。E-mail:wo119944595er@163.com

  • 中图分类号: TE328

A mathematical model for productivity calculation of shale-gas multilateral horizontal wells based on fractal theory

  • 摘要: 利用分形几何法,研究了页岩储集层内水平井分支数目和长度的分形维度,以及水平井各分支井段压降扰动的传播规律,得到压降漏斗边界随时间变化的关系,考虑了页岩储集层解吸吸附、扩散系数和启动压力梯度的影响,建立了页岩气分支水平井产能数学模型。采用拉普拉斯变换法,求解了水平井外边界定压条件下的不稳定渗流的产能方程。结合四川某海相页岩气藏储集层参数,计算分析了页岩气分支水平井产能特征及影响因素。研究表明:在页岩气生产过程中,同一分形维度下,水平井分支数目越大,气井的产量增大越快;分支井的井段长度越大,对主干井的流量影响越大。吸附页岩气的解吸使得压降漏斗的边界传播速度减慢,地层压力下降缓慢;产气量越大,井底压力随时间下降越快,但是下降的趋势随时间逐渐减缓。
  • 图  1  页岩气分支水平井平面径向流示意图

    Figure  1.  Schematic areal radial flow of shale-gas multilateral horizontal well

    图  2  页岩气水平井分支平面径向流示意图

    Figure  2.  Schematic areal radial flow of branch holes of shale-gas horizontal well

    图  3  不同分支数下压降漏斗边界与分支井半径比值随时间变化关系

    Figure  3.  Relationship of the ratio of pressure drop funeral boundary to branch hole radius vs. time for different branch hole amounts

    图  4  不同启动压力梯度下压降漏斗边界与主干井半径比值随时间变化关系

    Figure  4.  Relationship of the ratio of pressure drop funeral boundary to main hole radius vs. time for different startup pressure gradients

    图  5  不同分支长度对气井生产的影响

    Figure  5.  Effect of branch hole length on gas well production

    图  6  吸附解吸气体量对气井生产的影响

    Figure  6.  Effect of adsorption/desorption gas volume on gas well production

    图  7  基于分形模型的气井模拟值与实际产量对比

    Figure  7.  Comparison between actual production and simulated value of gas wells based on fractal model

  • [1] 郭小哲, 周长沙.基于滑脱的页岩气藏压裂水平井渗流模型及产能预测[J].石油钻采工艺, 2015, 37(3): 61-65. http://www.cqvip.com/main/zcps.aspx?c=1&id=664873205

    GUO Xiaozhe, ZHOU Changsha. Seepage model and productivity forecast based on slippage of fractured horizontal wells in shale gas pool[J]. Oil Drilling & Production Technology, 2015, 37(3): 61-65. http://www.cqvip.com/main/zcps.aspx?c=1&id=664873205
    [2] 邹才能, 董大忠, 王玉满, 李新景, 黄金亮, 王淑芳, 管全中, 张晨晨, 王红岩, 刘洪林, 拜文华, 梁峰, 吝文, 赵群, 刘德勋, 杨智, 梁萍萍, 孙莎莎, 邱振.中国页岩气特征、挑战及前景(二)[J].石油勘探与开发, 2016, 43(2): 166-178.

    ZOU Caineng, DONG Dazhong, WANG Yuman, LI Xinjing, HUANG Jinliang, WANG Shufang, GUAN Quanzhong, ZHANG Chenchen, WANG Hongyan, LIU Honglin, BAI Wenhua, LIANG Feng, LIN Wen, ZHAO Qun, LIU Dexun, YANG Zhi, LIANG Pingping, SUN Shasha, QIU Zhen. Shale gas in China: Characteristics, challenges and prospects(Ⅱ)[J]. Petroleum Exploration and Development, 2016, 43(2): 166-178.
    [3] 宁正福, 王波, 杨峰, 曾彦, 陈进娥, 张廉.页岩储集层微观渗流的微尺度效应[J].石油勘探与开发, 2014, 41(4): 445-452.

    NING Zhengfu, WANG Bo, YANG Feng, ZENG Yan, CHEN Jine, ZHANG Lian. Microscale effect of microvadose in shale reservoirs[J]. Petroleum Exploration and Development, 2014, 41(4): 445-452.
    [4] 曾传云, 谭勇志, 徐红国, 芦建江, 邹敏, 张斌成.多分支水平井钻井技术在玉平8井的应用[J].石油钻采工艺, 2013, 35(4): 22-25.

    ZENG Chuanyun, TAN Yongzhi, XU Hongguo, LU Jianjiang, ZOU Min, ZHANG Bincheng. Application of multi-lateral horizontal well drilling technique on Well Yuping-8[J]. Oil Drilling & Production Technology, 2013, 35(4): 22-25.
    [5] 陈作, 曾义金.深层页岩气分段压裂技术现状及发展建议[J].石油钻探技术, 2016, 44(1): 6-11.

    CHEN Zuo, ZENG Yijin. Present situations and prospects of multi-stage fracturing technology for deep shale gas development[J]. Petroleum Drilling Techniques, 2016, 44(1): 6-11.
    [6] 郭小哲, 周长沙.考虑应力敏感的页岩气藏压裂水平井渗流模型建立与分析[J].长江大学学报(自科版), 2015, 12(5): 44-50. http://www.cqvip.com/QK/88443X/201502/664074188.html

    GUO Xiaozhe, ZHOU Changsha. Establishment of a seepage model considering the stress sensitivity in fractured horizontal wells in shale gas reservoirs and itsanalysis[J]. Journal of Yangtze University(Natural Science Edition), 2015, 12(5): 44-50. http://www.cqvip.com/QK/88443X/201502/664074188.html
    [7] 朱维耀, 亓倩, 马千, 邓佳, 岳明, 刘玉章.页岩气不稳定渗流压力传播规律和数学模型[J].石油勘探与开发, 2016, 43(2): 261-267. http://d.wanfangdata.com.cn/Periodical/syktykf201602012

    ZHU Weiyao, QI Qian, MA Qian, DENG Jia, YUE Ming, LIU Yuzhang. Unstable seepage modeling and pressure propagation of shale gas reservoirs[J]. Petroleum Exploration and Development, 2016, 43(2): 261-267. http://d.wanfangdata.com.cn/Periodical/syktykf201602012
    [8] 张建军, 赵宇渊, 张南军, 王根常, 张一佳.不同类型油藏水平井开发适应性分析[J].石油钻采工艺, 2009, 31(6): 9-12.

    ZHANG Jianjun, ZHAO Yuyuan, ZHANG Nanjun, WANG Genchang, ZHANG Yijia. Adaptability analysis of horizontal well technology to oil reservoirs[J]. Oil Drilling & Production Technology, 2009, 31(6): 9-12.
    [9] PASCAL H. Non-steady flow through porous media in the presence of a threshold gradient[J]. Acta Mechanica, 1981, 39(3/4): 207-224.
    [10] 姜瑞忠, 滕文超, 徐建春.压裂改造复合页岩气藏不稳定压力与产量分析方法[J].天然气工业, 2015, 35(9): 42-47.

    JIANG Ruizhong, TENG Wenchao, XU Jianchun. Transient pressure and production analysis methods for composite shale gas reservoirs stimulated by fracturing[J]. Natural Gas Industry, 2015, 35(9): 42-47.
    [11] 谭晓华, 李晓平, 刘从领, 余燕, 宋昭杰, 袁淋.确定气水同产水平井流入动态关系的新方法[J].石油钻采工艺, 2014, 36(3): 59-64.

    TAN Xiaohua, LI Xiaoping, LIU Congling, YU Yan, SONG Zhaojie, YUAN Lin. A new method for ascertaining the influx dynamic relation in watergas producing horizontal wells[J]. Oil Drilling & Production Technology, 2014, 36(3): 59-64.
    [12] 刘晓华, 邹春梅, 姜艳东, 翟振宇.页岩气水平井动态评价方法[J].石油钻采工艺, 2013, 35(3): 55-58.

    LIU Xiaohua, ZOU Chunmei, JIANG Yandong, ZHAI Zhenyu. Horizontal well dynamic evaluation method in shale gas reservoirs[J]. Oil Drilling & Production Technology, 2013, 35(3): 55-58.
    [13] 李彦尊, 李相方, 徐敏, 武文涛, 张晓舟.页岩气藏产量递减预测模型研究及应用[J].石油钻采工艺, 2015, 37(6): 74-77. http://d.wanfangdata.com.cn/Periodical/dkyqt201504016

    LI Yanzun, LI Xiangfang, XU Min, WU Wentao, ZHANG Xiaozhou. Production decline prediction modelling and application in shale gas reservoir[J]. Oil Drilling & Production Technology, 2015, 37(6): 74-77. http://d.wanfangdata.com.cn/Periodical/dkyqt201504016
    [14] 范新欣, 王宝和, 于才渊, 徐朝辉, 齐永青, 陈登峰.基于分形理论的油页岩有效扩散系数研究[J].化学工程, 2010, 38(10): 238-242.

    FAN Xinxin, WANG Baolai, YU Caiyuan, XU Zhaohui, QI Yongqing, CHEN Dengfeng. Study on effective diffusion coefficient of oil shale based on fractal theory[J]. Chemical engineering, 2010, 38(10): 238-242.
    [15] 桑宇, 杨胜来, 郭小哲, 陈浩.一种已压裂页岩气水平井的产量预测新方法[J].西南石油大学学报(自然科学版), 2015, 37(3): 17-24. http://qikan.cqvip.com/Qikan/Article/Detail?id=664858056

    SANG Yu, YANG Shenglai, GUO Xiaozhe, CHEN Hao. A new productivity prediction method for fractured horizontal wells in shale gas reservoirs[J]. Journal of Southwest Petroleum University(Science & Technology Edition), 2015, 37(3): 17-24. http://qikan.cqvip.com/Qikan/Article/Detail?id=664858056
    [16] 朱维耀, 马千, 邓佳, 马旭东, 宋智勇, 岳明.纳微米级孔隙气体流动数学模型及应用[J].北京科技大学学报, 2014, 36(6): 709-715. http://www.cnki.com.cn/Article/CJFDTotal-BJKD201406001.htm

    ZHU Weiyao, MA Qian, DENG Jia, MA Dongxu, SONG Zhiyong, YUE Ming. Mathematical model and application of gas flow in nano-micron pores[J]. Journal of University of Science and Technology Beijing, 2014, 36(6): 709-715. http://www.cnki.com.cn/Article/CJFDTotal-BJKD201406001.htm
    [17] CLARKON C R. Production data analysis of unconventional gas wells: Review of theory and best practices[J]. International Journal of Coal Geology, 2013, 109: 101-146. http://www.sciencedirect.com/science/article/pii/S0166516213000050
    [18] CIPOLLA C L, LOLON E P, ERDLE J C, RUBIN B. Reservoir modeling in shale-gas reservoirs[J]. SPE Reservoir Evaluation & Engineering, 2010, 13(4): 638-653.
    [19] LANGMUIR I. The adsorption of gases on plane surfaces of glass, mica and platinum[J]. Journal of American Chemical Society, 1918, 40(9): 1403-1461. doi:  10.1021/ja02242a005
    [20] 郭为, 雄伟, 高树生, 胡志明.页岩气等温吸附/解吸特征[J].中南大学学报(自然科学版), 2013, 44(7): 2836-2840. http://d.wanfangdata.com.cn/Periodical/zngydxxb201307028

    GUO Wei, XIONG Wei, GAO Shusheng, HU Zhiming. Isothermal adsorption/desorption characteristics of shale gas[J]. Journal of Central South University(Science and Technology), 2013, 44(7): 2836-2840. http://d.wanfangdata.com.cn/Periodical/zngydxxb201307028
  • [1] 任千秋, 林然, 赵金洲, 吴建发, 宋毅.  基于初次压裂地质工程条件及改造效果的页岩气井重复压裂选井模型 . 石油钻采工艺, 2023, 45(2): 223-228, 243. doi: 10.13639/j.odpt.2023.02.014
    [2] 李冬梅, 柳志翔, 李林涛, 石鑫.  顺北超深断溶体油气藏完井技术 . 石油钻采工艺, 2020, 42(5): 600-605. doi: 10.13639/j.odpt.2020.05.013
    [3] 李准, 吴晓东, 韩国庆, 王瑞河, 郑磊.  考虑裂缝导流能力时效性的多级压裂水平井产能半解析模型 . 石油钻采工艺, 2019, 41(3): 354-362. doi: 10.13639/j.odpt.2019.03.015
    [4] 时贤, 程远方, 常鑫, 许洪星, 吴百烈, 蒋恕.  页岩气水平井段内多簇裂缝同步扩展模型建立与应用 . 石油钻采工艺, 2018, 40(2): 247-252. doi: 10.13639/j.odpt.2018.02.018
    [5] 周杨, 李莉, 吴忠宝, 甘俊奇, 王俊文, 刘翀.  低渗透油藏垂直裂缝井产能预测及分析 . 石油钻采工艺, 2017, 39(2): 146-150. doi: 10.13639/j.odpt.2017.02.004
    [6] 侯腾飞, 张士诚, 马新仿, 李栋, 孙延安.  支撑剂沉降规律对页岩气压裂水平井产能的影响 . 石油钻采工艺, 2017, 39(5): 638-645. doi: 10.13639/j.odpt.2017.05.019
    [7] 蒋贝贝, 李海涛, 王跃曾, 张春杰.  双层分支水平井注热海水开采海底 天然气水合物经济性评价 . 石油钻采工艺, 2015, 37(1): 87-91. doi: 10.13639/j.odpt.2015.01.022
    [8] 郭小哲, 周长沙.  基于滑脱的页岩气藏压裂水平井渗流模型及产能预测 . 石油钻采工艺, 2015, 37(3): 61-65. doi: 10.13639/j.odpt.2015.03.014
    [9] 刘晓华, 邹春梅, 姜艳东, 翟振宇.  页岩气水平井动态评价方法   . 石油钻采工艺, 2013, 35(3): 55-58.
    [10] 张睿, 韩国庆, 宁正福, 廖新维.  蛇曲井渗流特征与产能评价电模拟实验研究   . 石油钻采工艺, 2013, 35(6): 60-64.
    [11] 胡艾国, 熊佩, 姚昌宇, 李月丽.  红河油田长8油藏压裂水平井产能影响因素分析 . 石油钻采工艺, 2013, 35(3): 69-72.
    [12] 赵振峰, 吴晓东, 黄伟, 吕亿明, 周赤烽, 孔令琴.  特低渗透油藏压裂水平井产能影响因素分析 . 石油钻采工艺, 2011, 33(4): 81-83.
    [13] 陈冀嵋, 文浩, 徐乐, 杜岚, 蒋宜春.  酸压井酸蚀蚓孔区对产能的影响 . 石油钻采工艺, 2010, 32(4): 57-60.
    [14] 梁峰, 孟维龙, 刘晓辉.  分支水平井产能方程式研究进展 . 石油钻采工艺, 2010, 32(5): 64-68.
    [15] 王晓泉, 张守良, 吴奇, 刘玉章, 雷群, 丁云宏.  水平井分段压裂多段裂缝产能影响因素分析 . 石油钻采工艺, 2009, 31(1): 73-76.
    [16] 张磊, 王克雄, 程娣, 张航, 张锐, 程汝才.  高29支平1分支水平井钻井技术 . 石油钻采工艺, 2007, 29(1): 7-11. doi: 10.3969/j.issn.1000-7393.2007.01.003
    [17] 程林松, 张健琦, 李春兰.  水平井井网整体开发产能研究 . 石油钻采工艺, 2002, 24(2): 39-41. doi: 10.3969/j.issn.1000-7393.2002.02.013
    [18] 蒋廷学, 李安启, 姜东.  考虑井筒流动的垂直裂缝井稳态产能计算模型 . 石油钻采工艺, 2001, 23(4): 50-53,85. doi: 10.3969/j.issn.1000-7393.2001.04.017
    [19] 孙艾茵, 刘蜀知, 任书泉.  砂岩射孔井基质酸化数学模型的建立及应用 . 石油钻采工艺, 1997, 19(5): 70-75. doi: 10.3969/j.issn.1000-7393.1997.05.016
    [20] 李笑萍, 张丽囡, 翟云芳.  变井储气井试井的数学模型与解释方法 . 石油钻采工艺, 1993, 15(3): 71-78. doi: 10.3969/j.issn.1000-7393.1993.03.014
  • 加载中
图(7)
计量
  • 文章访问数:  590
  • HTML全文浏览量:  117
  • PDF下载量:  853
  • 被引次数: 0
出版历程
  • 刊出日期:  2017-02-10

基于分形理论的页岩气分支水平井产能数学模型

doi: 10.13639/j.odpt.2017.01.002
    基金项目:  国家科技重大专项“井震结合油藏精细结构表征技术”(编号:2011ZX05010-001);陕西省高水平大学建设专项资金资助项目(编号:2013SXTS03);延安大学专项科研基金项目(编号:YDK2012-1)
    作者简介:

    李忠厚(1982-),2008年毕业于中国石油大学(华东)油气田开发工程专业,现为延安大学讲师,从事物理法采油,非常规油气田开发方面的研究工作。通讯地址:(716000)中国陕西延安市圣地路580号延安大学石油工程与环境工程学院。电话:18329909987。E-mail:wo119944595er@163.com

  • 中图分类号: TE328

摘要: 利用分形几何法,研究了页岩储集层内水平井分支数目和长度的分形维度,以及水平井各分支井段压降扰动的传播规律,得到压降漏斗边界随时间变化的关系,考虑了页岩储集层解吸吸附、扩散系数和启动压力梯度的影响,建立了页岩气分支水平井产能数学模型。采用拉普拉斯变换法,求解了水平井外边界定压条件下的不稳定渗流的产能方程。结合四川某海相页岩气藏储集层参数,计算分析了页岩气分支水平井产能特征及影响因素。研究表明:在页岩气生产过程中,同一分形维度下,水平井分支数目越大,气井的产量增大越快;分支井的井段长度越大,对主干井的流量影响越大。吸附页岩气的解吸使得压降漏斗的边界传播速度减慢,地层压力下降缓慢;产气量越大,井底压力随时间下降越快,但是下降的趋势随时间逐渐减缓。

English Abstract

李忠厚, 吴小斌, DUZhongwei, 任茜莹. 基于分形理论的页岩气分支水平井产能数学模型[J]. 石油钻采工艺, 2017, 39(1): 7-13. doi: 10.13639/j.odpt.2017.01.002
引用本文: 李忠厚, 吴小斌, DUZhongwei, 任茜莹. 基于分形理论的页岩气分支水平井产能数学模型[J]. 石油钻采工艺, 2017, 39(1): 7-13. doi: 10.13639/j.odpt.2017.01.002
LI Zhonghou, WU Xiaobin, Zhongwei DU, REN Qianying. A mathematical model for productivity calculation of shale-gas multilateral horizontal wells based on fractal theory[J]. Oil Drilling & Production Technology, 2017, 39(1): 7-13. doi: 10.13639/j.odpt.2017.01.002
Citation: LI Zhonghou, WU Xiaobin, Zhongwei DU, REN Qianying. A mathematical model for productivity calculation of shale-gas multilateral horizontal wells based on fractal theory[J]. Oil Drilling & Production Technology, 2017, 39(1): 7-13. doi: 10.13639/j.odpt.2017.01.002
  • 页岩气水平井流动特征是气藏工程、产能分析、气藏数值模拟以及采气工程研究的基础。页岩气在纳米~微米级页岩孔隙及裂缝中流动时, 主要为解吸吸附、扩散、滑脱等, 为低速强非线性渗流[1-5]。郭小哲[6]等在实验的基础上认为, 页岩对内压具有很强的敏感效应; 朱维耀等[7-9]认为, 页岩气储集层渗流具有启动压力梯度; 分支水平井在生产时, 气体在地层中呈三维流向主干井和分支井, 其渗流规律就是一个页岩气平面径向流的压降漏斗动边界问题; 然后页岩气又沿着分支井与主干井汇合, 最后沿着主干井流向垂直井的井底。分支水平井的产量方程的建立和求解较难, 考虑压降漏斗动边界扰动问题则更难。中国页岩气整体处于工业起步阶段, 虽然目前国内外通过压裂手段对页岩气的开发已进入快速发展阶段, 然而由于页岩气储层的纳米级孔隙复杂, 对页岩孔隙分形特征和水平井产能预测的研究, 大多数研究成果给出的仅是页岩气压裂改造的渗流规律及其影响因素, 尽管也有部分文献提出了具体的渗流方程, 但考虑的因素尚少[10-15], 方程过于简化, 难于更好地反映分支水平井的产能变化规律。为此, 有必要揭示其页岩气的非线性流动规律, 研究页岩气分支水平井产能变化特征, 以便选取合理有效的开发方式和增产手段, 为页岩气的开发提供理论依据。

    基于分形理论对页岩储集层内水平井主干井和各分支井段长度和分支数目问题进行研究, 进而建立考虑了页岩储集层解吸吸附、水平井分支数目和长度影响的页岩气分支水平井不稳定渗流数学模型, 并推导和求解。结合中国四川某海相页岩储集层参数, 分析页岩气分支水平井产能变化特征及其影响因素。

    • 当页岩气分支水平井投入开发, 打开储集层后, 页岩气的流动不仅有各井段渗流过程, 还存在页岩气扩散、滑移、解吸流动, 流动阻力比常规天然气大。水平井主干井开采形成的压力降将逐渐向外传播, 地层各点压力分布是绕井轴旋转构成的曲面, 如图 1所示, 设某时刻t, 压力传播到地层r处, 在r范围内形成压降漏斗; r即为页岩气渗流过程中的压降漏斗边界。r随时间逐渐增大, 压降漏斗边界上的压力等于原始地层压力。

      图  1  页岩气分支水平井平面径向流示意图

      Figure 1.  Schematic areal radial flow of shale-gas multilateral horizontal well

    • 页岩气藏是自生自储式气藏, 是一种非均匀性的多孔储集层, 主要为纳米~微米级孔隙, 具有复杂的孔隙结构。页岩气赋存方式主要是分散在孔隙和裂缝中的游离态和吸附在干酪根和黏土矿物颗粒上的吸附态两大类[16], 以直井方式单纯依靠天然的孔-缝渗流不可能形成较大规模的产量。为了提高单井产量并获得较高的采收率, 分支水平井技术在国内外得到了广泛的应用。当以分支水平井开始生产时, 分支水平井的地层气体流量的变化不仅有以往水平井的地层气体流量的变化, 还存在水平井各分支之间地层气体流量的变化。因此, 在不稳定渗流过程中, 分支水平井地层气体流量的变化是一个扰动压降漏斗边界问题。

      对于气体流动已偏离达西定律的纳米~微米孔隙页岩储集层, 朱维耀等[16-20]建立了考虑扩散、滑脱系数的具有复杂孔隙结构的气体流动方程

      $$ v=-\frac{K_{0}}{\mu}\left(1+\frac{k_{\mathrm{B}} T}{\sqrt{2} \pi \delta^{2}} \frac{b}{p}\right) \frac{\mathrm{d} p}{\mathrm{d} x} $$ (1)

      则任一瞬间, 以水平井某分支为压降漏斗的中心井轴时, 附近地层压力为

      $$ \begin{aligned} & \left(p_{\mathrm{w} m}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}=\left(p_{\mathrm{e}}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}-\\ & \frac{\left(p_{\mathrm{e}}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}-\left(p_{\mathrm{w} m}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}}{\ln \frac{R_{1(m)}(t)}{r_{\mathrm{w}}}} \ln \frac{R_{1(m)}(t)}{r} \end{aligned} $$ (2)

      任一瞬间, 以水平井主干井为压降漏斗的中心井轴时, 附近地层压力为

      $$ \begin{array}{l} \left(p+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}=\left(p_{\mathrm{e}}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}- \\ \sum_{m=1}^{n} \frac{\left(p_{\mathrm{e}}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}-\left(p_{\mathrm{w} m}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}}{\ln \frac{R_{1(m)}(t)}{r_{\mathrm{w} m}} \ln \frac{R_{1(m)}(t)}{r}} \end{array} $$ (3)

      如图 2所示, 在地层中距离水平井某一分支为r处取微小圆环体, 水平长度为l, 厚度为dr, 此单元体中游离态气体的原始储量为2πrlφiρidr。某一任意时刻t, 该单元体中残留气体储量为2πrlφρdr, 因此, 从单元体孔隙中采出的游离态气体储量为2πrl(φiρiφρ)dr。考虑吸附态气体的解吸, 采出的总气体量为

      $$ Q_{\mathrm{sc}}=\sum_{m=1}^{n} \int_{r_{\mathrm{w} m}}^{R_{1(m)}(t)} 2 \pi r l\left[\left(\phi_{1} \rho_{\mathrm{i}}-\phi \rho\right)+q_{\mathrm{d}} t\right] \mathrm{d} r $$ (4)

      图  2  页岩气水平井分支平面径向流示意图

      Figure 2.  Schematic areal radial flow of branch holes of shale-gas horizontal well

      联立式(3)、式(4)代入气体状态方程, 进一步求得

      $$ \begin{aligned} Q_{\mathrm{sc}}=& \frac{\pi l \phi C T_{\mathrm{sc}} Z_{\mathrm{sc}} \rho_{\mathrm{gsc}}}{p_{\mathrm{sc}} T Z} \sum_{m=1}^{n}\left[\left(p_{\mathrm{e}}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}-\left(p_{\mathrm{w} m}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}\right] \times \\ & \frac{R_{1(m)}^{2}(t)-r_{\mathrm{w} m}^{2}}{2 \ln \frac{R_{1(m)}(t)}{r_{\mathrm{w} m}}}+\sum_{m=1}^{n} \pi l q_{\mathrm{d}} t\left[R_{1(m)}^{2}(t)-r_{\mathrm{w} m}^{2}\right] \end{aligned} $$ (5)

      按气体稳定渗流公式, 页岩气分支水平井流量可写为

      $$ q_{\mathrm{sc1}}=\frac{\pi K_{0} l T_{\mathrm{sc}} Z_{\mathrm{sc}} \rho_{\mathrm{gsc}}}{p_{\mathrm{sc}} T \mu Z} \sum_{m=1}^{n} \frac{\left(p_{\mathrm{e}}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}-\left(p_{\mathrm{w} m}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}\right)^{2}}{\ln \frac{R_{1(m)}(t)}{r_{\mathrm{w} m}}} $$ (6)

      假设水平井采用定产量生产, 即qsc1=const, 联立式(5)、式(6)可得

      $$ Q_{\mathrm{sc}}=\frac{\phi \mu C q_{\mathrm{sc1}}}{2 K_{0}} \sum_{m=1}^{n}\left(R_{1(m)}^{2}(t)-r_{\mathrm{w} m}^{2}\right)+\sum_{m=1}^{n} \pi l q_{\mathrm{d}} t\left(R_{1(m)}^{2}(t)-r_{\mathrm{w} m}^{2}\right) $$ (7)

      由此根据Qsc=qsc1t和式(7)得出页岩储集层任一时刻压降漏斗边界与水平井分支参数n变化的关系为

      $$ \sum_{m=1}^{n}\left(R_{1(m)}^{2}(t)-r_{\mathrm{w} m}^{2}\right)=\frac{q_{\mathrm{scl}} t}{\frac{\phi \mu C q_{\mathrm{sc} 1}}{2 K_{0}}+\pi l q_{\mathrm{d}} t} $$ (8)

      中国四川海相页岩气F气藏孔隙度为3.0%~6.0%, 绝对渗透率(0.3~7)×10-4mD, 真实气体压缩因子0.87, 气体黏度0.031 mPa·s, 平均解吸量4.611 7× 10-10kg/(m3·s), 储层压力为37.57 MPa, 井筒半径0.1 m, 水平井段长度1 182 m, 储层深度2 785 m, 气藏厚度80.5 m, 定产量18.2× 104m3/d。

      利用四川F气藏基本参数, 根据式(8)绘制水平井不同分支数目条件下页岩气储集层压降漏斗边界随时间变化关系图(图 3)。同一时刻, 分支井分支数目越多, 压降漏斗动用的区域越大; 在生产初期, 各个压降漏斗边界单独传播, 不产生扰动, 在生产进行到一定阶段后, 曲线上出现拐点, 表明不同分支产生的压降漏斗边界产生相互干扰, 拐点的参数表明产生干扰的次数和强度, 因此页岩气储集层中压力传播需要考虑水平井分支的影响。

      图  3  不同分支数下压降漏斗边界与分支井半径比值随时间变化关系

      Figure 3.  Relationship of the ratio of pressure drop funeral boundary to branch hole radius vs. time for different branch hole amounts

    • 在页岩气开采过程中, 由于页岩孔隙结构对页岩气赋存和运移的影响, 宏观上气体流动出现类似低渗油藏压力启动梯度现象。因此, 建立考虑启动压力梯度的稳态径向流气体常微分方程

      $$ \left\{\begin{array}{l} \frac{\mathrm{d}^{2} F}{\mathrm{d} r^{2}}+\frac{1}{r} \frac{\mathrm{d} F}{\mathrm{d} r}-\lambda_{\mathrm{t}} \frac{\mathrm{d} F}{\mathrm{d} r}=0 \\ r=r_{\mathrm{w} m}, r \frac{\partial F}{\partial r}=\frac{q_{\mathrm{sc} 2} T p_{\mathrm{sc}}}{\pi K_{0} l Z_{\mathrm{sc}} T_{\mathrm{sc}}} \\ r=r, F=F_{R_{2(m)}} \end{array}\right. $$ (9)

      其中

      $$ \begin{array}{c} \lambda_{\mathrm{t}}=C \lambda \\ F=2 \int_{p_{\mathrm{a}}}^{p} \frac{p}{\mu Z} \mathrm{d} p \end{array} $$

      求解得到考虑启动压力梯度的页岩气分支水平井渗流表达式

      $$ \begin{aligned} q_{\mathrm{sc} 2}=& \pi l \frac{K_{0} T_{\mathrm{sc}} Z_{\mathrm{sc}} \rho_{\mathrm{gsc}}}{p_{\mathrm{sc}} T \mu Z} \times \\ & \sum_{m=1}^{n} \frac{p_{\mathrm{e}}^{2}-p_{\mathrm{w} m}^{2}}{\exp \left(-\lambda_{\mathrm{t}} r_{\mathrm{w} m}\right)\left[\operatorname{Ei}\left(-\lambda_{\mathrm{t}} r_{\mathrm{e}}\right)-\operatorname{Ei}\left(-\lambda_{\mathrm{t}} r_{\mathrm{w} m}\right)\right]} \end{aligned} $$ (10)

      考虑解吸, 地层中采出的总气体量

      $$ \begin{array}{l} Q_{\mathrm{sc}}=\frac{\pi l \phi C T_{\mathrm{sc}} Z_{\mathrm{sc}} \rho_{\mathrm{gsc}}}{p_{\mathrm{sc}} T Z} \times \\ \sum_{m=1}^{n} \frac{\left(p_{\mathrm{e}}^{2}-p_{\mathrm{w} m}^{2}\right)\left[\mathrm{e}^{\lambda_{1} R_{2(m)}(t)}\left(\lambda_{\mathrm{t}} R_{2(m)}(t)-1\right)-\mathrm{e}^{\lambda_{1} r_{\mathrm{wm}}}\left(\lambda_{\mathrm{t}} r_{\mathrm{w} m}-1\right)\right]}{\lambda_{\mathrm{t}}^{2}\left[\operatorname{Ei}\left(-\lambda_{\mathrm{t}} R_{2(m)}(t)\right)-\operatorname{Ei}\left(-\lambda_{\mathrm{t}} r_{\mathrm{w} m}\right)\right]}+ \\ \sum_{m=1}^{n} \pi l q_{\mathrm{d}}\left[R_{2(m)}^{2}(t)-r_{\mathrm{w} m}^{2}\right] \end{array} $$ (11)

      将式(11)代入式(10), 得到页岩储集层考虑启动压力梯度压降边界与时间的关系为

      $$ \begin{aligned} t=& \frac{\phi \mu C}{K_{0}} \sum_{m=1}^{n} \frac{\left.\exp \left(\lambda_{\mathrm{t}} R_{2\left(r_{\mathrm{L}_{\mathrm{m}}, \theta_{m}}\right.}\right)(t)\right)}{\lambda_{\mathrm{t}}^{2}}\left(\lambda_{\mathrm{t}} R_{2(m)}(t)-\lambda_{\mathrm{t}} r_{\mathrm{w} m}\right)+\\ & \frac{\pi l q_{\mathrm{d}} t}{q_{\mathrm{sc} 2}} \sum_{m=1}^{n}\left[R_{2(m)}^{2}(t)-r_{\mathrm{w} m}^{2}\right] \end{aligned} $$ (12)

      图 4为不同启动压力梯度下压降边界随时间变化关系, 可以看出, 在同一时刻, 启动压力梯度越大, 压降漏斗传播边界越小。

      图  4  不同启动压力梯度下压降漏斗边界与主干井半径比值随时间变化关系

      Figure 4.  Relationship of the ratio of pressure drop funeral boundary to main hole radius vs. time for different startup pressure gradients

    • 分形理论由于其具有描述不规则物体的优越性和精确性, 是研究自然界中不规则性和不均匀性几何形态的有力工具。页岩气水平井分支在地层是一种不规则分布, 分支段可以认为是水平井段的自相似性的分形, 计算分支水平井的分形维数, 就可以近似计算分支水平井的产能模型。

      分支水平井的水平井段长度为L, 分支数目为n, 分支段长度为Lm; 以分支井与主干井的交点为分割点, 水平井段分为数目n+1与水平井相似的分形, 分支井为n个与水平井相似的分形, 则分支水平井的分形维数为

      $$ D_{n}=\frac{\ln (2 n+1)}{\ln \frac{L}{L_{m}}} $$ (13)

      在考虑页岩储集层解吸吸附和水平井分支数目、长度的基础上, 利用分形理论的分支井长度的分形维数法, 建立页岩气分支水平井不稳定渗流模型, 并进行推导和求解。

      基于页岩气渗流的连续性方程、运动方程和状态方程以及页岩气解吸吸附的Langmuir模型[17-18], 考虑不稳定渗流过程中水平井分支和启动压力梯度对压降漏斗传播的影响, 建立页岩气分支水平井不稳定渗流控制方程

      $$ \begin{aligned} & \frac{1}{r} \frac{\partial}{\partial r}\left[r \frac{1}{\mu Z}\left(p+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}-\lambda_{\mathrm{t}}\right) \frac{\partial p}{\partial r}\right]=\\ & \frac{\phi \mu C}{K_{0}} \frac{1}{\mu Z} \frac{\partial p}{\partial t}+\frac{p_{\mathrm{sc}} T p_{\mathrm{L}} V_{\mathrm{L}}}{T_{\mathrm{sc}} Z_{\mathrm{sc}} K_{0}\left(p+p_{\mathrm{La}}\right)^{2}} \frac{\partial p}{\partial t} \end{aligned} $$ (14)

      引入拟压力函数和无因次量, 并定义如下

      $$ \begin{array}{c} \Psi(p)=2 \int_{p_{\mathrm{a}}}^{p} \frac{1}{\mu Z}\left(p+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}-\lambda_{\mathrm{t}}\right) \mathrm{d} p \\ r_{\mathrm{D}}=\frac{r}{r_{\mathrm{w}}}, \quad r_{\mathrm{Dw}}=\frac{r_{\mathrm{w}}}{r_{\mathrm{e}}} \\ q_{\mathrm{D}}=\frac{p_{\mathrm{sc}} T}{\pi K_{0} l Z_{\mathrm{sc}} T_{\mathrm{sc}}\left(\Psi_{\mathrm{i}}-\Psi_{\mathrm{w}}\right)} q_{\mathrm{sc}} \\ \Psi_{\mathrm{D}}=\frac{\pi K_{0} l Z_{\mathrm{sc}} T_{\mathrm{sc}}}{p_{\mathrm{sc}} q_{\mathrm{sc}} T}\left(\Psi_{\mathrm{i}}-\Psi\right) \\ t_{\mathrm{D}}=\left(p+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}-\lambda_{\mathrm{t}}\right) \times \\ {\left[1+\frac{p_{\mathrm{sc}} T p_{\mathrm{L}} V_{\mathrm{L}}}{T_{\mathrm{sc}} Z_{\mathrm{sc}} K_{0}\left(p+p_{\mathrm{L}}\right)^{2}} \frac{Z K_{0}}{\phi C}\right]^{-1} \frac{K_{0} t}{\phi \mu C r_{\mathrm{w} m}^{2}}} \end{array} $$

      则(14)式可转换为

      $$ \frac{1}{r_{\mathrm{D}}} \frac{\partial}{\partial r_{\mathrm{D}}}\left(r_{\mathrm{D}} \frac{\partial \Psi_{\mathrm{D}}}{\partial r_{\mathrm{D}}}\right)=\frac{\partial \Psi_{\mathrm{D}}}{\partial t_{\mathrm{D}}} $$ (15)

      当页岩气水平井以某一恒定产量生产时, 外边界定压, 则其定解条件为

      $$r_{\mathrm{D}} \frac{\partial \Psi_{\mathrm{D}}}{\partial r_{\mathrm{D}}}=q_{\mathrm{D}}$$ (16)
      $$\left.\Psi_{\mathrm{D}}\right|_{r_{\mathrm{D}}=\infty}=0$$ (17)
      $$\left.\Psi_{\mathrm{D}}\right|_{t_{\mathrm{D}}=0}=0$$ (18)

      对于分支数为n的页岩气水平井, 运用拉普拉斯变换法求解其在气藏中的拟压力函数可表示为

      $$ \begin{aligned} & \psi_{\mathrm{w}}=\psi_{\mathrm{i}}-\frac{q_{\mathrm{sc}} p_{\mathrm{sc}} T}{\pi K_{0} l Z_{\mathrm{sc}} T_{\mathrm{sc}}} \times\\ & \left\{2 \ln \frac{R(m)}{r_{\mathrm{w}}}-4 \sum_{m=1}^{\infty}\left[\frac{1}{\beta_{\mathrm{n}}^{2}} \frac{\exp \left(-\beta_{\mathrm{n}}^{2} \frac{K_{0} t_{\mathrm{D}}}{\phi \mu C_{\mathrm{t}} r_{\mathrm{w}}^{2}}\right)}{\mathrm{~J}_{1}^{2}\left(\beta_{\mathrm{n}}\right)-\mathrm{J}_{0}^{2}\left(\frac{R(m)}{r_{\mathrm{w}}} \beta_{\mathrm{n}}\right)} \mathrm{J}_{0}^{2}\left(\frac{R(m)}{r_{\mathrm{w}}} \beta_{\mathrm{n}}\right)\right]\right\} \end{aligned} $$ (19)

      其中, βn为方程(20)的根。

      $$ \mathrm{J}_{1}\left(\frac{r_{\mathrm{e}}}{r_{\mathrm{w}}} \beta\right) Y_{1}(\beta)-Y\left(\frac{r_{\mathrm{e}}}{r_{\mathrm{w}}} \beta\right) J_{1}(\beta)=0 $$ (20)

      将拟压力函数代入式(19), 求得分支井水平井的产量为

      $$ \begin{aligned} & Q_{\mathrm{sc}}=q_{\mathrm{sc}}^{D_{n}} \frac{\pi K_{0} l Z_{\mathrm{sc}} T_{\mathrm{sc}}}{p_{\mathrm{sc}} T}\left[\left(p_{\mathrm{i}}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}-\lambda_{\mathrm{t}}\right)^{2}-\left(p_{\mathrm{w}}+\frac{k_{\mathrm{B}} T b}{\sqrt{2} \pi \delta^{2}}-\lambda_{\mathrm{t}}\right)^{2}\right] \times\\ & \left\{2 \ln \frac{R(m)}{r_{\mathrm{w}}}-4 \sum_{m=1}^{\infty}\left[\frac{1}{\beta_{\mathrm{n}}^{2}} \frac{\exp \left(-\beta_{\mathrm{n}}^{2} \frac{K_{0} t_{\mathrm{D}}}{\phi \mu C_{\mathrm{t}} r_{\mathrm{w}}^{2}}\right)}{\mathrm{~J}_{1}^{2}\left(\beta_{\mathrm{n}}\right)-J_{0}^{2}\left(\frac{R(m)}{r_{\mathrm{w}}} \beta_{\mathrm{n}}\right)} \mathrm{J}_{0}^{2}\left(\frac{R(m)}{r_{\mathrm{w}}} \beta_{\mathrm{n}}\right)\right]\right\}^{-1} \end{aligned} $$ (21)
    • 根据上文推导出的考虑水平井分支、启动压力梯度、解吸吸附的页岩气分支水平井产量变化规律, 结合四川F气藏基本参数, 对页岩气分支水平井不稳定渗流压力分布及其影响因素进行分析。

      图 5为不同分支长度对气井生产的影响曲线, 纵坐标表示分形维度计算下的气井产量QDn与单一水平井气井产量Q之比。由图可见, 同一分形维度下, 水平井分支数目越大, 气井的产量增大越快; 分支井的井段长度越大, 对主干井的流量影响越大。

      图  5  不同分支长度对气井生产的影响

      Figure 5.  Effect of branch hole length on gas well production

      图 6为吸附解吸气体量对气井生产的影响。由图可见, 考虑吸附页岩气解吸量时分支水平井日产量变化幅度不大, 分支水平井的累计产量变化较大, 具有明显的产量提升; 吸附页岩气的解吸使得页岩气的产量增加。

      图  6  吸附解吸气体量对气井生产的影响

      Figure 6.  Effect of adsorption/desorption gas volume on gas well production

      图 7为基于分形模型的气井模拟值与实际产量对比图。图中水平井的水平段长度1 182 m, 分支段长度450 m, 分支数目为1, 带入式(13)计算出分形维数Dn为1.14.根据同地层单一水平井的生产数据带入分形维数产能模型计算出分支井的模拟数据, 并与分支井的实际数据相比较。由图可见, 分支井的产量初期明显高于单一水平井, 模拟值与实际产量拟合较好; 由于水平分支井与主干井之间的相互干扰, 分支井产量递减较快; 随着生产时间的增加, 模拟值与实际产量变化逐渐平稳, 并且拟合较好。

      图  7  基于分形模型的气井模拟值与实际产量对比

      Figure 7.  Comparison between actual production and simulated value of gas wells based on fractal model

    • (1) 利用几何分形法得出了分支水平井井段长度的有效分形维数表达式。同一分形维度下, 水平井分支数目越大, 气井的产量增大越快; 分支井的井段长度越大, 对主干井的流量影响越大。

      (2) 考虑吸附页岩气解吸量时分支水平井日产量变化幅度不大, 分支水平井的累计产量变化较大, 具有明显的产量提升

      (3) 基于分形理论的分支井模型模拟值与实际产量的拟合曲线, 在生产初期和递减期拟合较好, 中间拟合度较差。

      符号说明:

      Symbol description:

      a为修正系数; b为滑脱系数, m2/s; C为气体的等温压缩系数, Pa-1; Dn为水平井段长度的分形维数; Ei为幂积分函数; F为拟压力函数, Pa/s; J0为第一零阶贝塞尔函数; J1为一阶贝塞尔函数; K0为储集层绝对渗透率, m2; kB为玻尔兹曼常数; L为气层水平井长度, m;Lm为气层水平分支长度, m;l为气层水平长度, m;n为水平井分支数(m=0, 1, 2, 3, …, n);p为储集层压力, MPa; pa为某一已知压力, MPa; pe为原始地层压力, MPa; pL为Langmuir压力, 表示吸附量为最大吸附量一半时的压力, Pa;psc为标准状态下压力, MPa; pwm为水平井分支井底压力, MPa; Qsc为地层中采出的总气体量, m3; qd为单位体积页岩单位时间的解吸量, kg/(m3·s); qsc1为标准条件下气井流量, m3/s; qsc2为启动压力梯度影响下气井流量, m3/s; R1(t)为水平井分支影响压降漏斗边界, m;R2(t)为启动压力梯度影响压降漏斗边界, m;r为压降漏斗边界距井筒距离, m;re为气井供给半径, m;rw为井筒半径, m;rwm为分支井井筒半径, m;rLm为水平井第m分支井径向距离, m;T为地层温度, K;Tsc为标准状态下温度, K;t为生产时间, d;VL为Langmuir体积, 表示最大吸附量, m3; v为气体渗流速度, m/s; Z为气体压缩因子, 无因次; Zsc为标准状态下气体压缩因子, 无因次; δ为气体分子的碰撞直径, m;θm为水平井第m分支与主干井夹角角度, rad; λ为启动压力梯度, Pa/m; μ为气体黏度, Pa·s; ρ为气体密度, kg/m3; ρi为地层原始气体密度, kg/m3; ρgsc为标准状态下气体密度, kg/ m3; φ为储集层孔隙度, 小数; φi为储集层原始孔隙度; Ψ为非线性拟压力函数, Pa/s; Ψw为井底非线性拟压力函数, Pa/s。下标D表示无因次。

参考文献 (20)

目录

    /

    返回文章
    返回