Timestep selection during streamline simulation via transverse flux correction
Streamline simulators have received increased attention because of their ability to effectively handle multimillion cell detailed geologic models and large simulation models. The efficiency of streamline simulation has relied primarily on their ability to take large timesteps with fewer pressure solutions within an IMPES formulation. However, unlike conventional finite-difference simulators, no clear guidelines are currently available for the choice of timestep for pressure and velocity updates. That is why we need largely an uncontrolled approximation, either managed by engineering judgment or by potentially time-consuming timestep size sensitivity studies early in a project. This will clearly lead us to the lack of understanding of numerical stability and error estimates during the solution. This research presents a novel approach for timestep selection during streamline simulation that is based on three elements. First, we reformulate the equations to be solved by a streamline simulator to include all of the three-dimensional flux terms - both aligned with and transverse to the flow directions. These transverse flux terms are totally neglected within the existing streamline simulation formulations. Second, we propose a simple grid-based corrector algorithm to update the saturation to account for the transverse flux. Third, we provide a discrete CFL (Courant-Friedrich-Levy) formulation for the corrector step that leads to a mechanism to ensure numerical stability via the choice of a stable timestep for pressure updates. This discrete CFL formulation now provides us with the same tools for timestep control as are available within conventional reservoir simulators. We demonstrate the validity and utility of our approach using a series of numerical experiments in homogeneous and heterogeneous ? five-spot patterns at various mobility ratios. For these numerical experiments, we pay particular attention to favorable mobility ratio displacements, as they are known to be challenging to streamline simulation. Our results clearly demonstrate the impact of the transverse flux correction on the accuracy of the solution and on the appropriate choice of timestep, across a range of mobility ratios. The proposed approach eliminates much of the subjectivity associated with streamline simulation, and provides a basis for automatic control of pressure timestep within full field streamline applications.