基于TLS技術的典型建筑物震害信息三維建模分析——以彭州市白鹿中學為例

焦其松， 張景發， 蔣洪波， 宿淵源， 王旭

(1.中國地震局工程力學研究所，哈爾濱　150080； 2.中國地震局地殼應力研究所地殼

動力學重點實驗室，北京　100085)

?

基于TLS技術的典型建筑物震害信息三維建模分析
——以彭州市白鹿中學為例

焦其松1，2， 張景發2， 蔣洪波2， 宿淵源2， 王旭2

(1.中國地震局工程力學研究所，哈爾濱150080； 2.中國地震局地殼應力研究所地殼

動力學重點實驗室，北京100085)

摘要：傳統光學及雷達遙感技術在地震災害損失評估中無法獲取建筑物立面破壞信息，對結構已遭破壞，但頂面基本完好的建筑物識別能力較差。三維激光掃描技術(terrestrial laser scanning,TLS)可以獲取建筑物結構破壞狀況的高精度三維信息。利用2013年10月23日RIEGL VZ-1000三維激光掃描儀采集到的汶川地震震后彭州市白鹿中學的精細點云數據，通過不規則三角網建模和紋理貼圖等處理進行三維建模分析，在此基礎上對建筑物的典型破壞特征及地表破裂情況進行了詳細建模研究，真實再現了白鹿中學破壞場景，實現了震害信息的永久保存。研究結果表明： 位于斷層上盤的勤學樓整體無明顯破損； 位于下盤的求知樓損毀相對嚴重，普遍出現“X”形的共軛剪切裂縫、構造柱斷裂及露筋等現象； 地震陡坎沉降量不大，高程剖面顯示現今該陡坎最低高度為205 cm，最高為231 cm，與GPS測量結果大致吻合。

關鍵詞：白鹿中學； 三維激光掃描(TLS)； 震害信息提?。?三維建模

0引言

隨著遙感數據多樣化和數據處理技術的發展，利用遙感技術探索與地震相關的信息已經逐漸成為地震研究中的一個熱點，覆蓋了地震研究的方方面面[1]。遙感技術可以大范圍快速獲取地震造成的建筑物破壞信息并進行損失評估，對于抗震救災快速響應、減少震害損失和災后重建等具有重要意義，尤其是在2008年汶川地震及后續多次嚴重地震災害中，多種遙感技術的應用發揮了重要作用[2-5]。地面三維激光掃描技術(terrestrial laser scanning,TLS)是繼GPS之后的又一項測繪新技術，它利用激光測距快速、非接觸、高精度地獲取目標對象表面高密度的空間坐標，實現目標三維信息提取與建模。其應用領域早期主要集中于文物保護及工業逆向工程檢測。近年來，TLS不斷發展并日漸成熟，國內對該技術的研究與應用已逐漸興起并取得了一定成果，其應用已經拓展到土木工程、地質調查、數字城市以及災害監測等領域[6-9]。傳統光學及雷達遙感技術獲取的建筑物破壞信息是宏觀的，雖然能較好地識別頂面破壞及完全倒塌的建筑物信息，但無法反映建筑物立面的破壞情況； 對于結構已經破壞，但頂面基本完好的建筑的識別能力更低，影響了損失評估精度。TLS則可以獲取建筑物結構破壞的高精度三維信息，實現任意角度建筑物的側面損毀狀況及地表破裂特征的定量量測。

2008年5月12日汶川地震使正好位于龍門山前山斷裂帶上的彭州市白鹿中學慘遭破壞。該斷裂從學校中間穿過，斷裂兩側分別是中學的初中部求知樓和小學部勤學樓，兩樓體均為砌體結構，承重墻體的構造特點有別于當地的普通砌體結構。其中求知樓破壞嚴重，勤學樓基本完好，其破壞程度與普通等剛度砌體結構破壞形式有很大不同，沒有出現底層破壞嚴重、上層結構破壞較輕的一般形態，極具研究價值[10]。2013年10月23日筆者利用RIEGL VZ-1000三維激光掃描儀采集了白鹿中學的精細點云數據，并利用不規則三角網建模和紋理貼圖等處理進行三維建模分析，在此基礎上選取建筑物典型破壞特征及地表破裂情況進行了詳細建模研究。

1研究區概況及數據采集

1.1研究區概況

龍門山斷裂帶是汶川8.0級地震的發震構造帶，其中前山斷裂帶穿過彭州市白鹿中學和湔江支流，并斷錯了河床、河漫灘和階地。地震發生后，在白鹿中學兩棟教學樓之間形成走向NE48°，右旋位移1 m，垂直位移1.8 m的地震陡坎。南東側的下盤距主斷面18 m處的求知樓地面破裂，教學樓墻體多處開裂，窗戶玻璃全部破損； 而上盤距離主斷面約7.6 m的勤學樓框架完好，且窗戶玻璃有部分保存完好，斷層下盤教學樓損壞大于上盤的教學樓[11-13]。

1.2數據采集

本次數據采集采用的儀器為VZ-1000，由奧地利RIEGL公司開發，整個系統由地面三維掃描儀(VZ-1000)、數碼相機(尼康D600)、后處理軟件(RISCAN Pro)、電源以及附屬設備構成，其最大掃描距離為1 200 m，單次掃描精度為5 mm(100 m距離處，一次單點掃描)，激光發射頻率為300 000點/s，掃描視場范圍為100°×360°(垂直×水平)，為快速、高效及高精度地獲取震害信息提供了保障。掃描重點區域為白鹿中學的兩棟教學樓及樓間的地震陡坎。圖1為白鹿中學地貌特征及掃描儀站點位置平面圖。

圖1　白鹿中學地貌特征及掃描站點位置平面圖

在進行站點勘查時先確定了站點的視角范圍，相鄰站點間保證有重疊區域，以利于后續數據拼接處理。根據現場實際情況合理制定了數據采集方案，共設置10個站點，掃描范圍覆蓋了研究區內各類特征地形及震害建筑物特征。

2數據處理及建模

2.1數據處理

對三維激光掃描儀采集到的點云數據的處理主要包括噪聲剔除、點云匹配、區域分割裁剪、坐標變換及點云抽稀等，其目的是為了將不同站點掃描獲得的點云進行正確匹配，統一到特定的相對坐標系中，以便對目標進行建模和定量化分析。

數據在采集過程中易受環境和系統本身的影響，如空氣中的水汽、儀器旋轉引起的抖動以及運動物體的干擾等，因這些影響所導致的噪聲點主要包括空中懸浮的點和三維物體附近的孤立點等。儀器在出廠時曾經過嚴格檢校，在非極端環境下使用時產生的誤差可以忽略不計。本次數據采集于秋季，當地氣候潮濕，空氣中氣溶膠含量較高，在點云數據中表現為低反射率、成團聚集的點簇，對建模影響較大，因此，需進行噪聲剔除工作。具體做法是將點云投影到正交平面上，以水平面為基準，基于高程信息，手動刪除建筑物周圍以及空中飄浮的散點。由于掃描目標范圍較大且地形復雜，激光掃描儀需要從不同位置、多視角進行若干幅掃描，而掃描所得到的點云數據都是以掃描儀位置為零點的局部坐標，因此必須對點云數據進行匹配。白鹿中學測區范圍較小，點云重疊區域較大，相鄰站點之間同名點較多，因此采用基于多測站獲取重疊區域點云數據的方法進行自動匹配。匹配前先將點云生成多面數據(polydata)，以降低匹配的復雜度和誤差，從而提高拼接效率； 然后在保證不影響曲面重構和保持精度的情況下對點云數據進行平滑，剔除冗余數據，降低原始點云數據密度以節省計算時間。匹配時首先導入一站點云數據作為起始匹配站點，再導入與之相鄰的點云數據，設置搜索距離，采用迭代最近點法(iterative closest point)搜索同名點，同名點之間采用最小二乘擬合法進行計算直至拼接的誤差最小，完成數據拼接。為保證多站拼接的精度，將標準差小于0.005 mm的點視為完全匹配，得到的匹配結果如圖2(a)所示。最后將經過幾何糾正的實景照片與點云數據精確配準，將照片上每個像素點的RGB彩色信息提取出來并賦到相應位置的點云上，實現點云數據的真彩色顯示。賦色之后的點云能夠較好地顯示地表信息，如圖2(b)所示。

(a) 點云匹配效果(b) 點云賦色效果

圖2白鹿中學局部點云

Fig.2Partial point clouds of Bailu middle school

2.2不規則三角網建模

利用激光三維掃描儀獲取的原始點云是不連續的離散點，對其進行不規則三角網建模，變成連續的曲面可以更加直觀地顯示震害信息。 RISCAN Pro三角網建模方法分為基于平面投影和基于極點投影2種，這2種方法在地形建模上有一定的優勢而在建筑物建模上優勢不大。因此，本文將建筑物及地震陡坎分割以后分別建模，再通過紋理貼圖得到點云重建的真三維模型，如圖3所示。

圖3　白鹿中學三維建模效果

3結果分析

在汶川地震中，白鹿中學教學樓破壞嚴重，其中砌體結構的破壞尤為嚴重。圖4為求知樓背面墻體的一組“X”形的共軛剪切裂縫三維建模圖，原始點云按0.01 m間距抽稀以后進行不規則三角網建模，然后在保證精度和模型曲率的條件下合并冗余三角形，對三角網模型進行簡化(圖4b)。墻體裂縫通過計算模型曲率顯示，最后利用照片對其進行紋理貼圖，得到該組裂縫的真三維模型，真實再現了墻體破壞情況。白鹿中學兩棟教學樓距斷層最近，位于斷層上盤的勤學樓隨上盤的逆沖整體上移而相對位置抬高，房屋框架完好，墻面有剝落的痕跡，整體無明顯破損； 位于下盤的求知樓在上盤擠壓作用下，墻體受到剪切破壞，并沿粘結強度最薄弱、應力最集中的地方(窗臺下方)裂開，普遍出現“X”形的共軛剪切裂縫、構造柱斷裂及露筋等現象。

地震發生時，求知樓與勤學樓之間的水泥地面經歷了瞬時性的逆沖和構造縮短形變過程，形成長度約65 m的地震陡坎。任俊杰等[11]震后地表調查發現白鹿中學坐落于T2階地面老陡坎坡面上(圖1)，GPS測量結果表明穿過白鹿中學水泥地面陡坎測線的垂直位錯量為(196±20) cm[13]。

(a) 點云渲染(b) 不規則三角網建模

(c) 模型曲率 (d) 真實紋理貼

圖4求知樓“X”形共軛剪切裂縫三維建模

Fig.43D model of “X”-shaped conjugate shear fracture in Qiuzhi building

圖5為三維激光掃描儀采集的點云數據不規則三角網建模結果。

圖5　地震陡坎不規則三角網建模及剖面圖

由于白鹿中學現今已劃為地震遺址進行保護，人為破壞較小，而且在震后遺留的水泥地板支撐保護下，陡坎沉降量不大，高程剖面顯示現今該陡坎最低高度為205 cm，最高為231 cm，與GPS測量結果大致吻合。

4結論

1)本文利用RIEGL VZ-1000 TLS采集了震后彭州市白鹿中學教學樓及地震陡坎的高精度點云數據，并利用不規則三角網建模和紋理貼圖等處理進行三維建模分析; 在此基礎上選取建筑物典型破壞特征及地表破裂的產狀信息進行了詳細建模研究，實現了任意角度建筑物的表面損毀狀況及地表破裂特征的定量量測，真實再現了現今白鹿中學遭受破壞的場景，實現了震害信息的永久保存，為科研人員進行進一步研究分析提供了可能。本次數據采集及建模精度可以達到mm級。

2)與傳統GPS基于單點的測量方法相比，TLS數據采集速度快、分辨率和精度高且屬于非接觸式測量，可以在很多復雜環境下應用，未來必將成為傳統的地震地質調查、房屋結構調查及分析的重要手段。

3)本文只對損壞程度較小且外形基本完整的房屋進行了模型重建。因此在后續工作中還要進一步研究多種復雜形狀與受損程度的建筑物，以及地表破裂的三維模擬重建技術，以實現對地震災害更為有效的評估。

參考文獻(References)：

[1]荊鳳,申旭輝,洪順英,等.遙感技術在地震科學研究中的應用[J].國土資源遙感,2008,20(2):5-8.doi:10.6046/gtzyyg.2008.02.02.

Jing F,Shen X H,Hong S Y,et al.The application of remote sensing technology to earthquake science research[J].Remote Sensing for Land and Resources,2008,20(2):5-8.doi:10.6046/gtzyyg.2008.02.02.

[2]王治華,周英杰,徐斌,等.“5·12”汶川大地震震中區映秀鎮地震災情及次生地質災害遙感初步調查[J].國土資源遙感,2008,20(2):1-4.doi:10.6046/gtzyyg.2008.02.01.

Wang Z H,Zhou Y J,Xu B,et al.Preliminary remote sensing investigation of damage caused by the“5·12”Wenchuan strong earthquake as well as secondary hazards at the quake center-Yingxiu town[J].Remote Sensing for Land and Resources,2008,20(2):1-4.doi:10.6046/gtzyyg.2008.02.01.

[3]郭兆成,童立強,鄭雄偉,等.四川蘆山地震次生地質災害遙感調查及災害特征初探[J].國土資源遙感,2014,26(3):99-105.doi:10.6046/gtzyyg.2014.03.16.

Guo Z C,Tong L Q,Zheng X W,et al.Remote sensing survey of secondary geological disasters triggered by Lushan earthquake in Sichuan Province and tentative discussion on disaster characteristics[J].Remote Sensing for Land and Resources,2014,26(3):99-105.doi:10.6046/gtzyyg.2014.03.16.

[4]于海洋,程鋼,張育民,等.基于LiDAR和航空影像的地震災害倒塌建筑物信息提取[J].國土資源遙感,2011,23(3):77-81.doi:10.6046/gtzyyg.2011.03.14.

Yu H Y,Cheng G,Zhang Y M,et al.The detection of earthquake-caused collapsed building information from LiDAR data and aerophotograph[J].Remote Sensing for Land and Resources,2011,23(3):77-81.doi:10.6046/gtzyyg.2011.03.14.

[5]魏永明,魏顯虎,陳玉.岷江流域映秀——茂縣段地震次生地質災害分布規律及發展趨勢分析[J].國土資源遙感,2014,26(4):179-186.doi:10.6046/gtzyyg.2014.04.28.

Wei Y M,Wei X H,Chen Y.Analysis of distribution regularity and development tendency of earthquake secondary geohazards in Yingxiu-Maoxian section along the Minjiang River[J].Remote Sensing for Land and Resources,2014,26(4):179-186.doi:10.6046/gtzyyg.2014.04.28.

[6]姚艷麗,蔣勝平,王紅平.基于地面三維激光掃描儀的滑坡整體變形監測方法[J].測繪地理信息,2014,39(1):50-53,56.

Yao Y L,Jiang S P,Wang H P.Overall deformation monitoring for landslide by using ground 3D laser scanner[J].Journal of Geomatics,2014,39(1):50-53,56.

[7]袁小祥,王曉青,竇愛霞,等.基于地面LiDAR玉樹地震地表破裂的三維建模分析[J].地震地質,2012,34(1):39-46.

Yuan X X,Wang X Q,Dou A X,et al.Terrestrial LiDAR-based 3D modeling analysis of surface rupture caused by Yushu earthquake[J].Seismology and Geology,2012,34(1):39-46.

[8]董秀軍,黃潤秋.三維激光掃描技術在高陡邊坡地質調查中的應用[J].巖石力學與工程學報,2006,25(s2):3629-3635.

Dong X J,Huang R Q.Application of 3D laser scanning technology to geologic survey of high and steep slope[J].Chinese Journal of Rock Mechanics and Engineering,2006,25(s2):3629-3635.

[9]梅文勝,周燕芳,周俊.基于地面三維激光掃描的精細地形測繪[J].測繪通報,2010(1):53-56.

Mei W S,Zhou Y F,Zhou J.Fine topographic mapping based on ground three-dimensional laser scanning[J].Bulletin of Surveying and Mapping,2010(1):53-56.

[10]孫柏濤,梅云弟.豎向變剛度約束砌體結構抗震能力研究[J].地震工程與工程振動,2012,32(2):1-13.

Sun B T,Mei Y D.A study on the aseismic capacity of the confined masonry structure with variable vertical stiffness [J].Earthquake Engineering and Engineering Vibration,2012,32(2):1-13.

[11]任俊杰,張世民,馬保起,等.龍門山斷裂帶中北段大震復發特征與復發間隔估計[J].地震學報,2009,31(2):160-171.

Ren J J,Zhang S M,Ma B Q,et al.Characteristics and recurrence intervals of large earthquakes along the middle-northern segment of the Longmenshan fault zone[J].Acta Seismologica Sinica,2009,31(2):160-171.

[12]陳浩,李勇,董順利,等.汶川8.0級地震地表破裂白鹿鎮段的變形特征[J].自然雜志,2009,31(5):268-271.

Chen H,Li Y,Dong S L,et al.Deformation features of surface rupture of the Wenchuan Earthquake(M8.0) at the Bailu Town[J].Chinese Journal of Nature,2009,31(5):268-271.

[13]陳立春,冉勇康,陳杰,等.龍門山前山斷裂大地震原地重復發生的地質地貌證據[J].第四紀研究,2009,29(3):472-483.

Chen L C,Ran Y K,Chen J,et al.Geological and geomorphic evidence for repeatedly occurring severe earthquakes on fronter-range fault of Longmenshan,China[J].Quaternary Sciences,2009,29(3):472-483.

[14]周強,孫柏濤.汶川地震中極震區砌體結構教學樓典型震害分析[J].震災防御技術,2010,5(4):467-476.

Zhou Q,Sun B T.Seismic damage analysis of typical masonry teaching buildings in meizoseismic area of Wenchuan earthquake[J].Technology for Earthquake Disaster Prevention,2010,5(4):467-476.

[15]謝瑞,胡敏捷,程效軍,等.三維激光HDS3000掃描儀點位精度分析與研究[J].遙感信息,2008(6):53-57,84.

Xie R,Hu M J,Cheng X J,et al.Study on point precision of 3D laser scanning based on HDS3000[J].Remote Sensing Information,2008(6):53-57,84.

(責任編輯： 刁淑娟)

Typical earthquake damage extraction and three-dimensional modeling analysis based on terrestrial laser scanning: A case study of Bailu middle school of Pengzhou city

JIAO Qisong1,2, ZHANG Jingfa2, JIANG Hongbo2， SU Yuanyuan2, WANG Xu2

(1.InstituteofEngineeringMechanics,CEA,Harbin150080,China; 2.InstituteofCrustalDynamics,CEA,Beijing100085,China)

Abstract:Traditional optical and Radar remote sensing technology fails to acquire the building facade damage in the case of earthquake damage assessment. The capability of identifying the buildings whose structure has been destroyed and whose roof top is basically complete is poor. This reduces the assessment accuracy. However, terrestrial laser scanning could acquire the high precision 3D information of the damaged structure. Using RIEGL VZ-1000 laser scanner, the authors collected laser point clouds data of Bailu middle school on October 23, 2013. The point clouds were used to conduct 3D modeling and analysis through triangular irregular network modeling and texture mapping. On such a basis, a detailed modeling study was carried out on the typical architectural damage characteristics and surface rupture destruction. The authors achieved damaged buildings surface observation at any angle and quantitative measurement of the surface rupture through point clouds 3D modeling. The results truly reproduced the current scene of Bailu middle school after earthquake and permanently preserved the damage information. The results showed that the Qinxue building located on the hanging wall had no significant damage, whereas the Qiuzhi building located on the heading wall was seriously damaged. “X”-shaped conjugate shear fracture, structural column fracture, rebar exposure and other phenomena were widespread. The earthquake scarp settlement was small and the scarp elevation profile showed that the minimum and maximum height was 205cm and 231cm respectively. This was consistent with the GPS measurements.

Keywords:Bailu middle school; terrestrial laser scanning(TLS); earthquake damage extraction; 3D modeling

作者簡介：第一 焦其松(1985-)，男，博士研究生，主要從事遙感地質方面的研究。Email： 1985jqs@163.com。

中圖法分類號：TP 79

文獻標志碼：A

文章編號：1001-070X(2016)01-0166-06

基金項目：中央級公益性科研院所基本科研業務專項(編號： J2213802)和國家高技術研究發展計劃(編號： 2012AA121304)共同資助。

收稿日期：2014-09-09；

修訂日期：2014-11-02

doi:10.6046/gtzyyg.2016.01.24

引用格式： 焦其松,張景發,蔣洪波,等.基于TLS技術的典型建筑物震害信息三維建模分析——以彭州市白鹿中學為例[J].國土資源遙感，2016，28(1)：166-171.(Jiao Q S,Zhang J F,Jiang H B,et al.Typical earthquake damage extraction and three-dimensional modeling analysis based on terrestrial laser scanning:A case study of Bailu middle school of Pengzhou City[J].Remote Sensing for Land and Resources,2016,28(1)：166-171.)