CN108614076A - A kind of carbon dioxide geologic sequestration evaluation method - Google Patents

A kind of carbon dioxide geologic sequestration evaluation method Download PDF

Info

Publication number
CN108614076A
CN108614076A CN201810393143.1A CN201810393143A CN108614076A CN 108614076 A CN108614076 A CN 108614076A CN 201810393143 A CN201810393143 A CN 201810393143A CN 108614076 A CN108614076 A CN 108614076A
Authority
CN
China
Prior art keywords
core sample
mineral
experiment
formula
reaction
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.)
Pending
Application number
CN201810393143.1A
Other languages
Chinese (zh)
Inventor
远光辉
操应长
昝念民
王艳忠
葸克来
杨田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201810393143.1A priority Critical patent/CN108614076A/en
Publication of CN108614076A publication Critical patent/CN108614076A/en
Pending legal-status Critical Current

Links

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
    • 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

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Food Science & Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Medicinal Chemistry (AREA)
  • Pathology (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Gas Separation By Absorption (AREA)

Abstract

The present invention relates to a kind of carbon dioxide geologic sequestration evaluation method, step is：S1, core sample is selected, determines core sample feature；S2, sample pretreatment：It cuts core sample and surface polishing is carried out to the core sample after cutting；S3, indoor water rock simulated experiment is carried out to core sample；S4, the core sample after experiment is handled, the solution after removal experiment reaction；The core sample that S5, analysis test experiments post-process；S6, carry out long-term Numerical Experiment.The present invention determines CO by short-term indoor water rock simulated experiment2The response path and reaction process occurred after filling oil-gas Layer carries out the amount of sealing up for safekeeping evaluation by long-term numerical simulation, makes CO on this basis2It fills simulation process and fills experiment closer to actual carbon dioxide scale, avoid and CO is speculated according to typical chemical equation2The process of sealing up for safekeeping after filling causes the uncertainty that the amount of sealing up for safekeeping is evaluated, hence it is evident that improves the accuracy for sealing Potential Evaluation up for safekeeping for actually sealing geologic body up for safekeeping.

Description

A kind of carbon dioxide geologic sequestration evaluation method
Technical field
The invention belongs to geochemical techniques fields, are related to the assessment technique of carbon dioxide geologic sequestration, specifically, relating to And a kind of oil-gas Layer carbon dioxide geologic sequestration evaluation method.
Background technology
It is sharply increased as human industry is movable, a large amount of CO2It is discharged into the Nature, the carbon dioxide in air is dense Degree continues to increase, and greenhouse effects constantly accumulate, and so as to cause global warming, glacier and frozen soil is caused to melt, on sea level The natural ecosystems disasters such as liter, seriously threaten human survival.CO2Capture with seal up for safekeeping as control CO2 emission most have One of technical method of potentiality, a large amount of CO both domestic and external2Geological storage project demonstrates its feasibility technically.CO2Catch Obtaining and sealing up for safekeeping is by carbon dioxide capture and to seal the long-term (example of progress in oil gas field, coal gas layer, the medium geologic body of brine layer up for safekeeping Such as：Thousands of years) it seals up for safekeeping, to prevent or substantially reduce greenhouse gas emission, to mitigate the influence to terrestrial climate.The whole world is discarded Carbon dioxide sequestration amount is up to 920GT in oil gas field, after carbon dioxide is filled to underground oil gas reservoir, it may occur that following four Seal process up for safekeeping：(1) construction is sealed up for safekeeping, is sealed up for safekeeping by stratum self construction physics with the property of free state；(2) hollow billet is sealed up for safekeeping, with ground Capillary force in plastid micro throat seals CO up for safekeeping2；(3) dissolving is sealed up for safekeeping, CO2It is soluble in water to be sealed up for safekeeping in the form of carbanion CO2；(4) mineral are sealed up for safekeeping, with mineral and dissolving CO2Water flooding reaction seal CO up for safekeeping2。
Domestic and foreign scholars mainly substantially evaluate CO by volumetric method at present2Seal potentiality up for safekeeping, we have discovered that current There is also following problems for carbon sequestration evaluation method：
(1) the Water-rock interaction process occurred during carbon dioxide geologic sequestration is all based at present traditional CO2 seals geochemical reaction process up for safekeeping, speculates CO according only to typical chemical equation2Process is sealed up for safekeeping after filling, this Sample can not represent CO in underground geologic bodies2The reaction process actually occurred after filling, this reaction process it is uncertain into And causes practical geologic body mineral to seal the Potential Evaluation process Minerals amount of sealing up for safekeeping estimated value accuracy up for safekeeping and remain to be discussed.
(2) CO now2In geological storage Potential Evaluation, volumetric method or Monte Carlo method are either used, these are commented Valence method is not only deviating from practically matter condition, and does not consider carbon dioxide occurrence status and Water-rock interaction in the different time stage Influence of the variation to geological storage amount, it is relatively low to be easy to cause evaluation result confidence level, can not accurately instruct practical CO2Envelope Inventory item, it is this to be directed to CO2The geological storage evaluation method of infilling time variation is not set up also.
Accurate CO2Geological storage evaluation method is to influence one of oil-gas Layer carbon dioxide sequestration potentiality key factor, certainly Determine to be directed to oil-gas Layer CO2Seal the design and planning of project up for safekeeping.Therefore.Need to develop it is more acurrate, more meet practical oil-gas Layer CO2Geological storage evaluation method.
Invention content
The present invention for the mineral amount of sealing up for safekeeping estimated value accuracy of the existing technology is poor, evaluation result confidence level is low etc. on Problem is stated, the carbon dioxide geologic sequestration evaluation method that a kind of accuracy is high, evaluation result is with a high credibility is provided.
In order to achieve the above object, the present invention provides a kind of carbon dioxide geologic sequestration evaluation method, specific steps For：
S1, core sample is selected, determines core sample feature；
S2, sample pretreatment：It cuts core sample and surface polishing is carried out to the core sample after cutting；
S3, indoor water rock simulated experiment is carried out to core sample；
S4, the core sample after experiment is handled, the solution after removal experiment reaction；
The core sample that S5, analysis test experiments post-process, step are：
(1) analysis test experiments reaction after core sample, determine experiment after reactive mineral and neogenic mineral type and at Point；
(2) solution ion concentration after the reaction of analysis test experiments；
(3) mineral saturation indices are calculated；
(4) reaction process of indoor water rock simulated experiment is determined；
S6, carry out long-term Numerical Experiment, step is：
(1) under practical geological condition constraints, analog temperature T' simulates the carbon dioxide partial pressure of fillingFor A length of Time when simulation, according to the mineralogical composition feature, hole ratio and formation fluid average chemical of objective area oil reservoir at Dtex is levied, with 1cm3The starting mineral and initial fluid data of unit area oil reservoir, establish numerical simulator；
(2) the long-term CO of simulation oil-gas Layer2Filling；
(3) quantitative assessment CO2Dissolve the amount of sealing up for safekeeping and CO2The mineral amount of sealing up for safekeeping.
Preferably, in step S1, determine core sample feature the specific steps are：
(1) core sample physical property characteristic is determined
The porosity and permeability of oil-gas Layer are determined by rock core physical property measurement；
(2) core sample Characters of Petrology And Mineralogy is determined
By the X-ray diffraction and thin section identification of rock core, oil-gas Layer rock type, mineralogical composition and content, diagenesis are determined Function Characteristics, the Diagenetic include diagenesis types and diagenesis evolution stage.
Preferably, it in step S2, cuts core sample and carries out the specific step of surface polishing to the core sample after cutting Suddenly it is：
(1) core sample is cut into a diameter of 25mm using rock sample cutting machine, the thin rounded flakes that thickness is 3mm；
(2) thick bodyguard rock is carried out to a surface of thin rounded flakes to cut, repeat to throw using 600# granularities abrasive pastes and polishing cloth Light 2-3 times, burnishing surface is rinsed with ultra-pure water 2-3 times repeatedly, until cleaning surface；
(3) surface after being cut to the thick bodyguard rock of thin rounded flakes using 8000# granularities abrasive pastes and polishing cloth carries out accurate brightness The abrasive pastes on thin rounded flakes surface are cleaned in polishing treatment with ultrasonic cleaner；
(4) drying box is utilized to be less than 40 DEG C of low temperature drying thin rounded flakes surfaces；
(5) metal spraying processing is carried out to the thin rounded flakes surface after drying using sputter；
(6) scanning electron microscopic observation thin rounded flakes surface is utilized；
(7) step (3)-step (6) 3-4 times is repeated, until thin rounded flakes surfacing is smooth, recycles 8000# granularities Abrasive pastes and polishing cloth remove the gold-plated film on thin rounded flakes surface.
Preferably, in step S3, to core sample carry out indoor water rock simulated experiment the specific steps are：
(1) indoor water rock simulated experiment parameter is designed, the experiment parameter includes carbon dioxide partial pressureIt is temperature, ultrapure Water volume.
(2) core sample after polishing treatment in step S2 is put into ultra-pure water in central reaction container, it will be central anti- It answers container to be put into heating mantle, opens simultaneously CO on central reaction container2The intake valve and CO at gas inlet end2Gas outlet side Air bleeding valve, be passed through CO2Air in gas 1-2min discharge central reaction containers, closes air bleeding valve, and intake valve is in opening state State is passed through CO2Reach preset carbon dioxide partial pressureIt is inserted into temperature/pressure sensor, temperature parameter is completed at temperature/pressure sensor Setting, real-time monitoring temperature and pressure keep temperature to be divided with carbon dioxideIt is constant.
Preferably, in step S4, to the core sample after experiment handled the specific steps are：After experiment, lead to Excess temperature control pressurer system closes the heating device in heating mantle, closes air bleeding valve and central reaction on dioxide bottle Intake valve on container；Heating mantle is opened, after cooled to room temperature, CO is discharged in the air bleeding valve opened on central reaction container2 Gas, takes out core sample and solution example after reaction, and the core sample after being reacted with ultrapure water removes residual solution In ion interference, then with drying box be less than 40 DEG C of low temperature drying core samples, with injection by the solution example after reaction from It is taken out in central reaction container, it is medium to be tested to be packed into reagent bottle.
Preferably, in step S5, determine reactive mineral and neogenic mineral type and ingredient after experiment the specific steps are：Profit Metal spraying processing is carried out to the core sample surface after drying with sputter, after the reaction of scanning electron microscopic observation core sample Minerals Surface morphology characteristics and newly-generated authigenic mineral morphological feature；Mineral chemical composition is carried out using energy disperse spectroscopy qualitative and quantitative Point analysis determines mineral reaction and newly-generated mineral type and ingredient.
Preferably, in step S5, analysis test experiments reaction after solution ion concentration the specific steps are：It is anti-from center Before core sample after answering container to take out reaction, using the PH of solution after acidometer test reaction, formula (1) and formula are utilized (2) it calculatesWithIon concentration utilizes the cation of solution after plasma emission spectrometer test reaction later Concentration and SiO2Concentration；The expression formula of formula (1) and formula (2) is：
In formula, K is carbonic acid first ionization constant；[H+] it is hydrogen ion concentration, unit：mol/L；For bicarbonate radical Ion concentration, unit：mol/L；K' is carbonic acid secondary ionization constant；For carbon acid ion concentration, unit：mol/L.
Preferably, in step 5, calculate mineral saturation indices the specific steps are：Various mineral are calculated using formula (3) Saturation index, the SI of mineral>0, indicate that mineral tend to precipitate；The SI of mineral<0, indicate that mineral tend to dissolve；The table of formula (3) It is up to formula：
SI=log10 (Q/Keq) n=1,2,3,4......i (3)
In formula, SI is expressed as the saturation index of different minerals, and Q is expressed as the ionic product of certain chemical reaction；Keq is indicated at certain Chemical reaction equilibrium constant；
Determine the reaction process of indoor water rock simulated experiment the specific steps are；Rock sample after comprehensive different simulated experiments Scanning electron microscope analysis, energy spectrum analysis, ion concentration analysis and saturation index analysis, determine short-term CO2Indoor water rock after filling Simulated experiment response path and reaction process.
Preferably, in step S6, the long-term CO of oil-gas Layer is simulated2Filling the specific steps are：Using numerical simulator as base Plinth carries out long-term CO under the dynamics data constraint of different minerals using numerical simulation software2Fill Water-rock interaction Numerical-Mode Draft experiment.
Preferably, quantitative assessment CO2Dissolve the amount of sealing up for safekeeping and CO2The step of mineral amount of sealing up for safekeeping is：
CO2After filling, gas CO2It is dissolved in solution, withWithIonic species is fixed in the solution, and public affairs are passed through Formula (4), formula (5), formula (6) calculate CO2After injection Time in solutionWithIon concentration Assessment of Changes CO2Dissolve the amount of sealing up for safekeeping；Formula (4), formula (5), formula (6) expression formula be：
In formula,It is expressed asMolal quantity knots modification, unit：mol；It is expressed as numerical simulation After experimentMolal quantity, unit：mol；Before being expressed as Numerical ExperimentMolal quantity, unit： mol；It is expressed asMolal quantity knots modification, unit：mol；After being expressed as Numerical Experiment's Molal quantity, unit：mol；Before being expressed as Numerical ExperimentMolal quantity, unit：mol；MsolubilityIt is expressed as CO2Dissolve the amount of sealing up for safekeeping, unit：mg；V is expressed as CO2The oil-gas Layer total volume of filling, unit：m3；
Gas CO2After being dissolved in solution, the Ca in carbanion, with solution is generated2+、Mg2+Ions binding forms carbonate Product calculates CO by formula (7), formula (8), formula (9)2Carbonate mineral mass change evaluates CO after injecting Time2 The mineral amount of sealing up for safekeeping；Formula (7), formula (8), formula (9) expression formula be：
ΔN(calcite)=N'(calcite)-N(calcite) (7)
ΔN(dolomite)=N'(dolomite)-N(dolomite) (8)
Mmineral=(Δ N(calcite)+ΔN(dolomite))×44×V (9)
In formula, Δ N(calcite)It is expressed as the molal quantity knots modification of calcite, unit：mol；N'(calcite)It is expressed as numerical value The molal quantity of calcite, unit after simulated experiment：mol；N(calcite)It is expressed as the molal quantity of calcite before Numerical Experiment, Unit：mol；ΔN(dolomite)It is expressed as the molal quantity knots modification of dolomite, unit：mol；N'(dolomite)It is expressed as numerical simulation The molal quantity of dolomite, unit after experiment：mol；N(dolomite)It is expressed as the molal quantity of dolomite before Numerical Experiment, it is single Position：mol；MmineralIt is expressed as CO2The mineral amount of sealing up for safekeeping, unit：mg.
Compared with prior art, the beneficial effects of the present invention are：
The present invention determines CO by short-term indoor water rock simulated experiment2The response path occurred after filling oil-gas Layer and reaction Process carries out the amount of sealing up for safekeeping evaluation by long-term numerical simulation, makes CO on this basis2Simulation process is filled closer to actual The filling experiment of carbon dioxide scale, avoids and speculates CO according to typical chemical equation2The process of sealing up for safekeeping after filling is made The uncertainty evaluated at the amount of sealing up for safekeeping, hence it is evident that improve the accuracy for sealing Potential Evaluation up for safekeeping for actually sealing geologic body up for safekeeping.
Description of the drawings
Fig. 1 is indoor water rock analogue experiment installation structural schematic diagram described in the embodiment of the present invention.
Fig. 2 a are after Shengli Oil Field of the embodiment of the present invention somewhere thick bodyguard rock of s_4 formation 2938.13m Sandstone Cores samples is cut Core sample surface scan Electronic Speculum characteristic pattern.
Fig. 2 b throw for the s_4 formation 2938.13m Sandstone Cores sample precision brightness of Shengli Oil Field of embodiment of the present invention somewhere Core sample surface scan Electronic Speculum characteristic pattern after light.
Fig. 3 a-c are core sample mineral surfaces Dissolution Characteristics figure after reaction of the embodiment of the present invention.
Fig. 4 is core sample cement corrosion of the embodiment of the present invention and neogenic mineral scanning electron microscope and characteristics of energy spectrum figure.
Fig. 5 is 150 DEG C of temperature of the embodiment of the present invention, PCO24MPa, infilling time Shengli Oil Field somewhere sand four in 10000 Section 2938.13m Sandstone Cores sample numerical value simulation drawings.
Fig. 6 is numerical simulation of embodiment of the present invention porosity and permeability variation figure.
In figure, 1, CO2Gas cylinder, 2, central reaction container, 3, pressure sensor, 4, temperature sensor, 5, heating mantle, 6, temperature Spend control pressurer system.
Specific implementation mode
In the following, the present invention is specifically described by illustrative embodiment.It should be appreciated, however, that not into one In the case of step narration, element, structure and features in an embodiment can also be advantageously incorporated into other embodiment In.
Present invention is disclosed a kind of carbon dioxide geologic sequestration evaluation method, the specific steps are：
S1, core sample is selected, determines core sample feature.Its specific steps are：
(1) core sample physical property characteristic is determined
The porosity and permeability of oil-gas Layer are determined by rock core physical property measurement, porosity is with permeability content referring to table 1.
(2) core sample Characters of Petrology And Mineralogy is determined
By the X-ray diffraction and thin section identification of rock core, oil-gas Layer rock type, mineralogical composition and content, diagenesis are determined Function Characteristics, the Diagenetic include diagenesis types and diagenesis evolution stage.Oil-gas Layer rock type, mine Object ingredient and content are referring to table 1.
Table 1
Content
Quartzy (Vol%) VQ
Potassium feldspar (Vol%) VKF
Albite (Vol%) VAL
Landwaste (Vol%) Vlf
Shale (Vol%) Vm
Cement (Vol%) Vc
Porosity (%) P
Permeability (× 10-3μm2) K
In table 1, VQIt is expressed as the volume fraction of quartz, %；VKFIt is expressed as the volume fraction of potassium feldspar, %；VALIt is expressed as The volume fraction of albite, %；VlfIt is expressed as the volume fraction of landwaste, %；VmIt is expressed as the volume fraction of shale, %；Vc tables It is shown as the volume fraction of cement, %；P is expressed as the porosity of core sample, %；K is expressed as the permeability of core sample, and 1 ×10-3μm2。
S2, sample pretreatment：It cuts core sample and surface polishing is carried out to the core sample after cutting；Its specific steps For：
(1) core sample is cut into a diameter of 25mm using rock sample cutting machine, the thin rounded flakes that thickness is 3mm；
(2) thick bodyguard rock is carried out to a surface of thin rounded flakes to cut, repeat to throw using 600# granularities abrasive pastes and polishing cloth Light 2-3 times, burnishing surface is rinsed with ultra-pure water 2-3 times repeatedly, until cleaning surface；
(3) surface after being cut to the thick bodyguard rock of thin rounded flakes using 8000# granularities abrasive pastes and polishing cloth carries out accurate brightness The abrasive pastes on thin rounded flakes surface are cleaned in polishing treatment with ultrasonic cleaner；
(4) drying box is utilized to be less than 40 DEG C of low temperature drying thin rounded flakes surfaces；
(5) metal spraying processing is carried out to the thin rounded flakes surface after drying using sputter；
(6) scanning electron microscopic observation thin rounded flakes surface is utilized；
(7) step (3)-step (6) 3-4 times is repeated, until thin rounded flakes surfacing is smooth, recycles 8000# granularities Abrasive pastes and polishing cloth remove the gold-plated film on thin rounded flakes surface.
S3, indoor water rock simulated experiment is carried out to core sample；Its specific steps are：
(1) indoor water rock simulated experiment parameter is designed, the experiment parameter includes carbon dioxide partial pressureIt is temperature, ultrapure Water volume.Indoor water rock simulated experiment parameter is referring to table 2.
Table 2
(2) core sample after polishing treatment in step S2 is put into ultra-pure water in central reaction container, it will be central anti- It answers container to be put into heating mantle, opens simultaneously CO on central reaction container2The intake valve and CO at gas inlet end2Gas outlet side Air bleeding valve, be passed through CO2Air in gas 1-2min discharge central reaction containers, closes air bleeding valve, and intake valve is in opening state State is passed through CO2Reach preset carbon dioxide partial pressureIt is inserted into temperature/pressure sensor, temperature parameter is completed at temperature/pressure sensor Setting, real-time monitoring temperature and pressure keep temperature to be divided with carbon dioxideIt is constant.
Device involved by above-mentioned simulating lab test mainly has referring to Fig. 1：CO2Gas cylinder 1, for providing CO2Gas；In Reaction vessel 2 is entreated, for accommodating rock sample, CO2And solution, and be rock sample, CO2Reacting environment is provided with solution；Pressure Sensor 3, the pressure for measuring central reaction container 2；Temperature sensor 4, the temperature for measuring central reaction container 2； Heating mantle 5 is heated for placing central reaction container 2 for central reaction container 2；Temperature, pressure control system 6, for controlling The temperature and pressure of central reaction container 2.Preferably, the temperature, pressure control system 6 uses temperature and pressure controller.
S4, the core sample after experiment is handled, the solution after removal experiment reaction；Its specific steps are：
After experiment, the heating device in heating mantle is closed by temperature, pressure control system, closes carbon dioxide gas The intake valve on air bleeding valve and central reaction container on bottle；Heating mantle is opened, after cooled to room temperature, opens central reaction CO is discharged in air bleeding valve on container2Gas takes out core sample and solution example after reaction, after being reacted with ultrapure water Core sample removes the ion interference in residual solution, then is less than 40 DEG C of low temperature drying core samples with drying box, uses injection Solution example after reaction is taken out from central reaction vessel, it is medium to be tested to be packed into reagent bottle.
The core sample that S5, analysis test experiments post-process, step are：
(1) analysis test experiments reaction after core sample, determine experiment after reactive mineral and neogenic mineral type and at Point；Its specific steps are：
Metal spraying processing, scanning electron microscopic observation core sample chats are carried out to the core sample surface after drying using sputter Surface morphology characteristics after object reaction and newly-generated authigenic mineral morphological feature；Mineral chemical composition is carried out using energy disperse spectroscopy Qualitative and quantitative point analysis, determines mineral reaction and newly-generated mineral type and ingredient.
(2) solution ion concentration after the reaction of analysis test experiments；Its specific steps are：
Before the core sample that central reaction vessel takes out after reaction, using the PH of solution after acidometer test reaction, It is calculated using formula (1) and formula (2)WithIon concentration utilizes plasma emission spectrometer to test anti-later The cation concn and SiO of solution after answering2Concentration, the cation concn and SiO of solution after reaction2Concentration is referring to table 3.Formula (1) Expression formula with formula (2) is：
In formula, K is carbonic acid first ionization constant；[H+] it is hydrogen ion concentration, unit：mol/L；For bicarbonate radical Ion concentration, unit：mol/L；K' is carbonic acid secondary ionization constant；For carbon acid ion concentration, unit：mol/L.
Table 3
(3) mineral saturation indices are calculated；Its specific steps are：
The saturation index of various mineral, the SI of mineral are calculated using formula (3)>0, indicate that mineral tend to precipitate；Mineral SI<0, indicate that mineral tend to dissolve；The expression formula of formula (3) is：
SI=log10 (Q/Keq) n=1,2,3,4......i (3)
In formula, SI is expressed as the saturation index of different minerals, and Q is expressed as the ionic product of certain chemical reaction；Keq is indicated at certain Chemical reaction equilibrium constant.
100 DEG C with 150 DEG C under the conditions of mineral equilibrium constant K eq referring to table 4.The saturation index of differential responses mineral referring to Table 5.
Table 4
Mineral type Mineral dissolution reacts 100℃ 150℃
Potassium feldspar KAlSi3O8+8H2O=K++Al(OH)4-+3H4SiO4 -18.54 -16.49
Albite NaAlSi3O8+8H2O=Na++Al(OH)4 -+3H4SiO4 -15.75 -14.3
Quartz SiO2+2H2O=H4SiO4 0 -2.98 -2.61
Kaolinite Al2Si2O5(OH)4+7H2O=2Al (OH)4 -+2H++2H4SiO4 0 -33.72 -30.47
Dolomite CaMg(CO3)2+2H+=Ca2++Mg2++2HCO3 - 0.76 -0.76
Calcite CaCO3+H+=Ca2++HCO3 - 0.83 0.15
Ferrocalcite CaFe(CO3)2+2H+=Ca2++Fe2++2HCO3 - -2.09 -3.31
Boehmite AlO2H+3H+=Al+3+2H2O 4.849 2.354
Table 5
(4) reaction process of indoor water rock simulated experiment is determined；Its specific steps are：
The scanning electron microscope analysis of rock sample, energy spectrum analysis, ion concentration analysis and saturation after comprehensive different simulated experiments Index analysis determines short-term CO2Indoor water rock simulated experiment response path and reaction process after filling.
S6, carry out long-term Numerical Experiment, step is：
(1) under practical geological condition constraints, analog temperature T' simulates the carbon dioxide partial pressure of fillingFor A length of Time when simulation, according to the mineralogical composition feature, hole ratio and formation fluid average chemical of objective area oil reservoir at Dtex is levied, with 1cm3The starting mineral and initial fluid data of unit area oil reservoir, establish numerical simulator.1cm3Unit plane The starting mineral and initial fluid data of long-pending oil layer are referring to table 6.
Table 6
(2) the long-term CO of simulation oil-gas Layer2Filling；Its specific steps are：
Based on numerical simulator, different minerals dynamics data constraint under, using numerical simulation software into The long-term CO of row2Fill Water-rock interaction Numerical Experiment.The dynamics data of different minerals is referring to table 7 in Numerical Experiment.
Table 7
(3) quantitative assessment CO2Dissolve the amount of sealing up for safekeeping and CO2The mineral amount of sealing up for safekeeping.Its specific steps are：
CO2After filling, gas CO2It is dissolved in solution, withWithIonic species is fixed in the solution, and public affairs are passed through Formula (4), formula (5), formula (6) calculate CO2After injection Time in solutionWithIon concentration Assessment of Changes CO2Dissolve the amount of sealing up for safekeeping；Formula (4), formula (5), formula (6) expression formula be：
In formula,It is expressed asMolal quantity knots modification, unit：mol；It is expressed as numerical simulation reality After testingMolal quantity, unit：mol；Before being expressed as Numerical ExperimentMolal quantity, unit：mol；It is expressed asMolal quantity knots modification, unit：mol；After being expressed as Numerical ExperimentMole Number, unit：mol；Before being expressed as Numerical ExperimentMolal quantity, unit：mol；MsolubilityIt is expressed as CO2 Dissolve the amount of sealing up for safekeeping, unit：mg；V is expressed as CO2The oil-gas Layer total volume of filling, unit：m3；
Gas CO2After being dissolved in solution, the Ca in carbanion, with solution is generated2+、Mg2+Ions binding forms carbonate Product calculates CO by formula (7), formula (8), formula (9)2Carbonate mineral mass change evaluates CO after injecting Time2 The mineral amount of sealing up for safekeeping；Formula (7), formula (8), formula (9) expression formula be：
ΔN(calcite)=N'(calcite)-N(calcite) (7)
ΔN(dolomite)=N'(dolomite)-N(dolomite) (8)
Mmineral=(Δ N(calcite)+ΔN(dolomite))×44×V (9)
In formula, Δ N(calcite)It is expressed as the molal quantity knots modification of calcite, unit：mol；N'(calcite)It is expressed as numerical value The molal quantity of calcite, unit after simulated experiment：mol；N(calcite)It is expressed as the molal quantity of calcite before Numerical Experiment, Unit：mol；ΔN(dolomite)It is expressed as the molal quantity knots modification of dolomite, unit：mol；N'(dolomite)It is expressed as numerical simulation The molal quantity of dolomite, unit after experiment：mol；N(dolomite)It is expressed as the molal quantity of dolomite before Numerical Experiment, it is single Position：mol；MmineralIt is expressed as CO2The mineral amount of sealing up for safekeeping, unit：mg.
It is with indoor Water-rock interaction and Numerical Experiment in the s_4 formation 2938.13m sandstone of Shengli Oil Field somewhere below Example further illustrates the present invention.
S1, core sample is selected, determines core sample feature.Its specific steps are：
(1) core sample physical property characteristic is determined
The porosity and permeability of Shengli Oil Field somewhere s_4 formation 2938.13m sandstone are determined by rock core physical property measurement, Porosity is with permeability content referring to table 1.
(2) core sample Characters of Petrology And Mineralogy is determined
By the X-ray diffraction and thin section identification of rock core, Shengli Oil Field somewhere s_4 formation 2938.13m sandstone are determined Rock type, mineralogical composition and content, Diagenetic, the Diagenetic include that diagenesis types and diagenesis are made Use the evolutionary phase.Oil-gas Layer rock type, mineralogical composition and content are referring to table 8.
Table 8
Content
Quartzy (Vol%) 38
Potassium feldspar (Vol%) 21
Albite (Vol%) 20
Landwaste (Vol%) 21.5
Shale (Vol%) 0.5
Cement (Vol%) 25
Porosity (%) 4.4
Permeability (× 10-3μm2) 0.165
S2, sample pretreatment：It cuts core sample and surface polishing is carried out to the core sample after cutting；Its specific steps For：
(1) core sample is cut into a diameter of 25mm using HBQ-4 type rock sample cutting machines, thickness is the round thin of 3mm Piece；
(2) thick bodyguard rock is carried out to a surface of thin rounded flakes using 600# granularity someiyoshine abrasive pastes and polishing cloth to cut (referring to Fig. 2 a) repeats polishing 2-3 times, rinses burnishing surface repeatedly with ultra-pure water 2-3 times, until cleaning surface；
(3) surface after being cut to the thick bodyguard rock of thin rounded flakes using 8000# granularity someiyoshine abrasive pastes and polishing cloth is carried out The abrasive pastes on thin rounded flakes surface are cleaned in accurate brightness polishing treatment with KH2200DB table type numerical control ultrasonic cleaners；
(4) one perseverance DHG-9145A air dry ovens of Shanghai are utilized to be less than 40 DEG C of low temperature drying thin rounded flakes surfaces；
(5) metal spraying processing is carried out to the thin rounded flakes surface after drying using SPT-20 small ions sputter；
(6) scanning electron microscopic observation thin rounded flakes surface is utilized；
(7) step (3)-step (6) 3-4 times is repeated, until thin rounded flakes surfacing is smooth (referring to Fig. 2 b), is recycled 8000# granularities abrasive pastes and polishing cloth remove the gold-plated film on thin rounded flakes surface.
S3, indoor water rock simulated experiment is carried out to core sample；Its specific steps are：
(1) indoor water rock simulated experiment parameter is designed, the experiment parameter includes carbon dioxide partial pressureIt is temperature, ultrapure Water volume, wherein ultra-pure water uses resistance>The water of 18.5M Ω.Indoor water rock simulated experiment parameter is referring to table 9.
Table 9
(2) core sample and the 200ml ultra-pure waters after polishing treatment in step S2 are put into central reaction container, will in Centre reaction vessel is put into heating mantle, opens simultaneously CO on central reaction container2The intake valve and CO at gas inlet end2Gas goes out The air bleeding valve at gas end, is passed through CO2Air in gas 1-2min discharge central reaction containers, closes air bleeding valve, and intake valve is in open State is opened, CO is passed through2Reach 4MPa；It is inserted into temperature/pressure sensor, temperature parameter setting is completed at temperature/pressure sensor, is monitored in real time Temperature and pressure keep temperature to be divided with carbon dioxideIt is constant.Wherein, central reaction container uses temperature range 0-200 ℃、PCO2The central reaction container of range 0-10MPa, heating mantle use the heating mantle of 0-200 DEG C of temperature range, temperature/pressure sensor Using 0-200 DEG C of temperature range, PCO2The temperature/pressure sensor of range 0-50MPa, CO2Gas cylinder uses PCO2For the CO of 10MPa2Gas cylinder.
S4, the core sample after experiment is handled, the solution after removal experiment reaction；Its specific steps are：
After experiment, the heating device in heating mantle is closed by temperature, pressure control system, closes carbon dioxide gas The intake valve on air bleeding valve and central reaction container on bottle；Heating mantle is opened, after cooled to room temperature, opens central reaction CO is discharged in air bleeding valve on container2Gas takes out core sample and solution example after reaction, after being reacted with ultrapure water Core sample removes the ion interference in residual solution, then is less than 40 DEG C of low temperature with one perseverance DHG-9145A air dry ovens of Shanghai Core sample is dried, the solution example after reaction is taken out from central reaction vessel with injection, is fitted into reagent bottle and waits for Test.
The core sample that S5, analysis test experiments post-process, step are：
(2) analysis test experiments reaction after core sample, determine experiment after reactive mineral and neogenic mineral type and at Point；Its specific steps are：
Metal spraying processing is carried out to the core sample surface after drying using SPT-20 small ions sputter, scanning electron microscope is seen Core sample is examined, at 100 DEG C, 4MPa, under the conditions of 2 days reaction time, it is seen that with potassium feldspar faint corrosion, quartz occur for albite Corrosion does not occur (referring to Fig. 3 a)；At 150 DEG C, 4MPa, 2 days reaction time condition, it is seen that albite surface generates a large amount of corrosions Hole, potassium feldspar surface see that faint corrosion (referring to Fig. 3 b) occurs for a small amount of pocket, quartz；At 150 DEG C, 4MPa, 4 days reaction time Under the conditions of, it is seen that albite stitches strong corrosion along joint, and potassium feldspar surface generates a small amount of pocket (referring to figure with quartz surfaces 3c).Strong corrosion (referring to Fig. 4 a, Fig. 4 b) all occurs for the dolomite cement in three groups of experiments.Utilize Brooker 430-M power spectrums Instrument carries out qualitative and quantitative point analysis to neogenic mineral chemical composition, under 150 DEG C, 4MPa, 2 days reaction time experiment condition, It can be seen that the small hexagon sheet kaolinite of feldspar Surface Creation (referring to Fig. 4 c)；With the increase in reaction time and temperature, 150 DEG C, 4MPa, under the conditions of 4 days reaction time, it is seen that larger hexagon sheet kaolinite (referring to Fig. 4 d), and in quartz and length The ferrocalcite (referring to Fig. 4 e) of the visible diamond shape of stone particle surface, the visible lath-shaped calcite of intergranular (referring to Fig. 4 f).
(2) solution ion concentration after the reaction of analysis test experiments；Its specific steps are：
Before the core sample that central reaction vessel takes out after reaction, using the PH of solution after acidometer test reaction, It is calculated using formula (1) and formula (2)WithIon concentration utilizes plasma emission spectrometer to test anti-later The cation concn and SiO of solution after answering2Concentration.The cation concn of solution, SiO after reaction2Concentration,Ion concentration WithIon concentration is referring to table 10.
Table 10
(3) mineral saturation indices are calculated；Its specific steps are：
The saturation index of various mineral, the SI of mineral are calculated using formula (3)>0, indicate that mineral tend to precipitate；Mineral SI<0, indicate that mineral tend to dissolve.The saturation index of differential responses mineral is referring to table 11.
Table 11
Experiments K-feldspar Albite Quartz Kaolinite Dolomite Calcite Ferrocalcite Boehmite
100℃-2days -5.5931 -8.2132 -1.0848 1.5589 -5.8016 -1.0640 -2.1399 1.0662
150℃-2days -3.9044 -6.3825 -0.4617 0.5251 -5.1885 -0.7397 -1.5551 4.4762
150℃-4days -1.1618 -3.5727 -0.0605 2.4860 -3.2730 0.5024 1.9030 7.0555
(4) reaction process of indoor water rock simulated experiment is determined；Its specific steps are：
The scanning electron microscope analysis of rock sample, energy spectrum analysis, ion concentration analysis and saturation after comprehensive different simulated experiments Index analysis determines short-term CO2Indoor water rock simulated experiment response path and reaction process after filling.
S6, it is had determined that in by the test analysis of indoor water rock simulated experiment core sample short-term S5Fixed numbers At a temperature of and, CO2Significance path and process carry out Numerical Experiment and inquire into long-term filling process.Carry out length The step of phase Numerical Experiment is：
(1) under practical geological condition constraints, analog temperature is 150 DEG C, and the carbon dioxide partial pressure for simulating filling is 4MPa, A length of 10000 when simulation, according to the mineralogical composition feature, hole ratio and formation fluid average chemical of objective area oil reservoir at Dtex is levied, with 1cm3The starting mineral and initial fluid data of unit area oil reservoir, establish numerical simulator.1cm3Unit plane The starting mineral and initial fluid data of long-pending oil layer are referring to table 12.
Table 12
(2) the long-term CO of simulation oil-gas Layer2Filling；Has CO establishing2The representative numerical simulator basis of geological storage body On, in conjunction with Petrogenetic Simulation CO2Injection front and back physical property, ion concentration, mineral quality, mineral saturation indices variation characteristic, are grown Phase CO2Fill numerical simulation.Its specific steps are：
Based on numerical simulator, under the dynamics data constraint of different minerals, Geochemist's is utilized Workbench (GWB) 9.0 carries out long-term CO2Fill Water-rock interaction Numerical Experiment.
Stage I (0-50)：To dissolve the CO based on sealing up for safekeeping2Filling process
CO2It is dissolved in solution after filling, water flooding PH shows as acid (PH=3.97) (referring to Fig. 5 f), dolomite corrosion Corrosion (referring to Fig. 5 d, 5e) occurs first, Ca in solution2+、Mg2+(referring to Fig. 5 a) increases.With K+、Na+Ionic activity is continuous Increase (referring to Fig. 5 c), corrosion (referring to Fig. 5 d, 5e) occurs for a large amount of potassium feldspars, albite.Kaolinite, quartz precipitation are (referring to figure 5d, 5e) lead to formation pore and permeability reduction (referring to Fig. 6), CO in this stage2It seals up for safekeeping mainly to be dissolved in water flooding generation HCO3 -And CO3 2-Based on, the dissolving amount of sealing up for safekeeping accounts for mineral and seals up for safekeeping and seal the 94.46% of summation up for safekeeping with dissolving up to 3860Mt.
Stage II (50-3960)：With mineral seal up for safekeeping based on CO2Filling process
In this stage, potassium feldspar continues that corrosion occurs with albite, and feldspar corrosion causes solution PH to be continuously increased (referring to figure 5f).With CO2Constantly Ca in filling and solution2+、Mg2+Ionic activity reduces (referring to Fig. 5 c), and calcite mineral largely precipitates, Ca2+、Mg2+Ion constantly reduces, this stage feldspar corrosion promotes carbonate mineral precipitation (referring to Fig. 5 d, 5e), formation pore It is further decreased (referring to Fig. 6) with permeability, the mineral amount of sealing up for safekeeping can account for 49.93% up to 3980Mt.
Stage III (3960-10000)：Carbon dioxide-water flooding-mineral react equilibrium process
Each ion concentration, PH, mineral quality no longer change in water flooding, and Water-rock interaction reaches balance.
(3) quantitative assessment CO2Dissolve the amount of sealing up for safekeeping and CO2The mineral amount of sealing up for safekeeping.Its specific steps are：
CO2After filling, gas CO2It is dissolved in solution, withWithIonic species is fixed in the solution, and public affairs are passed through Formula (4), formula (5), formula (6) calculate CO2Solution after injection 10000WithIon concentration Assessment of Changes CO2 The dissolving amount of sealing up for safekeeping is 4328.2Mt.
Gas CO2After being dissolved in solution, the Ca in carbanion, with solution is generated2+、Mg2+Ions binding forms carbonate Product calculates CO by formula (7), formula (8), formula (9)2Carbonate mineral mass change evaluation after injection 10000 CO2The mineral amount of sealing up for safekeeping is 4210Mt.
Existing CO2Mineral are sealed up for safekeeping in Potential Evaluation, with most widely sharp to North America geological conditions for U.S. Department of Energy The value that Monte Carlo method obtains, specific method are：The potassium feldspar amount of sealing up for safekeeping is calculated by formula (10), is counted by formula (11) Calculate the albite amount of sealing up for safekeeping, formula (10), formula (11) expression formula be：
In formula,It is expressed as the mineral amount of sealing up for safekeeping of potassium feldspar, unit：mg；ρK-feldsparIt is expressed as potassium feldspar Density, unit：kg/m3；VK-feldsparIt is expressed as the volume of potassium feldspar, unit：m3；MK-feldsparIt is expressed as rubbing for potassium feldspar That quality, unit：g/mol；It is expressed as CO2Molal weight, unit：g/mol；A is expressed as CO2Store coefficient of efficiency； It is expressed as CO2The mineral amount of sealing up for safekeeping, unit：mg；It is expressed as the mineral amount of sealing up for safekeeping of albite, unit：mg；ρAlbiteTable It is shown as the density of albite, unit：kg/m3；VAlbiteIt is expressed as the volume of albite, unit：cm3；MAlbiteIt is expressed as albite Molal weight, unit：g/mol.
In order to illustrate the effect of the present invention, it is ancient that Dongying Depression is calculated using existing above-mentioned Monte Carlo method and the method for the present invention Close is the reservoir Minerals amount of sealing up for safekeeping.The potassium feldspar and the albite mineral amount of sealing up for safekeeping being calculated using existing Monte Carlo EGS4 method are most Height only has 3814.3Mt, and potassium feldspar is calculated using method for numerical simulation of the present invention and the albite mineral amount of sealing up for safekeeping is reachable 4210Mt.Mineral are sealed up for safekeeping under the conditions of the mineral amount of sealing up for safekeeping that widely used Monte Carlo method is calculated is much smaller than practical geology Amount, the mineral amount of sealing up for safekeeping computational methods used today are not particularly suited for different geologic body CO2Seal evaluation up for safekeeping.Therefore, the present invention carries The method of confession significantly improves CO2Seal Potential Evaluation accuracy and confidence level up for safekeeping.
Embodiment provided above only with illustrating the present invention for convenience, and it is not intended to limit the protection scope of the present invention, Technical solution scope of the present invention, person of ordinary skill in the field make various simple deformations and modification, should all include In the above claim.

Claims (10)

1. a kind of carbon dioxide geologic sequestration evaluation method, which is characterized in that its specific steps are：
S1, core sample is selected, determines core sample feature；
S2, sample pretreatment：It cuts core sample and surface polishing is carried out to the core sample after cutting；
S3, indoor water rock simulated experiment is carried out to core sample；
S4, the core sample after experiment is handled, the solution after removal experiment reaction；
The core sample that S5, analysis test experiments post-process, step are：
(1) core sample after the reaction of analysis test experiments, determines reactive mineral and neogenic mineral type and ingredient after experiment；
(2) solution ion concentration after the reaction of analysis test experiments；
(3) mineral saturation indices are calculated；
(4) reaction process of indoor water rock simulated experiment is determined；
S6, carry out long-term Numerical Experiment, step is：
(1) under practical geological condition constraints, analog temperature T' simulates the carbon dioxide partial pressure of fillingForSimulation Shi Changwei Time, it is special according to the mineralogical composition feature of objective area oil reservoir, hole ratio and formation fluid average chemical constitution Sign, with 1cm3The starting mineral and initial fluid data of unit area oil reservoir, establish numerical simulator；
(2) the long-term CO of simulation oil-gas Layer2Filling；
(3) quantitative assessment CO2Dissolve the amount of sealing up for safekeeping and CO2The mineral amount of sealing up for safekeeping.
2. carbon dioxide geologic sequestration evaluation method as described in claim 1, which is characterized in that in step S1, determine rock core Sample characteristic the specific steps are：
(1) core sample physical property characteristic is determined
The porosity and permeability of oil-gas Layer are determined by rock core physical property measurement；
(2) core sample Characters of Petrology And Mineralogy is determined
By the X-ray diffraction and thin section identification of rock core, oil-gas Layer rock type, mineralogical composition and content, diagenesis are determined Feature, the Diagenetic include diagenesis types and diagenesis evolution stage.
3. carbon dioxide geologic sequestration evaluation method as claimed in claim 2, which is characterized in that in step S2, cut rock core Sample and to after cutting core sample carry out surface polishing the specific steps are：
(1) core sample is cut into a diameter of 25mm using rock sample cutting machine, the thin rounded flakes that thickness is 3mm；
(2) thick bodyguard rock is carried out to a surface of thin rounded flakes to cut, repeat to polish 2-3 using 600# granularities abrasive pastes and polishing cloth It is secondary, burnishing surface is rinsed repeatedly with ultra-pure water 2-3 times, until cleaning surface；
(3) surface after being cut to the thick bodyguard rock of thin rounded flakes using 8000# granularities abrasive pastes and polishing cloth carries out accurate brightness polishing Processing, the abrasive pastes on thin rounded flakes surface are cleaned with ultrasonic cleaner；
(4) drying box is utilized to be less than 40 DEG C of low temperature drying thin rounded flakes surfaces；
(5) metal spraying processing is carried out to the thin rounded flakes surface after drying using sputter；
(6) scanning electron microscopic observation thin rounded flakes surface is utilized；
(7) step (3)-step (6) 3-4 times is repeated, until thin rounded flakes surfacing is smooth, recycles the grinding of 8000# granularities Cream and polishing cloth remove the gold-plated film on thin rounded flakes surface.
4. carbon dioxide geologic sequestration evaluation method as claimed in claim 3, which is characterized in that in step S3, to core sample Product carry out indoor water rock simulated experiment the specific steps are：
(1) indoor water rock simulated experiment parameter is designed, the experiment parameter includes carbon dioxide partial pressureTemperature, ultrapure water body Product.
(2) core sample after polishing treatment in step S2 is put into ultra-pure water in central reaction container, central reaction is held Device is put into heating mantle, opens simultaneously CO on central reaction container2The intake valve and CO at gas inlet end2The row of gas outlet side Air valve is passed through CO2Air in gas 1-2min discharge central reaction containers, closes air bleeding valve, and intake valve is in open state, is led to Enter CO2Reach preset carbon dioxide partial pressureIt is inserted into temperature/pressure sensor, temperature parameter setting is completed at temperature/pressure sensor, Real-time monitoring temperature and pressure, keep temperature to be divided with carbon dioxideIt is constant.
5. carbon dioxide geologic sequestration evaluation method as claimed in claim 4, which is characterized in that in step S4, after experiment Core sample handled the specific steps are：After experiment, closed in heating mantle by temperature, pressure control system Heating device closes the air bleeding valve on dioxide bottle and the intake valve on central reaction container；Heating mantle is opened, it is naturally cold But CO is discharged to the air bleeding valve after room temperature, opened on central reaction container2Gas takes out core sample and solution sample after reaction Product, the core sample after being reacted with ultrapure water remove the ion interference in residual solution, then low less than 40 DEG C with drying box Temperature drying core sample, the solution example after reaction is taken out from central reaction vessel, it is medium to be packed into reagent bottle with injection It is to be tested.
6. carbon dioxide geologic sequestration evaluation method as claimed in claim 5, which is characterized in that in step S5, determine experiment Afterwards reactive mineral and neogenic mineral type and ingredient the specific steps are：Using sputter to the core sample surface after drying into The processing of row metal spraying, the surface morphology characteristics after the reaction of scanning electron microscopic observation core sample Minerals and newly-generated authigenic mineral shape State feature；Qualitative and quantitative point analysis is carried out to mineral chemical composition using energy disperse spectroscopy, determines mineral reaction and newly-generated mineral Type and ingredient.
7. carbon dioxide geologic sequestration evaluation method as claimed in claim 6, which is characterized in that in step S5, analysis test Experiment reaction after solution ion concentration the specific steps are：Before the core sample that central reaction vessel takes out after reaction, Using the PH of solution after acidometer test reaction, calculated using formula (1) and formula (2)WithIon concentration, it The cation concn and SiO of solution after plasma emission spectrometer test reaction are utilized afterwards2Concentration；Formula (1) and formula (2) Expression formula be：
In formula, K is carbonic acid first ionization constant；[H+] it is hydrogen ion concentration, unit：mol/L；For bicarbonate ion Concentration, unit：mol/L；K' is carbonic acid secondary ionization constant；For carbon acid ion concentration, unit：mol/L.
8. carbon dioxide geologic sequestration evaluation method as claimed in claim 7, which is characterized in that in step 5, it is full to calculate mineral With index the specific steps are：The saturation index of various mineral, the SI of mineral are calculated using formula (3)>0, indicate that mineral tend to Precipitation；The SI of mineral<0, indicate that mineral tend to dissolve；The expression formula of formula (3) is：
SI=log10 (Q/Keq) n=1,2,3,4......i (3)
In formula, SI is expressed as the saturation index of different minerals, and Q is expressed as the ionic product of certain chemical reaction；Keq is indicated in certain chemistry Reaction equilibrium constant；
Determine the reaction process of indoor water rock simulated experiment the specific steps are；Rock sample sweeps after comprehensive different simulated experiments Electron microscope analysis, energy spectrum analysis, ion concentration analysis and saturation index analysis are retouched, determines short-term CO2Indoor water rock is simulated after filling Test response path and reaction process.
9. carbon dioxide geologic sequestration evaluation method as claimed in claim 8, which is characterized in that in step S6, simulate oil gas The long-term CO of layer2Filling the specific steps are：Based on numerical simulator, under the dynamics data constraint of different minerals, Long-term CO is carried out using numerical simulation software2Fill Water-rock interaction Numerical Experiment.
10. carbon dioxide geologic sequestration evaluation method as claimed in claim 8, which is characterized in that in step S6, quantitative assessment CO2Dissolve the amount of sealing up for safekeeping and CO2The step of mineral amount of sealing up for safekeeping is：
CO2After filling, gas CO2It is dissolved in solution, withWithIonic species is fixed in the solution, and formula is passed through (4), formula (5), formula (6) calculate CO2After injection Time in solutionWithIon concentration Assessment of Changes CO2 Dissolve the amount of sealing up for safekeeping；Formula (4), formula (5), formula (6) expression formula be：
In formula,It is expressed asMolal quantity knots modification, unit：mol；After being expressed as Numerical ExperimentMolal quantity, unit：mol；Before being expressed as Numerical ExperimentMolal quantity, unit：mol；It is expressed asMolal quantity knots modification, unit：mol；After being expressed as Numerical ExperimentRub That number, unit：mol；Before being expressed as Numerical ExperimentMolal quantity, unit：mol；MsolubilityIt is expressed as CO2Dissolve the amount of sealing up for safekeeping, unit：mg；V is expressed as CO2The oil-gas Layer total volume of filling, unit：m3；
Gas CO2After being dissolved in solution, the Ca in carbanion, with solution is generated2+、Mg2+Ions binding forms carbonate product, CO is calculated by formula (7), formula (8), formula (9)2Carbonate mineral mass change evaluates CO after injecting Time2Mineral The amount of sealing up for safekeeping；Formula (7), formula (8), formula (9) expression formula be：
ΔN(calcite)=N'(calcite)-N(calcite) (7)
ΔN(dolomite)=N'(dolomite)-N(dolomite) (8)
Mmineral=(Δ N(calcite)+ΔN(dolomite))×44×V (9)
In formula, Δ N(calcite)It is expressed as the molal quantity knots modification of calcite, unit：mol；N'(calcite)It is expressed as numerical simulation The molal quantity of calcite, unit after experiment：mol；N(calcite)It is expressed as the molal quantity of calcite before Numerical Experiment, unit： mol；ΔN(dolomite)It is expressed as the molal quantity knots modification of dolomite, unit：mol；N'(dolomite)It is expressed as Numerical Experiment The molal quantity of dolomite afterwards, unit：mol；N(dolomite)It is expressed as the molal quantity of dolomite before Numerical Experiment, unit： mol；MmineralIt is expressed as CO2The mineral amount of sealing up for safekeeping, unit：mg.
CN201810393143.1A 2018-04-27 2018-04-27 A kind of carbon dioxide geologic sequestration evaluation method Pending CN108614076A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810393143.1A CN108614076A (en) 2018-04-27 2018-04-27 A kind of carbon dioxide geologic sequestration evaluation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810393143.1A CN108614076A (en) 2018-04-27 2018-04-27 A kind of carbon dioxide geologic sequestration evaluation method

Publications (1)

Publication Number Publication Date
CN108614076A true CN108614076A (en) 2018-10-02

Family

ID=63661057

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810393143.1A Pending CN108614076A (en) 2018-04-27 2018-04-27 A kind of carbon dioxide geologic sequestration evaluation method

Country Status (1)

Country Link
CN (1) CN108614076A (en)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109828099A (en) * 2019-01-28 2019-05-31 中国石油大学(华东) A kind of rock-water-hydrocarbon analogue experiment method
CN110133223A (en) * 2019-05-20 2019-08-16 中国石油大学(华东) A kind of dawsonite is in CO2Preservation condition evaluation system during geological reservoir
CN111855715A (en) * 2020-07-10 2020-10-30 中国石油天然气股份有限公司 Method and system for simulating and evaluating water-rock reaction of dolostone under salt
CN112924648A (en) * 2021-01-26 2021-06-08 大连理工大学 Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity
CN113610914A (en) * 2021-07-02 2021-11-05 中铁二十局集团有限公司 Sandstone pore structure analysis method, device, equipment and storage medium
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
CN114781287A (en) * 2022-06-13 2022-07-22 武汉恒信楚才科技有限公司 Method and system for calculating natural gas-water-rock interaction under high-temperature and high-pressure conditions
CN115165700A (en) * 2022-06-17 2022-10-11 中国石油大学(北京) Calculation method and experimental device for full-core-scale carbon dioxide adsorption and storage amount
CN115902160A (en) * 2022-11-28 2023-04-04 新疆敦华绿碳技术股份有限公司 Evaluation of CO in saline water reservoir geology 2 Method for sequestration of potential
WO2023071552A1 (en) * 2021-10-27 2023-05-04 中国华能集团清洁能源技术研究院有限公司 Gas solubility measuring apparatus and method, and storage potential predicting system and method
CN116255198A (en) * 2023-02-13 2023-06-13 大连理工大学 Carbon dioxide sealing method based on reservoir wettability optimization design and layered regulation
CN117250133A (en) * 2023-11-16 2023-12-19 国建大数据科技(辽宁)有限公司 Smart city large-scale data acquisition method and system
CN117405718A (en) * 2023-12-12 2024-01-16 山东石油化工学院 Quantitative evaluation method for carbon sequestration capacity of magma calcium element based on XRF scanning
CN117686689A (en) * 2024-02-02 2024-03-12 中国科学院地质与地球物理研究所 Basalt carbon sequestration efficiency measurement method based on basalt carbon sequestration potential evaluation

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102313790A (en) * 2011-07-19 2012-01-11 北京师范大学 Submarine geologic body carbon dioxide sequestration potential assessment method
CN103544361A (en) * 2013-11-04 2014-01-29 西北大学 Evaluation method of CO2 geological sequestration potential in oil-gas field development
CN104794329A (en) * 2015-04-02 2015-07-22 中国科学院地质与地球物理研究所 Method for evaluating storage potential of place scale CO<2> salt water layer

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102313790A (en) * 2011-07-19 2012-01-11 北京师范大学 Submarine geologic body carbon dioxide sequestration potential assessment method
CN103544361A (en) * 2013-11-04 2014-01-29 西北大学 Evaluation method of CO2 geological sequestration potential in oil-gas field development
CN104794329A (en) * 2015-04-02 2015-07-22 中国科学院地质与地球物理研究所 Method for evaluating storage potential of place scale CO<2> salt water layer

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
于淼: "CO2驱油后的流体-岩石相互作用的研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *
张冠儒: "咸水层碳封存中CO2迁移转化过程的数值模拟研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *
李海燕 等: "南堡凹陷CO_2在咸水层中的矿化封存机理", 《东北石油大学学报》 *
莫绍星 等: "咸水层CO_2矿物封存数值模拟研究进展", 《地质科技情报》 *
莫绍星 等: "基于TOUGHREACT_MP的_省略_组咸水层CO_2矿物封存数值模拟", 《吉林大学学报(地球科学版)》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109828099A (en) * 2019-01-28 2019-05-31 中国石油大学(华东) A kind of rock-water-hydrocarbon analogue experiment method
CN110133223B (en) * 2019-05-20 2022-04-15 中国石油大学(华东) Dawsonite in CO2Preservation condition evaluation system in geological storage process
CN110133223A (en) * 2019-05-20 2019-08-16 中国石油大学(华东) A kind of dawsonite is in CO2Preservation condition evaluation system during geological reservoir
CN111855715A (en) * 2020-07-10 2020-10-30 中国石油天然气股份有限公司 Method and system for simulating and evaluating water-rock reaction of dolostone under salt
CN112924648B (en) * 2021-01-26 2022-01-04 大连理工大学 Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity
CN112924648A (en) * 2021-01-26 2021-06-08 大连理工大学 Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity
CN113610914A (en) * 2021-07-02 2021-11-05 中铁二十局集团有限公司 Sandstone pore structure analysis method, device, equipment and storage medium
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
WO2023071552A1 (en) * 2021-10-27 2023-05-04 中国华能集团清洁能源技术研究院有限公司 Gas solubility measuring apparatus and method, and storage potential predicting system and method
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
CN114781287A (en) * 2022-06-13 2022-07-22 武汉恒信楚才科技有限公司 Method and system for calculating natural gas-water-rock interaction under high-temperature and high-pressure conditions
CN115165700B (en) * 2022-06-17 2023-08-25 中国石油大学(北京) Experimental device for full-core-scale carbon dioxide adsorption sealing quantity
CN115165700A (en) * 2022-06-17 2022-10-11 中国石油大学(北京) Calculation method and experimental device for full-core-scale carbon dioxide adsorption and storage amount
CN115902160A (en) * 2022-11-28 2023-04-04 新疆敦华绿碳技术股份有限公司 Evaluation of CO in saline water reservoir geology 2 Method for sequestration of potential
CN115902160B (en) * 2022-11-28 2024-04-09 新疆敦华绿碳技术股份有限公司 Evaluation of CO in brine layer geology 2 Method for storing potential
CN116255198B (en) * 2023-02-13 2023-10-13 大连理工大学 Carbon dioxide sealing method based on reservoir wettability optimization design and layered regulation
CN116255198A (en) * 2023-02-13 2023-06-13 大连理工大学 Carbon dioxide sealing method based on reservoir wettability optimization design and layered regulation
CN117250133A (en) * 2023-11-16 2023-12-19 国建大数据科技(辽宁)有限公司 Smart city large-scale data acquisition method and system
CN117250133B (en) * 2023-11-16 2024-02-20 国建大数据科技(辽宁)有限公司 Smart city large-scale data acquisition method and system
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
CN117686689A (en) * 2024-02-02 2024-03-12 中国科学院地质与地球物理研究所 Basalt carbon sequestration efficiency measurement method based on basalt carbon sequestration potential evaluation
CN117686689B (en) * 2024-02-02 2024-05-07 中国科学院地质与地球物理研究所 Basalt carbon sequestration efficiency measurement method based on basalt carbon sequestration potential evaluation

Similar Documents

Publication Publication Date Title
CN108614076A (en) A kind of carbon dioxide geologic sequestration evaluation method
Mason et al. Petrology of the metamorphic rocks
CN108956435B (en) Simulation experiment method and device for high-temperature high-pressure reservoir corrosion
CN103196807B (en) The analytical approach of a kind of sandstone diagenetic process and pore evolution
Schaef et al. Basalt reactivity variability with reservoir depth in supercritical CO2 and aqueous phases
Huet et al. A review of Portland cement carbonation mechanisms in CO2 rich environment
CN112924648B (en) Evaluation geological sequestration CO2Method for mineralizing evolution law and sealing storage quantity
Mangenot et al. Coupling Δ47 and fluid inclusion thermometry on carbonate cements to precisely reconstruct the temperature, salinity and δ18O of paleo-groundwater in sedimentary basins
CN108333100B (en) Test device and method for simulating soil corrosion
Han et al. The curved 14C vs. δ13C relationship in dissolved inorganic carbon: A useful tool for groundwater age-and geochemical interpretations
Gholami et al. CO2 sequestration in sandstone reservoirs: How does reactive flow alter trapping mechanisms?
Liu et al. Microscopic & macroscopic characterizations of Beijing marble as a building material for UNESCO heritage sites: New insights into physico-mechanical property estimation and weathering resistance
Yanzhong et al. Geologic CO2 storage in arkosic sandstones with CaCl2-rich formation water
CN104406904B (en) A kind of method predicting the Concrete Sulfate Attack degree of depth
Holm et al. Groundwater geochemistry of short‐term aquifer thermal energy storage test cycles
Wells et al. Spatially-variable carbonation reactions in polycrystalline olivine
Gong et al. Dissolution of sandstone powders in deionised water over the range 50–350 C
Wolff-Boenisch et al. A foray into false positive results in mineral dissolution and precipitation studies
Tishchenko et al. The contribution of organic matter to the alkaline reserve of natural waters
Kirste et al. Geochemical modelling of CO2-water-rock interaction in the Pretty Hill Formation, Otway Basin
Holzheid Dissolution kinetics of selected natural minerals relevant to potential CO2-injection sites–Part 2: Dissolution and alteration of carbonates and feldspars in CO2-bearing brines
Onuk et al. Formation of helictite in the cave Dragon Belly (Sardinia, Italy)—Microstructure and incorporation of Mg, Sr, and Ba
Ngwenya et al. Reservoir sensitivity to water flooding: An experimental study of seawater injection in a North Sea reservoir analog
Hurd et al. On the estimation of minimum mechanical loss during an in situ biogenic silica dissolution experiment
CN104850742B (en) A kind of method that calculating CO2 salt water layer mineral seal potentiality up for safekeeping

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