Pore-Scale Reconstruction and Simulation of Non-Darcy Flow in Synthetic Porous Rocks

Yixin Zhao, Guangpei Zhu, Cun Zhang, Shimin Liu, Derek Elsworth, Tong Zhang

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

We image synthetic porous rocks of varied porosity and pore size by micro–computed tomography with pore-scale finite element modeling representing the pore space for single-phase fluid flow. The simulations quantify the key features of microscale flow behavior in the synthetic cores. The smaller the permeability, the greater the critical pressure gradient required for the onset of non-Darcy fluid flows, and the easier the emergence of nonlinear seepage within the tested cores. The relationship between permeability and porosity from different methods shows a power law correlation with pressure. Structural heterogeneity and anisotropy of the pore systems are shown to have a significant impact on transport through the three-dimensional pore model—exhibiting irregular flow fields. The simulated permeabilities of the tested cores vary by a factor of 2–5 depending on the fluid flow directions. With the increase of flow seepage velocity, the flow regime deviates from Darcy linear relationship and non-Darcy behavior (inertia) leads to a reduction in the effective permeability of the core. Both experiments and modeling demonstrate that the larger the porosity and permeability, the stronger the non-linear phenomenon of the seepage within the pore space. A method is proposed to estimate the non-Darcy coefficient β based on simulation results, which provides a good prediction for all the tested cores. A new equation is established to predict the transition of flow patterns as a function of the apparent permeability K* and Reynolds number Re. The K*-Re model provides a theoretical approach to dynamically describe the transition from Darcy to Forchheimer flow.

Original languageEnglish (US)
Pages (from-to)2770-2786
Number of pages17
JournalJournal of Geophysical Research: Solid Earth
Volume123
Issue number4
DOIs
StatePublished - Apr 2018

Fingerprint

Seepage
Flow of fluids
Porosity
Rocks
rocks
permeability
porosity
rock
simulation
seepage
fluid flow
pore space
Pressure gradient
Flow patterns
Pore size
Tomography
Flow fields
Reynolds number
Anisotropy
single-phase flow

All Science Journal Classification (ASJC) codes

  • Geophysics
  • Geochemistry and Petrology
  • Earth and Planetary Sciences (miscellaneous)
  • Space and Planetary Science

Cite this

@article{16277511f4fc48649c8adc282d131f5e,
title = "Pore-Scale Reconstruction and Simulation of Non-Darcy Flow in Synthetic Porous Rocks",
abstract = "We image synthetic porous rocks of varied porosity and pore size by micro–computed tomography with pore-scale finite element modeling representing the pore space for single-phase fluid flow. The simulations quantify the key features of microscale flow behavior in the synthetic cores. The smaller the permeability, the greater the critical pressure gradient required for the onset of non-Darcy fluid flows, and the easier the emergence of nonlinear seepage within the tested cores. The relationship between permeability and porosity from different methods shows a power law correlation with pressure. Structural heterogeneity and anisotropy of the pore systems are shown to have a significant impact on transport through the three-dimensional pore model—exhibiting irregular flow fields. The simulated permeabilities of the tested cores vary by a factor of 2–5 depending on the fluid flow directions. With the increase of flow seepage velocity, the flow regime deviates from Darcy linear relationship and non-Darcy behavior (inertia) leads to a reduction in the effective permeability of the core. Both experiments and modeling demonstrate that the larger the porosity and permeability, the stronger the non-linear phenomenon of the seepage within the pore space. A method is proposed to estimate the non-Darcy coefficient β based on simulation results, which provides a good prediction for all the tested cores. A new equation is established to predict the transition of flow patterns as a function of the apparent permeability K* and Reynolds number Re. The K*-Re model provides a theoretical approach to dynamically describe the transition from Darcy to Forchheimer flow.",
author = "Yixin Zhao and Guangpei Zhu and Cun Zhang and Shimin Liu and Derek Elsworth and Tong Zhang",
year = "2018",
month = "4",
doi = "10.1002/2017JB015296",
language = "English (US)",
volume = "123",
pages = "2770--2786",
journal = "Journal of Geophysical Research",
issn = "0148-0227",
publisher = "American Geophysical Union",
number = "4",

}

Pore-Scale Reconstruction and Simulation of Non-Darcy Flow in Synthetic Porous Rocks. / Zhao, Yixin; Zhu, Guangpei; Zhang, Cun; Liu, Shimin; Elsworth, Derek; Zhang, Tong.

In: Journal of Geophysical Research: Solid Earth, Vol. 123, No. 4, 04.2018, p. 2770-2786.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Pore-Scale Reconstruction and Simulation of Non-Darcy Flow in Synthetic Porous Rocks

AU - Zhao, Yixin

AU - Zhu, Guangpei

AU - Zhang, Cun

AU - Liu, Shimin

AU - Elsworth, Derek

AU - Zhang, Tong

PY - 2018/4

Y1 - 2018/4

N2 - We image synthetic porous rocks of varied porosity and pore size by micro–computed tomography with pore-scale finite element modeling representing the pore space for single-phase fluid flow. The simulations quantify the key features of microscale flow behavior in the synthetic cores. The smaller the permeability, the greater the critical pressure gradient required for the onset of non-Darcy fluid flows, and the easier the emergence of nonlinear seepage within the tested cores. The relationship between permeability and porosity from different methods shows a power law correlation with pressure. Structural heterogeneity and anisotropy of the pore systems are shown to have a significant impact on transport through the three-dimensional pore model—exhibiting irregular flow fields. The simulated permeabilities of the tested cores vary by a factor of 2–5 depending on the fluid flow directions. With the increase of flow seepage velocity, the flow regime deviates from Darcy linear relationship and non-Darcy behavior (inertia) leads to a reduction in the effective permeability of the core. Both experiments and modeling demonstrate that the larger the porosity and permeability, the stronger the non-linear phenomenon of the seepage within the pore space. A method is proposed to estimate the non-Darcy coefficient β based on simulation results, which provides a good prediction for all the tested cores. A new equation is established to predict the transition of flow patterns as a function of the apparent permeability K* and Reynolds number Re. The K*-Re model provides a theoretical approach to dynamically describe the transition from Darcy to Forchheimer flow.

AB - We image synthetic porous rocks of varied porosity and pore size by micro–computed tomography with pore-scale finite element modeling representing the pore space for single-phase fluid flow. The simulations quantify the key features of microscale flow behavior in the synthetic cores. The smaller the permeability, the greater the critical pressure gradient required for the onset of non-Darcy fluid flows, and the easier the emergence of nonlinear seepage within the tested cores. The relationship between permeability and porosity from different methods shows a power law correlation with pressure. Structural heterogeneity and anisotropy of the pore systems are shown to have a significant impact on transport through the three-dimensional pore model—exhibiting irregular flow fields. The simulated permeabilities of the tested cores vary by a factor of 2–5 depending on the fluid flow directions. With the increase of flow seepage velocity, the flow regime deviates from Darcy linear relationship and non-Darcy behavior (inertia) leads to a reduction in the effective permeability of the core. Both experiments and modeling demonstrate that the larger the porosity and permeability, the stronger the non-linear phenomenon of the seepage within the pore space. A method is proposed to estimate the non-Darcy coefficient β based on simulation results, which provides a good prediction for all the tested cores. A new equation is established to predict the transition of flow patterns as a function of the apparent permeability K* and Reynolds number Re. The K*-Re model provides a theoretical approach to dynamically describe the transition from Darcy to Forchheimer flow.

UR - http://www.scopus.com/inward/record.url?scp=85045401925&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85045401925&partnerID=8YFLogxK

U2 - 10.1002/2017JB015296

DO - 10.1002/2017JB015296

M3 - Article

AN - SCOPUS:85045401925

VL - 123

SP - 2770

EP - 2786

JO - Journal of Geophysical Research

JF - Journal of Geophysical Research

SN - 0148-0227

IS - 4

ER -