費(fèi)畢剛，韓鵬飛，朱茂，萬小莉，王春偉，孫珍珍

(1.國質(zhì)(北京)建設(shè)工程檢測(cè)鑒定中心，北京 100081； 2.北京東方至遠(yuǎn)科技股份有限公司，北京 100081)

1 引 言

截至2017年，我國既有房屋建筑總面積超過500億平方米，每年新建的房屋面積占到世界總量的50%。房屋建筑的安全隱患可造成其產(chǎn)生傾斜或裂縫，甚或危及安全，從而導(dǎo)致資源浪費(fèi)、環(huán)境危害、生命財(cái)產(chǎn)損失等。確保房屋建筑安全正常使用，研究其形變機(jī)理，對(duì)房屋建筑形變進(jìn)行監(jiān)測(cè)、檢測(cè)、鑒定一直是政府部門、高校和科研機(jī)構(gòu)、企業(yè)和工程技術(shù)人員努力推進(jìn)的重要工作。

隨著星載合成孔徑雷達(dá)(InSAR)技術(shù)的迅猛發(fā)展，其在地震監(jiān)測(cè)、地表各類地物形變等監(jiān)測(cè)中得到了重視和應(yīng)用[1～4]。與傳統(tǒng)測(cè)量手段，比如水準(zhǔn)儀、近景攝影等測(cè)量技術(shù)相比，其綜合優(yōu)勢(shì)明顯：形變監(jiān)測(cè)精度可達(dá)毫米級(jí)[5，6]；監(jiān)測(cè)范圍一次可達(dá)成百上千平方公里；可定量獲得建筑物時(shí)間序列形變信息，數(shù)據(jù)更新快且豐富；受天氣影響小，可全天候、全天時(shí)地連續(xù)長期監(jiān)測(cè)[7]；監(jiān)測(cè)實(shí)施便利和安全，一般僅需衛(wèi)星重復(fù)獲取地表影像，無須直接接觸目標(biāo)；成本相對(duì)低，無須維護(hù)觀測(cè)網(wǎng)，對(duì)廣范圍建筑物監(jiān)測(cè)服務(wù)好。

近5年來，許多研究人員對(duì)InSAR監(jiān)測(cè)建筑物及基礎(chǔ)設(shè)施的形變進(jìn)行了研究和應(yīng)用。趙億[8]等以渤海大樓為例，與水準(zhǔn)測(cè)量值對(duì)比，提出星載InSAR技術(shù)可達(dá)到建筑物三級(jí)測(cè)量精度指標(biāo)。Ozden[9]等研究說明星載InSAR技術(shù)可用于路面和基礎(chǔ)設(shè)施監(jiān)測(cè)管理，形變量和形變速率可達(dá)毫米級(jí)；楊魁[10]等以天津市不同類型建筑物為例，從密度、點(diǎn)位分布、質(zhì)量等方面分析了星載PS-InSAR的識(shí)別優(yōu)勢(shì)，并驗(yàn)證其有效性。張躍[11]等研究了城市場(chǎng)景下高分辨率InSAR相干系數(shù)圖中的統(tǒng)計(jì)特性，并以建筑物、樹木等作為目標(biāo)，對(duì)比選取了較高精度的統(tǒng)計(jì)模型，利用該模型可以較好實(shí)現(xiàn)建筑物的檢測(cè)。Selvakumarana[12]等運(yùn)用星載SBAS-InSAR對(duì)橋梁遭受洪水沖刷坍塌風(fēng)險(xiǎn)進(jìn)行監(jiān)測(cè)，監(jiān)測(cè)結(jié)果與實(shí)際一致，InSAR可用于橋梁結(jié)構(gòu)健康監(jiān)測(cè)和風(fēng)險(xiǎn)預(yù)警預(yù)報(bào)。馬培峰[13]等指出多種星載InSAR可勝任對(duì)粵港澳大灣區(qū)實(shí)現(xiàn)多尺度地表(含建筑物)沉降監(jiān)測(cè)和精細(xì)化局部檢測(cè)；Jordan[14]等將InSAR成功應(yīng)用于監(jiān)測(cè)英國頁巖開采區(qū)域的地表變化情況，為基線監(jiān)測(cè)開創(chuàng)了先例。朱茂[6，15]等運(yùn)用星載InSAR技術(shù)對(duì)青島膠州灣跨海大橋的形變，以及深圳地鐵沿線的建筑物進(jìn)行了形變監(jiān)測(cè)，用案例分析證實(shí)了InSAR技術(shù)有能力在未來建筑物風(fēng)險(xiǎn)評(píng)估及綜合治理的過程中起到重要作用。

研究現(xiàn)狀表明：星載InSAR技術(shù)用于形變監(jiān)測(cè)，前景廣闊，適合房屋建筑目標(biāo)。InSAR技術(shù)的分析方法很多，需根據(jù)SAR影像數(shù)量、研究區(qū)域不同特點(diǎn)來選擇，其對(duì)安全監(jiān)測(cè)應(yīng)用研究有非常重要的影響。本文主要從形變測(cè)量算法、形變曲線時(shí)序分解、建筑物安全風(fēng)險(xiǎn)識(shí)別準(zhǔn)則以及大區(qū)域案例分析方面對(duì)城市房屋安全風(fēng)險(xiǎn)識(shí)別來進(jìn)行研究，使得星載InSAR技術(shù)應(yīng)用于房屋建筑安全監(jiān)測(cè)中，動(dòng)態(tài)掌握建筑物形變情況，從而實(shí)現(xiàn)大面積、廣范圍房屋建筑的安全監(jiān)測(cè)，及時(shí)主動(dòng)篩查結(jié)構(gòu)風(fēng)險(xiǎn)。

2 PSP形變測(cè)量算法

InSAR形變測(cè)量方法目前應(yīng)用最為廣泛的是二軌法。1993年Massonnet[16]等人把InSAR技術(shù)引入地震的形變監(jiān)測(cè)研究中，利用ERS-1/2的影像數(shù)據(jù)測(cè)量了地震的形變場(chǎng)，獲得了厘米級(jí)的精度。二軌法最基本思想是利用實(shí)驗(yàn)地區(qū)地表變化前后的兩幅SAR圖像生成干涉相位圖，再利用事先獲取的DEM數(shù)據(jù)模擬干涉相位圖，最后從干涉相位圖中減去模擬的干涉相位圖，就可以得到地表形變信息，那么目標(biāo)點(diǎn)在SAR視線方向的形變量△r與形變相位φdef的關(guān)系如式(1)所示，其中λ表示雷達(dá)信號(hào)的波長，由式(1)可以看出InSAR技術(shù)的形變測(cè)量精度與雷達(dá)波長相關(guān)。

(1)

隨著InSAR技術(shù)的發(fā)展，F(xiàn)erretti[17]等人提出了永久散射體合成孔徑雷達(dá)干涉測(cè)量(簡稱PS-InSAR)方法。該方法找出研究區(qū)域內(nèi)散射特性較穩(wěn)定、對(duì)雷達(dá)波反射較強(qiáng)的硬目標(biāo)，通過對(duì)它們的一系列觀測(cè)值進(jìn)行時(shí)間序列分析，對(duì)PS點(diǎn)相位的各個(gè)誤差項(xiàng)進(jìn)行建模，估計(jì)并補(bǔ)償這些誤差相位項(xiàng)，可有效提取這些目標(biāo)的形變信息，同時(shí)亦可獲得目標(biāo)的三維信息。為了進(jìn)一步提升PS InSAR算法的性能，Constatini[18]等提出了PSP(Persistent Scatterer Pair)方法，其核心思想是建立PS網(wǎng)格，對(duì)網(wǎng)格進(jìn)行不斷擴(kuò)建，在SAR數(shù)據(jù)中最終選出PS點(diǎn)集合，該算法是通過比較PS點(diǎn)對(duì)之間的相位特性來選擇PS點(diǎn)，能夠降低空間相關(guān)性誤差(如大氣相位誤差)對(duì)形變反演結(jié)果的影響。由于算法中形成PS點(diǎn)對(duì)的2個(gè)PS點(diǎn)距離較近，那么這2個(gè)PS點(diǎn)的相位差可以建立如下關(guān)系式：

(2)

式中，Ti表示第i幅SAR圖像相對(duì)于參考SAR圖像的時(shí)間，δva表示第a個(gè)PS對(duì)中兩個(gè)PS點(diǎn)的相對(duì)形變速率，Bi表示第i幅干涉圖的有效基線，δha表示第a個(gè)PS對(duì)中兩個(gè)PS點(diǎn)的相對(duì)高程差，r表示目標(biāo)點(diǎn)的斜距，θ表示下視角，εa，i表示噪聲和未建模的誤差。一般來說，如果εa，i較小，那么這個(gè)像素點(diǎn)對(duì)就可以被認(rèn)為是PS對(duì)，并稱之為PSP。為了從數(shù)學(xué)模型上判斷PSP，可以先定義第a個(gè)PSP的時(shí)間相關(guān)系數(shù)：

(3)

式中，ωa，i表示第a個(gè)像素點(diǎn)在第i幅干涉圖像中所對(duì)應(yīng)的權(quán)重值，簡單情況下可全設(shè)置為1，在處理過程中，先對(duì)未知參數(shù)δva和δha進(jìn)行最佳估計(jì)，然后計(jì)算γa。同時(shí)，設(shè)定時(shí)間相關(guān)系數(shù)閾值，當(dāng)γa大于該閾值時(shí)，這個(gè)像素點(diǎn)對(duì)可以被判定為PS對(duì)。

3 聚類點(diǎn)形變曲線時(shí)序分解

一般來說，建筑物形變是非線性的，不同時(shí)間段的形變速率可能存在差異。同時(shí)，針對(duì)剛性目標(biāo)，形變會(huì)隨溫度的變化而變化。因此，目標(biāo)點(diǎn)的形變能以周期型信號(hào)和分段線性信號(hào)建模，周期信號(hào)描述了建筑物的膨脹和收縮的規(guī)則運(yùn)動(dòng)，分段線性信號(hào)表示每個(gè)時(shí)間段內(nèi)的變形趨勢(shì)，用于識(shí)別目標(biāo)的安全風(fēng)險(xiǎn)，其數(shù)學(xué)表達(dá)式為：

(4)

式中，di表示第ti時(shí)刻InSAR測(cè)得的形變量，其中i=N0，N0+1，N0+2，…，Nm+1；vj和bj分別表示在第j分段區(qū)間內(nèi)的形變速率和恒定形變量，它們是線性形變分量的參數(shù)；A和φ分別是周期性信號(hào)的參數(shù)；εi為第ti時(shí)刻的隨機(jī)噪聲，N1，N2，…，Nm為m個(gè)斷點(diǎn)時(shí)刻。式(4)中周期性分量描述了膨脹和收縮運(yùn)動(dòng)，不同目標(biāo)的振幅不同；而分段線性分量表示每個(gè)時(shí)間分區(qū)的變形趨勢(shì)，可評(píng)估目標(biāo)的穩(wěn)定性。關(guān)鍵點(diǎn)形變曲線分解的目的是分析形變的趨勢(shì)變化，那么與溫度變化相關(guān)的形變信號(hào)可去除。

在數(shù)據(jù)處理過程中，vj，bj，A，φ及斷點(diǎn)時(shí)刻N(yùn)1，N2，…，Nm是未知參數(shù)，需要進(jìn)行估計(jì)。數(shù)據(jù)估計(jì)過程采用最小二乘方法，使得式(5)的值最小，斷點(diǎn)數(shù)采用式(6)進(jìn)行計(jì)算，處理過程可運(yùn)用動(dòng)態(tài)規(guī)劃理論提升處理效率[19]。

(5)

(6)

4 建筑物安全風(fēng)險(xiǎn)識(shí)別準(zhǔn)則

根據(jù)上述PSP形變測(cè)量算法，可以形成PS點(diǎn)集合，PS點(diǎn)集合中可以包括一個(gè)PS點(diǎn)或者多個(gè)PS點(diǎn)，空間位置臨近，且形變曲線相似的PS點(diǎn)集合可以聚合為一個(gè)聚類點(diǎn)，本文建筑物安全風(fēng)險(xiǎn)是根據(jù)建筑物上聚類點(diǎn)的形變指標(biāo)來識(shí)別。針對(duì)每一棟建筑物，當(dāng)建筑物內(nèi)無聚類點(diǎn)，則無法進(jìn)行安全風(fēng)險(xiǎn)識(shí)別；當(dāng)僅有一個(gè)聚類點(diǎn)時(shí)，可獲得建筑物的沉降量；至少含有兩個(gè)聚類點(diǎn)時(shí)，可獲得建筑物的沉降量和差異沉降量；當(dāng)至少含有一對(duì)高程相似的聚類點(diǎn)時(shí)，可獲得建筑物的沉降量、差異沉降量和傾斜率等。

基于PSP-InSAR的處理結(jié)果，可以同時(shí)提取目標(biāo)點(diǎn)的累積形變信息和近期形變信息，建筑物的異常形變可以基于這兩個(gè)方面展開。一般來說，如果累計(jì)形變和近期形變?nèi)疾怀^閾值，那么這棟建筑物被評(píng)估為穩(wěn)定；如果長期形變和近期形變有一項(xiàng)超過閾值，那么這棟建筑物應(yīng)現(xiàn)場(chǎng)踏勘；如果長期形變和近期形變?nèi)慷汲^閾值，那么這棟建筑就需要立即進(jìn)行現(xiàn)場(chǎng)監(jiān)測(cè)。目前建筑物安全風(fēng)險(xiǎn)識(shí)別相關(guān)評(píng)價(jià)標(biāo)準(zhǔn)主要有《危險(xiǎn)房屋鑒定標(biāo)準(zhǔn)》(JGJ 125-2016)、《民用建筑可靠性鑒定標(biāo)準(zhǔn)》(GB 50292-2015)、《工業(yè)建筑可靠性鑒定標(biāo)準(zhǔn)》(GB 50144-2008)和《建筑地基基礎(chǔ)設(shè)計(jì)規(guī)范》(GB 50007-2011)等。本文的建筑物形變閾值主要是綜合參考上述標(biāo)準(zhǔn)而選取。主要指標(biāo)是：當(dāng)房屋處于自然狀態(tài)時(shí)，地基沉降速率連續(xù)2個(gè)月是否大于 4 mm/月；基礎(chǔ)傾斜量是否大于2‰。本文的風(fēng)險(xiǎn)識(shí)別分級(jí)標(biāo)準(zhǔn)主要分為A、B、C、D四級(jí)，具體如表1所示。根據(jù)InSAR技術(shù)的特點(diǎn)，測(cè)得的形變量與監(jiān)測(cè)時(shí)長相關(guān)。因此，在設(shè)定閾值的過程中，也需要考慮InSAR數(shù)據(jù)處理的時(shí)間段，考慮監(jiān)測(cè)時(shí)長的閾值設(shè)定如表2所示。

分級(jí)標(biāo)準(zhǔn) 表1

形變指標(biāo)閾值 表2

5 大區(qū)域房屋安全風(fēng)險(xiǎn)識(shí)別

在北京市選取某區(qū)域，依據(jù)COSMO-SkyMed數(shù)據(jù)(32期)的InSAR形變測(cè)量結(jié)果，對(duì)該區(qū)域中的 12 272棟建筑物進(jìn)行了安全風(fēng)險(xiǎn)識(shí)別。由于數(shù)據(jù)的分辨率較高(3 m)，在超過90%的建筑物上都能獲取PS點(diǎn)。監(jiān)測(cè)時(shí)間為2011年9月～2018年8月，圖1給出了2015年6月～2018年8月目標(biāo)區(qū)域內(nèi)的建筑形變風(fēng)險(xiǎn)評(píng)級(jí)結(jié)果。結(jié)合監(jiān)測(cè)時(shí)間段內(nèi)的監(jiān)測(cè)結(jié)果，本文選取了部分風(fēng)險(xiǎn)評(píng)估等級(jí)為C級(jí)和D級(jí)的房屋，并對(duì)此進(jìn)行了實(shí)地踏勘，對(duì)有測(cè)量條件的房屋進(jìn)行了現(xiàn)場(chǎng)變形測(cè)量。下文將詳細(xì)介紹這些房屋的InSAR風(fēng)險(xiǎn)識(shí)別結(jié)果和實(shí)地踏勘驗(yàn)證結(jié)果。

圖1 目標(biāo)區(qū)域建筑形變風(fēng)險(xiǎn)評(píng)級(jí)結(jié)果

5.1 案例1

該案例房屋為磚混一層，高約4 m，房屋實(shí)景圖如圖2所示?；贗nSAR數(shù)據(jù)處理結(jié)果，在該建筑上共監(jiān)測(cè)到17個(gè)PS點(diǎn)。為了更好地揭示建筑的沉降演化歷史，對(duì)形變數(shù)據(jù)分析后，在該棟建筑上能選出反映沉降特征的2個(gè)關(guān)鍵點(diǎn)，如圖2所示，可獲取該棟建筑的沉降信息和傾斜信息，關(guān)鍵點(diǎn)的沉降曲線如圖3和圖4所示。

圖2 案例1房屋實(shí)景圖與關(guān)鍵點(diǎn)位置

圖3 案例1關(guān)鍵點(diǎn)原始沉降曲線

圖4 案例1關(guān)鍵點(diǎn)分段沉降曲線

根據(jù)測(cè)量結(jié)果，在2011年9月～2017年11月的時(shí)間范圍內(nèi)，該棟建筑最大累計(jì)傾斜量為-2.30‰，近期傾斜速率為-0.066‰/月，且以這個(gè)速率持續(xù)傾斜超過6個(gè)月，對(duì)照表2和表3的安全風(fēng)險(xiǎn)識(shí)別準(zhǔn)則，結(jié)合圖4的沉降曲線，經(jīng)過綜合分析，該棟建筑的風(fēng)險(xiǎn)評(píng)估結(jié)果為C級(jí)。根據(jù)實(shí)地踏勘結(jié)果，該院內(nèi)相鄰房屋連接處開裂，裂縫寬度為 15 mm，如圖5所示，而且存在地基不均勻沉降的現(xiàn)象。

圖5 案例1現(xiàn)場(chǎng)踏勘結(jié)果

5.2 案例2

該案例是二層磚混式小樓，實(shí)景圖如圖6所示?；贗nSAR數(shù)據(jù)處理結(jié)果，在該建筑上共監(jiān)測(cè)到26個(gè)PS點(diǎn)。為了分析建筑的沉降演化歷史，對(duì)形變數(shù)據(jù)分析后，在該棟建筑上可選出反映沉降特征的2個(gè)關(guān)鍵點(diǎn)，如圖6所示，可獲取房屋沉降信息和傾斜信息，關(guān)鍵點(diǎn)的沉降曲線如圖7和圖8所示。

圖6 案例2房屋實(shí)景圖與關(guān)鍵點(diǎn)位置

圖7 案例2關(guān)鍵點(diǎn)原始沉降曲線

圖8 案例2關(guān)鍵點(diǎn)分段沉降曲線

根據(jù)測(cè)量結(jié)果，在2011年9月～2017年11月的時(shí)間范圍內(nèi)，該棟房屋最大累計(jì)傾斜量為-3.70‰，近期傾斜速率為-0.21‰/月，且以這個(gè)速率持續(xù)傾斜為4.8個(gè)月。對(duì)照表2和表3的安全風(fēng)險(xiǎn)識(shí)別準(zhǔn)則，結(jié)合圖8的沉降曲線，經(jīng)過綜合分析，該棟建筑的風(fēng)險(xiǎn)評(píng)估結(jié)果為D級(jí)。如圖9所示，根據(jù)實(shí)地踏勘結(jié)果，院內(nèi)自建房屋傾斜嚴(yán)重，且存在多條縫寬大于 1.0 mm的斜向裂縫及豎向裂縫，依據(jù)《危險(xiǎn)房屋鑒定標(biāo)準(zhǔn)》(JGJ 125-2016)可評(píng)定為C級(jí)局部危房。

圖9 案例2現(xiàn)場(chǎng)踏勘結(jié)果

5.3 案例3

該棟建筑為磚混一層，高約3.5 m，實(shí)景圖如圖10所示?；贗nSAR數(shù)據(jù)處理結(jié)果，在該建筑上共監(jiān)測(cè)到19個(gè)PS點(diǎn)。為了分析沉降變形，對(duì)數(shù)據(jù)分析后，在該棟建筑上可選出反映沉降特征的2個(gè)關(guān)鍵點(diǎn)，如圖10所示，可獲取最大沉降信息和傾斜信息。關(guān)鍵點(diǎn)的沉降曲線如圖11和圖12所示。

圖10 案例3房屋實(shí)景圖與關(guān)鍵點(diǎn)位置

圖11 案例3關(guān)鍵點(diǎn)原始沉降曲線

圖12 案例3分段原始沉降曲線

根據(jù)測(cè)量結(jié)果，在2011年9月～2017年11月的時(shí)間范圍內(nèi)，該棟建筑最大累計(jì)傾斜量為-3.20‰，近期傾斜速率為-0.20‰/月，且以這個(gè)速率持續(xù)傾斜達(dá)4.8個(gè)月。經(jīng)過綜合分析，該棟建筑的風(fēng)險(xiǎn)評(píng)估結(jié)果為D級(jí)，需要立即現(xiàn)場(chǎng)監(jiān)測(cè)。如圖13所示，根據(jù)實(shí)地踏勘結(jié)果，該院內(nèi)自建房與正規(guī)房屋連接處開裂，院內(nèi)地基存在下沉現(xiàn)象，自建房屋墻體存在傾斜，應(yīng)加強(qiáng)監(jiān)測(cè)。

圖13 案例3現(xiàn)場(chǎng)踏勘結(jié)果

6 結(jié)論和展望

(1)聚類點(diǎn)形變曲線時(shí)序分解方法，可以在城市房屋安全風(fēng)險(xiǎn)識(shí)別中得到較好的運(yùn)用，可有效識(shí)別監(jiān)測(cè)時(shí)間段內(nèi)平均形變不大，近期形變明顯的房屋，從時(shí)間維度上對(duì)房屋異常形變進(jìn)行合理識(shí)別。

(2)對(duì)研究區(qū)域進(jìn)行現(xiàn)場(chǎng)踏勘，對(duì)比InSAR的風(fēng)險(xiǎn)識(shí)別結(jié)果與現(xiàn)場(chǎng)踏勘結(jié)果，兩者的結(jié)論總體上比較一致，從一定程度上驗(yàn)證了InSAR測(cè)量技術(shù)對(duì)房屋安全風(fēng)險(xiǎn)識(shí)別的正確性。從而說明本文給出的考慮監(jiān)測(cè)時(shí)長的閾值以及房屋分級(jí)準(zhǔn)則具有一定的合理性，可以從技術(shù)指標(biāo)層面來指導(dǎo)大范圍的城市房屋進(jìn)行安全風(fēng)險(xiǎn)識(shí)別。不過房屋風(fēng)險(xiǎn)識(shí)別準(zhǔn)則是需要考慮不同房屋種類而給出的，因此房屋安全風(fēng)險(xiǎn)識(shí)別準(zhǔn)則需要更加深入的研究，從而使得InSAR技術(shù)能夠合理、準(zhǔn)確地在大范圍房屋安全風(fēng)險(xiǎn)識(shí)別中運(yùn)用。

(3)基于COSMO—SkyMed數(shù)據(jù)的InSAR形變測(cè)量技術(shù)能在觀測(cè)區(qū)域獲得大量測(cè)量點(diǎn)，形變精度也能夠達(dá)到mm級(jí)，可以對(duì)大面積目標(biāo)區(qū)域內(nèi)中的房屋進(jìn)行長時(shí)間段、高精度的形變監(jiān)測(cè)。但是InSAR獲取數(shù)據(jù)的分辨率偏低，從而使得信息的獲取存在一定局限性。下階段的工作需要結(jié)合多重監(jiān)測(cè)技術(shù)對(duì)房屋的形變狀態(tài)開展深入分析和研究。

(4)基于星載InSAR技術(shù)對(duì)大區(qū)域房屋安全風(fēng)險(xiǎn)的識(shí)別是一種高效且低成本的方法。隨著該方法的不斷運(yùn)用和發(fā)展，可為房屋鑒定專家提供技術(shù)支持，也可合理地對(duì)城市房屋開展安全普查，實(shí)現(xiàn)風(fēng)險(xiǎn)評(píng)估和預(yù)警。