CN109655940A - 页岩各向异性岩石物理模型建模方法 - Google Patents
页岩各向异性岩石物理模型建模方法 Download PDFInfo
- Publication number
- CN109655940A CN109655940A CN201710949034.9A CN201710949034A CN109655940A CN 109655940 A CN109655940 A CN 109655940A CN 201710949034 A CN201710949034 A CN 201710949034A CN 109655940 A CN109655940 A CN 109655940A
- Authority
- CN
- China
- Prior art keywords
- shale
- clay
- anisotropic
- model
- elastic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 50
- 239000004927 clay Substances 0.000 claims abstract description 111
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 108
- 239000011707 mineral Substances 0.000 claims abstract description 108
- 239000000203 mixture Substances 0.000 claims abstract description 83
- 239000011148 porous material Substances 0.000 claims abstract description 82
- 239000005416 organic matter Substances 0.000 claims abstract description 73
- 239000011159 matrix material Substances 0.000 claims abstract description 49
- 239000002245 particle Substances 0.000 claims abstract description 49
- 239000012530 fluid Substances 0.000 claims abstract description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 28
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 20
- 239000002734 clay mineral Substances 0.000 claims description 17
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 15
- 229910052799 carbon Inorganic materials 0.000 claims description 15
- 238000002637 fluid replacement therapy Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 239000010433 feldspar Substances 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims 2
- 229910052742 iron Inorganic materials 0.000 claims 1
- 238000012512 characterization method Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 10
- 230000000704 physical effect Effects 0.000 description 10
- 238000012360 testing method Methods 0.000 description 10
- 238000011161 development Methods 0.000 description 9
- 238000004088 simulation Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 229910021532 Calcite Inorganic materials 0.000 description 4
- 229910000514 dolomite Inorganic materials 0.000 description 4
- 239000010459 dolomite Substances 0.000 description 4
- 239000011028 pyrite Substances 0.000 description 4
- 229910052683 pyrite Inorganic materials 0.000 description 4
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 4
- 239000010453 quartz Substances 0.000 description 4
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 4
- 235000021185 dessert Nutrition 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 239000002023 wood Substances 0.000 description 3
- CWYNVVGOOAEACU-UHFFFAOYSA-N Fe2+ Chemical compound [Fe+2] CWYNVVGOOAEACU-UHFFFAOYSA-N 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 229910001608 iron mineral Inorganic materials 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000008093 supporting effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000001564 chemical vapour infiltration Methods 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000011019 hematite Substances 0.000 description 1
- 229910052595 hematite Inorganic materials 0.000 description 1
- LIKBJVNGSGBSGK-UHFFFAOYSA-N iron(3+);oxygen(2-) Chemical compound [O-2].[O-2].[O-2].[Fe+3].[Fe+3] LIKBJVNGSGBSGK-UHFFFAOYSA-N 0.000 description 1
- 238000011545 laboratory measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000007779 soft material Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Toys (AREA)
Abstract
本发明公开了一种页岩各向异性岩石物理模型建模方法,包括将页岩基质视为脆性矿物、有机质和粘土组成的混合物;将粘土颗粒视为具有固定不变弹性刚度矩阵的各向异性元,引入粘土颗粒定向指数表征粘土颗粒定向排列程度;将页岩总孔隙划为脆性孔、粘土孔和有机质孔三种孔隙;获取含孔脆性矿物和含孔有机质的混合物的弹性参数;求取脆性矿物和有机质的混合物的弹性参数;求取含束缚水的孔隙粘土介质弹性参数;求取页岩基质和页岩干岩石骨架等效弹性参数;基于页岩基质和页岩干岩石骨架等效弹性参数得到流体饱和页岩等效弹性参数,完成流体饱和状态下页岩各向异性岩石物理模型的构建。本发明能较准确地预测页岩的纵横波速度和各向异性参数。
Description
技术领域
本发明属于岩石物理模型建模领域,更具体地,涉及一种页岩各向异性岩石物理模型建模方法。
背景技术
页岩的微观结构比较复杂,影响页岩弹性性质的因素多种多样,复杂的矿物组分、有机质、粘土颗粒的定向排列以及微裂缝的发育等对页岩弹性性质都会产生影响,再加上页岩具有典型的强各向异性特征,这都是常规储层岩石物理建模所没有遇到的难题。目前,还没有一种比较有效的描述泥页岩弹性性质的岩石物理模型建模方法。
基于对页岩微观结构和各向异性特征的研究表明,粘土矿物含量和其定向排列是其各向异性产生的主要因素,例如国内龙马溪组页岩。而国内页岩岩石物理模型的构建多没有考虑粘土矿物的定向性这一特征,把粘土矿物看作各向同性介质添加到模型中,也没有考虑页岩的微观结构特征。
因此,有必要提供一种页岩的各向异性岩石物理模型的建模方法,能够考虑页岩各向异性产生的主要因素,有效的描述页岩的弹性性质。
发明内容
本发明的目的是解决页岩岩石物理模型研究中存在的难题,旨在提供一种适用各向异性岩石物理模型建模方法,以期能够更加有效的表征页岩宏观岩石物理性质与微观结构之间的依赖关系。
为了实现上述目的,本发明提供了一种页岩各向异性岩石物理模型建模方法,该建模方法包括:
将页岩基质视为脆性矿物、有机质和粘土组成的混合物;
将粘土颗粒视为具有固定不变弹性刚度矩阵的各向异性元,引入粘土颗粒定向指数表征粘土颗粒定向排列程度;
将页岩总孔隙划分为脆性孔、粘土孔和有机质孔三种孔隙类型;
求取含孔隙脆性矿物混合物的弹性参数和含有有机质孔隙的有机质的弹性参数,进而建立含孔脆性矿物和含孔有机质的混合物的物理模型,获取含孔脆性矿物和含孔有机质的混合物的弹性参数;
求取脆性矿物和有机质的混合物的弹性参数;
求取含束缚水的孔隙粘土介质弹性参数;进而求取页岩基质等效弹性参数和求取页岩干岩石骨架等效弹性参数;
基于页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数得到流体饱和页岩等效弹性参数,从而完成流体饱和状态下页岩各向异性岩石物理模型的构建。
本发明基于页岩岩石物理实验测试和微观结构研究,在明确了页岩岩石物理性质、各向异性特征及其影响因素的基础上,研发了页岩各向异性岩石物理模型建模方法。通过本发明所建立的页岩各向异性岩石物理模型,能够更加有效的表征页岩宏观岩石物理性质与其微观结构之间的依赖关系,可较准确地预测页岩的纵波速度、横波速度和各向异性参数,从而为页岩气甜点预测、储层改造和监测等勘探开发技术研究和生产应用提供岩石物理基础支持。
本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。
附图说明
通过结合附图对本发明示例性实施方式进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了本发明的页岩各向异性岩石物理模型建模方法的流程图。
图2示出了根据本发明的一个实施例的页岩各向异性岩石物理模型建模方法的建模示意图。
图3示出了根据本发明的一个实施例的页岩各向异性岩石物理模型建模方法中求取页岩岩石基质参数的示意图。
图4示出了根据本发明的一个实施例的页岩各向异性岩石物理模型建模方法中求取页岩干岩石骨架参数的示意图。
图5示出了根据本发明的一个实施例的页岩各向异性岩石物理模型建模方法中求取页岩饱和岩石参数的示意图。
图6示出了根据本发明的一个实施例的页岩各向异性岩石物理模型模拟结果与页岩样品实验测试结果的对比图。
图7示出了根据本发明的一个实施例的页岩各向异性岩石物理模型模拟结果与页岩气井测井资料的对比图。
图8a-图8f分别示出了本发明的一个实施例的页岩各向异性岩石物理模型模拟有机质含量和有机质孔发育对页岩垂向纵波速度VPv、垂向SV横波速度VSV、Vp/Vs和各向异性参数ε、δ、γ影响的结果图。
图9a-图9f分别示出了本发明的一个实施例的页岩各向异性岩石物理模型模拟脆性矿物含量和粘土矿物定向性对页岩垂向纵波速度VPV、垂向SV横波速度VSV、Vp/Vs和各向异性参数ε、δ、γ影响的结果图。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了本发明的页岩各向异性岩石物理模型建模方法的流程图。图1中,步骤104、105和106中求取各混合物的弹性参数和含束缚水的孔隙粘土介质弹性参数的步骤不分先后,以下描述顺序仅为示例性的。
如图1所示,本发明提供了一种页岩各向异性岩石物理模型建模方法,该建模方法包括:步骤101,将页岩基质视为脆性矿物、有机质和粘土组成的混合物;步骤102,将粘土颗粒视为具有固定不变弹性刚度矩阵的各向异性元,引入粘土颗粒定向指数表征粘土颗粒定向排列程度;步骤103,将页岩总孔隙划分为脆性孔、粘土孔和有机质孔三种孔隙类型;步骤104,求取含孔隙脆性矿物混合物的弹性参数和含有有机质孔隙的有机质的弹性参数,进而建立含孔脆性矿物和含孔有机质的混合物的物理模型,获取含孔脆性矿物和含孔有机质的混合物的弹性参数;步骤105,求取脆性矿物和有机质的混合物的弹性参数;步骤106,求取含束缚水的孔隙粘土介质弹性参数;进而求取页岩基质等效弹性参数和求取页岩干岩石骨架等效弹性参数;步骤107,基于页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数得到流体饱和页岩等效弹性参数,从而完成流体饱和状态下页岩各向异性岩石物理模型的构建。
本发明基于页岩的岩性特征、微观结构和岩石物理特征,在保证理论模型物理机制和岩石结构相统一的基础上,将粘土颗粒视为各向异性元,引入粘土定向指数表征其定向排列程度,综合利用各向异性自洽模型和各向异性微分等效模型基础上,构建了页岩各向异性岩石物理模型,并利用岩石物理实验测试数据和测井数据验证了模型的准确性和适用性。
以下具体说明本发明的:
步骤101,将页岩基质视为脆性矿物、有机质和粘土组成的混合物。
在一个示例中,所述脆性矿物包括硅质、长石、钙质和铁质。
页岩脆性矿物包括硅质(石英)、长石、钙质(方解石和白云石)、铁质(黄铁矿和赤铁矿)等。
步骤102,将粘土颗粒视为具有固定不变弹性刚度矩阵的各向异性元,引入粘土颗粒定向指数表征粘土颗粒定向排列程度。
粘土含量和定向性是页岩各向异性特征产生的主要因素,因此将页岩中的粘土颗粒视为各向异性元,各向异性元具有固定不变的弹性刚度矩阵,通过引入各向异性元来实现对粘土矿物定向排列所表征的各向异性特征的模拟。各向异性元的弹性刚度矩阵为:
步骤103,将页岩总孔隙划分为脆性孔、粘土孔和有机质孔三种孔隙类型。建立物理模型时,将页岩总孔隙划分为脆性孔、粘土孔和有机质孔三种孔隙类型。
步骤104,求取含孔隙脆性矿物混合物的弹性参数和含有有机质孔隙的有机质的弹性参数,进而建立含孔脆性矿物和含孔有机质的混合物的物理模型,获取含孔脆性矿物和含孔有机质的混合物的弹性参数。
在一个示例中,所述求取含孔隙脆性矿物混合物的弹性参数包括:
采用Hashin-Shtrikman边界模型计算脆性矿物各向同性混合情况下的弹性参数上下边界,求上下边界的平均值作为脆性矿物混合物的弹性参数;
利用各向同性微分等效介质模型,将脆性孔加入到所述脆性矿物混合物中,得到含孔隙脆性矿物混合物的弹性参数。
Hashin-Shtrikman边界模型的计算公式由下式给出:
K1,K2为各组分的体积参数;μ1、μ2为各组分的剪切参数;f1和f2为各组分的体积含量。KHS+和KHS-为体积参数的上限和下限;μHS+和μHS-为剪切参数的上限和下限。上下限是通过交换哪个材料是1和2来求得的,当硬的材料定义为1时求得的是上限,当软的材料定义为1时求得的是下限。
在一个示例中,所述求取含有有机质孔隙的有机质弹性参数包括:
基于经验公式VKerogen(v%)=2.2466×TOC(wt%),将总有机碳质量百分含量转换为有机质的体积含量;
基于各向同性自洽模型和各向同性微分等效介模型的联合模型在有机质中添加有机质孔隙,得到含有有机质孔隙的有机质的弹性参数。
在一个示例中,所述建立含孔脆性矿物和含孔有机质的混合物的物理模型包括:
当总有机碳质量百分含量小于2%时,采用微分等效介质模型建立含孔脆性矿物和含孔有机质的混合物的物理模型;
当总有机碳质量百分含量大于2%时,采用自洽模型和微分等效介质模型的联合模型建立含孔脆性矿物和含孔有机质的混合物的物理模型。
脆性矿物和有机质的混合物,当总有机碳质量百分含量(TOC)含量小于2%时,有机质多分散分布在矿物颗粒间,不起骨架支撑作用,建模过程采用微分等效介质模型(DEM),将有机质当成包含物的形式添加;当总有机碳质量百分含量(TOC)含量大于2%时,有机质就会存在部分连续分布的状态,对岩石骨架将起支撑作用,有机质就视为岩石骨架的一部分,建模过程中使用自洽模型和微分等效介模型的联合模型(SCA-DEM结合的方式)实现有机质和脆性矿物的混合,从而得到脆性矿物和有机质的混合物的弹性参数。
步骤105,求取脆性矿物和有机质的混合物的弹性参数。
基于脆性矿物混合物和有机质,建立脆性矿物和有机质的混合物的物理模型,获取脆性矿物和有机质的混合物的弹性参数。
作为优选方案,所述建立脆性矿物和有机质的混合物的物理模型包括:
当总有机碳质量百分含量小于2%时,采用微分等效介质模型建立脆性矿物和有机质的混合物的物理模型;
当总有机碳质量百分含量大于2%时,采用自洽模型和微分等效介质模型的联合模型建立脆性矿物和有机质的混合物的物理模型。
步骤106,求取含束缚水的孔隙粘土介质弹性参数;进而求取页岩基质等效弹性参数和求取页岩干岩石骨架等效弹性参数。
在一个示例中,所述求取页岩基质等效弹性参数包括:
基于脆性矿物混合物和有机质,建立脆性矿物和有机质的混合物的物理模型,获取脆性矿物和有机质的混合物的弹性参数;
求取含束缚水的孔隙粘土介质弹性参数;
基于脆性矿物和有机质的混合物的弹性参数和含束缚水的孔隙粘土介质弹性参数,求取页岩基质等效参数。
在一个示例中,所述求取含束缚水的孔隙粘土介质弹性参数包括:
基于粘土颗粒各向异性元和粘土颗粒定向指数,求取粘土介质的弹性参数;
采用各向异性微分等效介质模型,在粘土中添加粘土孔隙和束缚水,得到含束缚水的孔隙粘土介质弹性参数;
将粘土颗粒视为具有固定不变弹性刚度矩阵的各向异性元,引入粘土颗粒定向指数表征粘土颗粒定向排列程度。将页岩中粘土孔隙视为不连通孔隙,采用各向异性微分等效介质模型(各向异性DEM模型),连同其中不可动的束缚水以包含物的形式向泥质中添加,从而得到含束缚水的孔隙粘土介质弹性参数。
在一个示例中,当cdi=1时,粘土介质的弹性参数等于粘土颗粒各向异性元的弹性参数;当cdi=0时,粘土介质的弹性参数为当0﹤cdi<1时,其中,cdi为粘土颗粒定向指数,为粘土颗粒各向异性元的弹性参数,Cv为粘土介质的弹性参数,为cdi=0粘土颗粒完全随机分布时粘土介质的弹性参数。
在一个示例中,所述求取页岩基质等效参数和页岩干岩石骨架等效参数包括:
当粘土矿物含量大于20%时,粘土作为页岩岩石的背景介质,脆性矿物和有机质的混合物作为包含物,采用各向异性微分等效介质模型分别计算页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数;
当粘土矿物含量小于20%时,以脆性矿物和有机质的混合物作为岩石的背景介质,粘土作为包含物,采用各向异性自洽模型和各向异性微分等效介质模型的联合模型分别计算页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数。
步骤107,基于页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数得到流体饱和页岩等效弹性参数,从而完成流体饱和状态下页岩各向异性岩石物理模型的构建。
在一个示例中,采用Brown-Korringa各向异性流体替换理论,基于页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数得到流体饱和页岩等效弹性参数,从而完成流体饱和状态下页岩各向异性岩石物理模型的构建;
Brown-Korringa各向异性流体替换理论由下式给出:
其中,为页岩干岩石骨架等效弹性参数,为页岩孔隙流体饱和状态下的等效弹性参数,页岩基质等效弹性参数,βfl为孔隙流体可压缩性参数,β0为矿物可压缩性参数,φ为孔隙度。
通过本发明的方法所建立的页岩各向异性岩石物理模型,能够更加有效的表征页岩宏观岩石物理性质与其微观结构之间的依赖关系,可较准确地预测页岩的纵波速度、横波速度和各向异性参数,从而为页岩气甜点预测、储层改造和监测等勘探开发技术研究和生产应用提供岩石物理基础支持。
实施例1
图2示出了根据本发明的一个实施例的页岩各向异性岩石物理模型建模方法的建模示意图。
如图2所示,1)首先,由Hashin-Shtrikman边界模型模拟脆性矿物各向同性混合情况下的弹性性质,由于石英、长石、方解石、白云石、黄铁矿等脆性矿物均可视为刚性矿物,弹性性质较为接近,可用该模型进行混合;
2)脆性矿物孔添加:龙马溪组页岩中的脆性矿物分布有粒间孔和粒内溶孔,且多为不连通孔隙,该类孔隙孔径和纵横比通常较大,但含量较少。因此可使用各向同性微分等效介质模型(各向同性DEM模型)加入到脆性矿物混合物中,得到含孔隙脆性矿物混合物;
3)有机质和有机质孔隙的添加:TOC为总有机碳质量百分含量(wt%),而在岩石物理建模等分析中用到的是干酪根的体积含量,需要把TOC转化为体积含量后再加入到模型中去,我们在岩石物理模型中使用由实验得到的经验公式VKerogen(v%)=2.2466×TOC(wt%),将TOC总有机碳质量百分含量(wt%)转换为有机质的体积含量。采用各向同性自洽模型和各向同性微分等效介模型的联合模型(各向同性SCA_DEM模型)在干酪根中添加有机质孔隙,确保有机质中孔隙的连通性。同时,有机质孔占总有机质体积的比重也反映了有机质的成熟度,从而使得模型可以模拟页岩有机质含量和成熟度对页岩岩石物理性质的影响;
4)含孔有机质和含孔脆性矿物的混合物:龙马溪组页岩中脆性矿物以生物成因为主,有机质含量和硅质含量成正相关,且有机质形态多为团块状分散分布。当TOC含量小于一定值时(2%),有机质多分散分布在矿物粒间,不起支撑作用,有机质不作为岩石的骨架,因此建模过程采用微分等效介质模型(DEM)当成包含物的形式添加;当TOC含量大于一定值时(2%),有机质就会存在部分连续分布的状态,对岩石骨架将起支撑作用,有机质就应视为岩石骨架的一部分,建模过程中使用自洽模型和微分等效介模型的联合模型(SCA_DEM联合模型)实现有机质和脆性矿物的混合;
5)粘土颗粒各向异性元:实验分析表明龙马溪组页岩具有较强的各向异性特征,粘土含量和定向性是其各向异性特征产生的主要因素,将页岩中的粘土颗粒视为具有固定不变弹性刚度矩阵的各向异性元,通过在模型中引入各向异性元来实现对粘土矿物的定向排列和颗粒形状所表征的各向异性特征的模拟。该各向异性元表示最微观的粘土结构单元,从各向异性元出发,利用各向异性等效介质理论,得到页岩中粘土整体所表现出的等效弹性性质和各向异性特征。
其中,Cs表示各向异性元的弹性刚度矩阵。
6)粘土定向排列指数:引入粘土颗粒定向指数cdi表征粘土颗粒定向排列程度。cdi=1表示粘土颗粒完全定向,粘土介质的弹性参数等于粘土颗粒各向异性元的弹性参数cdi=0表示粘土颗粒完全随机分布,不具有定向性,粘土介质宏观上表现为各向同性,其弹性参数为当cdi介于0-1之间时,粘土矿物为部分定向,其弹性参数CV由下式给出。
7)粘土矿物孔和孔隙束缚水的添加:粘土矿物的孔隙纵横比相对较小,且多为非连通孔隙,通常填充为不可动的束缚水。因此,将页岩储层中粘土孔隙视为不连通孔隙,采用各向异性DEM模型,连同其中不可动的束缚水以包含物的形式向泥质中添加。
8)粘土颗粒和含孔脆性矿物和含孔有机质的混合物的混合:当粘土矿物含量大于20%时,粘土矿物作为背景介质,脆性矿物和有机质的混合物1作为包含物,采用各向异性DEM模型计算等效参数;当粘土矿物含量小于20%时,以混合物1作为背景介质,粘土作为包含物,采用各向异性SCA_DEM联合模型计算等效参数。
9)采用各向同性微分等效介质模型(各向同性DEM模型)计算脆性矿物和有机质的混合物的弹性参数,进而采用各向异性自洽模型和各向异性微分等效介质模型的联合模型(各向异性SCA-DEM模型)计算页岩基质等效弹性参数。
10)根据储层中油、气、水的实际饱和度,利用Wood公式计算油、气、水混合时的流体弹性参数;根据页岩岩石基质参数、页岩干岩石参数和流体的弹性参数,利用Brown-Korringa各向异性流体替换理论,计算得到流体饱和状态下的页岩等效弹性参数;进而计算出流体饱和状态下页岩的纵、横波速度,以及各向异性参数。
本实施例,在明确了页岩岩石物理性质、各向异性特征及其影响因素的基础上,通过将粘土颗粒视为具有固定不变弹性刚度矩阵的各向异性元,引入粘土颗粒定向指数表征粘土颗粒定向排列程度,实现了对页岩各向异性特征的表征;考虑了页岩多种孔隙类型,有机质的成熟度反映在有机质孔的发育,可实现有机质孔隙发育以及有机质成熟度对页岩岩石物理性质的影响研究;模型中考虑了页岩复杂的矿物组分,通过矿物组分的替换分析,可实现脆性矿物类型及含量对页岩岩石物理性质的影响研究;考虑了页岩不同孔隙类型及形态,可实现孔隙结构对页岩岩石物理参数影响分析;采用BK各向异性流体替换,可实现各向异性介质情况下的流体替换分析,从而可获得流体饱和状态下的页岩岩石物理性质。
通过本实施例所建立的页岩各向异性岩石物理模型,能够更加有效的表征页岩宏观岩石物理性质与其微观结构之间的依赖关系,可较准确地预测页岩的纵波速度、横波速度和各向异性参数,从而为页岩气甜点预测、储层改造和监测等勘探开发技术研究和生产应用提供岩石物理基础支持。
实施例2
本实施例通过页岩岩石基质参数计算、页岩干岩石骨架参数计算和页岩饱和岩石参数计算三大步骤来具体说明本发明的页岩各向异性岩石物理模型建模方法。
图3示出了根据本发明的一个实施例的页岩各向异性岩石物理模型建模方法中求取页岩岩石基质参数的示意图。图4示出了根据本发明的一个实施例的页岩各向异性岩石物理模型建模方法中求取页岩干岩石骨架参数的示意图。图5示出了根据本发明的一个实施例的页岩各向异性岩石物理模型建模方法中求取页岩饱和岩石参数的示意图。其中,HS表示Hashin-Shtrikman边界理论模型,ani-SCA_DEM表示各向异性自洽模型和各向异性微分等效介模型的联合模型,ani-DEM表示各向异性微分等效介模型,TOC表示总有机碳质量百分含量,Kerogen表示有机质,Brown-Korringa公式表示Brown-Korringa各向异性流体替换理论公式,SCA_DEM表示各向同性自洽模型和各向同性微分等效介模型的联合模型。利用Wood公式求取油、气、水等流体混合物的等效弹性参数,流体混合物等效体积参数的计算公式为流体混合物密度的计算公式为fi、Ki和ρi分别是流体各组分的体积分量、体积参数和密度。
图3所示,本实施例提供的页岩各向异性岩石物理模型建模方法中求取页岩岩石基质参数的方法,具体包括以下步骤:
1)由Hashin-Shtrikman边界平均模型计算硅质(石英)、长石、钙质(方解石、白云石)、铁质(黄铁矿)等脆性矿物的脆性矿物混合物弹性参数;
2)使用由实验得到的经验公式VKerogen(v%)=2.2466×TOC(wt%),将TOC总有机碳质量百分含量(wt%)转换为有机质的体积百分含量;
3)计算混合物1(无孔的有机质和脆性矿物的混合物)的等效弹性参数,根据脆性矿物混合物弹性参数和含量,以及有机质弹性参数和含量,当TOC含量小于2%时,采用微分等效介质模型(DEM),将有机质当成包含物的形式添加;当TOC含量大于2%时,使用SCA_DEM结合的方式实现有机质和脆性矿物的混合,从而得到无孔的脆性矿物和有机质的混合物1的弹性参数;
4)将页岩中的粘土颗粒视为一个各向异性元,各向异性元具有固定不变的弹性刚度矩阵,并计算出不同定向指数情况下的页岩粘土弹性参数;
5)采用各向异性微分等效介质理论(aniDEM)模型向粘土中添加含有不可动束缚水的孔隙,得到含孔粘土的弹性参数;
6)由混合物1(即无孔的脆性矿物和有机质的混合物)的弹性参数和含孔粘土的弹性参数,利用各向异性自洽模型和各向异性微分等效介质模型的联合模型(ani-SCA_DEM)计算得到页岩岩石基质的弹性参数。
如图4所示,本实施例提供的页岩各向异性岩石物理模型建模方法中求取页岩干岩石骨架参数的方法,具体包括以下步骤:
1)由Hashin-Shtrikman边界平均模型计算硅质(石英)、长石、钙质(方解石、白云石)、铁质(黄铁矿)等脆性矿物的脆性矿物混合物弹性参数;
2)采用各向同性微分等效介质模型(DEM),将脆性矿物孔加入到脆性矿物混合物中,得到含孔隙脆性矿物混合物的弹性参数;
3)使用由实验得到的经验公式VKerogen(v%)=2.2466×TOC(wt%),将TOC总有机碳质量百分含量(wt%)转换为有机质的体积百分含量;
4)采用各向同性自洽模型和微分等效介质模型的联合模型(SCA_DEM),将有机质孔加入到有机质中,得到含孔隙有机质的弹性参数;
5)含孔混合物1(即含孔脆性矿物和含孔有机质的混合物)弹性等效弹性参数的计算:根据含孔脆性矿物混合物弹性参数和含量,以及含孔有机质弹性参数和含量,把含孔脆性矿物混合物当做背景介质,当TOC含量小于一定值时(2%),采用各向同性微分等效介质模型(DEM)把有机质当做成包含物的形式添加;当TOC含量大于一定值时(2%),建模过程中使用各向同性SCA_DEM结合的方式实现有机质和脆性矿物的混合;
6)将页岩中的粘土颗粒视为一个各向异性元,各向异性元具有固定不变的弹性刚度矩阵,并计算出不同定向指数情况下的页岩粘土弹性参数;
7)采用各向异性微分等效介质理论(ani-DEM)模型向粘土中添加含有不可动束缚水的孔隙,得到含孔粘土的弹性参数;
8)由含孔混合物1(即含孔脆性矿物和含孔有机质的混合物)的弹性参数和含孔粘土的弹性参数,利用各向异性自洽模型和各向异性微分等效介质模型的联合模型(ani-SCA_DEM)计算得到页岩干岩石的弹性参数。
如图5所示,本实施例提供的页岩各向异性岩石物理模型建模方法中求取页岩饱和流体状态下的岩石参数的方法,具体包括以下步骤:
1)根据储层中油、气、水的实际饱和度,利用Wood公式计算油、气、水混合时的流体弹性参数;
2)根据页岩岩石基质参数、页岩干岩石参数和流体的弹性参数,利用Brown-Korringa各向异性流体替换理论,计算得到流体饱和状态下的页岩等效弹性参数;
流体饱和状态下的弹性参数表示为:
3)进而根据弹性参数同速度和各向异性参数的关系式,计算出流体饱和状态下页岩的纵、横波速度,以及各向异性参数。
平行层理纵波速度VPH、垂直层理纵波速度VPV、SH横波速度VSH、SV横波速度VSV等与弹性参数各参数间的关系有下式给出:
ε、γ和δ三个各向异性参数与弹性参数各参数间的关系有下式给出:
测试例1
图6示出了根据本发明的一个实施例的页岩各向异性岩石物理模型模拟结果与页岩样品实验测试结果的对比图。其中,VP为纵波速度,VS为横波速度,ρ为密度,VPV表示垂直方向的纵波速度,VPH表示水平方向上的纵波速度,VSV表示垂直方向的SV横波速度,VSH表示水平方向上的SV横波速度,ε、γ、δ为各向异性参数。
根据本发明实施例1的页岩各向异性岩石物理模型建模方法进行仿真试验验证。利用研发的页岩各向异性岩石物理模型建模方法对页岩岩石物理测试样品进行了模拟试算,如图6所示,由图可知模型预测水平向、垂向的纵波和横波速度、各向异性参数等结果都与岩石物理测试数据有非常好的吻合,从而验证了本发明的页岩各向异性岩石物理模型建模方法的正确性和适用性。
测试例2
图7示出了根据本发明的一个实施例的页岩各向异性岩石物理模型模拟结果与页岩气井测井资料的对比图。
使用测井数据对根据本发明实施例1的页岩各向异性岩石物理模型建模方法进行验证。如图7所示,分析所建页岩各向异性岩石物理模型所预测的纵、横波速度和实测测井数据也都有非常好的吻合,从而也验证了本发明的页岩各向异性岩石模型建模方法的正确性和适用性。
测试例3
图8a-图8f分别示出了本发明的一个实施例的页岩各向异性岩石物理模型模拟有机质含量和有机质孔发育对页岩垂向纵波速度VPv、垂向SV横波速度VSV、Vp/Vs和各向异性参数ε、δ、γ影响的结果图。图9a-图9f分别示出了本发明的一个实施例的页岩各向异性岩石物理模型模拟脆性矿物含量和粘土矿物定向性对页岩垂向纵波速度VPV、垂向SV横波速度VSV、Vp/Vs和各向异性参数ε、δ、γ影响的结果图。其中,VPV为垂向纵波速度,VSV为垂向SV横波速度、Vp/Vs为垂向纵横波速度比,ε、δ、γ表示各向异性参数。
利用本发明实施例1的页岩各向异性岩石物理模型建模方法分析了页岩岩石速度、弹性参数、各向异性参数随有机质含量、有机质孔发育、脆性矿物含量及粘土定向性的变化情况。如图8a-9f所示,通过本发明的方法得出的结论与实验室测量结果有较好的一致性。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种页岩各向异性岩石物理模型建模方法,其特征在于,该建模方法包括:
将页岩基质视为脆性矿物、有机质和粘土组成的混合物;
将粘土颗粒视为具有固定不变弹性刚度矩阵的各向异性元,引入粘土颗粒定向指数表征粘土颗粒定向排列程度;
将页岩总孔隙划分为脆性孔、粘土孔和有机质孔三种孔隙类型;
求取含孔隙脆性矿物混合物的弹性参数和含有有机质孔隙的有机质的弹性参数,进而建立含孔脆性矿物和含孔有机质的混合物的物理模型,获取含孔脆性矿物和含孔有机质的混合物的弹性参数;
求取脆性矿物和有机质的混合物的弹性参数;
求取含束缚水的孔隙粘土介质弹性参数;进而求取页岩基质等效弹性参数和求取页岩干岩石骨架等效弹性参数;
基于页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数得到流体饱和页岩等效弹性参数,从而完成流体饱和状态下页岩各向异性岩石物理模型的构建。
2.根据权利要求1所述的页岩各向异性岩石物理模型建模方法,其中,所述求取页岩基质等效弹性参数包括:
基于脆性矿物混合物和有机质,建立脆性矿物和有机质的混合物的物理模型,获取脆性矿物和有机质的混合物的弹性参数;
求取含束缚水的孔隙粘土介质弹性参数;
基于脆性矿物和有机质的混合物的弹性参数和含束缚水的孔隙粘土介质弹性参数,求取页岩基质等效参数。
3.根据权利要求1所述的页岩各向异性岩石物理模型建模方法,其中,所述脆性矿物包括硅质、长石、钙质和铁质。
4.根据权利要求1所述的页岩各向异性岩石物理模型建模方法,其中,所述求取含孔隙脆性矿物混合物的弹性参数包括:
采用Hashin-Shtrikman边界模型计算脆性矿物各向同性混合情况下的弹性参数上下边界,求上下边界的平均值作为脆性矿物混合物的弹性参数;
利用各向同性微分等效介质模型,将脆性孔加入到所述脆性矿物混合物中,得到含孔隙脆性矿物混合物的弹性参数。
5.根据权利要求1所述的页岩各向异性岩石物理模型建模方法,其中,所述求取含有有机质孔隙的有机质弹性参数包括:
基于经验公式VKerogen(v%)=2.2466×TOC(wt%),将总有机碳质量百分含量转换为有机质的体积含量;
基于各向同性自洽模型和各向同性微分等效介模型的联合模型在有机质中添加有机质孔隙,得到含有有机质孔隙的有机质的弹性参数。
6.根据权利要求5所述的页岩各向异性岩石物理模型建模方法,其中,所述建立含孔脆性矿物和含孔有机质的混合物的物理模型包括:
当总有机碳质量百分含量小于2%时,采用微分等效介质模型建立含孔脆性矿物和含孔有机质的混合物的物理模型;
当总有机碳质量百分含量大于2%时,采用自洽模型和微分等效介质模型的联合模型建立含孔脆性矿物和含孔有机质的混合物的物理模型。
7.根据权利要求1所述的页岩各向异性岩石物理模型建模方法,其中,所述求取含束缚水的孔隙粘土介质弹性参数包括:
基于粘土颗粒各向异性元和粘土颗粒定向指数,求取粘土介质的弹性参数;
采用各向异性微分等效介质模型,在粘土介质中添加粘土孔隙和束缚水,得到含束缚水的孔隙粘土介质弹性参数。
8.根据权利要求7所述的页岩各向异性岩石物理模型建模方法,其中,
当cdi=1时,粘土介质的弹性参数等于粘土颗粒各向异性元的弹性参数;
当cdi=0时,粘土介质的弹性参数为
当0﹤cdi<1时,
其中,cdi为粘土颗粒定向指数,为粘土颗粒各向异性元的弹性参数,Cv为粘土介质的弹性参数,为cdi=0粘土颗粒完全随机分布时粘土介质的弹性参数。
9.根据权利要求1所述的页岩各向异性岩石物理模型建模方法,其中,所述求取页岩基质等效参数和页岩干岩石骨架等效参数包括:
当粘土矿物含量大于20%时,粘土作为页岩岩石的背景介质,脆性矿物和有机质的混合物作为包含物,采用各向异性微分等效介质模型分别计算页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数;
当粘土矿物含量小于20%时,以脆性矿物和有机质的混合物作为岩石的背景介质,粘土作为包含物,采用各向异性自洽模型和各向异性微分等效介质模型的联合模型分别计算页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数。
10.根据权利要求1所述的页岩各向异性岩石物理模型建模方法,其中,采用Brown-Korringa各向异性流体替换理论,基于页岩基质等效弹性参数和页岩干岩石骨架等效弹性参数得到流体饱和页岩等效弹性参数,从而完成流体饱和状态下页岩各向异性岩石物理模型的构建;
Brown-Korringa各向异性流体替换理论由下式给出:
其中,为页岩干岩石骨架等效弹性参数,为页岩孔隙流体饱和状态下的等效弹性参数,页岩基质等效弹性参数,βfl为孔隙流体可压缩性参数,β0为矿物可压缩性参数,φ为孔隙度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710949034.9A CN109655940B (zh) | 2017-10-12 | 2017-10-12 | 页岩各向异性岩石物理模型建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710949034.9A CN109655940B (zh) | 2017-10-12 | 2017-10-12 | 页岩各向异性岩石物理模型建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109655940A true CN109655940A (zh) | 2019-04-19 |
CN109655940B CN109655940B (zh) | 2020-09-29 |
Family
ID=66109833
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710949034.9A Active CN109655940B (zh) | 2017-10-12 | 2017-10-12 | 页岩各向异性岩石物理模型建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109655940B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110133725A (zh) * | 2019-05-08 | 2019-08-16 | 中国石油大学(华东) | 地震岩石横波速度预测方法及装置 |
CN110275206A (zh) * | 2019-08-12 | 2019-09-24 | 河海大学 | 一种裂隙-孔隙型岩石物理弹性模板 |
CN110320575A (zh) * | 2019-06-24 | 2019-10-11 | 中国石油大学(北京) | 基于岩石物理模型的页岩有机质含量确定方法及装置 |
CN110376347A (zh) * | 2019-07-12 | 2019-10-25 | 太原理工大学 | 一种煤的研究模型及其构建方法 |
CN110488386A (zh) * | 2019-09-20 | 2019-11-22 | 西南石油大学 | 一种基于页岩晶体几何因子取向函数的各向异性岩石物理标定方法 |
CN110515126A (zh) * | 2019-09-12 | 2019-11-29 | 中国石油大学(华东) | 一种含随机分布裂缝横向各向同性岩石的声速计算方法 |
CN110909486A (zh) * | 2019-12-06 | 2020-03-24 | 西南石油大学 | 一种正交各向异性页岩岩石物理模型的建立方法 |
CN110909487A (zh) * | 2019-12-06 | 2020-03-24 | 西南石油大学 | 一种精细各向异性页岩岩石物理模型的建立方法 |
CN111353237A (zh) * | 2020-03-23 | 2020-06-30 | 河海大学 | 一种基于矿物晶粒定向发育的各向异性岩石建模方法 |
CN112082918A (zh) * | 2020-08-04 | 2020-12-15 | 中国石油大学(北京) | 孔隙度的确定方法、装置及设备 |
CN112925022A (zh) * | 2021-01-28 | 2021-06-08 | 大庆油田有限责任公司 | 页岩vti介质各向异性参数的预测方法 |
CN113960659A (zh) * | 2021-10-14 | 2022-01-21 | 中国矿业大学 | 一种地震岩石物理驱动的煤层气储层含气量预测方法 |
WO2022198363A1 (zh) * | 2021-03-22 | 2022-09-29 | 中国石油大学(华东) | 一种预测页岩储层弹性参数的方法、装置及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103792581A (zh) * | 2013-12-31 | 2014-05-14 | 刘致水 | 基于岩石物理模型的页岩气储层地应力测井预测方法 |
CN105095631A (zh) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | 一种页岩各向异性岩石物理建模方法 |
CN105203739A (zh) * | 2015-10-29 | 2015-12-30 | 中国石油大学(北京) | 一种基于岩石物理模型的页岩地层各向异性参数预测方法 |
-
2017
- 2017-10-12 CN CN201710949034.9A patent/CN109655940B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103792581A (zh) * | 2013-12-31 | 2014-05-14 | 刘致水 | 基于岩石物理模型的页岩气储层地应力测井预测方法 |
CN105095631A (zh) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | 一种页岩各向异性岩石物理建模方法 |
CN105203739A (zh) * | 2015-10-29 | 2015-12-30 | 中国石油大学(北京) | 一种基于岩石物理模型的页岩地层各向异性参数预测方法 |
Non-Patent Citations (7)
Title |
---|
LIU XI-WU ET AL.: "Anisotropy rock physics model for the Longmaxi shale gas reservoir, Sichuan Basin, China", 《APPLIED GEOPHYSICS》 * |
YUCHONG PEI ET AL.: "Influence of water saturation on acoustic properties of two-phase medium based artificial cores", 《2016 SEG INTERNATIONAL EXPOSITION AND 86TH ANNUAL MEETING》 * |
ZHIQI GUO ET AL.: "A shale rock physics model for analysis of brittleness index, mineralogy and porosity in the Barnett Shale", 《JOURNAL OF GEOPHYSICS AND ENGINEERING》 * |
ZHIQI GUO ET AL.: "Rock physics modeling and anisotropy parameters inversion for shale reservoirs", 《2017 CGS/SEG INTERNATIONAL GEOPHYSICAL CONFERENCE》 * |
印兴耀等: "储层地震岩石物理建模研究现状与进展", 《石油物探》 * |
胡起等: "基于各向异性岩石物理模型的页岩气储层横波速度预测", 《石油物探》 * |
葛瑞 马沃可等: "《岩石物理手册 孔隙介质中地震分析工具》", 31 March 2008, 合肥:中国科学技术大学出版社 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110133725B (zh) * | 2019-05-08 | 2021-05-14 | 中国石油大学(华东) | 地震岩石横波速度预测方法及装置 |
CN110133725A (zh) * | 2019-05-08 | 2019-08-16 | 中国石油大学(华东) | 地震岩石横波速度预测方法及装置 |
CN110320575A (zh) * | 2019-06-24 | 2019-10-11 | 中国石油大学(北京) | 基于岩石物理模型的页岩有机质含量确定方法及装置 |
CN110320575B (zh) * | 2019-06-24 | 2020-05-15 | 中国石油大学(北京) | 基于岩石物理模型的页岩有机质含量确定方法及装置 |
CN110376347A (zh) * | 2019-07-12 | 2019-10-25 | 太原理工大学 | 一种煤的研究模型及其构建方法 |
CN110275206A (zh) * | 2019-08-12 | 2019-09-24 | 河海大学 | 一种裂隙-孔隙型岩石物理弹性模板 |
CN110515126A (zh) * | 2019-09-12 | 2019-11-29 | 中国石油大学(华东) | 一种含随机分布裂缝横向各向同性岩石的声速计算方法 |
CN110488386A (zh) * | 2019-09-20 | 2019-11-22 | 西南石油大学 | 一种基于页岩晶体几何因子取向函数的各向异性岩石物理标定方法 |
CN110488386B (zh) * | 2019-09-20 | 2022-03-25 | 西南石油大学 | 一种基于页岩晶体几何因子取向函数的各向异性岩石物理标定方法 |
CN110909486A (zh) * | 2019-12-06 | 2020-03-24 | 西南石油大学 | 一种正交各向异性页岩岩石物理模型的建立方法 |
CN110909487A (zh) * | 2019-12-06 | 2020-03-24 | 西南石油大学 | 一种精细各向异性页岩岩石物理模型的建立方法 |
CN110909487B (zh) * | 2019-12-06 | 2020-06-23 | 西南石油大学 | 一种精细各向异性页岩岩石物理模型的建立方法 |
CN111353237A (zh) * | 2020-03-23 | 2020-06-30 | 河海大学 | 一种基于矿物晶粒定向发育的各向异性岩石建模方法 |
CN111353237B (zh) * | 2020-03-23 | 2022-03-08 | 河海大学 | 一种基于矿物晶粒定向发育的各向异性岩石建模方法 |
CN112082918B (zh) * | 2020-08-04 | 2021-11-09 | 中国石油大学(北京) | 孔隙度的确定方法、装置及设备 |
CN112082918A (zh) * | 2020-08-04 | 2020-12-15 | 中国石油大学(北京) | 孔隙度的确定方法、装置及设备 |
CN112925022A (zh) * | 2021-01-28 | 2021-06-08 | 大庆油田有限责任公司 | 页岩vti介质各向异性参数的预测方法 |
CN112925022B (zh) * | 2021-01-28 | 2022-03-29 | 大庆油田有限责任公司 | 页岩vti介质各向异性参数的预测方法 |
WO2022198363A1 (zh) * | 2021-03-22 | 2022-09-29 | 中国石油大学(华东) | 一种预测页岩储层弹性参数的方法、装置及存储介质 |
CN113960659A (zh) * | 2021-10-14 | 2022-01-21 | 中国矿业大学 | 一种地震岩石物理驱动的煤层气储层含气量预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109655940B (zh) | 2020-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109655940B (zh) | 页岩各向异性岩石物理模型建模方法 | |
Chen et al. | Methane diffusion in shales with multiple pore sizes at supercritical conditions | |
Naveen et al. | Sorption kinetics of CH4 and CO2 diffusion in coal: theoretical and experimental study | |
Wei et al. | Geometrical, fractal and hydraulic properties of fractured reservoirs: A mini-review | |
Chen et al. | Permeability prediction of shale matrix reconstructed using the elementary building block model | |
CN110909487B (zh) | 一种精细各向异性页岩岩石物理模型的建立方法 | |
Feng et al. | Apparent permeability model for shale oil with multiple mechanisms | |
CN108399270B (zh) | 一种确定页岩地层中各向异性泥质比例的方法 | |
Niu et al. | Influencing factors and selection of CH4 and CO2 adsorption on Silurian shale in Yibin, Sichuan province of China | |
CN103510933B (zh) | 一种渗吸剂渗吸采油效果的定量评价方法 | |
Guo et al. | Adsorption and absorption of supercritical methane within shale kerogen slit | |
Czerw et al. | Application of the stretched exponential equation to sorption of mine gases and sorption induced swelling of bituminous coal | |
CN110909486B (zh) | 一种正交各向异性页岩岩石物理模型的建立方法 | |
Yang et al. | Recent advances in multiscale digital rock reconstruction, flow simulation, and experiments during shale gas production | |
CN109581490A (zh) | 一种致密储层的地震波速度预测方法及系统 | |
Fu et al. | Abrupt changes in reservoir properties of low-rank coal and its control factors for methane adsorbability | |
Dong et al. | Simulating coal permeability change as a function of effective stress using a microscale digital rock model | |
Chen et al. | A new approach to calculate gas saturation in shale reservoirs | |
Kumar et al. | Experimental and numerical evaluation of CBM potential in Jharia Coalfield India | |
Bai et al. | Controlling effect of shale pore connectivity on the adsorption state of shale gas: A case study of continental shale of Jurassic Ziliujing Formation, Sichuan Basin | |
Díaz-Curiel et al. | New granulometric expressions for estimating permeability of granular drainages | |
Miao et al. | Numerical investigation of methane seepage behaviour in coal with lattice Boltzmann approach: The synergistic effects of oxidizing acid and ionic liquid | |
El Amin et al. | Comparative study of shale-gas production using single-and dual-continuum approaches | |
Kim et al. | Improved capillary pressure model considering dual-pore geometry system in carbonate reservoirs | |
Mu et al. | Water distribution in marine shales: based on two-dimensional nuclear magnetic resonance and low-temperature nitrogen adsorption |
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 |