CN110909487B - 一种精细各向异性页岩岩石物理模型的建立方法 - Google Patents
一种精细各向异性页岩岩石物理模型的建立方法 Download PDFInfo
- Publication number
- CN110909487B CN110909487B CN201911239404.5A CN201911239404A CN110909487B CN 110909487 B CN110909487 B CN 110909487B CN 201911239404 A CN201911239404 A CN 201911239404A CN 110909487 B CN110909487 B CN 110909487B
- Authority
- CN
- China
- Prior art keywords
- shale
- model
- gpa
- dry
- tensor
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种精细各向异性页岩岩石物理模型的建立方法,包括以下步骤:通过测井、录井、室内测试、各向异性SCA和DEM模型、各向同性SCA和DEM模型、Brown‑Korringa模型,从而建立一种可以预测页岩刚度张量的岩石物理模型。本发明分别考虑无机质孔隙和有机质孔隙对建模结果的影响,提出利用吸附气和游离气的含量,结合总孔隙度分别计算有机质孔隙和无机质孔隙的含量,综合考虑了页岩各向异性的成因,与页岩地层的实际物性特征更加相符,并且该方法相对其他方法的建模精度更高,可以为后续的波速预测、各向异性评价、岩石弹性参数评价、地应力评价等提供基础数据。
Description
技术领域
本发明涉及一种精细各向异性页岩岩石物理模型的建立方法,属于岩石物理技术领域。
背景技术
页岩的微观性质会对宏观性质产生重要的影响,例如有机质矿物的定向排列、成层以及岩石中扁平状孔隙(横向发育微裂缝)的存在是导致页岩表现出横观各向同性的主要原因。页岩岩石在平行于层面方向和垂直于层面方向的弹性性质存在较大差异,最明显的就是两个方向上的杨氏模量和泊松比有所不同。常规油气藏(例如砂岩油气藏、碳酸盐岩油气藏)往往表现出各向同性的特性,其弹性参数仅需2个独立弹性刚度系数即可表征。而一般的横观各向同性的页岩弹性特征需要5个独立的弹性刚度系数才能表征。因此,从常规油气藏获取的一些结论和认识,不能照搬照抄至页岩地层。常规的地球物理解释方法对页岩不再适用。目前的测井仪器最多能获取5个独立弹性刚度系数中的3个,剩余2个仍然无能为力。岩石物理方法是架接岩石微观物性特征(岩石矿物组分、孔隙及其填充物、矿物颗粒形状、胶结程度、连通性等)和宏观弹性特征(波速、岩石弹性参数)的桥梁。岩石物理建模则是利用合理的基础岩石物理方法,针对岩石的矿物、孔隙、流体的特点进行描述,以达到表征岩石弹性特征的目的。页岩岩石物理建模为各向异性页岩地层的刚度系数计算、波速各向异性评价以及岩石的弹性特征评价提供了一种有力手段。
为模拟富含有机质页岩的等效弹性性质,以扫描电镜获得的黏土-干酪根矿物定向排列特征为基础,Vernik通过对Backus平均进行改进来模拟页岩的等效弹性性质,尽管修正后的Backus取得了良好的预测效果,但其计算公式为半经验公式,不具有普适性。Hornby等人从各向同性自洽模型(SCA)和微分等效介质模型(DEM)出发,推导了各向异性的SCA和DEM模型,基于扫描电镜,引入高斯分布模拟了页岩中黏土和干酪根的定向排列分布和成层特性。
Bandyopadhyay和Wu将将干酪根视为背景介质,利用V-R-H界限计算除干酪根以外的混合岩石矿物的等效弹性模量,并通过各向异性DEM模型将混合物添加到干酪根中,来模拟富含有机质页岩的弹性属性。通常情况下干酪根含量很少,不足以在地层中形成连续的背景,并且V-R-H界限计算的混合物的等效弹性模量有较大误差,此外该模型模拟出的岩石内部流体处于离散状态。考虑到流体的连通问题,Hu利用各向异性SCA和DEM模拟的黏土和流体混合物作为背景介质,并将其余矿物用V-R-H界限模拟干酪根、石英、方解石等混合物的等效弹性模量,最后用各向异性DEM将混合物添加到背景介质中,由于干酪根和石英等矿物的弹性模量差异较大,因此该模型也存在较大误差;刘子淳针对Wu和Hu模型的局限性,将干酪根作为背景介质,用SCA+DEM模型将孔隙添加到干酪根中以保证孔隙和干酪根的相互连通,利用H-S界限计算除干酪根以外的矿物等效弹性模量并作为包裹体,然后利用DEM添加到背景介质中,最后进行流体替换得到饱和的页岩岩石物理模型,该模型考虑了干酪根及孔隙的相互连通,但没有考虑不同孔隙对建模结果的影响。
发明内容
针对上述问题,本发明主要是克服现有技术中的不足之处,提出一种精细各向异性页岩岩石物理模型的建立方法。
本发明解决上述技术问题所提供的技术方案是:一种精细各向异性页岩岩石物理模型的建立方法,包括以下步骤:
步骤S1、通过测井、录井以及室内测试,得到建模所需的页岩物性参数,所述页岩物性参数包括矿物组成及含量、地层流体组成及含量、孔隙度、饱和度;
步骤S2、利用Hashin-Shtrikman界限计算基质矿物混合物的等效弹性模量;
步骤S3、利用各向同性SCA模型和各向同性DEM模型,将无机质孔隙加入到基质矿物混合物中得到干基质,建立干基质页岩岩石物理模型并计算干基质的等效弹性模量;
步骤S4、利用各向异性SCA模型和各向异性DEM模型计算有机质矿物混合物的等效弹性刚度张量;
步骤S5、利用CL系数描述有机质矿物的成层分布,结合Bond变换得到取向分布的有机质矿物混合物等效弹性刚度张量;
步骤S6、利用各向异性DEM模型,将有机质孔隙加入到取向分布的有机质矿物混合物中得到干有机质,建立干有机质页岩岩石物理模型并计算干有机质的等效弹性刚度张量;
步骤S7、利用各向异性DEM模型,将干基质添加到干有机质中,得到包含空孔隙的干页岩,建立干页岩岩石物理模型并计算干页岩的刚度张量;
步骤S8、根据含气饱和度、含水饱和度,利用Wood公式计算有机质孔隙内气-水混合物的体积模量,并转换成刚度张量;并结合干页岩呈VTI性质的情况,利用Brown-Korringa模型,将混合流体添加到干页岩中得到饱和流体页岩,建立饱和流体页岩岩石物理模型并计算饱和流体页岩的刚度张量。
进一步的技术方案是,所述步骤S2的具体过程是:利用Hashin-Shtrikman界限计算基质矿物混合物体积模量和剪切模量的上下限,并取上下限平均作为基质矿物混合物的体积模量和剪切模量。
进一步的技术方案是,所述Hashin-Shtrikman界限计算基质矿物混合物的公式如下:
式中:KHS+为基质矿物混合物的等效体积模量上限,GPa;KHS-为基质矿物混合物的等效体积模量下限,GPa;μHS+为基质矿物混合物的等效剪切积模量上限,GPa;μHS-为基质矿物混合物的等效剪切模量下限,GPa。
进一步的技术方案是,所述步骤S3中各向同性SCA模型的计算公式如下:
各向同性DEM模型的计算公式如下:
式中:K1、μ1分别为背景介质的体积模量和剪切模量,GPa;K2、μ2分别为包裹体的体积模量和剪切模量,GPa;v为包裹体的体积分数,无因次。
进一步的技术方案是,所述步骤S4中各向异性SCA模型的计算公式如下:
各向异性DEM模型的计算公式如下:
进一步的技术方案是,所述步骤S8中Brown-Korringa模型的计算公式如下:
本发明具有以下有益效果:本发明分别考虑无机质孔隙和有机质孔隙对建模结果的影响,提出利用吸附气和游离气的含量,结合总孔隙度分别计算有机质孔隙和无机质孔隙的含量,综合考虑了页岩各向异性的成因,与页岩地层的实际物性特征更加相符,并且该方法相对其他方法的建模精度更高,可以为后续的波速预测、各向异性评价、岩石弹性参数评价、地应力评价等提供基础数据。
附图说明
图1是本发明的建模流程图;
图2是x1井矿物及流体的测井解释图;
图3是有机质混合物的SCA和DEM联合等效刚度系数图;
图4是x1井岩石物理建模预测结果以及同测井测量结果的对比图。
具体实施方式
下面结合实施例和附图对本发明做更进一步的说明。
如图1所示,本发明的一种精细各向异性页岩岩石物理模型的建立方法,具体实施过程如下:
步骤S1、通过测井、录井以及室内测试,得到建模所需的页岩物性参数,其x1井矿物及流体的测井解释图如图2所示,所述页岩物性参数包括矿物组成及含量、地层流体组成及含量、孔隙度、饱和度;
利用测井、录井以及室内测试等手段,获取页岩的矿物组分及相对含量,矿物组分包括基质矿物和有机质矿物两大类,基质矿物主要包括石英、长石、方解石、白云石、黄铁矿等,有机质矿物主要黏土和干酪根;
再利用室内测试或测井,获取页岩孔隙度,以及孔隙中所含的页岩气及水对应的包含度;页岩孔隙中的气体分为吸附气和游离气,基于Xu-White模型的原理,将孔隙分为无机质孔隙和有机质孔隙两类,其中无机质孔隙包含了粒内孔、粒间孔以及微裂缝;利用游离气和吸附气相对含量计算无机质孔隙度和有机质孔隙度;
吸附气的含量可由兰格缪尔等温吸附方程进行计算,具体公式如下:
式中,ga为吸附气含量,m3/t;Vlc为在油藏温度下经过TOC校正后的兰格缪尔体积,m3/t;Plt为油藏温度下的兰格缪尔压力,Mpa;p为储层压力,Mpa。
游离气的含量可由以下方程进行计算:
式中,gf为游离气含量,m3/t;ψ为常数,无因次;ρ为体积密度,g/cm3;Bg为气体体积压缩系数,无因次;φe为有效孔隙度,无因次;Sw为有效孔隙度,无因次。
有机质孔隙度和无机质孔隙度的计算方式如下:
式中,φo为有机质孔隙度,无因次;φg为无机质孔隙度,无因次;φ为总孔隙度,无因次;ga为吸附气含量,m3/t;gf为游离气含量,m3/t。
步骤S2、利用Hashin-Shtrikman界限计算基质矿物混合物的等效弹性模量;
利用Hashin-Shtrikman界限计算基质矿物混合物体积模量和剪切模量的上下限,并取上下限平均作为基质矿物混合物的体积模量和剪切模量;基质矿物是指除了黏土和干酪根以外的矿物,主要包括:石英、长石、方解石、白云石、黄铁矿等;
所述Hashin-Shtrikman界限计算基质矿物混合物的公式如下:
式中:KHS+为基质矿物混合物的等效体积模量上限,GPa;KHS-为基质矿物混合物的等效体积模量下限,GPa;μHS+为基质矿物混合物的等效剪切积模量上限,GPa;μHS-为基质矿物混合物的等效剪切模量下限,GPa;
本实施例页岩岩石矿物及流体弹性模量如表1所示:
表1
步骤S3、利用各向同性SCA模型和各向同性DEM模型,将无机质孔隙加入到基质矿物混合物中得到干基质,建立干基质页岩岩石物理模型并计算干基质的等效弹性模量;
具体是给定页岩基质岩石的临界孔隙度,利用各向同性SCA模型计算干岩石在临界孔隙度时候的体积模量和剪切模量,然后再利用各向同性DEM模型将孔隙度逐步调整至真实无机质孔隙度φg,以计算真实无机质孔隙度时的干页岩基质的体积模量和剪切模量;
各向同性SCA模型的计算公式如下:
各向同性DEM模型的计算公式如下:
式中:K1、μ1分别为背景介质的体积模量和剪切模量,GPa;K2、μ2分别为包裹体的体积模量和剪切模量,GPa;v为包裹体的体积分数,无因次;
步骤S4、利用各向异性SCA模型和各向异性DEM模型计算有机质矿物混合物的等效弹性刚度张量,其结果图3所示;
具体的是首先利用各向异性SCA模型计算干酪根和黏土含量各占50%情况下的弹性特征,然后用各向异性DEM模型将两者的体积含量调整至对应的百分比,这既保证了干酪根和黏土的相互连通,又避免了由于干酪根和黏土添加顺序的不同而导致混合物弹性模量的不对称;
各向异性SCA模型的计算公式如下:
各向异性DEM模型的计算公式如下:
步骤S5、利用CL系数描述有机质矿物的成层分布,结合Bond变换得到取向分布的有机质矿物混合物等效弹性刚度张量;
S5.1利用扫描电镜资料,确定有机质颗粒随对称轴的偏转角度,给出表征有机质颗粒偏转角度的排列分布系数CL;
S5.2对于偏转到单个角度下的有机质颗粒,利用Bond变换将完全定向排列的有机质颗粒旋转到对应的偏转角度,得到该角度下对应的有机质颗粒的刚度张量;
S5.3将分布函数上所有排列角度的有机质颗粒的刚度矩阵利用V-R-H模型进行叠加,得到不同角度叠加的最终的有机质等效弹性刚度张量;
步骤S6、利用各向异性DEM模型,将有机质孔隙加入到取向分布的有机质矿物混合物中得到干有机质,建立干有机质页岩岩石物理模型并计算干有机质的等效弹性刚度张量;
步骤S7、利用各向异性DEM模型,将干基质添加到干有机质中,得到包含空孔隙的干页岩,建立干页岩岩石物理模型并计算干页岩的刚度张量;
步骤S8、根据含气饱和度、含水饱和度,利用Wood公式计算有机质孔隙内气-水混合物的体积模量,并转换成刚度张量;并结合干页岩呈VTI性质的情况,利用Brown-Korringa模型,将混合流体添加到干页岩中得到饱和流体页岩,建立饱和流体页岩岩石物理模型并计算饱和流体页岩的刚度张量,其结果如图4(c)所示;
Brown-Korringa模型的计算公式如下:
本实施例得到的岩石物理建模预测结果以及同测井测量结果的对比如图4所示;
图4(a)、(b)分别为测井实测的刚度系数和建模预测得到的刚度系数以及实测结果与预测结果之间的误差曲线。经统计,c33的实测结果与建模预测结果之间的误差介于-6.7%-7.1%之间(误差绝对值平均3.9%),c44的实测结果与建模预测结果之间的误差介于-5.1%-5.0%之间(误差绝对值平均2.5%),两个弹性刚度系数的预测结果与实测结果的差异较小,验证了模型的精确性。
以上所述,并非对本发明作任何形式上的限制,虽然本发明已通过上述实施例揭示,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些变动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (6)
1.一种精细各向异性页岩岩石物理模型的建立方法,其特征在于,包括以下步骤:
步骤S1、通过测井、录井以及室内测试,得到建模所需的页岩物性参数,所述页岩物性参数包括矿物组成及含量、地层流体组成及含量、孔隙度、饱和度;
步骤S2、利用Hashin-Shtrikman界限计算基质矿物混合物的等效弹性模量;
步骤S3、利用各向同性SCA模型和各向同性DEM模型,将无机质孔隙加入到基质矿物混合物中得到干基质,建立干基质页岩岩石物理模型并计算干基质的等效弹性模量;
步骤S4、利用各向异性SCA模型和各向异性DEM模型计算有机质矿物混合物的等效弹性刚度张量;
步骤S5、利用CL系数描述有机质矿物的成层分布,结合Bond变换得到取向分布的有机质矿物混合物等效弹性刚度张量;
步骤S6、利用各向异性DEM模型,将有机质孔隙加入到取向分布的有机质矿物混合物中得到干有机质,建立干有机质页岩岩石物理模型并计算干有机质的等效弹性刚度张量;
步骤S7、利用各向异性DEM模型,将干基质添加到干有机质中,得到包含空孔隙的干页岩,建立干页岩岩石物理模型并计算干页岩的刚度张量;
步骤S8、根据含气饱和度、含水饱和度,利用Wood公式计算有机质孔隙内气-水混合物的体积模量,并转换成刚度张量;并结合干页岩呈VTI性质的情况,利用Brown-Korringa模型,将混合流体添加到干页岩中得到饱和流体页岩,建立饱和流体页岩岩石物理模型并计算饱和流体页岩的刚度张量。
2.根据权利要求1所述的一种精细各向异性页岩岩石物理模型的建立方法,其特征在于,所述步骤S2的具体过程是:利用Hashin-Shtrikman界限计算基质矿物混合物体积模量和剪切模量的上下限,并取上下限平均作为基质矿物混合物的体积模量和剪切模量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911239404.5A CN110909487B (zh) | 2019-12-06 | 2019-12-06 | 一种精细各向异性页岩岩石物理模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911239404.5A CN110909487B (zh) | 2019-12-06 | 2019-12-06 | 一种精细各向异性页岩岩石物理模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110909487A CN110909487A (zh) | 2020-03-24 |
CN110909487B true CN110909487B (zh) | 2020-06-23 |
Family
ID=69823117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911239404.5A Active CN110909487B (zh) | 2019-12-06 | 2019-12-06 | 一种精细各向异性页岩岩石物理模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110909487B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111537345B (zh) * | 2020-05-12 | 2021-08-10 | 西南石油大学 | 一种横观各向同性岩样抗压剪强度参数的快速测定方法 |
CN111859632B (zh) * | 2020-07-01 | 2023-09-22 | 广州海洋地质调查局 | 一种水合物储层的岩石物理模型构建方法及处理终端 |
CN114594531A (zh) * | 2020-12-07 | 2022-06-07 | 中国石油化工股份有限公司 | 一种可动油的陆相页岩油地震岩石物理建模方法 |
CN112525803A (zh) * | 2021-02-10 | 2021-03-19 | 西南石油大学 | 基于黏土壳状模型的砂岩孔隙度有效应力系数计算方法 |
CN113642157B (zh) * | 2021-07-20 | 2024-05-31 | 西南石油大学 | 一种基于高斯过程回归的页岩孔隙结构参数预测方法 |
CN113960659B (zh) * | 2021-10-14 | 2023-02-03 | 中国矿业大学 | 一种地震岩石物理驱动的煤层气储层含气量预测方法 |
CN114371187B (zh) * | 2021-12-17 | 2023-11-03 | 成都理工大学 | 页岩孔隙定量结构复杂性测定方法、系统、介质、设备 |
CN118553325A (zh) * | 2024-07-30 | 2024-08-27 | 中国石油大学(华东) | 含裂隙水合物储层各向异性电学性质模拟方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095631A (zh) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | 一种页岩各向异性岩石物理建模方法 |
CN109655940A (zh) * | 2017-10-12 | 2019-04-19 | 中国石油化工股份有限公司 | 页岩各向异性岩石物理模型建模方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7676349B2 (en) * | 2004-12-06 | 2010-03-09 | Exxonmobil Upstream Research Co. | Integrated anisotropic rock physics model |
CN108983312B (zh) * | 2018-05-25 | 2019-09-24 | 中国海洋石油集团有限公司 | 一种裂缝型储层各向异性岩石模量估算方法 |
-
2019
- 2019-12-06 CN CN201911239404.5A patent/CN110909487B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095631A (zh) * | 2014-05-21 | 2015-11-25 | 中国石油化工股份有限公司 | 一种页岩各向异性岩石物理建模方法 |
CN109655940A (zh) * | 2017-10-12 | 2019-04-19 | 中国石油化工股份有限公司 | 页岩各向异性岩石物理模型建模方法 |
Non-Patent Citations (3)
Title |
---|
"Prediction of brittleness based on anisotropic rock physics model for kerogen-rich shale";Qian Ke-Ran et al.;《APPLIED GEOPHYSICS》;20171231;第14卷(第4期);第463-479页 * |
"包裹体理论在富有机质页岩岩石物理建模中的依据因素应用";刘子淳 等;《地球物理学进展》;20171231;第32卷(第6期);第2505-2512页 * |
"泥页岩岩石物理建模研究";董宁 等;《地球物理学报》;20140615;第57卷(第6期);第1990-1998页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110909487A (zh) | 2020-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110909487B (zh) | 一种精细各向异性页岩岩石物理模型的建立方法 | |
CN111460601B (zh) | 基于岩石物理建模的正交各向异性地层地应力预测方法 | |
CN109655940B (zh) | 页岩各向异性岩石物理模型建模方法 | |
Zhang et al. | Flow mechanism and simulation approaches for shale gas reservoirs: A review | |
Liu et al. | A new theoretical approach to model sorption-induced coal shrinkage or swelling | |
CN111460602B (zh) | 基于岩石物理建模的横观各向同性地层地应力预测方法 | |
CN105927218B (zh) | 一种陆相泥页岩储层含气量预测方法及装置 | |
Yang et al. | Review of shale gas sorption and its models | |
Wang et al. | Experimental and numerical investigation of dynamic gas adsorption/desorption–diffusion process in shale | |
CN110909486B (zh) | 一种正交各向异性页岩岩石物理模型的建立方法 | |
Fan et al. | Recent advances and perspectives of CO2-enhanced coalbed methane: experimental, modeling, and technological development | |
Pang et al. | Analysis of effective porosity and effective permeability in shale-gas reservoirs with consideration of gas adsorption and stress effects | |
Gou et al. | Quantitative calculated shale gas contents with different lithofacies: A case study of Fuling gas shale, Sichuan Basin, China | |
Luan et al. | Creation of synthetic samples for physical modelling of natural shale | |
Shafer | Recent advances in core analysis | |
CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
CN104007482A (zh) | 一种基于各向异性有效场的泥页岩岩石物理模型方法 | |
Liu et al. | Experimental study of the effective stress coefficient for coal anisotropic permeability | |
Zheng et al. | Dynamic fluid interactions during CO2-enhanced coalbed methane and CO2 sequestration in coal seams. Part 1: CO2–CH4 interactions | |
Yang et al. | Comparative investigations on wettability of typical marine, continental, and transitional shales in the middle Yangtze Platform (China) | |
Ding et al. | Competitive adsorption between CO2 and CH4 in tight sandstone and its influence on CO2-injection enhanced gas recovery (EGR) | |
Ren et al. | Control mechanism and parameter simulation of oil-water properties on spontaneous imbibition efficiency of tight sandstone reservoir | |
Berawala et al. | Evaluation of multicomponent adsorption kinetics for carbon dioxide enhanced gas recovery from tight shales | |
Chen et al. | Shale gas transport behavior considering dynamic changes in effective flow channels | |
Mo et al. | Influence of water on gas transport in shale nanopores: Pore-scale simulation study |
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 |