CN112347608A - Clastic rock reservoir diagenetic stage characterization method, system, device and storage medium - Google Patents

Clastic rock reservoir diagenetic stage characterization method, system, device and storage medium Download PDF

Info

Publication number
CN112347608A
CN112347608A CN202011004770.5A CN202011004770A CN112347608A CN 112347608 A CN112347608 A CN 112347608A CN 202011004770 A CN202011004770 A CN 202011004770A CN 112347608 A CN112347608 A CN 112347608A
Authority
CN
China
Prior art keywords
temperature
paleo
diagenetic
montmorillonite
illite
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
CN202011004770.5A
Other languages
Chinese (zh)
Other versions
CN112347608B (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.)
Yangtze University
Original Assignee
Yangtze University
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 Yangtze University filed Critical Yangtze University
Priority to CN202011004770.5A priority Critical patent/CN112347608B/en
Publication of CN112347608A publication Critical patent/CN112347608A/en
Application granted granted Critical
Publication of CN112347608B publication Critical patent/CN112347608B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/10Geothermal energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Silicates, Zeolites, And Molecular Sieves (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention relates to a clastic rock reservoir diagenetic stage characterization method, a clastic rock reservoir diagenetic stage characterization system, a clastic rock diagenetic stage characterization device and a storage medium, wherein the clastic rock diagenetic stage characterization method comprises the steps of obtaining a geothermal gradient, a stratum buried depth, a constant temperature zone buried depth and an ancient surface temperature, and obtaining a stratum ancient temperature according to the geothermal gradient, the stratum buried depth, the constant temperature zone buried depth and the ancient surface temperature; obtaining the reflectivity of a source rock vitrinite and the illite/montmorillonite ratio of the illite/montmorillonite mixed layer according to the formation paleo-temperature; and dividing the diagenetic stages of the clastic rock reservoir according to the paleo-temperature of the stratum, the vitrinite reflectivity of the source rock and the montmorillonite ratio of the illite/montmorillonite mixed layer. The method realizes more precise characterization of the diagenetic rock reservoir diagenesis stage and is beneficial to highlighting the difference of the reservoir evolution process.

Description

Clastic rock reservoir diagenetic stage characterization method, system, device and storage medium
Technical Field
The invention relates to the technical field of oil and gas exploration and development, in particular to a clastic rock reservoir diagenetic stage characterization method, system and device and a computer readable storage medium.
Background
In the sedimentary burying and diagenesis processes, a clastic rock diagenesis field, diagenesis substances, diagenesis and hydrocarbon generation capacity are closely related to diagenesis stages, and the clastic rock property is greatly influenced, so that the research on the diagenesis process has important significance on oil and gas exploration and development. In recent years, with the improvement of the oil and gas exploration degree, deep ultra-deep layer oil and gas exploration becomes the most main exploration target in the future, and deep ultra-deep layer reservoirs undergo a more complex diagenetic evolution process in the burying process, so that the research on diagenetic stages is enhanced, and important guidance can be provided for the effectiveness evaluation of the reservoirs in the oil and gas reservoir formation period.
A plurality of scholars conduct extensive and intensive research on diagenesis of clastic rock and put forward corresponding diagenesis stage division standards; in the early 90 s of the last century, domestic scholars need to be in a phoenix-eye and auspicious system to search and collect results of domestic and foreign scholars on diagenesis research and diagenesis stage division, provide clastic rock diagenesis stage division specifications in the petroleum and gas industry, and divide diagenesis of clastic rock into an early diagenesis stage and a late diagenesis stage; in the beginning of the twenty-first century, the early clastic rock diagenesis stage division standard is modified by adopting Fengxiang, clastic rock diagenesis evolution is divided into 4 stages, namely early, middle and late stages, a diagenesis stage is divided into A, B stages, and the A stage of the diagenesis stage is further divided into A1 and A2; menyanglin subdivides the diagenesis phase a of clastic rock into 4 diagenesis micro-phases on the basis of prior manual work.
The existing scheme is not beneficial to highlighting the difference of the reservoir evolution process and identifying the key geological elements causing the reservoir differential evolution.
Disclosure of Invention
In view of the above, there is a need to provide a clastic rock reservoir diagenetic stage characterization method for solving the technical problem that the prior art is not favorable for highlighting the reservoir evolution process difference.
The invention provides a clastic rock reservoir diagenetic stage characterization method, which comprises the following steps:
acquiring a ground temperature gradient, a stratum buried depth, a constant temperature zone buried depth and an ancient surface temperature, and acquiring a stratum ancient temperature according to the ground temperature gradient, the stratum buried depth, the constant temperature zone buried depth and the ancient surface temperature;
obtaining the reflectivity of a source rock vitrinite and the illite/montmorillonite ratio of the illite/montmorillonite mixed layer according to the formation paleo-temperature;
and dividing the diagenetic stages of the clastic rock reservoir according to the paleo-temperature of the stratum, the vitrinite reflectivity of the source rock and the montmorillonite ratio of the illite/montmorillonite mixed layer.
Further, acquiring the formation paleo-temperature according to the geothermal gradient, the formation burial depth, the constant temperature zone burial depth and the paleo-surface temperature, and concretely comprises the following steps of acquiring the paleo-temperature of the formation according to a formula
T=A*(H-H0)+T0
Obtaining the paleotemperature T of the stratum, wherein A is the ground temperature gradient, H is the buried depth of the stratum, and H0For constant temperature zone buried depth, T0Is the paleo-surface temperature.
Further, obtaining the reflectivity of the source rock vitrinite according to the paleo-temperature of the stratum specifically comprises obtaining the reflectivity of the source rock vitrinite according to a formula
Figure BDA0002695520420000021
Obtaining the reflectivity Ro% of the source rock vitrinite, wherein T is the formation paleo-temperature,
Figure BDA0002695520420000022
are coefficients.
Further, obtaining the illite/montmorillonite ratio according to the paleo-temperature of the stratum specifically comprises obtaining the illite/montmorillonite ratio according to a formula
I/S-S%=-3.293e-07*T^4+0.0001625*T^3-0.02466*T^2+0.5809*T+ 97.06
And obtaining the I/S-S% ratio of montmorillonite of the illite/montmorillonite mixed layer, wherein T is the paleotemperature of the stratum.
Further, dividing the clastic rock reservoir diagenetic stages according to the formation paleo-temperature, the source rock vitrinite reflectivity and the illite/montmorillonite ratio, wherein the method specifically comprises the following steps of dividing the formation paleo-temperature in a range of less than or equal to 200 ℃, the source rock vitrinite reflectivity in a range of more than or equal to 4 and the I/S-S% illite/montmorillonite ratio in a range of 0 to 1 into 36 diagenetic stages; and determining the division of the clastic rock reservoir diagenetic stages according to the formation paleo-temperature, the source rock vitrinite reflectivity and the illite/montmorillonite ratio and 36 diagenetic stages.
Further, dividing the paleo-temperature of the stratum in the range of less than or equal to 200 ℃, the vitrinite reflectivity of the source rock in the range of more than or equal to 4 and the I/S-S% illite/montmorillonite ratio in the range of 0 to 1 into 36 diagenesis stages, specifically comprising,
dividing a diagenesis stage in which the paleo-temperature of the stratum is within a range of less than or equal to 25 ℃, the vitrinite reflectance of the source rock is within a range of less than 0.24, and the I/S-S% illite/montmorillonite ratio is greater than 99 and within a range of less than or equal to 100 into a1 st diagenesis stage;
dividing a diagenesis stage into 2 nd to 35 th diagenesis stages, wherein the formation paleo-temperature is more than 25 ℃ and less than or equal to 195 ℃, the source rock vitrinite reflectivity is more than or equal to 0.24 and less than 3.45, and the I/S-S% illite/montmorillonite ratio is less than or equal to 99 and more than 1;
dividing a diagenesis stage into a 36 th diagenesis stage, wherein the formation paleo-temperature is more than 195 and less than or equal to 200 ℃, the source rock vitrinite reflectance is more than or equal to 3.45 and less than or equal to 4, and the I/S-S% illite/montmorillonite ratio is less than or equal to 1 and more than or equal to 0.
Furthermore, the clastic rock reservoir diagenetic stage characterization method further comprises the step of simulating time domain and space domain distribution of clastic rock reservoir diagenesis after the clastic rock reservoir diagenesis stage is divided.
The invention also provides a fine characterization system for the diagenetic stage of the clastic rock reservoir, which comprises a formation paleo-temperature acquisition module, a reflectivity and montmorillonite ratio acquisition module and a diagenetic stage division module;
the stratum paleo-temperature acquisition module is used for acquiring a ground temperature gradient, a stratum buried depth, a constant temperature zone buried depth and an paleo-surface temperature, and acquiring a stratum paleo-temperature according to the ground temperature gradient, the stratum buried depth, the constant temperature zone buried depth and the paleo-surface temperature;
the reflectivity and montmorillonite ratio acquisition module is used for acquiring the reflectivity of a source rock vitrinite and the montmorillonite ratio of an illite/montmorillonite mixed layer according to the formation paleo-temperature;
and the diagenetic stage division module is used for dividing the diagenetic stages of the clastic rock reservoir according to the paleo-temperature of the stratum, the reflectivity of the source rock vitrinite and the montmorillonite ratio of the illite/montmorillonite layer.
The invention also provides a clastic rock reservoir diagenetic phase fine characterization device which comprises a processor and a memory, wherein the memory is stored with a computer program, and when the computer program is executed by the processor, the clastic rock reservoir diagenetic phase characterization method is realized according to any technical scheme.
The invention also provides a computer-readable storage medium, on which a computer program is stored, which, when executed by a processor, implements the method for characterizing the diagenetic phase of a detritus reservoir as defined in any of the above claims.
Compared with the prior art, the invention has the beneficial effects that: acquiring the formation paleo-temperature according to the geothermal gradient, the formation burial depth, the constant temperature zone burial depth and the paleo-surface temperature by acquiring the geothermal gradient, the formation burial depth, the constant temperature zone burial depth and the paleo-surface temperature; obtaining the vitrinite reflectivity of the source rock and the illite/montmorillonite ratio of the illite/montmorillonite mixed layer according to the formation paleo-temperature; dividing the diagenetic stages of the clastic rock reservoir according to the paleo-temperature of the stratum, the vitrinite reflectivity of the source rock and the montmorillonite ratio of the illite/montmorillonite mixed layer; the method realizes more precise characterization of the diagenetic rock reservoir diagenesis stage and is beneficial to highlighting the difference of the reservoir evolution process.
Drawings
Fig. 1 is a schematic flow chart of a clastic rock reservoir diagenetic stage characterization method provided by the invention.
Detailed Description
The preferred embodiments of the present invention will now be described in detail with reference to the accompanying drawings, which form a part hereof, and which together with the embodiments of the invention serve to explain the principles of the invention and not to limit its scope.
Example 1
The embodiment of the invention provides a clastic rock reservoir diagenetic stage characterization method, which is a schematic flow chart shown in figure 1 and comprises the following steps:
s1, acquiring a ground temperature gradient, a stratum buried depth, a constant temperature zone buried depth and an ancient surface temperature, and acquiring a stratum ancient temperature according to the ground temperature gradient, the stratum buried depth, the constant temperature zone buried depth and the ancient surface temperature;
s2, obtaining the reflectivity of the source rock vitrinite and the illite/montmorillonite ratio according to the paleo-temperature of the stratum;
and S3, dividing the clastic rock reservoir diagenetic stages according to the paleo-temperature of the stratum, the vitrinite reflectivity of the source rock and the illite/montmorillonite ratio of the illite/montmorillonite layer.
It should be noted that sandstone diagenesis is a complex reaction process controlled by multiple geological factors, which mainly comprises a geologic structure background, a stratum burial speed, a deposition environment, hydrodynamic conditions, sediment components and structures, a diagenesis fluid environment, a water-rock reaction speed, a geothermal gradient, a stratum abnormal pressure and the like, in the process, diagenesis continuously changes the physical properties of a reservoir, and diagenesis types and physical properties characteristics show a certain rule; the influence of the mechanical compaction on the physical properties of the reservoir is mainly concentrated in the early diagenesis stage (early diagenesis stage A and early diagenesis stage B); the reservoir evolves to a medium diagenesis A1 stage, the rock is gradually consolidated into rock, and the rock has certain compressive resistance and mainly takes the cementing effect as the main point; the reservoir evolves to a medium-formation rock A2 stage, densification begins, mechanical compaction begins to convert to chemical compaction, particles are mostly in line-concave contact and concave-convex contact, and erosion action develops under the action of organic acid; after entering the diagenetic stage B, the compaction mainly takes chemical compaction, but the overall strength is not large, the diagenetic environment is converted to weak alkali, the erosion holes formed in the stage A2 are filled by carbonate cement and clay minerals, and the physical properties are further deteriorated; however, the degree of reservoir diagenesis evolution and the physical property change characteristics can have large differences;
it should be noted that the diagenetic characteristics of the clastic rock are influenced by geological factors such as the original texture of the sediment, temperature, pressure and fluid, wherein the temperature and pressure conditions and the fluid environment are the main control factors. The regular evolution and combination of different geological factors in the period of the geological history form a series of diagenetic stages and diagenetic orders, and simultaneously, various diagenetic parameters such as the reflectivity Ro of vitrinite, the content (I/S-S%) of smectite in an illite/montmorillonite mixed layer and the like are regularly changed, so that the diagenetic history of the sedimentary strata can be comprehensively researched and divided into diagenetic evolution stages by simulating the change rule of the diagenetic parameters along with time;
based on the burying history and the thermal history, and by means of a Sweeney chemical kinetics mode, the Ro% of different historical moments can be calculated. The process of converting the smectite into the illite is mainly influenced by temperature, time, pressure, deposition environment and source characteristics, and the conversion process can be expressed by an Elliot model formula;
preferably, the method for obtaining the formation paleo-temperature according to the geothermal gradient, the formation burial depth, the constant temperature zone burial depth and the paleo-surface temperature specifically comprises the following steps of obtaining the paleo-temperature according to a formula
T=A*(H-H0)+T0
Obtaining the paleotemperature T of the stratum, wherein A is the ground temperature gradient, H is the buried depth of the stratum, and H0For constant temperature zone buried depth, T0Is the paleo-surface temperature.
In specific implementation, the formation paleo-temperatures at different geological times are calculated based on geothermal gradient simulation;
T=A*(H-H0)+T0 (1)
wherein T is the paleo-temperature of the stratum, A is the gradient of the geothermal temperature, H is the buried depth of the stratum, H0For constant temperature with buried depth, T0Is the paleo-surface temperature;
preferably, the obtaining of the reflectivity of the source rock vitrinite according to the paleo-temperature of the stratum specifically comprises obtaining the reflectivity of the source rock vitrinite according to a formula
Figure BDA0002695520420000061
Obtaining the reflectivity Ro% of the source rock vitrinite, wherein T is the formation paleo-temperature,
Figure BDA0002695520420000062
are coefficients.
It should be noted that the reflectivity (Ro%) of the source rock vitrinite in the sedimentary basin is in positive correlation with the temperature, and in the burying process, the higher the formation temperature is, the larger the value is, and the process is irreversible and stable for ten minutes;
the vitrinite reflectivity method is extended to the determination of the maturity of the organic matter dispersed in sedimentary rock, is the most widely applied index of the maturity of the organic matter and can be used as a digital scale; because of the advantage of quantifiable organic matter maturity, the method is most widely applied to the aspect of sedimentary basin heat history research; in some complex temperature models, in order to simplify the establishment of the temperature models, the establishment of some complex temperature models often ignores the influence of time and pressure factors on the temperature; in the initial stage of burying, the value of the reflectivity of the vitrinite is lower, the reflectivity is increased faster along with the gradual increase of the buried depth, and the reflectivity of the vitrinite is continuously increased and the change rate is slowed down till the deep high-temperature gas generation stage, so that the whole change is S-shaped; the chemical structures of different types of kerogen are different, so that the maturation time is different, and when the vitrinite reflectivity is used for judging the maturity of organic matters, the kerogen of different types is different, and the coefficient can be used
Figure BDA0002695520420000063
Carrying out correction;
Figure BDA0002695520420000071
preferably, obtaining the illite/montmorillonite ratio according to the paleo-temperature of the stratum specifically comprises obtaining the illite/montmorillonite ratio according to a formula
I/S-S%=-3.293e-07*T^4+0.0001625*T^3-0.02466*T^2+0.5809*T+ 97.06
And obtaining the I/S-S% ratio of montmorillonite of the illite/montmorillonite mixed layer, wherein T is the paleotemperature of the stratum.
In one embodiment, the clay mineral combination and the transformation condition thereof, including the change of the mixed layer ratio in the mixed layer clay mineral, are important bases for dividing into diagenetic stages, and can lead diagenetic action research to be developed from qualitative to quantitative; the method is divided into rock stages by utilizing the change of clay minerals and the mixing layer ratio thereof, and has wide research and application abroad; the process of conversion of smectite to illite is mainly affected by temperature, time, pressure, deposition environment and source characteristics; this conversion process can be formulated using a model
I/S-S%=-3.293e-07*T^4+0.0001625*T^3-0.02466*T^2+0.5809*T+97.06(3)
The conversion of the montmorillonite is closely related to diagenesis, the deeper the diagenesis is, the higher the conversion rate of the montmorillonite is, the more violent the diagenesis is and the higher the conversion rate of the montmorillonite is; according to the research of diagenetic sequences of clastic rock reservoirs in different basin regions, the results show that the dissolution of feldspar and carbonate and the precipitation of quartz mainly occur in the diagenetic metaphase A1 stage (85-110 ℃) in an acid environment, and the carbonate and clay minerals are easy to be cemented in the early diagenetic stage A, B (<85 ℃) and the diagenetic mesodiagenetic stage A2 (110-140 ℃); the dissolution of carbonate minerals and aluminosilicate minerals is closely related to the formation of hydrocarbons by organic matters, because organic acids are formed by the thermal evolution of the organic matters, so that the property of the aqueous solution is changed, and the organic acids are very stable and can be stored at the temperature of 80-120 ℃, so that the quantity and the proportion of elements entering the solution of the aluminosilicate minerals and the carbonate cements are positively influenced;
the dissolution of feldspar is also influenced by the components, and when the content of plagioclase in sandstone is higher, the dissolution of feldspar is relatively strong, mainly because the content of calcium ions (Ca2+) in plagioclase feldspar is higher, and the plagioclase feldspar is easier to be combined with carbon dioxide;
with the increase of temperature and pressure, the dissolution rate increases sharply and slowly decreases after reaching the maximum value; feldspar dissolution is usually accompanied by precipitation of secondary minerals including kaolinite, illite, and quartz; potassium feldspar, albite and calcite feldspar can be spontaneously converted into kaolinite and illite in the diagenesis;
kaolinite formed by dissolving feldspar is mainly produced in the middle stage B of diagenesis, the cementation rate of illite is increased along with the increase of temperature, and the kaolinite is sharply reduced after the kaolinite reaches the maximum value; the illite cementation rate reaches the maximum value at a medium temperature (130-135 ℃);
macroscopically, the rate of quartz cementation is negligible below 80 ℃, but then exponentially increases with temperature, and oxygen isotope data of quartz overgrowth shows that cementation is easy to occur in the temperature range of 80-150 ℃;
diagenesis (erosion and cementation) is easier to occur in the phase A of the diagenesis stage, and at the moment, the mechanical compaction of sandstone is not compact enough in the early stage, enough residual interparticle pores are reserved, so that the flow of diagenesis fluid and the exchange of diagenesis substances are easy to realize, and mineral erosion and precipitation are more active;
in the diagenetic stage A, the hydrocarbon forming efficiency of organic matters is highest, and the aqueous solution contains enough organic acids, so that the corrosion of aluminosilicate minerals and carbonate minerals is facilitated. The model 1-3 montmorillonite conversion model is based on the theory that the strength of montmorillonite conversion is increased along with the increase of diagenesis stage and is sharply reduced after reaching the middle diagenesis stage A1;
in the early stage of burying, the content of authigenic illite is low, the conversion of montmorillonite is increased rapidly along with the gradual increase of the burying depth, and the conversion amount continues to increase and the change rate is reduced until the deep high-temperature gas generation stage, so that the whole change is in an inverse S shape;
preferably, the clastic rock reservoir diagenetic stages are divided according to the paleo-temperature of the stratum, the vitrinite reflectivity of the source rock and the illite/montmorillonite ratio, and the method specifically comprises the following steps of dividing the paleo-temperature of the stratum into 36 diagenetic stages within the range of less than or equal to 200 ℃, the vitrinite reflectivity of the source rock within the range of more than or equal to 4 and the I/S-S% illite/montmorillonite ratio within the range of 0 to 1; determining the division of the clastic rock reservoir diagenetic stages according to the formation paleo-temperature, the source rock vitrinite reflectivity and the illite/montmorillonite ratio and 36 diagenetic stages;
in a specific embodiment, selecting three parameters of the paleo-geothermal temperature T (DEG C), the vitrinite reflectance (Ro%) and the content (I/S-S%) of smectite in an illite/montmorillonite mixed layer, and combining a diagenesis evolution stage division scheme to divide the diagenesis stage into diagenesis stages;
comprehensively considering the influence of four geological factors of formation paleo-temperature (T), pressure (P), fluid and time (T) on diagenesis, selecting three parameters of formation paleo-temperature (DEG C), vitrinite reflectance (Ro%) and smectite content (I/S-S%) in clay mineral illite/montmorillonite mixed layer as main parameters of diagenesis evolution comprehensive simulation, dividing diagenesis stages, and simulating time domain and space domain distribution of diagenesis stages; dividing diagenesis stages into 36 levels according to three parameters of formation paleo-temperature (T), vitrinite reflectivity (Ro%) and illite/montmorillonite ratio (I/S-S%), wherein the levels are DS1-DS 36;
preferably, the paleo-formation temperature is less than or equal to 200 ℃, the source rock vitrinite reflectance is greater than or equal to 4, and the I/S-S% illite/montmorillonite ratio is in the range of 0 to 1, and the paleo-formation temperature is divided into 36 diagenesis stages, specifically comprising,
dividing a diagenesis stage in which the paleo-temperature of the stratum is within a range of less than or equal to 25 ℃, the vitrinite reflectance of the source rock is within a range of less than 0.24, and the I/S-S% illite/montmorillonite ratio is greater than 99 and within a range of less than or equal to 100 into a1 st diagenesis stage;
dividing a diagenesis stage into 2 nd to 35 th diagenesis stages, wherein the formation paleo-temperature is more than 25 ℃ and less than or equal to 195 ℃, the source rock vitrinite reflectivity is more than or equal to 0.24 and less than 3.45, and the I/S-S% illite/montmorillonite ratio is less than or equal to 99 and more than 1;
dividing a diagenesis stage into a 36 th diagenesis stage, wherein the formation paleo-temperature is more than 195 and less than or equal to 200 ℃, the source rock vitrinite reflectance is more than or equal to 3.45 and less than or equal to 4, and the I/S-S% illite/montmorillonite ratio is less than or equal to 1 and more than or equal to 0;
in one embodiment, the clastic reservoir diagenesis stage is finely divided, as shown in table 1,
TABLE 1
Figure BDA0002695520420000101
Diagenesis type, diagenesis strength, diagenesis evolution sequence, diagenesis fluid, hydrocarbon source rock hydrocarbon discharge time and strength, oil and gas reservoir formation and reservoir formation process in geological period are closely and internally connected with diagenesis stage of reservoir; deep and ultra-deep reservoirs with strong heterogeneity usually have higher diagenesis in the process of burying and undergo complex diagenesis conversion; by finely dividing the diagenesis stage and further finely depicting the diagenesis evolution process, the diagenesis environment, diagenesis events, diagenesis actions, diagenesis reservoir formation process and the difference of the time-space coupling relation thereof can be accurately revealed;
the reservoir diagenetic stage is divided on the basis of diagenetic environment parameters such as temperature (T), time (T), pressure (p) and the like, and geological variables such as particle contact relation, pore structure, authigenic mineral combination, vitrinite reflectivity (Ro%), organic matter peak-solving temperature (Tmax), illite-smectite mixed layer clay mineral conversion rate (I/S-S%), sterane isomerization index (SI), quartz secondary growth increasing index and the like are integrated, so that the reservoir diagenetic stage is divided from the angle of combination of qualitative and quantitative;
in the exploration and development process of oil and gas fields, the burying history and the thermal history of a sedimentary basin are required to be researched, the lithofacies paleogeography and sedimentary environment of the basin in different geological periods are modeled, the formation and development of a structure are analyzed, the maturity threshold depth of organic matters of hydrocarbon source rocks is determined, the oil and gas maturity horizon is judged, the oil and gas migration history is predicted and the like;
preferably, the clastic rock reservoir diagenetic phase characterization method further comprises the steps of simulating time domain and space domain distribution of clastic rock reservoir diagenesis after dividing the clastic rock reservoir diagenesis phase;
after the diagenetic stage is divided, diagenetic time domain and spatial domain distribution of the clastic reservoir can be simulated.
Example 2
The embodiment of the invention provides a fine characterization system for a diagenetic stage of a clastic rock reservoir, which comprises a stratum paleo-temperature acquisition module, a reflectivity and montmorillonite ratio acquisition module and a diagenetic stage division module;
the stratum paleo-temperature acquisition module is used for acquiring a ground temperature gradient, a stratum buried depth, a constant temperature zone buried depth and an paleo-surface temperature, and acquiring a stratum paleo-temperature according to the ground temperature gradient, the stratum buried depth, the constant temperature zone buried depth and the paleo-surface temperature;
the reflectivity and montmorillonite ratio acquisition module is used for acquiring the reflectivity of a source rock vitrinite and the montmorillonite ratio of an illite/montmorillonite mixed layer according to the formation paleo-temperature;
and the diagenetic stage division module is used for dividing the diagenetic stages of the clastic rock reservoir according to the paleo-temperature of the stratum, the reflectivity of the source rock vitrinite and the montmorillonite ratio of the illite/montmorillonite layer.
Example 3
The invention provides a clastic rock reservoir diagenetic phase fine characterization device which comprises a processor and a memory, wherein the memory is stored with a computer program, and the computer program is executed by the processor to realize the clastic rock reservoir diagenetic phase characterization method in the embodiment 1.
Example 4
An embodiment of the present invention provides a computer-readable storage medium, on which a computer program is stored, which, when executed by a processor, implements the clastic rock reservoir diagenetic stage characterization method as described in embodiment 1.
The invention discloses a clastic rock reservoir diagenetic stage characterization method, a clastic rock diagenetic stage characterization system, a clastic rock diagenetic stage characterization device and a computer readable storage medium, wherein the method comprises the steps of obtaining a geothermal gradient, a stratum buried depth, a constant temperature zone buried depth and an ancient surface temperature, and obtaining the ancient temperature of the stratum according to the geothermal gradient, the stratum buried depth, the constant temperature zone buried depth and the ancient surface temperature; obtaining the reflectivity of a source rock vitrinite and the illite/montmorillonite ratio of the illite/montmorillonite mixed layer according to the formation paleo-temperature; dividing the diagenetic stages of the clastic rock reservoir according to the paleo-temperature of the stratum, the vitrinite reflectivity of the source rock and the montmorillonite ratio of the illite/montmorillonite mixed layer; the method realizes more precise characterization of the diagenetic rock reservoir diagenesis stage and is beneficial to highlighting the difference of the reservoir evolution process.
According to the technical scheme, three parameters of the paleoterrestrial temperature T (DEG C), the vitrinite reflectivity (Ro%) and the content (I/S-S%) of the smectite in the clay mineral illite/montmorillonite mixed layer are comprehensively considered, the diagenetic stage is finely divided, and the fine division scheme is used for fine characterization of reservoir diagenetic evolution and diagenetic numerical simulation, can highlight the difference of the reservoir evolutionary process and is beneficial to identifying key geological factors causing the reservoir differential evolution.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are also within the scope of the present invention.

Claims (10)

1. A clastic rock reservoir diagenesis stage characterization method is characterized by comprising the following steps:
acquiring a ground temperature gradient, a stratum buried depth, a constant temperature zone buried depth and an ancient surface temperature, and acquiring a stratum ancient temperature according to the ground temperature gradient, the stratum buried depth, the constant temperature zone buried depth and the ancient surface temperature;
obtaining the reflectivity of a source rock vitrinite and the illite/montmorillonite ratio of the illite/montmorillonite mixed layer according to the formation paleo-temperature;
and dividing the diagenetic stages of the clastic rock reservoir according to the paleo-temperature of the stratum, the vitrinite reflectivity of the source rock and the montmorillonite ratio of the illite/montmorillonite mixed layer.
2. The clastic rock reservoir diagenetic stage characterization method of claim 1, wherein obtaining the formation paleo-temperature according to the geothermal gradient, the formation burial depth, the burial depth of the constant temperature zone and the paleo-surface temperature, specifically comprises obtaining the paleo-temperature of the formation according to a formula
T=A*(H-H0)+T0
Obtaining the paleotemperature T of the stratum, wherein A is the ground temperature gradient, H is the buried depth of the stratum, and H0For burying at depth of constant temperature zone, T0Is the paleo-surface temperature.
3. The clastic rock reservoir diagenetic stage characterization method of claim 1, wherein obtaining source rock vitrinite reflectivity from the formation paleo-temperatures comprises, in particular, obtaining source rock vitrinite reflectivity according to a formula
Figure FDA0002695520410000011
Obtaining the reflectivity Ro% of the source rock vitrinite, wherein T is the formation paleo-temperature,
Figure FDA0002695520410000012
are coefficients.
4. The clastic rock reservoir diagenetic stage characterization method of claim 1, wherein obtaining the illite/smectite ratio according to the formation paleo-temperature comprises, in particular, obtaining the illite/smectite ratio according to a formula
I/S-S%=-3.293e-07*T^4+0.0001625*T^3-0.02466*T^2+0.5809*T+97.06
And obtaining the I/S-S% ratio of montmorillonite of the illite/montmorillonite mixed layer, wherein T is the paleotemperature of the stratum.
5. The clastic rock reservoir diagenetic stage characterization method of claim 1, wherein the clastic rock reservoir diagenetic stages are divided according to the formation paleo-temperature, the source rock vitrinite reflectance and the illite/montmorillonite ratio, and the method specifically comprises the steps of dividing the formation paleo-temperature in a range of less than or equal to 200 ℃, the source rock vitrinite reflectance in a range of greater than or equal to 4 and the I/S-S% illite/montmorillonite ratio in a range of 0 to 1 into 36 diagenetic stages; and determining the division of the clastic rock reservoir diagenetic stages according to the formation paleo-temperature, the source rock vitrinite reflectivity and the illite/montmorillonite ratio and 36 diagenetic stages.
6. The clastic rock reservoir diagenetic stage characterization method of claim 5, wherein the formation paleo-temperature in the range of less than or equal to 200 ℃, the vitrinite reflectance of the source rock in the range of greater than or equal to 4, and the I/S-S% illite/montmorillonite ratio in the range of 0 to 1 are divided into 36 diagenetic stages, including in particular,
dividing a diagenesis stage in which the paleo-temperature of the stratum is within a range of less than or equal to 25 ℃, the vitrinite reflectance of the source rock is within a range of less than 0.24, and the I/S-S% illite/montmorillonite ratio is greater than 99 and within a range of less than or equal to 100 into a1 st diagenesis stage;
dividing a diagenesis stage into 2 to 35 diagenesis stages, wherein the formation paleo-temperature is more than 25 ℃ and less than or equal to 195 ℃, the source rock vitrinite reflectivity is more than or equal to 0.24 and less than or equal to 3.45, and the I/S-S% illite/montmorillonite ratio is less than or equal to 99 and more than 1;
dividing a diagenesis stage into a 36 th diagenesis stage, wherein the paleotemperature of the stratum is more than 195 and less than or equal to 200 ℃, the vitrinite reflectivity of the source rock is more than or equal to 3.45 and less than or equal to 4, and the I/S-S% illite/montmorillonite ratio is less than or equal to 1 and more than or equal to 0.
7. The clastic reservoir diagenetic phase characterization method of claim 1, further comprising, after dividing the clastic reservoir diagenetic phase, simulating a time domain and a space domain distribution of the clastic reservoir diagenetic.
8. The clastic rock reservoir diagenetic stage fine characterization system of claim 1, comprising a formation paleo-temperature acquisition module, a reflectivity and montmorillonite ratio acquisition module and a diagenetic stage division module;
the stratum paleo-temperature acquisition module is used for acquiring a ground temperature gradient, a stratum buried depth, a constant temperature zone buried depth and paleo-surface temperature, and acquiring the stratum paleo-temperature according to the ground temperature gradient, the stratum buried depth, the constant temperature zone buried depth and the paleo-surface temperature;
the reflectivity and montmorillonite ratio acquisition module is used for acquiring the reflectivity of a source rock vitrinite and the montmorillonite ratio of an illite/montmorillonite mixed layer according to the formation paleo-temperature;
and the diagenetic stage division module is used for dividing the diagenetic stages of the clastic rock reservoir according to the paleo-temperature of the stratum, the vitrinite reflectivity of the source rock and the montmorillonite ratio of the illite/montmorillonite mixed layer.
9. A clastic reservoir diagenetic phase fine characterization apparatus comprising a processor and a memory, the memory having stored thereon a computer program which, when executed by the processor, implements a clastic reservoir diagenetic phase characterization method as claimed in any one of claims 1 to 7.
10. A computer-readable storage medium, having stored thereon a computer program, wherein the computer program, when executed by a processor, implements a clastic reservoir diagenetic phase characterization method as claimed in any one of claims 1 to 7.
CN202011004770.5A 2020-09-22 2020-09-22 Clastic rock reservoir diagenetic stage characterization method, system, device and storage medium Active CN112347608B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011004770.5A CN112347608B (en) 2020-09-22 2020-09-22 Clastic rock reservoir diagenetic stage characterization method, system, device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011004770.5A CN112347608B (en) 2020-09-22 2020-09-22 Clastic rock reservoir diagenetic stage characterization method, system, device and storage medium

Publications (2)

Publication Number Publication Date
CN112347608A true CN112347608A (en) 2021-02-09
CN112347608B CN112347608B (en) 2023-03-17

Family

ID=74358055

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011004770.5A Active CN112347608B (en) 2020-09-22 2020-09-22 Clastic rock reservoir diagenetic stage characterization method, system, device and storage medium

Country Status (1)

Country Link
CN (1) CN112347608B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112285322A (en) * 2020-10-12 2021-01-29 湖南科技大学 Fine characterization method for diagenesis stage in diagenesis process of clastic rock reservoir

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103628867A (en) * 2013-11-26 2014-03-12 中国石油天然气股份有限公司 Method and system for simulating and analyzing shale reservoir diagenetic evolution process
CN108344854A (en) * 2018-02-06 2018-07-31 长江大学 The cementing quantitative forecasting technique for subtracting hole amount of clastic reservoir rock based on diagenetic process
CN108345046A (en) * 2018-02-06 2018-07-31 长江大学 Diagenetic Facies division methods during clastic reservoir rock is buried
CN110517794A (en) * 2019-08-23 2019-11-29 长安大学 A method of that establishes shale gas reservoir buries-thermal evolution history figure

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103628867A (en) * 2013-11-26 2014-03-12 中国石油天然气股份有限公司 Method and system for simulating and analyzing shale reservoir diagenetic evolution process
CN108344854A (en) * 2018-02-06 2018-07-31 长江大学 The cementing quantitative forecasting technique for subtracting hole amount of clastic reservoir rock based on diagenetic process
CN108345046A (en) * 2018-02-06 2018-07-31 长江大学 Diagenetic Facies division methods during clastic reservoir rock is buried
CN110517794A (en) * 2019-08-23 2019-11-29 长安大学 A method of that establishes shale gas reservoir buries-thermal evolution history figure

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
QIAN WENDAO, YIN TAIJU,AND ETC,.: "Identification of Diagenetic Facies Based on Diagenetic Process: An Example from Paleogene Lake Sediment Sandstone in Bozhong X Oilfield", 《EARTH SCIENCES》, 7 July 2018 (2018-07-07) *
QIAN WENDAO,YIN TAIJU AND ETC,.: "Forming Condition and Geology Prediction Techniques of Deep Clastic Reservoirs", 《ACTA GEOLOGICA SINICA (ENGLISH EDITION)》, 30 April 2017 (2017-04-30) *
WENDAO QIAN, TAIJU YIN,AND ETC,.: "Diagenesis and diagenetic stages prediction of Ed2 reservoir in the west of Bozhong sag", 《PETROLEUM》, 19 April 2019 (2019-04-19) *
何东博,应凤祥等: "碎屑岩成岩作用数值模拟及其应用", 《石油勘探与开发》, 23 December 2004 (2004-12-23) *
国家经济贸易委员会: "《中国人民共和国石油天然气行业标准5477-2003》", 18 March 2003 *
孟元林,张磊等: "碎屑岩成岩阶段细分与应用—以松辽盆地齐家—古龙地区为例", 《中国矿业大学学报》, 4 June 2020 (2020-06-04) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112285322A (en) * 2020-10-12 2021-01-29 湖南科技大学 Fine characterization method for diagenesis stage in diagenesis process of clastic rock reservoir

Also Published As

Publication number Publication date
CN112347608B (en) 2023-03-17

Similar Documents

Publication Publication Date Title
Tobin et al. Reservoir quality modeling of tight-gas sands in Wamsutter field: Integration of diagenesis, petroleum systems, and production data
Chaudhary et al. Shale oil production performance from a stimulated reservoir volume
CA2852050C (en) Methods and systems of modeling hydrocarbon flow from kerogens in a hydrocarbon bearing formation
Kontorovich et al. Historical-geological modeling of hydrocarbon generation in the Mesozoic–Cenozoic sedimentary basin of the Kara Sea (basin modeling)
US8359184B2 (en) Method, program and computer system for scaling hydrocarbon reservoir model data
Wang et al. Effect of burial history on diagenetic and reservoir-forming process of the Oligocene sandstone in Xihu sag, East China Sea Basin
AU2018203671A1 (en) Method of exploiting a sedimentary basin comprising hydrocarbons, using stratigraphic modelling
Erfani et al. Impact of microheterogeneity on upscaling reactive transport in geothermal energy
Jia et al. Some key issues on the unconventional petroleum systems
Dang Mechanistic modeling, design, and optimization of low salinity waterflooding
CN112347608B (en) Clastic rock reservoir diagenetic stage characterization method, system, device and storage medium
CN109187161B (en) Quantitative evaluation method for feldspar corrosion degree in clastic rock
Wang et al. Assessment of CO2 storage potential in high water-cut fractured volcanic gas reservoirs—Case study of China’s SN gas field
Shan et al. Study on hydrocarbon accumulation periods based on fluid inclusions and diagenetic sequence of the subsalt carbonate reservoirs in the Amu Darya right bank block
CN115964842A (en) Method and system for determining natural gas hydrate resource amount
Huang Assessment and management of sedimentary geothermal resources
CN112285322B (en) Fine characterization method for diagenetic stage in diagenetic rock reservoir diagenetic process
CN114139432A (en) Fractured reservoir CO using neural network technology2Flooding simulation method
Moldabayeva et al. Hydrodynamic modeling of field development using enhanced oil recovery methods
Chen et al. Pore-Scale Investigation on Dynamic Permeability Characterization of Hydrate-Bearing Sediments
CN111101930A (en) Method for evaluating single-well exploitation yield-increasing potential in gas reservoir development mode
Liu et al. A new method for recovering paleoporosity of sandstone: case study of middle Es3 member of Paleogene formation in Niuzhuang Sag, Dongying Depression, Bohai Bay Basin in China
Khasanov et al. Modeling the acid treatment of a polymictic reservoir
Agzamov et al. On the degree of influence of waterflooding on the oil recovery factor from productive formations of high-viscosity reservoirs X, represented by terrigenous reservoirs
Saeed Effect of temperature-dependent relative permeability functions on thermal recovery of heavy oil

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