CN113670960B - 一种基于分子模拟的真实页岩气吸附量预测方法 - Google Patents
一种基于分子模拟的真实页岩气吸附量预测方法 Download PDFInfo
- Publication number
- CN113670960B CN113670960B CN202110837284.XA CN202110837284A CN113670960B CN 113670960 B CN113670960 B CN 113670960B CN 202110837284 A CN202110837284 A CN 202110837284A CN 113670960 B CN113670960 B CN 113670960B
- Authority
- CN
- China
- Prior art keywords
- shale
- gas adsorption
- molecular
- shale gas
- adsorption capacity
- 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.)
- Active
Links
- 238000001179 sorption measurement Methods 0.000 title claims abstract description 180
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000012900 molecular simulation Methods 0.000 title claims abstract description 45
- 239000002734 clay mineral Substances 0.000 claims abstract description 74
- 238000002474 experimental method Methods 0.000 claims abstract description 59
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 claims abstract description 53
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 33
- 239000001569 carbon dioxide Substances 0.000 claims abstract description 26
- 229910002092 carbon dioxide Inorganic materials 0.000 claims abstract description 26
- 125000000524 functional group Chemical group 0.000 claims abstract description 18
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 18
- 239000011707 mineral Substances 0.000 claims abstract description 18
- 239000005416 organic matter Substances 0.000 claims abstract description 12
- 239000000203 mixture Substances 0.000 claims abstract description 9
- 239000007789 gas Substances 0.000 claims description 159
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 128
- 238000004088 simulation Methods 0.000 claims description 36
- 239000013078 crystal Substances 0.000 claims description 35
- 238000002441 X-ray diffraction Methods 0.000 claims description 26
- 239000011435 rock Substances 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 16
- 230000003993 interaction Effects 0.000 claims description 16
- 238000005457 optimization Methods 0.000 claims description 15
- 238000001035 drying Methods 0.000 claims description 12
- 238000005303 weighing Methods 0.000 claims description 12
- 238000004510 Lennard-Jones potential Methods 0.000 claims description 7
- 229910052799 carbon Inorganic materials 0.000 claims description 7
- 238000000227 grinding Methods 0.000 claims description 7
- 238000011049 filling Methods 0.000 claims description 6
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 239000001307 helium Substances 0.000 claims description 4
- 229910052734 helium Inorganic materials 0.000 claims description 4
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 claims description 4
- 238000011068 loading method Methods 0.000 claims description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 238000010812 external standard method Methods 0.000 claims description 3
- 238000000302 molecular modelling Methods 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 210000004027 cell Anatomy 0.000 claims 6
- 238000000342 Monte Carlo simulation Methods 0.000 claims 2
- 229910052573 porcelain Inorganic materials 0.000 claims 2
- 238000002447 crystallographic data Methods 0.000 claims 1
- 239000000126 substance Substances 0.000 abstract description 6
- 125000004429 atom Chemical group 0.000 description 17
- 239000010410 layer Substances 0.000 description 11
- 238000001228 spectrum Methods 0.000 description 9
- 239000000463 material Substances 0.000 description 8
- 230000009881 electrostatic interaction Effects 0.000 description 5
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 239000000919 ceramic Substances 0.000 description 4
- 229910001919 chlorite Inorganic materials 0.000 description 4
- 229910052619 chlorite group Inorganic materials 0.000 description 4
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 4
- 230000002860 competitive effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- GUJOJGAPFQRJSV-UHFFFAOYSA-N dialuminum;dioxosilane;oxygen(2-);hydrate Chemical group O.[O-2].[O-2].[O-2].[Al+3].[Al+3].O=[Si]=O.O=[Si]=O.O=[Si]=O.O=[Si]=O GUJOJGAPFQRJSV-UHFFFAOYSA-N 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 238000005984 hydrogenation reaction Methods 0.000 description 4
- NLYAJNPCOHFWQQ-UHFFFAOYSA-N kaolin Chemical group O.O.O=[Al]O[Si](=O)O[Si](=O)O[Al]=O NLYAJNPCOHFWQQ-UHFFFAOYSA-N 0.000 description 4
- VGIBGUSAECPPNB-UHFFFAOYSA-L nonaaluminum;magnesium;tripotassium;1,3-dioxido-2,4,5-trioxa-1,3-disilabicyclo[1.1.1]pentane;iron(2+);oxygen(2-);fluoride;hydroxide Chemical group [OH-].[O-2].[O-2].[O-2].[O-2].[O-2].[F-].[Mg+2].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[K+].[K+].[K+].[Fe+2].O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2 VGIBGUSAECPPNB-UHFFFAOYSA-L 0.000 description 4
- 125000004430 oxygen atom Chemical group O* 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- UPMLOUAZCHDJJD-UHFFFAOYSA-N 4,4'-Diphenylmethane Diisocyanate Chemical compound C1=CC(N=C=O)=CC=C1CC1=CC=C(N=C=O)C=C1 UPMLOUAZCHDJJD-UHFFFAOYSA-N 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 229910052900 illite Inorganic materials 0.000 description 3
- 229910052622 kaolinite Inorganic materials 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 229910052901 montmorillonite Inorganic materials 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 3
- 150000001721 carbon Chemical group 0.000 description 2
- 125000004432 carbon atom Chemical group C* 0.000 description 2
- 239000012153 distilled water Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011067 equilibration Methods 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 239000000843 powder Substances 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000000087 stabilizing effect Effects 0.000 description 2
- 239000000725 suspension Substances 0.000 description 2
- 238000001132 ultrasonic dispersion Methods 0.000 description 2
- 238000005406 washing Methods 0.000 description 2
- 150000001335 aliphatic alkanes Chemical class 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 229910002090 carbon oxide Inorganic materials 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000036571 hydration Effects 0.000 description 1
- 238000006703 hydration reaction Methods 0.000 description 1
- 125000002887 hydroxy group Chemical group [H]O* 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004898 kneading Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000004570 mortar (masonry) Substances 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000007873 sieving Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
- G01N23/20008—Constructional details of analysers, e.g. characterised by X-ray source, detector or optical system; Accessories therefor; Preparing specimens therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
- G01N23/2055—Analysing diffraction patterns
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
-
- 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/02—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by absorbing or adsorbing components of a material and determining change of weight of the adsorbent, e.g. determining moisture content
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Engineering & Computer Science (AREA)
- Crystallography & Structural Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于分子模拟的真实页岩气吸附量预测方法,包括通过XRD衍射实验获取页岩中矿物成分信息;通过TOC实验获取页岩中有机物组分信息;进行页岩气吸附实验获取吸附实验数据;利用获取的真实页岩成分信息进行多种物质混合分子模型构建,同时建立气体分子模型;模拟计算不同压力、温度下的真实页岩气吸附量;对比计算出的页岩气吸附量与实验测试的真实页岩气吸附量之间的误差并计算更高压力下页岩气吸附量;模拟计算不同含水量及不同二氧化碳含量对于页岩气吸附量的影响。本发明全面获取了页岩中粘土矿物的物相种类、晶胞结构信息以及干酪根官能团种类及TOC含量占比信息,使得建立的页岩分子模拟更符合页岩地层真实情况。
Description
技术领域
本发明涉及油气田开发领域,具体涉及一种基于分子模拟的真实页岩气吸附量预测方法。
背景技术
随着能源结构的调整,页岩气作为清洁能源发挥着越来越大的作用。然而页岩作为页岩气的存储基质,结构成分复杂,赋存方式多样,导致人们对页岩气微观吸附机理认识不清。因此,需要对页岩气赋存方式及吸附量进行研究。目前,有两种方法用于研究页岩基质中的页岩气(主要为CH4)吸附,分别为实验和模拟方法。通过实验方法获得的结果是高度可靠的,但只能获得吸附量,并不能展现出微观吸附机理,同时需要耗费大量测试时间,而且目前还没有高压吸附量的测试装置;分子模拟的方法是获得微观吸附机理的良好方法,但是页岩中含有多种矿物成分,许多学者认为难以表征矿物质中每种成分,因而大多数研究使用单一成分(例如干酪根或一种粘土矿物)来代表页岩表面,这些表面在矿物的异质性和多样性方面与实际的页岩表面没有相似之处。在使用一种以上组分的情况下,它们没有混合而是在纳米孔的不同末端,所建立的模型与真实页岩的差异性导致模拟吸附量与实验吸附量具有较大差异。因此,急需一种全新的真实页岩气吸附量预测方法,通过建立可表征真实页岩成分、含量和结构的分子模型,实现真实页岩气吸附量的直接准确预测,并且可用于目前实验无法达到的高温高压下的页岩气吸附量,为国内外各大页岩气田储量预测及后续勘探与开发提供依据。
目前涉及分子模拟计算页岩气吸附量的方法存在的问题有:
(1) 无法运用分子模型有效表征真实页岩多种矿物成分、含量和结构;
(2) 所建分子模型计算得出的吸附量与实验值差距过大;
(3) 没有一种可靠的方式证明分子模型的准确性。
发明内容
本发明的目的在于提供一种基于分子模拟的真实页岩气吸附量预测方法,以解决现有技术中无法运用分子模型有效表征真实页岩多种矿物成分、含量和结构,以及所建分子模型计算得出的吸附量与实验值差距过大的技术问题。
为实现上述目的,本发明提供了以下技术方案:
本发明提供的一种基于分子模拟的真实页岩气吸附量预测方法,包括以下步骤:
S1、通过XRD衍射实验获取页岩中矿物成分信息,所述矿物成分信息包括密度值;
S2、通过TOC实验获取页岩中有机物组分信息,所述页岩中有机物组分信息包括实验所得页岩气吸附量;
S3、进行页岩气吸附实验获取吸附实验数据,所述吸附实验数据包括不同压力点下的页岩气吸附量数据,用于与步骤S6所得模拟值进行对比;
S4、根据S1与S2中获取的真实页岩成分信息进行真实页岩分子模型构建,同时建立甲烷气体分子模型;
S5、优化S4中建立的真实页岩分子模型及甲烷气体分子模型;
S6、利用步骤S5中优化后的真实页岩分子模型和甲烷气体分子模型,计算不同压力、温度下的真实页岩气吸附量;
S7、对比6计算出的页岩气吸附量与3实验测试的真实页岩气吸附量,计算误差;
S8、计算比S6中设定的压力更高的压力下页岩气吸附量。
S9、利用步骤S5中优化后的真实页岩分子模型和甲烷气体分子模型,模拟计算不同含水量及不同二氧化碳含量对于页岩气吸附量的影响。
进一步的,所述步骤S1中,通过XRD衍射实验获取页岩中矿物成分信息,包含以下子步骤:
S101、对页岩样品进行粉末化处理;
S102、提取岩石样品中的粘土矿物并将其制片;
S103、对制片的样品进行X射线衍射,得到XRD图谱;
S104、运用MDI jade软件提取XRD图谱中的粘土矿物晶体信息,其中包括页岩中各矿物组分占比、粘土矿物各组分占比以及粘土矿物各组分晶胞参数。
进一步的,所述步骤S2中,通过TOC实验获取页岩中有机物组分信息,包含以下子步骤:
S201、对页岩样品进行粉末化处理;
S202、称取1g试样置于瓷坩埚中;
S203、反复向盛有试样的瓷坩埚内少量、多次地缓慢加入HCl,放入60~80℃的烘箱内烘干,直至反应完全为止;
S204、在100~120℃烘24h,直至除尽过量的HCl;
S205、通入温度为800℃以上的氧气流使页岩样品燃烧,使其中总有机碳全部转化成二氧化碳,用热导检测器(或红外检测器)检测二氧化碳的量,用外标法计算试样的总有机碳含量。
进一步的,所述步骤S3中,进行页岩气吸附实验获取吸附实验数据,包含以下子步骤:
S301、将样品装入样品池,打开温度控制系统和抽真空装置,使样品池温度稳定在实验预设温度;
S302、向样品池中充入高纯氮气,采集压力数据,确保系统气密性良好;
S303、调节并设置称量系统,逐点降低样品池压力,称量对应压力点下的样品桶重量;
S304、装入实验样品并抽真空,真空状态稳定后关闭抽真空装置,向系统通入氦气;
S305、调节并设置称量系统,逐点降低样品池压力,称量对应压力点下样品桶与样品总重量;
S306、对系统进行抽真空,真空状态稳定后关闭抽真空装置,向系统内通入甲烷,控制压力为第一个实验点;
S307、逐渐增大实验压力,直到最后一个压力点实验结束;
S308、记录每个压力点页岩气吸附量数据。
进一步的,所述步骤S4中,进行分子模拟建模包括以下子步骤:
S401、构建岩石中粘土矿物晶胞:
利用步骤S1中所得到的粘土矿物晶体信息,建立晶胞,设置晶胞信息,分别为晶胞空间群、基本晶胞参数,并在每种粘土矿物晶胞中输入每种原子位置信息;
S402、划分岩石中粘土矿物晶面:
通过晶面划分模块,利用步骤S1中得到的页岩中每种粘土矿物的晶面指数与面间距d值,确定出岩石样品中每种粘土矿物的层间距与分割面;
S403、建立岩石样品中每种粘土矿物混层结构:
通过步骤S1中得出的页岩中粘土矿物每种物相含量,并结合每种物相测量的峰强度比例,获得每种粘土矿物所占比例,进而得出分子模拟中每种粘土矿物最小层间占比与分布情况;并通过Materials Studio中建立晶体层模块,对每种划分晶面后的粘土矿物进行空间布局。
S404、建立页岩中有机物官能团,并随机插入所建立的粘土矿物混层结构中:
通过步骤S2中得出的TOC数据,建立页岩中干酪根官能团,并通过TOC含量占比换算为干酪根官能团含量占比,将不同种类的官能团随机插入所建立的粘土矿物混层结构中。
S405、建立甲烷气体分子模型:
(1)通过Materials Studio建立碳原子,然后加氢的方式建立出甲烷分子的初始模型。
(2)通过Materials Studio建立氧原子,然后加碳原子方式建立出二氧化碳分子的初始模型。
(3)通过Materials Studio建立氧原子,然后加氢的方式建立出水蒸气分子的初始模型。
进一步的,所述步骤S5中,优化真实页岩分子模型与甲烷气体分子模型包括以下子步骤:
S501、针对步骤S4中所构建的每种粘土矿物晶胞,通过Materials Studio分子模拟软件得到所建粘土矿物晶胞密度信息,并与步骤S1中得到的每种粘土矿物晶胞信息的密度值对比,通过微调晶体中原子位置、键长及增减原子个数的方式,使得通过所建模型得到的密度值与步骤S1中所测密度值相符。
S502、针对步骤S4中所构建的岩石样品中各粘土矿物最终的晶胞模型,通过MS分子模拟软件得到所建模型总的晶胞密度信息,并与步骤S1中得到的页岩中粘土矿物晶胞信息的密度值对比,通过调整每层粘土矿物层间距的方式,使得通过所建模型得到的密度值与步骤S1中所测密度值相符。
S503、运用MS中Forcite模块内的Geometry Optimization任务对步骤S4中真实页岩分子模型进行结构优化,使真实页岩分子模型处于能量最低化,保证后续计算结果更加精确。
S504、运用MS中Forcite模块内的Geometry Optimization任务对甲烷气体分子进行结构优化,使甲烷分子模型处于能量最低化,保证后续计算结果更加精确。
进一步的,所述步骤S6中,模拟计算步骤S5中优化后的真实页岩分子模型和甲烷气体分子模型在与实验值相同温度和压力条件下的页岩气吸附量,包含以下子步骤:
S601、页岩气吸附量的计算应当应用巨正则蒙特卡洛(GCMC)的模拟方法,选用MS中的中的Sorption模块进行模拟计算。
S602、运用Sorption模块中的Fixed pressure任务,并采取Metropolis方法开展不同压力点下的页岩气吸附模拟。力场选择为Clayff力场,原子之间的相互作用由LennardJones 势和库仑项描述。分子间的范德华相互作用和静电作用分别采用Atom和Ewald求和方法,非键截断半径设置为1.25nm。
S603、在分子模拟系统中设置与实验值相同温度与压力范围。对每个压力点开展1.5×107步模拟,其中前5×106步用于吸附平衡,后1×107步用于平衡吸附量的数据统计,得到不同压力点下真实页岩分子模型单位晶胞上的甲烷分子吸附个数。
进一步的,所述步骤S7中,对比步骤6中所得真实页岩分子模型的页岩气吸附量与步骤2中实验所得页岩气吸附量之间的误差值,包含以下子步骤:
S701、步骤S3中页岩气吸附实验测得的吸附量数据属于过剩吸附量,需通过三元Langmuir非线性拟合将过剩吸附量转化为绝对吸附量,并绘制实验所得的页岩气吸附等温线。
(1)页岩气绝对吸附量计算公式
(2)页岩气过剩吸附量计算公式
(3)页岩气过剩吸附量转化为绝对吸附量
S702、将步骤S6中真实页岩分子模型单位晶胞上的甲烷分子吸附个数转化为页岩中绝对吸附量数据,绘制基于真实页岩分子模型的页岩气吸附等温线。
S703、对比实验与分子模拟得到的页岩气吸附等温线,确定模拟值与实验值间的差异。
进一步的,所述步骤S8中,模拟计算在更高的压力条件下步骤S5中优化得到的真实页岩分子模型上的页岩气吸附量。包含以下子步骤:
S801、比较真实页岩分子模型计算所得页岩气吸附量与实验值间的误差值;
S802、目前由于页岩气吸附实验测试压力普遍在30MPa以下,若误差值小于1%,则通过步骤S6中巨正则蒙特卡洛(GCMC)的模拟方法模拟计算得到更高压力下页岩气吸附量数据。
进一步的,所述步骤S9中,利用步骤S5中优化后的真实页岩分子模型和甲烷气体分子模型,模拟计算不同含水量及不同二氧化碳含量对于页岩气吸附量的影响。包含以下子步骤:
S901、运用步骤S6中所采用的巨正则蒙特卡洛(GCMC)的模拟方法,选用MS中的中的Sorption模块中的Fixed pressure任务,并采取Metropolis方法开展不同压力点下的甲烷与水蒸气以及甲烷与二氧化碳的竞争吸附模拟。力场选择为Clayff力场,原子之间的相互作用由Lennard Jones 势和库仑项描述。分子间的范德华相互作用和静电作用分别采用Atom和Ewald求和方法,非键截断半径设置为1.25nm。
S902、由于此时模拟的气体中分别包含甲烷与水蒸气、甲烷与二氧化碳等物质,因此采用Peng-Robinson(PR)状态方程分别计算出每个压力点对应的甲烷与水蒸气及甲烷与二氧化碳混合气体中每种气体受到的有效压力(逸度)。
S903、在分子模拟系统中设置与实验值相同温度以及步骤S902中计算得出的每个压力点对应的甲烷与水蒸气及甲烷与二氧化碳混合气体中每种气体受到的有效压力。对每个压力点开展1.5×107步模拟,其中前5×106步用于吸附平衡,后1×107步用于平衡吸附量的数据统计,分别得到不同压力点下不同含水量与不同二氧化碳含量对于页岩气吸附量的影响。
本发明技术效果如下:
(1)本发明提供的基于分子模拟的真实页岩气吸附量预测方法,全面获取了页岩中粘土矿物的物相种类、晶胞结构信息,衍射角、面间距d值、测量的峰强度(峰高)、对应的物相及晶面指数;同时获取了页岩中干酪根官能团种类及TOC含量占比信息;进而模拟计算得出页岩中粘土矿物与干酪根官能团的占比与排列方式,该页岩分子模拟更符合页岩地层真实情况。
(2)本发明提供的基于分子模拟的真实页岩气吸附量预测方法,选用巨正则蒙特卡洛(GCMC)模拟方法,采用Clayff力场,通过设定与实验相同的温度与压力的方式计算得出的页岩气吸附量与实验吸附量的误差非常小,以往运用分子模拟计算所得的吸附量与实验值之间误差很大,进一步说明本发明所建真实页岩分子模型的准确性以及模拟方法的先进性。
(3)本发明提供的基于分子模拟的真实页岩气吸附量预测方法,可以对国内外不同地区不同油田的页岩分子尺度的空间构型进行建模,并能准确的计算得出该地区页岩的页岩气吸附量,能够对更高温度压力条件下页岩气吸附量与吸附趋势的计算起指导作用;同时模拟计算了不同压力与温度条件下页岩气与二氧化碳及水蒸气的竞争吸附量,可为后续页岩地层所进行的勘探开发提供理论依据。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。
图1是本发明实施例1的分子模拟方法流程图;
图2是本发明实施例1页岩样品中粘土矿物XRD曲线原始谱图;
图3是本发明实施例1对页岩样品粘土矿物XRD曲线原始谱图中第一条曲线寻峰后的XRD曲线谱图;
图4是本发明实施例1对页岩样品粘土矿物XRD曲线原始谱图中第二条曲线寻峰后的XRD曲线谱图;
图5是本发明实施例1对页岩样品粘土矿物XRD曲线原始谱图中第三条曲线寻峰后的XRD曲线谱图;
图6是本发明实施例1页岩样品中蒙脱石晶胞模型图;
图7是本发明实施例1页岩样品中伊利石晶胞模型图;
图8是本发明实施例1页岩样品中高岭石晶胞模型图;
图9是本发明实施例1页岩样品中绿泥石晶胞模型图;
图10是本发明实施例1结构优化后的页岩分子模型图;
图11是本发明实施例1结构优化后的甲烷分子模型图;
图12是本发明实施例1结构优化后的水分子模型图;
图13是本发明实施例1结构优化后的二氧化碳分子模型图;
图14是本发明实施例1真实页岩分子模型进行吸附模拟后甲烷分子分布概率图;
图15是本发明实施例1页岩分子模拟在不同压力下甲烷分子吸附个数趋势图;
图16是本发明实施例1公式转换前后在不同压力下实验测得的页岩气吸附量图;
图17是本发明实施例1实验与分子模拟得到的页岩气吸附量对比图;
图18是本发明实施例1分子模拟得到的更高压力下页岩气吸附量图;
图19是本发明实施例1分子模拟得到的不同含水量下页岩气吸附量图;
图20是本发明实施例1分子模拟得到的不同CO2含量下页岩气吸附量图。
具体实施方式
实施例1:
如图1所示,一种基于分子模拟的真实页岩气吸附量预测方法,包括下述步骤:
S1、通过XRD衍射实验获取页岩中矿物成分信息,包括以下子步骤:
S101、对页岩样品进行粉末化处理,包括以下处理方式:
(1)含油页岩样品应洗油处理至荧光四级以下;
(2)将页岩样品在低于60℃的温度下烘干,冷却至室温后粉碎待用;
(3)将粉碎后的岩石样品研磨至全部粒径小于40μm或手指搓捏无颗粒感。
S102、提取岩石样品中的粘土矿物并将其制片,包括以下处理方式:
(1)把粉末化处理后的样品放在40ml高型烧杯中,加10ml蒸馏水浸泡,采用CQX-25-06型(功率为250W,频率为25kHz)超声仪对样品进行超声分散。
(2)将超声分散后的样品倒入内径为64mm的沉淀量筒,加蒸馏水至 1000m1.每个沉淀量筒在 1000ml刻度往下10cm处有小孔,用注射器插入小孔吸取悬浮液样品,实验在20℃的恒温下进行,沉降时间为8-10小时。
(3)将吸取的悬浮液离心,使粘土颗粒沉降;
(4)将离心后的样品在低于60℃电热干燥箱中烘干;
(5)将烘干后的样品用玛瑙研钵磨至手搓捏无颗粒感,用纸包好,标明样品编号;
(6)将样品粉末置于样品载片的凹槽中,使样品被测量面保持粗糙,与样品载片的表面保持平整均匀,避免晶体粉末择优取向。
S103、对制片的四川省内江市威远WY-A井页岩样品进行X射线衍射,得到XRD图谱,如图2所示;
S104、运用MDI jade软件提取XRD图谱中的粘土矿物晶体信息,如图3-图5所示,其中包括页岩中各矿物组分占比、粘土矿物各组分占比以及粘土矿物各组分晶胞参数,具体提取信息如下表所示:
表1 页岩中各矿物组分占比
表2 页岩粘土矿物各组分占比
表3 页岩粘土矿物各组分晶胞参数
S2、通过TOC实验获取页岩中有机物组分信息,包括以下子步骤:
S201、对页岩样品进行粉末化处理,包括以下处理方式:
(1)页岩试样用水洗净,于阴凉通风处晾干;
(2)将页岩试样运用粉碎机粉碎,过100目筛,置于干燥器中备用。
S202、称取1g试样置于瓷坩埚中;
S203、反复向盛有试样的瓷坩埚内少量、多次地缓慢加入HCl,放入60~80℃的烘箱内烘干,直至反应完全为止;
S204、在100~120℃烘24h,直至除尽过量的HCl;
S205、通入温度为800℃以上的氧气流使页岩样品燃烧,使其中总有机碳全部转化成二氧化碳,用热导检测器(或红外检测器)检测二氧化碳的量,用外标法计算试样的总有机碳含量。具体所得的TOC含量数据如下表所示:
表4 页岩TOC含量
S3、进行页岩气吸附实验获取吸附实验数据,包括以下子步骤:
S301、将粉碎后样品装入样品池,打开温控系统和抽真空装置,使样品池温度稳定在实验预设温度,包括以下处理方式:
(1)将页岩试样运用粉碎机粉碎,筛选粒度为40-60目不少于20g的页岩样品备用;
(2)运用温控系统与抽真空装置,使得样品池温度稳定在135℃。
S302、向样品池中充入高纯氮气(纯度不低于99.99%),采集压力数据,确保系统气密性良好;
S303、调节并设置称量系统,逐点降低样品池压力,称量对应压力点下的样品桶重量;
S304、装入实验样品并抽真空,真空状态稳定后关闭抽真空装置,向系统冲入氦气,氦气纯度不低于99.99%;
S305、调节并设置称量系统,逐点降低样品池压力,称量对应压力点下样品桶与样品总重量;
S306、对系统进行抽真空,真空状态稳定后关闭抽真空装置,向系统内冲入甲烷,控制压力为第一个实验点,实验过程中每2min采集温度、压力与天平读数,选取实验平衡后5个记录数据平均值,为实验测点平衡后,温度、压力与天平读数;
S307、逐渐增大实验压力,实验压力点分别为0、0.5MPa、1.0MPa、2.0MPa、4.0MPa、8.0MPa、12.0MPa、18.0 MPa、24.0 MPa、30.0 MPa直到最后一个压力点实验结束;
S308、记录每个压力点页岩气吸附量数据。
S4、根据步骤S1与S2中获取的页岩成分信息进行分子模拟建模,同时建立气体分子模型,包含以下子步骤:
S401、构建页岩中粘土矿物晶胞:
利用步骤S1中所得到的粘土矿物晶体信息,首先,建立晶胞(Build Crystals),设置晶胞信息,分别为晶胞空间群、基本晶胞参数,并在每种粘土矿物晶胞中输入每种原子位置信息,通过调整晶体结构使其密度达到与测量值相符的密度,具体四川省内江市威远WY-A井页岩样品中粘土矿物晶胞信息参照步骤S1中通过XRD衍射获取的页岩晶体信息;所构建的页岩单个粘土矿物分子模型,如图6-图9所示。
S402、划分岩石中粘土矿物晶面:
通过晶面划分模块,利用步骤S1中得到的页岩中每种粘土矿物的晶面指数与面间距d值,确定出岩石样品中每种粘土矿物的层间距与分割面;具体每种粘土矿物晶面信息如下:
(1)该区块岩石样品中蒙脱石、伊利石、高岭石、绿泥石均是以(0 0 1)晶面为主;
(2)该区块岩石样品中蒙脱石的衍射角为8.883,面间距d=9.9280 Å;
(3)该区块岩石样品中伊利石的衍射角为10.1193,面间距d=10.1193Å;
(4)该区块岩石样品中高岭石的衍射角为12.089,面间距d=7.3152Å;
(5)该区块岩石样品中绿泥石的衍射角为26.68,面间距d=3.3396Å。
S403、建立岩石样品中每种粘土矿物混层结构:
通过步骤S1中得出的页岩中粘土矿物每种物相含量,并结合每种物相测量的峰强度比例,获得每种粘土矿物所占比例,进而得出分子模拟中每种粘土矿物最小层间占比与分布情况;并通过MS中建立晶体层模块(Build Layers),对每种划分晶面过后的粘土矿物进行空间布局。
S404、建立页岩中有机物官能团,并随机插入所建立的粘土矿物混层结构中:
通过步骤S2中得出的TOC数据,建立6种页岩中干酪根官能团,并通过TOC含量占比换算为干酪根官能团含量占比,将不同种类的官能团随机插入所建立的粘土矿物混层结构中,页岩中有机物官能团占比情况如下表所示:
表5 页岩中有机物官能团占比
S405、建立气体分子模型:
(1)通过Materials Studio建立碳原子,然后加氢的方式建立出甲烷分子的初始模型,此时甲烷分子的键长为1.140 Å,键角为109.580°。
(2)通过Materials Studio建立氧原子,然后加碳原子方式建立出二氧化碳分子的初始模型,此时二氧化碳分子的键长为1.51 Å,键角为179.952°。
(3)通过Materials Studio建立氧原子,然后加氢的方式建立出水蒸气分子的初始模型,此时水蒸气分子的键长为1.110 Å,键角为109.471°。
S5、优化步骤S4中建立的页岩分子模型与气体分子模型,包含以下子步骤:
S501、针对步骤S4中所构建的每种粘土矿物晶胞,通过MS分子模拟软件得到所建粘土矿物晶胞密度信息,并与步骤S1中得到的每种粘土矿物密度测量值对比,通过微调晶体中原子位置,键长及增减原子个数的方式,使得所建模型密度值与所测密度值相符;即,蒙脱石晶体通过调整晶体结构使其密度达到0.637g/cm3;伊利石晶体通过调整晶体结构使其密度达到1.44g/cm3;高岭石晶体通过调整晶体结构使其密度达到2.28g/cm3;绿泥石晶体通过调整晶体结构使其密度达到3.1g/cm3。
S502、针对步骤S4中所构建的岩石样品中各粘土矿物最终的晶胞模型,通过MS分子模拟软件得到所建模型总的晶胞密度信息,并与步骤S1中得到的页岩中粘土矿物晶胞信息的密度值对比,通过调整每层粘土矿物层间距的方式,使得通过所建模型得到的密度值与步骤S1中所测密度值相符。
S503、运用MS中Forcite模块内的Geometry Optimization任务对步骤页岩分子模型进行结构优化,利用Smart优化算法对结构进行最大5000次迭代,收敛精度设置为精细级别。力场选择为Clayff力场,原子之间的相互作用由Lennard Jones 势和库仑项描述。分子间的范德华相互作用和静电作用采用Atom求和方法,使页岩分子模型处于能量最低化,保证后续计算结果更加精确,最终结构优化后而得的页岩分子模型,如图10所示。
S504、运用MS中Forcite模块内的Geometry Optimization任务对气体分子进行结构优化,使气体分子模型处于能量最低化,优化后甲烷分子的键长为1.099 Å,键角为109.471°;二氧化碳分子的键长为1.16Å,键角为180°;水蒸气分子的键长为0.970 Å,键角为103.701°,如图11-图13所示。
S6,模拟计算步骤S5中优化后的真实页岩分子模型和气体分子模型在与实验值相同温度与压力条件下的页岩气吸附量,包含以下子步骤:
S601、页岩气吸附量的计算应当应用巨正则蒙特卡洛(GCMC)的模拟方法,选用MS中的中的Sorption模块进行模拟计算:
页岩气吸附模拟时,确定的宏观量有体系的体积 V,页岩气所处的环境温度T,以及所处压力条件下甲烷分子的化学势能μ,所以在模拟页岩气吸附时采取巨正则系综(μVT系综),因此使用巨正则蒙特卡罗(GCMC)方法进行模拟。
S602、运用Sorption模块中的Fixed pressure任务,并采取Metropolis方法开展不同压力点下的页岩气吸附模拟。力场选择为Clayff力场,原子之间的相互作用由LennardJones 势和库仑项描述。分子间的范德华相互作用和静电作用分别采用Atom和Ewald求和方法,非键截断半径设置为1.25nm。模拟计算生成页岩分子模型进行吸附模拟后甲烷分子分布概率图,如图14所示。分子模拟中力场选用对于最终模拟结果的准确性影响巨大,这里分子模拟选取Clayff 力场原因如下:
Clayff 力场是适用于模拟物质水合和多组分矿物体系及其与水溶液的界面,原子间相互作用势来源于各种简单水合化合物的结构和参数化的光谱数据。该力场使用基于SPC的键角可弯曲水模型来描述水和羟基的行为。原子之间的相互作用由Lennard Jones势和库仑项描述,力场中原子电荷来自于密度泛函理论计算所得结果。该力场将大多数原子间的相互作用使用非键形式进行描述。这种非键描述方式使得该力场在描述多相相互作用时具有广泛适用性,且可考虑流体相与固体之间的能量和动量交换,该力场也可保证足够少的参数数量,以允许构建较大且高度无序的系统。通过以上描述可知,Clayff 的特点适用于描述烷烃、水与粘土矿物之间的相互作用。
S603、在分子模拟系统中设置与实验值相同温度与压力范围,温度设置为135℃,压力分别设置为0、0.5MPa、1.0MPa、2.0MPa、4.0MPa、8.0MPa、12.0MPa、18.0 MPa、24.0 MPa、30.0 MPa。对每个压力点开展1.5×107步模拟,其中前5×106步用于吸附平衡,后1×107步用于平衡吸附量的数据统计,得到不同压力点下页岩分子模型单位晶胞上的甲烷分子吸附个数并绘制不同压力下甲烷分子吸附个数趋势图,如图15所示。
S7,对比步骤6中所得分子模型的页岩气吸附量与步骤2中实验所得页岩气吸附量之间的误差值,包含以下子步骤:
S701、步骤S3中页岩气吸附实验测得的吸附量数据属于过剩吸附量,需通过三元Langmuir非线性拟合将过剩吸附量转化为绝对吸附量,并绘制实验所得的页岩气吸附等温线,如图16所示。
(1)页岩气绝对吸附量计算公式
(2)页岩气绝对吸附量计算公式
(3)页岩气过剩吸附量转化为绝对吸附量
S702、将步骤S6中页岩分子模型单位晶胞上的甲烷分子吸附个数转化为页岩中绝对吸附量数据,绘制基于页岩分子模型的页岩气吸附等温线。
S703、对比实验与分子模拟得到的页岩气吸附等温线,如图17所示。确定分子模拟值与实验值间的差异。
S8、模拟计算在更高的压力条件下步骤S5中优化得到的页岩分子模型上的页岩气吸附量,包含以下子步骤:
S801、比较页岩分子模型计算所得页岩气吸附量与实验值间的误差值,在每个压力点上,分子模拟值与实验值的误差均小于1%。
S802、由于页岩气吸附实验测试压力普遍不超过30MPa,通过分子模拟的方法模拟计算得到30MPa-80MPa页岩气吸附量数据,并绘制得到更高压力下页岩气吸附量图,如图18所示。
S9、利用步骤S5中优化后的真实页岩分子模型和气体分子模型,模拟计算不同含水量及不同二氧化碳含量对于页岩气吸附量的影响,包含以下子步骤:
S901、运用步骤S6中所采用的巨正则蒙特卡洛(GCMC)的模拟方法,选用MS中的中的Sorption模块中的Fixed pressure任务,并采取Metropolis方法开展不同压力点下的甲烷与水蒸气以及甲烷与二氧化碳的竞争吸附模拟。力场选择为Clayff力场,原子之间的相互作用由Lennard Jones 势和库仑项描述。分子间的范德华相互作用和静电作用分别采用Atom和Ewald求和方法,非键截断半径设置为1.25nm。
S902、由于此时模拟的气体中分别包含甲烷与水蒸气、甲烷与二氧化碳等物质,因此采用Peng-Robinson(PR)状态方程分别计算出0、0.5MPa、1.0MPa、2.0MPa、4.0MPa、8.0MPa、12.0MPa、18.0 MPa、24.0 MPa、30.0 MPa每个压力点对应的甲烷与水蒸气及甲烷与二氧化碳混合气体中每种气体受到的有效压力(逸度)。
S903、在分子模拟系统中设置与实验值相同温度,温度设置为135℃以及步骤S902中计算得出的每个压力点对应的甲烷与水蒸气及甲烷与二氧化碳混合气体中每种气体受到的有效压力,设置了4种含有甲烷与水蒸气的混合气体,其中的含水量分别为0.5%、1.0%、1.5%、2.0%;另外设置了4种含有甲烷与CO2的混合气体,其中的CO2含量分别为5%、10%、15%、20%。对每个压力点开展1.5×107步模拟,其中前5×106步用于吸附平衡,后1×107步用于平衡吸附量的数据统计。得到不同压力点下不同含水量的混合气体中页岩气吸附量的变化趋势,如图19所示;以及不同压力点下不同CO2含量的混合气体中页岩气吸附量的变化趋势,如图20所示。
(1)由图19可知,页岩气混合气体在相同压力点下,甲烷吸附量随含水量的增加而降低,并且降低趋势逐渐放缓;
(2)由图20可知,页岩气混合气体在相同压力点下,甲烷吸附量随CO2含量的增加而降低,并且降低趋势逐渐放缓;
(3)运用该分子模拟方法计算得出的页岩气中甲烷与水以及甲烷与CO2竞争吸附的结果对国内外页岩气田勘探与注气开采有指导作用。
Claims (9)
1.一种基于分子模拟的真实页岩气吸附量预测方法,其特征在于,包括以下步骤:
S1、通过XRD衍射实验获取页岩中矿物成分信息,所述矿物成分信息包括矿物各组分占比、矿物各组分晶胞参数及密度值;
S2、通过TOC实验获取页岩中有机物组分信息,所述页岩中有机物组分信息包括有机物官能团信息及含量;
S3、进行页岩气吸附实验,获取吸附实验数据,所述吸附实验数据包括不同压力点下的页岩气吸附量数据,用于与步骤S6所得模拟值进行对比;
S4、根据S1与S2中获取的矿物成分信息和有机物组分信息,进行真实页岩分子模型构建,同时建立甲烷气体分子模型;
S5、优化S4中建立的真实页岩分子模型及甲烷气体分子模型;
S6、利用S5中优化后的真实页岩分子模型和甲烷气体分子模型,计算不同压力、温度下的真实页岩气吸附量;
S7、对比S6计算出的真实页岩气吸附量与S3中的吸附实验数据,计算误差;
S8、计算比S6中记载的压力更高的压力下的页岩气吸附量;
S9、利用S5中优化后的真实页岩分子模型和甲烷气体分子模型,模拟计算不同含水量及不同二氧化碳含量对于页岩气吸附量的影响;
所述步骤S5中,优化步骤S4中建立的真实页岩分子模型及甲烷气体分子模型,包括以下子步骤:
S501、通过微调每种粘土矿物晶体中原子位置、键长及增减原子个数的方式,使得通过所建模型得到的密度值与S1中所测密度值相符;
S502、基于步骤S1中XRD衍射数据中粘土矿物晶胞面间距与密度值,通过调整所构建的页岩中各粘土矿物最终的晶胞模型每层粘土矿物层间距的方式,此优化方式为分子建模常用方式,使得通过所建模型得到的密度值与S1中所测密度值相符;
S503、对真实页岩分子模型进行结构优化,使真实页岩分子模型处于能量最低化,保证后续计算结果更加精确;
S504、对甲烷气体分子模型进行结构优化,使甲烷气体分子模型处于能量最低化,保证后续计算结果更加精确。
2.根据权利要求1所述的一种基于分子模拟的真实页岩气吸附量预测方法,其特征在于,所述S1包括以下子步骤:
S101:对页岩样品进行制片;
S102:对制片的页岩样品进行X射线衍射得到XRD图谱;
S103:运用MDI jade软件提取得到XRD图谱中的页岩粘土矿物成分信息。
3.根据权利要求1所述的一种基于分子模拟的真实页岩气吸附量预测方法,其特征在于,所述S2中的TOC实验包括以下子步骤:
S201、对页岩样品进行粉末化处理;
S202、反复向盛有试样的瓷坩埚内少量、多次地缓慢加入HCl,放入60~80℃的烘箱内烘干,直至反应完全为止;
S203、100~120℃烘24h,直至除尽过量的HCl;
S204、通入温度为800℃以上的氧气流使页岩样品燃烧,使其中总有机碳全部转化成二氧化碳,用热导检测器检测二氧化碳的量,用外标法计算试样的有机物组分信息。
4.根据权利要求1所述的一种基于分子模拟的真实页岩气吸附量预测方法,其特征在于,所述步骤S3中,进行页岩气吸附实验获取吸附实验数据,包括以下子步骤:
S301、装入实验样品并抽真空,当真空状态稳定后,关闭抽真空装置,向系统冲入氦气;
S302、调节并设置称量系统,逐点降低样品池压力,称量对应压力点下样品桶与样品总重量;
S303、对系统进行抽真空,当真空状态稳定后,关闭抽真空装置,向系统内冲入甲烷,控制压力为第一个实验点;
S304、逐渐增大实验压力,直到最后一个压力点实验结束;
S305、记录每个压力点页岩气吸附量数据。
5.根据权利要求1所述的一种基于分子模拟的真实页岩气吸附量预测方法,其特征在于:所述S4中建立真实页岩分子模型与甲烷气体分子模型包括以下子步骤:
S401、构建岩石中粘土矿物晶胞;
S402、划分岩石中粘土矿物晶面;
S403、建立岩石样品中每种粘土矿物混层结构;
S404、建立页岩中有机物官能团,并随机插入所建立的粘土矿物混层结构中;
S405、建立真实页岩分子模型与甲烷气体分子模型。
6.根据权利要求1所述的一种基于分子模拟的真实页岩气吸附量预测方法,其特征在于,所述S6中,模拟计算S5中优化后的真实页岩分子模型和甲烷气体分子模型在与实验值相同温度和压力条件下的页岩气吸附量,应用巨正则蒙特卡洛的模拟方法,选用MS中的Sorption模块中的Fixed pressure任务,并采取Metropolis方法开展不同压力点下的页岩气吸附模拟;力场选择为Clayff力场,原子之间的相互作用由Lennard Jones 势和库仑项描述;得到不同压力点下真实页岩分子模型单位晶胞上的甲烷分子吸附量。
7.根据权利要求1所述的一种基于分子模拟的真实页岩气吸附量预测方法,其特征在于,所述S7中,对比S6中所得真实页岩分子模型的页岩气吸附量与S3中实验所得页岩气吸附量之间的误差值,包含以下子步骤:
S701、S3中页岩气吸附实验测得的吸附量数据属于过剩吸附量,需通过三元Langmuir非线性拟合将过剩吸附量转化为绝对吸附量,并绘制实验所得的页岩气吸附等温线;
S702、将S6中真实页岩分子模型中单位晶胞上的甲烷分子吸附个数转化为页岩中绝对吸附量数据,绘制基于真实页岩分子模型的页岩气吸附温线;
S703、对比实验与分子模拟得到的页岩气吸附温线,确定模拟值与实验值间的差异。
8.根据权利要求1所述的一种基于分子模拟的真实页岩气吸附量预测方法,其特征在于,所述S8中,模拟计算在更高的压力条件下步骤S5中优化得到的真实页岩分子模型上的页岩气吸附量。
9.根据权利要求6所述的一种基于分子模拟的真实页岩气吸附量预测方法,其特征在于,S9中,利用S5中优化后的真实页岩分子模型和甲烷气体分子模型,应用S6中所采用的巨正则蒙特卡洛的模拟方法模拟计算不同含水量对于页岩气吸附量的影响以及不同二氧化碳含量对于页岩气吸附量的影响。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110837284.XA CN113670960B (zh) | 2021-07-23 | 2021-07-23 | 一种基于分子模拟的真实页岩气吸附量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110837284.XA CN113670960B (zh) | 2021-07-23 | 2021-07-23 | 一种基于分子模拟的真实页岩气吸附量预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113670960A CN113670960A (zh) | 2021-11-19 |
CN113670960B true CN113670960B (zh) | 2021-12-28 |
Family
ID=78540125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110837284.XA Active CN113670960B (zh) | 2021-07-23 | 2021-07-23 | 一种基于分子模拟的真实页岩气吸附量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113670960B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114492089B (zh) * | 2022-04-06 | 2022-07-22 | 西南石油大学 | 一种基于分子模拟的页岩凝析气吸附量预测方法 |
CN115165700B (zh) * | 2022-06-17 | 2023-08-25 | 中国石油大学(北京) | 全岩心尺度二氧化碳吸附封存量的实验装置 |
CN117494452B (zh) * | 2023-11-14 | 2024-05-17 | 中国地质大学(北京) | 非常规储层吸附气量的确定方法、装置、终端及存储介质 |
Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102162785A (zh) * | 2011-01-11 | 2011-08-24 | 中国石油天然气股份有限公司 | 高温高压吸附测试仪 |
CN102607991A (zh) * | 2012-03-02 | 2012-07-25 | 中国石油天然气股份有限公司 | 煤/页岩吸附量测定装置 |
CN102944666A (zh) * | 2012-12-05 | 2013-02-27 | 西南石油大学 | 页岩气藏开采模拟实验装置 |
CN103994943A (zh) * | 2014-05-23 | 2014-08-20 | 中国石油大学(北京) | 一种煤/页岩等温吸附实验装置 |
CN104239602A (zh) * | 2014-07-15 | 2014-12-24 | 浙江科技学院 | 一种页岩气藏特性模拟方法 |
EP2928579A1 (en) * | 2012-12-06 | 2015-10-14 | ExxonMobil Research and Engineering Company | Selectivation of adsorbents for gas separation |
WO2016028564A1 (en) * | 2014-08-22 | 2016-02-25 | Schlumberger Canada Limited | Methods for monitoring fluid flow and transport in shale gas reservoirs |
AU2013403405A1 (en) * | 2013-10-23 | 2016-03-31 | Halliburton Energy Services, Inc. | Volatile surfactant treatment for subterranean formations |
CN106018167A (zh) * | 2016-07-27 | 2016-10-12 | 中国石油大学(华东) | 一种对等温吸附实验中吸附质气体的过剩吸附量进行校正的方法 |
CN106290061A (zh) * | 2016-07-27 | 2017-01-04 | 中国石油大学(华东) | 测定吸附质气体在吸附质上的吸附相密度的方法 |
CN106326519A (zh) * | 2015-07-09 | 2017-01-11 | 中国科学院微电子研究所 | 一种确定表面活性剂溶解吸附行为的方法及系统 |
CN106940279A (zh) * | 2017-05-18 | 2017-07-11 | 中国石油大学(华东) | 一种评价泥页岩储层赋存吸附气量的方法 |
CN106970000A (zh) * | 2017-04-21 | 2017-07-21 | 西南石油大学 | 煤/页岩超高压气体吸附和渗流实验评价页岩气吸附方法 |
CN108460219A (zh) * | 2018-03-16 | 2018-08-28 | 中石化华东石油工程有限公司测井分公司 | 基于超临界单分子层吸附模型计算页岩吸附气含气量方法 |
WO2018161116A1 (en) * | 2017-03-06 | 2018-09-13 | Commonwealth Scientific And Industrial Research Organisation | Permeable graphene and permeable graphene membranes |
CN109854236A (zh) * | 2019-02-25 | 2019-06-07 | 中国石油大学(华东) | 一种用于页岩油流动的数值模拟方法及装置 |
CN110489927A (zh) * | 2019-09-02 | 2019-11-22 | 西南石油大学 | 页岩吸附气吸附相密度模型的构建及绝对吸附量计算方法 |
CN111007233A (zh) * | 2019-12-25 | 2020-04-14 | 西南石油大学 | 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法 |
CN111307684A (zh) * | 2019-12-04 | 2020-06-19 | 西南石油大学 | 一种用于计算微纳米孔隙中气体渗透率的分子模拟方法 |
CN111912957A (zh) * | 2020-08-24 | 2020-11-10 | 东北石油大学 | 一种基于赋存状态的页岩油量检测方法及系统 |
CN112132337A (zh) * | 2020-09-22 | 2020-12-25 | 中国矿业大学(北京) | 一种预测煤层瓦斯参数的方法 |
CN112414891A (zh) * | 2020-10-27 | 2021-02-26 | 成都理工大学 | 一种页岩气吸附回滞现象的分子模拟方法 |
CN112613174A (zh) * | 2020-12-18 | 2021-04-06 | 西南石油大学 | 考虑多重吸附机制的页岩吸附甲烷能力评估方法 |
CN112986096A (zh) * | 2021-02-25 | 2021-06-18 | 中国地质大学(北京) | 一种地下含水泥页岩储层吸附气量快速估算方法 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070274867A1 (en) * | 2005-02-28 | 2007-11-29 | Honeywell International Inc. | Stationary phase for a micro fluid analyzer |
US8869889B2 (en) * | 2010-09-21 | 2014-10-28 | Palmer Labs, Llc | Method of using carbon dioxide in recovery of formation deposits |
US20140212982A1 (en) * | 2013-01-29 | 2014-07-31 | University Of Saskatchewan | Methods of selectively detecting the presence of a compound in a gaseous medium |
US10088597B2 (en) * | 2013-08-27 | 2018-10-02 | Schlumberger Technology Corporation | Determining phase behavior of a reservoir fluid |
US10160904B2 (en) * | 2013-10-23 | 2018-12-25 | Halliburton Energy Services, Inc. | Volatile surfactant treatment for subterranean formations |
US10151737B2 (en) * | 2015-04-08 | 2018-12-11 | King Fahd University Of Petroleum And Minerals | Method for permeability prediction of shale gas |
US10591399B2 (en) * | 2015-07-17 | 2020-03-17 | Saudi Arabian Oil Company | Methods for analyzing natural gas flow in subterranean reservoirs |
MX2016002733A (es) * | 2016-03-02 | 2017-09-01 | Univ Nacional Autónoma De México | Proceso de sintesis para la obtencion de sales cationicas cuaternarias utiles como trazadores quimicos en yacimientos carbonatados, y productos obtenidos con el mismo. |
CN106568922A (zh) * | 2016-10-19 | 2017-04-19 | 中国石油天然气股份有限公司 | 一种地层温度和压力条件下页岩吸附气含量的计算方法 |
CN107462508B (zh) * | 2017-08-16 | 2018-10-02 | 西南石油大学 | 一种多场耦合渗流多功能实验装置及测试方法 |
CN107422100B (zh) * | 2017-09-05 | 2019-09-03 | 中国石油大学(华东) | 一种计算页岩气储层赋存吸附气含量的方法 |
WO2019133191A1 (en) * | 2017-12-29 | 2019-07-04 | Exxonmobil Research And Engineering Company | Lubrication of oxygenated diamond-like carbon surfaces |
US11506815B2 (en) * | 2018-03-15 | 2022-11-22 | The Penn State Research Foundation | Method for improved recovery in ultra-tight reservoirs based on diffusion |
CA3042922A1 (en) * | 2018-05-15 | 2019-11-15 | Uti Limited Partnership | Method and system for measuring a parameter of a core sample |
CN112162000A (zh) * | 2020-08-24 | 2021-01-01 | 长江大学 | 一种基于页岩基质组分进行吸附气含量预测方法及系统 |
CN111912960B (zh) * | 2020-08-24 | 2022-11-08 | 东北石油大学 | 一种泥页岩中干酪根游离油定量评价方法及系统 |
-
2021
- 2021-07-23 CN CN202110837284.XA patent/CN113670960B/zh active Active
Patent Citations (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102162785A (zh) * | 2011-01-11 | 2011-08-24 | 中国石油天然气股份有限公司 | 高温高压吸附测试仪 |
CN102607991A (zh) * | 2012-03-02 | 2012-07-25 | 中国石油天然气股份有限公司 | 煤/页岩吸附量测定装置 |
CN102944666A (zh) * | 2012-12-05 | 2013-02-27 | 西南石油大学 | 页岩气藏开采模拟实验装置 |
EP2928579A1 (en) * | 2012-12-06 | 2015-10-14 | ExxonMobil Research and Engineering Company | Selectivation of adsorbents for gas separation |
AU2013403405A1 (en) * | 2013-10-23 | 2016-03-31 | Halliburton Energy Services, Inc. | Volatile surfactant treatment for subterranean formations |
CN103994943A (zh) * | 2014-05-23 | 2014-08-20 | 中国石油大学(北京) | 一种煤/页岩等温吸附实验装置 |
CN104239602A (zh) * | 2014-07-15 | 2014-12-24 | 浙江科技学院 | 一种页岩气藏特性模拟方法 |
WO2016028564A1 (en) * | 2014-08-22 | 2016-02-25 | Schlumberger Canada Limited | Methods for monitoring fluid flow and transport in shale gas reservoirs |
CN106326519A (zh) * | 2015-07-09 | 2017-01-11 | 中国科学院微电子研究所 | 一种确定表面活性剂溶解吸附行为的方法及系统 |
CN106018167A (zh) * | 2016-07-27 | 2016-10-12 | 中国石油大学(华东) | 一种对等温吸附实验中吸附质气体的过剩吸附量进行校正的方法 |
CN106290061A (zh) * | 2016-07-27 | 2017-01-04 | 中国石油大学(华东) | 测定吸附质气体在吸附质上的吸附相密度的方法 |
WO2018161116A1 (en) * | 2017-03-06 | 2018-09-13 | Commonwealth Scientific And Industrial Research Organisation | Permeable graphene and permeable graphene membranes |
CN106970000A (zh) * | 2017-04-21 | 2017-07-21 | 西南石油大学 | 煤/页岩超高压气体吸附和渗流实验评价页岩气吸附方法 |
CN106940279A (zh) * | 2017-05-18 | 2017-07-11 | 中国石油大学(华东) | 一种评价泥页岩储层赋存吸附气量的方法 |
CN108460219A (zh) * | 2018-03-16 | 2018-08-28 | 中石化华东石油工程有限公司测井分公司 | 基于超临界单分子层吸附模型计算页岩吸附气含气量方法 |
CN109854236A (zh) * | 2019-02-25 | 2019-06-07 | 中国石油大学(华东) | 一种用于页岩油流动的数值模拟方法及装置 |
CN110489927A (zh) * | 2019-09-02 | 2019-11-22 | 西南石油大学 | 页岩吸附气吸附相密度模型的构建及绝对吸附量计算方法 |
CN111307684A (zh) * | 2019-12-04 | 2020-06-19 | 西南石油大学 | 一种用于计算微纳米孔隙中气体渗透率的分子模拟方法 |
CN111007233A (zh) * | 2019-12-25 | 2020-04-14 | 西南石油大学 | 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法 |
CN111912957A (zh) * | 2020-08-24 | 2020-11-10 | 东北石油大学 | 一种基于赋存状态的页岩油量检测方法及系统 |
CN112132337A (zh) * | 2020-09-22 | 2020-12-25 | 中国矿业大学(北京) | 一种预测煤层瓦斯参数的方法 |
CN112414891A (zh) * | 2020-10-27 | 2021-02-26 | 成都理工大学 | 一种页岩气吸附回滞现象的分子模拟方法 |
CN112613174A (zh) * | 2020-12-18 | 2021-04-06 | 西南石油大学 | 考虑多重吸附机制的页岩吸附甲烷能力评估方法 |
CN112986096A (zh) * | 2021-02-25 | 2021-06-18 | 中国地质大学(北京) | 一种地下含水泥页岩储层吸附气量快速估算方法 |
Non-Patent Citations (8)
Title |
---|
A New Production Data Analysis Method of Shale Gas: Based on Flowing Material Balance Theory and Considering the Complex Flow Mechanisms of Multiple Pressure Systems;Xutao You.et;《Procedia》;20191231;第158卷;第3626-3632页 * |
Methane and CO2 Adsorption Capacities of Kerogen in the Eagle Ford Shale from Molecular Simulation;Peter Psarras.et;《Accounts》;20171231;第1-11页 * |
Quantification of CH4 adsorption capacity in kerogen-rich reservoir shales: An experimental investigation and molecular dynamic simulation;Yang Ju.et;《Energy》;20191231;第1-39页 * |
Review of Molecular Simulation Method for Gas Adsorption/desorption and Diffusion in Shale Matrix;WANG Hui.et;《Journal of Thermal Science》;20191231;第28卷(第1期);第1-16页 * |
有机质纳米孔隙吸附页岩气的分子模拟;张廷山等;《天然气地球科学》;20171231;第28卷(第01期);第146-155页 * |
页岩储层纳微米孔隙CO2/CH4吸附及驱替特性研究进展;邓佳等;《力学学报》;20201018;第53卷(第10期);第2880-2890页 * |
页岩气吸附的分子模拟及实验研究;孟伟杰等;《中国优秀硕士学位论文全文数据库 (基础科学辑)》;20170331;第A011-440页 * |
页岩纳米孔隙中小分子气体超临界吸附微观机理研究;代建伟;《中国优秀硕士学位论文全文数据库 (工程科技Ⅰ辑)》;20170331;第B019-315页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113670960A (zh) | 2021-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113670960B (zh) | 一种基于分子模拟的真实页岩气吸附量预测方法 | |
Zhai et al. | Adsorption and diffusion of shale gas reservoirs in modeled clay minerals at different geological depths | |
Wang et al. | Pore characterization of different clay minerals and its impact on methane adsorption capacity | |
Li et al. | Molecular simulation study on methane adsorption capacity and mechanism in clay minerals: Effect of clay type, pressure, and water saturation in shales | |
Yang et al. | Experimental investigation of water vapor adsorption isotherm on gas-producing Longmaxi shale: Mathematical modeling and implication for water distribution in shale reservoirs | |
Eliebid et al. | Adsorption role in shale gas recovery and the feasibility of CO2 in shale enhanced gas recovery: a study on shale gas from Saudi Arabia | |
Wang et al. | High-pressure hydrogen adsorption in clay minerals: Insights on natural hydrogen exploration | |
Cao et al. | Molecular Simulation of Methane Adsorption Behaviour in Kerogen Nanopores for Shale Gas Resource Assessment | |
CN105717202A (zh) | 表征页岩吸附解吸天然气过程中组分和同位素分馏作用的方法 | |
Zhang et al. | Molecular dynamics simulations of CH 4 diffusion in kaolinite: Influence of water content | |
Ortiz Cancino et al. | Effect of organic matter and thermal maturity on methane adsorption capacity on shales from the middle Magdalena Valley Basin in Colombia | |
Zhang et al. | Monte Carlo simulation of the adsorption and displacement of CH4 by CO2 Injection in shale organic carbon slit micropores for CO2 enhanced shale gas recovery | |
Song et al. | Molecular simulation on competitive adsorption mechanism of CH 4/CO 2 on shale kerogen | |
Gao et al. | Study on the adsorption and thermodynamic characteristics of methane under high temperature and pressure | |
Babaei et al. | Methane adsorption of nanocomposite shale in the presence of water: Insights from molecular simulations | |
Hu et al. | Insights into shale gas adsorption and an improved method for characterizing adsorption isotherm from molecular perspectives | |
Liu et al. | Molecular simulation calculation method for shale gas adsorption in realistic shale model | |
Liang et al. | Adsorption models for shale gas: A mini-review | |
Wang et al. | Carbon isotope fractionation during shale gas transport through organic and inorganic nanopores from molecular simulations | |
Wang et al. | Characteristics and influencing factors of supercritical methane adsorption in deep gas shale: A case study of marine Wufeng and Longmaxi formations from the Dongxi area, southeastern Sichuan basin (China) | |
Tao et al. | Displacement of shale gas confined in illite shale by flue gas: A molecular simulation study | |
Ukaomah et al. | Experimental and molecular investigation of water adsorption controls in marine and lacustrine shale reservoirs | |
Mu et al. | Sorption of deep shale gas on minerals and organic matter from molecular simulation | |
Qian et al. | Methane adsorption characteristics under in situ reservoir conditions of the Wufeng–Longmaxi shale in Southern Sichuan Basin, China: Implications for gas content evaluation | |
Ma et al. | Assessing impact of shale gas adsorption on free-gas permeability via a pore network flow model |
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 |