黃河河道冰層雷達波特征圖譜的現場實驗研究

張寶森+張防修+劉滋洋+韓紅衛+李志軍

摘要：2013年-2014年冬季，利用RIS K2型探地雷達在內蒙古頭道拐水文站黃河河道斷面開展冰層厚度探測試驗。試驗使用不同頻率天線探測冰厚，并與鉆孔實測冰厚對比。結果顯示200 MHz頻率天線可探測到冰下更深部位的狀況，并能夠識別出冰花層。利用雷達探測圖譜中電磁波傳播時間和實測冰厚得到盛冰期雷達電磁波在冰內的實際傳播速度為16.3 cm/ns；而融化期雷達電磁波在冰內的實際傳播速度小于15.0 cm/ns。雷達探測結果還顯示垂直于河道斷面方向的冰層厚度分布不均勻。主流區是以熱力學生長的顆粒冰和柱狀冰為主，最大冰厚約為60 cm；而非主流區是以冰花冰為主，最大冰厚約為70 cm。另外，雷達圖譜也可以確定冰層裂縫的具體位置、走向。

關鍵詞：冰層；探地雷達；雷達圖譜；實驗分析；黃河

中圖分類號：TV698.2；TN959.3文獻標志碼：A文章編號：1672-1683（2017）01-0121-05

Abstract：Radar measurement experiments were performed with a RIS K2 at the Toudaoguai Hydrological Station section in the winter of 2013-2014.The experiments detected ice thickness using radar with different frequency antenna and compared the detected thickness with measured thickness by drilling.The results showed the 200-MHz antenna can detect deeper depth below river ice and identify the frazil ice layer.The radar transmission speed in ice was determined based on the propagation time of radar electromagnetic wave and measured ice thickness.The radar transmission speed was about 16.3 cm/ns in the freezing period and less than 15.0 cm/ns in melting period.The GPR images showed that the thickness of the ice layer was not uniform along the section perpendicular to the course of river.In the main flow area between Pile No.400 and No.520，the maximum ice thickness was about 60 cm and the ice was mainly made up of granular ice and columnar ice created by thermodynamic growth.In the non-main flow area between Pile No.520 and No.740，the maximum ice thickness was about 70 cm and the ice was mainly made up of frazil ice.In addition，The GPR images can be used to locate the position and direction of the ice cracks in ice layer.

Key words： ice layer；ground-penetrating radar （GPR）；radar images；experimental analysis；Yellow River

我國東北、西北及青藏高原地區的冬季均存在河流湖泊結冰現象。河湖冰的存在對其包圍的結構物構成威脅，如河中流冰對橋墩的撞擊[1]，融冰期湖冰熱膨脹應力對護坡的破壞[2]；另外冰作為一種天然材料可以利用，如冰上汽車通行和冰上飛機起降[3-4]。也可以作為冰燈藝術的材料[5]。無論是冰的防治，還是冰的利用，冰厚度均是一個關鍵參數。因此，冰厚度在冰凌成因分析及預報[8]、冰工程[9]中是至關重要的因子之一。它對南水北調工程中的冰蓋膨脹力[6]和輸水能力[7]有直接貢獻。目前有多種方法獲得冰厚度，其中最直接的方式是原位鉆孔測量和熱電阻絲接觸式測量[10]。這兩種方法雖然簡單，但是效率偏低，很難實現大范圍冰厚度調查；冰下超聲測距儀（測量距離為150 cm，設計精度為1 mm，數據分辨率為0.01 mm）可以實現定點高精度連續測量冰厚，在北極海冰[11]、國內大慶水庫冰[10]、青藏高原北麓河熱融湖塘冰[12]厚度測量中得到應用。另外用磁致位移傳感器，在冰層的底面和表面均結合機械運動磁環，獲得過高精度定點接觸式測量冰厚[13]。而海冰厚度測量注重觀測范圍大，它利用激光高度計衛星（ICESat）和合成孔徑雷達（SAR）對海冰大區域觀測。其中ICESat原理為通過測量海水與海冰之間的高度差（干舷），然后根據冰在水中的浮力原理計算海冰厚度[14-17]，而這種方法無法消除冰上積雪的影響。當表面積雪越厚，觀測所得海冰厚度的誤差越大。SAR使用主動微波傳感器，與被動微波傳感器相比其精度有所提高，研究表明SAR 具有探測海冰厚度的潛力[18-20]。

目前黃河冰厚測量方式主要是人工測量，它只有在河流封凍且工作人員在冰上能夠自由活動時才能實現測量。另外，黃河部分位置的寬度均是百米級，給同量級的衛星遙感帶來挑戰。另外河冰下存在的冰花，發生的冰塞影響到定點測量技術的推廣。因而尋求非接觸式高效測量是必然趨勢，其中無人機懸掛探地雷達探測，能夠解決不穩定冰層厚度探測。對于初冰期、融冰期和冰塞冰壩期的冰厚度探測具有快速安全特點。本文使用的探地雷達技術是原位冰面測量，它具有無損、連續、快速、高效、高精度、實時成像探測等特點，在淡水冰[21]、海冰[22]、冰川冰[23]均有應用。是發展非接觸探測的雛形。通過定期對比同一河道斷面冰層探測圖像，可準確確定冰層厚度和結構變化[24]，能夠及時發現異常，為解譯冰層厚度變化和結構變化提供科學判據。

本文以黃河頭道拐為依托，于2013年-2014年度冬季對河道斷面開展反復雷達探測，獲得了一批雷達測試冰層厚度數據。同時取樣測試冰晶體結構和密度等參數。在此介紹冰層雷達圖譜，以期為淡水冰雷達測厚技術提供參考和經驗積累。

1 雷達探測河道冰層方法

由于冰、純水、底層沉積物介電常數分別為3～4，81和5～40，根據電磁波傳播理論，電磁波在材料內的傳播速度只與介電常數有關。因此探地雷達的電磁脈沖波在三種介質中的傳播速度不同，接收到的信號可反映出不同介質界面間的位置。探地雷達測量冰層厚度時，假設雷達波到達冰-水界面和返回接收端的時間相同，因此計算冰厚時只用到雙程走時的一半。

式中：H為冰層厚度（cm）；ΔT為雙程用時（ns）；ε為冰的介電常數；V為雷達波在冰內傳遞速度，其計算公式如下：

式中：c為電磁波在真空中的傳播速度（30 cm/ns）。

當實際冰的介電常數不確定時，可以根據實際冰厚反演得到冰的等效介電常數。本試驗所用的雷達是意大利IDS公司生產的RIS K2型探地雷達。它包括主機、天線和處理軟件三部分。

2 雷達探測黃河河道冰層結果

以黃河頭道拐水文站水文斷面附近為試驗場。將雷達天線緊貼冰面，以步行方式拖動雷達。邊行走，邊顯示雷達電磁波圖像，邊記錄。每條測線記錄均獨立編寫文件號，所有正式探測過程中如果遇隨機情況且在雷達波圖譜上出現異常時，則在現場重復觀測。測量模式采用實時測量和距離測量模式。設置完畢后還需在冰面上拖動雷達對雷達進行增益處理，之后開展斷面探測。在斷面探測中采用垂直于河道和平行于河道兩種方式，以便后期建立河道內冰層內部組構和厚度的三維分布狀態。

由于冰的介電常數與冰的晶體組構和冰溫有關。所以在每次探測冰厚時需要鉆孔實測冰厚。根據實測冰厚反演探測時冰層的有效介電常數，然后再依據該參數修正當時整個斷面的冰厚。另外，在典型河道位置，鉆取冰芯并分析冰的晶體組構。

影響雷達探測冰厚精度的因素除了雷達波在冰內的傳播速度外，還與探測天線的頻率有關[9]。前者由冰層的物理性質決定，后者由設備硬件決定。在2014年1月4日探測黃河河道冰層中，分別使用頻率200、400、900和1 500 MHz天線對頭道拐黃河垂直于河道和平行于河道斷面的河道冰層進行探測，表1給出了各頻率天線測量冰厚與鉆孔實測冰厚對比，可以看出各頻率天線都能夠滿足冰層探測要求，但使用200 MHz頻率天線可探測到更深部位的冰-水狀況，能夠識別出冰花厚度。因此，選用200 MHz頻率天線對黃河冬季冰層進行探測。

圖1是頭道拐水文站黃河水文斷面不同日期的冰厚斷面分布。該冰厚斷面在2014年2月中旬達到最大值。之后，冰層開始逐步消融，冰厚逐漸變薄。在同一時區，近岸冰較厚，非主流區次之，主流區冰厚較薄。在整個冰期斷面樁號400～520為主流區，斷面樁號520～740為非主流區，主流區的冰厚始終比非主流區的冰厚小10～20 cm。在主流區流速較快，下層冰花少，主要以熱力學生長的平整冰為主，最大冰厚約為60 cm；而在樁號520～740為非主流區流速慢，冰花容易在冰層下堆積，水內冰花容易凍結成為冰花冰，最大冰厚約為70 cm。冰晶體分析也證明樁號520～740的非主流區冰晶體類型為冰花冰，而樁號400～520的主流區冰晶體類型為粒狀冰和柱狀冰。

圖2給出了雷達電磁波在冰內傳播速度隨冰期變化規律。由圖可知，在2014年1月19日盛冰期時雷達電磁波在冰內傳播速度約為16.3 cm/ns；在2014年2月24日以后的融冰期時降低到15.0 cm/ns以下。雷達電磁波在冰內傳播速度減慢與冰內晶體組構和冰內溫度升高后的融水增加有關。不同時期，由于冰晶體組構不同，雷達波在冰內傳播速度也不同；特別是在進入融冰期時，冰溫度升高，晶體之間的融水增加，冰晶體的連接變得松軟，甚至有冰水混合存在。冰融化引起冰內含水量增加同時，也引起雷達電磁波在水內傳播速度降低；也可以認為是引起等效介電常數增加所產生的結果。

3 黃河河道現場典型冰層雷達圖譜特征描述

3.1 平整冰層

在平整冰層位置處，雷達圖像的波形特征表現為：能量團分布均勻或僅在局部存在強反射細亮條紋；電磁波能量衰減緩慢，探測距離遠且規律性較強；一般形成低幅反射波組，波形均勻，無雜亂反射，自動增益梯度相對較小。圖3給出了黃河頭道拐水文站平整冰層雷達圖譜。

3.2 含破碎帶和裂隙帶冰層

河道水位發生變化時，冰層將發生變形，但受到河道邊界約束，冰層內部會發生擠壓或斷裂，或在冰層內產生尺寸不同的裂縫。另外，河冰在動力作用下，冰層內也會產生破碎冰體。特別是在裂縫處會夾雜泥沙或水。這些破碎冰體、泥沙或水的介電常數同冰的介電常數差異大，導致破碎帶和裂縫帶的雷達圖譜和波形特征與平整冰層的有較大差異。一般在破碎帶和裂縫帶處雷達信號反射強烈，反射面的振幅顯著加強且變化劇烈；還會表現出能量團分布不均勻，在破碎帶和裂縫帶內常產生繞射和散射，同相軸出現錯斷；另外在深部甚至出現雷達圖像模糊不清現象。通常在裂縫帶會出現雙曲線強反射波，在破碎帶處會出現多條雙曲線強反射波。反射波同相軸的連線為破碎帶或裂縫帶的位置。圖 4給出含有破碎帶和裂縫帶的冰層雷達特征圖譜。

3.3 冰花集聚區的冰層

黃河中上游在冬季除了冰封外，其水流中還有冰花。由于冰花密度低于水的密度，冰花會上浮聚集于冰層下方，形成冰水混合狀態。對應的等效介電常數有別于液態水，所以探地雷達電磁波在冰花聚集區內產生繞射和散射；在穿過冰花聚集區時會產生一定規律的多次強反射。由冰花聚集區連續分布時，雷達電磁波穿過時的反射波連續性較好，波形相對均一。由于雷達探測冰層和冰花時雷達電磁波從高阻抗介質到低阻抗介質，因而反射電磁波與入射電磁波相位相反，以此可以判斷冰層和冰花的界面位置。圖 5給出含有冰花層的雷達特征圖譜。

4 結論

（1）雷達波在冰內傳播速度的變化是影響雷達探測冰厚準確性的決定因素。不同時期，由于冰內部組構的差異，凍結期的雷達波在冰內傳播速度為16.3 cm/ns，融冰期的對應速度小于15.0 cm/ns。

（2）雷達探測黃河頭道拐水文站不同時期的冰厚結果顯示，垂直于河道斷面的冰厚不均勻。在樁號400～520的主流區，以熱力學生長的顆粒冰和柱狀冰為主，最大冰厚約為60 cm；而在樁號520～740的非主流區，以冰花凍結形成的冰花冰為主，最大冰厚約為70 cm。

（3）根據雷達探測冰層所得雷達圖譜，可以定性判斷冰層的結構，冰層的連續性；冰層內部裂縫的具體位置和走向。

參考文獻（References）：

[1] 張宿峰.流冰與橋墩的相互作用[D].哈爾濱：東北林業大學，2014.（ZHANG Su-feng.Ice and the interaction of bridge piers[D].Harbin：Northeast Forestry University，2014.（in Chinese））

[2] 姜連杰，盧永超，韓紅衛，等.紅旗泡水庫冰層變形及護坡破壞現象成因分析[J].水利水電技術，2011，42（7）：85-88.（JIANG Lian-jie，LU Yong-chao，HAN Hong-wei，et al.Causal analysis on ice layer deformation and slope protection damage of Hongqipao Reservoir[J].Water Resources and Hydropower Engineering，2011，42（7）：85-88.（in Chinese）） DOI：10.13928/j.cnki.wrahe.2011.07.024

[3] Gold L W.Use of Ice Covers for Transportation[J].National Research Council of Canada，1971，8（2）：170-181.DOI：10.1139/t71-018

[4] 徐伯孟，于生清，盧興良.嫩江冰上交通[J].冰川凍土，1992，14（4）：374-380.（XU Bo-meng，YU Sheng-qing，LU Xing-laing.Traffic and transport on ice cover in Nenjiang River[J].Journal of Glaciology and Geocryology，1992，14（4）：374-380.（in Chinese））

[5] Dewar K，Meyer D，and Li W.Harbin，Lanterns of Ice，Sculptures of Snow[J].Tourism Management，2001，22：523-532.DOI：10.1016/S0261-5177（01）00007-3

[6] 周庭正，趙新，黃焱.南水北調中線工程渡槽結構冰蓋溫度膨脹力研究[J].南水北調與水利科技，2012，10（3）：7-11.（ZHOU Ting-zheng，ZHAO Xin，HUANG Yan.Ice thermal expansive pressure acting on aqueduct structure of Middle Route Scheme of SNWDP[J].South-to-North Water Transfers and Water Science and Technology，2012，10（3）：7-11.（in Chinese）） DOI：10.3724/SP.J.1201.2012.03007

[7] 趙新，張雅卓，王旭.加厚冰蓋滲流對渠道輸水能力的影響[J].南水北調與水利科技，2012，10（4）：7-9.（ZHAO Xin，ZHANG Ya-zhuo，WANG Xu.Impact analysis of seepage flow in thickened ice cover on the water diversion capacity in canal[J].South-to-North Water Transfers and Water Science and Technology，2012，10（3）：7-9.（in Chinese）） DOI：10.3724/SP.J.1201.2012.04007

[8] 霍世青，溫麗葉，范旻昊，等.2014-2015年度黃河寧蒙河段凌情及氣象成因[J].人民黃河，2016，38（2）：16-18，23.（HUO Shi-qing，WEN Li-ye，FAN Min-hao，et al.Analysis of ice flood characteristics and meteorological causes at Ningxia-Inner Mongolia reach of Yellow River in year 2014-2015[J].Yellow River，2016，38（2）：16-18，23.（in Chinese）） DOI：10.3969/j.issn.1000-1379.2016.02.005

[9] 劉詩華，侯樹強.寒區核電廠取水口防冰設施設計研究[J].人民黃河，2015，37（5）：112-115.（LIU Shi-hua，HOU Shu-qiang.Research on ice-preventing facility schemes of water-intake in a nuclear power plant in cold regions[J].Yellow River，2015，37（5）：112-115.（in Chinese）） DOI：10.3969/j.issn.1000-1379.2015.05.029

[10] 李志軍，楊宇，彭旭明，等.黑龍江紅旗泡水庫冰生長過程現場觀測數據的剖析[J].西安理工大學學報，2009，25（3）：270-27.LI Zhi-jun，YANG Yu，PENG Xu-ming，et al.The analysis of the field observation data of fresh ice growing process in Hongqipao Reservoir of Heilongjiang[J].Journal of Xi′an University of Technology，2009，25（3）：270-27.（in Chinese））

[11] Lei Rui-bo，Li Zhi-jun，Cheng Bin，et al.Investigation of the Thermodynamic Processes of a Floe-lead System in the Central Arctic during later Summer[J].Advances in Polar Science，2011，22（1）：42-48.DOI：10.3724/SP.J.1085.2011.00010

[12] Huang Wen-feng，Li Zhi-jun，Han Hong-wei，et al.Structural Analysis of Thermokarst Lake Ice in Beiluhe Basin，Qinghai-Tibet Plateau[J].Cold Region Science and Technology，2012，72：33-42.DOI：10.1016/j.coldregions.2011.11.005

[13] Lei Rui-bo，Li Zhi-jun，Cheng Yan-feng，et al.A New Apparatus for Monitoring Sea Ice Thickness Based on the Magnetostrictive-delay-line Principle[J].Journal of Atmospheric and Oceanic Technology，2009，26（4）：818-827.DOI：10.1175/2008JTECHO613.1

[14] Zwally H J，Yi D，Kwok R，et al.ICESat Measurements of Sea Ice Freeboard and Estimates of Sea Ice Thickness in the Weddell Sea[J].Journal of Geophysical Research，2008，113（C2）：C02S15.DOI：10.1029/2007JC004284，2008

[15] Kwok R，Cunningham G F.ICESat over Arctic Sea Ice：Estimation of Snow Depth and Ice Thickness[J].Journal of Geophysical Research，2008，113（C8）：C08010.DOI：10.1029/2008JC004753，2008

[16] Farrell S L，Laxon S W，McAdoo D C，et al.Five Years of Arctic Sea Ice Freeboard Measurements from the Ice，Cloud and Land Elevation Satellite[J].Journal of Geophysical Research，2009，114（C4）：C04008.DOI：10.1029/2008JC005074，2009

[17] Price D，Rack W，Haas C，et al.Sea Ice Freeboard in McMurdo Sound，Antarctica，Derived by Surface-validated ICESat Laser Altimeter Data[J].Journal of Geophysical Research，2013，118（3634-3650.DOI：10.1002/jgrc.20266，2013

[18] Nakamura K，Wakabayashi H，Naoki K，et al.Observation of Sea-ice Thickness in the Sea of Okhotsk by Using Dual-frequency and Fully Polarimetric Airborne SAR （Pi-SAR） Data[J].Geoscience and Remote Sensing，IEEE Transactions，2005，43（11）：2460-2469.DOI：10.1109/TGRS.2005.853928

[19] Nakamura K，Wakabayashi H，Uto S，et al.Observation of Sea-ice Thickness Using ENVISAT Data from Lützow-Holm Bay，East Antarctica[J].Geoscience and Remote Sensing Letters，IEEE，2009，6（2）：277-281.DOI：10.1109/LGRS.2008.2011061

[20] Kwok R，Nghiem S V，Yueh S H，et al.Retrieval of Thin Ice Thickness from Multifrequency Polarimetric SAR Data[J].Remote Sensing of Environment，1995，51（3）：361-374.DOI：10.1016/0034-4257（94）00017-H

[21] Li Zhi-jun，Jia Qing，Zhang Bao-sen，et al.Influences of Gas Bubble and Ice Density on Ice Thickness Measurement by GPR[J].Applied Geophysics，2010，7（2）：105-113.DOI：10.1007/s11770-010-0234-4

[22] 王幫兵，孫波，田鋼，等.利用三維GPR探測北極夏季海冰厚度及下表面形態特征分析[J].物探與化探，2007，31（3）：189-192.（WANG Bang-bing，SUN Bo，TIAN Gang，et al.The application of the 3D GPR method to measuring the depth of arctic summer sea ice and analysis its underside morphology[J].Geophysical and Geochemical Exploration，2007，31（3）：189-192.（in Chinese））

[23] 王寧練，蒲健辰.祁連山八一冰川雷達測厚與冰儲量分析[J].冰川凍土，2009，31（3）：431-435.（WANG Ning-lian，PU Jian-chen.Ice thickness，sounded by ground penetrating radar，on the Bayi Glacier in the Qilian Mountains，China[J].Journal of Glaciology and Geocryology，2009，31（3）：431-435.（in Chinese））

[24] 戴長雷，于成剛，廖厚初，等.寒區水科學及國際河流研究系列叢書3-冰情監測與預報[M].北京：中國水利水電出版社，2010.（DAI Chang-lei，YU Cheng-gang，LIAO Hou-chu，et al.Research on water science in frigid zone and international rivers 3-survey and forecast of river ice[M].Beijing：China Water and Power Press，2010.（in Chinese））