楊乾坤

（寧夏煤礦設(shè)計(jì)研究院有限責(zé)任公司，寧夏 銀川 750001）

地質(zhì)災(zāi)害中存在一種變化較慢且占地寬廣的災(zāi)害，即為地面沉降。在地面沉降中，可以采用雙軌D-InSAR 技術(shù)對(duì)地面沉降進(jìn)行檢測(cè)[1]。雙軌D-InSAR 技術(shù)發(fā)展至今，理論與數(shù)據(jù)處理技術(shù)已經(jīng)相對(duì)成熟。雙軌D-InSAR 技術(shù)可準(zhǔn)確檢測(cè)出地面沉降情況，減少資源浪費(fèi)，對(duì)采礦事業(yè)的發(fā)展有著重要意義。

1 礦山開(kāi)采區(qū)地面沉降采用雙軌D-InSAR技術(shù)實(shí)測(cè)驗(yàn)證

礦山開(kāi)采區(qū)地面出現(xiàn)的沉降現(xiàn)象與普通的地表形變有較大差異，礦區(qū)地面沉降具有時(shí)域非線性、沉降梯度大的特點(diǎn)。因此，需要選擇適合礦區(qū)地面沉降監(jiān)測(cè)技術(shù)。本文將對(duì)雙軌D-InSAR 技術(shù)在礦區(qū)地表沉降檢測(cè)中的應(yīng)用進(jìn)行研究。

在礦區(qū)地表變化信息與三維信息的獲取上，采用雙軌D-InSAR 技術(shù)，以合成孔徑雷達(dá)數(shù)據(jù)提取的相位信息為信息源[2]。雙軌D-InSAR 干涉圖由于兩副天線與地面某一目標(biāo)之間的距離不等而形成，根據(jù)兩次成像的相位差與地面目標(biāo)的三維空間位置之間存在的幾何關(guān)系。在探測(cè)礦區(qū)地表形變時(shí)，干涉圖中的相位值即為兩次成像的相位差測(cè)量值[3]，可測(cè)定地面目標(biāo)的三維坐標(biāo)，從而提供大范圍的高精度數(shù)字高程模型。雙軌D-InSAR 的處理過(guò)程就是通過(guò)外部數(shù)字高程模型或其它干涉影像對(duì)獲取的相位進(jìn)行差分處理，以便去除無(wú)用相位信息而獲取形變相位的一個(gè)過(guò)程。下圖為地表沉降信息提取流程。

2 基于雙軌D-InSAR技術(shù)的礦區(qū)地面沉降檢測(cè)方法

2.1 利用雙軌D-InSAR 技術(shù)精確定位地面控制點(diǎn)

礦區(qū)地形出現(xiàn)起伏情況下，不同高程點(diǎn)處干涉條紋顯現(xiàn)出的顏色有所差異，以像元值為依據(jù)將礦區(qū)地面變形情況進(jìn)行等級(jí)劃分，能夠清晰看出礦區(qū)地面的沉降情況，所以得出結(jié)論為礦區(qū)地面沉降高程位置的干涉環(huán)的顏色會(huì)發(fā)生變化。地面控制點(diǎn)作為零形變點(diǎn)進(jìn)行反向推算出沉陷區(qū)的形變量，控制點(diǎn)的選取為影像處理非常關(guān)鍵的一個(gè)步驟，控制點(diǎn)選取的質(zhì)量直接影響到地表形變反演的精度，因此在選取過(guò)程中，應(yīng)該避免控制點(diǎn)位于形變相位上。

2.2 模擬礦區(qū)地面沉降差分干涉相位

地面差分干涉相位值與雙軌D-InSAR 最大沉降量相同，相位模擬即對(duì)某個(gè)波段不同點(diǎn)號(hào)最大沉降量進(jìn)行求導(dǎo)，將數(shù)據(jù)匯總得到差分干涉相位圖與相應(yīng)數(shù)據(jù)。不同類型雷達(dá)干涉數(shù)據(jù)對(duì)地面沉降監(jiān)測(cè)的能力因雙軌D-InSAR 檢測(cè)的最大沉降量而決定，考慮礦區(qū)實(shí)際沉降、入射角、波長(zhǎng)，以及地面分辨率等因素，對(duì)礦區(qū)地面沉降的差分干涉相位進(jìn)行模擬。

由于L和C波段SAR 干涉數(shù)據(jù)的波長(zhǎng)、入射角和地面分辨率不同，因此相位入侵角度會(huì)發(fā)生一定的改變，對(duì)于同一礦區(qū)地面沉降而言，其在差分干涉圖距離向上，會(huì)均勻分布相應(yīng)的相位元數(shù)以及相位方向，當(dāng)轉(zhuǎn)換對(duì)應(yīng)的視線方向上，便可識(shí)別出最大沉降量相位以及不同數(shù)據(jù)下沉降相位。

圖1 礦區(qū)地表沉降信息提取流程

圖2 兩種方法沉降量對(duì)比圖

2.3 分析地面沉降量數(shù)據(jù)

在礦區(qū)地面沉降檢測(cè)過(guò)程中，模擬數(shù)據(jù)的變化是監(jiān)測(cè)的重點(diǎn)。根據(jù)最大沉降量相位以及不同數(shù)據(jù)下沉降相位，對(duì)地面沉降量數(shù)據(jù)進(jìn)行分析。設(shè)計(jì)的波段點(diǎn)號(hào)雷達(dá)入射角度為45°、雷達(dá)入射波波長(zhǎng)為130mm；規(guī)定的下沉盆地主要影響半徑為1530m，SAR 傳感器的地面分辨率為60%。查閱多方資料得到十個(gè)點(diǎn)號(hào)的十組數(shù)據(jù)，并對(duì)數(shù)據(jù)進(jìn)行參數(shù)矯正后，匯總得出它們的最大沉降量。令最大沉降量與礦區(qū)地面沉降參量作對(duì)比，得出沉降相位值如表1 所示。

表1 最大沉降量與礦區(qū)地面沉降參量分析表

通過(guò)以上方法步驟，完成了基于雙軌D-InSAR 技術(shù)對(duì)礦區(qū)地面沉降量的檢測(cè)。

3 對(duì)比實(shí)驗(yàn)

為研究基于雙軌D-InSAR 技術(shù)的礦區(qū)地面沉降檢測(cè)方法與傳統(tǒng)方法的區(qū)別，實(shí)驗(yàn)選取地面沉降情況較嚴(yán)峻的某地作為實(shí)驗(yàn)環(huán)境，礦區(qū)地屬盆地，地勢(shì)東北略高、西南稍低。實(shí)驗(yàn)環(huán)境基本數(shù)據(jù)為：礦井范圍南北長(zhǎng)3km～3.9km，東西寬3.5km～5.5km；地理坐標(biāo)120°31'～120°35'E，37°31'～339°35'N。然后采用本文方法與傳統(tǒng)方法分別對(duì)礦區(qū)地面沉降量與時(shí)間進(jìn)行趨勢(shì)匯總，在同樣時(shí)間內(nèi)，兩種方法下沉降量與實(shí)際沉降量的關(guān)系如圖2 所示。

圖2 中，a 圖顯示的是傳統(tǒng)方法下礦區(qū)地面沉降量變化情況，b 圖為基于雙軌D-InSAR 技術(shù)的礦區(qū)地面沉降檢測(cè)方法沉降量變化圖。根據(jù)圖中變化趨勢(shì)圖能夠發(fā)現(xiàn)：在相同的60s 時(shí)間內(nèi)，a 圖沉降量總體呈上升趨勢(shì)，b 圖雖然也為上升趨勢(shì)，但其上升坡度較小，相對(duì)圖a 在每秒鐘內(nèi)沉降量上升較慢；在60s 時(shí)，a 圖沉降量為289mm，b 圖沉降量為196mm，在同樣的上升趨勢(shì)下，b 圖沉降量較a 圖小了93mm。a 圖中沉降量變化更貼近實(shí)際沉降量，b 圖中沉降量與實(shí)際沉降量偏差較大。得出基于雙軌D-InSAR 技術(shù)的礦區(qū)地面沉降檢測(cè)方法檢測(cè)出的沉降量更貼近實(shí)際沉降量，使檢測(cè)結(jié)果更精確。

4 結(jié)語(yǔ)

通過(guò)最大沉降量監(jiān)測(cè)、相位模擬和礦區(qū)地面沉降數(shù)據(jù)分析，明確了L 波段與C 波段雷達(dá)干涉數(shù)據(jù)礦區(qū)地面沉降的監(jiān)測(cè)能力。再由基于雙軌D-InSAR 技術(shù)的礦區(qū)地面沉降檢測(cè)方法與傳統(tǒng)方法作對(duì)比，得出基于雙軌D-InSAR 技術(shù)的礦區(qū)地面沉降檢測(cè)方法可以獲取更接近實(shí)際沉降量的值，應(yīng)用雙軌D-InSAR 技術(shù)能夠得到更加精確的礦區(qū)地面沉降檢測(cè)結(jié)果。