CN110501758A - 一种砂砾岩储层纵向连续脆性指数预方法 - Google Patents

Abstract

本发明属于砂砾岩油气储层水力压裂技术领域，涉及一种砂砾岩储层纵向连续脆性指数预测方法。步骤如下：一、根据测井解释数据获取砂砾岩储层纵向连续的动态力学参数、孔隙度和砾石粒径及含量；二、建立弹性参数动静态转化关系；三、建立砂砾岩储层强度参数‑测井数据之间的解释模型；四、建立数值岩心模型，从而建立数值岩心宏细观参数转化关系；五、建立不同纵向深度的储层数值岩心并开展单轴压缩数值试验，获取应力、应变数据和全应力应变曲线；六、结合静态弹性模量计算脆性指数。本发明实现了对全井段砂砾岩储层岩心的数值化表征，可为砂砾岩储层可压性评价和射孔压裂层段优选提供指导，有助于促进油气藏经济高效开发。

Description

一种砂砾岩储层纵向连续脆性指数预方法

技术领域

本发明属于砂砾岩油气储层水力压裂技术领域，涉及一种砂砾岩储层纵向连续脆性指数预测方法。

背景技术

砂砾岩油气藏是近些年备受关注的一种重要的非常规油气资源。砂砾岩油气储层一般具有地质条件复杂、岩性变化大、物性复杂、非均质性强、渗透率低等特点，常规方法和技术手段难以实现有效改造和经济开采。水力压裂改造技术是砂砾岩储层增透开发的一项重要手段。水力压裂改造通过在井筒附近岩层中形成人工裂缝，并沟通天然弱面形成复杂裂缝，降低流体渗流阻力，改善渗流条件，从而实现增加油气开发产能。

实际压裂施工过程中，由于砂砾岩储层控制因素复杂，开采难度大，水力压裂改造后可能难以达到预期产能。导致压裂效果不理想的原因有很多，其中重要的一点在于对压裂射孔层段地质力学条件的把握不够准确，可压性评价不够准确。普遍认为储层可压性是表征储层压裂形成复杂裂缝并达到预期压裂效果的能力，受岩体脆性、断裂韧性、天然弱面和地应力条件等多种复杂因素影响；其中，脆性是影响可压性的最重要的指标。

砂砾岩一般由基质、随机分布的砾石和基质与砾石之间的胶结界面构成。胶结界面可视为一种天然弱面，砂砾岩力学性质主要受基质力学性质和胶结界面控制。界面总是围绕砾石存在，因此砾石的含量和粒径分别决定了界面的分布密度和面积，因此砂砾岩脆性主要受基质力学参数和砾石粒径及含量控制。脆性受砂砾岩储层脆性与压裂裂缝起裂和延伸以及复杂裂缝的形成等密切相关，脆性预测和评价直接影响到压裂改造过程中射孔层段选取和施工工艺参数选择等关键技术环节，对砂砾岩储层开发增产具有重要意义。

目前常用的脆性预测和评价方法多是针对页岩、煤岩和砂岩提出的，并没有针对砂砾岩的结构特点和力学特性提出的脆性评价方法。由于砾石颗粒的存在和天然裂隙不发育等结构特点，砂砾岩的力学特性与其他类型的岩石具有明显差异，常用的脆性评价方法对砂砾岩的适用性值得商榷。

目前常用的脆性评价方法中，基于矿物含量的脆性评价模型未考虑岩体内部胶结方式、成岩环境和地应力不同导致的脆性差异；基于室内岩心测试所得力学参数的脆性评价模型虽然能在一定程度上反映力学特性对脆性的影响，但只能得到取心层段的脆性特征，无法获得储层全井段纵向连续的脆性评价结果。基于测井曲线得到归一化弹性模量和泊松比定义的脆性评价模型，虽然能够得到纵向连续的脆性评价结果，但其是针对页岩储层提出的，对砂砾岩储层的适用性有待研究。因此目前常用的脆性评价方法对砂砾岩储层脆性预测均存在不适用性。需要针对砂砾岩储层结构特点和力学特征，结合岩石力学和地球物理方法，根据有限的物理岩心测试获取纵向上连续的脆性指数，提出适用于砂砾岩储层的脆性预测方法。

发明内容

本发明的目的在于提供一种砂砾岩储层纵向连续脆性指数预测方法，该方法结合测井解释数据、室内岩心物模测试结果和数值试验方法，对纵向上任意层位的砂砾岩岩心进行数值表征并开展数值试验，获取岩心残余强度特征预测砂砾岩储层纵向连续的脆性指数，为射孔层段优化及施工工艺参数优选提供参考，有助于提高储层纵向上的有效改造体积。

本发明采用的技术方案是：

一种砂砾岩储层纵向连续脆性指数预测方法，包括如下步骤：

第一步：根据测井解释数据获取砂砾岩储层纵向连续的动态岩石力学参数和孔隙度；

所述测井数据具体包括横波时差、纵波时差、自然伽马和密度测井曲线等，所述纵向连续的动态岩石力学参数具体为动态弹性模量和动态泊松比，以及由测井解释得到的单轴抗压强度和残余应力水平；

第二步：开展岩心单轴压缩测试，获取纵向非连续的静态力学参数，结合由测井数据得到的动态弹性参数和由室内测试得到静态力学参数中的静态弹性参数，采用最小二乘回归方法建立砂砾岩储层弹性参数的动静态转化关系，见公式(1)和(2)，由此获得纵向连续的静态弹性参数分布；

E_s＝f₁(E_d) (1)

ν_s＝f₂(ν_d) (2)

式中：E_s和E_d分别为静态和动态弹性模量，单位GPa；ν_s和ν_d别为静态和动态泊松比，无量纲；

所述岩心单轴压缩测试采用轴向位移加载，加载速率为2×10^-6m/s；所述静态力学参数具体为静态弹性模量、静态泊松比、单轴抗压强度和内摩擦角；弹性参数具体为弹性模量和泊松比；

第三步：利用砂砾岩岩心的静态强度参数对由测井解释得到的动态强度参数进行归位校核，建立砂砾岩储层强度参数-测井数据之间的解释模型，见公式 (3)和(4)；

σ_c＝F(C₁) (3)

θ＝F(C₂,Φ) (4)

式中：σ_c为校核后的单轴抗压强度；θ为校核后的残余强度系数，θ＝σ_r/σ_c，σ_r为残余强度；Φ为测井解释得到的岩心孔隙度；C₁和C₂分别为不同类型的测井数据；

利用砂砾岩储层弹性参数的动静态转化关系和强度参数-测井数据解释模型，对测井解释力学参数进行校核，得到砂砾岩储层纵向连续的精细的力学参数表征。

所述强度参数具体为单轴抗压强度和内摩擦角。

第四步：根据真实物理岩心的尺寸以及砾石含量和粒径，建立储层数值岩心有限元模型，对模型划分网格和施加轴向位移荷载并开展多次单轴压缩数值试验；利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定，保持砾石和界面材料的细观力学参数不变，从而建立数值岩心宏细观参数转化关系，见公式(6)和(7)；

式中：σ_c和σ_c0分别为岩心宏观和细观单轴抗压强度，MPa；E_s和E_s0分别为韦伯分布(Weibull分布)赋值时的岩心宏观和细观弹性模量均值，GPa；m为均值度系数，用以表征岩心力学参数的非均质性，m值越小则非均质性越强。

所建立的数值岩心一般为直径25mm、高50mm的圆柱体模型，其中，砾石建模时将其概化为球状或椭球状，按照测井数据显示的含量和粒径随机分布于数值岩心模型中；

所述利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定，保持砾石和界面材料的细观力学参数不变，具体做法为：在数值模型中给定基质细观单轴抗压强度均值和细观弹性模量均值，其他所需输入参数如泊松比和内摩擦角按相应的宏观参数输入，通过改变均值度系数分别拟合岩心宏观单轴抗压强度和宏观弹性模量与均值度系数的关系，进而根据物理岩心的真实宏观单轴抗压强度和宏观弹性模量反算与之对应的基质细观力学参数。最后将标定后的细观力学参数输入数值岩心模型进行单轴压缩试验，根据应力应变曲线形态对细观力学参数进行校验。要求物理岩心和数值岩心的应力应变曲线形态相同，且峰值强度和应变以及残余强度和应变相对大小均不超过5％，如超过5％则对细观力学参数进行多次标定和校验。

第五步：根据步骤第二步、第三步中得到的岩心静态宏观力学参数和第四步中得到的宏细观力学参数转化关系，可对任意层位储层对应的数值岩心细观力学参数进行标定，

将细观单轴抗压强度和E_s0细观弹性模量均值输入到真实破裂过程分析软件(RFPA3D)中，得到不同纵向深度的储层数值岩心模型，对储层数值岩心模型开展单轴压缩数值试验，获取应力、应变数据和全应力应变曲线。

所述通过对不同纵向深度对应的数值岩心开展细观力学参数标定和单轴压缩数值试验，实现了砂砾岩储层纵向上连续的岩心精细化数值表征，为储层纵向连续的脆性预测提供了技术支持。

第六步：根据不同纵向深度储层对应的数值岩心应力应变曲线，获取对应的脆性破坏残余强度特征参数σ_r和ε_r，进而根据式(8)计算脆性指数BI，由此得到纵向连续的脆性指数预测值。

式中，σ_r为残余强度，ε_r为残余应变；

所述结合岩心弹性模量和残余强度特征定义的脆性指数，不仅有效的综合考虑了峰前和峰后两个变形阶段的破坏特征，而且参数获取方便，计算简单，物理意义明确，有助于实际压裂现场应用和推广。

本发明具备的有益效果：

1.本发明结合岩石力学和地球物理方法，充分利用测井数据和室内物理岩心测试结果，借助砂砾岩岩心数值表征技术实现对任意纵向层位砂砾岩岩心进行数值建模，尤其是未能钻取岩心层段的储层力学参数表征，可有效克服常规脆性评价方法难以实现对全井段纵向连续脆性指数进行预测的不足。

2.针对砂砾岩的结构特点和力学特性，考虑了基质力学参数以及砾石粒径和含量对砂砾岩脆性的影响，且砂砾岩脆性评价方法。采用弹性模量和残余强度特征定义脆性指数，综合考虑了峰前和峰后两个变形阶段的破坏特征，脆性评价结果较其他方法更加准确可靠。

3.本发明脆性评价方法参数获取方便，仅需要常规测井解释数据和岩心室内测试结果，且计算简单，物理意义明确，有助于在砂砾岩油气储层压裂工程中应用和推广。

4.采用本发明提供的砂砾岩脆性评价方法，可为砂砾岩储层可压性评价、射孔压裂层段优选和施工工艺参数选取提供参考，提高压裂效果和储层改造体积。

附图说明

图1为本发明的砂砾岩岩心数值表征图。

图2为本发明的物理实验和数值模拟所得单轴压缩应力应变曲线图。

图3为本发明的数值岩心单轴压缩试验所得应力应变曲线示意图。

图4为本发明砂砾岩储层纵向连续的动态力学参数和孔隙度分布图。

具体实施方式

以下结合技术方案和附图详细叙述本发明的具体实施方式。

一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，包括如下步骤：

(a)根据测井解释数据获取砂砾岩储层纵向连续的动态力学参数和孔隙度；

(b)开展岩心单轴压缩测试，获取纵向上非连续的静态力学参数，结合由测井数据得到的动态弹性参数和由室内测试得到的静态弹性参数，采用最小二乘回归方法建立砂砾岩储层弹性参数的动静态转化关系，见公式(1)和(2)，由此获得纵向连续的静态弹性参数分布；

E_s＝f₁(E_d) (1)

ν_s＝f₂(ν_d) (2)

式中：E_s和E_d分别为静态和动态弹性模量，GPa，ν_s和ν_d别为静态和动态泊松比，无量纲；

(c)利用砂砾岩岩心静态强度参数对由测井解释得到的动态强度参数进行归位校核，建立砂砾岩储层强度参数-测井数据之间的解释模型，见公式(3)和(4)；

σ_c＝F(C₁) (3)

θ＝F(C₂,Φ) (4)

式中：σ_c为校核后的单轴抗压强度；θ为校核后的残余强度系数，θ＝σ_r/σ_c，σ_r为残余强度；Φ为测井解释得到的岩心孔隙度；C₁和C₂分别为不同类型的测井数据；

利用砂砾岩储层弹性参数的动静态转化关系和强度参数-测井数据解释模型，对测井解释力学参数进行校核，得到砂砾岩储层纵向连续的精细的力学参数表征。

(d)根据真实物理岩心的尺寸以及砾石含量和粒径，建立储层岩心数值有限元模型(如图1所示)，对模型划分网格和施加轴向位移荷载并开展多次单轴压缩数值试验；利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定(如图2所示)，保持砾石和界面材料的细观力学参数不变，从而建立数值岩心宏细观参数转化关系，见公式(6)和(7)；

式中：σ_c和σ_c0分别为岩心宏观和细观单轴抗压强度，MPa；E_s和E_s0分别为 Weibull分布赋值时的岩心宏观和细观弹性模量均值，GPa；m为均值度系数，用以表征岩心力学参数的非均质性，m值越小则非均质性越强；

(e)根据步骤b、c中得到的岩心静态宏观力学参数和步骤d中得到的宏细观力学参数转化关系，可对任意层位储层对应的数值岩心细观力学参数进行标定，从而建立不同纵向深度的储层数值岩心并开展单轴压缩数值试验，获取应力、应变数据和全应力应变曲线。

(f)根据不同纵向深度储层对应的数值岩心应力应变曲线，获取对应的脆性破坏残余强度特征参数σ_r和ε_r(如图3所示)，进而根据式(8)计算脆性指数BI，由此得到纵向连续的脆性指数预测值。

式中，σ_r为残余强度，ε_r为残余应变；

如图1所示，本发明采用有限元数值模型对任意纵向深度的物理岩心进行表征，建立的砂砾岩数值岩心一般为直径25mm、高50mm的圆柱体模型，其结构构成包括基质、砾石和基质与砾石之间的胶结界面。砾石建模时将其概化为球状或椭球状，按照测井数据显示的含量和粒径随机分布于数值岩心模型中；

如图2所示，本发明利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定，要求数值试验和物理实验得到应力应变曲线形态保持基本一致，重点关注峰值强度点和残余强度点是否基本重合，如不重合则对细观力学参数进行多次标定和校验。

如图3所示，本发明根据不同纵向深度储层对应的数值岩心应力应变曲线，获取对应的脆性破坏残余强度特征参数σ_r和ε_r，进而根据式(8)计算脆性指数。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

实例1

以胜利油田砂砾岩储层某区块-3382.4m至-3685.2m储层为例，根据钻取岩心室内物理实验结果和测井解释数据对其纵向连续脆性指数进行预测。

第一步，根据测井解释数据获取砂砾岩储层纵向连续的动态力学参数和孔隙度分布，如图4所示；

第二步，根据岩心室内单轴压缩测试结果，采用最小二乘回归方法建立砂砾岩储层弹性参数的动静态转化关系，如图4所示；

E_s＝3603.148-0.236E_d

ν_s＝-0.046+0.734ν_d

第三步，利用砂砾岩岩心静态强度参数对由测井解释得到的动态强度参数进行归位校核，建立如下式所示的砂砾岩储层强度参数-测井数据之间的解释模型，根据测井数据计算砂砾岩储层纵向连续的单轴抗压强度和内摩擦角，如图4 所示；

其中，V_p为纵波速度，单位km/s；Δt为声波时差，单位μs/m；ρ为密度，单位g/cm³。

第三步，以-3549.3m层段为例，根据测井解释数据获得砾石含量为38.5％，平均粒径为6.3mm，建立直径25mm、高50mm的圆柱体岩心数值有限元模型(如图1所示)，对模型划分网格和施加轴向位移荷载并开展多次单轴压缩数值试验；利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定(如图2所示)，保持砾石和界面材料的细观力学参数不变，从而建立数值岩心宏细观参数转化关系，见公式(6)和(7)；

第四步，取均值度m为3，根据物理岩心单轴抗压强度76.5MPa，弹性模量为29.8GPa，分别代入式公式(6)和(7)计算得到细观力学参数σ_c0为247.45MPa， E_s0为37.12Gpa，并输入数值岩心模型。开展单轴压缩数值试验得到数值岩心应力应变曲线并根据特征点应力应变数据对细观力学参数进行校核，要求物理岩心和数值岩心应力应变曲线形态一致且峰值强度、峰值应变和残余强度、残余应变相对大小均低于5％，若相对大小高于5％则调整细观力学进行多次校核，直至满足要求为止。最后残余强度特征σ_r为76.5Mpa，E_s为29.8Gpa，ε_r为0.064，代入式(8)计算得到该层段岩心脆性指数为1.278。

第五步，根据测井解释数据对任意层段岩心进行数值化表征，进而按照上述步骤计算脆性指数，获取砂砾岩储层纵向连续的脆性指数预测结果，如图4 所示。

Claims (8)

1.一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，包括如下步骤：

第一步：根据测井解释数据获取砂砾岩储层纵向连续的动态力学参数和孔隙度；

第二步：开展岩心单轴压缩测试，获取纵向上非连续的静态力学参数，结合由测井数据得到的动态弹性参数和由室内测试得到的静态弹性参数，采用最小二乘回归方法建立砂砾岩储层弹性参数的动静态转化关系，见公式(1)和(2)，由此获得纵向连续的静态弹性参数分布；

E_s＝f₁(E_d) (1)

ν_s＝f₂(ν_d) (2)

式中：E_s和E_d分别为静态和动态弹性模量，GPa，ν_s和ν_d别为静态和动态泊松比，无量纲；

第三步：利用砂砾岩岩心静态强度参数对由测井解释得到的动态强度参数进行归位校核，建立砂砾岩储层强度参数-测井数据之间的解释模型，见公式(3)和(4)；

σ_c＝F(C₁) (3)

θ＝F(C₂,Φ) (4)

式中：σ_c为校核后的单轴抗压强度；θ为校核后的残余强度系数，θ＝σ_r/σ_c，σ_r为残余强度；Φ为测井解释得到的岩心孔隙度；C₁和C₂分别为不同类型的测井数据；

利用砂砾岩储层弹性参数的动静态转化关系和强度参数-测井数据解释模型，对测井解释力学参数进行校核，得到砂砾岩储层纵向连续的精细的力学参数表征；

第四步：根据真实物理岩心的尺寸以及砾石含量和粒径，建立储层岩心数值有限元模型，对模型划分网格和施加轴向位移荷载并开展多次单轴压缩数值试验；利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定，保持砾石和界面材料的细观力学参数不变，从而建立数值岩心宏细观参数转化关系，见公式(6)和(7)；

式中：σ_c和σ_c0分别为岩心宏观和细观单轴抗压强度，MPa；E_s和E_s0分别为Weibull分布赋值时的岩心宏观和细观弹性模量均值，单位GPa；m为均值度系数，用以表征岩心力学参数的非均质性，m值越小则非均质性越强；

第五步：根据步骤第二步、第三步中得到的岩心静态宏观力学参数和第四步中得到的宏细观力学参数转化关系，可对任意层位储层对应的数值岩心细观力学参数进行标定，从而建立不同纵向深度的储层数值岩心并开展单轴压缩数值试验，获取应力、应变数据和全应力应变曲线；

第六步：根据不同纵向深度储层对应的数值岩心应力应变曲线，获取对应的脆性破坏残余强度特征参数σ_r和ε_r，进而根据式(8)计算脆性指数BI，由此得到纵向连续的脆性指数预测值；

式中，σ_r为残余强度，ε_r为残余应变。

2.根据权利要求1所述的一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，第一步中，所述测井数据具体包括横波时差、纵波时差、自然伽马和密度测井曲线，所述纵向连续的动态岩石力学参数具体为动态弹性模量和动态泊松比，以及由测井解释得到的单轴抗压强度和残余应力水平。

3.根据权利要求1或2所述的一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，第二步中，所述岩心单轴压缩测试采用轴向位移加载，加载速率为2×10^-6m/s；所述静态力学参数具体为静态弹性模量、静态泊松比、单轴抗压强度和内摩擦角；弹性参数具体为弹性模量和泊松比。

4.根据权利要求1或2所述的一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，第三步中，所述强度参数具体为单轴抗压强度和内摩擦角；第四步中，所建立的数值岩心一般为直径25mm、高50mm的圆柱体模型，砾石建模时将其概化为球状或椭球状，按照测井数据显示的含量和粒径随机分布于数值岩心模型中。

5.根据权利要求3所述的一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，第三步中，所述强度参数具体为单轴抗压强度和内摩擦角；第四步中，所建立的数值岩心一般为直径25mm、高50mm的圆柱体模型，砾石建模时将其概化为球状或椭球状，按照测井数据显示的含量和粒径随机分布于数值岩心模型中。

6.根据权利要求1或2或5所述的一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，第四步中，所述利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定，保持砾石和界面材料的细观力学参数不变，具体做法为：在数值模型中给定基质细观单轴抗压强度均值和细观弹性模量均值，其他所需输入参数包括泊松比和内摩擦角按相应的宏观参数输入，通过改变均值度系数分别拟合岩心宏观单轴抗压强度和宏观弹性模量与均值度系数的关系，进而根据物理岩心的真实宏观单轴抗压强度和宏观弹性模量反算与之对应的基质细观力学参数；最后将标定后的细观力学参数输入数值岩心模型进行单轴压缩试验，根据应力应变曲线形态对细观力学参数进行校验；要求物理岩心和数值岩心的应力应变曲线形态相同，且峰值强度和应变以及残余强度和应变相对大小均不超过5％，如超过5％则对细观力学参数进行多次标定和校验。

7.根据权利要求3所述的一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，第四步中，所述利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定，保持砾石和界面材料的细观力学参数不变，具体做法为：在数值模型中给定基质细观单轴抗压强度均值和细观弹性模量均值，其他所需输入参数包括泊松比和内摩擦角按相应的宏观参数输入，通过改变均值度系数分别拟合岩心宏观单轴抗压强度和宏观弹性模量与均值度系数的关系，进而根据物理岩心的真实宏观单轴抗压强度和宏观弹性模量反算与之对应的基质细观力学参数；最后将标定后的细观力学参数输入数值岩心模型进行单轴压缩试验，根据应力应变曲线形态对细观力学参数进行校验；要求物理岩心和数值岩心的应力应变曲线形态相同，且峰值强度和应变以及残余强度和应变相对大小均不超过5％，如超过5％则对细观力学参数进行多次标定和校验。

8.根据权利要求4所述的一种砂砾岩储层纵向连续脆性指数预测方法，其特征在于，第四步中，所述利用物理岩心的力学参数和应力应变曲线对相应数值岩心基质的细观力学参数进行标定，保持砾石和界面材料的细观力学参数不变，具体做法为：在数值模型中给定基质细观单轴抗压强度均值和细观弹性模量均值，其他所需输入参数包括泊松比和内摩擦角按相应的宏观参数输入，通过改变均值度系数分别拟合岩心宏观单轴抗压强度和宏观弹性模量与均值度系数的关系，进而根据物理岩心的真实宏观单轴抗压强度和宏观弹性模量反算与之对应的基质细观力学参数；最后将标定后的细观力学参数输入数值岩心模型进行单轴压缩试验，根据应力应变曲线形态对细观力学参数进行校验；要求物理岩心和数值岩心的应力应变曲线形态相同，且峰值强度和应变以及残余强度和应变相对大小均不超过5％，如超过5％则对细观力学参数进行多次标定和校验。