劉 哲

(漢中職業(yè)技術(shù)學(xué)院，陜西 漢中 723002)

0 引言

地面沉降可能會(huì)造成建筑物地基下沉、房屋開(kāi)裂、地下管道受損等危害。近年來(lái)，地面沉降問(wèn)題已成為制約城市發(fā)展的重大地質(zhì)災(zāi)害問(wèn)題之一。目前，我國(guó)華北平原、長(zhǎng)江三角洲以及汾渭地埑區(qū)域等均發(fā)生了不同程度的地面沉降[1]。西安位于渭河盆地，區(qū)內(nèi)構(gòu)造發(fā)育強(qiáng)烈，自20世紀(jì)60年代發(fā)現(xiàn)至今，西安已經(jīng)成為我國(guó)地面沉降最為嚴(yán)重的地區(qū)之一[2]。因此，為了控制西安市地面沉降及地裂縫的發(fā)展，急需開(kāi)展大范圍、高精度的地表形變監(jiān)測(cè)，從而為此項(xiàng)工作提供技術(shù)支撐。

目前，常用的地表形變監(jiān)測(cè)手段有水準(zhǔn)測(cè)量、全站儀測(cè)量、GNSS測(cè)量等，這些常規(guī)的監(jiān)測(cè)手段均為點(diǎn)狀測(cè)量，監(jiān)測(cè)時(shí)間分辨率高，但空間分辨率低，無(wú)法進(jìn)行大面積的監(jiān)測(cè)。InSAR技術(shù)作為近年來(lái)新興的地表形變監(jiān)測(cè)手段，具有全天時(shí)、全天候、大面積、高精度的優(yōu)點(diǎn)，目前已經(jīng)廣泛應(yīng)用到了地面沉降、滑坡、地震等地質(zhì)災(zāi)害監(jiān)測(cè)中[3]。但是，DInSAR技術(shù)容易受時(shí)間失相干、空間失相干以及大氣延遲的影響，近年來(lái)發(fā)展的時(shí)間序列InSAR技術(shù)(PSInSAR[4]和SBAS[5]等)克服了DInSAR技術(shù)的缺點(diǎn)，目前已經(jīng)廣泛應(yīng)用到城市地表形變監(jiān)測(cè)中[6]。例如，彭米米等[7]于2018年利用InSAR技術(shù)監(jiān)測(cè)西安市2015—2017年地面沉降分布特征，并分析地面沉降和地裂縫的相互關(guān)系。

本文利用2020年1月—2021年12月共52景覆蓋西安地區(qū)的Sentinel-1A雷達(dá)衛(wèi)星影像和SBAS-InSAR技術(shù)對(duì)西安市地表形變進(jìn)行了監(jiān)測(cè)，分析近2年來(lái)西安市地表形變規(guī)律，為后續(xù)政府部門研究治理西安市地面沉降及地裂縫提供一定參考依據(jù)。

1 研究區(qū)概況

西安市位于關(guān)中平原的中部(33.42°～34.45°N，107.40°～109.49°E)，北至渭河，南接秦嶺。西安市區(qū)域內(nèi)除河流區(qū)域外均被濕陷性黃土覆蓋，濕陷性黃土在受到水浸潤(rùn)后，土壤結(jié)構(gòu)容易發(fā)生破壞，從而導(dǎo)致地表發(fā)生沉降，甚至形成地裂縫，目前西安市已經(jīng)發(fā)現(xiàn)的地裂縫共有14條，嚴(yán)重制約了西安市的城市發(fā)展[8]?！段靼彩械刭|(zhì)災(zāi)害防治“十四五”規(guī)劃》顯示，全市地質(zhì)災(zāi)害隱患點(diǎn)共分為地質(zhì)災(zāi)害高易發(fā)區(qū)、中易發(fā)區(qū)、低易發(fā)區(qū)和不易發(fā)區(qū)，其中地質(zhì)災(zāi)害高易發(fā)區(qū)占全市總面積的16.18%，發(fā)育有崩塌、滑坡、地裂縫、泥石流等。

2 研究數(shù)據(jù)與方法

2.1 研究數(shù)據(jù)

為了研究西安市的地表形變信息，本次研究采用了Sentinel-1A雷達(dá)衛(wèi)星影像。本次研究區(qū)域地理范圍位于34°04′23″～34°24′26″N，108°37′39″～109°10′49″E，采用了2020年1月—2021年12月共52景覆蓋西安地區(qū)的Sentinel-1A雷達(dá)衛(wèi)星影像。圖像中心入射角34.243°，分辨率為5 m×20 m，SAR影像覆蓋區(qū)域約為50 km×30 km，總面積1 500 km2。此外，還獲取了對(duì)應(yīng)時(shí)間的精密軌道信息以及30 m分辨率的STRM數(shù)據(jù)。

2.2 研究方法

(1)

式中，d(tj,x,r)和d(ti,x,r)分別為相對(duì)于t0時(shí)刻LOS向的累積形變量。為方便理解，式中并未考慮大氣誤差、噪聲誤差以及殘余地形誤差。

因此，M幅干涉圖可根據(jù)式(2)得到：

δφ(x,r)=Aφ(x,r)

(2)

式中，A為(M×N)矩陣，由0,1,-1組成。

當(dāng)所有數(shù)據(jù)位于一個(gè)小基線集時(shí)，此時(shí)M≥N，利用最小二乘方法可以得到:

φ(x,r)=(ATA)-1ATδφ(x,r)

(3)

當(dāng)所有數(shù)據(jù)不在同一個(gè)子集時(shí)，此時(shí)方程秩虧，方程得到無(wú)限組解。因此，為了得到唯一解，Berardino等[5]提出了SVD分解的方法得到最小范數(shù)上的最小二乘解，從而得到了LOS向的累積形變量。

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

首先，選取第一景SAR影像為主影像，將其余的SAR影像與主影像進(jìn)行配準(zhǔn)；其次，設(shè)置時(shí)空基線閾值進(jìn)行構(gòu)網(wǎng)，為了避免大氣對(duì)干涉對(duì)的影響，人工篩選刪除了一些受大氣影響較為嚴(yán)重的干涉圖，最終得到188個(gè)干涉對(duì)進(jìn)行后續(xù)計(jì)算。利用DEM數(shù)據(jù)模擬地形相位進(jìn)行差分干涉，將得到的差分干涉相位進(jìn)行濾波，再選定穩(wěn)定參考點(diǎn)后進(jìn)行相位解纏；由于本次實(shí)驗(yàn)只有一個(gè)基線集，因此可直接通過(guò)最小二乘解算形變序列和平均形變速率，最后將結(jié)果投影到地理坐標(biāo)系下，得到WGS84坐標(biāo)系下的形變結(jié)果。

3 研究結(jié)果分析

3.1 LOS向平均速率結(jié)果

通過(guò)前文所述SBAS方法得到了西安市2020年1月—2021年12月的部分形變序列結(jié)果。通過(guò)形變序列結(jié)果，可以看出西安市近2年整體上形變保持穩(wěn)定，部分區(qū)域出現(xiàn)了沉降或抬升現(xiàn)象，最大累積形變量達(dá)-100～100 mm。根據(jù)形變序列結(jié)果，計(jì)算得到了西安市近2年LOS向平均形變速率結(jié)果：在2020年1月—2021年12月，西安市整體上保持穩(wěn)定狀態(tài)，形變速率在-60～55 mm/a，但局部地區(qū)出現(xiàn)了部分地面沉降及地面回彈現(xiàn)象。

3.2 西安市形變特征點(diǎn)分析

為了更好地分析西安市的形變信息，結(jié)合歷史光學(xué)遙感影像，選擇了具有明顯形變的2個(gè)形變特征點(diǎn)進(jìn)行形變序列分析。

(1)Ⅰ號(hào)形變區(qū)域位于灞河的東北方向，形變區(qū)域面積約1.16 km2，A點(diǎn)位于西安市灞橋區(qū)洪慶街道西侯村與東侯村交界處。特征點(diǎn)A 2020年1月—2021年12月LOS向最大累積形變量達(dá)-48 mm，其中在2020年3月、2021年3月形變速率加快，2021年5月出現(xiàn)地面回彈，最后趨于穩(wěn)定。結(jié)合光學(xué)影像數(shù)據(jù)發(fā)現(xiàn)，2019年7月后，該處區(qū)域城市建設(shè)進(jìn)程加快，修建了大量住宅，導(dǎo)致城市人口增多，地下水需求增大，從而大量抽取地下水，形成了地面沉降。

(2)Ⅱ號(hào)形變區(qū)域位于雁塔區(qū)西南部，形變區(qū)域面積約6.58 km2，B點(diǎn)位于西安市雁塔區(qū)電子城街道唐家村西側(cè)。特征點(diǎn)B 2020年1—12月出現(xiàn)下沉后，2020年12月—2021年12月出現(xiàn)大面積地面回彈，累積形變量達(dá)100 mm。結(jié)合前人研究結(jié)果，發(fā)現(xiàn)該區(qū)域發(fā)生地面回彈的原因主要有政府進(jìn)行城中村改造，同時(shí)關(guān)閉了抽水井，并將處理過(guò)后的地表水注入地下，從而造成地面回彈現(xiàn)象。

4 結(jié)語(yǔ)

本文采用2020年1月—2021年12月共52景覆蓋西安地區(qū)的Sentinel-1A雷達(dá)衛(wèi)星影像對(duì)西安市地表形變進(jìn)行監(jiān)測(cè)，得到了西安市年平均LOS向地表形變速率，通過(guò)結(jié)合歷史谷歌地球影像數(shù)據(jù)，得到以下結(jié)論：

(1)結(jié)合SBAS技術(shù)，本文對(duì)西安市SAR影像覆蓋的約50 km×30 km、總面積1 500 km2區(qū)域進(jìn)行計(jì)算，發(fā)現(xiàn)在2020年1月—2021年12月，西安市整體上保持穩(wěn)定狀態(tài)，但局部地區(qū)出現(xiàn)了地面沉降及回彈現(xiàn)象，其LOS向平均形變速率在-60～55 mm/a。

(2)本次實(shí)驗(yàn)共發(fā)現(xiàn)4個(gè)地表形變區(qū)域，其中Ⅰ區(qū)域?yàn)榈孛娉两祬^(qū)域，Ⅱ,Ⅲ,Ⅳ區(qū)域?yàn)榈乇硖齾^(qū)域。通過(guò)提?、?Ⅱ區(qū)域特征點(diǎn)發(fā)現(xiàn)，Ⅰ，Ⅱ區(qū)域LOS向最大累積形變量分別達(dá)-48 mm，100 mm。

(3)結(jié)合歷史谷歌地球影像和前人文獻(xiàn)分析可得，西安市地表形變主要與地下水開(kāi)采以及人口搬遷有關(guān)，過(guò)量開(kāi)采地下水以及建筑物的附加荷載是造成地面沉降的主要原因。

(4)針對(duì)近2年西安市的InSAR形變結(jié)果，西安市的地面沉降逐漸減弱，但多種因素(承壓水開(kāi)采、工程降水、潛水開(kāi)采等)誘發(fā)地面沉降的風(fēng)險(xiǎn)仍舊存在。因此，大力發(fā)展多源數(shù)據(jù)融合技術(shù)(天-空-地)一體化監(jiān)測(cè)技術(shù)，才能有效防范化解地質(zhì)災(zāi)害的威脅，保障人民群眾的生命財(cái)產(chǎn)安全。