Modeling complex transport processes in naturally fractured reservoirs (NFRs) using classical continuum models may not be practically possible, because using classical algorithms for the detailed structure of fracture-matrix system requires unreasonable computational time. Also, fractured reservoirs are highly heterogeneous, and finite-difference calculations for such models often cause convergence problems. In addition to these, an exact representation of a complex fracture network in classical continuum modeling algorithms is highly difficult. An alternative is to use a non-classical technique known as the Random Walk Particle Tracking (RWPT) algorithm.We showed earlier (Stalgorova and Babadagli, 2009) that the random walk (RW) technique can be adapted to model miscible flooding in a fractured porous medium at the lab scale. The unknown parameters used to match the model results were only diffusion coefficients for oil and solvent, as the diffusive/dispersive transport (effective if fracture and matrix) was coupled with viscous (effective in fracture) and gravity (effective in fracture and matrix) displacement. Advantages of this method over classical simulation are: (1) shorter computational time, which allows avoidance of simplifications, and (2) the ability to model the matrix-fracture diffusion process without any transfer function.In the present paper, we modified this lab scale RW model for field scale applications. For validation, a series of tracer test results from the Midale field in Canada was used. Fracture network model was constructed based on geological data, and then we used the RWPT model to calibrate the fracture network against tracer test results.We performed a sensitivity analysis to identify the importance of different parameters for the simulation results. The new model and observations can be used to validate and calibrate stochastically generated fracture network models and to estimate the EOR performance of NFRs.
IntroductionFluid flow in fractured geological media is very common and has been the subject of many research studies in different engineering disciplines including oil and gas production, enhanced oil recovery, groundwater contamination, nuclear waste disposal, and CO 2 sequestration in underground reservoirs. Yet its simulation is quite a challenging task. The main reason is the complexity of this type of heterogeneity and its representation in mathematical models.Fluid flow modeling starts from reservoir characterization and representation of a geological model. It is mainly the description of fracture network geometry, which also describes the matrix size distribution. Characteristically, there is not enough information available to describe the properties of each single fracture. Instead, only bigger fractures are defined directly, and a set of smaller fractures is assigned by a number of stochastic characteristics, such as distributions of fracture directions, lengths, widths, permeabilities, etc. Once reservoir characterization and fracture network description (static m...