石膏礦采空區穩定性研究進展

肖 超, 鄭懷昌, 王 棟, 李 明, 張曉君, 劉志河

(山東理工大學 資源與環境工程學院, 山東 淄博255091)

21世紀初，國內先后發生十余起石膏礦采空區大冒落，給人們人身安全帶來了巨大危害，對采空區進行穩定評價是安全管理采空區的必要條件[1].我國現有石膏礦山300余座，采空區1億多立方米，普遍采用房柱法開采，再加上石膏礦經濟效益低，很多礦山不能及時治理采空區，產生了大批未治理的采空區.隨著時間的延續，采空區面積不斷變大，給礦山生產帶來了極大安全隱患[2].目前缺乏可靠、普適、系統的穩定性評價方法，這也是采空區安全管理的重大障礙，及時進行石膏礦山穩定性分析具有重要意義.

采空區穩定性是一個包含采空區本身和周圍地質構造及所覆巖層巖性性質，同時又受到開采和爆破造成的擾動﹑地下水﹑時間引起的流變等復雜多因素影響的大尺度平衡問題.

1 經典理論的采空區穩定性研究

早期采空區理論研究大多基于支撐采空區的礦柱和頂板自持能力，運用巖體靜力學進行分析.M.M.普羅托吉亞科夫﹑哈元﹑舒爾茨﹑施托克﹑李特維尼申分別提出拱形冒落論﹑壓力拱假說﹑懸臂梁(板)冒落論﹑冒落巖塊碎脹充填論及開采沉陷的隨機介質理論.

我國學者發展成為了概率積分法冒落巖塊鉸結論以及砌體梁平衡說及采場薄板礦壓理論，王學濱[3]根據漸進破壞的礦柱，整理出礦柱剪切失穩判據的解析式.近期學者在采空區穩定性分析中引入突變理論，王連國[4]由分支曲線方程及定態曲面方程，利用突變理論建立了礦柱失穩的尖點突變模型，并基于此模型又分析預測了礦柱失穩機理；吳啟紅[5]建立了基于突變理論的強度折減法來分析多層采空區頂板穩定性.

鄭懷昌[6]引用界殼理論探討采空區穩定性支撐結構的演變對采空區穩定性的影響；趙忠明[7]在采空區遭受庫水的浸泡﹑沖掏和溶蝕條件下，利用相似模擬實驗研發治理；基于尖點突變理論，分析不同形式下的采空區的穩定性，頂板冒頂規律的準確預測在控制和消除相關問題是非常重要的，模糊邏輯被應用到煤礦預測冒頂規律[8].張曉君[9]認為礦柱支撐面積與礦柱面積之比和采深等因素對采空區的穩定性影響很大.陳建平[10]通過對采空區的現場調查礦巖取樣及圍巖節理裂隙的綜合分析，對現有空區穩定性進行了論證.鄧紅衛[11]通過分析采空區穩定性各影響因素，建立基于不確定型層次分析法(UAHP) 的采空區穩定性模糊評判方法.鄭茂興[12]運用綜合評判法，依據綜合得分對采空區穩定性進行合理劃分.

2 采空區穩定性靜態數值分析

理論分析給予石膏礦采空區穩定性分析定性指導；而數值模擬能定量﹑快速分析采空區穩定性，且能解決繁瑣的材料本構聯系，所以有廣泛的應用價值.根據不同特征的采空區評價采空區的穩定性，現有可利用數值模擬軟件有ANSYS、FLAC、RFPA、3DEC等.

李明等[13]采用FLAC3D軟件模擬分析了石膏礦山采空區的穩定性狀態,分析對比了長期強度對采空區穩定性的影響；楊宇江等[14]利用FLAC3D軟件和Mohr-Coulomb屈服準則，研究了在附加載荷的作用下，不同厚度的空區頂板的穩定性；張曉君[15]運用巖石破裂過程分析軟件RFPA2D,對采空區的破壞過程以及采空區巖性對采空區頂板大面積冒落破壞過程進行了相應數值模擬研究;張曉君[16]利用ANSYS軟件，數值模擬了采空區在不同深度不同傾角時的應力應變情況；運用FLAC3D軟件建立粘彈塑性本構模型(cvisc模型)，對石膏礦房礦柱的長期穩定性進行了數值模擬研究[17]；李長洪[18]利用3DEC軟件數值模擬分析了采空區形成后圍巖變形和應力分布情況.

3 石膏礦采空區穩定性動態研究

隨著時間推移，受巖石流變﹑爆破﹑地下水多因素影響，采空區巖石強度及礦柱-頂板-底板系統抵抗冒落能力都已大大折減，經典常規的靜態分析石膏礦采空區已不能滿足現實情況.

決定采空區穩定性的是頂板、礦柱、底板及支護共同組成的支撐系統的穩定性.王金安[19]認為漸進破壞的頂板和礦柱造成采空區堅硬頂板瞬間失穩及塌陷；運用力學分析采空區的模型中，J. A. Wang等[20]闡述了頂板力學狀態與礦柱總面積之間的聯系.分析了風化及剝落造成礦柱承載面積減小對采空區頂板穩定性狀態的影響，然而沒有總結怎樣將礦柱自身流變與采空區頂板的穩定性時間聯系起來.

在采空區穩定性長期強度方面，劉沐宇[21]采取室內純扭轉流變實驗的手段，探究了硬石膏的流變特性；Auvray等[22]分析了相對濕度造成礦柱強度的衰減；Biswas[23]針對不同時間的礦柱作出了應力測量,建立了風化厚度與礦柱時間之間的關系；Castellanza等[24]由水影響下石膏礦柱單軸壓縮實驗數據建立相應風化模型來預測礦柱穩定時間.

巖石流變導致長期強度下降，影響采空區長期穩定性，王金安[25]基于自建的礦柱-頂板體系流變力學模型，又依據推出的位移控制方程分析討論了不同階段下采空區頂板破壞.構成采空區支撐系統的礦柱和頂板的巖體受到多次爆破擾動作用造成了損傷，長久的累計損傷更近一步弱化了支撐系統[26-27]；而且，采空區應力場受到巖石的流變作用發生了變化[28]，加快了采空區失穩的演化歷程.

爆破對支撐采空區穩定性的礦柱及頂板上的巖石產生了損傷，趙延林等[29]提出采空區重疊頂板安全系數的概念，建立其豎向位移序列與折減系數的尖點突變模型；易平[30]運用GeoFBARV4.0數值模擬軟件，對采空區圍巖破壞發展的動態數值過程進行模擬分析.地震、爆破等動力作為影響采空區穩定性的主要因素之一，不但對采空區圍巖產生損傷和擾動，甚至引起采空區破壞和塌陷.楊金林等[31]通過改變擾動應力波峰值的大小，考查動力擾動強度的變化對采空區穩定性的影響.鄭懷昌[1]認為多因素耦合動態演化導致了采空區長期穩定性研究的復雜性.

4 存在的問題

石膏礦采空區應該看成一個系統，深入探討石膏礦采空區礦柱-頂板-底板系統多因素耦合下動態演變機制，為采空區長期穩定性探究新研究方法，可使石膏礦采空區穩定性分析研究上升到一個新的臺階.問題如下：

(1)石膏礦采空區分析周期長，加上石膏礦利潤低，不能單純追求效果，還要考慮經濟因素.

(2)對石膏礦采空區無論靜態數值分析或動態分析，大多基于單因素分析.如何運用系統觀點，基于數學模型﹑數值模擬等定量科學手法，缺乏石膏礦采空區巖石流變﹑爆破﹑地下水﹑長期強度等多因素耦合模型及動態演變機制.

(3)現階段的采空區穩定性理論方面研究距離實際有不少差距，指導實際工程作用不大，并且當前多集中在對礦柱、頂板的穩定性研究，缺少對頂板-礦柱-底板系統以及采空區整體穩定性的研究.

(4)當前大范圍利用數值模擬技術分析采空區穩定性，對采空區安全管理具有重要的指導意義.但是，分析中所運用的巖石力學參數和采空區應力分布多是根據某一時間的靜態值.這種分析方法的準確程度有待進一步提高.需要加深石膏礦采空區穩定性動態分析，建立基于常應力的石膏礦采空區動態破壞準則.

5 石膏礦采空區穩定性研究展望

以下為石膏礦采空區穩定研究的一些建議及新思路：

(1)采空區的長期穩定性研究是一個非常復雜的問題.根據巖石力學參數與靜態地壓分析得到的采空區穩定性與經受幾十年流變和開采爆破影響的礦柱-頂板-底板體系實際穩定性有較大差異，所以準確分析采空區長期穩定性需要將流變和爆破擾動對巖石力學性質的作用等量化到數值分析中.

(2)數值模擬可以比較系統的研究在動應力作用下采空區圍巖的動應力場分布規律及動靜應力場的疊加作用機理，可以對FLAC3D自定義本構模型進行二次開發，將動態數值模擬編入FLAC3D模擬中，也可運用虛擬現實等軟件模擬宏觀的應力場變化與微觀的巖石損傷累積關系.

(3)國內外分別從巖石流變、爆破擾動等單因素方面對采空區及地下空間工程的長期穩定性進行了相關研究，也有相關學者通過流變擾動效應試驗開展了巖石處于強度鄰域時的雙因素耦和作用的影響研究. 將石膏礦采空區巖石流變﹑爆破﹑地下水﹑長期強度等多因素進行耦合探究，可以為采空區這類大尺度地下空間工程長期穩定性研究開辟新的途徑.

(4)可以采用抗塌性衡量采空區“礦柱-頂板-底板”體系.作用于抗塌性體系的靜應力場與爆破沖擊能量場的耦合“力場”，隨開采推進，對于特定研究點而言也是處于衰減狀態.運用能量觀點，可以建立適應常應力狀態的動態破壞準則.

(5)可以把石膏礦采空區看成一個系統，運用科學數值方法建立量化的以時間﹑應力狀態﹑損傷特征﹑巖石含水率﹑擾動次數為變量和彈模﹑長期強度為結果的時域曲面和石膏礦采空區多因素耦合模型，建立動態演變機制.

[1]鄭懷昌，李明，劉志河，等.石膏礦采空區穩定性動態演化規律研究[J].化工礦物與加工，2012(12)：28-31.

[2]鄭懷昌，宋存義，胡龍，等.采空區頂板大面積冒落誘發沖擊氣浪模擬[J].北京科技大學學報，2010,32(3)：277-281.

[3]王學濱.屈服礦柱漸進破壞及應力分布數值模擬[J].中國地質災害與防治報,2006,17(2):50-56.

[4]王連國，繆協興.基于尖點突變模型的礦柱失穩機理研究[J].采礦與安全工程學報，2006,23(2):137-140.

[5]吳啟紅.礦山復雜多層采空區穩定性綜合分析及安全治理研究[D].長沙：中南大學，2010.

[6]鄭懷昌，李明.界殼理論在采空區失穩判定與危害控制研究中的應用探討[J].黃金，2005，26(12)：19-22.

[7]趙忠明.小浪底庫區采空區穩定性研究[D].焦作:河南理工大學,2011:1-90.

[8]Razani M , Yazdani-Chamzini A , Yakhchali S H.A novel fuzzy inference system for predicting roof fall rate in underground coal mines[J].Safety Science,2013,55:26-33.

[9]張曉君.影響采空區穩定性的因素敏感性分析[J].礦業研究與開發,2006(2):14-16.

[10]陳建平.采空區穩定性論證[J].采礦技術，2012(3)：52-54.

[11]鄧紅衛.基于UAHP 的采空區穩定性模糊綜合評判[J].中國安全科學學報， 2012 (4)：24-28.

[12]鄭茂興.石膏礦采空區穩定性分類的綜合評判法研究[J].化工礦物與加工，2012(5)：24-27.

[13]李明，鄭懷昌,劉志河，等.石膏長期強度對采空區穩定性的影響分析[J].化工礦物與加工，2010(2)：21-23.

[14]楊宇江，高晗，李元輝.附加載荷作用下空區頂板穩定性分析，東北大學學報：自然科學版，2011，32(9)：1332-1335.

[15]張曉君.采空區頂板大面積冒落的數值模擬[J].化工礦物與加工，2007，27(1)：17-18

[16]張曉君.礦柱及圍巖對采空區破壞影響的數值模擬研究[J]采礦與安全工程學報，2006(3)：123-126.

[17]楊布華.我國非煤礦山地下采空區穩定性研究進展[J].礦產與地質，2008(10) :473-478.

[18]李長宏，程巖.基于3DEC 的采空區穩定性分析[J].金屬礦山，2012 (4) :144-147.

[19]王金安，尚新春，劉紅等.采空區堅硬頂板破斷機理與災變塌陷研究[J].煤炭學報,2008,33(8):850-855.

[20]Wang J A，Shang X C，Ma H T. Investigation of catastrophic ground collapse in Xingtai gypsum mines in China[J]. InternationalJournal of Rock Mechanics and Mining Sciences，2008，45(8)：1480-1499.

[21]劉沐宇,徐長佑.硬石膏的流變特性及其長期強度的確定[J].中國礦業，2000,9(2): 53-55.

[22]Castellanza R, Gerolymatou E, Nova R. An attempt to predict the failure time of abandoned mine pillars [J]. Rock Mechanics and Rock Engineering,2008, 41( 3 ): 377-401.

[23]Auvray C，Homand F, Hoxha D. The influence of relative humidity on the rate of convergence in an underground gypsum mine[J]. International Journal of Rock Mechanics and MinngSciences, 2008,45(8) : 1454-1463.

[24]Sterpi D, Gioda G. Visco-plastic behaviour around advancing tunnels in squeezing rock[J]. Rock Mechanics and Rock Engineering, 2009, 42(2):319-339.

[25]王金安,李大鐘,馬海濤.采空區礦柱——頂板體系流變力學模型研究[J].巖石力學與工程學報,2010,29(3)：577-582.

[26]Nazarov L A, Nazarova L A.Estimate of the interchamber pillar stability based on the damage accumulation criterion[J]. Journal of Mining Science,2007，43(6)：575-584.

[27]Ma G W, Hao H，Wang F. Simulations of explosion-induced damage to underground rock chambers[J]. Journal of Rock Mechanics and Geotechnical Engineering,2011, 3 (1): 19-29.

[28]王波,高延法,王軍.流變擾動效應引起圍巖應力場演變規律分析[J].煤炭學報,2010，35(9)：1446-1450.

[29]趙延林，吳啟紅，王衛軍，等.基于突變理論的采空區重疊頂板穩定性強度折減法及應用[J].巖石力學與工程學報，2010，29(7)：1424-1434.

[30]易平.地下采空區圍巖破壞發展的有限元數值分析運用[J].有色設計與研究，2010，31(5)：5-7.

[31]楊金林，李夕兵，周子龍，等.動力擾動對采空區穩定性影響的離散元分析[J].科技導報,2010,28(4)：48-52.