### Abstract

The presence of cavities connected by fracture networks at multiple levels make the simulation of fluid flow in naturally fractured carbonate karst reservoirs a challenging problem. The challenge arises in properly treating the Darcy and non-Darcy flow in the different areas of fractured medium. In this paper, we present a single-phase transient flow model which is based on the Stokes-Brinkman equation and a generalized material balance equation. The generalized material balance equation proves to be exact in both cavities and porous media, and the Stokes-Brinkman equation mathematically combines Darcy and Stokes flow, thus allowing a seamless transition between the cavities and porous media with only minor amounts of perturbation introduced into the solutions. Finite differences are implemented for the solution of the proposed transient flow model. This solution method provides a smooth transition from standard multiple-porosity/permeability reservoir simulators and moreover, it is physically more straightforward, mathematically easier to derive and implement, and more apt to generalization from two-dimensional to three-dimensional cases than alternative techniques. Application of the derived transient flow model is shown by examples of three fine-scale 2-D geological models. The first two models, although simple, provide verification of the proposed transient flow model. The third example presents a more complex and realistic geological model derived from multiple-point statistics simulation technique with the second model used as the training image. The results of the third model form the foundation for future study of multi-phase and 3-D reservoir cases.

Original language | English |
---|---|

Title of host publication | Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2015 |

Publisher | Society of Petroleum Engineers |

Pages | 1229-1242 |

Number of pages | 14 |

Volume | 2 |

ISBN (Electronic) | 9781510800618 |

Publication status | Published - 2015 |

Externally published | Yes |

Event | SPE Reservoir Simulation Symposium 2015 - Houston, United States Duration: 23 Feb 2015 → 25 Feb 2015 |

### Other

Other | SPE Reservoir Simulation Symposium 2015 |
---|---|

Country | United States |

City | Houston |

Period | 23/2/15 → 25/2/15 |

### Fingerprint

### ASJC Scopus subject areas

- Modelling and Simulation
- Geochemistry and Petrology

### Cite this

*Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2015*(Vol. 2, pp. 1229-1242). Society of Petroleum Engineers.

**A unified finite difference model for the simulation of transient flow in naturally fractured carbonate karst reservoirs.** / He, Jie; Killough, John E.; Fadlelmula, Mohamed; Fraim, Michael.

Research output: Chapter in Book/Report/Conference proceeding › Conference contribution

*Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2015.*vol. 2, Society of Petroleum Engineers, pp. 1229-1242, SPE Reservoir Simulation Symposium 2015, Houston, United States, 23/2/15.

}

TY - GEN

T1 - A unified finite difference model for the simulation of transient flow in naturally fractured carbonate karst reservoirs

AU - He, Jie

AU - Killough, John E.

AU - Fadlelmula, Mohamed

AU - Fraim, Michael

PY - 2015

Y1 - 2015

N2 - The presence of cavities connected by fracture networks at multiple levels make the simulation of fluid flow in naturally fractured carbonate karst reservoirs a challenging problem. The challenge arises in properly treating the Darcy and non-Darcy flow in the different areas of fractured medium. In this paper, we present a single-phase transient flow model which is based on the Stokes-Brinkman equation and a generalized material balance equation. The generalized material balance equation proves to be exact in both cavities and porous media, and the Stokes-Brinkman equation mathematically combines Darcy and Stokes flow, thus allowing a seamless transition between the cavities and porous media with only minor amounts of perturbation introduced into the solutions. Finite differences are implemented for the solution of the proposed transient flow model. This solution method provides a smooth transition from standard multiple-porosity/permeability reservoir simulators and moreover, it is physically more straightforward, mathematically easier to derive and implement, and more apt to generalization from two-dimensional to three-dimensional cases than alternative techniques. Application of the derived transient flow model is shown by examples of three fine-scale 2-D geological models. The first two models, although simple, provide verification of the proposed transient flow model. The third example presents a more complex and realistic geological model derived from multiple-point statistics simulation technique with the second model used as the training image. The results of the third model form the foundation for future study of multi-phase and 3-D reservoir cases.

AB - The presence of cavities connected by fracture networks at multiple levels make the simulation of fluid flow in naturally fractured carbonate karst reservoirs a challenging problem. The challenge arises in properly treating the Darcy and non-Darcy flow in the different areas of fractured medium. In this paper, we present a single-phase transient flow model which is based on the Stokes-Brinkman equation and a generalized material balance equation. The generalized material balance equation proves to be exact in both cavities and porous media, and the Stokes-Brinkman equation mathematically combines Darcy and Stokes flow, thus allowing a seamless transition between the cavities and porous media with only minor amounts of perturbation introduced into the solutions. Finite differences are implemented for the solution of the proposed transient flow model. This solution method provides a smooth transition from standard multiple-porosity/permeability reservoir simulators and moreover, it is physically more straightforward, mathematically easier to derive and implement, and more apt to generalization from two-dimensional to three-dimensional cases than alternative techniques. Application of the derived transient flow model is shown by examples of three fine-scale 2-D geological models. The first two models, although simple, provide verification of the proposed transient flow model. The third example presents a more complex and realistic geological model derived from multiple-point statistics simulation technique with the second model used as the training image. The results of the third model form the foundation for future study of multi-phase and 3-D reservoir cases.

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

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

M3 - Conference contribution

AN - SCOPUS:84939427800

VL - 2

SP - 1229

EP - 1242

BT - Society of Petroleum Engineers - SPE Reservoir Simulation Symposium 2015

PB - Society of Petroleum Engineers

ER -