Real-Space Renormalization Estimates for Two-Phase Flow in Porous Media

We present a spatial renormalization group algorithm to handle immiscibletwo-phase flow in heterogeneous porous media. We call this algorithmFRACTAM-R, where FRACTAM is an acronym for Fast Renormalization Algorithmfor Correlated Transport in Anisotropic Media, and the R stands for relativepermeability. Originally, FRACTAM was an approximate iterative process thatreplaces the L × L lattice of grid blocks, representing the reservoir,by a (L/2) × (L/2) one. In fact, FRACTAM replaces the original L× L lattice by a hierarchical (fractal) lattice, in such a way thatfinding the solution of the two-phase flow equations becomes trivial. Thistriviality translates in practice into computer efficiency. For N=L ×L grid blocks we find that the computer time necessary to calculatefractional flow F(t) and pressure P(t) as a function of time scales as N^1.7 for FRACTAM-R. This should be contrasted with thecomputational time of a conventional grid simulator N^2.3. The solution we find in this way is an accurateapproximation to the direct solution of the original problem.