礦區(qū)沉降監(jiān)測的InSAR技術(shù)應(yīng)用研究分析

在礦區(qū)常規(guī)化管理工作中，地表沉降是較為常見的問題，要想有效減少沉降造成的安全隱患，就要在資源開采過程中配合相應(yīng)的監(jiān)測工序，提升大范圍監(jiān)測水平，InSAR技術(shù)以其全天時、全天候監(jiān)測的優(yōu)勢受到了廣泛應(yīng)用，能進(jìn)一步實(shí)現(xiàn)實(shí)時性動態(tài)監(jiān)測目標(biāo)。

隨著人們生活水平的提高、消費(fèi)支出結(jié)構(gòu)的變化，大眾旅游時代已經(jīng)到來，人們的生活方式與旅行方式也發(fā)生了變化，由于人均收入水平和閑暇時間的保證，我國居民人均出游率由100%上升到150%，旅游逐漸轉(zhuǎn)變?yōu)榫用癯B(tài)化的生活方式。可以預(yù)見，未來五年，旅游需求還將高速增長，到2020年，國內(nèi)旅游的人次將達(dá)到55億-60億人次。

1 InSAR技術(shù)概述

InSAR技術(shù)在應(yīng)用過程中，借助SAR實(shí)現(xiàn)平行軌道上同一個地區(qū)兩幅或者是兩幅以上單視復(fù)數(shù)影像的獲取，并形成干涉（表1），能讀取地區(qū)三維地表信息。近幾年，InSAR技術(shù)和GPS集成處理的模式逐漸得到推廣，能及時推算出電離層和對流層延遲信息，確保參數(shù)分析的合理性，最大程度上保證InSAR技術(shù)不受到外界因素的影響，實(shí)現(xiàn)集成化分析處理。

1.1 原理

另外，利用永久性反射體，配合雷達(dá)影像，就能直接獲取控制點(diǎn)雷達(dá)影像分析內(nèi)容。

1.2 流程

應(yīng)用InSAR技術(shù)完成相應(yīng)工作的過程中，要按照標(biāo)準(zhǔn)化流程進(jìn)行數(shù)據(jù)的處理，從而維持?jǐn)?shù)據(jù)應(yīng)用價值，提升數(shù)據(jù)管理的基本效果。

移取100 μg碲標(biāo)準(zhǔn)溶液于25 mL比色管中，加入4 mg鈣，以水為參比，按實(shí)驗(yàn)方法1.2進(jìn)行試驗(yàn)，結(jié)果見表5。試驗(yàn)結(jié)果表明，沉淀經(jīng)過放置30 min后過濾，能有效消除硫酸鈣的干擾，使吸光度值保持正常，本法選擇沉淀靜置30 min后過濾。

第一，影像配準(zhǔn)。主要是在參考影像和輸入影像之間，要利用技術(shù)處理手段完成同名點(diǎn)的匯總，以此確定相對幾何變換模型，并且配合模型輸入影像完成重新采樣工作，較為常見的方法包括相干系數(shù)法、最大干涉頻譜法、相位差影像平均波動函數(shù)法等

。

第二，干涉圖生成，將精確配準(zhǔn)的復(fù)數(shù)共軛相乘處理后，配合相關(guān)性測度處理方式，就能提升干涉圖像信噪比，發(fā)揮濾波處理的應(yīng)用優(yōu)勢，實(shí)現(xiàn)降噪處理等目標(biāo)。

第一，傳感器ASAR。①主圖像07/2019、輔助圖像08/2019，垂直基線-174m，時間基線35d。②主圖像12/2020、輔助圖像05/2021，垂直基線151m，時間基線140d。③主圖像05/2021、輔助圖像08/2021，垂直基線77m，時間基線105d。

Type: Vietnam. Transplant collected from Son La Province, Van Ho District, Tan Xuan Municipality, Cot Moc Village, at ca. 1000 m a.s.l. 20°40′33.3″N, 104°39′0.3″E, 10 Nov 2015, L. Averyanov, CPC 7158a,b /13279 (holotype, LE, not seen).

1.3 應(yīng)用局限性

盡管InSAR技術(shù)能實(shí)現(xiàn)相關(guān)數(shù)據(jù)信息的分析，但是其應(yīng)用也存在一定的局限性。一方面，針對兩次獲得的成像數(shù)據(jù)可知，隨機(jī)附加噪聲不同或者是不想管，就會造成相位差無法完全抵消的情況，使得信噪比不能滿足實(shí)際數(shù)據(jù)分析處理的應(yīng)用要求和標(biāo)準(zhǔn)，呈現(xiàn)出不明顯的干涉圖，必然會增加數(shù)據(jù)處理的難度。另一方面，大氣相位存在一定的延遲性，主要是因?yàn)榇髿鈼l件存在一定的動態(tài)變化特性，使得相位延遲的現(xiàn)象較為明顯，會出現(xiàn)時間延遲、空間延遲甚至是對有效信息的遮蓋。

綜上所述，在艾滋病抗體檢測中日常檢測可以應(yīng)用酶聯(lián)免疫法以及膠體金法，膠體金法速度更快，但是酶聯(lián)免疫法應(yīng)用價值較高。膠體金法檢測陽性的樣本可以再用酶聯(lián)免疫法復(fù)核。艾滋病傳播與高危險行為存在密切關(guān)聯(lián)，通常自高危險人群開始傳播并在一般人群中流行，因此，必須重視對高危人群加強(qiáng)性病艾滋病防治的健康教育，同時進(jìn)行行為學(xué)和血清學(xué)檢測，檢測過程中減少人為因素對檢出率造成的負(fù)性影響，盡量確保各個實(shí)驗(yàn)步驟的一致性。

因此，在InSAR技術(shù)應(yīng)用過程中，要明確其應(yīng)用要求和標(biāo)準(zhǔn)，踐行更加科學(xué)合理的控制模式，維持分析機(jī)制的規(guī)范性，最大程度上提高監(jiān)測水準(zhǔn)，保證技術(shù)應(yīng)用運(yùn)行的合理性

。

3.2.1 案例

2 InSAR技術(shù)監(jiān)測礦區(qū)沉陷適應(yīng)性分析

結(jié)合數(shù)據(jù)可知，選擇了2018到2021年的干涉對，最長時間基線控制在6個月以下，空間基線也在臨界基線以內(nèi)，相應(yīng)的數(shù)據(jù)具有一定的參考性，能作為后續(xù)進(jìn)行沉降評估的基本依據(jù)。值得一提的是，在外部DEM應(yīng)用的過程中，為了保證測試分析的合理性，對數(shù)據(jù)關(guān)聯(lián)性也進(jìn)行了評估控制。

由圖可知，TerraSAR-X影像在一定條件下時序數(shù)據(jù)會存在失相關(guān)現(xiàn)象。圖中工作面的北側(cè)與正斷層較為接近，南側(cè)為設(shè)計的工作面，西側(cè)與煤礦區(qū)域較為接近，工作面對應(yīng)上覆是采空區(qū)與，回撤通道外側(cè)500m左右就是省級公路區(qū)域。要依據(jù)工作面回采進(jìn)度完成相應(yīng)工作內(nèi)容，利用全部垮落法進(jìn)行采煤工序，且地表位置多數(shù)區(qū)域都是第四系松散沉積物覆蓋，區(qū)域有零星的偽頂區(qū)分布，易形成冒落問題，并且，頂板巖性是泥質(zhì)膠結(jié)材料，為了保證后續(xù)測試的合理性，應(yīng)用InSAR技術(shù)進(jìn)行數(shù)據(jù)分析，選取精準(zhǔn)主影像和輔助影像，對應(yīng)的差分干涉圖時序分析和累積沉降圖分析均符合要求，就能獲取最終的工作面低保沉降時序熱圖，明確時間失相關(guān)對相干性會產(chǎn)生較大的影響，甚至?xí)斐杀O(jiān)測結(jié)果的失效。

裝配式建筑結(jié)構(gòu)對住宅產(chǎn)業(yè)的發(fā)展和推動起到了很大的促進(jìn)作用，這與日本政府的干預(yù)和大力支持是密不可分的：日本政府方面制定了一系列的建設(shè)住宅工業(yè)化方針，并建立了相應(yīng)統(tǒng)一的標(biāo)準(zhǔn)。這樣就解決了標(biāo)準(zhǔn)化與住宅多樣化、復(fù)雜化之間的各種矛盾。通過在不斷的探索和研究中，繼續(xù)進(jìn)一步的快速進(jìn)步。

3 礦區(qū)沉降監(jiān)測的InSAR技術(shù)應(yīng)用內(nèi)容

明確InSAR技術(shù)應(yīng)用要點(diǎn)的基礎(chǔ)上，就要結(jié)合技術(shù)流程和關(guān)鍵落實(shí)相應(yīng)工作，并分析礦區(qū)沉降和開采的關(guān)系，從而確保InSAR技術(shù)能發(fā)揮其優(yōu)勢作用，實(shí)現(xiàn)全面處理的目標(biāo)。

3.1 礦山開采和礦區(qū)地表沉降的關(guān)系

在礦山采集過程中，由于采集工作的頻繁性，使得地下煤礦的開采會造成大范圍的采空區(qū)，并且，采空區(qū)受到上覆巖土層的重力作用，會逐漸出現(xiàn)變形問題，這種變形若是不能得到有效的控制和處理，就會演變?yōu)榈V區(qū)地表沉降現(xiàn)象。在開采工作不斷推進(jìn)的過程中，地表沉降中心也會出現(xiàn)相應(yīng)的移動。通常情況下，地表的最快沉降會發(fā)生在采空區(qū)形成后的3個月到6個月，隨后，相應(yīng)的沉降速度就會趨于穩(wěn)定。為了更好地維持礦區(qū)的安全質(zhì)量，要整合具體參數(shù)，有效維持分析規(guī)范的合理性和可控性。

因?yàn)槔走_(dá)是采取成像的方式完成地表成像結(jié)構(gòu)和數(shù)據(jù)的讀取工作，若是在理想狀態(tài)下，礦區(qū)地表形變在雷達(dá)干涉圖中會以非常清晰的形態(tài)反映出來，以便于技術(shù)部門能對地表沉降情況和地下開采活動的具體狀態(tài)予以實(shí)時性監(jiān)測分析

。

消渴方治療氣陰兩虛型2型糖尿病療效的Meta分析…………………………………………………… 張 蔚等（20）：2851

3.2 具體應(yīng)用

在討論礦區(qū)采空區(qū)沉降分析的過程中，依據(jù)案例的內(nèi)容開展相應(yīng)步驟的分析，明確評估內(nèi)容的要點(diǎn)，從而保證后續(xù)分析結(jié)構(gòu)的規(guī)范性。

綜上所述，不同基因及其產(chǎn)物會對AM發(fā)病具有不同的影響，根據(jù)對不同基因?qū)颊邫C(jī)體中所產(chǎn)生的作用以及途徑進(jìn)行研究，并對有效的阻斷方式進(jìn)行尋找，對治療AM疾病具有重要的影響。通過對不同基因在活化時的相關(guān)信號轉(zhuǎn)導(dǎo)系統(tǒng)的利用，可以有效實(shí)現(xiàn)對異常信號進(jìn)行抑制，并研究出相應(yīng)的信號轉(zhuǎn)導(dǎo)系統(tǒng)的具體方法，并將其與AM治療之間的關(guān)系進(jìn)行研究，以此將其作為治療AM的研究重點(diǎn)。

本文以河北FF礦山為例，其本身是老礦區(qū)，雷達(dá)圖像覆蓋了FF礦區(qū)的背部和中部地區(qū)，礦區(qū)中央有太行山山脈的余脈，海拔高度約為800m，最低位置僅為100m。在對其進(jìn)行InSAR技術(shù)處理分析的基礎(chǔ)上，獲取了覆蓋峰礦區(qū)的基礎(chǔ)數(shù)據(jù)。具體數(shù)據(jù)如下：

第三，平地效應(yīng)消除。要想有效獲取形變信息，就要結(jié)合趨勢面和地形起伏的相位參數(shù)，完成一次差分相位去除處理。主要采取的是相位解纏技術(shù)，結(jié)合路徑控制積分法在滿足最小二乘整體要求的同時完成算法求解。

第二，傳感器JERSI。①主圖像03/2018、輔助圖像04/2018，垂直基線1173m，時間基線44d。②主圖像12/2020、輔助圖像02/2021，垂直基線-1170m，時間基線68d。③主圖像02/2021、輔助圖像03/2021，垂直基線-325m，時間基線40d

。

正是因?yàn)槊旱V區(qū)地表?xiàng)l件較為復(fù)雜，塌陷區(qū)、積水區(qū)的產(chǎn)生，建筑物的破壞等等問題都會影響沉陷分析工作的進(jìn)程。因此，在應(yīng)用InSAR技術(shù)的過程中，要結(jié)合形變區(qū)域沉降信息完成對應(yīng)工作。并且，契合SAR數(shù)據(jù)成像分析工作內(nèi)容，最大程度上滿足干涉條紋評估的要求。比如，干涉相位時空相干對InSAR技術(shù)的應(yīng)用和結(jié)果會產(chǎn)生影響，影像對應(yīng)像元的相關(guān)性、高程精密度等都會形成不同的作用。

3.2.2 數(shù)據(jù)處理

教學(xué)目標(biāo)以傳授知識和技能向運(yùn)用已有技能和知識轉(zhuǎn)變；教學(xué)形式以教師教為主、學(xué)生被動學(xué)習(xí)向以學(xué)生為主體、在老師指導(dǎo)下主動學(xué)習(xí)轉(zhuǎn)變；交流方式以單方面、既有互動也有被動向雙向轉(zhuǎn)變；參與程度以學(xué)生要我學(xué)向我要學(xué)轉(zhuǎn)變；運(yùn)用激勵手段挖掘?qū)W生內(nèi)在潛力，提高學(xué)生的動手能力、合作精神和學(xué)習(xí)效果。

在完成數(shù)據(jù)羅列控制后，要結(jié)合數(shù)據(jù)情況進(jìn)行集中處理分析，從而維持?jǐn)?shù)據(jù)應(yīng)用的合理性。圖像大小選擇滿足圖像應(yīng)用標(biāo)準(zhǔn)。區(qū)域中心位置的經(jīng)緯度是東經(jīng)114.25°、北緯35.56°，整個數(shù)據(jù)規(guī)劃和圖像處理的區(qū)域內(nèi)，基礎(chǔ)尺寸內(nèi)容設(shè)定為20km＊30km，20km為方位向、30km為距離向

。

所以，在InSAR技術(shù)處理過程中，保證SAR數(shù)據(jù)充足、相關(guān)大氣噪聲處理方式合理，才能更好地滿足礦區(qū)沉降監(jiān)測的基本需求，匹配InSAR技術(shù)完成具體工序，提高技術(shù)的適應(yīng)性和可行性

。

與此同時，針對ENVISAT數(shù)據(jù)還進(jìn)行了DORIS處理，能對數(shù)據(jù)予以精準(zhǔn)軌道的分析和調(diào)控，確保數(shù)據(jù)處理后估算和應(yīng)用效果更加貼合實(shí)際需求。并且，能減少生成形成干涉相位圖像的噪音，保證形變干涉相位圖的效果，發(fā)揮濾波處理的有效價值，最終形成形變干涉相位示意圖。

3.3 結(jié)果分析

在完成數(shù)據(jù)處理工作后，就要對獲取的數(shù)據(jù)顯示的結(jié)果進(jìn)行統(tǒng)籌管理，從而獲取相匹配的依據(jù)和內(nèi)容

。

（1）相干特性。

對于各項(xiàng)績效考核指標(biāo)可以采取評分制度，根據(jù)工作人員各項(xiàng)指標(biāo)的表現(xiàn)進(jìn)行打分，滿分為100分。如果沒有完成考核或者出現(xiàn)違反有關(guān)規(guī)章的行為直接扣分。

在干涉對分析的基礎(chǔ)上，結(jié)合ENVISAT形變干涉相位圖，進(jìn)行相應(yīng)數(shù)據(jù)的分析和處理，結(jié)合數(shù)據(jù)顯示的結(jié)果可知，時間基線較短的干涉對其實(shí)際的相干擾效果較好，并且，干涉對只有一個衛(wèi)星回訪周期，相較于其他數(shù)據(jù)結(jié)果，其相干性能更加突出，從而維持分析的實(shí)效性。

第一，借助不同波段雷達(dá)圖像干涉對時間基線的分析可知，ENVISAT相位圖的相干性相較于JERSI要低很多，究其原因，主要是因?yàn)镋NVISAT借助C波段完成雷達(dá)成像處理的同時，其實(shí)際應(yīng)用效果也會受到周圍環(huán)境的影響，大氣環(huán)境、郊區(qū)地表目標(biāo)等因素，都會對數(shù)據(jù)最終結(jié)果的呈現(xiàn)過程形成不同程度的作用。并且，造成C波段去相干作用遠(yuǎn)遠(yuǎn)高于L波段，使得C波段形變干涉相位圖的相干性就要低于L波段

。

廳堂石，指天然形成的，石質(zhì)高雅、色澤細(xì)潤、造型獨(dú)特、個體適中、具有較高品位的觀賞石。如龍化石、貴州紅、晶體石等。

第二，結(jié)合ENVISAT具體結(jié)構(gòu)和形變相干相關(guān)數(shù)據(jù)關(guān)系可知，相干性也會受到時相的影響，盡管形變干涉之間空間基線較短，卻存在兩景雷達(dá)圖像獲取時間差異的問題。比如說，作物的生長以及收割月份，地表作物產(chǎn)生了很大的變化，此時，就會對相干性數(shù)據(jù)產(chǎn)生影響，也就是說，在地表變化大的月份，干涉對會呈現(xiàn)出模糊的條紋，而在秋冬季，地表變化較少，甚至基本無變化，此時，干涉對就會呈現(xiàn)出非常清晰的相變干涉條紋。

第三，結(jié)合JERSI形變干涉圖的具體結(jié)果可知，空間基線也會影響相干性，并且，幾何結(jié)構(gòu)、時間去相干作用都非常明顯，若是無法建立完整且統(tǒng)一的空間基線處理模式，就會制約最終數(shù)據(jù)結(jié)果的分析效率

。

（2）形變干涉相位。

利用InSAR技術(shù)進(jìn)行礦區(qū)地表沉降的測試時，也要關(guān)注形變干涉相位參數(shù)，其有效反映出地表的實(shí)時性動態(tài)變化，以便于能更好地了解相位結(jié)果。

第一，ENVISAT雷達(dá)形變干涉相位對地表形變的反映效果，要遠(yuǎn)遠(yuǎn)優(yōu)于JESRI的效果，其靈敏性更好，能打造更加科學(xué)且直觀的數(shù)據(jù)分析平臺。在相同采空區(qū)內(nèi)，會形成同等大小的地表沉降結(jié)構(gòu)，并且，ENVISAT雷達(dá)干涉也會對可觀測沉降范圍予以直觀的展示和處理，這就能更便于后續(xù)分析。也就是說，ENVISAT雷達(dá)干涉圖能有效劃分可觀測沉降范圍，配合沉降區(qū)域的評估，就能盡可能縮小區(qū)域，以保證后續(xù)操作內(nèi)容都能圍繞標(biāo)志區(qū)域開展。

第二，ENVISAT雷達(dá)干涉在相位圖分析的過程中，要結(jié)合沉降區(qū)域表現(xiàn)形式和狀態(tài)進(jìn)行評估，依據(jù)ENVISAT在沉降區(qū)域相對多的條紋數(shù)目就能更好地完成分析和對比。與此同時，ENVISAT雷達(dá)干涉相位圖還能反映出較大的形變范圍，盡管中心相干性較低，干涉條紋也不明顯，但是卻能維持沉降區(qū)域中心相干性的可控效果。主要是因?yàn)槌上衿骷乇沓两底鳂I(yè)使得地物存在一定的相干失效過程，相較于C波段，L波段能更好地展現(xiàn)出地表變化的范圍，并且將其和實(shí)際參數(shù)予以對比，有效提升相干特性的評估效果

。

（3）沉降監(jiān)測水準(zhǔn)。

改革開放 40 年 浙江自然資源領(lǐng)域改革中那些重要的“印記”（路雄英） .................................................11-14

在地下煤礦開采活動實(shí)時性監(jiān)督管理工作落實(shí)的過程中，為了保證后續(xù)相關(guān)決策體系的合理性，就要維持?jǐn)?shù)據(jù)應(yīng)用的規(guī)范水平，保證開采活動調(diào)研的最優(yōu)化。地表變化會直接反映在雷達(dá)形變干涉相位圖層面，因此，ENVISAT雷達(dá)形變干涉相位圖會借助綠色區(qū)域?qū)Τ两当容^明顯的區(qū)域進(jìn)行劃分和處理，此時，能讀取直觀干涉相位圖。而此時，依據(jù)干涉相位圖還能對礦區(qū)地表變化予以分析，從而有效推斷地下煤礦的開采情況和效果。比如，干涉對“主圖像12/2020、輔助圖像02/2021，垂直基線-1170m，時間基線68d”和“主圖像02/2021、輔助圖像03/2021，垂直基線-325m，時間基線40d”對應(yīng)位置就會形成變相位環(huán)，就能評估出2021年以后，區(qū)域內(nèi)存在明顯的采空區(qū)，并且，形變相位環(huán)說明采空區(qū)地表形變較小，無法被有效觀測，而在地表數(shù)據(jù)出現(xiàn)變動后，證明區(qū)域在形變干涉處理顯示出更加明確的地下采礦活動，這就能有效評估區(qū)域性地表沉降問題的發(fā)生，為全面提升地表沉降處理和后續(xù)決策控制工作水準(zhǔn)奠定基礎(chǔ)

。

另外，依據(jù)基礎(chǔ)干涉對的研究數(shù)據(jù)和對比結(jié)果，繼續(xù)進(jìn)行開采活動的評估和分析，以此為對象，完整分析不同區(qū)域不同時間段地表沉降現(xiàn)象的變化，就能逐漸推斷出變化的趨勢，以便于能減少不安全因素造成的隱患問題，實(shí)現(xiàn)經(jīng)濟(jì)效益和安全效益的和諧統(tǒng)一。

（4）定量分析。

在干涉對評估和數(shù)據(jù)處理的過程中，依據(jù)相關(guān)數(shù)據(jù)進(jìn)行地表沉降現(xiàn)象的分析，也要建立礦區(qū)地表沉降的定量估計。結(jié)合干涉對數(shù)據(jù)可知，視向最大沉降距離逐漸增加，配合水準(zhǔn)測量工序，就能最大程度上完成地表沉降數(shù)值分析的趨勢圖，在沉降量數(shù)量級評估的同時開展后續(xù)作業(yè)

。

4 結(jié)語

總而言之，應(yīng)用InSAR技術(shù)完成雷達(dá)形變干涉相位的評估，有效分析礦區(qū)采空區(qū)的沉降情況，實(shí)現(xiàn)地表沉降分析結(jié)果和監(jiān)測工序的處理目標(biāo)，以便于能開展更加可控且規(guī)范的處理機(jī)制，提升礦區(qū)安全管理的綜合水平，實(shí)現(xiàn)采礦行業(yè)可持續(xù)發(fā)展。

[1]陳洋,陶秋香,劉國林,等．InSAR與概率積分法聯(lián)合的礦區(qū)地表沉降精細(xì)化監(jiān)測方法[J]．地球物理學(xué)報,2021,64(10)：3554-3566．

[2]周定義,左小清．基于SBAS-InSAR和PSO-BP神經(jīng)網(wǎng)絡(luò)算法的礦區(qū)地表沉降監(jiān)測及預(yù)測[J]．云南大學(xué)學(xué)報（自然科學(xué)版）,2021,43(5)：895-905．

[3]劉海洋,吳文清,王振海．基于SBAS-InSAR技術(shù)的煤炭壩礦區(qū)地表沉降監(jiān)測與分析[J]．測繪與空間地理信息,2021,44(7)：69-72．

[4]王建,祁鑫博,王麗華,等．SBAS-InSAR技術(shù)在礦區(qū)地表沉降監(jiān)測中的應(yīng)用[J]．北京測繪,2021,35(4)：548-552．

[5]楊光欣,朱煜峰,張龍剛．Stacking-InSAR技術(shù)在礦區(qū)地表沉降監(jiān)測中的應(yīng)用[J]．礦山測量,2021,49(2)：98-102,119．

[6]徐郡,胡晉山,丁進(jìn),等．結(jié)合DInSAR和MAI-InSAR的礦區(qū)地表沉降監(jiān)測[J]．江西理工大學(xué)學(xué)報,2020,41(1)：15-20．

[7]李勇發(fā),左小清,麻源源,等．基于PS-InSAR技術(shù)和遺傳神經(jīng)網(wǎng)絡(luò)算法的礦區(qū)地表沉降監(jiān)測與預(yù)計[J]．地球物理學(xué)進(jìn)展,2020,35(3)：845-851．

[8]鐘曉春．基于InSAR技術(shù)的礦區(qū)沉降監(jiān)測研究進(jìn)展[J]．測繪與空間地理信息,2020,43(2)：178-181．

[9]王志勇,孟慶穎．InSAR相位解纏算法在礦區(qū)沉降監(jiān)測中的對比與分析[J]．遙感信息,2014(6)：47-51,57．

[10]任文靜,賈洪果,閆斌．SBAS-InSAR方法支持下的礦區(qū)地表沉降監(jiān)測及參數(shù)反演[J]．測繪通報,2021(3)：113-117,155．

[11]潘光永,陶秋香,陳洋,等．基于SBAS-InSAR的山東濟(jì)陽礦區(qū)沉降監(jiān)測與分析[J]．中國地質(zhì)災(zāi)害與防治學(xué)報,2020,31(4)：100-106,120．

[12]丁劉建,陶秋香,高騰飛,等．SBAS InSAR技術(shù)在礦區(qū)地面沉降監(jiān)測中的應(yīng)用[J]．中國科技論文,2019,14(3)：320-325．