Cyclic Rhythms and Macroheterogenity of the Carbonated Intersaliferous Reservoirs in the Pripiat Basin ( Byelorussia )

The intersaliferous complex of the Pripiat basin was deposited during a rift stage that occurred at Hercynian cycle. Reef formations are locally present within this complex, proof that the basin had undergone significant variations of tectonic regime during sedimentation processes. The natural reservoirs rocks of this complex consist essentially of biogenic, bioclastic limestones and dolomites which had suffered important diagenetics transformations. The vertical succession of these reservoirs rocks point up the existence of cycles quantified by the mathematical methods of autocorrelation and spectral analysis. The application of appointed mathematical methods for interpretation of neutron gamma logs of natural reservoirs intervals allowed the structuring of studied oil fields in linear hydrodynamic blocks proved by the water injections made with intention to maintain reservoirs pressure during fields’ exploitation.


Introduction
The actual importance of oil resources and the progressive exhaustion of proven World reserves as we exploit them swiftly, require the increasing of effectiveness of prospecting works as well as the development and optimization of already exploited oilfields.Complex reservoirs are generally underexploited because it's difficult to control their heterogeneity which must be well known to predict the best exploitation regime in intention to maximize the recovery processes.Consequently, new approaches allowing more detailed petrophysical study of complex reservoirs rocks are needed.In fact, a good knowledge of carbonated reservoirs heterogeneity characterized by microheterogeneity linked to post-diagenetics transformations and macroheterogeneity related to combination of reservoir, false seals and screens rocks; make easier the orientation of prospection and exploitation works.
Based on geological and geophysical materials, this paper presents the use of the mathematical methods of autocorrelation in purpose to study and appreciate the heterogeneity of carbonated reservoirs in three oilfields (Davidovskaya, Yuzno-Sosnovskaya and Yuzno-Alexandrovskaya) located in the North-East of the Pripiat basin.The study brings a solution to the problem of hydrodynamic discontinuity encountered at the time of water injection during exploitation.

Geology
Located in the South-East of Byelorussia, the Pripiat basin is a negative structure of 250 × 150 km, where prospection works actively began in 1940.The basin occupies a part of the Precambrian Eastern-European platform and belongs to the large Pripiat-Donetz grabens system which links Eastern Byelorussia to Ukraina (Figure 1).
The crystalline basement of the basin which is composed of Achaean and Proterozoïc granites, granodiorites and gneiss is sheared by several faults systems [1].Some ones are East-West listric faults, which limit the basin and the others are East-West regional faults which subdivide the basin into three various structural zones: the Southern, Central and Northern zones [2,3].
The sedimentary cover of the Pripiat basin includes Precambrian (Neoproterozoïc), Palaeozoic (Devonian, Carboniferous and Permian), Mesozoïc and Cenozoïc formations subdivided into great lithological three  The upper complex which consists essentially of Mesoview, the basin presents a gr

Form
that the reservoirs of rganisms had locally formed reef co tasomatose, recrystallization which have conduced to the Cenozoic layers of low thicknesses.The formation of this last complex is related to the basin immersion during the alpine cycle.From economic point of eat interest as part of the enormous Pripiat-Donetz petroleum basin.The main producing stratigraphics intervals are inter and under salt formations of Devonian age.

Intersaliferous Complex
er petrographic studies showed the intersaliferous complex consist of biogenic (algae), bioclastic and crystalline limestones and dolomites represented on Board 1.
The algae and other o nstructions [6][7][8][9].These carbonated rocks have undergone various catagenetic processes of dissolution, me-formation of new minerals explained by the nature of vacuums observed in the rocks [10][11][12].In objective to control the reservoirs characteristics for a better productivity, several diagraphics methods of exploration (electric, acoustic, gamma, neutron-gamma and cavernometric) were used to study the interval of carbonated intersaliferous complex.Analysis and interpretations were done in association with cores data.The neutron gamma which provides information relative to the quantity of hydrogen present in rocks and therefore total porosity of rocks, was largely employed to evaluate the qualitative types of porosity: pores, caves, cracks [13,14] and determinate the porosity of the reservoirs.The content of clay in carbonates induces a rising of the radioactivity of Gamma type (Gamma Ray) and a reduction of the porosity [15].The results of porosities and clay contents obtained by interpretation of geophysical data gave possibility to establish the limited value of porosity between the reservoirs rocks and the seals in interval 2.5% -7.5% ith an average of 4.8%, corresponding to a limited clay of wells logs, give possibility to reveal the masked cyclic rhythms resul w content of 15% [5,13].The subdivision of the layers in reservoirs rocks and seals rocks made from porosity values and clay content, without taking into account the permeability values and the decrease of the layers pressure, led to non conformity between considered reservoirs and the test flow results.Authors' works of 1993 and 2000, which established that open porosity within the intersalt reservoirs rocks varies with the nature and the extent of the secondary transformations in the carbonates, were used as basis for the subdivision of the massive reservoirs in effective reservoirs rocks and false seals intervals.Indeed, it come out from this work that the rocks of the intersalt complex have a density which varies from 2.4 to 2.78 kg/m 3 for the carbonated rocks; and reaches more than 2.93 kg/m 3 for the evaporates [11].The open porosity of these rocks varies from 0.7 to more than 12% with more than 5% for the limestones.High values of porosity (12% and more) were fixed in dolomites of the Eastern part of the basin.The reservoirs rocks are characterized by a density lower than 2.7 kg/m 3 , a porosity higher than 2% and a chaotic spatial distribution (as well vertical as horizontal), which had generated a macroscopic heterogeneity of these natural reservoirs.The study of this heterogeneity gave possibility to make a typisation of these reservoirs and found a perfect explication to the negative results of tests flow obtained in false seals intervals [11,16].The massive structure of the reservoirs, branded by a horizontal oil water contact must normally make easier maintain of layers pressure through water injection during exploitation.But it lets appear a serious problem of hydrodynamic discontinuity between several wells of the studied Yuzno-Sosnovskaya and Yuzno-Alexandrovskaya fields.
The detailed study of sedimentation cycles and postsedimentation transformations will be useful to determina anential and point out cy ting from the processes of sedimentation tion l whic its to detect repeated regularities of profiles ed by a large een two unspecified points te blocks with similar hydrodynamic continuity in intention to implant successfully the water injection wells.

Rhythmicity of the Natural Reservoirs
The rhythmicity in sedimentation of carbonated form tions is conditioned by several factors.In particular co dition of carbonated sedimentation, the size of the particles and the carbonate contents depends of the depth and shape of the basin.The sediments configuration varies in time and space with the changes of depotbelt centers and rhythm of vertical movements, in particular oscillations changes of the basin depth.
Variations of the basin depth just like variation of porosity and CaCO 3 contents are sequ clic rhythms of various orders in the sedimentary column.The autocorrelation function and spectral analysis and post sedimentation [17][18][19][20][21].

The Autocorrelation Func
The autocorrelation is a mathematical tool used in signa treatment.It is a cross correlation of a signal by itself, h perm in a signal, like a periodic signal troubl amount of noise or a fundamental frequency of a signal which does not contains effectively this fundamental frequency, but implies it with several of its harmonics [22,23].
If given N is a time series representing the log neutron gamma on the entire studied reservoir interval and contains n points isolated from each other of ∆; the di  stance betw is called lag of length  .Therefore, the autocorrelation of lag  is the correlation of all the points    The autocorrelation function describes the internal ruct re of a process and permit to evaluate up to which point the process prese es it correlation function gives n mmetry, closure and Schwarz inequality [21].

Spectral Analysis
Spectral analysis comes in complement to the autocorrelation and consists in decomposition of a random comprising this variable.T ken is regarded as being the sum of a multitude of short simpler diagrams, having each one a regular sinusoidal form, amplitude and wave length in a starting point.The analysis is used to detect the masked periodicities in layers formation, and to determine the wave lengths (T) corresponding to the thicknesses of the cycles with structure of repetition within the geological cuts [21].It describes the dispersion for each period.
The time series can be expressed as follow: where: k-number of cycles (harmonics); 2πx   an Φ k : angular phase; A : cycle amplitude As: If we consider that cos (3) For a time series, if we have n points distant one from other of Δ and j one of these po cos sin harmonic k, we can determine the amplitude A k and the angular phase Φ k of its periodic form (undulatory) as follow: In the case of a discrete series with regular sinusoidal structure, the dispersion of the harmonic component k can be expressed by: Since the Fourier theorem affirms that the dispersion of a variable is equal to the sum of dispersions of its components deleted by the peaks of the spectrogram; ones can determines the dominant components of the variable by calculating and putting the dispersions (energies) values on the y-axis.

Application and Results
Considering neutron gamma logs measured with an indepth lag of sampling of 0.8 m inside the productive (D 3 zd -D 3 el) interval of less than 200 m thick, we have made up for each well, a time series.
The obtained time series were subjected to the autocorrelation and spectral analysis using "Autocor" software elaborate in the PASCAL language.Therefore, we proceed to the presentation of the time series in the form of model and we calculate the values of dispersion and complete the process by construction of a series of correlograms and spectrograms.The analysis of the periodicity is in fact based on the correlograms and spectrograms (Figure 3).On the correlograms, it was possible to determine graphically the ray for which the correlation (R o ) is null and, on the spectrograms, the harmonics materialized by peaks.The values of periods and dispersions of these harmonics are calculated.The correlograms and the spectrograms are subjects of analyse which is based on their form and indexes of heterogeneity.Like quantitative indexes of heterogeneity, one considers for the autocorrelation functions for which the correlation (R o ) is null (R o expressed in meters and determined on the correlograms); and for the spectral analysis, the maximum period (T max expressed in meters) of the main wave.These two selected parameters will have values proportional to the degree of organization of layers within the studied reservoirs.

Cyclicity of the Cuts and
cuts are more or less homogeneous ease of the ization of ts show a cyclic repetition of petrocillatory movements is twice higher th the total absence of cyclic organization of the cuts.For this last type, we observe a swift decrease of the correlograms, a fluctuation of the intensity ln and his oscillation around the x-axis.This type consists of spectra reflecting a great differentiation of the ln diagraphic struc-ture One notes a polyperiodic composition of the graphs, showing that the vertical variation of petrophysical pro-perties of reservoirs was relied to high frequency oscilla-tory movements.The cuts of this type, with R o = 10 m and T max < 50 m, show a regular repetition of the petro-physical properties (see Figure 4).

Cyclicity of the cuts
For this field, all identified correlograms and spectrograms are characterized by a regular pe odicity of the cut (type II, III and IV).They are diffe nt from each degree of organization of the cut.The alues (T < 50 m) are observed in the Ea and non reservoirs ro r less homo-113, 47, 127, 121, 128, 8, 118) the homogeneity is

Macroheterogeneity of the Natural Reservoirs
other only by the R o and T max values.Consequently, they Cyclicity analysis based on correlograms and spectrograms diagrams (obtained after treatment of the neutron gamma logs) reveal the inner organisation of the cuts and show four different types of cuts cyclicity in the studied fields (Figure 4).
The 1st type shows a slow decrease of the correlograms values from +1 towards the x-axis which is, crossed for R o values  15 m.The correlogram is prolonged under the x-axis, by a rectilinear way.This variation is associated with a fluctuation ln γ With a T max  145 m, this type reflect the existence of a long periodicity in the cuts.The vertical variation of petrophysical properties of the rocks is ordered by oscillatory movements which wavelength is almost equal to the total thickness of the layers.The (see Figure 4).
The 2nd type is characterized by a regular decr correlograms and an apparent periodic organ the cuts.The cu physical properties which can be explained by a vertical succession of reservoirs and non reservoirs rocks.The cyclic configuration of type II is characterized by R o = 10 -15 m and T max = 69 -90 m.These spectra testify that the frequency of the os an in type I.The 3rd type is close to the 2nd type, but differs of it by a more fine structure which indicates a reduction, within the interval, of cyclic repetition of petrophysical properties of the rocks.The spectra often present two spades poorly matched with close values of dispersion.The reservoirs are characterized by a great heterogeneity and a greater subdivision of the cuts with R o = 5 -10 m and T max = 50 -68 m (see Figure 4).
The 4th type reflects the maximum heterogeneity and where the cut consists of non thick reservoirs layers (Figure 5).
Considering that the high value of R o indicates a relative homogeneity of the cut described by the autocorrelation, the Davidovskaya field has a vertical heterogeneity related to the succession of reservoirs cks in the Eastern part of the field.The reservoirs rocks thicknesses increase towards the South-eastern and the central parts of the field.

Cyclicity of the cuts
Four (4) types are found in this field: The 1st type with for period T max > 180 m, reflects the existence in the cuts of a long periodicity controlled by oscillatory movements which wavelength is almost equal to the total thickness of the layers.Cuts acquire from Eastern and South-eastern parts are more o geneous (wells 126, 124, 119, and 56).In the South (wells 10 not so pronounced. The 2nd type is characterized by a periodic organization  shows the existence of a zonation in distribution of these parameters (Figure 6).The charts of the ray R o and Tmax let observe several zones with minimal values (Figure 7).
They are the western zones (wells 32, 31, 33, 7), the eastern zones (wells 6, 53) and the heart of the structure (wells 5, 39) and south-east (wells 41, 50) of the field.The fields Yuzno-Alexandrovskaya and Yuzno-Sosnovskaya present a great similarity as well from the point of view of the space distribution of the heterogeneity as from the values of the examined parameters.

Structuring of the Oil Fields
The study of the sedimentation rhythms in the considered fields and the quantitative analysis of heterogeneity, as reflection of spatial distribution of the reservoirs properties, give possibility to propose charts of reservoirs structures.The zonation is made by juxta sition of the 8) and are examined as the , formed

Yuzno-Alexandrovskaya Field
charts of correlation ray (R o ) and maximum period (T max ).The various zones are blocks whose limits correspond to facies variations (Figure

Cyclicity of the cuts
The correlograms and the spectrograms gather are of 3 various types: The first type not structured, with a big ray (R ) and a result of the tectonic activities of low amplitude by po vertical movements.On the tectonic level, these zones are structurally and genetically independent elements.Within each block, the processes of secondary  For the Yuzno-Alexandrovskaya field, it is released from the linear blocks whose most significant a submeridian orientation, with a width of 400 -450 m.It recur all the structure and the fault of direction NW under an angle of 30˚ (see Figure

Conclusion
The carbonated intersaliferous reservoirs of Eliets and Zadon horizons consist of carbonates which had undergone catagenetic stage of transformations (dissolution, recrystallization).The irregular distribution of zones of reservoirs rocks characterized by a secondary porosity led to the creation of great spatial (vertical and horizontal) macroheterogeneity.The autocorrelation fun spectral analysis used to evaluate the n logs of the productive interval of the ca saliferous complex, have highlighted the existence of cyclicity in the sedimentation and secondary transformations of the reservoirs rocks.Taking into account the periods of various cycles, and the repartition of petrophysical characteristics of the reservoirs, it appears, on the studied fields, some isometric or linear blocks with common hydrodynamic regime.During controls of interferential continuity between the wells on the fields Yuzno-Sosnovskaya and Yuzno-Alexandrovskaya, the best interactions were observed between the wells of the same block.Consequently, to maintain the layers pressure during exploitation of the oilfields, we recommend consider the blocks which come into view for the implantation of water injection wells.
have a different type II is localise ri re d in South-East (wells 80, 60, 70, 79 and 59) and the types III and IV gather the cuts located in the West and North (zone of the wells P38, 53, 54, 33, 55, 105, 61, 63, 90, 65 and 67).The analysis of the charts of the R o parameter and its trend shows a variation of the ray of correlation (R o ) between 4 and 16 m.The highest values are observed in South-East and the smallest ones are observed in the East of the field.The charts of maximum (T max ) period are the same with the ray (R o ).The minimum v max st of the field.The maximum values (T max > 75 m) are localised in the South-eastern and central parts (well 81)
The zones with minimal values cover the North-West and are characterized by a southern most orientation.The zones of maximum values cover the central and southeastern parts of the field and are characterized by a north-western orientation.The zones of low values occupy the periphery of the structure.The charts of R o trend give possibility to subdivide the field in two (2) different parts: the western and the north-western peripheral parts which are characterized by a soft variation of R o .In each part the minimal values are in the form of submeridian linear zones.The central part of the field is dominated by the maximum values of R o with isometric distribution in the zone of wells 127, 126, 103, 125, 118.The analysis of the T max charts shows two zones entirely different from shape and curves of the same values.In the West, a zone of distribution West-east and another in the central and southern parts of field are characterized by high values of period T max (Figure6).The meeting of these announced zones is done in the area of the wells 117, 128, and 114 where the cuts present a minimal value of T max and more subdivision of the natural reservoir.It is a west-east directed linear zone.o period equivalent to the total thickness of the productive complex occupies the Western and Eastern slopes of the anticline.The 2nd and 3rd types describe a rhythmic succession of reservoirs layers and false covers with average (for the 2nd type) and high (for the 3rd type) frequencies.They respectively occupy the flanks and the heart of the anticline.

Figure 7 .
Figure 7. Reservoirs macroheterogeneity chart within Yuzno-Alexandrovskaya field.1: Site and well number; 2: Faults; 3: Isopacs of measured parameter.transformationsproceeded in an identical way.At the same time, they are different in intensity and character from a block to another.For Davidovskaya field (Figure8), a structuring in linear blocks with 3 principal zones.The Eastern zone is subdivided from the point of view of the alternation of reservoirs layers and false covers.The central zone is characterized by relative increase thicknesses of reservoirs layers and false covers.o the heart of e organization eservoirs layers and false covers.It co