CN110223736A - 一种用于计算页岩有机质不同孔径表面吸附量的方法 - Google Patents
一种用于计算页岩有机质不同孔径表面吸附量的方法 Download PDFInfo
- Publication number
- CN110223736A CN110223736A CN201910502651.3A CN201910502651A CN110223736A CN 110223736 A CN110223736 A CN 110223736A CN 201910502651 A CN201910502651 A CN 201910502651A CN 110223736 A CN110223736 A CN 110223736A
- Authority
- CN
- China
- Prior art keywords
- atom
- organic matter
- methane
- energy
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Chemical & Material Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种用于计算页岩有机质不同孔径表面吸附量的方法,该方法包含:(1),利用石墨构建石墨烯的层间结构建立有机质模型,将甲烷构型导入建好的有机质模型中;(2)进行吸附模拟,对分子体系进行能量最小化,通过调整原子位置变化获得初始的稳定构型;(3)计算有机质孔径中充填的甲烷分子数,优先吸附位的判断,计算全过程采用COMPASS力场,计算分子间范德华作,及甲烷与骨架的静电势能;(4)在与孔壁面平行方向上将纳米孔划分为K个单元,统计每个盒子内的原子种类和数目,计算其流体密度平均值;(5)计算Gibbs吸附量。本发明的方法对吸附相密度的计算与实际情况吻合更好,精度更高。
Description
技术领域
本发明属于页岩气勘探领域,具体涉及一种用于计算页岩有机质不同孔径表面吸附量的方法。
背景技术
国内外学者在对页岩的吸附性能进行研究时发现,在利用容量法对页岩吸附甲烷进行测试时,吸附相体积的存在将减小样品中自由气体的体积,进而对页岩吸附能力的测试结果产生影响。然而,容量法忽略了吸附相体积存在,这必将导致对于页岩吸附能力的评价存在较大偏差。
针对上述问题,国内外学者试图从不同的角度校正吸附相体积对吸附气能力测试所产生的影响。然而,无论以何种方式对吸附相体积所产生的影响进行校正,终将遵循物质守恒原理,因此吸附相密度的确定是无法回避的问题,且吸附相密度的准确性将直接影响校正结果。
目前,国内外学者基于不同方法对甲烷吸附相密度进行多种求解,主要可以归纳为以下3种常见的方法:
方法1:理论上认为气体吸附类似于气体液化,所以可以直接选取常压沸点时液体甲烷密度作为吸附相密度,取值为0.423g/cm3;
方法2:通过实验得到的Gibbs吸附量下降段线性化拟合求取吸附相密度。此法认为,Gibbs吸附量达到极大值拐点后随压力进一步增加呈线性下降趋势,说明吸附达到饱和,吸附饱和后吸附相密度和吸附相体积不在进一步增加而保持恒定。利用数据拟合线性下降段直线,直线斜率的绝对值即为吸附相体积,直线与横坐标交点值即为吸附相密度;
方法3:以数学优化思想根据实测Gibbs吸附量通过最小二乘法拟合求得吸附相密度。
上述三种吸附相密度求取的方法被广泛应用于校正Gibbs吸附量曲线,从而得到页岩甲烷表面吸附量,但哪一种方法得到的校正结果更加合理,更接近实际情况仍然存在争议。
为了校正页岩超临界吸附曲线,使校正获得的表面吸附量更加准确合理,对于以上3种方法进行了实际应用分析:
对于第1种方法,因为吸附相密度在不同样品、不同实验条件下会有较大变化,不能简单地认为吸附相的密度为一定值。并且,研究发现临界状态下甲烷吸附相的密度应处于临界密度(0.162g/cm3)与常压沸点液体甲烷密度(0.423g/cm3)之间差距较大,因此该方法的适用性较差;
对于第2种方法,部分样品拟合得到的吸附相密度大于常压沸点液体甲烷密度,所以该方法并不对所有样品都适用;
对于第3种方法,依然存在部分样品,拟合得到的甲烷吸附相密度大于常压沸点液体甲烷密度。这表明,该方法并不对所有样品均可靠。
综上所述,以上3种常用的计算吸附相密度的方法对于页岩中甲烷吸附相密度的计算并不适用,不具有普遍性。而且,以上三种计算方法都没有考虑孔径大小对于吸附相密度的影响,不能很好的反应真实的吸附情况,这必将对客观评价页岩吸附能力产生较大影响。
发明内容
本发明的目的是提供一种用于计算页岩有机质不同孔径表面吸附量的方法,该方法解决了现有方法不能很好的反应甲烷的真实吸附情况的问题,对吸附相密度的计算与实际情况吻合更好,精度更高。
为了达到上述目的,本发明提供了一种用于计算页岩有机质不同孔径表面吸附量的方法,其特征在于,该方法包含:
(1)建立有机质模型,利用石墨构建石墨烯的层间结构,相邻层间距为0.335nm,将甲烷构型导入建好的有机质模型中;
(2)进行吸附模拟,利用Smart算法对分子体系进行能量最小化,通过调整原子位置变化获得初始的稳定构型;
(3)利用巨正则系统蒙特卡洛法计算有机质孔径中充填的甲烷分子数,优先吸附位的判断利用巨正则系综蒙特卡洛法,设定平衡步数和化产步数;计算全过程采用COMPASS力场,其中分子间范德华作用采用Atom Based求和方法,甲烷与骨架的静电势能采用Ewald求和方法;
(4)在与孔壁面平行方向上将纳米孔划分为K个单元,导出充填流体分子后的骨架模型,获得吸附体系中各原子的坐标,统计每个盒子内的原子种类和数目,并计算甲烷吸附相密度ρa;
(5)以Gibbs吸附量作为表面吸附量,为:
式(2)中,Vex为平衡压力下甲烷的Gibbs吸附量;Vabs为平衡压力下甲烷表面吸附量;ρa为甲烷吸附相密度;ρg为平衡压力下气体相密度。
其中,利用Smart算法对分子体系进行能量最小化,依靠解牛顿运动方程的模拟计算,对体系中不同状态分子进行样本抽取,计算分子体系的积分,并以体系积分结果为初始状态,计算模拟体系的热力学性质、以及系统的势能、各原子所受的力、加速度、经过t时间后各分子的位置及速度,通过对这些数据运用统计平均的方法计算宏观性质以及空间分布,获得初始的稳定构型。
优选地,所述流体密度平均值ρa:
式(1)中,i为原子种类;Mi为i原子的相对分子质量;ki为小盒子内含有的i原子数量;NA为阿伏伽德罗常数;Vj为第j个小盒子的体积。
优选地,在步骤(2)中,体系的总能量由体系中分子的总动能和总势能构成,总势能通过体系中分子中各个原子的函数来表达,为:
U=UVDW+Uin (2);
式(2)中,U为总势能,UVDW为范德华作用,Uin为内部势能;
式(3)中,uij表示各i原子和j原子之间的范德华作用;rij表示各i原子和j原子之间的距离;
系统中任意一个原子i所受的力Fi为:
式(4)中,为在三维空间直角坐标系内与x轴、y轴、z轴方向相同的单位向量,表示空间向量的一组基底;xi、yi、zi分别表示原子i在三维空间直角坐标系中x轴、y轴、z轴方向的坐标值;
由牛顿运动定律,可得原子i所受的加速度ai为:
将牛顿运动定律方程对时间积分,可得原子i经过时间t后的速度vi与位置ri为:
式(8)中,为原子i的初始速度,为原子i的初始位置。
优选地,在步骤(3)中,设定平衡步数和化产步数均为1×106。
本发明的用于计算页岩有机质不同孔径表面吸附量的方法,解决了现有方法不能很好的反应甲烷的真实吸附情况的问题,具有以下优点:
(1)本发明的方法,利用分子动力学原理,模拟了甲烷气体在页岩有机质表面的吸附情况,本发明方法在吸附相密度的模拟中考虑到孔径,更加合理,明显优于传统方法;
(2)本发明的方法,在计算结果上,对吸附相密度的计算与实际情况吻合更好,精度更高,明显优于传统方法;
(3)本发明的方法,利用分子动力学模拟从微观的角度考虑了宏观实验没有办法进行的不同孔径大小对于吸附相密度的影响;
(4)本发明方法简单易行,准确度高,实用性较强,且有较好的推广应用价值。
附图说明
图1为分子动力学模拟方法流程图。
图2为本发明方法的流程图。
图3为本发明在某油田某井吸附气量测井评价成果图。
图4为未修正模型计算结果与实验测试结果对比图。
图5为本发明修正模型计算结果与实验测试结果对比图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
一种用于计算页岩有机质不同孔径表面吸附量的方法,如图2所示,为本发明方法的流程图,该方法在某油田进行测试时,具体步骤如下:
(1)建立模型
利用石墨构建石墨烯的层间结构来代替有机质的层间结构,采用3层石墨稀片层作为有机质纳米孔的固体壁,且保证相邻两层之间平行距离为0.335nm。然后,将MaterialsStudio软件包中已经构建完成的甲烷构型导入建好的有机质模型中,进行后续的吸附模拟。
(2)模拟处理
如图1所示,为分子动力学模拟方法流程图,利用Smart算法对分子体系进行能量最小化,通过调整原子位置变化获得初始的稳定构型。然后,利用巨正则系统蒙特卡洛法(GCMC,Grand Canonical Monte Carlo)计算有机质孔径中充填的甲烷分子数,优先吸附位的判断利用巨正则系综蒙特卡洛法,平衡步数和化产步数均为1×106。
计算全过程采用COMPASS力场,其中分子间范德华作用采用Atom Based求和方法,甲烷与骨架的静电势能采用Ewald求和方法。
其中,体系的总能量由体系中分子的总动能和总势能构成,总势能通过体系中分子中各个原子的函数来表达,为:
U=UVDW+Uin (2);
式(2)中,U为总势能,UVDW为范德华作用,Uin为内部势能;
式(3)中,uij表示各i原子和j原子之间的范德华作用;;rij表示各i原子和j原子之间的距离;
系统中任意一个原子i所受的力Fi为:
式(4)中,为在三维空间直角坐标系内与x轴、y轴、z轴方向相同的单位向量,表示空间向量的一组基底;xi、yi、zi分别表示原子i在三维空间直角坐标系中x轴、y轴、z轴方向的坐标值;
由牛顿运动定律,可得原子i所受的加速度ai为:
将牛顿运动定律方程对时间积分,可得原子i经过时间t后的速度vi与位置ri为:
式(8)中,为原子i的初始速度,为原子i的初始位置。
(3)吸附相密度求取
为直接模拟流体在狭缝纳米孔的密度分布情况,首先需在与孔壁面平行方向上将纳米孔划分为K个单元(即K个小盒子,一般厚度取0.04~0.1nm)。然后,通过MaterialStudios软件导出充填流体分子后的骨架模型,获得吸附体系中各原子的坐标,统计每个盒子内的原子种类和数目,并计算其流体密度平均值,其函数表达式见式(1):
式(1)中,i为原子种类;Mi为i原子的相对分子质量;ki为小盒子内含有的i原子数量,单位为个;NA为阿伏伽德罗常数,单位为mol-1;Vj为第j个小盒子的体积,单位为m3。
(4)计算表面吸附量
页岩等温吸附应用的基本理论为单分子层吸附理论,Langmiur方程模型就基于该理论建立。等温吸附实验测得的气体吸附量为Gibbs吸附量,而理论吸附模型中的气体吸附量为表面吸附量。Gibbs对吸附量定义为吸附相中超过体相密度的过剩量,其关系为公式(2)所示。
式(9)中,Vex为平衡压力下甲烷的Gibbs吸附量;Vabs为平衡压力下甲烷表面吸附量(m3/t);ρa为甲烷吸附相密度(g/cm3);ρg为平衡压力下气体相密度(g/cm3)。
由于页岩等温吸附需要区分Gibbs吸附量和表面吸附量,这就需要知道吸附相密度ρa和气体相密度ρg,其中ρg可以通过查找表1得到。
表1甲烷气体密度与压力关系
如图3所示,为本发明在某油田某井吸附气量测井评价成果图(1315~1440m),修正模型的计算结果和实验测试结果吻合效果较好,并且与未修正模型计算结果差值达到了100%以上。因此为了正确的估算页岩储层的吸附能力必须使用修正后的模型进行精确计算。
本发明的方法在某油田某井进行使用,如图4结果显示,未修正模型的定量计算结果与实验室测试吸附气量结果相差较大,误差为41.7%,说明未修改模型已经不能够用于精确计算页岩吸附气量。如图5结果显示,修正模型的定量计算结果在整体趋势上与实验室测试吸附气量结果基本一致且误差为7.5%,说明修正模型的定量计算比较精确,能够用于精确计算页岩吸附气量。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (4)
1.一种用于计算页岩有机质不同孔径表面吸附量的方法,其特征在于,该方法包含:
(1)建立有机质模型,利用石墨构建石墨烯的层间结构,相邻层间距为0.335nm,将甲烷构型导入建好的有机质模型中;
(2)进行吸附模拟,利用Smart算法对分子体系进行能量最小化,通过调整原子位置变化获得初始的稳定构型;
(3)利用巨正则系统蒙特卡洛法计算有机质孔径中充填的甲烷分子数,优先吸附位的判断利用巨正则系综蒙特卡洛法,设定平衡步数和化产步数;计算全过程采用COMPASS力场,其中分子间范德华作用采用AtomBased求和方法,甲烷与骨架的静电势能采用Ewald求和方法;
(4)在与孔壁面平行方向上将纳米孔划分为K个单元,导出充填流体分子后的骨架模型,获得吸附体系中各原子的坐标,统计每个盒子内的原子种类和数目,并计算甲烷吸附相密度ρa;
(5)以Gibbs吸附量作为表面吸附量,为:
式中,Vex为平衡压力下甲烷的Gibbs吸附量;Vabs为平衡压力下甲烷表面吸附量;ρa为甲烷吸附相密度;ρg为平衡压力下气体相密度;
其中,利用Smart算法对分子体系进行能量最小化,依靠解牛顿运动方程的模拟计算,对体系中不同状态分子进行样本抽取,计算分子体系的积分,并以体系积分结果为初始状态,计算模拟体系的热力学性质、以及系统的势能、各原子所受的力、加速度、经过t时间后各分子的位置及速度,通过对这些数据运用统计平均的方法计算宏观性质以及空间分布,获得初始的稳定构型。
2.根据权利要求1所述的用于计算页岩有机质不同孔径表面吸附量的方法,其特征在于,所述流体密度平均值ρa:
式中,i为原子种类;Mi为i原子的相对分子质量;ki为小盒子内含有的i原子数量;NA为阿伏伽德罗常数;Vj为第j个小盒子的体积。
3.根据权利要求1所述的用于计算页岩有机质不同孔径表面吸附量的方法,其特征在于,在步骤(2)中,体系的总能量由体系中分子的总动能和总势能构成,总势能通过体系中分子中各个原子的函数来表达,为:
U=UVDW+Uin (2);
式(2)中,U为总势能,UVDW为范德华作用能,Uin为内部势能;
范德华作用能UVDW通过体系中各原子之间的相互作用近似加成得到,为:
式(3)中,uij表示各i原子和j原子之间的范德华作用能;rij表示各i原子和j原子之间的距离;
系统中任意一个原子i所受的力Fi为:
式(4)中,为在三维空间直角坐标系内与x轴、y轴、z轴方向相同的单位向量,表示空间向量的一组基底;xi、yi、zi分别表示原子i在三维空间直角坐标系中x轴、y轴、z轴方向的坐标值;
由牛顿运动定律,可得原子i所受的加速度ai为:
将牛顿运动定律方程对时间积分,可得原子i经过时间t后的速度vi与位置ri为:
式(8)中,为原子i的初始速度,为原子i的初始位置。
4.根据权利要求1所述的用于计算页岩有机质不同孔径表面吸附量的方法,其特征在于,在步骤(3)中,设定平衡步数和化产步数均为1×106。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910502651.3A CN110223736A (zh) | 2019-06-11 | 2019-06-11 | 一种用于计算页岩有机质不同孔径表面吸附量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910502651.3A CN110223736A (zh) | 2019-06-11 | 2019-06-11 | 一种用于计算页岩有机质不同孔径表面吸附量的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110223736A true CN110223736A (zh) | 2019-09-10 |
Family
ID=67816339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910502651.3A Pending CN110223736A (zh) | 2019-06-11 | 2019-06-11 | 一种用于计算页岩有机质不同孔径表面吸附量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110223736A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111104735A (zh) * | 2019-12-12 | 2020-05-05 | 中国人民解放军96901部队23分队 | 一种用于计算氡在多孔材料中吸附量的方法 |
CN111664825A (zh) * | 2020-06-10 | 2020-09-15 | 西安石油大学 | 一种提高孔结构参数测定准确性的方法 |
CN112229754A (zh) * | 2020-09-29 | 2021-01-15 | 西南石油大学 | 一种迭代平均法计算甲烷在页岩中吸附量的方法 |
CN112798462A (zh) * | 2021-01-04 | 2021-05-14 | 西南石油大学 | 一种考虑页岩孔径的气体吸附量计算模型构建方法 |
CN113591408A (zh) * | 2021-08-02 | 2021-11-02 | 西南石油大学 | 液态硫在碳酸盐岩纳米孔隙中微观吸附特征参数的确定方法 |
CN113779830A (zh) * | 2021-08-31 | 2021-12-10 | 西安交通大学 | 一种纳米材料微孔结构与气体吸附扩散性质间构效关系的获取方法 |
CN114166723A (zh) * | 2021-12-21 | 2022-03-11 | 中国石油大学(华东) | 纳米多孔介质中气体的量子物理吸附行为预测方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108446526A (zh) * | 2018-05-22 | 2018-08-24 | 中国地质大学(武汉) | 基于吸附特征曲线与特征方程的页岩吸附量的预测方法 |
CN108595778A (zh) * | 2018-04-02 | 2018-09-28 | 中国石油大学(北京) | 一种页岩复合有机质分子模型的构建方法 |
-
2019
- 2019-06-11 CN CN201910502651.3A patent/CN110223736A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108595778A (zh) * | 2018-04-02 | 2018-09-28 | 中国石油大学(北京) | 一种页岩复合有机质分子模型的构建方法 |
CN108446526A (zh) * | 2018-05-22 | 2018-08-24 | 中国地质大学(武汉) | 基于吸附特征曲线与特征方程的页岩吸附量的预测方法 |
Non-Patent Citations (3)
Title |
---|
代建伟: ""页岩纳米孔隙中小分子气体超临界吸附微观机理研究"", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
张亮仁: "《常用计算机辅助设计软件教程》", 30 April 2017, 中国医药科技出版社 * |
熊健 等: ""甲烷在官能团化石墨中吸附行为的影响因素研究"", 《中国矿业大学学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111104735B (zh) * | 2019-12-12 | 2023-05-09 | 中国人民解放军96901部队23分队 | 一种用于计算氡在多孔材料中吸附量的方法 |
CN111104735A (zh) * | 2019-12-12 | 2020-05-05 | 中国人民解放军96901部队23分队 | 一种用于计算氡在多孔材料中吸附量的方法 |
CN111664825A (zh) * | 2020-06-10 | 2020-09-15 | 西安石油大学 | 一种提高孔结构参数测定准确性的方法 |
CN111664825B (zh) * | 2020-06-10 | 2021-07-09 | 西安石油大学 | 一种提高孔结构参数测定准确性的方法 |
CN112229754B (zh) * | 2020-09-29 | 2022-01-28 | 西南石油大学 | 一种迭代平均法计算甲烷在页岩中吸附量的方法 |
CN112229754A (zh) * | 2020-09-29 | 2021-01-15 | 西南石油大学 | 一种迭代平均法计算甲烷在页岩中吸附量的方法 |
CN112798462B (zh) * | 2021-01-04 | 2022-06-24 | 西南石油大学 | 一种考虑页岩孔径的气体吸附量计算模型构建方法 |
CN112798462A (zh) * | 2021-01-04 | 2021-05-14 | 西南石油大学 | 一种考虑页岩孔径的气体吸附量计算模型构建方法 |
CN113591408A (zh) * | 2021-08-02 | 2021-11-02 | 西南石油大学 | 液态硫在碳酸盐岩纳米孔隙中微观吸附特征参数的确定方法 |
CN113779830A (zh) * | 2021-08-31 | 2021-12-10 | 西安交通大学 | 一种纳米材料微孔结构与气体吸附扩散性质间构效关系的获取方法 |
CN113779830B (zh) * | 2021-08-31 | 2022-12-09 | 西安交通大学 | 一种纳米材料微孔结构与气体吸附扩散性质间构效关系的获取方法 |
CN114166723A (zh) * | 2021-12-21 | 2022-03-11 | 中国石油大学(华东) | 纳米多孔介质中气体的量子物理吸附行为预测方法及系统 |
CN114166723B (zh) * | 2021-12-21 | 2023-06-16 | 中国石油大学(华东) | 纳米多孔介质中气体的量子物理吸附行为预测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110223736A (zh) | 一种用于计算页岩有机质不同孔径表面吸附量的方法 | |
Vasileiadis et al. | Transport properties of shale gas in relation to kerogen porosity | |
Zeng et al. | Molecular simulation of carbon dioxide and methane adsorption in shale organic nanopores | |
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 | |
Wolfsberg | Isotope effects | |
CN106198297A (zh) | 一种精确计算甲烷在页岩上真实吸附量的方法 | |
Imberti et al. | Formic and acetic acid aggregation in the liquid state | |
CN110672813A (zh) | 一种页岩含气量计算方法 | |
Brugman et al. | Calcite (104) surface–electrolyte structure: A 3D comparison of surface X-ray diffraction and simulations | |
CN106442268A (zh) | 一种页岩介孔孔径分布的检测方法 | |
CN106767671B (zh) | 基于三维电子罗盘的地质结构面产状计算方法 | |
Pethes | A comparison of classical interatomic potentials applied to highly concentrated aqueous lithium chloride solutions | |
CN107818379A (zh) | 吸附气量确定方法和装置 | |
Zhao et al. | Microscopic simulation of methane adsorption in organic matter | |
Thakre et al. | Computing anisotropic cavity potential for clathrate hydrates | |
Dürholt et al. | The impact of mesopores on the mechanical stability of HKUST‐1: A multiscale investigation | |
Fortini et al. | Phase behavior and structure of model colloid-polymer mixtures confined between two parallel planar walls | |
Zhan et al. | Fractal Characteristics of Deep Shales in Southern China by Small-Angle Neutron Scattering and Low-Pressure Nitrogen Adsorption | |
Wang et al. | The characteristic development of micropores in deep coal and its relationship with adsorption capacity on the eastern margin of the Ordos Basin, China | |
Wang et al. | Determining the geometric surface area of mesoporous materials | |
CN104200025B (zh) | 一种定量分析沸石咪唑酯骨架材料路易斯酸性的方法 | |
Soper et al. | Adsorption of simple gases into the porous glass MCM-41 | |
Dimitrova et al. | A spherical cavity model for quadrupolar dielectrics | |
You et al. | Impacts of pore structure on the occurrence of free oil in lacustrine shale pore networks | |
Gatta et al. | Single-crystal neutron diffraction and Mössbauer spectroscopic study of hureaulite,(Mn, Fe) 5 (PO4) 2 (HPO4) 2 (H2O) 4 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190910 |
|
RJ01 | Rejection of invention patent application after publication |