Method for inverting sedimentary basin ancient pore fluid property evolution process by comprehensively applying multiple geological data

Abstract

The invention relates to the technical field of oil exploration, in particular to a method for inverting a sedimentary basin ancient pore fluid property evolution process by comprehensively applying various geological data. According to the method, the type of water-rock reaction occurring in a reservoir and the property of ancient pore fluid are judged according to the depth range, formation time and comprehensive geological background of different diagenetic products in a research area, the property of formation water in the sandstone reservoir, the content and change of skeleton particles, the type and isotope characteristics of the diagenetic products and the O isotope characteristics of ancient pore water; and establishing a corresponding relationship among the diagenetic product, the formation water property, the skeleton particle change rule and the paleo-pore fluid property, and inverting the evolution process of the paleo-pore fluid in the research area according to the corresponding relationship. According to the method, the evolution process of the strata fluid in the petroliferous basin can be determined, and the determination of related chemical reaction types in the reservoir in different evolution stages has important theoretical significance for determining the formation mechanism of the deep reservoir and predicting the distribution mode of the deep reservoir.

Description

Technical Field

[0001] The present invention relates to the technical field of petroleum exploration, and in particular to a method for inverting the evolution process of paleo-pore fluid properties in a sedimentary basin by comprehensively applying a variety of geological data. Background Art

[0002] During the evolution of sedimentary basins, matter and energy migrate and redistribute using pore fluids as carriers. Pore fluids are the core and key to this process. The evolution of the physical and chemical properties of pore fluids leads to the dissolution and precipitation of minerals within the reservoir. This process affects the preservation and destruction of primary pores in sedimentary basins and controls the formation and evolution of secondary pores. In other words, the evolution of fluids within oil and gas basins is the essence and core of their diagenesis, and determines the formation, evolution, preservation, and distribution of deep favorable reservoirs within the basins. Therefore, determining the evolution process of pore fluids within oil and gas basins and determining the dissolution-precipitation process of various minerals under different types of pore fluid environments have important theoretical significance and practical application value for the inversion of reservoir evolution processes and distribution prediction within oil and gas basins.

[0003] Currently, researchers generally treat formation fluids and diagenesis as two independent scientific issues. Their origins and evolution are often determined based on their isotopic or trace element signatures.

[0004] Diagenesis is primarily studied based on sedimentological theories and petrological testing methods, focusing on diagenetic facies, diagenetic sequences, and diagenetic evolution processes to evaluate and predict secondary porosity development zones. Separating the different materials within a diagenetic system for separate studies makes it difficult to accurately invert the reservoir-fluid interaction process and the genesis of secondary porosity under different formation fluid environments.

[0005] Patent CN105572131A discloses a comprehensive paleofluid geochemical analysis method, primarily addressing existing methods' shortcomings: a single tracing method, subjective interpretation, loose integration with the geological background of reservoir formation, and low credibility of the obtained results. This invention effectively addresses the shortcomings of existing isotope geochemical tracing techniques in the study of formation fluid activity and oil and gas accumulation, including a single tracing method, low credibility of evidence, subjective interpretation, loose integration with the geological background, and unreliable results. This invention is of great significance for promoting the application of fluid geochemical tracing in oil and gas accumulation research.

[0006] Patent CN116401851A discloses a method for recovering the paleo-porosity and compaction evolution history of reservoirs. This method establishes a vertical profile of measured porosity and permeability using statistical methods; recovers paleo-porosity and determines the reservoir porosity evolution history using thin-section backstripping; forward models the compaction evolution history of the reservoir based on physical simulation experiments; and recovers reservoir physical property evolution curves based on thin-section backstripping inversion and physical simulation forward modeling, and establishes a reservoir physical property parameter recovery chart. By combining thin-section backstripping inversion with physical simulation forward modeling, the dynamic evolution process of the reservoir is analyzed, resolving the problem that conventional physical property analysis and testing only yield static results. By recovering reservoir physical property evolution curves and establishing a reservoir physical property parameter recovery chart, the problem that conventional analysis of the diagenetic evolution of high-quality reservoirs can only be qualitative is resolved. The method also quantitatively recovers reservoir physical property parameters at critical moments of oil and gas accumulation and clarifies the compaction evolution history of the reservoir. This patent only focuses on the physical property evolution process of sandstone at the thin section scale, and does not focus on the evolution process of the core element of the process, the paleo-pore fluid.

[0007] Patent CN108388709B discloses a quantitative reservoir porosity simulation method based on diagenetic facies prediction. This method, based on research into the sedimentary facies, lithology, fluid properties, sedimentary cycles, diagenetic stages, and fault development of clastic reservoirs, establishes a "geological parameter-diagenesis-porosity" prediction model under different geological environments and diagenetic processes. This model then predicts the porosity of clastic reservoirs and determines the spatial distribution of pores within them, providing a basis for reservoir evaluation. This patent primarily identifies sandstone physical properties based on well logging data and does not address the evolution of paleo-pore fluids.

[0008] The methods used in the above patents are different from those of the present invention, and their accuracy in the evolution of the properties of paleo-pore fluids in sedimentary basins is not as good as that of the present invention.

[0009] In short, the technical solutions, technical problems to be solved and beneficial effects produced by the above-mentioned disclosed technologies are all different from the present invention. Regarding more technical features, technical problems to be solved and beneficial effects of the present invention, the above-mentioned disclosed technical documents do not provide any technical inspiration. Summary of the Invention

[0010] The main purpose of the present invention is to provide a method for inverting the evolution process of paleopore fluid properties in sedimentary basins by comprehensively applying a variety of geological data. The present invention is based on the fact that the changes in formation fluids, sandstone skeleton particles, and diagenetic products are the result of the coordinated evolution within a water-rock reaction system. The current formation fluid properties, skeleton particle evolution, and diagenetic product properties can all reflect the paleopore fluid properties and evolution process in the region. The method of the present invention utilizes formation water constant ion analysis to classify formation water and determine the depletion-enrichment characteristics of constant ions of different types of formation water; utilizes cast thin section observation to determine the evolution process of the skeleton particle components of the sandstone reservoir during burial; utilizes cast thin section observation, scanning electron microscopy, elemental analysis, isotope analysis, and inclusion analysis as means to determine the genetic mechanism and material source of the diagenetic products; and finally determines the paleopore fluid properties and their evolution process in the study area by integrating the research results of these three methods and combining them with the regional geological overview. The present invention forms a method for inverting the evolution process of paleopore fluid properties in sedimentary basins by comprehensively applying a variety of geological data.

[0011] To achieve the above object, the present invention adopts the following technical solutions:

[0012] In one aspect, the present invention provides a method for inverting the evolution of paleopore fluid properties in a sedimentary basin by comprehensively applying multiple geological data, comprising the following steps:

[0013] S1. Determine the formation water properties within the sandstone reservoirs in the study area;

[0014] S2. Sandstone samples from different burial depths in the study area were selected to grind thin sections of casts and determine the variation of the content of various types of skeleton particles with depth;

[0015] S3. Determine the types and isotopic characteristics of diagenetic products, determine the oxygen isotopic characteristics of pore water in the reservoir during the formation of different carbonate cements, and determine the formation temperatures of carbonate cements of different occurrences and authigenic quartz of different occurrences;

[0016] S4. Determine the depth range and formation time of different occurrences and types of diagenetic products;

[0017] S5. Determine the type of water-rock reaction and paleopore fluid properties within the reservoir based on the depth range and formation time of different diagenetic products in the study area, the geological background, the properties of formation water within the sandstone reservoir, the content and variation of skeletal grains, the types and isotopic characteristics of diagenetic products, and the O isotopic characteristics of paleopore water;

[0018] Establish the corresponding relationship between diagenetic products, formation water properties, skeleton particle change patterns and paleo-pore fluid properties; and invert the evolution process of paleo-pore fluid in the study area based on the corresponding relationship.

[0019] Furthermore, in S1:

[0020] S101. Determine the hydrochemical characteristics of the formation water in the study area, the content of the major ions it contains, and its variation trend with burial depth;

[0021] S102. Classify the formation water in the reservoir and clarify the spatial distribution characteristics of different types of formation water;

[0022] S103. Utilize the evaporation curves of different water bodies to determine the depletion-enrichment characteristics of constant particles in different types of formation water, and thus determine the properties of formation water in sandstone reservoirs.

[0023] Furthermore, in S101, the hydrochemical characteristics of the formation water are the total mineralization of the formation water and its variation trend with burial depth;

[0024] In S102, the SPSS cluster analysis module is used to classify the formation water in the reservoir based on the constant ion content of the formation water. The specific method is as follows:

[0025] Taking the milligram equivalents of various constant ions in different formation water samples as the numerical values, the formation water samples were analyzed using the cluster analysis module in SPSS20 software, and a dendrogram was made and classified.

[0026] Furthermore, in S103, the specific method for determining the depletion-enrichment characteristics of constant particles in different types of formation water using the evaporation curves of different water bodies is as follows:

[0027] Using chloride ions as a standard, a relationship diagram between chloride ion concentration and the constant ions was drawn, and evaporation-concentration curves of different water bodies were drawn at the same time to determine the loss-enrichment characteristics of constant particles in different types of formation water.

[0028] Furthermore, the constant ions include Na + , K + , Ca 2+ Mg 2+ 、SO4 2- 、HCO3 - .

[0029] Furthermore, in S2, sandstone samples from different burial depths in the study area were selected to grind thin sections of castings and observed under a polarizing microscope to determine the variation pattern of the content of various types of skeleton particles with depth and clarify their spatial variation pattern;

[0030] The rock skeleton particles in the sandstone reservoir mainly include: the proportion of quartz particles in the total particles, the proportion of feldspar particles in the total particles, the proportion of rock debris particles in the total particles, the proportion of sedimentary rock debris particles in the total rock debris particles, the proportion of igneous rock debris particles in the total rock debris particles, the proportion of igneous rock debris particles in the total rock debris particles, the proportion of potassium feldspar particles in the total feldspar, and the proportion of plagioclase particles in the total feldspar.

[0031] Furthermore, in S3:

[0032] S301. Based on polarizing microscope observation of cast thin sections, determine the type, content, occurrence, contact and replacement relationship, and formation sequence of diagenetic products within the sandstone reservoirs in the study area; and calculate the vertical distribution and spatial distribution characteristics of diagenetic products of different occurrences;

[0033] S302. Select samples of carbonate cements of different occurrences and authigenic quartz cements of different occurrences, prepare inclusion thin sections, and determine the homogenization temperature of the inclusions, i.e., the formation temperature of the corresponding minerals;

[0034] S303. Based on the observation results of S301 and S302, select typical samples for isotope thin section observation and determine the δ 13 C PDB and δ 18 O PDB ;

[0035] S304. Based on S302 and S303, determine the O isotope characteristics of the pore water in the reservoir during the formation period of different carbonate cements, and determine the formation temperature of carbonate cements of different occurrences and the formation temperature of authigenic quartz of different occurrences.

[0036] Furthermore, in S303, the typical samples are all samples of carbonate cements observed and temperature measured in S302, as well as samples in S301 in which the content of carbonate cements of a certain occurrence accounts for more than 70% of the total content of carbonate cements in the sample, and the single crystals or aggregates of the carbonate cements of the occurrence are larger than 50um.

[0037] Furthermore, in S304, the method for determining the formation temperature of carbonate cements of different occurrences and the formation temperature of authigenic quartz of different occurrences is as follows:

[0038] The δ 13 C PDB and δ 18 O PDB At the same time, the homogenization temperature of inclusions of carbonate cements of different occurrences and authigenic quartz of different occurrences was tested using the hot and cold stage programmed temperature method;

[0039] Based on the δ18 O PDB and δ 18 O PDB , comprehensively considering the homogenization temperature of brine inclusions in carbonate cements of different occurrences, determine the formation temperature of carbonate cements of different occurrences; based on the homogenization temperature of brine inclusions in authigenic quartz of different occurrences, determine the homogenization temperature of authigenic quartz of different occurrences.

[0040] Furthermore, in S4, a burial history-thermal history map of the study area is established. Based on the homogenization temperatures of carbonate cements of different occurrences and authigenic quartz of different occurrences determined in S304 and the contact and replacement relationship of different diagenetic products determined in S301, the diagenetic sequence of the study area is established; the formation depth range and formation time of diagenetic products of different occurrences and types are determined.

[0041] In another aspect, the present invention provides a system for inverting the evolution process of the ancient pore fluid properties in a sedimentary basin by comprehensively applying multiple geological data, including a processor and a memory storing a computer program, and the processor is configured to execute the computer program to implement a method for inverting the evolution process of the ancient pore fluid properties in a sedimentary basin by comprehensively applying multiple geological data.

[0042] In a third aspect, the present invention provides a computer storage medium having a computer program stored thereon, which, when executed by a processor, implements a method for inverting the evolution process of paleo-pore fluid properties in a sedimentary basin by comprehensively applying a variety of geological data.

[0043] Compared with the prior art, the present invention has the following advantages:

[0044] The method proposed in this paper utilizes constant ion analysis of formation water to determine the source and evolution of formation fluids. It also uses thin-section observation of casts, isotope testing, and inclusion observation to determine the evolution of the sandstone reservoir's skeletal particle composition and diagenetic products. Finally, a method is proposed that comprehensively utilizes the paleo-formation fluid properties of formation fluids and diagenetic products to determine the water-rock reaction processes experienced by the reservoir under different formation fluid environments and the evolution of reservoir physical properties. This method can clarify the evolution of formation fluids within oil and gas basins and identify the relevant chemical reaction types within reservoirs at different evolutionary stages, which has important theoretical significance for determining the genetic mechanism of deep reservoirs and predicting their distribution patterns. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] Figure 1 The main constant ion composition of the current formation water in Bonan Sag;

[0046] Figure 2 Cluster analysis of current formation water chemistry in the Bonan Subsag;

[0047] Figure 3Conventional ion ratios and salinity characteristics of different types of formation water in the Bonan Sag;

[0048] Figure 4 Depletion-enrichment characteristics of different types of formation water in the Bonan Sag;

[0049] Figure 5 Variation of skeleton grain content in sandstone samples from different strata in the Bonan Sag;

[0050] Figure 6 Typical diagenetic phenomena and contact-replacement relationships of sandstone reservoirs in the Bonan Sag;

[0051] Figure 7 Vertical distribution characteristics of carbonate cements of different occurrences in the deep sag area of ​​Bonan Sag;

[0052] Figure 8 Vertical distribution characteristics of authigenic kaolinite of different occurrences in the deep sag area of ​​Bonan Sag;

[0053] Figure 9 Homogenization temperature of brine inclusions in carbonate cements of different occurrences in sandstone reservoirs of Bonan Sag;

[0054] Figure 10 Homogenization temperature of brine inclusions in authigenic quartz of different occurrences in sandstone reservoirs in the Bonan Sag;

[0055] Figure 11 Isotopic characteristics of carbonate cements of different occurrences in sandstone reservoirs of Bonan Sag;

[0056] Figure 12 Oxygen isotope evolution history of paleopore fluids;

[0057] Figure 13 Diagenetic sequence of sandstone reservoirs in the Bonan Sag;

[0058] Figure 14 Evolution of paleopore fluids in different tectonic zones of the Bonan Sag. DETAILED DESCRIPTION

[0059] It should be noted that the following detailed descriptions are exemplary and intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the art to which the present invention belongs.

[0060] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is intended to include the plural form. In addition, it should be understood that when the terms "comprise" and "include" are used in this specification, they indicate the presence of features, steps, operations and combinations thereof.

[0061] In order to enable those skilled in the art to more clearly understand the technical solution of the present invention, the technical solution of the present invention will be described in detail below with reference to specific embodiments.

[0062] Example 1

[0063] A method for inverting the evolution process of paleopore fluid properties in a sedimentary basin by comprehensively applying multiple geological data, the method comprising the following steps:

[0064] S1. Statistical analysis of the hydrochemical characteristics of formation water in the study area, the content of constant ions contained therein, and their changing trends with burial depth; the hydrochemical characteristics of formation water are the total mineralization of formation water, the content of constant ions, and their changing trends with burial depth. The constant ions include but are not limited to Na + , K + , Ca 2+ Mg 2+ 、SO4 2- 、HCO3 - .

[0065] Classify the formation water in the reservoir. Using the cluster analysis module in SPSS, the formation water in the reservoir is classified based on the constant ion content and mineralization of the formation water: determine the mineralization (TDS (g / L)), sodium chloride coefficient (rNa / rCl), metamorphic coefficient (r(Cl-Na) / rMg), desulfurization coefficient (rSO4 2- *100 / rCl), carbonate equilibrium coefficient (r(HCO3 - +CO3 2- ) / rCa), if one or more hydrochemical parameters of different types of formation water have obvious differences, then the classification method is verified to be scientific and effective, and can serve as the basis for the next step of research; at the same time, - Ions are used as calibration to make Cl - The relationship diagram between ion concentration and the constant ions is drawn to draw the evaporation-concentration curves of different types of formation water to determine the depletion-enrichment characteristics of constant particles in different types of formation water.

[0066] S2. Based on cast thin section observation, the types and contents of skeleton particles in sandstone reservoir samples were statistically analyzed; the observation results were statistically analyzed to determine the composition of rock skeleton particles and their longitudinal variation patterns in the sandstone reservoirs in the study area; the skeleton particles mainly include: the proportion of quartz particles in the total particles (quartz particles / total particles (%)); the proportion of feldspar particles in the total particles (feldspar particles / total particles (%)); the proportion of rock fragments in the total particles (rock fragments / total particles (%)); the proportion of sedimentary rock fragments in the total rock fragments (sedimentary rock fragments / total rock fragments (%)); the proportion of igneous rock fragments in the total rock fragments (igneous rock fragments / total rock fragments (%)); the proportion of igneous rock fragments in the total rock fragments (igneous rock fragments / total rock fragments (%)); the proportion of potassium feldspar particles in the total feldspar (potassium feldspar particles / feldspar particles (%)); the proportion of plagioclase particles in the total feldspar (plagioclase particles / feldspar particles (%));

[0067] S301. Based on polarizing microscope observation of cast thin sections, determine the type, content, occurrence, and contact and replacement relationship of the diagenetic products in the study area; and calculate the vertical distribution and spatial characteristics of diagenetic products with different occurrences.

[0068] S302. Select the samples of carbonate cements and authigenic quartz with good crystal form and large crystals observed in S301 to grind into inclusion thin sections (to ensure that the carbonates and authigenic quartz cements of different occurrences observed in S301 are basically covered in the samples and have a certain number of observation samples); observe under a microscope, select brine inclusions of appropriate size and shape (round or elliptical, >5um) in carbonate cements of different occurrences and authigenic quartz, and use the hot and cold stage programmed temperature technique to test the uniform temperature of brine inclusions in different diagenetic products.

[0069] During the programmed temperature increase, when the temperature is below 70°C, the temperature is increased in 10°C increments; when the temperature is above 70°C, the temperature is increased in 5°C increments. After each temperature increase, the sample is stopped for 20 minutes to observe whether the inclusions have reached a homogeneous state. The homogenization temperature of the salt water inclusions in the mineral is the formation temperature of the mineral.

[0070] S303. Grind the remaining samples of the samples selected in S302 to prepare isotope test slices; using the results observed in S301, select samples with a certain occurrence of carbonate cement accounting for more than 70% of the total carbonate cement content and with large single crystals or aggregates of carbonate cement (>50um), and grind isotope test slices. Observe the isotope test slices under a microscope, and select carbonate cements of different occurrences. Use laser ablation technology to vaporize them, and the generated gas enters the mass spectrometer to measure the δ 13 C PDB (‰) and δ 18 OPDB (‰).

[0071] S304: Homogenization temperature and δ of brine inclusions of carbonate cements of the same sample and occurrence in S302 and S303 18 O PDB (‰) was tested. The homogenization temperature and δ 18 O PDB (‰); the homogenization temperature of the brine inclusions was taken as the formation temperature of the carbonate cements in this period, and the δ 18 O PDB (‰) is used as the O isotope characteristic value of the mineral. Using the carbonate cement O isotope characteristic differentiation chart proposed by Friedman I et al., 1977, the O isotope characteristics of the paleopore fluids during the formation period of carbonate cements of different occurrences are calculated.

[0072] Using the calculated O isotope characteristic values ​​of paleopore fluids during the formation period of carbonate cements of different occurrences and the carbonate cement O isotope characteristic differentiation chart proposed by Friedman I et al., 1977, the formation temperature of carbonate cement samples without brine inclusion testing was calculated. Note that the calculation process must ensure that the occurrence of the calculated carbonate cement and the O isotope characteristics of the paleopore fluid selected for use must correspond one to one.

[0073] The formation temperatures of carbonate cements of different occurrences calculated based on O isotope characteristics and the homogenization temperatures of brine inclusions within carbonate cements of different occurrences are used to determine the formation temperatures of carbonate cements of different occurrences (direct superposition is sufficient).

[0074] The uniform temperature of brine inclusions in authigenic quartz of different occurrences is taken as the formation temperature of authigenic quartz of that occurrence.

[0075] S4: Collect statistics on single-well stratification data, the time of start and end of stratum deposition, the start and end time of erosion, erosion thickness, thermal history data, etc. of key wells in the study area, and use the single-well burial history map production module in the petrel software to establish a burial history-thermal history (geothermal and evolution) map of the study area.

[0076] Based on the homogenization temperature of carbonate cements of different occurrences and authigenic quartz of different occurrences determined by S304, they are stuck on the burial history-thermal history diagram of typical wells;

[0077] Based on the contact and metasomatism relationships of different diagenetic products determined by S301, other diagenetic products of different occurrences (diagenetic products other than carbonate cements and quartz secondary enlargement, such as clay minerals and authigenic pyrite) are placed on the burial history-thermal history map to establish the diagenetic sequence of the study area; the formation depth range and formation time of different occurrences and types of diagenetic products are determined;

[0078] It can be read directly from the diagenetic sequence map: the horizontal axis of the map (bottom) is the time of mineral formation, in Ma (million years); the vertical axis of the map (left) is the temperature and depth of mineral formation.

[0079] S5: Based on the depth range and formation time of different diagenetic products extracted in S4, the geological background of this stage, the properties of formation water in the sandstone reservoir buried in this stage (S1), the content and changes of skeleton particles in the sandstone reservoir buried in this stage (S2), the types and isotopic characteristics of diagenetic products in the sandstone reservoir buried in this stage (S301, S303), and the O isotopic characteristics of the ancient pore water in this stage (S304) are comprehensively considered to comprehensively judge the type of water-rock reaction occurring in the reservoir at this stage and the properties of the ancient pore fluid, and establish the corresponding relationship between diagenetic products, current formation water properties, the change law of skeleton particles and the properties of the ancient pore fluid; based on this corresponding relationship, the evolution process of the ancient pore fluid in the study area is inverted.

[0080] Example 2

[0081] Taking the Bonan Subsag as an example, the method for inverting the evolution process of paleo-pore fluid properties in a sedimentary basin by comprehensively applying multiple geological data includes the following steps:

[0082] S1: Statistical analysis of the hydrochemical characteristics of formation water in the study area, the content of constant ions contained in it, and its changing trend with burial depth ( Figure 1 ); mainly including TDS (total mineralization) of formation water, Na + +K + , Ca 2+ Mg 2+ 、SO4 2- 、HCO3 - The concentration of constant ions.

[0083] The cluster analysis module of SPSS, a mathematical statistics software, was used to classify the formation water in the study area. This classification is mainly based on the similarity of the formation water chemistry (TDS and the concentration of constant ions). Figure 2 ). Figure 2The horizontal axis represents the chemical composition of hydrochemical samples, while the vertical axis represents cluster distance. As cluster distance increases, hydrochemical samples of varying similarity gradually cluster together, depending on the degree of similarity. Within a certain range of cluster distance variations, the classification results remain stable, which is the classification standard selected by the study (in this study, the classification distance was 12, dividing the formation water samples into three categories: Class I, Class II, and Class III).

[0084] Determine the mineralization degree (TDS (g / L)), sodium chloride coefficient (rNa / rCl), metamorphic coefficient (r(Cl-Na) / rMg), desulfurization coefficient (rSO4 2- *100 / rCl), carbonate equilibrium coefficient (r(HCO3 - +CO3 2- ) / rCa)( Figure 3 ), Figure 3 It shows that there are obvious differences in the salinity and various constant ion parameters of different types of formation water, indicating that the classification has certain scientificity and accuracy and can serve as the basis for the next step of research.

[0085] Using the evaporation curves of different water bodies (lake water, sea water, early formation water), the loss-enrichment characteristics of constant particles in different types of formation water are determined. + +K + ) with Cl - The concentration change curve is basically located near the evaporation line of lake water. This type of formation water increases with the increase of Cl- concentration (Na + +K + ) increases slightly. + +K + ) with Cl - The concentration evolution trend is close to the evaporation line of Class I formation water. Similar to Class I formation water, with the increase of Cl- concentration, (Na + +K + The increase rate of Class III formation water has slightly decreased. Class III formation water as a whole has greatly deviated from the concentration line of lake water, sea water and Class II formation water, but is closer to the evaporation line of Class II formation water. In addition, Class III formation water obviously shows a stronger + +K + ) loss or Cl - Enrichment ( Figure 4 -a). Ca in Type I formation water 2+ With Cl - The concentration change curve is located near the evaporation line of the lake water, with a slight Ca 2+ Loss; the evaporation line Ca of type II formation water compared with seawater, lake water and type I formation water 2+Although the type III formation water shows a significant Ca 2+ enriched, but compared with type II formation water, it shows Ca 2+ Loss or Cl - The enrichment ( Figure 4 -b). Compared with the concentration lines of seawater and lake water, both Type I and Type II formation water have Mg 2+ The phenomenon of loss, Class III formation water is Mg 2+ Loss or Cl - ion enrichment; relative to the concentration line of type I formation water, type II formation water partially shows Mg 2+ Enrichment characteristics; relative to the concentration line of type II formation water, type III formation water partially shows Mg 2+ Enrichment features ( Figure 4 -c). Class I formation water SO4 2- The content of Cl increases with the degree of concentration (Cl - content increased) decreased significantly, indicating that SO4 2- Compared with the concentration lines of seawater and lake water, Class II formation water also shows obvious SO4 2- Loss phenomenon; Class III formation water current SO4 2- Compared with the concentration line of Class II formation water, the content characteristics show SO4 2- The phenomenon of relative enrichment ( Figure 4 -d). Type I formation water HCO3 - The content of Cl increases with the concentration - The content of HCO3 in the evolution of Type I formation water decreased significantly. - Compared with the evaporation line of seawater and lake water, Class II formation water has obvious HCO3 - The enrichment characteristics of type III formation water are HCO3 - The loss phenomenon is obvious ( Figure 4 -e). HCO3 - Concentration and CO3 2- The concentration is obviously negatively correlated, indicating that HCO3 - To CO3 2- The transformation ( Figure 4 -f).

[0086] S2: Based on the observation of cast thin sections, the types and contents of skeleton particles in sandstone reservoir samples were counted; the observation results were statistically analyzed to determine the composition of rock skeleton particles and their longitudinal variation patterns in the sandstone reservoirs in the study area ( Figure 5); the skeleton particles mainly include: the proportion of quartz particles in the total particles (quartz particles / total particles (%)); the proportion of feldspar particles in the total particles (feldspar particles / total particles (%)); the proportion of rock debris particles in the total particles (rock debris particles / total particles (%)); the proportion of sedimentary rock debris particles in the total rock debris particles (sedimentary rock debris particles / total rock debris particles (%)); the proportion of igneous rock debris particles in the total rock debris particles (igneous rock debris particles / total rock debris particles (%)); the proportion of igneous rock debris particles in the total rock debris particles (igneous rock debris particles / total rock debris particles (%)); the proportion of potassium feldspar particles in the total feldspar (potassium feldspar particles / feldspar particles (%)); the proportion of plagioclase particles in the total feldspar (plagioclase particles / feldspar particles (%)).

[0087] The composition of the skeleton particles in the sandstone reservoir changes with increasing burial depth. When the burial depth is less than 3000m, the proportion of quartz particles in the sandstone in the study area increases with increasing burial depth; the proportion of feldspar particles in the total particles decreases with increasing burial depth; and the proportion of rock debris particles in the total particles remains basically unchanged with burial depth. Figure 5 When the burial depth is greater than 3000m, the proportion of quartz particles and rock debris particles in the total particles decreases with increasing burial depth; the proportion of feldspar particles in the total particles decreases with increasing burial depth. During the burial diagenesis process, the changes in the content of different types of feldspar particles also vary ( Figure 5 ). When the burial depth is less than 3000m, the ratio of potassium feldspar particles to feldspar particles decreases with increasing burial depth; the ratio of plagioclase particles to feldspar particles increases with increasing burial depth. When the burial depth is greater than 3000m, the ratio of plagioclase particles to feldspar particles decreases with increasing burial depth; the ratio of potassium feldspar particles to feldspar particles increases with increasing burial depth. During the burial diagenesis process, there are also differences in the changes in the content of different types of rock debris particles. When the burial depth is less than 3000m, the ratio of various types of rock debris particles to the total rock debris particles is relatively stable and does not change significantly with burial depth. When the burial depth is greater than 3000m, the ratio of sedimentary rock debris to the total rock debris increases with increasing burial depth; the ratio of igneous rock debris and metamorphic rock debris to the total rock debris increases with increasing burial depth ( Figure 5 ).

[0088] S301: Based on polarizing microscope observation of cast thin sections, determine the type, content, occurrence, and contact and replacement relationship of the diagenetic products in the study area; and calculate the vertical distribution of diagenetic products with different occurrences.

[0089] Based on the differences in the optical properties of different minerals, the types of diagenetic products developed within the rock sample are determined (for example, when observing thin sections of casts using an optical microscope, quartz is a first-order grayish white without twinning; feldspar is a first-order grayish white with slight twinning structures such as kappa twinning or polysynthetic twinning; dolomite is generally rhombic, has two sets of cleavage, and has a high-order white interference color). Based on the morphological characteristics of different types of diagenetic products, the stages of different diagenetic products, the contact relationships between diagenetic products of different stages, and the order of diagenetic products of different stages are determined (for example, calcite filling uncompacted pores and calcite filling feldspar dissolution pores are morphologically distinct stages; calcite filling the outer side of quartz secondary overgrowths formed later than the quartz secondary overgrowths). The probability method is used to calculate the vertical distribution range of different diagenetic products. For example, there are 50 thin section samples in the range of 2000m-2100m, of which 20 contain a certain type of diagenetic product with a certain occurrence. Then the content of this type of diagenetic product in the range of 2000m-2100m is 40%. Based on this method, the vertical distribution characteristics of different diagenetic products can be established.

[0090] The study area has discovered multiple phases and types of diagenetic products of different occurrences. These mainly include four different occurrences of carbonate cements, two different occurrences of authigenic kaolinite, four different occurrences of authigenic quartz, two different occurrences of chlorite and illite ( Figure 6 The first-stage carbonate cement (Cc1) of the sandstone reservoir in the study area is micritic calcite, which mainly appears in the form of equal-thickness ring edges at the edges of some sandstone particles ( Figure 6 -a). The content of carbonate cements in this period is low (the surface porosity is less than 2%); carbonate cements in this period are developed in sandstone reservoirs at a burial depth of 2000m-4200m, and the trend of their content changing with burial depth is not obvious ( Figure 7 The second-stage carbonate cement fills the primary pores and is located outside the first-stage carbonate cement ( Figure 6 -b). The content of carbonate cement (Cc2) in the sandstone reservoirs at a depth of 2000m-2300m increases with depth; the content of Cc2 in the sandstone reservoirs at a depth of 2300m-3000m does not change significantly with depth; the content of Cc2 in the sandstone reservoirs at a depth greater than 3000m decreases significantly with depth ( Figure 7 The third stage carbonate cement (Cc3) is filled in the feldspar dissolution pores and residual primary pores ( Figure 6 -c). The content of carbonate cement (Cc3) in this period (Cc2) increases with the depth in the sandstone reservoir at 2500m-2900m; the content in the sandstone reservoir at 2900m-4200m decreases with the depth. Figure 7 The fourth-stage carbonate cement (Cc4) is mainly developed outside the third-stage carbonate cement or in its dissolved pores ( Figure 6 -d); it is not developed in sandstones with a burial depth of less than 3400m, but its content in sandstones with a burial depth of more than 3400m increases with the increase of burial depth ( Figure 7 ).

[0091] Four types of authigenic quartz with different occurrences are developed in the study area. They are mainly quartz secondary enlargement edges; euhedral quartz grains; intragranular cracks in quartz grains; and quartz healing cracks that cut through the quartz secondary enlargement edges. The tertiary carbonate cement ( Figure 6 -e, f).

[0092] Two types of authigenic kaolinite with different occurrences are developed in the study area. The aggregates of the first-stage authigenic kaolinite are worm-like, filling primary pores or feldspar dissolution pores, and showing the characteristics of secondary enlargement of quartz ( Figure 6 -g). The content of kaolinite in this period increases with the increase of burial depth in sandstones with a burial depth of less than 2800m; the content in sandstones with a burial depth of 2800m-4000m remains basically stable; in sandstones with a burial depth of more than 4000m, the content of this type of kaolinite decreases with the increase of burial depth ( Figure 8 The second stage of authigenic kaolinite is characterized by the kaolinization of the entire feldspar particles, and the kaolinite aggregates are combined into the form of feldspar particles ( Figure 6 -h). This type of kaolinite coexists with euhedral quartz. This type of kaolinite is generally not developed in sandstones with a burial depth of less than 2800m; the content of this type of kaolinite in sandstones with a burial depth of 2800m-4000m increases with increasing burial depth; the content of this type of kaolinite in sandstones with a burial depth of more than 4000m decreases with increasing burial depth ( Figure 8 ).

[0093] Two types of pyrite with different occurrences are developed in the study area. The first-stage pyrite is located inside the first-stage carbonate cement and is dispersed fine-grained ( Figure 6 -i). The content of this type of pyrite is relatively low, making it difficult to accurately calculate its content and distribution. The second-stage pyrite is generally associated with the fourth-stage carbonate cement phase ( Figure 6 -j). This type of pyrite is generally developed in sandstone reservoirs with a burial depth greater than 3500m. Within the depth range of 3500m-4000m, the content of pyrite in this period increases with increasing burial depth; when the burial depth is greater than 4500m, the content of pyrite in this period decreases with increasing burial depth.

[0094] Three different types of chlorite and illite occur in the study area. The first type of chlorite and illite develops in a rim-like pattern between the first and second carbonate cements. The second type of chlorite and illite primarily manifests as illitization and chloritization of authigenic kaolinite from various phases. The third type of chlorite and illite is primarily the product of authigenic precipitation from pore fluids. The third type of illite and chlorite minerals are primarily located outside or coexist with the fourth type of carbonate minerals. The chlorite and illite content in sandstone reservoirs increases rapidly at depths greater than 3500m, and their content increases with depth.

[0095] S302: Select the samples of carbonate cements and authigenic quartz with good crystal form and large crystals observed in S301 to grind inclusion thin sections (to ensure that the carbonates and authigenic quartz cements of different occurrences observed in S301 are basically covered in the samples and have a certain number of observation samples); observe under a microscope, select salt water inclusions of appropriate size and shape (round or elliptical, >5um) in carbonate cements of different occurrences and authigenic quartz, and use the hot and cold stage programmed temperature technology to test the uniform temperature of salt water inclusions in different diagenetic products.

[0096] During the programmed temperature increase, when the temperature is below 70°C, the temperature is increased in 10°C increments; when the temperature is above 70°C, the temperature is increased in 5°C increments. After each temperature increase, the sample is stopped for 20 minutes to observe whether the inclusions have reached a homogeneous state. The homogenization temperature of the salt water inclusions in the mineral is the formation temperature of the mineral.

[0097] The homogenization temperature of the brine inclusions in the first-stage authigenic quartz cement (Q1, quartz secondary enlargement) is 75-120°C; the homogenization temperature of the brine inclusions in the second-stage authigenic quartz cement (Q2, euhedral quartz) is 110-145°C.

[0098] The first-stage carbonate cements are micrite and do not contain brine inclusions whose temperature can be measured; the second-stage carbonate cements are 60-90°C; the third-stage carbonate cements are 80-130°C; the fourth-stage carbonate cements are 120-150°C;

[0099] S303: Grind the remaining samples of the samples selected in S302 to prepare isotope test slices; using the results observed in S301, select samples with a certain occurrence of carbonate cement accounting for more than 70% of the total carbonate cement content and with large single crystals or aggregates of carbonate cement (>50um), and grind isotope test slices. Observe the isotope test slices under a microscope, and select carbonate cements of different occurrences, gasify them using laser ablation technology, and the generated gas enters the mass spectrometer to test the δ 13 C PDB (‰) and δ18 O PDB (‰). The results are shown in ( Figure 9 )

[0100] S304: Select samples that have been tested for saltwater inclusions (S302) and isotopes (S303) in carbonate cements of different occurrences, select the test results of homogenization temperature and isotopes of carbonate cement inclusions of the same sample and the same occurrence, and calculate the values ​​of the O isotopic characteristics of the paleopore fluid during the formation period of different carbonate cements:

[0101] The δ 18 O V-PDB / ‰ is about -4.8‰. The δ 18 O V-PDB / ‰ is similar to the first-stage carbonate cement, which is also about -4.8‰. The δ 18 O V-PDB / ‰ is obviously positive, about -2.0‰. The δ 18 O V-PDB / ‰ is about -1.5‰( Figure 12 ).

[0102] The δ 18 O V-PDB The calculated results of / ‰ are used in the O isotope characteristic differentiation formula proposed by Riedman I et al., 1977 to calculate the formation temperature of carbonate cements of different occurrences:

[0103] Research shows that the first phase of carbonate cement was formed at 16℃-45℃; the second phase was formed at 20℃-85℃; the third phase was formed at 80℃-120℃; and the fourth phase was formed at 95℃-160℃.

[0104] Based on the homogenization temperature of brine inclusions and isotope test data, the formation temperatures of carbonate cements and authigenic quartz of different occurrences were calculated:

[0105] The first phase of carbonate cement formed at 16-45°C; the second phase formed at 20-85°C; the third phase formed at 80-130°C; and the fourth phase formed at 95-160°C. The homogenization temperature of the brine inclusions in the first phase of authigenic quartz cement (Q1, quartz secondary enlargement) is 75-120°C; the homogenization temperature of the brine inclusions in the second phase of authigenic quartz cement (Q2, euhedral quartz) is 110-145°C.

[0106] S4: Based on the burial history, thermal history, formation temperature of typical diagenetic products (S304), and contact and replacement relationship of diagenetic products (S301), a diagenetic sequence diagram of the study area is established ( Figure 13 ).

[0107] According to the diagenetic sequence, read the diagenetic combination and the formation time and depth of different diagenetic combinations;

[0108] The first carbonate cement (Cc1) was formed at 47-40 Ma; at this time, the reservoir depth was <1000 m;

[0109] The second phase of carbonate cement (Cc2) was formed between 40 and 32 Ma; at this time, the reservoir was buried at a depth of 1000 to 2000 m.

[0110] Plagioclase dissolution, the first stage of authigenic quartz and the first stage of authigenic kaolinite were formed at 32-23 Ma; at this time, the reservoir was buried at a depth of 2000-3200 m;

[0111] The third stage carbonate cement (Cc3) was formed at 23-8 Ma; at this time, the reservoir was buried at a depth of 2500-3200 m (uplift);

[0112] The second stage of authigenic quartz and the second stage of kaolinite and large-scale kaolinization of potassium feldspar were formed at 8-4 Ma; at this time, the reservoir was buried at a depth of 3200-4000 m;

[0113] The fourth phase of carbonate cement (Cc4), the dissolution of large-scale rock fragments, the dissolution of anhydrite, and the second phase of kaolinite and pyrite filling were formed at 4-2 Ma; at this time, the reservoir was buried at a depth of 4000-4500 m.

[0114] Large-scale chloritization of kaolinite; authigenic chlorite and illite precipitation formed between 2 Ma and the present; at this time, the reservoir depth was >4500 m;

[0115] S5: Based on the depth range and formation time of different diagenetic products extracted in S4, the geological background of this stage, the properties of the formation water in the sandstone reservoir buried at this stage (S1), the content and changes of the skeleton grains in the sandstone reservoir buried at this stage (S2), the types and isotopic characteristics of the diagenetic products in the sandstone reservoir buried at this stage (S301, S303), and the O isotopic characteristics of the paleopore water at this stage (S304) are comprehensively considered to comprehensively determine the type of water-rock reaction that occurred in the reservoir at this stage and the properties of the paleopore fluid, and establish the corresponding relationship between the diagenetic products, the properties of the current formation water, the change pattern of the skeleton grains and the properties of the paleopore fluid;

[0116] For example, the second-stage carbonate cements were mainly formed within the burial range of 1000m-2000m; the formation water in this process was Class I formation water, mainly Cl - Concentration and Na+ +K + , Ca 2+ Mg 2+ The plasma concentration increases at a ratio of 1:1, which is mainly due to the evaporation and concentration of formation water. At the same time, the content and composition of various skeleton particles in this burial depth section do not change significantly, that is, there is no phenomenon of skeleton particles dissolving to provide materials for formation water and diagenetic products. At the same time, the isotopic characteristics of the second-stage carbonate cement show obvious characteristics of lake precipitation carbonate, and the δ 18 O PDB / ‰ is -4.8‰, which also falls within the typical lake water isotope range. In other words, the pore fluid in this process was lake water that underwent evaporation and concentration. However, simple evaporation and concentration cannot form carbonate cements with a surface porosity exceeding 7%. Considering that the compaction effect during the formation of this period of secondary carbonate cements caused the mudstone porosity to drop rapidly (from 60% to 10%-20%), a large amount of sedimentary water in the mudstone was discharged into the sandstone. At the same time, the solubility of carbonate minerals decreased due to the increase in temperature, and a large amount of calcite precipitated from the pore water. The concentration effect during the burial process and the mudstone compaction and drainage effect led to the formation of Cc2 and controlled its distribution characteristics. That is, the paleopore fluid represented by the second-stage carbonate (Cc2) was the lake water that had undergone concentration and was discharged by mudstone compaction.

[0117] Based on the same analytical approach, the corresponding relationships between the diagenetic mineral assemblage at different stages, the properties of current formation water, the variation patterns of skeleton particles and the properties of paleopore fluids were determined (Table 1); and the evolution process of paleopore fluids in the study area was restored ( Figure 14 ).

[0118] Table 1

[0119]

[0120]

[0121] The above embodiments are preferred implementation modes of the present invention, but the implementation modes of the present invention are not limited to the above embodiments. Any other changes, modifications, substitutions, combinations, and simplifications that do not deviate from the spirit and principles of the present invention should be considered as equivalent replacement methods and are included in the scope of protection of the present invention.

Claims

1. A method for inverting the evolution of paleopore fluid properties in a sedimentary basin by comprehensively applying multiple geological data, characterized in that: The following steps are involved: According to the depth range and formation time of different diagenetic products in the study area, the properties of formation water in the reservoir, the content and changes of skeleton particles, the types and isotopic characteristics of diagenetic products, the isotopic characteristics of paleopore water, the types of water-rock reactions occurring in the reservoir and the properties of paleopore fluids; Establish the corresponding relationship between diagenetic products, formation water properties, skeleton particle change patterns and paleo-pore fluid properties; and invert the evolution process of paleo-pore fluid in the study area based on the corresponding relationship.

2. The method according to claim 1, characterized in that The specific steps are: S1. Determine the formation water properties within the sandstone reservoirs in the study area; S2. Sandstone samples from different burial depths in the study area were selected to grind thin sections of casts and determine the variation of the content of various types of skeleton particles with depth; S3. Determine the types and isotopic characteristics of diagenetic products, determine the oxygen isotopic characteristics of pore water in the reservoir during the formation of different carbonate cements, and determine the formation temperatures of carbonate cements of different occurrences and authigenic quartz of different occurrences; S4. Determine the depth range and formation time of different occurrences and types of diagenetic products; S5. Determine the type of water-rock reaction and paleopore fluid properties within the reservoir based on the depth range and formation time of different diagenetic products in the study area, the geological background, the properties of formation water within the sandstone reservoir, the content and variation of skeletal grains, the types and isotopic characteristics of diagenetic products, and the O isotopic characteristics of paleopore water; Establish the corresponding relationship between diagenetic products, formation water properties, skeleton particle change patterns and paleo-pore fluid properties; and invert the evolution process of paleo-pore fluid in the study area based on the corresponding relationship.

3. The method according to claim 2, characterized in that In S1: S101. Determine the hydrochemical characteristics of the formation water in the study area, the content of the major ions it contains, and its variation trend with burial depth; S102. Classify the formation water in the reservoir and clarify the spatial distribution characteristics of different types of formation water; S103. Utilize the evaporation curves of different water bodies to determine the depletion-enrichment characteristics of constant particles in different types of formation water, and thus determine the properties of formation water in sandstone reservoirs.

4. The method according to claim 3, characterized in that In S101, the hydrochemical characteristics of the formation water are the total mineralization of the formation water and its variation trend with burial depth; In S102, the SPSS cluster analysis module is used to classify the formation water in the reservoir based on the constant ion content of the formation water. The specific method is as follows: Taking the milligram equivalents of various constant ions in different formation water samples as the numerical values, the formation water samples were analyzed using the cluster analysis module in SPSS20 software, and a dendrogram was made and classified.

5. The method according to claim 4, characterized in that: In S103, the specific method for determining the depletion-enrichment characteristics of constant particles in different types of formation water using the evaporation curves of different water bodies is as follows: Using chloride ions as a standard, a relationship diagram between chloride ion concentration and the constant ions was drawn, and evaporation-concentration curves of different water bodies were drawn at the same time to determine the loss-enrichment characteristics of constant particles in different types of formation water.

6. The method according to any one of claims 3 to 5, characterized in that: The constant ions include Na + , K + , Ca 2+ Mg 2 + 、SO4 2- 、HCO3 - .

7. The method according to claim 2, characterized in that: In S2, sandstone samples from different burial depths in the study area were selected to grind thin sections of castings and observed under a polarizing microscope to determine the variation pattern of the content of various types of skeleton particles with depth and clarify its spatial variation pattern; The rock skeleton particles in the sandstone reservoir mainly include: the proportion of quartz particles in the total particles, the proportion of feldspar particles in the total particles, the proportion of rock debris particles in the total particles, the proportion of sedimentary rock debris particles in the total rock debris particles, the proportion of igneous rock debris particles in the total rock debris particles, the proportion of igneous rock debris particles in the total rock debris particles, the proportion of potassium feldspar particles in the total feldspar, and the proportion of plagioclase particles in the total feldspar.

8. The method according to claim 2, characterized in that: In S3: S301. Based on polarizing microscope observation of cast thin sections, determine the type, content, occurrence, contact and replacement relationship, and formation sequence of diagenetic products within the sandstone reservoirs in the study area; and calculate the vertical distribution and spatial distribution characteristics of diagenetic products of different occurrences; S302. Select samples of carbonate cements of different occurrences and authigenic quartz cements of different occurrences, prepare inclusion thin sections, and determine the homogenization temperature of the inclusions, i.e., the formation temperature of the corresponding minerals; S303. Based on the observation results of S301 and S302, select typical samples for isotope thin section observation and determine the δ 13 C PDB and δ 18 O PDB ; S304. Based on S302 and S303, determine the O isotope characteristics of the pore water in the reservoir during the formation period of different carbonate cements, and determine the formation temperature of carbonate cements of different occurrences and the formation temperature of authigenic quartz of different occurrences.

9. The method according to claim 8, characterized in that In S303, the typical samples are all samples of carbonate cements observed and measured in S302, as well as samples in S301 in which the content of carbonate cements in a certain occurrence accounts for more than 70% of the total content of carbonate cements in the sample, and the single crystals or aggregates of the carbonate cements in the occurrence are larger than 50um.

10. The method according to claim 8, characterized in that: In S304, the method for determining the formation temperature of carbonate cements of different occurrences and the formation temperature of authigenic quartz of different occurrences is as follows: The δ 13 C PDB and δ 18 O PDB At the same time, the homogenization temperature of inclusions of carbonate cements of different occurrences and authigenic quartz of different occurrences was tested using the hot and cold stage programmed temperature method; Based on the δ 18 O PDB and δ 18 O PDB , comprehensively considering the homogenization temperature of brine inclusions in carbonate cements of different occurrences, determine the formation temperature of carbonate cements of different occurrences; based on the homogenization temperature of brine inclusions in authigenic quartz of different occurrences, determine the homogenization temperature of authigenic quartz of different occurrences.

11. The method according to claim 8, characterized in that: In S4, a burial history-thermal history map of the study area was established. Based on the homogenization temperatures of carbonate cements of different occurrences and authigenic quartz of different occurrences determined in S304 and the contact and replacement relationship of different diagenetic products determined in S301, the diagenetic sequence of the study area was established; the formation depth range and formation time of diagenetic products of different occurrences and types were determined.

12. A system for inverting the evolution of paleopore fluid properties in sedimentary basins by comprehensively applying multiple geological data, comprising a processor and a memory storing a computer program. It is characterized by: The processor is configured to execute the computer program to implement a method for inverting the evolution process of paleo-pore fluid properties in a sedimentary basin by comprehensively applying multiple geological data.

13. A computer storage medium having a computer program stored thereon, characterized in that ； When the computer program is executed by a processor, a method for inverting the evolution process of paleo-pore fluid properties in a sedimentary basin by comprehensively applying multiple geological data is implemented.

Images