CN112528226A - Deep carbonate reservoir pore evolution recovery method - Google Patents
Deep carbonate reservoir pore evolution recovery method Download PDFInfo
- Publication number
- CN112528226A CN112528226A CN202011227968.XA CN202011227968A CN112528226A CN 112528226 A CN112528226 A CN 112528226A CN 202011227968 A CN202011227968 A CN 202011227968A CN 112528226 A CN112528226 A CN 112528226A
- Authority
- CN
- China
- Prior art keywords
- stage
- cement
- carbonate
- burial
- pore
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 title claims abstract description 159
- 239000011148 porous material Substances 0.000 title claims abstract description 142
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000011084 recovery Methods 0.000 title claims abstract description 31
- 239000004568 cement Substances 0.000 claims abstract description 304
- 238000009933 burial Methods 0.000 claims abstract description 191
- 239000011435 rock Substances 0.000 claims abstract description 126
- 125000005587 carbonate group Chemical group 0.000 claims abstract description 42
- 230000003628 erosive effect Effects 0.000 claims abstract description 40
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 36
- 238000012360 testing method Methods 0.000 claims abstract description 12
- 238000012937 correction Methods 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 15
- 239000000843 powder Substances 0.000 claims description 15
- 230000009467 reduction Effects 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 11
- 238000005520 cutting process Methods 0.000 claims description 10
- 238000011161 development Methods 0.000 claims description 10
- 230000001973 epigenetic effect Effects 0.000 claims description 8
- 238000011065 in-situ storage Methods 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 4
- 238000005260 corrosion Methods 0.000 description 24
- 230000007797 corrosion Effects 0.000 description 24
- 230000000694 effects Effects 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 8
- 238000009529 body temperature measurement Methods 0.000 description 7
- 239000010459 dolomite Substances 0.000 description 7
- 229910000514 dolomite Inorganic materials 0.000 description 7
- 238000005266 casting Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 238000004090 dissolution Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 229910021532 Calcite Inorganic materials 0.000 description 2
- 229910001748 carbonate mineral Inorganic materials 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- INJRKJPEYSAMPD-UHFFFAOYSA-N aluminum;silicic acid;hydrate Chemical compound O.[Al].[Al].O[Si](O)(O)O INJRKJPEYSAMPD-UHFFFAOYSA-N 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000010443 kyanite Substances 0.000 description 1
- 229910052850 kyanite Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
- G01K11/30—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of the effect of a material on X-radiation, gamma radiation or particle radiation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/84—Systems specially adapted for particular applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Food Science & Technology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Evolutionary Biology (AREA)
- Medicinal Chemistry (AREA)
- Dispersion Chemistry (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Computational Biology (AREA)
- Algebra (AREA)
- Remote Sensing (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Geology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention provides a deep carbonate reservoir pore evolution recovery method. The method comprises the following steps: obtaining a work area representative rock sample; determining the stage of the carbonate cement, and carrying out isotope dating and cluster isotope testing on the carbonate cement in each stage to obtain the absolute age and the formation temperature of the carbonate cement in each stage; acquiring an ancient earth temperature model and a burial history model of a work area; correcting the burial history model by using the absolute age and the forming temperature of the subcarbonate cement of each stage and combining the ancient geothermal model of the work area to obtain a corrected burial history model, and determining the burial depth formed by the subcarbonate cement of each stage by combining the absolute age of the subcarbonate cement of each stage on the basis; and (3) counting the porosity of the existing pore surface of each rock sample, the area percentage of each stage of the secondary carbonate cement and the amount of each stage of the secondary erosion pore-increasing quantity based on the rock sample, and reconstructing the pore evolution history changing along with the buried depth by combining the buried depth formed by each stage of the secondary carbonate cement.
Description
Technical Field
The invention belongs to the technical field of carbonate rock oil and gas reservoir evaluation methods in petroleum and gas geological exploration, and particularly relates to a deep carbonate rock reservoir pore evolution recovery method.
Background
Deep carbonate rock has become a very important field of succession for oil and gas exploration in China, but due to the complex geological conditions of deep high temperature, high pressure and the like and the high chemical activity of deep carbonate rock, the deep carbonate rock is not known in the pore evolution process along with the depth change, the reservoir distribution rule is not known, and the evaluation of the deep carbonate rock reservoir has become one of the bottlenecks of deep oil and gas exploration. Therefore, deep carbonate rock pore evolution history recovery needs to be carried out, and a basis is provided for deep carbonate rock reservoir evaluation.
The pore evolution of carbonate reservoirs is mainly affected by erosion and cementationWith control, the predecessors made a great deal of research work on the aspects of diagenetic-pore evolution of carbonate rocks, but the diagenetic-pore evolution history is established qualitatively based on diagenetic action (mainly corrosion action and cementation action) and diagenetic sequence research under a microscope. Due to the absence of absolute age data and cluster isotopes (Δ) for the U-Pb isotope47) Mutual evidence of temperature data cannot ensure the reliability of a burial history model, cannot determine the burial depth of cementation and dissolution, and is difficult to realize the reconstruction of the pore evolution process of the carbonate rock along with the change of the burial depth. Moreover, not all cements and erosion cavities have definite mutual cutting relations, sometimes the cement, the erosion period and the diagenesis sequence are difficult to establish, and even the qualitative diagenesis-pore evolution history has great uncertainty.
Disclosure of Invention
The invention aims to provide a method capable of solving deep carbonate reservoir pore evolution history recovery, which realizes the reconstruction of the deep carbonate reservoir pore quantitative evolution history along with the change of the buried depth and provides a basis for evaluating the deep carbonate reservoir.
In order to achieve the above object, the present invention provides a deep carbonate reservoir pore evolution recovery method, wherein the method comprises:
obtaining a work area representative rock sample, wherein the characteristics of the work area representative rock sample comprise: the development of the holes of the rock sample, the filling of carbonate cement in the holes and the mutual cutting of the carbonate cement in the rock sample;
determining the stage of the carbonate cement in the rock sample, and carrying out isotope year measurement on the carbonate cement of each stage to obtain the absolute age of the carbonate cement of each stage; performing a cluster isotope (such as delta 47 temperature) test on the carbonate cement of each stage to obtain the formation temperature of the carbonate cement of each stage;
acquiring an ancient earth temperature model and a burial history model of a work area; correcting the burial history model of the work area by using the absolute age of the subcarbonate cement of each stage and the formation temperature of the subcarbonate cement of each stage and combining the ancient geothermal temperature model of the work area to obtain a burial history model after correction of the work area;
determining the buried depth formed by each stage of the sub-carbonate cement by using the absolute age of the sub-carbonate cement of each stage based on a buried history model after work area correction;
on the basis of the representative rock sample of the work area, counting the porosity of the existing pore surface of each rock sample, the area percentage of the subcarbonate cement at each stage (corresponding to the cementation effect to reduce the pore effect) and the amount of the subcarbonate cement at each stage (corresponding to the pore effect of the dissolution effect at each stage); and (3) reconstructing the evolution history of the diagenetic porosity (namely a curve of the diagenetic porosity along with the change of the burial depth) along with the burial depth of the subcarbonate cement at each stage based on the statistical result.
In the deep carbonate rock reservoir pore evolution recovery method, the rock sample with the characteristics of hole development, filling of carbonate cement in the hole and mutual intersection of the carbonate cement in the rock sample is easy to establish a complete and reliable diagenetic sequence and is suitable for year measurement and temperature measurement; preferably, the carbonate cement bond characteristics and the reciprocal intersection relationship of the representative rock sample of the work area are clear.
In the above deep carbonate reservoir pore evolution restoration method, preferably, the determining of the age of the carbonate cement intersected with each other in the rock sample is performed using a sample slice a made of a rock sample representative of the work area. More preferably, the thickness of the sample sheet a is 30 ± 3 μm.
In the above deep carbonate reservoir pore evolution restoration method, preferably, the isotope year measurement is performed by using a sample slice B made of a rock sample representative of the work area. More preferably, the thickness of the sample sheet B is 80 to 100 μm.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the performing of the cluster isotope test is performed by using a powder sample of each stage of the sub-carbonate cement.
In a specific embodiment, the method for recovering pore evolution of a deep carbonate reservoir comprises the following steps:
obtaining a work area representative rock sample, wherein the characteristics of the work area representative rock sample comprise: the development of the holes of the rock sample, the filling of carbonate cement in the holes and the mutual cutting of the carbonate cement in the rock sample;
respectively preparing at least 2 parallel samples corresponding to the representative rock samples of the obtained work area, preparing sample slices A and sample slices B of the representative rock samples by using the parallel samples, and reserving the residual parts of the parallel samples;
observing the carbonate cement of the sample slice A, and determining the period of the carbonate cement in the rock sample;
in the corresponding sample slice B, the carbonate cements of each stage corresponding to the carbonate cements of each stage in the sample slice A are defined, and the absolute age of the carbonate cements of each stage is obtained by isotope year measurement;
in the corresponding parallel sample residual part, obtaining a powder sample of carbonate cement of each stage corresponding to the carbonate cement of each stage in the sample slice A, and carrying out cluster isotope test to obtain the formation temperature of the carbonate cement of each stage;
acquiring an ancient earth temperature model and a burial history model of a work area; correcting the burial history model of the work area by using the absolute age of the subcarbonate cement of each stage and the formation temperature of the subcarbonate cement of each stage and combining the ancient geothermal temperature model of the work area to obtain a burial history model after correction of the work area;
determining the burial depth of each stage of the sub-carbonate cement by using the absolute age of the sub-carbonate cement of each stage based on a buried history model after work area correction;
based on the representative rock samples of each work area, counting the porosity of the existing pore surface of each rock sample, the area percentage of the sub-carbonate cement at each stage and the sub-erosion pore-increasing amount at each stage; based on the statistical result, combining the burial depth formed by the secondary carbonate cement at each stage, and reconstructing the pore evolution history along with the change of the burial depth;
preferably, the thickness of the sample sheet a is 30 ± 3 μm;
preferably, the sample slice a has a diameter of 1.5-2.5 cm;
preferably, the thickness of the sample sheet B is 80 to 100 μm;
preferably, the sample slice B has a diameter of 1.5-2.5 cm;
preferably, the separately preparing at least 2 parallel samples corresponding to each representative rock sample is performed by: cutting each representative rock sample into a cylinder with the diameter of 1.5-2.5cm and the thickness of 0.8cm, and preparing 2 parallel samples along two sides of a section;
preferably, the mirror image similarity of the sample slice A and the sample slice B is not lower than 90%; the consistency of the carbonate cement stages for year measurement and temperature measurement, and the one-to-one correspondence of age data and temperature data are ensured.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the determining the age of the carbonate cement in the rock sample includes:
establishing a complete and reliable diagenetic sequence according to the mutual intersection relationship of the carbonate cements, and determining the period of the carbonate cements with the mutual intersection relationship;
carbonate cements that do not have an inter-relationship, as a single stage.
In the deep carbonate reservoir pore evolution recovery method, preferably, the isotope dating is performed by using a laser in-situ U-Pb isotope dating mode.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the mass of the powder sample is not less than 10 mg.
In the above deep carbonate reservoir pore evolution recovery method, the powder sample may be obtained using a microdriller.
In the deep carbonate reservoir pore evolution recovery method, the work area burial history model can be obtained according to the conventional technical means in the field, for example, the work area ancient geothermal model and the burial history model are established according to the regional geological background, the drilling and seismic data.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the correcting the work area burial history model by using the absolute age of each stage of the sub-carbonate cement and the formation temperature of each stage of the sub-carbonate cement in combination with the work area ancient geothermal model includes:
the absolute age of each stage of the secondary carbonate cement is put into a burial history model to obtain the first burial depth of each stage of the secondary carbonate cement, and the formation temperature of each stage of the secondary carbonate cement is used for calculating the second burial depth of each stage of the secondary carbonate cement according to an ancient geothermal model;
if the first burial depth of the secondary carbonate cement at each stage is inconsistent with the second burial depth, the burial history model is unreliable, and the burial history curve is modified to ensure that the burial depth of the secondary carbonate cement at each stage is positioned at the second burial depth, so that the burial history curve after the work area correction is obtained;
if the first burial depth of the subcarbonate cement at each stage is consistent with the second burial depth, the absolute age of the subcarbonate cement at each stage and the forming temperature of the subcarbonate cement at each stage are considered to form a mutual evidence-based relationship, the burial history curve is reliable, and the burial history curve is used as a burial curve model after correction of a work area.
In the deep carbonate reservoir pore evolution recovery method, the determination of the burial depth of each stage of the sub-carbonate cement by using the absolute age of each stage of the sub-carbonate cement based on the burial history model after work area correction can be realized by the following method: and acquiring the burial depth of the subcarbonate cement at each stage based on the absolute age point of the subcarbonate cement at each stage to the corrected burial history model.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the reconstructing a diagenetic pore evolution history varying with the burial depth based on the statistical results and combined with the burial depth formed by the sub-carbonate cements at each stage includes: the initial porosity for the epigenetic stage is determined based on the porosity of the existing pore face and the area percentage of the subcarbonate cement at each stage, and the sum of the porosity of the existing pore face and the area percentage of the subcarbonate cement at each stage is the initial porosity for the epigenetic stage. To improve the representativeness and accuracy of the statistics, the initial porosity of the epigenetic period can be averaged based on the initial porosities of the plurality of rock samples.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the reconstructing a diagenetic pore evolution history varying with the burial depth based on the statistical results and combined with the burial depth formed by the sub-carbonate cements at each stage includes: the amount of porosity reduction for each stage of cementation was determined based on the area percentage of the sub-carbonate cement at each stage. To improve the representativeness and accuracy of the statistics, the cementing hole reducing amount of each phase of the multiple rock samples can be averaged to be used as the cementing hole reducing amount of the phase.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the reconstructing a diagenetic pore evolution history varying with the burial depth based on the statistical results and combined with the burial depth of each stage of the sub-carbonate cement comprises: and determining the corresponding buried depth of each stage of cementation subtractive hole based on the buried depth of each stage of secondary carbonate cement, wherein the buried depth of each stage of secondary carbonate cement is regarded as the buried depth when the pore is filled, namely the buried depth corresponding to the stage of cementation subtractive hole.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the reconstructing a diagenetic pore evolution history varying with the burial depth based on the statistical results and combined with the burial depth formed by the sub-carbonate cements at each stage includes: the porosity after each stage of cementation debulking is determined based on the existing porosity of the porosity face, the area percentage of each stage of sub-carbonate cement, the porosity after each stage of cementation debulking being equal to the sum of the existing porosity of the porosity face and the area percentage of each stage of sub-carbonate cement minus the total area percentage of the stage and the previous stage of sub-carbonate cement. To improve the representativeness and accuracy of the statistics, the porosity after cementing and pore reducing of a plurality of rock samples can be averaged to be used as the porosity after cementing and pore reducing of the period.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the reconstructing a diagenetic pore evolution history varying with the burial depth based on the statistical results and combined with the burial depth formed by the sub-carbonate cements at each stage includes: and determining the erosion pore-increasing amount of each stage based on the erosion pore-increasing amounts of the stages. To improve the representativeness and accuracy of the statistics, the erosion pore volume of each phase of the multiple rock samples can be averaged to be used as the erosion pore volume of the phase.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the reconstructing a diagenetic pore evolution history varying with the burial depth based on the statistical results and combined with the burial depth of each stage of the sub-carbonate cement comprises: and determining the depth corresponding to each stage of erosion increase hole based on the buried depth of the sub-carbonate cement of each stage, wherein the depth corresponding to each stage of erosion increase hole is between the buried depth of the eroded carbonate cement and the buried depth of the next stage of carbonate cement.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the reconstructing a diagenetic pore evolution history varying with the burial depth based on the statistical results and combined with the burial depth formed by the sub-carbonate cements at each stage includes: the porosity after each stage of erosion and pore enlargement is determined based on the porosity of the existing pore face, the area percentage of the subcarbonate cement at each stage and the erosion and pore enlargement amount of each stage, and the porosity after each stage of cementation and pore reduction is the sum of the porosity of the existing pore face and the area percentage of the subcarbonate cement at each stage-the total area percentage of the subcarbonate cement at the stage and the previous stage + the erosion and pore enlargement amount of the subcarbonate cement at the stage and the previous stage. In order to improve the representativeness and the accuracy of statistics, the porosity after the erosion and the pore growth of a plurality of rock samples can be averaged to be used as the porosity after the erosion and the pore growth of the period.
In the above deep carbonate reservoir pore evolution recovery method, preferably, the reconstructing a diagenetic pore evolution history varying with the burial depth based on the statistical results and combined with the burial depth formed by the sub-carbonate cements at each stage includes: the final porosity is determined based on the existing porosity of the pore face. To improve statistical representation and accuracy, the final porosity may be averaged based on the existing porosity of the porosity face of the plurality of rock samples.
In one embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of sub-carbonate cement formation, comprises:
determining the first point of the curve, namely the initial porosity (P1), which is the sum of the existing pore face porosity and the area percentage of each stage of the secondary carbonate cement, and the depth is the apparent stage depth (D1 ═ 0); preferably, the initial porosity is finally determined based on the average of the initial porosities of the respective rock samples;
determining a second point of the curve, namely a first-stage bond post-porosity-reduction (P2), wherein the calculation method is that the first-stage bond post-porosity-reduction (P2) is equal to the initial porosity (P1) -the first diagenetic-stage bond post-porosity-reduction amount (A1), and the depth is obtained by casting points in the burial history model according to the age of the first-stage bond (D2);
determining a third point of the curve, namely the porosity after corrosion and hole increase in the first stage (P3), wherein the calculation method is that the porosity after corrosion and hole increase in the first stage (P3) is equal to the porosity after cementation and hole decrease in the first stage (P2) + the corrosion and hole increase amount of the first diagenetic stage (B1), and the depth is between the burial depth of the cement in the first stage and the burial depth of the cement in the second stage (D3);
determining the fourth point of the curve, namely the porosity after second-stage cementation and hole reduction (P4), wherein the calculation method is that the porosity after second-stage cementation and hole reduction (P4) is equal to the porosity after first-stage corrosion and hole increase (P3) -the second-diagenetic-stage cementation hole reduction amount (A2), and the age of the first-stage cementation is the depth obtained by casting points in the burial history model (D4);
determining the fifth point of the curve, namely the porosity after the second-stage corrosion hole increase (P5), wherein the calculation method is that the porosity after the second-stage corrosion hole increase (P5) is equal to the porosity after the second-stage cementation hole decrease (P4) + the corrosion hole increase amount of the second diagenetic stage (B2), and the depth is between the second-stage cement buried depth and the third-stage cement buried depth (D5);
according to the stage of cementing material in the rock sample, if the cementing material in the third stage and the fourth stage also exists, the control points of the curve can be sequentially increased according to the method.
In the technical scheme provided by the invention, the absolute age data of the carbonate cement isotope and the temperature data of the cluster isotope (delta 47) are obtained by utilizing the isotope dating technology in combination with the cluster isotope (delta 47) temperature measuring technology, and a reliable burying history model is established based on the two data, so that the establishment of the carbonate reservoir pore evolution history changing along with the depth becomes possible. Specifically, according to the technical scheme provided by the invention, on the basis of measuring the absolute age of the isotope of the carbonate cement and measuring the temperature of the cluster isotope (delta 47), an ancient geothermal model and a burial history model are corrected, a reliable tectonic-burial history and basin thermal history geological background is provided for diagenesis-pore evolution history recovery, the absolute age of the carbonate cement is put on the burial history model, the depth of cementation and corrosion is determined, the quantitative pore evolution history reconstruction of the deep carbonate rock along with the change of the burial depth is realized, and a basis is provided for evaluating a deep carbonate rock reservoir. The technical scheme provided by the invention provides a means for recovering the quantitative pore evolution history of the deep carbonate reservoir and a technical means for evaluating the effectiveness of the deep carbonate reservoir.
Drawings
Fig. 1 is a flowchart of a deep carbonate reservoir pore evolution recovery method provided in embodiment 1 of the present invention.
FIG. 2A is a graph of characteristics and diagenesis sequence of a Takenorthwest seismic denier system Qigonglake dolomitic reservoir sample Q-56-1 in example 1 of the present invention.
FIG. 2B is a graph of characteristics and diagenesis sequence of a Takenorthwest seismic denier system Qigonglake dolomitic reservoir sample Q-56-1 in example 1 of the present invention.
FIG. 3A is a graph of characteristics and diagenesis sequence of a Takenorthwest seismic denier system Qigonglake dolomitic reservoir sample Q-58-1-1 in example 1 of the present invention.
FIG. 3B is a graph of characteristics and diagenesis sequence of a Takenorthwest seismic denier system Qigonglake dolomitic reservoir sample Q-58-1-1 in example 1 of the present invention.
FIG. 4 is a graph of characteristics and diagenesis sequence of a Takenorthwest seismic denier system Qigonglake dolomitic reservoir sample Q-58-1-2 in example 1 of the present invention.
FIG. 5A is a graph of characteristics and diagenesis sequence of a Takenorthwest seismic denier system Qigonglake dolomitic reservoir sample Q-76-1 in example 1 of the present invention.
FIG. 5B is a graph of characteristics and diagenesis sequence of a Takenorthwest seismic denier system Qigonglake dolomitic reservoir sample Q-76-1 in example 1 of the present invention.
FIG. 6A is a graph of characteristics and diagenesis sequence of a Takenorthwest seismic denier system Qigonglake dolomitic reservoir sample Q-151-1 in example 1 of the present invention.
FIG. 6B is a graph of characteristics and diagenesis sequence of a Takenorthwest seismic denier system Qigonglake dolomitic reservoir sample Q-151-1 in example 1 of the present invention.
FIG. 7 is a graph of the King Gebraker set burial history of the tower in the northwest seismic system reconstructed in example 1 of the present invention.
Fig. 8 is a reservoir pore evolution curve diagram of the qigbulake group carbonate rock of the northwest seismic denier system of the tower reconstructed in the embodiment 1 of the invention along with the change of the depth.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be described in detail and completely with reference to the drawings in the embodiments of the present invention. It is to be understood that the embodiments described are only a few embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
An embodiment of the invention provides a deep carbonate reservoir pore evolution recovery method, wherein the method comprises the following steps:
obtaining a work area representative rock sample, wherein the characteristics of the work area representative rock sample comprise: the development of the holes of the rock sample, the filling of carbonate cement in the holes and the mutual cutting of the carbonate cement in the rock sample;
determining the stage of the carbonate cement in the rock sample, and carrying out isotope year measurement on the carbonate cement of each stage to obtain the absolute age of the carbonate cement of each stage; performing a cluster isotope (such as delta 47 temperature) test on the carbonate cement of each stage to obtain the formation temperature of the carbonate cement of each stage;
acquiring an ancient earth temperature model and a burial history model of a work area; correcting the burial history model of the work area by using the absolute age of the subcarbonate cement of each stage and the formation temperature of the subcarbonate cement of each stage and combining the ancient geothermal temperature model of the work area to obtain a burial history model after correction of the work area;
determining the buried depth formed by each stage of the sub-carbonate cement by using the absolute age of the sub-carbonate cement of each stage based on a buried history model after work area correction;
on the basis of the representative rock sample of the work area, counting the porosity of the existing pore surface of each rock sample, the area percentage of the subcarbonate cement at each stage (corresponding to the cementation effect to reduce the pore effect) and the amount of the subcarbonate cement at each stage (corresponding to the pore effect of the dissolution effect at each stage); and (3) reconstructing the evolution history of the diagenetic porosity (namely a curve of the diagenetic porosity along with the change of the burial depth) along with the burial depth of the subcarbonate cement at each stage based on the statistical result.
The rock sample with the characteristics of hole development, filling of carbonate cement in the hole and mutual intersection of the carbonate cement in the rock sample is easy to establish a complete and reliable diagenetic sequence and is suitable for year measurement and temperature measurement; further, the carbonate cement bond characteristics and the mutual intersection relationship of representative rock samples of the work area are clear.
In a preferred embodiment, the determination of the age of the intercleaved carbonate cement in the rock sample is performed using a sample slice a made from a representative rock sample for the work area. Further, the thickness of the sample sheet A was 30. + -. 3. mu.m.
In a preferred embodiment, the isotope year measurement is performed using a sample slice B made from a rock sample representative of the work area. Further, the thickness of the sample sheet B was 80 to 100. mu.m.
In a preferred embodiment, the cluster isotope test is performed using powder samples of each stage of a subcarbonate cement.
In a preferred embodiment, the method for recovering pore evolution of a deep carbonate reservoir comprises:
obtaining a work area representative rock sample, wherein the characteristics of the work area representative rock sample comprise: the development of the holes of the rock sample, the filling of carbonate cement in the holes and the mutual cutting of the carbonate cement in the rock sample;
respectively preparing at least 2 parallel samples corresponding to the representative rock samples of the obtained work area, preparing sample slices A and sample slices B of the representative rock samples by using the parallel samples, and reserving the residual parts of the parallel samples;
observing the carbonate cement of the sample slice A, and determining the period of the carbonate cement in the rock sample;
in the corresponding sample slice B, the carbonate cements of each stage corresponding to the carbonate cements of each stage in the sample slice A are defined, and the absolute age of the carbonate cements of each stage is obtained by isotope year measurement;
in the corresponding parallel sample residual part, obtaining a powder sample of carbonate cement of each stage corresponding to the carbonate cement of each stage in the sample slice A, and carrying out cluster isotope test to obtain the formation temperature of the carbonate cement of each stage;
acquiring an ancient earth temperature model and a burial history model of a work area; correcting the burial history model of the work area by using the absolute age of the subcarbonate cement of each stage and the formation temperature of the subcarbonate cement of each stage and combining the ancient geothermal temperature model of the work area to obtain a burial history model after correction of the work area;
determining the burial depth of each stage of the sub-carbonate cement by using the absolute age of the sub-carbonate cement of each stage based on a buried history model after work area correction;
based on the representative rock samples of each work area, counting the porosity of the existing pore surface of each rock sample, the area percentage of the sub-carbonate cement at each stage and the sub-erosion pore-increasing amount at each stage; based on the statistical result, combining the burial depth formed by the secondary carbonate cement at each stage, and reconstructing the pore evolution history along with the change of the burial depth;
further, the thickness of the sample sheet a was 30 ± 3 μm;
further, the diameter of the sample sheet A is 1.5-2.5 cm;
further, the thickness of the sample sheet B is 80 to 100 μm;
further, the diameter of the sample sheet B is 1.5-2.5 cm;
further, the preparation of at least 2 parallel samples corresponding to each representative rock sample, respectively, was carried out by: cutting each representative rock sample into a cylinder with the diameter of 1.5-2.5cm and the thickness of 0.8cm, and preparing 2 parallel samples along two sides of a section;
further, the mirror image similarity of the sample slice A and the sample slice B is not lower than 90%; the consistency of the carbonate cement stages for year measurement and temperature measurement, and the one-to-one correspondence of age data and temperature data are ensured.
In a preferred embodiment, defining the stage of carbonate cement in the rock sample comprises:
establishing a complete and reliable diagenetic sequence according to the mutual intersection relationship of the carbonate cements, and determining the period of the carbonate cements with the mutual intersection relationship;
carbonate cements that do not have an inter-relationship, as a single stage.
In the deep carbonate reservoir pore evolution recovery method, preferably, the isotope dating is performed by using a laser in-situ U-Pb isotope dating mode.
In a preferred embodiment, the mass of the powder-like is not less than 10 mg.
In a preferred embodiment, the powder sample may be obtained using a microdriller.
In a preferred embodiment, the acquisition of the burial history model of the work area can be achieved according to the conventional technical means in the field, for example, an ancient earth temperature model and a burial history model of the work area are established according to the geological background, well drilling and seismic data of the area.
In a preferred embodiment, the correcting the model of the burial history of the work area using the absolute age of each stage of the sub-carbonate cement and the formation temperature of each stage of the sub-carbonate cement in combination with the paleogeothermal model of the work area comprises:
the absolute age of each stage of the secondary carbonate cement is put into a burial history model to obtain the first burial depth of each stage of the secondary carbonate cement, and the formation temperature of each stage of the secondary carbonate cement is used for calculating the second burial depth of each stage of the secondary carbonate cement according to an ancient geothermal model;
if the first burial depth of the secondary carbonate cement at each stage is inconsistent with the second burial depth, the burial history model is unreliable, and the burial history curve is modified to ensure that the burial depth of the secondary carbonate cement at each stage is positioned at the second burial depth, so that the burial history curve after the work area correction is obtained;
if the first burial depth of the subcarbonate cement at each stage is consistent with the second burial depth, the absolute age of the subcarbonate cement at each stage and the forming temperature of the subcarbonate cement at each stage are considered to form a mutual evidence-based relationship, the burial history curve is reliable, and the burial history curve is used as a burial curve model after correction of a work area.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of sub-carbonate cement formation, comprises: the depth of the carbonate cement at which the pore is filled is considered to be the depth of the burial at which the pore is filled, and the depth at which the dissolution pore is formed is considered to be between the depth of the burial of the eroded carbonate cement and the burial depth of the next phase of the carbonate cement.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of sub-carbonate cement formation, comprises: the initial porosity for the epigenetic stage is determined based on the porosity of the existing pore face and the area percentage of the subcarbonate cement at each stage, and the sum of the porosity of the existing pore face and the area percentage of the subcarbonate cement at each stage is the initial porosity for the epigenetic stage. To improve the representativeness and accuracy of the statistics, the initial porosity of the epigenetic period can be averaged based on the initial porosities of the plurality of rock samples.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of sub-carbonate cement formation, comprises: the amount of porosity reduction for each stage of cementation was determined based on the area percentage of the sub-carbonate cement at each stage. To improve the representativeness and accuracy of the statistics, the cementing hole reducing amount of each phase of the multiple rock samples can be averaged to be used as the cementing hole reducing amount of the phase.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of the secondary carbonate cement, comprises: and determining the corresponding buried depth of each stage of cementation subtractive hole based on the buried depth of each stage of secondary carbonate cement, wherein the buried depth of each stage of secondary carbonate cement is regarded as the buried depth when the pore is filled, namely the buried depth corresponding to the stage of cementation subtractive hole.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of sub-carbonate cement formation, comprises: the porosity after each stage of cementation debulking is determined based on the existing porosity of the porosity face, the area percentage of each stage of sub-carbonate cement, the porosity after each stage of cementation debulking being equal to the sum of the existing porosity of the porosity face and the area percentage of each stage of sub-carbonate cement minus the total area percentage of the stage and the previous stage of sub-carbonate cement. To improve the representativeness and accuracy of the statistics, the porosity after cementing and pore reducing of a plurality of rock samples can be averaged to be used as the porosity after cementing and pore reducing of the period.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of sub-carbonate cement formation, comprises: and determining the erosion pore-increasing amount of each stage based on the erosion pore-increasing amounts of the stages. To improve the representativeness and accuracy of the statistics, the erosion pore volume of each phase of the multiple rock samples can be averaged to be used as the erosion pore volume of the phase.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of the secondary carbonate cement, comprises: and determining the depth corresponding to each stage of erosion increase hole based on the buried depth of the sub-carbonate cement of each stage, wherein the depth corresponding to each stage of erosion increase hole is between the buried depth of the eroded carbonate cement and the buried depth of the next stage of carbonate cement.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of sub-carbonate cement formation, comprises: the porosity after each stage of erosion and pore enlargement is determined based on the porosity of the existing pore face, the area percentage of the subcarbonate cement at each stage and the erosion and pore enlargement amount of each stage, and the porosity after each stage of cementation and pore reduction is the sum of the porosity of the existing pore face and the area percentage of the subcarbonate cement at each stage-the total area percentage of the subcarbonate cement at the stage and the previous stage + the erosion and pore enlargement amount of the subcarbonate cement at the stage and the previous stage. In order to improve the representativeness and the accuracy of statistics, the porosity after the erosion and the pore growth of a plurality of rock samples can be averaged to be used as the porosity after the erosion and the pore growth of the period.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of sub-carbonate cement formation, comprises: the final porosity is determined based on the existing porosity of the pore face. To improve statistical representation and accuracy, the final porosity may be averaged based on the existing porosity of the porosity face of the plurality of rock samples.
In a preferred embodiment, reconstructing the history of the evolution of diagenetic pores as a function of the depth of burial, based on the statistical results, in combination with the depth of burial for each stage of sub-carbonate cement formation, comprises:
determining the first point of the curve, namely the initial porosity (P1), which is the sum of the existing pore face porosity and the area percentage of each stage of the sub-carbonate cement, and the depth is the apparent stage depth (D1 is 0); preferably, the initial porosity is finally determined based on the average of the initial porosities of the respective rock samples;
determining a second point of the curve, namely a first-stage bond post-porosity-reduction (P2), wherein the calculation method is that the first-stage bond post-porosity-reduction (P2) is equal to the initial porosity (P1) -the first diagenetic-stage bond post-porosity-reduction amount (A1), and the depth is obtained by casting points in the burial history model according to the age of the first-stage bond (D2);
determining a third point of the curve, namely the porosity after corrosion and hole increase in the first stage (P3), wherein the calculation method is that the porosity after corrosion and hole increase in the first stage (P3) is equal to the porosity after cementation and hole decrease in the first stage (P2) + the corrosion and hole increase amount of the first diagenetic stage (B1), and the depth is between the burial depth of the cement in the first stage and the burial depth of the cement in the second stage (D3);
determining the fourth point of the curve, namely the porosity after second-stage cementation and hole reduction (P4), wherein the calculation method is that the porosity after second-stage cementation and hole reduction (P4) is equal to the porosity after first-stage corrosion and hole increase (P3) -the second-diagenetic-stage cementation hole reduction amount (A2), and the age of the first-stage cementation is the depth obtained by casting points in the burial history model (D4);
determining the fifth point of the curve, namely the porosity after the second-stage corrosion hole increase (P5), wherein the calculation method is that the porosity after the second-stage corrosion hole increase (P5) is equal to the porosity after the second-stage cementation hole decrease (P4) + the corrosion hole increase amount of the second diagenetic stage (B2), and the depth is between the second-stage cement buried depth and the third-stage cement buried depth (D5);
according to the stage of cementing material in the rock sample, if the cementing material in the third stage and the fourth stage also exists, the control points of the curve can be sequentially increased according to the method.
Example 1
The embodiment provides a deep carbonate reservoir pore evolution recovery method, which is used for carrying out pore evolution curve reconstruction on a tower northwest seismic denier system odd-Gebraker group dolomite reservoir and providing a technical means for evaluating the effectiveness of the deep carbonate reservoir, and as shown in fig. 1, the method specifically comprises the following steps:
step S1: obtaining a work area representative rock sample, wherein the characteristics of the work area representative rock sample comprise: the development of the holes of the rock sample, the filling of carbonate cement in the holes and the mutual cutting of the carbonate cement in the rock sample;
the rock sample with the characteristics of hole development, filling of carbonate cement in the hole and mutual intersection of the carbonate cement in the rock sample is easy to establish a complete and reliable diagenetic sequence and is suitable for year measurement and temperature measurement.
Step S2: respectively preparing 2 parallel samples, namely a parallel sample A and a parallel sample B, corresponding to each representative rock sample of the obtained work area;
specifically, each representative rock sample was cut into a cylinder having a diameter of 1.5 to 2.5cm and a thickness of 0.8cm, and 2 parallel samples, one parallel sample a and the other parallel sample B, were prepared along both sides of the cut surface.
Step S3: preparing a sample slice A of each representative rock sample from the parallel sample A corresponding to each representative rock sample in the work area, preparing a sample slice B of each representative rock sample from the parallel sample B corresponding to each representative rock sample, and reserving the residual part of each parallel sample; wherein the thickness of the sample sheet A was 30 μm and the thickness of the sample sheet B was 100. mu.m.
Step S4: carrying out mirror image relation consistency screening on the sample slice A and the sample slice B of each representative rock sample;
specifically, a microscope is used for respectively observing the sample slice A and the sample slice B of each representative rock sample in a mirror image relationship, if the similarity of the mirror image relationship of the sample slice A and the sample slice B of a certain sample is not lower than 90%, the sample slice A and the sample slice B are reserved, otherwise, the sample is removed;
the mirror image corresponding relation research of the thin slice A and the thin slice B ensures the consistency of the carbonate cement period for year measurement and temperature measurement and the one-to-one correspondence of age data and temperature data.
Step S5: observing the carbonate cement of the sample slice A, and determining the period of the carbonate cement in the rock sample;
specifically, the sample sheet A is subjected to carbonate cement observation, and the type, characteristics, duration and the like of the carbonate cement are mainly observed; establishing a complete and reliable diagenetic sequence according to the mutual intersection relationship of the carbonate cements, and determining the period of the carbonate cements with the mutual intersection relationship; carbonate cements that do not have a mutual cross-cut relationship, as a single stage;
the structural components of the dolomite of the odd-Gebraker group of the northwest seismic denier system of the tower are as follows from early to late in sequence: firstly, the surrounding rock → fibrillar annular edge dolomite → foliated dolomite → fine powder grain dolomite → middle grain dolomite → hydrothermal dolomite and quartz (as shown in fig. 2A-fig. 6B). There is also a period of carbonate cement, in which the calcite structural component is not in an intersecting relationship with other structural components, as a single period, namely calcite filled in the fracture.
Step S6: in the corresponding sample slice B, the carbonate cements of each stage corresponding to the carbonate cements of each stage in the sample slice A are defined, and laser in-situ U-Pb isotope dating is carried out to obtain the absolute age of the carbonate cements of each stage; the results are shown in Table 1;
and (3) laser in-situ U-Pb isotope dating according to the specifications and requirements of the carbonate mineral laser in-situ U-Pb isotope dating technology.
TABLE 1 temperatures of different structural component cluster isotopes (Delta 47) of the Kigelac dolomites of the northwest seismic denier system of the tower
Step S7: in the corresponding parallel sample residual part, powder samples of carbonate cements of each stage corresponding to the carbonate cements of each stage in the sample slice A are drilled by a micro drill (the mass of each powder sample is 10mg), and cluster isotope (delta 47 temperature) test is carried out to obtain the formation temperature of the carbonate cements of each stage; the results are shown in Table 1;
if the residual part of the parallel sample corresponding to one rock sample cannot drill enough powder samples, the synchronous powder samples can be drilled through the residual parts of the parallel samples corresponding to a plurality of rock samples, so that the problem of insufficient powder sample amount is solved;
the cluster isotope (delta 47 temperature) test is carried out according to the specification and the requirement of the carbonate mineral cluster isotope (delta 47) temperature measurement technology.
Step S8: acquiring an ancient earth temperature model and a burial history model of a work area; correcting the burying history model by using the absolute age of each stage of secondary carbonate cement and the forming temperature of each stage of secondary carbonate cement to obtain a burying history model after the work area is corrected (as shown by a curve B in figure 7);
specifically, the method comprises the following steps:
8.1, according to the geological background, the drilling well and the seismic data of the region, primarily establishing a burial history model of the work area (shown as a curve A in figure 7), and acquiring an ancient geothermal model of the work area (in the embodiment, an ancient geothermal gradient);
8.2, casting the absolute age of the subcarbonate cement of each stage to a burial history model to obtain a first burial depth of the subcarbonate cement of each stage, and calculating the formation temperature of the subcarbonate cement of each stage according to the paleoterrestrial temperature gradient to obtain a second burial depth of the subcarbonate cement of each stage;
if the first burial depth of the secondary carbonate cement at each stage is inconsistent with the second burial depth, the burial history model is unreliable, and the burial history curve is modified to ensure that the burial depth of the secondary carbonate cement at each stage is positioned at the second burial depth, so that the burial history curve after the work area correction is obtained;
if the first burial depth of the subcarbonate cement at each stage is consistent with the second burial depth, the absolute age of the subcarbonate cement at each stage and the forming temperature of the subcarbonate cement at each stage are considered to form a mutual evidence-based relationship, the burial history curve is reliable, and the burial history curve is used as a burial curve model after correction of a work area.
The temperature gradient of the cambrian-early Orotan earth is 3.2 ℃/3.5 ℃/100m, the temperature gradient of the prime mover-mud basin earth is 3.0 ℃/100m, the temperature gradient of the carbolite-eclipse earth is 3.0 ℃/100m, the temperature gradient of the triage-chalk end earth is 2.5 ℃/100m, and the temperature gradient of the new generation earth is 2.0 ℃/100 m. By continuously correcting the absolute age of the U-Pb isotope and the temperature of the cluster isotope (. DELTA.47), a Kirgbrak group burial history curve (shown as curve B in FIG. 7) of the northwest seismic denier system of the tower was established.
Step S9: based on the burying history model after the correction of the work area, the burying depth of each stage of subcarbonate cement is obtained by using the absolute age of each stage of subcarbonate cement to cast points in the corrected burying history curve (as shown in table 2).
TABLE 2 buried depth of different structural components of kyanite of Qigrelac group in northwest seismic system of tower
Step S10: and (3) counting the porosity of the existing pore surface of each rock sample, the area percentage of the sub-carbonate cement at each stage and the corrosion pore-increasing amount of each stage based on the representative rock sample of the work area, and providing statistical data for quantitative recovery of the evolution history of the diagenetic pores.
Step S11: reconstructing a diagenetic pore evolution history curve which changes along with the depth based on the statistical result of the step S10 and by combining the burial depth of the secondary carbonate cement at each stage, wherein the result is shown in a figure 8;
determining a first point of a curve, wherein the initial porosity P1 is determined to be about 35% according to the sum of the face porosity of all pores of each rock sample and the area percentage of each stage of sub-carbonate diagenetic minerals, the time is a stage of the region suffering from surface-induced erosion, and the near-surface depth is D1 ═ 0;
determining a second point of the curve, wherein according to the cement below the mirror (7%), the porosity P2 after the cement is equal to the initial porosity (P1) - (. 35% -7%) 28% of the cementing pore-reducing amount (A1), and the depth D2 of the cement (600) plus 620 m;
determining a third point of the curve, counting about 2% of the corrosion hole increasing amount of the cement (II) according to the statistics under a mirror, wherein the porosity P3 after corrosion hole increasing is the porosity (P2) of the second point and the corrosion hole increasing amount (B1) in the period is 28% + 2% + 30%, and the depth is that D3 between the depths of the cement (II) and the cement (III) is between 620m and 1300 m;
determining the fourth point of the curve, counting about 5% of the cement (c) under the mirror, determining the porosity P4 (third point (P3)) after cementing and hole reducing, determining the cementing and hole reducing amount (A2) as 30% -5% -25%, and determining the depth as the depth D3 of the cement (c) as 1300 m;
determining the fifth point of the curve, counting about 1% of the erosion pore-increasing amount of the cement (c) under a mirror, wherein the porosity P5 after the erosion pore-increasing is 25% of the porosity (P4) + the erosion pore-increasing amount (B2) in the period of 25% and 26%, and the depth is that D5 is between 1300m and 1750m between the depths of the cement (c) and the cement (c);
determining the sixth point of the curve, counting the cement (r) by about 5% under a mirror, determining the porosity P6 after cementing and pore reducing (P5) as the fifth point, determining the cementing and pore reducing amount (A3) as 26% -5% as 21%, and determining the depth D6 as 1750m of the cement (r);
determining the seventh point of the curve, counting the corrosion pore-increasing amount of the cement (r) under a mirror to be about 1%, wherein the porosity after corrosion pore-increasing P7 is the sixth point porosity (P6) + the corrosion pore-increasing amount (B3) in the period is 21% + 1% + 22%, and the depth is that D7 is between 1750-;
determining the eighth point of the curve, counting according to the number of cement (c) under the mirror, wherein the porosity after cementing and hole reducing is P8 which is the seventh point (P7), the cementing and hole reducing amount (A4) which is 22-8% which is 14%, and the depth which is the depth D8 of cement (c) which is 4100 m;
determining the ninth point of the curve, counting according to the position under the mirror that the erosion pore-increasing amount of the cement is about 1 percent, the porosity P9 after the erosion pore-increasing is the eighth point porosity (P8) + the erosion pore-increasing amount of the stage (B4) is 14 percent and 1 percent is 15 percent, and the depth is that D9 is between 4100m and 3000m between the depths of the cement;
determining a tenth point of the curve, counting the cementing material under the mirror to reach 2%, and after cementing and hole reducing, determining a porosity P10 as a ninth point (P9) and a cementing and hole reducing amount (A5) as 15-2 as 13%, and determining a depth D9 as 3000 m;
determining the tenth point of the curve, wherein the corrosion hole increasing amount is not obvious and is about 0% according to the statistics of the cement under the mirror;
determining the twelfth point of the curve, counting about 2% of the cement according to the space under the mirror, determining the porosity P12 after cementing and reducing the pores as the tenth point (P11), determining the cementing and reducing the pores as the A6 as 13% -2% as 11%, and determining the depth D9 as 4000m of the cement.
Claims (19)
1. A deep carbonate reservoir pore evolution recovery method, wherein the method comprises the following steps:
obtaining a work area representative rock sample, wherein the characteristics of the work area representative rock sample comprise: the development of the holes of the rock sample, the filling of carbonate cement in the holes and the mutual cutting of the carbonate cement in the rock sample;
determining the stage of the carbonate cement in the rock sample, and carrying out isotope year measurement on the carbonate cement of each stage to obtain the absolute age of the carbonate cement of each stage; carrying out cluster isotope test on the carbonate cement of each stage to obtain the formation temperature of the carbonate cement of each stage;
acquiring an ancient earth temperature model and a burial history model of a work area; correcting the burial history model of the work area by using the absolute age of the subcarbonate cement of each stage and the formation temperature of the subcarbonate cement of each stage and combining the ancient geothermal temperature model of the work area to obtain a burial history model after correction of the work area;
determining the buried depth formed by each stage of the sub-carbonate cement by using the absolute age of the sub-carbonate cement of each stage based on a buried history model after work area correction;
on the basis of the representative rock sample of the work area, counting the porosity of the existing pore surface of each rock sample, the area percentage of the sub-carbonate cement of each stage and the sub-erosion pore-increasing amount of each stage; and based on the statistical result, reconstructing the evolution history of the diagenetic porosity along with the change of the burial depth by combining the burial depth formed by the secondary carbonate cement at each stage.
2. The method of claim 1, wherein the determining the age of the intercrossed carbonate cement in the rock sample is performed using a sample slice a made from a representative rock sample of the work area, the sample slice a having a thickness of 30 ± 3 μ ι η.
3. The method of claim 1, wherein the performing isotope year measurements is performed using a sample slice B made from a rock sample representative of the work area, the sample slice B having a thickness of 80-100 μ ι η.
4. The method of claim 1, wherein the performing a cluster isotope test is performed using a powder sample of each stage of a subcarbonate cement.
5. The method of claim 1, wherein the method comprises:
obtaining a work area representative rock sample, wherein the characteristics of the work area representative rock sample comprise: the development of the holes of the rock sample, the filling of carbonate cement in the holes and the mutual cutting of the carbonate cement in the rock sample;
respectively preparing at least 2 parallel samples corresponding to the representative rock samples of the obtained work area, preparing sample slices A and sample slices B of the representative rock samples by using the parallel samples, and reserving the residual parts of the parallel samples;
observing the carbonate cement of the sample slice A, and determining the period of the carbonate cement in the rock sample;
in the corresponding sample slice B, the carbonate cements of each stage corresponding to the carbonate cements of each stage in the sample slice A are defined, and the absolute age of the carbonate cements of each stage is obtained by isotope year measurement;
in the corresponding parallel sample residual part, obtaining a powder sample of carbonate cement of each stage corresponding to the carbonate cement of each stage in the sample slice A, and carrying out cluster isotope test to obtain the formation temperature of the carbonate cement of each stage;
acquiring an ancient earth temperature model and a burial history model of a work area; correcting the burial history model of the work area by using the absolute age of the subcarbonate cement of each stage and the formation temperature of the subcarbonate cement of each stage and combining the ancient geothermal temperature model of the work area to obtain a burial history model after correction of the work area;
determining the burial depth of each stage of sub-carbonate cement by using the absolute age of each stage of sub-carbonate cement based on a buried history model after work area correction;
based on the representative rock samples of each work area, counting the porosity of the existing pore surface of each rock sample, the area percentage of the sub-carbonate cement at each stage and the sub-erosion pore-increasing amount at each stage; based on the statistical result, combining the burial depth formed by the secondary carbonate cement at each stage, and reconstructing the pore evolution history along with the change of the burial depth;
preferably, the thickness of the sample sheet a is 30 ± 3 μm;
preferably, the sample slice a has a diameter of 1.5-2.5 cm;
preferably, the thickness of the sample sheet B is 80 to 100 μm;
preferably, the sample slice B has a diameter of 1.5-2.5 cm.
6. The method of claim 5, wherein the separately preparing at least 2 parallel samples corresponding to each representative rock sample is performed by: each representative rock sample was cut into cylinders with a diameter of 1.5-2.5cm and a thickness of 0.8cm, and 2 parallel samples were prepared along both sides of the cut surface.
7. The method according to claim 5, wherein the sample sheet A and the sample sheet B preferably have a mirror image similarity of not less than 90%.
8. The method of any one of claims 1-7, wherein the stage of carbonate cement in the definitive rock sample comprises:
establishing a complete and reliable diagenetic sequence according to the mutual intersection relationship of the carbonate cements, and determining the period of the carbonate cements with the mutual intersection relationship;
carbonate cements that do not have an inter-relationship, as a single stage.
9. The method of any one of claims 1-7, wherein the isotope dating is performed using a laser in situ U-Pb isotope dating modality.
10. The method as claimed in any one of claims 1 to 7, wherein the obtaining of the work area burial history model is performed as follows: and establishing an ancient earth temperature model and a burial history model of the work area according to the geological background, the drilling well and the seismic data of the area.
11. The method of any one of claims 1-7, wherein the correcting the site burial history model using the absolute age of each stage of the sub-carbonate cement and the formation temperature of each stage of the sub-carbonate cement, in conjunction with a site paleogeothermal model, comprises:
the absolute age of each stage of the secondary carbonate cement is put into a burial history model to obtain the first burial depth of each stage of the secondary carbonate cement, and the formation temperature of each stage of the secondary carbonate cement is used for calculating the second burial depth of each stage of the secondary carbonate cement according to an ancient geothermal model;
if the first burial depth of the secondary carbonate cement at each stage is inconsistent with the second burial depth, the burial history model is unreliable, and the burial history curve is modified to ensure that the burial depth of the secondary carbonate cement at each stage is positioned at the second burial depth, so that the burial history curve after the work area correction is obtained;
if the first burial depth of the subcarbonate cement at each stage is consistent with the second burial depth, the absolute age of the subcarbonate cement at each stage and the forming temperature of the subcarbonate cement at each stage are considered to form a mutual evidence-based relationship, the burial history curve is reliable, and the burial history curve is used as a burial curve model after correction of a work area.
12. The method of any one of claims 1-7, wherein reconstructing the history of diagenetic pore evolution as a function of burial depth based on statistical results in combination with the burial depth for each stage of sub-carbonate cement formation comprises:
the initial porosity for the epigenetic stage is determined based on the porosity of the existing pore face and the area percentage of the sub-carbonate cement at each stage, and the sum of the porosity of the existing pore face and the area percentage of the sub-cement at each stage is the initial porosity for the epigenetic stage.
13. The method of any one of claims 1-7, wherein reconstructing the history of diagenetic pore evolution as a function of burial depth based on statistical results in combination with the burial depth for each stage of sub-carbonate cement formation comprises:
the amount of porosity reduction for each stage of cementation was determined based on the area percentage of the sub-carbonate cement at each stage.
14. The method of any one of claims 1-7, wherein reconstructing the history of diagenetic pore evolution as a function of burial depth based on statistical results in combination with the burial depth for each stage of the secondary carbonate cement comprises: and determining the corresponding burial depth of each stage of cementation subtractive hole based on the burial depth of each stage of secondary carbonate cement, wherein the burial depth of each stage of secondary carbonate cement is the corresponding burial depth of the stage of cementation subtractive hole.
15. The method of any one of claims 1-7, wherein reconstructing the history of diagenetic pore evolution as a function of burial depth based on statistical results in combination with the burial depth for each stage of sub-carbonate cement formation comprises:
the porosity after each stage of cementation debulking is determined based on the existing porosity of the porosity face, the area percentage of each stage of sub-carbonate cement, the porosity after each stage of cementation debulking being equal to the sum of the existing porosity of the porosity face and the area percentage of each stage of sub-carbonate cement minus the total area percentage of the stage and the previous stage of sub-carbonate cement.
16. The method of any one of claims 1-7, wherein reconstructing the history of diagenetic pore evolution as a function of burial depth based on statistical results in combination with the burial depth for each stage of sub-carbonate cement formation comprises:
and determining the erosion pore-increasing amount of each stage based on the erosion pore-increasing amounts of the stages.
17. The method of any one of claims 1-7, wherein reconstructing the history of diagenetic pore evolution as a function of burial depth based on statistical results in combination with the burial depth for each stage of the secondary carbonate cement comprises:
and determining the depth corresponding to each stage of erosion increase hole based on the buried depth of the sub-carbonate cement of each stage, wherein the depth corresponding to each stage of erosion increase hole is between the buried depth of the eroded carbonate cement and the buried depth of the next stage of carbonate cement.
18. The method of any one of claims 1-7, wherein reconstructing the history of diagenetic pore evolution as a function of burial depth based on statistical results in combination with the burial depth for each stage of sub-carbonate cement formation comprises:
the porosity after each stage of erosion and pore enlargement is determined based on the porosity of the existing pore face, the area percentage of the secondary carbonate cement at each stage and the erosion and pore enlargement amount of each stage, and the porosity after each stage of cementation and pore reduction is the sum of the porosity of the existing pore face and the area percentage of the cement at each stage-the total area percentage of the secondary carbonate cement at the stage and the previous stage + the erosion and pore enlargement amount of the stage and the previous stage.
19. The method of any one of claims 1-7, wherein reconstructing the history of diagenetic pore evolution as a function of burial depth based on statistical results in combination with the burial depth for each stage of sub-carbonate cement formation comprises:
the final porosity is determined based on the existing porosity of the pore face.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011227968.XA CN112528226B (en) | 2020-11-06 | 2020-11-06 | Deep carbonate reservoir pore evolution recovery method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011227968.XA CN112528226B (en) | 2020-11-06 | 2020-11-06 | Deep carbonate reservoir pore evolution recovery method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112528226A true CN112528226A (en) | 2021-03-19 |
| CN112528226B CN112528226B (en) | 2022-10-04 |
Family
ID=74979775
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011227968.XA Active CN112528226B (en) | 2020-11-06 | 2020-11-06 | Deep carbonate reservoir pore evolution recovery method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112528226B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113049471A (en) * | 2021-03-23 | 2021-06-29 | 中国石油大学(北京) | Recovery method for porosity evolution process of carbonate rock sequence stratum |
| CN114371272A (en) * | 2022-01-11 | 2022-04-19 | 中国石油大学(北京) | Shale reservoir diagenesis numerical simulation method, device, equipment and medium |
| CN115128048A (en) * | 2021-03-25 | 2022-09-30 | 中国石油天然气股份有限公司 | Method and device for determining formation time of multi-stage carbonate cement |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102748016A (en) * | 2012-07-10 | 2012-10-24 | 中国石油大学(华东) | Geologic history period sandstone reservoir porosity evolution recovery method |
| CN102778421A (en) * | 2012-07-10 | 2012-11-14 | 中国石油大学(华东) | Permeability evolution recovery method for sandstone reservoir in geological history period |
| CN109324345A (en) * | 2018-12-03 | 2019-02-12 | 中国石油化工股份有限公司 | Superimposed Basins oil and gas entrapment timing rock porosity restoration methods |
-
2020
- 2020-11-06 CN CN202011227968.XA patent/CN112528226B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102748016A (en) * | 2012-07-10 | 2012-10-24 | 中国石油大学(华东) | Geologic history period sandstone reservoir porosity evolution recovery method |
| CN102778421A (en) * | 2012-07-10 | 2012-11-14 | 中国石油大学(华东) | Permeability evolution recovery method for sandstone reservoir in geological history period |
| CN109324345A (en) * | 2018-12-03 | 2019-02-12 | 中国石油化工股份有限公司 | Superimposed Basins oil and gas entrapment timing rock porosity restoration methods |
Non-Patent Citations (2)
| Title |
|---|
| 潘荣等: "库车坳陷克拉苏冲断带深部巴什基奇克组致密储层孔隙演化定量研究", 《地学前缘》 * |
| 陆江等: "孔隙度反演回剥法在储层物性定量预测中的应用:以珠Ⅲ坳陷文昌区为例", 《地质科技情报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113049471A (en) * | 2021-03-23 | 2021-06-29 | 中国石油大学(北京) | Recovery method for porosity evolution process of carbonate rock sequence stratum |
| CN113049471B (en) * | 2021-03-23 | 2021-10-08 | 中国石油大学(北京) | Recovery method for porosity evolution process of carbonate rock sequence stratum |
| US11487045B2 (en) | 2021-03-23 | 2022-11-01 | China University Of Petroleum-Beijing | Method for recovering porosity evolution process of sequence stratigraphy of carbonate rocks |
| CN115128048A (en) * | 2021-03-25 | 2022-09-30 | 中国石油天然气股份有限公司 | Method and device for determining formation time of multi-stage carbonate cement |
| CN115128048B (en) * | 2021-03-25 | 2024-05-28 | 中国石油天然气股份有限公司 | Method and device for determining formation time of multi-stage carbonate cement |
| CN114371272A (en) * | 2022-01-11 | 2022-04-19 | 中国石油大学(北京) | Shale reservoir diagenesis numerical simulation method, device, equipment and medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112528226B (en) | 2022-10-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112528226B (en) | Deep carbonate reservoir pore evolution recovery method | |
| Houben et al. | Pore morphology and distribution in the Shaly facies of Opalinus Clay (Mont Terri, Switzerland): Insights from representative 2D BIB–SEM investigations on mm to nm scale | |
| US9121808B2 (en) | Procedure for the determination of effective and total porosity of carbonated sedimentary rocks, and morphology characterization of their micro and nanopores | |
| Zhang et al. | Major factors controlling fracture development in the Middle Permian Lucaogou Formation tight oil reservoir, Junggar Basin, NW China | |
| Liu et al. | Impact of thermal maturity on the diagenesis and porosity of lacustrine oil-prone shales: Insights from natural shale samples with thermal maturation in the oil generation window | |
| CN112528463B (en) | Old carbonate reservoir diagenesis-pore evolution recovery method | |
| CN112305637B (en) | Reconstruction method of buried history of old sea carbonate rock | |
| Yuan et al. | Genetic mechanism of high-quality reservoirs in Permian tight fan delta conglomerates at the northwestern margin of the Junggar Basin, northwestern China | |
| Chen et al. | Characterization of the paleocrusts of weathered Carboniferous volcanics from the Junggar Basin, western China: Significance as gas reservoirs | |
| Roverato et al. | Characterization of the Cubilche debris avalanche deposit, a controversial case from the northern Andes, Ecuador | |
| Hobbs et al. | Engineering geology of British rocks and soils: Mudstones of the Mercia Mudstone Group | |
| Zhang et al. | Hydrogeomorphologic architecture of epikarst reservoirs in the middle-lower ordovician, Tazhong uplift, Tarim Basin, China | |
| Fattahpour et al. | Effect of grain characteristics and cement content on the unconfined compressive strength of artificial sandstones | |
| Nader et al. | Diagenetic effects of compaction on reservoir properties: the case of early callovian “Dalle Nacrée” formation (Paris basin, France) | |
| CN111610126A (en) | Method and system for identifying and evaluating anti-compaction effect of dolomite reservoir pores | |
| CN112946780B (en) | Method and device for determining running and sliding fracture activity period | |
| CN112881385B (en) | Carbonate oil-gas reservoir geological process reconstruction method based on dating technology | |
| CN108663719A (en) | The method and system that reservoir historical relics restore | |
| CN109187161A (en) | A kind of quantitative evaluation method of clastic rock andesine dissolution extent | |
| Guerra et al. | A comparative study of the effects of clay content on the fracture conductivity of the Eagle Ford Shale and Marcellus Shale Formations | |
| Martins et al. | Interaction between geological and geotechnical investigations of a sandstone residual soil | |
| CN112304839B (en) | Timed and quantitative carbonate rock formation-pore evolution recovery method | |
| Wang et al. | Soil–water retention behavior of a loess–paleosol sequence and its significance for hydrology and paleoclimate: a case study from the Luochuan profile of the Loess Plateau, China | |
| Shi et al. | Difference of microfeatures among diagenetic facies in tight sandstone reservoirs of the Triassic Yanchang Formation in the Midwestern region, Ordos Basin | |
| CN112859199B (en) | Method for recovering evolution history of carbonate rock diagenetic environment |
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 |