CN104007485B - 一种建立复杂孔隙泥页岩岩石物理模型方法 - Google Patents
一种建立复杂孔隙泥页岩岩石物理模型方法 Download PDFInfo
- Publication number
- CN104007485B CN104007485B CN201410153372.8A CN201410153372A CN104007485B CN 104007485 B CN104007485 B CN 104007485B CN 201410153372 A CN201410153372 A CN 201410153372A CN 104007485 B CN104007485 B CN 104007485B
- Authority
- CN
- China
- Prior art keywords
- kerogen
- hole
- modulus
- mud shale
- model
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 239000011148 porous material Substances 0.000 claims abstract description 44
- 239000011159 matrix material Substances 0.000 claims abstract description 25
- 239000000203 mixture Substances 0.000 claims abstract description 13
- 239000007787 solid Substances 0.000 claims abstract description 9
- 239000011435 rock Substances 0.000 claims description 28
- 238000010008 shearing Methods 0.000 claims description 16
- 239000012530 fluid Substances 0.000 claims description 9
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 7
- 239000011707 mineral Substances 0.000 claims description 7
- 239000000470 constituent Substances 0.000 claims description 5
- 229920006395 saturated elastomer Polymers 0.000 claims description 5
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 239000005416 organic matter Substances 0.000 description 4
- 238000001615 p wave Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005211 surface analysis Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明一种建立复杂孔隙泥页岩岩石物理模型方法，基于改进的Xu-Payne理论的复杂孔隙泥页岩岩石物理模型，该模型通过改进Xu-Payne模型针对碳酸盐储层复杂孔隙特征分布，引入有机质分布，并将孔隙类型分为基质孔隙和干酪根孔隙，基质孔隙分为粒内孔、粒间孔和裂缝，分别用孔隙纵横比0.8、0.15、0.01来表征，干酪根孔隙为微裂隙且充填油气，利用孔隙纵横比0.01来表征，干酪根与油气混合物通过Brown-Korringa方程来实现固体替代，并且假设其纵横比为0.8，从而建立考虑有机质与复杂孔隙类型影响的泥页岩岩石物理模型。
Description
技术领域
本发明属于岩石物理技术领域，特别涉及一种建立复杂孔隙泥页岩岩石物理模型的方法，是一种基于改进Xu-Payne理论的考虑了复杂基质孔隙结构及干酪根微裂隙影响建立的可以预测含有机质饱和流体岩石弹性张量的岩石物理模型。
背景技术
泥页岩由于其复杂的孔隙特征，以及有机质的分布形态及特征，传统的岩石物理模型既没有考虑复杂的基质孔隙形态，更没有考虑有机质的存在及其分布特征，实验表明有机质的存在会显著影响岩石的弹性性质，且有机质的分布形态及特征也是评价页岩气气源条件的重要因素，而且镜下观测表明有机质可能以充填物的方式存在于孔隙中，且干酪根中也存在充填油气的微裂隙，迫切需要研发一种能够模拟有机质及复杂孔隙响应的岩石物理模型。
发明内容
本发明的目的在于提供一种建立复杂孔隙泥页岩岩石物理模型方法，通过改进Xu-Payne模型针对碳酸盐储层复杂孔隙特征分布，引入有机质分布，并将孔隙类型分为基质孔隙和干酪根孔隙，基质孔隙分为粒内孔、粒间孔和裂缝，分别用孔隙纵横比0.8、0.15、0.01来表征，干酪根孔隙为微裂隙且充填油气，利用孔隙纵横比0.01来表征，干酪根与油气混合物通过Brown-Korringa方程来实现固体替代，并且假设其纵横比为0.8，从而建立考虑有机质与复杂孔隙类型影响的泥页岩岩石物理模型。
为了实现上述目的，本发明一种建立复杂孔隙泥页岩岩石物理模型方法，包括以下步骤：
步骤1，结合测井曲线，岩心数据及地质地化报告，得到组成泥页岩各组分的含量、弹性参数、孔隙度、有机质分布特征；
步骤2，根据步骤1得到的矿物组分的体积分数及弹性参数，通过Voigt-Reuss-Hill平均公式（1952）得到岩石骨架的体积模量和剪切模量；
步骤3，将步骤1得到的孔隙度分为基质孔隙和干酪根孔隙，通过DEM（DifferentialEffectiveMedium）模型得到干酪根与干酪根孔隙的混合物，设定干酪根孔隙的纵横比为0.01，并通过Gassmann（1951）方程在微裂隙中充满油气；
步骤4，通过DEM模型向岩石骨架中加入与干酪根与油气混合物所占体积相等的孔隙，设定其纵横比为0.8，并通过Brown-Korringa（1975）方程实现固体替代；
步骤5，利用DEM模型向前面形成的混合物中加入基质孔隙，分为粒内孔、粒间孔、裂缝，得到干岩石的体积模量与剪切模量；
步骤6，通过Gassmann（1951）方程进行流体替代，计算饱和流体岩石的体积模量与剪切模量；
步骤7，利用计算的体积模量与剪切模量，结合计算的总密度和岩石速度之间的关系，得到泥页岩的纵波速度和横波速度，并与实测数据进行对比，判断是否吻合，如果吻合，结束，如果不吻合，返回步骤5，调整基质孔隙中不同孔隙形态所占的体积分数，直至吻合。
本发明与已有技术相比达到的有益效果是：
基于改进的Xu-Payne（2009）理论的复杂孔隙泥页岩岩石物理模型，该模型通过改进Xu-Payne模型针对碳酸盐储层复杂孔隙特征分布，引入有机质分布，并将孔隙类型分为基质孔隙和干酪根孔隙，基质孔隙分为粒内孔、粒间孔和裂缝，分别用孔隙纵横比0.8、0.15、0.01来表征，干酪根孔隙为微裂隙且充填油气，利用孔隙纵横比0.01来表征，干酪根与油气混合物通过Brown-Korringa方程来实现固体替代，并且假设其纵横比为0.8，从而建立考虑有机质与复杂孔隙类型影响的泥页岩岩石物理模型。
附图说明
图1是基于改进Xu-Payne理论的复杂孔隙泥页岩岩石物理模型建模流程图。
图2是根据测井曲线计算得到的干酪根与地化分析得到的干酪根数据准确度分析图。
图3是根据本模型预测得到的纵波速度与实测纵波速度准确度分析图。
图4（a）是含有机质岩石纵波阻抗随干酪根体积分数的变化示意图。
图4（b）是含有机质岩石纵波阻抗与横波阻抗随干酪根体积分数的变化示意图。
图5（a）是含有机质岩石纵波阻抗与纵横波速度比随干酪根体积分数的变化示意图。
图5（b）是含有机质岩石杨氏模量与脆性因子随干酪根体积分数的变化示意图。
图6基于改进Xu-Payne理论的复杂孔隙形态泥页岩岩石物理模型得到的页岩气井综合分析评价图。
具体实施方式
以下结合实例与附图说明本发明具体实施方式。
一种建立复杂孔隙泥页岩岩石物理模型方法，包括以下步骤：
步骤1，结合测井曲线，岩心数据及地质地化报告，得到组成泥页岩各组分的含量、弹性参数、孔隙度、有机质分布特征；
步骤2，根据步骤1得到的矿物组分的体积分数及弹性参数，通过Voigt-Reuss-Hill平均公式（1952）得到岩石骨架的体积模量和剪切模量；
步骤3，将步骤1得到的孔隙度分为基质孔隙和干酪根孔隙，通过DEM（DifferentialEffectiveMedium）模型得到干酪根与干酪根孔隙的混合物，设定干酪根孔隙的纵横比为0.01，并通过Gassmann（1951）方程在微裂隙中充满油气；
步骤4，通过DEM模型向岩石骨架中加入与干酪根与油气混合物所占体积相等的孔隙，设定其纵横比为0.8，并通过Brown-Korringa（1975）方程实现固体替代；
步骤5，利用DEM模型向前面形成的混合物中加入基质孔隙，分为粒内孔、粒间孔、裂缝，得到干岩石的体积模量与剪切模量；
步骤6，通过Gassmann（1951）方程进行流体替代，计算饱和流体岩石的体积模量与剪切模量；
步骤7，利用计算的体积模量与剪切模量，结合计算的总密度和岩石速度之间的关系，得到泥页岩的纵波速度和横波速度，并与实测数据进行对比，判断是否吻合，如果吻合，结束，如果不吻合，返回步骤5，调整基质孔隙中不同孔隙形态所占的体积分数，直至吻合。
实施例中，所述吻合的误差范围为10%。
实施例中，当实际测井曲线中只有纵波测井曲线时，用纵波速度进行约束。
实施例中，步骤3所述的孔隙度分为基质孔隙与干酪根孔隙，基质孔隙分为粒内孔、粒间孔、裂缝，分别用纵横比0.8、0.15、0.01来表征，干酪根孔隙为充填其中的微裂隙，设定纵横比为0.01，通过DEM模型实现干酪根与微裂隙的混合，并通过Gassmann实现微裂隙中充满油气。
实施例中，步骤4所述干酪根与油气混合物与基质矿物的混合，通过DEM模型加入干孔隙，设定纵横比为0.8，并通过Brown-Korringa方程来实现固体替代，加入干酪根与油气混合物。
以上具体实施方式仅用于说明本发明，而非用于限定本发明。
下面通过一口页岩气井纵横波速度预测与优势层段特征分析来说明：
图1是基于改进Xu-Payne理论的复杂孔隙泥页岩岩石物理模型建模流程图。
步骤1，结合测井曲线，岩心数据及地质地化报告，得到组成泥页岩各组分的含量、弹性参数、孔隙度、有机质分布特征；
步骤2，根据步骤1得到的矿物组分的体积分数及弹性参数，通过Voigt-Reuss-Hill平均公式（1952）得到岩石骨架的体积模量和剪切模量；
步骤3，将步骤1得到的孔隙度分为基质孔隙和干酪根孔隙，通过DEM（DifferentialEffectiveMedium）模型得到干酪根与干酪根孔隙的混合物，设定干酪根孔隙的纵横比为0.01，并通过Gassmann（1951）方程在微裂隙中充满油气；
步骤4，通过DEM模型向岩石骨架中加入与干酪根与油气混合物所占体积相等的孔隙，设定其纵横比为0.8，并通过Brown-Korringa（1975）方程实现固体替代；
步骤5，利用DEM模型向前面形成的混合物中加入基质孔隙，分为粒内孔、粒间孔、裂缝，得到干岩石的体积模量与剪切模量；
步骤6，通过Gassmann方程进行流体替代，计算饱和流体岩石的体积模量与剪切模量；
步骤7，利用计算的体积模量与剪切模量，结合计算的总密度和岩石速度之间的关系，得到泥页岩的纵波速度和横波速度，并与实测数据进行对比，判断是否吻合（允许误差范围10%），如果吻合，结束，如果不吻合，返回步骤5，调整基质孔隙中不同孔隙形态所占的体积分数，直至吻合（注意，实际测井曲线中往往只有纵波测井曲线，此时可以用纵波速度进行约束）。
图2为通过测井曲线计算，结合地化分析干酪根数据进行约束得到的干酪根体积分数分布特征，其与地化分析数据的相关系数达到0.86，证明了用来评价的干酪根分布特征的可靠性及准确性。
图3为依据测井曲线得到的各种矿物组分分布及干酪根计算数据，根据本文的模型预测得到的纵波速度与实测测井纵波速度的对比，其相关系数达到0.99，误差率在1%，证明了该模型的准确性。
图4与图5为根据计算的干酪根数据与基于本文模型得到的纵横波速度得到的含有机质泥页岩的弹性参数随干酪根体积分数的变化特征，这与前人的研究结果是一致的，侧面说明了本文模型建模思路的正确性。
图6为依据本模型得到的综合数据以及上面分析的干酪根对含有机质岩石弹性特征影响的分析，优选的页岩气井发育较好层段，从微电阻率扫描成像上看裂缝很发育，这与气测结果及实际压裂结果是一致的。
Claims (5)
1.一种建立复杂孔隙泥页岩岩石物理模型方法，其特征在于，包括以下步骤：
步骤1，结合测井曲线，岩心数据及地质地化报告，得到组成泥页岩各组分的含量、弹性参数、孔隙度、有机质分布特征；
步骤2，根据步骤1得到的矿物组分的体积分数及弹性参数，通过Voigt-Reuss-Hill平均公式-1952得到岩石骨架的体积模量和剪切模量；
步骤3，将步骤1得到的孔隙度分为基质孔隙和干酪根孔隙，通过DEM微分等效介质模型得到干酪根与干酪根孔隙的混合物，设定干酪根孔隙的纵横比为0.01，并通过Gassmann-1951方程在微裂隙中充满油气；
步骤4，通过DEM模型向岩石骨架中加入与干酪根与油气混合物所占体积相等的孔隙，设定其纵横比为0.8，并通过Brown-Korringa-1975方程实现固体替代；
步骤5，利用DEM模型向前面形成的混合物中加入基质孔隙，分为粒内孔、粒间孔、裂缝，得到干岩石的体积模量与剪切模量；
步骤6，通过Gassmann-1951方程进行流体替代，计算饱和流体岩石的体积模量与剪切模量；
步骤7，利用计算的体积模量与剪切模量，结合计算的总密度和岩石速度之间的关系，得到泥页岩的纵波速度和横波速度，并与实测数据进行对比，判断是否吻合，如果吻合，结束，如果不吻合，返回步骤5，调整基质孔隙中不同孔隙形态所占的体积分数，直至吻合。
2.根据权利要求1所述的方法，其特征在于，所述吻合的误差范围为10%。
3.根据权利要求1所述的方法，其特征在于，当实际测井曲线中只有纵波测井曲线时，用纵波速度进行约束。
4.根据权利要求1所述的方法，其特征在于，步骤3所述的孔隙度分为基质孔隙与干酪根孔隙，基质孔隙分为粒内孔、粒间孔、裂缝，分别用纵横比0.8、0.15、0.01来表征，干酪根孔隙为充填其中的微裂隙，设定纵横比为0.01，通过DEM模型实现干酪根与微裂隙的混合，并通过Gassmann实现微裂隙中充满油气。
5.根据权利要求1所述的方法，其特征在于，步骤4所述干酪根与油气混合物与基质矿物的混合，通过DEM模型加入干孔隙，设定纵横比为0.8，并通过Brown-Korringa方程来实现固体替代，加入干酪根与油气混合物。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410153372.8A CN104007485B (zh) | 2014-04-16 | 2014-04-16 | 一种建立复杂孔隙泥页岩岩石物理模型方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410153372.8A CN104007485B (zh) | 2014-04-16 | 2014-04-16 | 一种建立复杂孔隙泥页岩岩石物理模型方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104007485A CN104007485A (zh) | 2014-08-27 |
| CN104007485B true CN104007485B (zh) | 2016-01-20 |
Family
ID=51368228
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410153372.8A Expired - Fee Related CN104007485B (zh) | 2014-04-16 | 2014-04-16 | 一种建立复杂孔隙泥页岩岩石物理模型方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104007485B (zh) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104316977B (zh) * | 2014-11-12 | 2017-02-15 | 中国石油大学(华东) | 一种页岩气储层的岩石模量计算方法 |
| CN106285653A (zh) * | 2015-06-05 | 2017-01-04 | 中国石油化工股份有限公司 | 确定页岩气储层的有机孔隙度的方法 |
| CN105203739A (zh) * | 2015-10-29 | 2015-12-30 | 中国石油大学(北京) | 一种基于岩石物理模型的页岩地层各向异性参数预测方法 |
| CN105301642B (zh) * | 2015-11-02 | 2017-11-07 | 中国石油天然气股份有限公司 | 非均匀孔隙岩石及其固态有机质体积含量确定方法及装置 |
| CN106290105B (zh) * | 2016-07-20 | 2017-08-11 | 中国石油大学(华东) | 一种碳酸盐岩储层溶蚀孔隙体积含量预测方法 |
| CN108957580B (zh) * | 2017-05-22 | 2020-07-03 | 中国石油化工股份有限公司 | 一种反演页岩地层孔隙结构的方法 |
| CN109143397B (zh) * | 2017-06-28 | 2020-05-19 | 中国石油化工股份有限公司 | 碳酸盐岩储层缝洞充填物识别方法及系统 |
| CN109655903B (zh) * | 2017-10-11 | 2020-09-15 | 中国石油化工股份有限公司 | 页岩层横波速度预测方法及系统 |
| CN110552690A (zh) * | 2018-05-30 | 2019-12-10 | 中国石油化工股份有限公司 | 一种页岩储层脆性评价的方法 |
| CN110836836B (zh) * | 2018-08-17 | 2022-08-16 | 中国石油化工股份有限公司 | 干酪根弹性模量测试方法 |
| CN110189353B (zh) * | 2019-06-10 | 2021-01-19 | 中国石油大学(华东) | 一种泥页岩能谱矿物分布图校准方法及系统 |
| CN110471129A (zh) * | 2019-09-20 | 2019-11-19 | 西南石油大学 | 一种深层页岩高温高压下的各向异性岩石物理建模方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110292767A1 (en) * | 2010-05-27 | 2011-12-01 | Jianchun Dai | Determining a quantity of a given material in a subterranean structure |
| CN103713327A (zh) * | 2013-12-31 | 2014-04-09 | 孙赞东 | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 |
| CN103713320A (zh) * | 2013-12-31 | 2014-04-09 | 孙赞东 | 一种富有机质泥页岩岩石物理模型的建立方法 |
-
2014
- 2014-04-16 CN CN201410153372.8A patent/CN104007485B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110292767A1 (en) * | 2010-05-27 | 2011-12-01 | Jianchun Dai | Determining a quantity of a given material in a subterranean structure |
| CN103713327A (zh) * | 2013-12-31 | 2014-04-09 | 孙赞东 | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 |
| CN103713320A (zh) * | 2013-12-31 | 2014-04-09 | 孙赞东 | 一种富有机质泥页岩岩石物理模型的建立方法 |
Non-Patent Citations (2)
| Title |
|---|
| 江苏油田泥灰岩裂缝性储层测井解释方法研究;施振飞;《测井技术》;20040831;第28卷(第1期);第41-44页 * |
| 适用复杂孔隙结构地层消除背景导电饱和度模型;章海宁等;《测井技术》;20110228;第35卷(第1期);第301-305页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104007485A (zh) | 2014-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104007485B (zh) | 一种建立复杂孔隙泥页岩岩石物理模型方法 | |
| CN102445708B (zh) | 三维等效富泥质砂岩速度预测模型 | |
| CN104007482B (zh) | 一种基于各向异性有效场的泥页岩岩石物理模型方法 | |
| CN103713320B (zh) | 一种富有机质泥页岩岩石物理模型的建立方法 | |
| CN105095631B (zh) | 一种页岩各向异性岩石物理建模方法 | |
| CN110824556B (zh) | 一种非常规致密砂岩储层的岩石物理模型建立方法及应用 | |
| CN105089615B (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
| CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
| CN105974474A (zh) | 一种基于孔隙替换的储层孔隙类型评价的方法 | |
| CN102156297B (zh) | 基于砂岩油藏叠后地震数据的流体替换方法 | |
| CN104181585A (zh) | 一种地球物理勘探横波估算方法和系统 | |
| CN105840189A (zh) | 一种页岩气岩石物理与测井方法 | |
| Sun et al. | Advanced petrophysical, geological, geophysical and geomechanical reservoir characterization–key to the successful implementation of a geo-engineered completion optimization program in the eagle ford shale | |
| CN105317429A (zh) | 一种确定页岩地层干酪根含量的方法 | |
| CN109211745A (zh) | 一种富含有机质泥页岩孔隙度的演化过程的恢复方法 | |
| Li et al. | A rock-physical modeling method for carbonate reservoirs at seismic scale | |
| Al Ibrahim et al. | Thermal maturation effects on the elastic properties of organic-rich mudrocks | |
| CN108875115B (zh) | 一种确定岩石强度的方法 | |
| CN110954942A (zh) | 粘土含量识别的图版制作方法及系统 | |
| Montaut et al. | Detection and quantification of rock physics properties for improved hydraulic fracturing in hydrocarbon-bearing shale | |
| Saberi | Rock physics integration: From petrophysics to simulation | |
| CN105866003A (zh) | 一种页岩气储层总孔隙度获取方法及装置 | |
| CN107167844B (zh) | 一种砂泥岩薄互层流体替换方法 | |
| Muhammad et al. | Challenges in Developing Part Field Models to Assess Waterflood Performance for the Upper Burgan Reservoir, Greater Burgan Field, Kuwait | |
| CN107764697A (zh) | 基于孔隙介质渐进方程非线性反演的含气性检测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| ASS | Succession or assignment of patent right |
Owner name: SUN ZANDONG Effective date: 20150626 Owner name: SUN YONGYANG Free format text: FORMER OWNER: SUN ZANDONG Effective date: 20150626 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20150626 Address after: 102249, China University of Petroleum, 18, Xuefu Road, Beijing, Changping District (Beijing) Applicant after: Sun Yongyang Applicant after: Sun Zandong Address before: Beijing city Changping District Road No. 18, China University of Petroleum (Beijing) Applicant before: Sun Zandong |
|
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190314 Address after: 102249 18 Fu Xue Road, Changping District, Beijing Patentee after: China University of Petroleum (Beijing) Address before: 102249 Xuefu Road, Changping District, Beijing, China University of Petroleum (Beijing) Co-patentee before: Sun Zandong Patentee before: Sun Yongyang |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 Termination date: 20210416 |