曹曉衛(wèi)李春江顏小飛吳一帆李志軍

摘要：黃河河道凍結(jié)結(jié)冰過程是黃河防凌汛研究的一個重點，進行彎道和橋墩周圍冰層厚度的探測對防治冰凌災(zāi)害具有重要的意義。探地雷達對冰厚探測具有便攜、高效、連續(xù)、快速、實時等優(yōu)勢。雷達天線的頻率決定雷達波的穿透深度和分辨率。從雷達圖像上可以清楚的識別空氣冰界面和冰水界面。通過現(xiàn)場打孔測量冰厚數(shù)據(jù)反算雷達波在冰層中的傳播速度進一步提高了測量結(jié)果的準(zhǔn)確性。測量結(jié)果清楚的顯示了彎道及橋墩處冰厚分布不均勻。橋墩北側(cè)冰層較厚，彎道主河道處冰層較厚。200 MHz天線的探地雷達可以穿透黃河冰層，雷達圖像可以清楚地顯示冰水界面。冰層厚度測量為分析黃河冰的凍結(jié)、融化過程以及冬春開河期冰塞、冰壩的治理提供了基礎(chǔ)數(shù)據(jù)。

關(guān)鍵詞：黃河；冰厚；探地雷達；彎道；橋墩

中圖分類號：TN959.3文獻標(biāo)志碼：A文章編號：

16721683（2016）06009105

Measuring ice thickness around the curve and piers in the Yellow River with ground penetrating radar

CAO Xiaowei1，LI Chunjiang2，YAN Xiaofei2，WU Yifan3，LI Zhijun1

（1.Stake Key Laboratory of Coastal and Offshore Engineering，Dalian University of Technology，Dalian116024，China；2.Yellow River Institute of Hydraulic Research，Yellow River Conservancy Commission，Zhengzhou 450000，China；3.School of Civil Engineering，Hefei University of Technology，Hefei 230009，China）

Abstract：The freezing process of yellow river is a key point in the study of Yellow River ice flood prevention，detection of the ice thickness around the curve and piers has important significance for prevention of the ice flood.The ground penetrating radar has the advantages of portable，highly efficient，continuous，fast and real time for ice thickness detection.The frequency of the radar antenna determines the penetration depth and resolution of the radar wave.From the radar images，the airwater interface and icewater interface could be clearly identified.The accuracy of the measurement results was further improved by using ice thickness data from the field drilling to calculate radar propagation velocity in the ice.The measurement results clearly showed that the distribution of ice thickness at the bend and pier was uneven.The ice on the north side of piers and the main channel of the curve was thicker.The ground penetrating radar with 200 Mhz antenna could penetrate the Yellow River ice，and the radar images could clearly show the icewater interface.This study could provide the basic data for the analysis of the yellow river ice freezing and melting process and governance of the ice jam in breakup period.

Key words：Yellow River；ice thickness；ground penetrating radar；curve；pier

[JP+1]黃河是我國的第二大河，河道全長5 400 km[1]。黃河河道凍結(jié)結(jié)冰過程是黃河防凌汛研究的一個重點，黃河河道結(jié)冰過程與氣溫，水流，河道邊界條件等因素有關(guān)。每到冬春黃河開河時期，融冰洪水中的冰塊受到冰蓋結(jié)構(gòu)的影響失去流動的動力，極易在冰蓋結(jié)構(gòu)上堆積、堵塞形成冰塞、冰壩結(jié)構(gòu)體，嚴重堵塞河道，造成冰凌洪水災(zāi)害[23]。彎道由于其本身地形的特殊性，浮冰及冰花易在彎道處堆積形成冰塞、冰壩。橋墩在冰期對河道中冰塞、冰壩的形成和演變會產(chǎn)生不可忽視的影響作用[4]。橋墩冰環(huán)境中受到的作用主要包括兩個方面：一方面是冰和冰蓋對橋墩的靜冰作用力[5]；另一方面是流冰對橋墩的撞擊和摩擦力[6]，這些都可能造成橋墩的破壞。彎道和橋墩周圍冰層厚度作為分析的基礎(chǔ)數(shù)據(jù)，對防治冰凌災(zāi)害和橋墩破壞具有重要的意義。

冰厚是較難監(jiān)測的物理指標(biāo)[7]。目前在大中尺度上的冰厚觀測上，國內(nèi)外已發(fā)展了許多自動化觀測技術(shù)，衛(wèi)星遙感[8]、雷達探測[9]、電磁感應(yīng)和激光測距技術(shù)的組合[10]、艦載聲吶[11]等。在定點冰厚觀測方面，系泊聲吶和超聲探測技術(shù)實現(xiàn)了連續(xù)監(jiān)測，但是測量精度只能達到±1 cm[1213]。雷瑞波等[14]研究了三種冰厚定點觀測方法：基于磁致伸縮原理的冰厚測量儀、電阻式冰厚測量儀、熱電阻絲冰厚測量裝置。在小尺度冰厚觀測方面，鉆孔測量雖然是測量冰厚最可靠的手段，但其效率較低，難以滿足一定范圍的冰厚觀測。探地雷達攜帶方便，可以高效、連續(xù)、快速、實時的對冰厚進行測量，同時已經(jīng)成功的應(yīng)用在海冰和淡水冰冰層厚度探測[1516]。

[JP+1]本文主要研究內(nèi)容是使用探地雷達探測黃河彎道和橋墩冰層厚度，探究探地雷達對黃河冰厚探測的可能性，通過現(xiàn)場打孔測量冰厚數(shù)據(jù)反算雷達波在冰層中的傳播速度來提高冰厚探測的準(zhǔn)確性，并根據(jù)現(xiàn)場水深數(shù)據(jù)對冰層厚度分布進行相關(guān)分析。

1測量區(qū)域和工作方法

本次測量分別于2016年3月1日和2016年3月3日對黃河巨合灘大橋（40°13.151′ N，111°11075′ E）和黃河什四份村彎道（40°17714′ N，111°2738′ E）進行冰層厚度的探測。彎道探測斷面分布見圖1，從進彎道到出彎道共選取了7個斷面，每個斷面間距不等的布置了9～13個測點進行冰層厚度和水深的測量。通過7個斷面冰層厚度的測量來對整個彎道冰層厚度分布進行分析。

巨合灘黃河公路大橋由西向東跨越黃河，大橋全長937 m，主橋跨徑為5×100 m[17]，共五個橋墩，從此次測量的現(xiàn)場圖（2（a））中可以看出，只有2個橋墩周圍凍結(jié)有冰層，定義為1號橋墩和2號橋墩，此次測量就圍繞這2個橋墩展開。為了研究橋墩的存在對于河道冰層厚度的影響，如圖（2（b））所示，本次測量沿水流方向設(shè)置2個縱斷面，垂直水流方向設(shè)置3個橫斷面。用探地雷達對這5個斷面的冰厚進行測量，來分析橋墩對冰層厚度分布的影響。

探測設(shè)備選用的是意大利IDS公司生產(chǎn)的RISK2型探地雷達系統(tǒng)，天線的頻率決定雷達波的穿透深度和分辨率。孫波等[18]在第一次中國北極科學(xué)考察期間進行有關(guān)雷達天線頻率與探測深度以及分辨率關(guān)系的研究，研究結(jié)果表明200 MHz的天線具有較高的分辨率，但穿透能力較弱。根據(jù)現(xiàn)場打孔測量冰厚情況，本次測量冰厚均在1 m左右，不需要較強穿透能力，為了得到較高的分辨率，故采用200 MHz的天線進行測量。

2雷達測冰厚的原理

探地雷達發(fā)射機向冰層下發(fā)射超高頻寬帶短脈沖電磁波，遇到不同介質(zhì)和界面時，部分脈沖波將反射回地面，被天線接收，然后進行數(shù)據(jù)處理[19]。在探測冰層厚度時，空氣冰界面和冰水界面處介質(zhì)發(fā)生改變，有反射波被天線接收，從圖3可以看出兩條顏色變化明顯的分界線，由此判斷出空氣冰界面和冰水界面的位置。

雷達波在冰中的傳播速度是決定雷達測量冰層厚度準(zhǔn)確性的重要因素，雷達波在冰中的傳播速度可以用下式計算：

[WTB1X]v[WT]=C/[KF（]ε[KF）][JY]（1）

式中：C為雷達波在空氣中的傳播速度（30 cm/ns），ε為物質(zhì)的相對介電常數(shù)，其中純冰的相對介電常數(shù)為317。天然冰相比純冰結(jié)構(gòu)成分復(fù)雜，雷達波的傳播速度受到冰晶體結(jié)構(gòu)，氣泡含量，泥沙含量等影響。李志軍等[20]對冰內(nèi)氣泡對雷達波在冰中傳播速度的影響進行了相關(guān)研究，研究結(jié)果表明：雷達波在粒狀冰中的傳播速度 1702 cm/ns，在柱狀冰中的傳播速度為1698 cm/ns。

本次測量雷達波在冰層中的傳播速度采用反算方法。如圖3所示，從雷達圖像中可以清楚的分辨出空氣冰界面和冰水界面，快速的提取出雷達波在冰層中的雙程傳播時間。在彎道處選取10個位置進行打孔測量冰厚，利用其中5個位置的實測冰厚，從對應(yīng)雷達波圖像提取雷達的雙程傳播時間，利用公式（2）推算得到雷達在冰中的傳播速度，如表1所示。通過計算得到雷達波的平均傳播速率是1609 cm/ns，由于此次測量處于融冰期，冰內(nèi)含水量較高，較盛冰期的雷達傳播速率慢。通過計算出來的平均速度，得到冰厚的計算公式（3）。

[WTB1X]v[WT]=2D/t[JY]（2）

式中：[WTB1X]v[WT]為雷達波在冰中傳播速率；D為人工測量冰厚；t為雷達波雙程傳播時間。

i=[WTB1X][AKv-D][WT]·t/2=7.86t[JY]（3）

式中：i為雷達測量冰厚（cm），[WTB1X][AKv-D]為雷達波[WT]在冰層中平均傳播速度（1609 cm/ns），t為雷達波在冰層中的雙程傳播時間（ns）。

利用其余5個位置冰厚對公式（3）進行驗證，見圖4。從圖中可以看出根據(jù)公式（3）計算得到雷達測量冰厚與實測冰厚相關(guān)系數(shù)為0994 1，后續(xù)冰厚計算均按公式（3）進行計算。

3結(jié)果分析

3.1彎道冰層厚度

[JP3]對圖（1（b））中虛線區(qū)域進行冰層厚度分析，得到彎道水深分布等值線圖，彎道冰厚等值線圖[21]，見圖5。從水深分布等值線圖（5（a））中可以清晰的分辨出主河道位置。

從圖（5（b））中可以明顯看出主河道冰層較主河道外側(cè)淺灘冰層厚，造成這種現(xiàn)象的原因是外側(cè)淺灘主要是自然凍結(jié)，而上游水流帶來的浮冰和冰花會在主河道處堆積，同時由于彎道的特性，浮冰和冰花會堆積在主河道外側(cè)，造成主河道外側(cè)冰層厚度較厚。從冰厚等值線圖（5（b））中可以看出，1號斷面中間部分冰厚較薄，主要是因為下游存在清溝造成1號斷面的冰厚變薄。比較4號斷面和7號斷面可以發(fā)現(xiàn)，上游進入彎道處的7號斷面冰厚較均勻，而處在彎道中間處4號斷面冰厚分布不均勻，特別是在主河道位置冰層厚度明顯增加。彎道的存在對彎道處的冰層厚度分布影響較大，對上下游的冰層厚度分布影響較小。

從圖（5（b））中可以看出彎道處冰層厚度分布不均勻，冰底形態(tài)高低起伏，極不平整，造成這種形態(tài)的主要原因如下。

（1）黃河主河道位置會發(fā)生改變，隨著主河道位置的變化，水流會對冰底不同位置產(chǎn)生沖刷，造成了冰層底部的高低不平；

（2）在黃河冰蓋凍結(jié)過程中，上游水流攜帶的浮冰和冰花會附著在冰蓋的不同位置，造成不同位置的冰厚差異。

（3）黃河本身泥沙含量較高，泥沙對黃河的凍結(jié)過程也會產(chǎn)生一定影響。

3.2橋墩周圍冰層厚度

使用探地雷達對橋墩縱橫斷面冰厚進行測量，結(jié)果見圖6、圖7。本次測量的南北橫斷面均在距離橋墩50 m遠的位置，從圖6可以看出橋墩南北兩側(cè)橫斷面冰厚分布無明顯差別，橋墩的存在對于南北兩側(cè)橫斷面冰厚分布并沒有產(chǎn)生較大的影響。從圖7可以看出，在縱斷面方向上，橋墩北側(cè)冰層厚度較厚，通過分析造成這種差異的原因如下。

（1） 從凍結(jié)過程分析，橋墩北側(cè)位于河流上游，水流攜帶的浮冰和冰花是冰蓋凍結(jié)的重要組成部分[22]，浮冰和冰花會在橋墩北側(cè)產(chǎn)生堆積，使橋墩北側(cè)的冰層凍結(jié)較厚。

（2） 從熱力學(xué)方面分析，橋墩北側(cè)受橋墩陰影的影響，受太陽光照較少，氣溫較橋墩南側(cè)低，凍結(jié)厚度較厚。

4結(jié)果和討論

結(jié)果表明，200 MHz天線的雷達波可以穿透黃河冰層，在選擇適合頻率的天線條件下，雷達圖像可以清楚地顯示冰水界面的位置。由于黃河冰是一種復(fù)雜介質(zhì)，又受到冰內(nèi)泥沙、氣泡以及冰溫等因素的影響，導(dǎo)致雷達波在冰層中的傳播速度出現(xiàn)一定的離散性。提高探地雷達探測精度最理想的方法是與冰鉆測厚數(shù)據(jù)相結(jié)合，得到雷達波在冰中傳播的實時速度，來提高雷達探測冰厚的準(zhǔn)確性。本次研究，對于未來衛(wèi)星遙感和航空雷達探測來講，有重要的參考價值。

利用雷達可以快速探測出冰層厚度，對于黃河冰情的監(jiān)測具有很好的應(yīng)用前景，對于分析黃河冰的凍結(jié)、融化過程以及冬春開河期冰塞、冰壩的治理具有重要的意義。未來希望可以加強對雷達波傳播速度的認識，來進一步提高探地雷達對黃河冰厚探測的準(zhǔn)確性。

參考文獻（References）：

[1]可素娟，呂光圻，任志遠.黃河巴彥高勒河段冰塞機理研究[J].水利學(xué)報，2000（7）：6669.（ KE Sujuan，LV Guangqi，REN Zhiyuan.Study on mechanism of ice jam formation in Bayangaole Section of Yellow River[J].Journal of Hydraulic Engineering，2000（7）：6669.（in Chinese）） DOI：10.3321/j.issn：05599350.2000.07.012

[2]劉東常，孟聞遠，謝巍，等.黃河冰塞冰壩防治研究新思路[J].華北水利水電學(xué)院學(xué)報，2011，32（2）：14.（LIU Dongchang，MENG Wenyuan，XIE Wei，et al.New ideas of prevention and control of ice jam and ice dam in Yellow River[J].Journal of North China Institute of Water Conservancy and Hydroelectric Power，2011，32（2）：14.（in Chinese）） DOI：10.3969/j.issn.10025634.2011.02.001

[3]王富強，王雷.近10年黃河寧蒙河段凌情特征分析[J].南水北調(diào)與水利科技，2014，12（4）：2124.（WANG Fuqiang，WANG Lei.Analysis of ice regime characteristics in the NingxiaInner Mongolia Reach of Yellow River in the recent ten years[J].South to North Water Transfers and Water Science & Technology，2014，12（4）：2124.（in Chinese））.DOI：10.13476/j.cnki.nsbdqk.2014.03.005

[4]Xiong F，Xu G.Numerical Investigation of River IceBridge Pier Interaction[C].Structures Congress.2009：110.

[5]陸欽年，段忠東，歐進萍，等.河冰對橋墩作用的冰荷載計算方法（Ⅱ）——冰壓力計算公式[J].自然災(zāi)害學(xué)報，2002，11（4）：112118.（LU Qinnian，DUAN Zhongdong，OU Jinping，et al.Calculational method of river ice loads on piers（Ⅱ）：The formula for ice pressure[J].Journal of Natural Disasters，2002，11（4）：112118.（in Chinese）） DOI：10.3969/j.issn.10044574.2002.04.019

[6]張麗敏.冰單軸壓縮強度與影響因素試驗研究[D].大連理工大學(xué)，2012.（ZHANG Limin.Experimental Study on uniaxial compressive strength of ice and influence factors[D].Dalian University of Technology，2012.（in Chinese））

[7]Laxon S，Peacock N，Smith D.High Interannual Variability of Sea Ice Thickness In The Arctic Region.[J].Nature，2003，425（6961）：94750.DOI：10.1038/nature02050

[8]Tamura T，Ohshima K I，Markus T，et al.Estimation of Thin Ice Thickness and Detection of Fast Ice from SSM/I Data in the Antarctic Ocean[J].Journal of Atmospheric & Oceanic Technology，2007，24（10）：431443.DOI：http：//dx.doi.org/10.1175/JTECH2113.1

[9]趙寶剛.渤海遼東灣冰區(qū)工程點雷達海冰監(jiān)測和預(yù)報技術(shù)研究[D].大連：大連海事大學(xué)，2008.（ZHAO Baogang.Study on ice monitoring and forecasting with radar of ice area projects point in Liaodong Bay Bohai Sea[D].Dalian：Dalian Maritime University，2008.（in Chinese））

[ZK）]

[10][ZK（#]Haas C.Evaluation of ShipBased ElectromagneticInductive Thickness Measurements of Summer SeaIce in the Bellingshausen And Amundsen Seas，Antarctica[J].Cold Regions Science & Technology，1998，27（1）：116.DOI：10.1016/S0165232X（97）000190

[11]D A Rothrock，Yu Y，G A Maykut.Thinning of the Arctic SeaIce Cover[J].Geophysical Research Letters，1999，26（23）：34693472.DOI：10.1029/1999GL010863

[12]崔華義，郭紀(jì)捷.相關(guān)技術(shù)在冰水界面測量中的應(yīng)用[J].海洋技術(shù)學(xué)報，2004，23（1）：3537.（CUI Huayi，GUO Jijie.The application of correlation in the measurement of icewater interface[J].Ocean Technology，2004，23（1）：3537.（in Chinese））.DOI：10.3969/j.issn.10032029.2004.01.008

[13]Richtermenge J A，Perovich D K，Elder B C，et al.Ice MassBalance Buoys：A Tool for Measuring and Attributing Changes in the Thickness of the Arctic SeaIce Cover[J].Annals of Glaciology，2006，44（1）：205210（6）.DOI：10.3189/172756406781811727

[14]雷瑞波，李志軍，秦建敏，等.定點冰厚觀測新技術(shù)研究[J].水科學(xué)進展，2009，20（2）：287292.（ LEI Ruibo，LI Zhijun，QIN Jianmin，et al.Investigation of new technologies for insitu ice thickness observation[J].Advances In Water Science，2009，20（2）：287292.（in Chinese））.DOI：10.3321/j.issn：10016791.2009.02.021

[15]王幫兵，孫波，田鋼，等.利用三維GPR探測北極夏季海冰厚度及下表面形態(tài)特征分析[J].物探與化探，2007，31（3）：189192.（WANG Bangbing，SUN Bo，TIAN Gang，et al.The application of the 3D GPR method to measuring the depth of arctic summer sea ice and analyzing its underside morphology[J].Geophysical&Geochemical Exploration，2007，31（3）：189192.（in Chinese）） DOI：10.3969/j.issn.10008918.2007.03.001

[16]王寧練，蒲健辰.祁連山八一冰川雷達測厚與冰儲量分析[J].冰川凍土，2009，31（3）：431435.（ WANG Ninglian，PU Jianchen.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）：431435.（in Chinese））

[17]張寧，侯慶軍，李衛(wèi)華，等.巨合灘黃河大橋混凝土配合比設(shè)計[J].內(nèi)蒙古公路與運輸，2011（2）：3032.doi：10.3969/j.issn.10050574.2011.02.012.（ZHANG Ning，HOU Qingjun，LI Weihua，at el.Design of the concrete mixture ratio of the Juhetan Yellow River Rigde[J].Highways & Transportation in Inner Mongolia，2011（2）：3032.（in Chinese）） DOI：10.3969/j.issn.10050574.2011.02.012

[18]SUN Bo，WEN Jiahong，HE Maobing，et al.Sea ice thickness measurement and its underside morphol-ogy analysis using radar penetration in the arctic ocean[J].Science in China，2003，46（11）：11511160.DOI：10.3969/j.issn.16747313.2003.11.007

[19]吳凱濱.淺地層探地雷達合成孔徑成像應(yīng)用研究[D].華南理工大學(xué)，2012.（WU Kaibin.Application research on synthetic aperture imaging of shallow subsurface ground penetrating radar[D].South China University of Technology，2012 （in Chinese））

[20]Li Z J，Jia Q，Zhang B S，et al.Influences of Gas Bubble and Ice Density on Ice Thickness Measurement By GPR[J].Applied Geophysics，2010，7（7）：10513.DOI：10.1007/s1177001002344

[21]張弛.基于Surfer軟件的地形圖快速繪制方法[J].水利技術(shù)監(jiān)督，2009，17（3）：4143.（ZHANG Chi.Fast rendering method of topographic map based on Surfer Software[J].[JP2]Technical Supervision in Water Resources，2009，17（3）：4143.（in Chinese）） DOI：10.3969/j.issn.10081305.2009.03.017

[22]張紅彪，韓紅衛(wèi)，張寶森，等.黃河小型清溝封凍過程現(xiàn)場模擬實驗與實測數(shù)據(jù)分析[J].數(shù)學(xué)的實踐與認識，2014（23）：134139.（ZHANG Hongbiao，HAN Hongwei，ZHANG Baosen，et al.Site simulation test and data analysis of freezing process of open water surface in The Yellow River[J].Mathematics in Practice and Theory，2014（23）：134139.（in Chinese））