| 返回中国矿井物探网 | 返回华虹智能 |
 
Email:用户名 @ kjwt.cn 密码
 
新闻检索
往期回顾
您现在的位置:第60期 -> 第四版 -> 新闻内容
并行电法技术与数值模拟在采煤面覆岩破坏发育规律中的应用
新闻作者:张周鑫  发布时间:2013-04-24  查看次数:  放大 缩小 默认
安徽理工大学 张周鑫
  摘要:为实现工作面安全开采,有必要对其上覆岩层运动特征进行研究,确定“两带”发育特征。本文结合工程背景,综合运用并行电法测试技术和UDEC数值模拟计算的方法,得出覆岩受采动影响破坏变形特征。电法测试技术基于岩层电性差异,记录不同时间段覆岩电阻率参数变化,并根据计算及图形模拟结果,反映覆岩破坏的变化状况。UDEC模拟基于离散元思想,通过建立地质模型并赋予相应实测参数,计算得出的结果反映采场上覆岩体随工作面推进变形破坏特征,确定“两带”范围与发育形态。综合利用物理模拟和数值模拟的技术手段能更好地为生产服务,并希望广泛应用于工程实践。
  关键词:覆岩破坏;并行电法技术;UDEC;安全生产
  1 引言
  采煤面覆岩变形与破坏发育规律探查与研究一直是矿井安全生产十分关注的问题,对矿井顶板管理及水害防治具有重要的指导意义,在这方面我国学者已经进行了较深入的研究并取了重要成果,代表性的有“上三带”理论[1]、关键层理论[2]、“三图-双预测”法[3]等。这些理论、方法为研究煤矿开采引起覆岩破坏“两带”发育特征[4-5]提供了很好的指导作用。然而,采动影响下的覆岩变形破坏是一个极其复杂的问题,对此还需深入研究;近年来,利用地球物理方法和计算机技术研究分析覆岩破坏产生裂隙的发育状况,初步取得了比较明显效果。本文针对煤层开采前后覆岩变形和破坏引起岩层电阻率变化这一特征,在煤层顶板布置电极,观测电阻率变化,从而反演顶板覆岩破坏范围与形态[6-8];同时,将数值模拟方法应用到预测覆岩“两带”发育形态特征上。UDEC[9-11]数值模拟软件是在离散单元法基础上编制的数值计算程序,基于的原理是将岩体视为由离散的岩块和岩块间的节理面所组成,允许岩块平移、转动和变形,而节理面可被压缩、分离或滑动,较适用于煤炭开采引起上覆岩体滑移、变形、破裂等非线性问题的研究分析。
  2 并行电法测试技术在覆岩破坏“两带”发育规律研究中应用
  2.1 并行电法基本原理
  并行电法探测原理基于静电场理论,采用阵列布极方法,以目标体电性差异为前提,集电剖面与电测深一体,进行二维地电断面测量;根据记录的电阻率值,进行计算、处理、分析,从而用于划分地层、圈闭异常、确定“两带”发育等方面。并行电法是改进了的高密度电法,并行采集数据,一般采用单极供电(AM法)和偶极子供电(ABM法)两种方式来进行数据采集和处理,还可以远程控制及在线观测,较传统方法优势在于很大程度提高采集效率。
  2.2 岩层电阻率解释
  岩体的结构特征是影响电阻率的主要因素之一,两者有着显著的相关性。通常不同岩性电阻率值有一定差别,同一岩层,由于其结构特征发生变化,使得电阻率值也会发生改变。对于煤层顶板岩层来说,从电性上分析煤层电阻率值相对较高,砂岩次之,黏土岩类最低。当岩层发生变形与破坏时,若岩层不含水,则其导电性变差,局部电阻率值增高。采矿过程中岩层电性在纵向和横向上发生变化,代表了其破坏和裂隙发育特征;因此,通过测取顶板不同深度岩层电阻率变化来分析其变形与破坏规律,这是电法测试覆岩破坏特征的地质基础。“两带”发育仅为一定高度以下岩层电阻率所发生的变化,每次测试时以岩层背景电阻率分布作为基础,通过动态测试可以从时空规律上直观分析岩层的破坏过程和规律。
  2.3 岩层背景电阻率分布
  覆岩破坏探测中,背景场电阻率分布是在工作面距测孔较远,一般距离120m以上时的岩层电阻率分布,即未受采动影响岩层电阻率分布。图一为某矿工作面开采前一段时间测试电阻率成像结果剖面,测试时回采工作面距监测范围较远,离孔口130m以上,电阻率值分布集中在<2000Ω•m范围内,其中泥岩段电阻率最低,细砂岩和粉砂岩段电阻率值稍高。其地层岩性变化特征较为明晰,为正常岩层电性反应,局部岩层电阻率的高低变化反映出岩性变化或岩体完整性的不同,可作为背景值,为后续采动影响时岩层变形与破坏电阻率值对比提供基础。  

图一 电阻率背景值
  2.4 覆岩变形破坏“两带”发育特征分析
  据于师建等通过相似材料实验模拟表明:垮落带内岩层破碎垮落,岩石碎胀造成的视电阻率变化是正常情况的4-6倍,导水裂缝带范围内岩层裂隙普遍发育,横向、垂向裂隙交错,视电阻率是正常情况的2-3倍。冒落带和导水裂缝带合称为“两带”,其上部的弯曲下沉带内,由于下部岩层垮落,上部岩层会逐个向下弯曲下沉,在岩层界面附近常有裂隙或真空离层发育,随着上部岩层的下沉,先形成的裂缝或离层会有所闭合,岩层少有裂隙发育。根据钻孔电法监测结果,对于探测区内的“两带”及离层的发育过程及形态进行有效判断。图二为随工作面向前推进过程中不同时段电阻率剖面的节选图。  

(a)   

(b)

(c)
  图二 覆岩变形与破坏电阻率分布
  从节选的图上可以看出:随工作面的不断推进,监测的电阻率剖面变化规律明显,煤层顶板煤岩体结构受到应力破坏后产生裂隙或垮落,表现为电阻率升高的特点,随着工作面向监测孔口推进,高阻变化也动态的向孔口移动。据于师建老师判定“两带”发育高度和形态特征的总结,图上所示不同时间的测试结果可连续动态分析采场覆岩变形与破坏基本规律,为矿井安全生产提供有效技术参数。
  3 煤层采动顶板UDEC离散元数值模拟
  3.1 离散元模拟基本思想
  UDEC模拟软件最早是由P.A.Cundall于1971年提出的一种不连续数值方法模拟,这种方法明显的优点是适用于模拟节理系统在准静态或动态条件下的变形过程。离散单元法以受节理裂隙切割或分离的块体为出发点,将研究区域划分为单元。单元因受节理等不连续面控制,在运动过程中,单元结点之间可以接触,也可分开,任一块体单元总会受到相邻单元的力和力矩的作用,正是在其合力和合力矩的作用下产生运动。以每个单元刚体运动方程为基础,在建立描述整个破坏状态的显式方程组之后,根据牛顿第二运动定律并结合不同的本构关系,以动态松弛法进行迭代计算。单元之间相互作用力的大小可以根据力和位移的关系求出,而个别单元的运动则完全按该单元所受的不平衡力和不平衡力矩的大小,按牛顿运动定律确定。
  在UDEC中,为完整块体和不连续面开发了几种材料特性模型,用来模拟不连续地质界面可能显现的典型特性,是基于拉格朗日算法能够很好地模拟块体系统的变形和大位移,因而离散元在矿山岩体力学和矿山压力研究中得以较广泛的应用。
  3.2 数值模型设计原则
  地质模型的建立是模拟计算的前提,为保证较真实地反应现场实际,在离散元模型设计过程中,必须遵循一定的原则:
  (1)影响采动覆岩移动的因素很多,模型的设计必须突出影响采动覆岩移动的主要因素,并尽可能多地考虑其它重要因素。
  (2)模型是由实体简化,体现一般性,模型的设计必须能很好地反映材料的物理力学特性与结构面特性。
  (3)地下煤层开挖是半无限域范围,数值模拟只能是在有限的范围内进行。因此,模型的设计必须考虑其边界效应和初始地应力作用。
  (4)模型的设计,必须能体现工作面的推进与接续,能体现出覆岩“两带”随工作面向前推进变化过程,需具有适合的尺寸。
  模型的边界条件为:左右边界约束水平位移,下部边界约束竖向位移,模型的上部条件为自由边界,允许表面及内部所有单元体产生变形和移动,垂向上施加自重应力,见图三。  

图三 初始模型及边界
  3.3 数值模型建立
  本次离散元计算模型以现场钻孔勘探资料为基础,更好地揭示工作面开采时覆岩层移动和地应力变化规律,其中以某矿13-1煤层为具体研究背景,定义一个工作面为模型框架,建立计算模型。
  根据开采沉陷学相关公式及其数值模型试算,确定数值计算区域尺寸为:走向模型400m*120m,定义开挖步距20m,共开挖140m,煤岩层物理力学参数具体见下表一、二。  
表一 材料物理力学参数

表二 节理面力学参数
  3.4 模拟计算结果
  设计模型开挖7步,步距20m,切眼留设煤岩柱100m,通过计算机135761步运算,输出结果图,图件综合反应了模型开挖、破坏、变形、稳定的整个动态过程。为便于分析,从中选取几个代表时间段(40m、80m、140m)进行解释说明,具体见图四-图九。  

图四 开挖40m屈服破坏形态  

图五 开挖80m屈服破坏形态  

图六 开挖140m屈服破坏形态
  图中所示棕色﹡代表岩体屈服面,红色○区域代表岩体正受到拉张断裂,绿色×代表塑性破坏区。从数值模拟图可以看出,天然岩体受采动影响,采场覆岩表现出屈服和破坏的特征。
  图四为模型开挖40m的形态,由此可以看出此时煤层直接顶岩层内部出现了很多的裂隙,有小范围的破坏,顶板有所下沉,顶板的上部岩层基本未受影响。随着开挖进行到80m(图五),顶板下沉量加大,裂隙进一步向上扩展,原先开挖垮落部分顶板裂隙发育稳定,局部有压实、节理闭合的情况。随着工作面继续推进采场覆岩破坏形态未有较大的变化,趋于稳定形态,工作面上方及前方裂隙最为发育,高度也最大,整个冒落裂隙范围大致呈两端高中间低的“马鞍”形分布(图六)。  

图七 开挖40m覆岩应力等值线图  

图八 开挖80m覆岩应力等值线图  

图九 开挖140m覆岩应力等值线图
  从图七-图九覆岩应力变化等值线图上可以看出:煤层的开采破坏了岩土体内的天然应力平衡状态,使采空区周围岩土体的应力发生重新分布,上覆岩体最大主应力的方向由开采前的铅直方向改变为向采空区两侧的偏斜,形成在采空区两侧的拱型支撑压力增高带;采空区中部受上覆岩体的垮落、充填和压实,形成稳定的支撑拱,应力分布均匀,总体应力分布形态呈现出均匀对称的“之”字型特征。
  由上述分析可知,在开采过程中,随着工作面的推进,每一个岩层均经历了变形、离层、断裂和冒落的过程,离层伴随着下部岩层的冒落发展由下向上逐步发展。每一个岩层的冒落是随着开采范围扩大分段形成的,具有周期性,每一个冒落形成后,均会达到一个新的平衡状态。在达到充分开采前,离层的起裂、扩展、闭合和分布是一个动态的时空演化过程。最终达到稳定状态,覆岩“两带”形态上呈现出“马鞍形”特征;应力在切眼和收作线附近达到最大值,表现出“之”字型特征。
  4 结论
  (1)并行电法测试技术对煤层采动引起的顶板岩层破坏规律认识具有针对性,能够较准确的确定覆岩“两带”发育特征,且其数据远程采集和动态监测等特点,使其具有显著地推广价值。
  (2)UDEC数值模拟软件可以很好地应用于煤层开采覆岩破坏特征的研究中,该软件能较好与地下工程特点结合,真实反应其动态变化。
  (3)两种技术方法同时应用于顶板预测,来验证和对比分析两种不同方法得出的规律性与统一性认识,可以丰富该类问题研究的内容。
  参考文献:
  [1]钱鸣高,缪协兴.岩层控制的关键层理论[J].煤炭学报,1996.
  [2]武强,黄晓玲,董东林.评价煤层顶板涌(突)水条件的“三图-双预测法”[J].煤炭学报,2000.
  [3]侯忠杰.组合关键层理论的应用研究及其参数确定[J].煤炭学报,2001,26(6).
  [4]许延春,李俊成,刘世奇等.综放开采覆岩“两带”高度的计算公式及适用性分析[J].煤矿开采,2011.
  [5]邹海,桂和荣,王桂梁等.综放开采导水裂隙带高度预测方法[J].煤田地质与勘探,1998.
  [6]张平松,刘盛东,吴荣新等.采煤面覆岩变形与破坏立体电法动态测试[J].岩石力学与工程学报,2009.
  [7]张平松,刘盛东,舒玉峰.煤层开采覆岩破坏发育规律动态测试分析[J].煤炭学报,2011.
  [8]吴荣新,方良成,周继生.采用网络并行电法仪探测采煤工作面无煤区[J].安徽理工大学学报(自然科学版),2007.
  [9]陈陆望,桂和荣,李一帆.UDEC模拟厚松散层及超薄覆岩条件下开采防水煤柱覆岩突水可能性[J].水文地质工程地质,2007.
  [10]郝志勇,林柏泉,张家山等.基于UDEC的保护层开采中覆岩移动规律的数值模拟与分析[J].中国矿业,2007.
  [11]李忠建,魏久传,施龙青等.浅埋煤层开采数值模拟及顶板突水危险性分析[J].煤矿安全,2011.
  作者简介:
  张周鑫,男,安徽理工大学在读研究生,主要从事工程地质研究。
 
上一篇 下一篇
  | 联系我们 | 内部员工系统 | Copyright? 2012-2015,All Rights Reserved 闽ICP备11010421号  
主办单位:中国矿业大学地球物理研究所 安徽理工大学物探研究中心 安徽惠洲地下灾害研究设计院 福州华虹智能科技股份有限公司
域名:www.kjwt.cn◆◆◆中国矿井物探网.cn ◆◆◆ www.kjwt.org