周 振，李 青，池金谷，李 雄，杜通波

（中國計量學(xué)院 機(jī)電工程學(xué)院，浙江 杭州 310018）

滑坡一般指由于受到地質(zhì)結(jié)構(gòu)、地表或地下水活動、人工切坡等因素的影響，斜坡上的部分巖體和土體沿地層軟弱面（或軟弱帶）在重力作用下向下滑動的不良地質(zhì)現(xiàn)象［1-2］.國際上采用 Varnes提出的五種運(yùn)動類型來對滑坡進(jìn)行分類：崩塌，傾倒，滑動，側(cè)向擴(kuò)展和流動［3-4］.美國地質(zhì)調(diào)查局（USGS）于2004年進(jìn)一步細(xì)分，包括：旋轉(zhuǎn)型滑坡，平移型滑坡，塊體滑坡，巖崩，彎折傾倒，碎屑流，碎屑崩塌，泥土流動，蠕滑和側(cè)向擴(kuò)展.

滑坡的監(jiān)測內(nèi)容包括地表位移、地下位移、地下水位、地下推力、土壤濕度、位移速率、地下水壓力、地應(yīng)力等，其中地下位移的監(jiān)測是重點(diǎn)也是難點(diǎn).國內(nèi)外對于坡體深部位移監(jiān)測主要利用埋設(shè)在地質(zhì)鉆孔中的測斜儀［5-6］、多 點(diǎn) 位 移 計［7-9］和TDR同軸電纜［10-13］來獲取數(shù)據(jù).鉆孔測斜儀不能測量滑坡體不同深度的垂直位移；多點(diǎn)位移計只能測量一維位移；TDR監(jiān)測技術(shù)很難確定滑坡滑動的方向，且同軸電纜的伸縮性很差；三種主流測量方法都不能測量出地下位移的三維變化，而只有三維變化的測量才能真正反應(yīng)出地下巖土的實際變化情況.鑒于此，筆者提出一種基于霍爾效應(yīng)的地下位移三維測量方法.

1 工作原理

旋轉(zhuǎn)型滑坡、平移型滑坡和塊體滑移是常見的滑坡類型，旋轉(zhuǎn)型滑坡示意圖如下：

圖1 旋轉(zhuǎn)型滑坡Figure 1 Rotational landslide

旋轉(zhuǎn)型滑坡是指沿著向上凹的弧形破裂面滑動的滑坡，滑坡移動只沿內(nèi)部滑面發(fā)生，坡體裂縫在平面上呈環(huán)狀且凹向移動方向，其下伏破裂面連同外露破裂壁呈勺形，滑坡的坡體移動時近似繞平行于斜坡的軸轉(zhuǎn)動［14］.

平移型滑坡和塊體滑移的示意圖如圖2.

圖2 平移型滑坡和塊體滑移Figure 2 Translational landslide and block slide

對旋轉(zhuǎn)滑坡的地下位移測量需要同時測量出水平及垂直位移，這就要求測量方法能實現(xiàn)三維測量.本文提出的基于霍爾效應(yīng)的地下位移三維測量方法，其工作原理如圖3.

圖3中的測量單元均為圓柱體，內(nèi)置傳感器及電路.各個測量單元通過電源線、信號線連成一串與主機(jī)相連，且各個測量單元之間用高彈性連接物相連.在實際安裝過程中，保持各個測量單元的上下平面平行，并找到地下巖層的不動層，將一個測量單元埋入其中作為基準(zhǔn)，建立三維坐標(biāo)系.

當(dāng)最復(fù)雜的多級旋轉(zhuǎn)型滑坡發(fā)生時，每一個滑坡坡體移動近似繞平行于斜坡的軸轉(zhuǎn)動，且每一個滑坡坡體基本保持原樣.這樣，坡體的運(yùn)動將改變各個測量單元之間的相對位置，且相鄰測量單元之間的上下表面基本保持平行.如圖4.

只要能測量出各個測量單元之間的相對位置變化，并以埋入不動層中的測量單元為基準(zhǔn)，就可以測量出整個地下巖土的三維變化情況.類似的測量思想可以在很多文獻(xiàn)中得到佐證，并有大量的實際應(yīng)用項目支持，筆者所在的科研團(tuán)隊對此相關(guān)領(lǐng)域已有多年的研究［15-16］.能否完成三維測量，關(guān)鍵在于如何測量出各個測量單元之間的相對位置變化.基于霍爾效應(yīng)設(shè)計集成式傳感器并配合相應(yīng)的算法，就可實現(xiàn)這一目的.

圖4 滑坡示意圖Figure 4 Schematic plan of landslide

霍爾效應(yīng)：當(dāng)電流垂直于外磁場通過導(dǎo)體（或半導(dǎo)體）時，在導(dǎo)體（或半導(dǎo)體）的垂直于磁場和電流方向的兩個端面之間會出現(xiàn)電勢差.可用如下公式表示：

式（1）中，RH為霍爾系數(shù)；I為通過的電流；B為垂直于I的磁感應(yīng)強(qiáng)度；d為導(dǎo)體的厚度.當(dāng)電流I和霍爾片的厚度d保持不變時，霍爾電勢差UH與磁感應(yīng)強(qiáng)度B成正比.通過電勢差UH即可反映磁感應(yīng)強(qiáng)度的變化［17］.

在本方案中，磁場來源為圓柱形釹鐵硼永磁體.當(dāng)磁鋼與霍爾傳感器的相對位置發(fā)生變化時，霍爾傳感器的輸出電壓就能指示出具體的變化情況.

2 設(shè)計方案

本文提出的測量方法，核心內(nèi)容是傳感器的設(shè)計以及所采用的測量算法.

什么是發(fā)展？列寧曾認(rèn)為，有兩種基本的(或可能的？或兩種在歷史上常見的？)發(fā)展(進(jìn)化)觀點(diǎn)：一種觀點(diǎn)認(rèn)為發(fā)展就是減少和增加，是重復(fù)，另一種觀點(diǎn)認(rèn)為發(fā)展就是對立面的統(tǒng)一。也有學(xué)者從哲學(xué)視角分析，認(rèn)為發(fā)展是指事物由小到大、由簡單到復(fù)雜、由低級到高級的變化。還有人認(rèn)為發(fā)展指變化的趨勢。一般說來，發(fā)展的內(nèi)涵包括自然的演化和人類社會的進(jìn)步，是指事物自身的矛盾運(yùn)動和空間相互作用和時間序列的永恒歷史運(yùn)動。發(fā)展是人類社會的永恒主題，也是人類文明演進(jìn)的邏輯起點(diǎn)。

圖5 傳感器結(jié)構(gòu)Figure 5 Sensor structure

圖5中，傳感器分為兩個部分.第一部分是空心圓柱體V1，一塊圓柱形釹鐵硼永久磁鋼M（φ5×h18，單位mm）固定在其下端面面板S1上.第二部分是空心圓柱體V2，在距離其上端h2（13.5mm）的面板S2上，是由三個水平緊密貼合的霍爾傳感器組成的陣列，如圖6.霍爾傳感器陶瓷片的中心點(diǎn)可感受磁感應(yīng)強(qiáng)度，三個霍爾傳感器陶瓷片的中心點(diǎn)構(gòu)成一個等邊三角形的三個頂點(diǎn)A、B、C，邊長L為12mm.

圖6 面板S2俯視圖Figure 6 Plan view of panel S2

這里將空心圓柱體V1與V2相隔h1，是為了便于觀察.在實際測量中，只需將兩者貼合即可.在進(jìn)行巖土地下三維位移測量時，將空心圓柱體V1埋入巖土滑動層，V2則埋入巖土不動層.當(dāng)?shù)叵峦翆影l(fā)生相對位移時，就可以測量土層的三維位移情況.如果增加圓柱體V1與V2的數(shù)量，連成葫蘆串，就可測量大厚度土層的三維位移情況.

進(jìn)行如圖7的磁鋼位移標(biāo)定實驗.為配合實驗，搭建如圖8的測量系統(tǒng).

圖7 位移實驗Figure 7 Displacement experiment

圖8 測量系統(tǒng)Figure 8 Measurement system

本文設(shè)計的標(biāo)定和測試實驗是在7SC405高精度五軸運(yùn)動臺上進(jìn)行的.將霍爾傳感器的位置固定，并將磁鋼固定在7SC405運(yùn)動臺的運(yùn)動臂（非金屬材質(zhì)）上.通過上位機(jī)控制運(yùn)動臺的動作，精確改變磁鋼與傳感器之間的相對位置.圖7中，以霍爾傳感器陶瓷片中心點(diǎn)正上方h4（13.5mm）處為坐標(biāo)原點(diǎn)，建立xz直角坐標(biāo)系.z軸垂直于霍爾傳感器所在平面，且經(jīng)過陶瓷片中心點(diǎn).x軸平行于霍爾傳感器所在平面.

首先，讓圓柱形磁鋼下端面的中心點(diǎn)與坐標(biāo)系原點(diǎn)重合，然后使之水平向右移動.每次移動0.2mm，采集霍爾傳感器輸出的電壓值，一共移動100次.然后，讓磁鋼回到坐標(biāo)系原點(diǎn)，并增加高度0.2mm，重復(fù)同樣的實驗.接著，重復(fù)前兩步的實驗.整個實驗將得到一個100行100列的二維數(shù)組.利用MATLAB進(jìn)行三維繪圖，如圖9.

圖9 電壓采集分布圖Figure 9 Distribution map of sampled voltage

當(dāng)磁鋼逐漸遠(yuǎn)離霍爾傳感器，無論是沿著x軸方向或是z軸方向，傳感器輸出電壓都將相應(yīng)減小.在x軸及z軸坐標(biāo)最小處，傳感器輸出電壓達(dá)到最大值.相反，在x軸及z軸坐標(biāo)最大處傳感器輸出電壓最小，因為此時霍爾傳感器與釹鐵硼永磁體之間的相對距離最大.對圖9進(jìn)一步處理，利用MATLAB繪制等高線，每一條等高線代表相同的電壓值，且相鄰兩條等高線之間的差值是相同的，如圖10.圖10中可以看到，等高線基本上是單調(diào)遞減的.

圖10 15條等高線分布圖Figure 10 Distribution map of 15contour lines

回顧之前的位移標(biāo)定實驗，可以得出這樣的結(jié)論：當(dāng)磁鋼在一條等高線上移動時，霍爾傳感器輸出的電壓將不會改變.由于磁鋼是圓柱體，且沿著中軸線被磁化，那么其磁場應(yīng)該是關(guān)于它的中軸線對稱分布的.將等高線繞著z軸旋轉(zhuǎn)一圈，可以得到一個三維等值曲面.當(dāng)磁鋼在這個等值曲面上移動時，霍爾傳感器輸出的電壓值不變.再看之前的傳感器俯視圖（圖6），當(dāng)磁鋼處在霍爾傳感器陣列所在平面的上方時，三個霍爾傳感器將輸出三個電壓值.根據(jù)之前所述的磁鋼位移標(biāo)定實驗，依據(jù)每一個電壓值可以得到相應(yīng)的等值曲面.兩個等值曲面相交得到一條曲線，這條曲線與第三個等值曲面相交就能得到一個唯一交點(diǎn)，交點(diǎn)的三維坐標(biāo)就是磁鋼下端面中心點(diǎn)的坐標(biāo).

以上闡述了此測量方法在理論上的可行性.理想的交點(diǎn)計算方法是利用三條空間曲線直接計算求取交點(diǎn)，這三條空間曲線分別對應(yīng)于三個等值曲面.但從MATLAB中能提取的僅僅是等高線上離散點(diǎn)的坐標(biāo)，并不是等高線的函數(shù)表達(dá)式.由于對離散點(diǎn)采用曲線擬合的方法獲得函數(shù)表達(dá)式，其誤差較大，故不妨換一個角度思考問題：假設(shè)已經(jīng)得到能夠準(zhǔn)確表達(dá)等高線的函數(shù)表達(dá)式，用一個平行于xy平面的切面橫切上文中提到的三個等值曲面，那么一定能得到三個圓.設(shè)三個等值曲面交點(diǎn)的坐標(biāo)是（x0，y0，z0），則必定會有一個切面z＝z0，使得此切面橫切出來的三個圓能交與一點(diǎn)，這個點(diǎn)就是所要求的交點(diǎn).前文述及構(gòu)成一條等高線的是一簇離散點(diǎn)，那么將相鄰的兩個離散點(diǎn)全部用直線相連就能得到一條折線.將此折線同樣繞著z軸旋轉(zhuǎn)一圈，也能得到一個等值曲面，只是這個等值曲面并不光滑.用一個切面z＝z1去橫切三個這樣的等值曲面，也能得到三個圓.當(dāng)z1足夠接近甚至是等于z0時，能確保橫切出來的三個圓當(dāng)中，有兩個圓相交于兩點(diǎn).當(dāng)然，無論這兩個點(diǎn)當(dāng)中的哪一個點(diǎn)與第三個圓圓心之間的距離都可能不等于第三個圓的半徑，但只要兩者相差的值小于某個設(shè)定值，就可以認(rèn)為三個圓能交與一點(diǎn)，由此也就計算出了三個等值曲面的交點(diǎn).實際操作中，利用 MATLAB得到40條等高線對應(yīng)的離散點(diǎn)的坐標(biāo)，存儲在軟件中.當(dāng)需要用到某個特定電壓對應(yīng)的等高線時，就利用軟件插值去求取相應(yīng)的離散點(diǎn)的坐標(biāo).

圖10所描繪的是20mm×20mm磁鋼位移標(biāo)定實驗中，霍爾傳感器輸出電壓的等高線簇.由于是軟件插值，可以無限增加等高線簇的密度，即近似地認(rèn)為等高線簇構(gòu)成了一個面.將這個面繞著z軸旋轉(zhuǎn)一圈，就能得到一個圓柱體（φ40×h20，單位mm）.如圖6，三個霍爾傳感器的排列呈一個等邊三角形.三個傳感器對應(yīng)三個這樣的圓柱體，三個圓柱體在三維空間中相交得到一個立體圖形，如圖11.圖中的柱體描繪了測量的空間范圍，高為20mm，其上下底面如圖12（圖中灰區(qū)域）.也就是說，只要磁鋼下端面的中心點(diǎn)處在這個柱體中，就能測量出其所處位置的三維坐標(biāo).

圖11 測量范圍的三維視圖Figure 11 3Dview of measurement range

圖12 柱體的底面Figure 12 Underside of the cylinder

將測量裝置埋入巖土深部以后，測量出此時磁鋼下端面中心點(diǎn)處的三維坐標(biāo)（x1，y1，z1），也就是起始點(diǎn)坐標(biāo).當(dāng)土層發(fā)生形變，可測量出形變之后的坐標(biāo)值（x2，y2，z2）.通過簡單的幾何運(yùn)算可以得到垂直及水平的位移大小值，也可以采用公式（2）得出總位移的大小.

3 數(shù)據(jù)分析

建立圖13的系統(tǒng)三維笛卡爾坐標(biāo)系：

圖13 測量裝置系統(tǒng)坐標(biāo)系Figure 13 Coordinate system of the measurement equipment

三個霍爾傳感器陶瓷片的中心點(diǎn)A、B、C在V2上端面上的投影為A0，B0，C0.以B0點(diǎn)為三維笛卡爾坐標(biāo)系的原點(diǎn)，x軸沿著線段B0C0水平向右，y軸同處在V2上端面且垂直于x軸，z軸垂直xy平面向上.實際制作的測量裝置，B點(diǎn)的坐標(biāo)為（0，0，－13.5），C 點(diǎn)的坐標(biāo)為（12.325，0，－13.5），A 點(diǎn)的坐標(biāo)為（5.95，11，－13.5），單位mm.選取一個8mm×8mm×8mm區(qū)域進(jìn)行測量，一共216個點(diǎn).其中x軸坐標(biāo)分別取為2.1625、3.7625、5.3625、6.9625、8.5625、10.1625，y軸和z軸坐標(biāo)分別取為0、1.6、3.2、4.8、6.4、8.

測量結(jié)果顯示，有些點(diǎn)的x軸坐標(biāo)的測量誤差能達(dá)到0.5mm，y軸、z軸坐標(biāo)誤差較小.造成這種現(xiàn)象的一個原因是：圖10中的等高線是基本單調(diào)遞減的，而且等高線的起始段是趨向于水平的.起始段上的兩個點(diǎn)（x1，y1）、（x2，y2），差值y2－y1遠(yuǎn)小于差值x2－x1.實際測量時，由于時漂以及溫漂等因素的影響，電路測得的霍爾傳感器的輸出電壓會有微小的偏差.這個偏差的影響相當(dāng)于對等高線造成了y2－y1的差值，反應(yīng)在測量結(jié)果上就是x坐標(biāo)的測量誤差較大.進(jìn)行多次測量，對誤差進(jìn)行分析，得出非常相似的規(guī)律，由此表明，算法仍需改進(jìn).利用這些規(guī)律，結(jié)合實驗數(shù)據(jù)，采用優(yōu)化算法中的麥夸特法＋通用全局優(yōu)化法進(jìn)行擬合，得到三個坐標(biāo)各自的補(bǔ)償值的計算公式.擬合模型均采用如公式（3）的形式：

式（3）其中d為補(bǔ)償值，分別有dx、dy以及dz，對應(yīng)三個坐標(biāo)軸.p為常量，x、y、z均為初步測得的未經(jīng)補(bǔ)償?shù)淖鴺?biāo)值.將初步測得的值減去補(bǔ)償值就得到了最終的測量值.選取一次實驗，代入216個點(diǎn)的數(shù)據(jù)進(jìn)行擬合.三條擬合公式的殘差平方和均小于0.21，相關(guān)系數(shù)均大于0.95.x坐標(biāo)的補(bǔ)償值較大，最大能接近0.4mm，y、z坐標(biāo)的補(bǔ)償值較小.這里只將補(bǔ)償值dx計算公式的常量p羅列如下：35.38，－1.28，－1.86，－3.52，0.718，－0.436，2.64，－0.057，0.036，－0.182，72.142，6.07，－8.1，－5.155，－0.013，－0.328，3.491，－0.0288，0.0324，－0.265.將得到的三條補(bǔ)償公式去修正另幾次實驗的測量結(jié)果，效果較好.進(jìn)行一次新的實驗，并加入補(bǔ)償公式進(jìn)行校準(zhǔn)，一共測量216個點(diǎn)，x軸坐標(biāo)分別取為2.1625、3.7625、5.3625、6.9625、8.5625、10.1625，y軸和z軸坐標(biāo)分別取為0、1.6、3.2、4.8、6.4、8.圖14表示的是x坐標(biāo)測量的情況，216個點(diǎn)中x坐標(biāo)為2.1625的一共有36個點(diǎn)，其他五種取值同理.

圖14 x坐標(biāo)測量Figure 14 x－coordinate measurement

圖中橫坐標(biāo)表示第1到第36個點(diǎn)，縱坐標(biāo)表示x坐標(biāo)值.星點(diǎn)表示最終測量的坐標(biāo)值，黑線表示實際的坐標(biāo)值，兩者的差值就是實際測量的誤差.同樣地，y軸及z軸坐標(biāo)的測量誤差如圖15、圖16.

由于磁鋼產(chǎn)生的磁場分布非線性，空間中不同區(qū)域的點(diǎn)的測量精度是不同的，這點(diǎn)可以在圖形中反應(yīng)出來，即誤差有大有小.

4 結(jié) 語

土體的運(yùn)動將改變傳感器與磁鋼之間的相對位置，采用上述方法可計算出土體變形前后磁鋼的三維坐標(biāo)，通過簡單的幾何運(yùn)算就可以得到位移的大小值.本方案設(shè)計的三維位移測量方法能較全面地測量地下土層的形變情況，并達(dá)到一定的測量精度，從而為地質(zhì)災(zāi)害地下位移監(jiān)測提供了一種新的解決方案.

［1］仝達(dá)偉，張平之，吳重慶，等.滑坡監(jiān)測研究及其最新進(jìn)展［J］.傳感器世界，2005，11（6）：10－14.Tong Dawei，Zhang Pingzhi，Wu Chongqin，et al.Research and new development of landslide activity monitoring［J］.Sensor World，2005，11（6）：10－14.

［2］董勖壯.淺談滑坡成因及防治措施［J］.礦業(yè)工程，2009，7（1）：8－11.Dong Xunzhuang.Analysis of reasons for hill creep and its prevention［J］.Mining Engineering，2009，7（1）：8－11.

［3］Varnes D J.Slope movement types and processes［C］／／National Research Council.Washington，DC：［s.n］，1978：11－33.

［4］Schuster R L，Krized R J.Landslide Analysis and Control［M］.Washington，DC：National Academy of Sciences，1978：5－10.

［5］王翠珀，李繼濤.鉆孔測斜儀在某露天礦深層巖移監(jiān)測中的應(yīng)用［J］.山西建筑，2010，36（14）：91.Wang Cuipo，Li Jitao.The clinograph application in an open pit mine deep rock movement monitoring［J］.Shanxi Architecture，2010，36（14）：91.

［6］張麗芬，姚運(yùn)生，曾夏生，等.鉆孔測斜儀在高臺滑坡深部位移監(jiān)測中的應(yīng)用［J］.地質(zhì)災(zāi)害與環(huán)境保護(hù)，2007，18（4）：91－94.Zhang Lifen，Yao Yunsheng，Zeng Xiasheng，et al.Application of CX－03borehole inclinometer to enshi slope monitoring［J］.Journal of Genological Hazards and Environment Preservation，2007，18（4）：91－94.

［7］陳義軍，劉長武，徐 進(jìn)，等.新型套管式同軸多點(diǎn)位移計及其在隧道圍巖變形監(jiān)測中的應(yīng)用［J］.巖土工程學(xué)報，2008，30（7）：1085.Chen Yijun，Liu Changwu，Xu jin，et al.Application of new spigot coaxial rod extensometer in monitoring deformation of surrounding rock［J］.Chinese Journal Geotechnical Engineering，2008，30（7）：1085.

［8］張乾兵，朱維申，李 勇，等.洞群模型試驗中微型多點(diǎn)位移計的設(shè)計及應(yīng)用［J］.巖土力學(xué)，2011，32（2）：623－628.Sun Qianbin，Zhu Weishen，Li Yong，et al.Design of mini multipoint extensometer in geomechanical model test of cavern group and its application［J］.Rock and Soil Mechanics，2011，32（2）：623－628.

［9］劉文慶，趙 飛，董建輝，等.多點(diǎn)位移計在高速公路高邊坡穩(wěn)定性監(jiān)測中的應(yīng)用［J］.地質(zhì)災(zāi)害與環(huán)境保護(hù)，2010，21（4）：104－107.Liu Wenqin，Zhao Fei，Dong Jianhui，et al.Application of multiple－point rod extensometer to high slope of high slope of expressway stability monitoring［J］.Journal of Genological Hazards and Environment Preservation，2010，21（4）：104－107.

［10］譚捍華，傅鶴林.TDR技術(shù)在公路邊坡監(jiān)測中的應(yīng)用試驗［J］.巖土力學(xué)，2010，31（4）：1331－1336.Tan Hanhua，F(xiàn)u Helin.Testing study of application of time domain reflectometry to highway slope monitoring［J］.Rock and Soil Mechanics，2010，31（4）：1331－1336.

［11］張 青，史彥新.基于TDR的滑坡監(jiān)測系統(tǒng)［J］.儀器儀表學(xué)報，2005，26（11）：1199－1202.Zhang Qing，Shi Yanxin.The landslide monitoring system based on TDR［J］.Chinese Journal of Scientific Instrument，2005，26（11）：1199－1202.

［12］陳云敏，陳 赟，陳仁朋，等.滑坡監(jiān)測TDR技術(shù)的試驗研究［J］.巖石力學(xué)與工程學(xué)報，2004，23（16）：2748－2755.Chen Yunmin，Chen Yun，Chen Renpeng，et al.Testing study on applications of time domain reflectometry to slope monitoring［J］.Chinese Journal of Rock Mechanics and Engineering，2004，23（16）：2748－2755.

［13］Singer J.Development and testing of a time domain reflectometry（TDR）monitoring system for subsurface deformations［C］／／Proceedings of the European Rock Mechanics Symposium.Eurock：［s.n］，2010：613－616.

［14］李 濱，殷躍平，吳樹仁，等.多級旋轉(zhuǎn)黃土滑坡基本類型及特征分析［J］.工程地質(zhì)學(xué)報，2011，19（5）：704.Li Bin，Yin Yuepin，Wu Shuren，et al.Basic types and characteristics of multiple rotational landslides in loess［J］.Journal of Engineering Geology，2011，19（5）：704.

［15］楊 杰，李 青.地下位移實時監(jiān)測系統(tǒng)研究［J］.中國計量學(xué)院學(xué)報，2008，19（1）：41－46.Yang Jie，Li Qing.On real－time underground displacement monitor systems［J］.Journal of China Jiliang University，2008，19（1）：41－46.

［16］Shentu Nanying，Zhang Hongjian.Research on an electromagnetic induction－based deep displacement sensor［J］.IEEE Sensors Journal，2011，11（6）：1504－1515.

［17］池金谷，李 青.一種基于霍爾效應(yīng)的巖土地下位移精密測量方法［J］.工業(yè)控制計算機(jī)，2010，23（6）：63－64.Chi Jingu，Li Qing.Method of precise measurement of underground displacement based on hall effect［J］.Industrial Control Computer，2010，23（6）：63－64.