本发明涉及岩石力学技术领域，特别涉及一种岩石力学与储层工程参数评估方法。
背景技术：
充分了解深部储层岩石的力学特性对深部资源高效开发具有重要意义。工程上，通常采用室内试验测试岩芯的基本力学参数。随着储层逐渐向深部发展，地质构造环境复杂，取芯作业变得十分困难，也就难以通过室内试验获取岩石力学参数，造成钻井和压裂作业具有一定的盲目性，导致储层产能不理想等后果。但是，钻井过程中排出的岩屑包含了储层微观尺度的矿物成分、矿物颗粒结构和矿物微观力学参数等特征基因，且这些特征基因决定了岩石的宏观力学特性。因此，应充分利用随钻排出的岩屑，在准确测定微观特征基因的基础上有效评估储层岩石的宏观力学性能，并在深井勘探开发中应用，解决深部岩石力学参数难以获取的难题，为深部资源高效开发提供科学依据。
现有技术中多采用地表露头岩芯开展室内试验，对实际井下岩石力学参数进行估算，或通过测井信息进行评估，均存在较大的不确定性、估算结果精度差等问题。
技术实现要素：
本发明提供一种岩石力学与储层工程参数评估方法，解决了或部分解决了现有技术中深部储层取芯费时费力、无法准确获取储层工程参数的技术问题。
为解决上述技术问题，本发明提供了一种岩石力学与储层工程参数评估方法包括：获取深井随钻排出的岩屑；对所述岩屑进行x射线衍射，获取所述岩屑的矿物成分及含量；对所述岩屑进行ct扫描，并在显微镜下进行观察，获取所述岩屑的各矿物颗粒的微观结构特征；对所述岩屑进行微/纳米压痕试验，测定所述岩屑的各矿物成分的微观力学参数；由所述岩屑的矿物成分及含量、所述岩屑的各矿物颗粒的微观结构特征及所述岩屑的各矿物成分的微观力学参数获取岩屑微观特征基因；根据所述岩屑微观特征基因构建三维计算模型，取代表性体积单元得到岩石宏观力学参数；根据岩石宏观力学参数及工程参数评估模型对现场需求的工程参数进行评估。
进一步地，所述对岩屑进行x射线衍射包括：将所述岩屑研磨成200目。
进一步地，所述对岩屑进行ct扫描，并在显微镜下进行观察包括：将所述岩屑进行微米ct扫描，将ct扫描切片进行三维重构，形成三维数字岩芯，获取所述岩屑的矿物颗粒的微观结构特征；将所述岩屑制成薄片，放置于偏光显微镜下观察，确定所述岩屑的矿物颗粒的微观结构特征，与ct扫描结果互为补充。
进一步地，所述薄片的厚度为28-32μm。
进一步地，所述对岩屑进行微/纳米压痕试验包括：在压痕试验机的显微镜下区分所述岩屑的各矿物的分布区域，对每种矿物进行单独测试，获取每种矿物成分的微观力学参数。
进一步地，所述岩屑的矿物成分的微观力学参数包括：弹性模量、硬度及断裂韧性。
进一步地，
进一步地，所述根据岩石宏观力学参数及工程参数评估模型对现场需求的工程参数进行评估包括：结合测井资料，对脆性指标、断裂韧性和可压裂性指标进行评估。
进一步地，所述测井资料包括：储层埋深及岩性。
进一步地，所述根据岩石宏观力学参数及工程参数评估模型对现场需求的工程参数进行评估包括：bi＝σcσt/2；ki＝0.313 0.027e′；
本申请实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：
由于获取深井随钻排出的岩屑，对岩屑进行x射线衍射，获取岩屑的矿物成分及含量，对岩屑进行ct扫描，并在显微镜下进行观察，获取岩屑的各矿物颗粒的微观结构特征，对岩屑进行微/纳米压痕试验，测定所述岩屑的各矿物成分的微观力学参数，由岩屑的矿物成分及含量、岩屑的各矿物颗粒的微观结构特征及岩屑的各矿物成分的微观力学参数获取岩屑微观特征基因，根据岩屑微观特征基因构建三维计算模型，取代表性体积单元得到岩石宏观力学参数，根据岩石宏观力学参数及工程参数评估模型可以准确地评估现场需求的工程参数，对现场钻井和压裂工程设计与优化具有重要意义，不需要深部储层取芯，准确性高，估算精度高。
附图说明
图1为本发明实施例提供的岩石力学与储层工程参数评估方法的流程示意图；
图2为图1中岩屑微观特征基因测定中压痕载荷-位移曲线。
具体实施方式
参见图1，本发明实施例提供的一种岩石力学与储层工程参数评估方法，其特征在于，包括：
获取深井随钻排出的岩屑。
对岩屑进行x射线衍射，获取岩屑的矿物成分及含量。
对岩屑进行ct扫描，并在显微镜下进行观察，获取岩屑的各矿物颗粒的微观结构特征。
对岩屑进行微/纳米压痕试验，测定岩屑的各矿物成分的微观力学参数，获取岩屑微观特征基因。
根据岩屑微观特征基因构建三维计算模型，将岩屑的各矿物成分的微观力学参数赋予模型中离散的材料点，进而选取合适尺寸的代表性体积单元，用于评估储层岩石的宏观力学性能。
本申请具体实施方式由于获取深井随钻排出的岩屑，对岩屑进行x射线衍射，获取岩屑的矿物成分及含量，对岩屑进行ct扫描，并在显微镜下进行观察，获取岩屑的各矿物颗粒的微观结构特征，对岩屑进行微/纳米压痕试验，测定所述岩屑的各矿物成分的微观力学参数，由岩屑的矿物成分及含量、岩屑的各矿物颗粒的微观结构特征及岩屑的各矿物成分的微观力学参数获取岩屑微观特征基因，根据岩屑微观特征基因构建三维计算模型，取代表性体积单元得到岩石宏观力学参数，根据岩石宏观力学参数及工程参数评估模型可以准确地评估现场需求的工程参数，对现场钻井和压裂工程设计与优化具有重要意义，不需要深部储层取芯，准确性高，估算精度高。
其中，构建三维计算模型，将岩屑的各矿物成分的微观力学参数赋予模型中离散的材料点，进而选取合适尺寸的代表性体积单元，初步得到储层岩石的宏观力学性能，包括抗压/拉强度、弹性模量和泊松比。
具体地，对岩屑进行x射线衍射包括：
将岩屑研磨成200目，将200目的岩屑，置于x射线衍射仪中进行主要矿物成分定性和定量分析。
x射线衍射试验，是指从随钻排出的岩屑中筛选出代表性岩屑，并将其辗磨成200目细粉，采用x-射线衍射仪进行定性和定量分析，获取岩屑的主要矿物成分及含量。
具体地，对岩屑进行ct扫描，并在显微镜下进行观察包括：
将岩屑进行微米ct扫描，将ct扫描切片进行三维重构，形成三维数字岩芯，获取所述岩屑的矿物颗粒的微观结构特征。
将岩屑制成薄片，放置于偏光显微镜下观察，确定岩屑的矿物颗粒的微观结构特征，与ct扫描结果互为补充。
薄片的厚度为28-32μm。
ct扫描，是指将岩屑放置于微米ct机内进行高精度扫描，获取岩屑矿物颗粒的微观结构信息。
具体地，对岩屑进行微/纳米压痕试验包括：
在压痕试验机的显微镜下区分所述岩屑的各矿物的分布区域，对每种矿物进行单独测试，获取每种矿物成分的微观力学参数。
参见图2，所述微/纳米压痕试验，是指将岩屑表面抛光后在压痕试验机上进行微观力学参数测试，应根据矿物颗粒结构和尺寸选择合适的压头及加载条件。在压痕试验机的显微镜下区分各矿物的分布区域，对每种矿物进行单独测试。压头，纳米尺度以berkovich(金刚石玻氏压头)压头为主，微米尺度以vicker压头为主。加载条件包括最大力、加载速率、压入深度和饱载时间。微观力学参数主要包括弹性模量、硬度和断裂韧性，均可由压痕载荷—位移曲线计算得出。
岩屑的矿物成分的微观力学参数包括：弹性模量、硬度及断裂韧性。
式中，e为对应矿物的弹性模量；νs为待测样品的泊松比；ei和νi分别为压头的弹性模量和泊松比；er为压痕模量；a(hc)为压头与试样的接触投影面积，与压入深度hc有关；s为接触刚度；β和ε是与压头形状有关的常数；fmax和hmax分别为最大力和最大压入深度；h为对应矿物的硬度。
ut＝ue up＝ue upp ufrac；
ufrac＝up-upp；
式中，ut表示压痕过程中的总能量，主要包括弹性能ue和塑性能up；fl和ful分别表示加载力和卸载力，力-位移曲线满足幂函数关系，k，m，n，α均为数据拟合参数；hh为饱载开始时的压入深度；hr为压痕残余深度；ufrac为断裂能；upp为纯塑性能；gc为能量释放率；kc为对应矿物的断裂韧性。
具体地，根据岩石宏观力学参数及工程参数评估模型对现场需求的工程参数进行评估包括：
结合测井资料，对脆性指标、断裂韧性和可压裂性指标进行评估。
测井资料包括：储层埋深及岩性。
bi＝σcσt/2；
ki＝0.313 0.027e′；
式中，bi为储层岩石的脆性指标；σc和σt分别为岩石单轴抗压强度和抗拉强度，由数值计算得出；ki为储层岩石的断裂韧性；e’为岩石宏观弹性模量，由数值计算得出；bi_n和ki_n分别为归一化脆性指数和归一化断裂韧性；bmax和bmin为目标储层的最大和最小脆性指标；kmax和kmin为目标储层的最大和最小断裂韧性。
为了更清楚介绍本发明实施例，下面从本发明实施例的使用方法上予以介绍。
将深井随钻排出的代表性岩屑烘干24h，然后将其分成若干份用于不同的测试项目，便于比较。
取出其中一份岩屑碾磨成200目细粉，置于x射线衍射仪中进行主要矿物成分定性和定量分析。
取出其中一份岩屑用于微米ct扫描，将ct扫描切片进行三维重构，形成三维数字岩芯，获取矿物颗粒的微观结构特征；进一步将岩屑制成厚度为30μm的薄片放置于偏光显微镜下观察，确定矿物颗粒的微观结构，与ct扫描结果互为补充。
将其中一份岩屑表面抛光后在压痕试验机上进行微观力学参数测试，根据矿物颗粒结构和尺寸选择合适的压头及加载条件。在压痕试验机的显微镜下区分各矿物的分布区域，对每种矿物进行单独测试，获取每种矿物成分的微观力学参数。微观力学参数主要包括弹性模量e、硬度h和断裂韧性kc，均可以通过荷载-位移曲线计算得出，公式如下：
式中，e为对应矿物的弹性模量；νs为待测样品的泊松比；ei和νi分别为压头的弹性模量和泊松比；er为压痕模量；a(hc)为压头与试样的接触投影面积，与压入深度hc有关；s为接触刚度；β和ε是与压头形状有关的常数；fmax和hmax分别为最大力和最大压入深度；h为对应矿物的硬度。
对于断裂韧性的求解，计算过程如下：
ut＝ue up＝ue upp ufrac
ufrac＝up-upp
式中，ut表示压痕过程中的总能量，主要包括弹性能ue和塑性能up；fl和ful分别表示加载力和卸载力，力-位移曲线满足幂函数关系，k，m，n，α均为数据拟合参数；hh为饱载开始时的压入深度；hr为压痕残余深度；ufrac为断裂能；upp为纯塑性能；gc为能量释放率；kc为对应矿物的断裂韧性。
基于上述岩屑的测试分析，获取岩石特征基因，包括矿物成分、矿物结构特征和矿物微观力学参数。采用数值计算方法，构建三维计算模型，将微观力学参数赋予模型中离散的材料点(各矿物成分)，选取合适尺寸的代表性体积单元(rve)，评估岩石的宏观力学性能，包括抗压/拉强度、弹性模量和泊松比。
采用工程参数评价模型对储层宏观力学参数进行有效评估，如下所示：
bi＝σcσt/2
ki＝0.313 0.027e′
式中，bi为储层岩石的脆性指标；σc和σt分别为岩石单轴抗压强度和抗拉强度，由数值计算得出；ki为储层岩石的断裂韧性；e’为岩石宏观弹性模量，由数值计算得出；bi_n和ki_n分别为归一化脆性指数和归一化断裂韧性；bmax和bmin为目标储层的最大和最小脆性指标；kmax和kmin为目标储层的最大和最小断裂韧性。
上述工程参数可直接应用于现场钻井和压裂工程设计与优化，更加准确和便捷。关于脆性、断裂韧性和可压裂性的评价方法有多种，本发明仅列举了现场应用广泛且准确度高的方法。基于岩石微观特征基因，也可以通过矿物成分含量、压痕载荷-位移关系以及数值计算中的应力-应变关系对工程参数进行评估，从而指导压裂等相关工程设计和实施。
本申请提供了一种深井岩屑微观特征基因测定以及岩石力学与储层工程参数评估方法，通过岩屑微观特征基因测定，可解决深部储层取芯费时费力、无法准确获取储层工程参数的难题，基于岩屑微观特征基因可以更准确地获取现场需求的工程参数，对钻井和压裂工程设计与优化具有重要意义。
最后所应说明的是，以上具体实施方式仅用以说明本发明的技术方案而非限制，尽管参照实例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
1.一种岩石力学与储层工程参数评估方法，其特征在于，包括：
获取深井随钻排出的岩屑；
对所述岩屑进行x射线衍射，获取所述岩屑的矿物成分及含量；
对所述岩屑进行ct扫描，并在显微镜下进行观察，获取所述岩屑的各矿物颗粒的微观结构特征；
对所述岩屑进行微/纳米压痕试验，测定所述岩屑的各矿物成分的微观力学参数；
由所述岩屑的矿物成分及含量、所述岩屑的各矿物颗粒的微观结构特征及所述岩屑的各矿物成分的微观力学参数获取岩屑微观特征基因；
根据所述岩屑微观特征基因构建三维计算模型，取代表性体积单元得到岩石宏观力学参数；
根据岩石宏观力学参数及工程参数评估模型对现场需求的工程参数进行评估。
2.根据权利要求1所述的岩石力学与储层工程参数评估方法，其特征在于，所述对岩屑进行x射线衍射包括：
将所述岩屑研磨成200目。
3.根据权利要求1所述的岩石力学与储层工程参数评估方法，其特征在于，所述对岩屑进行ct扫描，并在显微镜下进行观察包括：
将所述岩屑进行微米ct扫描，将ct扫描切片进行三维重构，形成三维数字岩芯，获取所述岩屑的矿物颗粒的微观结构特征；
将所述岩屑制成薄片，放置于偏光显微镜下观察，确定所述岩屑的矿物颗粒的微观结构特征，与ct扫描结果互为补充。
4.根据权利要求3所述的岩石力学与储层工程参数评估方法，其特征在于：
所述薄片的厚度为28-32μm。
5.根据权利要求1所述的岩石力学与储层工程参数评估方法，其特征在于，所述对岩屑进行微/纳米压痕试验包括：
在压痕试验机的显微镜下区分所述岩屑的各矿物的分布区域，对每种矿物进行单独测试，获取每种矿物成分的微观力学参数。
6.根据权利要求5所述的岩石力学与储层工程参数评估方法，其特征在于：
所述岩屑的矿物成分的微观力学参数包括：弹性模量、硬度及断裂韧性。
7.根据权利要求6所述的岩石力学与储层工程参数评估方法，其特征在于：
式中，e为对应矿物的弹性模量；νs为待测样品的泊松比；ei和νi分别为压头的弹性模量和泊松比；er为压痕模量；a(hc)为压头与试样的接触投影面积，与压入深度hc有关；s为接触刚度；β和ε是与压头形状有关的常数；fmax和hmax分别为最大力和最大压入深度；h为对应矿物的硬度；
ut＝ue up＝ue upp ufrac；
ufrac＝up-upp；
式中，ut表示压痕过程中的总能量，主要包括弹性能ue和塑性能up；fl和ful分别表示加载力和卸载力，力-位移曲线满足幂函数关系，k，m，n，α均为数据拟合参数；hh为饱载开始时的压入深度；hr为压痕残余深度；ufrac为断裂能；upp为纯塑性能；gc为能量释放率；kc为对应矿物的断裂韧性。
8.根据权利要求1所述的岩石力学与储层工程参数评估方法，其特征在于，所述根据岩石宏观力学参数及工程参数评估模型对现场需求的工程参数进行评估包括：
结合测井资料，对脆性指标、断裂韧性和可压裂性指标进行评估。
9.根据权利要求8所述的岩石力学与储层工程参数评估方法，其特征在于：
所述测井资料包括：储层埋深及岩性。
10.根据权利要求1所述的岩石力学与储层工程参数评估方法，其特征在于，所述根据岩石宏观力学参数及工程参数评估模型对现场需求的工程参数进行评估包括：
bi＝σcσt/2；
ki＝0.313 0.027e′；
式中，bi为储层岩石的脆性指标；σc和σt分别为岩石单轴抗压强度和抗拉强度，由数值计算得出；ki为储层岩石的断裂韧性；e’为岩石宏观弹性模量，由数值计算得出；bi_n和ki_n分别为归一化脆性指数和归一化断裂韧性；bmax和bmin为目标储层的最大和最小脆性指标；kmax和kmin为目标储层的最大和最小断裂韧性。
技术总结