礦區(qū)地表觀測站輔助設計與數(shù)據(jù)處理程序實現(xiàn)

方 齊 郭廣禮 朱曉峻

(1.中國礦業(yè)大學環(huán)境與測繪學院，江蘇 徐州 221116；2.國土環(huán)境與災害監(jiān)測國家測繪地理信息局重點實驗室，江蘇 徐州 221116；3.江蘇省資源環(huán)境信息工程重點實驗室，江蘇 徐州 221116)

礦區(qū)地表觀測站輔助設計與數(shù)據(jù)處理程序實現(xiàn)

方 齊1,2,3郭廣禮1,2,3朱曉峻1,2,3

(1.中國礦業(yè)大學環(huán)境與測繪學院，江蘇 徐州 221116；2.國土環(huán)境與災害監(jiān)測國家測繪地理信息局重點實驗室，江蘇 徐州 221116；3.江蘇省資源環(huán)境信息工程重點實驗室，江蘇 徐州 221116)

為了提高礦區(qū)地表移動觀測站設計的自動化程度、數(shù)據(jù)處理速度以及繪制移動變形曲線圖的效率，基于Visual Basic(VB)及Visual Basic for Application (VBA)編程語言，并融合Access數(shù)據(jù)庫、AutoCAD、Excel軟件技術，實現(xiàn)了地表移動觀測站設計中的參數(shù)查詢、觀測線長度計算、觀測線布設、數(shù)據(jù)輸出、后期觀測數(shù)據(jù)處理及繪制曲線等一系列功能。首先在AutoCAD軟件中快速進行工作面上方地表、線狀路線觀測站的布設；然后采用Excel軟件自動化處理觀測數(shù)據(jù)并繪制變形曲線圖，并改進了Excel軟件中繪制帶平滑線散點圖坐標軸的標注功能；最后將上述2個功能模塊集成到同一程序中，從而克服以往程序功能單一、自動化程度差、修改不便等缺陷。結果表明，該程序能夠適應多樣化的布站需求，滿足繪制地表移動變形曲線圖的特殊要求，有助于提高礦區(qū)地表觀測站設計的自動化程度，為建立地表移動觀測系統(tǒng)提供參考。

地表移動觀測站 數(shù)據(jù)處理 移動變形曲線 VB VBA 地表移動觀測系統(tǒng)

目前對于礦區(qū)地表移動觀測應用程序的開發(fā)得到了學者們的關注，如朱群康基于VB.NET開發(fā)了地表移動觀測站設計系統(tǒng)[1]，張姣姣對礦區(qū)地表移動觀測成果的自動化求解方法進行了研究[2]。但現(xiàn)有的相當一部分地表移動觀測程序僅是針對某一特定功能開發(fā)而成的，觀測線布設的便利性與曲線圖的顯示效果欠佳，地表移動曲線圖與觀測數(shù)據(jù)無法進行同步修改，程序系統(tǒng)的適用性有待進一步完善。

為此，兼顧程序的實用性及礦區(qū)變形監(jiān)測設計和數(shù)據(jù)處理的特殊性，采用VB、VBA語言進行程序開發(fā)，程序功能包括：①將地表移動觀測站設計及數(shù)據(jù)處理模塊作為2個獨立的功能插件，提高本程序對其他編程語言的兼容性，并將該2個功能插件集成到同一可執(zhí)行程序中，解決插件的加載問題，為后續(xù)地表移動觀測系統(tǒng)建立框架；②建立礦區(qū)角值參數(shù)數(shù)據(jù)庫，能夠進行模糊查詢和快速檢索，豐富觀測站布設的快捷功能，如連續(xù)或間斷標注點號、不同符號標注控制點及觀測點、快速創(chuàng)建不規(guī)則觀測線和公路、鐵路等線狀路線觀測線、常用繪圖屬性的二次更改等，提高了設計的便捷性，簡化了觀測線布設的步驟；③數(shù)據(jù)處理繪圖模塊中，觀測數(shù)據(jù)可在Excel軟件中進行輸入和批量計算，擴展了Excel軟件中帶平滑線散點圖橫軸的標注功能，實現(xiàn)了橫軸標注為文本類型、非等間隔標注及采空區(qū)位置對照等功能，曲線樣式調整方便，便于數(shù)據(jù)添加和圖形維護，數(shù)據(jù)與曲線能夠同步調整，從而提高了工作效率。

1 地表移動觀測站設計與數(shù)據(jù)處理方法

1.1 地表移動觀測站設計

地表移動觀測站按布設形式分為網狀和剖面線狀2類，本研究采用剖面線狀觀測站布設形式進行程序設計[3-4]。地表移動觀測站設計主要包括觀測線位置的確定和觀測線長度的計算2部分內容，如圖1所示。

1.1.1 走向線位置及長度確定

走向觀測線應布置于移動盆地的走向主斷面上，走向觀測線位置偏移工作面中心點的距離S可通過下式計算：

(1)

式中，H1、H2為分別為采區(qū)下邊界和上邊界的開采深度，m；θ為最大下沉角，(°)。

任意一條走向觀測線的長度通過下式計算：

LZ=2hcotφ+2(H0-h)cot(δ-Δδ)+l，

(2)

式中，φ為松散層移動角，(°)；H0為回采工作面平均開采深度，m；h為松散層厚度，m；δ為走向移動角，(°)；Δδ為走向移動角修正值，(°)；l為工作面走向長度，m。

圖1 觀測站設計示意

1.1.2 傾向觀測線位置及長度確定

傾向觀測線應布設于充分采動區(qū)內，考慮鄰近開采的影響，傾向觀測線到開切眼的距離D1及到停采線的距離D3應分別滿足如下條件：

(3)

式中各變量意義同上。

程序設計中采用式(3)作為傾向觀測線位置選擇的約束依據(jù)，當?shù)乇砦催_到充分采動時，傾向觀測線的位置由采區(qū)工作面的中心位置確定。任意一條傾向觀測線的長度可由下式計算：

LQ=2hcotφ+(H1-h)cot(β-Δβ)+

(H2-h)cot(γ-Δγ)+Lcosα，

(4)

式中，β為下山移動角，(°)；Δβ為下山移動角修正值，(°)；γ為上山移動角，(°)；Δγ為上山移動角修正值；L為工作面傾斜長度，m；α為工作面傾角，(°)；其余各變量意義同上。

1.2 觀測數(shù)據(jù)處理方法

觀測數(shù)據(jù)前期預處理采用文獻[5]所介紹的坐標規(guī)劃方法，觀測數(shù)據(jù)后期處理主要包括：①各測點的下沉和水平移動值計算；②相鄰兩測點間的傾斜和水平變形值計算；③相鄰兩線段(或相鄰三點)的曲率變形值[4]計算。

1.3 地表移動變形曲線圖

地表移動變形曲線圖包括下沉、傾斜、曲率、水平移動和水平變形曲線。繪制移動變形曲線圖應滿足如下要求：①根據(jù)下沉、水平移動數(shù)據(jù)的正、負號分別將其展在測點的上下兩側；②依據(jù)傾斜和水平變形數(shù)據(jù)的正、負號將其展在兩測點間中點的正上方或正下方；③曲率數(shù)據(jù)應展在相鄰兩線段之和的中點的正上方或正下方。展點完畢后，將各曲線點分別采用平滑曲線連接，完成曲線圖的繪制工作[6-7]。

2 程序實現(xiàn)

2.1 程序簡介

作為一種通用的自動化語言，VBA語言開發(fā)具有高效性并且語法簡單[8-10]，適合地表移動觀測站設計及數(shù)據(jù)處理功能模塊的開發(fā)。AutoCAD作為專業(yè)繪圖軟件，數(shù)據(jù)處理便捷性不夠，且曲線樣式調整不便，數(shù)據(jù)更改后圖形維護欠佳；Excel軟件中圖形與數(shù)據(jù)的關聯(lián)性較強，數(shù)據(jù)與圖形能夠同步更新，便于后期數(shù)據(jù)維護，程序運行速度較快。因此，利用VBA語言開發(fā)觀測站設計插件及數(shù)據(jù)處理繪圖插件，采用VB語言生成能夠自動加載的可執(zhí)行程序，將上述2個功能模塊作為插件調用，自動化加載VBA語言程序，使程序集成一體化，便于利用其他編程語言豐富后續(xù)功能，增加程序的可移植性。程序調用Excel軟件時則自動加載開發(fā)程序菜單，結果如圖2所示。

圖2 Excel軟件菜單自動加載

2.2 功能模塊

2.2.1 觀測站設計模塊

根據(jù)地表移動觀測站設計原理開發(fā)相應功能模塊，在AutoCAD軟件中輸入命令可直接執(zhí)行程序，并顯示功能模塊界面，如圖3所示。

2.2.1.1 參數(shù)查詢

程序建立了角值參數(shù)數(shù)據(jù)庫，實現(xiàn)數(shù)據(jù)共享[11]，數(shù)據(jù)庫文件中內置《建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規(guī)程》[12]內所列的各礦區(qū)的角值參數(shù)，根據(jù)所提供的關鍵字可模糊查詢到包含關鍵字在內的礦區(qū)名稱及所涉及的角值參數(shù)，便于使用類比法選取角值參數(shù)設計觀測線，界面如圖4所示。

圖3 AutoCAD功能模塊界面

圖4 參數(shù)查詢界面

程序部分核心代碼如下：

Set adoCon=New ADODB.Connection

adoCon.ConnectionString="Provider=Microsoft.jet.OLEDB.4.0;data source=" & pathexe & "parater.mdb;"adoCon.Open

Set daoRs = New ADODB.Recordset

daoRs.Open "gczcs"，adoCon，adOpenForwardOnly，adLockOptimistic

2.2.1.2 參數(shù)設置

參數(shù)設置界面主要包括基本參數(shù)、走向參數(shù)、傾向參數(shù)及繪圖參數(shù)等子界面，觀測線位置、長度的確定，均按照文中第1節(jié)所介紹的方法進行計算。繪圖參數(shù)界面功能主要用于作圖設置，批量修改，快速對字體類型、字體高度、測點尺寸、控制點及監(jiān)測點圖層顏色進行二次修改。在作圖之前，除了需要指定上述參數(shù)外，還需要輸入監(jiān)測點間斷標注間隔、控制點號前綴、監(jiān)測點號前綴等參數(shù)。

2.2.1.3 功能區(qū)

功能區(qū)主要是軟件各項功能按鈕，包括加密測點、標注點號、參數(shù)查詢、創(chuàng)建觀測線等。

2.2.1.4 功能特點

觀測線設計類型主要有規(guī)則觀測線和非規(guī)則觀測線2種形式，本程序兼顧該2種形式實現(xiàn)觀測線設計功能，以第1節(jié)介紹的設計原理為依據(jù)，對Auto CAD軟件進行二次開發(fā)[13-14]，通過軟件定義的順序選取工作面角點坐標，將工作面近似為矩形工作面，并計算其中心點位置及方位角等參數(shù)，依次輸入?yún)?shù)進行觀測線的設計。對于觀測線的位置可以通過計算得到或鼠標選取，相關命令提示顯示在AutoCAD軟件中的命令行中，并在工作面自動標出傾向觀測線位置到開切眼、停采線的設計最小距離，滿足傾向線觀測位置選擇的設計要求。傾向觀測線位置選擇部分核心代碼如下：

getpointstyle = 1

getpointstyle = ThisDrawing.Utility.GetInteger("默認基點(1)或[自定義基點(2)]:")

If getpointstyle = 2 Then

basePick = ThisDrawing.Utility.GetPoint("選取點:")

basePnt1 = basePick(0)

basePnt2 = basePick(1)

對于復雜形狀的觀測線或在采動區(qū)內布設公路、鐵路等線狀路線觀測線時，測點位置受地形影響較大，多為不規(guī)則觀測線[15-16]，通過對類似程序分析[17-20]，本程序添加了輔助設計功能，利用軟件的快捷功能最大限度提高自動化程度，利用建立的規(guī)則觀測線的影響范圍，采用少許點繪制出觀測線的大致趨勢線，之后通過創(chuàng)建測點命令，建立觀測線，并可一鍵選擇監(jiān)測點、控制點及觀測線層，整體移動測點位置及點號。由于某些重要位置需要加密觀測點，可以選取“加密觀測點”命令進行加密。觀測線中的監(jiān)測點、控制點布設完畢后，點擊“創(chuàng)建點號”，即可完成觀測線點號標記、間斷標注或全部標注，同時提取設計點坐標存入 “.txt”格式文件中。布置多條觀測線時，重復上述步驟即可。

2.2.2 數(shù)據(jù)處理及繪圖模塊

地表移動觀測站周期觀測之后，需要計算移動變形值，將每一次觀測數(shù)據(jù)繪制在同一觀測線的剖面圖上，用以研究地表移動變形規(guī)律。Excel軟件提供了強大的數(shù)據(jù)處理及繪圖功能，但是帶平滑線的散點圖不支持對于橫坐標軸標注為文本類型以及不等間隔標注等功能，而實際觀測站任意2點之間的距離因地形的影響不一定完全相等，如果按等間隔標注會對研究地表移動變形規(guī)律產生影響；另外，點號通常是文本型，并非單純的數(shù)字。

通過程序開發(fā)實現(xiàn)移動變形值計算以及繪制移動變形曲線圖功能，功能特點為：通過相關命令建立格式文件，輸入原始數(shù)據(jù)，進行單期或多期數(shù)據(jù)批量計算，輸出可視化計算結果，并可隨時調整數(shù)據(jù)，其中下沉、傾斜、曲率在同一表中，水平變形和移動在同一表中。變形值計算的部分代碼如下：

If Selection.Cells(j，i)= "" Then

ActiveSheet.Cells(j + rownum + 3，i)= 0

Else

ActiveSheet.Cells(j + rownum + 3，i)= Format((Selection.Cells(j，4)- Selection.Cells(j，i))* 1000，"0")

End If

通過以上計算，程序自動將計算結果按規(guī)定格式整理，之后選擇繪圖源數(shù)據(jù)，點擊按鈕，程序將橫軸標簽及標簽坐標分別作為橫坐標源數(shù)據(jù)列，自動繪制出橫坐標軸，之后將標簽及標簽坐標展點在橫坐標軸上，插入標簽標記，即可完成橫坐標文本標簽及非等間隔標注。另外，選擇的觀測數(shù)據(jù)將分別繪制擬合曲線，通過“采空區(qū)設置”功能，利用線條填充繪制采空區(qū)位置示意圖，方便與點位信息對照。橫軸繪制的部分核心代碼如下：

i = mchtchart.SeriesCollection.Count

For n = 1 To cellnum

PositionY(n)= 0

PositionX(n)= Selection.Cells(n，1).Value

Next n

mchtchart.SeriesCollection.NewSeries

mchtchart.SeriesCollection(i + 1).XValues = PositionX

mchtchart.SeriesCollection(i + 1).Values = PositionY

mchtchart.SeriesCollection(i + 1).ChartType = xlXYScatter

3 應用實例

內蒙古自治區(qū)某礦首采充填工作面上方有公路、礦區(qū)專用鐵路等重要設施，地下開采可能會對其產生影響，故需在首采工作面上方布設地表移動觀測站及專用鐵路、公路觀測線，以便及時維護公路、鐵路等設施并積累該地質采礦條件下的觀測資料，根據(jù)地質采礦條件及設計參數(shù)，進行觀測線輔助設計，效果如圖5所示。

鐵路、公路的專用觀測線分別沿鐵路、公路鋪設，并超出沉陷影響范圍外一段距離。工作面走向上方地勢平坦，觀測條件較好，故沿直線布設走向觀測線，用于研究工作面走向主斷面上地表移動變形規(guī)律。

圖5 觀測站布設效果

安徽省某礦充填工作面長約310 m，寬約110 m。工作面對應上方地表有村莊、學校以及大型工業(yè)建筑等，地面情況復雜，為了控制地表沉陷及提高煤炭采出率，采用固體充填采煤技術進行壓煤開采，并建立了地表移動觀測站及巖層移動觀測站。2012年8月—2014年5月經過了29期觀測，積累了豐富的觀測數(shù)據(jù)，后期可用觀測點共24個。利用本研究所編寫的程序將其中6期的巖層移動觀測數(shù)據(jù)進行處理并繪制出下沉曲線圖，如圖6所示。

圖6 下沉曲線

通過與工作面位置對照，圖6所示的下沉規(guī)律與非充分采動時主斷面內的移動變形規(guī)律基本符合。

4 結 語

所編寫的程序能夠完成地表移動觀測站的輔助設計及后期觀測數(shù)據(jù)的處理繪圖任務，滿足規(guī)則觀測線設計要求，并可輔助設計非規(guī)則觀測線，增強了地表移動觀測數(shù)據(jù)處理繪圖自動化程度，減輕了設計處理的工作量，避免了大量重復性的工作內容，具有簡單實用、便于調用的特點，對于構建功能齊全的地表移動觀測設計處理系統(tǒng)具有一定的參考價值。

[1] 朱群康，趙興旺.基于VB.NET的地表移動觀測站設計系統(tǒng)開發(fā)[J].煤炭技術，2008(2)：127-130. Zhu Qunkang，Zhao Xingwang.Design of surface ground movement observation station based on VB.NET[J].Coal Technology，2008(2):127-130.

[2] 張姣姣，鄭禮權，徐良驥，等.礦區(qū)地表移動觀測成果的自動化求解[J].礦山測量，2013(2)：87-89. Zhang Jiaojiao，Zheng Liquan，Xu Liangji，et al.Automatic computation of observation results of surface movement in mining area[J].Mine Surveying，2013(2):87-89.

[3] 董合祥，朱子清.地表移動觀測站的設計研究[J].煤礦現(xiàn)代化，2010(1)：40-41. Dong Hexiang，Zhu Ziqing.The design and study on observation station of ground movement's foundation[J].Coal Mine Modernization，2010(1):40-41.

[4] 趙建紅，郭志磊.地表移動觀測站設計[J].煤炭技術，2010(4)：156-158. Zhao Jianhong，Guo Zhilei.The design on observation station of ground movement's foundation[J].Coal Technology，2010(4):156-158.

[5] 王明柱，郭廣禮，王 磊.礦區(qū)地表移動變形監(jiān)測數(shù)據(jù)處理[J].測繪科學，2012(4)：67-69. Wang Mingzhu，Guo Guangli，Wang Lei.Surface movement observation data processing in mining area[J].Science of Surveying and Mapping，2012(4):67-69.

[6] 何國清，楊 倫，凌賡娣，等.礦山開采沉陷學[M].徐州：中國礦業(yè)大學出版社，1991. He Guoqing，Yang Lun，Ling Gengdi，et al.Mining Subsidence Science[M].Xu Zhou:China University of Mining and Technology Press，1991.

[7] 郭廣禮，查劍鋒.礦山開采沉陷及其防治[M].徐州：中國礦業(yè)大學出版社， 2012. Guo Guangli，Zha Jianfeng.Mining Subsidence and Its Prevention[M].Xuzhou:China University of Mining and Technology Press,2012.

[8] 王漢雄，喬景順.Excel VBA在測量數(shù)據(jù)處理中的應用[J].測繪科學，2008(2)：210-212. Wang Hanxiong，Qiao Jingshun.Application of Excel VBA in data processing of engineering surveying[J].Science of Surveying and Mapping，2008(2):210-212.

[9] 劉曉云，張世娟，程傳錄.精密水準測量數(shù)據(jù)處理自動化系統(tǒng)的研究與實現(xiàn)[J].測繪通報，2013(10)：67-69. Liu Xiaoyun，Zhang Shijuan，Cheng Chuanlu.Research and realization of precise leveling data processing auto-system[J].Bulletin of Surveying and Mapping，2013(10):67-69.

[10] 羅剛君，章蘭新，黃朝陽.Excel 2010 VBA編程與實踐[M].北京：電子工業(yè)出版社，2010. Luo Gangjun，Zhang Lanxin，Huang Chaoyang.Excel 2010 VBA Programming and Practice[M].Beijing:Electronics Industry Press,2010.

[11] 劉橋喜，毛善君，馬藹乃，等.煤礦安全地理信息系統(tǒng)的設計與實現(xiàn)[J].測繪通報，2004(2)：60-62. Liu Qiaoxi，Mao Shanjun，Ma Ainai，et al.The design and realization of coal mine safety geographical information system[J].Bulletin of Surveying and Mapping，2004(2):60-62.

[12] 國家煤炭工業(yè)局.建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規(guī)程[M].北京：煤炭工業(yè)出版社，2004. State Bureau of Coal Industry.Buildings，Water，Railway and Coal Pillar Design and Coal Mining Regulations[M].Beijing:China Coal Industry Press,2004.

[13] 胡祝敏，劉 星，艾鴻敏.CAD二次開發(fā)在宗地數(shù)據(jù)處理中的應用[J].測繪工程，2014(5)：62-65. Hu Zhumin，Liu Xing，Ai Hongmin.Application of CAD secondary development to parcel data processing[J].Engineering of Surveying and Mapping,2014(5):62-65.

[14] 王新生，馬海濤，寧化展.基于Java 3D的場地平整土方計算小軟件的研究與實現(xiàn)[J].金屬礦山，2013(6)：49-52. Wang Xinsheng，Ma Haitao，Ning Huazhan.Research and realization of small software for calculation of earthwork in site leveling based on Java 3D[J].Metal Mine，2013(6):49-52.

[15] 王 剛，郭廣禮，李 伶.山區(qū)移動觀測站布設模式與數(shù)據(jù)處理[J].金屬礦山，2011(6)：116-119. Wang Gang，Guo Guangli，Li Ling.Layout mode and data processing of mountain surface movement observation station[J].Metal Mine，2011(6):116-119.

[16] 王 剛.地表移動觀測站數(shù)據(jù)處理模式[J].煤礦安全，2011(11)：129-133. Wang Gang.The data processing model of surface movement observation station[J].Safety in Coal Mine，2011(11):129-133.

[17] 朱黎明，吳 捷，于雪濤，等.礦山開采沉陷地表移動觀測站的設計與實現(xiàn)[J].測繪與空間地理信息，2014(3)：96-99. Zhu Liming，Wu Jie，Yu Xuetao，et al.Design and realization of the observation station of surface movement caused by mining subsidence[J].Geomatics & Spatial Information Technology，2014(3):96-99.

[18] 賈小敏，余學祥，趙 勇.基于VB與CAD的礦區(qū)地表移動變形求解及圖形繪制[J].測繪信息與工程，2009(1)：38-39. Jia Xiaomin,Yu Xuexiang，Zhao Yong.Surface movement deformation calculation and graph drawing based on VB and CAD[J].Journal of Geomatics，2009(1):38-39.

[19] 袁佳佳，戴佳琪，姚光虎.煤首采面地表移動觀測線與觀測方案設計[J].測繪地理信息，2013(5)：34-37. Yuan Jiajia，Dai Jiaqi，Yao Guanghu.Surface movement observation line and the observation design of first working face[J].Journal of Geomatics，2013(5):34-37.

[20] 趙興旺，余學祥，劉 輝.基于MO的地表移動觀測站設計系統(tǒng)開發(fā)[J].遼寧工程技術大學學報，2007(S)：43-46. Zhao Xingwang，Yu Xuexiang，Liu Hui.Design of surface ground movement observation station based on MO[J].Journal of Liaoning Technical University，2007(S):43-46.

(責任編輯 王小兵)

Realization of Aided Design of Surface Movement Observation Station and Data Processing in Mining Area

Fang Qi1,2,3Guo Guangli1,2,3Zhu Xiaojun1,2,3

(1.SchoolofEnvironmentScienceandSpatialInformatics，ChinaUniversityofMiningandTechnology,Xuzhou221116，China;2.NationalAdministrationofSurveying，MappingandGeo-informationKeyLaboratoryforLandEnvironmentandDisasterMonitoring，Xuzhou221116，China；3.KeyLaboratoryofResourceandEnvironmentInformationEngineeringofJiangsuProvince，Xuzhou221116，China)

In order to improve the automation of the surface movement observation station in design and the efficiency of data processing and mapping of movement and deformation curve，based on the programming language of Visual Basic(VB) and Visual Basic for Application (VBA)，the functions such as query of parameters，computation of observation line，observation station design，output of data，data processing and automatic mapping of surface movement observation station are realized by adopting the technology of Access，AutoCAD and Excel.Firstly，the surface and linear route movement observation stations at working face are arranged quickly by using the AutoCAD software；then the observation data are processed automatically and the deformation and movement graphs are mapped based on the Excel software.In addition，the function of axis label of scatter chart with a smooth line is optimized；Finally，both of the functions are integrated into one program to realize the function module integration so as to overcome the insufficiency of single function，low degree of automation and difficulty in modification.The research results show that，the program written in this paper can satisfy the needs for variable observation station design and the special requirement of mapping the movement and deformation curve.Besides that，this program is conducive to improve the automation degree of mining surface movement observation and it can also provide some reference for establishing the surface movement observation system.

Surface movement observation station，Data processing，Movement deformation curve，VB，VBA，Surface movement observation system

2015-03-04

“十二五”國家科技支撐計劃項目(編號：2012BAB13B03)，國家自然科學基金青年基金項目(編號：41104011)，江蘇高校優(yōu)勢學科建設工程項目(編號：SZBF2011-6-B35)，江蘇省資源環(huán)境信息工程重點實驗室基金項目(編號：JS201309)。

方 齊(1989—)，男，碩士研究生。 通訊作者 郭廣禮(1965—)，男，教授，博士，博士研究生導師。

TD173

A

1001-1250(2015)-05-129-06