考虑温度效应的裂缝型碳酸盐岩酸化数值模拟

杨俊超; 程远方; 韩忠英; 韩松财; 闫传梁; 赵珠宇; 孙波

doi:10.3969/j.issn.1006-6535.2025.04.018

特种油气藏 >

2025 , Vol. 32 >Issue 4: 149 - 157

DOI: https://doi.org/10.3969/j.issn.1006-6535.2025.04.018

钻采工程

考虑温度效应的裂缝型碳酸盐岩酸化数值模拟

杨俊超 ,
程远方 ,
韩忠英 ,
韩松财 ,
闫传梁 ,
赵珠宇 ,
孙波

展开

1.中国石油大学(华东),山东青岛 266580;
2.中国石油大学(华东)非常规油气开发教育部重点实验室,山东青岛 266580;
3.成都理工大学,四川成都 610059

杨俊超(1996—),男,2019年毕业于中国石油大学(华东)石油工程专业,现为该校油气井工程专业在读博士研究生,主要从事油气井酸化压裂方面的研究。

收稿日期: 2024-03-11

修回日期: 2025-05-19

网络出版日期: 2025-09-03

基金资助

国家自然科学基金“超深特深层油气钻采流动调控”(52288101);四川省自然科学基金“液氮压裂页岩储层多重致裂机制耦合分析及复杂缝网动态演化模拟研究”(2022NSFSC1227);山东省高等学校“青创团队计划”“天然气水合物多场耦合创新团队”(2022KJ072)

收起

Numerical simulation of acidizing in fractured carbonate rocks considering temperature effects

YANG Junchao ,
CHENG Yuanfang ,
HAN Zhongying ,
HAN Songcai ,
YAN Chuanliang ,
ZHAO Zhuyu ,
SUN Bo

Expand

1. China University of Petroleum (East China), Qingdao, Shandong 266580, China;
2. Key Laboratory of Unconventional Oil and Gas Development, Ministry of Education, China University of Petroleum (East China), Qingdao, Shandong 266580, China;
3. Chengdu University of Technology, Chengdu, Sichuan 610059, China

Received date: 2024-03-11

Revised date: 2025-05-19

Online published: 2025-09-03

Fold

摘要

针对深部裂缝型碳酸盐岩储层酸化蚓孔扩展机理不清的问题,基于双尺度连续模型和裂缝离散处理方法,建立了渗流场-化学场-温度场耦合多矿物成分酸化的MFCAC数学模型,分析了考虑温度效应的裂缝型储层蚓孔扩展机理和酸蚀特征,研究了CaMg(CO₃)₂与CaCO₃初始含量比、酸液温度和反应放热对蚓孔扩展的影响。研究表明:主蚓孔扩展方向主要受周围距离较近的裂缝控制域主导;温度效应导致酸岩反应速率增加,H⁺有效传质距离缩短,导致主蚓孔沟通的天然裂缝数量减少;反应放热导致的H⁺扩散速率增加使酸液更容易向各方向孔隙扩散,蚓孔形态更复杂;酸液突破体积比随CaMg(CO₃)₂初始含量的增加先减小后增大;CaMg(CO₃)₂与CaCO₃初始含量比为1.5时,对应的酸液突破体积比最小;酸液突破体积比随酸液温度的增加而变大;与不考虑反应放热相比,考虑反应放热的酸液突破体积比更小。该研究为深部裂缝型碳酸盐岩储层酸化施工中动态蚓孔扩展的精准调控及施工参数优化提供了理论依据。

关键词： 蚓孔扩展; 含量比; 酸液突破体积比; 酸液温度; 反应放热

本文引用格式

杨俊超 , 程远方 , 韩忠英 , 韩松财 , 闫传梁 , 赵珠宇 , 孙波 . 考虑温度效应的裂缝型碳酸盐岩酸化数值模拟[J]. 特种油气藏, 2025 , 32(4) : 149 -157 . DOI: 10.3969/j.issn.1006-6535.2025.04.018

Abstract

To address the unclear mechanism of wormhole propagation during acidizing in deep fractured carbonate rocks,an MFCAC mathematical model coupling seepage field,chemical field,and temperature field for multi-mineral component acidizing was established based on a dual-scale continuum model and discrete fracture treatment method.The mechanism of wormhole propagation and acid-etching characteristics in fractured reservoirs considering temperature effects were analyzed,and the influences of the initial content ratio of CaMg(CO₃)₂ to CaCO₃,acid temperature,and reaction exotherm on wormhole propagation were studied.The study shows that the propagation direction of the main wormhole is mainly dominated by the control domain of nearby fractures;temperature effects increase the acid-rock reaction rate but decrease the H⁺ mass transfer rate,resulting in fewer natural fractures connected by the main wormhole;the increased H⁺ diffusion rate caused by temperature effects makes acid easier to diffuse into pores in all directions,leading to more complex wormhole morphology;the pore volume to breakthrough (PVBT) first decreases and then increases with increasing initial CaMg(CO₃)₂ content;when the initial content ratio of CaMg(CO₃)₂ to CaCO₃ is 1.5,the corresponding PVBT is the smallest;PVBT increases with increasing acid temperature;compared to ignoring reaction exotherm,considering reaction exotherm results in a smaller PVBT.This research provides a theoretical basis for the precise control of dynamic wormhole propagation and optimization of operational parameters during acidizing operations in deep fractured carbonate reservoirs.

Key words： wormhole propagation; content ratio; pore volume to breakthrough (PVBT); acid temperature; reaction exotherm

参考文献

[1] 马乃拜,马新平,杜伟维,等.顺北一区奥陶系碳酸盐岩储层发育特征与预测[J].内蒙古石油化工,2019,45(7):106-112.
MA Naibai,MA Xinping,DU Weiwei,et al.Developmental characteristics for the carbonate reservoir in Ordovician in first area of Shunbei[J].Inner Mongolia Petrochemical Industry,2019,45(7):106-112.
[2] 姚军,黄朝琴,王子胜,等.缝洞型油藏的离散缝洞网络流动数学模型[J].石油学报,2010,31(5):815-809.
YAO Jun,HUANG Chaoqin,WANG Zisheng,et al.Mathematical model of fluid flow in fractured vuggy reservoirs based on discrete fracture-vug network[J].Acta Petrolei Sinica,2010,31(5):815-809.
[3] 康志江.缝洞型碳酸盐岩油藏耦合数值模拟新方法[J].新疆石油地质,2010,31(5):514-516.
KANG Zhijiang.New method of coupling numerical simulation and application to fracture-cavern carbonate reservoir[J].Xinjiang Petroleum Geology,2010,31(5):514-516.
[4] 黄朝琴,姚军,王月英,等.基于离散裂缝模型的裂缝性油藏注水开发数值模拟[J].计算物理,2011,28(1):41-49.
HUANG Chaoqin,YAO Jun,WANG Yueying,et al.Numerical simulation on water flooding development of fractured reservoirs in a discrete-fracture model[J].Chinese Journal of Computational Physics,2011,28(1):41-49.
[5] 高艳霞,李小波,彭小龙,等.缝洞型油藏大尺度缝洞体等效模拟方法研究[J].长江大学学报(自然科学版),2016,13(14):66-69.
GAO Yanxia,LI Xiaobo,PENG Xiaolong,et al.The equivalent simulation of big size cavity-fracture in cavity-fractured reservoirs[J].Journal of Yangtze University(Natural Science Edition),2016,13(14):66-69.
[6] 赵立强,王润宇,刘平礼,等.考虑天然微裂缝的酸蚀蚓孔扩展数值模拟研究[J].油气藏评价与开发,2020,10(2):76-82.
ZHAO Liqiang,WANG Runyu,LIU Pingli,et al.Numerical simulation of wormhole propagation considering natural micro-fractures[J].Reservoir Evaluation and Development,2020,10(2):76-82.
[7] UPADHYAY V K,SZYMCZAK P,LADD A J C.Initial conditions or emergence:what determines dissolution patterns in rough fractures?[J].Journal of Geophysical Research:Solid Earth,2015,120(9):6102-6121.
[8] HANNA R B,RAJARAM H.Influence of aperture variability on dissolutional growth of fissures in karst formations[J].Water Resources Research,1998,34(11):2843-2853.
[9] HILL A D,ZHU D,DONG C,et al.Computer model predicts matrix acidizing of naturally fractured carbonate[J].Journal of Petroleum Technology,2001,53(10):20-25.
[10] DONG C,ZHU D,HILL A D.Acid penetration in natural fracture networks[J].SPE Production & Facilities,2002,17(3):160-170.
[11] UGURSAL A,ZHU D,HILL A D.Development of acid fracturing model for naturally fractured reservoirs[J].SPE Production & Operations,2019,34(4):735-748.
[12] 李勇明,廖毅,赵金洲,等.基于双尺度等效渗流模型的复杂碳酸盐岩蚓孔扩展形态研究[J].天然气地球科学,2016,27(1):121-127.
LI Yongming,LIAO Yi,ZHAO Jinzhou,et al.Wormhole dissolution pattern study in complicated carbonate rock based on two-scale continuum model and equivalent seepage theory[J].Natural Gas Geoscience,2016,27(1):121-127.
[13] KALIA N,BALAKOTAIAH V.Effect of medium heterogeneities on reactive dissolution of carbonates[J].Chemical Engineering Science,2009,64(2):376-390.
[14] GAROPOV T T,KARIMI-FARD M,TCHELEPI H A.Discrete fracture model for coupled flow and geomechanics[J].Computational Geosciences,2016,20(1):149-160.
[15] SANDVE T H,BERRE I,NORDBOTTEN J M.An efficient multi-point flux approximation method for discrete fracture-matrix simulations[J].Journal of Computational Physics,2012,231(9):3784-3800.
[16] THERRIEN R,SUDICKY E A.Three-dimensional analysis of variably-saturated flow and solute transport in discretely-fractured porous media[J].Journal of Contaminant Hydrology,1996,23(1/2):1-12.
[17] PANG A M K R,ZIAUDDIN M,BALAKOTAIAH V.Two-scale contineem model for simulation of wormholes in carbonate acidization[J].American Institute of Chemical Engineers Journal,2005,51(12):3231-3248.
[18] 华青,王昱珩,张娜,等.碳酸盐岩酸蚀蚓孔影响因素模拟分析研究[J].石油地质与工程,2023,37(6):97-102.
HUA Qing,WANG Yuheng,ZHANG Na,et al.Simulation analysis of influencing factors of carbonate acidizing wormholes[J].Petroleum Geology & Engineering,2023,37(6):97-102.
[19] 牟建业,李双明,赵鑫,等.基于真实孔隙空间分布的酸蚀蚓孔扩展规律数值模拟研究[J].科学技术与工程,2014,14(35):40-46.
MU Jianye,LI Shuangming,ZHAO Xin,et al.Modeling wormhole propagation behavior based on real pore spatial distributions[J].Science Technology and Engineering,2014,14(35):40-46.
[20] LIU P,YAO J,COUPLES G D,et al.Numerical modelling and analysis of reactive flow and wormhole formation in fractured carbonate rocks[J].Chemical Engineering Science,2017,67(172):143-157.
[21] 齐宁,陈国彬,李振亮,等.基于分步算法的裂缝性碳酸盐岩油藏大尺度酸化数值模拟[J].石油学报,2020,41(3):348-362.
QI Ning,CHEN Guobin,LI Zhenliang,et al.Numerical simulation of large-scale acidification in fractured carbonate reservoirs based on a step-by-step algorithm[J].Acta Petrolei Sinica,2020,41(3):348-362.
[22] 潘林.白云岩储层酸压非均匀刻蚀机理研究[D].青岛:中国石油大学(华东),2021.
PAN Lin.Mechanism of heterogeneous etching in acid fracturing of dolomite reservoirs[D].Qingdao:China University of Petroleum(East China),2021.
[23] 齐宁,李柏杨,方明君,等.基于碳酸盐岩酸化溶蚀形态的酸液最优注入速度界限[J].中国石油大学学报(自然科学版),2017,41(5):117-122.
QI Ning,LI Baiyang,FANG Mingjun,et al.Injection rate optimization for acidizing process of carbonate rocks based on dissolution morphology[J].Journal of China University of Petroleum(Edition of Natural Sciences),2017,41(5):117-122.
[24] FREDD C N,FOGLER H S.Alternative stimulation fluids and their impact on carbonate acidizing[J].SPE Journal,1998,3(1):34-41.
[25] FREDD C N,FOGLER H S.Influence of transport and reaction on wormhole formation in porous media[J].AIChE Journal,1998,44(9):1933-1949.
[26] FREDD C N,FOGLER H S.Optimum conditions for wormhole formation in carbonate porous media:influence of transport and reaction[J].SPE Journal,1999,4(3):196-205.
[27] QI Ning,CHEN Guobin,LIANG Chong,et al.Numerical simulation and analysis of the influence of fracture geometry on wormhole propagation in carbonate reservoirs[J].Chemical Engineering Science,2019,198(12):124-143.
[28] DONG Kai,ZHU Ding,HILL A D.The role of temperature on optimal conditions in dolomite acidizing:an experimental study and its applications[J].Journal of Petroleum Science and Engineering,2018,165:736-742.
[29] BAZIN B.From matrix acidizing to acid fracturing:a laboratory evaluation of acid/rock interactions[J].SPE Production & Facilities,2001,16(1):22-29.

Options

文章导航

模态框（Modal）标题

摘要

本文引用格式

Abstract

参考文献