趙世振,翟 迪
(西安科技大學 能源學院,陜西 西安 710054)
條帶開采法是一種非充分的部分開采方法[1-3],條帶開采的沉陷控制效果較好,但采出率低、生產效率低,顯然不應成為主要的沉陷控制方法。充填開采技術不僅可以解放“三下”壓煤,提高采煤率,還可以有效保護環境,但其實際應用效果、可靠性受到多種因素影響,與預期還存在較大的差異,全采全充布置模式因覆巖整體性擾動,還存在沉陷失控的風險,并且充填成本較高,材料來源難以保障,充填可靠性還有待提高,這些問題給充填開采方法的推廣應用造成了障礙[4-6]。
協調開采是一種寬泛的巖層移動時空控制理念[7-9],是指通過采區或工作面優化布置和開采過程調整來控制巖層移動的程度、均勻性和影響過程的技術手段。朱曉俊等[10]對帶狀充填開采巖層移動機理進行了研究,并得出了充填體與煤柱復合支撐體的強度規律;馬立強等[11]研究了淺埋煤層保水開采技術工藝并在適宜地質條件下得到了成功應用,戴華陽等[8]通過應用設計,提出了村莊下“采-充-留”協調開采方案設計,并驗證了該方法對地表移動變形的良好控制效果。結合已有研究,提出了一種網狀帶式充填開采技術,該技術結合條帶開采及充填開采的優點,在條帶開采的同時,對上一個采空條帶即時進行充填;再以充填體為保護條帶,對條帶煤柱采用部分充填開采,最后充填大巷,形成一個整體的網狀結構充填體,充填體支撐覆巖達到減沉目的,最終實現煤炭全采。
網狀帶式充填開采技術是一種第Ⅰ階段先部分開采,再即時充填,第Ⅱ階段開采條帶保護煤柱,再部分充填的全采方法。依據開采沉陷基本原理,在非充分開采條件下,地表的移動變形可以得到較好的控制;而極不充分開采[12]能更大幅度地減小地表移動變形,從而達到控制地表沉陷與保護地表建(構)筑物的目的。一般來說,在采礦地質條件一定的情況下,工作面尺寸越小,采動程度越小,對地表的影響程度越弱。結合協調開采的基本思想,在煤層中構建多個非充分采動工作面,是網狀帶式充填開采的必要條件,也是其基本原理。
第Ⅰ階段為即時充填的條帶開采。根據具體地質采礦條件,遵循開采寬度小于頂板破斷距的原則[13]設計條帶開采參數,在進行下一個條帶開采的同時,對上一個開采條帶進行充填如圖1(a)所示。第Ⅱ階段為置換條帶煤柱的部分充填開采。在工作面完成預先設定的條帶開采,假設完成即時充填措施后,巷道圍巖受開采擾動較小,巷道可再次使用。在工作面形成回采條帶煤柱的生產系統后,對保留條帶煤柱進行開采,每開采一個條帶,隨即進行充填,但只充填回采巷道的一部分,當整個工作面充填完成后,充填回風與運輸順槽,整個采區會形成一個整體連接的網狀結構充填體以減緩上覆巖層的垮落和地表沉陷,如圖1(b)所示。

圖1 開采流程示意
根據我國開采實踐經驗表明,考慮到現有的連采機連采技術的熟練使用,為提高采煤效率支巷寬度應該盡可能的大于5 m。對于充填條帶開采而言,根據保留煤柱最小強度規定[14],寬高比應大于2,因此保留條帶寬度應該大于4 m。
某礦采區走向長約3.3~4 km,傾斜寬約1.9 km,面積約6.8 km2。煤層開采上限標高+380 m,開采下限標高+140 m。本區主采煤層平均埋深308 m,煤層厚度2 m,為近水平煤層。
2.2.1 計算模型的建立
考慮到該地區的地表變形移動角和邊界角等因素,為了保證地表的充分采動,模擬開采長度定為500 m,約1.5倍采深,向開采邊界以外延伸600 m,約2倍采深,建立對稱模型;模型高度從地表到煤層底板以下30 m。煤巖層頂底板組成及物理力學參數見表1。

表1 煤巖層頂底板及物理力學參數
2.2.2 地表移動觀測線的布置
沿采煤工作面走向和傾向布設兩條觀測線,如圖2所示。記錄下沉值,走向測線從切眼前50 m開始每隔10 m設一測點,考慮到模型邊界效應,測點布置到距采空區右邊界50 m處,共60個頂板測點,編號從“2”到“62”。傾向測線從下到上多出工作面寬度各50 m,共53個,編號從“63”到“115”。每個開采循環,記錄所有測點數據,每開挖10步記一次,由此可得出各測點在整個開采過程中的動態下沉值,繪制出各點動態變形曲線。

圖2 地表移動觀測線布置
充填體力學參數采用密度為1 370 kg/m3、體積模量為2.83 GPa、剪切模量為1.31 GPa、內聚力為3.8 MPa、抗拉強度0.9 MPa及內摩擦角為20°的高強度膏體充填材料。模型如圖3所示。

圖3 FLAC3D模型
將采用傳統條帶開采采煤方法與采用帶狀充填法開采時地表變形的豎向位移提取出來,并繪制其等值線下沉圖。從圖4~圖7可以看出無論采用何種充填工藝進行充填開采時地表都會出現中間下沉盆地,這與開采沉陷理論相符。當采用采10留5條帶開采時地表變形明顯,盆地中心下沉值最大值為120 mm,如圖4所示;當采用這種帶式充填采煤方法進行開采時,地表變形比采用條帶開采時明顯有所減小,此時盆地中心的最大下沉值為55 mm,比圖4中最大下沉值減小了一半以上,如圖5所示;當采用采10留10的條帶開采采煤方法進行開采時,相比于采用采10留5條帶開采的地表沉陷三維圖,此時的地表沉陷下沉值縮小了約60%,如圖6所示;當采用采10留10帶式充填開采方法進行開采,盆地中心的最大下沉值與采10留5網狀充填開采方法相比相差15 mm,但回采率低了17%,如圖7所示。

圖4 采10留5條帶開采

圖5 采10留5條帶充填協調開采

圖6 采10留10條帶開采

圖7 采10留10條帶充填協調開采
不同方法開采時的最終下沉值如圖8所示(TD-條帶開采;DC-帶式充填),對應地表下沉系數見表2。地表下沉系數計算式為

表2 不同開采方法的地表下沉最大值與下沉系數

圖8 地表下沉值
qc=wmax/Mcosα
(1)
式中,qc為下沉系數;wmax為地表最大下沉值;M為煤層厚度;α為煤層傾角。
由表2中對比可以看出對于條帶開采來說采留比與下沉系數之間有著一定的聯系,當保留煤柱為10 m時下沉系數qc僅為保留煤柱為5 m時的42%左右,所以保留煤柱是值得關注并研究的問題,既能保證采出率達到效益最大化而又不影響地表建構筑物的安全。而對于這種帶式充填開采,由表3可知地表下沉主要發生在條采充填的第Ⅰ階段,第Ⅱ階段置換保留煤柱的開采方式對地表的影響微乎其微,地表最大下沉值變化對于兩種采留寬方式的不同,差別僅在10~20 mm左右,但回采率相差17%,為進一步研究這種帶式充填開采地表下沉系數與采留比之間的關系,在不改變地質條件的前提下,以上述高強度充填體為原型通過改變其采留寬進行數值模擬,得到了開采地表下沉系數和采留比的數據,見表3。
從表3可看出,隨著地表最大下沉值與下沉系數qc的增加,隨開采條帶寬度與保留煤柱寬度的增加而顯著增大。當采留寬分別為40時,第Ⅱ階段的地表下沉值是第Ⅰ階段的2倍,因此可以推測了這種充填方式充填完成后空頂破壞的臨界采留寬為40:40。

表3 不同采留比條件下的地表下沉系數
為進一步探究該網狀帶式充填采煤方法的影響因素,以FLAC3D建模采20留10為背景,以上述充填體為原型,通過改變充填體彈性模量或改變采高的方式進行數值模擬,得到數據見表4、表5。

表4 不同彈性模量充填體探究

表5 不同采高對下沉值的影響
從表4中可以看出,下沉值W隨充填體的強度增加而減小,當強度增大到一定程度,下沉值不會有所改變,而此時的下沉值近似為未充填區域的平均下沉值。
而表5中對于不同采高的變化,下沉值的增長呈現線性增長,對比條帶開采的臨界采高值[16-18],這種網狀帶式充填開采,臨界采高相對更大,也適用于大采高。
基于FLAC3D建立數值模型,得出各因素對最大下沉值及下沉系數關聯度排序,各因素對地表最大下沉值影響程度排序為:采留寬>采出率>采高>充填體強度。
在半無限開采情況下,除下沉外,主要的移動和變形值均發生在x=-r~+r的范圍之內。故稱r為主要影響半徑如圖10所示。

圖10 參數r與tanβ的幾何意義
將x=±r的地表點與煤壁(煤層水平上x=0的點)相連,其連線與水平線之間所夾的銳角β稱為主要影響角,其正切tanβ稱為主要影響角正切。其公式為
(2)
規程中列出了我國各大礦區的實測tanβ值,一般tanβ在1.2~2.6。根據實測資料求取tanβ的方法是在傾向充分采動、走向半無限開采的走向實測下沉曲線W(x)上,分別確定下沉值為0.16W0和0.84W0的點,它們和下沉值為0.5W0的點之平距均應為0.4r,從中求出r(若兩個平距求得的r值稍有不同,可取均值);然后再根據采深計算出tanβ。按照覆巖性質區分的概率積分參數的經驗值,此處取tanβ為1.5。
拐點偏移距就是拐點偏離采空區邊界的距離。實測曲線上拐點的位置可按其性質確定。具體為,拐點處的下沉值為0.5W0,拐點處的傾斜和水平移動值最大,距拐點左右0.4r處的曲率和水平變形值最大。
由于實測誤差,求定的拐點偏距常不重合。可根據實際情況選定最適當的拐點位置,拐點偏距一般為(0.05~0.3)H之間。
通過地表移動預計表明地表建(構)筑物受開采影響輕微,且均在Ⅰ級以內,下沉量較小,能有效地控制地表建(構)筑物的損害。
(1)針對條帶開采采出率較低、全采全充成本較高的問題,提出了一種條帶開采與充填開采相結合的網狀帶式充填開采方法,該方法實現了煤層全采,并闡述了其基本原理及技術步驟。
(2)通過數值模擬分析條帶開采與帶式充填開采控制地表變形沉陷,相比較帶式充填開采,地表下沉值減少50%以上,并給出了充填區域的極限垮距為40:40,當采留寬為40:40時,第Ⅱ階段置換保留煤柱地表最大下沉值會翻倍。
(3)地表下沉值隨充填體的強度增加而減小,當強度增大到一定程度,下沉值不會有所改變,而對于不同采高的變化,下沉值的增長呈線性增長,對比條帶開采的臨界采高值,網狀帶式充填開采,臨界采高相對更大,因此這種方法同樣適用于大采高。
(4)通過模擬,得出各個因素對地表下沉影響程度的排列順序:采留寬>采出率>采高>充填體強度。并通過對臨界采留寬之內的不同采留寬布置對地表下沉的預計,結果顯示地表下沉對建(構)筑物的影響都在Ⅰ級影響之內。