### Abstract

The modeling and numerical simulation of fluid flow in naturally fractured carbonate karst reservoirs are extremely challenging due to non-Darcy flow in vugs and caves connected by fracture networks. The momentum balance of such flow has been shown to be better described by the Brinkman equation both physically and mathematically, and many methods have been proposed in the literature dealing with the steady-state Brinkman model. We carry Brinkman's idea one step further and propose a transient flow model which consists of the Brinkman equation and a generalized material balance equation, and the latter has proven to be exact in the fractured carbonate karst reservoir. Finite differences are implemented for the solution of the proposed transient flow model. This solution method is more straightforward, easier to derive and implement, and more apt to generalization from 2D to 3D cases than alternative techniques. Numerical simulation of the transient Brinkman model requires explicit solution of not only pressure at the center of each grid block, but also velocities at the interfaces between the blocks, which exaggerates the computational cost and makes the computational process more difficult and less stable. In this paper, we propose a simplified finite difference formulation of the transient Brinkman model, which significantly reduces the computational time of the simulation process, and improves accuracy and stability of the simulation results. We update our reservoir simulator with this new formulation and illustrate it with a complex 3D naturally fractured carbonate karst reservoir model. The results of this study form the foundation for future 3D multi-phase reservoir cases.

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

Title of host publication | Society of Petroleum Engineers - SPE Annual Technical Conference and Exhibition, ATCE 2015 |

Publisher | Society of Petroleum Engineers (SPE) |

Pages | 5420-5431 |

Number of pages | 12 |

Volume | 2015-January |

ISBN (Electronic) | 9781510813229 |

Publication status | Published - 2015 |

Externally published | Yes |

Event | SPE Annual Technical Conference and Exhibition, ATCE 2015 - Houston, United States Duration: 28 Sep 2015 → 30 Sep 2015 |

### Other

Other | SPE Annual Technical Conference and Exhibition, ATCE 2015 |
---|---|

Country | United States |

City | Houston |

Period | 28/9/15 → 30/9/15 |

### Fingerprint

### ASJC Scopus subject areas

- Fuel Technology
- Energy Engineering and Power Technology

### Cite this

*Society of Petroleum Engineers - SPE Annual Technical Conference and Exhibition, ATCE 2015*(Vol. 2015-January, pp. 5420-5431). Society of Petroleum Engineers (SPE).

**Unified finite difference modeling of transient flow in naturally fractured carbonate karst reservoirs - A 3D case study.** / He, Jie; Killough, John E.; Fadlelmula, Mohamed; Fraim, Michael.

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

*Society of Petroleum Engineers - SPE Annual Technical Conference and Exhibition, ATCE 2015.*vol. 2015-January, Society of Petroleum Engineers (SPE), pp. 5420-5431, SPE Annual Technical Conference and Exhibition, ATCE 2015, Houston, United States, 28/9/15.

}

TY - GEN

T1 - Unified finite difference modeling of transient flow in naturally fractured carbonate karst reservoirs - A 3D case study

AU - He, Jie

AU - Killough, John E.

AU - Fadlelmula, Mohamed

AU - Fraim, Michael

PY - 2015

Y1 - 2015

N2 - The modeling and numerical simulation of fluid flow in naturally fractured carbonate karst reservoirs are extremely challenging due to non-Darcy flow in vugs and caves connected by fracture networks. The momentum balance of such flow has been shown to be better described by the Brinkman equation both physically and mathematically, and many methods have been proposed in the literature dealing with the steady-state Brinkman model. We carry Brinkman's idea one step further and propose a transient flow model which consists of the Brinkman equation and a generalized material balance equation, and the latter has proven to be exact in the fractured carbonate karst reservoir. Finite differences are implemented for the solution of the proposed transient flow model. This solution method is more straightforward, easier to derive and implement, and more apt to generalization from 2D to 3D cases than alternative techniques. Numerical simulation of the transient Brinkman model requires explicit solution of not only pressure at the center of each grid block, but also velocities at the interfaces between the blocks, which exaggerates the computational cost and makes the computational process more difficult and less stable. In this paper, we propose a simplified finite difference formulation of the transient Brinkman model, which significantly reduces the computational time of the simulation process, and improves accuracy and stability of the simulation results. We update our reservoir simulator with this new formulation and illustrate it with a complex 3D naturally fractured carbonate karst reservoir model. The results of this study form the foundation for future 3D multi-phase reservoir cases.

AB - The modeling and numerical simulation of fluid flow in naturally fractured carbonate karst reservoirs are extremely challenging due to non-Darcy flow in vugs and caves connected by fracture networks. The momentum balance of such flow has been shown to be better described by the Brinkman equation both physically and mathematically, and many methods have been proposed in the literature dealing with the steady-state Brinkman model. We carry Brinkman's idea one step further and propose a transient flow model which consists of the Brinkman equation and a generalized material balance equation, and the latter has proven to be exact in the fractured carbonate karst reservoir. Finite differences are implemented for the solution of the proposed transient flow model. This solution method is more straightforward, easier to derive and implement, and more apt to generalization from 2D to 3D cases than alternative techniques. Numerical simulation of the transient Brinkman model requires explicit solution of not only pressure at the center of each grid block, but also velocities at the interfaces between the blocks, which exaggerates the computational cost and makes the computational process more difficult and less stable. In this paper, we propose a simplified finite difference formulation of the transient Brinkman model, which significantly reduces the computational time of the simulation process, and improves accuracy and stability of the simulation results. We update our reservoir simulator with this new formulation and illustrate it with a complex 3D naturally fractured carbonate karst reservoir model. The results of this study form the foundation for future 3D multi-phase reservoir cases.

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

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

M3 - Conference contribution

VL - 2015-January

SP - 5420

EP - 5431

BT - Society of Petroleum Engineers - SPE Annual Technical Conference and Exhibition, ATCE 2015

PB - Society of Petroleum Engineers (SPE)

ER -