馮 晅，梁帥帥，恩和得力海，張明賀， 董澤君，周皓秋，齊嘉慧，趙瑋昌

1.吉林大學(xué)地球信息探測儀器教育部重點實驗室，長春 130026 2.吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院，長春 130026 3.河北省地礦局第一地質(zhì)大隊，河北 邯鄲 056001

0 引言

地下管道是城市的生命線，是城市中最重要的基礎(chǔ)設(shè)施，它擔(dān)負(fù)著人們正常工作和生活的責(zé)任，具有傳輸能量、傳遞信息，水源供應(yīng)和排放廢液等一系列重要的職能。但絕大多數(shù)城市地下管道分布并不是很清楚，相關(guān)部門對于管道分布圖的資料管理也不夠規(guī)范，甚至一些廢舊居民區(qū)和廢舊工業(yè)區(qū)沒有管道圖。這些問題都給城市擴張及建設(shè)造成了不小的困難，有時給人們的生產(chǎn)和生活造成困難(如停水停電等)或?qū)е率鹿?，甚至人員傷亡。管道探測不僅對城市建設(shè)、改造等起到了極為重要的作用，甚至對于整個城市的發(fā)展進(jìn)步都至關(guān)重要。因此對地下管道快速準(zhǔn)確地識別，是我們目前很有必要做的工作[1]。

探地雷達(dá)(GPR)是一種無損探測技術(shù)，通過天線向地下發(fā)射高頻電磁波，接收到的電磁波由于地下媒介電磁特性和物理性質(zhì)的差別，其路徑、強度、波形都會發(fā)生改變，以此探知地下物質(zhì)的物性差別[2]。

20世紀(jì)90年代出現(xiàn)了研究探地雷達(dá)的熱潮，這是我國探地雷達(dá)對地下管道探測技術(shù)發(fā)展的起步階段。在此之后，我國科研人員在探地雷達(dá)方面投入了大量的人力物力，我國的探測技術(shù)趨于成熟。在90年代中后葉，國防科技大學(xué)研究出一種對道路分辨率為4 mm左右、使系統(tǒng)分辨率達(dá)到3～5 cm的探地雷達(dá)[3]。此后，探地雷達(dá)在工程方面的應(yīng)用日漸增多，對管道探測的研究亦日趨成熟[4]。然而常規(guī)探地雷達(dá)大多數(shù)是單極化雷達(dá)，單極化雷達(dá)是采用固定極化方式的天線來發(fā)射和接收電磁波信號，只接收共極化信息[5]，故得到的是單極化數(shù)據(jù)。但單極化數(shù)據(jù)信息量遠(yuǎn)不足以滿足實際需求，對管道目標(biāo)準(zhǔn)確快速的識別造成了困難，于是全極化探地雷達(dá)在單極化雷達(dá)之后得到了快速的發(fā)展[6-12]。

全極化探地雷達(dá)是通過接收來自不同極化方向的電磁波，利用不同目標(biāo)體反射回來的電磁波信號中矢量特性的不同以及各矢量分量之間的關(guān)系，更好地獲得目標(biāo)體的位置及走向等各方面信息[12]。其優(yōu)點是能夠分析目標(biāo)體的極化特性，由全極化散射矩陣表述。通過該矩陣可以了解目標(biāo)體在每個測點上的極化特征，這些極化特征既包括共極化信息也包括交叉極化信息，從而能夠獲取目標(biāo)體較為全面的極化信息。極化分解技術(shù)主要就是通過分析目標(biāo)體的全極化信息，提取目標(biāo)體的極化屬性，有效地識別目標(biāo)體[7]。

1 全極化探地雷達(dá)地下管道探測原理

1.1 全極化探地雷達(dá)系統(tǒng)

全極化步進(jìn)頻率探地雷達(dá)系統(tǒng)(圖1)由矢量網(wǎng)絡(luò)分析儀、三維直角坐標(biāo)機器人、PC控制單元和全極化天線構(gòu)成。將發(fā)射天線、接收天線分別接在發(fā)射端口和接收端口，即可構(gòu)成完整的全極化探地雷達(dá)系統(tǒng)。該系統(tǒng)具有精確定位、測量、自動化存儲數(shù)據(jù)等功能，可獲得高質(zhì)量的全極化數(shù)據(jù)。

圖1 全極化探地雷達(dá)測量系統(tǒng)[7]Fig.1 Full polarimetric GPR measurement system[7]

在本文中用到的矢量網(wǎng)絡(luò)分析儀的型號是E5071C，頻帶寬度為：9 kHz～8.5 GHz，掃描點數(shù)可達(dá)1 601，輸出功率為-15～10 dBm，最大輸入功率為20 dBm。通過網(wǎng)絡(luò)分析儀，可以發(fā)射信號和接收信號，能夠進(jìn)行簡單的數(shù)據(jù)處理。在使用網(wǎng)絡(luò)分析儀之前首先要進(jìn)行校準(zhǔn)，校準(zhǔn)主要是去除電纜及儀器本身的干擾。我們一般用雙端口校準(zhǔn)，校準(zhǔn)件的型號是E5052C，包含4個校準(zhǔn)件，分別是短路校準(zhǔn)件、開路校準(zhǔn)件、負(fù)載校準(zhǔn)件和直通校準(zhǔn)件。

三維直角坐標(biāo)機器人用于空間位置精確定位，能在x,y,z3個方向上移動，測量精度可達(dá)0.1 mm。通過PC控制單元控制三維直角坐標(biāo)機器人的移動。

在該全極化探地雷達(dá)系統(tǒng)中，所用天線是Vivaldi天線。天線在探地雷達(dá)系統(tǒng)中至關(guān)重要，極化天線分為水平極化(H極化)天線和垂直極化(V極化)天線兩種。電磁波的傳播方向與極化方向一致時，我們稱之為H極化；波的傳播方向與極化方向垂直稱為V極化，由H極化和V極化可以構(gòu)成交叉極化(VH極化或者HV極化)[7]。

假設(shè)在散射和散射體全部確定的情況下，完全散射就能夠由散射矩陣描述。在單站后向散射情況下，發(fā)射天線和接收天線在同一位置時，測線中每個測點的散射矩陣S可表示為

(1)

式中,SHH，SHV，SVH和SVV是散射系數(shù)。SHH和SVV代表入射場和散射場的極化方式相同，稱之為“同極化”項；SVH和SHV代表入射場和散射場極化方式相互正交，稱之為“交叉極化”項。由于發(fā)射天線和接收天線是可交換的，且傳播的介質(zhì)也滿足互易性，因此后向散射矩陣也滿足互易性，即有SHV=SVH[14]。那么目標(biāo)矢量Ω可寫成

(2)

那么協(xié)方差矩陣C就可以寫成

C=Ω·Ω*T。

(3)

式中：*代表共軛變換；上角標(biāo)T表示轉(zhuǎn)置。

將目標(biāo)矢量Ω代入式(3)可得

(4)

計算協(xié)方差矩陣C的特征值和特征向量，可將公式(4)化為標(biāo)準(zhǔn)化形式[17-19]：

(5)

其中，

(6)

(7)

(8)

(9)

Z1--Z9 含義見表1。圖的二維特征分類空間[30]Fig.2 Two dimensional feature classification space of the

2 地下管道分類識別實驗

我們在實驗室均勻的沙槽中分別對單一管道、多個管道彼此之間相互影響以及管道受其他屬性目標(biāo)體影響3種不同類型的實驗進(jìn)行全極化測量。

2.1 單一管道極化屬性分析

將直徑為5 cm的金屬圓柱體(圖3)埋在2.5 m×2.5 m均勻的白沙沙槽中，目標(biāo)體大致埋在沙槽中間位置，埋深15 cm左右。用金屬圓柱體模擬實際環(huán)境中的管道，測線布設(shè)在目標(biāo)體正上方，用全極化探地雷達(dá)系統(tǒng)進(jìn)行測量，分別得到了HH，HV，VV 3種極化方式的測量剖面(圖4)。

圖3 兩個實驗?zāi)繕?biāo)體Fig.3 Two experimental objects

2.2 多個管道極化屬性分析

沙槽中均勻沙的相對介電常數(shù)約為4，沙中埋置4個相距很近的管道，其分布如圖6所示，從左到右管道的埋深依次為45,45,40,20 cm，橫向間隔為60，45，55 cm。中間兩個管道(2#、3#)是金屬管，兩側(cè)分別為PVC管，測量間隔為2 cm。用全極化探地雷達(dá)系統(tǒng)進(jìn)行全極化測量，分別得到了HH，HV，VV 3種極化方式的測量剖面(圖7)。

表1 Z1--Z9 代表的散射機制

a.HH，水平發(fā)射水平接收；b.HV，水平發(fā)射垂直接收；c.VV，垂直發(fā)射垂直接收。圖4 單一管道目標(biāo)體剖面圖Fig.4 Profiles of a single pipeline

圖5 單一管道目標(biāo)的的分布圖Fig.5 Distribution of of a single pipeline

圖6 管道的三維立體圖Fig.6 Three dimensional view of four pipes

a.HH，水平發(fā)射水平接收；b.HV，水平發(fā)射垂直接收；c.VV，垂直發(fā)射垂直接收。圖7 4根管道的剖面圖Fig.7 Profiles of the four pipelines

a.第一根管道；b.第二根管道；c.第三根管道；d.第四根管道。圖8 4根金屬管道的的分布圖Fig.8 Distribution of of the four metal pipelines

2.3 管道和其他非線性目標(biāo)體相互影響的極化屬性分析

將直徑為5 cm的金屬圓柱體和直徑為15 cm的實心金屬球(圖3)埋在各向同性的白沙中，2個目標(biāo)體大致在沙槽中間位置，目標(biāo)體相距60 cm左右。用全極化探地雷達(dá)系統(tǒng)進(jìn)行測量，測線依然布設(shè)在目標(biāo)體正上方，分別得到了HH，HV，VV 3種極化方式的測量剖面(圖9)。

為了驗證目標(biāo)體之間由于極化屬性相互影響致使數(shù)據(jù)點有一定的偏移想法的正確性，對此非線性目標(biāo)體金屬球進(jìn)行了極化屬性的分析，得到圖11。從圖11看出，球的數(shù)據(jù)點集中在Z9低熵表面散射區(qū)。但認(rèn)真觀察會發(fā)現(xiàn)，由于線性管道極化屬性的影響，位于Z9區(qū)的數(shù)據(jù)點卻向Z8線性目標(biāo)體區(qū)有一定的偏移，但主要的數(shù)據(jù)點依然集中在Z9區(qū)。

a.HH，水平發(fā)射水平接收；b.HV，水平發(fā)射垂直接收；c.VV，垂直發(fā)射垂直接收。圖9 管道和金屬球在全極化GPR系統(tǒng)測量下得到的剖面圖Fig.9 Profiles of pipeline and ball with full polarimetric GPR system

圖10 管道在非線性目標(biāo)體球影響下的的分布圖Fig.10 Distribution of of the pipe under the influence of nonlinear target ball

圖11 金屬球在管道的影響下的的分布圖Fig.11 Distribution of of the ball under the influence of pipeline

3 結(jié)論與展望

[1] 張鵬，董韜，馬彬，等. 基于探地雷達(dá)的地下管線管徑探測與判識方法[J]. 地下空間與工程學(xué)報，2015, 11(4): 1023-1032.

Zhang Peng, Dong Tao, Ma Bin, et al. Research on Interpreting the Information of Underground Pipeline’s Diameter Detected by GPR[J]. Chinese Journal of Underground Space and Engineering, 2015, 11(4): 1023-1032.

[2] 陳義群,肖柏勛. 論探地雷達(dá)現(xiàn)狀與發(fā)展[J]. 工程地球物理學(xué)報，2004, 2(2): 149-155.

Chen Yiqun, Xiao Boxun.On the Status Quo and Development of Ground Penetrating Radar[J]. Chinese Journal of Engineering Geophysics, 2004, 2(2): 149-155.

[3] 張文軍. 基于雷達(dá)的地下管道探測關(guān)鍵技術(shù)研究[D]. 長春：長春理工大學(xué)，2016.

Zhang Wenjun. Research on Key Technology of Underground Pipeline Detection Based on Radar[D]. Changchun: Changchun University of Science and Technology, 2016.

[4] 張文波，魏文博，景建恩，等. 利用探地雷達(dá)的極化特性檢測建筑物結(jié)構(gòu)[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版)，2008, 38(1): 156-160.

Zhang Wenbo, Wei Wenbo, Jing Jian’en, et al. Application of Ground Penetrating Radar Polarization in the Concrete Structure Detection[J]. Journal of Jilin University(Earth Science Edition), 2008, 38(1): 156-160.

[5] 汪夢奇. 分時極化雷達(dá)系統(tǒng)及其實現(xiàn)技術(shù)研究[D]. 北京：北京理工大學(xué)，2015.

Wang Mengqi. The Studies on Time-Division Polarimetric Radar System and Its Implementation Technology[D]. Beijing: Beijing Institute of Technology, 2015.

[6] 李麗麗，馮晅，鹿琪，等. 極化步頻探地雷達(dá)系統(tǒng)初步研究[J]. 吉林大學(xué)學(xué)報(地球科學(xué)版)，2008, 38(1): 150-152.

Li Lili, Feng Xuan, Lu Qi, et al. Preliminary Study on Polarization Step Frequency Radar Systems[J]. Journal of Jilin University(Earth Science Edition), 2008, 38(1): 150-152.

[7] 于月. 全極化探地雷達(dá)H-α特征分解技術(shù)研究[D]. 長春：吉林大學(xué)，2016.

Yu Yue. Development ofH-Alpha Feature Decomposition Technology For Full-Polarimetric Ground Penetrating Radar[D]. Changchun: Jilin University, 2016.

[8] 梁文婧. 全極化步進(jìn)頻率探地雷達(dá)硬件系統(tǒng)研究[D]. 長春：吉林大學(xué)，2012.

Liang Wenjing. Full-Polarimetric Step-Frequency GPR System Research[D]. Changchun: Jilin University, 2012.

[9] 李麗麗. 全極化探地雷達(dá)采集系統(tǒng)及校準(zhǔn)技術(shù)初探[D]. 長春：吉林大學(xué)，2010.

Li Lili. Fully-Polarimetric GPR Acquisition System and Calibration Technique Preliminary Study[D]. Changchun: Jilin University, 2010.

[10] Xuan F, Yue Y, Cai L, et al. Subsurface Polarimetric Migration Imaging for Full Polarimetric Ground-Penetrating Radar[J]. Geophysical Journal International, 2015, 202(2): 1324-1338.

[11] 鄒立龍. 全極化探地雷達(dá)正演模擬及極化校準(zhǔn)技術(shù)[D]. 長春：吉林大學(xué)，2012.

Zou Lilong. Forward Modeling and Calibration Technique of Full-Polarimetric GPR[D]. Changchun: Jilin University， 2012.

[12] 尤志鑫. 探地雷達(dá)地下淺層砂體和管道探測技術(shù)研究[D]. 長春：吉林大學(xué)，2016.

You Zhixin. Study of Underground Shallow Sandstone and Pipeline with Ground Penetrating Radar[D]. Changchun: Jilin University, 2016.

[13] Huynen J R. A Revisitation of the Phenomenological Approach with Applications to Radar Target Decomposition[R]. University of Illinois at Chicago: Department of Electrical Engineering and Computer Sciencesrr, 1982.

[14] Lee J S, Pottier E. 極化雷達(dá)成像基礎(chǔ)與應(yīng)用[M]. 北京：電子工業(yè)出版社，2013.

Lee J S, Pottier E. Fundamentals and Applications of Polarimetric Radar Imaging[M]. Beijing: Publishing House of Electronics Industry, 2013.

[15] Cloude S R, Pottier E. An Entropy Based Classifi-cation Scheme for Land Applications of Polarimetric SAR[J]. IEEE Trans Geosci Remote Sens, 1997, 35(1): 68-78.

[16]Zhao J G, Sato M. Radar Polarimetry Analysis Applied to Single-Hole Fully Polarimetric Borehole Radar[J]. IEEE Trans Geosci Remote Sens, 2006, 44(12): 3547-3554.

[17] Yilmaz ?, Doherty S M. Seismic Data Analysis: Processing, Inversion, and Interpretation of Seismic Data[M]// 2nd ed. Tulsa: Soc Exploration Geophys, 2001.

[18]Lauren F F, Pottier E, Lee J S. Unsupervised Classification of Multifrequency and Fully Polarimetric SAR Images Based on theH/A/Alpha-Wishart Classifier[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(11): 2332-2342.

[19]Lee J S, Grunes M R, Ainsworth T L, et al. Unsupervised Classification of Polarimetric SAR Images by Applying Target Decomposition and Complex Wishart Distribution[J]. IEEE Trans Geosci Remote Sensing, 1999, 37: 2249-2258.

[20] Lee J S, Schuler D L, Ainsworth T L. Polarimetric SAR Data Compensation for Terrain Azimuth Slope Variation[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(5): 2153-2163.

[21] Lee J S, Schuler D L, Ainsworth T L, et al. On the Estimation of Radar Polarization Orientation Shifts Induced by Terrain Slopes[J]. IEEE Transactions on Geoscience and Remote Sensing, 2002, 40(1): 30-41.

[22]Pottier E. Unsupervised Classification Scheme and Topography Derivation of POLSAR Data on theH/A/αPolarimetric Decomposition[C] //Theorem Proceedings of the 4th International Workshop on Radar Polarimetry. Nantes: Taylor & Francis Online, 1998: 535-548.

[23] Pottier E, Boerner W M, Schuler D L. Estimation of Terrain Surface Azimuthal /Range Slopes Using Polarimetric Decomposition of POLSAR Data[C]// Proceedings of IGARSS 1999. Hambourg: IEEE, 1999.

[24] Schuler D L, Lee J S, Grandi G D. Measurement of Topography Using Polarimetric SAR Images[J]. IEEE Transactions on Geoscience and Remote Sensing, 1996, 34(5): 1266-1277.

[25] Schuler D L, Lee J S, Ainsworth T L, et al. Terrain Slope Measurement Accuracy Using Polarimetric SAR data[C]// Proceedings of IGARSS 1999. Hambourg: IEEE, 1999.

[26] Schuler D L, Lee J S, Ainsworth T L, et al. Pola-rimetric DEM Generation from POLSAR Image Information[C]// Proceedings of URSI-XXVIth General Assembly. Toronto: University of Toronto, 1999.

[27] Schuler D L, Lee J S, Ainsworth T L, et al. Terrain Topography Measurement Using Multipass Polarimetric Synthetic Aperture Radar Data[J]. Radio Science, 2000, 35(3): 813-832.

[28] Fisher E, McMechan G A, Annan A P, et al. Exa-mples of Reverse-Time Migration of Single-Channel, Ground-Penetrating Radar Profiles[J]. Geophysics, 1992, 57(4): 577-586.

[29]Lee J S, Pottier E. Polarimetric Radar Imaging: From Basics to Applications[M]. Boca Raton: University of Rochester, 2009.

[30]Cassidy N J. Evaluating LNAPL Contamination Using GPR Signal Attenuation Analysis and Dielectric Property Measurements: Practical Implications for Hydrological Studies[J]. J. Contaminant Hydrol, 2007, 94(1/2): 49-75.