典型區域氣候下的D-InSAR大氣形變誤差分析

孫方 朱慶林 劉琨 趙振維 康士峰 董翔

(中國電波傳播研究所,青島 266107)

?

典型區域氣候下的D-InSAR大氣形變誤差分析

孫方 朱慶林 劉琨 趙振維 康士峰 董翔

(中國電波傳播研究所,青島 266107)

大氣結構常數;大氣湍流;對流層延遲;差分干涉合成孔徑雷達;形變誤差

DOI 10.13443/j.cjors.2016032802

引 言

合成孔徑雷達差分干涉(Differential Interferometry Synthetic Aperture Radar,D-InSAR)技術是在合成孔徑雷達干涉(Interferometry Synthetic Aperture Radar,InSAR)技術的基礎上發展起來的．其基本原理是:利用重復軌觀測獲取的干涉相位,通過差分處理去除2次觀測相位中的共有量(平地效應、地形相位和大氣延遲等),可以得到形變相位,進而獲得高程形變量．由于干涉相位對微小形變極其敏感,毫米級的形變在干涉相位中都會有所反映,因此,D-InSAR技術可以監測地表的微小形變,被廣泛應用于地震、泥石流、火山爆發等自然災害的預報與監測研究中,具有廣闊的應用前景．而隨著經濟建設的需要,如高鐵建設沿線地面沉降災害的調查、滑坡及危巖體監測等,對D-InSAR的監測精度提出了更高的要求,但是由于大氣延遲相位的存在,制約了D-InSAR監測精度的進一步提高,因此準確實時地預報監測地表處的大氣環境特性,評估大氣延遲相位的影響對提高D-InSAR監測精度至關重要．

雷達信號穿越大氣層時,不同密度的大氣介質折射率的變化對信號傳播會產生傳播路徑彎曲和速度延遲,從而產生大氣延遲相位．大氣延遲相位分為兩種：一種是由于大氣湍流造成的,它與地形不相關;另一種是由對流層的垂直分層變化造成的,它與地形相關．由對流層的垂直分層變化造成的大氣延遲可以通過準確測量地面的氣象參數消除絕大部分;而湍流效應是由大氣水汽隨時間和區域擾動變化造成的,具有很大的隨機性和不均一性,因此很難預測．研究表明,大氣的相對濕度在空間和時間上20%的變化可以導致10～14 cm的形變測量誤差[1-2],大氣湍流效應的不確定性已經成為限制D-InSAR測量精度進一步提高的瓶頸問題之一,只有準確獲取D-InSAR成像時間內的大氣湍流強度,才能進一步提高大氣延遲相位的修正精度．

1 三軌法形變誤差

圖1 D-InSAR三軌法干涉測量示意圖

令兩次觀測時所經歷的大氣延遲量分別為δρ、δρ′,則大氣干涉相位可表示為[2]:

(1)

φatm的隨機變化對D-InSAR提取地形和探測地表形變有很大干擾,甚至可能掩蓋真實的形變信號．假定D-InSAR系統中每幅干涉圖中的大氣延遲標準差δφatm相同,則大氣延遲相位對于三軌法差分干涉測量的形變量標準差可表示為[3]:

(2)

2 大氣結構常數與干涉相位

圖2 北京地區不同對應的天頂延遲隨時間的變化

根據D-InSAR在不同時間測得的干涉相位結構函數與大氣結構常數之間的關系,與時間相關的相位結構函數可以寫成[9]:

=〈|ΔS(t)-ΔS(t+Δt)|2〉.

(3)

(4)

(5)

式中：d為基線距離;

(6)

(7)

3 計算分析結果

表1 典型區域在不同時間的對比

圖3 不同InSAR衛星形變誤差隨的變化

由圖3可知,SAR的基線長度越大,大氣湍流導致的形變誤差越大．對于ERS-1/2衛星,上述典型區域采樣樣本中的大氣湍流強度都能滿足形變誤差小于厘米級的高精度指標;對于TerraSAR衛星,北京夏季、上海春季的大氣環境導致的形變誤差很有可能大于厘米級;而對于ALOS衛星,除了北京和烏魯木齊的冬季,絕大部分區域在采樣時間的大氣湍流強度都會導致較大的形變誤差,普遍高于厘米級,基本都在1～2 cm之間,上海的春季甚至達到了4 cm．以上只是分析了正常大氣環境下的形變誤差,當遇到極端的天氣條件,如暴風雨等強烈的對流運動時,大氣結構常數可達到10-12m-2/3的量級,此時大氣導致的形變誤差甚至能達到十幾厘米,嚴重影響D-InSAR的測量性能．

4 結 論

[1]李陶.重復軌道星載SAR差分干涉監測地表形變研究[M].武漢：武漢大學出版社, 2004.

[2]余景波, 劉國林, 曹振坦, 等.大氣延遲誤差對In-SAR數據處理影響的定量分析[J].空間科學學報, 2012, 32(5)：720-729.

YUJB,LIUGL,CAOZT,etal.quantitativeanalysisofeffectofatmosphericdelayerroroninSARdataprocessing[J].Chinesejournalofspacescience, 2012, 32(5)：720-729.(inChinese)

[3]CHARLESW,URSW,ANDREASW.Interferometricpointtargetanalysiswithjers-1l-bandSARdata[C]//2003IEEEInternationalGeoscienceandRemoteSensingSymposium.Toulouse,France, 21-25July, 2003.

[4]BIANCOL,WILCZAKJM.Convectiveboundarylayerdepth：improvedmeasurementbyDopplerradarwindprofilerusingfuzzylogicmethods[J].Journalofatmospheric&oceanictechnology, 2002, 19(11)：1745-1758.

[5]NILSSONT,HAASR.Impactofatmosphericturbulenceongeodeticverylongbaselineinterferometry[J].Journalofgeophysicalresearch, 115(10), 2010:1-11.

[6]朱慶林, 趙振維, 吳振森.精密單點定位方法測量對流層天頂延遲的精度改善[J].武漢大學學報(信息科學版), 2009, 34(9)：1098-1101.

ZHUQL,ZHAOZW,WUZS.Precisionimprovementoftroposphericzenithpathdelayestimationbyprecisepointpositioning[J].GeomaticsandInformationScienceofWuhanUniversity, 2009,34(9),1098-1101.(inChinese)

[7]林樂科.利用GNSS信號的地基大氣折射率剖面反演技術研究[D].南京：南京郵電大學, 2011.

LINLK.Profilingtroposphericrefractivitybasedonsingleground-basedGPSreceiver[D].Nanjing：NanjingUniversityofPostsandTelecommunications, 2011.(inChinese)

[8]朱慶林, 吳振森, 趙振維, 等.單臺地基衛星導航接收機測量對流層斜延遲[J].電波科學學報, 2010, 25(1)：37-40.

ZHUQL,WUZS,ZHAOZW,etal.Troposphericslantdelaymeasuredbysingularground-basedsatellitenavigationreceiver[J].Chinesejournalofradioscience,2010, 25(1)：37-40.(inChinese)

[9]NESSELJA,MANNINGRM.DerivationofmicrowaverefractiveindexstructureconstantoftheatmospherefromK-bandinterferometricphasemeasurements[J].IEEEtransactionsonantennasandpropagation,2014, 62(11)：5590-5598.

[10]BERARDINOP,FORNAROG,LANARIR,etal.AnewalgorithmforsurfacedeformationmonitoringbasedonsmallbaselinedifferentialSARinterferograms[J].IEEEtransactionsongeoscienceandremotesensing, 2002, 40(11)：2375-2383.

[11]LIZ,MULLERJP,CROSSP.Comparisonofprecipitablewatervaporderivedfromradiosonde,GPS,andmoderate-resolutionimagingspectroradiometermeasurements[J].Journalofgeophysicalresearch, 2003, 108 (20)：4651.

[12]程顯海, 張玉生.基于氣象探空數據計算的大氣折射率誤差分析[J].電波科學學報, 2014, 29(5)：951-956.

CHENGXH,ZHANGYS.Analysisoftheerrorsinestimatingatmosphericrefractiveindexfromradiosondemeasurements[J].Chinesejournalofradioscience, 2014, 29(5)：951-956.(inChinese)

[13]熊皓.電磁波傳播與空間環境[M].北京:電子工業出版社,2004.

[14]黃捷.電波大折射誤差修正[M].北京：國防工業出版社, 2000.

[15]潘中偉, 劉成國.對流層折射率起伏的基本規律[J].電波科學學報, 1997, 12(1)：89-95.

PANZW,LIUCG.Thecharacteristicsoftroposphericrefractivityfluctuation[J].Chinesejournalofradioscience, 1997, 12(1)：89-95.(inChinese)

[16]BONFORTEA,FERRETTIA,PRATIC,etal.CalibrationofatmosphericeffectsonSARinterferogramsbyGPSandlocalatmospheremodels：firstresults[J].Journalofatmosphericandsolar-terrestrialphysics, 2001, 63(12)：1343-1357.

[17]楊成生.差分干涉雷達測量技術中水汽延遲改正方法研究[D].西安：長安大學, 2011.

YANGCS.Studyonatmosphericdelaycorrectionmethodsindifferentialinterferometricsyntheticapertureradar[D].Xi’an：Chang’anUniversity, 2011.(inChinese)

D-InSAR atmospheric deformation error analysis in representative areas with different climate

SUN Fang ZHU Qinglin LIU Kun ZHAO Zhenwei KANG Shifeng DONG Xiang

(ChinaResearchInstituteofRadiowavePropagation,Qingdao266107,China)

atmospheric structure constant; atmospheric turbulence; troposphere delay; differential SAR interferometry; deformation error

孫方, 朱慶林, 劉琨, 等.典型區域氣候下的D-InSAR大氣形變誤差分析[J].電波科學學報,2016,31(5):1023-1027.

10.13443/j.cjors.2016032802

SUN F, ZHU Q L, LIU K, et al.D-InSAR atmospheric deformation error analysis in representative areas with different climate [J].Chinese journal of radio science,2016,31(5):1023-1027.(in Chinese).DOI：10.13443/j.cjors.2016032802

2016-03-28

國家自然科學基金(批準號：41305024)

TN

A

孫方 (1982-),女,山東人,高級工程師．現在中國電波傳播研究所從事電波傳播、電波折射修正等方面的工作．

朱慶林 (1981-),男,河南人,博士,高級工程師．現在中國電波傳播研究所從事大氣遙感、GNSS技術以及環境監測等工作．

劉琨 (1979-),男,山東人,博士生,高級工程師．現在中國電波傳播研究所從事GPS測量、電離層探測與反演以及電波折射修正等工作．

趙振維 (1965-),男,河北人,博士,研究員,中國電子科技集團公司首席專家,中國電波傳播研究所總工程師．現在中國電波傳播研究所從事電波傳播、無線電氣象、大氣環境遙感等電波環境及其傳播特性的理論和應用基礎研究．

康士峰 (1966-),男,河南人,博士,研究員,中國電波傳播研究所副總工程師．現在中國電波傳播研究所從事電磁散射理論、電波環境機理與預報、微波遙感等方面的工作．

董翔 (1987-),男,寧夏人,碩士,現在中國電波傳播研究所從事GNSS散射信號傳播、大氣波導監測等方面的工作.

聯系人：孫方 E-mail：sf427@163.com