A HIGHER-ORDER GODUNOV SCHEME COUPLED WITH DYNAMIC LOCAL GRID REFINEMENT FOR FLOW IN A POROUS-MEDIUM

Today's reservoir simulators are generally inefficient due to the use of low-order accurate schemes in space and time on regular grids. The numerical diffusion inherent in the single-point upstream weighting sc heme can cause large errors in the numerical solutions, which can only be removed by the use of very fine uniform grids with many grid block s. The aim of this work is to develop a method which can produce physi cally reliable and accurate results efficiently. A higher-order Goduno v scheme is coupled with dynamic grid adaptivity, where grid blocks ar e inserted in highly active regions of the how field and removed from regions of inactivity, thereby concentrating the computational effort where it is most needed. The quality of results computed by the adapti ve higher-order scheme are comparable with those computed by the highe r-order scheme on a uniform grid, globally refined to the level of the finest adaptive grid zones, while great savings in computer time are obtained. The adaptive high-order scheme is vastly superior compared t o the first-order scheme on a uniform or adaptive grid. The results pr esented demonstrate the benefits of the method in reservoir simulation .