›› 2015, Vol. 36 ›› Issue (11): 3315-3321.doi: 10.16285/j.rsm.2015.11.037

• 数值分析 • 上一篇    下一篇

考虑非达西效应的酸蚀裂缝流场数值模拟

郭建春1,庄 园2,刘 超3   

  1. 1. 西南石油大学 油气藏地质及开发工程国家重点实验室,四川 成都 610500;2. 中石化西南油气分公司川东北采气厂,四川 阆中 637400; 3. 中石油长庆油田分公司油气工艺研究院,陕西 西安,710018
  • 收稿日期:2015-04-12 出版日期:2015-11-11 发布日期:2018-06-14
  • 通讯作者: 庄园,男,1988年生,硕士,主要从事油开采和增产技术方面的研究工作。E-mail: 88194754@qq.com E-mail:guojianchun@vip.163.com
  • 作者简介:郭建春,男,1970年生,博士,教授,博士生导师,主要从事油开采和增产技术方面的研究工作。
  • 基金资助:

    国家科技重大专项(No.2011ZX05045-03);四川省青年科技创新研究团队资助计划(No.2011JTD0018)。

Numerical simulation of flow field of acid etched fractures considering non-Darcy effect

GUO Jian-chun 1, ZHUANG Yuan 2, LIU Chao 3   

  1. 1. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu, Sichuan 610500, China; 2. Northeast Sichuan Gas Recovery Plant, Southwest Oil & Gas Company, Sinopec, Langzhong, Sichuan 637400, China; 3. Research Institute of Oil/Gas Technology, Changqing Oilfield Branch Company, Sinopec, Xi'an, Shaanxi 710018, China
  • Received:2015-04-12 Online:2015-11-11 Published:2018-06-14
  • Supported by:

    Project supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China (Grant No.2011ZX05045-03) and the Subsidy Scheme for Youth Science and Technology Innovation Team of Sichuan Provincial(Grant No.2011JTD0018).

摘要: 酸压改造后的酸蚀裂缝是储层流体的主要流动通道,研究流体在酸蚀裂缝中的流动规律是十分必要的。本文介绍了目前裂隙流动研究中的4种常用控制方程及各自的适用条件,确定了Navier-Stokes方程为酸蚀裂缝流场数值模拟的控制方程。为了分析酸蚀裂缝流场分布规律以及非达西效应对裂缝导流能力的影响,利用逆向工程技术对两种不同刻蚀形态的酸蚀岩样壁面进行了实体重构,并采用有限元数值模拟方法进行了不同流量下的酸蚀裂缝流动试验。结果表明:酸刻蚀形态和裂缝接触关系对流场分布的影响较大。隙宽分布平滑的均匀刻蚀裂缝,其流态稳定、曲折度低,但由于裂缝开度较窄且比表面较高,边界层效应明显,流体流动时产生的黏滞阻力较高;沟槽裂缝曲面粗糙,隙宽分布复杂,流体流动时流态不稳定,曲折度高,在高流量下会产生明显的涡流流动,增大惯性阻力。裂缝通道存在缩颈现象会引起流体绕流和多次加减速,产生额外的压力损耗。随着模拟流量的增大,压降与流量会逐渐偏离线性关系而呈现出非达西流动现象。酸蚀裂缝壁面越粗糙,产生非达西效应时的临界流量和临界雷诺数就越小,在相同流量下非达西效应就越强烈,导流能力下降速率更快。

关键词: 酸蚀裂缝, 数值模拟, 流场分布, 非达西效应, 导流能力

Abstract: Acid etched fracture caused by acid fracturing is the main flow channel for fluid, thus, it is necessary to study the flowing rules of fluid in acid etched fracture. Four governing equations commonly used in fracture flow research and their applied conditions are introduced, among which Navier-Stokes equation is chosen as the governing equation for numerical simulation on acid etched fracture. In order to analyze the behavior of fluid flow through acid etched fractures and high-speed non-Darcy effect on fracture conductivity capacity, two kinds of fractures with different acid etching morphologies are reconstructed to solid models by using reverse engineering techniques, and flow experiments are carried out under different flow rates by using finite element numerical simulation. It indicates that, the acidized surface morphology and contact area of acid etched fracture have a great impact on the flow field distribution. The uniform etching fracture with smooth aperture distribution leads to stable flow pattern and low tortuosity, however, with the presence of obvious boundary layer effect caused by narrow fracture width and wide surface area, high viscous resistance will be generated while flowing. The channeling fracture with rough surface and complex aperture distribution results in unstable flow pattern and high tortuosity, and there is an obvious vortex under high flow rate, hence, increasing the inertial resistance loss. The necking phenomenon in contact area leads to circumfluence, and multiple acceleration-decelerations cause additional pressure loss. Meanwhile, with the increase in simulating flow rate, pressure drop and flow rate will gradually deviate from the linear relationship and present the non-Darcy flow phenomenon. The rougher the fracture surface is, the smaller the critical flow and critical Reynolds number are, the stronger the non-Darcy effect of fluid is under the same flow rate, the faster the conductivity capacity decreases.

Key words: acid etched fracture, numerical simulation, flow field distribution, non-Darcy effect, conductivity capacity

中图分类号: 

  • O 241
[1] 李翻翻, 陈卫忠, 雷江, 于洪丹, 马永尚, . 基于塑性损伤的黏土岩力学特性研究[J]. 岩土力学, 2020, 41(1): 132-140.
[2] 夏 坤, 董林, 蒲小武, 李璐, . 黄土塬地震动响应特征分析[J]. 岩土力学, 2020, 41(1): 295-304.
[3] 修乃岭, 严玉忠, 胥云, 王欣, 管保山, 王臻, 梁天成, 付海峰, 田国荣, 蒙传幼, . 基于非达西流动的自支撑剪切裂缝导流 能力实验研究[J]. 岩土力学, 2019, 40(S1): 135-142.
[4] 郭院成, 李明宇, 张艳伟, . 预应力锚杆复合土钉墙支护体系增量解析方法[J]. 岩土力学, 2019, 40(S1): 253-258.
[5] 闫国强, 殷跃平, 黄波林, 张枝华, 代贞伟, . 三峡库区巫山金鸡岭滑坡成因机制与变形特征[J]. 岩土力学, 2019, 40(S1): 329-340.
[6] 刘红岩. 宏细观缺陷对岩体力学特性及边坡稳定影响研究[J]. 岩土力学, 2019, 40(S1): 431-439.
[7] 金爱兵, 刘佳伟, 赵怡晴, 王本鑫, 孙浩, 魏余栋, . 卸荷条件下花岗岩力学特性分析[J]. 岩土力学, 2019, 40(S1): 459-467.
[8] 韩征, 粟滨, 李艳鸽, 王伟, 王卫东, 黄健陵, 陈光齐, . 基于HBP本构模型的泥石流动力过程SPH数值模拟[J]. 岩土力学, 2019, 40(S1): 477-485.
[9] 吴锦亮, 何吉, . 岩质边坡动态开挖模拟的复合单元模型[J]. 岩土力学, 2019, 40(S1): 535-540.
[10] 吴凤元, 樊赟赟, 陈剑平, 李军, . 基于不同侵蚀模型的高速崩滑碎屑 流动力过程模拟分析[J]. 岩土力学, 2019, 40(8): 3236-3246.
[11] 孙峰, 薛世峰, 逄铭玉, 唐梅荣, 张翔, 李川, . 基于连续损伤的水平井射孔-近井筒三维破裂模拟[J]. 岩土力学, 2019, 40(8): 3255-3261.
[12] 穆锐, 浦少云, 黄质宏, 李永辉, 郑培鑫, 刘 旸, 刘 泽, 郑红超, . 土岩组合岩体中抗拔桩极限承载力的确定[J]. 岩土力学, 2019, 40(7): 2825-2837.
[13] 金俊超, 佘成学, 尚朋阳. 基于应变软化指标的岩石非线性蠕变模型[J]. 岩土力学, 2019, 40(6): 2239-2246.
[14] 张 聪, 梁经纬, 阳军生, 曹 磊, 谢亦朋, 张贵金, . 堤坝脉动注浆浆液扩散机制及应用研究[J]. 岩土力学, 2019, 40(4): 1507-1514.
[15] 严 健, 何 川, 汪 波, 蒙 伟, . 高地温对隧道岩爆发生的影响性研究[J]. 岩土力学, 2019, 40(4): 1543-1550.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 黄建华,宋二祥. 大型锚碇基础围护工程冻结帷幕力学性态研究[J]. , 2009, 30(11): 3372 -3378 .
[2] 王观石,李长洪,陈保君,李世海. 应力波在非线性结构面介质中的传播规律[J]. , 2009, 30(12): 3747 -3752 .
[3] 王朝阳,许 强,倪万魁. 原状黄土CT试验中应力-应变关系的研究[J]. , 2010, 31(2): 387 -391 .
[4] 邓 琴,郭明伟,李春光,葛修润. 基于边界元法的边坡矢量和稳定分析[J]. , 2010, 31(6): 1971 -1976 .
[5] 万少石,年廷凯,蒋景彩,栾茂田. 边坡稳定强度折减有限元分析中的若干问题讨论[J]. , 2010, 31(7): 2283 -2288 .
[6] 闫 铁,李 玮,毕雪亮. 基于分形方法的多孔介质有效应力模型研究[J]. , 2010, 31(8): 2625 -2629 .
[7] 刘 嘉,王 栋. 正常固结黏土中平板锚基础的吸力和抗拉力[J]. , 2009, 30(3): 735 -740 .
[8] 赵尚毅,郑颖人,李安洪,邱文平,唐晓松,徐 俊. 多排埋入式抗滑桩在武隆县政府滑坡中的应用[J]. , 2009, 30(S1): 160 -164 .
[9] 刘振平,贺怀建,朱发华. 基于钻孔数据的三维可视化快速建模技术的研究[J]. , 2009, 30(S1): 260 -266 .
[10] 赵跃堂,林家炜,石 磊. 冲击荷载作用下层裂问题研究[J]. , 2011, 32(S2): 122 -126 .