CN112924648A - Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity - Google Patents
Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity Download PDFInfo
- Publication number
- CN112924648A CN112924648A CN202110101034.XA CN202110101034A CN112924648A CN 112924648 A CN112924648 A CN 112924648A CN 202110101034 A CN202110101034 A CN 202110101034A CN 112924648 A CN112924648 A CN 112924648A
- Authority
- CN
- China
- Prior art keywords
- core
- reaction
- mineralization
- rock
- mineral
- 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
Images
Classifications
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N1/00—Sampling; Preparing specimens for investigation
- G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention belongs to the technical field of geochemistry, and provides a method for evaluating geological sequestration CO2Method for mineralizing evolution law and sequestration amount, exploring spatial and temporal evolution law of true core mineralization degree from multiple angles, and quantitatively calculating CO2Mineralization of sequestration rate and amount, and establishment of CO2Model for evaluating mineralization and sequestration potential for developing CO2Geological sequestration provides technical support for the assessment of mineralization sequestration potential. The analysis method is clear and easy to realize, can synchronously carry out a plurality of water rock reaction experiments, can form batch analysis, and is convenient for industry implementation.
Description
Technical Field
The invention belongs to the technical field of geochemistry, and relates to CO2An assessment technique for mineralization and sequestration, in particular to quantitative assessment of CO in formation cores2Method for mineralizing evolution law and sealing storage quantity
Background
Since the first industrial revolution, the widespread use of fossil energy has released a great deal of energy as CO2Is a representative greenhouse gas, causing global climate crisis. CO 22Geological sequestration is an efficient, clean carbon treatment technique, particularly with CO2The mineralization is most effective and safe for sealing and can nearly permanently seal CO2And has gained wide attention from the academic world and governments of all countries. Basalt sequestration project Carbfix for iceland, CO injection started in 20122Is the first CO based on the mineralization and sequestration mechanism2And sealing the item. A related study (Science,2016,352,1312-2The effectiveness of mineralization and storage has the advantages of early mineralization and storage starting time, large storage amount, small storage risk and the like.
However, CO is compared to other sequestration mechanisms2Mineral sequestration lacks intensive research. Due to the experimental lack of the corresponding CO2The research method of the mineralization evolution law has long mineralization reaction time scale, so that the previous research is mostly based on actual stratum simplified models to develop CO2The reaction transport numerical simulation is carried out, but the numerical model is limited by parameter uncertainty, and the geochemical model of most stratum scales cannot be effectively verified, so that the reliability of the result is not high. For CO2The mineralization and sequestration potential is evaluated, and the existing volume method is presumed according to a typical chemical reaction equation and does not consider the actually generated reaction process and the reaction progress degree, namely CO2The accuracy of the occlusion quantity estimation remains to be questioned. In summary, the complexity of the formation structure, lithofacies composition, and CO are considered to be truly reactive2Unsteady state characteristics of the mineralization process, difficulty in exploring the mineralization reaction process in real time and quasi-steady state characteristicsDetermination of CO2And (4) mineralizing and storing the water.
Aiming at the problem that the spatio-temporal evolution rule of mineralization degree and CO are difficult to be explored at present2The accuracy of the evaluation of the mineralization sequestration amount is not high, and provides an evaluation method for geological sequestration CO2The method for mineralizing evolution law and sealing stock can effectively evaluate the long-term CO of the target stratum2And (5) mineralizing and sealing.
Disclosure of Invention
Aiming at the current CO2The invention integrates various test characterization data, provides a method for evaluating geological sequestration CO based on thermodynamic equation of state2Method for mineralizing evolution law and sequestration amount, exploring spatial and temporal evolution law of true core mineralization degree from multiple angles, and quantitatively calculating CO2Mineralization of sequestration rate and amount, and establishment of CO2Model for evaluating mineralization and sequestration potential for developing CO2Geological sequestration provides technical support for the assessment of mineralization sequestration potential.
The technical scheme of the invention is as follows:
evaluation geological sequestration CO2The method for mineralizing evolution law and sealing storage amount comprises the following steps:
step one, processing stratum core and lithofacies slice and pretreating
(1.1) selecting a stratum characteristic rock sample, processing a rock core and representing a rock stratum with corresponding real composition and structure; preferably, the core machining specific dimensions are typically 25mm in diameter and 50mm in length.
(1.2) at the end of the core, a slice of lithofacies is cut out, with a thickness of 2-3mm, for initial characterization of the sample before reaction.
Step two, initial characterization of the core
(2.1) carrying out a scanning electron microscope and energy spectrum analysis (SEM + EDS) test on the core slice, analyzing the micro appearance and the composition elements of the rock sample mineral and preliminarily determining the mineral type;
(2.2) lightly polishing the rock core slice, removing the gold spraying layer, and then crushing and grinding; performing XRD test, and obtaining the types and proportions of rock minerals based on the quantitative analyte phases of the constituent elements tested by SEM and EDS; performing a total carbon analysis (TC) test to obtain the initial carbon content of the rock sample as a basic reference value; preferably, in order to reduce the influence of particle size and the preferential orientation of X-ray diffraction, the screen should be 300-400 meshes, and a ground glass sheet is used for sample preparation;
and (2.3) carrying out CT scanning on the core, constructing a digital core, and obtaining the pore distribution characteristics of the digital core.
Coating the outer surfaces of the core except one radial end face with epoxy resin, and then normally air-drying; preferably, the epoxy is used to coat the outer surface of the core because it reacts with CO2Water does not react and can effectively block CO2Migration and diffusion are carried out, and meanwhile, the experimental result is not influenced.
Step four, carrying out a simulated mineralization experiment and calculating CO2Mineralization sequestration rate;
and (3) carrying out a static constant-pressure mineralization reaction test by considering the characteristics of long mineralization reaction time scale, low stratum fluid flow rate, diffusion limitation on mass transfer and the like in a real stratum.
(4.1) putting the rock core into a reaction kettle, adding ultrapure water to completely submerge, closing the system, vacuumizing, and simultaneously heating the air bath to the reaction temperature;
(4.2) isothermal CO was injected from the upper end of the reactor2Maintaining constant pressure until the target constant pressure reaction time; preferably, the constant pressure reaction time is dependent on the overall experimental plan and the mineral composition of the characteristic core, and is at least 10 days (the general principle is that the core reactive mineral content is high and the minimum constant pressure experimental time can be suitably shortened).
(4.3) CO stopping2Injecting at constant pressure, sealing the reaction kettle, and recording the pressure data change for a period of time. Preferably, to efficiently calculate CO of reactive rocks2The mineralization reaction rate and the closed mineralization experiment last for at least 5 days, but the content of the reactive minerals of the core needs to be comprehensively considered (the general principle is that the content of the reactive minerals is high, and the minimum closed experiment time can be properly shortened).
(4.4) slow mineralization in CO2After the initial injection is stabilized, holdThe continuous mineralization reaction can be regarded as a constant speed reaction in a long time later, and CO can be obtained based on the pressure data of a closed experiment2The rate of mineralization reaction sequestration. Calculating CO in the closed experiment based on the following formula2Mineralization and storage amount:
wherein V is CO in the reaction kettle2Volume, unit m3(ii) a R is the gas equation of state constant, 8.314J/(mol K); t is the temperature of the reaction kettle system in K; p is the pressure of the reaction kettle, the unit is pa, the subscript I represents an initial system of a closed experiment, and the subscript t represents a real-time system at the moment of the closed experiment; z is CO2The compression factor of the gas state equation is calculated or inquired in an NIST database by an empirical formula; n is CO2The amount of substance(s) in mol.
Further, the CO per mass of the core can be calculated2The mineralization rate, the formula is as follows:
wherein k is the core CO per unit mass2The mineralization and storage rate is expressed in mol/(kg. day); m isRock (A. B. EThe dry mass of the core before reaction is unit kg; Δ n and Δ t are CO in the blocking experiment, respectively2The change in the amount of substance and the corresponding change in the reaction time in mol and day, respectively;
and (4.5) after the experimental time reaches the target time, closing the air bath, and taking out the rock core after gradient pressure relief. And drying and removing the epoxy resin layer. Preferably, drying is carried out for 48h at 60 ℃.
Step five, testing and analyzing the rock core after reaction, exploring the spatial and temporal evolution rule of mineralization reaction, and establishing CO2And (3) a mineralization sequestration potential evaluation model.
(5.1) carrying out CT scanning after reaction on the rock core, constructing the digital rock core after reaction, acquiring microscopic pore distribution of the digital rock core, comparing CT images before reaction, determining a mineral dissolution and precipitation area, and acquiring a spatial evolution law of a pore three-dimensional structure and size distribution, particularly porosity change.
(5.2) the core is cut axially into two halves, a and b.
(5.3) cutting the a into thin slices at equal intervals along the axial direction, and grinding the thin slices into powder respectively. A part of samples of each slice pass through a 300-mesh screen, an XRD (X-ray diffraction) phase quantitative analysis test is carried out, a secondary precipitated carbonate phase is determined based on the analysis of the initial mineral phase composition, and the spatial evolution rule of the precipitated carbonate phase composition and the mass fraction thereof is obtained; performing TC test on the other part of the sample powder, obtaining the spatial distribution of the C sequestration amount, and calculating the total CO of the core2And (4) mineralizing and storing the water. The detailed calculation is according to the following formula:
mCO2＝∫mCO2,i
wherein m isC,IIs the carbon content of the pre-reaction core in grams/kg for the TC test; m isC,iThe carbon content of the ith slice of the rock mass a after reaction in the TC test is g/kg; mCO2、MCIs respectively CO2The molecular and atomic mass of C in g/mol; m isCO2,i、mCO2Respectively the ith layer CO of the converted rock mass a2Mineralized occluded amount and rock mass a total CO2The amount of mineralized seal stock is g/kg. Preferably, the sheet thickness is 3-4 mm.
Further, based on the proportion of different carbonate phases analyzed by XRD data of each slice and the combination of C sealed storage of corresponding rock slices, the amount of substances of different carbonates can be calculated, and then core CO is obtained2The spatial distribution rule of carbonate phase precipitation in mineralization and storage.
And (5.4) processing the b sheets along the axial direction, and performing Map test on the central axial plane by using Raman to obtain a two-dimensional Raman spectrum of the central axial plane of the core. Based on the core mineral phase and the carbonate precipitation phase obtained by the analysis, inquiring a database to obtain a standard spectrum of the core mineral phase and the carbonate precipitation phase, obtaining a mineral content distribution map of a central axial surface through data processing, and comparing the obtained spatial distribution of the carbonate phase precipitation to obtain the spatial distribution and the content of the original mineral corrosion and the carbonate precipitation of the axial surface.
(5.5) combining with the dissolution and precipitation characteristic region in the gray-scale image of the central axial plane of the CT digital core after the corresponding reaction, the front depth of dissolution and precipitation of the minerals of the digital core can be obtained, particularly whether carbonate exists in the deep part of the core, and whether the pores are blocked by the dissolution and precipitation of the minerals is researched.
(5.6) observing the mineral dissolution, carbonate precipitation microscopic morphology and mineral distribution of the feature region of the central axial plane of the slice along the reaction depth direction by adopting SEM + EDS based on the dissolution and precipitation feature region in the central axial plane gray scale image of the CT digital core, and obtaining the local mineral dissolution and precipitation space-time law.
(5.7) CO per unit mass of core calculated according to step four2Rate of mineralization, core CO was estimated according to the following formula2And (3) sealing potential:
N＝k·T·a
wherein N is CO per unit mass of the core2Mineralized and sealed amount, unit is mol/kg, k is core CO of unit mass2The mineralization and storage rate is expressed in mol/(kg. day); t is evaluation time in day; a is a correction coefficient, mainly comprehensively considers the mineral proportion, the pore structure and the like of the core, and can calculate the core CO based on the TC test data of the core after reaction2Mineralized sealing amount mCO2And the linear fitting is carried out by utilizing the formula. Preferably, the fit for a particular formation is best obtained from a plurality of sets of experimental data of mineralization reactions at different times.
Preferably, for the analysis of a certain characteristic stratum in practical application, mineralization reaction experiments of different times should be carried out based on the method so as to find out geological sequestration CO of the stratum2Mineralisation evolution law and establishment of CO2And (3) a mineralization sequestration potential evaluation model.
The invention has the beneficial effects that:
(1) the method of the invention is based on characteristic rock core and adoptsMultiple testing techniques, multi-angle analysis of CO2Mineralization processes, including lithofacies change, CO2The mineralization rate and the like can be used for exploring the space-time evolution law of the dissolution and the precipitation of minerals and CO of the whole and local parts of the characteristic rock core2Spatial distribution of mineralized occluded amount to build up CO2The mineralization and sequestration potential evaluation model is finally helpful for evaluating CO of a real target stratum2Mineralization and sequestration potential. At the same time, for subsequent CO2And (4) carrying out mineralization reaction transport simulation, and providing a standard for verifying the correctness of basic parameters and models.
(2) The analysis method is clear and easy to realize, can synchronously develop a plurality of water rock reaction experiments (mineralization experiments with different time, temperature and pressure are developed according to actual evaluation requirements), can form batch analysis, and is convenient for industry implementation.
Drawings
FIG. 1 is a flow chart of the experiment and analysis of the method of example 1 of the present invention.
Detailed Description
The technical solutions of the present invention will be described in detail below for the purpose of clearly understanding the technical features, objects and advantages of the present invention, but the present invention is not limited thereto.
Example 1
In view of the above problems, the present invention provides a method for evaluating geological sequestration of CO2The specific flow of the method for mineralizing evolution law and sealing storage is shown in figure 1. The method comprises the following steps:
step one, processing stratum rock cores and rock phase slices and preprocessing.
And (1.1) selecting a stratum characteristic rock sample aiming at a target reactive stratum, and processing a rock core to represent a rock stratum with corresponding real composition and structure. Preferably, the core machining specific dimensions are typically 25mm in diameter and 50mm in length.
(1.2) cutting out a petrographic thin slice with the thickness of 2-3mm at the end of the core by using a miniature core cutting machine, and using the petrographic thin slice for initial characterization of a sample before reaction.
And step two, initial characterization of the core.
(2.1) carrying out a scanning electron microscope and energy spectrum analysis (SEM + EDS) test on the core slice, analyzing the microscopic morphology and the constituent elements of the rock sample mineral and preliminarily determining the mineral type.
And (2.2) slightly grinding and polishing the thin sheet by using a standard rock sample polishing machine, removing the gold spraying layer, and then crushing and grinding by using an agate lapping body. One part is subjected to XRD test, and based on the quantitative analysis phase of the constituent elements tested by SEM and EDS, the types and proportions of rock minerals are obtained. In another part, a total carbon analysis (TC) test is performed to obtain the initial carbon content of the rock sample, typically organic matter or native carbonate, as a basic reference value. Preferably, to reduce the particle size effect and the preferred orientation of X-ray diffraction, the screen should be 300-400 mesh, and the sample is prepared using a ground glass plate.
And (2.3) carrying out CT scanning on the core, constructing a digital core, and obtaining the pore distribution characteristics of the digital core.
And step three, coating the rest outer surfaces of the core except one radial end face with epoxy resin, and then normally air-drying. Preferably, the epoxy is used to coat the outer surface of the core because it reacts with CO2Water does not react and can effectively block CO2Migration and diffusion are carried out, and meanwhile, the experimental result is not influenced.
Step four, carrying out a simulated mineralization experiment and calculating CO2Mineralization sequestration rate.
And (3) carrying out a static constant-pressure mineralization reaction test by considering the characteristics of long mineralization reaction time scale, low stratum fluid flow rate, diffusion limitation on mass transfer and the like in a real stratum.
(4.1) putting the rock core into a reaction kettle, adding ultrapure water to completely submerge, closing the system, vacuumizing for 24h, and simultaneously heating the air bath to the reaction temperature. Preferably, the diameter of the reaction kettle is selected according to the actual water-rock ratio planned by the experiment, and is generally selected to be 30-40mm, which is closer to the real formation water-rock ratio.
(4.2) injecting isothermal CO from the upper end of the reaction kettle by adopting a plunger pump2And (4) reaching the target pressure, and maintaining the constant pressure until the target constant pressure reaction time. Preferably, the constant pressure reaction time is dependent on the overall experimental plan and the mineral composition of the characteristic core for a minimum of not less than 10 days (The general principle is that the core reactive mineral content is high, and the minimum constant pressure experiment time can be properly shortened).
(4.3) CO stopping2Injecting at constant pressure, sealing the reaction kettle, and recording the pressure data change for a period of time. Preferably, to efficiently calculate CO of reactive rocks2The mineralization reaction rate and the closed mineralization experiment last for at least 5 days, but the content of the reactive minerals of the core needs to be comprehensively considered (the general principle is that the content of the reactive minerals is high, and the minimum closed experiment time can be properly shortened).
(4.4) slow mineralization in CO2After the initial injection is stable, the continuous mineralization reaction can be regarded as a constant-speed reaction for a long time later, and CO can be obtained based on the pressure data of a closed experiment2The rate of mineralization reaction sequestration. Calculating CO in the closed experiment based on the following formula2Mineralization and storage amount:
wherein V is CO in the reaction kettle2Volume, unit m3(ii) a R is the gas equation of state constant, 8.314J/(mol K); t is the temperature of the reaction kettle system in K; t is reaction time, in units of s; p is the pressure of the reaction kettle, the unit is pa, the subscript I represents an initial system of a closed experiment, and the subscript t represents a real-time system at the moment of the closed experiment; z is CO2The compression factor of the gas state equation can be calculated by an empirical formula or inquired into a NIST database; n is CO2The amount of substance(s) in mol.
Further, the CO per mass of the core can be calculated2The mineralization rate, the formula is as follows:
wherein k is the core CO per unit mass2The mineralization and storage rate is expressed in mol/(kg. day); m isRock (A. B. EThe dry mass of the core before reaction is unit kg; Δ n and Δ t are respectivelyIs CO in a closed experiment2The change in the amount of substance corresponds to the change in the reaction time in mol and day.
And (4.5) after the experimental time reaches the target time, closing the air bath, and taking out the rock core after gradient pressure relief. And drying and removing the epoxy resin layer. Preferably, drying is carried out for 48h at 60 ℃.
Step five, testing and analyzing the rock core after reaction, exploring the spatial and temporal evolution rule of mineralization reaction, and establishing CO2And (3) a mineralization sequestration potential evaluation model.
(5.1) carrying out CT scanning after reaction on the rock core, constructing the digital rock core after reaction, obtaining the micro-pore distribution of the digital rock core, comparing CT images before reaction, and determining CO2And (3) mineralizing the mineral dissolution and precipitation area, and acquiring the space evolution law of the three-dimensional structure and size distribution of the pores, particularly the change of the porosity.
(5.2) the core is cut axially into two halves, a and b.
(5.3) cutting the a into thin slices at equal intervals along the axial direction, and grinding the thin slices into powder respectively. A part of samples of each slice pass through a 300-mesh screen, an XRD (X-ray diffraction) phase quantitative analysis test is carried out, a secondary precipitated carbonate phase is determined based on the analysis of the initial mineral phase composition, and the spatial evolution rule of the precipitated carbonate phase composition and the mass fraction thereof is obtained; performing TC test on the other part of the sample powder, obtaining the spatial distribution of the C sequestration amount, and calculating the total CO of the core2And (4) mineralizing and storing the water. The detailed calculation is according to the following formula:
mCO2＝∫mCO2,i
wherein m isC,IIs the carbon content of the pre-reaction core in grams/kg for the TC test; m isC,iThe carbon content of the ith slice of the rock mass a after reaction in the TC test is g/kg; mCO2、MCIs respectively CO2The molecular and atomic mass of C in g/mol; m isCO2,i、mCO2Respectively the ith layer C of the converted rock mass aO2Mineralized occluded amount and rock mass a total CO2The amount of mineralized seal stock is g/kg. Preferably, the sheet thickness is 3-4 mm.
Further, based on the proportion of different carbonate phases analyzed by XRD data of each slice and the combination of C sealed storage of corresponding rock slices, the amount of substances of different carbonates can be calculated, and then core CO is obtained2The spatial distribution rule of carbonate phase precipitation in mineralization and storage.
And (5.4) processing the b sheets along the axial direction, polishing the central axial plane, and performing Map test on the central axial plane by using Raman to obtain a two-dimensional Raman spectrum of the central axial plane of the core. Based on the core mineral phase and the carbonate precipitation phase obtained by the analysis, inquiring a database to obtain a standard spectrum of the core mineral phase and the carbonate precipitation phase, obtaining a mineral content distribution map of a central axial surface through data processing, and comparing the obtained spatial distribution of the carbonate phase precipitation to obtain the spatial distribution and the content of the original mineral corrosion and the carbonate precipitation of the axial surface.
(5.5) combining with the dissolution and precipitation characteristic region in the gray-scale image of the central axial plane of the CT digital core after the corresponding reaction, the front depth of dissolution and precipitation of the minerals of the digital core can be obtained, particularly whether carbonate exists in the deep part of the core, and whether the pores are blocked by the dissolution and precipitation of the minerals is researched.
(5.6) carrying out gold spraying on the central axial surface of the lithofacies slice, observing the mineral dissolution, the carbonate precipitation micro-morphology and the mineral distribution of the central axial surface characteristic area of the slice along the reaction depth direction by adopting SEM + EDS based on the dissolution and precipitation characteristic area in the central axial surface gray scale image of the CT digital core, and obtaining the local space-time law of mineral dissolution and precipitation.
(5.7) CO per unit mass of core calculated according to step four2Rate of mineralization, core CO was estimated according to the following formula2And (3) sealing potential:
N＝k·T·a
wherein N is CO per unit mass of the core2Mineralized and sealed amount, unit is mol/kg, k is core CO of unit mass2The mineralization and storage rate is expressed in mol/(kg. day); t is evaluation time in day; a is a correction coefficient, mainlyThe core CO can be calculated based on the TC test data of the reacted core by comprehensively considering the mineral proportion, the pore structure and the like of the core2Mineralized sealing amount mCO2And the linear fitting is carried out by utilizing the formula. Preferably, the fit for a particular formation is best obtained from a plurality of sets of experimental data of mineralization reactions at different times.
Preferably, the CO is aimed at the characteristic stratum in practical application2Mineralization and sequestration analysis, based on the method, mineralization reaction experiments should be carried out at different times so as to find out geological sequestration CO of the stratum2Mineralisation evolution law and establishment of CO2And (3) a mineralization sequestration potential evaluation model.
The above-described embodiment is only one embodiment of the present invention, and it will be apparent to those skilled in the art that various modifications or variations can be made based on the principle disclosed in the present invention, and not limited to the description of the above-described specific embodiment of the present invention, so that the above-described embodiment is only for describing the method in detail, but not in a limiting sense.
Claims (1)
1. Evaluation geological sequestration CO2The method for mineralizing evolution law and sealing storage amount is characterized by comprising the following steps:
step one, processing stratum core and lithofacies slice and pretreating
(1.1) selecting a stratum characteristic rock sample, processing a rock core and representing a rock stratum with corresponding real composition and structure;
(1.2) cutting out a lithofacies slice with the thickness of 2-3mm at the end part of the core for initial characterization of a sample before reaction;
step two, initial characterization of the core
(2.1) carrying out a scanning electron microscope and energy spectrum analysis test on the core slice, analyzing the micro appearance and the composition elements of the rock sample mineral and preliminarily determining the mineral type;
(2.2) lightly polishing the rock core slice, removing the gold spraying layer, and then crushing and grinding into powder; performing XRD test on a part of powder, and obtaining the types and proportions of rock minerals based on a scanning electron microscope and a quantitative analysis phase of constituent elements tested by energy spectrum analysis; the other part of the powder is subjected to total carbon analysis and test to obtain the initial carbon content of the rock sample as a basic reference value;
(2.3) carrying out CT scanning on the rock core, constructing a digital rock core, and obtaining pore distribution characteristics of the digital rock core;
coating the outer surfaces of the core except one radial end face with epoxy resin, and then normally air-drying;
step four, carrying out a simulated mineralization experiment and calculating CO2Mineralization sequestration rate;
(4.1) putting the rock core into a reaction kettle, adding ultrapure water to completely submerge, closing the system, vacuumizing, and simultaneously heating the air bath to the reaction temperature;
(4.2) isothermal CO was injected from the upper end of the reactor2Maintaining constant pressure until the target constant pressure reaction time;
(4.3) CO stopping2Injecting at constant pressure, sealing the reaction kettle, and recording the pressure data change of the reaction kettle for a period of time;
(4.4) slow mineralization in CO2After the initial injection is stable, the continuous mineralization reaction is regarded as a constant-speed reaction for a long time later, and based on the pressure data of the closed experiment, CO is obtained2The mineralization reaction sealing rate; calculating CO in the closed experiment based on the following formula2Mineralization and storage amount:
wherein V is CO in the reaction kettle2Volume, unit m3(ii) a R is the gas equation of state constant, 8.314J/(mol K); t is the temperature of the reaction kettle system in K; p is the pressure of the reaction kettle, the unit is pa, the subscript I represents an initial system of a closed experiment, and the subscript t represents a real-time system at the moment of the closed experiment; z is CO2The compression factor of the gas state equation is calculated or inquired in an NIST database by an empirical formula; n is CO2The amount of substance(s) in mol;
further, calculating the core mass per unitCO of2The mineralization rate, the formula is as follows:
wherein k is the core CO per unit mass2The mineralization and storage rate is expressed in mol/(kg. day); m isRock (A. B. EThe dry mass of the core before reaction is unit kg; Δ n and Δ t are CO in the blocking experiment, respectively2The change in the amount of substance and the corresponding change in the reaction time in mol and day, respectively;
(4.5) after the experimental time reaches the target time, closing the air bath, and taking out the rock core after gradient pressure relief; then drying and removing the epoxy resin layer;
step five, testing and analyzing the rock core after reaction, exploring the spatial and temporal evolution rule of mineralization reaction, and establishing CO2A mineralization sequestration potential assessment model;
(5.1) carrying out CT scanning after reaction on the rock core, constructing a digital rock core after reaction, acquiring microscopic pore distribution of the digital rock core, comparing CT images before reaction, determining mineral dissolution and precipitation influence areas, and acquiring a spatial evolution law of a pore three-dimensional structure and size distribution, particularly porosity change;
(5.2) cutting the core into two equal parts, a and b, along the axial direction;
(5.3) cutting the a into slices at equal intervals along the axial direction, and grinding the slices into powder respectively; a part of samples of each slice pass through a 300-mesh screen, an XRD (X-ray diffraction) phase quantitative analysis test is carried out, a secondary precipitated carbonate phase is determined based on the analysis of the initial mineral phase composition, and the spatial evolution rule of the precipitated carbonate phase composition and the mass fraction thereof is obtained; the other part of the sample is subjected to total carbon analysis and test to obtain the spatial distribution of the C sequestration, and the total CO of the core is calculated2Mineralizing the sealing stock; calculated according to the following formula:
mCO2＝∫mCO2,i
wherein m isC,IIs the carbon content of the pre-reaction core as measured by total carbon analysis in g/kg; m isC,iThe carbon content of the ith slice of the rock mass a after reaction is measured by total carbon analysis, and the unit is g/kg; mCO2、MCAre respectively CO2The molecular and atomic mass of C in g/mol; m isCO2,i、mCO2Respectively the ith layer CO of the converted rock mass a2Mineralized occluded amount and rock mass a total CO2Mineralizing and storing amount, wherein the unit is g/kg;
(5.4) processing the b sheets along the axial direction, and performing Map test on the central axial plane by using Raman to obtain a two-dimensional Raman spectrum of the central axial plane of the core; based on the core mineral phase and the carbonate precipitation phase obtained by the analysis, inquiring a database to obtain a standard spectrum of the core mineral phase and the carbonate precipitation phase, obtaining a mineral content distribution map of a central axial surface through data processing, and comparing the obtained spatial distribution of the carbonate phase precipitation to obtain the spatial distribution and content of original mineral corrosion and carbonate precipitation of the axial surface;
(5.5) combining with a dissolution and precipitation characteristic region in a gray scale image of a central axial plane of the CT digital core after the corresponding reaction, acquiring the mineral dissolution and precipitation front depth of the digital core, particularly whether carbonate exists in the deep part of the core, and researching whether the mineral dissolution and precipitation blocks a pore throat;
(5.6) observing the mineral dissolution, carbonate precipitation micro-morphology and mineral distribution of the feature region of the central axial plane of the slice along the reaction depth direction by adopting SEM + EDS based on the dissolution and precipitation feature region in the central axial plane gray scale image of the CT digital core, and obtaining the local mineral dissolution and precipitation space-time law;
(5.7) CO per unit mass of core calculated according to step four2Rate of mineralization, core CO was estimated according to the following formula2And (3) sealing potential:
N＝k·T·a
wherein N is CO per unit mass of the core2Mineralized and sealed amount, unit is mol/kg, k is core CO of unit mass2The mineralization and storage rate is expressed in mol/(kg. day); t is evaluation time in day;a is a correction coefficient, mineral proportion and pore structure of the core are comprehensively considered, and core CO is calculated based on TC test data of the core after reaction2Mineralized sealing amount mCO2And the linear fitting is carried out by utilizing the formula.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110101034.XA CN112924648B (en) | 2021-01-26 | 2021-01-26 | Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110101034.XA CN112924648B (en) | 2021-01-26 | 2021-01-26 | Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112924648A true CN112924648A (en) | 2021-06-08 |
| CN112924648B CN112924648B (en) | 2022-01-04 |
Family
ID=76166045
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110101034.XA Active CN112924648B (en) | 2021-01-26 | 2021-01-26 | Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112924648B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113740196A (en) * | 2021-09-02 | 2021-12-03 | 太原理工大学 | Device and method for measuring carbon sequestration amount of in-situ pyrolysis organic mineral layer |
| CN114544463A (en) * | 2022-03-24 | 2022-05-27 | 中国石油大学(北京) | Fractured reservoir geologic body CO2Buried storage evaluation method and device |
| CN117233253A (en) * | 2023-09-19 | 2023-12-15 | 中国矿业大学 | CO (carbon monoxide) 2 Mineralization enhancement and monitoring evaluation simulation device and method |
| CN117405718A (en) * | 2023-12-12 | 2024-01-16 | 山东石油化工学院 | Quantitative evaluation method for carbon sequestration capacity of magma calcium element based on XRF scanning |
| CN117665262A (en) * | 2024-02-02 | 2024-03-08 | 中国科学院地质与地球物理研究所 | Basalt CO 2 Quantitative assessment method, system and equipment for carbon sequestration potential |
| CN118348121A (en) * | 2024-06-14 | 2024-07-16 | 中国矿业大学 | Imaging and carbon fixation effect evaluation method and system in carbon dioxide mineralization process |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103544361A (en) * | 2013-11-04 | 2014-01-29 | 西北大学 | Evaluation method of CO2 geological sequestration potential in oil-gas field development |
| US20140205520A1 (en) * | 2011-06-17 | 2014-07-24 | Kjeoy Research & Education Center | Leaching of minerals and sequestration of co2 |
| CN104850742A (en) * | 2015-05-08 | 2015-08-19 | 河北工程大学 | Method for calculating CO2 salty water layer mineral sequestration potential |
| CN106383133A (en) * | 2016-10-21 | 2017-02-08 | 大连理工大学 | Method for measuring capillary pressure of CO2-saline-core system based on X-ray computed tomography (CT) technology |
| CN108614076A (en) * | 2018-04-27 | 2018-10-02 | 中国石油大学(华东) | A kind of carbon dioxide geologic sequestration evaluation method |
| CN111537018A (en) * | 2019-02-06 | 2020-08-14 | 印度电子信息技术部 | Estimating sequestered CO2Automatic carbon sequestration estimation system and method for quantities of |
-
2021
- 2021-01-26 CN CN202110101034.XA patent/CN112924648B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140205520A1 (en) * | 2011-06-17 | 2014-07-24 | Kjeoy Research & Education Center | Leaching of minerals and sequestration of co2 |
| CN103544361A (en) * | 2013-11-04 | 2014-01-29 | 西北大学 | Evaluation method of CO2 geological sequestration potential in oil-gas field development |
| CN104850742A (en) * | 2015-05-08 | 2015-08-19 | 河北工程大学 | Method for calculating CO2 salty water layer mineral sequestration potential |
| CN106383133A (en) * | 2016-10-21 | 2017-02-08 | 大连理工大学 | Method for measuring capillary pressure of CO2-saline-core system based on X-ray computed tomography (CT) technology |
| CN108614076A (en) * | 2018-04-27 | 2018-10-02 | 中国石油大学(华东) | A kind of carbon dioxide geologic sequestration evaluation method |
| CN111537018A (en) * | 2019-02-06 | 2020-08-14 | 印度电子信息技术部 | Estimating sequestered CO2Automatic carbon sequestration estimation system and method for quantities of |
Non-Patent Citations (6)
| Title |
|---|
| LANLAN JIANG ET AL: "Behavior of CO2/water flow in porous media for CO2 geological storage", 《MAGNETIC RESONANCE IMAGING》 * |
| YONGCHEN SONG ET AL.: "Experimental study on the CO2-decane displacement front behavior in", 《ENERGY》 * |
| YONGCHEN SONG ET AL: "Pore-scale visualization study on CO2 displacement of brine in micromodels", 《INTERNATIONAL JOURNAL OF GREENHOUSE GAS CONTROL》 * |
| 刘树阳: "CO2提高天然气采收率的重力效应与", 《中国石油大学学报( 自然科学版)》 * |
| 张 倩: "CO2 地质封存过程中传热特性研究", 《大连理工大学学报》 * |
| 张毅 等: "地质封存中CO2盐水溶液的表观摩尔体积模型", 《环境科学与技术》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113740196A (en) * | 2021-09-02 | 2021-12-03 | 太原理工大学 | Device and method for measuring carbon sequestration amount of in-situ pyrolysis organic mineral layer |
| CN113740196B (en) * | 2021-09-02 | 2024-04-16 | 太原理工大学 | Device and method for measuring carbon sequestration amount of in-situ pyrolysis organic ore layer |
| CN114544463A (en) * | 2022-03-24 | 2022-05-27 | 中国石油大学(北京) | Fractured reservoir geologic body CO2Buried storage evaluation method and device |
| CN114544463B (en) * | 2022-03-24 | 2023-02-03 | 中国石油大学(北京) | Fractured reservoir geologic body CO 2 Method and device for estimating burial |
| CN117233253A (en) * | 2023-09-19 | 2023-12-15 | 中国矿业大学 | CO (carbon monoxide) 2 Mineralization enhancement and monitoring evaluation simulation device and method |
| CN117233253B (en) * | 2023-09-19 | 2024-04-05 | 中国矿业大学 | CO (carbon monoxide) 2 Mineralization enhancement and monitoring evaluation simulation device and method |
| CN117405718A (en) * | 2023-12-12 | 2024-01-16 | 山东石油化工学院 | Quantitative evaluation method for carbon sequestration capacity of magma calcium element based on XRF scanning |
| CN117405718B (en) * | 2023-12-12 | 2024-02-13 | 山东石油化工学院 | Quantitative evaluation method for carbon sequestration capacity of magma calcium element based on XRF scanning |
| CN117665262A (en) * | 2024-02-02 | 2024-03-08 | 中国科学院地质与地球物理研究所 | Basalt CO 2 Quantitative assessment method, system and equipment for carbon sequestration potential |
| CN118348121A (en) * | 2024-06-14 | 2024-07-16 | 中国矿业大学 | Imaging and carbon fixation effect evaluation method and system in carbon dioxide mineralization process |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112924648B (en) | 2022-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112924648B (en) | Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity | |
| Meckler et al. | Glacial–interglacial temperature change in the tropical West Pacific: A comparison of stalagmite-based paleo-thermometers | |
| Rankin et al. | A year‐long record of size‐segregated aerosol composition at Halley, Antarctica | |
| Aeschbach-Hertig et al. | The physical structure and dynamics of a deep, meromictic crater lake (Lac Pavin, France) | |
| CN108614076A (en) | A kind of carbon dioxide geologic sequestration evaluation method | |
| Mix et al. | A hot and high Eocene Sierra Nevada | |
| Dingwell et al. | Eruptive shearing of tube pumice: pure and simple | |
| Osborne et al. | Calcification of the planktonic foraminifera Globigerina bulloides and carbonate ion concentration: Results from the Santa Barbara Basin | |
| Brenner et al. | Influence of water on clumped-isotope bond reordering kinetics in calcite | |
| Stern et al. | Processes controlling the oxygen isotope ratio of soil CO2: Analytic and numerical modeling | |
| Jayasekara et al. | Effect of salinity on supercritical CO2 permeability of caprock in deep saline aquifers: An experimental study | |
| Bénézeth et al. | Mineral solubility and aqueous speciation under hydrothermal conditions to 300 C–the carbonate system as an example | |
| Ma et al. | Role of CO2-water-rock interactions and implications for CO2 sequestration in Eocene deeply buried sandstones in the Bonan Sag, eastern Bohai Bay Basin, China | |
| Patterson et al. | Atmospheric history of H2 over the past century reconstructed from South Pole firn air | |
| Busch et al. | Predicting effective diffusion coefficients in mudrocks using a fractal model and small‐angle neutron scattering measurements | |
| Misyura et al. | Improving the efficiency of storage of natural and artificial methane hydrates | |
| CN112485282B (en) | Measuring system and method for soil-water characteristic curve of gas hydrate-containing sediment | |
| Watson et al. | Nitrogen diffusion in silicate minerals, with implications for nitrogen transport and cycling in the lithosphere | |
| CN109030775A (en) | The analytical equipment and analysis method that water flow loses solvable karst under a kind of closed system | |
| Erickson III et al. | The seasonal cycle of atmospheric CO2: A study based on the NCAR Community Climate Model (CCM2) | |
| Liu et al. | Improving the Estimation of Salt Distribution during Evaporation in Saline Soil by HP1 Model | |
| Honda et al. | Utility of an automatic water sampler to observe seasonal variability in nutrients and DIC in the Northwestern North Pacific | |
| Ding et al. | Exploring the role of the structural heterogeneity of fractured carbonate reservoirs in contact with dissolved CO2 based on fracture-water-rock simulation experiments | |
| Randi et al. | Experimental and numerical simulation of the injection of a CO2 saturated solution in a carbonate reservoir: application to the CO2-DISSOLVED concept combining CO2 geological storage and geothermal heat recovery | |
| Anovitz et al. | A quantitative approach to the analysis of reactive mineralogy and surface area |
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 |