張耀平 楊雙有 盧志剛 蘇曉萍 夏弋江

(1.江西理工大學(xué)應(yīng)用科學(xué)學(xué)院;2.江西理工大學(xué)資源與環(huán)境工程學(xué)院)

礦區(qū)地表沉陷對礦區(qū)生產(chǎn)安全有重大影響，并會破壞周邊自然環(huán)境，同時造成耕地?zé)o法耕種，房屋開裂毀壞等［1］。因此，對地表監(jiān)測，并對其沉降規(guī)律進行分析研究將具有重要意義。目前，常規(guī)研究方法通常是通過地表監(jiān)測，獲得相關(guān)數(shù)據(jù)，建立二維坐標曲線［2］。這種方法分析不夠直觀形象，可視化太低，分析片面，準確性不高，預(yù)測性能差，不利于地表沉降防治措施的制定和實施［3］。本研究利用SURFER軟件對地表監(jiān)測點數(shù)據(jù)進行可視化處理，繪制出地表沉降三維表面圖，地表沉降三維切片圖，地表監(jiān)測點張貼圖，大大提高了研究的直觀性，可靠性，準確性，預(yù)測性，為礦山地表沉陷防治措施的制定和實施，提供可靠準確依據(jù)［4-5］。

1 采空區(qū)地表沉降規(guī)律三維可視化研究

1.1 地表的移動破壞形式

在礦體的開挖過程中，逐漸形成空區(qū)。同時，由于采動的影響，地表緩慢下沉，從而在空區(qū)上方地表形成一個較大的下沉區(qū)域，即地表移動盆地。該移動趨勢蔓延到地表后，隨采空區(qū)再次開挖，暴露時間的增長，地表下沉范圍及最大沉降累積量也逐漸增大，逐漸形成“碗狀”，即非充分采動。隨開采工作向前進行，空區(qū)繼續(xù)下沉，直到最大沉降累積量不再隨采空區(qū)增大，但盆地形態(tài)仍呈“碗狀”，即充分采動。開采活動繼續(xù)進行，同時移動盆地的沉降累積量不變，但盆地呈平底的“盆狀”，此時，地表的采動影響稱超充分采動［6-9］。

此外，常見的地表的移動破壞形式還有塌陷坑裂縫及臺階。當(dāng)?shù)V層為急傾斜形態(tài)時，開采過程容易出現(xiàn)塌陷坑。當(dāng)該礦層松散層較薄時，則容易引起裂縫或臺階。

1.2 地表變形監(jiān)測

(1)地表水平位移監(jiān)測。采用萊卡全站儀和GPS，監(jiān)測空區(qū)地表水平位移。GPS針對工作基點實施監(jiān)測，全站儀針對監(jiān)測點實施測量。監(jiān)測點的水平位移使用全站儀TC402觀測，主要測量監(jiān)測點的點位相對水平位移變化情況。在有條件的情況下，布設(shè)3～5個GPS基準點，與礦山已有的GPS地面控制點進行聯(lián)測。水平位移觀測周期，對于地表比較穩(wěn)固區(qū)域的觀測，亦可與沉降觀測協(xié)調(diào)考慮確定，或變形情況及工程進展情況而定。

(2)地表沉降位移監(jiān)測。水準觀測由于地勢起伏大，導(dǎo)致精度降低，特別是對于坡度大的地區(qū)，受水準尺和三腳架高度影響，無法進行觀測。GPS由于垂直位移精度尚不理想，無法滿足要求，監(jiān)測范圍比較小，因此考慮采用三角高程方法進行垂直位移監(jiān)測，使用全站儀TC402觀測，測量變形監(jiān)測點的點位垂直位移變化情況。沉降監(jiān)測周期

式中，m為2沉降觀測點之間的高差誤差;v為沉降速度，一般取平均沉降量與間隔時間的比值;K為變形值與其誤差之比，在5～10之間選擇。

1.3 基于SURFER的地表沉降可視化處理與分析

SURFER軟件功能強大，可以輕松制作線框圖及三維表面圖等［10-12］。研究通過對空區(qū)地表進行水平位移和沉降位移監(jiān)測，獲得地表移動數(shù)據(jù)，將該數(shù)據(jù)利用導(dǎo)入SURFER軟件中，對地表移動數(shù)據(jù)進行可視化處理，分析地表沉降規(guī)律。其操作過程見圖1。

2 工程實例

安徽某鐵礦+100 m以上水平分層已采空，形成多個空區(qū)。目前正進行+100 m以下水平分層的地下開采，井下施工中沿及西沿均掘進至0線與1線中間部位，2條巷道均為沿脈巷道，在礦層中掘進。礦層結(jié)構(gòu)比較松軟，礦石較破碎，穩(wěn)定性差，巷道掘進維護較困難。礦層產(chǎn)狀不明顯，但總體向西傾，厚度大，礦層傾角總體5°～15°。礦山的地質(zhì)和地貌情況比較復(fù)雜:地勢高低不平，起伏大，測區(qū)范圍不大，礦區(qū)植被分布廣泛。為保證地表生產(chǎn)安全，礦山對空區(qū)地表進行定時監(jiān)測。

圖1 地表下沉可視化處理分析過程

借助SURFER軟件對監(jiān)測數(shù)據(jù)進行可視化處理，繪制出地表沉降三維表面圖見圖2，地表監(jiān)測點張貼圖見圖3，地表沉降切片圖見圖4，自下向上為4—5月和7—9月，并利用Excel生成地表沉降量如圖5。最后，結(jié)合地表的移動破壞形式，從地表沉陷形狀，沉陷蔓延趨勢，最大沉降累積量等角度，分析了空區(qū)地表沉降規(guī)律，并基于礦山安全考慮給出合理建議。

由圖2～圖4可知:由于礦體開挖，空區(qū)的形成，地表部分地區(qū)逐漸下沉，在點X4，X5左側(cè)形成碗狀形態(tài)，X9，X3點附近向上凸起。隨著空區(qū)暴露時間的增長，地表沉降范圍進一步向X6，X7，X8點附近擴大，沉降累積量逐漸增加，碗狀形態(tài)逐漸趨于盆狀，即形成下沉盆地。但由于地表沉降范圍向X4，X5點擴大不明顯，則碗狀向沉陷盆地轉(zhuǎn)變時，無法形成規(guī)則的下沉盆地。XN附近點剛開始沉降非常明顯，后面基本沒有變化。

由圖3～圖5可知:第2次監(jiān)測最大沉降累積量為－15 mm，位于XN點，其次為－8 mm，位于X4點，最近一次監(jiān)測最大沉降累積量為－13 mm，位于X7點。說明沉降范圍從 X4，X5點左側(cè)向 X6，X7，X8點附近擴大的同時，最大沉降點也隨之轉(zhuǎn)移。XN附近點剛開始沉降非常明顯，為－15 mm，后面略有回升，基本上保持在－13 mm左右。隨開挖的持續(xù)進行，暴露面積進一步增大，地表變形經(jīng)歷了碗狀形態(tài)和下沉盆地，并在繼續(xù)擴大，但最大沉降累積量趨于穩(wěn)定。

圖2 地表沉降三維表面

圖3 地表監(jiān)測點張貼示意

圖4 地表沉降切片示意

圖5 地表沉降量

綜上所述，該空區(qū)地表經(jīng)歷了碗狀變形，并最終形成了不規(guī)則的地表移動盆地，同時，沉降范圍仍在繼續(xù)擴大。XN點沉降累積量幾乎保持不變，不再下沉，但該點附近下沉后的形態(tài)及范圍不明朗。出于礦山地表生產(chǎn)安全考慮，應(yīng)在原有監(jiān)測點的基礎(chǔ)上增設(shè)監(jiān)測點定期監(jiān)測，即應(yīng)在XN附近增設(shè)監(jiān)測點，以進一步確定其沉降規(guī)律，同時在 X4，X5，X6，X7，X8點左側(cè)增設(shè)至少1排監(jiān)測點以確定沉降范圍、下沉值等。

3 結(jié)論

(1)利用SURFER軟件對地表監(jiān)測點數(shù)據(jù)進行可視化處理，繪制出地表沉降三維表面圖、地表沉降三維切片圖、地表監(jiān)測點張貼圖。大大提高了研究的直觀性、可靠性、準確性、預(yù)測性，為礦山地表沉陷防治措施的制定和實施提供了可靠依據(jù)。

(2)通過綜合分析，該空區(qū)地表沉降規(guī)律為空區(qū)地表經(jīng)歷了碗狀變形，并最終形成了不規(guī)則的地表移動盆地，同時，沉降范圍仍在繼續(xù)擴大。XN點沉降累積量幾乎保持不變，不再下沉，但該點附近下沉后的形態(tài)及范圍不明朗。

(3)基于分析結(jié)果，本研究認為應(yīng)在XN附近增設(shè)監(jiān)測點，以進一步確定其沉降規(guī)律，同時在X4，X5，X6，X7，X8左側(cè)增設(shè)至少1排監(jiān)測點以確定沉降范圍、下沉值等。

［1］ 孫 超，薄景山，劉紅帥，等.采空區(qū)地表沉降影響因素研究［J］.吉林大學(xué)學(xué)報:地球科學(xué)版，2009(3):498-502.

［2］ 趙德深，范學(xué)理.礦區(qū)地面塌陷控制技術(shù)研究現(xiàn)狀與發(fā)展方向［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，2001(2):89-92.

［3］ 羅澤華.凡口鉛鋅礦地表沉降監(jiān)測與研究［J］.礦山測量，2003(2):3-4，8.

［4］ 楊慶坤，郭福生，王繼堯，等.Surfer軟件在煤層儲量計算和可視化中的應(yīng)用［J］.煤炭技術(shù)，2010(1):181-184.

［5］ 楊宇宏，許 勝.基于Surfer和ANSYS復(fù)雜地質(zhì)體的FLAC3D實體建模［J］.山西建筑，2011(5):63-64.

［6］ 張忠堯，張國駿，張忠華.地面沉降塌陷對區(qū)域環(huán)境經(jīng)濟的影響及防治對策研究［J］.北方環(huán)境，1999(4):20-21.

［7］ 劉 毅.地面沉降研究的新進展與面臨的新問題［J］.地學(xué)前緣，2001(2):273-278.

［8］ 羅 敏，劉巨保.開采區(qū)塊地表沉降位移的有限元分析［J］.大慶石油學(xué)院學(xué)報，2005(3):55-58，124.

［9］ 徐福喜.萬家寨引黃工程北干線煤礦采空區(qū)變形監(jiān)測分析［J］.水利水電技術(shù)，2007(5):62-64.

［10］ 鄭睿博，周榮福，徐海林，等.基于Surfer及AutoCAD軟件的等厚線圖繪制［J］.中州煤炭，2010(10):92-93.

［11］ 邰軍良.用VB和Surfer繪制彩色平面等值線圖［J］.西部探礦工程，2011(12):122-124.

［12］ 馬勝利，賈東方，肖新成.AutoCAD及Surfer軟件在煤層等厚線圖繪制中的應(yīng)用［J］.中州煤炭，2009(9):91-92.