摘要：2023年12月18日，甘肅省臨夏州積石山縣發(fā)生MS6.2地震，造成了大量的人員傷亡和經(jīng)濟損失。為了研究該地震的發(fā)震構(gòu)造，以多源觀測為約束，利用GNSS和InSAR數(shù)據(jù)分別獨立和聯(lián)合反演其同震破裂模型;基于反演得到的模型，計算積石山地震在周邊斷裂上產(chǎn)生的同震庫侖應(yīng)力變化，并進一步分析其對周圍地區(qū)的應(yīng)力影響。結(jié)果顯示：積石山地震的發(fā)震斷層走向310°，傾角51°，是一次以逆沖為主兼有少量右旋走滑分量的地震，最大滑動量為56.02 cm，大致位于5 km深度處;該地震為調(diào)節(jié)拉脊山構(gòu)造轉(zhuǎn)換帶兩側(cè)區(qū)域的應(yīng)力-應(yīng)變關(guān)系發(fā)揮著重要作用。地震發(fā)生后，積石山東緣斷裂南段、拉脊山南緣斷裂中南段、倒淌河—臨夏斷裂中段，以及西秦嶺北緣斷裂西段受到的庫侖應(yīng)力加載較高，這些區(qū)域未來的地震危險性需要重點關(guān)注。

關(guān)鍵詞：積石山地震; GNSS; InSAR; 同震庫侖應(yīng)力變化; 地震危險性

中圖分類號： P315文獻標志碼：A文章編號： 1000-0844（2024）04-0867-13

DOI：10.20000/j.1000-0844.20240326004

GNSS and InSAR derived coseismic slip distribution

of the 2023 Jishishan MS6.2 earthquakeSUN Yunqiang， QIU Xinpeng， CHEN Changyong， ZHAO Zhuo， LIN Junqi， GONG Weicheng

（College of Transportation and Civil Engineering，F(xiàn)ujian Agriculture and Forestry University，F(xiàn)uzhou 350002，F(xiàn)ujian，China）Abstract：

On December 18， 2023， an MS6.2 earthquake struck Jishishan County， Linxia， Gansu， resulting in significant casualties and property losses. This study aims to investigate the seismogenic structure of the earthquake using multisource observational data. The coseismic rupture model was inverted independently and jointly using GNSS and InSAR data. Based on this model， the coseismic Coulomb stress changes on surrounding faults were calculated to analyze the stress influence caused by the MS6.2 earthquake in the surrounding areas. The seismogenic fault of the Jishishan earthquake was found to have a strike of 310° and a dip angle of 51°. The fault movement was primarily characterized by thrust with a small right-lateral strike-slip component. The maximum slip occurred at a depth of approximately 5 km， with a maximum displacement of 56.02 cm. This earthquake played an important role in regulating the stress-strain relationship across the Lajishan tectonic transition zone. Post-earthquake analysis revealed that high-stress loading was observed in several fault sections： the southern section of the eastern margin fault of Jishishan， the south-central section of the southern margin fault of Lajishan， the middle section of the Daotanghe-Linxia fault， and the western section of the northern margin of the western Qinling fault. These findings indicate that the seismic hazards of these faults require close monitoring in the future.

Keywords：Jishishan earthquake; GNSS; InSAR; coseismic Coulomb stress changes; seismic hazards

0引言

根據(jù)中國地震臺網(wǎng)中心（https：//news.ceic.ac.cn/）測定，北京時間2023年12月18日23時59分，中國甘肅省臨夏州積石山縣發(fā)生MS6.2地震（35.7°N，102.79°E），震源深度10 km，最高烈度達Ⅷ度（圖1）。全球?qū)Ш蕉ㄎ幌到y(tǒng)（Global Navigation Satellite System，GNSS）捕捉到此次地震的同震形變覆蓋震中約50 km的區(qū)域范圍，總體與地震烈度結(jié)果保持一致［1］。在主震發(fā)生后，該地區(qū)還發(fā)生了多次余震，截至2023年12月22日，共發(fā)生3.0級以上的余震20次［2］。積石山地震是一次典型的“主震-余震”型地震［3］，造成了大量的人員傷亡（148人死亡，979人受傷，3人失蹤）和基礎(chǔ)設(shè)施毀壞，給社會經(jīng)濟造成巨大損失［4-5］。

青藏高原東北緣作為青藏高原朝NE向擴展的活動前緣［6-7］，長期以來受到印度板塊和歐亞板塊的碰撞擠壓作用，構(gòu)造活動強烈，而強烈的構(gòu)造活動在控制著該區(qū)域構(gòu)造演化的同時也深刻影響其地震活動［7-8］。本次地震序列就發(fā)生于青藏高原東北緣的拉脊山斷裂帶上。拉脊山斷裂帶形成于加里東期，主要由拉脊山北緣斷裂、拉脊山南緣斷裂和積石山東緣斷裂組成［9］，起到調(diào)節(jié)日月山右旋走滑和西秦嶺北緣左旋走滑運動之間應(yīng)力-應(yīng)變關(guān)系的作用。拉脊山斷裂帶經(jīng)歷過多次的擠壓變形與地殼運動，最終形成了現(xiàn)今較為復雜的朝NE向凸出的弧形構(gòu)造形態(tài)［10］。據(jù)歷史記載，拉脊山南北兩側(cè)共發(fā)生過20余次5級以上的中強地震［10］，但歷史強震的離逝時間較為久遠，導致該區(qū)域經(jīng)歷了較長時間的中強地震平靜期。此次積石山地震打破了該區(qū)域中強地震平靜期，因此，其未來的地震危險性需要進一步評估。

地震破裂過程對于深入了解地震的發(fā)震構(gòu)造和破裂過程，以及進一步分析區(qū)域地震危險性具有十分重要的意義。積石山地震發(fā)生后的短時間內(nèi)，大量學者投入到對其發(fā)震斷層的破裂反演研究中［2，9，11-14］。例如，王衛(wèi)民等（https：//www.itpcas.cas.cn/new_kycg/new_kyjz/202312/t20231219_6945583.html）利用遠場波形數(shù)據(jù)反演了積石山地震的破裂過程，發(fā)現(xiàn)其由一個主破裂區(qū)和一個次級破裂區(qū)組成，主破裂長約15 km，最大滑動量20 cm。張勇等（https：//pku-geophysics-source.group/htmls/20231218155930Jishishan.html）利用近場的強震記錄對該地震的破裂過程進行了反演，揭示其為NNW向的單側(cè)破裂形式，主破裂長約12 km，最大滑動量37 cm。為避免遠震波形數(shù)據(jù)在傳播過程中受到各種因素的影響，羅艷等［13］利用區(qū)域的寬頻帶地震波形記錄反演了積石山地震的破裂過程，發(fā)現(xiàn)其最大錯動量約為35 cm，主體破裂尺度約10 km。楊九元等［14］基于合成孔徑雷達干涉測量（Interferometric Synthetic Aperture Radar，InSAR）同震形變場反演了積石山地震的斷層幾何參數(shù)和同震滑動分布，結(jié)果顯示，其同震滑動主要分布于5.6～16.1 km的深度范圍，最大滑動量可達66 cm。目前，關(guān)于積石山MS6.2地震的同震破裂模型很多，但大多是通過單一的觀測數(shù)據(jù)反演得到的。僅用單一的觀測數(shù)據(jù)難以保證對震中區(qū)域的空間全覆蓋，采用多種數(shù)據(jù)進行聯(lián)合反演不僅能使數(shù)據(jù)覆蓋得更全面，而且可以實現(xiàn)不同數(shù)據(jù)之間的相互約束［15-16］。如今，已有大量研究采用以GNSS和InSAR為代表的多種數(shù)據(jù)進行聯(lián)合反演［15，17-19］，這種方式能夠在一定程度上彌補單一形變觀測數(shù)據(jù)的缺陷［15，19-20］。

基于此，本文首先采用Okada［21］的彈性位錯模型來計算矩形斷層滑動引起的地表位移，并使用混沌麻雀搜索優(yōu)化算法［22］（Chaos Sparrow Search Optimization Algorithm，CSSOA）求解斷層的最優(yōu)滑動分布;然后，以GNSS和InSAR數(shù)據(jù)為約束，分別獨立和聯(lián)合反演了積石山MS6.2地震的同震破裂模型;最后，基于反演的同震破裂模型計算了該地震在周邊斷裂上產(chǎn)生的同震庫侖應(yīng)力變化，以期進一步理解積石山地震的發(fā)震機制，為該區(qū)域的地震危險性分析提供參考。

1GNSS和InSAR數(shù)據(jù)

目前，包括GNSS和InSAR在內(nèi)的大地測量數(shù)據(jù)已經(jīng)被廣泛地應(yīng)用于地震的同震形變研究中［23］。大地測量數(shù)據(jù)，特別是近場位移，比遠震波對斷層幾何和滑動分布更敏感［24］，更有利于區(qū)域斷層的發(fā)震構(gòu)造和破裂行為研究。

本次研究所使用的GNSS數(shù)據(jù)來自李志才等［1］對32個衛(wèi)星導航定位基準站進行高精度動態(tài)單歷元數(shù)據(jù)處理后的結(jié)果，覆蓋了震中附近50 km范圍內(nèi)的同震形變情況，共包含5個站點的數(shù)據(jù)，即圖1中的GUTI、LXYJ、XUNH、LXJS和gsl9。各站點均較為完整、清晰地記錄到了地震引起的瞬時變形。然而，積石山地震周邊的GNSS觀測站點較少，獲得的GNSS觀測數(shù)據(jù)相對稀疏，所以同震形變誤差也相對較大（置信度35%）。其中，LXJS站點距震中最近，觀測到的同震形變也最高，東西向形變達（13.6SymbolqB@

2.6） mm，南北向形變達（10.9SymbolqB@

5.1） mm［1］，遠高于其余站點的觀測值。

本次研究采用的InSAR數(shù)據(jù)來源于歐洲空間局（https：//dataspace.copernicus.eu/）Sentinel-1衛(wèi)星（表1），通過解算提取了如圖2所示的積石山地震前后的InSAR同震干涉圖和同震視線向（Line of Sight，LOS）形變圖。結(jié)果顯示，InSAR升、降軌觀測數(shù)據(jù)獲得的同震干涉圖［圖2（a）、（b）］均為一個近似橢圓狀的形變條紋，圖像內(nèi)還有少量的失相干區(qū)域，或許與地震引發(fā)的滑坡、泥流等災害有關(guān)［4-5］。InSAR升、降軌的LOS向形變圖［圖2（c）、（d）］顯示的隆升與沉降分界線大致沿NW—SE方向，可大致確定發(fā)震斷層位于拉脊山南緣斷裂和積石山東緣斷裂的交界區(qū)域，且主要沿NW—SE方向破裂。發(fā)震斷層北側(cè)的同震形變以朝向于衛(wèi)星的抬升運動為主，結(jié)合區(qū)域構(gòu)造形態(tài)特征［4-5］，認為積石山地震符合逆沖型地震的形變特征，且暗示了其發(fā)震斷層可能為拉脊山南緣斷裂未出露地表的隱伏逆沖分支斷層。

為了提高InSAR反演的效率并降低遠場噪聲對反演結(jié)果的影響，采用均勻采樣方法對積石山地震的InSAR升、降軌同震形變場進行降采樣處理，在提高震中周邊區(qū)域采樣點密度的同時降低遠離震中區(qū)域的采樣點密度，最終獲取了升、降軌各1 500個形變數(shù)據(jù)點。

2反演方法

研究者們曾使用模擬退火算法（Simulated Annealing Algorithm，SAA）［25-26］、粒子群優(yōu)化（Particle Swarm Optimization，PSO）算法［12，27］、灰狼優(yōu)化（Grey Wolf Optimization，GWO）算法［28］等方法進行斷層滑動的反演。這些算法雖然在一定程度上提高了搜索能力，卻存在易陷入局部最優(yōu)、收斂速度慢等缺陷。麻雀搜索算法（Sparrow Search Algorithm，SSA）［29］是近年來新興的一種群體智能優(yōu)化算法，能夠在一定程度上解決收斂速度慢的問題，同時具有搜索精度高、求解效率高、實現(xiàn)簡單等優(yōu)點。但是，當搜索接近全局最優(yōu)、種群多樣性減少時，該算法還是容易陷入局部最優(yōu)狀態(tài)。而CSSOA在SSA的基礎(chǔ)上引入了Tent混沌和高斯變異，增強了個體在搜索過程中跳出局部最優(yōu)的能力，并進一步提高了收斂速度［22］。

因此，本文采用Okada［21］的彈性位錯模型計算斷層滑動引起的地表位移，并使用CSSOA［22］反演斷層的最優(yōu)滑動分布。同時，為了提高反演精度，CSSOA反演策略中的種群規(guī)模N=100，最大迭代次數(shù)T=50 000，收斂誤差1×10-10，發(fā)現(xiàn)者和偵察預警麻雀數(shù)占種群規(guī)模的20%［22］。

3結(jié)果分析

本節(jié)將以GNSS和InSAR數(shù)據(jù)的不同聯(lián)合為約束，結(jié)合Okada彈性位錯模型和CSSOA，反演得到積石山地震的斷層幾何參數(shù)和同震滑動分布。之后，基于反演得到的同震破裂模型計算積石山地震產(chǎn)生的同震庫侖應(yīng)力變化，以進一步分析周邊斷層未來的地震危險性。

3.1斷層幾何

對積石山地震發(fā)震斷層幾何的研究中，不同研究者從不同的角度給出了自己的判斷［2，9，11-14］。他們對于走向的判斷較為一致，認為主要沿NW—SE方向，這與本文獲取的InSAR同震LOS向形變結(jié)果一致［圖2（c）、（d）］。對于發(fā)震斷層的傾向，研究結(jié)果卻有較大差異［9，11-14］。其中，劉振江等［12］認為，積石山地震的發(fā)震斷層傾向為NE，這與InSAR反演結(jié)果以及余震和同震滑坡的空間分布均較為一致［4，12］。因此，本文選擇NE方向作為積石山地震的發(fā)震斷層傾向。

基于積石山地震NE傾的假設(shè)，在反演時設(shè)置搜索區(qū)間為走向［270°，360°］、傾角［40°，90°］。結(jié)果顯示，當斷層走向310°、傾角51°時，能夠與InSAR數(shù)據(jù)實現(xiàn)最佳擬合［圖4（a）］。該結(jié)果與楊九元等［14］、張喆等（https：//www.cea-igp.ac.cn/cxdt/280419.html）的反演結(jié)果較為接近。

根據(jù)擬合得到的斷層幾何參數(shù)，反演了積石山地震的同震滑動分布。通過固定斷層的走向、傾角和位置，沿走向和傾向?qū)鄬用娴拈L度和寬度均擴大至20 km，等間距地將斷層面離散為30 × 20個小矩形斷層單元。為了權(quán)衡數(shù)據(jù)擬合度和滑動分布的平滑度，根據(jù)測試結(jié)果，設(shè)置最優(yōu)平滑因子為0.005［圖4 （b）］。

3.2同震滑動分布反演

3.2.1GNSS獨立反演

以GNSS同震觀測數(shù)據(jù)為約束，獨立反演了積石山地震的同震滑動分布。結(jié)果顯示（圖5），同震滑動分布幾乎遍布整個斷層面，呈現(xiàn)出由斷層面中心向四周遞減的形式，在中部（5～15 km）呈現(xiàn)逆沖兼右旋的特征，而在～3 km深度處和模型深部則更傾向于純右旋走滑的趨勢。中心的高滑移區(qū)近似圓形，最大滑動集中于7.5～10 km深度，最大滑動量為25.64 cm，在地表還有接近15 cm的滑動聚集。平均滑動角130.07°，釋放的地震矩約1.11×1018 N·m，相當于一次MW6.0地震。

3.2.2InSAR獨立反演

以InSAR數(shù)據(jù)為約束，獨立反演了積石山地震的同震滑動分布。結(jié)果顯示（圖6），同震滑動的分布特征相對單一，近似呈半橢圓形，以逆沖為主，伴隨有少量的左旋走滑特征。同震滑動分布于2.5～20 km深度范圍，主要由5 km深度處向兩側(cè)及更深處遞減;最大滑動量為57.69 cm，位于～5 km深度。地表淺層0～2.5 km深度范圍內(nèi)未有滑動分布存在，暗示積石山地震是一次隱伏的逆沖型地震。破裂長度大致為15 km，平均滑動角77.30°，釋放的地震矩約2.47×1018 N·m，相當于一次MW6.2地震。

3.2.3GNSS和InSAR聯(lián)合反演

圖5、6的結(jié)果顯示，使用GNSS和InSAR數(shù)據(jù)進行獨立反演時存在很大差異，不僅同震滑動的運動性質(zhì)不一致，滑動量也相差了近兩倍?？梢?，不同觀測數(shù)據(jù)的單一反演無法準確、真實地表達地震破裂形態(tài)。

基于衛(wèi)星獲取的InSAR圖像可以完整記錄不規(guī)則地表的同震一維LOS向形變［33］，但InSAR技術(shù)可能出現(xiàn)相位失相干和平行衛(wèi)星軌道跡線的斷層形變精度缺失［14，34］，因此，很難僅依靠InSAR觀測數(shù)據(jù)獲得真實可靠的三維同震形變場［35］。GNSS可以觀測到同震的三維形變場，實現(xiàn)長期的時序形變監(jiān)測，但積石山地震周邊的GNSS站點稀疏，無法捕捉到大范圍的同震形變場的細節(jié)。很明顯，這兩種數(shù)據(jù)的觀測特性能夠?qū)崿F(xiàn)互補，且為斷層的滑動分布反演帶來更豐富的約束。因此，本研究結(jié)合GNSS和InSAR數(shù)據(jù)進行聯(lián)合反演，使兩種數(shù)據(jù)相互約束。

由于GNSS數(shù)據(jù)較為稀疏，且置信率較低，因此在聯(lián)合反演中給予它較小的權(quán)重，并相應(yīng)地提高InSAR數(shù)據(jù)的權(quán)重，最終將GNSS和InSAR數(shù)據(jù)聯(lián)合反演的相對權(quán)重設(shè)置為0.2∶0.8。

GNSS和InSAR數(shù)據(jù)聯(lián)合反演的結(jié)果（圖7）顯示，同震滑動分布主要呈現(xiàn)出逆沖兼右旋的特征，越接近地表，右旋走滑分量逐漸減少，進而轉(zhuǎn)變成以逆沖為主兼有少量右旋走滑分量的形式。主破裂長度大致為15 km;同震滑動主要分布于2.5～7.5 km的深度范圍，最大滑動位于5 km深度，最大滑動量為56.02 cm;10～20 km深度范圍內(nèi)還散布有多處滑動量為20～30 cm的次級破裂。平均滑動角102.61°，釋放的地震矩約1.79×1018 N·m，相當于一次MW6.1地震。

3.3同震庫侖應(yīng)力變化

地震是在一定構(gòu)造應(yīng)力條件下斷裂達到極限狀態(tài)突然破裂失穩(wěn)的結(jié)果，一次地震的發(fā)生會導致區(qū)域斷裂的應(yīng)力狀態(tài)發(fā)生改變［36］。計算大地震產(chǎn)生的庫侖應(yīng)力變化，對于評估區(qū)域未來的地震危險性具有重要意義［37-38］。

為了分析積石山地震對周邊主要斷裂的同震庫侖應(yīng)力影響，將庫侖應(yīng)力投影至各接收斷層上（接收斷層參數(shù)列于表2，深度為5 km）。結(jié)果顯示（圖8），地震產(chǎn)生的庫侖應(yīng)力變化對震中50 km范圍內(nèi)的影響較大，即庫侖應(yīng)力變化超過千帕量級的斷裂段主要集中在此范圍內(nèi)。積石山東緣斷裂中段、拉脊山南緣斷裂東南段和倒淌河—臨夏斷裂中東段的庫侖應(yīng)力變化為負，其庫侖應(yīng)力卸載最高分別為72.61、15.41和9.92 kPa。莊浪河斷裂、馬銜山北緣斷裂和西秦嶺北緣斷裂中段也產(chǎn)生了一定的應(yīng)力卸載作用，但影響程度較低，僅有百帕量級。

如圖8所示，積石山地震發(fā)生后，周邊的多條斷裂都受到了應(yīng)力加載效應(yīng)的影響，加載效應(yīng)的傳遞主要沿NW—SE方向，并隨著震中距的增加而減小。其中，臨潭—宕昌斷裂受到的應(yīng)力加載影響最小，應(yīng)力加載量不到百帕;拉脊山北緣斷裂和日月山斷裂受到的影響也相對較低，分別為0.20 kPa和0.22 kPa;西秦嶺北緣斷裂西段和倒淌河—臨夏斷裂中段也受到了一定的應(yīng)力加載作用，庫侖應(yīng)力最大增加量分別為1.60 kPa和2.10 kPa。積石山東緣斷裂和拉脊山南緣斷裂與震中較為接近，因此所受應(yīng)力加載影響較大。其中，積石山東緣斷裂受到的應(yīng)力加載最高達15.58 kPa，加載區(qū)域主要位于南段;拉脊山南緣斷裂受到的應(yīng)力加載最高達15.81 kPa，加載區(qū)域主要位于中南段；兩者均超過了10 kPa的應(yīng)力觸發(fā)閾值［36］，地震危險性增大。雖然積石山地震的發(fā)生對周邊斷裂的影響僅有千帕，但考慮到該區(qū)域歷史強震的離逝時間較為長久以及構(gòu)造的持續(xù)加載，這些斷裂未來的地震危險性仍需重點關(guān)注。

4討論

4.1不同數(shù)據(jù)反演結(jié)果差異

上文基于GNSS和InSAR數(shù)據(jù)對積石山地震的同震破裂模型進行了獨立反演，在不同數(shù)據(jù)約束下，最終的反演結(jié)果呈現(xiàn)一定的差異。其中，單一GNSS數(shù)據(jù)反演的滑動分布在地表呈現(xiàn)大量純右旋走滑特征，斷層中部轉(zhuǎn)變?yōu)槟鏇_兼右旋的形式（圖5），而單一InSAR數(shù)據(jù)反演的滑動分布在斷層中則呈現(xiàn)出集中逆沖兼少量左旋的形式（圖6）。由于InSAR技術(shù)可以獲取到高精度的一維LOS向形變，特別是地表垂向形變，而GNSS獲取的是地表的三維形變場，具有更高的水平分量精度，若僅依靠單一數(shù)據(jù)進行反演，滑動分布形態(tài)會更傾向于所使用數(shù)據(jù)的自身優(yōu)勢，難以全面反映真實的地震破裂過程［35］。

GNSS和InSAR數(shù)據(jù)的聯(lián)合反演能夠使多方互為約束，同時保證水平分量和垂直分量的精度，從而提高反演結(jié)果的可靠性。聯(lián)合反演的結(jié)果顯示，積石山地震的同震滑動分布從斷層底部向地表逐漸由右旋走滑向逆沖轉(zhuǎn)變，呈現(xiàn)出以逆沖為主兼有少量右旋走滑分量的形式（圖7），且最大滑動分布靠近地表（～5 km）。這可能表明了此次積石山地震更傾向于淺層釋放［27］，也在一定程度上解釋了其震級不高卻震害嚴重的原因。

同時，在進行聯(lián)合反演時，數(shù)據(jù)的權(quán)重主要根據(jù)其置信度進行分配，反演結(jié)果自然會受到權(quán)重比的影響。未來在嘗試聯(lián)合反演時，除了需要考慮多種數(shù)據(jù)的最優(yōu)組合外，還需要進一步考慮權(quán)重的合理分配，以期最大程度地吻合現(xiàn)場觀測值。

針對積石山地震，有學者采用波形數(shù)據(jù)和大地測量數(shù)據(jù)進行了一系列的反演工作［2，9，11-14］。由于對地震的傾向尚未取得共識，而傾向的差異又會對結(jié)果產(chǎn)生較大影響，因此本文主要與其他NE傾的同震破裂模型進行對比。不同研究反演得到的積石山地震的同震滑動分布形態(tài)各不相同，但都顯示其是一次以逆沖為主兼少量右旋走滑分量的地震，與本文的聯(lián)合反演結(jié)果一致。不同研究者得出的最大滑動量結(jié)果存在較大差異：最小的反演結(jié)果僅為24 cm［2］，而最大的反演結(jié)果達66 cm［14］，總體呈現(xiàn)出InSAR數(shù)據(jù)結(jié)果略高于波形數(shù)據(jù)的現(xiàn)象，這或許與InSAR數(shù)據(jù)對近場形變更為敏感有關(guān)［24］。大部分研究顯示同震最大滑動在10 km深度上下浮動，而李雨森等［9］的反演結(jié)果則顯示最大滑動位于地表淺層，呈現(xiàn)出由地表淺層向四周破裂的發(fā)震形式，這可能與滑動限制有關(guān)（考慮到積石山地震淺層破裂的存在［27］）。

4.2最大剪切應(yīng)變率

最大剪切應(yīng)變率能夠反映區(qū)域構(gòu)造變形的強烈程度和能量積累的狀態(tài)，從而揭示構(gòu)造變形與地震活動之間的關(guān)系［35，44］。大多數(shù)時候，最大剪切應(yīng)變率的分布與潛在的地震危險區(qū)直接相關(guān)，其高梯度區(qū)及邊緣區(qū)域很有可能就是未來幾年的主要地震危險區(qū)［45］。前人的研究［7］也曾指出區(qū)域的最大剪應(yīng)變率分布可以反映各斷裂潛在的地震危險性，并為區(qū)域未來的地震活動性提供科學依據(jù)。

將研究區(qū)域劃分為0.1°×0.1°的網(wǎng)格，基于1991—2016年間觀測到的全球定位系統(tǒng)速度場數(shù)據(jù)［46］，采用Shen等［47］ 開發(fā)的高斯函數(shù)和Voronoi單元加權(quán)函數(shù)的應(yīng)變率插值代碼，計算了研究區(qū)域的最大剪切應(yīng)變率分布，并進一步分析其構(gòu)造變形和斷層應(yīng)變積累的特征［35］。在計算時，假設(shè)區(qū)域的構(gòu)造運動是剛性塊體連續(xù)變形，忽略地表破裂或蠕變產(chǎn)生的非彈性變形［48-49］（即地殼介質(zhì)是彈性的，應(yīng)力率和應(yīng)變率呈線性關(guān)系）。結(jié)果顯示（圖9），研究區(qū)域最大剪切應(yīng)變率的低值區(qū)位于拉脊山斷裂和積石山斷裂的交界部位，而高值區(qū)則主要位于拉脊山斷裂帶的NE、NW和SW三個方向，總體呈現(xiàn)由拉脊山構(gòu)造轉(zhuǎn)換帶向周緣遞增的趨勢。值得關(guān)注的是，積石山地震就發(fā)生在最大剪切應(yīng)變率高低值的過渡區(qū)域，這表明該地震可能是日月山右旋走滑運動和西秦嶺北緣左旋走滑運動之間的擠壓作用所致［10］，并為調(diào)節(jié)拉脊山構(gòu)造轉(zhuǎn)換帶兩側(cè)區(qū)域的應(yīng)力-應(yīng)變關(guān)系發(fā)揮著重要作用。相關(guān)研究指出，2008年于田地震［50］、2013年松原地震［51］、2014年于田地震［50］和2019年瀘定地震［52］等均發(fā)生于最大剪切應(yīng)變率高低值的過渡區(qū)域。此次積石山地震的發(fā)生也再次驗證了這一觀點，并為最大剪切應(yīng)變率高低值過渡區(qū)域的地震危險性研究提供了又一可供參考的震例。

5結(jié)論

本文以GNSS和InSAR數(shù)據(jù)為約束，反演了2023年積石山MS6.2地震的同震破裂模型，并以此計算了地震產(chǎn)生的同震庫侖應(yīng)力變化及其對周邊主要斷裂的影響。主要結(jié)論如下：

（1） 聯(lián)合反演結(jié)果顯示，積石山地震的發(fā)震斷層走向為310°，傾角為51°，是一次以逆沖為主兼有少量右旋走滑分量的地震;最大滑動大致位于5 km深度處，最大滑動量為56.02 cm，釋放的地震矩約為1.79×1018 N·m，相當于一次MW6.1地震。

（2） 積石山地震產(chǎn)生的同震庫侖應(yīng)力變化結(jié)果顯示，積石山東緣斷裂中段、拉脊山南緣斷裂東南段以及倒淌河—臨夏斷裂中東段產(chǎn)生了大幅度的應(yīng)力卸載作用，地震危險性降低;而積石山東緣斷裂南段、拉脊山南緣斷裂中南段、倒淌河—臨夏斷裂中段和西秦嶺北緣斷裂西段受到的應(yīng)力加載較高，這些斷裂的潛在地震危險性需要重點關(guān)注。

（3） 積石山地震正好發(fā)生于最大剪切應(yīng)變率的過渡區(qū)域，很有可能是日月山右旋走滑運動和西秦嶺北緣左旋走滑運動之間的擠壓作用導致的，其為調(diào)節(jié)拉脊山構(gòu)造轉(zhuǎn)換帶兩側(cè)區(qū)域的應(yīng)力-應(yīng)變關(guān)系發(fā)揮著重要作用。

參考文獻（References）

［1］李志才，陳智，武軍酈，等.基于高頻GNSS觀測的甘肅積石山6.2級地震同震形變［J/OL］.武漢大學學報（信息科學版），2024：1-13 （2024-01-15）［2024-03-26］.https：//doi.org/10.13203/j.whugis20240004.LI Zhicai，CHEN Zhi，WU Junli，et al.Co-seismic deformation of the Jishishan 6.2 earthquake in Gansu Province based on high-frequency GNSS observation［J/OL］.Geomatics and Information Science of Wuhan University，20241：1-13 （2024-01-15）［2024-03-26］.https：//doi.org/10.13203/j.whugis20240004.

［2］王安簡，高原.利用遠場記錄反演2023年甘肅積石山M6.2地震破裂過程［J］.地震，2024，44（1）：195-203.WANG Anjian，GAO Yuan.Inversion of the rupture process of the 2023 Jishishan M6.2 earthquake in Gansu by teleseismic data［J］.Earthquake，2024，44（1）：195-203.

［3］周明月，萬永革，關(guān)兆萱，等.2023年甘肅積石山MS6.2地震對周圍斷層的應(yīng)力影響及對余震的觸發(fā)作用［J/OL］.地震工程學報，2024：1-8 （2024-02-03）［2024-03-26］.https：//doi.org/10.20000/j.1000-0844.20240101001.ZHOU Mingyue，WAN Yongge，GUAN Zhaoxuan，et al.Stress change on surrounding faults and triggering of after shocks by the 2023 Jishishan MS6.2 earthquake in Gansu，China［J/OL］.China Earthquake Engineering Journal，2024：1-8 （2024-02-03）［2024-03-26］.https：//doi.org/10.20000/j.1000-0844.20240101001.

［4］陳博，宋闖，陳毅，等.2023年甘肅積石山MS6.2地震同震滑坡和建筑物損毀情況應(yīng)急識別與影響因素研究［J/OL］.武漢大學學報（信息科學版），2024：1-16 （2024-01-02）［2024-03-26］.https：//doi.org/10.13203/J.whugis20230497.CHEN Bo，SONG Chuang，CHEN Yi，et al.Emergency identification and influencing factor analysis of coseismic landslides and building damages induced by the 2023 MS6.2 Jishishan （Gansu，China） earthquake［J/OL］.Geomatics and Information Science of Wuhan University，2024：1-16 （2024-01-02）［2024-03-26］.https：//doi.org/10.13203/J.whugis20230497.

［5］白卓立，季靈運，朱良玉，等.2023年甘肅積石山MS6.2地震中川鄉(xiāng)泥流成因及破壞性分析［J/OL］.地震工程學報，2024：1-10 （2024-02-27）［2024-03-26］.https：//doi.org/10.20000/j.1000-0844.20240110001.BAI Zhuoli，JI Lingyun，ZHU Liangyu，et al.Cause and destructive analysis of the mudflow in Zhongchuan town，triggered by the 2023 Jishishan，Gansu，MS6.2 earthquake［J/OL］.China Earthquake Engineering Journal，2024：1-10 （2024-02-27）［2024-03-26］.https：//doi.org/10.20000/j.1000-0844.20240110001.

［6］張培震，鄭德文，尹功明，等.有關(guān)青藏高原東北緣晚新生代擴展與隆升的討論［J］.第四紀研究，2006，26（1）：5-13.ZHANG Peizhen，ZHENG Dewen，YIN Gongming，et al.Discussion on Late Cenozoic growth and rise of northeastern margin of the Tibetan Plateau［J］.Quaternary Sciences，2006，26（1）：5-13.

［7］周書紅，張懷，孫云強，等.青藏高原東北緣與主要斷裂帶相關(guān)的構(gòu)造應(yīng)力率和應(yīng)變率場的數(shù)值模擬［J］.中國科學（地球科學），2023，53（10）：2392-2406.ZHOU Shuhong，ZHANG Huai，SUN Yunqiang，et al.Numerical simulation of tectonic stress and strain rate fields in the northeastern margin of the Tibetan Plateau in association with major fault zones［J］.Scientia Sinica （Terrae），2023，53（10）：2392-2406.

［8］張培震，鄧起東，張國民，等.中國大陸的強震活動與活動地塊［J］.中國科學（D輯：地球科學），2003，33（增刊1）：12-20.

ZHANG Peizhen，DENG Qidong，ZHANG Guomin，et al.Strong earthquake activity and active blocks in Chinese mainland［J］.Science in China （Series D：Earth Sciences），2003，33（Suppl01）：12-20.

［9］李雨森，李為樂，許強，等.2023年積石山MS6.2級地震InSAR同震形變探測與斷層滑動分布反演［J］.成都理工大學學報（自然科學版），2024，51（1）：22-32，75.LI Yusen，LI Weile，XU Qiang，et al.Coseismic deformation and slip distribution of the 2023 Jishishan MS6.2 earthquake revealed by InSAR observations［J］.Journal of Chengdu University of Technology （Science amp; Technology Edition），2024，51（1）：22-32，75.

［10］袁道陽，張培震，雷中生，等.青海拉脊山斷裂帶新活動特征的初步研究［J］.中國地震，2005，21（1）：93-102.YUAN Daoyang，ZHANG Peizhen，LEI Zhongsheng，et al.A preliminary study on the new activity features of the Lajishan Mountain fault zone in Qinghai Province［J］.Earthquake Research in China，2005，21（1）：93-102.

［11］TANG X W，GUO R M，ZHANG Y J，et al.Rapid rupture characterization for the 2023 MS6.2 Jishishan earthquake［J］.Earthquake Research Advances，2024（2）：22-26.

［12］劉振江，韓炳權(quán)，能懿菡，等.InSAR觀測約束下的2023年甘肅積石山地震震源參數(shù)及其滑動分布［J/OL］.武漢大學學報（信息科學版），2024：1-13 （2024-01-11）［2024-03-26］.https：//doi.org/10.13203/j.whugis20240008.LIU Zhenjiang，HANG Bingquan，NAI Yihan，et al.Source parameters and slip distribution of the 2023 MW6.0 Jishishan （Gansu，China） earthquake constrained by InSAR observations［J/OL］.Geomatics and Information Science of Wuhan University，2024：1-13 （2024-01-11）［2024-03-26］.https：//doi.org/10.13203/j.whugis20240008.

［13］羅艷，朱音杰，高原.區(qū)域地震記錄揭示的2023年甘肅積石山6.2級地震震源破裂過程［J］.地震，2024，44（1）：189-194.LUO Yan，ZHU Yinjie，GAO Yuan.Source rupture process of the 2023 M6.2 Jishishan，Gansu earthquake using the regional seismic data［J］.Earthquake，2024，44（1）：189-194.

［14］楊九元，溫揚茂，許才軍.InSAR觀測揭示的2023年甘肅積石山MS6.2地震發(fā)震構(gòu)造［J/OL］.武漢大學學報（信息科學版），2024：1-11 （2024-01-11）［2024-03-26］.https：//doi.org/10.13203/j.whugis20230501.YANG Jiuyuan，WEN Yangmao，XU Caijun.Seismogenic fault structure of the 2023 Jishishan （Gansu） MS6.2 earthquake revealed by InSAR observations［J/OL］.Geomatics and Information Science of Wuhan University，2024：1-11 （2024-01-11）［2024-03-26］.https：//doi.org/10.13203/j.whugis20230501.

［15］黎朕靈，金明培，繆素秋.2021年云南漾濞MS6.4地震同震位移場和震源滑動模型反演［J］.地震研究，2021，44（3）：330-337.LI Zhenling，JIN Mingpei，MIAO Suqiu.Slip model and co-seismic displacement field of the 2021 Yangbi，Yunnan MS6.4 earthquake［J］.Journal of Seismological Research，2021，44（3）：330-337.

［16］張勇.地震破裂過程反演及其在大地震應(yīng)急中的應(yīng)用研究［D］.北京：中國地震局地球物理研究所，2010.ZHANG Yong.Rupture process inversion and its application to the earthquake emergency of large earthquakes［D］.Beijing：Institute of Geophysics，China Earthquake Administration，2010.

［17］陳威，喬學軍，劉剛，等.基于GNSS與InSAR約束的九寨溝MS7.0地震滑動模型及其庫侖應(yīng)力研究［J］.地球物理學報，2018，61（5）：2122-2132.CHEN Wei，QIAO Xuejun，LIU Gang，et al.Study on the coseismic slip model and Coulomb stress of the 2017 Jiuzhaigou MS7.0 earthquake constrained by GNSS and InSAR measurements［J］.Chinese Journal of Geophysics，2018，61（5）：2122-2132.

［18］余鵬飛，熊維，陳威，等.基于GNSS和InSAR約束的2021年瑪多MS7.4地震同震滑動分布及應(yīng)用［J］.地球物理學報，2022，65（2）：509-522.YU Pengfei，XIONG Wei，CHEN Wei，et al.Slip model of the 2021 MS7.4 Maduo earthquake constrained by GNSS and InSAR coseismic deformation［J］.Chinese Journal of Geophysics，2022，65（2）：509-522.

［19］王永哲.利用GPS和InSAR數(shù)據(jù)反演2011年日本東北MW9.0地震斷層的同震滑動分布［J］.地震學報，2015，37（5）：796-805，885.WANG Yongzhe.Coseismic slip distribution of the 2011 Tohoku MW9.0 earthquake inferred from GPS and InSAR data［J］.Acta Seismologica Sinica，2015，37（5）：796-805，885.

［20］許雙安.2022年1月8日青海門源MS6.9地震GNSS同震形變場及其斷層滑動分布［J］.地震工程學報，2023，45（2）：401-409.XU Shuang'an.Co-seismic deformation filed and fault-slip distribution of the Menyuan MS6.9 earthquake on January 8，2022，constrained by GNSS observations［J］.China Earthquake Engineering Journal，2023，45（2）：401-409.

［21］OKADA Y.Internal deformation due to shear and tensile faults in a half-space［J］.Bulletin of the Seismological Society of America，1992，82（2）：1018-1040.

［22］呂鑫，慕曉冬，張鈞，等.混沌麻雀搜索優(yōu)化算法［J］.北京航空航天大學學報，2021，47（8）：1712-1720.L Xin，MU Xiaodong，ZHANG Jun，et al.Chaos sparrow search optimization algorithm［J］.Journal of Beijing University of Aeronautics and Astronautics，2021，47（8）：1712-1720.

［23］張文朋，張春麗，高武平，等.基于PS-InSAR的薊運河斷裂震間形變監(jiān)測研究［J］.地震工程學報，2023，45（1）：120-129.ZHANG Wenpeng，ZHANG Chunli，GAO Wuping，et al.Interseismic deformation monitoring of the Jiyunhe fault based on PS-InSAR［J］.China Earthquake Engineering Journal，2023，45（1）：120-129.

［24］WANG Q，QIAO X J，LAN Q G，et al.Rupture of deep faults in the 2008 Wenchuan earthquake and uplift of the Longmen Shan［J］.Nature Geoscience，2011，4（9）：634-640.

［25］王小怡，徐克科，劉新奇，等.InSAR觀測與小震分布聯(lián)合分析2021年云南漾濞MS6.4地震應(yīng)力卸載對周圍區(qū)域的影響［J］.地震學報，2023，45（5）：849-862.WANG Xiaoyi，XU Keke，LIU Xinqi，et al.Impact of stress unloading induced by the 2021 MS6.4 Yangbi，Yunnan earthquake on the surrounding areas from the joint analysis of InSAR observations and small earthquake distribution［J］.Acta Seismologica Sinica，2023，45（5）：849-862.

［26］字城岱，滕興發(fā)，關(guān)舒丹，等.2023年塔吉克斯坦MW7.2地震InSAR同震形變場與滑動分布［J］.大地測量與地球動力學，2024，44（3）：234-239.ZI Chengdai，TENG Xingfa，GUAN Shudan，et al.Coseismic deformation field and slip distribution of the 2023 Tajikistan MW7.2 earthquake revealed by InSAR observations［J］.Journal of Geodesy and Geodynamics，2024，44（3）：234-239.

［27］楊彥明，蘇淑娟.2023年甘肅積石山MS6.2地震：一次逆沖為主的淺源強震［J］.地震，2024，44（1）：167-174.YANG Yanming，SU Shujuan.The 2023 Jishishan MS6.2 earthquake in Gansu Province，China：a shallow strong earthquake with thrust-dominated components［J］.Earthquake，2024，44（1）：167-174.

［28］王樂洋，孫龍翔，許光煜.利用GPS觀測數(shù)據(jù)反演震源參數(shù)的單純形組合加權(quán)距離灰狼新算法［J/OL］.武漢大學學報（信息科學版），2022：1-20 （2022-01-06）［2024-03-26］.https：//doi.org/10.13203/j.whugis20210114.WANG Leyang，SUN Longxiang，XU Guangyu.Combinations of the simplex and weighted distance-based grey wolf algorithms for the seismic source parameter inversion with GPS measurements［J/OL］.Geomatics and Information Science of Wuhan University，2022：1-20 （2022-01-06）［2024-03-26］.https：//doi.org/10.13203/j.whugis20210114.

［29］XUE J K，SHEN B.A novel swarm intelligence optimization approach：sparrow search algorithm［J］.Systems Science amp; Control Engineering，2020，8（1）：22-34.

［30］馬悅，宋國治，張大坤.基于改進模擬退火的三維片上網(wǎng)絡(luò)映射算法研究［J］.鄭州大學學報（理學版），2017，49（3）：9-13.MA Yue，SONG Guozhi，ZHANG Dakun.Research on mapping for three-dimensional network-on-chip based on improved simulated annealing algorithm［J］.Journal of Zhengzhou University （Natural Science Edition），2017，49（3）：9-13.

［31］杜茂康，楊單，陸聰.震后多時段應(yīng)急物資聯(lián)合配送路徑優(yōu)化研究［J］.河北科技大學學報，2023，44（6）：621-631.DU Maokang，YANG Dan，LU Cong.Research on optimization of joint delivery path for emergency materials in multiple time periods after earthquake［J］.Journal of Hebei University of Science and Technology，2023，44（6）：621-631.

［32］牛東曉，趙磊，張博，等.粒子群優(yōu)化灰色模型在負荷預測中的應(yīng)用［J］.中國管理科學，2007，15（1）：69-73.NIU Dongxiao，ZHAO Lei，ZHANG Bo，et al.The application of particle swarm optimization based grey model to power load forecasting［J］.Chinese Journal of Management Science，2007，15（1）：69-73.

［33］XIAO L H，ZHENG R，ZOU R.Coseismic slip distribution of the 2021 MW7.4 Maduo，Qinghai earthquake estimated from InSAR and GPS measurements［J］.Journal of Earth Science，2022，33（4）：885-891.

［34］蔣凱，楊瑩輝，許強，等.基于InSAR形變的滑面摩擦殘余估計與滑坡危險性評價［J/OL］.工程地質(zhì)學報，2024：1-16 （2024-03-15）［2024-03-26］.https：//doi.org/10.13544/j.cnki.jeg.2023-0245.JIANG Kai，YANG Yinghui，XU Qiang，et al.Friction residual estimation and landslide risk assessment based on InSAR deformation［J/OL］.Journal of Engineering Geology，2024：1-16 （2024-03-15）［2024-03-26］.https：//doi.org/10.13544/j.cnki.jeg.2023-0245.

［35］蔣鋒云，季靈運，朱良玉，等.聯(lián)合GPS和InSAR研究海原—六盤山斷裂現(xiàn)今的地殼變形特征［J］.地震地質(zhì)，2023，45（2）：377-400.JIANG Fengyun，JI Lingyun，ZHU Liangyu，et al.The present crustal deformation characteristics of the Haiyuan—Liupanshan fault zone from InSAR and GPS observations［J］.Seismology and Geology，2023，45（2）：377-400.

［36］KING G C P，STEIN R S，LIN J.Static stress changes and the triggering of earthquakes［J］.Bulletin of the Seismological Society of America，1994，84（3）：935-953.

［37］龔煒程，孫云強，邱鑫鵬，等.青藏高原東北緣主要斷裂帶的庫侖應(yīng)力演化及其強震危險性［J］.地震工程學報，2023，45（3）：585-597.GONG Weicheng，SUN Yunqiang，QIU Xinpeng，et al.Coulomb stress evolution and seismic hazards along major faults in the northeastern margin of the Tibetan Plateau［J］.China Earthquake Engineering Journal，2023，45（3）：585-597.

［38］朱俊文，蔡藝萌，張波，等.昌馬MS7.6地震三維同震形變場數(shù)值模擬及庫侖應(yīng)力觸發(fā)研究［J］.地震工程學報，2023，45（3）：598-605.ZHU Junwen，CAI Yimeng，ZHANG Bo，et al.Numerical simulation of 3D coseismic deformation field and Coulomb stress triggering Changma MS7.6 earthquake［J］.China Earthquake Engineering Journal，2023，45（3）：598-605.

［39］MELOSH H J，RAEFSKY A.A simple and efficient method for introducing faults into finite element computations［J］.Bulletin of the Seismological Society of America，1981，71（5）：1391-1400.

［40］孫云強，羅綱，黃祿淵.基于分裂節(jié)點法的地震同震和震后形變數(shù)值模擬及其在汶川大地震中的應(yīng)用［J］.地質(zhì)力學學報，2021，27（2）：241-253.SUN Yunqiang，LUO Gang，HUANG Luyuan.Numerical simulation of coseismic and postseismic deformation through a node-splitting algorithm：a case study of the Wenchuan earthquake［J］.Journal of Geomechanics，2021，27（2）：241-253.

［41］岳沖，屈春燕，牛安福，等.瑪多MS7.4地震對周邊斷層的應(yīng)力影響分析［J］.地震地質(zhì)，2021，43（5）：1041-1059.YUE Chong，QU Chunyan，NIU Anfu，et al.Analysis of stress influence of Qinghai Maduo MS7.4 earthquake on surrounding faults［J］.Seismology and Geology，2021，43（5）：1041-1059.

［42］郭曉玉，高銳，高建榮，等.青藏高原東北緣馬銜山斷裂帶構(gòu)造屬性的綜合研究［J］.地球物理學報，2018，61（2）：560-569.GUO Xiaoyu，GAO Rui，GAO Jianrong，et al.Integrated analysis on the tectonic features of the Maxian Shan fault zone in the northeastern Tibetan Plateau［J］.Chinese Journal of Geophysics，2018，61（2）：560-569.

［43］李光濤，程理，吳昊，等.臨潭—宕昌主干斷裂南東段晚第四紀活動的地質(zhì)地貌證據(jù)［J］.地震工程學報，2020，42（2）：376-383.LI Guangtao，CHENG Li，WU Hao，et al.Geological and geomorphological evidence of Late Quaternary activity along the southeastern segment of the Lintan—Tanchang major fault，China［J］.China Earthquake Engineering Journal，2020，42（2）：376-383.

［44］金陽，金紅林，高原，等.基于GNSS數(shù)據(jù)的三種應(yīng)變計算方法在青藏高原東南緣地區(qū)的應(yīng)用［J］.地震，2021，41（4）：57-67.JIN Yang，JIN Honglin，GAO Yuan，et al.Application of three strain calculation methods based on GNSS data in the southeastern margin of Tibet Plateau［J］.Earthquake，2021，41（4）：57-67.

［45］王伶俐，洪敏，張勇，等.云南地區(qū)GNSS應(yīng)變場時空演化與地震事件關(guān)聯(lián)分析［J］.中國地震，2020，36（1）：91-104.WANG Lingli，HONG Min，ZHANG Yong，et al.Analysis of the relationship between the spatiotemporal evolution of GNSS strain field and earthquake events in Yunnan region［J］.Earthquake Research in China，2020，36（1）：91-104.

［46］WANG M，SHEN Z K.Present-day crustal deformation of continental China derived from GPS and its tectonic implications［J］.Journal of Geophysical Research：Solid Earth，2020，125（2）：e2019JB018774.

［47］SHEN Z K，WANG M，ZENG Y H，et al.Optimal interpolation of spatially discretized geodetic data［J］.The Bulletin of the Seismological Society of America，2015，105（4）：2117-2127.

［48］TAPPONNIER P，PELTZER G，LE DAIN A Y，et al.Propagating extrusion tectonics in Asia：new insights from simple experiments with plasticine［J］.Geology，1982，10（12）：611-616.

［49］PELTZER G，TAPPONNIER P.Formation and evolution of strike-slip faults，rifts，and basins during the India—Asia Collision：an experimental approach［J］.Journal of Geophysical Research：Solid Earth，1988，93（B12）：15085-15117.

［50］魏文薪，武艷強，江在森，等.2014年新疆于田MS7.3地震前地殼變形特征研究［J］.地震學報，2014，36（4）：595-601.WEI Wenxin，WU Yanqiang，JIANG Zaisen，et al.Study on crustal deformation characteristics before Xinjiang Yutian MS7.3 earthquake in 2014［J］.Acta Seismologica Sinica，2014，36（4）：595-601.

［51］楊博，朱爽，周海濤，等.2013年松原MS5.8地震的地殼水平形變背景［J］.地震，2015，35（2）：72-79.YANG Bo，ZHU Shuang，ZHOU Haitao，et al.Horizontal crust deformation field before the 2013 Songyuan MS5.8 earthquake in Jilin Province［J］.Earthquake，2015，35（2）：72-79.

［52］李瑩，田建慧，李心怡，等.2022年9月5日四川瀘定MS6.8地震深部構(gòu)造特征［J］.地球物理學報，2023，66（4）：1385-1396.LI Ying，TIAN Jianhui，LI Xinyi，et al.Deep tectonic pattern of the Luding MS6.8 earthquake on 5th September 2022 in Sichuan Province，China［J］.Chinese Journal of Geophysics，2023，66（4）：1385-1396.

（本文編輯：趙乘程）