樟村坪國家礦山公園采空區穩定性研究

郭 峰, 王水華, 蔣新洪, 渠 婧

(1.湖北省地質局 水文地質工程地質大隊,湖北 宜昌 443000; 2.湖北省地質局 第八地質大隊,湖北 襄陽 441002)

樟村坪國家礦山公園采空區穩定性研究

郭 峰1, 王水華1, 蔣新洪2, 渠 婧2

(1.湖北省地質局 水文地質工程地質大隊,湖北 宜昌 443000; 2.湖北省地質局 第八地質大隊,湖北 襄陽 441002)

樟村坪國家礦山公園為全國首家磷礦國家礦山公園,園區已形成大面積采空區,采空區冒落形成多處地質災害。通過研究分析區內采空區及地質災害發育分布特征與形成條件,運用FLAC/FLAC3D數值模擬方法對采空區穩定性進行研究,并有針對性地提出防治措施建議,對樟村坪國家礦山公園建設及類似礦山防災減災具有現實意義。

采空區;穩定性;礦山公園

樟村坪國家礦山公園所在的湖北昌達化工有限責任公司樟村坪磷礦位于樟村坪礦區Ⅲ礦段,毗鄰樟村坪集鎮,礦區面積2.107 9 km2,是宜昌磷礦開發最早的礦區之一。2005年以來,因面臨資源枯竭,主要以回采為主,同時對原采空區進行部分礦柱置換,截至2011年12月底,累計出礦5 869.6 kt。2013年,湖北昌達化工有限責任公司宜昌樟村坪磷礦獲批國家礦山公園。

湖北昌達化工有限責任公司宜昌樟村坪磷礦區已形成大面積采空區,面積約129.7×104m2[1],占礦山總面積62%。采空區冒落導致礦區內地質災害活動日漸頻繁,因該礦山毗鄰樟村坪集鎮,采空冒落型地質災害不僅危及礦山安全生產,還影響國家礦山公園旅游產業規模的提升,對樟村坪集鎮的生態、環境退化也有加劇趨勢,危害性與潛在威脅巨大。因此,通過研究分析區內采空區及地質災害發育分布特征與形成條件,提出有針對性的防治措施建議,對國家礦山公園開發利用及同類礦山防災減災具有現實意義。

1 區域地質背景及開采歷史

1.1 地形地貌

區內地貌類型屬黃陵背斜東翼單斜構造侵蝕溶蝕中低山區。區內地勢北西向南東逐步降低,制高點老君坐高程1 252 m,最低點樟村坪集鎮高程894 m,最大高差258 m。本區為黃柏河東支的源頭——黃馬河發源地之一。西側坡陡東側斜坡稍緩,坡體中、上部震旦系上統燈影組白云巖出露區多為高陡絕壁、奇峰林立,坡體中下部陡山沱組炭質微晶白云巖、炭質頁巖出露區多構成緩坡。具體地形參見圖1。

圖1 樟村坪礦段地形立體圖Fig.1 Topographic stereogram of Zhangcunping ore section

1.2 地層巖性

礦區出露地層由老至新為:震旦系陡山沱組(Z1d)巖性為微晶白云巖、含碳質微晶白云巖,夾薄層泥質白云巖、含磷礦層,礦體厚度平均3～4 m。震旦系燈影組(Z2dn)巖性為白云巖、硅質白云巖、白云質磷塊巖、泥質白云巖及角礫狀白云巖組成。

圖2 樟村坪礦山公園采空區及地質災害分布圖Fig.2 Distribution map of goafs and geological disasters in Zhangcunping mine park1.礦山公園范圍;2.礦柱;3.礦層露頭線;4.地下井巷及編號;5.采空區礦柱、頂板、底鼓集中區;6.采空區界線;7.地面塌陷及編號;8.危巖及編號;9.地裂縫及編號。

1.3 地質構造

礦區為一平緩的單斜構造,走向近南北,傾向東,傾角5°～15°,平均8°,褶皺不發育。本區段主要斷層有8條,其中落差>20 m的有2條,按走向大體可分為北西向及北東向兩組,以北西向斷層為主。對礦區有影響的主要為樟村坪斷裂(F1),為正斷層,通過礦區南部邊緣,構成礦段南界。斷層走向110°～290°,延伸1 000 m,傾向20°,傾角52°～62°;斷層下盤為崆嶺群變質巖,上盤為陡山沱組;該斷層破碎帶一般為8～32 m,向深部變寬,斷裂帶發育碎裂巖、破碎巖塊,并見混合花崗巖的角礫等。

1.4 礦床地質特征

1.5 水文地質條件

區內地下水類型主要為碳酸鹽巖裂隙溶隙水。磷礦層賦存于當地侵蝕基準面以上,由于礦段四周被溝溪切割較深,地表水和礦段地下水無水力聯系。地下水主要接受大氣降水補給,丁家河為區內主要供水水源地。

1.6 礦山開采歷史及現狀

1966—1968年原湖北省地質局第九地質大隊即對樟村坪礦區進行了普查和勘探,上世紀70年代有多家鄉鎮礦山開采;在80年代,由原化工部與湖北省政府共同投資6 648萬元興建湖北省宜昌樟村坪磷礦,2003年企業改制轉為民營,更名為湖北昌達化工有限責任公司樟村坪磷礦。

樟村坪磷礦采用地下開采方式,為房柱法和全面法回采,整個礦段布設有1 028、1 010、 987和904四個中段及1 050副中段。該礦山已形成大面積采空區,總面積約129.7×104m2,此外鄰近磷礦企業在礦體北側和東側附近也形成約50×104m2的采空區面積。這些采空區基本都未加處理,依賴預留的礦柱來維持采空區的穩定,但采空區內預留的保安礦柱、點柱和條形礦柱曾遭受嚴重的盜采,使采空區礦柱系統的支撐能力削弱。大面積采空區分布對井下工作人員、設備和規劃中的國家礦山公園安全構成了嚴重威脅。

2 采空區及地質災害分布

2.1 采空區分布及變形特征

湖北昌達化工有限責任公司樟村坪磷礦已分別對1 028 m、1 010 m、987 m、984 m、982 m、981 m、980 m等標高中段進行開采,其中1 028中段、1 010中段、987中段已經閉坑,已在礦山西部形成一個約71×104m2的采空區,部分已塌陷;東部及東南部形成一個長條狀約41.3×104m2的采空區,地表已開裂,局部已塌陷。全礦主巷道已深達1 600 m,采空區總面積達112.3×104m2。2008年以后,礦山生產以規?；夭蔀橹?同時對原采空區進行部分礦柱置換,生產中段主要為984中段、983中段、982中段、981中段、980中段。西南角原設計預留的陡巖山坡及民房保安礦柱亦已部分開采。截至2011年12月底,礦山已形成總面積約129.7×104m2采空區(圖2),全礦備采礦量已形成,已至生產尾期。

井下采空區調查顯示,井下的頂板穩定狀況、礦柱完整情況整體上呈現出從北到南趨好的現象。在1 028中段280～284、1 010中段100～105、1 010中段106～109與987中段875～877之間的雙包至T3#古裂縫區域頂板局部已冒落,冒落區邊緣部分尚未接頂,冒落體塊度大小不一,冒落區附近礦柱片落嚴重,如照片1所示。109、877以北(雙包一帶)及南側984采場1#、2#溜井附近分別出現礦柱壓裂、頂板顯著下沉,底鼓現象,礦柱受力變形破壞以受壓劈裂方式為主(照片2),破壞礦柱表層出現縱向裂縫并側鼓、剝離,此外部分礦柱變形破壞呈剪切方式(照片3),絕大多數破壞的礦柱并未完全壓潰或剪斷,仍有部分支撐能力,少量破壞礦柱由于留設尺寸偏小,在重力作用下已完全失去承載能力;空區內出現底鼓的地點較少,底鼓開裂延伸長度也僅為2 m左右(照片4);頂板下沉現象較普遍,空區內很多信號柱被壓彎、折斷,所設置的木滑尺也紀錄了一定的變形量(照片5),表明頂板圍巖由于采場開挖的卸荷回彈或受力產生了彎曲變形。

照片2 礦柱受壓劈裂破壞Photo 2 Splitting failure of pillar under compression

以上變形現象表明該區域采空區有顯著的地壓活動。從282中部、104、874～870往南胡家屋場一帶,采空區的整體穩定性較好,空區內礦柱分布不均勻,局部最大跨度達80～100 m,礦柱尺寸普遍為3 m×3.5 m、4 m×4 m、4 m×5 m點柱,連續礦壁很少,也有1～2 m寬的小礦柱,礦柱間最大間距一般為11～14.5 m,雖有個別礦柱破壞現象,但所占比例較小。在1 028中段預留的保安礦柱目前已被部分回采,邊緣的連續礦柱厚度顯著減小,回采區段多為點柱,但穩定性較好,礦柱完好。

照片3 礦柱受壓剪切破壞Photo 3 Shear failure of pillar under compression

照片4 采場底板底鼓現象Photo 4 Floor heave of stope floor

照片5 采空區頂板下沉Photo 5 Roof subsidence in goafs

2.2 地質災害分布特征

長期采礦使樟村坪國家礦山公園區內形成129.7×104m2的采空區,引起地表山體變形開裂,形成地面塌陷及崩塌(危巖、滾石),破壞森林植被和耕地,嚴重威脅采空區地表及樟村坪集鎮人民生命財產安全。區內共分布各類地質災害7處[2],其中地面塌陷(地裂縫)5處、危巖2處(表1)。

區內地質災害形成的直接原因為采礦形成大面積的采空區,主要動力來源于礦層頂板重力及雨水入滲。受毗鄰樟村坪斷裂影響,礦區巖層中高角度構造裂隙發育,礦層頂板較薄,頂板由薄層狀的脆性巖體構成,頂板的巖層被切割呈塊狀,完整性較差。采礦后采空區上部山體應力重新調整及雨水的入滲加載,使頂板逐漸破壞從而產生塌落,因頂板薄而在地表形成地裂縫、塌陷坑及危巖體,地表變形區大致沿著淺層采空區或淺層井巷發育。如雙包地裂縫(危巖)、T3#古裂縫及羅家屋場危巖即位于1 010中段106～109與987中段875～877之間的早期采空區上部山體,據井下采空區調查資料,該區段采空區普遍存在礦柱壓裂、頂板顯著下沉、底鼓現象,該段地表羅家屋場、邊家屋場一帶也產生了地面塌陷(照片6、照片7)。

表1 樟村坪國家礦山公園地質災害基本特征一覽表Table 1 List of characteristics of geological disasters in Zhangcunping national mine park

照片6 羅家屋場地面塌陷導致房屋倒塌Photo 6 Building collapses caused by ground collapse in Luojiawuchang

照片7 邊家屋場地面塌陷導致地坪開裂Photo 7 Floor cracks caused by ground collapse in Bianjiawuchang

3 采空區穩定性數值模擬分析

3.1 分析方法

本次運用數值模擬的方法對昌達樟村坪磷礦采空區進行研究。FLAC/FLAC3D分析軟件適合于連續巖土介質和巖土結構工程問題的求解。該程序建立在拉格朗日算法基礎上,適合模擬大變形和扭曲。FLAC/FLAC3D采用顯式算法來獲得模型全部運動方程(包括內變量)的時間步長解,從而可以追蹤材料的漸進破壞和垮落。

3.2 計算幾何模型及邊界條件

計算模型[3]選取垂直于走向的一段山體進行三維模擬計算,計算軟件采用FLAC3D。三維計算模型在水平方向的長×寬尺寸為:408 m×152 m。模型的寬度方向為礦體的走向方向,長度方向為礦體傾向方向,模型高度范圍為標高900 m至地表。

計算模型共劃分為92 550個六面體單元,98 952個節點。計算模型的幾何形態及單元劃分如圖3所示。計算模型中的礦柱及間距,有實測數據的按實際尺寸處理,對于沒有實測數據的則按礦山標準采礦方法設計參數確定。

計算模型中,上部邊界自由;左、右邊界的水平位移固定約束,垂直位移不約束;底部邊界的垂直位移固定約束。本構模型采用Mohr-Columb彈塑性模型。

圖3 三維數值模型Fig.3 Three-dimensional numerical model

3.3 巖體力學參數與初始應力場

力學參數依據長沙礦冶研究院提供的湖北宜昌磷礦巖石力學參數測定報告,并采用工程地質類比法,借鑒類似礦區所得到的巖石力學參數來獲得,初始應力場按自重應力場計算。具體計算參數參見表2,表中的體積模量和剪切模型按照下式得出:

K= (1)

數值模擬計算中,以壓應力為負,拉應力為正。

3.4 計算結果

數值模擬計算根據采空區條件,分別設定多種工況計算,現分述如下:

3.4.1 當前工況

在當前開挖情況下,巖體最大變形區位于采空區頂面,隨高度增加位移減小,山體西側位移量大于東側,在西側坡面下部位移較大。應力變化分布較均勻,在最大壓應力出現在不規則礦柱分布的采空區區域,最大拉應力基本上沿著開挖層的底板分布(圖4)。底板和部分礦柱產生受拉和剪切屈服,屈服礦柱占總礦柱的10%左右,主要集中在山峰的下方,與該位置承受的巖體自重壓力較大有關,頂板也出現有零星的剪切屈服區,但分布較分散。地表水平最大位移出現在山峰處,地表水平位移規律基本上是山峰左側往左位移,山峰右側往右位移。

圖4 當前工況下山體位移等值線圖Fig.4 Contour map of the mountain displacement at the present condition

由上可見,在當前工況下,采空區及山體整體是穩定的,但頂板、礦柱和底板均出現屈服破壞區,表明其穩定狀態正逐步向臨界狀態轉化。

3.4.2 礦柱破壞工況

為了預測隨著采空區礦柱的破壞,采空區圍巖力學狀態狀況和變化趨勢,分別計算了在礦柱有10%、30%破壞和局部礦柱全部垮塌等情況下的變化。

隨著礦柱的減少,山體整體位移值呈增加趨勢,開采空間的屈服面積在逐漸增大,具體表現為礦柱部分主要是剪切破壞,底板和頂板表現為張拉破壞。當目前礦柱破壞的數量達到30%后,頂板的張拉屈服區顯著增加(圖5),同時底板和礦柱會出現連片的大面積屈服破壞,因此當礦柱破壞數量超過20%后采空區穩定性狀態將可能嚴重惡化。在現有條件下,若采空區出現局部一定范圍內礦柱全部破壞失效的情況,該區域的頂、底板位移量呈扇形分布形式增加,顯著增加的范圍集中在附近有限的區域內,頂板產生應力集中。四周的礦柱和底板出現屈服。

圖5 礦柱破壞30%時整體位移分布圖Fig.5 Overall displacement distribution map when the pillar is destroyed by 30%

3.4.3 礦柱力學性能削弱工況

為了研究礦柱經過長期的風化、蠕變、爆破震動、人為破壞等自然和人為因素的作用,力學特性弱化引起的采空區穩定性狀態的變化,計算了該工況下,采空區圍巖位移、應力變化特征。計算結果顯示,在礦柱的力學性能逐漸削弱下降過程中,如采空區圍巖屈服破壞區域顯著增加(圖6),礦柱幾乎完全屈服,位移量顯著增加,可能導致頂板的冒頂、垮落。說明大規模采空區經過長期的弱化,必然導致破壞,所以在采空區治理中,采用加固支撐方法長久維護空區穩定是不可取的。

圖6 礦柱力學性能降低后整體位移圖Fig.6 Overall displacement map after the reduction of mechanical properties of the pillar

4 采空區治理方法建議

采空區處理通常采用封閉、崩落、加固和充填四種方法,有時也聯合使用這幾種方法。

樟村坪磷礦現有空區體積約260×104m3,采用充填治理經濟性差、組織實施復雜、周期長,在現有條件下,實施難度較大。采用加固支撐措施,工程量較大,只能起到減緩推遲頂板冒落時間的作用,安全性較差。

采用強制崩落法處理采空區,其結果對礦山公園而言具有不可控性。故均在樟村坪礦不適用。

樟村坪磷礦已批準建設成為礦山公園,在保證采空區不影響到旅游和周邊環境安全的前提下,采用自然崩落法與封堵相結合的方法比較符合礦山當前的實際情況。

5 結語

(1) 區內地質災害發育、分布受井下采空區頂板穩定狀況、礦柱完整情況影響明顯。地面塌陷區大致沿著淺層采空區或淺層井巷發育,變形區段采空區普遍存在礦柱壓裂、頂板顯著下沉、底鼓現象,其地表山體形成多處開裂、危巖。

(2) 利用地形圖、采空區及地質災害的分布資料,結合FLAC3D分析技術,建立了三維數值模型。通過模擬顯示,在當前工況下,采空區及山體整體是穩定的;若采空區礦柱破壞加劇及礦柱力學性能削弱,采空區穩定性狀態將可能惡化。

(3) 樟村坪國家礦山公園目前已啟動實施一期建設——礦山地質環境治理。采空區穩定性研究為下一步劃分治理區域,結合礦山公園規劃因區施策具有指導意義。

[1] 呂小虎.湖北省宜昌磷礦樟村坪礦區Ⅲ礦段2011年度資源儲量報告[R].宜昌:湖北昌達化工有限責任公司,2012.

[2] 郭峰,馬振,閻巍.宜昌市夷陵區地質災害詳細調查報告[R].武漢:湖北省水文地質工程地質勘察院,2016.

[3] 汪曉霖.湖北昌達化工有限責任公司宜昌樟村坪磷礦井下采空區治理方案[R].武漢:中鋼集團武漢安全環保研究院,2006.

Study on the Stability of Goaf in Zhangcunping National Mine Park

GUO Feng1, WANG Shuihua1, JIANG Xinhong2, QU Jing2

(1.HubeiInstituteofHydrogeologyandEngineeringGeology,Yichang,Hubei443000; 2.EighthGeologicalBrigadeofHubeiGeologicalBureau,Xiangyang,Hubei441002)

Zhangcunping national mine park is the first phosphate park in China with widespread goafs.There are many geological disasters formed in collapsed goafs.In the study of goafs,through research on development and distribution characteristics and formation conditions of goafs and geological disasters in the study area,some targeted prevention measures are put forward using numerical simulation,which have important practical significance to development and utilization of national mine parks and disaster prevention and reduction in related mine.

mine park; goaf; stability

TD853.391+.2

A

1671-1211(2017)06-0746-06

10.16536/j.cnki.issn.1671-1211.2017.06.015

2017-08-21;改回日期2017-09-25

郭峰(1971-),男,高級工程師,水文地質與工程地質專業,從事水文地質、工程地質、環境地質勘查、調查、設計工作。E-mail:582404490@qq.com

數字出版網址:http://www.cnki.net/kcms/detail/42.1736.P.20171027.0937.004.html數字出版日期2017-10-27 09:37

李雯)