摘要

Modeling viscous fingering perturbations in stochastic fields have always been a formidable task for geoscientists. Although fingering and channeling are often used synonymously to describe uneven front displacements, in fact, both phenomena represent different bypassing effects on hydrocarbon ultimate recovery. However, fingering perturbations in stochastic fields to invoke them during simulation often lead to severe channeling. Hence, the focus of this paper is to develop a finite volume numerical simulator with a level of certainty where fingering perturbations in stochastic fields do not result into severe channeling. The simulator development involves discretization of convection-diffusion and pressure equations with total variation diminishing and central differencing schemes, respectively. Correlated stochastic fields for simulation runs were created using convolution with a Gaussian filter achieved through 2D fast and inverse Fourier transform, Dykstra-Parsons coefficient, random porosity generator, and autocorrelation lengths. The obtained pressure and tracer concentration profiles were solved by iterative implicit pressure and explicit concentration approach, respectively. Numerical results of concentration profiles indicated that, unlike flow channeling, viscous fingering occurred in both homogeneous and heterogeneous fields. The obtained flow pattern map captures these flow regimes, namely fluid mobility fingering, channeling, and heterogeneous fingering. Sensitivity analyses of diffusion coefficient contributions to channeling and fingering phenomena are of extreme importance. The simulator shows a good validation with analytical solution while insensitive to grid refinements and spurious oscillations. The simulator is an improvement on first-order scheme simulators with the capability of computationally tracking fingering perturbations in stochastic fields which often result into channeling anomalies.