段寶林

摘? 要：在地震地殼形變監(jiān)測(cè)工作中，InSAR測(cè)量技術(shù)發(fā)揮出了重要作用。其測(cè)量結(jié)果可以廣泛應(yīng)用于研究地震動(dòng)力學(xué)機(jī)制及地震斷層幾何學(xué)特征方面，為進(jìn)一步認(rèn)識(shí)到斷層活動(dòng)、巖石圈動(dòng)力規(guī)律等奠定下了堅(jiān)實(shí)的基礎(chǔ)。本文主要分析了InSAR測(cè)量在地震地殼形變研究中的關(guān)鍵技術(shù)，使InSAR技術(shù)能夠更加合理的應(yīng)用到地震地殼形變檢測(cè)工作中。

關(guān)鍵詞：地震地殼形變? InSAR? 關(guān)鍵技術(shù)? 時(shí)序分析

中圖分類號(hào)：P237? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A? ? ? ? ? ? ? ? ? ? 文章編號(hào)：1674-098X（2020）12（b）-0045-03

Abstract： InSAR has played an important role in crustal deformation monitoring. The measurement results can be widely used in the study of seismic dynamic mechanism and seismic fault geometry characteristics， and lay a solid foundation for further understanding of fault activity and lithosphere dynamic law. This paper mainly analyzes the key technologies of InSAR measurement in the study of seismic crustal deformation， so that InSAR technology can be more reasonable applied to the detection of seismic crustal deformation.

Key Words： Earthquake crustal deformation; InSAR; Key technology; Time series analysis

InSAR是一種用以測(cè)量地面高度及地表形變特征的技術(shù)，具有極高的測(cè)量精度。InSAR時(shí)序分析方法選取了某些SAR數(shù)據(jù)干涉對(duì)組合，通過建模、DEM誤差等，能夠觀察出地面形變?cè)陂L時(shí)間內(nèi)的累積情況。InSAR技術(shù)因使用方法簡單、作業(yè)效率極高，可適應(yīng)極端環(huán)境下的測(cè)量任務(wù)，通過將InSAR所獲取到的地面形變數(shù)據(jù)與相關(guān)資料相結(jié)合，可以實(shí)現(xiàn)更為精確的地震預(yù)報(bào)，對(duì)于當(dāng)今的地震原理研究有著深刻影響。

1? 地震地殼形變InSAR測(cè)量中的關(guān)鍵技術(shù)

1.1 大氣校正及其對(duì)時(shí)序分析流程改進(jìn)

在應(yīng)用InSAR測(cè)量地震地殼形變時(shí)，往往會(huì)因較強(qiáng)的大氣效應(yīng)影響到測(cè)量結(jié)果的精確度。距離震中區(qū)域的形變量較大，干涉條紋不易受到大氣影響，使得各項(xiàng)數(shù)據(jù)較為穩(wěn)定，因此往往會(huì)在此種狀況掩蓋效應(yīng)下，忽視了大氣效應(yīng)的干擾。但在研究中，可以明顯觀察到大氣效應(yīng)所產(chǎn)生的影響，今后需在此方面加強(qiáng)相應(yīng)的研究力度。

InSAR時(shí)序分離出的大氣相位與GNSS數(shù)據(jù)等不盡相同，使得在判定大氣狀況是不能僅依靠此類方法。在使用InSAR時(shí)序分析方法用以分析大氣狀況較為復(fù)雜且地表波動(dòng)較大的區(qū)域時(shí)，因受到傳統(tǒng)時(shí)序分析方法較為狹隘的影響，使得大氣狀況判斷出現(xiàn)偏差，形變測(cè)量誤差較大。造成此現(xiàn)象的因素主要有如下兩方面。

首先，由于大氣的周期性運(yùn)動(dòng)，使得大氣相位不再具有“時(shí)域高頻”的相關(guān)特征;在大氣較厚并且分布不均時(shí)，往往會(huì)導(dǎo)致時(shí)空濾波處理無法進(jìn)行。現(xiàn)代研究表明，在大氣垂直與DEM較為密切時(shí)，其時(shí)空濾波雖已經(jīng)過相應(yīng)處理，但會(huì)殘留一定的大氣相位，極易對(duì)監(jiān)測(cè)造成干擾。其次，InSAR時(shí)序分析需借助于高相干點(diǎn)獲取各類信息，大氣相位分離主要依靠PS點(diǎn)，在分析范圍擴(kuò)大的情況下，使得此方式的精準(zhǔn)度無法保障。

InSAR時(shí)序分析獲取相位信息主要依靠高相干點(diǎn)，PS點(diǎn)具有的密度質(zhì)量、質(zhì)量可保障大氣相位分離的精準(zhǔn)度。在測(cè)量范圍擴(kuò)大的情況下，往往會(huì)受到相干性的影響，導(dǎo)致自身精準(zhǔn)度下降。如僅采取濾波的方式對(duì)大氣相位進(jìn)行處理具有較大的局限性。這主要表現(xiàn)在以下兩方面。

首先，會(huì)導(dǎo)致出現(xiàn)誤判的情況。在進(jìn)行濾波處理時(shí)，由于大氣相位處理不徹底，將會(huì)產(chǎn)生誤判，將某些非形變問題作為地表形變處理。另外如處理過多，使得某些關(guān)鍵信息丟失，無法完整解讀地質(zhì)運(yùn)動(dòng)的基本情況。其次，因大氣條紋與軌道條紋具有重合性特征，降低了基線的重估精準(zhǔn)度，在流程處理環(huán)節(jié)出現(xiàn)較為嚴(yán)重的誤差，無法獲取到精確的結(jié)果。為進(jìn)一步提升InSAR時(shí)序分析形變測(cè)量應(yīng)對(duì)多種環(huán)境的能力，應(yīng)將大氣校正作為關(guān)鍵流程予以重視。

由于大氣校正具有諸多優(yōu)勢(shì)，使其在InSAR時(shí)序分析形變測(cè)量中的應(yīng)用愈發(fā)廣泛。隨著技術(shù)的發(fā)展，現(xiàn)今已出現(xiàn)了采取模擬大氣相位的方法優(yōu)化PS-InSAR，從而進(jìn)一步強(qiáng)化了大氣相位估計(jì)的準(zhǔn)確性，使形變測(cè)量更加合理。在PS-InSAR中加入ECMWF大氣數(shù)值模式后，將地形相位去除，利用ECMWF模擬大氣垂直層流分量，從干涉圖中將其去除，然后采取時(shí)空濾波去除水平湍流大氣相位。除上述方法外，還可以在InSAR時(shí)序分析中采用基于GNSS數(shù)據(jù)的大氣校正的方法，能夠顯著提升監(jiān)測(cè)能力。在此過程中，應(yīng)借助于專業(yè)化軟件，使大氣校正法更為有效。GACOS處理軟件是使用效果較好的軟件之一，其采用了ECMWF大氣產(chǎn)品及GNSS數(shù)據(jù)，可快速生成高精度大氣對(duì)流延遲相位，實(shí)現(xiàn)InSAR大氣校正目的。GACOS因覆蓋范圍極為廣泛、不受時(shí)間及空間的限制、生成SAR數(shù)據(jù)及時(shí)、性能評(píng)價(jià)便捷等優(yōu)勢(shì)，在業(yè)內(nèi)有著極高的應(yīng)用前景。

1.2 三維形變信息提取

在實(shí)際應(yīng)用中，InSAR可與GNSS、Offset Tracking等內(nèi)容相互配合，從而推導(dǎo)出地表三維形變信息。但由于其精確度較低，往往不能適用于精度要求較高的場景中。造成此種狀況的原因在于以下3方面。

（1）Offset Tracking等技術(shù)自身精確度難以滿足實(shí)際使用要求;（2）GNSS的數(shù)據(jù)相對(duì)較少，不具備整體覆蓋的能力;（3）SAR系統(tǒng)的主要工作內(nèi)容不是地表三維形變信息的獲取。從上述內(nèi)容中可見，三維形變信息提取的關(guān)鍵因素在于獲取到足夠且精確的數(shù)據(jù)。為實(shí)現(xiàn)此目標(biāo)，要加大大軌道傾角SAR衛(wèi)星系統(tǒng)的研究力度，并且將獲取地表三維形變信息作為主要功能，扭轉(zhuǎn)SAR衛(wèi)星系統(tǒng)對(duì)三維形變信息提取工作不力的局面。

1.3 大范圍形變監(jiān)測(cè)

（1）ScanSAR時(shí)序分析。

在SAR系統(tǒng)中，ScanSAR的應(yīng)用最為廣泛。雖然ScanSAR有著極強(qiáng)的數(shù)據(jù)覆蓋能力，但由于數(shù)據(jù)干涉效果不佳，加之其自身分辨率不能得到保障，使得其應(yīng)用具有一定的局限性。產(chǎn)生此問題的根源在于，星載SAR系統(tǒng)在設(shè)計(jì)之初即存在Burst不能同步的現(xiàn)象，造成數(shù)據(jù)干涉相干性不強(qiáng)。另外由于ScanSAR模式分辨率較低，使高相干目標(biāo)明顯不足。

（2）InSAR形變結(jié)果拼接。

在地殼形變監(jiān)測(cè)中，因觀測(cè)區(qū)域的差異，需對(duì)干涉測(cè)量結(jié)果進(jìn)行拼接處理?？赏ㄟ^以下策略開展此項(xiàng)工作：首先，在InSAR處理時(shí)需采取參考點(diǎn)才能夠完成，因此應(yīng)考慮參考點(diǎn)對(duì)監(jiān)測(cè)結(jié)果的影響，制定出選擇參考點(diǎn)的標(biāo)準(zhǔn)。其次，需厘清InSAR處理參考點(diǎn)的關(guān)系，通過建立起傳遞理論模型，使其形成統(tǒng)一性。由于空間坐標(biāo)系存在偏差，當(dāng)SAR影像重疊部分的高程差不明顯時(shí)，可將其中一幅SAR影像作為主影像，采取多項(xiàng)式法統(tǒng)一其它SAR影像。如高程差較為顯著，可將待拼接的SAR影像放置于同一坐標(biāo)系中。

1.4 高相干目標(biāo)選取

在InSAR時(shí)序分析中，通常采取多幅SAR影像用來獲取高相干目標(biāo)的形變信息，保證形變結(jié)果的準(zhǔn)確性。高相干目標(biāo)在大氣相位分離、相位解纏精度方面有著較大作用，其數(shù)量直接影響到精確度。

InSAR時(shí)序分析受到多種因素的共同作用，要重點(diǎn)關(guān)注高相干點(diǎn)的數(shù)量級(jí)密度問題。如在大范圍區(qū)域且缺少足夠植被的情況下，則在選取高相干點(diǎn)時(shí)出現(xiàn)較多難點(diǎn)。由于當(dāng)前InSAR時(shí)序分析在選取高相干點(diǎn)時(shí)采取了統(tǒng)一的參數(shù)閾值方法，如將其應(yīng)用于較大范圍時(shí)，極易出現(xiàn)PS點(diǎn)分布不均的現(xiàn)象，使得某些區(qū)域內(nèi)的計(jì)算結(jié)果失真。現(xiàn)今出現(xiàn)了SqueeSAR技術(shù)及StaMPS方法。經(jīng)長期觀察可見，上述兩種方法在地表形變監(jiān)測(cè)有著較好的應(yīng)用空間。其中，SqueeSAR技術(shù)采用了雙樣本KS檢驗(yàn)法，選擇受時(shí)間失相干影響較小的統(tǒng)計(jì)像元，使得高相干目標(biāo)的質(zhì)量與密度均得到了顯著提升。StaMPS方法將PS按照散射體進(jìn)行對(duì)待，以相位穩(wěn)定性作為重要指標(biāo)，在提取高相干目標(biāo)時(shí)，利用了相位空間與振幅離差之間的關(guān)聯(lián)性。

某些區(qū)域中的自然地表相干性條件較差，使得穩(wěn)定高相干像元數(shù)量不足。針對(duì)此種狀況，應(yīng)采用InSAR時(shí)序分析選點(diǎn)的自適應(yīng)法，及時(shí)調(diào)整自身閾值，達(dá)到穩(wěn)定識(shí)別高相干像元的目的。要制定出相應(yīng)的策略，依據(jù)所處環(huán)境，將閾值參數(shù)設(shè)定在合理范圍內(nèi)。要及時(shí)評(píng)價(jià)結(jié)果的差異性，使形變測(cè)量范圍得以增大，最后還應(yīng)繪制出相應(yīng)的分布圖，使結(jié)果更為形象直觀。

1.5 大數(shù)據(jù)背景下數(shù)據(jù)處理方法

隨著信息科技、大數(shù)據(jù)技術(shù)、衛(wèi)星發(fā)射等行業(yè)取得巨大進(jìn)步，在地震地殼形變監(jiān)測(cè)中應(yīng)及時(shí)引入多種方法，充分利用各類檔案數(shù)據(jù)，使數(shù)據(jù)處理更為高效、準(zhǔn)確，不斷提升應(yīng)急響應(yīng)的效率。由于地震地殼形變監(jiān)測(cè)過程中會(huì)產(chǎn)生大量的數(shù)據(jù)資料，為確保觀測(cè)過程的連貫性，其InSAR時(shí)序分析數(shù)據(jù)處理方式需加以改進(jìn)。針對(duì)新獲取到的數(shù)據(jù)資料，InSAR時(shí)序分析數(shù)法可將其及時(shí)加入到數(shù)據(jù)庫中，而無需從初始端進(jìn)行操作。另外由于數(shù)據(jù)累積量及監(jiān)測(cè)任務(wù)逐漸增大，會(huì)影響到計(jì)算機(jī)的運(yùn)行速度，從長遠(yuǎn)角度出發(fā)，應(yīng)及時(shí)將繁多的數(shù)據(jù)遷移到云端，提升計(jì)算機(jī)的效能。

在星載SAR技術(shù)日益成熟的背景下，加之衛(wèi)星重訪周期的時(shí)間越來越短，使得SAR的空間分辨率得以大幅提升。基于衛(wèi)星重訪周期時(shí)間被縮短的狀況，從而有機(jī)會(huì)能夠在極短的時(shí)間內(nèi)獲取到足夠多的資料，但大量的數(shù)據(jù)使得計(jì)算機(jī)處理速度緩慢。在InSAR時(shí)序分析方法中，往往僅依靠對(duì)已存儲(chǔ)的信息進(jìn)行處理，由于地殼形變監(jiān)測(cè)具有顯著的連貫性，使得其數(shù)據(jù)處理滯后性嚴(yán)重，無法適應(yīng)現(xiàn)今對(duì)數(shù)據(jù)資料的大量需求。在今后的工作中，應(yīng)重點(diǎn)研究利用較短重訪周期下SAR數(shù)據(jù)的方法，將其運(yùn)用到InSAR時(shí)序分析中，從而在地震地殼形變監(jiān)測(cè)中發(fā)揮出重要作用，實(shí)現(xiàn)實(shí)時(shí)動(dòng)態(tài)式的監(jiān)測(cè)。正是由于存在此類問題，因此需在InSAR時(shí)序分析技術(shù)采取漸進(jìn)方法，使傳統(tǒng)模式下的數(shù)據(jù)集中處理轉(zhuǎn)變成為增量式處理模式，發(fā)揮出大數(shù)據(jù)技術(shù)的作用。在處理技術(shù)不斷發(fā)展的狀況下，現(xiàn)已出現(xiàn)了漸進(jìn)式SBAS處理方法，可提升Sentine-lA星載SAR干涉數(shù)據(jù)的收集能力，從而獲取到SAR數(shù)據(jù)。此方法能夠依據(jù)現(xiàn)有的形變監(jiān)測(cè)數(shù)據(jù)，推導(dǎo)出當(dāng)前的形變結(jié)果，摒棄了傳統(tǒng)模式下需再次從初始端進(jìn)行操作的弊端，極大提升了時(shí)序分析的準(zhǔn)確性及效率。

2? 結(jié)語

通過上文分析可見，在地殼形變監(jiān)測(cè)中，InSAR技術(shù)在其中發(fā)揮出了重要作用，在地震全過程中均可以監(jiān)測(cè)地殼形變的發(fā)生規(guī)律。我國的空間立體觀測(cè)技術(shù)在衛(wèi)星發(fā)射技術(shù)的引領(lǐng)下，現(xiàn)已取得了較為滿意的成績。尤其隨著SAR測(cè)繪衛(wèi)星和形變監(jiān)測(cè)衛(wèi)星陸續(xù)投入使用，使得InSAR技術(shù)已成為當(dāng)前的熱點(diǎn)研究領(lǐng)域。SAR系統(tǒng)不同于光學(xué)衛(wèi)星，其自身功率有限，不能實(shí)現(xiàn)實(shí)時(shí)且長期的對(duì)地成像。要按照既定計(jì)劃，選取重點(diǎn)觀測(cè)區(qū)域，從而獲得某一區(qū)域內(nèi)充足的數(shù)據(jù)資料，為后續(xù)進(jìn)一步研究提供了數(shù)據(jù)支持。基于實(shí)際應(yīng)用的考慮，將來可在SAR系統(tǒng)中設(shè)置不同任務(wù)的優(yōu)先級(jí)別、成像時(shí)間于運(yùn)行周期之間的時(shí)長關(guān)系等。

參考文獻(xiàn)

[1] 曾琪明，朱猛，焦健.地震地殼形變InSAR測(cè)量中的關(guān)鍵技術(shù)分析[J].遙感學(xué)報(bào)，2018，22（S1）：17-31.

[2] 王慶良.地殼形變觀測(cè)在地震監(jiān)測(cè)預(yù)報(bào)中的發(fā)展與應(yīng)用研究[J].地震研究，2018，41（3）：343.

[3] 董曉燕，趙柯柯，張貴豪.遙感影像匹配技術(shù)在地震形變監(jiān)測(cè)中的應(yīng)用研究[J].黃河水利職業(yè)技術(shù)學(xué)院學(xué)報(bào)，2019，31（2）：31-35.

[4] 鞏丹丹，尚俊斌.地殼形變觀測(cè)在地震監(jiān)測(cè)預(yù)報(bào)中的發(fā)展與應(yīng)用研究[J].甘肅科技縱橫，2019，48（4）：12-14.

[5] Y.G.Wan，Z.-K.Shen，R.Bürgmann，J.B.Sun，M.Wang，靳志同.根據(jù)GPS和InSAR測(cè)量推斷2008年M_W7.9汶川地震的斷層形狀和破裂分布[J].世界地震譯叢，2018，49（3）：195-221.

[6] 李寧，王卓識(shí).地殼形變觀測(cè)在地震監(jiān)測(cè)預(yù)報(bào)中的發(fā)展與應(yīng)用[J].河南科技，2019（34）：138-140.