CN115165700A - 全岩心尺度二氧化碳吸附封存量的计算方法及实验装置 - Google Patents
全岩心尺度二氧化碳吸附封存量的计算方法及实验装置 Download PDFInfo
- Publication number
- CN115165700A CN115165700A CN202210692457.8A CN202210692457A CN115165700A CN 115165700 A CN115165700 A CN 115165700A CN 202210692457 A CN202210692457 A CN 202210692457A CN 115165700 A CN115165700 A CN 115165700A
- Authority
- CN
- China
- Prior art keywords
- carbon dioxide
- core
- adsorption
- calculating
- pore throat
- 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
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 title claims abstract description 398
- 238000001179 sorption measurement Methods 0.000 title claims abstract description 212
- 229910002092 carbon dioxide Inorganic materials 0.000 title claims abstract description 199
- 239000001569 carbon dioxide Substances 0.000 title claims abstract description 199
- 238000004364 calculation method Methods 0.000 title claims abstract description 26
- 238000003860 storage Methods 0.000 title claims abstract description 13
- 239000011148 porous material Substances 0.000 claims abstract description 71
- 238000000034 method Methods 0.000 claims abstract description 54
- 239000011435 rock Substances 0.000 claims abstract description 47
- 238000007789 sealing Methods 0.000 claims abstract description 28
- 238000002474 experimental method Methods 0.000 claims abstract description 26
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims abstract description 19
- 238000012900 molecular simulation Methods 0.000 claims abstract description 16
- 229910052757 nitrogen Inorganic materials 0.000 claims abstract description 9
- 238000002336 sorption--desorption measurement Methods 0.000 claims abstract description 6
- 238000012856 packing Methods 0.000 claims description 8
- 239000002699 waste material Substances 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 6
- 239000002912 waste gas Substances 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000011065 in-situ storage Methods 0.000 claims description 3
- 230000009919 sequestration Effects 0.000 description 11
- 239000007789 gas Substances 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000000342 Monte Carlo simulation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 239000003245 coal Substances 0.000 description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 2
- 239000002156 adsorbate Substances 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000004587 chromatography analysis Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000012634 fragment Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000002411 thermogravimetry Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- ICXADQHBWHLSCI-UHFFFAOYSA-N dubinine Natural products C1=CC=C2C(OC)=C(CC(O3)C(C)(O)COC(C)=O)C3=NC2=C1 ICXADQHBWHLSCI-UHFFFAOYSA-N 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000009830 intercalation Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000329 molecular dynamics simulation Methods 0.000 description 1
- JCXJVPUVTGWSNB-UHFFFAOYSA-N nitrogen dioxide Inorganic materials O=[N]=O JCXJVPUVTGWSNB-UHFFFAOYSA-N 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
- G01N5/04—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by removing a component, e.g. by evaporation, and weighing the remainder
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02C—CAPTURE, STORAGE, SEQUESTRATION OR DISPOSAL OF GREENHOUSE GASES [GHG]
- Y02C20/00—Capture or disposal of greenhouse gases
- Y02C20/40—Capture or disposal of greenhouse gases of CO2
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pure & Applied Mathematics (AREA)
- Biochemistry (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Dispersion Chemistry (AREA)
- Geometry (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Computational Biology (AREA)
- Algebra (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Carbon And Carbon Compounds (AREA)
Abstract
本发明提供一种全岩心尺度二氧化碳吸附封存量的计算方法及实验装置，该方法包括通过氮气吸附/解吸实验获取岩心的孔喉半径分布；基于热重法二氧化碳吸附实验分别获取不同温度下的二氧化碳的过剩吸附量曲线；采用蒙特卡罗分子模拟方法对过剩吸附量曲线进行校正；根据孔喉半径分布和校正后的过剩吸附量曲线计算岩心的二氧化碳的绝对吸附封存量。本发明的计算方法在考虑了二氧化碳吸附相的体积后，并修正了孔喉尺寸的影响，因此，通过该方法计算得到的结果更加精确。
Description
技术领域
本发明属于油田二氧化碳地质埋存技术领域，尤其涉及一种全岩心尺度二氧化碳吸附封存量的计算方法及实验装置。
背景技术
随着二氧化碳排放量的日益增加，温室效应导致全球气候环境日趋恶劣。为了减少二氧化碳排放量，可以将二氧化碳封存至地层中，目前，二氧化碳封存原理主要有地质构造封存、残余气封存、溶解封存以及矿物封存四种机理。而吸附封存二氧化碳是物理形式的封存，是二氧化碳分子附在岩心内部孔喉表面的界面现象。而目前国内外学者对于二氧化碳封存量的计算，没有考虑孔喉吸附作用，并且前人对吸附作用的研究存在一定的误差，没有考虑吸附相体积的影响。因此，为了更准确的得到二氧化碳的封存量，吸附作用是不容忽视的。而前人对吸附相密度的研究尚不能反映岩石的真实吸附能力，这必将对客观评价岩石二氧化碳封存量产生影响。
发明内容
本发明的主要目的是提出一种全岩心尺度二氧化碳吸附封存量的计算方法及实验装置，旨在解决现有技术中的全岩心尺度二氧化碳吸附封存量量的计算不准确的技术问题。
为了实现上述目的，本发明提供一种全岩心尺度二氧化碳吸附封存量的计算方法，包括：
通过氮气吸附/解吸实验获取岩心的孔喉半径分布；
基于热重法二氧化碳吸附实验分别获取不同温度下的二氧化碳的过剩吸附量曲线；
采用蒙特卡罗分子模拟方法对所述过剩吸附量曲线进行校正；
根据所述孔喉半径分布和校正后的过剩吸附量曲线计算所述岩心的二氧化碳的绝对吸附封存量。
在本发明的实施例中，所述基于热重法二氧化碳吸附实验分别获取不同温度下的二氧化碳的过剩吸附量曲线的步骤包括：
在预设温度下，对所述岩心的密封腔内充入二氧化碳至不同的体系压力；
采用热重法进行二氧化碳吸附测试实验并获取每个体系压力下的二氧化碳的过剩吸附量；
根据不同体系压力下的二氧化碳的过剩吸附量数据绘制预设温度下的二氧化碳的过剩吸附量曲线；
改变预设温度并依次重复上述步骤，得到多条在不同预设温度下的二氧化碳的过剩吸附量曲线。
在本发明的实施例中，所述采用蒙特卡罗分子模拟方法对所述过剩吸附量曲线进行校正的步骤包括：
采用蒙特卡罗分子模拟方法计算二氧化碳在预设温度和不同体系压力下，不同孔喉半径所对应的二氧化碳吸附层的平均密度；
根据不同孔喉半径所对应的所述二氧化碳吸附层的平均密度校正所述岩心的过剩吸附量曲线。
在本发明的实施例中，所述采用蒙特卡罗分子模拟方法计算二氧化碳在预设温度和不同体系压力下，不同孔喉半径所对应的二氧化碳吸附层的平均密度的步骤包括：
计算在预设温度和不同体系压力下，孔喉半径为R1所对应的二氧化碳密度并绘制二氧化碳密度分布曲线；
根据二氧化碳分子的直径划分吸附层，得到孔喉半径为R1下的二氧化碳吸附层；
计算预设温度和不同体系压力下，孔喉半径为R1所对应的二氧化碳吸附层的平均密度；
改变选取的孔喉半径并依次重复上述步骤，并获取在预设温度和不同体系压力下，不同孔喉半径所对应的二氧化碳吸附层的平均密度曲线。
在本发明的实施例中，所述根据所述孔喉半径分布和校正后的过剩吸附量曲线计算所述岩心的二氧化碳的绝对吸附封存量的步骤包括：
根据岩心孔喉半径分布和不同孔喉半径所对应的二氧化碳吸附层的平均密度计算整个岩心的二氧化碳吸附层的平均密度；
根据过剩吸附量和整个岩心的二氧化碳吸附层的平均密度计算预设温度且不同体系压力状态下所述岩心的绝对吸附量。
在本发明的实施例中，岩心中单个孔喉半径下的二氧化碳吸附层的平均密度采用如下计算公式计算得到：
式中，ρave为单个二氧化碳吸附层的平均密度，kg/m3；ρads为二氧化碳吸附层的原位密度，kg/m3；a是二氧化碳吸附层的起点，b是二氧化碳吸附层的终点，zab是a与b之间的距离，为二氧化碳单个分子的直径长度，nm。
在本发明的实施例中，整个岩心的二氧化碳吸附层的平均密度采用如下计算公式得到：
式中，
为整个岩心的二氧化碳吸附层的平均密度，kg/m3；α1,α2,α3……αn分别是孔喉半径为R1、R2、R3、…Rn所占的频数；
分别是孔喉半径为R1、R2、R3、…Rn的吸附层的平均密度，kg/m3。
在本发明的实施例中，所述岩心的绝对吸附量采用如下计算公式能够得到：
式中，Mads为岩心的绝对吸附量，mmol/g；Mex为岩心的过剩吸附量，mmol/g；ρ为岩心的二氧化碳的体相密度，kg/m3；
为整个岩心的二氧化碳吸附层的平均密度，kg/m3。
在本发明的实施例中，还提出一种计算全岩心尺度二氧化碳吸附封存量的实验装置，应用于如上所述的全岩心尺度二氧化碳吸附封存量的计算方法中，所述计算全岩心尺度二氧化碳吸附封存量的实验装置包括依次连接的恒温箱、岩心粉碎机和岩心密封箱，所述岩心粉碎机用于将恒温箱内输送的岩心粉碎，所述岩心密封箱内形成有用于盛装所述岩心的收容腔，所述岩心密封箱选择连接有氮气瓶和二氧化碳瓶，所述岩心密封箱上设置有温度计和压力计。
在本发明的实施例中，所述计算全岩心尺度二氧化碳吸附封存量的实验装置还包括废弃瓶，所述废弃瓶和所述岩心密封箱管道连接并用于收集实验中产生的废气。
通过上述技术方案，本发明实施例所提供的全岩心尺度二氧化碳吸附封存量的计算方法具有如下的有益效果：
在计算全岩心尺度的二氧化碳吸附封存量时，首先通过氮气吸附/解吸实验获取岩心的孔喉半径分布；然后基于热重法二氧化碳吸附实验分别获取不同温度下的二氧化碳的过剩吸附量曲线；并采用蒙特卡罗分子模拟方法对过剩吸附量曲线进行校正；最后根据孔喉半径分布和校正后的过剩吸附量曲线计算岩心的二氧化碳的绝对吸附封存量；本发明的计算方法在考虑了二氧化碳吸附相的体积后，并修正了孔喉尺寸的影响，因此，通过该方法计算得到的结果更加精确。该方法可为技术人员提供准确的岩心尺度下的二氧化碳吸附封存量，进而可以计算油田尺度下二氧化碳封存量。
本发明的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明的理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明，但并不构成对本发明的限制。在附图中：
图1是根据本发明一实施例中全岩心尺度二氧化碳吸附封存量的计算方法的流程示意图；
图2是根据本发明一实施例中计算全岩心尺度二氧化碳吸附封存量的实验装置的结构示意图；
图3是根据本发明一实施例中岩心的孔喉半径分布图；
图4是根据本发明一实施例中在预设温度下二氧化碳的过剩吸附量曲线；
图5是根据本发明一实施例中在预设温度下，孔喉半径为R1所对应的二氧化碳密度分布图；
图6是根据本发明一实施例中在预设温度下，孔喉半径为R1所对应的二氧化碳吸附层的示意图；
图7是根据本发明一实施例中在预设温度下，孔喉半径为R1所对应的二氧化碳吸附层的平均密度分布图；
图8是根据本发明一实施例中在预设温度下，孔喉半径分别为R1、R2、R3和R4所对应的二氧化碳吸附层的平均密度分布图；
图9是根据本发明一实施例中在预设温度下，整个岩心的吸附层平均密度分布图；
图10是根据本发明一实施例中在预设温度下，整个岩心的绝对吸附量曲线图。
附图标记说明
具体实施方式
以下结合附图对本发明的具体实施例进行详细说明。应当理解的是，此处所描述的具体实施例仅用于说明和解释本发明，并不用于限制本发明。
下面参考附图描述根据本发明的全岩心尺度二氧化碳吸附封存量的计算方法及实验装置。
如图1所示，在本发明的实施例中，提供一种全岩心尺度二氧化碳吸附封存量的计算方法，包括：
步骤S10：通过氮气吸附/解吸实验获取岩心的孔喉半径分布；
首先选取一定尺寸的岩心样品进行实验，将直径为3.0cm、长6.0cm的岩心放入温度为T的恒温箱中恒温烘干24小时直至质量不发生变化，并通过粉碎设备将岩心机械打碎。将岩样碎块放入岩心密封箱的密封腔内，通过N2气吸附/解吸测试获得岩心的孔喉半径分布，孔喉半径分别为R1、R2、R3、…Rn，如图3所示，可以明显看出不同的孔喉半径所对应的频数，孔喉半径所对应的频数越高时，说明该孔喉半径在岩样中所占的比例越多，根据本实验所选取的岩样实验结果来看，本实施例中岩样的孔喉半径主要集中分布在14nm～28nm之间。
步骤S20：基于热重法二氧化碳吸附实验分别获取不同温度下的二氧化碳的过剩吸附量曲线；
步骤S30：采用蒙特卡罗分子模拟方法对过剩吸附量曲线进行校正；
步骤S40：根据孔喉半径分布和校正后的过剩吸附量曲线计算岩心的二氧化碳的绝对吸附封存量。
该方法首先通过N2吸附实验获得岩心孔喉半径分布，在此基础上基于热重法实验，获得二氧化碳在不同温度压力条件下的过剩吸附量曲线，最后，结合蒙特卡罗分子模拟方法，计算得到全岩心尺度下的二氧化碳吸附封存量，本发明的计算方法在考虑了二氧化碳吸附相的体积后，并修正了孔喉尺寸的影响，因此，通过该方法计算得到的结果更加精确。该方法可为技术人员提供准确的岩心尺度下的二氧化碳吸附封存量，进而可以计算油田尺度下二氧化碳封存量。
在本发明的实施例中，步骤S20还包括：
步骤S21：在预设温度下，对岩心的密封腔内充入二氧化碳至不同的体系压力；
步骤S22：采用热重法进行二氧化碳吸附测试实验并获取每个体系压力下的二氧化碳的过剩吸附量；
步骤S23：根据不同体系压力下的二氧化碳的过剩吸附量数据绘制预设温度下的二氧化碳的过剩吸附量曲线；
步骤S24：改变预设温度并依次重复上述步骤，得到多条在不同预设温度下的二氧化碳的过剩吸附量曲线。
在其中一个实施例中，在预设温度条件下的岩心密封腔内充入二氧化碳至体系压力为P1、P2、P3……Pn，在每个压力采用热重法进行二氧化碳吸附测试实验，得到预设温度，压力为P1、P2、P3……Pn下的二氧化碳过剩吸附量并绘制预设温度，不同压力下的二氧化碳的过剩吸附量曲线，如图4所示，本例中预设温度条件下，随着压力的增加，二氧化碳的过剩吸附量逐渐增加，当压力超过10MPa时，随着压力增加，二氧化碳的过剩吸附量逐渐降低。改变预设温度，重复上述步骤，便可获得在不同温度下，不同压力下的二氧化碳的过剩吸附量曲线。
在本发明的实施例中，步骤S30还包括：
步骤S31：采用蒙特卡罗分子模拟方法计算二氧化碳在预设温度和不同体系压力下，不同孔喉半径所对应的二氧化碳吸附层的平均密度；
步骤S32：根据不同孔喉半径所对应的二氧化碳吸附层的平均密度校正岩心的过剩吸附量曲线。
国内外学者对于吸附的研究主要在煤层气及页岩领域，吸附实验基本主要包括色谱分析法，常压流动法、容量法及热重法。色谱分析法实验时间短，但实验精度低，常压流动法的测试压力较低，具有一定的局限性，容量法主要针对于煤层气吸附实验，而热重法的实验精度更高，目前被广泛使用。国内外学者对煤层气和页岩等温吸附发现，等温吸附实验获得的吸附量为过剩吸附量，忽略了吸附相体积的影响，该实验结果无法反映真正的吸附能力，存在一定误差。目前，常用的方法是通过实验或者数值模拟的方法对过剩吸附量进行校正。Dubinin(1960)提出吸附相的密度为一个恒定值，与范德华常数b有关。后来，人Tsai(1985)认为吸附相的密度等于液体吸附质的密度。Li et al.(2002)对上述方法进行了比较，认为吸附相的密度是系统温度的函数。Ambrose(2012)通过分子模拟表明，吸附相的密度与体系温度、压力和孔径有关。
近些年来，分子模拟方法随着计算机技术的快速发展被广泛运用到气体吸附特性的研究中。该方法以统计力学的研究为基础，结合分子力场研究理论，在分子、原子水平上进一步进行微观研究，建立分子模型来模拟不同条件下其相关结构与参数，进而研究其物理化学性质。分子模拟方法根据研究体系的微观状态量分不同，可以将其划分为两种方法：蒙特卡罗方法和分子动力学方法。本发明基于蒙特卡罗模拟方法对二氧化碳吸附封存进行计算。蒙特卡洛方法的原理是基于通过随机抽取样本的方式来分析研究目标。它的基本步骤是先构建一个相应的概率模型去代替实际系统，然后在此概率模型的基础上随机抽样模拟，就能计算得到实际系统的一些有关统计方面的近似解。其中，蒙特卡罗分子模拟方法的详细步骤如下：
(1)通过随机数产生器产生一个相应的分子构象，然后将要研究的吸附质分子以任意方向随机加入该构象中的任意位置。
(2)改变此分子构象内的相关粒子发生相应的位置变化，从而产生另一种新构象。每次变化均应包括三种扰动：插入分子；删除分子；移动分子。
(3)计算分析获取新构象与原构象之间的能量改变量ΔE，根据ΔE的值确定新构象是否加入到统计平均中。
(4)若不能加入到统计平均中，则放弃此构型，重新进行计算，直至能够选取得到相应最稳定的分子构象。
(5)当体系达到统计平均后，当前体系计算的吸附量为最终的二氧化碳吸附量。
进一步地，步骤S31包括：
步骤S311：计算在预设温度和不同体系压力下，孔喉半径为R1所对应的二氧化碳密度并绘制二氧化碳密度分布曲线；
步骤S312：根据二氧化碳分子的直径划分吸附层，得到孔喉半径为R1下的二氧化碳吸附层；
如图5和图6所示，为在预设温度下，当孔喉半径为R1下的二氧化碳吸附层的密度分布和吸附层所在的示意图；本实施例以二氧化碳分子的直径划分为一个吸附层，从图6可以看出，当孔喉半径为为R1时，存在两个二氧化碳吸附层：a至b之间、b’至a’之间，每个二氧化碳吸附层的密度也可以从图上可以直接得到。
步骤S313：计算预设温度和不同体系压力下，孔喉半径为R1所对应的二氧化碳吸附层的平均密度；对于如图6所示的实施例中，该孔喉半径下存在多个二氧化碳吸附层时，每个吸附层的不同位置所对应的密度也不同，因此需要计算该孔喉半径下的每个二氧化碳吸附层的平均密度；具体地，岩心中单个孔喉半径下的二氧化碳吸附层的平均密度采用如下计算公式计算得到：
式中，ρave为单个二氧化碳吸附层的平均密度，kg/m3；ρads为二氧化碳吸附层的原位密度，kg/m3；a是二氧化碳吸附层的起点，b是二氧化碳吸附层的终点，zab是a与b之间的距离，为单个二氧化碳分子的直径长度，nm。
步骤S314：改变选取的孔喉半径并依次重复上述步骤，并获取在预设温度和不同体系压力下，不同孔喉半径所对应的二氧化碳吸附层的平均密度曲线。
如图8所示，选取了孔喉半径为R1、R2、R3以及R4所对应的二氧化碳吸附层的平均密度曲线图，其中，R1＞R2＞R3＞R4，对于不同的孔喉半径来说，相应的二氧化碳吸附层的平均密度也不同，并且，从图8上来看，对于不同的孔喉半径下二氧化碳吸附层的平均密度的增长曲线趋势一致，在相同的体系压力和温度下，孔喉半径越大，所对应的二氧化碳吸附层的平均密度越小。
在本发明的实施例中，步骤S40还包括：
步骤S41：根据岩心孔喉半径分布和不同孔喉半径所对应的二氧化碳吸附层的平均密度计算整个岩心的二氧化碳吸附层的平均密度；其中，整个岩心的二氧化碳吸附层的平均密度采用如下计算公式得到：
式中，
为整个岩心的二氧化碳吸附层的平均密度，kg/m3；α1,α2,α3……αn分别是孔喉半径为R1、R2、R3、…Rn所占的频数；
分别是孔喉半径为R1、R2、R3、…Rn的吸附层的平均密度，kg/m3。
如图9所示，即为本实施例在预设温度下，随着体系压力的增大，整个岩心的二氧化碳吸附层的平均密度曲线图，在该实施例中，选取了0～30MPa压力范围内的岩心的二氧化碳吸附层的平均密度数据，并绘制成曲线；从图9上来看，随着体系压力的增大，岩心内二氧化碳的吸附量随着增大，整个岩心二氧化碳吸附层的平均密度成逐渐增大的趋势，当压力大于15MPa时，整个岩心的二氧化碳吸附层的平均密度增大速度降低；并且可以看出，随着体系压力的增大，整个岩心的二氧化碳吸附层的平均密度分布曲线和岩心中单个孔喉半径下的二氧化碳吸附层的平均密度分布曲线趋势接近。由于对于不同的孔喉半径所在的吸附层来说，二氧化碳的吸附情况也不同，而本实施例在计算整个岩心的二氧化碳吸附层的平均密度时考虑了孔喉半径的分布情况，使得对整个岩心的二氧化碳吸附层的平均密度的计算更加准确。
步骤S42：根据过剩吸附量和整个岩心的二氧化碳吸附层的平均密度计算预设温度且不同体系压力状态下岩心的绝对吸附量；而该岩心的绝对吸附量即为全岩心尺度二氧化碳的吸附封存量。其中，岩心的绝对吸附量采用如下计算公式能够得到：
式中，Mads为岩心的绝对吸附量，mmol/g；Mex为岩心的过剩吸附量，mmol/g；ρ为岩心的二氧化碳的体相密度，kg/m3，可从美国国家标准与技术研究院化学数据库(NIST)中直接查询获取；
为整个岩心内二氧化碳吸附层的平均密度，kg/m3。绝对吸附量表示的页岩中二氧化碳的实际吸附量，当压力增加到一定程度，吸附必然会达到饱和，表现为绝对吸附量的不再增加。具体地的实验结果如图10所示，随着横坐标体系压力的增加，该实施例纵坐标中的岩心的绝对吸附量增大，但当压力到达一定值时，岩心的绝对吸附量增加缓慢，当体系压力达到23～25MPa左右时，岩心的绝对吸附量不再变化，说明此时二氧化碳的吸附已经达到饱和。本发明通过考虑岩心的过剩吸附量和二氧化碳吸附层的密度，也就是说考虑了二氧化碳吸附相的体积，修正了孔喉半径的影响，使得岩心的绝对吸附量的计算结果更加精确。
此外，在本发明的实施例中，如图2所示，还提出一种计算全岩心尺度二氧化碳吸附封存量的实验装置，应用于如上的全岩心尺度二氧化碳吸附封存量的计算方法中，计算全岩心尺度二氧化碳吸附封存量的实验装置包括依次连接的恒温箱1、岩心粉碎机2和岩心密封箱3，岩心粉碎机2用于将恒温箱1内输送的岩心粉碎，岩心密封箱3内形成有用于盛装岩心的收容腔，岩心密封箱3选择连接有氮气瓶8和二氧化碳瓶9，岩心密封箱3上设置有用于检测收容腔内温度的温度计和检测压力的压力计。其中，岩心粉碎机2为机械粉碎的结构原理形式，可以采用现有技术中常规的岩心粉碎机2的结构形式进行，但并不仅限于本发明的实施例。岩心密封箱3的下端通过连通管道上的支管道分别连接氮气瓶8和二氧化碳瓶9，在连通管道的主管道和两个支管道上均设有阀门6，以选择氮气或二氧化碳通入岩心密封箱3内。在连通管道的主管道上设有检测气体流量或流速的第一流量计4，通过采用第一流量计4实时监控对岩心密封箱3内充入气体的流量，以达到调节整个岩心密封箱3内的体系压力的目的。
在本发明的实施例中，该测定CO2在孔隙中不同封存状态封存量的实验装置还包括废弃瓶5，废弃瓶5和岩心密封箱3管道连接并用于收集实验中产生的废气，且在废弃瓶5和岩心密封箱3之间的连接管道上设置有检测废气流量的第二流量计7。
在本发明的描述中，需要理解的是，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
在本发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接或彼此可通讯；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系，除非另有明确的限定。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (10)
1.一种全岩心尺度二氧化碳吸附封存量的计算方法，其特征在于，包括：
通过氮气吸附/解吸实验获取岩心的孔喉半径分布；
基于热重法二氧化碳吸附实验分别获取不同温度下的二氧化碳的过剩吸附量曲线；
采用蒙特卡罗分子模拟方法对所述过剩吸附量曲线进行校正；
根据所述孔喉半径分布和校正后的过剩吸附量曲线计算所述岩心的二氧化碳的绝对吸附封存量。
2.根据权利要求1所述的全岩心尺度二氧化碳吸附封存量的计算方法，其特征在于，所述基于热重法二氧化碳吸附实验分别获取不同温度下的二氧化碳的过剩吸附量曲线的步骤包括：
在预设温度下，对所述岩心的密封腔内充入二氧化碳至不同的体系压力；
采用热重法进行二氧化碳吸附测试实验并获取每个体系压力下的二氧化碳的过剩吸附量；
根据不同体系压力下的二氧化碳的过剩吸附量数据绘制预设温度下的二氧化碳的过剩吸附量曲线；
改变预设温度并依次重复上述步骤，得到多条在不同预设温度下的二氧化碳的过剩吸附量曲线。
3.根据权利要求1所述的全岩心尺度二氧化碳吸附封存量的计算方法，其特征在于，所述采用蒙特卡罗分子模拟方法对所述过剩吸附量曲线进行校正的步骤包括：
采用蒙特卡罗分子模拟方法计算二氧化碳在预设温度和不同体系压力下，不同孔喉半径所对应的二氧化碳吸附层的平均密度；
根据不同孔喉半径所对应的所述二氧化碳吸附层的平均密度校正所述岩心的过剩吸附量曲线。
4.根据权利要求3所述的全岩心尺度二氧化碳吸附封存量的计算方法，其特征在于，所述采用蒙特卡罗分子模拟方法计算二氧化碳在预设温度和不同体系压力下，不同孔喉半径所对应的二氧化碳吸附层的平均密度的步骤包括：
计算在预设温度和不同体系压力下，孔喉半径为R1所对应的二氧化碳密度并绘制二氧化碳密度分布曲线；
根据二氧化碳分子的直径划分吸附层，得到孔喉半径为R1下的二氧化碳吸附层；
计算预设温度和不同体系压力下，孔喉半径为R1所对应的二氧化碳吸附层的平均密度；
改变选取的孔喉半径并依次重复上述步骤，并获取在预设温度和不同体系压力下，不同孔喉半径所对应的二氧化碳吸附层的平均密度曲线。
5.根据权利要求4所述的全岩心尺度二氧化碳吸附封存量的计算方法，其特征在于，所述根据所述孔喉半径分布和校正后的过剩吸附量曲线计算所述岩心的二氧化碳的绝对吸附封存量的步骤包括：
根据岩心孔喉半径分布和不同孔喉半径所对应的二氧化碳吸附层的平均密度计算整个岩心的二氧化碳吸附层的平均密度；
根据过剩吸附量和整个岩心的二氧化碳吸附层的平均密度计算预设温度且不同体系压力状态下所述岩心的绝对吸附量。
6.根据权利要求5所述的全岩心尺度二氧化碳吸附封存量的计算方法，其特征在于，岩心中单个孔喉半径下的二氧化碳吸附层的平均密度采用如下计算公式计算得到：
式中，ρave为单个二氧化碳吸附层的平均密度，kg/m3；ρads为二氧化碳吸附层的原位密度，kg/m3；a是二氧化碳吸附层的起点，b是二氧化碳吸附层的终点，zab是a与b之间的距离，为二氧化碳单个分子的直径长度，nm。
7.根据权利要求6所述的全岩心尺度二氧化碳吸附封存量的计算方法，其特征在于，整个岩心的二氧化碳吸附层的平均密度采用如下计算公式得到：
式中，
为整个岩心的二氧化碳吸附层的平均密度，kg/m3；α1,α2,α3……αn分别是孔喉半径为R1、R2、R3、…Rn所占的频数；
分别是孔喉半径为R1、R2、R3、…Rn的吸附层的平均密度，kg/m3。
8.根据权利要求7所述的全岩心尺度二氧化碳吸附封存量的计算方法，其特征在于，所述岩心的绝对吸附量采用如下计算公式能够得到：
式中，Mads为岩心的绝对吸附量，mmol/g；Mex为岩心的过剩吸附量，mmol/g；ρ为岩心的二氧化碳的体相密度，kg/m3；
为整个岩心的二氧化碳吸附层的平均密度，kg/m3。
9.一种计算全岩心尺度二氧化碳吸附封存量的实验装置，其特征在于，应用于如权利要求1至8中任意一项所述的全岩心尺度二氧化碳吸附封存量的计算方法中，所述计算全岩心尺度二氧化碳吸附封存量的实验装置包括依次连接的(1)、岩心粉碎机(2)和岩心密封箱(3)，所述岩心粉碎机(2)用于将恒温箱内输送的岩心粉碎，所述岩心密封箱(3)内形成有用于盛装所述岩心的收容腔，所述岩心密封箱(3)选择连接有氮气瓶(8)和二氧化碳瓶(9)，所述岩心密封箱(3)上设置有温度计和压力计。
10.根据权利要求9所述的计算全岩心尺度二氧化碳吸附封存量的实验装置，其特征在于，所述计算全岩心尺度二氧化碳吸附封存量的实验装置还包括废弃瓶(5)，所述废弃瓶(5)和所述岩心密封箱(3)管道连接并用于收集实验中产生的废气。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210692457.8A CN115165700B (zh) | 2022-06-17 | 2022-06-17 | 全岩心尺度二氧化碳吸附封存量的实验装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210692457.8A CN115165700B (zh) | 2022-06-17 | 2022-06-17 | 全岩心尺度二氧化碳吸附封存量的实验装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115165700A true CN115165700A (zh) | 2022-10-11 |
| CN115165700B CN115165700B (zh) | 2023-08-25 |
Family
ID=83484433
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210692457.8A Active CN115165700B (zh) | 2022-06-17 | 2022-06-17 | 全岩心尺度二氧化碳吸附封存量的实验装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115165700B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104568702A (zh) * | 2015-01-27 | 2015-04-29 | 中国石油大学(北京) | 聚合物的不可及孔隙体积和不可入孔隙半径的测定方法 |
| CN108181200A (zh) * | 2017-12-27 | 2018-06-19 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 泥页岩等温吸附实验曲线校正方法及装置 |
| CN108614076A (zh) * | 2018-04-27 | 2018-10-02 | 中国石油大学(华东) | 一种二氧化碳地质封存评价方法 |
| CN111189758A (zh) * | 2020-01-13 | 2020-05-22 | 闽江学院 | 一种致密砂岩或页岩油储层孔喉全尺度表征的方法 |
| CN113670960A (zh) * | 2021-07-23 | 2021-11-19 | 西南石油大学 | 一种基于分子模拟的真实页岩气吸附量预测方法 |
| CN114112856A (zh) * | 2021-12-07 | 2022-03-01 | 核工业北京化工冶金研究院 | 一种含铀砂岩全尺寸孔径结构表征方法及系统 |
-
2022
- 2022-06-17 CN CN202210692457.8A patent/CN115165700B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104568702A (zh) * | 2015-01-27 | 2015-04-29 | 中国石油大学(北京) | 聚合物的不可及孔隙体积和不可入孔隙半径的测定方法 |
| CN108181200A (zh) * | 2017-12-27 | 2018-06-19 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 泥页岩等温吸附实验曲线校正方法及装置 |
| CN108614076A (zh) * | 2018-04-27 | 2018-10-02 | 中国石油大学(华东) | 一种二氧化碳地质封存评价方法 |
| CN111189758A (zh) * | 2020-01-13 | 2020-05-22 | 闽江学院 | 一种致密砂岩或页岩油储层孔喉全尺度表征的方法 |
| CN113670960A (zh) * | 2021-07-23 | 2021-11-19 | 西南石油大学 | 一种基于分子模拟的真实页岩气吸附量预测方法 |
| CN114112856A (zh) * | 2021-12-07 | 2022-03-01 | 核工业北京化工冶金研究院 | 一种含铀砂岩全尺寸孔径结构表征方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| YUELIANG LIUA等: "Determination of the absolute adsorption/desorption isotherms of CH4 and n-C4H10 on shale from a nano-scale perspective", 《FUEL》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115165700B (zh) | 2023-08-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107807143B (zh) | 水合物专用低场核磁共振多探头定量测试系统及方法 | |
| CN102252918B (zh) | 含天然气水合物沉积物三轴试验装置及其试验方法 | |
| Lu et al. | Adsorption measurements in Devonian shales | |
| Liu et al. | CO2 sequestration in depleted methane hydrate sandy reservoirs | |
| CN101762435B (zh) | 宽温区储氢材料低温储氢性能试验装置 | |
| Zhu et al. | CO 2 capture with activated carbons prepared by petroleum coke and KOH at low pressure | |
| CN112710554B (zh) | 含水合物沉积物三轴剪切、渗流一体化实验方法 | |
| CN103226089B (zh) | 一种页岩气体渗透率测定方法 | |
| CN110534161B (zh) | 吸附质气体吸附相密度模型构建及绝对吸附量计算方法 | |
| CN110345904B (zh) | 水合物分解过程中沉积物变形及渗透率测试装置及方法 | |
| CN115144298A (zh) | 测定co2在孔隙中不同封存状态封存量的实验方法及装置 | |
| Belandria et al. | Compositional analysis of the gas phase for the CO2+ N2+ tetra-n-butylammonium bromide aqueous solution systems under hydrate stability conditions | |
| WO2011133885A1 (en) | Total storage capacity and total porosity of porous media | |
| CN106872328A (zh) | 一种低渗透岩心孔隙度和渗透率的测试装置及测试方法 | |
| US20160290907A1 (en) | Determination of free volume of a rock sample using high pressure adsorption data | |
| WO2022161137A1 (zh) | 岩石动态物性测量系统 | |
| CN109655372A (zh) | 一种PdO吸气剂吸气量测试装置及方法 | |
| Lee et al. | Novel Sieverts’ type volumetric measurements of hydrogen storage properties for very small sample quantities | |
| Li et al. | Effects of particle size and adsorption pressure on methane gas desorption and diffusion in coal | |
| Fan et al. | A novel experimental system for accurate gas sorption and its application to various shale rocks | |
| Feng et al. | Effects of temperature and pressure on gas desorption in coal in an enclosed system: a theoretical and experimental study | |
| CA3174170A1 (en) | Hydrogen production, storage and recovery | |
| CN115165700B (zh) | 全岩心尺度二氧化碳吸附封存量的实验装置 | |
| CN111664825B (zh) | 一种提高孔结构参数测定准确性的方法 | |
| CN105474011A (zh) | 确定油组合物与流体的最小混溶压力 |
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 |