地表移動(dòng)觀測(cè)站數(shù)據(jù)處理新方法研究＊

王 剛 郭廣禮 王 磊

(1)國(guó)土環(huán)境與災(zāi)害監(jiān)測(cè)國(guó)家測(cè)繪局重點(diǎn)實(shí)驗(yàn)室(中國(guó)礦業(yè)大學(xué)),徐州 221116 2)江蘇省資源環(huán)境信息工程重點(diǎn)實(shí)驗(yàn)室 (中國(guó)礦業(yè)大學(xué)),徐州 221116)

地表移動(dòng)觀測(cè)站數(shù)據(jù)處理新方法研究＊

王 剛1,2)郭廣禮1,2)王 磊1,2)

(1)國(guó)土環(huán)境與災(zāi)害監(jiān)測(cè)國(guó)家測(cè)繪局重點(diǎn)實(shí)驗(yàn)室(中國(guó)礦業(yè)大學(xué)),徐州 221116 2)江蘇省資源環(huán)境信息工程重點(diǎn)實(shí)驗(yàn)室 (中國(guó)礦業(yè)大學(xué)),徐州 221116)

建立地表移動(dòng)觀測(cè)站是研究開(kāi)采沉陷和建筑物破壞規(guī)律的最可靠手段。目前,實(shí)地觀測(cè)數(shù)據(jù)在后期的處理過(guò)程中大多都是采用傳統(tǒng)的支距法等方法,將觀測(cè)值改正到工作面的走向和傾向方向上,使得觀測(cè)站布設(shè)與數(shù)據(jù)處理具有較大的局限性,地表移動(dòng)和變形計(jì)算精度也受到較大的影響。針對(duì)傳統(tǒng)的地表移動(dòng)觀測(cè)站數(shù)據(jù)處理存在的問(wèn)題,結(jié)合工程實(shí)例介紹一種地表移動(dòng)觀測(cè)站數(shù)據(jù)處理新方法。

地表移動(dòng);觀測(cè)站;數(shù)據(jù)處理;坐標(biāo)轉(zhuǎn)換;精度

1 引言

研究開(kāi)采引起地表移動(dòng)變形的方法很多,如實(shí)地觀測(cè)法、理論模擬法和相似材料模型法等。但目前公認(rèn)最可靠最有效的認(rèn)識(shí)巖層與地表移動(dòng)規(guī)律的方法還是實(shí)地觀測(cè)法,即按一定要求在受開(kāi)采影響的地表、巖層內(nèi)部、建筑物上設(shè)置一系列的觀測(cè)點(diǎn),在采動(dòng)過(guò)程中,定期獲取這些點(diǎn)的坐標(biāo),以確定各觀測(cè)點(diǎn)位置的變化和觀測(cè)點(diǎn)間的相對(duì)移動(dòng),從而找出地表和巖層移動(dòng)變形的規(guī)律。因此,除了觀測(cè)站的設(shè)計(jì)方案必須合理外,實(shí)測(cè)資料數(shù)據(jù)處理的準(zhǔn)確性也十分重要,其結(jié)果準(zhǔn)確與否直接關(guān)系到所得規(guī)律的可信度。如果處理方法恰當(dāng),所得結(jié)論就能正確反映某種地質(zhì)采礦條件下的地表移動(dòng)和變形規(guī)律;反之,如果處理方法不當(dāng)或者出現(xiàn)錯(cuò)誤,研究就會(huì)失去意義。

因此,在地表移動(dòng)觀測(cè)站的布設(shè)過(guò)程中,其正確的設(shè)計(jì)方案和合理的數(shù)據(jù)處理模式是保證研究意義能夠?qū)崿F(xiàn)的關(guān)鍵所在。鑒于此,本文結(jié)合具體的工程實(shí)例,研究了一種地表移動(dòng)觀測(cè)站的數(shù)據(jù)處理新方法。

2 傳統(tǒng)觀測(cè)站布置及數(shù)據(jù)處理存在的問(wèn)題

地表移動(dòng)觀測(cè)站可分為網(wǎng)狀觀測(cè)站和斷面線狀觀測(cè)站。但有時(shí)限于地形條件,在布設(shè)斷面線狀觀測(cè)站不便時(shí),在山區(qū)也可以布置網(wǎng)狀觀測(cè)站[1],一般情況下主要還是布設(shè)為斷面線狀觀測(cè)站。傳統(tǒng)的地表移動(dòng)觀測(cè)站要求設(shè)計(jì)的觀測(cè)線必須嚴(yán)格位于地表下沉盆地主斷面線上,利用經(jīng)緯儀、水準(zhǔn)儀、鋼尺和支距尺等進(jìn)行測(cè)量,其優(yōu)點(diǎn)是所用儀器價(jià)格低廉,使用也較方便,一般能滿足測(cè)量精度的要求。但由于水平移動(dòng)觀測(cè)數(shù)據(jù)只有各測(cè)點(diǎn)沿觀測(cè)線方向的水平移動(dòng)值,垂直觀測(cè)線方向的水平移動(dòng)值一般不進(jìn)行觀測(cè),從而造成了以下問(wèn)題：

1)外業(yè)作業(yè)勞動(dòng)強(qiáng)度大,內(nèi)業(yè)計(jì)算復(fù)雜,精度較低。特別對(duì)于水平移動(dòng)測(cè)量,由于鋼尺量邊的誤差隨著邊長(zhǎng)的增大而逐漸積累,在測(cè)線較長(zhǎng)時(shí),水平移動(dòng)測(cè)量的精度就會(huì)很低,甚至達(dá)不到要求的精度。

2)由于受各種因素的干擾,地表移動(dòng)觀測(cè)站不可能按理論上布置成一條直線而是沿某一直線左右擺動(dòng)[2,3],從而使得水平移動(dòng)值存在一定誤差 (圖1)。

圖1 傳統(tǒng)觀測(cè)站布置Fig.1 Layout of traditional observatory

各觀測(cè)點(diǎn)的水平移動(dòng)值一般是按以下方法求出的：假設(shè)A為控制點(diǎn),其末次觀測(cè)與首次觀測(cè)相比, AB伸長(zhǎng)量為ΔLAB,BC伸長(zhǎng)量為ΔLBC,CD伸長(zhǎng)量為ΔLCD等。則 B、C、D3點(diǎn)沿傾斜方向的水平移動(dòng)值分別為ΔLAB、ΔLAB+ΔLBC、ΔLAB+ΔLBC+ΔLCD。

這樣求出的水平移動(dòng)值存在以下誤差：

首先,ΔLAB、ΔLBC、ΔLCD是 3個(gè)不同方向的矢量,不能求代數(shù)和,相鄰線段方向變化越大其代數(shù)和的準(zhǔn)確性越差,而且觀測(cè)點(diǎn)距離控制點(diǎn)越遠(yuǎn),求出的位移值累計(jì)誤差值越大。

其次,不在主斷面線的各點(diǎn)的實(shí)際位移值與主斷面線方向存在一個(gè)夾角,因此實(shí)際位移必須分解為沿主斷面線方向的一個(gè)分量和垂直主斷面線方向的一個(gè)分量,而在計(jì)算時(shí)經(jīng)常把非主斷面線方向的位移看作沿主斷面線方向的位移。

3)由于地形、地貌等原因,觀測(cè)線可能無(wú)法沿走向和傾向布置,而必須布置成折線 (圖 2)。在傳統(tǒng)觀測(cè)方法中,對(duì)于AB線上的各觀測(cè)點(diǎn),只觀測(cè)順著AB方向位移,而不觀測(cè)垂直于 AB方向的位移,對(duì)于BC線上的各觀測(cè)點(diǎn)同樣如此。即在實(shí)際數(shù)據(jù)的處理過(guò)程中,常把沿觀測(cè)線方向的值當(dāng)作總位移,從而無(wú)法知曉沿?cái)嗝婢€方向的實(shí)際位移值。由于AB線上的觀測(cè)點(diǎn)與 BC線上的觀測(cè)點(diǎn)的觀測(cè)值方向不同,并且無(wú)法分解,而傳統(tǒng)的求參程序是采用主斷面的移動(dòng)進(jìn)行擬合的,因此,無(wú)論按哪個(gè)方向擬合都會(huì)產(chǎn)生誤差。

圖2 特殊情況下觀測(cè)站的布置Fig.2 Layout of observatory under special circumstance

4)測(cè)量工作量大,數(shù)據(jù)處理復(fù)雜,尤其在進(jìn)行支距測(cè)量和支距改正時(shí)更是如此。支距改正的目的是將相鄰觀測(cè)點(diǎn)間的距離改正到成觀測(cè)線方向上,如圖 3,A、B、C為觀測(cè)線上相鄰的 3個(gè)點(diǎn),點(diǎn)間距分別為 L1、L2,設(shè) 3個(gè)點(diǎn)在觀測(cè)線方向的距離為 L11、L22,支距改正ΔL按式 (1)計(jì)算：

圖3 支距改正Fig.3 Distance correct

3 地表移動(dòng)觀測(cè)數(shù)據(jù)處理新方法

為了克服傳統(tǒng)觀測(cè)站數(shù)據(jù)在處理過(guò)程中的缺陷,許多學(xué)者曾進(jìn)行了大量的研究。李春意等[4]改進(jìn)了觀測(cè)站的數(shù)據(jù)處理方法：先求出各觀測(cè)點(diǎn)的坐標(biāo)變化量,再將坐標(biāo)變化量投影到觀測(cè)線上(圖 4)。周萬(wàn)茂[5]提出兩種方法來(lái)改進(jìn)數(shù)據(jù)處理的復(fù)雜程度：一是將傳統(tǒng)的相鄰觀測(cè)點(diǎn)間距的觀測(cè)法改為按大坐標(biāo)或礦區(qū)坐標(biāo)或其他相對(duì)坐標(biāo)觀測(cè)各觀測(cè)點(diǎn)坐標(biāo)的變化,二是把沿觀測(cè)線方向的位移改為沿兩個(gè)坐標(biāo)方向的位移。

圖4 觀測(cè)點(diǎn)數(shù)據(jù)處理Fig.4 Data processing of observation point

但這種模式對(duì)于地表移動(dòng)觀測(cè)站的數(shù)據(jù)處理計(jì)算量較大。為此,本文試圖在實(shí)際應(yīng)用中探索出另一種較為簡(jiǎn)潔方便的數(shù)據(jù)處理方法。

觀測(cè)數(shù)據(jù)經(jīng)過(guò)整理改正后,便可計(jì)算觀測(cè)線上各測(cè)點(diǎn)和各測(cè)點(diǎn)間的移動(dòng)和變形。移動(dòng)和變形計(jì)算主要包括：各測(cè)點(diǎn)的下沉和水平移動(dòng),相鄰兩側(cè)點(diǎn)間的傾斜和水平變形,相鄰兩線段 (或相鄰 3點(diǎn))的曲率變形,觀測(cè)點(diǎn)的下沉速度等。在計(jì)算之前,先將用GPS測(cè)得各點(diǎn)的WGS84坐標(biāo)轉(zhuǎn)換為BJ54坐標(biāo)系坐標(biāo),再將BJ54坐標(biāo)系坐標(biāo)轉(zhuǎn)換為礦區(qū)工作面坐標(biāo)系統(tǒng)坐標(biāo)。具體轉(zhuǎn)換過(guò)程如下(圖 5)。

圖5 BJ54坐標(biāo)系和礦區(qū)工作面坐標(biāo)系的轉(zhuǎn)換Fig.5 Conversion bet ween BJ54 coordinate system and mining face coordinate system

用 GPS和全站儀分別精確測(cè)出控制點(diǎn)A、B、C、D的WGS84坐標(biāo)和高斯坐標(biāo),利用坐標(biāo)轉(zhuǎn)換工具求出兩套坐標(biāo)系統(tǒng)的轉(zhuǎn)換參數(shù),進(jìn)而實(shí)現(xiàn)其他各點(diǎn)的坐標(biāo)轉(zhuǎn)換。在圖上量出控制點(diǎn) A、B、C、D的礦區(qū)工作面坐標(biāo)系統(tǒng)坐標(biāo),利用坐標(biāo)轉(zhuǎn)換工具求出BJ54坐標(biāo)系和礦區(qū)工作面坐標(biāo)系的轉(zhuǎn)換參數(shù),進(jìn)而將其他各點(diǎn)的坐標(biāo)轉(zhuǎn)換為礦區(qū)工作面坐標(biāo)系統(tǒng)坐標(biāo)。最后進(jìn)行各移動(dòng)和變形的計(jì)算,公式如下：

1)m次觀測(cè)時(shí) n點(diǎn)的下沉值：

式中,hn0、hnm為首次和m次觀測(cè)時(shí) n點(diǎn)的高程。

也可用下列公式計(jì)算m次觀測(cè)時(shí) n點(diǎn)的下沉：

式中,zn0、znm為首次和 m次觀測(cè)時(shí) n點(diǎn)在礦區(qū)坐標(biāo)系下的 z坐標(biāo)。

2)相鄰兩點(diǎn)的傾斜

x軸(傾向)方向傾斜分量：

y軸 (走向)方向傾斜分量：

總的傾斜量：

式中,Δxn～n+1為 n號(hào)點(diǎn)和 n+1號(hào)點(diǎn)的 x坐標(biāo)差, Δyn～n+1為 n號(hào)點(diǎn)和 n+1號(hào)點(diǎn)的 y坐標(biāo)差,wn+1、wn為分別表示 n+1號(hào)點(diǎn)和 n號(hào)點(diǎn)的下沉值。

3)n號(hào)點(diǎn)附近的曲率

式中,in+1～n、in～n-1為 n+1號(hào)點(diǎn)至 n號(hào)點(diǎn)和 n號(hào)點(diǎn)至 n-1號(hào)點(diǎn)的傾斜;Δxn+1～n、Δxn～n-1為 n+1號(hào)點(diǎn)至n號(hào)點(diǎn)和 n號(hào)點(diǎn)至 n-1號(hào)點(diǎn)的 x坐標(biāo)差; Δyn+1～n、Δyn～n-1為 n+1號(hào)點(diǎn)至 n號(hào)點(diǎn)和 n號(hào)點(diǎn)至 n -1號(hào)點(diǎn)的 y坐標(biāo)差。

4)n號(hào)點(diǎn)的水平移動(dòng)

x軸方向水平移動(dòng)量：

y軸方向水平移動(dòng)量：

總的水平移動(dòng)量：

式中,xnm、xn0為 m次觀測(cè)時(shí)和首次觀測(cè)時(shí) n號(hào)點(diǎn)的x坐標(biāo);ynm、yn0為 m次觀測(cè)時(shí)和首次觀測(cè)時(shí) n號(hào)點(diǎn)的y坐標(biāo)。

5)n號(hào)點(diǎn)至 n+1號(hào)點(diǎn)間的水平變形

x軸方向水平變形量：

y軸方向水平變形量：

總的水平變形量：

式中,(xn+1～n)m、(xn+1～n)0為 n+1號(hào)點(diǎn)和 n號(hào)點(diǎn)在m次觀測(cè)時(shí)和首次觀測(cè)時(shí)的 x坐標(biāo);(yn+1～n)m、 (yn+1～n)0為 n+1號(hào)點(diǎn)和 n號(hào)點(diǎn)在m次觀測(cè)時(shí)和首次觀測(cè)時(shí)的 y坐標(biāo)。

4 工程實(shí)例

4.1 礦區(qū)概況

濟(jì)寧三號(hào)煤礦位于山東省西南部濟(jì)寧市南郊,距濟(jì)寧市 14 km。礦井設(shè)計(jì)年產(chǎn)量 500萬(wàn)噸,是我國(guó)特大型現(xiàn)代化礦井。井田范圍內(nèi)地表面以農(nóng)田和村莊為主。63下 04矸石充填試采區(qū)位于濟(jì)三礦工業(yè)廣場(chǎng)西北部,辛店村以東,南陽(yáng)湖堤從工作面南部穿過(guò),地表大部分為農(nóng)田,工作面開(kāi)采南陽(yáng)湖堤和辛店村的保護(hù)煤柱。地面標(biāo)高 +33.0～+33.3 m,平均 33.14 m,工作面標(biāo)高 -639.0～-663.0 m,平均-650 m。該面東臨 63下 03工作面采空區(qū),西臨尚未回采 63下 05工作面,北臨 63下 04工作面采空區(qū),初步設(shè)計(jì)停采線向南距六采區(qū)輔助運(yùn)巷巷中107.0 m。

矸石充填試采區(qū)域煤巖層總體趨勢(shì)呈現(xiàn)南高北低的單斜構(gòu)造,局部伴生寬緩的波狀起伏。共發(fā)育斷層 10條,落差大于 3.0 m的斷層兩條,最大落差6.5 m,斷層在采區(qū)內(nèi)延展長(zhǎng)度 8.0～523.0 m。根據(jù)井巷工程揭露情況,走向以南北為主,延展長(zhǎng)度一般較長(zhǎng),落差變化較小,以張性正斷層為主。

矸石充填試采區(qū)域位于辛店村保護(hù)煤柱東北部,其正上方地表為辛店村農(nóng)田和積水塌陷區(qū),西距辛店村主體建筑群約 480 m。

4.2 應(yīng)用成果及分析

為了使 GPS的測(cè)量成果轉(zhuǎn)換為礦區(qū)坐標(biāo)系統(tǒng),必須先計(jì)算出坐標(biāo)轉(zhuǎn)換參數(shù),而求取坐標(biāo)轉(zhuǎn)換參數(shù)必須己知至少 3個(gè)高精度控制點(diǎn)的礦區(qū)坐標(biāo)值 x、y、h。根據(jù)已知控制點(diǎn)均勻分布在測(cè)區(qū)的原則,本次工程中礦區(qū)內(nèi)選擇 4個(gè)己知的靜態(tài) GPS控制點(diǎn),分別為A、B、C、D(圖 6)。

圖6 濟(jì)三煤礦地表移動(dòng)觀測(cè)站沉降觀測(cè)點(diǎn)布設(shè)圖Fig.6 Layout of subsidence observation spots at surface movement observation station in Jisan mine

利用隨機(jī)軟件求出坐標(biāo)轉(zhuǎn)換參數(shù)為：比例系數(shù)S=-1 000 016,旋轉(zhuǎn)度 JR=-0.101 449,高程參數(shù)a= -0.000 047 331,b = 0.000 023 634,c= 1.373 809 179;然后利用此參數(shù)將礦區(qū)其他測(cè)點(diǎn)轉(zhuǎn)換至BJ54坐標(biāo)系;最后利用控制點(diǎn) A、B、C、D求出BJ54坐標(biāo)系和該礦區(qū)新建坐標(biāo)系的轉(zhuǎn)換參數(shù) (表1),進(jìn)而將其他測(cè)點(diǎn)的 BJ54坐標(biāo)轉(zhuǎn)換為礦區(qū)坐標(biāo)。最后利用式(1)～(12)分別求出各測(cè)點(diǎn)的移動(dòng)變形值。該礦區(qū)坐標(biāo)系以 83201工作面上順槽上一點(diǎn) o為坐標(biāo)原點(diǎn),傾向方向?yàn)?x軸,走向方向?yàn)?y軸 (圖7)。

圖7 礦區(qū)坐標(biāo)系Fig.7 Coordinate system ofmining area

限于篇幅,在此只隨機(jī)抽取其中的 3個(gè)觀測(cè)點(diǎn)進(jìn)行驗(yàn)算。表 1表示控制點(diǎn)A、B、C、D4點(diǎn)的 BJ54坐標(biāo)系坐標(biāo)和礦區(qū)坐標(biāo)系坐標(biāo)以及通過(guò)坐標(biāo)轉(zhuǎn)換參數(shù)軟件利用四參數(shù)法計(jì)算出來(lái)的兩坐標(biāo)系統(tǒng)的轉(zhuǎn)換參數(shù),四參數(shù)法的轉(zhuǎn)換模型如式(14)所示。

為了研究開(kāi)采后觀測(cè)點(diǎn) a、b、c的實(shí)測(cè)變化和用本文所述新方法計(jì)算所得變化間的關(guān)系,將二者進(jìn)行了對(duì)比。表 3表示地表觀測(cè)點(diǎn) a、b、c受開(kāi)采影響后的BJ54坐標(biāo)和由表 1中轉(zhuǎn)換參數(shù)計(jì)算得到的相應(yīng)礦區(qū)坐標(biāo),表 4表示由實(shí)測(cè)所得各觀測(cè)點(diǎn)的位移和由本文所述新方法計(jì)算所得各觀測(cè)點(diǎn)的位移,實(shí)測(cè)位移值由圖7量得。

表1 控制點(diǎn)的 BJ54坐標(biāo)和礦區(qū)坐標(biāo)(單位:m)Tab.1 BJ54 coordinates and m ine surface coordinates of control poi nts(un it:m)

表2 開(kāi)采前各觀測(cè)點(diǎn)的 BJ54坐標(biāo)和礦區(qū)坐標(biāo)(單位:m)Tab.2 BJ54 coordinates and m ine surface coordinates of observation points before m in i ng(un it:m)

表3 開(kāi)采后各觀測(cè)點(diǎn)的 BJ54坐標(biāo)和礦區(qū)坐標(biāo)(單位:m)Tab.3 BJ54 coordinates and m ine surface coordi nates of observation points after m i n ing(un it:m)

表4 實(shí)測(cè)變化值和計(jì)算變化值(單位:m)Tab 4 M easured and calculated changes(un it:m)

5 結(jié)論

1)該新方法不需支距測(cè)量、鋼尺量邊,克服了傳統(tǒng)地表移動(dòng)觀測(cè)站數(shù)據(jù)處理過(guò)程中存在的弊端,便于用計(jì)算機(jī)進(jìn)行數(shù)據(jù)處理,從而可大大提高數(shù)據(jù)處理的精度和工作效率。

2)GPS的測(cè)量成果只需要進(jìn)行兩次坐標(biāo)轉(zhuǎn)換即可直接計(jì)算出地表各移動(dòng)變形指標(biāo),過(guò)程簡(jiǎn)單。

3)雖然在兩次坐標(biāo)轉(zhuǎn)化中精度會(huì)有所降低,但從實(shí)例分析可知,坐標(biāo)轉(zhuǎn)換中精度損失很少,計(jì)算精度達(dá)到毫米級(jí),完全可以滿足對(duì)開(kāi)采沉陷分析的精度要求。

4)由于觀測(cè)手段和計(jì)算手段的局限性,傳統(tǒng)地表移動(dòng)觀測(cè)站在數(shù)據(jù)處理方法上存在種種問(wèn)題,但隨著先進(jìn)觀測(cè)儀器的使用和計(jì)算機(jī)的普及,本文提出的方法將具有可推廣性。

1 原東方.山區(qū)地表移動(dòng)觀測(cè)站數(shù)據(jù)處理方法[J].焦作工學(xué)院學(xué)報(bào),1996,15(2)：14-19.(Yuan Dongfang.Resarch on data processingmethod ofmountain surfacemovementobservation[J].Journal of Henan Polytechnic University(Natural Science),1996,15(2)：14-19)

2 李萬(wàn)武.地表移動(dòng)觀測(cè)站的建立 [J].煤炭技術(shù),2002 (6)：69-70.(LiWanwu.Building of the surface migrating observation station[J].Coal Technology,2002(6)：67-70)

3 趙立武,等.非正規(guī)地表移動(dòng)觀測(cè)站數(shù)據(jù)處理的一種新方法[J].煤炭科技,1999,(2)：29-30.(Zhao Liwu,et al. A new kind ofmethod of infor mal surface movement observatory data processing[J].Coal Technology,1999,(2)：29-30)

4 李春意,郭增長(zhǎng),楊福芹.地表移動(dòng)觀測(cè)站數(shù)據(jù)處理方法探討[J].河南理工大學(xué)學(xué)報(bào),2006,25(6)：481-483.(Li Chunyi,Guo Zengchang and Yang Fuqin.Data procession methods discussion of land surfacemovementobservation station[J].Journal of Henan Polytechnic University,2006,25 (6)：481-483)

5 周萬(wàn)茂.地表移動(dòng)觀測(cè)站布設(shè)及數(shù)據(jù)處理方法簡(jiǎn)析[J].煤礦開(kāi)采,2000,1.(Zhou Wanmao.A kind new method of surface movement observatory arrangement and data processing[J].Coal Mining Technology,2006,25(6)：481-483.)

A NEW M ETHOD FOR DATA PROCESSING OF LAND SURFACE MOVEM ENT OBSERVATI ON STATI ON

Wang Gang1,2),Guo Guangli1,2)andWangLei1,2)

(1)Key Laboratory for Resources Environm ent and D isasterM onitoring of SBSM(CUM T),Xuzhou 221116 2)M ain Laboratory of Resource Environm ent Infor m ation of Jiangsu(CUM T),Xuzhou 221116)

Establishingmobile stations is the most reliable way for researching surface subsidence and building damage.At present,the processing in later stage of field observations is mostly with offset method or other traditionalmethods to converse the observations to the strike direction and tendency direction,thus station’s layout and data processing are subject to greater li mitations,the calculation accuraciesofmovement and deformation of the surface are significantly affected.For the Shortcomings of traditional data processing of observation station,a new method of data processing of observation station is introduced with practical engineering example.

land surface movement;observation station;data processing;coordinate transformation;accuracy

1671-5942(2011)04-0127-05

2010-11-23

國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(50834004);國(guó)家公益性行業(yè)基金(200811050)

王剛,男,1986年生,碩士研究生,主要研究方向?yàn)殚_(kāi)采沉陷與“三下采煤”.E-mail：wanggangcumt@163.com

P207

A