Modeling of oil phase displacement processes in heterogeneous anisotropic reservoirs

Methods of computer modeling of slightly permeable anisotropic oil-bearing layers are needed nowadays because they give us a possibility to obtain a concept on filtration processes near producing and forcing wells in different practical situations and in this way to raise the level of exploitation of such layers significantly. On the other hand, they allow evaluation and taking into account some uncertainties which appear as a result of inefficient information on the structure and properties of the layer outside the wells. In order to investigate the practical aspects of supporting efficiency of oil production in anisotropic heterogeneous low permeable reservoirs on the base of combined finite-element-difference method for solving the non-stationary anisotropic piezoconductivity problem, modeling of distribution of layer pressure was carried out in the vicinity of the production and forcing wells taking into account the anisotropy of the permeability and conditions of the oil phase infiltration on the margins of the examined layer. It has been found that the intensity of filtration process between producing and forcing wells depends essentially on their spacing in both shear-isotropic and anisotropic oil-bearing layers. In addition the effect of oil phase permeability in shear direction dominates over the effect of permeability in axes directions. Starting from the obtained information for the effective exploitation of anisotropic slightly permeable layers we need to locate producing and forcing wells in the areas with relatively low permeability of the layer and especially to avoid the places with presence of shear permeability. It is important to locate the wells in such a way that blocking the oil in the direction of reduced permeability and fast depletion of the layer in the direction of increased permeability would not happen as well as mutual exchange between producing and forcing wells would not stop. While locating the system of specified wells within anisotropic layers of oil deposit it is necessary to conduct a systemic analysis of environmental anisotropy of layers aimed at such a location of these wells which would guarantee the effective dynamics of filtration processes around them. Application of quadratic isoparametric approximation of finite-elemental net of examined area of oil-bearing layer and implicit differential time approximation brings to increase of precision and stability of numerical solution of the problem.