3-采动断层活化分形界面效应的数值模拟研究.docVIP

采动断层活化分形界面效应的数值模拟研究 摘要　采动断层活化是采矿科学中的重要研究课题,以往的研究没有揭示断层面几何形态对它的影响.针对地质断层面的分形性质,构造了自仿射意义下的分形曲线及含有这种分形性质断层面的地质采矿模型,用数值方法模拟地下开挖引起分形断层面活化的现象,研究和总结了分形断层面对采动岩体沉陷的影响规律.研究表明,采动断层活化具有明显的分形界面效应;采动断层剪切滑移力学行为与其分形维数是相关的,不同分形维数的断层面采动后,导致其附近岩体中的应力场和位移场也不同. 关键词　采动断层活化　分形界面效应　分形维数 采动断层(节理)活化是矿山开采沉陷学科中典型的非线性力学问题,是指岩体中的断层(节理)在采动影响下产生滑移的现象.断层(节理)广泛存在于含煤地层中,一些断层中有充填物,一些断层中无充填物或很少.对于后者,其采动剪切滑移力学特征明显受控于断层(节理)面的几何结构形态,而以往的研究是把断层(节理)面看做光滑表面,显然不符合实际.研究表明[1~5],地质断裂面具有明显的统计分形性质,且表现为各向异性特征;进一步研究表明:地质断裂面是自仿射分形[6~8].实验证明,岩石(体)断裂面的剪切力学行为与其分形维数具有明显的相关性[9].所以,我们借助数值模拟的方法研究和验证断层(节理)采动剪切滑移力学特征与其分形维数的相关性,总结了分形断层面对采动岩体沉陷的影响规律. 1　数值模拟基本工作 采用Weierstrass-Mandelbrot函数构造自仿射分形曲线,函数表达式为 取=0, b=1.5并考虑其正弦函数,得 式中,D为分维,其值为D∈[1, 2]. 分别以D=1·1, 1·2, 1·3, 1·4, 1·5构造了5条自仿射意义下的分形曲线,并预设于5个地质采矿模拟原型中.模拟原型:煤层上履岩层厚为50 m,煤层底板岩层厚为10 m,地层总厚为60 m,走向长200 m,煤层为水平煤层,采厚1·6 m,整个地层岩性单一,全部为砂岩,断层倾角为65°,采空区在断层上盘一侧,在整个计算中,所有参数均不发生变化,只研究不同维数(D=1·1, 1·2, 1·3, 1·4, 1·5)的分形断层面受地下开采影响时,地表下沉、水平移动及采空区上方应力的不同分布规律. 开采分两步进行,从右向左首先采出40 m,然后继续推进到60 m (图1、图2)[10]. 图1　Y方向的拉应力分布 Fig·1　Distribution of tensile stresses in theYdirection 图中应力单位为105Pa (图2同) 图2　Y方向的压应力分布 Fig·2　Distribution of compressive stresses inYdirection 2　采动断层活化分形界面效应数值模拟结果分析 2·1　分形断层面对采动岩体应力分布和破坏规律的影响 本次模拟获得了含不同分形维数断层面的岩体在采动影响条件下煤层上覆岩层内的应力分布状态,覆岩破坏性态和地表移动破坏状况.当采宽为40 m时,采动影响没有诱使断层面活化, 5个模型的应力分布状态和覆岩的破坏性态基本相同,均表现为:在靠近采空区的上覆层岩中产生双向拉伸破坏,而形成冒落区,并逐渐向上发展为单向(竖向)拉伸,易产生离层,而在采空区两侧煤体的上方产生横向拉伸、竖向压缩的支承压力区,采空区两侧竖向压应力分布大体对称,形态基本相同,在地表移动边界附近的岩体内部产生了横向受拉、竖向受压的边界破坏区,在此区域内易产生地表裂缝.地表上各点的位移维持正常的分布状态,断层两侧的节点没有发生不连续的变形,只产生连续变形. 当采宽达到60 m时,采动影响波及了整个分形断层面, 5个模型上的分形断层面附近均发生了应力集中现象,但因断层面分形维数的不同而有所差异,为了清楚起见,我们分离了应力分布,提取X方向的拉、压应力分布,Y方向的拉、压应力分布图(图1、图2).从5个模型的应力分布状态来看,断层面的分形几何形态对X方向的拉应力影响不大,但受拉区域的应力值较采宽为40 m时有所增大,X方向的压应力分布则随着断层面分形几何形态的不同而有所差异:随着分形维数的增加,靠近地表附近的X方向压应力集中程度逐渐减弱.说明断层活化量越来越小,从而造成断层面活化后其两侧岩体的接触程度越来越弱.Y方向的拉应力在分形断层面的下部集中程度大体一致,而上部随分形维数的增大,集中程度越来越小,说明断层面的起伏形态对岩体Y方向的拉应力有一定影响.由于开采的影响,致使断层面活化,随着分形维数的增加,断层面两侧岩体的互相咬合程度越来越强,断层活化时的阻力越来越大,因此下部岩体移动值及向上传递的量值越来越小,从而造成断层面上部的拉应力集中程度越来越小.采空区正上方岩体的受拉区域随着分形维数的增加,范围有所增大,而且受拉位置