張 迪 李家存 吳中海 劉紹堂 盧 燕

1. 首都師范大學(xué)資源環(huán)境與旅游學(xué)院, 北京 100048;2. 河南工程學(xué)院土木工程學(xué)院, 河南 鄭州 451191;3. 中國地質(zhì)科學(xué)院地質(zhì)力學(xué)研究所, 北京 100081

0 引言

典型第四紀(jì)地貌體的高精度測(cè)量是活斷層定量研究的基礎(chǔ), 斷層活動(dòng)的最新證據(jù)往往是通過這些微地貌形態(tài)刻畫出來的。 通過分析活斷層的典型微地貌形態(tài), 可估算斷層位移量, 重構(gòu)斷層破裂歷史, 揭示斷層的活動(dòng)習(xí)性, 這對(duì)量化斷層的活動(dòng)強(qiáng)度、 評(píng)價(jià)區(qū)域內(nèi)地震的危險(xiǎn)性及復(fù)發(fā)性具有重要意義 (劉靜等, 2013; He et al. , 2016;Wei et al. , 2019; 唐毅, 2019; 肖坤澤和童亨茂,2020)。 傳統(tǒng)的測(cè)量方法受地質(zhì)環(huán)境、 數(shù)據(jù)采集密度和效率、 數(shù)據(jù)精度等因素影響, 已無法滿足活斷層研究中對(duì)精細(xì)化空間數(shù)據(jù)的要求。 隨著遙感技術(shù)的發(fā)展和遙感影像分辨率的不斷提高, 提供了從宏觀角度研究活斷層形態(tài)分布與空間展布的方法; 然而受遙感影像分辨率的限制, 難以刻畫出活斷層上構(gòu)造微小的地貌形態(tài)且無法獲取垂直方向上的地表變形數(shù)據(jù) ( 付碧宏等, 2008;Arrowsmith and Zielke, 2009; 張景發(fā)等, 2016)。

激光雷達(dá)掃描技術(shù) ( Light Detection And Ranging, LiDAR) 具有數(shù)據(jù)采集效率高、 高分辨率和非接觸測(cè)量等優(yōu)點(diǎn), 可深入復(fù)雜、 危險(xiǎn)環(huán)境中快速、 高效地獲取斷層微地貌形態(tài), 實(shí)現(xiàn)復(fù)雜地質(zhì)條件下斷層微地貌形態(tài)的高精度三維重現(xiàn), 并為活斷層研究提供多樣化的測(cè)繪數(shù)據(jù)產(chǎn)品, 如數(shù)字表面模型 (Digital Surface Model, DSM)、 數(shù)字高程模型 (Digital Elevation Model, DEM) 和等高線等, 這為活斷層微地貌形態(tài)識(shí)別與量化提供新的技術(shù)方法 (馬洪超, 2011; Oskin et al. , 2012;任治坤等, 2014; 馬曉雪等, 2016; Ren et al. ,2018)。 目前, 基于LiDAR (主要是機(jī)載LiDAR 和地面LiDAR) 的活斷層應(yīng)用研究主要集中在以下四個(gè)方面: ①地表破裂三維重建及斷錯(cuò)微地貌的精細(xì)化測(cè)量 (Kayen et al. , 2006; Gold et al. ,2012; 袁小祥, 2012; 劉靜等, 2013; 魏占玉等,2014; 譚錫斌等, 2015; Wilkinson et al. , 2015;劉宇平等, 2016; 康帥等, 2017); ②震間變形監(jiān)測(cè) (Wilkinson et al. , 2012; Gold et al. , 2013);③古地震探槽剖面數(shù)字化測(cè)量 (鄭文俊等, 2015;李占飛等, 2016); ④斷層面形貌變化的高精度測(cè)量 ( Renard et al. , 2006; 魏 占 玉 等, 2010;Candela and Renard, 2012)。 相對(duì)于地面LiDAR 技術(shù), 機(jī)載LiDAR 具有數(shù)據(jù)獲取成本昂貴、 數(shù)據(jù)采集及處理專業(yè)性較強(qiáng)等特點(diǎn), 比較適合開展大范圍內(nèi)、 大尺度上的活動(dòng)構(gòu)造形態(tài)識(shí)別與量化研究(馬洪超等, 2008; 馬洪超, 2011; Oskin et al. ,2012; 佘金星等, 2018), 國內(nèi)的大多數(shù)研究機(jī)構(gòu)一般不具備開展機(jī)載LiDAR 測(cè)量的條件; 而地面LiDAR 具有操作簡單、 數(shù)據(jù)采樣密度高和精度高等優(yōu)點(diǎn), 可實(shí)現(xiàn)一定范圍內(nèi)復(fù)雜地質(zhì)條件下活斷層典型地貌體的精細(xì)化測(cè)量, 重建斷層微地貌真三維模型, 進(jìn)而劃分和識(shí)別地貌單元, 并對(duì)斷層微地貌的重要特征參數(shù)進(jìn)行測(cè)量分析。

文章以川西理塘毛埡壩盆地北緣正斷層崖為研究對(duì)象, 詳細(xì)介紹了利用地面LiDAR 的斷層崖微地貌的數(shù)據(jù)采集及處理流程, 建立了斷層崖高精度的數(shù)字高程模型和真三維地表模型, 并分析了正斷層崖的地貌特征, 計(jì)算出正斷層錯(cuò)動(dòng)兩期最新地貌面的精確垂直位錯(cuò)量, 為利用地面LiDAR技術(shù)精細(xì)測(cè)量活斷層典型微地貌體形態(tài)提供了方法參考。

1 研究區(qū)概況

毛埡壩斷裂位于理塘左旋走滑活動(dòng)斷裂帶北西段, 順無量河上游, 沿大、 小毛埡壩盆地北緣發(fā)育, 整體上近西東向展布, 是控制著大、 小毛埡壩斷陷盆地發(fā)育的主邊界斷裂 (圖1; 周春景等, 2015; 白永健等, 2019)。 在遙感影像上沿毛埡壩盆地北緣山前地帶展布, 并依次形成斷層崖、斷層三角面和沖-洪積扇等斷裂展布的線性標(biāo)志,均指示該斷裂活動(dòng)性顯著。 另外, 地質(zhì)調(diào)查中在318 國道上二四九工班的西側(cè)發(fā)現(xiàn)巨型崩塌滑坡體, 其水平距離可達(dá)4 km, 前后緣高差約為750 m, 這顯然指示此斷裂帶上發(fā)生過大地震活動(dòng)。 目前學(xué)界對(duì)毛埡壩斷裂的活動(dòng)性質(zhì)仍有不同認(rèn)識(shí), 主要爭議在該斷裂為逆斷層性質(zhì)還是典型的正斷性質(zhì) (馬丹等; 2014; 周榮軍等, 2015;徐錫偉等, 2015; 吳中海等; 2015; Chevalier et al. , 2016)。 從衛(wèi)星遙感影像上可清晰分辨出a 處為斷層垂直錯(cuò)斷山前沖積扇形成的斷層崖 (圖1),流經(jīng)斷層的沖溝在斷層下盤呈羽狀展布 (圖2a),b 處為整齊排列斷面均傾向盆地的斷層三角面和斷層垂直切割山前沖積扇的形成斷層崖 (圖1, 圖2b)。 c 處斷層上升盤處為兩級(jí)冰水臺(tái)地, 斷層下盤處存在明顯的地塹和地壘構(gòu)造 (圖1, 圖2c)。在數(shù)據(jù)采集點(diǎn)洞玉溝溝口處 (圖1 中的紅色區(qū)域,圖2d) 存在階梯狀的斷層崖發(fā)育, 這些典型的構(gòu)造地貌均表明毛埡壩斷裂具有正斷性質(zhì)。

圖1 毛埡壩斷裂幾何展布Fig.1 Geometric distribution of the Maoyaba fault

實(shí)驗(yàn)點(diǎn)位于毛埡壩盆地北緣禾尼鄉(xiāng)洞玉溝溝口處 (30.23°N; 99.86°E), 海拔約4100 m。 此處存在兩級(jí)階梯狀的斷層崖, 整體地貌全景圖如圖2d。 斷層崖地貌形態(tài)清晰, 兩階梯斷層崖T1和T2高差明顯, 適合開展地面LiDAR 的數(shù)據(jù)采集實(shí)驗(yàn)。

2 數(shù)據(jù)采集與處理方法

2.1 數(shù)據(jù)采集

研究選擇FARO 公司生產(chǎn)的Focus3D X 330 型地面三維激光掃描儀進(jìn)行數(shù)據(jù)采集, 系統(tǒng)主要由激光測(cè)距系統(tǒng)、 掃描儀旋轉(zhuǎn)平臺(tái)、 內(nèi)置高分辨率相機(jī)、 軟件控制平臺(tái)、 附屬配套設(shè)備及后處理軟件等組成, 此類型的激光掃描儀高度集成、 體積小、重量輕且操作簡單, 其最大掃描距離為330 m, 點(diǎn)云密度較高且掃描精度可達(dá)到毫米級(jí) (測(cè)距10 m和25 m 時(shí)誤差為±2 mm), 滿足復(fù)雜地質(zhì)環(huán)境下高精度地形測(cè)量的要求 (圖3)。 數(shù)據(jù)采集流程主要包括工作路徑規(guī)劃、 掃描站點(diǎn)布置、 掃描站點(diǎn)之間的靶球放置、 掃描參數(shù)設(shè)置和控制點(diǎn)絕對(duì)坐標(biāo)的獲取 (圖4)。 根據(jù)斷層崖的分布范圍和地形條件, 現(xiàn)場(chǎng)勘察后規(guī)劃出最佳掃描路徑并布置掃描站點(diǎn)。 受采集區(qū)域地面植被反射率的影響, 地面三維激光掃描儀的實(shí)際測(cè)量距離值往往小于理論值, 為保證掃描范圍涵蓋整個(gè)采集區(qū)內(nèi)詳盡的地形特征, 此次數(shù)據(jù)采集共設(shè)置28 個(gè)站, 通過站與站之間放置的靶球 (3 個(gè)或3 個(gè)以上) 來配準(zhǔn)點(diǎn)云, 進(jìn)而實(shí)現(xiàn)所有站點(diǎn)云的整體配準(zhǔn)。 地面三維激光掃描儀在獲取精確空間形態(tài)時(shí), 內(nèi)置的高分辨相機(jī)也同步獲取彩色紋理信息。 為便于和其他遙感數(shù)據(jù)進(jìn)行協(xié)同分析, 需將整體配準(zhǔn)后的點(diǎn)云從相對(duì)坐標(biāo)系轉(zhuǎn)換到大地坐標(biāo)系下, 文中采用高精度GPS 靜態(tài)同步觀測(cè)3 個(gè)或3 個(gè)以上控制點(diǎn)的方法實(shí)現(xiàn)點(diǎn)云從相對(duì)坐標(biāo)系到絕對(duì)大地坐標(biāo)系的轉(zhuǎn)換。

圖2 毛埡壩斷裂地貌特征與實(shí)驗(yàn)點(diǎn)地貌Fig.2 Geological features of the Maoyaba fault and geomorphologic environment of the study area

圖3 地面三維激光掃描儀數(shù)據(jù)采集方法Fig.3 Data acquisition using TLS

圖4 地面三維激光數(shù)據(jù)采集與預(yù)處理流程圖Fig.4 Flow chart of the TLS data acquiring and processing

2.2 數(shù)據(jù)預(yù)處理

點(diǎn)云數(shù)據(jù)預(yù)處理過程主要包括點(diǎn)云配準(zhǔn)、 彩色點(diǎn)云生成、 坐標(biāo)轉(zhuǎn)換、 點(diǎn)云濾波和點(diǎn)云重采樣(圖4), 其中點(diǎn)云配準(zhǔn)、 彩色點(diǎn)云生成和坐標(biāo)轉(zhuǎn)換可通過激光掃描后處理軟件Faro scene 6.2 完成,而點(diǎn)云濾波、 重采樣及不規(guī)則三角構(gòu)網(wǎng)建模則是在軟件Geomagic Studio 中完成。 由于受研究區(qū)范圍、 地形遮擋影響以及掃描儀最大測(cè)量距離限制,單站式掃描無法實(shí)現(xiàn)斷層崖全方位地表信息的測(cè)量, 需采用多測(cè)站掃描通過后期點(diǎn)云拼接才能獲取完整的三維空間數(shù)據(jù)。 以首站點(diǎn)云為起始站,借助站與站之間靶球在軟件Faro scene 6.2 中可自動(dòng)完成兩站點(diǎn)云之間的配準(zhǔn), 進(jìn)而將28 站點(diǎn)云都拼接到同一相對(duì)坐標(biāo)系下實(shí)現(xiàn)點(diǎn)云整體配準(zhǔn) (圖5)。 由于掃描視角的限制, 圖上缺失點(diǎn)云的白色圓形區(qū)域?yàn)榈孛嫒S激光掃描儀測(cè)站所在位置。原始點(diǎn)云僅包含距離、 角度和反射強(qiáng)度信息, 將糾正后的數(shù)字圖像與點(diǎn)云進(jìn)行融合, 即將圖像上彩色信息精確賦予到對(duì)應(yīng)的點(diǎn)云位置以生成真彩色的三維點(diǎn)云, 從而豐富點(diǎn)云特征, 增加反映斷層崖的地表信息量。 由于受設(shè)備系統(tǒng)內(nèi)部因素和研究區(qū)內(nèi)植被、 較大石塊和車輛等環(huán)境影響 (圖2, 圖3), 地面三維激光掃描儀在數(shù)據(jù)采集過程中不可避免地會(huì)產(chǎn)生噪聲點(diǎn), 需進(jìn)行點(diǎn)云濾波。 針對(duì)斷層崖周圍比較明顯的散點(diǎn)及空中漂浮點(diǎn), 可采用手動(dòng)方法直接刪除; 而對(duì)于地面上的植被和較大石塊等噪聲點(diǎn), 選擇最大局部坡度濾波法通過不斷設(shè)定坡度閾值來實(shí)現(xiàn)噪聲點(diǎn)消除, 濾波前后效果如圖6。 多站拼接后的點(diǎn)云中包含重疊區(qū)域, 若直接進(jìn)行三角構(gòu)網(wǎng)建模, 既占用大量的計(jì)算機(jī)資源, 也會(huì)增加構(gòu)網(wǎng)擬合算法的難度, 同時(shí)也降低了模型構(gòu)建效率。 因此, 在保證曲面重構(gòu)精度的前提下, 需對(duì)點(diǎn)云進(jìn)行重采樣處理, 以剔除冗余數(shù)據(jù), 降低點(diǎn)云密度, 提高三維模型重建效率。

圖5 配準(zhǔn)后點(diǎn)云圖Fig.5 Point cloud image after registration

圖6 點(diǎn)云濾波效果圖Fig.6 Point cloud image after filtering

3 三維模型及定量分析

3.1 三維建模效果

數(shù)據(jù)預(yù)處理后的點(diǎn)云可通過軟件顯示 (圖7a、7b), 點(diǎn)云中空白部分為測(cè)站點(diǎn)和地面植被、 石塊濾波后的區(qū)域。 由于點(diǎn)云是不連續(xù)的離散點(diǎn), 從圖7a、 7b 上難以直觀識(shí)別出斷層崖的地貌特征及細(xì)微地形變化, 也無法進(jìn)行定量分析。 因此, 為更加形象、 直觀地分析斷層崖的地形地貌特征,在軟件Geomagic Studio 中對(duì)點(diǎn)云進(jìn)行不規(guī)則三角形構(gòu)網(wǎng) (圖7c、 7d), 其中圖7c 為分辨率為0.05 m 的DEM, 而圖7d 為紋理貼圖后的真彩色三維模型。 相對(duì)點(diǎn)云圖, DEM 和真彩色三維模型圖為連續(xù)曲面, 隨著地形的起伏變化, 三角面形狀、大小和密度也會(huì)隨著采樣點(diǎn)的位置和密度而發(fā)生變化, 由空間采樣點(diǎn)和斷線組成的不規(guī)則三角面可將復(fù)雜的地形特征更精確真實(shí)地表達(dá)出來, 并對(duì)地貌、 地形和地表變形精細(xì)特征進(jìn)行多視角分析。

圖7 斷層崖三維建模效果Fig.7 3D modeling of the fault scarp

3.2 斷層崖的地貌特征

地面LiDAR 可實(shí)現(xiàn)斷層崖微地貌形態(tài)的高精度三維重現(xiàn), 構(gòu)建了分辨率為0.05 m 的DEM 和真彩色三維模型, 這為活斷層微地貌地表變形及斷層空間分布的準(zhǔn)確識(shí)別與量化提供了高精度基礎(chǔ)數(shù)據(jù)。 活斷層經(jīng)過區(qū)域往往伴隨著地貌形態(tài)變化,如地表破裂、 斷層崖、 斷層三角面和水系形態(tài)變化等, 通過地貌形態(tài)特征參數(shù)的提取可確定出斷層的位置及走向。 斷層崖的坡度圖 (圖8a) 是將斷層崖上坡度變化用不同顏色繪制, 紅色區(qū)域表示坡度較大區(qū)域 (坡度值一般在22°以上), 即地貌形態(tài)變化較大; 而藍(lán)色區(qū)域則表示坡度較小區(qū)域 (坡度值在5°以內(nèi)), 即地貌形態(tài)變化較小, 通過坡度圖可清晰識(shí)別出斷層崖微地貌變形區(qū)域。斷層崖的等高線圖中等高距為0.5 m (圖8b), 斷層崖微地貌地形變化可通過等高線的密集程度清晰表示 (圖8b)。 將坡度圖與等高線圖的疊加 (圖8c), 綜合坡度和等高線變化可完成對(duì)斷層崖地貌形態(tài)的解譯, 此斷層崖近西東向發(fā)育, 存在階梯狀斷層崖 (T1和T2), T1和上盤交界區(qū)存在明顯地表破裂 (圖8d), 這清晰指示出區(qū)域內(nèi)斷層位置和走向。 由于斷層活動(dòng)的非均勻性, 兩斷崖在垂向地貌特征差別明顯, 斷崖T1地表變形較小, 與上盤交界處存在明顯地表破裂 (圖8c 中白色虛線區(qū)域), 而斷崖T2地表變形較大且分布范圍較大,斷崖T2西部形態(tài)保留較好, 東部可能受自然或人為破壞, 形態(tài)發(fā)生一定程度變化。 斷層近西東向依次錯(cuò)斷沖積扇形成斷崖T1和斷崖T2, 并在斷層上盤存在地勢(shì)較低的凹槽型區(qū)域, 初步判定為地塹區(qū) (圖8d)。

圖8 斷層崖地貌特征分析Fig.8 Slope maps and contour maps of the fault scarp

3.3 垂直位錯(cuò)量確定

由于受地殼構(gòu)造活動(dòng)的持續(xù)作用以及外界環(huán)境因素影響, 如長期風(fēng)化、 侵蝕和沉積作用, 使斷層上下盤的地形發(fā)生變化, 導(dǎo)致斷層上下盤的分界難以識(shí)別。 為最大程度上減少斷層崖垂直位錯(cuò)量的計(jì)算誤差, 在高精度DEM 上橫跨活斷層方向依次提取6 條地形剖面進(jìn)行分析 (圖8d、 圖9), 其中剖面1 和2 為T1斷崖處的垂直位錯(cuò)剖面,剖面3 和4 為T2斷崖左側(cè)的垂直位錯(cuò)剖面, 剖面5和6 為T2斷崖右側(cè)的垂直位錯(cuò)剖面。 在剖面1 和剖面2 上水平距離12 ～22 m 范圍內(nèi)高程值發(fā)生明顯變化, 且兩側(cè)坡度基本一致, 以12 m 和22 m 作為斷層崖上盤與下盤分界, 剖面1 和剖面2 擬合后的垂直斷距分別為1.1 m 和1.16 m, 因此T1斷崖處的平均垂直抬升量取兩剖面的平均值為1.13 m。同理, 剖面3 和4 在水平距離12 ～45 m 范圍內(nèi)高程值發(fā)生明顯變化, 以12 m 和45 m 作為此處斷層上下盤分界, 擬合后的垂直斷距分別為6.3 m 和6.4 m, 以平均值6.35 m 作為T2斷崖左側(cè)處的垂直斷距。 剖面5 和6 為T2斷崖右側(cè)處的二維地形剖面, 以水平距離10 m 和48 m 作為此處上下盤分界, 擬合后的垂直斷距分別為 6.1 m 和5.55 m,取平均值5.83 m 作為T2斷崖右側(cè)處的垂直斷距。綜合T2斷崖左側(cè)和右側(cè)的垂直斷距, 以兩處的平均值6.09 m 作為T2斷崖的垂直斷距, 與T1斷崖處的垂直斷距相比, T2斷崖左右兩側(cè)的垂直斷距擬合結(jié)果不一致, 主要是由于T2斷崖的分布范圍較大, 地表地貌形態(tài)分布不均勻。

4 討論

地面LiDAR 具有操作簡單、 數(shù)據(jù)采集效率高、高分辨率和非接觸測(cè)量等優(yōu)點(diǎn), 可快速、 高效地實(shí)現(xiàn)復(fù)雜地質(zhì)條件下活斷層典型地貌體的高精度三維測(cè)量, 更加直觀、 形象地識(shí)別厘米級(jí)甚至是亞厘米級(jí)的地貌構(gòu)造形態(tài), 為活斷層微地貌形態(tài)精細(xì)化測(cè)量及地貌特征參數(shù)量化提供有效的技術(shù)手段。 與無人機(jī)攝影測(cè)量技術(shù)相比, 地面LiDAR技術(shù)可去除植被等噪聲信息, 揭示普通光學(xué)影像無法獲取的活斷層地貌構(gòu)造信息, 特別適合復(fù)雜地質(zhì)條件下有限范圍內(nèi)的活斷層典型地貌體精細(xì)測(cè)量與分析。 但受采集區(qū)域地面植被等因素影響,地面LiDAR 的實(shí)際測(cè)量距離值往往小于理論值,為避免測(cè)量過程中遺漏信息, 數(shù)據(jù)采集過程中要充分考慮研究區(qū)范圍、 地形遮擋以及掃描儀最大測(cè)量距離限制等。 另外, 點(diǎn)云濾波中非地面點(diǎn)的去除直接關(guān)系到后期三維模型建立的精度, 應(yīng)根據(jù)研究區(qū)地物實(shí)際情況選擇合適的濾波算法。

綜合分析三維模型、 坡度圖和等高線圖, 可明顯識(shí)別出此探測(cè)點(diǎn)為階梯狀斷層崖, 近西東向發(fā)育, 在斷層上盤存在地塹區(qū), 兩斷層崖地貌特征差別明顯, 斷崖T1地表變形較小, 存在明顯地表破裂; 而斷崖T2地表變形較大且分布范圍較大,斷崖T2左側(cè)形態(tài)保留較好, 右側(cè)可能受自然或人為破壞, 形態(tài)發(fā)生一定程度變化。 根據(jù)沿兩斷層崖依次提取的6 條地形剖面, 計(jì)算出了兩斷層崖的垂直位錯(cuò)量分別為1.13 m 和6.09 m。 斷崖T1分布比較均勻且地表變形較小, 利用剖面1 和2 計(jì)算的垂直位錯(cuò)量基本一致。 斷層T2分布范圍較大且地貌形態(tài)變化較大, 斷層面坡度較大, 斷崖左側(cè)上下盤分界明顯, 可明顯識(shí)別出主斷層位置; 由于受自然或人類活動(dòng)等因素影響, 斷崖右側(cè)地貌形態(tài)變化明顯, 斷層面分為兩部分, 坡度較小, 上下盤分界不清晰, 導(dǎo)致剖面3, 4 與剖面5, 6 垂直位錯(cuò)量結(jié)果不一致。 對(duì)于地貌形態(tài)變化較大斷層崖垂直位錯(cuò)量獲取, 可選取多組二維地形剖面進(jìn)行擬合計(jì)算。

5 結(jié)論

利用地面LiDAR 技術(shù)重建了川西理塘毛埡壩斷裂禾尼處斷層崖的高精度三維模型, 詳細(xì)介紹了數(shù)據(jù)采集及處理方法, 驗(yàn)證了地面LiDAR 在活斷層典型地貌體的精細(xì)化測(cè)量上的有效性和實(shí)用性, 并獲得以下初步認(rèn)識(shí):

(1) 經(jīng)過點(diǎn)云配準(zhǔn)、 濾波、 重采樣和三角形構(gòu)網(wǎng)處理后, 建立了0.05 m 分辨率的數(shù)字高程模型和真彩色三維模型, 綜合坡度圖和等高線圖分析了斷層崖的地貌特征, 并獲取了川西理塘毛埡壩斷裂禾尼處正斷層錯(cuò)動(dòng)兩期最新地貌面的精確垂直位錯(cuò)量分別為1.13 m 和6.09 m。

(2) 地面LiDAR 技術(shù)可為研究活斷層的地貌和形變演化提供精細(xì)的基礎(chǔ)數(shù)據(jù), 將有助于進(jìn)一步定量或半定量分析斷層崖的幾何學(xué)、 運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)特征, 為活斷層微地貌高精度三維重現(xiàn)提供重要方法。