一种水作用下碳酸盐岩力学性能时效劣化数值模拟方法技术

本发明专利技术公开了一种水作用下碳酸盐岩力学性能时效劣化数值模拟方法，包括以下步骤：步骤1、建立碳酸盐岩在水作用下的力学性能劣化机制概化模型；步骤2、根据步骤1得到的机制概化模型建立不同浸水条件下碳酸盐岩力学性能时效劣化离散元模型，得到不同浸水条件下多边形块体间的块体间粘结的力学性能随时间的量化相关性；步骤3、标定离散元模型的力学参数；进行室内不同浸水条件和浸水时间下碳酸盐岩力学试验，对步骤2中离散元模型的力学参数进行标定；步骤4、运用已标定参数的离散元模型，开展不同浸水条件下碳酸盐岩力学性能离散元模拟。本发明专利技术真实反映了碳酸盐岩‑水相互作用机理，提出了其粘结二阶段劣化模型及其离散元模拟方法。

【技术实现步骤摘要】
一种水作用下碳酸盐岩力学性能时效劣化数值模拟方法
本专利技术涉及水作用下碳酸盐岩力学性能衰减数值模拟领域，具体是一种水作用下碳酸盐岩力学性能时效劣化数值模拟方法。
技术介绍
碳酸盐岩在我国分布十分普遍，典型岩石类型包括石灰岩、大理岩、白云岩、泥灰岩等。碳酸盐岩主要化学成分为方解石(主要成份碳酸钙)。水对碳酸盐岩的作用不仅仅限于软化作用，同时与碳酸钙在水长期作用下发生溶解密切相关。我国西南许多大型水电站库区的碳酸盐岩边坡(如三峡库区存在大量的石灰岩边坡，小湾水电站坝址区主要为大理岩)，另外我国石灰岩地层含有丰富的有机物残骸，使得该地层内存在大量的煤炭(奥陶系石灰岩地层及石炭二叠系)、石油及天然气资源。碳酸盐岩在地表水(水库、河流)及地下水等短期及长期浸泡作用下，具有不同的损伤破坏机制，其强度劣化具有明显的时效性。这种强度时效劣化效应可能造成岩体工程灾变，造成滑坡及地下工程围岩破坏。三峡水库自2003年开始蓄水起，造成了大量的石灰岩、泥灰岩边坡的破坏，如2008年发生的巫峡龚家坊滑坡。奥陶纪及二叠纪碳酸盐岩地层采煤突水灾害也经常发生，其浸水损伤破坏造成突水灾害是煤矿灾害第二灾害事件(仅次于瓦斯)，如开滦范各庄煤矿(石灰岩)2171工作面发生特大陷落柱突水，造成全矿被淹。碳酸盐岩长期与水接触，导致岩石内裂隙、孔洞中的矿物或者胶结物(如碳酸盐、铁的氧化物等)在水中发生溶解、侵蚀，长期蠕变效应使岩石损伤、缺陷累进性扩展，产生突发破坏，形成滑坡、地下工程围岩破坏及突水灾害。基于地表水、地下水作用下岩石强度弱化常直接采用强度折减法对工程岩体的破坏过程进行数值模拟。例如周桂云等运用有限元强度折减法，计算十里铺水电站边坡水位下降时边坡的稳定安全系数。但是强度折减方法无法体现碳酸盐岩-水相互作用机理，而且折减多少以及如何折减没有明确的理论依据。
技术实现思路
针对现有技术的不足，本专利技术拟解决的技术问题是，提供一种水作用下碳酸盐岩力学性能时效劣化数值模拟方法。本专利技术解决所述技术问题的技术方案是，提供一种水作用下碳酸盐岩力学性能时效劣化数值模拟方法，其特征在于该方法包括以下步骤：步骤1、根据水对碳酸盐岩中粘结的作用机理，建立碳酸盐岩在水作用下的力学性能劣化机制概化模型，该机制概化模型由矿物颗粒、成岩胶结物、沉积胶结物和水组成；运用离散单元法将碳酸盐岩内的矿物颗粒用多边形块体表征，胶结物用沉积类粘结和成岩类粘结共同表征；步骤2、根据步骤1得到的机制概化模型建立不同浸水条件下碳酸盐岩力学性能时效劣化离散元模型，得到不同浸水条件下多边形块体间的块体间粘结的力学性能随时间的量化相关性；步骤3、标定离散元模型的力学参数；进行室内不同浸水条件和浸水时间下碳酸盐岩力学试验，对步骤2中离散元模型的力学参数进行标定；步骤4、运用已标定参数的离散元模型，开展不同浸水条件下碳酸盐岩力学性能离散元模拟。与现有技术相比，本专利技术有益效果在于：(1)本专利技术的突出效益是真实反映了碳酸盐岩-水相互作用机理，提出了其粘结二阶段劣化模型(沉积胶结物瞬时破坏，而成岩胶结物则随浸水时间及酸性浓度渐进性时效劣化)及其离散元模拟方法。(2)数值方法成本低且具有较好的拓展性，本专利技术通过离散元方法实现了水作用下碳酸盐岩力学性能时效劣化的数值模拟，这样可以通过少量的室内试验标定力学性能参数后，可开展大量的不同压力及浸水时间的碳酸盐岩力学行为的数值模拟研究，而数值模拟不产生除计算机运行及人工费之外其它费用。(3)可应用工程岩体灾变模拟，具有很好的工程应用前景。本专利技术所建立的碳酸盐岩力学性能时效劣化数值模拟方法，可以有效地模拟浸水时间、水酸性环境等对碳酸盐岩力学性能的劣化规律，可以应用于浸水边坡及地下工程等灾变预测及演化过程的模拟，有利于工程岩体稳定性评价及灾害机制分析。附图说明图1为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法一种实施例的方法流程框图。图2为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法一种实施例的机制概化模型的初始状态图。图3为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法一种实施例的机制概化模型的短期浸水图。图4为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法一种实施例的机制概化模型的长期浸水图。图5为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法一种实施例的离散元模型的初始状态图。图6为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法一种实施例的离散元模型的短期浸水图。图7为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法一种实施例的离散元模型的长期浸水图。图8为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法一种实施例的室内不同浸水条件下碳酸盐岩单轴压缩试验图。图9为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法一种实施例的室内不同浸水条件下碳酸盐岩劈裂试验图。图10为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法实施例1的碳酸盐岩单轴压缩强度随浸泡时间变化的试验与模拟结果对比图。图11为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法实施例1的碳酸盐岩抗拉强度随浸泡时间变化的试验与模拟结果对比图。图12为本专利技术水作用下碳酸盐岩力学性能时效劣化数值模拟方法实施例1的碳酸盐岩弹性模量随浸泡时间变化的试验与模拟结果对比图。(图中：1、成岩胶结物；2、矿物颗粒；3、沉积胶结物；4、水；5、多边形块体；6、成岩类粘结；7、沉积类粘结；8、块体间粘结)具体实施方式下面给出本专利技术的具体实施例。具体实施例仅用于进一步详细说明本专利技术，不限制本申请权利要求的保护范围。本专利技术提供了一种水作用下碳酸盐岩力学性能时效劣化数值模拟方法(简称方法)，其特征在于该方法包括以下步骤：步骤1、根据水对碳酸盐岩中粘结的作用机理，建立碳酸盐岩在水作用下的力学性能劣化机制概化模型，该机制概化模型由矿物颗粒2、两种胶结物(成岩胶结物1和沉积胶结物3)和水4组成(参见图2-4)；运用离散单元法将碳酸盐岩内的矿物颗粒2用多边形块体5(block)表征，胶结物用沉积类粘结7和成岩类粘结6共同表征；水对沉积胶结物3的作用机制是快速破坏，对成岩胶结物1则为渐进性时效劣化的作用机制；步骤2、根据步骤1得到的机制概化模型建立不同浸水条件下碳酸盐岩力学性能时效劣化离散元模型，得到不同浸水条件下多边形块体5间的块体间粘结8的力学性能(强度和刚度)随时间的量化相关性；步骤3、标定离散元模型的力学参数；进行室内不同浸水条件和浸水时间下碳酸盐岩力学试验，对步骤2中离散元模型的力学参数进行标定；步骤4、运用已标定参数的离散元模型，开展不同浸水条件下碳酸盐岩力学性能离散元模拟。进一步，所述步骤1中，该机制概化模型包括：1)在短期浸水作用下，浸水仅破坏碳酸盐岩中沉积胶结物3；2)在长期浸水作用下，将逐渐破坏碳酸盐岩中成岩胶结物1且这类破坏具有显著的时效性。在机制概化模型中，多边形块体5间通过块体间粘结8(bond)产生相互作用，此块体间粘结8用沉积类粘结7(bond1)和成岩类粘结6(bond2)共同表征。在机制概化模型中，沉积类粘结7具有碳酸盐岩中沉积胶结物3的水理性质，成岩类粘结6具有碳酸盐岩中成岩胶结物1的水理性质。进一步，所述步骤2中在干燥状态下碳酸盐岩力学性能时效劣化离散本文档来自技高网...

【技术保护点】
1.一种水作用下碳酸盐岩力学性能时效劣化数值模拟方法，其特征在于该方法包括以下步骤：步骤1、根据水对碳酸盐岩中粘结的作用机理，建立碳酸盐岩在水作用下的力学性能劣化机制概化模型，该机制概化模型由矿物颗粒、成岩胶结物、沉积胶结物和水组成；运用离散单元法将碳酸盐岩内的矿物颗粒用多边形块体表征，胶结物用沉积类粘结和成岩类粘结共同表征；步骤2、根据步骤1得到的机制概化模型建立不同浸水条件下碳酸盐岩力学性能时效劣化离散元模型，得到不同浸水条件下多边形块体间的块体间粘结的力学性能随时间的量化相关性；步骤3、标定离散元模型的力学参数；进行室内不同浸水条件和浸水时间下碳酸盐岩力学试验，对步骤2中离散元模型的力学参数进行标定；步骤4、运用已标定参数的离散元模型，开展不同浸水条件下碳酸盐岩力学性能离散元模拟。

【技术特征摘要】
1.一种水作用下碳酸盐岩力学性能时效劣化数值模拟方法，其特征在于该方法包括以下步骤：步骤1、根据水对碳酸盐岩中粘结的作用机理，建立碳酸盐岩在水作用下的力学性能劣化机制概化模型，该机制概化模型由矿物颗粒、成岩胶结物、沉积胶结物和水组成；运用离散单元法将碳酸盐岩内的矿物颗粒用多边形块体表征，胶结物用沉积类粘结和成岩类粘结共同表征；步骤2、根据步骤1得到的机制概化模型建立不同浸水条件下碳酸盐岩力学性能时效劣化离散元模型，得到不同浸水条件下多边形块体间的块体间粘结的力学性能随时间的量化相关性；步骤3、标定离散元模型的力学参数；进行室内不同浸水条件和浸水时间下碳酸盐岩力学试验，对步骤2中离散元模型的力学参数进行标定；步骤4、运用已标定参数的离散元模型，开展不同浸水条件下碳酸盐岩力学性能离散元模拟。2.根据权利要求1所述的水作用下碳酸盐岩力学性能时效劣化数值模拟方法，其特征在于所述步骤2中在干燥状态下碳酸盐岩力学性能时效劣化离散元模型中，多边形块体间的块体间粘结的强度和刚度分别为：式(1)中，下标dry表示干燥状态；jcoh为与单轴抗压强度相关的微观参数，jten为与抗拉强度相关的微观参数，jkn为与弹性模量相关的微观参数；ηjcoh、ηjten、ηjkn分别对应微观上单轴抗压强度、抗拉强度、弹性模量的软化系数；等号右边第一项为沉积类粘结所提供的强度或刚度，第二项为成岩类粘结所提供的强度或刚度；在短期浸水作用下，沉积类粘结所提供的强度和刚度消失，等号右边第一项为零；在长期浸水作用下，成岩类粘结所对应的强度和刚度与时...

【专利技术属性】
技术研发人员：黄达，顾东明，杨芸芸，钟助，宋宜祥，
申请(专利权)人：河北工业大学，
类型：发明
国别省市：天津,12