UDK 622.276.5

Y. S. Oryngozhin, Dr. of Tech. Sci., Professor

A. A. Shukmanova, doctoral PhD

 

K. I. Satpayev Kazakh National Technical University

 

THE TRANSFORMATION PARAMETERS OF THE GEOLOGICAL MODEL PARAMETERS TO FILTER GRID MODEL

 

For a correct simulation of the formation region of interest is paramount consideration the variability of filtration and reservoir properties of the rocks , as well as the inclusion of spatial geometry of the structure of the productive formations and their discontinuity .

In a non-uniform thickness , permeability, and other physical properties of the formation, but on the area exposed for the process of oil displacement occurs in all areas, but at different rates depending on the changes in physical properties. Another character in the process acquires a mercurial formation characterized by discontinuities in one form or another . Real formations may have lenticular structure , flaking collectors on individual streaks , their thinning , etc. This ultimately leads to losses in the extraction of oil .

The discontinuity of the formation may occur as a result of disjunctive dislocations. In many cases, the discontinuity of the formation is the result of combined effects of various factors on the process of sedimentation.

Even the detailed geological model consisting of millions of cells that can not display the actual structure of the layers in the space between the holes . Created probabilistic models are in most cases non-uniform , but continuous between wells. These models do not reflect discontinuous reservoirs , thereby deducing the CIN weakly dependent on well spacing . In the transition from a three-dimensional geological model , which can consist of tens of millions of cells to the hydrodynamic model consisting of tens or hundreds of thousands of cells , there is desensitization variability of reservoir properties .

In general, the inhomogeneity during the transition from the formation geological model to hydrodynamic produced, at least two methods:

1 ) construct a modified phase curves , taking into account the heterogeneity of sand members of permeability ;

2) by calculating the sweep efficiency of the selected pattern of wells , taking into account the discontinuous layers .

Currently , when averaged geological models , the established practice of their preparation for reservoir simulation development process is to divide the model to approximately equal intervals , the amount of which is determined by the capabilities of the software and hardware , and averaging the characteristics of reservoir within these intervals. With this averaging often blurs and fades discontinuous layers .

To improve the quality produced by geological and engineering models should be significantly strengthened the role of geological experts to conduct detailed correlation of simulated objects , making structural and lithologic patterns. Detailed correlation of well sections is the main method for studying the internal structure of the subsoil , built on the principle of comparing the well sections to highlight in cuts and tracking in area reservoirs and impermeable layers between them.

One way to study the detailed geological structure of the reservoir is to compare the well sections of the zonal method of correlation. Zonal correlation is to compare a section of the reservoir and the allocation within individual zones , which are characterized by common lithologic and physical properties spread out over the square, and known differences vertically. Detail or zonal correlation is carried out in relation to one or more selected Tetrad horizons. The more the horizons of reference can be used in the correlation , the more accurate the results. In wells and reservoirs collectors tight partitions between them isolated using conditionally limits collectors .

In separate wells of a partial fusion occurs designated zones , and sometimes all of the intervals , are connected together to form one monolithic permeable formation . The provision of such zones within the reservoir provides a more detailed understanding of the geological structure of it and build the zonal maps, which are the basis for the design and adjustment of the accepted system design, as well as measures used to influence the formation.

For the design of field development , as well as for calculating the reserves should be during the detailed correlation allocate zonal intervals thin , and not limited to the dismemberment of the cut deposits only a few large packs . Practice shows that the minimum thickness of the layers available for that purpose, is approximately 0.5-1.0 m. These values can be considered decisive for the selection of individual zonal intervals .

In any case, the choice of the number of intervals productive zone of the geological formations in the filter and the corresponding models can not be carried out such consolidation , in which some sections of clay with thickness of 0.5-1.0 m, and more neighboring would be absorbed zonal intervals.

To characterize the discontinuity layer is necessary to determine the spread over the area of each selected interval zone or layer. Under existing rules mining geometry , the border of the zone interval on the map is held at half the distance between two adjacent wells , one in the context of which this zone interval is present , and the other is wedged out . As already mentioned, some individual wells zonal intervals merge together to form an area of their distribution zones merge. The boundaries of these zones on the card are set in the same manner as during the boundaries of zone intervals .

Areas dedicated zone intervals are filled to their respective grid cells of the geological model with preservation of all allocated impermeable boundaries .

Grid geological model usually is using in modern software , has a large dimension of the order of millions of grid blocks. The direct use of such a network, a digitized appropriate filtration-capacity parameters in the filter model almost makes it impossible to perform multiple calculations forecasting technological parameters . Hence the dimension of the grid geological model is reduced by 1-2 orders of magnitude with the appropriate procedures of averaging .

At coincidence grid geological and model input data is done directly from the geological and mathematical models . If the geological data grid is not built, its construction and filling is made using cards of relevant parameters by various means (graphic, software , manual, including funds that are part of the filter model).

A common method of averaging parameters detailed geological grid on the cell enlarged geologic grid technology is a rescaling - Upscaling. In this case, the average porosity values are calculated by arithmetic averaging , averaging is performed by weighting the saturation parameter pore volume. To obtain the values of permeability units in the technology used Upscaling method diagonal permeability tensor , built on the numerical solution of incompressible fluid -phase filtration for each block filtration model . For each block of the dynamic grid solves the problem of filtering a single-phase incompressible flow in three directions. This method allows the anisotropic permeability field that reflects the presence of a cross-sectional thin streaks - non-reservoir and collectors with different properties.

It will be appreciated that averaging methods can be used only within the areas of permeable reservoir rocks and should be averaged properties packs comprising impermeable sections with thicknesses over 0.5-1.0 m.

At present the problem of averaging given much attention , it is devoted a lot of publications. However, the correct solution to the problem of averaging the equations of multiphase flow is absent. Therefore averaging algorithms used in modern software packages , in varying degrees, are approximate , and the procedure for converting the grid geological model of the grid filter models require in each case an individual approach and conduct special test calculations . As a first test calculation is checked coincidence of geological reserves of geological and models.

The recent emergence of modern simulation software , high-performance computer with a clock speed of over three gigahertz and parallel cluster computing systems to create and effectively calculate filtration models , numbering several million cells , thereby minimizing the averaging procedure for converting the parameters of a geological model parameters filtration grid model .