吳詩(shī)勇,時(shí)帥印,黃 楷
(安徽理工大學(xué)地球與環(huán)境學(xué)院,安徽 淮南 232001)
隨著國(guó)家經(jīng)濟(jì)快速發(fā)展,煤炭等礦產(chǎn)資源需求日益緊張[1-2],而地下煤層賦存條件較為復(fù)雜,不易開(kāi)采,事故頻發(fā)[3-4]。在煤炭開(kāi)采過(guò)程中絕大多數(shù)突水事故是由隱伏斷層引起的,斷層作為重要的突水通道,許多專家學(xué)者進(jìn)行了長(zhǎng)期的大量的研究[5-8]。文獻(xiàn)[9]通過(guò)建立不同傾角斷層的三維地質(zhì)模型,研究得出小傾角正斷層比大傾角正斷層對(duì)底板突水影響大等結(jié)論。文獻(xiàn)[10]建立主、次斷層的三維地質(zhì)模型,模擬采動(dòng)過(guò)程中對(duì)斷層活化的影響,研究得出次生斷層的存在易造成斷層活化等結(jié)論。因此,在礦井生產(chǎn)安全的前提下,如何合理留設(shè)防水煤巖柱的尺寸,提高開(kāi)采工作面煤炭回收率,顯得愈益重要。
本文根據(jù)楊柳礦10412工作面實(shí)際地質(zhì)條件,采用FLAC3D來(lái)模擬采動(dòng)條件下斷層活化對(duì)煤層頂、底板及采空區(qū)的影響,確定留設(shè)斷層保護(hù)煤柱的合理寬度,為礦區(qū)安全生產(chǎn)及防突水提供理論依據(jù)和參考。
10412工作面位于楊柳煤礦井田中部,介于戴廟斷層與小沈家斷層之間(見(jiàn)圖1)。10412是10煤采區(qū)工作面。 10煤層原始沉積較好,以單一煤層出現(xiàn),煤層厚度0~7.97m,平均厚2.95m。該工作面煤巖呈近S-N向延伸,傾向東,傾角1°~10°。煤層頂、底板以中砂巖、粉砂巖為主。工作面采用自由垮落的方式管理頂板。
圖1 楊柳煤礦構(gòu)造簡(jiǎn)圖
南部小沈家斷層呈NE-SW向延伸,傾向SE,傾角在50°~70°之間,延伸長(zhǎng)度約10km,落差20~220m。斷層上盤以細(xì)砂巖為主,下盤以泥巖為主;斷層面有擦痕,有少量淋水。
模擬區(qū)位于10412工作面南部,包含小沈家斷層。選取模型長(zhǎng)300m,沿巖層走向;寬250m,沿?cái)鄬友由旆较?;?20m,沿豎直方向。
根據(jù)小沈家斷層展布情況,本次地質(zhì)模型選取的斷層傾角為60°,落差30m,斷層帶寬5m。各巖層按水平地層處理,煤層厚度選取3m。劃分模型網(wǎng)格總數(shù)49 410,網(wǎng)格節(jié)點(diǎn)數(shù)53 816(見(jiàn)圖2)。
圖2 三維地質(zhì)模型
模擬采用摩爾-庫(kù)倫本構(gòu)模型。四周及底面固定,頂部邊界為活動(dòng)面。模型埋深440m,上覆巖層平均密度取2 280kg/m3,則上部巖層等效荷載為9.8MPa。頂、底面只受上覆巖層荷載的垂向作用,四周無(wú)剪切力。
相關(guān)巖石力學(xué)參數(shù)如表1所示。
表1 巖石力學(xué)參數(shù)
工作面從煤層走向100m處開(kāi)挖,向斷層方向推進(jìn)。分五次開(kāi)挖,每次開(kāi)挖20m。
最大主應(yīng)力云圖(見(jiàn)圖3a~圖3d)顯示,應(yīng)力變化主要表現(xiàn)在豎直方向,從采空區(qū)到工作面,大致形成卸壓區(qū)、應(yīng)力集中區(qū)、初始狀態(tài)區(qū)3個(gè)區(qū)域。卸壓區(qū)主要分布在采空區(qū)四周,其中頂板應(yīng)力明顯高于底板;應(yīng)力集中區(qū)出現(xiàn)在煤壁附近,向外圍逐漸過(guò)渡到初始應(yīng)力區(qū)。隨著煤層開(kāi)采向前推進(jìn),應(yīng)力集中區(qū)隨之前移,卸壓區(qū)逐漸擴(kuò)大,工作面與斷層之間的初始應(yīng)力區(qū)逐漸縮小。當(dāng)開(kāi)挖90m時(shí),應(yīng)力影響范圍靠近斷層帶,頂板應(yīng)力集中區(qū)最先抵達(dá)斷層帶(見(jiàn)圖3e),斷層下盤近斷層處應(yīng)力未超過(guò)煤巖強(qiáng)度。繼續(xù)開(kāi)挖5m后,斷層附近應(yīng)力增大,導(dǎo)致水平方向應(yīng)力超過(guò)巖石的抗拉強(qiáng)度(見(jiàn)圖3f),出現(xiàn)明顯的導(dǎo)水裂隙面,造成工作面附近巖石破碎區(qū)與斷層破壞區(qū)相連通,形成導(dǎo)水通道。
(a)開(kāi)挖20m (b)開(kāi)挖40m
圖4為煤層頂板2m處工作面開(kāi)挖方向應(yīng)力分布,坐標(biāo)原點(diǎn)處表示斷層位置。工作面推進(jìn)過(guò)程中,卸壓區(qū)向斷層方向逐漸擴(kuò)展, 煤壁前后集中應(yīng)力逐漸增大, 從初始應(yīng)力11.20MPa達(dá)到20.4MPa, 且煤柱寬度逐漸減小, 直至出現(xiàn)溝通斷層帶的導(dǎo)水裂隙,此時(shí)垂向應(yīng)力值達(dá)到20.61MPa,此后則逐漸減小。
圖4 頂板應(yīng)力分布圖
工作面的開(kāi)挖,造成采空區(qū)圍巖擾動(dòng)明顯,開(kāi)始出現(xiàn)塑性變形。從模擬變形云圖(見(jiàn)圖5)來(lái)看,巖石變形區(qū)呈“馬鞍”型,相對(duì)于開(kāi)挖面大致呈對(duì)稱分布。當(dāng)工作面推進(jìn)20m時(shí)(見(jiàn)圖3a),頂、底板圍巖破壞區(qū)域小,主要集中在采空區(qū)中部,其前后煤壁均發(fā)生塑性變形,且頂板變形大于底板。隨著工作面不斷推進(jìn)(見(jiàn)圖5b~圖5d),頂、底板變形加劇,采煤工作面附近塑性區(qū)不斷擴(kuò)大,頂板塑性區(qū)與斷層帶之間的有效隔水寬度越來(lái)越小。當(dāng)工作面推進(jìn)90m時(shí)(見(jiàn)圖5e),采空區(qū)頂、底板破壞嚴(yán)重,頂板塑性區(qū)抵達(dá)斷層;在此基礎(chǔ)上,再往前推進(jìn)5m(見(jiàn)圖5f),頂板破壞帶與斷層帶導(dǎo)通,斷層導(dǎo)水,易引起煤層突水事故,此時(shí)工作面距離斷層帶36.53m。
(a)開(kāi)挖20m (b)開(kāi)挖40m
根據(jù)《煤礦防治水細(xì)則》[11],斷層區(qū)采煤,煤柱留設(shè)寬度計(jì)算公式為
(1)
式中:K為安全系數(shù),一般取2~5;M為煤層的厚度,m;P為斷層煤柱承受的水壓,MPa;FP為煤層的抗拉強(qiáng)度,MPa。
將K=4,M=2.95m,P=3.0MPa,F(xiàn)P=0.26MPa代入式(1)得L=34.12m。
從以上模擬結(jié)果來(lái)看,此臨界寬度必然會(huì)造成采動(dòng)引起斷層活化,導(dǎo)致開(kāi)采引起的頂板破壞帶與斷層帶導(dǎo)通,形成導(dǎo)水通道,容易引發(fā)突水事故。因此,10412工作面小沈家斷層防水煤柱留設(shè)寬度應(yīng)不小于36.53m。
本文主要采用數(shù)值模擬方法,對(duì)楊柳礦10412工作面小沈家斷層區(qū)在采動(dòng)條件下煤層頂、底板變形及應(yīng)力分布進(jìn)行了研究,并確定了小沈家斷層防水煤柱留設(shè)的合理寬度。
(1)沿工作面推進(jìn)方向,從采空區(qū)到斷層帶,應(yīng)力大致形成卸壓區(qū)、應(yīng)力集中區(qū)、初始狀態(tài)區(qū)三個(gè)區(qū)域;
(2)巖石變形區(qū)主要集中在采空區(qū)中部,呈“馬鞍”型,隨著工作面的推進(jìn),頂、底板圍巖破壞區(qū)域擴(kuò)大,塑性變形區(qū)逐漸靠近斷層帶;
(3)當(dāng)工作面推進(jìn)95m時(shí),集中應(yīng)力超過(guò)斷層帶圍巖強(qiáng)度導(dǎo)致下盤巖石出現(xiàn)裂縫,造成頂板破壞帶與斷層帶導(dǎo)通,形成導(dǎo)水通道,易引發(fā)突水事故;
(4)根據(jù)采動(dòng)過(guò)程中頂、底板變形及應(yīng)力分布特征,建議10412工作面小沈家斷層防水煤柱留設(shè)寬度應(yīng)不小于36.53m。