Lithologic Reservoirs ›› 2026, Vol. 38 ›› Issue (3): 173-181.doi: 10.12108/yxyqc.20260315

• PETROLEUM ENGINEERING AND OIL & GAS FIELD DEVELOPMENT • Previous Articles     Next Articles

Pressure dynamic characteristics of CO2 flooding in dual porosity media reservoirs

CAO Xiutai1(), ZHONG Huiying1,2(), SUN Yuxin1, ZHOU Hongliang3, FU Jing4   

  1. 1 Key Laboratory of Enhanced Recovery of Ministry of Education, Northeast Petroleum University, Daqing 163318, Heilongjiang, China
    2 Key Laboratory of Reservoir Stimulation, China National Petroleum Corporation, Daqing 163318, Heilongjiang, China
    3 Yushulin Company Geological Technology Research InstituteDaqing Oilfield Co., Ltd., Daqing 163318, Heilongjiang, China
    4 Department of Petroleum Engineering, Colorado School of Mines, Golden 80401, Colorado, USA
  • Received:2025-09-23 Revised:2025-11-06 Online:2026-05-01 Published:2026-01-22

Abstract:

Shale oil reservoirs subjected to complex hydraulic fracturing are prone to evolve into dual porosity media. Due to differences in matrix-fracture seepage behavior and the areal heterogeneity of “concentration-viscosity” after CO2 injection, the pressure-response mechanism of CO2 flooding becomes highly complex, resulting in insufficient accuracy in pressure-transient characterization and parameter inversion. Fick’s law was employed to investigate the CO2 distribution characteristics. A dual porosity flow model was established by simultaneously accounting for concentration-viscosity-pressure coupling and the matrix threshold pressure gradient. The governing flow equations are solved numerically to generate pressure-transient curves for CO2 flooding in dual porosity reservoirs, and effects of the matrix threshold pressure gradient, elastic storativity ratio, interporosity-flow coefficient, injection rate, and diffusion coefficient on the pressure-transient behavior were analyzed. The results show that: (1) A larger matrix threshold pressure gradient leads to a more pronounced late-time upturn of both pressure and pressure-derivative curves on the well testing curve. (2) A smaller elastic storativity ratio produces a wider and deeper “trough” in the pressure-derivative curve during the interporosity-flow stage. (3) With decreasing interporosity-flow coefficient, matrix-to-fracture crossflow slows down,the “trough” in the pressure-derivative curve moves to the right and becomes deeper. A larger interporosity-flow coefficient makes the late-time upturn in the pressure curve more evident. (4) A higher injection rate increases both pressure and pressure-derivative levels over the entire testing period and narrows the “trough” during the interporosity-flow stage. (5) The diffusion coefficient mainly influences the radial-flow stage of the composite system in the late stage: a larger diffusion coefficient causes this stage to occur earlier, reduces flow resistance, and shifts the pressure and pressure-derivative curves downward. (6) The proposed model enables quantitative characterization of pressure-transient features during CO2 injection in fractured wells of Daqing Oilfield.

Key words: dual porosity media, CO2 flooding, pressure-transient curve, threshold pressure gradient, elastic storativity ratio, interporosity-flow coefficient, injection rate, diffusion coefficient, seepage mathematical model

CLC Number: 

  • TE353

Fig. 1

Physical model of CO2 flooding in a dual porosity media reservoir"

Fig. 2

Schematic of the configuration of fracture-pore dual porosity media reservoir"

Table 1

Dimensionless variables"

无因次变量 定义 无因次变量 定义
无因次渗透率 ${K}_{D}=\frac{{K}_{m}}{{K}_{f}}$ 无因次半径 ${r}_{D}=\frac{r}{{r}_{w}}$
无因次压力 ${{p}_{D}}_{{}_{j}}\left({r}_{D},{t}_{D}\right)=\frac{2\pi {K}_{f}h}{q{\mu }_{g}B}\left({p}_{0}-{p}_{j}\right)$
j=m,f
无因次启动压力梯度 ${G}_{D}=\frac{2\pi {K}_{f}h{r}_{w}}{q{\mu }_{g}B}G$
无因次井筒储集系数 ${C}_{D}=\frac{1}{2\pi {\varphi }_{f}{C}_{tf}h{r}_{w}^{2}}C$ 无因次黏度 ${\mu }_{D}=\frac{{\mu }_{0}}{{\mu }_{mix}}$

Fig. 3

Comparison of the pressure dynamic curves calculated by the established model and tNavigator software"

Table 2

Reservoir parameters from S151 block of Daqing Oilfield"

参数 数值 参数 数值
CO2注入速度/
(m3·s-1)
4.63×10-4 油藏厚度/m 5.0
扩散系数/(m2·s-1) 2.46×10-6 井筒半径/m 0.1
外边界半径/m 100 体积系数 1.12
基质渗透率/mD 0.500 裂缝渗透率/mD 1 000
启动压力梯度/
(Pa·m-1)
4×104 井筒存储系数/(m3·Pa-1) 1×10-8
表皮系数 0.1 原油初始黏度/(Pa·s) 0.35×10-3
CO2黏度/(Pa·s) 0.05×10-3 窜流系数 1×10-6
弹性储能比 0.01 综合压缩系数/
Pa-1
4×10-8

Fig. 4

CO2 concentration variation curves in dual porosity media reservoirs"

Fig. 5

Mixed viscosity change curves in dual porosity media reservoirs"

Fig. 6

Pressure dynamic curves of CO2 flooding in dual porosity media reservoirs"

Fig. 7

Pressure and pressure-derivative curves under different dimensionless threshold gradient"

Fig. 8

Pressure and pressure-derivative curves of w under different elastic storativity ratio"

Fig. 9

Pressure and pressure-derivative curves of λ under different interporosity-flow coefficient"

Fig. 10

Pressure and pressure-derivative curves of Q under different CO2 injection rate"

Fig. 11

Mixed viscosity curves of H under different diffusion coefficients with identical time (180 days)"

Fig. 12

Pressure and pressure-derivative curves of H under different diffusion coefficients"

Fig. 13

CO2 flooding well testing fitting curves of dual porosity media reservoirs in typical measured wells"

Table 3

Well testing interpretation results of well X11 and well X302 in S151 block and S11 block of Daqing Oilfield"

井号 井筒储集系数/(m3·Pa-1) 表皮
系数
基质
渗透率/mD
裂缝
渗透率/mD
窜流系数 弹性
储能比
扩散系数/
(m2·s-1)
启动压力梯度/(Pa·m-1)
X11井 2.91×10-8 0.15 0.132 1 120 2.314×10-5 0.112 2.46×10-6 2.4×104
X302井 2.48×10-8 0.20 0.793 1 460 1.056×10-6 0.009 2.12×10-6
[1] 郭平, 许清华, 孙振, 等. 天然气藏CO2驱及地质埋存技术研究进展[J]. 岩性油气藏, 2016, 28(3):6-11.
GUO Ping, XU Qinghua, SUN Zhen, et al. Research progress of CO2 flooding and geological storage in gas reservoirs[J]. Lithologic Reservoirs, 2016, 28(3):6-11.
[2] 宋新民, 王峰, 马德胜, 等. 中国石油二氧化碳捕集、驱油与埋存技术进展及展望[J]. 石油勘探与开发, 2023, 50(1):206-218.
SONG Xinmin, WANG Feng, MA Desheng, et al. Progress and prospect of carbon dioxide capture,utilization and storage in CNPC oilfields[J]. Petroleum Exploration and Development, 2023, 50(1):206-218.
[3] ZHANG Yiqi, YANG Shenglai, BI Lufei, et al. A technical review of CO2 flooding sweep-characteristics research advance and sweep-extend technology[J]. Petroleum Science, 2024, 22(1),225-276.
[4] 何希鹏, 张培先, 高玉巧, 等. 中国非常规油气资源效益开发面临的挑战与对策[J]. 中国石油勘探, 2025, 30(1):28-43.
HE Xipeng, ZHANG Peixian, GAO Yuqiao, et al. Challenges and countermeasures for beneficial development of unconventional oil and gas resources in China[J]. China Petroleum Exploration, 2025, 30(1):28-43.
[5] 贾承造. 论非常规油气对经典石油天然气地质学理论的突破及意义[J]. 石油勘探与开发, 2017, 44(1):1-11.
JIA Chengzao. Breakthrough and significance of unconventional oil and gas to classical petroleum geological theory[J]. Petroleum Exploration and Development, 2017, 44(1):1-11.
[6] 张衍君, 刘拯君, 徐豪, 等. 页岩油储层前置CO2压裂液体滞留效应研究进展[J]. 岩性油气藏, 2026, 38(1):180-190.
ZHANG Yanjun, LIU Zhengjun, XU Hao, et al. Research pro-gress on retention effects of pre-CO2 fracturing fluid of shale oil reservoirs[J]. Lithologic Reservoirs, 2026, 38(1):180-190.
[7] 徐田录, 吴承美, 张金凤, 等. 吉木萨尔凹陷二叠系芦草沟组页岩油储层天然裂缝特征与压裂模拟[J]. 岩性油气藏, 2024, 36(4):35-43.
XU Tianlu, WU Chengmei, ZHANG Jinfeng, et al. Natural fracture characteristics and fracture network simulation in shale reservoirs of Permian Lucaogou Formation in Jimsar Sag[J]. Lithologic Reservoirs, 2024, 36(4):35-43.
[8] TANG R W, AMBASTHA A K. Analysis of CO2 pressure transient data with two- and three-region analytical radial compo-site models[R]. Houston,SPE Annual Technical Conference and Exhibition, 1988.
[9] 姜瑞忠, 张海涛, 张伟, 等. CO2驱三区复合油藏水平井压力动态分析[J]. 油气地质与采收率, 2018, 25(6):63-70.
JIANG Ruizhong, ZHANG Haitao, ZHANG Wei, et al. Dynamic pressure analysis of three-zone composite horizontal well in oil reservoirs for CO2 flooding[J]. Petroleum Geology and Reco-very Efficiency, 2018, 25(6):63-70.
[10] SU Kun, LIAO Xinwei, ZHAO Xiaoliang. Transient pressure analysis and interpretation for analytical composite model of CO2 flooding[J]. Journal of Petroleum Science and Enginee-ring, 2015, 125:128-135.
[11] 吴明录, 李涛, 赵高龙, 等. 双重孔隙介质三区复合油藏水平井试井模型[J]. 油气井测试, 2022, 31(4):6-12.
WU Minglu, LI Tao, ZHAO Gaolong, et al. Well test model for horizontal well in dual-porosity three-zone composite reservoir[J]. Well Testing, 2022, 31(4):6-12.
[12] 吴明录, 赵高龙, 姚军, 等. 稠油热采三区复合油藏水平井试井解释模型及压力响应特征[J]. 大庆石油地质与开发, 2016, 35(6):117-122.
WU Minglu, ZHAO Gaolong, YAO Jun, et al. Well test interpreting model and pressure response characteristics of the horizontal well in thermally recovered three-block composite heavy oil reservoirs[J]. Petroleum Geology & Oilfield Development in Daqing, 2016, 35(6):117-122.
[13] 魏明强, 段永刚, 方全堂, 等. 页岩气藏多级压裂水平井压力动态分析[J]. 中南大学学报(自然科学版), 2016, 47(12):4141-4147.
WEI Mingqiang, DUAN Yonggang, FANG Quantang, et al. Pressure dynamic analysis of multistage fractured horizontal well in shale gas reservoirs[J]. Journal of Central South University (Science and Technology), 2016, 47(12):4141-4147.
[14] 卢婷, 王鸣川, 马文礼, 等. 考虑多重应力敏感效应的页岩气藏压裂水平井试井模型[J]. 新疆石油地质, 2021, 42(6):741-748.
LU Ting, WANG Mingchuan, MA Wenli, et al. Fractured horizontal well test model for shale gas reservoirs with considering multiple stress sensitive factors[J]. Xinjiang Petroleum Geo-logy, 2021, 42(6):741-748.
[15] 陈志明, YU Wei, 石璐铭, 等. 压裂水平井的多模式裂缝网络试井模型及参数评价:以吉木萨尔页岩油为例[J]. 石油学报, 2023, 44(10):1706-1726.
CHEN Zhiming, YU Wei, SHI Luming, et al. Well test model and parameter evaluation of multi-mode fracture network in fractured horizontal well:A case study of Jimsar shale oil[J]. Acta Petrolei Sinica, 2023, 44(10):1706-1726.
[16] 徐阳. 低渗油藏CO2近混相驱提高采收率机理研究[D]. 青岛: 中国石油大学(华东), 2011.
XU Yang. Mechanisms of near-miscible CO2 flooding for EOR in low permeability reservoirs[D]. Qingdao: China University of Petroleum (East China), 2011.
[17] ZHANG Na, YIN Mingfei, WEI Mingzhen, et al. Identification of CO2 sequestration opportunities:CO2 miscible flooding guidelines[J]. Fuel, 2019, 241:459-467.
[18] GONG Rundong, LIU Junrong, LIU Shuyang, et al. Impact of interlayer mass transfer on CO2 flooding performance in heterogeneous tight reservoirs[J]. Physics of Fluids, 2025, 37(2):023134.
[19] 侯艳红. 低渗透油藏CO2混相驱油机理研究[D]. 青岛: 中国石油大学(华东), 2010.
HOU Yanhong. Mechanism study of CO2 miscible displacement in low permeability reservoirs[D]. Qingdao: China University of Petroleum(East China), 2010.
[20] 苏玉亮, 吴晓东, 侯艳红, 等. 低渗透油藏CO2混相驱油机制及影响因素[J]. 中国石油大学学报(自然科学版), 2011, 35(3):99-102.
SU Yuliang, WU Xiaodong, HOU Yanhong, et al. Mechanism of CO2 miscible displacement in low permeability reservoir and influencing factor[J]. Journal of China University of Petroleum (Edition of Natural Science), 2011, 35(3):99-102.
[21] 崔传智, 李静, 吴忠维. 扩散吸附作用下CO2非混相驱微观渗流特征模拟[J]. 岩性油气藏, 2024, 36(6) :181-188.
CUI Chuanzhi, LI Jing, WU Zhongwei. Simulation of microscopic seepage characteristics of CO2 immiscible flooding under the effect of diffusion and adsorption[J]. Lithologic Reservoirs, 2024, 36(6):181-188.
[22] 李友全, 韩秀虹, 阎燕, 等. 低渗透油藏CO2吞吐压力响应曲线分析[J]. 岩性油气藏, 2017, 29(6):119-127.
LI Youquan, HAN Xiuhong, YAN Yan, et al. Pressure transient analysis curves on CO2 huff and puff in low permeability reservoir[J]. Lithologic Reservoirs, 2017, 29(6):119-127.
[23] 苏玉亮, 吴春新. 倾斜油藏二氧化碳驱油中的黏性指进分析[J]. 石油钻探技术, 2009, 37(5):93-96.
SU Yuliang, WU Chunxin. Analysis of viscous fingering of CO2 flooding in tilted reservoirs[J]. Petroleum Drilling Techniques, 2009, 37(5):93-96.
[24] 钟会影, 余承挚, 沈文霞, 等. 考虑启动压力梯度的致密油藏水平井裂缝干扰渗流特征[J]. 岩性油气藏, 2024, 36(3):172-179.
ZHONG Huiying, YU Chengzhi, SHEN Wenxia, et al. Characteristics of fracture interference between horizontal wells in tight reservoirs considering threshold pressure gradient[J]. Lithologic Reservoirs, 2024, 36(3):172-179.
[25] 何吉祥, 徐有杰, 高阳, 等. 裂缝性致密油藏多级压裂水平井试井模型[J]. 断块油气田, 2021, 28(2):241-246.
HE Jixiang, XU Youjie, GAO Yang, et al. Well test model of multi-stage fractured horizontal well in fractured tight reservoirs[J]. Fault-Block Oil & Gas Field, 2021, 28(2):241-246.
[1] XIONG Yu, ZHANG Weicen, LI Yamei, GENG Wenshuang, WU Daoming, MU Dan, LIU Tong. Prediction method of water-to-gas gravity displacement efficiency and interface stability of high-dip angle reservoirs:A case study of Paleogene Shahejie Formation reservoir in Liuzan area of Bohai Bay Basin [J]. Lithologic Reservoirs, 2026, 38(2): 178-193.
[2] LI Jiamin, ZHANG Yizhong, ZHANG Maolin, QIN Bowen, YANG Yuxin. Minimum miscible pressure model of CO2 flooding based on kernel ridgeregression multi-parameter optimization [J]. Lithologic Reservoirs, 2025, 37(6): 172-179.
[3] XU Youjie, REN Zongxiao, XIANG Zuping, FAN Xiaohui, YU Mengnan. Numerical well testing model of fractured well with complex fractures multi-well interference in heterogeneous tight gas reservoirs [J]. Lithologic Reservoirs, 2025, 37(3): 194-200.
[4] ZHONG Huiying, YU Chengzhi, SHEN Wenxia, BI Yongbin, YI Ran, NI Haoming. Characteristics of fracture interference between horizontal wells in tight reservoirs considering threshold pressure gradient [J]. Lithologic Reservoirs, 2024, 36(3): 172-179.
[5] GUO Yongwei, YAN Fangping, WANG Jing, CHU Huili, YANG Jianlei, CHEN Yingchao, ZHANG Xiaoyang. Characteristics of solid deposition and reservoir damage of CO2 flooding in tight sandstone reservoirs [J]. Lithologic Reservoirs, 2021, 33(3): 153-161.
[6] SUN Huizhu, ZHU Yushuang, WEI Yong, GAO Yuan. Influence mechanism of acidification on oil recovery during CO2 flooding [J]. Lithologic Reservoirs, 2020, 32(4): 136-142.
[7] CUI Yongzheng, JIANG Ruizhong, GAO Yihua, QIAO Xin, WANG Qiong. Pressure transient analysis of hydraulic fractured vertical wells with variable conductivity for CO2 flooding [J]. Lithologic Reservoirs, 2020, 32(4): 172-180.
[8] TANG Meirong, ZHANG Tongwu, BAI Xiaohu, WANG Xuanyi, LI Chuan. Influence of pore throat structure on reservoir damage with CO2 flooding [J]. Lithologic Reservoirs, 2019, 31(3): 113-119.
[9] JIANG Ruizhong, SHEN Zeyang, CUI Yongzheng, ZHANG Fulei, ZHANG Chunguang, YUAN Jianwei. Dynamical characteristics of inclined well in dual medium low permeability reservoir [J]. Lithologic Reservoirs, 2018, 30(6): 131-137.
[10] SHANG Qinghua, WANG Yuxia, HUANG Chunxia, CHEN Longlong. Supercritical and non-supercritical CO2 flooding characteristics in tight sandstone reservoir [J]. Lithologic Reservoirs, 2018, 30(3): 153-158.
[11] MA Li, OUYANG Chuanxiang, TAN Zhengyang, WANG Changquan, SONG Yan, LIN Fei. Efficiency improvement of CO2 flooding in middle and later stage for low permeability reservoirs [J]. Lithologic Reservoirs, 2018, 30(2): 139-145.
[12] YANG Hong, WANG Hong, NAN Yufeng, QU Yaning, LIANG Kaiqiang, JIANG Shaojing. Suitability evaluation of enhanced oil recovery by CO2 flooding [J]. Lithologic Reservoirs, 2017, 29(3): 140-146.
[13] DENG Xuefeng. Optimization of reasonable production pressure difference of fractured horizontal well in low permeability tight reservoirs [J]. Lithologic Reservoirs, 2017, 29(1): 135-139.
[14] Li Youquan,Meng Fankun,Yan Yan,Han Fengrui,Yu Weijie,Zhou Shiyu. Pressure transient analysis on CO2 flooding in low permeability reservoirs considering fluid heterogeneity [J]. Lithologic Reservoirs, 2016, 28(4): 106-112.
[15] Guo Ping,Xu Qinghua,Sun Zhen,Du Jianfen,Wang Zhouhua. Research progress of CO2 flooding and geological storage in gas reservoirs [J]. Lithologic Reservoirs, 2016, 28(3): 6-11.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!
TRENDMD: