基于背景噪聲成像的沂沭斷裂帶及鄰區波速變化

中圖分類號：P631.4 文獻標志碼：A

Abstract： The Yishu fault zone is located in the middle section of the Tanlu fault zone，consisting of multiple nearly north-northeast trending faults. Significant variations in crustal structure are present on both sides of the fault zone. The seismic velocity and velocity variations are important parameters for understanding the deep structural characteristics of the earth and the potential damage of the medium. Based on the continuous seismic waveforms from permanent stations in the Yishu fault zone and surrounding regions from January 20ll to December 20l4，we perform ambient noise tomography to monitor the phase velocity variations of Rayleigh wave. The results show that the velocity variations correlate well with the small earthquake activities. The earthquakes （ ??Mgt;2 ）are mainly concentrated near the velocity reduction and velocity transition zones.The velocity variations near the Yishu fault zone are not significant，indicating that the fault activity is relatively weak recently. The velocity variations in the Jiaobei uplift are notable between 2Ol3 and 2014，which is related to surface deformation.It is speculated that the compression effect of regional stress is the main factor for the velocity increase in the Jiaobei uplift.

Key words： Yishu fault zone； ambient noise tomography； velocity variation； small earthquakeactivity

0 引言

地震波速測量是監測地殼介質變化的重要技術，在探究地球內部構造與動力學方面有著十分重要的作用[1-2]。大量研究表明[3-7]，火山和地震等構造因素以及降水、氣溫和大氣壓強等非構造因素均會導致地下介質地震波速變化。因此，監測地震波速變化有助于認識地下介質的應力變化狀態，為區域地震活動規律和危險性判斷提供重要的地震學證據。

目前，地下介質地震波速變化研究主要涉及天然重復地震、人工震源和背景噪聲監測等方面。天然重復地震主要利用發生位置幾乎相同、震源機制和波形相近的重復地震事件，通過移動窗互譜或移動窗互相關等方法分析重復地震事件之間的波形差異[8-11]，進而獲取地震波速度的變化特征。人工震源具有地震波形重復性較強的優點，結合尾波干涉方法可以獲得精度較高的波速變化結果[3，12]。然而，天然重復地震監測經常受到發震時間不連續和時空分布不均勻等因素的影響；而人工震源監測工作量較大且觀測成本較高，主要用于監測地震波速的時序變化特征。為了獲得地震波速變化的空間特征，Pei等[13利用2000—2014年期間的直達P波走時信息對龍門山斷裂帶地區開展了區域性波速變化研究，揭示了汶川地震前的波速降低以及震后波速的恢復過程；但是這一方法受區域震源分布不均的影響較大。隨著背景噪聲處理技術的發展，許多學者[14-20]通過背景噪聲互相關或自相關來重建經驗格林函數，并將不同時段疊加的經驗格林函數等效為重復震源來監測地下介質波速變化，已在火山、強震和工程地質等領域開展了波速變化時序監測研究。背景噪聲成像具有可重復性和時間穩定性等特點[21]，可以不依賴于地震或人工爆破源就能獲取精細的地下速度結構及其變化特征，有助于了解地殼介質的應力變化狀態。

沂沭斷裂帶位于郯廬斷裂帶山東段，北起萊州灣，向南沿沂河、沭河延伸至蘇北駱馬湖地區，由多條近于平行的北北東向斷裂組成[22]。斷裂帶兩側區域構造復雜，發育魯西隆起、濟陽凹陷、濟寧凹陷、膠北隆起和膠萊盆地等多個構造單元。歷史上，受太平洋板塊俯沖作用影響，沂述斷裂帶發生過強烈的左旋剪切作用[23-25]。近年來，該區域的地殼形變和構造活動仍十分頻繁，被認為是郯廬斷裂帶上地震危險性最高的區域之一[26-28]。因此，沂述斷裂帶及鄰區地震波速變化研究對于了解該區域的應力變化及地震風險性至關重要。本研究采用背景噪聲成像技術構建不同時段沂述斷裂帶及鄰區的瑞利面波相速度結構，通過分析波速變化的時空特征，探討地震波速變化與不同因素的相關性以及可能的影響機制，以期揭示沂述斷裂帶深部介質的演化特征。

1 數據與方法

選取2011年1月—2014年12月布設于山東和河南臺網固定臺站的垂直分量連續波形數據，其中山東臺網39個臺站，河南臺網4個臺站，共計43個臺站（圖1）。理論和實際應用均表明，長時間噪聲互相關函數疊加可以獲得臺站對的經驗格林函數[29-32]。本文數據預處理參照 Bensen 等[33]的處理流程。首先，將各臺站的原始數據截成時間長度為1d的波形數據；然后，對數據進行降采樣、去儀器響應、去均值、去傾斜分量、帶通濾波 （2～40s） ，以及基于滑動絕對平均方法的時域歸一化和頻譜白噪化等預處理；最后，對臺站對波形進行互相關計算和疊加，獲得臺站對的互相關函數。

利用圖形轉換技術提取瑞利面波相速度頻散曲線[34-36]。考慮到實際數據的信噪比及頻散曲線的連續性，選取臺間距大于2倍波長且信噪比大于8的波形用于頻散曲線的提取。圖2展示了同一臺站對（新鄉-安丘）不同年度的相速度頻散曲線。拾取過程中，本文均采用完全相同的參數和篩選原則。值得注意的是，盡管這些數據來自不同時期的觀測，但各年度的頻散曲線呈現出非常高的相似性。最終，我們獲取了不同年度 5～35 s周期范圍內、圖1所示所有地震臺站兩兩成對臺站對的頻散曲線（圖3）。從圖3中可以看出，不同年度間的頻散曲線呈現出較高的相似性，表明研究區域的地下介質結構在觀測期間保持了較好的穩定性。

本研究采用基于射線理論的相速度空間反演方法[37-38]反演區域瑞利面波相速度結構（圖4）。該方法采用最小二乘準則，通過不斷迭代使目標函數值達到最小。目標函數表達式為

Φ（s）=（t-t_obs^T）C_D^-1（t-t_obs）+

（s-s_p）^TC_M^-1（s-s_p）_e

式中： Ψ_tΨ_t 和 t_obs 分別為預測和觀測走時； s 和 s_p 分別代

圖1沂述斷裂帶區域地質構造與臺站分布圖

XX：新鄉臺;TIA：泰安臺；ANQ：安丘臺;YTA：煙臺臺;SDTA：泰安GPS觀測點;SDCY：昌邑GPS觀測點;SDYT：煙臺GPS觀測點。

圖2新鄉（XX）-安丘（ANQ）臺站對不同年度的相速度頻散曲線

Fig.2Phase velocity dispersion curves of the Xinxiang （XX）-Anqiu （ANQ） station pair in different years

圖3所有臺站對不同年度周期 的瑞利面波相速度頻散曲線

Fig.3Rayleigh wave phase velocity dispersion curves of allstation pairs with a period of 5-35 s in different year

表慢度模型和先驗慢度模型； c_D 為觀測數據的協方差矩陣； C_M 為反演模型參數的先驗協方差矩陣。網格點 r₁ 和 r₂ 之間的協方差滿足以下公式：

式中： σ_s 為先驗慢度的不確定度； L 表示模型的相關長度。

2結果

不同周期的瑞利面波相速度可以反映不同深度范圍的剪切波速度結構特征。圖4展示了基于不同年度頻散曲線反演獲得的沂沐斷裂帶及其鄰近區域的相速度結構分布。由于在頻散曲線拾取和相速度結構反演過程中均采用了相同的篩選原則和參數，因此不同年度的瑞利面波相速度結果呈現出非常高的相似性。周期為10s的相速度結果（圖 4a，d，g 、j）主要反映了上地殼的速度結構特征，與地表地質單元具有較好的一致性：濟陽凹陷、濟寧凹陷和膠萊盆地等區域表現為低波速異常，而魯西隆起和膠北隆起等區域主要表現為高波速異常。周期為18s瑞利面波相速度結果（圖 4b，e，h，k） 主要反映了中下地殼的結構特征，魯西隆起和膠北隆起區域主要呈現出高波速異常。周期為 25s 瑞利面波相速度結果（圖4c、f、i、1）主要反映了上地幔頂部的結構特征，沂沭斷裂帶附近主要呈現出低波速異常，而濟陽凹陷區域則表現為高波速異常。

為了分析沂述斷裂帶及鄰區不同年度之間的波速差異，我們通過差分計算獲得相鄰年份瑞利面波相速度隨時間的變化特征（圖5）。本文主要考慮不確定度小于 1.5% 的區域（圖5d）。為了更清晰地對比小震活動（ 2～5 級）與波速變化的關系，采用雙差定位方法對區域的小震進行了重定位。由于該區.瑞利面波相速度均值； Δv .瑞利面波相速度變化量。

圖4不同年度周期為 的瑞利面波相速度結構

Fig.4Phase velocity structures of Rayleigh wave at periods of 10，18，and 25s in different yea

域的小震發震深度主要集中于 10km 左右，與周期10s的瑞利面波相速度敏感核深度比較相近，而且周期10s的瑞利面波相速度結構與地表構造也具有較好的相關性，可靠性較高（圖5d）；因此，本文主要探討周期10s的瑞利面波相速度結構變化特征與小震活動的相關性及可能的影響因素。由圖5可以看出，波速變化呈現出強烈的不均勻性，而小震主要發生于波速減小以及波速變化的梯度帶附近。2011—2012年期間，膠北隆起和魯西隆起西部主要表現為波速降低，而膠萊盆地和魯西隆起東部表現為波速升高（圖5a），2012—2013年期間的波速變化（圖5b）整體與 2011—2012年波速變化比較相似，但2013—2014年期間膠萊盆地和膠北隆起波速變化幅度較大，形成鮮明對比（圖5c）。

3 討論

周期10s的相速度結構與區域地表構造密切相關。魯西隆起和膠北隆起區表現出相對高的波速異常，而膠萊盆地區域則表現為低波速異常，與前人噪聲成像結果基本一致[39]。不同時段的相速度結構表現出的高度相似性（圖4）也印證了背景噪聲成像的穩定性和可重復性。在波速差分變化結果中，膠北隆起區域在2013—2014年間存在較明顯的波速升高現象（圖5c）。為了驗證該波速變化異常的可靠性，選取泰安（TIA）和煙臺（YTA）兩臺站間不同年度的經驗格林函數進行對比，結果如圖6所示。可以看出，與2011—2013年相比，2014年的經驗格林函數存在明顯的相位提前，表明兩臺站之間存在波速升高的區域，這可能與我們觀測到的膠北隆起波速升高的現象相對應。另外，沂述斷裂帶上的波速差分變化幅度整體較低（圖5a—c），表明近期沂沐斷裂帶活動性較弱[40]

圖5周期10s的瑞利面波相速度差分變化（a、b、c）與不確定度分析（d）

Fig.5 Variation （a，b，c） and uncertainty analysis （d）of phase velocity difference of Rayleigh wave at period of

周期10s的瑞利面波相速度結構年度差分變化特征與2級以上小震活動具有明顯的相關性。小震主要發生于波速變化減小以及波速變化梯度帶附近（圖5）。實驗表明，巖石應力與彈性波速度存在相關性，巖石會隨應力增加而出現裂隙、擴容和壓實現象[41]。人工源地震[42]也證實了地震波速度與地應力場的相關性，主要表現為地震波速隨地應力的增大而增大。研究區內波速變化減小的區域可能代表了介質的壓實作用減弱，不穩定性上升，較易發生小震活動[43-45]；而在波速變化梯度帶附近，應力的差異性易在斷層內產生剪切作用，進而導致小震活動增加。

地震波速常常受區域降水和大氣壓強等因素的影響[12，15，46-48]。區域降水會影響地下水的補給條件，改變地下介質的孔隙水壓力，減緩淺層地質體的孔隙彈性效應，從而引起地震波速度的下降[49-50]而大氣壓強會在地表施加荷載，改變介質內部的裂隙狀態，進而引起地震波速度的變化。但是，大量研究表明，降水和大氣壓強主要對地表淺層波速產生影響，所導致的波速變化大概為0.015%[4，15，48，51]，遠小于2013—2014年期間膠東地區所觀測到的波速變化幅度（圖5c）。

研究表明[52-56]，地表形變可以反映區域的構造活動特征，與地震波速變化存在相關性。我們分析了泰安（SDTA）、煙臺（SDYT）和昌邑（SDCY）GPS觀測點的GPS基線解算結果。圖7為最小二乘法擬合的東西方向上的基線值及其基線變化速率，可以看出：在東西方向上，2011—2014年期間泰安-昌邑臺站基線變化速率由 3. 19mm/a 轉變為-0.61mm/a ，表明兩臺站之間的拉張趨勢逐漸減緩；而泰安-煙臺臺站基線在2011年處于拉張狀態，2012—2014年期間變化速率轉為負值，且擠壓速率逐漸增大，表明泰安-煙臺之間的擠壓強度遂漸增加，暗示著應力不斷積累，這與我們觀測到的膠北隆起波速升高的現象一致。

圖6泰安（TIA）-煙臺（YTA）臺站對經驗格林函數變化 Fig.6 Empirical Green'sfunctionvariations ofTaian（TIA） andYantai（YTA）stationpair

圖7不同GPS觀測點之間的東西方向基線時間序列Fig.7 Timeseriesofeast-westbaselinebetweendifferentGPSstations

4結論

本研究基于2011—2014年期間沂沐斷裂帶及鄰區43個地震臺站的垂直分量連續波形數據，利用背景噪聲成像方法反演得到了不同年度的瑞利面波相速度變化結果，主要獲得以下認識：

1）短周期瑞利面波相速度結構與地表構造密切相關，隆起區主要表現為高波速異常，而凹陷區則表現為低波速異常。

2）波速變化與小震活動具有明顯的相關性，小震主要發生于波速變化減小以及波速變化的梯度帶附近，可能與區域的應力調整有關。

3）膠北隆起2013—2014年期間波速變化尤為顯著，且與GPS基線變化結果相一致，表明區域的擠壓應力增強可能是導致膠北隆起地區地震波速顯著升高的主要因素。

參考文獻（References）：

[1]Wang Q Y，Yao H J. Monitoring of Velocity ChangesBased on Seismic Ambient Noise：A Brief Review andPerspective[J]. Earth and Planetary Physics，2020，4（5）：532-542.

[2]Hutapea F L，Tsuji T，Ikeda T. Real-Time CrustalMonitoring System of Japanese Islands Based onSpatio-Temporal Seismic Velocity Variation[J]. EarthPlanets Space，2020，72（1）：1-16.

[3] Niu FL，Silver PG，Daley T M，et al. PreseismicVelocity Changes Observed from Active SourceMonitoring at the Parkfield SAFOD Drill Site[J].Nature，2008，454：204-208.

[4]Wang Q Y，Brenguier F，Campillo M，et al. SeasonalCrustal Seismic Velocity Changes Throughout Japan[J].Journal of Geophysical Research： Solid Earth，2017，122（10）：7987-8002.

[5]Cloetingh S，Sternai P，Koptev K，et al. CoupledSurface to Deep Earth Processes：Perspectives fromTOPO-EUROPE with an Emphasis on Climate- andEnergy-Related Societal Challenges[J]. Global andPlanetary Change，2023，226：104140.

[6]Ermert L A，Cano E C，Chaussard E，et al. ProbingEnvironmental and Tectonic Changes UnderneathCiudad de México with the Urban Seismic Field[J].Solid Earth，2023，14（5）：529-549.

[7]Fan L F， Sun H Y. Seismic Wave PropagationThrough an In-Situ Stressed Rock Mass[J]. Journal ofApplied Geophysics，2015，121：13-20.

[8]Ratdomopurbo A，Poupinet G. Monitoring a TemporalChange of Seismic Velocity in a Volcano： Applicationto the 1992 Eruption of Mt Merapi （Indonesia）[J].Geophysical Research Letters，1995，22（7）：775-778.

L9」Pandolti D， Bean C J， Saccorotti G. Coda WaveInterferometric Detection of Seismic Velocity ChangesAssociated with the 1999 M=3 .6 Event at MtVesuvius[J]. Geophysical Research Letters， 2006，33（6）：272-288.

[10]Martini F，Bean C J，Saccorotti G，et al. SeasonalCycles of Seismic Velocity Variations Detected UsingCoda Wave Interferometry at Fogo Volcano， SaoMiguel，Azores，During 2003 - 2004[J]. Journal ofVolcanology and Geothermal Research， 2009，181（3/4）：231-246.

[11]Battaglia J，Metaxian JP，Garaebiti E. Earthquake-Volcano InteractionImagedbyCoda WaveInterferometry[J]. Geophysical Research Letters，2012，39（11）：L11309.

[12]Silver PG，Daley T M，Niu F L，et al. Active SourceMonitoring of Cross-Well Seismic Travel Time forStress-InducedChanges [J]. BulletinoftheSeismological Society of America，2007，97（1B）：281-293.

[13]Pei SP，Niu F L，Ben-Zion Y，et al. Seismic VelocityReduction andAcceleratedRecovery duetoEarthquakes on the Longmenshan Fault [J]. NatureGeoscience，2019，12：387-392.

[14]Campillo M，Paul A.Long-Range Correlations in theDiffuse Seismic Coda[J]. Science，2003，299：547 -549.

[15]Shapiro N M，Campillo M，Stehly L，et al. High-Resolution Surface Wave Tomography from AmbientSeismic Noise[J]. Science，2005，307：1615 - 1618.

[16]Sens-Schonfelder C， Wegler U. Passve ImageInterferometry and Seasonal Variations of SeismicVelocities atMerapiVolcano，Indonesia [J].Geophysical Research Letters，2006，33（21）：1-5.

[17]Brenguier F，Shapiro N M，Campillo M，et al.TowardsForecasting Volcanic EruptionsUsingSeismic Noise[J].Nature Geoscience，20o8，1（2）：126 -130.

[18]劉志坤，黃金莉.利用背景噪聲互相關研究汶川地震震源區地震波速度變化[J].地球物理學報，2010，53（4）：853-863.Liu Zhikun， Huang Jinli. Temporal Changes ofSeismic Velocity Around the Wenchuan EarthquakeFault Zone from Ambient Seismic Noise Correlation[J].Chinese Journal of Geophysics，2O1o，53（4）：853-863.

「10］安茄 王佳濤 楊微 笙 利田背昱噪聲研空 2021年An Yanru，Wang Weitao，Yang Wei，et al.UsingAmbient Noise to Study the Co-Seismic and Post-Seismic Velocity Changes of the 2O21 Yangbi M_S6.4 （20Earthquake in Yunnan [J]. Chinese Journal ofGeophysics，2023，66（8）：3185-3201.

[20]孔慶翰，呂子強，張廣偉.基于密集臺陣的遼寧老虎臺礦區微震事件重定位［J].吉林大學學報（地球科學版），2024，54（5）：1685-1695.Kong Qinghan， Lu Ziqiang， Zhang Guangwei.Microseismic Events Relocation of Laohutai CoalMine in Liaoning Province Based on Dense SeismicArray[J]. Journal of Jilin University （Earth ScienceEdition），2024，54（5）：1685-1695.

[21]王智雄，馮晅，侯賀晟，等.基于背景噪聲成像方法的蘿北地區淺地表速度結構［J]．世界地質，2022，41（2）：365-372.Wang Zhixiong，Feng Xuan，Hou Hesheng，et al.Shallow-Surface Velocity Structure in Luobei AreaBased on Ambient Noise Tomography Method[J].World Geology，2022，41（2）： 365-372.

[22]劉洪波，吳治國，鄭孝誠，等.沂沭斷裂帶萊州灣海域重力場特征[J].吉林大學學報（地球科學版），2019，49（5）：1438 - 1447.Liu Hongbo，Wu Zhiguo， Zheng Xiaocheng，et al.Gravity Field Characteristics of Yishu Fault Zone inLaizhou Bay[J]. Journal of Jilin University （EarthScience Edition），2019，49（5）：1438-1447.

[23]張鵬，王良書，鐘鍇，等.郯廬斷裂帶的分段性研究[J].地質論評，2007，53（5）：586－591，721，722.Zhang Peng，Wang Liangshu， Zhong Kai， et al.Research on the Segmentation of Tancheng-LujiangFault Zone[J]. Geological Review，2007，53（5）：586-591，721，722.

[24]朱成林，甘衛軍，李杰，等．日本 M_w9.0 地震后沂述斷裂帶兩側塊體相對運動及對地震活動的影響［J].地球物理學報，2018，61（3）：988－999.Zhu Chenglin，Gan Weijun，Li Jie，et al.RelativeMotion Between the Two Blocks on Either Side of theYishu Fault Zone After the 2Oll Japan M_w9.0 Earthquake and Its Effect on Seismic Activity[J].Chinese Journal of Geophysics， 2018，61（3）： 988 -999.

[25]晁洪太，李家靈，崔昭文，等.郯廬斷裂帶中段全新世活斷層的特征滑動行為與特征地震［J]．內陸地震，1994，8（4）：297- 304.Characteristic Earthquakes [J]. Inland Earthquake，1994，8（4）： 297-304.

[26]嚴樂佳，朱光，林少澤，等.沂沭斷裂帶新構造活動規律與機制[J]．中國科學：地球科學，2014，44（7）：1452-1467.Yan Lejia， Zhu Guang，Lin Shaoze，et al. The Lawand Mechanism of Neotectonic Activity in the YishuFault Zone[J]. Scientia Sinica （Terrae），2O14，44（7）：1452-1467.

[27]李兵，謝富仁，黃金水，等.龍門山斷裂帶大邑地震空區地應力狀態與地震危險性[J]．中國科學：地球科學，2022，52（7）：1409-1418.Li Bing，Xie Furen，Huang Jinshui，et al. In-SituStress State and Seismic Hazard of Dayi SeismicEmpty Area in Longmenshan Fault Zone[J]. ScientiaSinica（Terrae），2022，52（7）：1409-1418.

[28]鐘大賚，李理，丁增勇，等.魯西北西向斷裂系與沂述斷裂帶晚中生代演化關系及其動力學背景探討[J].地學前緣，2010，17（3）：166－190.Zhong Dalai， Li Li， Ding Zengyong， etal.Relationship Between NW Faults of West Shandongand Yi-Shu Fault Zone in Late Mesozoic and TheirGeotectonic Setting Operations [J]. Earth ScienceFrontiers，2010，17（3）：166-190.

[29]呂子強，雷建設.2015年尼泊爾 Ms8.1地震震源區S波三維速度結構與強震發生機理研究[J].地球物理學報，2016，59（12）：4529-4543.Lu Ziqiang，Lei Jianshe. 3 - D S- Wave VelocityStructure Around the 2O15 Ms8.1 Nepal EarthquakeSource Areas and Strong Earthquake Mechanism[J].Chinese Journal of Geophysics，2016，59（12）： 4529-4543.

[30]呂子強，趙俐紅，李鉑，等.天山造山帶地區瑞利面波相速度與方位各向異性[J].地球物理學報，2019，62（9）：3354 - 3364.Lu Ziqiang，Zhao Lihong，Li Bo，et al. RayleighWave Phase Velocity and Azimuthal Anisotropy ofTien Shan Orogenic Belt from Ambient NoiseTomography[J]. Chinese Journal of Geophysics，2019，62（9）：3354-3364.

[31]夏心茹，江國明，趙大鵬.郯廬斷裂中南段地殼速度結構與地震活動性研究［J].CT理論與應用研究，2025，34（2）：163-174.Xia Xinru， Jiang Guoming， Zhao Dapeng. CrustalVelocity Structure of the Middle-Southern Segment ofthe Tanlu Fault Zone（China） and Its Correlationwith the Seismic Activity in the Region[J].Computerized Tomography Theory and Applications，2025，34（2）：163-174.

[32]王娟娟，姚華建，王偉濤，等.基于背景噪聲成像方法的新疆呼圖壁儲氣庫地區近地表速度結構研究［J].地球物理學報，2018，61（11）：4436－4447.Wang Juanjuan，Yao Huajian，Wang Weitao，et al.Study of the Near-Surface Velocity Structure of theHutubi Gas Storage Area in Xinjiang from AmbientNoise Tomography [J]. ChineseJournalofGeophysics，2018，61（11）： 4436-4447.

[33]Bensen G D，Ritzwoller M H，Barmin M P，et al.Processing Seismic Ambient Noise Data to ObtainReliable Broad-BandSurfaceWaveDispersionMeasurements[J ].GeophysicalJournalInternational，2007，169（3）：1239－1260.

[34]Yao H J，Hilst R D，Hoop M V. Surface-WaveArray Tomography in SE Tibet from AmbientSeismic Noise and Two-Station Analysis： I： PhaseVelocity Maps[J]. Geophysical Journal International，2006，166（2）： 732-744.

[35]呂子強，雷建設，周智剛，等.環渤海地區Pn波速度結構與各向異性[J].地球物理學報，2016，59（6）：2047-2055.Lü Ziqiang，Lei Jianshe， Zhou Zhigang，et al. Pn-Wave Velocity and Anisotropy Around the Bohai SeaAreas[J]. Chinese Journal of Geophysics， 2O16，59（6）：2047-2055.

[36]張浩，雷建設，宋曉燕，等．山西斷陷帶及其鄰區背景噪聲面波直接反演成像[J].CT 理論與應用研究，2025，34（2）：175 -189.Zhang Hao，Lei Jianshe，Song Xiaoyan，et al. DirectSurface-Wave Tomography from Ambient Noise inthe Shanxi Rift Zone and Adjacent Areas[J].Computerized Tomography Theory and Applications，2025，34（2）：175-189.

[37]Tarantola A， Valette B. Generalized’ NonlinearInverse Problem Solved Using the Least SquaresCriterion[J]. Reviews of Geophysics，1982，20（2）：219-232.

[38]Tarantola A， Nercessian A. Three-DimensionalInversion Without Blocks[J]. Geophysical JournalInternational，1984，76（2）： 299-306.

[39]顧勤平，丁志峰，康清清，等.郯廬斷裂帶中南段及鄰區其干背昱噪聲的瑞利波群速度層析成像「I地球Gu Qinping，Ding Zhifeng，Kang Qingqing，et al.Group Velocity Tomography of Rayleigh Wave in theMiddle-Southern Segment of the Tan-Lu Fault Zoneand Adjacent Regions Using Ambient Seismic Noise[J]. Chinese Journal of Geophysics，2020，63（4）：1505-1522.

[40]杜存鵬，殷海濤，于勝文，等.基于GNSS 及地震活動的沂述斷裂帶危險性分析[J].大地測量與地球動力學，2023，43（1）：38-41.Du Cunpeng，Yin Haitao，Yu Shengwen，et al. RiskAnalysis of Yishu Fault Zone Based on GNSS andSeismicity[J]. Journal of Geodesy and Geodynamics，2023，43（1）： 38 -41.

[41]張逸凌，劉雷，杜建國.高溫高壓下巖石波速研究進展[J].地球物理學進展，2023，38（5）：1999－2022.Zhang Yiling， Liu Lei， Du Jianguo.Research Progressof Rock Wave Velocity Under High Temperature andHigh Pressure[J]. Progress in Geophysics， 2023，38（5）：1999 -2022.

[42]Furumoto M，Ichimori Y，Hayashi N. Seismic WaveVelocity Changes and Stress Build-up in the Ccrust ofthe Kanto-Tokai Region[J]. Geophysical ResearchLetters，2001，28（19）：3737-3740.

[43]祁玉萍，孟令媛，龍鋒，等.錦屏一級水電站庫區中小地震震源機制解與局部應力場特征［J].地球物理學報，2024，67（1）：189-204.Qi Yuping， Meng Lingyuan， Long Feng， et al.Characteristics of Focal Mechanism Solution andLocal Stress Field of Small to Moderate Earthquakesin Jinping Hydropower Station Reservoir Area[J].Chinese Journal of Geophysics，2024，67（1）：189 -204.

[44]郭云飛，雷建設，關鵬虎.海南島及鄰區中小地震重定位[J].CT理論與應用研究，2025，34（2）：191－204.Guo Yunfei，Lei Jianshe，Guan Penghu. Relocation ofModerate-to-Small Earthquakesin and AroundHainan Island[J]. Computerized Tomography Theoryand Applications，2025，34（2）：191-204.

[45]王天哲，張萬佶，祁善博，等.長短期記憶網絡在P波初至震相識別中的實驗研究[J].CT理論與應用研究，2025，34（2）：205-215.Wang Tianzhe， Zhang Wanji，Qi Shanbo，et al.Experimental Study on Long Short-Term MemoryNetworks for Identifying P-Wave Primary Phase[J].Computerized Tomography Theory and Applications ，2025，34（2）：205-215.Variations at TCDP are Controlled by MJO DrivenPrecipitation Pattern and High Fluid DischargeProperties[J]. Earth and Planetary Science Letters，2014，391（2）：121-127.

[47]Delouche E， Stehly L. Seasonal Seismic VelocityVariationsMeasuredUsingSeismicNoiseAutocorrelations to Monitor the Dynamic of Aquifersin Greece[J]. Journal of Geophysical Research： SolidEarth，2023，128（12）：B026759.

[48]Gradon C，Brenguier F，Stammeijer J，et al. SeismicVelocity Response to Atmospheric Pressure UsingTime-Lapse Passive Seismic Interferometry [J].Bulletin of the Seismological Society of America，2021，111（6）：3451-3458.

[49]Gret A，Snieder R， Ozbay U. Monitoringin SitustressChanges in a Mining Environment with Coda WaveInterferometry[J]. Geophysical Journal International，2006，167（2）：504-508.

[50]顧悅，裴爍瑾，梁姍姍，等.基于地震背景噪聲監測山東臨沂地區地下水位變化[J].中國地震，2021，37（4） ：857 -867.Gu Yue， Pei Shuojin， Liang Shanshan， et al.Monitoring of Groundwater Level Change by SeismicAmbient Method in Linyi Area of Shandong Province[J].Earthquake Research in China，2021，37（4）：857 -867.

[51]靳宇蓉，魯克新，李鵬，等.基于穩定同位素的土壤水分運動特征[J].土壤學報，2015，52（4）：792－801.Jin Yurong，Lu Kexin，Li Peng，et al. Research onSoil Water Movement Based on Stable Isotopes[J].Acta Pedologica Sinica，2015，52（4）： 792-801.

[52]張培震，王琪，馬宗晉.中國大陸現今構造運動的GPS速度場與活動地塊[J].地學前緣，2002，9（2）：430-441.Zhang Peizhen，Wang Qi，Ma Zongjin，et al.GPSVelocity Field andActive Crustal BlocksofContemporary Tectonic Deformation in ContinentalChina[J]. Earth Science Frontiers，2002，9（2）：430 -441.

[53］薄萬舉.GPS展示的中國大陸主要相對變形特征及強震活動研究[J]．地球物理學進展，2013，28（2）：599-606.Bo Wanju. Mainly Relative Deformation Features onChina Continent Revealed by GPS and Researches onStrong Earthquake Activities [J]. ProgressinGeophysics，2013，28（2）：599-606.

[54]Castaldo R， Nardis R， DeNovelis V， et al.CoseismicStress and Strain Field ChangesInvestigation Through 3-D Finite Element Modelingof DInSAR and GPS Measurements and Geological/Seismological Data：The L’Aquila （Italy） 2009Earthquake Case Study[J]. Journal of GeophysicalResearch：Solid Earth，2018，123（5）：4193-4222.

[55] Kang Y L， Shi Y R， Anderson J L. TectonicMechanism and Evolution of Eastern China During

the Early Cretaceous：A View from Magmatism in the Middle to Southern Tan - Lu Fault Zone[J]. International GeologyReview，2019，13（1）：21-46. Cubuk-Sabuncu Y，Jonsdottir K，Caudron C，et al. Temporal Seismic Velocity Changes During the 2020 Rapid Inflation at Mt borbjorn-Svartsengi，Iceland， Using Seismic Ambient Noise[J]. Geophysical ResearchLetters，2021，48（11）：GL092265.