CN109540764B - 评价泥页岩储层贡献孔隙内吸附态甲烷厚度和密度的方法 - Google Patents
评价泥页岩储层贡献孔隙内吸附态甲烷厚度和密度的方法 Download PDFInfo
- Publication number
- CN109540764B CN109540764B CN201811521652.4A CN201811521652A CN109540764B CN 109540764 B CN109540764 B CN 109540764B CN 201811521652 A CN201811521652 A CN 201811521652A CN 109540764 B CN109540764 B CN 109540764B
- Authority
- CN
- China
- Prior art keywords
- methane
- clay
- organic matter
- minerals
- density
- 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
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
-
- 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
Abstract
一种评价泥页岩储层有机质、粘土和其它矿物所贡献孔隙赋存吸附态甲烷厚度和密度的方法。该方法步骤：粉碎一块泥页岩样品筛选3个以上不同目数子样品进行TOC、干酪根元素含量、全岩分析、低温氮气吸附‑解吸和甲烷等温吸附实验；利用子样品TOC和干酪根元素含量计算有机质含量，并归一化有机质、粘土和其它矿物含量；根据子样品有机质、粘土、其它矿物含量和低温氮气吸附‑解吸结果评价单位质量有机质、粘土和其它矿物中各孔径孔隙体积；根据子样品有机质、粘土、其它矿物含量和甲烷等温吸附实验结果评价单位质量有机质、粘土和其它矿物赋存吸附态甲烷含量；建立模型计算有机质、粘土和其它矿物所贡献不同孔径孔隙内吸附态甲烷厚度和密度。
Description
技术领域
一种评价泥页岩储层有机质、粘土和其它矿物所贡献孔隙赋存吸附态甲烷厚度和密度的方法。
背景技术
页岩气是主体以吸附和游离状态赋存于具有生烃能力的泥岩及页岩等地层中的天然气聚集，游离气和吸附气是泥页岩储层中页岩气的两种主要赋存方式。游离气主要存在于泥页岩储层的裂缝、大孔和小孔中心空间，吸附气主要通过物理吸附作用赋存于纳米级孔隙的内表面。影响泥页岩储层吸附甲烷含量的主要因素是孔隙壁成分、孔隙体积、孔径分布、温度、压力和含水量等。不同种类多孔介质的比表面对赋存吸附气量起主导作用，吸附态甲烷密度除了受温度影响外，还受压力影响。如果能够确定泥页岩储层中有机质、粘土矿物和其它矿物所贡献不同孔径孔隙内吸附态甲烷厚度和密度随温度、压力的变化特征，即可确定泥页岩储层赋存吸附态甲烷的绝对吸附量。为评价页岩气资源、优选页岩气开发有利区和制定页岩气井开发方案具有重要指导作用。
目前，常利用甲烷液相密度或固定的密度计算绝对吸附量是不适用的。通过实验方法界定孔隙中游离气和吸附气的边界，以及分析吸附态密度是难以实现的。目前评价吸附态页岩气的密度和吸附层厚度的方法主要为分子模拟。该方法在评价吸附气密度和厚度时存在以下不足：构建的固体表面是平板型而非实际的圆弧型，即构建的狭缝型孔隙空间，忽略了圆弧形孔隙壁对气体分子的吸附势能叠加作用；所构建的固体结构比较简单，无法搭建结构复杂的有机质等模型，进而影响分子模拟评价泥页岩中有机质对甲烷的吸附；由于受到计算量的限制，分子模拟的体系尺寸和分子数量均有限；在分子模拟预测吸附作用时，密度是通过逸度而不是压力的函数来计算的，逸度的宽泛定义是真实气体所施加的蒸汽压与相应理想气体的偏差；分子模拟的结果缺少实验依据支持或者与实验分析结果误差较大，例如吸附气密度和厚度等缺少实验支持，分子模拟与实验测量的单位质量吸附质所赋存吸附气量存在较大差异。
为此，本发明通过有机碳含量、干酪根元素含量、全岩分析、低温氮气吸附-解吸和甲烷等温吸附实验结果，先建立有机质、粘土矿物和其它矿物对各孔径孔隙体积的贡献评价模型，确定单位质量有机质、粘土矿物和其它矿物对各孔径孔隙体积的贡献；然后建立有机质、粘土矿物和其它矿物赋存吸附态甲烷含量评价模型，确定单位质量有机质、粘土矿物和其它矿物赋存吸附态甲烷含量；最后建立有机质、粘土矿物和其它矿物所贡献的不同孔径孔隙赋存吸附态甲烷厚度和密度评价模型，定量评价有机质、粘土矿物和其它矿物所贡献不同孔径孔隙赋存吸附态甲烷厚度和密度随温度、压力的变化。
发明内容
本发明的目的是：提供一种评价泥页岩储层有机质、粘土和其它矿物所贡献孔隙赋存吸附态甲烷厚度和密度的方法，实现对一种评价泥页岩储层有机质、粘土矿物和其它矿物所贡献不同孔径孔隙内吸附态甲烷厚度和密度的定量评价。克服现有技术、方法无法有效对泥页岩储层有机质、粘土矿物和其它矿物所贡献不同孔径孔隙内吸附态甲烷厚度和密度开展评价的缺点。
本发明采用的技术方案是：评价泥页岩储层有机质、粘土和其它矿物所贡献不同孔径孔隙内吸附态甲烷厚度和密度的方法，其特征在于：
步骤1：粉碎一块泥页岩储层样品，筛选3个以上不同目数子样品进行有机碳含量、干酪根元素含量、全岩分析、低温氮气吸附-解吸和甲烷等温吸附实验，获得泥页岩储层各子样品中有机碳质量百分比分别为wTOC1 0、wTOC2 0、…和wTOCn 0，干酪根中碳元素质量百分比分别为wC1、wC2、…和wCn，粘土矿物质量百分比分别为w粘土1 0、w粘土2 0、…和w粘土n 0，其它矿物质量百分比分别为w其它1 0、w其它2 0、…和w其它n 0，单位质量的泥页岩储层各子样品中孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙体积为Vij，单位质量的泥页岩储层各子样品在温度为Tx、压力为Py条件下赋存吸附态甲烷含量为Qixy，其中i=1、2、…、n，是泥页岩储层子样品的编号，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，泥页岩储层子样品中有机碳、粘土矿物和其它矿物质量百分比的单位均为%，干酪根中碳元素质量百分比的单位为%，单位质量的泥页岩储层子样品各孔径范围的孔隙体积单位均为cm3/g，单位质量的泥页岩储层各子样品赋存吸附态甲烷含量的单位均为m3/t；
步骤2：利用步骤1中实验测量泥页岩储层各子样品有机碳质量百分比wTOC1 0、wTOC2 0、…和wTOCn 0，以及干酪根中碳元素质量百分比wC1、wC2、…和wCn，按照下列公式计算泥页岩储层各子样品中有机质质量百分比w有机质1 0、w有机质2 0、…和w有机质n 0，
w有机质i 0=wTOCi 0/ wCi×100%
式中，w有机质i 0是未做归一化处理的各子样品有机质质量百分比，wTOCi 0是实验测量各子样品有机碳质量百分比，wCi是实验测量各子样品干酪根中碳元素质量百分比，i=1、2、…、n，是泥页岩储层子样品的编号,泥页岩储层子样品有机碳质量百分比、干酪根中碳元素质量百分比和泥页岩储层子样品中有机质质量百分比的单位均为%；
然后按照下列公式对每个子样品有机质、粘土矿物和其它矿物的质量百分比进行归一化处理，让每个子样品有机质、粘土矿物和其它矿物的质量百分比之和等于100%，获得各子样品归一化之后有机质质量百分比w有机质i、粘土矿物质量百分比w粘土i，和其它矿物质量百分比w其它i,
w有机质i=w有机质i 0×100%
w粘土i=w粘土i 0×(100-w有机质i 0)/100%
w其它i=w其它i 0×(100-w有机质i 0)/100%
式中，w有机质i、w粘土i和w其它i是做归一化之后各子样品有机质、粘土矿物和其它矿物质量百分比，w有机质i 0、w粘土i 0和w其它i 0是做归一化之前各子样品有机质、粘土矿物和其它矿物质量百分比，i=1、2、…、n，是泥页岩储层子样品的编号，有机质、粘土矿物和其它矿物质量百分比的单位均为%；
步骤3：根据步骤2中获得归一化之后的泥页岩各子样品中有机质质量百分比w有机质1、w有机质2、…和w有机质n，粘土矿物质量百分比w粘土1、w粘土2、…和w粘土n，其它矿物质量百分比w其它1、w其它2、…和w其它n，以及步骤1中获得单位质量各子样品中孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙体积Vij，建立如下方程组和目标函数，当目标函数f(V有机质j，V粘土j，V其它j)的值最小时，即可获得单位质量有机质、粘土矿物和其它矿物所贡献孔径范围编号为j的孔隙体积V有机质j、V粘土j和V其它j，
式中，V有机质j、V粘土j和V其它j分别是单位质量有机质、粘土矿物和其它矿物所贡献孔径范围编号为j的孔隙体积，j=1、2、…、7，是孔径范围由小到大的编号，i=1、2、…、n，是泥页岩储层子样品的编号，泥页岩储层子样品中有机质、粘土矿物和其它矿物的质量百分比单位均为%，单位质量泥页岩子样品中各孔径范围的孔隙体积的单位均为cm3/g，单位质量有机质、粘土矿物和其它矿物所贡献各孔径范围的孔隙体积的单位均为cm3/g；
步骤4：根据步骤2中获得归一化之后的泥页岩各子样品中有机质质量百分比w有机质1、w有机质2、…和w有机质n，粘土矿物质量百分比w粘土1、w粘土2、…和w粘土n，其它矿物质量百分比w其它1、w其它2、…和w其它n，以及步骤1中获得单位质量的泥页岩各子样品在温度为Tx、压力为Py条件下赋存吸附态甲烷含量Qixy，建立如下方程组和目标函数，当目标函数f(Q有机质xy，Q粘土xy，Q其它xy)的值最小时，即可获得单位质量有机质、粘土矿物和其它矿物在温度为Tx、压力为Py条件下赋存吸附态甲烷含量Q有机质xy、Q粘土xy和Q其它xy，
式中，Q有机质xy、Q粘土xy和Q其它xy分别是在温度为Tx、压力为Py条件下单位质量的有机质、粘土矿物和其它矿物赋存吸附态甲烷含量，Qixy是在温度为Tx、压力为Py条件下单位质量且编号为i的子样品赋存吸附态甲烷含量， i=1、2、…、n，是泥页岩储层子样品的编号， x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，单位质量的泥页岩储层子样品、有机质、粘土矿物和其它矿物赋存吸附态甲烷含量的单位均为m3/t，温度Tx的单位为℃，压力Py的单位为MPa；
步骤5：在步骤3的基础上，将单位质量有机质所贡献的孔隙近似为相应孔径的圆柱体，在温度为Tx、压力为Py条件下单位质量有机质贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积V吸有机质jxy由下列公式表达，该公式显示有机质贡献孔径范围编号为j的孔隙孔径D有机质j在小于0.38nm, D有机质j不大于2倍的吸附态甲烷厚度h吸有机质jxy且D有机质j不小于0.38nm,和 D有机质j大于2倍的吸附态甲烷厚度h吸有机质jxy且D有机质j不小于0.38nm三种情况下V吸有机质jxy的值，
式中，V吸有机质jxy为温度为Tx、压力为Py条件下单位质量有机质贡献孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积，V有机质j为单位质量有机质贡献孔径范围编号为j的孔隙体积，D有机质j为有机质所贡献孔径范围编号为j的孔隙孔径，h吸有机质jxy为有机质贡献孔径范围编号为j的孔隙内吸附态甲烷厚度，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位为cm3/g，单位质量有机质所贡献孔隙体积的单位为cm3/g，孔隙孔径和吸附态甲烷厚度的单位均为nm；
在步骤3的基础上，将单位质量粘土矿物所贡献的孔隙近似为相应孔径的圆柱体，在温度为Tx、压力为Py条件下单位质量粘土矿物贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积V吸粘土jxy由下列公式表达，该公式显示粘土矿物所贡献的孔径范围编号为j的孔隙孔径D粘土j在小于0.38nm, D粘土j不大于2倍的吸附态甲烷厚度h吸粘土jxy且D粘土j不小于0.38nm,和 D粘土j大于2倍的吸附态甲烷厚度h吸粘土jxy且D粘土j不小于0.38nm三种情况下V吸粘土jxy的值，
式中，V吸粘土jxy为温度为Tx、压力为Py条件下单位质量粘土矿物贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积，V粘土j为单位质量粘土矿物所贡献的孔径范围编号为j的孔隙体积，D粘土j为粘土矿物所贡献的孔径范围编号为j的孔隙孔径，h吸粘土jxy为粘土矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位为cm3/g，单位质量粘土矿物所贡献孔隙体积的单位为cm3/g，孔隙孔径和吸附态甲烷厚度的单位均为nm；
在步骤3的基础上，将单位质量其它矿物所贡献的孔隙近似为相应孔径的圆柱体，在温度为Tx、压力为Py条件下单位质量其它矿物贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积V吸其它jxy由下列公式表达，该公式显示其它矿物所贡献的孔径范围编号为j的孔隙孔径D其它j在小于0.38nm, D其它j不大于2倍吸附态甲烷厚度h吸其它jxy且D其它j不小于0.38nm,和 D其它j大于2倍吸附态甲烷厚度h吸其它jxy且D其它j不小于0.38nm三种情况下V吸其它jxy的值，
式中，V吸其它jxy为温度为Tx、压力为Py条件下单位质量其它矿物贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，V其它j为单位质量其它矿物所贡献的孔径范围编号为j的孔隙体积，D其它j为其它矿物所贡献的孔径范围编号为j的孔隙孔径，h吸其它jxy为其它矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位为cm3/g，单位质量其它矿物所贡献孔隙体积的单位为cm3/g，孔隙孔径和吸附态甲烷厚度的单位均为nm；
步骤6：在步骤4和步骤5的基础上，根据吸附态甲烷密度小于固态甲烷密度，吸附态甲烷密度大于游离态甲烷密度，有机质所贡献孔隙内吸附态甲烷密度随孔隙孔径的增大而减小，吸附态甲烷密度随温度的升高而减小，吸附态甲烷密度随压力增大而增大的规律，建立如下方程组和目标函数，当目标函数f(ρ吸有机质jxy，h吸有机质jxy)的值最小时，即可获得温度为Tx、压力为Py条件下由有机质贡献孔径范围编号为j的孔隙内吸附态甲烷的密度ρ吸有机质jxy和厚度h吸有机质jxy，
式中， V吸有机质jxy为温度为Tx、压力为Py条件下孔径为有机质贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，ρ吸有机质jxy为温度为Tx、压力为Py条件下有机质贡献孔径范围编号为j的孔隙内吸附态甲烷密度，ρ游xy为温度为Tx、压力为Py条件下游离态甲烷密度，Q有机质xy为单位质量有机质赋存吸附态甲烷含量，ρ固为固体甲烷密度，h吸有机质jxy为温度为Tx、压力为Py条件下有机质所贡献孔径范围编号为j的孔隙内吸附态甲烷厚度，M为甲烷的摩尔质量，其值为16.0425 g/mol，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位是cm3/g，吸附态甲烷密度、游离态甲烷密度和固体甲烷密度的单位均为kg/m3，单位质量有机质赋存吸附态甲烷含量的单位是m3/t，吸附态甲烷厚度的单位是nm;
在步骤4和步骤5的基础上，根据吸附态甲烷密度小于固态甲烷密度，吸附态甲烷密度大于游离态甲烷密度，粘土矿物所贡献孔隙内吸附态甲烷密度随孔径的增大而减小，吸附态甲烷密度随温度的升高而减小，吸附态甲烷密度随压力增大而增大的规律，建立如下方程组和目标函数，当目标函数f(ρ吸粘土jxy，h吸粘土jxy)的值最小时，即可获得温度为Tx、压力为Py条件下由粘土矿物贡献孔径范围编号为j的孔隙内吸附态甲烷的密度ρ吸粘土jxy和厚度h吸粘土jxy，
式中， V吸粘土jxy为温度为Tx、压力为Py条件下粘土矿物贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，ρ吸粘土jxy为温度为Tx、压力为Py条件下粘土矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷密度，ρ游xy为温度为Tx、压力为Py条件下游离态甲烷密度，Q粘土xy为单位质量粘土矿物赋存吸附态甲烷含量，ρ固为固体甲烷密度，h吸粘土jxy为温度为Tx、压力为Py条件下粘土矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，M为甲烷的摩尔质量，其值为16.0425 g/mol，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位是cm3/g，吸附态甲烷密度、游离态甲烷密度和固体甲烷密度的单位均为kg/m3，单位质量粘土矿物赋存吸附态甲烷含量的单位是m3/t，吸附态甲烷厚度的单位是nm;
在步骤4和步骤5的基础上，根据吸附态甲烷密度小于固态甲烷密度，吸附态甲烷密度大于游离态甲烷密度，其它矿物所贡献孔隙内吸附态甲烷密度随孔径增大而减小，吸附态甲烷密度随温度升高而减小，吸附态甲烷密度随压力增大而增大的规律，建立如下方程组和目标函数，当目标函数f(ρ吸其它jxy，h吸其它jxy)的值最小时，即可获得温度为Tx、压力为Py条件下由其它矿物贡献孔径范围编号为j的孔隙内吸附态甲烷的密度ρ吸其它jxy和厚度h吸其它jxy，
式中， V吸其它jxy为温度为Tx、压力为Py条件下其它矿物贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，ρ吸其它jxy为温度为Tx、压力为Py条件下其它矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷密度，ρ游xy为温度为Tx、压力为Py条件下游离态甲烷密度，Q其它xy为单位质量其它矿物赋存吸附态甲烷含量，ρ固为固体甲烷密度，h吸其它jxy为温度为Tx、压力为Py条件下其它矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，M为甲烷的摩尔质量，其值为16.0425 g/mol，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位是cm3/g，吸附态甲烷密度、游离态甲烷密度和固体甲烷密度的单位均为kg/m3，单位质量其它矿物赋存吸附态甲烷含量的单位是m3/t，吸附态甲烷厚度的单位是nm。
附图说明
图1是本发明的流程图。
具体实施方式：
实施例1：如图1所述，一种评价泥页岩储层有机质、粘土和其它矿物所贡献孔隙赋存吸附态甲烷厚度和密度的方法，含有以下步骤；
步骤1：粉碎一块泥页岩储层样品，筛选20-40目、40-60目、60-80目、80-100目和100-120目5种目数子样品进行TOC含量、干酪根元素含量、全岩分析、低温氮气吸附-解吸和甲烷等温吸附实验，获得泥页岩储层各子样品中TOC质量百分比分别为1.28%、1.10%、2.07%、2.22%和2.94%，干酪根中碳元素质量百分比分别为86.12%、86.72%、87.01%、85.57%和87.98%，粘土矿物质量百分比分别为41.6%、42.2%、23.0%、25.7%和30.3%，其它矿物质量百分比分别为58.4%、57.8%、77.0%、74.3%和69.7%，低温氮气吸附-解吸实验获得单位质量的5个子样品中孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙体积Vij，结果见表1，甲烷等温吸附实验获得单位质量的泥页岩储层各子样品在温度为30℃、压力为1MPa、2MPa、3MPa、4MPa、5MPa、6MPa、7MPa、8MPa、9MPa和10MPa条件下赋存吸附态甲烷含量为Qixy，结果见表2，泥页岩储层子样品中有机碳、粘土矿物和其它矿物质量百分比的单位均为%，干酪根中碳元素质量百分比的单位为%，单位质量的泥页岩储层子样品各孔径范围的孔隙体积单位均为cm3/g，单位质量的泥页岩储层各子样品赋存吸附态甲烷含量的单位均为m3/t。
表1
表2
步骤2：利用5个子样品TOC质量百分比1.28%、1.10%、2.07%、2.22%和2.94%，以及干酪根中碳元素质量百分比86.12%、86.72%、87.01%、85.57%和87.98%，按照下列公式计算泥页岩储层5个子样品中未做归一化处理的有机质质量百分比分别为1.49%、1.27%、2.38%、2.59%和3.34%。
w有机质i 0=wTOCi 0/ wCi×100%
式中，w有机质i 0是未做归一化处理的各子样品有机质质量百分比，wTOCi 0是实验测量各子样品有机碳质量百分比，wCi是实验测量各子样品干酪根中碳元素质量百分比，i=1、2、…、n，是泥页岩储层子样品的编号,泥页岩储层子样品有机碳质量百分比、干酪根中碳元素质量百分比和泥页岩储层子样品中有机质质量百分比的单位均为%。
然后按照下列公式对每个子样品有机质、粘土矿物和其它矿物的质量百分比进行归一化处理，让每个子样品有机质、粘土矿物和其它矿物的质量百分比之和等于100%，获得5个子样品归一化之后有机质质量百分比分别为1.49%、1.27%、2.38%、2.59%和3.34%，粘土矿物质量百分比分别为40.98%、41.67%、22.47%、25.05%和29.32%，其它矿物质量百分比分别为57.53%、57.08%、75.21%、72.42%和67.45%。
w有机质i=w有机质i 0×100%
w粘土i=w粘土i 0×(100-w有机质i 0)/100%
w其它i=w其它i 0×(100-w有机质i 0)/100%
式中，w有机质i、w粘土i和w其它i是做归一化之后各子样品有机质、粘土矿物和其它矿物质量百分比，w有机质i 0、w粘土i 0和w其它i 0是做归一化之前各子样品有机质、粘土矿物和其它矿物质量百分比，i=1、2、…、n，是泥页岩储层子样品的编号，有机质、粘土矿物和其它矿物质量百分比的单位均为%。
步骤3：根据步骤2中获得归一化之后的泥页岩各子样品中有机质质量百分比分别为1.49%、1.27%、2.38%、2.59%和3.34%，粘土矿物质量百分比分别为40.98%、41.67%、22.47%、25.05%和29.32%，其它矿物质量百分比分别为57.53%、57.08%、75.21%、72.42%和67.45%，以及步骤1中获得单位质量各子样品中孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙体积Vij（见表1），建立如下方程组和目标函数，当目标函数f(V有机质j，V粘土j，V其它j)的值最小时，即可获得单位质量有机质、粘土矿物和其它矿物所贡献孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙体积，结果见表3。
式中，V有机质j、V粘土j和V其它j分别是单位质量有机质、粘土矿物和其它矿物所贡献孔径范围编号为j的孔隙体积，j=1、2、…、7，是孔径范围由小到大的编号，i=1、2、…、5，是泥页岩储层子样品的编号，泥页岩储层子样品中有机质、粘土矿物和其它矿物的质量百分比单位均为%，单位质量泥页岩子样品中各孔径范围的孔隙体积的单位均为cm3/g，单位质量有机质、粘土矿物和其它矿物所贡献各孔径范围的孔隙体积的单位均为cm3/g。
表3
步骤4：根据步骤2中获得归一化之后的泥页岩各子样品中有机质质量百分比分别为1.49%、1.27%、2.38%、2.59%和3.34%，粘土矿物质量百分比分别为40.98%、41.67%、22.47%、25.05%和29.32%，其它矿物质量百分比分别为57.53%、57.08%、75.21%、72.42%和67.45%，以及步骤1中获得单位质量的泥页岩各子样品在温度为30℃、压力为1MPa、2MPa、3MPa、4MPa、5MPa、6MPa、7MPa、8MPa、9MPa和10MPa条件下赋存吸附态甲烷含量Qixy(见表2)，建立如下方程组和目标函数，当目标函数f(Q有机质xy，Q粘土xy，Q其它xy)的值最小时，获得单位质量有机质、粘土矿物和其它矿物在温度为30℃、压力为1MPa、2MPa、3MPa、4MPa、5MPa、6MPa、7MPa、8MPa、9MPa和10MPa条件下赋存吸附态甲烷含量Q有机质xy、Q粘土xy和Q其它xy，结果见表4。
式中，Q有机质xy、Q粘土xy和Q其它xy分别是在温度为Tx、压力为Py条件下单位质量的有机质、粘土矿物和其它矿物赋存吸附态甲烷含量，Qixy是在温度为Tx、压力为Py条件下单位质量且编号为i的子样品赋存吸附态甲烷含量， i=1、2、…、n，n=5，是泥页岩储层子样品的编号， x=1，是温度的编号，y=1、2、…、10，是压力由低至高的编号，单位质量的泥页岩储层子样品、有机质、粘土矿物和其它矿物赋存吸附态甲烷含量的单位均为m3/t，温度Tx的单位为℃，压力Py的单位为MPa。
表4
步骤5：在步骤3的基础上，将单位质量有机质所贡献的孔隙近似为相应孔径的圆柱体，在温度为Tx、压力为Py条件下单位质量有机质贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积V吸有机质jxy由下列公式表达，该公式显示有机质贡献孔径范围编号为j的孔隙孔径D有机质j在小于0.38nm, D有机质j不大于2倍的吸附态甲烷厚度h吸有机质jxy且D有机质j不小于0.38nm,和 D有机质j大于2倍的吸附态甲烷厚度h吸有机质jxy且D有机质j不小于0.38nm三种情况下V吸有机质jxy的值。
式中，V吸有机质jxy为温度为Tx、压力为Py条件下单位质量有机质贡献孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积，V有机质j为单位质量有机质贡献孔径范围编号为j的孔隙体积，D有机质j为有机质所贡献孔径范围编号为j的孔隙孔径，h吸有机质jxy为有机质贡献孔径范围编号为j的孔隙内吸附态甲烷厚度，j=1、2、…、7，是孔径范围由小到大的编号，x=1，是温度的编号，y=1、2、…、10，是压力由低至高的编号，吸附态甲烷所占孔体积的单位为cm3/g，单位质量有机质所贡献孔隙体积的单位为cm3/g，孔隙孔径和吸附态甲烷厚度的单位均为nm。
在步骤3的基础上，将单位质量粘土矿物所贡献的孔隙近似为相应孔径的圆柱体，在温度为Tx、压力为Py条件下单位质量粘土矿物贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积V吸粘土jxy由下列公式表达，该公式显示粘土矿物所贡献的孔径范围编号为j的孔隙孔径D粘土j在小于0.38nm, D粘土j不大于2倍的吸附态甲烷厚度h吸粘土jxy且D粘土j不小于0.38nm,和 D粘土j大于2倍的吸附态甲烷厚度h吸粘土jxy且D粘土j不小于0.38nm三种情况下V吸粘土jxy的值。
式中，V吸粘土jxy为温度为Tx、压力为Py条件下单位质量粘土矿物贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积，V粘土j为单位质量粘土矿物所贡献的孔径范围编号为j的孔隙体积，D粘土j为粘土矿物所贡献的孔径范围编号为j的孔隙孔径，h吸粘土jxy为粘土矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，j=1、2、…、7，是孔径范围由小到大的编号，x=1，是温度的编号，y=1、2、…、10，是压力由低至高的编号，吸附态甲烷所占孔体积的单位为cm3/g，单位质量粘土矿物所贡献孔隙体积的单位为cm3/g，孔隙孔径和吸附态甲烷厚度的单位均为nm。
在步骤3的基础上，将单位质量其它矿物所贡献的孔隙近似为相应孔径的圆柱体，在温度为Tx、压力为Py条件下单位质量其它矿物贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积V吸其它jxy由下列公式表达，该公式显示其它矿物所贡献的孔径范围编号为j的孔隙孔径D其它j在小于0.38nm, D其它j不大于2倍吸附态甲烷厚度h吸其它jxy且D其它j不小于0.38nm,和 D其它j大于2倍吸附态甲烷厚度h吸其它jxy且D其它j不小于0.38nm三种情况下V吸其它jxy的值。
式中，V吸其它jxy为温度为Tx、压力为Py条件下单位质量其它矿物贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，V其它j为单位质量其它矿物所贡献的孔径范围编号为j的孔隙体积，D其它j为其它矿物所贡献的孔径范围编号为j的孔隙孔径，h吸其它jxy为其它矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，j=1、2、…、7，是孔径范围由小到大的编号，x=1是温度的编号，y=1、2、…、10，是压力由低至高的编号，吸附态甲烷所占孔体积的单位为cm3/g，单位质量其它矿物所贡献孔隙体积的单位为cm3/g，孔隙孔径和吸附态甲烷厚度的单位均为nm。
步骤6：在步骤4和步骤5的基础上，根据吸附态甲烷密度小于固态甲烷密度，吸附态甲烷密度大于游离态甲烷密度，有机质所贡献孔隙内吸附态甲烷密度随孔隙孔径的增大而减小，吸附态甲烷密度随温度的升高而减小，吸附态甲烷密度随压力增大而增大的规律，建立如下方程组和目标函数，当目标函数f(ρ吸有机质jxy，h吸有机质jxy)的值最小时，获得在温度为30℃、压力为1MPa、2MPa、3MPa、4MPa、5MPa、6MPa、7MPa、8MPa、9MPa和10MPa条件下由有机质贡献孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙内吸附态甲烷的密度和厚度，结果见表5和表6。
式中， V吸有机质jxy为温度为Tx、压力为Py条件下孔径为有机质贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，ρ吸有机质jxy为温度为Tx、压力为Py条件下有机质贡献孔径范围编号为j的孔隙内吸附态甲烷密度，ρ游xy为温度为Tx、压力为Py条件下游离态甲烷密度，Q有机质xy为单位质量有机质赋存吸附态甲烷含量，ρ固为固体甲烷密度，h吸有机质jxy为温度为Tx、压力为Py条件下有机质所贡献孔径范围编号为j的孔隙内吸附态甲烷厚度，M为甲烷的摩尔质量，其值为16.0425 g/mol，j=1、2、…、7，是孔径范围由小到大的编号，x=1，是温度的编号，y=1、2、…、10，是压力由低至高的编号，吸附态甲烷所占孔体积的单位是cm3/g，吸附态甲烷密度、游离态甲烷密度和固体甲烷密度的单位均为kg/m3，单位质量有机质赋存吸附态甲烷含量的单位是m3/t，吸附态甲烷厚度的单位是nm。
在步骤4和步骤5的基础上，根据吸附态甲烷密度小于固态甲烷密度，吸附态甲烷密度大于游离态甲烷密度，粘土矿物所贡献孔隙内吸附态甲烷密度随孔径的增大而减小，吸附态甲烷密度随温度的升高而减小，吸附态甲烷密度随压力增大而增大的规律，建立如下方程组和目标函数，当目标函数f(ρ吸粘土jxy，h吸粘土jxy)的值最小时，即可获得在温度为30℃、压力为1MPa、2MPa、3MPa、4MPa、5MPa、6MPa、7MPa、8MPa、9MPa和10MPa条件下由粘土矿物贡献孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙内吸附态甲烷的密度和厚度，结果见表5和表6。
式中， V吸粘土jxy为温度为Tx、压力为Py条件下粘土矿物贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，ρ吸粘土jxy为温度为Tx、压力为Py条件下粘土矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷密度，ρ游xy为温度为Tx、压力为Py条件下游离态甲烷密度，Q粘土xy为单位质量粘土矿物赋存吸附态甲烷含量，ρ固为固体甲烷密度，h吸粘土jxy为温度为Tx、压力为Py条件下粘土矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，M为甲烷的摩尔质量，其值为16.0425 g/mol，j=1、2、…、7，是孔径范围由小到大的编号，x=1，是温度的编号，y=1、2、…、10，是压力由低至高的编号，吸附态甲烷所占孔体积的单位是cm3/g，吸附态甲烷密度、游离态甲烷密度和固体甲烷密度的单位均为kg/m3，单位质量粘土矿物赋存吸附态甲烷含量的单位是m3/t，吸附态甲烷厚度的单位是nm。
在步骤4和步骤5的基础上，根据吸附态甲烷密度小于固态甲烷密度，吸附态甲烷密度大于游离态甲烷密度，其它矿物所贡献孔隙内吸附态甲烷密度随孔径增大而减小，吸附态甲烷密度随温度升高而减小，吸附态甲烷密度随压力增大而增大的规律，建立如下方程组和目标函数，当目标函数f(ρ吸其它jxy，h吸其它jxy)的值最小时，即可获得在温度为30℃、压力为1MPa、2MPa、3MPa、4MPa、5MPa、6MPa、7MPa、8MPa、9MPa和10MPa条件下由其它矿物贡献孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙内吸附态甲烷的密度和厚度，结果见表5和表6。
式中， V吸其它jxy为温度为Tx、压力为Py条件下其它矿物贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，ρ吸其它jxy为温度为Tx、压力为Py条件下其它矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷密度，ρ游xy为温度为Tx、压力为Py条件下游离态甲烷密度，Q其它xy为单位质量其它矿物赋存吸附态甲烷含量，ρ固为固体甲烷密度，h吸其它jxy为温度为Tx、压力为Py条件下其它矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，M为甲烷的摩尔质量，其值为16.0425 g/mol，j=1、2、…、7，是孔径范围由小到大的编号，x=1，是温度的编号，y=1、2、…、10，是压力由低至高的编号，吸附态甲烷所占孔体积的单位是cm3/g，吸附态甲烷密度、游离态甲烷密度和固体甲烷密度的单位均为kg/m3，单位质量其它矿物赋存吸附态甲烷含量的单位是m3/t，吸附态甲烷厚度的单位是nm。
表5
表6
Claims (1)
1.一种评价泥页岩储层有机质、粘土和其它矿物所贡献孔隙赋存吸附态甲烷厚度和密度的方法，其特征在于：
步骤1：粉碎一块泥页岩储层样品，筛选3个以上不同目数子样品进行有机碳含量、干酪根元素含量、全岩分析、低温氮气吸附-解吸和甲烷等温吸附实验，获得泥页岩储层各子样品中有机碳质量百分比分别为wTOC1 0、wTOC2 0、…和wTOCn 0，干酪根中碳元素质量百分比分别为wC1、wC2、…和wCn，粘土矿物质量百分比分别为w粘土1 0、w粘土2 0、…和w粘土n 0，其它矿物质量百分比分别为w其它1 0、w其它2 0、…和w其它n 0，单位质量的泥页岩储层各子样品中孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙体积为Vij，单位质量的泥页岩储层各子样品在温度为Tx、压力为Py条件下赋存吸附态甲烷含量为Qixy，其中i=1、2、…、n，是泥页岩储层子样品的编号，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，泥页岩储层子样品中有机碳、粘土矿物和其它矿物质量百分比的单位均为%，干酪根中碳元素质量百分比的单位为%，单位质量的泥页岩储层子样品各孔径范围的孔隙体积单位均为cm3/g，单位质量的泥页岩储层各子样品赋存吸附态甲烷含量的单位均为m3/t；
步骤2：利用步骤1中实验测量泥页岩储层各子样品有机碳质量百分比wTOC1 0、wTOC2 0、…和wTOCn 0，以及干酪根中碳元素质量百分比wC1、wC2、…和wCn，按照下列公式计算泥页岩储层各子样品中有机质质量百分比w有机质1 0、w有机质2 0、…和w有机质n 0，
w有机质i 0=wTOCi 0/ wCi×100%
式中，w有机质i 0是未做归一化处理的各子样品有机质质量百分比，wTOCi 0是实验测量各子样品有机碳质量百分比，wCi是实验测量各子样品干酪根中碳元素质量百分比，i=1、2、…、n，是泥页岩储层子样品的编号,泥页岩储层子样品有机碳质量百分比、干酪根中碳元素质量百分比和泥页岩储层子样品中有机质质量百分比的单位均为%；
然后按照下列公式对每个子样品有机质、粘土矿物和其它矿物的质量百分比进行归一化处理，让每个子样品有机质、粘土矿物和其它矿物的质量百分比之和等于100%，获得各子样品归一化之后有机质质量百分比w有机质i、粘土矿物质量百分比w粘土i，和其它矿物质量百分比w其它i,
w有机质i=w有机质i 0×100%
w粘土i=w粘土i 0×(100-w有机质i 0)/100%
w其它i=w其它i 0×(100-w有机质i 0)/100%
式中，w有机质i、w粘土i和w其它i是做归一化之后各子样品有机质、粘土矿物和其它矿物质量百分比，w有机质i 0、w粘土i 0和w其它i 0是做归一化之前各子样品有机质、粘土矿物和其它矿物质量百分比，i=1、2、…、n，是泥页岩储层子样品的编号，有机质、粘土矿物和其它矿物质量百分比的单位均为%；
步骤3：根据步骤2中获得归一化之后的泥页岩各子样品中有机质质量百分比w有机质1、w有机质2、…和w有机质n，粘土矿物质量百分比w粘土1、w粘土2、…和w粘土n，其它矿物质量百分比w其它1、w其它2、…和w其它n，以及步骤1中获得单位质量各子样品中孔径范围分别为<2nm、2-5 nm、5-10nm、10-20nm、20-50nm、50-100nm和100-200nm的孔隙体积Vij，建立如下方程组和目标函数，当目标函数f(V有机质j，V粘土j，V其它j)的值最小时，即可获得单位质量有机质、粘土矿物和其它矿物所贡献孔径范围编号为j的孔隙体积V有机质j、V粘土j和V其它j，
式中，V有机质j、V粘土j和V其它j分别是单位质量有机质、粘土矿物和其它矿物所贡献孔径范围编号为j的孔隙体积，j=1、2、…、7，是孔径范围由小到大的编号，i=1、2、…、n，是泥页岩储层子样品的编号，泥页岩储层子样品中有机质、粘土矿物和其它矿物的质量百分比单位均为%，单位质量泥页岩子样品中各孔径范围的孔隙体积的单位均为cm3/g，单位质量有机质、粘土矿物和其它矿物所贡献各孔径范围的孔隙体积的单位均为cm3/g；
步骤4：根据步骤2中获得归一化之后的泥页岩各子样品中有机质质量百分比w有机质1、w有机质2、…和w有机质n，粘土矿物质量百分比w粘土1、w粘土2、…和w粘土n，其它矿物质量百分比w其它1、w其它2、…和w其它n，以及步骤1中获得单位质量的泥页岩各子样品在温度为Tx、压力为Py条件下赋存吸附态甲烷含量Qixy，建立如下方程组和目标函数，当目标函数f(Q有机质xy，Q粘土xy，Q其它xy)的值最小时，即可获得单位质量有机质、粘土矿物和其它矿物在温度为Tx、压力为Py条件下赋存吸附态甲烷含量Q有机质xy、Q粘土xy和Q其它xy，
式中，Q有机质xy、Q粘土xy和Q其它xy分别是在温度为Tx、压力为Py条件下单位质量的有机质、粘土矿物和其它矿物赋存吸附态甲烷含量，Qixy是在温度为Tx、压力为Py条件下单位质量且编号为i的子样品赋存吸附态甲烷含量， i=1、2、…、n，是泥页岩储层子样品的编号， x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，单位质量的泥页岩储层子样品、有机质、粘土矿物和其它矿物赋存吸附态甲烷含量的单位均为m3/t，温度Tx的单位为℃，压力Py的单位为MPa；
步骤5：在步骤3的基础上，将单位质量有机质所贡献的孔隙近似为相应孔径的圆柱体，在温度为Tx、压力为Py条件下单位质量有机质贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积V吸有机质jxy由下列公式表达，该公式显示有机质贡献孔径范围编号为j的孔隙孔径D有机质j在小于0.38nm, D有机质j不大于2倍的吸附态甲烷厚度h吸有机质jxy且D有机质j不小于0.38nm,和 D有机质j大于2倍的吸附态甲烷厚度h吸有机质jxy且D有机质j不小于0.38nm三种情况下V吸有机质jxy的值，
式中，V吸有机质jxy为温度为Tx、压力为Py条件下单位质量有机质贡献孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积，V有机质j为单位质量有机质贡献孔径范围编号为j的孔隙体积，D有机质j为有机质所贡献孔径范围编号为j的孔隙孔径，h吸有机质jxy为有机质贡献孔径范围编号为j的孔隙内吸附态甲烷厚度，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位为cm3/g，单位质量有机质所贡献孔隙体积的单位为cm3/g，孔隙孔径和吸附态甲烷厚度的单位均为nm；
在步骤3的基础上，将单位质量粘土矿物所贡献的孔隙近似为相应孔径的圆柱体，在温度为Tx、压力为Py条件下单位质量粘土矿物贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积V吸粘土jxy由下列公式表达，该公式显示粘土矿物所贡献的孔径范围编号为j的孔隙孔径D粘土j在小于0.38nm, D粘土j不大于2倍的吸附态甲烷厚度h吸粘土jxy且D粘土j不小于0.38nm,和 D粘土j大于2倍的吸附态甲烷厚度h吸粘土jxy且D粘土j不小于0.38nm三种情况下V吸粘土jxy的值，
式中，V吸粘土jxy为温度为Tx、压力为Py条件下单位质量粘土矿物贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积，V粘土j为单位质量粘土矿物所贡献的孔径范围编号为j的孔隙体积，D粘土j为粘土矿物所贡献的孔径范围编号为j的孔隙孔径，h吸粘土jxy为粘土矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位为cm3/g，单位质量粘土矿物所贡献孔隙体积的单位为cm3/g，孔隙孔径和吸附态甲烷厚度的单位均为nm；
在步骤3的基础上，将单位质量其它矿物所贡献的孔隙近似为相应孔径的圆柱体，在温度为Tx、压力为Py条件下单位质量其它矿物贡献的孔径范围编号为j的孔隙内赋存吸附态甲烷所占孔体积V吸其它jxy由下列公式表达，该公式显示其它矿物所贡献的孔径范围编号为j的孔隙孔径D其它j在小于0.38nm, D其它j不大于2倍吸附态甲烷厚度h吸其它jxy且D其它j不小于0.38nm,和 D其它j大于2倍吸附态甲烷厚度h吸其它jxy且D其它j不小于0.38nm三种情况下V吸其它jxy的值，
式中，V吸其它jxy为温度为Tx、压力为Py条件下单位质量其它矿物贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，V其它j为单位质量其它矿物所贡献的孔径范围编号为j的孔隙体积，D其它j为其它矿物所贡献的孔径范围编号为j的孔隙孔径，h吸其它jxy为其它矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位为cm3/g，单位质量其它矿物所贡献孔隙体积的单位为cm3/g，孔隙孔径和吸附态甲烷厚度的单位均为nm；
步骤6：在步骤4和步骤5的基础上，根据吸附态甲烷密度小于固态甲烷密度，吸附态甲烷密度大于游离态甲烷密度，有机质所贡献孔隙内吸附态甲烷密度随孔隙孔径的增大而减小，吸附态甲烷密度随温度的升高而减小，吸附态甲烷密度随压力增大而增大的规律，建立如下方程组和目标函数，当目标函数f(ρ吸有机质jxy，h吸有机质jxy)的值最小时，即可获得温度为Tx、压力为Py条件下由有机质贡献孔径范围编号为j的孔隙内吸附态甲烷的密度ρ吸有机质jxy和厚度h吸有机质jxy，
式中， V吸有机质jxy为温度为Tx、压力为Py条件下孔径为有机质贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，ρ吸有机质jxy为温度为Tx、压力为Py条件下有机质贡献孔径范围编号为j的孔隙内吸附态甲烷密度，ρ游xy为温度为Tx、压力为Py条件下游离态甲烷密度，Q有机质xy为单位质量有机质赋存吸附态甲烷含量，ρ固为固体甲烷密度，h吸有机质jxy为温度为Tx、压力为Py条件下有机质所贡献孔径范围编号为j的孔隙内吸附态甲烷厚度，M为甲烷的摩尔质量，其值为16.0425 g/mol，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位是cm3/g，吸附态甲烷密度、游离态甲烷密度和固体甲烷密度的单位均为kg/m3，单位质量有机质赋存吸附态甲烷含量的单位是m3/t，吸附态甲烷厚度的单位是nm;
在步骤4和步骤5的基础上，根据吸附态甲烷密度小于固态甲烷密度，吸附态甲烷密度大于游离态甲烷密度，粘土矿物所贡献孔隙内吸附态甲烷密度随孔径的增大而减小，吸附态甲烷密度随温度的升高而减小，吸附态甲烷密度随压力增大而增大的规律，建立如下方程组和目标函数，当目标函数f(ρ吸粘土jxy，h吸粘土jxy)的值最小时，即可获得温度为Tx、压力为Py条件下由粘土矿物贡献孔径范围编号为j的孔隙内吸附态甲烷的密度ρ吸粘土jxy和厚度h吸粘土jxy，
式中， V吸粘土jxy为温度为Tx、压力为Py条件下粘土矿物贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，ρ吸粘土jxy为温度为Tx、压力为Py条件下粘土矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷密度，ρ游xy为温度为Tx、压力为Py条件下游离态甲烷密度，Q粘土xy为单位质量粘土矿物赋存吸附态甲烷含量，ρ固为固体甲烷密度，h吸粘土jxy为温度为Tx、压力为Py条件下粘土矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，M为甲烷的摩尔质量，其值为16.0425 g/mol，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位是cm3/g，吸附态甲烷密度、游离态甲烷密度和固体甲烷密度的单位均为kg/m3，单位质量粘土矿物赋存吸附态甲烷含量的单位是m3/t，吸附态甲烷厚度的单位是nm;
在步骤4和步骤5的基础上，根据吸附态甲烷密度小于固态甲烷密度，吸附态甲烷密度大于游离态甲烷密度，其它矿物所贡献孔隙内吸附态甲烷密度随孔径增大而减小，吸附态甲烷密度随温度升高而减小，吸附态甲烷密度随压力增大而增大的规律，建立如下方程组和目标函数，当目标函数f(ρ吸其它jxy，h吸其它jxy)的值最小时，即可获得温度为Tx、压力为Py条件下由其它矿物贡献孔径范围编号为j的孔隙内吸附态甲烷的密度ρ吸其它jxy和厚度h吸其它jxy，
式中， V吸其它jxy为温度为Tx、压力为Py条件下其它矿物贡献的孔径范围编号为j的孔隙内吸附态甲烷所占孔体积，ρ吸其它jxy为温度为Tx、压力为Py条件下其它矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷密度，ρ游xy为温度为Tx、压力为Py条件下游离态甲烷密度，Q其它xy为单位质量其它矿物赋存吸附态甲烷含量，ρ固为固体甲烷密度，h吸其它jxy为温度为Tx、压力为Py条件下其它矿物所贡献的孔径范围编号为j的孔隙内吸附态甲烷厚度，M为甲烷的摩尔质量，其值为16.0425 g/mol，j=1、2、…、7，是孔径范围由小到大的编号，x=1、2、…、m，是温度由低至高的编号，y=1、2、…、z，是压力由低至高的编号，吸附态甲烷所占孔体积的单位是cm3/g，吸附态甲烷密度、游离态甲烷密度和固体甲烷密度的单位均为kg/m3，单位质量其它矿物赋存吸附态甲烷含量的单位是m3/t，吸附态甲烷厚度的单位是nm。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811521652.4A CN109540764B (zh) | 2018-12-13 | 2018-12-13 | 评价泥页岩储层贡献孔隙内吸附态甲烷厚度和密度的方法 |
| PCT/CN2019/087062 WO2020119021A1 (zh) | 2018-12-13 | 2019-05-15 | 一种评价泥页岩储层有机质、粘土和其它矿物所贡献孔隙赋存吸附态甲烷厚度和密度的方法 |
| US16/686,129 US20200191697A1 (en) | 2018-12-13 | 2019-11-16 | Method for evaluating thickness and density of adsorbed methane in pores contributed by organic matter, clay and other minerals in mud shale reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811521652.4A CN109540764B (zh) | 2018-12-13 | 2018-12-13 | 评价泥页岩储层贡献孔隙内吸附态甲烷厚度和密度的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109540764A CN109540764A (zh) | 2019-03-29 |
| CN109540764B true CN109540764B (zh) | 2020-02-07 |
Family
ID=65854885
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811521652.4A Active CN109540764B (zh) | 2018-12-13 | 2018-12-13 | 评价泥页岩储层贡献孔隙内吸附态甲烷厚度和密度的方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN109540764B (zh) |
| WO (1) | WO2020119021A1 (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109540764B (zh) * | 2018-12-13 | 2020-02-07 | 中国石油大学(华东) | 评价泥页岩储层贡献孔隙内吸附态甲烷厚度和密度的方法 |
| CN110018293B (zh) * | 2019-05-05 | 2021-09-21 | 西安石油大学 | 一种考虑水敏性的页岩黏土多相耦合甲烷含气量计算方法 |
| CN110849844B (zh) * | 2019-11-21 | 2022-03-11 | 中国石油大学(华东) | 一种测量纯矿物纳米级圆柱管内吸附态甲烷厚度的方法 |
| CN112304799B (zh) * | 2020-06-24 | 2021-06-22 | 成都理工大学 | 一种页岩油储层不同赋存状态有机质定量分析的方法 |
| CN112304837B (zh) * | 2020-06-24 | 2021-06-25 | 成都理工大学 | 一种判断页岩油储层有机质微观赋存结构的方法 |
| CN114441403A (zh) * | 2020-11-03 | 2022-05-06 | 中国石油化工股份有限公司 | 一种基于页岩孔隙分布的游离油含量测定方法 |
| CN113049785A (zh) * | 2021-04-08 | 2021-06-29 | 中国石油大学(华东) | 一种评价泥页岩储层各组分赋存吸附态甲烷量方法 |
| CN115615366B (zh) * | 2022-11-21 | 2023-03-10 | 武汉普锐赛斯科技有限公司 | 一种页岩孔隙吸附层厚度检测装置及方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106568922A (zh) * | 2016-10-19 | 2017-04-19 | 中国石油天然气股份有限公司 | 一种地层温度和压力条件下页岩吸附气含量的计算方法 |
| CN107422100A (zh) * | 2017-09-05 | 2017-12-01 | 中国石油大学(华东) | 一种计算页岩气储层赋存吸附气含量的方法 |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4686364A (en) * | 1985-07-19 | 1987-08-11 | Schlumberger Technology Corporation | In situ determination of total carbon and evaluation of source rock therefrom |
| US7663011B2 (en) * | 1999-09-07 | 2010-02-16 | Lummus Technology Inc. | Mesoporous material with active metals |
| CN105842416B (zh) * | 2015-01-16 | 2018-06-29 | 中石化西南石油工程有限公司 | 一种基于兰格缪尔模型的页岩吸附气含量动态计算方法 |
| CN105043957B (zh) * | 2015-07-06 | 2018-05-11 | 成都理工大学 | 通过泥页岩等温吸附曲线分类判断孔隙结构的方法 |
| CN105158816B (zh) * | 2015-10-19 | 2017-11-28 | 中国石油大学(华东) | 预测页岩不同类型吸附气非均质性分布的方法 |
| CN105445159B (zh) * | 2015-11-11 | 2021-08-13 | 成都理工大学 | 一种得到孔径分布曲线及样品比表面积的方法 |
| CN106940279B (zh) * | 2017-05-18 | 2019-05-21 | 中国石油大学(华东) | 一种评价泥页岩储层赋存吸附气量的方法 |
| CN107271322B (zh) * | 2017-07-27 | 2019-09-06 | 中国石油大学(华东) | 一种定量评价泥页岩页岩气资源量及表征吸附气和游离气转化规律的方法 |
| CN107560994B (zh) * | 2017-09-05 | 2019-08-02 | 中国石油大学(华东) | 一种评价泥页岩有机质、粘土和其它矿物中孔径分布方法 |
| CN108458960B (zh) * | 2018-03-27 | 2019-10-29 | 中国石油大学(华东) | 富有机质泥页岩含氢组分、孔隙度及孔径的评价方法 |
| CN109540764B (zh) * | 2018-12-13 | 2020-02-07 | 中国石油大学(华东) | 评价泥页岩储层贡献孔隙内吸附态甲烷厚度和密度的方法 |
-
2018
- 2018-12-13 CN CN201811521652.4A patent/CN109540764B/zh active Active
-
2019
- 2019-05-15 WO PCT/CN2019/087062 patent/WO2020119021A1/zh active Application Filing
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106568922A (zh) * | 2016-10-19 | 2017-04-19 | 中国石油天然气股份有限公司 | 一种地层温度和压力条件下页岩吸附气含量的计算方法 |
| CN107422100A (zh) * | 2017-09-05 | 2017-12-01 | 中国石油大学(华东) | 一种计算页岩气储层赋存吸附气含量的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2020119021A1 (zh) | 2020-06-18 |
| CN109540764A (zh) | 2019-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109540764B (zh) | 评价泥页岩储层贡献孔隙内吸附态甲烷厚度和密度的方法 | |
| Klewiah et al. | Review of experimental sorption studies of CO2 and CH4 in shales | |
| Cai et al. | Shale gas transport model in 3D fractal porous media with variable pore sizes | |
| Fu et al. | Characteristics of pore structure and fractal dimension of low-rank coal: A case study of Lower Jurassic Xishanyao coal in the southern Junggar Basin, NW China | |
| Liu et al. | Determination of the absolute adsorption/desorption isotherms of CH4 and n-C4H10 on shale from a nano-scale perspective | |
| Sun et al. | Transport capacity of gas confined in nanoporous ultra-tight gas reservoirs with real gas effect and water storage mechanisms coupling | |
| Yu et al. | Multiscale gas transport behavior in heterogeneous shale matrix consisting of organic and inorganic nanopores | |
| Wang et al. | Measurement of dynamic adsorption–diffusion process of methane in shale | |
| Aljamaan et al. | Experimental investigation and Grand Canonical Monte Carlo simulation of gas shale adsorption from the macro to the nano scale | |
| Song et al. | Methane surface diffusion capacity in carbon-based capillary with application to organic-rich shale gas reservoir | |
| Wang et al. | Organic shale micropore and mesopore structure characterization by ultra-low pressure N2 physisorption: Experimental procedure and interpretation model | |
| Sun et al. | Mechanistic insight into the displacement of CH 4 by CO 2 in calcite slit nanopores: the effect of competitive adsorption | |
| Pang et al. | Experimental measurement and analytical estimation of methane absorption in shale kerogen | |
| 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 | |
| Liu et al. | Selective adsorption of CO2/CH4 mixture on clay-rich shale using molecular simulations | |
| Liu et al. | Absolute adsorption of CH4 on shale with the simplified local-density theory | |
| Zhang et al. | Monte Carlo simulations of methane adsorption on kaolinite as a function of pore size | |
| CN107387068B (zh) | 一种用于确定页岩气储层游离气含量的方法及系统 | |
| CN105842416B (zh) | 一种基于兰格缪尔模型的页岩吸附气含量动态计算方法 | |
| CN106547966A (zh) | 一种页岩油吸附量与可动量评价模型及其建立、应用方法 | |
| Aljamaan et al. | In-depth experimental investigation of shale physical and transport properties | |
| Ren et al. | Modeling of mixed-gas adsorption on shale using hPC-SAFT-MPTA | |
| Zeng et al. | Characterization of gas transport in shale: A multi-mechanism permeability modeling approach | |
| CN106932323A (zh) | 一种页岩气藏气体有效孔隙度反演方法 | |
| Zhang et al. | Fractal pore structure model and multilayer fractal adsorption in shale |
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 |