探地雷達(dá)在監(jiān)測(cè)地表土壤水分中的研究進(jìn)展

1、    探地雷達(dá)在監(jiān)測(cè)地表土壤水分中的研究進(jìn)展    高勝國(guó)+翁海騰+朱忠禮摘要：探地雷達(dá)（ground penetrating radar，gpr）作為監(jiān)測(cè)中小尺度或田塊尺度地表土壤含水量的一種極有潛力的技術(shù)，在近幾年來得到了快速的發(fā)展。在介紹探地雷達(dá)測(cè)量地表土壤含水量基本原理的基礎(chǔ)上，總結(jié)探地雷達(dá)在地表土壤含水量監(jiān)測(cè)中的應(yīng)用進(jìn)展，主要包括反射波法、地面波法和地表反射系數(shù)法。同時(shí)，結(jié)合國(guó)內(nèi)外最新的研究成果，分析每種方法的特點(diǎn)，討論每種方法的測(cè)量精度以及測(cè)量的代表深度。最后，對(duì)探地雷達(dá)在地表土壤含水量監(jiān)測(cè)方面的應(yīng)用進(jìn)行展望。關(guān)鍵詞：探地雷達(dá)；土壤水分；反

2、射波法；地面波法；地表反射系數(shù)法；測(cè)量深度；衛(wèi)星遙感；地面驗(yàn)證數(shù)據(jù)： s127文獻(xiàn)標(biāo)志碼： a：1002-1302（2017）12-0001-05土壤水作為水資源的重要組成部分，是陸地生態(tài)系統(tǒng)尤其是農(nóng)田生態(tài)系統(tǒng)最重要的因素之一。掌握農(nóng)田尺度土壤水分在空間上的分布，對(duì)于農(nóng)田灌溉管理等農(nóng)業(yè)生產(chǎn)實(shí)踐活動(dòng)意義重大。長(zhǎng)期以來，對(duì)農(nóng)田地表土壤含水量測(cè)定方法的研究一直被業(yè)內(nèi)人士所重視。目前為止，國(guó)內(nèi)外科研人員提出的土壤含水量測(cè)定方法多達(dá)幾十種，比較傳統(tǒng)的方法主要包括土鉆取土稱質(zhì)量法、中子儀法、電容法、時(shí)域反射法（tdr）、頻域反射儀（fdr）法等1-3。雖然以上方法在探測(cè)精度上基本可以滿足應(yīng)用的需求，但它們

3、只能探測(cè)小尺度的土壤含水量，準(zhǔn)確地講均屬于單點(diǎn)測(cè)量。由于土壤含水量空間變異性很大，這幾種方法難以直接得到土壤含水量在面上的分布數(shù)據(jù)，也不能反映土壤含水量在空間上的連續(xù)變化。如果要實(shí)施大面積土壤含水量的監(jiān)測(cè)，這些技術(shù)需要大量的觀測(cè)數(shù)據(jù)以滿足一定的采樣密度，而后空間插值，均耗時(shí)耗力。另外，這些方法均屬于有損測(cè)量，測(cè)量的同時(shí)破壞了土層結(jié)構(gòu)，難以實(shí)現(xiàn)重復(fù)測(cè)量4。隨著遙感技術(shù)的發(fā)展，應(yīng)用遙感手段觀測(cè)地表土壤含水量已廣泛開展，并且已經(jīng)取得了一些實(shí)質(zhì)性的進(jìn)展5，這使得區(qū)域尺度或大尺度土壤含水量實(shí)時(shí)動(dòng)態(tài)監(jiān)測(cè)成為可能。作為傳統(tǒng)野外土壤含水量觀測(cè)（點(diǎn)觀測(cè)）和大尺度遙感觀測(cè)的一種補(bǔ)充，探地雷達(dá)（gpr）是一種快速、

4、無損測(cè)量技術(shù)，可以用來填補(bǔ)在中小尺度或是田塊尺度上測(cè)量土壤含水量技術(shù)上的空白2。大量的研究應(yīng)用工作證明，探地雷達(dá)可以實(shí)現(xiàn)田塊尺度上土壤含水量的觀測(cè)，并且在保證原位測(cè)定精度的同時(shí)，快速、便捷地獲取較大區(qū)域的土壤含水量信息，大大提高傳統(tǒng)區(qū)域范圍土壤含水量觀測(cè)的工作效率，也為衛(wèi)星遙感土壤監(jiān)測(cè)提供地面驗(yàn)證數(shù)據(jù)。1探地雷達(dá)監(jiān)測(cè)土壤水分理論基礎(chǔ)探地雷達(dá)是一種地球物理無損探測(cè)技術(shù)，從20世紀(jì)70年代開始逐漸被應(yīng)用到水文和農(nóng)業(yè)領(lǐng)域當(dāng)中。從原理上來講，探地雷達(dá)技術(shù)和地震波法、聲納法相似，由主機(jī)產(chǎn)生一定頻率（1 mhz1 ghz）的電磁波，由發(fā)射天線向地下介質(zhì)發(fā)射脈沖信號(hào)，然后由接收天線接收和處理回波信號(hào)，這樣回

5、波信號(hào)的雙程傳播時(shí)間、電磁場(chǎng)強(qiáng)度和波形等特征會(huì)反映脈沖信號(hào)所通過介質(zhì)的電磁特性以及幾何形態(tài)變化等特征6-7。在探地雷達(dá)的回波信號(hào)中可以識(shí)別及利用的回波類型主要有空氣直達(dá)波、地面直達(dá)波和反射回波3種（圖1）?？諝庵边_(dá)波為發(fā)射天線發(fā)出的脈沖信號(hào)通過空氣直接被接收天線所接收，由于電磁波在空氣中衰減較小，因此該回波強(qiáng)度較大，在應(yīng)用中可以用其代表發(fā)射脈沖信號(hào)的強(qiáng)度。地面直達(dá)波為在接觸地表測(cè)量過程中，由發(fā)射天線端沿地表直線傳播到接收天線端的脈沖信號(hào)，反射回波是接收到的介質(zhì)中反射界面的反射波，利用這些回波可以分析傳播介質(zhì)的電磁特性。利用探地雷達(dá)測(cè)量地表土壤含水量的主要理論基礎(chǔ)可以概括為：（1）探地雷達(dá)信號(hào)傳

6、播速度或反射率可以反映土壤介電常數(shù)特征。（2）常溫下，干土壤的相對(duì)介電常數(shù)是4，空氣的相對(duì)介電常數(shù)是1，淡水的相對(duì)介電常數(shù)是81，淡水的相對(duì)介電常數(shù)遠(yuǎn)大于空氣和干土壤的相對(duì)介電常數(shù)，則土壤含水量可以對(duì)土壤的相對(duì)介電常數(shù)產(chǎn)生很大的影響，那么可以通過測(cè)定介電常數(shù)來反映土壤含水量9-11。（3）在土壤含水量研究及應(yīng)用領(lǐng)域最經(jīng)典的土壤體積含水量（）和介電常數(shù)（）的關(guān)系是topp等提出的topp公式12：=-0.053+0.029 3-0.000 552+0.000 004 33。（1）2探地雷達(dá)監(jiān)測(cè)土壤含水量應(yīng)用進(jìn)展按照天線端接收回波類型，將現(xiàn)階段常用的探地雷達(dá)監(jiān)測(cè)土壤含水量的方法分為反射波法、地面波

7、法和地表反射系數(shù)法3種。2.1反射波法反射波法是依據(jù)土壤介質(zhì)中反射波的傳播距離和雙程走時(shí)來計(jì)算電磁波在土壤中的傳播速度，進(jìn)而確定土壤的相對(duì)介電常數(shù)，最后根據(jù)topp公式確定土壤含水量13。根據(jù)測(cè)量過程中發(fā)射天線和接收天線相對(duì)位置的差異，可以將這種方法分為固定天線距模式和變天線距模式。2.1.1固定天線距模式固定天線距模式測(cè)量土壤含水量一般適用于土壤中有明顯的反射界面，而且可以知道或估計(jì)反射界面的深度。測(cè)量過程中收發(fā)天線在一定的天線距下同時(shí)向同一方向移動(dòng)，每次移動(dòng)間隔進(jìn)行1次測(cè)量，保證連續(xù)測(cè)量。如圖2所示，假設(shè)土壤中存在明顯的反射面，反射層深度為h，以天線間距為x進(jìn)行測(cè)量，測(cè)得的雷達(dá)回波的雙程走

8、時(shí)為t，這樣可以計(jì)算得到雷達(dá)信號(hào)波在土壤中的平均傳播速度v：v=h2+（x/2）2/t。（2）則土壤的相對(duì)介電常數(shù)r可表達(dá)為：r=c2v2=4c2t24h2+x2。（3）stoffregen等基于此方法利用1 ghz天線測(cè)定蒸滲儀箱體中沙質(zhì)土壤的含水量，經(jīng)過對(duì)比驗(yàn)證，發(fā)現(xiàn)此方法測(cè)定的土壤體積含水量和蒸滲儀測(cè)量數(shù)據(jù)的標(biāo)準(zhǔn)差為0.01 m3/m314。lunt等利用此測(cè)量方法，并且利用試驗(yàn)站定標(biāo)后的介電常數(shù)和土壤體積含水量的關(guān)系計(jì)算得到gpr土壤體積含水量，同時(shí)和中子儀的測(cè)量結(jié)果相比較得到的均方根誤差為 0.018 m3/m315。shengguo等比較了一定采樣體積下，gpr測(cè)量土壤含水量數(shù)據(jù)和

9、稱質(zhì)量法的測(cè)量結(jié)果，得到均方根誤差為0.05 m3/m316。endprint固定天線距模式操作方便、計(jì)算簡(jiǎn)單，它得到的土壤相對(duì)介電常數(shù)為地面到反射面路徑上的相對(duì)介電常數(shù)的平均值，因此該方法得到的土壤含水量可以代表測(cè)點(diǎn)周圍的局部區(qū)域情況，這種方法的使用前提是反射層面存在且已知深度，如可以利用農(nóng)田的犁底層作為反射面。2.1.2變天線距模式在土壤層中存在連續(xù)反射面的情況下，可通過變天線距模式來估計(jì)土壤體積含水量，變天線距離模式又可以分為共中心點(diǎn)模式（common-midpoint，cmp）和廣角模式（wide angle reflection and refraction，warr）。共中心點(diǎn)模式

10、即雷達(dá)系統(tǒng)收發(fā)天線分離，測(cè)量時(shí)，收發(fā)天線分別向兩邊等距移動(dòng)，在每次移動(dòng)后進(jìn)行數(shù)據(jù)獲取，探測(cè)示意圖如圖3-a所示。廣角測(cè)量模式也采用收發(fā)天線分離，測(cè)量時(shí)與共中心點(diǎn)測(cè)量相類似，只是其中的一個(gè)天線固定，另一個(gè)天線向一邊等間距移動(dòng)，其測(cè)量示意圖如圖3-b。假設(shè)2次測(cè)量中發(fā)射、接收天線的水平間距分別為x1和x2（x1x2），測(cè)量相應(yīng)的雙程進(jìn)行時(shí)間，分別為t1和t2。土壤中電磁波速可表示為：v=x22-x21t21-t20。（4）那么可以得到土壤媒質(zhì)的相對(duì)介電常數(shù)為：r=c2v2=c2（t21-t20）x22-x21。（5）grote等比較了稱質(zhì)量法和cmp法土壤含水量測(cè)量結(jié)果，發(fā)現(xiàn) 900 mhz 天線

11、的數(shù)據(jù)結(jié)果和稱質(zhì)量法相關(guān)性（0.98）比450 mhz天線結(jié)果的線性相關(guān)性（0.92）強(qiáng)，并且它們都有很小的均方根誤差（分別為0.017、0.011 m3/m3）17。此外，huisman等利用地面回波的斜率估計(jì)土壤體積含水量，225 mhz 天線的精度為±0.024 m3/m3，450 mhz天線為±0.025 m3/m318。利用變天線距模式測(cè)量土壤含水量的優(yōu)點(diǎn)是它不須要知道地下反射層的深度信息，只須接收到雷達(dá)的反射數(shù)據(jù)即可，但是這種方法也存在一些局限，即它和固定天線距模式一樣，必須是在地下反射層存在的條件下才可以進(jìn)行；變天線距模式在每一個(gè)位置需要多次測(cè)量才可以估算土壤

12、含水量；在非均勻土壤介質(zhì)中，共中心點(diǎn)模式所得測(cè)量結(jié)果只能反映靠近共中心點(diǎn)的區(qū)域土壤含水量，廣角測(cè)量模式的測(cè)量結(jié)果反映靠近固定的天線附近區(qū)域土壤含水量19。2.2地面波法在雷達(dá)天線接觸地面測(cè)量模式下，地面波什發(fā)射脈沖在地面散射后剩余的一部分能量沿著土壤表層由發(fā)射天線傳播到接收天線，它可以反映表層土壤的電磁特性。這種方法在土壤下沒有明顯反射層的情況下，仍然可以使用，地面波的傳播如圖1所示。由于地面波具有能量容易耗散的特點(diǎn)，為了使測(cè)量精度更高，測(cè)量過程中須將接收天線和發(fā)射天線盡量緊挨土壤表層。在利用地面波估計(jì)土壤體積含水量時(shí)，可以利用廣角模式、共中心點(diǎn)模式和固定天線距模式3種測(cè)量模式，其中廣角模式（

13、warr）和共中心點(diǎn)模式（cmp）測(cè)量數(shù)據(jù)空間分辨率低而且在每個(gè)位置上測(cè)量所需時(shí)間較長(zhǎng)，固定天線距模式（fo）具有更高空間分辨率且節(jié)省測(cè)量時(shí)間，更適合大區(qū)域土壤含水量快速監(jiān)測(cè)。du等提出了利用地面波法測(cè)量土壤體積含水量的程序：（1）利用變天線距法估計(jì)不同的天線間距下地面波的傳播時(shí)間；（2）選擇一個(gè)能夠使地面波清晰的區(qū)別于空氣波和反射波的天線間距；（3）在這種天線距離下利用固定天線距法建立地面波傳播時(shí)間和介電常數(shù)之間的關(guān)系8，20。在利用地面波測(cè)量地表土壤含水量過程中，須要估算電磁脈沖在地表的傳播時(shí)間，因此地面波傳播時(shí)間的精度對(duì)于準(zhǔn)確估計(jì)土壤含水量非常重要。在地面波傳播時(shí)間的估計(jì)中一般會(huì)有12

14、ns的誤差，其主要來源有：（1）波的前緣難以被識(shí)別；（2）零時(shí)刻的定標(biāo)不一定準(zhǔn)確；（3）地面波傳播距離的不確定。galagedara等中關(guān)于零時(shí)刻的定標(biāo)如下：假設(shè)在一定天線距離x（m）下，計(jì)算得到的空氣波傳播時(shí)間為tair，觀測(cè)得到的空氣波傳播時(shí)間為taw21，這樣零時(shí)刻校準(zhǔn)量t0為：t0=tair-taw。（6）假設(shè)觀測(cè)得到地波的傳播時(shí)間為tgw，則校準(zhǔn)后的地波傳播時(shí)間tab為：tab=tgw-t0。（7）地面波在地表土壤中傳播的速度v可以表示為：v=xtab。（8）土壤相對(duì)介電常數(shù)可以表達(dá)為：r=t2abc2x2。（9）另外，sperl等通過公式推導(dǎo)直接給出了比較簡(jiǎn)單的介電常數(shù)和空氣波、地

15、面波傳播時(shí)間之間的關(guān)系8：=cv2=c（tgw-taw）+xx2。（10）huisman等利用225 mhz的天線以及固定天線距法中的地面波估計(jì)地表土壤含水量，并且和tdr數(shù)據(jù)作對(duì)比，指出gpr的誤差為±0.036 m3/m317。grote等通過試驗(yàn)比較固定天線距法地面波估計(jì)的土壤含水量和稱質(zhì)量法估計(jì)的土壤含水量，對(duì)于900 mhz的天線均方根誤差為0.11 m3/m316，對(duì)于450 mhz天線均方根誤差為0.17 m3/m3，并且最大的誤差出現(xiàn)在土壤非常干燥的情況下。galagedara等通過和tdr測(cè)量土壤含水量數(shù)據(jù)相比較，指出100 mhz天線的地面波所估計(jì)的土壤含水量應(yīng)該

16、主要是地表10 cm以內(nèi)的平均含水量，并且用空氣波進(jìn)行零時(shí)刻定標(biāo)比較精確，這樣得到土壤含水量誤差在0.01 m3/m3以內(nèi)21。另外，探地雷達(dá)的天線頻率和探測(cè)深度始終是1對(duì)矛盾關(guān)系，天線頻率越高，探測(cè)深度越淺。但采用地面波法測(cè)定土壤含水量時(shí)，由于測(cè)定的是土壤表層含水量，不需要很深的探測(cè)深度，所以應(yīng)該優(yōu)先采用較高的天線頻率，天線頻率越高，土壤含水量的測(cè)定精度相對(duì)越高。目前，利用地面波法估計(jì)土壤含水量也存在如下缺陷：（1）將地面波區(qū)別于回波中的反射波和折射波比較困難，尤其是固定天線距法；（2）難以確定合適的天線間距；（3）由于地面波的衰減比較快，這樣就限制了天線分離距離的可變化范圍22-23。en

17、dprint2.3地表反射系數(shù)法地表反射系數(shù)法是將天線置于空氣-土壤分界面上方（圖4），每次測(cè)量可以估計(jì)雷達(dá)足跡范圍內(nèi)的土壤含水量。為了測(cè)量方便，可以將天線架設(shè)在可移動(dòng)裝置或是低空飛行平臺(tái)上實(shí)現(xiàn)區(qū)域內(nèi)的快速覆蓋。這種方法利用了空氣-土壤分界面的振幅反射系數(shù)r（r可以表示為反射回波信號(hào)的振幅與理想完全反射面的反射回波信號(hào)振幅的比值），這里的理想界面可以用比雷達(dá)足跡面積大的鐵板來代替，也可以利用空氣直達(dá)波振幅代替鐵板的回波信號(hào)振幅使用，而振幅反射系數(shù)和土壤介電常數(shù)soil之間的關(guān)系可表示為24-25：r=1-soil1+soil。（11）如果應(yīng)用空氣直達(dá)波作為理想完全回波，土壤介電常數(shù)可以利用反射

18、波振幅和空氣直達(dá)波振幅表示為25-26：soil=1+ar/am1-ar/am2。（12）式中：ar為反射波振幅；am為空氣直達(dá)波振幅。由于回波信號(hào)受直流偏移的或雜波的影響而難以消除，在一些研究中認(rèn)為將回波中最大波峰值和最小波谷值之間的差定義為振幅a更為合理（圖5）8，27。利用此方法估計(jì)得到土壤相對(duì)介電常數(shù)后，同理可以根據(jù)topp公式計(jì)算得到土壤體積含水量。huisman等指出，利用這種方法得到的地表土壤含水量和0.20 m長(zhǎng)度探針的tdr數(shù)據(jù)非常相似，但土壤含水量測(cè)量結(jié)果在短距離上的差異非常大；同時(shí)分析了可能導(dǎo)致這種差異的3種解釋：（1）土壤含水量隨土層深度的變化會(huì)影響反射系數(shù)；（2）地表

19、粗糙度會(huì)對(duì)反射系數(shù)產(chǎn)生影響；（3）振幅測(cè)量過程中存在誤差26。redman等又指出，現(xiàn)在的gpr設(shè)備對(duì)振幅的測(cè)量完全可以達(dá)到要求的精度，所以振幅的測(cè)量誤差不是導(dǎo)致土壤含水量觀測(cè)誤差的主要原因，而地表粗糙度和土壤含水量廓線（也可以認(rèn)為是探測(cè)深度問題）對(duì)反射系數(shù)的影響成為當(dāng)前研究的重點(diǎn)27。weihermüller等指出，面反射法所反映的地表土壤含水量空間上的差異要比實(shí)際小，最可能的原因是面反射法對(duì)表面的土壤更加敏感，而表層干燥土壤的相對(duì)均一分布會(huì)導(dǎo)致土壤含水量測(cè)量結(jié)果空間差異變小28。redman等數(shù)字模擬土壤含水量的層狀分布（上層干燥下層濕潤(rùn)、上層濕潤(rùn)下層干燥）對(duì)gpr測(cè)量的影響，結(jié)

20、果表明，土壤含水量層狀結(jié)構(gòu)會(huì)嚴(yán)重影響探測(cè)到的含水量的大小，并且表層的濕度會(huì)影響探測(cè)的深度25。要應(yīng)用面反射法檢測(cè)區(qū)域土壤含水量，則須要嚴(yán)格考慮地表粗糙度、表層干燥土壤的影響。在植被覆蓋區(qū)域，還須要考慮植被覆蓋的影響。目前這種方法還沒有被大量地實(shí)踐應(yīng)用。3探地雷達(dá)監(jiān)測(cè)土壤含水量的代表深度探地雷達(dá)監(jiān)測(cè)地表土壤含水量所代表的深度一直是此應(yīng)用研究的一個(gè)瓶頸，也是一個(gè)研究熱點(diǎn)。在常用的幾種方法中，反射波法所代表的深度由反射層的深度決定，目前大部分的研究是針對(duì)地面波法和地表反射系數(shù)法。探地雷達(dá)測(cè)量土壤含水量所代表的深度受土壤基質(zhì)、土壤含水量廓線以及所使用的天線頻率的影響，由于土壤基質(zhì)及土壤含水量的空間異質(zhì)

21、性，探地雷達(dá)測(cè)量土壤含水量的深度仍處于數(shù)字模擬以及條件試驗(yàn)的探究階段29-30。redman等通過數(shù)字模擬的方法探究地表反射系數(shù)法探測(cè)的深度，指出土壤含水量的層狀分布會(huì)嚴(yán)重影響反射系數(shù)法所測(cè)得的土壤含水量，須要進(jìn)一步改進(jìn)模型反映地表散射和采樣體積內(nèi)土壤含水量空間分布差異的影響25。galagedara等通過試驗(yàn)的方法將灌溉過程中g(shù)pr地面波法得到的土壤含水量數(shù)據(jù)和tdr數(shù)據(jù)相比較發(fā)現(xiàn)，對(duì)于100 mhz的天線，gpr所得到的土壤含水量數(shù)據(jù)和探針長(zhǎng)為0.1 m的tdr土壤含水量數(shù)據(jù)最接近，并且回歸分析后沒有明顯的差異；對(duì)于450 mhz的天線而言，gpr所測(cè)量的土壤含水量大于試驗(yàn)所使用幾種長(zhǎng)度t

22、dr探針測(cè)量值，但是和最淺層（20 cm）的tdr測(cè)量值最接近21。galagedara等通過gprmax2d軟件模擬了在只有2個(gè)不同含水量土壤層的土壤剖面條件下地面波測(cè)量土壤含水量的試驗(yàn)，通過改變上層土壤層的厚度來測(cè)試gpr地面波的探測(cè)深度，結(jié)果表明gpr地面波在上層土壤干燥下層土壤濕潤(rùn)的模式下的探測(cè)深度比上層土壤濕潤(rùn)下層土壤干燥模式的探測(cè)深度大；gpr地面波的探測(cè)深度和所使用頻率呈負(fù)相關(guān)，并且給出了通過試驗(yàn)擬合的上層干燥下層濕潤(rùn)和上層濕潤(rùn)下層干燥2種模式下gpr探測(cè)深度和頻率的線性方程，且具有很高的相關(guān)性；同時(shí)發(fā)現(xiàn)在450 mhz天線條件下，不同的電導(dǎo)率并沒有對(duì)gpr地面波的探測(cè)深度產(chǎn)生明

23、顯的影響31。galagedara等通過gpr數(shù)據(jù)和一系列探針長(zhǎng)度的tdr數(shù)據(jù)作比較，指出在灌溉過程中，地面直達(dá)波的采樣深度在距離地表0.10.5 m之間；在同一試驗(yàn)中還利用土壤濕潤(rùn)峰移動(dòng)的理論來研究gpr在灌溉過程中的探測(cè)深度，并指出450 mhz天線在灌溉條件下的探測(cè)深度小于0.56 m32。grote等也通過試驗(yàn)的方法探究250、500、1 000 mhz的天線在沙土中的探測(cè)深度，試驗(yàn)結(jié)果表明，地面波的探測(cè)深度和使用的頻率成正比，在濕潤(rùn)土壤中的探測(cè)深度要比干燥土壤小一些，但是探測(cè)深度對(duì)土壤濕度的依賴程度遠(yuǎn)比通過數(shù)字模擬預(yù)期的?。辉谠撛囼?yàn)中探測(cè)的最小深度出現(xiàn)在利用1 000 mhz天線在濕

24、潤(rùn)土壤中測(cè)量的情況下，深度為12 cm，探測(cè)的最大深度出現(xiàn)在利用250 mhz在干燥土壤測(cè)量的情況下，深度為30 cm33。通過控制條件試驗(yàn)以及大量的農(nóng)田試驗(yàn)給出探地雷達(dá)測(cè)量地表土壤含水量參考深度的可行途徑。4總結(jié)與展望理論研究和大量的應(yīng)用實(shí)踐表明，探地雷達(dá)（gpr）是測(cè)量區(qū)域土壤含水量的潛有力工具?，F(xiàn)階段基于探地雷達(dá)測(cè)量地表土壤含水量的方法中，反射波法是一種發(fā)展相對(duì)成熟且應(yīng)用較多的方法。而在農(nóng)業(yè)和水文領(lǐng)域地面波法的研究和應(yīng)用更多，尤其在農(nóng)田尺度上開展土壤含水量監(jiān)測(cè)，因?yàn)榈孛娌ǚú恍枰叵掠忻黠@的反射層作為輔助條件，地面波法比反射波法更容易推廣使用。地表反射系數(shù)法測(cè)量地表土壤含水量的研究和應(yīng)用

25、還相對(duì)較少，但是這種方法具有更加獨(dú)特的優(yōu)勢(shì)。因?yàn)檫@種方法是和地表非接觸測(cè)量，可以憑借地面移動(dòng)裝置或飛行平臺(tái)更加容易實(shí)現(xiàn)區(qū)域土壤含水量監(jiān)測(cè)，并且可以應(yīng)用到地表有低矮植被覆蓋的環(huán)境條件下，適合于混合下墊面的測(cè)量。同時(shí)，鑒于探地雷達(dá)地表反射系數(shù)法可以更加快速獲取中小尺度區(qū)域土壤含水量空間數(shù)據(jù)，且地表反射系數(shù)法的機(jī)理和星（機(jī)）載主動(dòng)微波遙感有相似之處，這種方法可以發(fā)展成為目前遙感反演土壤含水量的理想地面驗(yàn)證手段14，34。endprint近幾年來，探地雷達(dá)監(jiān)測(cè)地表土壤含水量的研究應(yīng)用工作取得了實(shí)質(zhì)性的進(jìn)展，但在農(nóng)業(yè)、水文等領(lǐng)域的大范圍應(yīng)用推廣仍有較大的改進(jìn)空間，有多個(gè)技術(shù)瓶頸問題須要解決如地表粗糙度的

26、影響，非均勻土壤介質(zhì)條件下獲取的土壤含水量所代表的深度、空間分辨率以及植被覆蓋的影響等。另外，為了滿足實(shí)際應(yīng)用的需求，還須要在包括topp公式在內(nèi)的經(jīng)驗(yàn)或半經(jīng)驗(yàn)的-關(guān)系基礎(chǔ)上，針對(duì)不同探地雷達(dá)天線頻率、不同土壤基質(zhì)等建立相應(yīng)的經(jīng)驗(yàn)或物理-關(guān)系模型35。傳統(tǒng)的小尺度（點(diǎn)測(cè)量）測(cè)量數(shù)據(jù)和探地雷達(dá)測(cè)量的土壤含水量數(shù)據(jù)之間，由于觀測(cè)尺度的差異會(huì)導(dǎo)致明顯的驗(yàn)證誤差，因此也有必要設(shè)計(jì)更多的控制條件試驗(yàn)為gpr技術(shù)提供更加可信的驗(yàn)證途徑。隨著技術(shù)的不斷進(jìn)步、試驗(yàn)條件的改善以及國(guó)內(nèi)外研究應(yīng)用經(jīng)驗(yàn)的不斷積累，探地雷達(dá)可以成為中小尺度測(cè)量地表土壤含水量的一種更加有效的手段。參考文獻(xiàn)：1邵明安，王全九，黃明斌. 土

27、壤物理學(xué)m. 北京：高等教育出版社，2006：57-61.2evett s r，parkin g w. advances in soil water content sensingj. vadose zone journal，2005，4（4）：986-991.3topp g p. measuring soil water contentj. encyclopedia of hydrological sciences，2006，5（2）：1077-1089.4robinson d a，campbell c s，hopmans j w，et al. soil moisture measureme

28、nt for ecological and hydrological watershed-scale observatories：a reviewj. vadose zone journal，2008，7（1）：358-389.5仝兆遠(yuǎn)，張萬昌. 土壤水分遙感監(jiān)測(cè)的研究進(jìn)展j. 水土保持通報(bào)，2007，27（4）：107-113.6雷少剛，卞正富. 探地雷達(dá)測(cè)定土壤含水率研究綜述j. 土壤通報(bào)，2008，39（5）：1179-1183.7schmalz b，lennartz b. analyses of soil water content variations and gpr attribu

29、te distributionsj. journal of hydrology，2002，267（3/4）：217-226.8sperl c. determination of spatial and temporal variation of the soil water content in an agro-ecosystem with ground-penetrating radard. munich：technische universitat münchen，1999.9benedetto a. water content evaluation in unsaturated

30、 soil using gpr signal analysis in the frequency domainj. journal of applied geophysics，2010，71（1）：26-35.10曾昭發(fā)，劉四新，王者江，等. 探地雷達(dá)方法原理及應(yīng)用m. 北京：科學(xué)出版社，2006：36-49.11粟毅，黃春林，雷文太. 探地雷達(dá)理論語(yǔ)應(yīng)用m. 北京：科學(xué)出版社，2006：23-29.12topp g c，davis j l，annan a p. electromagnetic determination of soil water content：measurements i

31、n coaxial transmission linesj. water resources research，1980，16（3）：574-582.13causse e，senechal p. model-based automatic dense velocity analysis of gpr field data for the estimation of soil propertiesj. journal of geophysics and engineering，2006，3（2）：169-176.14stoffregen h，zenker t. wessolek g. accur

32、acy of soil water content measurements using ground penetrating radar：comparison of ground penetrating radar and lysimeter dataj. journal of hydrology，2002，267（3/4）：201-206.15lunt i a，hubbard s s，rubin y.soil moisture content estimation using ground-penetrating radar reflection dataj. journal of hyd

33、rology，2005，307（1/2/3/4）：254-269.16gao s g，zhu z l，jiao q s，et al. application of gpr surface reflection method for estimating soil water content of bare soilc.14th international conference on ground penetrating radar，2012：639-645.endprint17grote k，hubbard s，rubin y. field-scale estimation of volume

34、tric water content using ground-penetrating radar ground wave techniquesj. water resources research，2003，39（11）：1321.18huisman j a，sperl c，bouten w，et al. soil water content measurements at different scales：accuracy of time domain reflectometry and ground-penetrating radarj. journal of hydrology，200

35、1，245（1/2/3/4）：48-58.19huisman j a，bouten w. accuracy and reproducibility of measuring soil water content with the ground wave of ground penetrating radarj. geophys，2003，8：65-73.20du s. determination of water content in the subsurface with the ground wave of ground penetrating radard. munich：ludwig

36、maximilians universitat，1996.21galagedara l w，parkin g d. an analysis of the ground-penetrating radar direct ground wave method for soil water content measurementj. hydrological processes，2003，17（18）：3615-3628.22朱安寧，吉麗青. 基于探地雷達(dá)的土壤水分測(cè)定方法研究進(jìn)展j. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2009，17（5）：1039-1044.23oden c p，olhoeft g r，wrig

37、ht d l，et al. measuring the electrical properties of soil using a calibrated ground-coupled gpr systemj. vadose zone journal，2008，7（1）：171-183.24charlton m b. principles of ground-penetrationg radar for soil moisture assessment m/ol. 2008：http：/25redman d g l，parkin g. measuring soil water content w

38、ith the ground penetrating radar surface reflectivity method：effects of spatial variabilityj. asae meeting paper no.032276，2003.26lambot s，weihermüller l，huisman j a，et al. analysis of air-launched ground-penetrating radar techniques to measure the soil surface water contentj. water resources research，2006，42（11）：2526-2528.27huisman j a，hubbard s s，redman j d，et al. measuring soil water content with ground penetrating radar：a reviewj. vadose zone journal，2003，2（4）：476-491.28weihermueller l，huisman j a，lambot s，