CN111007233B - 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法 - Google Patents

一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法 Download PDF

Info

Publication number
CN111007233B
CN111007233B CN201911361351.4A CN201911361351A CN111007233B CN 111007233 B CN111007233 B CN 111007233B CN 201911361351 A CN201911361351 A CN 201911361351A CN 111007233 B CN111007233 B CN 111007233B
Authority
CN
China
Prior art keywords
shale
adsorption
reservoir
pore
gas
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
Application number
CN201911361351.4A
Other languages
English (en)
Other versions
CN111007233A (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.)
Petrochina Co ltd Zhejiang Oil Field Branch
Southwest Petroleum University
Original Assignee
Petrochina Co ltd Zhejiang Oil Field Branch
Southwest Petroleum University
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 Petrochina Co ltd Zhejiang Oil Field Branch, Southwest Petroleum University filed Critical Petrochina Co ltd Zhejiang Oil Field Branch
Priority to CN201911361351.4A priority Critical patent/CN111007233B/zh
Publication of CN111007233A publication Critical patent/CN111007233A/zh
Application granted granted Critical
Publication of CN111007233B publication Critical patent/CN111007233B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01N33/24—Earth materials
    • 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
    • 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
    • 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/207—Diffractometry using detectors, e.g. using a probe in a central position and one or more displaceable detectors in circumferential positions
    • 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/22—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 measuring secondary emission from the material
    • G01N23/225—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 measuring secondary emission from the material using electron or ion
    • G01N23/2251—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 measuring secondary emission from the material using electron or ion using incident electron beams, e.g. scanning electron microscopy [SEM]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01Q—SCANNING-PROBE TECHNIQUES OR APPARATUS; APPLICATIONS OF SCANNING-PROBE TECHNIQUES, e.g. SCANNING PROBE MICROSCOPY [SPM]
    • G01Q60/00—Particular types of SPM [Scanning Probe Microscopy] or microscopes; Essential components thereof
    • G01Q60/24—AFM [Atomic Force Microscopy] or apparatus therefor, e.g. AFM probes
    • 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
    • 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
    • G01N2015/0866—Sorption
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/05—Investigating materials by wave or particle radiation by diffraction, scatter or reflection
    • G01N2223/056—Investigating materials by wave or particle radiation by diffraction, scatter or reflection diffraction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/07—Investigating materials by wave or particle radiation secondary emission
    • G01N2223/09—Investigating materials by wave or particle radiation secondary emission exo-electron emission
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/10—Different kinds of radiation or particles
    • G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
    • G01N2223/1016—X-ray
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/10—Different kinds of radiation or particles
    • G01N2223/102—Different kinds of radiation or particles beta or electrons
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/401—Imaging image processing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/60—Specific applications or type of materials
    • G01N2223/649—Specific applications or type of materials porosity

Landscapes

  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Dispersion Chemistry (AREA)
  • Theoretical Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Remote Sensing (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Geology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了一种提供一种分析页岩微观孔隙中甲烷‑二氧化碳运动行为的方法，利用数值模拟研究在多场耦合条件下页岩中CO2/CH4的动力学特点以及在微观多孔介质中运动的特性，建立页岩中CO2/CH4吸附与扩散模型，为CO2高效置换页岩中的CH4，并为CO2在页岩中封存研究提供理论依据。通过且仅仅通过将基于CCSIM算法和信息熵的三维建模方法和基于AFM的三维建模方法两种方法结合起来对页岩非均质储层三维孔隙结构特征进行表征，只需要选取1‑3个特征横切面与CCSIM算法和信息熵三维建模法相结合，仅仅通过这两种三维建模方法，即可真实、经济、准确地表征页岩非均质储层三维孔隙结构特征，并且对于大尺寸的样品分析准确性大大提高，在实际生产中具有很好的应用前景。

Description

一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法
技术领域
本发明涉及气体运动行为分析领域，具体为一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法。
背景技术
近年来，国内外学者利用不同方法进行了粘土矿物对 CH4、CH4/CO2混合气体的吸附行为研究。Liu等(2015)用分子动力学的方法模拟了不同比例的 CH4/CO2 混合物在碳纳米管中吸附、解吸的异同, 曹伟等(2014)通过 GCMC 和 MD 方法研究纯组分和等物质的量的二组分 CO2/CH4混合物在CNT bundle中的吸附和扩散，探索了在膜孔道内受限的生物沼气的分离机理。He等(2015)也证实了用碳纳米管模型模拟页岩储层有机质纳米孔隙对页岩气的吸附和扩散的动力学特征的可行性。
田园媛（2018）通过建立伊利石、蒙脱石及高岭石粘土矿物平板型孔隙模型，运用力场叠加相关准则，模拟不同温度、压力、粘土矿物类型及纳米孔隙大小条件下的甲烷吸附以及考虑不同体相密度摩尔比的二氧化碳/甲烷竞争吸附。同时开展惰性气体（氦气）在不同粘土矿物模型中的吸附模拟，以确定考虑分子间斥力的有效孔隙体积。
左骁遥(2019)运用Materials Studio软件构建了高岭石、蒙脱石、伊利石以及1:1伊蒙混层几种常见的粘土矿物的晶体结构模型，采用量子力学和分子力学方法进行了优化，并建立这几种矿物宽度为2nm、4nm、6nm孔隙结构模型，然后利用蒙特卡罗模拟方法研究了进行二氧化碳地质封存的主要埋深下几种粘土矿物孔隙结构对二氧化碳的吸附特点，再运用分子动力学模拟探讨了几种粘土矿物孔隙结构对二氧化碳吸附机理。最后分析比较同一条件下几种粘土矿物对二氧化碳吸附特征。
自 2011 年来，页岩气在页岩中吸附的分子模拟逐渐成为研究热点。而对CH4/CO2吸附的分子模拟研究确相对不多，且主要集中于对干酪根、单一粘土矿物或煤的模拟研究，都未结合页岩实际地质埋深条件、CH4/CO2在页岩微观孔隙中的分布及吸附扩散机理进行分析解释，对页岩微观孔隙的研究手段还不够丰富。
发明内容
本发明的目的在于克服现有技术的不足，提供一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法，利用数值模拟研究在多场耦合条件下页岩中 CO2/CH4的动力学特点以及在微观多孔介质中运动的特性，建立页岩中 CO2/CH4 吸附与扩散模型，为 CO2高效置换页岩中的 CH4，并为 CO2在页岩中封存研究提供理论依据。
本发明的目的是通过以下技术方案来实现的：
一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法，其特征在于包括如下步骤：
（1） 现场资料收集：收集野外剖面样品和岩心资料、随钻测试、压裂排采数据，掌握开发页岩气层在地层条件下的压力、温度、区域应力特征参数，用以页岩储层特征的表述和模拟；
（2）页岩特征分析：在区域地质演化研究的基础上，通过薄片鉴定、X衍射、扫描电镜，结合测井资料等对研究区页岩岩芯样品进行矿物成分和含量分析，研究粘土矿物中蒙脱石向伊利石转化过程；进行有机碳测定、干酪根显微组分及类型、有机质成熟度、微观古生物实验分析鉴定，确定泥页岩的有机质性质，为页岩储层微观特征表征提供依据；
（3）页岩气非均质储层微观孔隙多尺度定量、定性表征：采用直接观测和气体吸附间接数值测定两种实验方法来表征页岩微观储集空间特征，1)微观储集空间二维形态、结构与大小的表征：直观地观察页岩表面孔隙结构特征，观察纳米级孔隙，获得精细的孔隙结构图像，测量其参数；2）微观储集空间的间接定量表征：采用高压压汞法+低压氮气吸附法+低压 CO2吸附法联合表征页岩孔隙结构；
（4）页岩气赋存富集机理分析：通过等温吸附实验、脱附气测定、损失气测定、残余气测定、储层敏感性评价实验分析技术，从页岩储层微孔裂隙（特别是纳米级孔隙）特征入手，分析页岩有机质、粘土矿物纳米级孔隙结构对页岩气的赋存状态的影响机制，掌握纳米微孔中流体的行为与性质的变化规律，研究水、烃类与矿物的反应及对储集物性的影响，评价页岩的含气性，更好的评价页岩微观孔隙的吸附和扩散能力，为气体分子在孔隙中的吸附和扩散提供模拟依据；
（5）页岩储层三维建模方法：将基于 CCSIM 算法和扫描电镜图像进行三维建模和基于 AFM的三维建模方法结合起来表征页岩非均质储层三维孔隙结构特征；
（6）页岩粘土矿物表面及单元层层间结构分析：采用蒙特卡洛方法研究气体分子在粘土矿物基面和边缘面上的吸附行为，探讨表面能量特征和吸附相的结构；通过分子动力学计算和实验测量，研究有机插层粘土的层间微观结构和动力学特征，以长链三甲基烷基铵离子 HDTMA+为柱撑剂，通过 X 射线衍射、电子探针实验数据综合分析和计算，对不同条件下 HDTMA+在粘土矿物层间的排列模式的演化过程进行详细而系统的研究，揭示HDTMA+在粘土矿物层间的排列方式，确定层间孔径、孔道率、比表面积及对气体的吸附性，并采用分子动力学（Molecular Dynamics，MD）方法研究粘土矿物的层间水化和膨胀，利用分子模拟软件 materials studio为实现模拟的操作平台，根据原子的排列结构构建分子模型，基于页岩相关数据，按照粘土矿物组分分别构建粘土矿物层间结构的三维模型，将粘土矿物含量分析得到的各粘土矿物含量比例建立混合模型，使模型更趋近于现实情况的页岩粘土矿物孔隙；
（7）CO2/CH4分子在页岩微观孔隙中的吸附和扩散模拟：以页岩孔隙系统的三维等效模型为载体，以不同参数作为模拟的初始条件和约束条件，基于淬火密度泛函（QSDFT）分子模拟方法及耗散粒子动力学（DPD）方法，对不同页岩样本孔隙系统中CH4吸附、运移的动力学特征进行模拟，分析不同条件下二氧化碳对甲烷的替代效果；以各种模型建立为基础，充分运用蒙特卡洛、分子力学以及分子动力学的理论依据对流体分子进行吸附和扩散模拟，即分析页岩储层不同组分的孔隙模型在页岩储层不同埋深条件下对甲烷、二氧化碳的吸附和扩散行为，计算页岩储层不同组分的孔隙模型在不同埋深下对甲烷、二氧化碳的吸附量及吸附热，以及甲烷、二氧化碳在页岩储层不同组分的孔隙模型中的相对密度和自扩散系数；
（8）对结果的验证：对上述模拟结果，进行实验分析验证，采用含气量数据、X 衍射检验数据以及二氧化碳甲烷实验数据进行综合分析，将实验数据与模拟结果进行比对和验证。
进一步的，所述步骤（3）中微观储集空间二维形态、结构与大小的表征是利用场发射扫描电镜、环境扫描电镜、原子力显微镜方法直观地观察黏土矿物微孔孔喉二维图像特征。
进一步的，所述步骤（5）中基于CCSIM 算法和信息熵的三维建模方法，是在仔细甄别扫描电镜图像后，选出能反应页岩强非均质性的典型图像，对选出的图像以信息熵为根据进行网格化处理，将每个网格区域四等分，信息熵值也被四等分，把这个运算过程会持续到分解为足够小的信息熵值；选出的初始图像作为3D网格的第一层，前后左右四个方位用CCSIM 算法生成图像；第二层则在第一层网格化后的图像基础上，选取一部分进行运算生成；以下各层将选能反映该层非均质性的数据和图像来生成；经过多次这样的迭代运算以后，得到所有的2D图像；之后将这些图像进行堆叠便可以生成三维模型；选择多个不同倍数的页岩 SEM 图像，逐一进行建模，来表征页岩三维非均质储层。
进一步的，所述步骤（7）中，是以不同孔隙结构、储层压力、温度、气体浓度参数作为模拟的初始条件，加以孔隙大小、连通性、液态水含量及分布为约束条件。
本发明的有益效果是：
1、微纳米孔隙空间作为非常规油气重要的储集空间，其孔隙结构特征是影响储集层物性的重要因素，准确全面地表征微纳米孔隙结构已成为非常规储集层研究的重要内容。从页岩微观孔隙、矿物类型、TOC含量、渗流力学、化学热力学和化学动力学原理等相关基础理论出发，结合数值模拟研究埋藏条件下页岩中CO2/CH4的特点以及在微观多孔介质中运动的特性，建立页岩中CO2/CH4吸附与扩散模型，最终为CO2高效置换页岩中的CH4并为CO2在页岩中封存的研究提供理论基础是重要的研究方向，本发明所提供的方法综合考察了地质作用及页岩成分组成对页岩储层的物理化学条件、孔隙系统特征及 CH4/CO2分子传递过程的影响，能够构建客观实际的页岩气非均质储层微观孔隙模型。
2、在本发明所提供的方法的基础上，可获得页岩中 CO2/CH4分子传递过程模拟的初始和边界条件、约束性因素，通过分子模拟的手段可以揭示页岩储层中 CO2/CH4分子吸附和扩散机理，并对整个传递过程进行定量刻画。
3、基于CCSIM 算法和信息熵的三维建模方法能够全面、迅速并且经济地反映页岩非均质储层三维孔隙结构特征，但和页岩实际孔隙结构分析对比，发现存在准确性不足的问题；基于AFM的三维建模方法能够对样品某一个横切面表面的起伏形态、单孔深度进行定量描述，同时也可以判别空隙类型，表征准确度高，但通过基于AFM三维建模方法定量描述样品多个横切面表面孔隙结构特征的方法的分析成本高，在实际应用中受到限制。本发明通过将基于CCSIM 算法和信息熵的三维建模方法和基于 AFM的三维建模方法两种方法结合起来对页岩非均质储层三维孔隙结构特征进行表征，能够充分发挥两种方法各自的优势，不需要通过AFM三维建模法描述样品若干个横切面特征以构建三维立体模型，只需要选取1-3个特征横切面与CCSIM 算法和信息熵三维建模法相结合，通过且仅仅通过这两种三维建模方法，即可真实、经济、准确地表征页岩非均质储层三维孔隙结构特征，并且对于大尺寸的样品分析准确性大大提高，这两种三维建模方法的结合有效降低了分析成本且提高了分析准确性，在实际生产中具有很好的应用前景。
4、压汞法（MIP）的优势表征范围在 20nm 以上，主要为介孔和宏孔；N2吸附在部分微孔到 50nm 之间，位于介孔范围，N2测定数据不仅能定量表征孔径分布，还能从等温线与滞后回线类型推测孔隙形状；CO2吸附在 0.35-1.5nm，优势表征范围为微孔，采用高压压汞法+低压氮气吸附法+低压 CO2吸附法联合表征页岩孔隙结构，用气体吸附技术与高压压汞技术相结合的方式来分析页岩，能够解决气体吸附技术无法测量大孔和高压压汞技术不能测量微孔的问题。
具体实施方式
下面进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。
（1）现场资料收集。
收集野外剖面样品和昭通页岩气国家级示范区及威远-长宁页岩气国家级示范区岩心资料、随钻测试、压裂排采数据，掌握开发页岩气层在地层条件下的压力、温度、区域应力特征参数，用以页岩储层特征的表述和模拟。
（2）页岩特征分析。
在区域地质演化研究的基础上，通过薄片鉴定、X衍射、扫描电镜，结合测井资料对研究区页岩岩芯样品进行矿物成分和含量分析，研究粘土矿物中蒙脱石向伊利石转化过程；进行有机碳测定、干酪根显微组分及类型、有机质成熟度、微观古生物实验分析鉴定，确定泥页岩的有机质性质。页岩特征是分析页岩的基础，以上参数可以为页岩储层微观特征表征提供依据。
（3）页岩气非均质储层微观孔隙多尺度定量、定性表征。
页岩基质孔隙的结构复杂、类型多样、孔径分布宽，并且孔隙形成与页岩有机质和黏土有关系，因此采用单一手段确定页岩微观非均质孔隙结构比较困难，采用直接观测和气体吸附间接数值测定两种实验方法来表征页岩微观储集空间特征。
1)微观储集空间二维形态、结构与大小的表征：
利用扫描电镜（场发射扫描电镜、环境扫描电镜）、原子力显微镜方法直观地观察黏土矿物微孔孔喉二维图像特征，直观地观察页岩表面孔隙结构特征，观察纳米级孔隙，获得精细的孔隙结构图像，测量其参数。
2）微观储集空间的间接定量表征：
用气体吸附技术与高压压汞技术相结合的方式来分析页岩，能够解决气体吸附技术无法测量大孔和高压压汞技术不能测量微孔的问题，N2吸附能较为准确的反映泥页岩微孔、介孔的分布情况，而压汞法则能弥补 N2 吸附对大孔分析的不足。在运用传统孔隙结构研究方法的同时，需要借鉴物理化学、分析化学以及材料科学多项测试方法，以提高致密储层微观孔隙特征表征的精度与准确性。
由于各种表征方法都有其优越性及局限性，压汞法（MIP）的优势表征范围在 20nm以上，主要为介孔和宏孔；N2吸附在部分微孔到 50nm 之间（DFT 法可致微孔，BJH 法下限仅为2nm），位于介孔范围，N2测定数据不仅能定量表征孔径分布，还能从等温线与滞后回线类型推测孔隙形状；CO2吸附在 0.35-1.5nm，优势表征范围为微孔。因此，准确定量表征同时具微孔、介孔和宏孔的泥页岩孔隙结构时，优选各个测试结果优势范围，采用高压压汞法+低压氮气吸附法+低压 CO2吸附法联合表征页岩孔隙结构。
（4）页岩气赋存富集机理分析。
通过等温吸附实验、脱附气测定、损失气测定、残余气测定、储层敏感性评价实验分析技术，从页岩储层微孔裂隙（特别是纳米级孔隙）特征入手，分析页岩有机质、粘土矿物纳米级孔隙结构对页岩气的赋存状态的影响机制，掌握纳米微孔中流体的行为与性质的变化规律，研究水、烃类与矿物的反应及对储集物性的影响，评价页岩的含气性，更好的评价页岩微观孔隙的吸附和扩散能力，为气体分子在孔隙中的吸附和扩散提供模拟依据。
（5）页岩储层三维建模方法。
①基于 CCSIM 算法和信息熵的三维建模方法
基于 CCSIM 算法（cross-correlation-based simulation）和扫描电镜图像，在仔细甄别扫描电镜图像后，选出能反应页岩强非均质性的典型图像，对选出的图像以信息熵为根据进行网格化处理，将每个网格区域四等分，信息熵值也被四等分，把这个运算过程会持续到分解为足够小的信息熵值。
选出的初始图像作为 3D 网格的第一层，前后左右四个方位用 CCSIM 算法生成图像。第二层则在第一层网格化后的图像基础上，选取一部分进行运算生成。以下各层将选能反映该层非均质性的数据和图像来生成。经过多次这样的迭代运算以后，得到所有的 2D图像。之后将这些图像进行堆叠便可以生成三维模型。选择多个不同倍数的页岩 SEM 图像，逐一进行建模，来表征页岩三维非均质储层。
②基于 AFM的三维建模方法
AFM 不需要进行更多的其他制样处理就可以得到样品表面的三维形貌图像，并且探针与样品表面相互作用力在10-8N 以下，不存在诸如 SEM 的电子束对样品损伤的问题。基于AFM（原子力显微镜）图像在接触模式下可以有多种方法测量样品的表面性质，因此，运用横切面分析，绘制样品某一个方向的横切面，即可得到此横切面的起伏趋势，不仅可对该样品表面的起伏形态、单孔深度进行定量描述，同时，也可判别孔隙类型。
（6） 页岩粘土矿物表面及单元层层间结构分析。
采用蒙特卡洛（Monte Carlo，MC）方法研究气体分子在粘土矿物基面和边缘面上的吸附行为，探讨表面能量特征和吸附相的结构；通过分子动力学计算和实验测量，研究有机插层粘土的层间微观结构和动力学特征，以长链三甲基烷基铵离子 HDTMA+为柱撑剂，通过 X 射线衍射、电子探针实验数据综合分析和计算，对不同条件下 HDTMA+在粘土矿物层间的排列模式的演化过程进行详细而系统的研究，揭示 HDTMA+在粘土矿物层间的排列方式，确定层间孔径、孔道率、比表面积及对气体的吸附性，并采用分子动力学（MolecularDynamics，MD）方法研究粘土矿物的层间水化和膨胀。利用分子模拟软件 materialsstudio为实现模拟的操作平台，根据原子的排列结构构建分子模型，基于川南地区五峰组、龙马溪组页岩相关数据，按照粘土矿物组分分别构建粘土矿物层间结构的三维模型，将粘土矿物含量分析得到的各粘土矿物含量比例建立混合模型，使模型更趋近于现实情况的页岩粘土矿物孔隙。
（7）CO2/CH4分子在页岩微观孔隙中的吸附和扩散模拟。
以页岩孔隙系统的三维等效模型为载体，以不同孔隙结构、储层压力、温度、气体浓度（包括了 CH4、CO2、水蒸气等）参数作为模拟的初始条件，加以孔隙大小、连通性、液态水含量及分布为约束条件，基于淬火密度泛函（QSDFT）分子模拟方法及耗散粒子动力学（DPD）方法，分别对昭通国家级页岩气示范区五峰组、龙马溪组页岩微-纳米级孔隙系统中 CH4吸附、运移的动力学特征进行模拟，分析不同条件下二氧化碳对甲烷的替代效果。以各种模型建立为基础，充分运用蒙特卡洛、分子力学以及分子动力学的理论依据对流体分子（CO2分子、CH4分子）进行吸附和扩散模拟，即分析页岩储层不同组分的孔隙模型在页岩储层不同埋深条件下对甲烷、二氧化碳的吸附和扩散行为，计算页岩储层不同组分的孔隙模型在不同埋深下对甲烷、二氧化碳的吸附量及吸附热，以及甲烷、二氧化碳在页岩储层不同组分的孔隙模型中的相对密度和自扩散系数。
（8）对结果的验证
对上述模拟结果，进行实验分析验证，采用含气量数据、X 衍射检验数据以及二氧化碳甲烷实验数据进行综合分析，将实验数据与模拟结果进行比对和验证。
以上所述仅是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。

Claims (3)

1.一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法，其特征在于，包括如下步骤：(1)现场资料收集：收集野外剖面样品和岩心资料、随钻测试和压裂排采数据，掌握开发页岩气层在地层条件下的压力、温度和区域应力特征参数，用以页岩储层特征的表述和模拟；
(2)页岩特征分析：在区域地质演化研究的基础上，通过薄片鉴定、X衍射和扫描电镜，结合测井资料等对研究区页岩岩芯样品进行矿物成分和含量分析，研究粘土矿物中蒙脱石向伊利石转化过程；进行有机碳测定、干酪根显微组分及类型、有机质成熟度和微观古生物实验分析鉴定，确定泥页岩的有机质性质，为页岩储层微观特征表征提供依据；
(3)页岩气非均质储层微观孔隙多尺度定量和定性表征：采用直接观测和气体吸附间接数值测定两种实验方法来表征页岩微观储集空间特征，1)微观储集空间二维形态、结构与大小的表征：直观地观察页岩表面孔隙结构特征，观察纳米级孔隙，获得精细的孔隙结构图像，测量其参数；2)微观储集空间的间接定量表征：采用高压压汞法+低压氮气吸附法+低压CO2吸附法联合表征页岩孔隙结构；
(4)页岩气赋存富集机理分析：通过等温吸附实验、脱附气测定、损失气测定、残余气测定和储层敏感性评价实验分析技术，从页岩储层微孔裂隙特征入手，分析页岩有机质和粘土矿物纳米级孔隙结构对页岩气的赋存状态的影响机制，掌握纳米微孔中流体的行为与性质的变化规律，研究水、烃类与矿物的反应及对储集物性的影响，评价页岩的含气性，更好的评价页岩微观孔隙的吸附和扩散能力，为气体分子在孔隙中的吸附和扩散提供模拟依据；
(5)页岩储层三维建模方法：通过且仅通过将基于CCSIM算法和扫描电镜图像进行三维建模和基于AFM的三维建模方法结合起来表征页岩非均质储层三维孔隙结构特征，基于CCSIM算法和信息熵的三维建模方法，是在仔细甄别扫描电镜图像后，选出能反应页岩强非均质性的典型图像，对选出的图像以信息熵为根据进行网格化处理，将每个网格区域四等分，信息熵值也被四等分，把这个运算过程会持续到分解为足够小的信息熵值；选出的初始图像作为3D网格的第一层，前后左右四个方位用CCSIM算法生成图像；第二层则在第一层网格化后的图像基础上，选取一部分进行运算生成；以下各层将选能反映该层非均质性的数据和图像来生成；经过多次这样的迭代运算以后，得到所有的2D图像；之后将这些图像进行堆叠便可以生成三维模型；选择多个不同倍数的页岩SEM图像，逐一进行建模，来表征页岩三维非均质储层；
(6)页岩粘土矿物表面及单元层层间结构分析：采用蒙特卡洛方法研究气体分子在粘土矿物基面和边缘面上的吸附行为，探讨表面能量特征和吸附相的结构；通过分子动力学计算和实验测量，研究有机插层粘土的层间微观结构和动力学特征，以长链三甲基烷基铵离子HDTMA+为柱撑剂，通过X射线衍射、电子探针实验数据综合分析和计算，对不同条件下HDTMA+在粘土矿物层间的排列模式的演化过程进行详细而系统的研究，揭示HDTMA+在粘土矿物层间的排列方式，确定层间孔径、孔道率、比表面积及对气体的吸附性，并采用分子动力学方法研究粘土矿物的层间水化和膨胀，利用分子模拟软件materials studio为实现模拟的操作平台，根据原子的排列结构构建分子模型，基于页岩相关数据，按照粘土矿物组分分别构建粘土矿物层间结构的三维模型，将粘土矿物含量分析得到的各粘土矿物含量比例建立混合模型，使模型更趋近于现实情况的页岩粘土矿物孔隙；
(7)CO2/CH4分子在页岩微观孔隙中的吸附和扩散模拟：以页岩孔隙系统的三维等效模型为载体，以不同参数作为模拟的初始条件和约束条件，基于淬火密度泛函分子模拟方法及耗散粒子动力学方法，对不同页岩样本孔隙系统中CH4吸附和运移的动力学特征进行模拟，分析不同条件下二氧化碳对甲烷的替代效果；以各种模型建立为基础，充分运用蒙特卡洛、分子力学以及分子动力学的理论依据对流体分子进行吸附和扩散模拟，即分析页岩储层不同组分的孔隙模型在页岩储层不同埋深条件下对甲烷、二氧化碳的吸附和扩散行为，并得到相应的模拟结果：计算页岩储层不同组分的孔隙模型在不同埋深下对甲烷、二氧化碳的吸附量及吸附热，以及甲烷、二氧化碳在页岩储层不同组分的孔隙模型中的相对密度和自扩散系数；
(8)对结果的验证：对上述模拟结果，进行实验分析验证，采用含气量数据、X衍射检验数据以及二氧化碳甲烷实验数据进行综合分析，将实验数据与模拟结果进行比对和验证。
2.根据权利要求1所述的一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法，其特征在于，所述步骤(3)中微观储集空间二维形态、结构与大小的表征是利用场发射扫描电镜、环境扫描电镜和原子力显微镜方法直观地观察黏土矿物微孔孔喉二维图像特征。
3.根据权利要求1所述的一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法，其特征在于，所述步骤(7)中，是以不同孔隙结构、储层压力、温度和气体浓度参数作为模拟的初始条件，加以孔隙大小、连通性和液态水含量及分布为约束条件。
CN201911361351.4A 2019-12-25 2019-12-25 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法 Active CN111007233B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911361351.4A CN111007233B (zh) 2019-12-25 2019-12-25 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911361351.4A CN111007233B (zh) 2019-12-25 2019-12-25 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法

Publications (2)

Publication Number Publication Date
CN111007233A CN111007233A (zh) 2020-04-14
CN111007233B true CN111007233B (zh) 2022-03-11

Family

ID=70118921

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911361351.4A Active CN111007233B (zh) 2019-12-25 2019-12-25 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法

Country Status (1)

Country Link
CN (1) CN111007233B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111912958B (zh) * 2020-08-24 2021-04-16 东北石油大学 一种泥页岩无机矿物赋存油中吸附及游离油量检测方法
CN112414891B (zh) * 2020-10-27 2021-08-24 成都理工大学 一种页岩气吸附回滞现象的分子模拟方法
CN113670960B (zh) * 2021-07-23 2021-12-28 西南石油大学 一种基于分子模拟的真实页岩气吸附量预测方法
CN113820472B (zh) * 2021-09-06 2022-09-06 中国地质大学（武汉） 一种评价成烃生物对页岩气储集能力影响的方法
CN113821937B (zh) * 2021-10-29 2022-12-20 中国华能集团清洁能源技术研究院有限公司 一种定量预测二氧化碳强化气藏开采和封存的方法
CN117388070A (zh) * 2023-12-11 2024-01-12 中国矿业大学(北京) 一种气固耦合微观动态原位试验系统及试验方法
CN118150510A (zh) * 2024-05-11 2024-06-07 中国地质大学（北京） 一种有机质气体吸附能力的确定方法

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011133885A1 (en) * 2010-04-23 2011-10-27 The Board Of Regents Of The University Of Oklahoma Total storage capacity and total porosity of porous media
CA2822696A1 (en) * 2011-01-27 2012-08-02 Schlumberger Canada Limited Gas sorption analysis of unconventional rock samples
CA2897689A1 (en) * 2013-02-08 2014-08-14 Schlumberger Canada Limited Apparatus and methodology for measuring properties of microporous material at multiple scales
CN104569344A (zh) * 2014-12-31 2015-04-29 长江大学 页岩储层脆性矿物地震定量表征方法
CN104569343A (zh) * 2014-12-31 2015-04-29 长江大学 页岩储层脆性矿物测井定量表征方法
CN104698152A (zh) * 2015-02-27 2015-06-10 中国石油大学(华东) 一种研究泥页岩吸附天然气主控因素的方法
CN105651966A (zh) * 2016-01-18 2016-06-08 山东科技大学 一种页岩油气优质储层评价方法及参数确定方法
CN105717202A (zh) * 2016-01-29 2016-06-29 中国石油大学(华东) 表征页岩吸附解吸天然气过程中组分和同位素分馏作用的方法
CN106940279A (zh) * 2017-05-18 2017-07-11 中国石油大学(华东) 一种评价泥页岩储层赋存吸附气量的方法
CN107271322A (zh) * 2017-07-27 2017-10-20 中国石油大学(华东) 一种定量评价泥页岩页岩气资源量及表征吸附气和游离气转化规律的方法
CN107449707A (zh) * 2017-07-03 2017-12-08 中国石油天然气股份有限公司 页岩储层中不同尺度孔隙定量的三维表征确定方法和装置
CN108458960A (zh) * 2018-03-27 2018-08-28 中国石油大学(华东) 一种富有机质泥页岩含氢组分、孔隙度及孔径的评价方法
CN108595778A (zh) * 2018-04-02 2018-09-28 中国石油大学(北京) 一种页岩复合有机质分子模型的构建方法
CN108956424A (zh) * 2018-07-19 2018-12-07 湖南科技大学 一种页岩中孔隙定量表征的方法
CN109164501A (zh) * 2018-09-07 2019-01-08 中国石油化工股份有限公司 一种基于分形维数评价页岩油可动性的方法
US10190998B1 (en) * 2018-08-29 2019-01-29 Research Institute Of Petroleum Exploration & Development, Dagang Oil Field Of Cnpc Method and device for evaluating and predicting a shale oil enrichment areas of fault lacustrine basins
CN109444015A (zh) * 2018-10-31 2019-03-08 成都理工大学 多视域、多尺度下的页岩储层微孔系统鉴定方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9279897B2 (en) * 2005-12-12 2016-03-08 Hannes Zuercher Locating oil or gas passively by observing a porous oil and gas saturated system giving off its characteristic resonance response to ambient background noise, including optional differentiation of oil, locatinggas and water
AR097834A1 (es) * 2013-10-01 2016-04-20 Shell Int Research Métodos para estimar la densidad del recurso utilizando espectroscopia raman de inclusiones en los plays de recursos de esquisto
US10386529B2 (en) * 2014-11-19 2019-08-20 Schlumberger Technology Corporation Subsurface estimation of level of organic maturity
US10151737B2 (en) * 2015-04-08 2018-12-11 King Fahd University Of Petroleum And Minerals Method for permeability prediction of shale gas
US10151715B2 (en) * 2015-09-03 2018-12-11 Saudi Arabian Oil Company Nano-level evaluation of kerogen-rich reservoir rock
US9869649B2 (en) * 2015-09-03 2018-01-16 Saudi Arabian Oil Company Nano-level evaluation of kerogen-rich reservoir rock

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011133885A1 (en) * 2010-04-23 2011-10-27 The Board Of Regents Of The University Of Oklahoma Total storage capacity and total porosity of porous media
CA2822696A1 (en) * 2011-01-27 2012-08-02 Schlumberger Canada Limited Gas sorption analysis of unconventional rock samples
CA2897689A1 (en) * 2013-02-08 2014-08-14 Schlumberger Canada Limited Apparatus and methodology for measuring properties of microporous material at multiple scales
CN104569344A (zh) * 2014-12-31 2015-04-29 长江大学 页岩储层脆性矿物地震定量表征方法
CN104569343A (zh) * 2014-12-31 2015-04-29 长江大学 页岩储层脆性矿物测井定量表征方法
CN104698152A (zh) * 2015-02-27 2015-06-10 中国石油大学(华东) 一种研究泥页岩吸附天然气主控因素的方法
CN105651966A (zh) * 2016-01-18 2016-06-08 山东科技大学 一种页岩油气优质储层评价方法及参数确定方法
CN105717202A (zh) * 2016-01-29 2016-06-29 中国石油大学(华东) 表征页岩吸附解吸天然气过程中组分和同位素分馏作用的方法
CN106940279A (zh) * 2017-05-18 2017-07-11 中国石油大学(华东) 一种评价泥页岩储层赋存吸附气量的方法
CN107449707A (zh) * 2017-07-03 2017-12-08 中国石油天然气股份有限公司 页岩储层中不同尺度孔隙定量的三维表征确定方法和装置
CN107271322A (zh) * 2017-07-27 2017-10-20 中国石油大学(华东) 一种定量评价泥页岩页岩气资源量及表征吸附气和游离气转化规律的方法
CN108458960A (zh) * 2018-03-27 2018-08-28 中国石油大学(华东) 一种富有机质泥页岩含氢组分、孔隙度及孔径的评价方法
WO2019184429A1 (zh) * 2018-03-27 2019-10-03 中国石油大学(华东) 一种富有机质泥页岩含氢组分、孔隙度及孔径的评价方法
CN108595778A (zh) * 2018-04-02 2018-09-28 中国石油大学(北京) 一种页岩复合有机质分子模型的构建方法
CN108956424A (zh) * 2018-07-19 2018-12-07 湖南科技大学 一种页岩中孔隙定量表征的方法
US10190998B1 (en) * 2018-08-29 2019-01-29 Research Institute Of Petroleum Exploration & Development, Dagang Oil Field Of Cnpc Method and device for evaluating and predicting a shale oil enrichment areas of fault lacustrine basins
CN109164501A (zh) * 2018-09-07 2019-01-08 中国石油化工股份有限公司 一种基于分形维数评价页岩油可动性的方法
CN109444015A (zh) * 2018-10-31 2019-03-08 成都理工大学 多视域、多尺度下的页岩储层微孔系统鉴定方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A Pore Scale Gas Flow Model for Shale Gas Reservoir;Vivek Swami等;《SPE Americas Unconventional Resources Conference》;20120630;第1-16页 *
Experimental investigation and Grand Canonical Monte Carlo simulation of gas shale adsorption from the macro to the nano scale;Hamza Aljamaan等;《Journal of Natural Gas Science and Engineering》;20171231;第48卷;第119-137页 *
Three-Dimensional Stochastic Characterization of Shale SEM Images;Pejman Tahmasebi等;《Transport in Porous Media》;20150904;第110卷;第521-531页 *
致密储层孔隙结构表征技术及发展趋势;蒋裕强 等;《地质科技情报》;20140531;第33卷(第3期);第63-70页 *
页岩微观孔隙模型构建——以滇黔北地区志留系龙马溪组页岩为例;陈科洛 等;《石油勘探与开发》;20181231;第45卷(第3期);第396-405页 *

Also Published As

Publication number Publication date
CN111007233A (zh) 2020-04-14

Similar Documents

Publication Publication Date Title
CN111007233B (zh) 一种分析页岩微观孔隙中甲烷-二氧化碳运动行为的方法
Yu et al. Transport of shale gas in microporous/nanoporous media: molecular to pore-scale simulations
Wang et al. Molecular simulation of gas adsorption in shale nanopores: A critical review
Wang et al. Supercritical methane diffusion in shale nanopores: effects of pressure, mineral types, and moisture content
Wang et al. Experimental investigation of gas mass transport and diffusion coefficients in porous media with nanopores
Sun et al. Characterization of gas transport behaviors in shale gas and tight gas reservoirs by digital rock analysis
Chen et al. Pore structure characterization of different rank coals using N2 and CO2 adsorption and its effect on CH4 adsorption capacity: A case in Panguan syncline, western Guizhou, China
Tang et al. Effect of organic matter and maturity on pore size distribution and gas storage capacity in high-mature to post-mature shales
Yang et al. Review of shale gas sorption and its models
Wang et al. Pore characterization of different clay minerals and its impact on methane adsorption capacity
Li et al. Adsorbed and free oil in lacustrine nanoporous shale: A theoretical model and a case study
Babatunde et al. Molecular simulation study of CO2/CH4 adsorption on realistic heterogeneous shale surfaces
Zhou et al. Confinement effects and CO2/CH4 competitive adsorption in realistic shale kerogen nanopores
Jia et al. Variations of pore structure in organic-rich shales with different lithofacies from the Jiangdong block, fuling shale gas field, SW China: Insights into gas storage and pore evolution
Sun et al. A review on shale oil and gas characteristics and molecular dynamics simulation for the fluid behavior in shale pore
Li et al. The general form of transport diffusivity of shale gas in organic-rich nano-slits—A molecular simulation study using Darken approximation
Meng et al. Rock fabric of lacustrine shale and its influence on residual oil distribution in the Upper Cretaceous Qingshankou Formation, Songliao Basin
Rother et al. Sorption phase of supercritical CO2 in silica aerogel: experiments and mesoscale computer simulations
Wu et al. Pore characterization and inner adsorption mechanism investigation for methane in organic and inorganic matters of shale
Yang et al. Recent advances in multiscale digital rock reconstruction, flow simulation, and experiments during shale gas production
Pang et al. Evaluation of gas adsorption in nanoporous shale by simplified local density model integrated with pore structure and pore size distribution
Pan et al. Mini review on wettability in the methane–liquid–rock system at reservoir conditions: implications for gas recovery and geo-storage
Qian et al. Study on the pore structure and fractal characteristics of different lithofacies of Wufeng–Longmaxi Formation shale in southern Sichuan Basin, China
Wang et al. New coupled apparent permeability models for gas transport in inorganic nanopores of shale reservoirs considering multiple effects
Jiang et al. Review of shale gas transport prediction: Basic theory, numerical simulation, application of ai methods, and perspectives

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