熊黎黎,賈 璐

(1. 南昌航空大學(xué)土木建筑學(xué)院,江西 南昌 330069;2. 南昌大學(xué)建筑工程學(xué)院,江西 南昌 330031)

1 引言

工人在建筑樓房時(shí),地基的穩(wěn)固是必不可少的一部分,但露天開采會(huì)破壞地基地表平衡,導(dǎo)致地基邊緣位置穩(wěn)定性差,極易出現(xiàn)建筑坍塌、建筑質(zhì)量低、穩(wěn)固性差的問(wèn)題,這種現(xiàn)象對(duì)工作人員及設(shè)備都會(huì)產(chǎn)生安全影響。而我國(guó)地質(zhì)結(jié)構(gòu)復(fù)雜、地形地貌起伏較大,南方地區(qū)容易出現(xiàn)地質(zhì)災(zāi)害,導(dǎo)致建筑地基出現(xiàn)二次變形的現(xiàn)象[1-2],地基巖土邊坡的失穩(wěn)或滑動(dòng)給工程建設(shè)及人民生命財(cái)產(chǎn)帶來(lái)嚴(yán)重?fù)p失。為了保證人民群眾的生命安全,需要對(duì)建筑地基二次形變檢測(cè)方法展開研究,從而避免地質(zhì)災(zāi)害造成的危害。

王軍飛等人[3]提出地基合成孔徑雷達(dá)的邊坡形變監(jiān)測(cè)應(yīng)用研究方法,利用粗差探測(cè)技術(shù)對(duì)PS點(diǎn)解纏處理,通過(guò)對(duì)邊坡的監(jiān)測(cè)數(shù)據(jù),取得PS點(diǎn)在監(jiān)測(cè)期間的時(shí)間序列形變,以此實(shí)現(xiàn)地基二次形變檢測(cè),該方法的確立結(jié)果存有誤差,存在檢測(cè)效果差的問(wèn)題。宋強(qiáng)等人[4]提出基于匹配點(diǎn)變形估計(jì)的布標(biāo)缺陷檢測(cè)算法,通過(guò)獲取的地基形變圖像,利用光流法計(jì)算出圖像與待檢測(cè)圖像之間的匹配點(diǎn)運(yùn)動(dòng)適量,通過(guò)逆變換去除待檢測(cè)圖像的受力形變,從而實(shí)現(xiàn)二次形變檢測(cè),該方法獲取的圖像不完善,存在檢測(cè)精度低的問(wèn)題。劉俊博等人[5]提出基于激光點(diǎn)云的鐵路邊坡表面形變檢測(cè)方法,通過(guò)檢測(cè)雷達(dá)裝置與三維點(diǎn)云數(shù)據(jù)采集系統(tǒng),獲得鐵路周邊環(huán)境的三維點(diǎn)云。與鐵路場(chǎng)景先驗(yàn)知識(shí)配準(zhǔn)算法相結(jié)合,通過(guò)獲取基準(zhǔn)點(diǎn)云數(shù)據(jù),計(jì)算出地基邊坡表面變化量,實(shí)現(xiàn)形變檢測(cè)。該方法的處理結(jié)果不妥當(dāng),存在檢測(cè)可行性低的問(wèn)題。

為了解決上述方法中存在的問(wèn)題,提出基于神經(jīng)網(wǎng)絡(luò)的建筑地基二次形變檢測(cè)方法。通過(guò)采集與處理建筑地基二次形變數(shù)據(jù),獲得二次形變差異圖,利用NSCT方法融合差異圖像,并提取其特征向量,通過(guò)BP神經(jīng)網(wǎng)絡(luò)實(shí)現(xiàn)建筑地基二次形變檢測(cè)。

2 數(shù)據(jù)采集及處理

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

由于建筑地形內(nèi)部結(jié)構(gòu)復(fù)雜,為此以激光雷達(dá)、慣性測(cè)量單元、GPS地形掃描系統(tǒng)作為采集設(shè)備[6]。根據(jù)選取的采集設(shè)備,分別對(duì)激光雷達(dá)、慣性測(cè)量單元及GPS接收機(jī)的基準(zhǔn)時(shí)鐘表述:

T=ts=tc=tg

(1)

式中,T表述基準(zhǔn)時(shí)鐘,ts表述激光雷達(dá)時(shí)間參量,tc記作慣性測(cè)量單位時(shí)間參量,tg表述GPS時(shí)間參量。

根據(jù)設(shè)定結(jié)果,選取1個(gè)控制器調(diào)節(jié)采集系統(tǒng)的多線程運(yùn)行模式,從而確保數(shù)據(jù)采集有效性。分別設(shè)置各個(gè)采集設(shè)備、采集目標(biāo),利用下式表述激光雷達(dá)的采集目標(biāo),定義如下:

s={t,sx,sy,sz}

(2)

式中,s表述采集的地基數(shù)據(jù)目標(biāo),t表述采集時(shí)間,sx、sy、sz均表述x、y、z三個(gè)軸向的地形地基位置信息。

慣性測(cè)量單元采集的地基數(shù)據(jù)目標(biāo)定義:

c={t,cy,ce,ch}

(3)

式中,c表述采集的目標(biāo)數(shù)據(jù),cy表述仰角信息,ce表述側(cè)角信息、ch均表述航向信息。

GPS的建筑地基采集數(shù)據(jù)即可定義為:

g={t,gj,gw,gc}

(4)

式中,g表述GPS采集數(shù)據(jù),gj表述地基經(jīng)度信息,gw表述緯度信息,gc表述高程信息。

選取的地形掃描裝置設(shè)備的實(shí)際輸出值會(huì)影響建筑地基形變數(shù)據(jù)采集效果,為確保數(shù)據(jù)采集可靠性,需要利用控制系統(tǒng)改進(jìn)采集設(shè)備,因而設(shè)立的控制邏輯表示如下:

K=(Ti,Pis,Pic,Pig)

(5)

式中,K表述多線程控制指令,Ti表述參量設(shè)置,Pis表述激光雷達(dá)采樣效率,Pic表述慣性測(cè)量單元采樣效率,Pig表述GPS接收機(jī)采樣效率。

通過(guò)單一指令控制采集設(shè)備對(duì)建筑地形形變數(shù)據(jù)的不同單元采集,以此避免滯后影響。

針對(duì)以上設(shè)置的控制邏輯,對(duì)采集的地基形變數(shù)據(jù)波動(dòng)情況分析,為此綜合激光雷達(dá)、慣性測(cè)量單位、GPS數(shù)據(jù)信息采集后的波動(dòng)情況,利用下式計(jì)算表示:

(6)

式中,ε表述波動(dòng)強(qiáng)度,st、st-1表述相鄰采集時(shí)刻所采集的數(shù)據(jù),ct、ct-1表述相鄰時(shí)刻慣性測(cè)量單元采集結(jié)果,gt、gt-1表述相鄰時(shí)刻GPS采集的數(shù)據(jù),T表述基準(zhǔn)時(shí)刻。

從方程(6)可以看出,當(dāng)ε的取值越高時(shí),說(shuō)明所采集的建筑地基二次形變數(shù)據(jù)波動(dòng)強(qiáng),證明該地形采集結(jié)果復(fù)雜;反之說(shuō)明形變穩(wěn)定。

將上述過(guò)程作為基礎(chǔ),設(shè)定數(shù)據(jù)采集過(guò)程中的參數(shù),定義如下:Ki/ε≤d。式中,Ki表述目標(biāo)指令,d表述地基形變采集數(shù)據(jù)精度。

根據(jù)設(shè)定結(jié)果,從而實(shí)現(xiàn)建筑地基二次形變數(shù)據(jù)采集。

2.2 變差異圖獲取

基于采集的建筑地基數(shù)據(jù),采用數(shù)字成像技術(shù)對(duì)采集的數(shù)據(jù)處理,獲取建筑地基二次形變差異圖?？紤]到圖像像素點(diǎn)周圍特征,利用均值對(duì)數(shù)比表述獲取的差異圖,定義如下:

(7)

式中,(i,j)表述像素點(diǎn),Xd表述差異圖,m1(i,j)表述X1圖像位于像素點(diǎn)的均值,m2(i,j)表述X2圖像的均值。

由于圖像之間的相干系數(shù)會(huì)對(duì)圖像的相似度產(chǎn)生影響,所以采用相干系數(shù)獲取成像處理后變化區(qū)域較明顯、背景信息較平滑,因而兩個(gè)圖像的相干系數(shù)定義如下:

(8)

式中,*表述共軛復(fù)數(shù),Xm(i,j)表述相關(guān)系數(shù)取值。

2.3 差異圖處理

由于獲取的差異圖中子帶圖像與原始圖像內(nèi)部細(xì)節(jié)一致,所以根據(jù)NSCT的特點(diǎn)及優(yōu)勢(shì),采用NSCT處理獲取的差異圖[7-8],達(dá)到去噪目的后融合差異圖像[9],有利于更好的捕捉圖像細(xì)節(jié),便于檢測(cè)出建筑地基二次形變。

采用拉普斯金字塔濾波器、方向組濾波器作為濾波設(shè)備,分別對(duì)差異圖Xd、Xm采取NSCT變換,取得兩種差異圖像之間的低頻系數(shù),即Ld、Lm,高頻系數(shù)Hd、Hm,從獲取的系數(shù)中反映出高低頻圖像信息。

為了使圖像的二次形變區(qū)域更加明顯,采用加權(quán)平均法進(jìn)一步處理圖像低頻系數(shù),令圖像與系數(shù)融合,那么低頻系數(shù)計(jì)算如下:

Lf=αLm+(1-α)+Ld

(9)

式中,Lf表述融合后的低頻系數(shù),α表述加權(quán)平均融合系數(shù)。

經(jīng)上述融合處理后,采用NSCT方法對(duì)Lf與差異圖Xd、Xm反變換,從中取得融合后的差異圖像Fd、Fm。采用局部能量法對(duì)背景圖像噪聲數(shù)據(jù)抑制,最終融合后的差異圖像定義如下:

(10)

式中,Ed(i,j)表述差異圖Fd的局部能量,Em(i,j)表述Fm的局部能量,Li,j表述像素中心,F(m,d)(q)表述第q個(gè)像素值,Fh表述合成后的差異圖。

2.4 差異圖特征提取

通過(guò)對(duì)差異圖像的濾波處理及融合,取得建筑地基二次形變?nèi)诤喜町悎D。采用PCA方法提取融合差異圖像的二次形變特征[10-11],根據(jù)分塊思想劃分融合差異圖,成為n個(gè)3×3的相鄰小塊,同時(shí)不具備重疊性。為此將各個(gè)小塊看成矩陣,對(duì)其轉(zhuǎn)換成列向量,利用下式計(jì)算出全部小塊均值向量。

(11)

式中,Pt表述列向量,m表述均值向量,t表述分塊數(shù)量。

重構(gòu)所有列向量,將其重構(gòu)為尺寸大小為32×n的矩陣,定義為S,因而利用下式計(jì)算S的協(xié)方差矩陣,表示如下:

(12)

式中,T表述矩陣轉(zhuǎn)置。

由于協(xié)方差矩陣C的尺寸為32×32,對(duì)其特征分解后獲取特征值、特征向量,從中挑選出相應(yīng)特征向量后構(gòu)成正交基。因而向量構(gòu)成的矩陣由A,該矩陣A的大小尺寸在32×r。

依據(jù)融合后的差異圖Fh,挑選出圖像各個(gè)像素的鄰域小塊,對(duì)小塊轉(zhuǎn)換后取得列向量,即Yη,其中η=1,2,…,k,k表述圖像總數(shù)量。

將列向量映射到矩陣A中,其定義如下:

Yη=AT(Vη-m)

(13)

式中,Yη表述提取的特征向量,AT表述轉(zhuǎn)置矩陣。

融合差異圖中的各個(gè)像素都可以通過(guò)R維特征向量表述,為此根據(jù)圖像像素點(diǎn)[Y1,Y2,…,Yk]組建成特征向量空間Q,以此實(shí)現(xiàn)建筑地基差異圖形變特征向量的提取。

3 建筑地基二次形變檢測(cè)

根據(jù)提取的建筑地基差異圖二次形變特征向量,將其輸入到建立的神經(jīng)網(wǎng)絡(luò)中訓(xùn)練[12-13],從而實(shí)現(xiàn)建筑地基二次形變檢測(cè)。

結(jié)合小波神經(jīng)網(wǎng)絡(luò)思想,改進(jìn)激勵(lì)函數(shù),為此提取出圖像相位信息:

h(x)=cos(5x)·e-x2/2

(14)

式中,h(x)表述相位信息,e表述誤差。

設(shè)置網(wǎng)絡(luò)輸入數(shù)量由R表述,節(jié)點(diǎn)輸出由N表述,那么神經(jīng)網(wǎng)絡(luò)中第p個(gè)樣本的第n個(gè)節(jié)點(diǎn)的輸入用下式定義:

(15)

式中,k表述隱含層單元數(shù)量,M表述輸入層單元數(shù)量,ωn,k表述各個(gè)單元之間的連接權(quán)值,ak表述伸縮因子,bk表述平移因子,SP(xm)表述輸入信號(hào),fn,p表述輸入期望。

為此輸出層函數(shù)通過(guò)下式定義:

(16)

式中,fz表述輸出層,e-z表述輸出層函數(shù)。

根據(jù)設(shè)置結(jié)果,將形變特征輸入到網(wǎng)絡(luò)內(nèi)訓(xùn)練從而實(shí)現(xiàn)建筑地基二次形變檢測(cè),其具體流程如下所示:

1)對(duì)提取的特征樣本數(shù)據(jù)預(yù)處理[14-15],用下式定義:

(17)

2)采用主成分分析法分析提取的特征向量,并確定神經(jīng)網(wǎng)絡(luò)輸入節(jié)點(diǎn)。

3)初始化網(wǎng)絡(luò)設(shè)置相關(guān)參數(shù),即學(xué)習(xí)效率a、動(dòng)量系數(shù)b、迭代誤差e。

4)向網(wǎng)絡(luò)中輸入一組學(xué)習(xí)樣本特征向量,設(shè)置輸入節(jié)點(diǎn)為輸入向量,并將fn,p輸入其中。

5)根據(jù)特征向量在網(wǎng)絡(luò)中的自學(xué)習(xí)流程,基于神經(jīng)網(wǎng)絡(luò)當(dāng)前參數(shù),進(jìn)一步計(jì)算網(wǎng)絡(luò)的實(shí)際輸出。當(dāng)輸出誤差小于e時(shí),則訓(xùn)練結(jié)束。

6)若不滿足上述條件,進(jìn)一步實(shí)施誤差反向傳播,進(jìn)一步修正神經(jīng)網(wǎng)絡(luò)參數(shù)。

7)利用修正后的神經(jīng)網(wǎng)絡(luò)再次訓(xùn)練形變特征向量,其輸出結(jié)果就是最終檢測(cè)結(jié)果。

基于上述流程,實(shí)現(xiàn)基于神經(jīng)網(wǎng)絡(luò)的建筑地基二次形變檢測(cè)。

4 實(shí)驗(yàn)與分析

為了驗(yàn)證基于神經(jīng)網(wǎng)絡(luò)的建筑地基二次形變檢測(cè)方法的整體有效性,需要對(duì)該方法開展實(shí)驗(yàn)對(duì)比測(cè)試。

采用基于神經(jīng)網(wǎng)絡(luò)的建筑地基二次形變檢測(cè)方法(方法1)、地基合成孔徑雷達(dá)的邊坡形變監(jiān)測(cè)應(yīng)用研究方法(方法2)和基于激光點(diǎn)云的鐵路邊坡表面形變檢測(cè)方法(方法3)實(shí)施實(shí)驗(yàn)測(cè)試。

1)濾波效果

開展建筑地基二次形變檢測(cè)前,為保證建筑地基二次形變檢測(cè)可靠性,將某城市施工場(chǎng)地建筑地基作為本次實(shí)驗(yàn)場(chǎng)景,詳情如圖1所示。

圖1 建筑地基場(chǎng)景示意圖

采用地基SAR系統(tǒng)工作,并按照每20分鐘采集一張建筑地基黑白圖像,且相鄰的圖像數(shù)據(jù)采集間隔時(shí)間為5分鐘,連續(xù)采集5天即可獲取建筑地基序列圖像,有利于檢測(cè)建筑地基二次變形。那么設(shè)定地基SAR系統(tǒng)工作參數(shù),詳情如表1所示。

表1 地基SAR系統(tǒng)工作參數(shù)

由于采集的建筑地基圖像內(nèi)會(huì)因?yàn)樵肼曈绊懼率箞D像相位點(diǎn)出現(xiàn)不連續(xù)的斷點(diǎn),導(dǎo)致相位解纏失敗,影響建筑地基二次形變檢測(cè)效果,所以為了可以獲取高質(zhì)量圖像,降低殘差點(diǎn)數(shù)目,需要對(duì)采集的圖像濾波處理,從而達(dá)到抑制噪聲的目的。

那么利用方法1、方法2和方法3分別對(duì)具有噪聲的建筑地基圖像濾波處理,三種方法的濾波效果測(cè)試如圖2所示。

圖2 三種方法的濾波效果測(cè)試

從圖2的實(shí)驗(yàn)結(jié)果可知,方法1對(duì)含噪聲的建筑地基圖像濾波處理后,圖像更能體現(xiàn)出內(nèi)部細(xì)節(jié)信息,驗(yàn)證了該方法的濾波性能強(qiáng)。而方法2、方法3在濾波處理時(shí),處理后圖像內(nèi)部依然存有較多噪聲,同時(shí)方法2的處理結(jié)果與原始圖像相同,說(shuō)明方法2的濾波效果不佳,也因此證明方法2和方法3的濾波效果要低于方法1。

2)形變位移檢測(cè)

在上述研究的基礎(chǔ)上,將Visual Studio.NET用作開發(fā)平臺(tái),C++用作開發(fā)語(yǔ)言,SQL Server作為數(shù)據(jù)庫(kù)平臺(tái),根據(jù)連續(xù)五天采集的建筑地基圖像數(shù)據(jù),結(jié)合貝葉斯估計(jì)方法,分別利用方法1、方法2和方法3對(duì)建筑地基二次形變開展詳細(xì)檢測(cè),測(cè)試結(jié)果如圖3所示。

圖3 建筑地基二次形變檢測(cè)測(cè)試

分析圖3中的數(shù)據(jù)發(fā)現(xiàn),選取的采集時(shí)間分別為10月11～10月15日,根據(jù)檢測(cè)結(jié)果可知,方法1的位移測(cè)試結(jié)果與實(shí)際位移結(jié)果基本一致,且運(yùn)動(dòng)軌跡相同,表明方法1的檢測(cè)精度最佳。而方法2的在10月11日～10月13日時(shí),運(yùn)動(dòng)軌跡與實(shí)際結(jié)果相同,但兩者之間存有誤差,同時(shí)方法2在10月13日～10月14日時(shí)運(yùn)動(dòng)軌跡發(fā)生轉(zhuǎn)變,偏離實(shí)際軌跡,可見(jiàn)方法2的檢測(cè)效果不佳。方法3的整體檢測(cè)結(jié)果與實(shí)際軌跡相差大,說(shuō)明方法3與實(shí)際結(jié)果偏差高。

綜上所述,方法1的檢測(cè)結(jié)果明顯優(yōu)于其余兩種方法,這是因?yàn)榉椒?對(duì)獲取的圖像實(shí)施濾波處理,以此消除圖像內(nèi)的噪聲,更有利于捕捉圖像細(xì)節(jié),從而提升了該方法的建筑地基二次形變檢測(cè)效果,使其具備更高的檢測(cè)精度。

5 結(jié)束語(yǔ)

建筑地基二次形變導(dǎo)致建筑邊緣位置穩(wěn)定性差,為有效避免這種問(wèn)題,提出基于神經(jīng)網(wǎng)絡(luò)的建筑地基二次形變檢測(cè)方法。

該方法采集了建筑地基二次形變數(shù)據(jù),對(duì)采集結(jié)果成像處理后,獲取建筑地基二次形變差異圖,通過(guò)對(duì)圖像的濾波處理及特征向量的提取,將其輸入到建立的BP神經(jīng)網(wǎng)絡(luò)內(nèi)訓(xùn)練,其輸出結(jié)果就是最終檢測(cè)結(jié)果。實(shí)驗(yàn)表明所提方法的濾波性能強(qiáng),且位移檢測(cè)結(jié)果與實(shí)際值基本一致,檢測(cè)精度高,該方法在建筑地基二次形變檢測(cè)方法中發(fā)揮著重要作用。