關(guān)鍵詞雙峰寺水庫(kù)；傾斜觀測(cè)；蓄水載荷；地殼形變

中圖分類號(hào)：P315.63 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：2096-7780（2025）08-0463-09

doi：10.19987/j.dzkxjz.2024-134

Impact of the Shuangfengsi reservoir on VP broadband tiltmeter observations at Chengde station

Guo Yaya1）， Zhang Keliang2，Zhang Yulin1，Ren Junfeng1，Yang Donghui1）

1） Chengde Earthquake Monitoring Center Station，Hebei Chengde O670oo， China

2） Institute of Geology， China Earthquake Administration，Beijing looo29， China

AbstractSinceits installation，theVPbroadband tiltmeter at Chengde stationhas exhibited along-term trend of northward tilt along the north-south （NS）axis and westward ilt along the east-west （EW）axis. However，during July to August 2021，both components showedamarked acceleration in tilting，withrates approximatelyfour times higherthan the long-term average. During the same period， the nearby Shuangfengsi reservoir， located about 10km northeast of the station，wasbeing filedandreached itsmaximumwaterlevelonAugust 24，2O21.Thisstudyanalyzes thecorrelation betweenthe inclinometer dataand changes in reservoir capacityand further quantifies theresulting crustal deformation using a non-uniform surface load model. Correlation analysis revealed a strong relationship between variation in the VP broadband tiltmeterdata andchanges inreservoir capacity.Model calculations indicated that thereservoir filingresulted in a maximum crustal subsidence of 0.4421mm ，with corresponding northward tilt of 0.1179mm and eastward tilt of （20 0.0882mm . These values are consistent with the observations recorded by the VP broadband tiltmeter.A similarly strong corelation was alsoobserved between thetiltmeterdataandreservoircapacitychanges in 2O22.The findingssuggest that recent annual variation trends observed by the VP broadband tiltmeter at Chengde station are primarily atributable to loading efects from the Shuangfengsi reservoir.By quantitatively comparing modeled ground tilt induced by reservoir impoundmentwith thetiltmeter observations，this studyeffectivelyrulesout tectonic originssuch asearthquake precursors.These results help prevent potential misinterpretationofthe dataand provide a valuable exampleof how nontectonic influences should be considered in the anomaly analysis used for earthquake forecasting.

KeywordsShuangfengsi reservoir; tilt observation; water storage load; crustal deformation

0 引言

地傾斜觀測(cè)通過測(cè)量地表或近地表某處垂向位移的水平梯度或水平位移的垂向梯度獲取地殼形變信息，其連續(xù)觀測(cè)能夠反映觀測(cè)點(diǎn)附近區(qū)域的地殼應(yīng)力變化1，常作為地震前兆的重要觀測(cè)手段之一。在監(jiān)測(cè)預(yù)報(bào)實(shí)踐中，傾斜觀測(cè)出現(xiàn)的傾斜速率加快，趨勢(shì)轉(zhuǎn)折，大幅跳升等破年變異常都是地震預(yù)報(bào)的重要信息[2-4]。例如，劉仲全[5]利用云南省19套傾斜儀的異常變化對(duì)2001年1月15日姚安M6.6地震做出預(yù)測(cè)預(yù)報(bào)。然而，大量觀測(cè)和研究表明，傾斜形變觀測(cè)還受到自然環(huán)境、場(chǎng)地環(huán)境等因素的干擾，當(dāng)這些非構(gòu)造干擾的空間范圍足夠大、幅度足夠強(qiáng)或者持續(xù)時(shí)間足夠長(zhǎng)時(shí)，也會(huì)引起階躍變化或顯著的加速趨勢(shì)。例如，常熟臺(tái)附近尚湖因強(qiáng)降雨水位上升導(dǎo)致臺(tái)站3套傾斜儀出現(xiàn)傾斜速度加快的現(xiàn)象；魯?shù)?.5級(jí)地震前昭通地震臺(tái)傾斜儀記錄到強(qiáng)降雨導(dǎo)致的異常變化。此外，水庫(kù)蓄放水、隧道施工卸載等人類活動(dòng)也會(huì)引起附近臺(tái)站傾斜記錄的異常變化[。例如，三峽水庫(kù)水體荷載變化引起宜昌地震臺(tái)洞體應(yīng)變和水管傾斜儀轉(zhuǎn)折變化，黃柏河上游水庫(kù)蓄放水造成宜昌臺(tái)地傾斜觀測(cè)水管、垂直擺傾斜儀NS分量出現(xiàn)趨勢(shì)性改變[9-10]。這些研究表明，傾斜儀能夠記錄降雨、水庫(kù)蓄放水等非構(gòu)造因素引起的大幅變化，通過相關(guān)性分析、物理模型模擬等方法確定非構(gòu)造因素的影響，有助于提高觀測(cè)異常識(shí)別的準(zhǔn)確性。

承德地震臺(tái)VP寬頻帶傾斜儀2021年7—8月傾斜速率突然加快，變化幅度超過 150×10^-3′′ ，約為正常年變速率的4倍。同一時(shí)期，承德臺(tái)東北方向約

10km 處的雙峰寺水庫(kù)竣工蓄水，水庫(kù)蓄水過程中的載荷變化會(huì)引起地殼垂直及水平方向的形變[11-13]。本文利用該地震臺(tái)VP傾斜儀觀測(cè)數(shù)據(jù)，通過與同時(shí)段水庫(kù)水文資料的相關(guān)性分析以及水庫(kù)載荷引起的傾斜變化，計(jì)算分析雙峰寺水庫(kù)蓄水對(duì)承德臺(tái)形變觀測(cè)的影響。

1數(shù)據(jù)分析

1.1承德臺(tái)VP寬頻帶傾斜儀觀測(cè)

承德地震臺(tái)位于承德市雙橋區(qū)，地質(zhì)構(gòu)造上屬于華北板塊北緣燕山造山帶東段。形變觀測(cè)山洞位于臺(tái)站院內(nèi)東側(cè)，洞深約 66m ，覆蓋層厚度約 45m 臺(tái)基巖性為片麻巖，洞室密封性好，年溫差小于 0.1^°C VP寬頻帶傾斜儀安裝在山洞最里側(cè)的洞室中。臺(tái)站2km 范圍內(nèi)無(wú)采礦、采石場(chǎng)和機(jī)械廣等干擾源，觀測(cè)環(huán)境穩(wěn)定。2013年12月12日正式入網(wǎng)運(yùn)行，記錄固體潮和同震曲線能力較好，數(shù)據(jù)精度高[14-16]，連續(xù)多年在國(guó)家局資料評(píng)比中獲獎(jiǎng)。

承德臺(tái)VP寬頻帶傾斜儀運(yùn)行以來(lái)，年變曲線變化趨勢(shì)比較穩(wěn)定，除校準(zhǔn)、維修等短期影響外，總體上一直處于向北（NS向斜率為正）、向西（EW向斜率為負(fù)）傾斜的狀態(tài)。2021年7月13日開始傾斜儀兩分向均發(fā)生傾斜加速，尤其8月1一9日傾斜速率達(dá)到正常4倍；然而，8月9日發(fā)生轉(zhuǎn)折變化，NS向快速轉(zhuǎn)S向傾斜，EW向則快速轉(zhuǎn)E向傾斜；此后，逐漸恢復(fù)正常。其中，NS向在8月11日結(jié)束轉(zhuǎn)折，恢復(fù)N向傾斜，19日再次出現(xiàn)S向傾斜，但變化幅度較小，28日恢復(fù)正常；EW向8月11日后傾斜速率減緩，22日轉(zhuǎn)平，9月4日恢復(fù)W向傾斜，變化速率恢復(fù)正常。因2019年開擺維修數(shù)據(jù)漂移較大，選用

2018年和2020年數(shù)據(jù)與2021年數(shù)據(jù)做年變曲線對(duì)比（圖1），可見2021年7—8月數(shù)據(jù)有明顯“異?！?。

1.2 降雨量

承德地區(qū)一般從每年5月下旬起，雨量逐漸充沛，進(jìn)入10月降雨減少，主汛期在7一8月。查閱歷史天氣資料，2021年為近10年降雨最多的年份，年降雨量 848mm ，遠(yuǎn)高于均值 504mm 。排除儀器故障和人為干擾后，首先考慮“異常”為降雨干擾。持續(xù)長(zhǎng)時(shí)間降雨干擾，雨水滲入會(huì)改變介質(zhì)密度，引起局部區(qū)域地殼變形，觀測(cè)曲線表現(xiàn)為連續(xù)緩慢的上升或下降形態(tài)，降雨結(jié)束后，逐漸恢復(fù)；而短暫強(qiáng)降雨干擾，若降雨量超過驅(qū)動(dòng)降雨量，觀測(cè)曲線多表現(xiàn)為同步噪聲增大、數(shù)據(jù)突跳，變化時(shí)間與降雨時(shí)間同步[17-18]

1.3雙峰寺水庫(kù)水位及蓄水量

雙峰寺水庫(kù)坐落于武烈河干流，水流北入南出，水庫(kù)總庫(kù)容1.37億 m³ ，正常蓄水位 389m ，死水位 382m 攔河壩高 50.1m ，壩頂高程 396.1m ，壩長(zhǎng) 533m 2013年8月開工，2020年12月竣工驗(yàn)收。2019年5月下閘蓄水，2021年8月25日，蓄水水位最高。該水庫(kù)位于承德臺(tái)東北方向，長(zhǎng) 5.1km ，最大寬度 1.4km 近端距承德臺(tái) 7.8km ，遠(yuǎn)端距承德臺(tái) 10.9km（ 圖2）。

雙峰寺水庫(kù)2020—2021年處于蓄水階段，7—10月是主要蓄水期，水庫(kù)水位上升較快，2020年水位升高 7.3m ，庫(kù)容增加 3.586×10⁷m³ 。2021年雨水較多，8月下旬水位升高 6.0m ，達(dá)到正常蓄水位，庫(kù)容增加 4.156×10⁷m³ ，期間有多次泄洪。

從庫(kù)容變化曲線（圖3）看，2020年水庫(kù)蓄水對(duì)VP寬頻帶傾斜儀并無(wú)影響，2021年7月13日VP寬頻帶傾斜儀出現(xiàn)“異常”，此時(shí)水庫(kù)蓄水量為

5.421×10⁷m³ 。2021年有兩次快速蓄水，第一次在7月1—14日，庫(kù)容從 3.818×10⁷m³ 增加到 6.217×10⁷m³ 增加 2.399×10⁷m³ ，第二次在8月9—25日，庫(kù)容從6.097×10⁷m³ 增加到 8.486×10⁷m³ ，增加 2.389×10⁷m³ 庫(kù)容進(jìn)入第一個(gè)穩(wěn)定階段后，開始泄洪（圖4）。

2雙峰寺水庫(kù)蓄水對(duì)傾斜觀測(cè)的影響

2.1 相關(guān)性分析

已有研究表明，降雨對(duì)地形變有較強(qiáng)的調(diào)制作用，降雨影響一般不是在降雨之后立即表現(xiàn)出來(lái)，有時(shí)若干天后才有反應(yīng)，是一個(gè)既調(diào)幅又調(diào)相的過程[19]。通過對(duì)承德臺(tái)VP寬頻帶傾斜儀數(shù)據(jù)與降雨數(shù)據(jù)對(duì)比，發(fā)現(xiàn)在個(gè)別強(qiáng)對(duì)流天氣時(shí)，有出現(xiàn)與氣壓同步的噪聲變大，卻沒有明顯的持續(xù)加載或卸載過程，如2021年6月28日—7月1日累計(jì)降雨 127mm 7月12日單日降雨量 72mm ，但VP寬頻帶傾斜儀兩分向均運(yùn)行穩(wěn)定，未出現(xiàn)延遲的上升或下降形態(tài)（圖5）。8月降雨較7月明顯減少，時(shí)間上與“異?！睂?duì)應(yīng)性差，受降雨干擾的可能性較小。

相比而言，VP寬頻帶傾斜儀兩分向時(shí)間序列與雙峰寺水庫(kù)單天庫(kù)容變化呈現(xiàn)出非常強(qiáng)的相關(guān)性，2021年NS向與庫(kù)容量的變化呈正相關(guān)，相關(guān)系數(shù)為0.83，EW向與庫(kù)容量的變化呈負(fù)相關(guān)，相關(guān)系數(shù)為 -0.62 不過，這兩個(gè)分向時(shí)間序列與泄洪數(shù)據(jù)的相關(guān)性較弱。

第一蓄水階段庫(kù)容相對(duì)穩(wěn)定后，傾斜儀數(shù)據(jù)有兩次變化（圖4， ①② ，其中，NS分量呈現(xiàn)出兩次顯著上升（北向傾斜加速），而EW分量在這兩次變化過程中表現(xiàn)出加速下降（西向傾斜加速），但兩次的速率不同。時(shí)段 ① ：NS向與庫(kù)容量的變化呈負(fù)相關(guān)，相關(guān)系數(shù)為-0.62，EW向與庫(kù)容量的變化呈正相關(guān)，相關(guān)系數(shù)為0.74；時(shí)段 ② ：NS向與庫(kù)容量的變化呈正相關(guān)，相關(guān)系數(shù)為0.81，EW向與庫(kù)容量的變化呈負(fù)相關(guān)，相關(guān)系數(shù)為 -0.65 第二蓄水階段，時(shí)段③： NS向表現(xiàn)出快速下降，與庫(kù)容變化呈負(fù)相關(guān)，相關(guān)系數(shù)為-0.82，而EW向則表現(xiàn)出快速上升，與庫(kù)容變化呈正相關(guān)，相關(guān)系數(shù)為 0.95 蓄水期間的相關(guān)系數(shù)明顯高于長(zhǎng)期變化的相關(guān)性。

2.2 模擬分析

2.2.1 載荷模型

Boussinesq[20解決了集中力作用在半無(wú)限空間彈性體表面時(shí)，任意一點(diǎn)附加應(yīng)力和位移解，在此基礎(chǔ)上Love[21]和Becker、Bevis[22]給出了半空間內(nèi)一個(gè)垂直均勻壓力作用于矩形區(qū)域上時(shí)，半空間任意點(diǎn)附加應(yīng)力和位移解。DUrso和Marmo[23用MATLB代碼加實(shí)例驗(yàn)證了Love和Becker、Bevis的理論在彈性半空間表面施加垂直均勻壓力作用任意多邊形區(qū)域的應(yīng)用。

按照Love和Becker、Bevis彈性半空間載荷理論，建立一個(gè)正交直角坐標(biāo)系 （o，x，y，z） ！ σ_o 為作用力原點(diǎn)， 為直角坐標(biāo)系三分量， z 軸向下，P_z 表示均勻的垂直壓力， p（x，y，z） 表示半空間的任意點(diǎn)， 為被施加載荷的平面上的任意點(diǎn)。 p 點(diǎn)的豎向應(yīng)力及豎向位移為：

公式（1）中 d 是位移向量；公式（2）的中 表示載荷作用面積，向量 rp 的模 ? 表示受力半徑，而U 和 V 是Boussinesq勢(shì)函數(shù)；公式（3）表示兩個(gè)Lame常數(shù) λ 和 μ 與楊氏模量 E 和泊松比 ν 的關(guān)系。

聯(lián)合公式（1）、（2）和（3）可以推導(dǎo)出水平分量d_x，d_y 和垂直分量 d_z 的計(jì)算公式（4）（5）：

本文采用D'Urso和Marmo[23]的程序計(jì)算水庫(kù)蓄水載荷引起的位移。

2.2.2 模擬分析

雙峰寺水庫(kù)整個(gè)庫(kù)區(qū)的水深并不一致，靠近南面壩體附近的水深，越往北水越淺。為了提高計(jì)算精度，把雙峰寺水庫(kù)按水位從深到淺劃分3個(gè)區(qū)域，3個(gè)區(qū)域的面積分別占總面積的 30% 、 45% / 25% ，參照庫(kù)容量與平均水深設(shè)置3個(gè)區(qū)域的水位，并將3個(gè)加載區(qū)域細(xì)分 0.5^′×0.5 的矩形網(wǎng)格，計(jì)算網(wǎng)格中每個(gè)單元的載荷產(chǎn)生的位移，由于水庫(kù)區(qū)域巖漿巖和變質(zhì)巖分布較廣，計(jì)算過程中取泊松比 ν=0.25 ，楊氏模量 E=4×10¹⁰ ，根據(jù)計(jì)算結(jié)果在建模區(qū)域內(nèi)繪制半空間自由表面水平和垂直位移的等高線圖。

用非均勻載荷模型模擬水庫(kù)蓄水對(duì)地殼的影響，由圖6可知，雙峰寺水庫(kù)蓄水引起庫(kù)區(qū)附近地殼西側(cè)向東，東側(cè)向西，南側(cè)向北，北側(cè)向南的傾斜。2021年7—8月，庫(kù)容從 3.74×10⁷m³ 上升到 8.48×10⁷m³ 第一階段蓄水結(jié)束（7月13日）時(shí)，庫(kù)容 5.421×10⁷m³ 水位為 385.31m ，庫(kù)區(qū) N 向形變量 0.0833mm ，E向形變量 0.0623mm ，沉降量 0.3121mm 第二階段蓄水，達(dá)到正常蓄水位 389.76m 時(shí)，庫(kù)容 8.486×10⁷m³ 庫(kù)區(qū)N向形變量 0.1179mm ，E向形變量 0.0882mm 沉降量 0.4421mm 水庫(kù)在蓄水過程中，隨著水位上升，庫(kù)區(qū)形變量增大，影響范圍也在擴(kuò)大，離庫(kù)區(qū)距離越遠(yuǎn)形變量越小。

2.2.3 實(shí)際位移變化量

VP寬頻帶傾斜儀為擺式儀器，記錄地面傾斜變化時(shí)鉛錘擺的變化角度，標(biāo)準(zhǔn)擺錘長(zhǎng)度 L=0.1m

△φ為地傾斜角， Δh 為垂直高差，角度轉(zhuǎn)換公式：

式中， Δφ 單位：\"， Δh 單位： μm

第一階段蓄水結(jié)束，VP寬頻帶傾斜儀NS向上升 226.560×10^-3′′ ，EW向下降 111.554×10^-3′′ ，根據(jù)角度轉(zhuǎn)換公式，NS向垂向高差 0.1098mm ，EW向垂向高差 0.0541mm 第二階段蓄水結(jié)束，NS向下降110.156×10^-3′′ ，EW向上升 184.739×10^-3′′ ，根據(jù)角度轉(zhuǎn)換公式，NS向垂向高差 0.0530mm ，EW向垂向高差 0.0896mm NS向、EW向最大變幅與水庫(kù)達(dá)正常蓄水位時(shí)，載荷模型計(jì)算的 N 向、E向形變量相差不大（表1）。

3分析與討論

雙峰寺水庫(kù)庫(kù)容自2022年5—7月庫(kù)容減少約

2.046×10⁷m³ ，隨著雨季水量的補(bǔ)充，7月中旬至10月期間庫(kù)容逐漸增加約 1.645×10⁷m³ ，全年無(wú)泄洪（圖7）。2022年VP寬頻帶傾斜儀EW向仍為W向單向傾斜，但NS向年變曲線上半年N向傾斜下半年轉(zhuǎn)S向傾斜，不再是N向單向傾斜的狀態(tài)。

對(duì)2022年VP寬頻帶傾斜儀與水庫(kù)庫(kù)容數(shù)據(jù)做相關(guān)性分析，NS向與庫(kù)容量的變化呈負(fù)相關(guān)，相關(guān)系數(shù)為-0.79，EW向與庫(kù)容量變化的相關(guān)性不大。5—10月庫(kù)容明顯變化期間，NS向與庫(kù)容量變化的相關(guān)系數(shù)-0.92，相關(guān)系數(shù)明顯高于年變化的相關(guān)性。2022年VP寬頻帶傾斜儀NS向數(shù)據(jù)仍然受水庫(kù)蓄水影響。

水庫(kù)蓄水與傾斜儀數(shù)據(jù)變化表現(xiàn)出一定的相關(guān)性，但在一些時(shí)段又并非完全相關(guān)。表明水庫(kù)蓄水對(duì)傾斜儀變化有較強(qiáng)的調(diào)節(jié)作用，但不是唯一因素。作為靈敏度較高的形變手段之一，VP寬頻帶傾斜儀能夠記錄到不同距離、多空間尺度載荷的變化，本文僅通過水庫(kù)水位和降雨數(shù)據(jù)進(jìn)行分析可能并不全面。例如，本文所采用的降雨量數(shù)據(jù)來(lái)自單個(gè)監(jiān)測(cè)站，僅代表該站附近或周圍區(qū)域的變化。實(shí)際上，臺(tái)站所在區(qū)域東西兩側(cè)均有河流通過，這些河流的水位會(huì)隨著流域內(nèi)降雨量的變化而變化。此外，除了雙峰寺水庫(kù)之外，臺(tái)站東南方向有人工湖。這些水體距離臺(tái)站更近，其等量水位變化也會(huì)造成一定程度的地殼傾斜變化。需要根據(jù)臺(tái)站的地理位置、氣候條件、地質(zhì)結(jié)構(gòu)、周圍環(huán)境變化等客觀條件，收集更多相關(guān)監(jiān)測(cè)資料，建立更完善的物理模擬模型才能確保異常核實(shí)結(jié)論的準(zhǔn)確性。

4結(jié)論

通過以上相關(guān)性分析和非均勻載荷模型初步分析，認(rèn)為近年承德臺(tái)VP寬頻帶傾斜儀年變趨勢(shì)改變主要由雙峰寺水庫(kù)蓄水引起。

（1）對(duì)地殼形變的影響：彈性半空間載荷模型計(jì)算表明，雙峰寺水庫(kù)蓄水庫(kù)區(qū)地殼垂直形變量隨水庫(kù)蓄水量的增加而增大，最大達(dá) 0.4421mm ，從水庫(kù)中心向外形變逐漸減小，水平形變影響較小，E向最大形變量 0.0882mm ，N向最大形變量 0.1179mm 。

（2）對(duì)VP寬頻帶傾斜儀的影響：雙峰寺水庫(kù)蓄水對(duì)VP寬頻帶傾斜儀日變化無(wú)明顯影響，當(dāng)庫(kù)容超過 5.421×10⁷m³ 后，VP寬頻帶傾斜儀觀測(cè)數(shù)據(jù)對(duì)庫(kù)區(qū)周圍地殼形變有響應(yīng)，NS向N向傾斜，最大變幅226.560×10^-3′′ ，EW向W向傾斜，最大變幅 111.554× 10^-3′′ ，NS向響應(yīng)程度大于EW向。NS向、EW向最大變幅與水庫(kù)達(dá)正常蓄水位時(shí)，載荷模型計(jì)算的N向、E向形變量相差不大。

綜上所述，本研究通過對(duì)雙峰寺水庫(kù)蓄水影響的分析，為排除地震前兆誤判提供了依據(jù)，同時(shí)也為在地震預(yù)測(cè)異常分析中綜合考慮非構(gòu)造因素的影響提供了實(shí)例參考。

致謝

本文所用水庫(kù)數(shù)據(jù)由雙峰寺水庫(kù)工程建設(shè)管理局提供。

承德臺(tái)與雙峰寺水庫(kù)位置圖從https：//zh-cntopographic-map.com/獲取。

參考文獻(xiàn)

[1]JohnstonMJS，MortensenCE.Tilt precursors beforeearthquakesontheSan AndreasFault，Califormia[J].Science，1974，186（4168）： 1031-1034

[2]HosoyaaK.CharacteristiciltofthgroundthatprecededtheoccueneofthestrongearthquakeofMarch7952J].alfPhysics of the Earth，1952，1（2）：75-81

[3］楊德賀，袁靜，王秀英，等.形變觀測(cè)數(shù)據(jù)的多異常形態(tài)統(tǒng)一識(shí)別[J].地球物理學(xué)報(bào)，2017，60（12）：4623-4632YangDH，YuanJ，Wang XY，etal.Identificationof multi-anomaliesof precursorydeformationdata[J].Chinese JoualofGeophysics，2017， 60（12）： 4623-4632

[4]楊玲英，毛先進(jìn)，徐金.2014年景谷6.6級(jí)地震前通海地震臺(tái)傾斜形變異常響應(yīng)[J].地震地磁觀測(cè)與研究，2018，39（5）： 124-130YangLY，MaoXJ，XuJ.Tiltdeformationresponsesandanomalous analysesatTonghai Seismic StationforJingguM6.6earthquakeon Oct.7， 2014[J]. Seismological and Geomagnetic Observation and Research， 2018， 39（5）： 124-130

[5］劉仲全.姚安6.5級(jí)地震前云南傾斜場(chǎng)變化特征[J].地震研究，2001，24（4）：301-306Liu ZQ.Thevariationcharacteristicsofthetilt field in Yunnan beforethe Yaoan earthquakeM6.5[J].Journalof SeismologicalResearch，2001，24（4）：301-306

[6]丁建國(guó)，陸德明，狄梁，等.常熟臺(tái)傾斜儀降雨干擾定量分析[J].華南地震，2011，31（3）：83-88DingJG，LuDM，DiL，etal. QuantitativelydeterminationanalysisofraiinfluenceatChangshuseismic station[J].South hinaJournal of Seismology，2011， 31（3）： 83-88

[7］李智蓉，付虹，張中旭.魯?shù)?.5級(jí)地震前昭通地震臺(tái)形變異常與降雨關(guān)系分析[J].云南大學(xué)學(xué)報(bào)（自然科學(xué)版），2018，40（5）： 919-928Li ZR，F(xiàn)u H，ZhangZX.Analysis ontherelationshipbetweendeformationanomalyandrainfalin ZhaotongstationbeforeLudian6.5 earthquake[J]. Journal of Yunnan University （Natural Sciences Edition），2018，40（5）： 919-928

[8］何應(yīng)文，岳沖，李智蓉，等.愛華隧道施工對(duì)云南云縣臺(tái)水管傾斜觀測(cè)影響分析[J].地震工程學(xué)報(bào)，2018，40（增刊1）：117-122HeY W，YueC，Li ZR，etal. Influenceof AiuatunnelconstructionontheobservationofDSQ water tubetiltmeterat Yunxianseismic station，Yunnan Province[J]. China Earthquake Engineering Journal， 2018，40（S1）：117-122

[9］孫伶俐，姚運(yùn)生，蔣玲霞，等.三峽水庫(kù)水體荷載變化對(duì)地殼形變的影響[J].長(zhǎng)江科學(xué)院院報(bào)，2015，32（12）：46-50SunL，YaoYS，JangLX，talEfectofwaterloadvarationintheThe GorgeResevoirotecrustaldefomationJ].JualofYangtze River Scientific Research Institute，2015，32（12）：46-50

[10］袁曲，黃仲，王慧，等.黃柏河水位變化對(duì)宜昌臺(tái)地傾斜觀測(cè)的影響[J].大地測(cè)量與地球動(dòng)力學(xué)，2009，29（增刊1）：79-83YuanQ，HuangZ，WangH，etal.EfetoflevelchangeofHangbairiverongroundtiltobservationsatichangstation[J]Joualof Geodesy and Geodynamics，2009，29（S1）： 79-83

[11］尹財(cái)，章傳銀，王偉，等.三峽水庫(kù)蓄放水對(duì)地殼形變影響的計(jì)算分析[J].測(cè)繪地理信息，2017，42（3）：102-105YinC，Zhang CY，Wang W，etal.Calculation and analysisoftheinfluenceof waterstorage and drainage inthe ThreeGorgesReservoir on crustal deformation[J]. Journal of Geomatics，2017，42（3）：102-105

[12］楊志榮，李士柱，黃宗英，等.新疆克孜爾水庫(kù)形變前兆異常的初步研究[J].內(nèi)陸地震，1998，12（2）：152-160Yang ZR，LiSZ，Huang ZY，etalrimarystudyondeformationanomalyatKizilreservoir[J].InlandEarthquake，998，12（2）：152-160

[13］杜瑞林，喬學(xué)軍，王琪，等.長(zhǎng)江三峽水庫(kù)蓄水荷載地殼形變：GPS觀測(cè)研究[J].自然科學(xué)進(jìn)展，2004，14（9）：1006-1011DuRL，QiaoXJ，WangQ，etal.Crustal deformatation byGPSin Three Gorges Resservoirafterimpoundment[J].ProgressinNatural Science，2004，14（9）：1006-1011

[14］張玉林，宋曉冰，郭亞亞，等.承德地震臺(tái)VP型寬頻帶垂直擺運(yùn)行效能分析[J].地震地磁觀測(cè)與研究，2018，39（2）：108-115Zhang YL ， Song X B，Guo Y Y，etal. Theanalysis of operation eficiencyof the VP type vertical pendulum at Chengde SeismicStation[J]. Seismological and Geomagnetic Observation and Research， 2018， 39（2）： 108-115

[15］王嘉琦，楊東輝，常玉柱，等.河北區(qū)域VP型寬頻帶垂直擺傾斜儀同震響應(yīng)分析：以2019年河北唐山 M4.5 地震為

例[J].山西地震，2022（1）：28-30，56 Wang JQ，YangDH， Chang Y Z，etal.Co-seismic response analysis of VPbroadband vertical pendulum tiltmeter in Hebei region： Takingthe2019TangshanM4.5earthquakeinHebeiProvinceasanexample[J].EarthquakeResearchinShanxi2O2（1）：28-3056 [16］袁國(guó)旭，王嘉琦，李萬(wàn)里.承德雙峰寺水庫(kù)臺(tái)網(wǎng)建設(shè)及監(jiān)測(cè)能力評(píng)估[J].地震科學(xué)進(jìn)展，2025，55（6）：340-347 Yuan GX，WangJQ，LiWL.Networkconstructionand monitoringcapabiltyevaluationoftheShuangfengsireservoirinChengde[J]. Progress in Earthquake Sciences，，2025，55（6）：340-347 [17］楊學(xué)慧，常玉巧，李燕玲，等.用同相疊加方法識(shí)別定點(diǎn)形變數(shù)據(jù)中的降雨影響[J].大地測(cè)量與地球動(dòng)力學(xué)，2016， 36（10）： 936-940 Yang XH，Chang YQ，LiYL，etal.Distinguishg theainfalldisturbing pattincrustaldeformationrecordingbyi-asstack of data[J]. Journal of Geodesy and Geodynamics，2016，36（10）： 936-940

[18]呂琳，楊艷芳，舒德亮.武漢臺(tái)形變觀測(cè)與降雨參數(shù)之間的定量分析與研究[J].地震工程學(xué)報(bào)，2020，42（1）：136-142 LuL，YangYF，ShuDL.Quantitativeanalysisandmechanismoftherelationbetweendeformationobservationandrainfal parameters based on data from Wuhan station[J]. China Earthquake Engineering Journal，2020，42（1）：136-142 [19］劉序儼.應(yīng)用褶積同態(tài)濾波排除降雨對(duì)地形變觀測(cè)的干擾[J].地震，1985（6）：48-51 LiuXY.Eliminatingprecipitatioefectsonbsevationdataofcustaldeformationbymeansofhomomophicconvolutiofilter[J]. Earthquake， 1985（6）： 48-51 20]BouinesqJApplicationdespotentielsaetuedelequlibreetdumouventdessoides élastiques：ricipalemntaualul des deformations etdespresionsque produisent，danscessolides，des efortsquelconquesexeres surunepetepartiedeleursurfaceou de leur interieur：memoiresuivi de notes étendues sur divers points de physique，mathematiqueet danalyse[M].Paris：GauthierVillars，1885： 231 21]LoveAEIXetrsprodediitesodyprueartofthoudyJ]losoaacsfte Royal Society A，1929，228（659-669）：377-420 [22]Becker JM， Bevis M. Love's problem[J]. Geophysical Jourmal International， 2004，156（2）： 171-178 [23]D'Urso MG，Marmo F.On a generalized Love’s problem[J]. Computers amp; Geosciences，2013，61： 144-151