曹 健,高 斌,黃慶享
(1.內蒙古科技大學礦業(yè)與煤炭學院,內蒙古包頭 014010;2.西安科技大學能源學院,陜西西安 710054)
長壁工作面開采的上行裂隙是導水導氣的重要通道,研究上行裂隙的發(fā)育規(guī)律及其機理,對于工作面的安全開采具有重要意義[1-3]。隨工作面推進,開切側邊界上行裂隙和工作面?zhèn)葎討B(tài)上行裂隙不斷發(fā)育,其實質是在集中拉應力和垂向應力的作用下,超過巖層的極限強度時,在巖層頂端產生裂紋,裂紋擴展直至貫通該巖層,上行裂隙發(fā)育??梢?,正是微觀上裂紋的擴展導致了宏觀中上行裂隙的發(fā)育,因此,有必要從斷裂力學的角度來研究上行裂隙的發(fā)育機理。
目前,關于采動上行裂隙發(fā)育規(guī)律及其機理的研究,主要集中在2 個方面:①覆巖裂隙發(fā)育規(guī)律的研究;②覆巖裂隙發(fā)育高度的研究。黃慶享[4-5]揭示了淺埋煤層開采的上行裂隙發(fā)育規(guī)律;張文忠等[6]研究了局部充填開采充填參數(shù)對上行裂隙發(fā)育規(guī)律的影響;張玉軍等[7]得出采動巖體裂隙以高角度縱向上行裂縫為主;賈后省等[8]研究了淺埋薄基巖采煤工作面縱向裂隙的“張開-閉合”規(guī)律;黃慶享等[9]得出了覆巖裂隙主要為采空區(qū)邊緣向上發(fā)育的上行裂隙和垮落區(qū)內的離層裂隙;謝曉鋒等[10]認為煤層厚度對導水裂隙帶高度的影響最大;張勇等[11]分析了厚煤層上覆巖層采動裂隙擴展原理;王曉振等[12]得出了導水裂隙的發(fā)育高度同時受到采厚和覆巖關鍵層結構的影響;曹祖寶等[13]研究了黃隴煤田煤層開采覆巖結構對裂隙發(fā)育高度的影響規(guī)律;黃慶享等[14]提出了采空區(qū)上覆巖層的“冒落帶、塊體鉸接帶、似連續(xù)帶和彎曲下沉帶”的“四帶”劃分;車曉陽等[15]基于工作面寬度、巖性及巖層破斷角等因素,給出了裂隙帶發(fā)育高度的計算公式;白建平等[16]研究了工作面回采后覆巖裂隙發(fā)育高度及演化形態(tài);王創(chuàng)業(yè)等[17]研究了淺埋煤層開采的裂隙網絡演化規(guī)律;黃健豐等[18]分析了綜放開采巖層裂隙發(fā)育特征;陳輝等[19]研究了厚松散層薄基巖淺埋煤層開采的導水斷裂帶高度。曹健等[20-21]揭示了集中應力導致的裂隙發(fā)育擴展與煤柱布置之間的關系。以上研究多采用鉆孔實測、物理模擬等方法從宏觀上得到開采覆巖上行裂隙發(fā)育規(guī)律,而從微觀裂紋的產生、擴展直至貫通巖層自由面的全過程揭示上行裂隙發(fā)育機理有待于進一步研究。
1.1.1 物理相似模擬設計
霍洛灣煤礦22103 工作面開采2-2上煤,工作面長249 m,推進長度1 265 m,煤層傾角1°~3°。工作面采高2.5 m,地面標高1 180~1 254.5 m,煤層底板標高1 084~1 098.07 m,平均埋深107.6 m。
建立物理相似模擬模型,采用1:100 的相似比,模型尺寸為:3 m(長)×0.2 m(寬)×1.2 m(高),模型采用河沙作骨料,石膏、大白粉作為膠結材料,物理相似模型全景如圖1。
圖1 物理相似模型全景Fig.1 Physical simulation model
1.1.2 上行裂隙發(fā)育特征
隨采煤工作面推進,上覆巖層在集中拉應力和重力的作用下發(fā)生撓曲、垮落,上行裂隙隨之沿一定垮落角不斷向上發(fā)育。
走向方向上行裂隙的發(fā)育過程可以分為上行裂隙產生段、上行裂隙發(fā)育延伸段和上行裂隙穩(wěn)定段3 個階段,上行裂隙的發(fā)育過程分段如圖2。
圖2 上行裂隙的發(fā)育過程分段Fig.2 Development stages of upward fractures
上行裂隙發(fā)育是垂直作用力與集中拉應力綜合作用的結果,為研究上行裂隙發(fā)育機理,運用斷裂力學對裂紋的受力和擴展過程進行分析,巖層裂紋端部受力分析如圖3。
圖3 巖層裂紋端部受力分析Fig.3 End force analysis of strata fractures
裂紋端部垂直方向上受巖層自身重力G 和上覆載荷q 的作用,可導致面內剪切型裂紋(Ⅱ型裂紋);垂直于上行裂隙延伸方向受邊界集中拉應力σθ的作用,可導致拉伸型裂紋(Ⅰ型裂紋),2 種裂紋的受力和擴展形式如圖4。
圖4 2 種不同的裂紋型式Fig.4 Two different fracture types
可見,采動巖層裂紋的擴展是以上2 種裂紋綜合作用的結果,因此,在斷裂力學中應當作為混合型裂紋(Ⅰ、Ⅱ型裂紋)進行分析[22]。
工作面開采,垮落帶頂板呈無序垮落狀態(tài),該范圍內的上行裂隙巖層完全破斷。其上部的斷裂帶頂板呈有序垮落,在自重、上覆載荷和邊界集中拉應力的作用下回轉下沉,上行裂隙頂板受力為“上拉下壓”狀態(tài)。由于巖層隨采動呈分層垮落,為分析斷裂帶上行裂隙發(fā)育機理,確定其發(fā)育高度,將裂隙帶巖層自下而上進行編號(1,2,…,n),上行裂隙發(fā)育機理如圖5。
圖5 上行裂隙發(fā)育機理Fig.5 Development mechanism of upward fractures
裂隙帶巖層的破斷擴展從下至上逐層計算,受拉應力和垂直方向的作用力,超過其臨界強度時,首先在巖層的頂端產生裂紋,隨后在拉應力的作用下沿著裂隙帶方向向下擴展直至貫通該巖層,那么該巖層破斷,成為導水導氣的通道,同時上行裂隙向上發(fā)育延伸。
上行裂隙發(fā)育高度的確定如圖6。
圖6 上行裂隙發(fā)育高度的確定Fig.6 Determination of development height of upward fracture
根據(jù)上行裂隙發(fā)育機理,隨著計算層位增高,斷裂帶第n-1 層巖層裂紋仍能夠貫通巖層自由面時,該巖層完全斷裂,上行裂隙繼續(xù)向上發(fā)育;而計算到其上部的第n 層時,由于集中拉應力和垂直方向作用力減小,裂紋不足以擴展,在該巖層不會貫通,則該巖層則不會完全斷裂,此時,上行裂隙不再沿破斷角向上延伸。由此可以得到在充分采動條件下,基于微觀裂紋擴展理論的采動上行裂隙發(fā)育高度的確定方法:煤層頂部直至第n 層巖層底部的厚度即為煤層開采后上行裂隙的發(fā)育高度。
根據(jù)Erdogan 和Sih 提出的σ(θ)max理論[16],控制巖層裂紋斷裂的參數(shù)是裂紋端部的最大環(huán)向拉應力σ(θ)max,據(jù)此,建立的上行裂隙巖層裂紋擴展理論模型如7。
圖7 裂紋擴展理論模型Fig.7 Theoretical model of fracture extension
裂紋端部的應力狀態(tài)可以用式(1)來表示。
式中:τ 為A 點處的剪應力,MPa;c 為裂紋的半長,m。
裂紋在其端部沿徑向方向擴展,在σ(θ)max達到巖層的臨界強度因子時,裂紋開始擴展。式(1)可用式(4)和式(5)來表示:
根據(jù)式(4)和式(5),得到的σ(θ)max理論的啟裂跡線如圖8。
圖8 σ(θ)max 理論的啟裂跡線Fig.8 Beginning broken curve of σ(θ)max theory
裂紋擴展的判據(jù)如下:
1)根據(jù)式(2)、式(3)和式(5)求得裂紋的擴展啟裂角θ0。
2)根據(jù)式(2)、式(3)、式(4)和所求θ0值,結合圖8 可以判斷裂紋是否發(fā)生擴展(上行裂隙是否繼續(xù)發(fā)育)。若(KⅠ/KⅠc,KⅡ/KⅠc)處于啟裂跡線內側,則裂紋不發(fā)生擴展,采動上行裂隙不向上發(fā)育;若處于啟裂跡線外側,則裂紋擴展直至達到自由面,上行裂隙繼續(xù)向上發(fā)育。
1)隨采煤工作面推進,上行裂隙發(fā)育可以分為上行裂隙產生段、發(fā)育延伸段和穩(wěn)定段,在集中拉應力和垂直應力的綜合作用下,上行裂隙發(fā)育頂端巖層裂紋的擴展表現(xiàn)為混合型裂紋。
2)基于裂紋擴展理論,揭示了上行裂隙發(fā)育機理。巖層頂端裂紋在拉應力的作用下沿裂隙帶方向向下擴展直至貫通,導致該巖層完全破斷,上行裂隙向上發(fā)育;反之,若裂紋不發(fā)生擴展,則上行裂隙不再發(fā)育,據(jù)此得到了基于微觀裂紋擴展理論的采動上行裂隙發(fā)育高度的確定方法。
3)建立了裂紋擴展的σ(θ)max理論,提出了裂紋擴展的判據(jù)。若(KⅠ/KⅠc,KⅡ/KⅠc)處于σ(θ)max理論啟裂跡線內,則裂紋不發(fā)生擴展,采動上行裂隙不向上發(fā)育延伸;若該點處于啟裂跡線外側,則裂紋發(fā)生擴展直至達到巖層自由面,上行裂隙繼續(xù)向上發(fā)育。