蘇小寧,朱慶,孟國杰,孫海麗,閆浩文
1 蘭州交通大學測繪與地理信息學院,蘭州 730070 2 西南交通大學地球科學與環(huán)境工程學院,成都 611756 3 中國地震局地震預(yù)測研究所,北京 100036 4 首都師范大學資源環(huán)境與旅游學院,北京 100048
GNSS(Global Navigation Satellite System)觀測技術(shù)在研究全球地殼運動和構(gòu)造活動等方面發(fā)揮了重要的作用,“中國地殼運動觀測網(wǎng)絡(luò)”、“中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)”、首都圈GNSS觀測網(wǎng)和省屬地震氣象測繪部門的CORS(Continuously Operating Reference Stations)網(wǎng)等覆蓋全國和區(qū)域的多源觀測數(shù)據(jù),為獲取中國大陸較高空間分辨率的變形場提供了豐富的資料(Wang and Shen,2020).由于GNSS速度場相對于某個參考基準,因此選擇不同參考基準導(dǎo)致速度場的表現(xiàn)形式差異巨大.而應(yīng)變率場與參考基準選擇無關(guān),通過對比應(yīng)變積累和地震矩釋放的關(guān)系,在地震危險性分析和長期預(yù)報方面發(fā)揮著重要作用(Ojo et al., 2021),這其中的關(guān)鍵是解算出能反映測站空間不規(guī)則分布特征且精準的應(yīng)變率場.
目前解算應(yīng)變場的方法可總結(jié)為兩類:一類是建立塊體-斷層模型,計算塊體內(nèi)部均勻應(yīng)變,稱為物理模型(申重陽等,2002;Tape et al., 2009);另一種是利用GNSS速度場,構(gòu)建空間連續(xù)的速度場模型,進而解算應(yīng)變率場,稱為數(shù)學模型(Tape et al., 2009).可見,在數(shù)學模型中關(guān)鍵是如何精準構(gòu)建空間連續(xù)的速度場模型,常用的方法有距離加權(quán)(Shen et al,1996)、多面函數(shù)(武艷強等,2009)、Kriging插值(劉曉霞等,2014)、最小二乘配置(江在森和劉經(jīng)南,2010)和多尺度球面小波(Tape et al., 2009;程鵬飛等,2015;蘇小寧等,2016)等,這些方法在確定觀測值權(quán)重時或?qū)NSS速度場看作為等精度的觀測值,或利用擬合的速度場形式誤差表征其不等精度的特征.但實際上經(jīng)常存在由于觀測數(shù)據(jù)質(zhì)量較差導(dǎo)致的誤差較大或者土層觀測點標墩不穩(wěn)定產(chǎn)生的局部速度差異,也就是在觀測數(shù)據(jù)中存在少量粗差.此時利用現(xiàn)有方法構(gòu)建速度場模型時,少量粗差將會對解算的應(yīng)變場產(chǎn)生重要的影響,即上述方法均不具備抗差功能.
本文基于能夠體現(xiàn)GNSS測站不規(guī)則分布特征的多尺度球面小波方法構(gòu)建GNSS速度場,獲取不同空間尺度的地殼運動場,在此基礎(chǔ)上通過抗差迭代算法顧及速度場中存在的少量粗差對建立速度場模型的影響,通過在實測數(shù)據(jù)中人為加入粗差的模擬計算驗證方法的抗差性,結(jié)合川滇地區(qū)最新的GNSS速度場數(shù)據(jù)給出了區(qū)域的速度場和應(yīng)變率場模型.
球面上觀測點的三維速度可以以南北向、東西向和垂直向3個分量來表示,本文只考慮水平向速度場,因此只涉及到南北向和東西向2個分量,考慮垂直向時的方法與本文所述相同.地面觀測站(λ,φ)的地殼運動速度ν(λ,φ)可表示為
(1)
基于球面小波算法構(gòu)建多尺度地殼運動場模型的本質(zhì)是采用小波基函數(shù)作為估計速度場的框架,利用對小波母函數(shù)進行平移和伸縮等仿射變換后得到的小波基函數(shù)來表示不同分量的速度場,因此公式(1)可以變換為
(2)
(3)
以高斯函數(shù)差建立的球面小波基函數(shù)(Difference of Gaussians,DOG)在構(gòu)建多尺度地殼運動場模型時得到了廣泛應(yīng)用(蘇小寧等,2016),DOG小波的函數(shù)模型為(包伯成等,2009)
(4)
空間球面小波分解的原則是根據(jù)觀測站的分布建立一系列球面小波基函數(shù).通過對球面進行三角形分解,可以獲得在球面上均勻分布的球面三角形網(wǎng)格,分解的層次與球面小波基函數(shù)的尺度相對應(yīng),也對應(yīng)于不同邊長的球面三角形,具體形式可參考 Wang和Dahlen(1995).首先根據(jù)分解尺度建立空間規(guī)則分布的小波基函數(shù),再根據(jù)地表觀測點的空間密度對規(guī)則分布的小波基函數(shù)進行篩選,使得最終的小波基函數(shù)集在空間分布上與地表觀測站的空間密度密切相關(guān),進而構(gòu)建的模型能夠凸顯測站不規(guī)則分布的特征.分解的具體步驟如下:
(1)選取研究區(qū)經(jīng)緯度的覆蓋范圍形成球面矩形,計算球面矩形對角線的長度,搜索小于該長度且距離最近的三角形邊長,確定最小分解尺度;
(2)根據(jù)不同分解尺度,建立多層次且規(guī)則分布的空間小波基函數(shù),標注對應(yīng)的球面三角形邊長;
(3)根據(jù)地表觀測點分布,對所有小波基函數(shù)逐個搜索,以小波基函數(shù)所在位置為圓心和分解尺度對應(yīng)的球面三角形邊長為半徑,確定小波基函數(shù)所對應(yīng)的“空間支撐”范圍,統(tǒng)計“空間支撐”內(nèi)地面觀測站的個數(shù),篩選出“空間支撐”內(nèi)部多于3個站點的小波基函數(shù),獲得模型可用的最大分解尺度、空間不均勻分布的小波基函數(shù)及其數(shù)量;
(4)統(tǒng)計分析不同分解尺度小波基函數(shù)的有效覆蓋范圍,選擇有效覆蓋范圍超過70%的最大分解尺度為最終的球面小波分解尺度,形成在空間不規(guī)則分布且不同分解尺度的小波基函數(shù)集.
根據(jù)地表觀測站的空間分布建立的球面小波基函數(shù)個數(shù)為N個,因此與其對應(yīng)的系數(shù),也就是待求參數(shù)的個數(shù)也是N個,而觀測點的個數(shù)為M個,也就是觀測量個數(shù)為M.在實際建立速度場模型時,觀測點個數(shù)往往小于待求參數(shù)個數(shù),導(dǎo)致觀測方程的參數(shù)求解成為對欠定方程的求解問題.同時,由母函數(shù)經(jīng)過仿射變換得到的DOG小波基函數(shù)為非正交基,導(dǎo)致法方程秩虧,涉及到不適定方程的求解問題.
在建立速度場模型時,可以通過最優(yōu)化算法中的正則化方法,求解出待求參數(shù)向量,目標函數(shù)和參數(shù)的最優(yōu)解分別為
(5)
m=(GTPDG+λ2I)-1GTPDd,
(6)
式中S(m)為目標函數(shù),通過以目標函數(shù)為最小準則獲取參數(shù)的最優(yōu)解.G為由小波基函數(shù)組成的系數(shù)矩陣,m為待求參數(shù)向量,d為地表GNSS觀測點觀測的速度場向量,CD為觀測數(shù)據(jù)的方差-協(xié)方差矩陣,其逆矩陣為權(quán)矩陣PD,考慮到不同位置測站計算的地殼運動速度為獨立觀測量,則CD和PD均為對角矩陣.CM為待求參數(shù)的方差-協(xié)方差矩陣,本文采用Tikhonov正則化方法(徐天河和楊元喜,2003),CM的具體形式為λ2I,I為單位對角陣,λ為正則化參數(shù).利用公式(6),通過多次迭代之后,則可以求解出參數(shù)m的平差值.可以看出,由于對法方程矩陣的對角元素進行了修正,顯然該平差值為有偏估計值.
從公式(6)可以看出,在參數(shù)最優(yōu)解求解的過程中,如果已知基函數(shù)分布和最優(yōu)化解算方法,則參數(shù)的解只依賴于如何確定觀測數(shù)據(jù)的權(quán)矩陣.以往普遍采用擬合的地表測站GNSS速度場的形式誤差給定,但如果在觀測數(shù)據(jù)中存在異常值或者粗差值,利用形式誤差確定權(quán)矩陣的方法并不能精確地表達其隨機特征,此時公式(6)不具備抗差功能,則無法獲得精確的參數(shù)解,最終建立的速度場模型值得商榷,因此,有必要對公式(6)進行修正,使其具備抗差功能.
一個好的估計方法不僅要具有“穩(wěn)健性”,還得具有“抗干擾性”.“穩(wěn)健性”是指當假定模型與實際模型有微小差異時,其估計方法的性能只受到微弱的影響;“抗干擾性”是指當觀測樣本中混入少量粗差時,估計值受到的影響不大(楊元喜等,2002).地表不規(guī)則分布的地殼運動速度場屬于獨立觀測,且不等精度,適合采用按等價權(quán)原理作抗差估計,目的是解決不等精度獨立觀測值權(quán)函數(shù)的確定問題(周江文,1989;楊元喜,1994).
在公式(5)所示最優(yōu)化算法的目標函數(shù)中,第一項為觀測數(shù)據(jù)殘差的加權(quán)平方和,顯然個別異常大殘差的出現(xiàn)會導(dǎo)致平方和顯著增大.為了達到平方和最小的目標函數(shù),平差值必然會遷就那些異常觀測值或者粗差值,使得平差值偏離最優(yōu)值,因此個別異常值或粗差會對整個估值產(chǎn)生較大的影響.抗差估計的本質(zhì)是利用增長較慢的函數(shù)代替平方和函數(shù),得到具有較好抗粗差性的估值.
(7)
(1)確定待估參數(shù)的初值,確定權(quán)矩陣的初值,給出觀測數(shù)據(jù)的方差-協(xié)方差矩陣;
(2)按照公式(5)和公式(6)所示的傳統(tǒng)最優(yōu)化算法估計參數(shù)的平差值,計算觀測數(shù)據(jù)與模型值的殘差;
(5)對參數(shù)平差值與初值的差值進行評價,如果其符合限差要求,則結(jié)束求解;否則重復(fù)步驟(2)—(4)進行迭代求解,直至待求參數(shù)平差值的差值符合要求為止.其中第k+1次的迭代解為
(8)
中國大陸觀測地殼運動的GNSS測站主要來自于中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò),該網(wǎng)絡(luò)包括260個連續(xù)觀測站和2000個流動觀測站,其中32個連續(xù)觀測站和1056個流動觀測站觀測起始于1998年,其余測站觀測起始于2009年,所有的數(shù)據(jù)均已經(jīng)獲得了超過10年的觀測資料.除此之外,在一些地震活動性較強的區(qū)域,例如首都圈地區(qū)、山西地塹帶和川滇地區(qū)等,也進行了GNSS流動觀測網(wǎng)的加密觀測;省屬的地震、氣象和測繪部門為了滿足區(qū)域性CORS網(wǎng)的需求,在所屬范圍也布設(shè)了連續(xù)觀測站,這些測站中的一部分能夠滿足地殼運動觀測的需求.Wang和Shen(2020)對以上數(shù)據(jù)進行了收集和整理,同時也收集了中國大陸周邊的GNSS觀測資料,通過高精度數(shù)據(jù)處理與震間形變場提取,給出了中國大陸高空間分辨率的GNSS震間速度場.數(shù)據(jù)處理采用GAMIT/GLOBK軟件(Herring et al., 2015, 2018),所有數(shù)據(jù)解算采用統(tǒng)一的策略,確保解算的坐標參數(shù)的高精度和強自洽性.在GNSS原始觀測數(shù)據(jù)解算階段,將地面測站坐標、地球定向參數(shù)(極移和協(xié)調(diào)時間時)和衛(wèi)星軌道一并估算,有效地消除了由于地球物理模型或者軌道參數(shù)變化所帶來的系統(tǒng)性偏差(Wang and Shen, 2020).在時間序列建模階段,同時考慮了大地震同震和震后位移的影響;考慮到不同天線類型和可能的天線設(shè)置誤差,對GAMIT給出的測站坐標誤差進行了重新分配,并利用有色噪聲隨機模型構(gòu)造觀測數(shù)據(jù)的方差-協(xié)方差矩陣,體現(xiàn)出時間序列中的時間相關(guān)性特征(Tao et al., 2021).
在Wang和Shen(2020)給定的GNSS速度場中,中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)的數(shù)據(jù)截止于2015年,本文采用與其一致的觀測數(shù)據(jù)解算策略、時間序列建模策略和地殼運動解算策略,補充處理了該觀測網(wǎng)2017、2018、2019和2020年的觀測數(shù)據(jù),增加了觀測數(shù)據(jù)的長度,也提高了速度場的穩(wěn)定性,更新的川滇地區(qū)GNSS速度場如圖1所示.
圖1 川滇地區(qū)相對于歐亞參考框架的GNSS速度場,誤差橢圓置信區(qū)間為95%Fig.1 GNSS velocity field with respect to Eurasian reference frame in the Sichuan and Yunnan regions. The confidence interval of error ellipse is 95%
圖2 川滇地區(qū)GNSS測站分布的球面小波分解尺度Fig.2 Decomposition scale of spherical wavelet based on the distribution of GNSS stations in the Sichuan and Yunnan regions
根據(jù)GNSS測站分布,依據(jù)空間球面小波分解原則,獲得川滇地區(qū)的球面小波分解尺度,如圖2所示.研究區(qū)內(nèi)最大分解尺度的最小值為6,最大分解尺度的最大值為9.除西南角和東南角的邊緣區(qū)域,其他區(qū)域的最大分解尺度均在7以上.最大分解尺度能達到9的區(qū)域主要沿著主要構(gòu)造帶分布,例如小江斷裂帶南端、紅河斷裂與麗江—小金河斷裂交匯的區(qū)域和鮮水河斷裂與龍門山斷裂南端交匯的三岔口區(qū)域.川滇菱形塊體覆蓋的主要區(qū)域最大分解尺度都能達到8,因此本文最終選擇的最大分解尺度為8,分解尺度的選擇范圍為3~8.
正則化參數(shù)的求解采用留一交叉驗證法(Leave-one-out Cross Validation,LOOCV),首先去除一個觀測點數(shù)據(jù),計算待求參數(shù),然后計算該點觀測值與模型估值之間的差值,通過全局搜索,選擇出最優(yōu)的正則化參數(shù),進而求解出待求參數(shù)的最優(yōu)解,正則化參數(shù)的搜索過程如圖3所示.選擇1~1010作為搜索范圍,可以看出該過程搜索到了函數(shù)的全局最小值,進而提取出正則化參數(shù)λ的取值,圖3中橫坐標取正則化參數(shù)λ的常用對數(shù),其在東西分量和南北分量分別為533和849.
圖3 基于LOOCV算法的最優(yōu)化算法正則化參數(shù)Fig.3 Selection of regulation parameters based on optimized LOOCV method
確定球面小波分解尺度和正則化參數(shù)之后,通過多次迭代即可求解出待估參數(shù)的最優(yōu)化值,進而正演出研究范圍內(nèi)空間任意一點的地殼運動速度,實現(xiàn)地殼形變場建模,給出空間連續(xù)的地殼運動圖像,如圖4所示,圖中冷色(藍色)代表小速率值,暖色(紅色)代表大速率值.可以看出,安寧河—則木河—小江斷裂組成的構(gòu)造帶是川滇菱形塊體地殼運動速率的分界帶.
圖4 川滇地區(qū)地殼運動場模型(a) 東西分量; (b) 南北分量.Fig.4 Crustal deformation model in the Sichuan and Yunnan regions(a) The EW component; (b) The NS component.
抗差迭代算法可有效削弱異常觀測值和粗差對地殼形變場建模的影響,進而獲得準確的應(yīng)變率場.通過多次抗差迭代計算,在所有參與建模的367個GNSS測站中,檢測出74個降權(quán)觀測值,占比為20%,檢測出7個除權(quán)觀測值,占比為2%,圖5中藍色三角形站點為正常觀測站點,黑色三角形站點為降權(quán)觀測站點,紅色三角形站點為除權(quán)觀測站點.
圖5 抗差迭代算法給出的川滇地區(qū)三類GNSS測站分布圖Fig.5 The distribution of three types of GNSS stations in Sichuan and Yunnan regions based on robust iterative algorithm
圖6給出了利用抗差迭代算法對部分測站進行降權(quán)和除權(quán)之后,基于多尺度球面小波算法構(gòu)建的川滇地區(qū)地殼運動場模型,可以看出考慮抗差性前后計算的速度場空間分布形態(tài)基本一致,安寧河—則木河—小江斷裂是川滇菱形塊體的東邊界這一基本特征未發(fā)生明顯改變.為了進一步分析抗差迭代算法在構(gòu)建地殼形變場時的作用,對本文顧及抗差性的方法與傳統(tǒng)多尺度球面小波方法獲得的地殼運動場模型進行對比,結(jié)果如圖7所示,圖中也給出了降權(quán)站點和除權(quán)站點的分布.可以看出,二者差異的量值在1.0 mm·a-1以內(nèi),說明雖然在構(gòu)建速度場模型時考慮了抗差性,但由于原有速度場本身數(shù)據(jù)質(zhì)量較好,因此得到的結(jié)果是穩(wěn)定的,二者的差異較小,而差異性較大的區(qū)域總是覆蓋著被降權(quán)或者除權(quán)的測站.
利用目視法對基于GNSS觀測資料解算的中國大陸地殼運動場進行分析,發(fā)現(xiàn)不可避免地存在一些與周圍測站的運動特征明顯不協(xié)調(diào)的測站.其產(chǎn)生的原因主要包括:(1)因觀測時間跨度太短而無法擬合出真正的震間運動;(2)觀測標墩建于土層之上,由于季節(jié)變化等因素導(dǎo)致標墩本身存在非線性的變化;(3)測站觀測環(huán)境的改變導(dǎo)致觀測噪聲增加或運動特征偏離真值.
本文所采用的觀測數(shù)據(jù)均超過了10年,超過解算震間變形時對時間跨度4.5年的要求(Blewitt et al., 2013).觀測環(huán)境的改變在測站的日志文件中都會有明確記載,在數(shù)據(jù)預(yù)處理階段已經(jīng)被考慮到.因此,為了分析測站被降權(quán)處理和被除權(quán)處理的原因,本文重點對這些測站的觀測標墩和數(shù)據(jù)質(zhì)量進行了進一步分析.其中測站標墩分為基巖點和土層點,通常認為基巖點能夠精準代表地殼的運動特征,而土層點會受到土壤隨著季節(jié)變化的熱脹冷縮和土層本身不穩(wěn)定等因素的影響(譚偉杰等,2017),產(chǎn)生局部效應(yīng)而可能無法精準代表地殼的運動特征.被降權(quán)和除權(quán)處理的81個站點全部為土層站點,可見觀測結(jié)果中存在明顯的局部影響.通過觀測數(shù)據(jù)形式誤差的大小體現(xiàn)出數(shù)據(jù)質(zhì)量的差異性,分析權(quán)不變、降權(quán)和除權(quán)三類測站形式誤差的差異性.計算三類數(shù)據(jù)東西向和南北向形式誤差的均值,權(quán)不變、降權(quán)和除權(quán)三類數(shù)據(jù)東西向形式誤差的均值分別為0.61、0.70和0.93 mm·a-1,南北向形式誤差的均值分別為0.52、0.60和0.83 mm·a-1,可以看出三類數(shù)據(jù)的形式誤差存在差異,表明其時間序列的離散度逐漸增大,數(shù)據(jù)質(zhì)量逐漸降低.
圖6 顧及抗差迭代算法的川滇地區(qū)地殼運動場模型(a) 東西分量; (b) 南北分量.Fig.6 Crustal deformation model in Sichuan and Yunnan regions considering robust iterative algorithm(a) The EW component; (b) The NS component.
圖7 顧及抗差迭代算法與傳統(tǒng)方法構(gòu)建的地殼運動場在東西分量的差值Fig.7 The difference of crustal deformation model in the EW component between our and traditional method
為了進一步檢驗抗差迭代算法的抗差性能,在速度場空間變化較大的川滇地區(qū)人為地引入粗差,分析粗差引入前后速度場模型的差異性特征.引入粗差的方式有在單個測站引入粗差和在相鄰的幾個測站引入粗差兩種,具體如圖8所示,其中測站H335、H387和H380為3個在區(qū)域上獨立的含有粗差的測站,測站H346、H111、H112、SCYY和H103為在區(qū)域上具有相關(guān)性的5個含有粗差的測站.用含有粗差的測站構(gòu)建運動場模型,并與不含粗差的數(shù)據(jù)構(gòu)建的結(jié)果進行對比,分析本文算法的抗差性.
圖8 速度場中單個測站和區(qū)域測站加入人為誤差的兩種方案Fig.8 Two schemes of adding gross errors artificially of individual station and regional stations in velocity field
圖9 加入人為粗差前后東西向速度場模型的差異分布Fig.9 The difference of crustal deformation model in EW component before and after adding gross errors artificially
圖9給出了對二者速度場求差的結(jié)果,圖中-0.5~0.5 mm·a-1以內(nèi)的值用同一個色標表示,可以看出,其中測站H335和H387兩個單獨的存在粗差的測站經(jīng)過抗差迭代算法之后,對建模的結(jié)果并未產(chǎn)生影響,而另一個獨立的H380測站的粗差卻對建模結(jié)果產(chǎn)生了明顯的影響,進一步對測站分布分析可以看出,在H335和H387兩個測站10 km以內(nèi)的范圍還有其他測站,而H380測站在10 km以內(nèi)的范圍內(nèi)沒有其他測站,也就是說抗差迭代算法從測站H335和H387檢測出了粗差,而從H380測站沒有檢測到粗差,把含有粗差的觀測值當作好的觀測值對待,這是由于該測站周邊沒有其他觀測站作為對比,即便其本身含有粗差也被當做了有效值.由測站SCYY、H346、H103、H111和H112等5個測站組成的含有粗差的區(qū)域?qū)=Y(jié)果也產(chǎn)生了明顯的影響,表明區(qū)域性存在的粗差會對建模結(jié)果產(chǎn)生影響.
從圖10可以看出,川滇地區(qū)重要應(yīng)變積累的區(qū)域沿著主要的活動構(gòu)造帶分布,其中主要的一個條帶是安寧河—則木河—小江斷裂帶,這也是川滇菱形塊體的東邊界,紅河斷裂與麗江—小金河斷裂交匯的區(qū)域是另外一個重要的應(yīng)變積累區(qū)域,該分布特征與已有的研究成果保持一致(Qu et al., 2018),但在應(yīng)變集中區(qū)的量值上普遍高于已發(fā)表的結(jié)果(Pan and Shen, 2017;李長軍等,2019),這主要是由于本文采用的多尺度球面小波的應(yīng)變計算方法能夠凸顯出測站不均勻分布的特征,獲取不同空間分辨率的應(yīng)變率場,進而提取出應(yīng)變積累的集中區(qū)域.圖10也給出本文方法與傳統(tǒng)多尺度球面小波方法解算的應(yīng)變率場的差異,從第二不變量的量值來看,應(yīng)變差異值的量值遠低于應(yīng)變的量值,表明兩個結(jié)果均是穩(wěn)健的,對比圖10和圖7可以看出,應(yīng)變率差異較大的區(qū)域也是速度場差異較大的區(qū)域,同時也是存在被降權(quán)或除權(quán)測站的區(qū)域.
基于多尺度球面小波算法解算應(yīng)變率場的優(yōu)點之一是能夠獲取不同空間尺度的應(yīng)變率場,從而體現(xiàn)出不同規(guī)模構(gòu)造運動的應(yīng)變積累特征.小的分解尺度能夠獲取大規(guī)模構(gòu)造運動的應(yīng)變積累,例如塊體邊界的應(yīng)變積累;大的分解尺度能夠獲取小規(guī)模構(gòu)造運動的應(yīng)變積累,例如塊體邊界應(yīng)變積累的差異性分布或塊體內(nèi)部次級斷層的應(yīng)變積累.圖11分別給出了分解尺度為3~7和8的應(yīng)變率第二不變量,二者之和為區(qū)域的總應(yīng)變,總體特征與已有結(jié)果保持一致(Jin et al, 2019; Zhang et al., 2019).從分解尺度為3~7的應(yīng)變率場可以看出,應(yīng)變積累的主要區(qū)域為塊體運動的邊界區(qū)域,作為川滇菱形塊體東邊界的安寧河—則木河—小江斷裂帶,其是川滇地區(qū)主要的應(yīng)變積累區(qū)域,其中以安寧河斷裂和則木河斷裂為東邊界的大涼山次級塊體內(nèi)部存在應(yīng)變積累,沿著小江斷裂的應(yīng)變積累穿過紅河斷裂之后一直向西南方向延伸.其他主要的應(yīng)變積累區(qū)域是位于紅河斷裂與麗江—小金河交匯的區(qū)域和中緬邊界的幾條西南向斷裂分布區(qū)域,主要包括畹町斷裂和南汀河斷裂等.分解尺度為8的應(yīng)變積累分布顯示區(qū)域內(nèi)除了沿著主要斷裂的應(yīng)變積累外,還存在局部小區(qū)域的應(yīng)變積累,這些應(yīng)變積累主要體現(xiàn)了沿著斷裂應(yīng)變積累的差異性特征.
圖10 川滇地區(qū)應(yīng)變率第二不變量(a) 本文方法解算的應(yīng)變率場; (b) 本文方法與傳統(tǒng)多尺度球面小波方法解算結(jié)果的差異.Fig.10 The second invariant of strain rates in Sichuan and Yunnan regions(a) The strain rate field calculated by our method; (b) The difference between our and traditional method.
圖11 川滇地區(qū)應(yīng)變率第二不變量(a) 分解尺度為3~7; (b) 分解尺度為8.Fig.11 The second invariant of strain rates in Sichuan and Yunnan regions(a) The decomposition scale is from 3 to 7; (b) The decomposition scale is 8.
利用長期積累的GNSS觀測資料經(jīng)過高精度統(tǒng)一數(shù)據(jù)處理提取地殼運動圖像,然后解算GNSS應(yīng)變率場,進而對應(yīng)變積累與地震矩釋放進行對比分析,這些是開展區(qū)域地震危險性分析的主要工作,并在實際地震潛勢研究中得到了廣泛應(yīng)用和驗證.由于GNSS時間序列中存在與長期構(gòu)造運動無關(guān)的信息致使擬合的速度場無法真實表達實際的震間地殼運動特征,進而影響到應(yīng)變率場的結(jié)果.
本文提出了顧及抗差性的速度場模型構(gòu)建方法,考慮到了測站因局部運動和數(shù)據(jù)質(zhì)量較差在速度場中產(chǎn)生的粗差.以解算的川滇地區(qū)最新GNSS觀測資料為例,在觀測的速度場中檢測出了20%的降權(quán)測站和2%的除權(quán)測站,對這些測站進一步分析發(fā)現(xiàn),粗差產(chǎn)生的原因是由于觀測數(shù)據(jù)質(zhì)量較差和測站為土層觀測點而存在局部非構(gòu)造變形,進而獲得穩(wěn)健性較好的速度場模型.通過在觀測數(shù)據(jù)人為加入粗差檢驗本文方法的抗差性,結(jié)果表明單個存在的粗差對速度場模型不會產(chǎn)生影響,而區(qū)域性存在的粗差會對結(jié)果產(chǎn)生影響,表明本文的方法具有好的抗差性.本文進一步分析了顧及抗差性算法對解算應(yīng)變率場的貢獻,進而給出了川滇地區(qū)精準的應(yīng)變率場.
致謝感謝中國地震臺網(wǎng)中心提供的中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡(luò)GNSS觀測數(shù)據(jù).