Incorporation of Mixed Finite Element Methods in Compositional Simulation for Reduction of Numerical Dispersion

R.E. Ewing, R.F. Heinemann
1983 SPE Reservoir Simulation Symposium   unpublished
Previous studies have shown that standard finite difference techniques cause numerical dispersion and grid orientation problems when used to simulate enhanced recovery processes with adverse mobility ratios. In compositional simulation, numerical dispersion can diffuse sharp fluid interfaces yielding erroneous predictions of fluid compositions and corresponding errors in the velocities of the miscible frontal advance. Numerical dispersion can also effect the computed locations of the boundaries
more » ... s of the boundaries of the regions of single-phase and twophase flow. Inaccurate fluid velocities and suboptimal use of upstream weighting of transport terms combine to cause many aspects of the numerical dispersion and grid orientation problems. A mixed finite element method has been developed to obtain more accurate approximations to the fluid velocities. In this method the Darcy velocities are considered as primary variables together with the total fluid pressure. Although finite element techniques are used to compute the more accurate fluid velocities, these velocities ar.e then incorporated into a more standard finite difference method for the bulk of the simulation process. This paper presents the use of mixed methods in a two-dimensional finite difference compositional simulator to reduce problems caused by numerical dispersion. Comparisons are made with a standard finite difference simulator on problems involving immiscible displacement and multiple contact miscibility phenomena.
doi:10.2118/12267-ms fatcat:cojiowkjh5axzjxe7kjocscnwq