考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度预测模型.pdf

年 月第 卷第 期西 北 工 业 大 学 学 报 ：收稿日期：基金项目：安徽省高等学校自然科学研究重点项目（）、矿山建设工程安徽省高校重点实验室开放课题（）、安徽理工大学校青年基金（）、安徽理工大学高层次引进人才科研启动基金（）、国家自然科学基金面上项目（，）、中央高校基本科研业务费（）与安徽省大学生创新创业训练计划（）资助作者简介：石浩（），安徽理工大学讲师，主要从事岩石力学及锚杆支护研究。通信作者：张后全（），中国矿业大学副教授，主要从事岩石力学及锚杆支护研究。：考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度预测模型石浩，张后全，吴疆宇，宋雷，李明，荣传新，陆鹏举安徽理工大学 矿山建设工程安徽省高校重点实验，安徽 淮南；中国矿业大学 深部岩土力学与地下工程国家重点实验室，江苏 徐州；中国矿业大学 力学与土木工程学院，江苏 徐州；安徽理工大学 土木建筑学院，安徽 淮南；山东科技大学 省部共建矿山岩层智能控制与绿色开采国家重点实验室培育基地，山东 青岛；中煤矿山建设集团有限责任公司，安徽 合肥 摘 要：构建统筹考虑岩石内部宏细微观缺陷影响的强度预测模型对于保障智慧矿山的建设以及井下人员设备的安全都具有重要意义。综合考虑岩石微元强度的韦伯分布与裂纹长度的幂律排布规律，采用数值模拟与理论推导方法分别建立损伤岩石试件的离散元计算模型与强度预测模型，并利用数值计算结果对理论模型的合理性进行验证。结果表明：通过编程同时实现了 计算模型中微元强度的韦伯分布与裂纹长度的幂律排布，分析了岩石中的宏、细微观损伤与相应分布参数间的定量对应关系。数值建立了 个同时考虑宏、细微观损伤的砂岩试件并对其进行了模拟加载，实现了对微元强度、预制裂纹信息及岩石强度间关系的统计分析。根据模拟结果构建了砂岩试件单轴抗压强度四维空间散点，得到了多损伤参量影响下的试件抗压强度流动规律。联合 方法及岩石损伤概率分布理论，推导建立了 参数的岩石强度预测模型。该模型能够同时描述细微观微元强度与宏观裂纹信息对岩石强度的影响，并且理论模型计算结果与数值模拟结果高度吻合，相关系数达。关 键 词：微元损伤统计；韦伯分布；预制裂纹长度；幂律排布；岩石强度预测模型中图分类号：文献标志码：文章编号：（）随着煤炭开采工艺的进步、机械化水平的提高以及信息技术的发展，以数字化、信息化为前提和基础的智慧矿山建设已经取得了阶段性进展；相信在不远的未来，我国必然会实现深部煤炭的无人化开采。智慧矿山的建设对岩石力学特别是强度的准确判识要求极高，而岩石内部的各类损伤对其强度特性的影响占重要地位，是引起岩石强度离散的关键。因此，为了保障智慧矿山的建设以及井下人员、设备的安全，有必要掌握岩石内部的各类损伤对其强度特性的影响规律。目前，关于损伤缺陷影响岩石强度特性的研究已取得了十分丰硕的成果，并且根据研究对象的不同可分为 类。第一类研究主要针对岩石微元体强度劣化等微细观损伤缺陷。此类研究基于自然岩石为非均质材料，认为其宏观破坏对应微细观破坏由量变到质变的过程，。专家学者多从唯第 期石浩，等：考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度预测模型象学角度出发，假定微元体损伤参量服从韦伯统计分布规律，基于 准则或 准则建立与温度、蠕变等特定工况相关的岩石损伤本构模型。相关成果对于丰富岩石损伤描述方法及揭示岩石损伤内在机制具有重要意义，但众多研究主要是基于理论方法开展的，重在分析韦伯分布参数对应力应变关系（曲线）的影响，。在本质上所研究岩石试件仍以黑箱的形式存在，无法设定微元体强度或损伤参量的分布规律，导致理论分布假设的合理性得不到充分验证。第二类研究主要针对裂纹、裂隙等宏观损伤缺陷。此类研究预制缺陷数量差异较大，多数研究预制缺陷数量小于，也有少量研究预制缺陷数量成百上千。对于预制缺陷数量小于 的研究，研究者主要分析单个裂纹、裂隙的倾角、数量、位置、张开程度以及多个裂纹、裂隙的组合形式等对岩石强度的影响。相关成果对于了解缺陷岩石的力学性能具有重要意义，但对缺陷信息（如倾角等）的设置过于理想，也导致研究结果的普适性较弱。对于预制缺陷数量成百上千的研究，研究者通常注重对缺陷生成方法（包括概率分布或利用地层扫描和计算机信息提取相结合）的提出和改善，而对缺陷如何影响岩石强度特性的分析相对不足。另外，裂纹长度服从幂律排布的规律也未被充分考虑。赵洪宝等回顾了关于裂隙（纹）岩石的研究，指出应当考虑分布裂纹数量及裂纹间的相互作用并建立与裂纹分布相关的强度理论，才能真正提高研究成果的工程应用价值。以赵洪宝等人的观点作为指引，本课题组根据 理论进行推导，首先建立了与预制裂纹分布信息相关的损伤砂岩强度模型。鉴于自然界的岩石试件其内部宏细微观缺陷并存，仅对宏观或细微观的缺陷进行单独研究难以全面地掌握岩石承载失效特性，而目前统筹考虑岩石宏细微观缺陷对其强度特性影响的研究较少。因此，本文综合考虑岩石微元强度的韦伯分布与裂纹长度的幂律排布规律，采用数值模拟与理论推导方法分别构建损伤岩石试件的离散元计算模型与强度预测模型，并利用数值计算结果对理论模型的合理性进行验证。岩石微元强度韦伯分布与裂纹长度幂律排布理论 岩石微元强度韦伯分布理论岩石可视为由众多基元介质组成的有机整体，岩石破坏过程对应基元体损伤由量变到质变，基元介质的力学性质描述适用于统计的方法。和众多专家学者的研究一致，本文假设岩石微元破坏的概率随微元强度 （）的分布密度 为 （）（）式中，为试件的许用强度。根据（）式可得损伤度 与微元破坏的概率密度的关系为（）（）（）根据（）式可知，确定岩石的微元强度与破坏概率是建立岩石损伤模型抑或岩石强度预测模型的前提。目前，韦伯函数已被证明适合用于描述岩石微元强度或损伤的统计分布规律，对应的岩石微元破坏概率密度方程为 （）式中：，均为韦伯分布参数；为均质系数，用以反映微元体强度等力学参数分布的密集程度；代表尺度参数。岩石裂纹长度幂律排布理论已有研究表明，岩石中的裂纹长度服从幂律排布，对应的分布方程为 （）式中：表示裂纹长度；表示有关裂纹长度分布的自变量；另外，为大于 的常数且 （），而 为大于 的幂律指数。在实际应用中，需要预先给出预制裂纹的最大长度、最小长度 和幂律指数，从而确定特定的幂律排布形式。然后，根据（）式的增减性可以建立（）式所示关系 （）式中，分别为 的上界和下界。另外，考虑到（）式为超静定方程，及 均为未知参数，需要提前给一个参数赋值，故假定（）式成立西 北 工 业 大 学 学 报第 卷 或 或 （）根据（）（）式可以确定裂纹长度的幂律排布形式。在此基础上，通过给定预制裂纹数量，具体的裂纹长度可以利用线性插值（理论）或概率分配（数值模拟）的方式获得。考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度模拟研究无论是岩石微元强度的韦伯分布抑或裂纹长度的幂律排布，对其研究都难以通过室内实验开展，但理论模型的建立需有数据支撑。考虑到韦伯分布及幂律排布在本质上都是数学运算，完全可以通过数值方法编程实现，故在此采用数值模拟方法分析微元强度韦伯分布与裂纹长度幂律排布对岩石强度的影响。颗粒黏结模型（）作为离散元颗粒流软件，被广泛应用于岩土力学方面的研究。该软件优势契合本文研究主题，能够较好地模拟微元强度及裂纹长度等分布对介质模型强度的影响，故被本文研究选用。软件对岩土介质的表征主要采用颗粒及黏结两部分，其中颗粒为刚体，通常用圆盘表示；而对于颗粒间的黏结，适用于岩石材料的模型主要有 种，分别为接触黏结模型和平行黏结模型，两模型对应的颗粒及黏结的力学行为如图 所示。接触黏结模型（见图）相对简单，仅能表征由颗粒接触、碰撞产生的法向力和切向力；而平行黏结模型（见图）既能反映颗粒间的接触作用，亦能反映颗粒间的弯矩作用。通常认为，此 种模型在岩石内部均存在，故本文研究采用包含这 种模型的颗粒黏结模型（）。图 黏结模型及其微观力学行为示意图，无初始损伤砂岩数值模型及细观参数的确定研究岩石微元强度韦伯分布及裂纹长度幂律排布对其强度的影响，理想状态下需获取无初始损伤试件的强度特性，。但在现实条件下，所有可获得的岩石试件都不可避免地存在各类损伤且难以量化，专家学者一般选择无明显损伤的试件作为无初始损伤试件开展研究。按照此种思路，笔者基于标准尺寸青砂岩的单轴压缩室内实验，在先前的研究中已经建立了相应的 数值模型并标定了模型细观参数。模型参数已在表 中列出，本文研究均基于该数值模型及标定的细观参数。需要特别指出的是，所提青砂岩试件其单轴压缩实验强度约为 。表 介质细观参数参数数值颗粒最小尺寸 颗粒最大尺寸 密度（）孔隙率接触黏结模量 接触黏结刚度比摩擦因数平行黏结拉伸强度 平行黏结黏聚力 平行黏结摩擦角（）平行黏结模量 平行黏结刚度比 微元强度韦伯分布与裂纹长度幂律排布的实现及数值计算结果 微元强度韦伯分布的实现对于表 中列出的众多模型细观参数，用以控制微元强度的主要为黏结黏聚力及黏结拉伸强度。所以对于数值模型微元强度韦伯分布的实现，仅需利用 软件自带的 语言编程生成特定形式的韦伯分布随机数，随机数数量与黏结数量相同，然后将韦伯分布随机数分别与黏结的黏聚力和拉伸强度同时结合（取乘积），生成新的平行黏结黏聚力及平行黏结拉伸强度数值。本文综合考虑均质系数 对 分布覆盖区间范围及区间内分布梯度的影响，设置其值分别为，。通第 期石浩，等：考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度预测模型过对模型颗粒间的黏结信息进行统计，得到的颗粒黏结强度 模拟与理论分布之间的对比情况如图 所示。可以看出，理论分布与模拟分布高度一致，颗粒黏结强度均呈减速降低趋势。并且，均质系数越小，颗粒黏结强度参数极大值越大；当均质系数为，时，对应的黏结强度参数极大值分别为，。另外，各均质系数下的曲线交汇点对应的黏结强度值为，与表 中黏结强度值一致。图 颗粒黏结强度模拟分布与理论分布对比 裂纹长度幂律排布的实现在 软件中，裂纹以构件的形式存在。长度服从幂律排布的裂纹是利用 语言编译软件中的（）模块生成的，具体步骤为：设置裂纹的长度范围、长度分布形式、倾角范围、位置分布形式；给定裂纹的数量以及在试件中的分布范围，由此生成裂纹；赋值裂纹力学属性参数。本文设置的裂纹参数与文献一致，裂纹长度范围为，服从幂律排布；倾角范围为；位置分布随机，分布范围与试件面积范围一致。预制裂纹力学参数赋值基于平滑节理模型（），模型参数已在表 中列出，。可以看出，裂纹存在的主要作用是破坏裂纹面两侧黏结的黏聚力和拉伸强度。表 平滑节理模型的力学参数，参数数值单位面积法向刚度 单位面积切向刚度 摩擦因数拉伸强度 黏结强度 温度应力、风化、开采扰动、海水侵蚀以及地质作用程度的不同都会对岩石损伤产生较大影响，导致其内部裂纹信息存在较大差异。基于此种考量，本文建立的数值模型根据预制裂纹数量及裂纹分布幂律指数可以分为 类。具体而言，预制裂纹数量分别设置为，而幂律指数 分别设置为，该参数同样是在考虑幂律排布的覆盖区间范围及区间内分布梯度变化的情况下确定的。当预制裂纹数量为，预制裂纹长度 在不同幂律指数下的理论分布与模拟分布对比如图 所示。可以看出，模拟与理论分布形式基本一致，均呈减速降低趋势。并且，幂律指数越小，裂纹长度降低速率越平缓（见图 中幂律指数为 对应的分布曲线），对应长裂纹占比越高；幂律指数越大，裂纹长度降低更接近于两段式，先急速降低再平稳变化（见图 中幂律指数为 对应的分布曲线），对应长裂纹占比越低。图 理论与实际模拟得到的预制裂纹长度在不同幂律指数下的分布 数值计算结果根据 节及 节可知，本文共计建立了 均质系数 裂纹数量 幂律指数 个参数组合。为了尽可能降低试件强度离散性的影响，每个参数组合包含 个试件。需要指出的是，每个参数组合 个试件的所有参数信息都是一致的，但 个试件对应的微元强度韦伯分布模块及 模块调用次数分别为。当微元强度韦伯分布模块及 模块调用次数一定时，试件内部的微元强度和裂纹分布信息是完全一致的，即试件的模拟计算是可重复的；而当调用次数变化时，微元强度和单个裂纹的实际存在位置及形式将再次发生随机分布。西 北 工 业 大 学 学 报第 卷由上述介绍可知，本文针对损伤岩石共计建立并计算了 个数值模型试件，基本能够实现对微元强度、预 制 裂 纹 信 息 及 岩 石 强 度 间 关 系 的 统 计研究。以预制裂纹数量 为 轴（第一维），均值系数 为 轴（第二维），幂律指数 为 轴（第三维），单轴抗压强度 为颜色轴（第四维），根据模拟结果构建得到的试件单轴强度四维散点（，），如图所示。可以看出，各参数组合试件 的强度分布规律基本一致，抗压强度最大值超过 ，最小值接近 ；另外，抗压强度具有明晰的定向流动规律，其最小值均对应预制裂纹数量较多、均质系数较小且幂律指数较小处，抗压强度最大值均对应预制裂纹数量较少、均质系数较大且幂律指数较大处。图 个试件抗压强度模拟结果四维空间散点 考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度预测模型 考虑微元强度韦伯分布的砂岩强度预测模型及验证（）（）式分别给出了损伤变量与微元破坏概率密度的关系以及岩石微元强度与破坏概率的影响因素。考虑到岩石最终破坏是损伤不断增加的连续过程所引发的结果，故将（）式代入（）式进行积分，得到了损伤度与微元强度概率间的对应关系为 （）试件的有效应力 与损伤程度的关系为（）（）式中，为无初始损伤试件的最大承载应力，数值上取为无初始损伤试件的单轴压缩强度，根据节知该值为 。根据 的定义，其数学期望值（）、标准差 均可表示为，的函数，为 （）（）（）（）（）（）式中，符号 为伽马函数，其在实数域上的定义为第 期石浩，等：考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度预测模型（）（）（）越大，越接近材料均值。将（）、（）式代入（）式，则有效应力的期望值 可表示为 （）考虑试件强度的离散性及（）式的增减特性，将实际应力矩阵以区间范围的形式表示为 （）实际应用中需要给出应力矩阵与均质系数之间的确切关系，如此便需确定方差系数。考虑到试件（非）均质离散性随均质系数的变化而变化，假定方差系数服从表达形式覆盖广泛的 函数，由此可得实际应力矩阵与均质系数之间的关系 （）式中，和 均为方差系数表达式的待定常数。利用（）式拟合模拟结果，如图 所示，可以看出将方差系数假定为均质系数的 函数能够较好地描述有效应力随均质系数的变化规律，相关系数 超，说明该假定是合理的。并且计算可得，在 的变化区间内，和 的组合数值介于 之间，满足要求。图 理论与模拟得到的岩石抗压强度随均质系数的变化规律 微元强度韦伯分布与裂纹长度幂律排布耦合影响下的砂岩强度预测模型及验证当仅考虑裂纹对试件强度的影响时，笔者已根据 方法建立了有效应力与裂纹密度之间的关系 （）式中：表示单位面积的裂纹数量；为预制裂纹数量；为二维标准试件的面积；和 分别为裂纹密度系数和裂纹平均半长，均与幂律系数密切相关。同样利用 函数对文献 的数据进行回归，可得 和 与幂律指数 的关系分别为 ，相关系数 均超。根据（）式中单位面积裂纹数量 与预制裂纹数量 的关系以及裂纹平均半长 与幂律指数 之间的关系，可得有效应力 与预制裂纹数量 及幂律指数 之间的关系为 （）西 北 工 业 大 学 学 报第 卷 裂纹作为结构单元，其存在直接导致所在位置黏结的拉伸强度及黏结强度降低为（见表）。虽然裂纹对试件强度的作用在本质上也是通过影响（“杀死”）微元强度实现的，但其对黏结的破坏程度决定了受裂纹影响的黏结将可能不再受微元强度韦伯分布的影响，即裂纹与微元强度韦伯分布对黏结的影响互不干扰。据此，结合（）及（）式，叠加考虑微元强度韦伯分布及裂纹长度幂律排布对试件强度的影响，可建立 参数的砂岩强度预测模型为 （）结合图 的拟合结果以及文献中的数据可得模型涉及的 个参数数值，在表 中列出。表 砂岩强度预测模型参数数值参数数值参数数值理论模型建立后，需要对其合理性进行验证。本文首先利用（）式计算砂岩试件在不同预制裂纹数量、裂纹分布幂律指数及微元强度均值系数影响下的单轴抗压强度值，然后构建得到了砂岩单轴抗压强度理论值的四维空间散点，其与模拟结果平均值的对比如图 所示。可以直观看出，抗压强度理论分布（见图）与模拟分布（见图）基本一致，即抗压强度最小值均对应预制裂纹数量较多、均质系数较小且幂律指数较小处；抗压强度最大值均对应预制裂纹数量较少、均质系数较大且幂律指数较大处。上述是对理论结果与数值结果分布相似性的主观描述，为了定量描述理论结果与模拟结果的吻合程度，依照数值模拟设置的参数组合，利用（）式计算获得了砂岩试件在相应预制裂纹数量、裂纹分布幂律指数及微元强度均值系数影响下的单轴抗压强度，其与模拟结果平均值对比如图 所示。图中，图 抗压强度模拟结果平均值与理论计算值四维空间散点横坐标表示 个参数组合按照所得砂岩强度模拟均值自小至大依次排列时对应的序号。可以看出，理论结果与模拟结果高度一致，并且计算得相关系数，说明利用所建强度预测模型能够较好地描述微元强度韦伯分布与裂纹长度幂律排布对试件强度的影响。第 期石浩，等：考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度预测模型图 理论计算值与模拟数据的对比 结 论本文综合考虑岩石微元强度的韦伯分布与裂纹长度的幂律排布规律，构建了损伤砂岩试件的计算模型及强度预测理论模型，并利用数值计算结果对理论模型的合理性进行了验证。主要结论为：）通过编程同时实现了 计算模型中微元强度的韦伯分布与裂纹长度的幂律排布，分析了岩石中的宏、细微观损伤与相应分布参数间的定量对应关系。）数值建立了 个同时考虑宏、细微观损伤的砂岩试件并对其进行了模拟加载，实现了对微元强度、预制裂纹信息及岩石强度间关系的统计分析。根据模拟结果构建了砂岩试件单轴抗压强度四维空间散点（，），得到了多损伤参量影响下的试件抗压强度流动规律。）联合 方法及岩石损伤概率分布理论，推导建立了 参数的岩石强度预测模型，经与数值模拟结果对比分析，该模型被证明能够有效地描述细微观微元强度及宏观裂纹信息对岩石强度的影响。参考文献：丁恩杰，俞啸，夏冰，等 矿山信息化发展及以数字孪生为核心的智慧矿山关键技术 煤炭学报，（）：，（）：（）王国法，赵国瑞，任怀伟 智慧煤矿与智能化开采关键核心技术分析 煤炭学报，（）：，（）：（）付蒙，李江红，吴亚锋，等 钻柱黏滑振动特性仿真与产生机理分析 西北工业大学学报，（）：，（）：（）尤明庆 岩石的损伤、黏结和摩擦特性研究 岩土工程学报，（）：，（）：（）秦庆词，李克钢，杨宝威，等 岩石全应力应变过程关键特征点损伤特征分析 岩土力学，（增刊）：，（）：（）尤明庆 围压对岩石试样强度的影响及离散性 岩石力学与工程学报，（）：，（）：（）曹文贵，赵明华，刘成学 基于韦伯分布的岩石损伤软化模型及其修正方法研究 岩石力学与工程学报，（）：，（）：（）蒋浩鹏，姜谙男，杨秀荣 基于韦伯分布的高温岩石统计损伤本构模型及其验证 岩土力学，（）：西 北 工 业 大 学 学 报第 卷，（）：（）张慧梅，雷利娜，杨更社 基于韦伯统计分布的岩石损伤模型 湖南科技大学学报，（）：，（）：（）张庆贺，袁亮，方致远，等 基于多元监测信息融合的含裂隙类岩石材料破坏规律和预警方法 采矿与安全工程学报，（）：，（）：（）李天斌，高美奔，陈国庆，等 基于热力损伤本构参数的硬岩脆性评价方法 岩石力学与工程学报，（增刊）：，（）：（）刘建，赵国彦，彭府华 岩石介质弹塑性应变软化本构细观力学参数统计分布模型 煤炭学报，（增刊）：，（）：（），：，：，（）：宋义敏，李肖飞，许海亮，等 基于裂纹界面摩擦滑动演化的岩石断裂研究 中国矿业大学学报，（）：，（）：（）汪杰，付建新，宋卫东，等 岩石充填体组合模型力学特性及微裂纹演化特征颗粒流模拟 中国矿业大学学报，（）：，（）：（），：，（）：，（）：，：，：，（）：赵洪宝，胡桂林，李伟，等 预制裂隙岩石裂纹扩展规律的研究进展与思考 地下空间与工程学报，（增刊）：，（）：（），第 期石浩，等：考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度预测模型 ，（）：杨忠民，高永涛，吴顺川，等 基于等效岩体技术的节理参数对岩体强度影响性研究 中国矿业大学学报，（）：，（）：（），（）：，：石浩，宋雷，王敦显，等 基于裂纹长度幂律排布的损伤砂岩强度模型中国矿业大学学报，（）：，（）：（），（）：石浩，张后全，宋雷，等 动载下含预制裂纹砂岩的力学特性及破裂过程研究 振动与冲击，（）：，（）：（）石浩 采动岩体强度演化及锚岩相互作用研究 徐州：中国矿业大学，：，（），：张慧梅，李焕容，刘小宁，等 基于 理论的岩石损伤本构模型及参数探讨 武汉理工大学学报，（）：，（）：（）董陇军，李夕兵，马举，等 未知波速系统中声发射与微震震源三维解析综合定位方法及工程应用 岩石力学与工程学报，（）：，（）：（）徐飞，徐卫亚，温森，等 基于 的围岩稳定性评价 岩土力学，（）：，（）：（）西 北 工 业 大 学 学 报第 卷 ，；，；，；，；，；，：，：），），），：；引用格式：石浩，张后全，吴疆宇，等 考虑微元强度韦伯分布与裂纹长度幂律排布的砂岩强度预测模型 西北工业大学学报，（）：，（）：（）（：），