CN113466957A - Reservoir body stratifying method and device based on astronomical gyrus - Google Patents

Reservoir body stratifying method and device based on astronomical gyrus Download PDF

Info

Publication number
CN113466957A
CN113466957A CN202110789813.3A CN202110789813A CN113466957A CN 113466957 A CN113466957 A CN 113466957A CN 202110789813 A CN202110789813 A CN 202110789813A CN 113466957 A CN113466957 A CN 113466957A
Authority
CN
China
Prior art keywords
eccentricity
reservoir
target block
short
astronomical
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110789813.3A
Other languages
Chinese (zh)
Other versions
CN113466957B (en
Inventor
张建国
姜在兴
张元福
江建辉
申静静
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN202110789813.3A priority Critical patent/CN113466957B/en
Publication of CN113466957A publication Critical patent/CN113466957A/en
Application granted granted Critical
Publication of CN113466957B publication Critical patent/CN113466957B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/04Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a reservoir body stratifying method and a reservoir body stratifying device based on astronomical gyrus, wherein the reservoir body stratifying method based on astronomical gyrus comprises the following steps: determining the long eccentricity, the short eccentricity and the time offset of the astronomical cycle of the target block according to the logging data of the target block; establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset; and under the stratum trellis, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference. The reservoir body stratum-fixing method and device based on the astronomical cycle, provided by the invention, have the advantages that the astronomical cycle is applied to oil-gas exploration for the first time, and the problems of high-precision stratum contrast of deep water reservoir bodies and the division of the internal period cycle of large sets of gravel rock bodies are well solved.

Description

Reservoir body stratifying method and device based on astronomical gyrus
Technical Field
The invention relates to the technical field of petroleum and natural gas exploration, in particular to a reservoir body stratifying method and device based on astronomical gyrus.
Background
At present, the reservoir is mainly subjected to isochronous grid division and horizon determination at home and abroad by using a sequence stratigraphy method. However, the stratigraphic method is based on the change of relative water depth as stratigraphic division (fig. 1 and fig. 2, in fig. 2, the water depth of deep lakes and deep sea is from hundreds of meters to thousands of meters, so that the general scale change (from a few meters to a few tens of meters) of the sea (lake) plane has little influence on the deep water environment and hardly causes the lithology change of the deep water reservoir body; 2) the inner curtain structure of the sand (gravel) rock sleeve reservoir body cannot be largely divided.
This is because deep water environments (such as deep lakes and deep sea in geological history) develop a large number of oil and gas reservoirs, and high-precision stratigraphic division of deep water reservoirs is the basis and precondition for accurate oil and gas exploration. Stratigraphic stratigraphy born in the 80's last century, marks of water depth changes in rock records, as a means of deep water reservoir compartmentalization. However, since the deep water environment itself has a large water depth, it is not sensitive to the water depth change (especially, it is not sensitive to the medium or small scale water depth change), and thus it is difficult for the deep water reservoir lithology change to be caused by the water depth change, and further it is difficult to perform the deep water sequence stratigraphy division with high precision. In addition, the evolution of deep water fine-grained sedimentary rocks is also controlled by other factors than water depth. For example, changes in water reducibility, water salinity, water circulation, organic matter productivity and storage conditions, biological disturbance degree, flood strength and frequency, wind field strength, etc., can affect the deposition and evolution of deepwater reservoirs. Therefore, the water depth is used as the only standard for stratigraphic division, and the method has one-sidedness.
In addition, in the technical field of oil and gas exploration, the internal period subdivision of a large set of conglomerate reservoir bodies is also a difficult point and a weak link. Because no obvious lithology or lithofacies change exists in the large-sleeve glutenite reservoir body, the conventional stratum division method adopts sequence stratigraphy, and the problem of secondary division of the internal period of the large-sleeve glutenite reservoir body cannot be solved (figure 3). Because the dividing basis of the sequence stratigraphy is lithology and lithofacies change, the sequence stratigraphy can divide each set of conglomerate reservoir bodies and shale reservoir bodies according to the lithology and lithofacies change. Each set of conglomerate reservoir bodies are actually formed by multi-period deposition, and because obvious lithologic and lithofacies interfaces do not exist in the set of conglomerate reservoir bodies formed by multi-period deposition, the sequence stratigraphy cannot perform periodic subdivision on the interiors of the conglomerate reservoir bodies. As shown in fig. 3, the accuracy of the stratigraphic divisions is equal to the thickness of a set of conglomerate reservoirs, about 10 meters.
Disclosure of Invention
Aiming at the problems in the prior art, the astronomical convolution-based reservoir body stratifying method and device provided by the invention apply astronomical convolution to oil and gas exploration for the first time, and well solve the problems of high-precision stratum contrast of water storage bodies and division of secondary convolution in large sets of gravel rock bodies.
In a first aspect, the present invention provides an astronomical cycle based reservoir stratifying method comprising:
determining the long eccentricity, the short eccentricity and the time offset of the astronomical cycle of the target block according to the logging data of the target block;
establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset;
and under the stratum trellis, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference.
In one embodiment, the determining long eccentricity, short eccentricity and time offset of astronomical convolution of a target block according to log data of the target block comprises:
performing Milnacular convolution comparison on the target block according to single well logging data to determine each level convolution of the target block;
and extracting the long eccentricity, the short eccentricity and the time difference from the respective levels of the gyrations.
In one embodiment, the creating a stratigraphic framework of the target block based on the long eccentricity, the short eccentricity and the time offset includes:
carrying out stratum calibration according to the age of the anchor point and/or the marker layer to establish a cycle stratum lattice of each well of the target block;
comparing the multiple wells in the target block according to the long eccentricity, the short eccentricity and the time offset under the cycle stratigraphic framework to generate a comparison result;
determining the deposition cycle corresponding to the long eccentricity, the short eccentricity and the offset according to the seismic data of the target block;
and establishing the stratum lattice according to the comparison result and the deposition cycle.
In one embodiment, the determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the age under the stratigraphic framework comprises:
predicting the size of the reservoir based on the long eccentricity;
predicting a vertical position of the reservoir from the long eccentricity and the short eccentricity;
according to the age identification the interior curtain structure and the argillaceous interlayer position of the reservoir body.
In a second aspect, the present invention provides an astronomical convolution based reservoir stratifying device, comprising:
the time scale determining module is used for determining the long eccentricity, the short eccentricity and the time offset of the astronomical cycle of the target block according to the logging data of the target block;
the stratum lattice building module is used for building the stratum lattice of the target block based on the long eccentricity, the short eccentricity and the time offset;
and the reservoir layer fixing module is used for determining the reservoir layer of the target block under the stratum framework according to the long eccentricity, the short eccentricity and the time difference.
In one embodiment, the time scale determination module comprises:
the cycle determining unit is used for carrying out Milnackov cycle comparison on the target block according to single well logging data so as to determine each level cycle of the target block;
and the time scale extraction unit is used for extracting the long eccentricity, the short eccentricity and the time difference from each level of convolution.
In one embodiment, the stratigraphic grid building module comprises:
the stratum calibration unit is used for calibrating the stratum according to the age of the anchor point and/or the mark layer so as to establish a gyrating stratum lattice of each well of the target block;
the comparison result generating unit is used for comparing the multiple wells in the target block according to the long eccentricity, the short eccentricity and the time difference under the cycle stratigraphic framework to generate a comparison result;
the deposition cycle determining unit is used for respectively determining deposition cycles corresponding to the long eccentricity, the short eccentricity and the years according to the seismic data of the target block;
and the stratigraphic framework establishing unit is used for establishing the stratigraphic framework according to the comparison result and the deposition cycle.
In one embodiment, the reservoir stratifying module comprises:
a reservoir size prediction unit for predicting the size of the reservoir based on the long eccentricity;
a reservoir vertical position prediction unit for predicting a vertical position of the reservoir from the long eccentricity and the short eccentricity;
interior curtain structure recognition unit is used for according to the age discernment the interior curtain structure and the mud interlayer position of reservoir body.
In a third aspect, the present invention provides an electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor when executing the program implementing the steps of the astronomical convolution based reservoir stratifying method.
In a fourth aspect, the present invention provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of an astronomical convolution based reservoir stratification method.
As can be seen from the above description, the astronomical cycle-based reservoir stratifying method and apparatus provided by the embodiments of the present invention determine the long eccentricity, short eccentricity and offset of the astronomical cycle of the target block according to the well log data of the target block; secondly, establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset; and finally, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference under the stratum framework. According to the invention, the vertical position and the inner tenths of the reservoir are predicted from astronomical geological parameters (long eccentricity, short eccentricity and time difference), so that the divided inner tenths are constrained by geological parameters and have high reliability.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a schematic diagram of the relative positions of sea level, relative sea level and water depth in the world according to the background art of the present invention;
FIG. 2 is a schematic diagram of a deep water environment in the background art of the present invention;
FIG. 3 is a schematic diagram of reservoir period subdivision based on sequence stratigraphy in the background of the present invention;
FIG. 4 is a schematic flow diagram of an astronomical convolution based reservoir stratifying method in an embodiment of the present invention;
FIG. 5 is a schematic view of an astronomical convolution deposition coupling in an embodiment of the present invention;
FIG. 6 is a flowchart illustrating a step 100 according to an embodiment of the present invention;
FIG. 7 is a flowchart illustrating step 200 according to an embodiment of the present invention;
FIG. 8 is a flowchart illustrating step 300 according to an embodiment of the present invention;
FIG. 9 is a schematic flow diagram of an astronomical convolution based reservoir stratifying method in a specific application example of the present invention;
FIG. 10 is a schematic diagram of the prediction of a high quality reservoir using long eccentricity in a specific application example of the present invention;
FIG. 11 is a schematic diagram of the prediction of a high quality reservoir using short eccentricity in a specific application example of the present invention;
FIG. 12 is a diagram showing the prediction of a high-quality reservoir using the age in a specific application example of the present invention;
FIG. 13 is a flow chart of a Chinese gyrus premium reservoir stratifying technique according to an example embodiment of the present invention;
FIG. 14 is a display of the results of an astronomical convolution based reservoir stratifying method in a specific application example of the present invention;
FIG. 15 is a schematic structural diagram of an astronomical convolution based reservoir stratifying device in an embodiment of the present invention;
FIG. 16 is a schematic diagram of the structure of the time scale determination module 10 in the astronomical convolution based reservoir stratifying device in an embodiment of the present invention;
FIG. 17 is a schematic diagram illustrating the structure of the stratigraphic framework building blocks 20 in the astronomical convolution based reservoir stratifying apparatus according to an embodiment of the present invention;
FIG. 18 is a schematic diagram of the module 30 composition of reservoir zonation in an astronomical convolution based reservoir zonation apparatus in an embodiment of the present invention;
fig. 19 is a schematic structural diagram of an electronic device in an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
It should be noted that the terms "comprises" and "comprising," and any variations thereof, in the description and claims of this application and the above-described drawings, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the embodiments with reference to the attached drawings.
Embodiments of the present invention provide an embodiment of an astronomical cycle based reservoir stratifying method, see fig. 4, which specifically includes the following:
step 100: and determining the long eccentricity, the short eccentricity and the time offset of the astronomical cycle of the target block according to the logging data of the target block.
In the prior art, the vertical reservoir stratigraphy method generally depends on sequence stratigraphy: the stratum division of the sequence stratigraphy is carried out according to the water depth change, namely, the stratum division is carried out on the basis of the theoretical basis of controlling the lithology and lithofacies change based on the water depth change. This is applicable to a shallow water environment, because the shallow water environment (water depth of several meters to ten and several meters) has a small water depth, and is therefore sensitive to water depth changes, and meter-level water depth changes can completely cause changes in lithology and lithofacies, so that these lithology and lithofacies changed interfaces can be used as stratigraphic division interfaces, and further, the distribution of reservoirs is constrained in a sequence stratigraphic framework.
As global oil and gas exploration enters a deep water environment from a shallow water environment in a geological historical period, the success or failure of the oil and gas exploration is directly determined by the high-precision vertical layer-fixing of a deep water reservoir body. For deep water reservoirs, the variation of water depth on a general scale has little effect on the depth of the reservoir itself, which is large. This results in the characteristic of large-scale homogeneity, insignificant lithology and lithofacies changes, and difficulty in achieving high-precision stratigraphic comparison with sequence stratigraphy.
On the other hand, in sedimentary rocks, there are constantly developing thick-bed conglomerate reservoirs that may be formed from multiple stages of deposition. However, because the internal lithology changes insignificantly, it is difficult to achieve subtyping of the inside tenths of a conglomerate reservoir using sequence stratigraphy.
It is understood that in step 100, referring to fig. 5, astronomical cycle controls the solar energy cycle, solar excess cycle controls the climate cycle, climate cycle controls glacier cycle, seasonal cycle, evaporation cycle, nutrient input cycle, atmospheric cycle, productivity cycle, horizontal plane cycle, it should be noted that the long eccentricity in step 100 means a cycle of 40.5 ten thousand years, the short eccentricity means a cycle of 10 ten thousand years, and the time difference means a cycle of 1.9 ten thousand years.
Step 200: establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset;
the stratum lattice is a regional space-time ordered arrangement pattern widely applied to various stratums or rock units in stratum sequences. Refers to the geometry of stratigraphic and lithologic units in the basin and their configuration (conybearec.e.b, 1979), which is a three-dimensional concept. The progress of the seismic exploration technology and the occurrence of the sequence stratigraphic method enable non-integration discontinuities and corresponding integration surfaces thereof to be rapidly identified in basin research, sequence stratigraphic units of different levels to be divided and compared, and an isochronous stratigraphic framework to be established. On the basis, the sedimentary system domain, the sedimentary system and the phase can be further researched, and the ancient geographic environment of the basin in each period and the distribution of the sedimentary system can be rebuilt. In contrast to general paleogeographic analysis, in hydrocarbon-bearing basin analysis, megasequence (megasequence), supersequence (sequence), sequence (sequence), system domain, depositional system, and facies are all considered as a kind of geologic body — namely, building units (building blocks) filling different levels of depositional basins. In the oil-gas-containing basin, the research can clarify the spatial configuration relationship of the source rock, the reservoir and the cover.
Step 300: under the stratum framework, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference;
specifically, vertical accurate stratification and inner curtain structure division is realized according to the evolution rule of the high-quality reservoir body in a three-dimensional high-accuracy astronomical cycle framework.
As can be seen from the above description, in the astronomical convolution-based reservoir stratifying method provided by the embodiment of the present invention, the long eccentricity, the short eccentricity and the time difference of the astronomical convolution of the target block are determined according to the well log data of the target block; secondly, establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset; and finally, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference under the stratum framework. The invention applies astronomical convolution to oil-gas exploration for the first time, and solves the problems of high-precision stratum contrast of deep water reservoir bodies and division of the secondary convolution of the internal period of large sets of gravel rock bodies.
In one embodiment, referring to fig. 6, step 100 further comprises:
step 101: performing Milnacular convolution comparison on the target block according to single well logging data to determine each level convolution of the target block;
the comparison of the yland vickers gyrus of the milanand is also called yland vickers gyrus analysis, which is a theory of climate change related to astronomical factors. It is believed that the climate on earth depends on the solar radiation energy received at the surface of the earth. The variation of solar radiation energy is related to the periodic variation of three parameters of orbital eccentricity, ecliptic inclination and offset of the earth around the sun. The alternation of ice periods and inter-ice periods on the earth is the result of this change. The orbit of the earth around the sun is elliptical. The eccentricity of the earth orbit is the ellipticity of the orbit of the earth, and the ellipticity varies periodically between 0.005 and 0.06, and the period of the ellipticity varies by 10 ten thousand years. The modern earth orbit has an ellipse ratio of 0.0167. The variation in eccentricity changes the distance between the earth and the sun, thus affecting the intensity of the solar short wave radiation reaching the earth. When the ellipticity of the orbit of the earth reaches a maximum, the difference of the solar radiation energy between the near-day point and the far-day point reaches 20-30%. At this time, the temperature on the earth is the lowest, and a cold period appears. Conversely, when the ellipse of the earth orbit reaches a minimum, a warm period occurs. Ecliptic slope refers to the inclination of the earth's axis relative to the earth's orbital plane. It also varies periodically from 21.5 ° to 24.5 °. The cycle time was 4.1 ten thousand years. The ecliptic slope of modern earth is 23.5 °. When the ecliptic slope is small, the temperature difference between the earth equator and the two poles is small, the temperature difference between the earth in winter and summer is small, the earth is relatively warm in winter, the air humidity is high, the snow fall amount is increased, the air temperature in summer is low, the ice melting of the two poles is reduced, and the expansion of the ice covers of the two poles is facilitated. And vice versa. The rotation axis of the earth and the vertical axis of the ecliptic plane form an included angle of about 23.5 degrees. The time difference is the time required for the earth to move around the vertical axis of the ecliptic plane for one circle when the earth rotates. The period is 2.3 to 1.9 ten thousand years. According to the theory of the Milankovichi, due to the periodic changes of the earth orbit eccentricity, ecliptic tilt and precision, the radiation energy received by the earth from the sun also changes correspondingly periodically. It is this change in radiant energy that causes alternating changes in temperature and elevation in the earth's climate and periodic elevation of the sea level. The formation of ice stages and inter-ice stages is the common result of these variations in the earth orbit parameters. The resulting sea level change is referred to as glacier sea level change. This variation exists over a period of 10, 4.1 and 2.3 million years.
Step 102: and extracting the long eccentricity, the short eccentricity and the time difference from the respective levels of the gyrations.
In steps 101 and 102, Miridae Virgilla gyrus analysis of curves such as single well GR logging is firstly carried out, gyrus signals of all levels are extracted, three parameters with a duration ratio of 40.5:10:1.9 are selected, and the three parameters comprise long eccentricity, short eccentricity and offset in sequence.
In one embodiment, referring to fig. 7, step 200 further comprises:
step 201: carrying out stratum calibration according to the age of the anchor point and/or the marker layer to establish a cycle stratum lattice of each well of the target block;
step 202: comparing the multiple wells in the target block according to the long eccentricity, the short eccentricity and the time offset under the cycle stratigraphic framework to generate a comparison result;
step 203: determining the deposition cycle corresponding to the long eccentricity, the short eccentricity and the offset according to the seismic data of the target block;
step 204: and establishing the stratum lattice according to the comparison result and the deposition cycle.
In steps 201 to 204, firstly, the anchor point age or the mark layer is used for calibration, and the cycle stratigraphic framework of each well is connected, so that each long eccentricity, short eccentricity and time cycle can be compared among a plurality of wells; and then, carrying out high-precision time-frequency processing on the seismic data, screening out the deposition cycle corresponding to each astronomical period, and establishing a three-dimensional high-precision astronomical cycle framework by means of the horizontal and longitudinal continuous distribution attributes of the seismic data. Finally, high-precision isochronous stratigraphic comparison between fine-grained and medium-grained sediment of the whole basin is achieved.
In one embodiment, referring to fig. 8, step 300 further comprises:
step 301: predicting the size of the reservoir based on the long eccentricity;
step 302: predicting a vertical position of the reservoir from the long eccentricity and the short eccentricity;
step 303: according to the age identification the interior curtain structure and the argillaceous interlayer position of the reservoir body.
In steps 301-303, the strong land-source debris and current input background in the middle of the long-eccentricity gyrus corresponds to a large set of sand (gravel) rock reservoirs, while the top and bottom of the long-eccentricity gyrus are thin interbed sand (gravel) rock reservoirs under weak land-source input and strong wind field conditions, so the long-eccentricity gyrus can be used to predict reservoir size.
The input of the land source in the middle of the short eccentricity ratio is relatively strong, mainly takes the deposition of sand (gravel) as the main stage, and is the development stage of a high-quality reservoir body; the short eccentricity takes argillaceous deposition as the main part, and the reservoir body does not develop relatively, so that the vertical position of the reservoir body can be predicted by utilizing the secondary cycle 'short eccentricity' under the long eccentricity.
The middle part of the time difference is sand (conglomerate) rock with a faster deposition rate, while the bottom part of the time difference is a deposition 'intermission' period, and only a small amount of argillaceous interlayer is deposited; one time cycle is exactly corresponding to one interior time cycle, and then utilizes the time top-bottom interface to realize the recognition of reservoir body interior structure and argillaceous interlayer, so can utilize the time to discern reservoir body interior structure and argillaceous interlayer position.
It will be appreciated that astronomical gyrus is not a novelty, and that astronomical gyrus analysis methods exist in the prior art, but are primarily used in ancient climate studies to determine the age of the formation. The invention introduces astronomical convolution for the first time and applies the astronomical convolution to oil and gas exploration. The accurate horizon of the high-quality reservoir is determined by utilizing the long eccentricity, the short eccentricity and the time-lapse gyrus for the first time.
To further illustrate the present solution, the present invention also provides a specific application example of the astronomical convolution-based reservoir stratifying method, which specifically includes the following, see fig. 9.
The existing stratigraphic technology of the sequence relies on lithology and lithofacies change caused by water depth change excessively, which leads to deep water reservoir bodies insensitive to water depth change and thick-layer conglomerate reservoir bodies with inconspicuous internal lithology change, and the stratigraphic technology of the sequence is difficult to divide the deep water reservoir bodies and restricts oil and gas exploration.
In addition, besides sequence stratigraphy, a few scholars propose to solve the difficult problem of the interior-canopy subdivision of the thick-layer reservoir by utilizing geophysical inversion technology, such as wavelet change analysis. However, the wavelet transformation technology lacks geological theory as a guide and belongs to a purely mathematical digital game. Whether the geological parameters of the corresponding level are actually used for controlling the wavelet transform is verified in a non-subordinate mode, namely, the result is verified in a non-subordinate mode in a wrong mode, and therefore the reliability of the wavelet transform technology is not high. However, the invention applies astronomical convolution to oil and gas exploration for the first time, and solves the problems of high-precision stratum contrast of deep water reservoir bodies and division of the secondary convolution of the internal period of large sets of gravel rock bodies.
The astronomical cycle based reservoir stratifying method provided by this specific application example is illustrated here by way of a lump 167 well.
S1: and establishing a high-precision astronomical convolution grid.
Specifically, a high-precision stratum framework with three parameters of long eccentricity (cycle 40.5 ten thousand years), short eccentricity (cycle 10 ten thousand years) and time difference (cycle 1.9 ten thousand years) of deep and shallow water astronomy cycle is established. The method specifically comprises 3 steps:
carrying out Miridae Virgiz cycle analysis of curves such as single well GR logging and the like, extracting cycle signals of all levels, and selecting three parameters with a duration ratio of 40.5:10:1.9, namely a long eccentricity, a short eccentricity and a time difference in sequence;
secondly, calibrating by utilizing the age of the anchor point or the marker layer, and connecting the gyrating stratum grillage of each well, so that each long eccentricity, short eccentricity and precession can be compared among multiple wells;
and thirdly, high-precision time-frequency processing is carried out on the seismic data, the deposition gyrus corresponding to each astronomical period are screened out, and a three-dimensional high-precision astronomical gyrus framework is built by combining well and earthquake with the aid of the property of transverse and longitudinal continuous distribution of the seismic data.
Through the 3 steps, high-precision isochronous stratigraphic comparison between fine grain deposition and medium coarse grain deposition of the whole basin is achieved.
S2: a premium reservoir is identified under astronomical cyclic trellis control.
And realizing vertical accurate layer-fixing and inner curtain structure division according to the evolution rule of the high-quality reservoir body in a three-dimensional high-accuracy astronomical cycle grid. The method comprises three dimensions from large to small:
predicting reservoir size using long eccentricity convolution. The strong land source debris and current input background at the middle of the long eccentricity gyre corresponds to a large set of sand (gravel) rock reservoirs, while the top and bottom of the long eccentricity gyre are thin interbed sand (gravel) rock reservoirs at weak land source input and strong wind field conditions (fig. 10).
In fig. 10, the premium reservoir is located in the middle of a long eccentricity (the ashen white circle filled offshore underwater fan in the figure is the premium reservoir) within one long eccentricity.
And secondly, predicting the vertical position of the reservoir body by using the short eccentricity of the secondary cycle under the long eccentricity. The input of the land source in the middle of the short eccentricity ratio is relatively strong, mainly takes the deposition of sand (gravel) as the main stage, and is the development stage of a high-quality reservoir body; while the short eccentricity top and bottom predominate in argillaceous deposits, the reservoir relatively does not develop (fig. 11).
In fig. 11, a long eccentricity consists of 4 short eccentricities, and the premium reservoir develops within 2 short eccentricities with the long eccentricity revolving back to the middle. For these two short eccentricities, the premium reservoir develops further in the middle of each short eccentricity revolution (the ashen white circle filled offshore underwater fan in fig. 11 is the premium reservoir).
Identifying the positions of the inner curtain structure and the argillaceous interlayer of the reservoir body by using the time difference. The middle part of the time difference is sand (conglomerate) rock with a faster deposition rate, while the bottom part of the time difference is a deposition 'intermission' period, and only a small amount of argillaceous interlayer is deposited; one offset cycle corresponds exactly to one insider period, and reservoir inside structure and argillaceous interlayer identification is realized by utilizing the offset apical-basal interface (fig. 12 and 13).
In fig. 12, one short eccentricity includes 5 years, and for the high-quality reservoir in the middle of the short eccentricity gyrus, the internal stages of the high-quality reservoir can be divided using the gyrus. As shown in fig. 12, two near-shore underwater fan reservoirs in the lump 167 wells were marked out 3 time laps, each time lap being the next to 1 introversion.
From the above description, it can be seen that the reservoir body stratifying method based on astronomical recursion provided by the embodiment of the present invention predicts the vertical position and the inside-canopy period of the reservoir body from astronomical geological parameters (long eccentricity, short eccentricity, and time difference), so that the divided inside-canopy period is constrained by geological parameters and has high reliability.
The method effectively guides drilling deployment and construction. For example, in the victory oil field 2019, the technology is used for deploying the 664 wells at positions which are ten kilometers away from the 167 wells, and the various levels of the 664 wells are predicted according to the various levels of the gyrating depths and the stratum inclination conditions of the 167 wells before drilling. In the shore 664 well drilling process, a large set of high-quality reservoir bodies are drilled in the middle of long eccentricity gyrus (short eccentricity gyrus two and three), small-scale poor reservoir bodies are drilled at the top and bottom of the long eccentricity gyrus (short eccentricity gyrus one and four), the reservoir bodies are vertically positioned in the middle of each short eccentricity gyrus, argillaceous interlayers are drilled at the top and bottom interfaces of the reservoir bodies, the error between the actual drilling condition and the prediction result is smaller than the one-time gyrus thickness (2-3 meters) (fig. 14, the drilling proves that the vertical prediction error of the high-quality reservoir bodies is reduced from 10 meters to 2.5 meters, the prediction precision is improved by 4 times), and the drilling rate of the high-quality reservoir bodies is effectively improved.
As can be seen from the above description, in the astronomical convolution-based reservoir stratifying method provided by the embodiment of the present invention, the long eccentricity, the short eccentricity and the time difference of the astronomical convolution of the target block are determined according to the well log data of the target block; secondly, establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset; and finally, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference under the stratum framework. According to the invention, the vertical position and the inner tenths of the reservoir are predicted from astronomical geological parameters (long eccentricity, short eccentricity and time difference), so that the divided inner tenths are constrained by geological parameters and have high reliability.
Based on the same inventive concept, the present application also provides an astronomical cycle based reservoir stratifying device, which can be used for implementing the methods described in the above embodiments, as described in the following embodiments. Since the principle of solving the problem of the astronomical cycle based reservoir stratifying device is similar to that of the astronomical cycle based reservoir stratifying method, the implementation of the astronomical cycle based reservoir stratifying device can be referred to the implementation of the astronomical cycle based reservoir stratifying method, and repeated parts are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. While the system described in the embodiments below is preferably implemented in software, implementations in hardware, or a combination of software and hardware are also possible and contemplated.
Embodiments of the present invention provide a specific implementation of an astronomical loop based reservoir stratifying apparatus capable of implementing an astronomical loop based reservoir stratifying method, which specifically includes the following contents, see fig. 15:
the time scale determining module 10 is used for determining the long eccentricity, the short eccentricity and the time offset of the astronomical cycle of the target block according to the logging data of the target block;
a stratigraphic framework building module 20 for building stratigraphic framework of the target block based on the long eccentricity, the short eccentricity and the time offset;
and a reservoir stratifying module 30 for determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the age under the stratum trellis.
In one embodiment, referring to fig. 16, the time scale determining module 10 includes:
a convolution determining unit 101, configured to perform a milnaci convolution comparison on the target block according to single-well logging data, so as to determine each level convolution of the target block;
a time scale extracting unit 102, configured to extract the long eccentricity, the short eccentricity, and the time offset from the respective level revolutions.
In one embodiment, referring to fig. 17, the stratigraphic grid building module 20 comprises:
the stratum calibration unit 201 is configured to perform stratum calibration according to the age of the anchor point and/or the marker layer to establish a rotation stratum lattice of each well of the target block;
a comparison result generating unit 202, configured to compare the multiple wells in the target block according to the long eccentricity, the short eccentricity, and the time offset under the recurrent stratigraphic framework to generate a comparison result;
a deposition cycle determining unit 203, configured to determine deposition cycles corresponding to the long eccentricity, the short eccentricity, and the age according to the seismic data of the target block;
and the stratigraphic framework establishing unit 204 is used for establishing the stratigraphic framework according to the comparison result and the deposition cycle.
In one embodiment, referring to fig. 18, the reservoir zonal module 30 comprises:
a reservoir scale prediction unit 301 for predicting the scale of the reservoir from the long eccentricity;
a reservoir vertical position prediction unit 302 for predicting a vertical position of the reservoir from the long eccentricity and the short eccentricity;
and the inner curtain structure identification unit 303 is used for identifying the inner curtain structure and the argillaceous interlayer position of the reservoir body according to the years.
As can be seen from the above description, the reservoir stratifying device based on astronomical gyrus provided by the embodiment of the present invention firstly determines the long eccentricity, the short eccentricity and the age of the astronomical gyrus of the target block according to the well logging data of the target block; secondly, establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset; and finally, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference under the stratum framework. According to the invention, the vertical position and the inner tenths of the reservoir are predicted from astronomical geological parameters (long eccentricity, short eccentricity and time difference), so that the divided inner tenths are constrained by geological parameters and have high reliability.
Embodiments of the present application further provide a specific implementation of an electronic device capable of implementing all steps in the astronomical convolution-based reservoir stratifying method in the foregoing embodiments, and referring to fig. 19, the electronic device specifically includes the following contents:
a processor (processor)1201, a memory (memory)1202, a communication Interface 1203, and a bus 1204;
the processor 1201, the memory 1202 and the communication interface 1203 complete communication with each other through the bus 1204; the communication interface 1203 is configured to implement information transmission between related devices, such as a server-side device, a measurement device, and a client device.
The processor 1201 is used to invoke a computer program in the memory 1202, the processor when executing the computer program implementing all the steps in the astronomical convolution based reservoir stratifying method of the above embodiments, e.g. the processor when executing the computer program implementing the following steps:
step 100: determining the long eccentricity, the short eccentricity and the time offset of the astronomical cycle of the target block according to the logging data of the target block;
step 200: establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset;
step 300: and under the stratum trellis, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference.
Embodiments of the present application also provide a computer-readable storage medium capable of implementing all the steps in the astronomical convolution based reservoir stratifying method of the above embodiments, the computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements all the steps of the astronomical convolution based reservoir stratifying method of the above embodiments, for example, the processor implements the following steps when executing the computer program:
step 100: determining the long eccentricity, the short eccentricity and the time offset of the astronomical cycle of the target block according to the logging data of the target block;
step 200: establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset;
step 300: and under the stratum trellis, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, for the hardware + program class embodiment, since it is substantially similar to the method embodiment, the description is simple, and the relevant points can be referred to the partial description of the method embodiment.
The foregoing description has been directed to specific embodiments of this disclosure. Other embodiments are within the scope of the following claims. In some cases, the actions or steps recited in the claims may be performed in a different order than in the embodiments and still achieve desirable results. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing may also be possible or may be advantageous.
Although embodiments of the present description provide method steps as described in embodiments or flowcharts, more or fewer steps may be included based on conventional or non-inventive means. The order of steps recited in the embodiments is merely one manner of performing the steps in a multitude of orders and does not represent the only order of execution. When an actual apparatus or end product executes, it may execute sequentially or in parallel (e.g., parallel processors or multi-threaded environments, or even distributed data processing environments) according to the method shown in the embodiment or the figures. The terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, the presence of additional identical or equivalent elements in a process, method, article, or apparatus that comprises the recited elements is not excluded.
The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.

Claims (10)

1. An astronomical convolution-based reservoir stratifying method, comprising:
determining the long eccentricity, the short eccentricity and the time offset of the astronomical cycle of the target block according to the logging data of the target block;
establishing a stratum framework of the target block based on the long eccentricity, the short eccentricity and the time offset;
and under the stratum trellis, determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the time difference.
2. The astronomical convolution-based reservoir stratifying method of claim 1, wherein the determining long eccentricity, short eccentricity and years of astronomical convolution of a target block from log data of the target block comprises:
performing Milnacular convolution comparison on the target block according to single well logging data to determine each level convolution of the target block;
and extracting the long eccentricity, the short eccentricity and the time difference from the respective levels of the gyrations.
3. The astronomical convolution-based reservoir stratifying method of claim 2, wherein the establishing a stratigraphic framework of the target block based on the long eccentricity, short eccentricity and years comprises:
carrying out stratum calibration according to the age of the anchor point and/or the marker layer to establish a cycle stratum lattice of each well of the target block;
comparing the multiple wells in the target block according to the long eccentricity, the short eccentricity and the time offset under the cycle stratigraphic framework to generate a comparison result;
determining the deposition cycle corresponding to the long eccentricity, the short eccentricity and the offset according to the seismic data of the target block;
and establishing the stratum lattice according to the comparison result and the deposition cycle.
4. The astronomical convolution-based reservoir stratifying method of claim 1, wherein the determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the years under the stratigraphic framework comprises:
predicting the size of the reservoir based on the long eccentricity;
predicting a vertical position of the reservoir from the long eccentricity and the short eccentricity;
according to the age identification the interior curtain structure and the argillaceous interlayer position of the reservoir body.
5. An astronomical convolution-based reservoir stratifying device, comprising:
the time scale determining module is used for determining the long eccentricity, the short eccentricity and the time offset of the astronomical cycle of the target block according to the logging data of the target block;
the stratum lattice building module is used for building the stratum lattice of the target block based on the long eccentricity, the short eccentricity and the time offset;
and the reservoir layer fixing module is used for determining the reservoir layer of the target block under the stratum framework according to the long eccentricity, the short eccentricity and the time difference.
6. The astronomical convolution-based reservoir stratifying device of claim 5, wherein the time scale determination module comprises:
the cycle determining unit is used for carrying out Milnackov cycle comparison on the target block according to single well logging data so as to determine each level cycle of the target block;
and the time scale extraction unit is used for extracting the long eccentricity, the short eccentricity and the time difference from each level of convolution.
7. The astronomical convolution-based reservoir stratifying device of claim 6, wherein the stratigraphic framework building module comprises:
the stratum calibration unit is used for calibrating the stratum according to the age of the anchor point and/or the mark layer so as to establish a gyrating stratum lattice of each well of the target block;
the comparison result generating unit is used for comparing the multiple wells in the target block according to the long eccentricity, the short eccentricity and the time difference under the cycle stratigraphic framework to generate a comparison result;
the deposition cycle determining unit is used for respectively determining deposition cycles corresponding to the long eccentricity, the short eccentricity and the years according to the seismic data of the target block;
and the stratigraphic framework establishing unit is used for establishing the stratigraphic framework according to the comparison result and the deposition cycle.
8. The astronomical convolution-based reservoir stratifying device of claim 5, wherein the determining the reservoir of the target block according to the long eccentricity, the short eccentricity and the years under the stratigraphic framework comprises:
a reservoir size prediction unit for predicting the size of the reservoir based on the long eccentricity;
a reservoir vertical position prediction unit for predicting a vertical position of the reservoir from the long eccentricity and the short eccentricity;
interior curtain structure recognition unit is used for according to the age discernment the interior curtain structure and the mud interlayer position of reservoir body.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the steps of the astronomical convolution based reservoir stratification method of any of the claims 1 to 4.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the astronomical convolution based reservoir stratification method of any of the claims 1 to 4.
CN202110789813.3A 2021-07-13 2021-07-13 Reservoir body stratifying method and device based on astronomical gyrus Active CN113466957B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110789813.3A CN113466957B (en) 2021-07-13 2021-07-13 Reservoir body stratifying method and device based on astronomical gyrus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110789813.3A CN113466957B (en) 2021-07-13 2021-07-13 Reservoir body stratifying method and device based on astronomical gyrus

Publications (2)

Publication Number Publication Date
CN113466957A true CN113466957A (en) 2021-10-01
CN113466957B CN113466957B (en) 2022-07-01

Family

ID=77880009

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110789813.3A Active CN113466957B (en) 2021-07-13 2021-07-13 Reservoir body stratifying method and device based on astronomical gyrus

Country Status (1)

Country Link
CN (1) CN113466957B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115217470A (en) * 2022-07-19 2022-10-21 中国石油大学(华东) Centimeter-micron scale cycloidal division and driving factor identification method in shale

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109001840A (en) * 2018-07-18 2018-12-14 中国地质大学(北京) A kind of mud shale genetic analysis method and device
CN111381277A (en) * 2018-12-27 2020-07-07 中国石油天然气股份有限公司 Method and device for establishing carbonate rock isochronous stratum framework
US20200393586A1 (en) * 2019-06-12 2020-12-17 John A. BABCOCK Lowstand erosional seismic stratigraphy

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109001840A (en) * 2018-07-18 2018-12-14 中国地质大学(北京) A kind of mud shale genetic analysis method and device
CN111381277A (en) * 2018-12-27 2020-07-07 中国石油天然气股份有限公司 Method and device for establishing carbonate rock isochronous stratum framework
US20200393586A1 (en) * 2019-06-12 2020-12-17 John A. BABCOCK Lowstand erosional seismic stratigraphy

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
闫建平: "湖相泥页岩天文地层旋回测井识别在沾化凹陷沙三下亚段的应用", 《测井技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115217470A (en) * 2022-07-19 2022-10-21 中国石油大学(华东) Centimeter-micron scale cycloidal division and driving factor identification method in shale

Also Published As

Publication number Publication date
CN113466957B (en) 2022-07-01

Similar Documents

Publication Publication Date Title
Feng et al. Sequence stratigraphy and importance of syndepositional structural slope-break for architecture of Paleogene syn-rift lacustrine strata, Bohai Bay Basin, E. China
CN104765065B (en) Method and device for determining sedimentary stratum ancient landform
Sullivan et al. Application of new seismic attributes to collapse chimneys in the Fort Worth Basin
Grammer et al. Integration of outcrop and modern analogs in reservoir modeling: Overview with examples from the Bahamas
Sayago et al. Characterization of a deeply buried paleokarst terrain in the Loppa High using core data and multiattribute seismic facies classification
CN107589469B (en) The determination method and apparatus of oil-water interfaces
CN105527653A (en) Virtual well construction method based on geological information
Yin et al. A three-dimensional high-resolution reservoir model of the Eocene Shahejie Formation in Bohai Bay Basin, integrating stratigraphic forward modeling and geostatistics
CN113466957B (en) Reservoir body stratifying method and device based on astronomical gyrus
Sinha et al. Real-Time 3D imaging of complex turbiditic reservoir architecture
Qi et al. Hierarchy and subsurface correlation of muddy baffles in lacustrine delta fronts: a case study in the X Oilfield, Subei Basin, China
CN110471107A (en) The depicting method and device of the disconnected solution of carbonate rock
CN110320574A (en) The method portrayed based on gentle slope delta Thin Sandbody
Yang et al. Gas accumulation conditions and key technologies for exploration & development of Sinian and Cambrian gas reservoirs in Anyue gasfield
Bond et al. Evidence of astronomical forcing of the Earth's climate in Cretaceous and Cambrian times
Naseer Application of instantaneous spectral decomposition-based porosity simulations for imaging shallow-marine stratigraphic traps of Lower-Eocene carbonates sequences of Indus Basin, Onshore Pakistan
CN107507095A (en) From fast conversion method from depth of stratum section to stratum time section and system
Naseer Seismic attributes and quantitative reverse reservoir simulations of shallow-marine eocene sedimentary fairways, extensional rifting Indus onshore petroleum system, southern Pakistan: Implications for hydrocarbon exploration
Rey et al. Stratigraphy: terminology and practice
CN113885096A (en) High-frequency sequence partitioning and small layer comparison method and device, electronic equipment and medium
SZABÓ Pore Pressure Prediction In Pannonian Hydrocarbon Reservoir Systems Using An Integrated Interpretation Approach
Wu et al. Characteristics of Milankovitch cycles recorded in Eocene strata in the eastern depression of North Yellow Sea Basin, North China
CN116203625A (en) Integrated prediction method and device for conventional and unconventional oil and gas reservoirs
Pipping et al. K18-Golf Field seismic and reservoir modelling challenges
Malin Improved Constraint on Salt Geometry in the Southern Nordkapp Basin-Modeling salt geometry using potential field data integrated with traditional reflection seismic and well data.

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant