， ，，，

(1.山東省第七地質(zhì)礦產(chǎn)勘查院，山東 臨沂 276000； 2.貴州省質(zhì)安交通工程監(jiān)控檢測(cè)中心有限責(zé)任公司，貴陽(yáng) 550081； 3.貴州大學(xué) 資源與環(huán)境工程學(xué)院，貴陽(yáng) 550025)

0 引 言

探地雷達(dá)(GPR)是研究高山積雪厚度和雪水當(dāng)量(SWE)的有效工具，相比于其他技術(shù)，GPR能在有限的時(shí)間里探查更廣大的區(qū)域[1-3]。與地震數(shù)據(jù)處理類似，要將雷達(dá)波的雙程旅行時(shí)間轉(zhuǎn)化為深度，必須要知道雷達(dá)波在介質(zhì)中傳播的速度，因此波速的推算就顯得尤為重要。在實(shí)際應(yīng)用中，可以通過(guò)設(shè)計(jì)雙天線探地雷達(dá)的觀測(cè)方式，并運(yùn)用相應(yīng)的地球物理數(shù)據(jù)的處理方法，達(dá)到測(cè)量介質(zhì)中電磁波波速的目的，如共中心點(diǎn)法(CMP)[4]、共深度點(diǎn)法(CDP)[1]、偏移分析法[5]、寬角反射/折射分析法(WARR)[6]以及反射波振幅反演等[7]。另外，還有其他電磁波方法和技術(shù)，如時(shí)域反射法(TDR)和雪特性分析儀(Snow Fork)，它們可以測(cè)量探針附近一小塊積雪區(qū)域的電磁波速度[8-10]，一般情況下，可以將其作為整個(gè)調(diào)查區(qū)域的電磁波波速。

除此之外，還有根據(jù)雪的密度求取波速的方法，即通過(guò)模型方程，將雪的密度與雪的介電常數(shù)聯(lián)系起來(lái)[1，3，6，8-9，11]。這種方法需要明確區(qū)分雪的干、濕狀態(tài)，以便準(zhǔn)確地選擇相應(yīng)的方程。

在最簡(jiǎn)單的干雪狀態(tài)下(僅僅是冰與空氣的混合)，表征雪的介電常數(shù)和密度關(guān)系的方程也有許多種，因此選擇合適的方程就有一定難度。倘若雪中還摻雜了部分水，這種三相混合的情況將會(huì)更加復(fù)雜，誤差也會(huì)大大增加。此外，大多數(shù)模型方程都沒(méi)有給出明確的不確定度，這會(huì)影響到對(duì)計(jì)算結(jié)果的可信度分析。

本文所研究的重點(diǎn)是在干雪狀態(tài)下，對(duì)表征雪的介電常數(shù)和密度關(guān)系的模型方程進(jìn)行比較，選出最合適的模型方程，并加以驗(yàn)證。

1 密度-介電常數(shù)的關(guān)系

幾十年來(lái)，許多學(xué)者使用時(shí)域反射法(TDR)或頻率域分析法測(cè)算干雪/濕雪的介電常數(shù)，根據(jù)雪的密度(ρ)和其介電常數(shù)(ε)的測(cè)量結(jié)果，提煉出兩者之間的關(guān)系。本文將一些較為經(jīng)典的ε(ρ)關(guān)系方程歸納見(jiàn)表1，表1中部分方程適用于濕雪狀態(tài)，可以將其模型中的含水率定為零，使其適用于干雪狀態(tài)[12]。

表1 雪的密度-介電常數(shù)關(guān)系

表1中雪的密度單位是g/cm3。根據(jù)表1，關(guān)系方程大致可以分為以下3類：

①線性關(guān)系：

ε=a1+a2ρ

②二次非線性關(guān)系：

ε=a1+a2ρ+a3ρ2或ε=(1+a4ρ)2

③混合規(guī)則：如文獻(xiàn)[13]、文獻(xiàn)[14]

相應(yīng)文章的作者參考自己測(cè)量得到的數(shù)據(jù)及其他文獻(xiàn)的數(shù)據(jù)，得到這些線性和二次非線性關(guān)系方程。一些作者通過(guò)更新測(cè)量數(shù)據(jù)，簡(jiǎn)單地修正了方程，因此表1中有些方程的常數(shù)項(xiàng)有許多相似之處。

大多數(shù)的文章作者都沒(méi)有給出相應(yīng)方程或者參數(shù)的不確定度，這限制了方程的實(shí)際應(yīng)用，并使得方程之間的對(duì)比變得困難[15]。只有個(gè)別作者給出了不確定度，具體如下：

文獻(xiàn)[17]估計(jì)ε的不確定度為±2%。

文獻(xiàn)[18]和文獻(xiàn)[19]的擬合曲線的標(biāo)準(zhǔn)差為0.031。

文獻(xiàn)[20]的擬合曲線的標(biāo)準(zhǔn)差為0.006。

缺少明確的不確定度，選擇哪種模型方程計(jì)算雪的介電常數(shù)就沒(méi)有對(duì)比依據(jù)。

表1中的部分文章作者給出的不確定度，已在圖1(b)中表示。在模型[13]、模型[14]、模型[21]和模型[22]中，冰的介電常數(shù)定為3.19±0.04[12]，冰的密度定為920±10 kg/m3[23]。由于模型[16]和模型[24]的擬合曲線基本一致，所以進(jìn)行合并處理。

圖1為雪的密度-介電常數(shù)關(guān)系擬合曲線。圖1(a)為表1方程的擬合曲線，圖1(b)是對(duì)圖1(a)矩形框中的圖形進(jìn)行放大，以更好地展示模型[16-20]的不確定度。

從圖1(b)中可以看出，Robin等的模型(黑色線，模型編號(hào)[16])可以作為最合理的選擇。因?yàn)樵撃Ｐ偷牟淮_定度區(qū)間幾乎囊括了其他模型的范圍，只有模型[32]在較高的密度區(qū)間上，稍微偏離了模型[16]的不確定度區(qū)間。

圖1 雪的密度-介電常數(shù)關(guān)系擬合曲線

因此，可以選擇模型[16]的方程來(lái)求取雪中的電磁波波速。

2 基于干雪的密度求取波速

使用擬合模型[16]，可以通過(guò)測(cè)得的密度值得到雪的介電常數(shù)及其不確定度，并得到雪中的電磁波波速：

(1)

式(1)中的c是真空中的光速，根據(jù)不確定度的傳播率[15]，波速υ的不確定度uυ取決于以下兩個(gè)因素：

①根據(jù)模型[16]得到的介電常數(shù)的不確定度：

②密度測(cè)量的不確定度：uρ

因此，可以得到波速υ的不確定度uυ為：

(2)

干雪的密度-電磁波波速關(guān)系曲線見(jiàn)圖2，數(shù)值的計(jì)算是基于模型[16]。在圖2中，左邊軸表征在干雪的密度范圍內(nèi)，同時(shí)在模型不確定度的影響下，其電磁波波速的最大值和最小值的范圍(以虛線表示)；右邊軸表征波速υ的相對(duì)不確定度，可由式(2)得到。

需要注意的是，在圖2中將密度測(cè)量的不確定度(uρ)定為零，因此圖2中的波速υ的相對(duì)不確定度僅與擬合模型中的介電常數(shù)的不確定度(uε)有關(guān)。若代入實(shí)驗(yàn)數(shù)據(jù)，就需要考慮雪密度測(cè)量的不確定度，所以實(shí)測(cè)條件下，uυ將會(huì)比圖2中的值更大。

圖2 干雪的密度-電磁波波速關(guān)系曲線

另外，由模型[16]得到的波速υ的相對(duì)不確定度，其量級(jí)與其他測(cè)量方法得到的波速υ的相對(duì)不確定度一致，如時(shí)域反射法(TDR)[8-9]、雪特性分析儀(Snow Fork)[10]、共中心點(diǎn)法(CMP)[4]、偏移分析法[5]等。

3 實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證

本文所使用的數(shù)據(jù)來(lái)源于Paolo等測(cè)量的數(shù)據(jù)[9]，其數(shù)據(jù)采集自意大利阿爾卑斯山東部的Solda冰川附近，其3個(gè)采樣位置(雪坑)中的雪均為干雪[9]。時(shí)域反射法(TDR)測(cè)量使用的儀器為美國(guó)泰克 Tektronix 1502C 金屬電纜測(cè)試儀，所連接的3根導(dǎo)線探針長(zhǎng)度為l=(30.0±2.9)mm，探針在空氣中的理論特性阻抗為(165±2)Ω(計(jì)算方法見(jiàn)文獻(xiàn)[34])。

實(shí)際測(cè)量時(shí)，在雪坑坑壁的不同高度，將探針?biāo)讲迦脒M(jìn)行探測(cè)。由于儀器連接的問(wèn)題，在采集數(shù)據(jù)的第一天，在第一個(gè)采樣坑中，只有坑壁上半部分的數(shù)據(jù)，見(jiàn)圖3(a)。

在傳輸線中，信號(hào)傳輸?shù)膯纬套邥r(shí)Δt可以使用導(dǎo)數(shù)法進(jìn)行估算[20]，其波速的計(jì)算公式為：

(3)

此外，還在雪坑邊緣沿垂直方向進(jìn)行取芯采樣，以測(cè)量積雪的密度，并得到相對(duì)連續(xù)的隨深度變化的密度剖面。采樣使用的是圓柱形的取樣器，其底面直徑為(47.00±0.29) mm，長(zhǎng)度為(300.0±2.9) mm，每次取樣的體積大約為0.5 L。測(cè)定樣本質(zhì)量時(shí)使用的儀器為科恩HDB 5K5N型電子鉤秤，其標(biāo)稱精度為5 g?；谀Ｐ蚚16]，可以由雪坑密度隨深度變化的情況，得到波速的垂直分布廓線。

圖3為3個(gè)采樣位置的波速垂直分布廓線。其中，藍(lán)色點(diǎn)表示基于時(shí)域反射法(TDR)測(cè)量，使用式(1)計(jì)算得到的結(jié)果，黑色線表示基于取芯密度實(shí)測(cè)得到的結(jié)果。虛線表示根據(jù)式(2)得到的波速不確定度，需要指出的是，這個(gè)不確定度包含模型[16]的介電常數(shù)的不確定度和取芯密度實(shí)際測(cè)量的不確定度的影響。

圖3 3個(gè)采樣位置的波速垂直分布廓線

從圖3(a)-圖3(c)可以知道，所有采樣位置都有類似的特征，即地表附近的波速值稍高，且在深度0～150 cm的范圍內(nèi)，均有波速值隨深度增加而下降的趨勢(shì)，當(dāng)深度大于150 cm時(shí)，波速值近似不變。之前的其他研究者也得出了類似的波速垂直分布廓線(如文獻(xiàn)[1]和文獻(xiàn)[35])。

在每個(gè)采樣位置的兩組數(shù)據(jù)(即TDR結(jié)果和取芯實(shí)測(cè)結(jié)果)中，它們之間差異產(chǎn)生可以歸因于測(cè)量方法的差異。因?yàn)門(mén)DR探針在測(cè)量時(shí)是水平插入坑壁的，而對(duì)雪進(jìn)行取芯采樣時(shí)是沿著垂直方向的。因此，這兩種方法不能同時(shí)測(cè)量同一部分雪，故而會(huì)有差異。此外，兩者之間的差異主要集中于采樣位置(雪坑)的表面，因?yàn)樵谘┛颖砻妗Ｑ┑目紫抖容^大，所以在進(jìn)行采樣和測(cè)量時(shí)會(huì)出現(xiàn)較明顯的差異。

4 結(jié) 論

通過(guò)比較不同的密度(ρ)-介電常數(shù)(ε)關(guān)系模型，基于如下兩個(gè)原因，證明了Robin等給出的模型方程是其中最為可靠的一個(gè)：

①Robin模型(即模型[16])的擬合曲線位于兩個(gè)極值曲線的中間(圖1)。

②Robin模型(即模型[16])的不確定度區(qū)間囊括了其他模型的范圍(圖1)。

此外，Robin模型的不確定度的大小，與不同測(cè)量手段得到的不確定度的大小，是具有可比性的(如TDR或CMP)。在干雪上所采集到的實(shí)驗(yàn)數(shù)據(jù)，其結(jié)果也證明了Robin等建立的模型是合適的，最后的比較結(jié)果也有較好的一致性。