風景名勝區下銅礦開采地表沉降安全影響研究

李同鵬 胡 崴 方勝勇 李何林 婁廣文 1

（1.中鋼集團馬鞍山礦山研究總院股份有限公司，安徽 馬鞍山 243000；2.金屬礦山安全與健康國家重點實驗室，安徽 馬鞍山 243000；3.華唯金屬礦產資源高效循環利用國家工程研究中心有限公司，安徽 馬鞍山 243000；4.銅陵有色股份鳳凰山礦業有限公司，安徽 銅陵 244000）

礦山地下開采對地表的安全影響是地下礦山開采過程中普遍存在的問題，開采過程中產生的地表沉降不僅會造成礦區生態環境破壞，還可能導致地表建構筑物破壞，嚴重時還將造成地表塌陷等地質災害。隨著礦產資源的持續開發，大部分易采礦體已經基本開采殆盡，未來越來越多的地下礦山面臨對建構筑物、水體等壓覆的礦產資源進行回采。且近年來，國家全面推進綠色礦山建設，在新時代生態文明建設和綠水青山就是金山銀山的發展理念下，為保護地表生態環境和建構筑物安全，對控制地表沉降造成的影響也提出了更高的要求［1-3］。

鳳凰山銅礦位于鳳凰山風景名勝區下方，礦區及周邊地表分布有景區已有的和規劃建設的各類建構筑物，礦山開采可能對景區地表建構筑物安全造成影響。本研究在對礦山現狀和地質條件進行全面調查和分析的基礎上，采用理論計算和數值模擬分析手段，研究礦山地下開采對風景名勝區的安全影響，為礦山后續開采和景區建設安全提供理論依據。

1 工程概況

鳳凰山銅礦始建于1965年，自1971年正式投產至今，已有近50年的開采歷史，井下共形成了8個中段，其中-440 m水平以上礦體已全部開采結束，前期形成的采空區均已采用尾砂和廢石進行了充填或已冒落塌實。目前主要開采范圍為-500 m中段和-560 m中段礦體，設計采用上向分層水平充填法進行開采，礦塊沿走向布置，長度為50 m，寬度為礦體水平厚度，間柱寬4～5 m，頂柱3～5 m，底柱7 m，礦房中間留有?5 m不規則的礦柱，采場高60 m。礦體與地表相對位置見圖1。

2 臨時采空區冒落安全影響分析

根據設計和生產實際，礦山-440 m以下全部采用上向分層充填采礦法開采，分層回采高度3 m，采場分層礦石出完，形成5～6 m高的臨時采空區后進行充填。故空頂高度最大為6 m。若形成的采空區未及時充填，則可能造成臨時采空區垮塌。本次研究借鑒“三帶”理論計算經驗公式，通過計算冒落帶、裂隙帶高度分析地下開采對地表的安全影響。

冒落帶高度［4-5］

裂隙帶高度［4-5］

式中，M為空頂高度。

根據計算，冒落帶高度為10.5～14.9 m，裂隙帶高度為39.85～51.05 m，則臨時采空區冒落影響高度最大為14.9 m+51.05 m=65.95 m。即鳳凰山銅礦-440 m水平以下按設計采礦方法和參數進行開采時，采空區垮塌造成的最大影響高度為-440 m水平以上65.95 m，不會對地表安全造成影響。

3 開采沉陷數值模擬

基于礦山開采現狀和設計情況，利用FLAC3D軟件對礦山不同開采階段下的開采及充填過程進行模擬，分析地表沉降和變形的總體情況和對風景區建構筑物的影響。

（1）模型建立。為合理模擬鳳凰山銅礦地下開采對風景區建構筑物的安全影響，在充分研究已有地質資料與多次試算基礎上，確定了數值模型的平面范圍，模型平面范圍東西向長1 200 m，南北向寬1 000 m。礦區地表分布有排土場和露天采坑，露天采坑底部采用膠結尾砂充填，根據礦區實測地形圖和各中段平面圖建立了地表和礦體模型，見圖2。

（2）巖體力學參數取值。依據礦山提供的巖石力學強度指標及礦區工程地質條件，參考采礦設計手冊中巖石力學參數的經驗數據及工程巖體質量分級標準，通過折減，確定本次數值計算模型中圍巖及礦體材料物理力學參數，詳見表1。

（3）模擬過程。為盡可能真實地反映鳳凰山銅礦礦體開采及覆巖的變形特征，數值計算按以下幾個步驟進行：①形成初始應力場；②模擬鳳凰山銅礦-440m水平以上已采礦體回采和充填過程；③模擬-440 m水平以下（-500 m中段和-560 m中段）礦體回采和充填過程。每個步驟的計算均在前一步計算基礎上連續進行，并且逐步推進和計算，從而客觀地反映多個過程的疊加效應［6-8］。

（4）-440 m以上已開采礦體開采數值模擬。-440 m以上為礦山前期已開采的礦體，包括-40 m、-80m、-120 m、-180 m、-240 m、-300 m、-360 m和-440 m共8個中段。-440 m以上各中段開采后的地表垂直位移和水平位移云圖見圖3。通過模擬分析可知，-440 m以上開采結束后，地表產生的最大垂直位移約41.05 cm，主要集中在露天采坑底部，其四周主要表現為整體均勻沉降。最大水平位移約8.47 cm，水平變形影響區主要位于露天采坑底部及其西側區域，數值模擬結果與礦山現狀情況較為吻合。

（5）-440 m以下礦體開采數值模擬。在-440 m以上礦體回采后地表產生的應力和位移場的基礎上，對-440 m以下礦體開采過程進行模擬分析。-440 m以下（-500 m中段和-560 m中段）回采后的地表垂直位移和水平位移見圖4。通過對-440 m水平以下礦體開采過程模擬分析可知，隨著開采過程的進行，地表變形進一步增大，地表最大垂直位移約43.75 cm，最大水平位移約8.96 cm，比-440 m以上開采結束時垂直位移和水平位移分別增大2.70 cm和0.39 cm。

4 地表沉降對建構筑物安全影響分析

鳳凰山銅礦地表需保護的建構筑物主要有采礦遺址、戶外露營基地、道路等景區設施，根據《GB 50771—2012有色金屬采礦設計規范》中的等級劃分，其保護等級分別為Ⅰ級和Ⅱ級，對建構筑物的影響程度主要考慮傾斜、曲率和水平變形值，變形允許值見表2。

根據鳳凰山銅礦開采結束后地表垂直位移云圖（見圖5），在建構筑物附近各剖面線（A線～I線）上每隔2 m提取1個垂直位移數據，每條剖面線提取了30個垂直位移數據，根據公式［9-12］計算得出了地表建構筑物附近的最大傾斜為0.75 mm/m；最大曲率為0.17 mm/m，滿足Ⅰ級建筑物變形允許值要求。各剖面線建構筑物附近的垂直位移曲線見圖6。

根據鳳凰山銅礦開采結束后地表水平位移云圖（見圖7），在建構筑物附近各剖面線（A線～I線）上每隔2 m提取1個水平位移數據，每條剖面線提取了30個水平位移數據，計算得出了地表建構筑物附近的最大傾斜為0.44 mm/m，滿足Ⅰ級建筑物變形允許值要求，各剖面線建構筑物附近的水平位移曲線見圖8。

5 結 論

（1）根據理論計算分析，鳳凰山銅礦深部-440 m以下采用上向水平充填采礦法開采，僅有單層5～6 m高的臨時采空區，采空區坍塌最大影響高度約為65.95 m，不會造成地表出現明顯沉降或塌陷，可以保障風景區地表安全。

（2）根據數值模擬計算分析，鳳凰山銅礦深部-440 m以下開采所產生的地表最大垂直位移增量為2.69 cm，最大水平位移增量為0.39 cm；地表建構筑物所處位置產生的傾斜、曲率和水平變形值均在《GB 50771—2012有色金屬采礦設計規范》所規定的Ⅰ級保護建構筑物變形允許值范圍內，可以保障風景區地表建構筑物安全。

（3）本次研究為鳳凰山銅礦深部開采和鳳凰山風景名勝區建設安全提供了理論依據，推動了礦山開采與風景名勝區建設共生共存、協調發展，同時也為類似礦山地表沉降安全影響分析提供了借鑒。