CN112347608B - 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
CN112347608B
CN112347608B CN202011004770.5A CN202011004770A CN112347608B CN 112347608 B CN112347608 B CN 112347608B CN 202011004770 A CN202011004770 A CN 202011004770A CN 112347608 B CN112347608 B CN 112347608B
Authority
CN
China
Prior art keywords
temperature
paleo
equal
illite
diagenetic
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.)
Active
Application number
CN202011004770.5A
Other languages
Chinese (zh)
Other versions
CN112347608A (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

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/10—Geothermal 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 gas exploration degree, deep ultra-deep layer oil gas exploration becomes the most main exploration target in the future, and deep ultra-deep layer reservoir undergoes a more complex diagenetic evolution process in the burying process, so that the research on the diagenetic stage is enhanced, and important guidance can be provided for the effectiveness evaluation of the reservoir in the oil gas reservoir 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, a domestic student Ying Fengxiang systematically searches and collects the results of the domestic and foreign students on diagenesis research and diagenesis stage division, provides the specification of clastic rock diagenesis stage division in the petroleum and gas industry, and divides diagenesis evolution of clastic rock into an early diagenesis stage and a late diagenesis stage; in the beginning of the twenty-first century, ying Fengxiang modifies the earlier clastic rock diagenesis stage division standard, clastic rock diagenesis is divided into 4 stages of early, middle and late, and a diagenesis stage is divided into A, B, wherein the A stage of the diagenesis stage is further divided into two subphases of A1 and A2; meng Yuanlin on the basis of former manual work, the diagenesis stage A of clastic rock is subdivided into 4 diagenesis micro-stages.
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 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 stratum paleo-temperature according to the geothermal gradient, stratum burial depth, constant temperature zone burial depth and paleo-surface temperature, and concretely comprising the steps of
T＝A*(H-H 0 )+T 0
Obtaining the paleotemperature T of the stratum, wherein A is the ground temperature gradient, H is the buried depth of the stratum, and H 0 For constant temperature zone buried depth, T 0 Is the paleo-surface temperature.
Further, obtaining the reflectivity of the source rock vitrinite according to the paleo-temperature of the stratum, specifically comprising obtaining the reflectivity of the source rock vitrinite according to a formula
Figure SMS_1
Obtaining the reflectivity Ro% of the source rock vitrinite, wherein T is the formation paleo-temperature,
Figure SMS_2
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 clastic rock reservoir as described in any of the above technical solutions.
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 a source rock vitrinite and the illite/montmorillonite ratio according to the paleo-temperature of the stratum;
and S3, 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.
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-formation rock 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; when the reservoir evolves to a medium-formation rock A2 stage, densification is started, mechanical compaction is started to be converted into chemical compaction, particles are mostly in line-concave contact and concave-convex contact, and the particles are corroded and developed under the action of organic acid; after entering the B stage of diagenetic formation, the compaction mainly takes chemical compaction, but the overall strength is not large, the diagenetic environment is changed to weak alkali, the erosion holes formed in the A2 stage 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 smectite content (I/S-S%) 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 each diagenetic parameter 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-H 0 )+T 0
The paleo-temperature T of the stratum is obtained,wherein A is the earth temperature gradient, H is the stratum buried depth, H 0 For constant temperature zone buried depth, T 0 Is 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-H 0 )+T 0 (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, H 0 For constant temperature with buried depth, T 0 Is 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 SMS_3
Obtaining the reflectivity Ro% of the source rock vitrinite, wherein T is the formation paleo-temperature,
Figure SMS_4
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 burying 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,therefore, when the vitrinite reflectivity is used for judging the maturity of the organic matters, different types of kerogen are distinguished, and the coefficient can be used
Figure SMS_5
Carrying out correction;
Figure SMS_6
preferably, the obtaining of 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 montmorillonite ratio I/S-S% of the illite/montmorillonite mixed layer, wherein T is the formation paleo-temperature.
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 as 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 middle diagenetic stage A1 (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 middle diagenetic stage A2 (110-140 ℃); the dissolution of carbonate minerals and aluminosilicate minerals is closely related to the hydrocarbon formation of organic matters, because organic acid is formed by the thermal evolution of the organic matters, so that the property of the aqueous solution is changed, and the organic acid is very stable and can be stored at the temperature of 80-120 ℃, so that the quantity and the proportion of elements of the aluminosilicate minerals and the carbonate cement entering the solution are positively influenced;
the dissolution of feldspar is also influenced by the components, when the content of plagioclase in sandstone is higher, the dissolution of feldspar is relatively strong, mainly because the content of calcium ions (Ca & lt 2+ & gt) in plagioclase is higher, and the plagioclase is easier to combine 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 Li Danhe 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 grows exponentially with temperature, and oxygen isotope data of quartz overgrowth indicates that cementation easily occurs 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, 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 is a montmorillonite conversion model, and the theoretical basis is that the montmorillonite conversion strength is increased along with the increase of a diagenetic stage and is sharply reduced after reaching a middle diagenetic 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 formation paleo-temperature, the source rock vitrinite reflectivity and the illite/montmorillonite layer Dan Bi, 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 reflectivity in a range of more than or equal to 4 and the I/S-S% illite/montmorillonite layer montmorillonite ratio in a range of 0 to 1 into 36 diagenetic stages; 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-DS36;
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 SMS_7
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 from the angle of qualitative and quantitative combination by integrating the geological variables such as particle contact relation, pore structure, authigenic mineral combination, vitrinite group reflectivity (Ro%), organic matter Jie Fengwen (Tmax), illite smectite mixed layer clay mineral conversion rate (I/S-S%), sterane hopane isomerization index (SI), quartz secondary growth increasing index and the like on the basis of diagenetic environment parameters such as temperature (T), time (T), pressure (p) and the like;
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 the time domain and space domain distribution of clastic rock reservoir diagenetic after the clastic rock reservoir diagenetic phase is divided;
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
Embodiments of the present invention provide a computer readable storage medium having stored thereon a computer program 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 paleo-geothermal 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 evolution process and is beneficial to identifying key geological factors causing the reservoir difference 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 (8)

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;
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;
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 illite/montmorillonite ratio, and concretely comprising 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 diagenetic stages of the clastic rock reservoir according to the formation paleo-temperature, the source rock vitrinite reflectivity and the illite/montmorillonite ratio and 36 diagenetic stages;
dividing the formation paleo-temperature in the range of less than or equal to 200 ℃, the source rock vitrinite reflectivity 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 diagenetic stage into a1 st diagenetic stage, wherein the paleo-temperature of the stratum is within the range of less than or equal to 25 ℃, the vitrinite reflectivity of the source rock is within the range of less than 0.24, and the I/S-S% illite/montmorillonite ratio of the montmorillonite mixed layer is greater than 99 and within the range of less than or equal to 100;
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.
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-H 0 )+T 0
Obtaining the ancient temperature T of the stratum, wherein A is the ground temperature gradient, and H is the stratum burialDeep, H 0 For burying at depth of constant temperature zone, T 0 Is 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
Ro%=∂*(3.684e-10*T^5-1.822e-07*T^4+3.325e-05*T^3-0.002658*T^2+
0.09441*T-0.8632)
Obtaining reflectivity of source rock vitriniteRo%Wherein T is the formation paleo-temperature, and ∂ is the coefficient.
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 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.
6. A clastic rock reservoir diagenetic stage characterization system is characterized by comprising a stratum paleo-temperature acquisition module, a reflectivity and montmorillonite ratio acquisition module and a diagenetic stage division module;
the formation paleo-temperature acquisition module is used for acquiring a geothermal gradient, a formation burial depth, a constant temperature zone burial depth and a paleo-surface temperature, and acquiring a formation paleo-temperature according to the geothermal gradient, the formation burial depth, the constant temperature zone burial 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;
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 layer;
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 illite/montmorillonite ratio, and concretely comprising 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 value of the illite/montmorillonite ratio and 36 diagenetic stages;
dividing the formation paleo-temperature in the range of less than or equal to 200 ℃, the source rock vitrinite reflectivity 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 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. A clastic reservoir diagenesis phase characterization apparatus comprising a processor and a memory, the memory having stored thereon a computer program which, when executed by the processor, carries out a clastic reservoir diagenesis phase characterization method according to any one of claims 1 to 5.
8. 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 set forth in any one of claims 1 to 5.
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 CN112347608A (en) 2021-02-09
CN112347608B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112285322B (en) * 2020-10-12 2022-10-11 湖南科技大学 Fine characterization method for diagenetic stage in diagenetic rock reservoir diagenetic process

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110517794A (en) * 2019-08-23 2019-11-29 长安大学 A method of that establishes shale gas reservoir buries-thermal evolution history figure

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103628867B (en) * 2013-11-26 2016-08-31 中国石油天然气股份有限公司 The simulation of a kind of shale reservoir diagenesis evolution process and analysis method and system
CN108344854B (en) * 2018-02-06 2021-01-08 长江大学 Quantitative prediction method for clastic rock reservoir cementation pore-reducing amount based on diagenesis process
CN108345046B (en) * 2018-02-06 2020-07-03 长江大学 Lithomorphic phase dividing method in clastic rock reservoir burying process

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110517794A (en) * 2019-08-23 2019-11-29 长安大学 A method of that establishes shale gas reservoir buries-thermal evolution history figure

Also Published As

Publication number Publication date
CN112347608A (en) 2021-02-09

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
US8868392B2 (en) Method and system for predicting the effect of microbes injected into an oil-bearing reservoir
US8359184B2 (en) Method, program and computer system for scaling hydrocarbon reservoir model data
Dong et al. Hydrocarbon migration and accumulation in the Upper Cretaceous Qingshankou Formation, Changling Sag, southern Songliao Basin: Insights from integrated analyses of fluid inclusion, oil source correlation and basin modelling
Li et al. Mathematical modeling of secondary precipitation from sandstone acidizing
Dang Mechanistic modeling, design, and optimization of low salinity waterflooding
CN112347608B (en) Clastic rock reservoir diagenetic stage characterization method, system, device and storage medium
He et al. Evolution of prototype basin and change of tectonic-sedimentary pattern in Paleozoic, Ordos Basin
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
CN111101930B (en) Single-well exploitation production-increasing potential evaluation method in gas reservoir development mode
Carr et al. Some thoughts on the influence of pressure and thermal history assumptions on petroleum systems modelling
CN114139432A (en) Fractured reservoir CO using neural network technology2Flooding simulation method
Watson An investigation of the geothermal potential of the Upper Devonian sandstones beneath eastern Glasgow
Moldabayeva et al. Hydrodynamic modeling of field development using enhanced oil recovery methods
Gaoxuan et al. Assessment of the hofsstadir geothermal field, W-Iceland, by lumped parameter modelling, Monte Carlo simulation and tracer test analysis
Hantschel et al. Introduction to Basin modeling
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
Bartels et al. Reactive flow and permeability prediction—numerical simulation of complex hydrogeothermal problems
Asperger Rapid, high-temperature, field test method for evaluation of geothermal calcium carbonate scale inhibitors
Wang et al. The Study of Optimizing Reservoir Model Using Experimental Design in Stochastic Models
Gao et al. Characterization of Lower Paleozoic Water and Gas Distribution in Yan’an Gas Field, Ordos Basin
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