王創(chuàng)業(yè)，張 琪，李俊鵬，劉 偉

（內蒙古科技大學 礦業(yè)研究院，內蒙古 包頭 014010）

礦床的工作面長度是影響上覆巖層應力和移動變化的主要因素之一，遭到破壞的覆巖應力需要向新的平衡發(fā)展[1-3]，近些年相關專家圍繞工作面長度對覆巖破壞和移動進行了研究，文獻[4-5]研究表明：覆巖內部應力隨著工作面長度的增大而呈指數(shù)性增大，液壓支架的載荷也顯著增加，就目前現(xiàn)代化的采礦工藝來看，布置長距離的工作面有利于減少巷道回采的工作量，提高回采效率，是現(xiàn)代化礦床開采的需要，如今已經(jīng)受到越來越多企業(yè)的重視。本文以內蒙古北山盤陀山鎢礦南部礦段為研究背景，通過FLAC3D有限差分軟件分析不同工作面長度對礦床上部覆巖及地表的移動和破壞的影響[6-9]，旨在有效的保障礦區(qū)安全生產(chǎn)、提高工作面生產(chǎn)效率和保護地表環(huán)境。

1 工程概況及數(shù)值模型

盤陀山鎢礦床位于塔里木板塊與哈薩克斯坦板塊縫合帶南側，內蒙古北山地區(qū)緣內接觸帶及巖體中部，隆起帶中出露地層為中元古界長城系古硐井群，為一套淺海相陸源碎屑巖組合，主要由泥質粉砂巖、粉砂質板巖、花崗砂巖等組成，主要分為北部、中北部、中東部及南部4個鎢礦段，其中南部礦床呈緩傾斜板狀，傾角小于10°，主要產(chǎn)于石英細脈帶中，礦床距地表標高為180～220m，平均厚度為3m，礦體走向長度為55～535m，傾向長度為20～300m，研究內容及模型采用的巖石力學參數(shù)見表1。

利用FLAC3D數(shù)值模擬軟件，本構模型采用庫倫摩爾準則，建立走向長度為500m，傾向長度為250m，厚度為180m的模型。x為傾向方向，y為礦層推進方向，z為巖層方向，在距離頂板不同高度的上覆巖層和地表分別布置測點，監(jiān)測其應力和位移的變化。對模型的底面和四周面的x、y、z面進行約束，上表面為自由面，采用FLAC3D分別對礦層工作面長度為90m、120m、150m、180m工作面推進長度為300m的采空區(qū)進行模擬，進而分析開采不同工作面長度對上覆巖層應力和移動變化的影響。

表1 巖石力學參數(shù)Tab.1 Rock mechanics parameters

2 模擬結果分析

2.1 應力分析

上覆巖層在自重應力的作用下裂斷、破壞和下沉并產(chǎn)生礦壓顯現(xiàn)，巖塊間重新咬合相互擠壓鉸接形成砌體梁并對上部巖體提供支撐和保護的作用，上覆巖層的應力變化程度直接影響著巖層的破壞程度，對工作面、巷道的安全和回采效率具有非常重要的影響，所以隨時掌握巖層間的應力變化進而判斷地表沉陷和巖層移動的影響具有重大意義。

盤陀山鎢礦上覆巖層巖質堅硬，通過公式（1）和公式（2）分別算出垮落帶和導水裂隙帶距工作面頂板的距離。

式中：Hd為礦層頂板上部垮落帶高度，m；Hk為礦層頂板導水裂隙帶高度，m；M為采高，m。

為了更好地了解礦層開采對上覆巖層移動破壞的應力變化，在覆巖的垮落帶、裂隙帶和彎曲下沉帶中沿著工作面推進方向分別布置監(jiān)測點，得到syy方向的應力進行分析，結果如圖1所示。

根據(jù)距離礦層頂板不同高度syy方向水平應力曲線，選取在距離頂板13m、65m、100m的上覆巖層所布置的測點，利用測點1、測點2和測點3、測點4的差值的平均值來表示覆巖內部應力變化的情況，進而得出覆巖內部擾動破壞情況。

圖1（a）為距離工作面頂板13m的syy方向水平應力曲線，從圖中看出應力在工作面推進的過程中發(fā)生應力轉移，由于四周礦層對上覆巖層的支撐保護作用，在工作面開切眼至推進到停采線處應力最大產(chǎn)生明顯的應力集中現(xiàn)象，中部頂板巖層下沉應力變化明顯。

圖1（b）（c）為距離工作面頂板65m和100m的syy方向水平應力曲線，曲線呈鋸齒狀而且波動非常大，在開切眼與工作面停采線處的應力集中現(xiàn)象相比于距離頂板13m處有所增高，當開采深度相同時，隨著工作面長度的變化，回采對上覆巖層的擾動呈現(xiàn)出不同的特征，在距離頂板的距離相同時，隨著工作面長度的增加應力變化的波動程度也在增大，巖層內部受擾動影響強烈。

圖1 距離頂板不同高度應力曲線Fig.1 Stress curve of different heights from the roof

表1 syy應力差值Tab.1 syy stress difference

根據(jù)在距離頂板不同高度測得應力差值，得出測點在離頂板不同距離應力差值的平均值，結果如表1所示。根據(jù)表1的結果顯示：隨著礦層工作面長度的增加，應力差值的均值是增大的，但是應力差值增大的幅度在減小，這說明礦層開采隨著工作面長度的增加對地表的影響在增加，但是增加的幅度是降低的。由此繪制出工作面長度與應力平均差值關系，如圖2所示。

圖2 距離頂板不同高度的應力差值曲線Fig.2 Stress difference curve of different heights from the roof

通過Origin擬合曲線，用公式表達出工作面長度和應力差值平均值的關系，見式（3）～式（5）。

式中：P為各層的平均應力差值，Pa；L為工作面長度，m。

2.2 地表位移分析

地下工作面開采造成了上覆巖層及地表沉陷，引起地表植被、建筑物、耕地的移動和破壞，通過整理出布置在沿走向方向的地表監(jiān)測點所得到的地表下沉值，可繪制出不同工作面長度對應的地表最大下沉值如圖3。

圖3 不同工作面長度對應的最大下沉值Fig.3 Maximum sinking value corresponding to different working face lengths

模擬結果表明：工作面長度是影響地表下沉的重要因素，隨著工作面長度的增加，地表最大下沉量增加的幅度呈現(xiàn)降低的趨勢，當沉陷區(qū)越深時對植被和道路的影響范圍就越廣，地表積水增加了對工作面發(fā)生滲水的可能性。為了更好地了解其關系，現(xiàn)將其擬合，如式（6）所示。

式中：W為地表最大下沉量，mm；L為工作面長度，m。

曲率是地表移動變形的主要指標，對建筑、道路和植被的破壞影響較大，地表拉伸壓縮狀態(tài)對應著曲率的正負值，而拉伸對地表植被的影響較大，這里研究地表的拉伸區(qū)狀態(tài)，所以可以通過曲率變化反映走向斷面上的彎曲程度，進而研究地表的移動變形情況。由于在工作面上方地表沉陷近似盆狀，選取下沉的一側為研究對象，根據(jù)公式得出在走向方向上，不同工作面長度距開切眼不同距離（280m，300m，320m，340m，360m）對應的地表下沉的正曲率變化如表2，繪制出工作面長度與所得曲率平均值的關系如圖4。

研究結果表明：地表曲率變形的動態(tài)發(fā)展對建筑物的影響也是一種動態(tài)的影響過程，隨著工作面長度的增加曲率平均值增加，但是曲率增加的幅度在降低。通過工作面長度可以推斷出地表曲率的變化程度，分析出地表拉伸區(qū)的范圍，進而合理的選擇栽種植被的品種和范圍，增加植被的覆蓋率，減少損失，同時對可能產(chǎn)生裂縫的地區(qū)要提前加強防治，避免地表水滲透影響工作面的生產(chǎn)效率及發(fā)生災害，對于地表植物栽種、建筑物范圍的選擇及保護起到了重要作用。

表2 不同工作面長度走向曲率值對比表Tab.2 Curvature comparison chart corresponding to different working face lengths

圖4 工作面長度-曲率平均值關系曲線Fig.4 Curve of working surface length and mean curvature

3 結論

（1）在相同的地質和開采條件下，隨著工作面長度的增加，應力波動現(xiàn)象明顯，上覆巖層所受的擾動和地表下沉量隨之增大的幅度在減小，通過覆巖應力的變化的大小，預測覆巖移動和破壞程度，對于降低工作面的損傷避免發(fā)生礦山災害并提高生產(chǎn)效率、安全生產(chǎn)起到重要作用。

（2）建立工作面長度與地表下沉量和曲率的關系曲線，通過曲率的變化分析地表沉陷區(qū)的拉伸區(qū)和壓縮區(qū)的范圍，對合理選擇抗拉強度大的植被，提高植被覆蓋率保護生態(tài)環(huán)境，合理的布置建筑物、農田的范圍提供了重要的參考價值，同時對預測地表裂縫并加強防治，避免地表水滲透降低生產(chǎn)效率具有重要的指導意義。

[1] 郭麒麟，喬世范，劉寶?。_采影響下的巖土體移動與變形規(guī)律[J]．采礦與安全學報，2011，28（1）：110-115.GUO Qilin，QIAO Shifan，LIU Baochen．The movement and deformation of soil and rock mass resulted from mining activity[J]．Journal of Mining&Safety Engineering，2011，28（1）：110-115.

[2] 謝廣祥，楊 科．采場圍巖宏觀應力殼演化特征[J]．巖石力學與工程學報，2010，29（增刊 1）：2676-2680.XIE Guangxiang，YANG Ke．Study of macro stress shell evolving characteristics of rock surrounding face[J]．Chinese Journal of Rock Mechanics and Engineering，2010，29（supply1）：2676-2680.

[3] 趙 康，陳斯妮，趙 奎．采動過程金屬礦山空區(qū)覆巖材料細觀破裂特征[J]．地下空間與工程學報，2016，12（4）：920-925．ZHAO Kang，CHEN Shini，ZHAO Kui．Meso-cracking characteristic of overburden strata material over metalmine goaf under mining process[J]．Chinese Journal of Underground Space and Engineering，2016，12（4）：920-925．

[4] 謝廣祥，王 磊．采場圍巖應力殼力學特征的工作面長度效應[J]．煤炭學報，2008，33（12）：1337-1340．XIE Guangxiang，WANG Lei．Effect of long wall length on mechanical characteristics of surrounding rock stress shell in mining face[J]．Journal of China Coal Society，2008，33（12）：1337-1340．

[5] 宋選民，顧鐵風，閆志海．淺埋煤層大采高工作面長度對礦壓顯現(xiàn)的影響規(guī)律研究[J]．巖石力學與工程學報，2007，26（增刊1）：4007-4013．SONG Xuanmin，GU Tiefeng，YAN Zhihai．Effects of increasing working faces length on underground pressure behaviors of mining super-high faces under shallow coal seam[J]．Chinese Journal of Rock Mechanics and Engineering，2007，26（supply1）：4007-4013．

[6] 劉書賢，魏曉剛，王 偉．深部采動覆巖破斷的力學模型及沉陷致災分析[J]．中國安全科學學報，2013，23（12）：71-76.LIU Shuxian，WEI Xiangang，WANG Wei．Mechanical model of overlying strata rupture caused by deep coal mining and analysis of strata subsidence disaster[J]．China Safety Science Journal，2013，23（12）：71-76.

[7] 許國勝，李德海，侯得峰，等．厚松散層下開采地表動態(tài)移動變形規(guī)律實測及預測研究[J]．巖石力學，2016，37（7）：2056-2062．XU Guosheng，LI Dehai，HOU Defeng，et al．Measurement and prediction of the transient surface movement and deformation induced by mining under thick alluvium[J]．Rock and Soil Mechanics，2016，37（7）：2056-2062．

[8] 郭廣禮，馬悅漢，馬占國．煤礦開采沉陷有效控制的新途徑[J]．中國礦業(yè)大學學報，2004，33（2）：151-154．GUO Guangli，MA Yuehan，MA Zhanguo．A new method for ground subsidence control in coal mining[J]．Journal of China University of Mining&Technology，2004，33（2）：151-154．

[9] 梁 冰，石占山，孫維吉．構造應力作用下覆巖卸壓演化規(guī)律試驗研究[J]．安全與環(huán)境學報，2014，14（1）：33-37.LIANG Bing，SHI Zanshan，SHUN Weiji．Experiment study of relief evolution of the overlying strata under the effect of the tectonic stress[J]．Journal of Safety and Environment，2014，14（1）：33-37.