潘 越， 馬 寧， 吳向上， 孟慶良

(1.上海無線電設(shè)備研究所，上海201109;2.上海機(jī)電工程研究所，上海201109)

0 引言

隨著橋梁行業(yè)的飛速發(fā)展，微形變監(jiān)測技術(shù)對橋梁的健康普查和災(zāi)害預(yù)防具有重要意義。與傳統(tǒng)微形變監(jiān)測方法相比，微波干涉測量法具有高精度，非接觸式遠(yuǎn)距離測量和實時測量等無可比擬的優(yōu)勢。微形變檢測是對變形體進(jìn)行測量以確定其空間位置及內(nèi)部形態(tài)的變化特征，目前微形變檢測已經(jīng)有很多相對成熟的測量方法，包括常規(guī)大地測量方法、專門測量法、攝影法等[1]。常規(guī)大地測量法并不是一種新興的測量技術(shù)，無論是測量儀器還是理論都經(jīng)過了時間的考驗，所以測量結(jié)果也相對安全可靠[2]。這種測量方法觀測費用低、適用性廣、精度高，能夠適用于不同工作環(huán)境和檢測物。但是同時也有一定的局限性，測量過程耗時久，需要大量的勞動力，導(dǎo)致測量效率低下，而且遇到雨雪等惡劣天氣還會影響測量結(jié)果的準(zhǔn)確性[3]。專門的測量方法(例如加速度傳感器)的優(yōu)點是相對精度高，而且專用儀器種類豐富、選擇空間大，其缺點是測量范圍小，安裝費時費力。攝影測量的方法可以利用少量的觀測點數(shù)來獲取觀測目標(biāo)的局部或者完整的樣貌[4]。攝影測量法的優(yōu)點是可以完整地獲取觀測目標(biāo)的變形信息，具有快速直觀全面的特點，可以減少戶外工作量，還可以遠(yuǎn)距離測量[5]。它的缺點是測量精度低，對光線要求高，需要拍攝大量照片，不適于定點測量[6]。常規(guī)的檢測方法雖然理論成熟，但大都存在一定的局限性，不能進(jìn)行大范圍短時間地測量，效率不高需要勞動量大，且檢測設(shè)備無法在惡劣天氣下正常工作等。本文所研究的微波干涉測量法，能夠?qū)h(yuǎn)處的目標(biāo)點進(jìn)行非接觸式測量，測量周期短，測量精度高且不易受天氣、溫度和環(huán)境的影響。此外，還具有便攜、安裝簡單等優(yōu)點，因此在微形變測量領(lǐng)域具有廣泛的應(yīng)用前景。

1 理論基礎(chǔ)

1.1 干涉測量技術(shù)

干涉測量技術(shù)是通過對比雷達(dá)不同時刻的兩次回波信號的相位差來獲取目標(biāo)的微形變信息，其原理如圖1所示。利用干涉測量技術(shù)，即使目標(biāo)產(chǎn)生毫米級別微小的變形也能被發(fā)現(xiàn)并測量出來。干涉測量法測量微形變的公式為[7]

式中：λ 為發(fā)射電磁波波長;φ1為被測目標(biāo)第一次回波信號的相位;φ2為被測目標(biāo)第二次回波信號的相位;d 為雷達(dá)視線方向的微形變。

圖1 干涉測量原理圖

1.2 合成孔徑雷達(dá)理論

合成孔徑雷達(dá)(Synthetic Aperture Radar，SAR)是選用小尺寸的天線向空間輻射電磁波信號，雷達(dá)沿著軌道方向運動，在運動過程中不斷地發(fā)射接收信號，并將接收到的信號按順序存儲起來[8]。真實孔徑雷達(dá)是在同一位置接收目標(biāo)的回波信號，合成孔徑雷達(dá)則是在不同位置接收同一目標(biāo)的回波信號。合成孔徑雷達(dá)天線的尺寸不需要很大，很小的天線就可以通過相對運動產(chǎn)生一個等效的大孔徑的線性陣列天線[9]。合成孔徑雷達(dá)接受到的信號不是同時采集的，所以需要將天線在各個位置發(fā)射的信號和接收到的目標(biāo)回波按順序儲存，并對一定時間間隔內(nèi)的接收信號進(jìn)行處理，完成雷達(dá)成像[10]。合成孔徑長度Ls是目標(biāo)仍在天線波瓣寬度之內(nèi)雷達(dá)運動的最大距離，合成孔徑雷達(dá)的最大合成孔徑長度Ls可以表示為

式中：λ 為天線輻射電磁波的波長;D 為天線的真實孔徑;R 為雷達(dá)與目標(biāo)之間的斜距。

1.3 后向投影算法

后向投影(Back Projection，BP)算法是一種典型的SAR 時域成像算法。延遲-求和為后向投影算法的核心思想，通過對回波信號延遲時間的計算實現(xiàn)回波信號的相干疊加，將接收到的回波信號反向投影到所需測量區(qū)域的每個像素點位置處，累加不同角度的像素值，進(jìn)而獲取對應(yīng)點的散射強(qiáng)度信息完成目標(biāo)高分辨成像[11]。

BP算法需要對成像區(qū)域中的各像素點執(zhí)行逐一匹配操作，算法第一步需要對成像區(qū)域進(jìn)行網(wǎng)格劃分，對成像空間進(jìn)行離散化操作得到M 個像素點;算法第二步對信號從發(fā)射天線到各像素點以及從各像素點返回至接收天線的時間和進(jìn)行計算，根據(jù)相關(guān)往返耗時獲取對應(yīng)位置的幅值信息，將軌道上各位置所有回波信號的幅值進(jìn)行累加得到對應(yīng)點幅值完成成像。BP 算法成像的工作流程圖如圖2所示。

后向投影法成像公式的離散形式為

式中：I(m)為第m 個像素點的散射系數(shù)，m=1，2，…，M;Sr[τ(k，m)，k]為第m 個像素點在第k 個位置處回波信號的幅度值;τ(k，m)為第m 個像素點到第k 個位置的雙程時延。

2 模型搭建

2.1 測量模型

橋梁的微形變測量主要是測量橋面垂直于地面方向上產(chǎn)生的微小位移，實際工程中很難直接測得，而是通過測量雷達(dá)與目標(biāo)連線方向上的位移變化量計算出橋面垂直方向上位移。測量模型如圖3所示，雷達(dá)與地面之間的夾角是α，目標(biāo)與

圖2 BP算法流程圖

雷達(dá)的相對距離為R，雷達(dá)與目標(biāo)連線方向上產(chǎn)生的變化量是ΔR，進(jìn)而計算出橋面垂直方向產(chǎn)生的位移為

圖3 微波干涉測量法測量模型

2.2 信號模型

步進(jìn)頻率信號是一系列窄帶脈沖組成的合成寬帶信號，能夠獲得高分辨率，在降低雷達(dá)體積、重量和成本方面具有顯著作用[12]。步進(jìn)頻率信號的實質(zhì)是線性調(diào)頻信號的離散化形式，是以一定步長連續(xù)變化的N 個正弦(或余弦)信號序列，主要分為脈沖和連續(xù)波兩種信號形式。脈沖信號是在每個頻率上發(fā)射一個窄帶脈沖信號，脈沖信號體制下信號源的結(jié)構(gòu)相對復(fù)雜，但其脈沖本身具有解算距離的能力，能夠求解距離模糊，所以可以測量更遠(yuǎn)的距離。連續(xù)波信號是在頻率變化之間信號持續(xù)發(fā)射，信號能量大，能夠降低對發(fā)射機(jī)峰值功率的要求，但不能計算距離模糊，作用距離相對較近[13]。步進(jìn)頻率連續(xù)波的示意圖如圖4所示。

圖4 步進(jìn)頻率連續(xù)波示意圖

步進(jìn)頻率信號的一般形式為

式中：Tr為脈沖重復(fù)周期;τ 為發(fā)射脈沖寬度;Δf為步進(jìn)頻率量;N 為頻率步進(jìn)數(shù)。

雷達(dá)在長度為L 的軌道上勻速運動，并按一定時間間隔向目標(biāo)區(qū)域發(fā)射N 個脈沖的步進(jìn)頻率連續(xù)波，由于信號持續(xù)時間很短，假設(shè)雷達(dá)在發(fā)射和接收信號期間是靜止的，即雷達(dá)與目標(biāo)之間未發(fā)生相對運動。雷達(dá)每移動長度Δl 的距離發(fā)射和接收一次信號，雷達(dá)在導(dǎo)軌上走完單程發(fā)射和接收M 幀信號完成一次數(shù)據(jù)的采集，然后雷達(dá)重新返回起點重復(fù)上述運動，完成多次數(shù)據(jù)的收集和存儲經(jīng)過信號處理后按照一定的成像算法分別成像。發(fā)射的步進(jìn)頻率連續(xù)波波形圖如圖5所示。

3 微波干涉測量法

基于步進(jìn)頻率連續(xù)波的微波干涉測量法具有高精度，非接觸式測量和實時測量的優(yōu)勢。雷達(dá)在固定的水平軌道上勻速運動，通過持續(xù)向外輻射電磁波信號獲取目標(biāo)點的位置信息。距離向利用電磁波從雷達(dá)到目標(biāo)的用時來計算二者的相對距離。方位向是利用合成孔徑獲取目標(biāo)高分辨率的方位向信息。確定距離信息與方位信息后，對已知目標(biāo)點進(jìn)行微形變的測量，利用干涉測量原理計算多次測量的回波信號的相位差來獲取目標(biāo)點的微形變。

圖5 發(fā)射信號模型

圖6 后向投影法

3.1 目標(biāo)定位

對目標(biāo)點微形變的測量首先要對目標(biāo)的定位，本文選用BP 算法對目標(biāo)進(jìn)行成像確定目標(biāo)的位置信息。首先對成像目標(biāo)區(qū)域進(jìn)行離散化處理，將它分成P×Q 個均勻的像素點，P、Q 分別代表雷達(dá)視線方向和沿軌道運動方向像素點的個數(shù)。設(shè)合成孔徑雷達(dá)的角度分辨率為Δθ，距離分辨率為Δr，軌道長度為L，在軌道上一共有m 個位置進(jìn)行數(shù)據(jù)的采集。雷達(dá)在每個位置上都進(jìn)行一次信號的發(fā)射和接收，采集到的數(shù)據(jù)可以看作是時間陣列，每個陣元的位置記為m，其中目標(biāo)區(qū)域中任意像素點(p，q)與第m 個位置之間的相對距離記作Rm。步進(jìn)頻率連續(xù)波以固定的重復(fù)周期沿方位向發(fā)射接收信號，在第m 個位置發(fā)射的信號形式為

從上述成像過程可以看出先對回波信號進(jìn)行傅立葉逆變換完成距離向聚焦，然后再利用合成孔徑后向投影算法完成方位向聚焦，最終完成目標(biāo)點的二維成像。

3.2 微形變測量

完成對目標(biāo)點的成像后，我們可以通過二維圖像上確定測量點的位置，然后利用干涉測量法同時對多點進(jìn)行微形變檢測，目標(biāo)點的微形變測量過程如圖7所示。設(shè)兩次接收到的回波信號表示為

接收信號的相位信息由往返路徑確定，可得

第一次回波信號與第二次回波信號的共軛相乘

提取干涉相位，得

從而得到微形變量為

最后利用公式(4)計算出橋面垂直方向的微小位移d。

圖7 二維點目標(biāo)微形變測量過程

4 仿真分析

發(fā)射一串含有256個單頻脈沖的步進(jìn)頻率連續(xù)波信號，每個脈沖信號持續(xù)時間為1μs，起始頻率為17.1 GHz，其單頻脈沖間的步進(jìn)頻率為1 MHz。設(shè)雷達(dá)在長度為2 m 的軌道上作勻速直線運動，每隔0.005 m 完成一次信號的發(fā)射和接收，2 min走完一次全程。下面在信噪比為20 dB，三個目標(biāo)點的位置坐標(biāo)分別為(100，4)，(100，10)，(100，16)的條件下，對多目標(biāo)點的二維成像進(jìn)行仿真。得到目標(biāo)點成像的仿真結(jié)果如圖8所示。

圖8 目標(biāo)點的成像結(jié)果

由仿真圖8可知，在停-走-停假設(shè)條件下，合成孔徑雷達(dá)在固定導(dǎo)軌上勻速運動，以固定頻率發(fā)射和接收步進(jìn)頻率連續(xù)波，完成數(shù)據(jù)的采集，然后對二維回波信號進(jìn)行信號處理，進(jìn)而能夠很好地完成合成孔徑雷達(dá)二維成像，能夠很好達(dá)到對目標(biāo)點定位的目的。

4.1 微形變測量仿真分析

在目標(biāo)位置和雷達(dá)仰角確定的前提下，可以利用干涉測量對目標(biāo)點進(jìn)行長時間實時測量。假設(shè)雷達(dá)的一次測量周期為3 min，目標(biāo)點在距離向每隔3 min發(fā)生0.2 mm 的微形變，下面對目標(biāo)點60 min內(nèi)形變量的測量進(jìn)行仿真。將測量的微形變值與理論微形變值進(jìn)行比較，得到的測量誤差如圖9所示。

仿真圖9反映的是目標(biāo)點微形變測量值與實際微形變值的差值，從圖中可知絕大部分點的測量誤差范圍在(-0.1 mm，0.1 mm)，證明了干涉測量的準(zhǔn)確性和高精度，精度能夠達(dá)到0.1毫米級。

4.2 信噪比對測量精度的影響

圖9 二維干涉測量微形變測量誤差

通過對微波干涉測量算法的仿真，已經(jīng)證實了設(shè)計的基于步進(jìn)頻率連續(xù)波的微波干涉測量法能夠有效準(zhǔn)確地測量出多目標(biāo)點的微形變，在20 d B信噪比的情況下，測量精度能夠達(dá)到0.1毫米級。下面就針對信噪比對多目標(biāo)點微形變測量精度的影響進(jìn)行討論。仿真參數(shù)設(shè)置與表1相同，分別在信噪比為5，10，15，20，25，30，35，40，45，50 d B，的情況下，對每一個點目標(biāo)重復(fù)300次微形變的測量，然后求出其均方根誤差。圖10給出了不同信噪比情況下，三個目標(biāo)點微形變測量值的均方根誤差。

圖10 信噪比對點目標(biāo)測量精度的影響

由仿真結(jié)果可知，隨著信噪比的增大，不同目標(biāo)點微形變測量的誤差逐漸減小，且不同點處的測量精度不同;當(dāng)信噪比大于30 dB時，三個目標(biāo)點的測量誤差相差不大;綜上所述，信噪比雖然在一定范圍內(nèi)能夠影響微形變檢測算法的測量精度，且對不同位置處目標(biāo)點測量精度的影響各不相同，但只要在信噪比大于15 dB 的環(huán)境下仍然能夠保證對多目標(biāo)點微形變測量的精度在亞毫米級，證明了該算法具有很高的測量精度，能夠有效測量目標(biāo)的微形變。

5 結(jié)論

本文分析了微形變檢測在橋梁微形變監(jiān)測領(lǐng)域的優(yōu)勢和應(yīng)用價值，利用微波干涉測量法對目標(biāo)微形變進(jìn)行高精度和遠(yuǎn)距離的測量，并通過Matlab仿真驗證了算法的可行性與準(zhǔn)確性。首先利用后向投影算法對目標(biāo)點進(jìn)行定位，然后利用干涉測量技術(shù)計算出微小位移，最后通過比較分析不同信噪比對測量精度的影響，當(dāng)信噪比大于15 dB時該算法測量精度能達(dá)到0.1 mm，說明本文提出的基于步進(jìn)頻率連續(xù)波的微波干涉測量法能夠有效的測量計算出目標(biāo)點微小位移。