发明名称 一种采煤沉陷灾变发生的预警方法
摘要 本发明公开了一种采煤沉陷灾变发生的预警方法,利用关键层和软岩的极限破断距及其下部的随工作面推进距离的增大自由空间高度的发育情况来判断导水裂隙带发育高度,然后利用导水裂隙带高度评估采煤沉陷灾变发生的可能性,随着工作面推进距离的增大导水裂隙带高度持续发育,当导水裂隙带高度接近导水裂隙带临界高度时发出预警:工作面距离如果进一步增大,就会发生采煤沉陷灾变,否则是安全的。本发明通过采煤沉陷灾变预警,预防地下水沿着采动裂隙加速向采空区或深部岩体渗透,造成地下水位降低,进而形成井下突水或淹井灾害事故。一方面最大限度的生产出宝贵的煤炭资源,另一方面有效保护了采煤工作的生产安全。
申请公布号 CN105332738A 申请公布日期 2016.02.17
申请号 CN201510532625.7 申请日期 2015.08.24
申请人 西安科技大学 发明人 孙学阳;夏玉成
分类号 E21F17/18(2006.01)I 主分类号 E21F17/18(2006.01)I
代理机构 代理人
主权项 一种采煤沉陷灾变发生的预警方法,其特征在于,包括如下步骤:S1、确定关键层在覆岩中的位置;S11、由下往上确定覆岩中的坚硬岩层位置,假设第1层岩层为坚硬岩层,其上直至第m层岩层与之协调变形,第m+1层岩层不与之协调变形,则第m+1层岩层是第2层坚硬岩层,由组合梁原理可导出作用在第1层硬岩层上的载荷为:<img file="91793dest_path_FSB0000148576590000011.GIF" wi="1232" he="256" />式中:q<sub>1</sub>(x)|<sub>m</sub>为考虑到第m层岩层对第1层坚硬岩层形成的载荷;h<sub>i</sub>、γ<sub>i</sub>、E<sub>i</sub>分别为第i岩层的厚度、容重、弹性模量(i=1,2,...,m);考虑到第m+1层,对第1层坚硬岩层形成的载荷为:<img file="8934dest_path_FSB0000148576590000012.GIF" wi="1366" he="263" />由于第m+1层为坚硬岩层,其挠度小于下部岩层的挠度,第m+1层以上岩层已不再需要其下部岩层去承担它所承受的载荷,则必然有:q<sub>1</sub>(x)|<sub>m+1</sub><q<sub>1</sub>(x)|<sub>m</sub>            (3) 将式(1)、(2)代入式(3)并化简可得判别坚硬岩层位置的公式:<img file="788671dest_path_FSB0000148576590000013.GIF" wi="1357" he="121" />判别时,从煤层上方第1层岩层开始往上逐层计算,当满足式(4)时,则不再往上计算,此时从第1层岩层往上,第m层岩层为第1层硬岩层,从第1层硬岩层开始,按上述方法确定第2层硬岩层的位置,以此类推,直至确定出 最上一层硬岩层(设为第n层硬岩层),得覆岩中硬岩层的位置及其所控软岩层组;S12、计算各硬岩层的破断距,由固支梁力学模型,根据材料力学理论分析可知,梁内任一点的正应力为:<img file="183880dest_path_FSB0000148576590000021.GIF" wi="996" he="113" />式中:M为任一点所在截面弯矩,(KN m);y为任一点与截面中性轴的距离,(m);h为岩梁厚度,(m);由对固支梁的分析可知,固支梁最大弯矩发生在梁的两端,即:<img file="998252dest_path_FSB0000148576590000022.GIF" wi="1352" he="122" />所对应的最大拉应力为:<img file="820715dest_path_FSB0000148576590000023.GIF" wi="1369" he="130" />当σ<sub>max</sub>=σ<sub>t</sub>时,岩梁断裂,其极限跨距为:<img file="150065dest_path_FSB0000148576590000024.GIF" wi="868" he="147" />式中:h<sub>k</sub>为第k层硬岩层的厚度(m);σ<sub>k</sub>为第k层硬岩层的抗拉强度(MPa);q<sub>k</sub>为第k层硬岩层承受的载荷(kN/m<sup>2</sup>);由(1)式可知,q<sub>k</sub>可按下式确定:<img file="348965dest_path_FSB0000148576590000025.GIF" wi="1172" he="275" />式中,下标k代表第k层硬岩层;下标j代表第k层硬岩层所控软岩层组的分层号;m<sub>k</sub>为第k层硬岩层所控软岩层的层数;E<sub>k</sub><sub>,</sub><sub>j</sub>,h<sub>k</sub><sub>,</sub><sub>j</sub>,γ<sub>k</sub><sub>,</sub><sub>j</sub>分别为第k层硬岩层所控软岩层组中第j层岩层弹性模量、分层厚度及容重;由于表土层的弹性模量可视为0,设表土层厚度为H,容重为γ,则最上 一层硬岩层即第n层硬岩层上的载荷可按下式计算:<img file="17844dest_path_FSB0000148576590000031.GIF" wi="1357" he="292" />S13、按以下原则对各硬岩层的破断距进行比较,确定关键层位置:①第k层硬岩层若为关键层,其破断距应小于其上部所有硬岩层的破断距,即满足L<sub>k</sub><L<sub>k+1</sub>        (11) ②若第k层硬岩层破断距L<sub>k</sub>大于其上方第k+1层硬岩层破断距,则将第k+1层硬岩层承受的载荷加到第k层硬岩层上,重新计算第k层硬岩层的破断距;③从最下一层硬岩层开始逐层往上判别L<sub>k</sub><L<sub>k+1</sub>是否成立,及当L<sub>k</sub>>L<sub>k+1</sub>时重新计算第k层硬岩层破断距;S14、通过下式计算关键层破断距:<img file="11208dest_path_FSB0000148576590000032.GIF" wi="1006" he="147" />其中ψ=PT;S2、确定所夹的软岩极限破断距;S21、设其挠曲方程为:<img file="765537dest_path_FSB0000148576590000033.GIF" wi="1664" he="145" />则其满足边界条件:ω|<sub>x</sub><sub>=</sub><sub>0</sub>=0,ω|<sub>x</sub><sub>=</sub><sub>l</sub>=0,<img file="502549dest_path_FSB0000148576590000034.GIF" wi="244" he="133" /><img file="25934dest_path_FSB0000148576590000035.GIF" wi="326" he="188" />利用伽辽金法可得:<img file="455778dest_path_FSB0000148576590000036.GIF" wi="1100" he="141" />通过解算,可得固支梁的最大挠度为:<img file="697404dest_path_FSB0000148576590000041.GIF" wi="954" he="127" />式(13)可用通项公式表示为:<img file="972528dest_path_FSB0000148576590000042.GIF" wi="1305" he="130" />式中:E为梁的弹性模量,I为惯性矩,<img file="615998dest_path_FSB0000148576590000043.GIF" wi="177" he="120" />由式(16)可得弯曲梁转角方程为:<img file="951165dest_path_FSB0000148576590000044.GIF" wi="1385" he="132" />进而,由式(17)得到梁的曲率方程为:<img file="680086dest_path_FSB0000148576590000045.GIF" wi="1369" he="123" />固支梁弯曲后其水平拉伸变形为:<img file="758901dest_path_FSB0000148576590000046.GIF" wi="621" he="123" />其中,y为横截面上某点到中性层的距离,所以固支梁弯曲后其截面上某点的水平拉伸变形为:<img file="256878dest_path_FSB0000148576590000047.GIF" wi="1293" he="127" />分析式(19)可知,当<img file="90842dest_path_FSB0000148576590000048.GIF" wi="406" he="113" />时,水平拉伸变形最大,最大值为:<img file="41481dest_path_FSB0000148576590000049.GIF" wi="1363" he="132" />式中,当n→∞时,<img file="923986dest_path_FSB00001485765900000410.GIF" wi="388" he="132" />所以水平拉伸变形的最大值为:<img file="276470dest_path_FSB00001485765900000411.GIF" wi="776" he="113" />由于在弯曲梁的中性轴面下端面产生水平拉伸变形,并且在梁的下端即y=h/2的端面上水平拉应变值最大;因此,只要此处的ε<sub>max</sub>不超过岩层的临界 水平拉伸应变值,那么该岩层就不会产生导水裂隙;所以:<img file="219018dest_path_FSB0000148576590000051.GIF" wi="774" he="114" />取泥岩等较软弱类岩层的临界水平拉伸变形值为1.0mm/m,由式(23)可以得到岩梁受力弯曲产生最大水平拉伸应变值时的跨距:<img file="656953dest_path_FSB0000148576590000052.GIF" wi="727" he="130" />此时所对应的该岩层下端面到煤层顶板的距离H就是导水裂隙带发育的高度,所对应的工作面推进距离为:<img file="77570dest_path_FSB0000148576590000053.GIF" wi="1033" he="151" />式中:φ<sub>q</sub>、φ<sub>h</sub>分别为岩层的前、后方断裂角;S3、通过下式计算每层岩层下的自由空间高度:<img file="550139dest_path_FSB0000148576590000054.GIF" wi="950" he="144" />式中:Δ<sub>i</sub>为第i层岩层下的自由空间高度;M为煤层开采厚度;h<sub>j</sub>为第j层岩层的厚度;k<sub>j</sub>为第j层岩石的残余碎胀系数;当工作面推进到其可产生最大拉伸的距离时,岩层的最大挠度为:<img file="663589dest_path_FSB0000148576590000055.GIF" wi="780" he="120" />此时,如果软岩层的最大挠度大于其下部自由空间高度,由于自由空间的限制,软岩层将保持塑性状态而不会破坏,导水裂隙带至此不在向上发展,于是有:ω<sub>i</sub><sub>,</sub><sub>max</sub>>Δ<sub>i</sub>        (28) 即:<img file="588820dest_path_FSB0000148576590000056.GIF" wi="921" he="132" />反之,当软岩层的最大挠度小于其下部自由空间高度,将发生破断并导水。于是有:ω<sub>i</sub><sub>,</sub><sub>max</sub><Δ<sub>i</sub>           (30) 即:<img file="813128dest_path_FSB0000148576590000061.GIF" wi="1008" he="135" />S4、利用关键层和软岩的极限破断距及其下部的自由空间高度来判断导水裂隙带高度:①当关键层的悬露小于其极限破断距时,导水裂隙带不会向上发展;此时结构关键层起到隔水作用,是隔水关键层;②当关键层的悬露大于其极限破断距时,如果不存在自由空间高度,导水裂隙带将终止,否则,继续向上发展,结构关键层失去隔水功能;③当软岩层的水平拉伸应变小于其水平极限拉伸应变时,导水裂隙带不会向上发展,该软岩就是隔水关键层;④当软岩层的水平拉伸应变大于其水平极限拉伸应变时,如果岩层的最大挠度大于其下部自由空间高度,导水裂隙带终止,否则,继续向上发展;S5、利用导水裂隙带高度评估采煤沉陷灾变发生的可能性,当导水裂隙带高度接近导水裂隙带临界高度时发出预警:工作面距离如果进一步增大,就会发生采煤沉陷灾变,否则是安全的,即开采不会引起采煤沉陷的发生。
地址 710054 陕西省西安市雁塔中路58号