Pore-scale Competitive Migration and Selective Regulation of Multi-component Hydrocarbons in Continental Shale Oil: A Phase-field Modeling Study
DOI:
https://doi.org/10.15377/2409-787X.2026.13.3Keywords:
Phase-field model, Continental shale oil, Competitive seepage, Pore-scale simulation, Multi-component hydrocarbon migrationAbstract
Continental shale oil is a vital alternative resource for oil and gas reserve and production growth in China. However, strong reservoir heterogeneity, widely developed micro-nano pore throats and complex multi-component crude oil properties lead to low crude oil producing degree and restrict efficient large-scale development. Most existing seepage studies are based on single-component assumptions, failing to fully reveal multi-component differential migration, competitive seepage mechanisms and the threshold regulation effect of pore scale. To address this issue, this study builds a pore-scale multi-component seepage model through numerical simulation, divides hydrocarbons into light, medium and heavy components, and systematically analyzes migration characteristics under varying injection flow rates and particle radii. The results show that component properties dominate differential migration and competitive seepage: light components with strong convection and diffusion capacities preferentially occupy dominant seepage channels and gradually compress the flow space of medium and heavy components. Within the range of 1.0-3.0 g/s, injection flow rate has a weak influence on the multi-component distribution pattern and does not alter the relative migration capacity of each component. Pore scale imposes selective regulation with a distinct threshold effect: heavier components are more sensitive to pore tightening, and the inter-component migration capacity gap widens sharply when the throat width shrinks below a critical value of δcr ∈ (5, 10) mm (corresponding to a critical particle radius rcr ∈ (20.0, 22.5) mm and a constriction ratio αcr ∈ (0.10, 0.20)). This study deepens the understanding of multi-component microscopic migration mechanisms of continental shale oil, and provides theoretical support for reservoir producing degree evaluation, development parameter optimization and hierarchical reservoir stimulation
References
[1] Zhao WZ, Hu SY, Hou LH, Yang T, Li X, Guo BC, et al. Types and resource potential of continental shale oil in China and its boundary with tight oil. Pet Explor Dev. 2020; 47(1): 1‑11. https://doi.org/10.1016/S1876-3804(20)60001-5
[2] Hu SY, Bai B, Tao SZ, Bian CS, Zhang TS, Chen YY, et al. Heterogeneous geological conditions and differential enrichment of medium and high maturity continental shale oil in China. Pet Explor Dev. 2022; 49(2): 257‑271. https://doi.org/10.1016/S1876-3804(22)60022-3
[3] Li H, Xu ZQ, Gao XD, Xie JT, Qin QR, Wang HJ, et al. Pore structure evolution and geological controls in lacustrine shale systems with implications for marine shale reservoir characterization. Sci Rep. 2025; 15: 17702. https://doi.org/10.1038/s41598-025-02415-y
[4] Zhao W, Jia CZ, Song Y, Li XF, Hou LH, Jiang L, et al. Hydrocarbon migration and accumulation simulation: A review and a novel multi‑scale quantitative numerical simulation method. Adv Colloid Interface Sci. 2025; 342: 103523. https://doi.org/10.1016/j.cis.2025.103523
[5] Xu F, Li Y, Xu X, Yang J, Xiao F, Shi Y. Thermal‑hydrological‑mechanical coupling simulation of matrix‑type shale oil: Analysis of oil‑gas‑water multiphase flow mechanism under various parameter conditions. Geoenergy Sci Eng. 2025; 246: 213656. https://doi.org/10.1016/j.geoen.2025.213656
[6] Wang YM, Wang ZL, Zhang ZC, Yao SS, Zhang H, Zheng GQ, et al. Recent techniques on analyses and characterizations of shale gas and oil reservoir. Energy Rev. 2024; 3(2): 100067. https://doi.org/10.1016/j.enrev.2023.100067
[7] Will J, Eckardt S. Optimization of Hydrocarbon Production from Unconventional Shale Reservoirs using Numerical Modelling. Int J Pet Technol. 2017; 4: 14‑23. https://doi.org/10.15377/2409-787X.2017.04.01.3
[8] Yu CL, Zhang SM, Chen T, Sun ZG, Xu JX. Seepage mode in lamina‑developed shale oil reservoirs under strong heterogeneous and strong fluid-solid coupling in Jiyang depression of China. Sci Rep. 2025; 15: 936. https://doi.org/10.1038/s41598-024-84882-3
[9] Li JH, Li YJ, Sang Q, Gong HJ, Xu L, Zhang HY, et al. Quantification of seepage characteristics in shale oil reservoirs: A triple medium model‑driven approach. Geoenergy Sci Eng. 2024; 242: 213225. https://doi.org/10.1016/j.geoen.2024.213225
[10] Xu S, Wen J, Gou QY, Li QQ, Liu BC. Research Progress and Significance of Shale Oil Micro‑Migration. J Earth Sci. 2024; 35: 1765‑9. https://doi.org/10.1007/s12583-024-0071-7
[11] Alpak FO, Riviere B, Frank F. A phase‑field method for the direct simulation of two‑phase flows in pore‑scale media using a non‑equilibrium wetting boundary condition. Comput Geosci. 2016; 20: 881‑908. https://doi.org/10.1007/s10596-015-9551-2
[12] Cai JC, Qin XJ, Xia XZ, Jiao XH, Chen H, Wang H, et al. Numerical modeling of multiphase flow in porous media considering micro‑ and nanoscale effects: A comprehensive review. Gas Sci Eng. 2024; 131: 205441. https://doi.org/10.1016/j.jgsce.2024.205441
[13] Ramstad T, Berg CF, Thompson K. Pore‑Scale Simulations of Single‑ and Two‑Phase Flow in Porous Media: Approaches and Applications. Transp Porous Med. 2019; 130: 77‑104. https://doi.org/10.1007/s11242-019-01289-9
[14] Li YH, Xu FQ, Xu XY, Bi CQ, Shi Y, Han RB. Investigation on differentiation behaviors of multiple hydrocarbon components in intra‑source migration within matrix‑type shale oil reservoirs. Fuel 2027; 429(Pt E): 141063. https://doi.org/10.1016/j.fuel.2026.141063
[15] Guan Z, Meng Q, Huang Q, Tang X, Shan Y, Liu G, et al. Hydrocarbon micro‑migration and differential enrichment mechanism of shale oil in upper Cretaceous Qingshankou Formation in Songliao Basin. Geol J. 2024; 59(5): 1503‑19. https://doi.org/10.1002/gj.4952
[16] Anonymous. Molecular dynamics study of shale oil adsorption and diffusion behavior in reservoir nanopores: Impact of hydrocarbon composition and surface type. J Mol Liq. 2024; 414: 126110. https://doi.org/10.1016/j.molliq.2024.126110
[17] Xue CL, Ji DL, Cheng D, Wen YT, Luo HH, Li Y. Adsorption behaviors of different components of shale oil in quartz slits studied by molecular simulation. ACS Omega. 2022; 7(45): 41189‑200. https://doi.org/10.1021/acsomega.2c04845
[18] Xiao L, Ju YW, Jin J, Liu J, Tan JQ, Jiang S, et al. Micro‑nanoscale pore structure characteristics and enrichment mechanisms of shale oil in mixed sedimentary reservoirs of the Jimusar Sag, Junggar Basin, northwestern China. Pet Sci. 2026; 23(3): 1089‑104. https://doi.org/10.1016/j.petsci.2025.10.014
[19] Xia DB, Yang ZM, Gao TN, Li HB, Lin W. Characteristics of micro‑ and nano‑pores in shale oil reservoirs. J Pet Explor Prod Technol. 2021; 11: 157‑69. https://doi.org/10.1007/s13202-020-01012-1
[20] Li C, Liu XP, Meng LJ, Wang JW, Gao YL, Sun B, et al. Fractal heterogeneity characteristics and influenced factors of lacustrine shale pore structure based on N₂ adsorption and mercury injection in the Paleogene Shahejie Formation, Bohai Bay Basin, China. Pet Res. 2026; 9(1): 1‑12. https://doi.org/10.1016/j.ptlrs.2026.01.019
[21] Chen ZM, Zhou B, Olorode O, Yu W, Sepehrnoori K. Evaluation of CO₂ sequestration potential in shale reservoirs using an improved embedded discrete fracture model. Int J Rock Mech Min Sci. 2026; 204: 106583. https://doi.org/10.1016/j.ijrmms.2026.106583
[22] Yan GG, Yan WD, Yuan YZ, Gong XJ, Tang ZQ, Bai XY. Research on the method of evaluating the pore throat structure of rock microscopically based on the 3d pore network model of digital core. Int J Pet Technol. 2022; 44: 44‑53. https://doi.org/10.54653/2409-787X.2022.09.6
[23] Qu YQ, Ouyang SQ, Gao JW, Shi J, Wu YY, Cheng YT, et al. Pore Space Characteristics and Migration Changes in Hydrocarbons in Shale Reservoir. Fractal Fract. 2024; 8(10): 588. https://doi.org/10.3390/fractalfract8100588
[24] Wang WD, Zhang Q, Xu JL, Zheng D, Liu LF, Su YL. Hydrodynamic resistance of pore-throat structures and its effect on shale oil apparent permeability. Int J Coal Sci Technol. 2024; 11: 22. https://doi.org/10.1007/s40789-024-00671-3
[25] Li JL, Ji LL, Jiang WB, Lin M, Mao M, Wei XL, et al. Investigation of pore structures and flow transport characteristics of continental shale in Dongying Sag. J Pet Explor Prod Technol. 2025; 15: 163. https://doi.org/10.1007/s13202-025-02079-4
[26] Xu FQ, Song XZ, Song GF, Ji JY, Song ZH, Shi Y, et al. Numerical studies on heat extraction evaluation and multi‑objective optimization of abandoned oil well patterns in intermittent operation mode. Energy 2023; 269: 126777. https://doi.org/10.1016/j.energy.2023.126777
[27] Gao XY, Yang SL, Zhang YQ, He BY, Wang MY, Hu JT, et al. Integrated modeling and economic assessment of hydrogen storage with wellbore‑reservoir‑thermo‑hydro‑mechanical‑diffusion coupling. Energy 2026; 353: 141065. https://doi.org/10.1016/j.energy.2026.141065
[28] Pati S, Borah A, Boruah MP, Randive PR. Critical review on local thermal equilibrium and local thermal non‑equilibrium approaches for the analysis of forced convective flow through porous media. Int Commun Heat Mass Transf. 2022; 132: 105889. https://doi.org/10.1016/j.icheatmasstransfer.2022.105889
[29] Hamidi S, Heinze T, Galvan B, Miller S. Critical review of the local thermal equilibrium assumption in heterogeneous porous media: Dependence on permeability and porosity contrasts. Appl Therm Eng. 2019; 147: 962‑71. https://doi.org/10.1016/j.applthermaleng.2018.10.130
[30] Xu FQ, Song XZ, Shi Y, Zhou MM, Wang GS. Erosion damage mechanism of fractured rock mass during cyclic compressed air storage and release in depleted sandstone oil‑gas fields. Geoenergy Sci Eng. 2026; 260: 214389. https://doi.org/10.1016/j.geoen.2026.214389
[31] Wang YX, Li HY, Wang ZZ, Xu FQ, Wang TR, Wang C, et al. Particle tracing‑based numerical investigation of thermal‑hydraulic coupling mechanisms and particle transport characteristics in multi‑branch well underground thermal energy storage (UTES). Int J Therm Sci. 2026; 230: 111175. https://doi.org/10.1016/j.ijthermalsci.2026.111175
[32] Shi Y, Cui QL, Song XZ, Li GS, Wang GS. Performances of a novel compressed CO₂ energy storage and heat storage integration system using multi‑stage hydraulic fractures of horizontal well in subsurface reservoirs. Renew Energy 2025; 251: 123396. https://doi.org/10.1016/j.renene.2025.123396
[33] Li S, Wang GS, Zhou MM, Song XZ, Shi Y, Yi JL, et al. Thermal performance of an aquifer thermal energy storage system: Insights from novel multilateral wells. Energy 2024; 294: 130915. https://doi.org/10.1016/j.energy.2024.130915
[34] Feng L, Xiao X, Wang E, Lu C, Li C, Li G, et al. Oil mobility evaluation and light hydrocarbon restoration in shale reservoirs: a review. Gondwana Res. 2026; 154: 122‑54. https://doi.org/10.1016/j.gr.2025.12.010
[35] Bai Y, Zhao J, Wu W. Methods to determine the upper limits of petrophysical properties in tight oil reservoirs: examples from the Ordos and Songliao Basins. J Pet Sci Eng. 2021; 196: 107983. https://doi.org/10.1016/j.petrol.2020.107983
[36] Xiao Y, Zeng Z, Zhang LQ, Wang JZ, Wang YW, Huang CG. A reduction‑consistent phase field model for non‑isothermal multiphase flows of N immiscible incompressible fluids. Int J Heat Mass Transf. 2024; 228: 125657. https://doi.org/10.1016/j.ijheatmasstransfer.2024.125657
[37] Lee D, Huh JY, Jeong D, Shin J, Yun A, Kim J. Physical, mathematical, and numerical derivations of the Cahn-Hilliard equation. Comput Mater Sci. 2014; 81: 216‑25. https://doi.org/10.1016/j.commatsci.2013.08.027
[38] Liu H, Cheng AJ, Wang H, Zhao J. Time‑fractional Allen-Cahn and Cahn-Hilliard phase‑field models and their numerical investigation. Comput Math Appl. 2018; 76(8): 1876‑92. https://doi.org/10.1016/j.camwa.2018.07.036
[39] Boyer F, Lapuerta C, Minjeaud S, Piar B, Quintard M. Cahn‑Hilliard/Navier‑Stokes model for the simulation of three‑phase flows. Transp Porous Med. 2010; 88: 463‑83. https://doi.org/10.1007/s11242-009-9408-z
[40] Djoko J, Koko J. Numerical methods for the Stokes and Navier‑Stokes equations driven by threshold slip boundary conditions. Comput Methods Appl Mech Eng. 2016; 305: 936‑58. https://doi.org/10.1016/j.cma.2016.03.026
[41] Chen Z, Huan G, Ma Y. Computational Methods for Multiphase Flows in Porous Media. Philadelphia: SIAM; 2006. https://doi.org/10.1137/1.9780898718942
[42] Xu F, Song X, Song Z, Wang G, Shi Y. Heat extraction and interlayer interference in heterogeneous multi‑layer commingled production oil reservoir. 56th US Rock Mechanics/Geomechanics Symposium. 2022; Santa Fe, New Mexico, USA, June 2022. Paper Number: ARMA-2022-0504. https://doi.org/10.56952/ARMA-2022-0504
Downloads
Published
Issue
Section
License
Copyright (c) 2026 Fuqiang Xu, Zhigang Liu, Can Ji, Rubing Han, Yan Liu, Aimin Chen

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
All the published articles are licensed under the terms of the Creative Commons Attribution Non-Commercial License (CC BY-NC 4.0) which permits unrestricted, non-commercial use, distribution and reproduction in any medium, provided the work is properly cited.


