TM50在跨斷層水平位移(短程)測量中的應用

易天陽, 王偉力, 溫軍軍, 卓如峰

(四川省地震局測繪工程院,四川 雅安 625000)

技術交流

TM50在跨斷層水平位移(短程)測量中的應用

易天陽, 王偉力, 溫軍軍, 卓如峰

(四川省地震局測繪工程院,四川 雅安 625000)

對于短程跨斷層水平位移測量，采用高精度光電測距儀是較好的選擇，就目前的測量技術與作業方法而言，徠卡測量機器人TM50無論在測量精度還是作業速度上都具有很好的表現。通過兩個跨斷層場地的實例展現了TM50在跨斷層水平位移測量中的應用，并對多周期實例數據進行了比較分析。

斷層； 位移測量； 氣象修正； 周期誤差修正； 儀器常數修正

0 引言

跨斷層微量位移測量的目的在于了解斷層活動的方式，以便探索斷層活動與地震的關系[1]。對于跨斷層位移測量，通常使用測量學的方法長期連續地精密測量斷層兩側的垂直位移和水平位移，以確定斷層的活動情況。由于觀測技術條件的局限，多年以來水平位移觀測常采用因瓦基線尺進行數據獲取。該方法測量精度高、受環境溫度和氣壓的影響小，但往、返測交換基線尺的方向有時會受到場地的限制而導致操作難度大、測量距離短、作業人員需求多、勞動強度大、場地跨越公路的情況存在安全隱患以及跨溝壑的斷層場地難以進行數據采集等缺點。隨著觀測技術的進步，GPS、光電測距成為目前跨斷層水平位移測量的主要觀測技術[2]，其中高精度GPS技術由于其自身特點主要用來監測火山地震、構造地震、全球板塊運動，尤其是板塊邊界地區的重要手段。文獻[3]提到了基于GPS技術等大地形變監測獲取的現今地殼運動的微動態定量信息對地震預測是十分重要的。而對于短程距離(3 km以內)的跨斷層場地，采用高精度光電測距儀是目前跨斷層位移測量的較好選擇(表1)。對于既有跨斷層的短程測距觀測場地，當點間通視存在困難時，為了能將多年累積下來的觀測資料予以利用和比較分析，可以采用自由設站的方式進行數據采集，但觀測限差和成果精度需滿足規范的指標要求，對于新建場地應考慮點間采用對向觀測的測量方式進行點位的選埋。

表 1 測距方式比較表

1 技術要求

依據《中、短程光電測距規范》及《地震地殼形變觀測方法 跨斷層位移測量》規范，儀器設備主要技術要求列于表2。

表 2 儀器設備技術指標[2，4]

2 數據采集

測量過程采用機載程序軟件進行測回數、正倒鏡、自動照準、觀測限差等的自動控制。為了保證測回間的測距較差限差，要求一次設站(可能會有多個觀測目標)的所有測回應在短時間內完成，從而保證在一次設站的觀測時間段內氣象要素不會有太大的變化，因此一次設站的測回數(規范要求不少于4個測回)和方向數都不宜太多，同時對儀器的觀測精度和觀測速度提出了較為苛刻的要求。在目前的技術條件下，徠卡測量機器人TM50能夠自動精確照準目標并快速完成測量，而且其觀測精度能夠滿足規范的要求。儀器測量時采用標準氣象條件下的氣象值進行數據采集，以輔助觀測的方式讀取并記錄相應站點觀測過程的儀器高、鏡高和氣象元素，內業數據處理時匹配其實際觀測時的氣象信息來進行觀測斜距的氣象改正，并輔以觀測標志間的高差、儀器高、棱鏡高進行斜距改正(或采用天頂距改平)。數據采集技術指標如表3、4。

表 3 距離觀測技術要求[2,4]

表 4 距離測量較差限值[2，4](單位:mm)

當采用三角高差(天頂距)進行斜距改平時，垂直角的觀測限差應滿足表5的要求。

表 5 垂直角觀測技術要求[5] (單位：s)

3 數據處理

儀器獲取數據為標準氣象條件下的斜距、天頂距、水平角，為了使周期間獲取的數據具有可比性，需要對標準氣象條件下獲得的數據進行各項修正計算。

3.1 氣象修正

這是電磁波測距的重要改正，因為電磁波在大氣中傳輸時受氣象條件的影響很大。此項改正的實質是大氣折射率對距離的改正,因折射率與氣壓、溫度、濕度有關，所以習慣上稱為氣象改正[6]。在作業時實際氣象條件下的群折射率n的計算公式為[6-7]：

(1)

式中:λ為真空中光波的有效波長，單位為μm。

氣象修正按式(3)計算[4]：

(3)

式中:D′(m)為觀測距離，下同。

3.2 周期誤差修正

3.3 儀器常數修正

儀器加、乘常數屬于系統誤差，可通過儀器檢定機構對儀器的年檢獲取，測量使用配套的徠卡圓棱鏡。加常數修正值按式(4)計算，乘常數修正值按式(5)計算[4]：

(4)

(5)

3.4 斜距改平

邊長觀測時在觀測墩上的每次設站儀器高和目標架設高度不能夠相同，使得期次間的觀測斜距值不能直接進行比較。當觀測距離D′經過各項修正后的測量斜距D化算為水平距離且將水平距離化算到同一參考面上才能將對應邊的數據進行周期比較，從而判定斷層隨時間的相對位移變化情況。

(1) 三角高差改平

當對向觀測三角高差成果滿足式(10)要求時利用觀測天頂距和斜距計算水平距離。對于短程距離，將測量距離(斜距)進行傾斜改正，化為水平距離按式(6)計算。

(6)

(2) 水準高差改平

觀測距離D′經過各項修正以后，對于短程距離將測量距離化算為水平距離按式(7)計算。

(7)

需要注意的是，ΔH并不是直接采用水準高差的測量值，因為儀器測量的是儀器中心與棱鏡中心間的斜距，計算采用的ΔH除包含水準高差的測量值外，還應包含儀器高與目標高度的差值。

圖1 水準高差、儀器高、鏡高示意圖Fig.1 Sketch of leveling height different,instrument height, and target height

依據圖1所示，i為儀器高，v為目標高，hAB是水準點A、B間的高差，ΔH是斜距改平采用高差，D為修正后的斜距。則IM間的高差ΔH有如下關系式：ΔH+i=hAB+v，從而有ΔH=hAB+v-i。另外，如果知道A、B點的高程HA、HB的話，則IM間的高差ΔH有如下關系式：ΔH=HB+v-(HA+i)。

4 精度評定

精度評定是對測量成果可靠性進行的評估。在跨斷層位移測量中似乎應該忽略測量角度的影響。對于場地環境能夠方便等級水準測量高差而進行傾斜改正的話，對測量角度加以關注是不必要的；但是對于端點間不方便等級水準觀測，但滿足采用三角高差進行傾斜改正的場地，除了對對向觀測距離進行精度評估外，還需對三角高差成果的符合性(對向高差之差)進行評估。

4.1 測量距離精度評定

距離測量結果的精度評定參照《中、短程光電測距規范》對向觀測距離的精度評定要求進行中誤差計算,并參照《跨斷層測量規范》關于測距相對中誤差的限值進行合格評定。

(1) 對向觀測的平均值中誤差[4]：

(8)

(2) 相對中誤差[4]：

(9)

4.2 三角高差成果符合性

采用三角測量確定高程差并進行傾斜修正。對向觀測時，高差之差應滿足式(10)的要求。

(10)

5 應用實例

(1) 點間高差滿足三角高差改平的符合性判定(表6)

(2) 計算實例

實例涉及4個周期的數據。由于數據表較多，不予詳細羅列，兩個場地分別以最近一期成果數據為例(表7，表8)。

依據測前點間高差符合性判定，寶興縣五龍鄉場地的A-D、B-D邊和西昌市大箐鄉A-C邊滿足三角高差進行改平，首先對三角高差成果的符合性進行判定(表9)。在滿足規范要求的前提下進行改平計算(表10)，并將高差改平的結果與水準高差改平的結果進行比較(表11)。

表 8 2015年9月西昌市大箐鄉場地距離測量成果表

表 9 三角高差成果用于傾斜改平的符合性判定表

表 10 三角高差改平的成果表

(3) 周期數據的測量成果比較分析

獨立的一周期測量只能反映當時周期的測量成果自身的精度情況。為了能夠跟蹤監測跨斷層場地端點間的水平位移相對變化情況，需要對對應邊的各期次進行期次間的數據比較，并對各期次的成果數據較初始期次數據成果的累積變化予以比較和分析(圖2)。

表 11 不同改平方式的成果較差表

從圖3、表12可以看出，寶興縣五龍鄉場地各監測邊在周期間的最大較差在B-C(斷層同側點)邊(2015.9)，達-2.49 mm。該場地測量邊相鄰期的距離相對變化值最大為6.44×10-6，其他較差均在±1.3 mm以內，期間變化較為平穩；較初始期的累積變化量均在±2 mm以內，最大為-1.76 mm，距離相對變化值最大為4.55×10-6，數據對場地的表現為平穩。

從圖4、表13可以看出，西昌市大箐鄉場地斷層同側的點間距離相對變化無論是期間較差還是累積變化，都表現得比較穩定，沒有急劇變化；斷層兩側監測點間距離的相對變化較為活躍，起初體現為斷層拉開(2015年3月—2015年5月)，接下來的期間變化體現斷層繼續拉開的幅度有所放緩(2015年5月—2015年7月)，在2015年9月時部分點間距的變化表現為急劇縮短，變化量在(-5.2，-2.9)mm，該場地測量邊相鄰期的距離相對變化值最大為5.77×10-6；從累積變化圖也不難看出，較初始期的累積變化量逐漸增大(2015年3月—2015年7月)，后又急劇減小(2015年7月—2015年9月)，向初始測量值靠近并呈壓縮態勢(平→張→壓)，距離相對變化值最大為4.42×10-6，數據對場地的表現為平穩。

圖2 兩綜合場地Fig.2 Two cross-fault sites

圖3 寶興縣五龍鄉場地Fig.3 Site in Wulong town，Baoxing county

邊名稱A?C月份邊長值/m期間較差累積較差A?D邊長值/m期間較差累積較差B?C邊長值/m期間較差累積較差2015．5769．24888相對變化相對變化545．67762相對變化相對變化386．71460相對變化相對變化2015．3769．247950．930．931．21E-061．21E-06545．677560．060．061．10E-071．10E-07386．71516-0．56-0．56-1．45E-06-1．45E-062015．7769．24814-0．740．19-9．62E-072．47E-07545．67611-1．51-1．45-2．77E-06-2．66E-06386．715891．290．733．34E-061．89E-062015．9769．24687-1．27-1．08-1．65E-06-1．40E-06545．67577-0．34-1．79-6．23E-07-3．28E-06386．71340-2．49-1．76-6．44E-06-4．55E-06邊名稱B?D月份邊長值/m期間較差累積較差C?D邊長值/m期間較差累積較差2015．5645．43370相對變化相對變化802．36973相對變化相對變化2015．3645．433810．110．111．70E-071．70E-07802．370160．430．435．36E-075．36E-072015．7645．43338-0．43-0．32-6．66E-07-4．96E-07802．36904-1．12-0．69-1．40E-06-8．60E-072015．9645．433490．11-0．211．70E-07-3．25E-07802．36810-0．94-1．63-1．17E-06-2．03E-06

6 結語

通過TM50在跨斷層位移測量中4期觀測的應用實例，在觀測指標滿足相關規范要求下進行數據采集，并對觀測數據進行相關修正和不同方式的傾斜改平后發現：

(1) 高差改平的測量邊精度滿足規范要求；

(2) 滿足三角高差改平的測量邊精度滿足規范要求；

表 13 西昌市大箐鄉場地測距成果分析表(較差值單位數為mm)

圖4 西昌市大箐鄉場地Fig.4 Site in Daqing town，Xichang city

(3) 滿足三角高差改平的測量邊與水準高差改平的測量邊的結果基本一致，其較差滿足規范要求的較差限值。

另外，采用TM50進行跨斷層水平位移的監測方法便于內業數據資料整理，其作業速度快，人員需求少，較常規因瓦基線尺測距能夠大大降低野外作業人員的勞動強度。

注：應用及數據來源于四川省地震局測繪工程院與中國地震局地震預測研究所協作關于“龍門山斷裂帶南段形變綜合觀測場及安寧河斷裂帶形變綜合觀測場”的動態跟蹤研究項目。

References)

[1] 龔鴻慶.關于跨斷層位移測量的資料處理問題[J].地殼形變與地震,1983(1):48. GONG Hong-qing.Data Processing Related to Displacement Measurements across Faults[J].Crustal Deformation and Earthquake,1983(1):48.(in Chinese)

[2] 中國地震局.DB/T 47-2012,地震地殼形變觀測方法 跨斷層位移測量[S]北京：地震出版社，2012. China Earthquake Administration.DB/T 47-2012,The Method of Earthquake-related Crust Monitoring——Fault-crossing Displacement Measurement[S]:Beijing:Seismological Press,2012.(in Chinese).

[3] 王雙緒,丁平,張希,等.利用GPS資料研究印尼8.7級地震對我國西部地震趨勢的影響[J].西北地震學報,2005,27(4):322. WANG Shuang-xu,DING Ping,ZHANG Xi, et al.Research on Influence of IndonesiaMS8.7 Earthquake to Western China Using GPS Data[J].Northwestern Seismological Journal,2005,27(4):322.(in Chinese).

[4] 國家測繪局測繪標準化研究所.GB/T 16818-2008,中、短程光電測距規范[S].北京：中國標準出版社，2008：4-10. Institute of Suveying and Mapping Standardization of State Bureau of Suveying and Maping.GB/T 16818-2008,Specifications for Short to Medium RangesElectro-optical Distance Measurements[S].Beijing:Standards Press of China,2008:4-10.(in Chinese).

[5] 國家測繪局測繪標準化研究所.GB/T 17942-2000,國家三角測量規范[S].北京：中國標準出版版，2000:13. Institute of Suveying and Mapping Standardization of State Bureau of Suveying and Maping.GB/T 17942-2000, Specifications for National Triangulation[S].Beijing:Standards Press of China,2000：13.(in Chinese).

[6] 孔祥元,郭際明,劉宗泉.大地測量學基礎[M].武漢:武漢大學出版社,2008:304-314. KONG Xiang-yuan,GUO Ji-ming,LIU Zong-quan.Foundation of Geodesy[M].Wuhan:Wuhan University Press,2008:304-314.(in Chinese).

[7] 姚輝,陳夙穎.全站儀氣象改正公式及氣象元素測量精度對距離的影響[J].測繪通報,2008(4):14. YAO Hui,CHEN Su-ying.Total Station’s Meteorological Correction Formula and the Influence of Measuring Precision of Meteorological Element on Distance[J].Bulletin of Surveying and Mapping,2008(4):14.(in Chinese).

[8] 國家地震局.跨斷層測量規范[S]:北京：地震出版社，1991:4,44. National Seismological Bureau.Specifications of Cross-fault Geodetic Measurements[S].Beijing:Seismlogical Press,1991:4,44.(in Chinese)

Application of Leica TM50 in Horizontal Displacement (Short Distance) Measurements across Faults

YI Tian-yang, WAMG Wei-li, WEN Jun-jun, ZHUO Ru-feng

(SurveyingandEngineeringInstitute,EarthquakeAdministrationofSichuanProvince,Ya’an625000,Sichuan,China)

Based on the application of the Leica TM50 robotic total station to short-distance measurement of cross-fault displacements, the aim of this study is to show the feasibility of an intelligent measuring robot with nominal precision in short-distance measurement. This paper compares the advantages and disadvantages of the traditional invar-wire baseline method, GPS, and measuring robots in relation to aspects such as measurement precision, operational mode, and personnel demand. It finds that the use of a high-precision electro-optical distance meter (EDM) is the best choice at present for relatively short (<3 km) cross-fault displacement measurements. The study then clarifies the operational basis and technical requirements of an EDM and the operational hazards in data collection. In order to ensure repeatable and comparable measurement data, the observational field data need to be corrected with data processing. Meteorological modifications, periodic error correction, and instrumental correction are all discussed. Results are presented from an application example, compared graphically, and briefly analyzed in relation to fault displacements that are contained therein.

fault; displacement measurement; meteorologic modification; periodic error correction; instrument constant correction

2015-10-27 基金項目：中國地震局地震預測研究所基本科研業務專項(2014IES010201,2014IES010201)

易天陽(1976-)，男(漢族)，四川遂寧人，工程師，碩士，從事地震監測工作。E-mail:sunsky-yi@163.com。

TD178

A

1000-0844(2016)06-1010-07

10.3969/j.issn.1000-0844.2016.06.1010