山東日照地區不同類型含水層地下水位多年動態變化特征及其影響因素

楊鵬，張壽川，郭本力，等.山東日照地區不同類型含水層地下水位多年動態變化特征及其影響因素.吉林大學學報（地球科學版）， 2024，54（3）：967979. doi：10.13278/j.cnki.jjuese.20220322.

Yang Peng，Zhang Shouchuan，Guo Benli，et al. Characteristic of Groundwater Level Dynamic Variation and the Potential Influence Factors in Different Aquifer Categories in Rizhao Area， Shandong Province. Journal of Jilin University （Earth Science Edition） ， 2024，54（3）：967979. doi：10.13278/j.cnki.jjuese.20220322.

摘要：自然條件下沿海地區地下水與海水呈現動態平衡，但地下水超采會破壞二者間平衡關系，進而引發一系列地質環境問題。地下水動態變化受多種因素影響呈現非平穩信號變化特征，傳統時間序列分析方法無法從多角度準確描述時間序列動態變化特征，基于時間域和頻率域分解的小波分析方法為探究地下水動態變化特征及影響因素提供了便利條件。利用山東日照地區不同類型含水層3口監測井1996—2015年間地下水位和降雨量監測數據，應用連續小波變換法和小波相干分析法開展研究，結果表明：不同類型含水層地下水位動態變化周期特征存在差異，

荻水監測井（第四系松散巖類孔隙水含水層）地下水位呈現年際、4～5 a和11 a周期變化，莒縣地震辦監測井（基巖裂隙含水層）和

東莞中學監測井（碳酸鹽巖類裂隙巖溶水含水層）地下水位均呈現11 a周期和14 a強趨勢項變化，總體而言，該地區降雨量呈現年際周期、4～5 a周期、11 a周期變化以及14 a強趨勢項變化；降雨量周期變化是導致不同類型含水層地下水位呈現年際周期變化和多年周期變化的主要因素；不同含水層地下水位對降雨量響應時間存在差異是造成其呈現不同周期變化特征的主要原因。

關鍵詞：地下水動態; 連續小波變換; 小波相干; 時滯性；日照地區

doi：10.13278/j.cnki.jjuese.20220322

中圖分類號：P641.8

文獻標志碼：A

收稿日期：20221128

作者簡介：楊鵬（1986-），男，高級工程師，主要從事水文地質、環境地質方面的研究，E-mail： 9804179@qq.com

通信作者：郭本力（1983-），男，高級工程師，主要從事水文地質、環境地質方面的研究，E-mail： 286486736@qq.com

基金項目：日照市城市地質調查項目（SDGP371100202102000475）；中國地質調查局項目（DD202216772）；中央財政基本科研業務費資助項目（JKY202004）

Supported by the Project of Rizhao Geological Survey （SDGP371100202102000475）， the Project of China Geological Survey （DD202216772） and the Fundamental Research Funds of CGS Research （JKY202004）

Characteristic of Groundwater Level Dynamic Variation and the Potential Influence Factors in Different Aquifer Categories in Rizhao Area， Shandong Province

Yang Peng1，2，3，4， Zhang Shouchuan5， Guo Benli1，2，3，4， Yang Chongjing1，2，3，4，

Feng Qiyuan1，2，3，4， Lü Youcheng1，2，3，4， Xie Luxi1，2，3，4

1. The Eighth Institute of Geology and Mineral Resources Exploration of Shandong Province， Rizhao 276826， Shandong， China

2. Key Laboratory of Nonferrous Metal Ore Exploration and Resource Evaluation of Shandong Provincial Bureau of Geology ""and Mineral Resources， Rizhao 276826， Shandong， China

3. Rizhao Big Data Research Institute of Geology and Geographic Information， Rizhao 276826， Shandong， China

4. Rizhao Key Laboratory of Land Quality Evaluation and Pollution Remediation， Rizhao 276826， Shandong， China

5. Chinese Academy of Geological Sciences， Beijing 100037， China

Abstract：

Under the natural conditions， there is the dynamic balance between the groundwater and sea in the coastal region. However， the overexploitation of groundwater causes the imbalance between them， which also induces the various of environmental problems. The analysis methods of traditional time series cannot accurately reveal the dynamic variation characteristic in both time and frequency domains， while wavelet analysis methods provide convenient way to explore the dynamic variation characteristics and external influence factors of groundwater. The continuous wavelet transform and wavelet coherence analysis are used to analyze the time series data of groundwater level in three monitoring wells and precipitation in Rizhao area during the interval of 19962015. The results show that： The periodic characteristics of groundwater level variation in different aquifer categories are different. The groundwater level in Dishui monitoring well （pore water of shallow loose rocks） shows periodic variation characteristics which are one year， 45 a and 11 a. The groundwater level in Earthquake Bureau of Juxian well （bedrock fissure water） and Dongguan High School well （carbonate fissure karst water） show the periodic variation of 11 a" and the strong tendency of 14 a. The precipitation in Rizhao area shows several different periodic variations， including one year， 45 a， 11 a and the strong tendency of 14 a." The periodic variation of precipitation in Rizhao area is the important influence factor on the variation of groundwater level which show the periodic variation of one year and several years." The lag time between the groundwater change and the variation of precipitation in different aquifer type is different， which is the main reason for the diversity characteristic of groundwater level periodic variation in different aquifer type.

Key words：

groundwater regime; continuous wavelet transform; wavelet coherence; lag time; Rizhao area

0" 引言

地下水是水資源的重要組成部分，其不僅是重要的地質信息載體［15］，也在飲水保障、農田灌溉、工業生產和生態環境保護等方面發揮著重要作用［68］。但由于人類活動加劇和極端氣象事件頻發，地下水資源管理面臨巨大挑戰。近年來，國家地下水監測工程的實施，不僅有助于深入了解地下水動態變化規律，也為地下水資源開發、利用提供了科學依據和有力保障［911］。

沿海地區地下水在自然條件下與海水達到動態平衡，但過量開采地下水則會造成地下水位下降，極易引發海水入侵、水質惡化、生態環境退化等問題［1213］。山東省日照市位于黃海之濱，該地區地下水監測工作已開展20余年，前人［1415］在不同類型含水層中均布設了監測井，已積累海量監測數據；但對于該地區地下水位動態變化規律的研究十分匱乏，僅有的研究只局限在基于時間域角度的地下水位變化趨勢分析方面，導致對該地區地下水動態變化規律的認識十分有限，缺乏對地下水監測數據所蘊含地質信息的深入挖掘；同時，該地區工農業生產對地下水的不合理開采，破壞了地下水和海水間的動態平衡，引發了海水入侵、地下咸水面積增加等一系列生態環境問題，嚴重威脅該地區地下水安全。因此，開展日照地區不同類型含水層地下水動態變化特征及影響因素研究迫在眉睫。

前人對時間序列的研究主要基于滑動平均法、濾波分析法和傅里葉變換法等方法［1618］，但這些方法具有一定的局限性，不僅無法準確獲取非平穩信號在局部時間域內的頻域特征，也無法對時間序列中的突變點進行嚴謹的數學判斷。地下水受構造因素、氣象因素以及人類生產生活的影響［1923］，呈現非平穩信號特征，因此無法利用上述分析方法開展研究。近些年，從傅里葉變換法演變而來的小波分析法有效彌補了上述幾種分析方法的不足，為識別地下水動態變化特征開辟了新途徑。該方法為結合樣條分析、調和分析、傅里葉分析、數值分析和泛函數分析的有效時頻分析工具［2425］，主要優勢為：其可以從時間域和頻率域兩個角度分析數據動態變化特征，并通過在時間域和頻率域內搜索不同時間序列共同特征進而分析二者間的相關性和滯后性。已有基于小波分析法對地下水動態的研究大多集中于識別單一含水層地下水動態變化特征，而對同一地區不同類型含水層地下水動態變化特征的對比分析關注較少，地下水位變化影響因素及機制尚未得到深入揭示。

受分析方法的制約，已有研究對日照地區地下水動態變化的分析僅停留在時間域角度層面，尚未開展基于時間域和頻率域的多角度地下水動態變化特征及影響因素的研究。為此，本文選取山東省日照地區不同類型含水層監測井——荻水監測井（第四系松散巖類孔隙水）、莒縣地震辦監測井（基巖裂隙水）和東莞中學監測井（碳酸鹽巖類裂隙巖溶水）為研究對象，利用多年地下水位和降雨量監測數據，基于連續小波變換法分析時間序列數據周期特征，采用小波相干分析法研究不同時間序列數據間的相關性和時滯性；以期準確識別自然條件下日照地區不同類型含水層地下水位動態變化特征及影響因素，為合理開發利用地下水資源，保護生態環境提供科學依據。

1 "區域概況

日照市位于山東省東南部黃海之濱，陸域面積5 317 km2（圖1），海域面積5 630 km2，大陸海岸線

總長226 km，地形地貌以山地、丘陵為主，二者約占總面積的69%，屬于典型魯東山地丘陵地區，域內最高點為五蓮縣馬耳山，海拔706 m。

1.1" 水文地質概況

日照地處一級構造單元魯東斷塊內的二級構造單元膠南隆起，沂沭斷裂帶末端和五蓮—成斷裂帶交會處，區域內構造復雜。結合地下水賦存、埋藏條件以及含水層性質，區域內地下水可劃分為4種不同類型［1415］：第四系松散巖類孔隙水、碎屑巖類孔隙裂隙水、基巖裂隙水和碳酸鹽巖類裂隙巖溶水（圖1）。其中：第四系松散巖類孔隙水主要分布于河流兩側沖積、沖洪積砂石礫層中，水化學類型為HCO3·ClCa·Na，該類型地下水分布范圍較廣，供水能力較強；碎屑巖類孔隙裂隙水分布于莒縣東部、五蓮縣西部和北部區域，含水層為白堊系礫巖、粉砂巖和碎屑巖，水化學類型為HCO3Ca·Na；基巖裂隙水分布于莒縣東部、五蓮縣東部和南部，以及東港區南部區域，含水層巖性主要為燕山晚期花崗巖，水化學類型為HCO3Ca·Mg，受控于構造條件等因素影響，基巖裂隙水分布范圍最廣，但水量較小，供水能力不足；碳酸鹽巖類裂隙巖溶水主要分布于莒縣北部區域，分布較少。

區域內地下水主要補給來源為大氣降水［15］，由于地勢呈現北高南低、西高東低的特點，因此地下水總體呈現自西北向東南流動的趨勢。該區域第四系沉積物組成物質較細，主要包括亞黏土、含礫亞黏土和粉砂土等，加之地形起伏較小、水力坡度小等特點，故大氣降水補給過程及地下水徑流速度較為緩慢。第四系沉積物之下為古風化裂隙水，由于與上層第四系松散巖類孔隙水聯系較差，因此該層地下水具有微承壓性。在基巖露頭區，大氣降水沿風化裂隙、斷層和其他導水裂隙入滲、補給，形成風化裂隙和構造裂隙水。

荻水監測井，位于日照市嵐山區荻水村，井深6.1 m，含水層類型為第四系松散巖類孔隙水，鉆孔柱狀圖如圖2a所示：0～1.1 m為耕植土；1.1～2.2 m為粉質黏土，粉土與黏土體積比約為3∶7；2.2～6.1 m為強風化花崗巖。莒縣地震辦監測井，位于日照市莒縣城陽街道，井深209.0 m，含水層類型為基巖裂隙水，鉆孔柱狀圖如圖2b所示：0～0.5 m為耕植土；0.5～9.0 m為強風化花崗巖，含水量較少；9.0～16.0 m為中風化花崗巖；16.0～40.0 m為石英巖；40.0～45.0 m為花崗巖，較為破碎，含水量較大；45.0～80.0 m為石英巖；80.0～115.0 m為花崗巖；115.0～118.0 m為石英巖，含水量較大；118.0～209.0 m為花崗巖，含水量較少。東莞中學監測井，位于日照市莒縣東莞鎮，井深5.9 m，含水層為碳酸鹽巖類裂隙巖溶水，鉆孔柱狀圖如圖2c所示，0～0.8 m為耕植土，0.8～5.9 m為強風化砂頁巖。

1.2" 區域地下水動態

日照地區地下水動態監測工作始于1993年，積

累了豐富的監測數據。監測站點的布設考慮不同地貌單元和地下水補徑排條件，主要分布于河流兩岸、供水水源地、海水入侵和地下水污染嚴重地段［15］。

地下水動態變化主要受氣象因素和人類生產生活用水影響［2628］。降雨入滲是日照地區地下水資源的主要補給來源［15， 22］。該地區屬暖溫帶濕潤季風大陸性氣候，雨量充沛，多年平均降雨量為817.6 mm，降雨多集中在每年6—8月，該時間段內降雨量約占全年降雨量的72.7%（圖3）。計算3口監測井水位以及降雨量多年月均值發現，每年6—9月井水位呈現明顯上升趨勢，每年10月—次年3月井水位呈現下降趨勢。

2" 研究方法

時間序列根據特征不同可分為兩種類型，分別為具有平穩信號特征的時間序列和具有非平穩信號特征的時間序列。其中：平穩信號是指參數分布規律不隨時間變化的信號；非平穩信號是指需要分析時間域和頻率域特性的時變信號或頻變信號，其參數特征（如均值、方差等）隨時間變化。

平穩信號特征時間序列常采用傅里葉變換方法進行分析;但該方法無法全面刻畫突變信號和非平穩信號的局部特征，并不適用于分析具有非平穩信號特征的時間序列。由于非平穩信號的頻率特性會隨時間變化，為了捕獲這一時變特性，需要對信號進行時頻分析，常用的時頻分析方法包括小波變換、希爾伯特變換、短時傅里葉變換等。地下水動態由于受到多種外部因素疊加影響而呈現非平穩動態變化的特征，在分析時需要應用時頻分析工具對其同時進行時間域和頻率域信息提取。

本研究中主要采用小波變換開展研究工作，該方法發展至今已40余年，在自然科學領域得到廣泛應用［2931］。該方法可以同時將時間序列數據在時間域和頻率域內進行分解并分析其固有特征（圖4），主要包括連續小波變換（continuous wavelet transform，CWT）和離散小波變換（disperse wavelet transform，DWT）兩種方法。在本次研究中應用連續小波變換方法開展時間序列特征分析［3233］。

圖中彩色曲線為動態監測示例數據。

2.1" 連續小波變換

連續小波變換通過利用已知頻率和振幅的小波基函數，將時間序列數據投影到相應的時間域和頻率域空間內，分析其動態變化特征。對時間序列進行連續小波變換時，選擇恰當的小波基函數尤為重要，常用小波基函數主要包括Morlet函數、Haar函數、Mexican Hat函數、Gaus函數等。在本次研究中，采用Morlet函數作為小波基函數，該函數是高斯包絡下的單頻率復正弦函數，具有非正交分解特性，可用于連續小波變換。

對于連續函數f（x），其連續小波變換可表示為

Wf（a，b）=（1/a）∫

f（t）ψ（（t－b）/a）dt。（1）

式中：W為小波變換系數；a為伸縮尺度；b為平移因子；a，b均為實數且agt;0；ψ（x）通過ψ（x）=（23π－1/4）（1－x2）e－x2/2計算；t為時間間隔。

對于長度為N的時間序列而言，其連續小波變換可表示為

W（a，b）=（1/a）∑Nn=1x（n）φ（（n－b）/a）。 （2）

式中：φ（（n-b）/a）為子小波；x（n）為變換系數。

小波方差可以確定時間序列中的主周期，其計算方式為將b域內所有小波系數的平方積分，可表示為

Var（a）=∫Wp（a，b）2db。（3）

式中：Var（a）為小波方差；Wp（a， b）為信號函數。

2.2" 小波相干

小波相干方法通過搜索兩個時間序列在時間域和頻率域內的共同特征，從而有效識別兩個時間序列的相關性。對于兩個具有周期特征的時間序列，常采用小波相干研究二者相關性。與傳統相關系數定義相似，兩個時間序列x、y小波相干函數的平方可定義［34］為

R2（x，y）=S（s-1W（x，y））2S（s-1W（x））·S（s-1W（y）） 。（4）

式中：R2（x，y）為小波相干函數的平方；S表示平滑操作；s為小波尺度。當只有一個時間序列時，W表示連續小波變換，當有兩個不同時間序列時，W表示兩個時間序列的小波相干變換。“相干”是指小波相干的平方，其取值為0～1之間。

2.3" 水位數據收集及處理

本次研究涉及數據主要包括日照地區不同類型含水層地下水監測數據和降雨數據。為準確掌握該地區地下水時空動態變化特征，在該地區不同類型含水層均布設監測井，監測時間間隔為5 d。由于部分監測井監測數據出現缺失，故選取數據連續性較強的監測井為研究對象。圖5為荻水、莒縣地震辦和東莞中學3口監測井1996—2015年水位動態變化圖。如圖5所示：荻水監測井在1996—2009年間水位呈現動態變化趨勢，水位變幅小于1 m，2009—2015年間水位呈現下降趨勢，水位變幅為3 m；莒縣地震辦監測井在1996—2010年間水位呈現平穩波動變化趨勢，在2010—2015年間水位呈現下降趨勢；東莞中學監測井在1996—2015年間水位未呈現明顯上升或下降趨勢，水位變幅小于2.5 m。降雨數據來自于中國氣象數據網（http：//data.cma.cn/）。

利用連續小波變換和小波相干分析法進行多時間尺度分析的步驟為：1）數據預處理。小波分析需要等時間間隔的數據，需將水位監測數據取月均值、降雨數據取月累積降雨量。2）連續小波變換。將預處理數據代入式（2）中，并選取Morlet小波基函數進行計算，得到小波系數。小波系數為復數，通過計算可以得到小波系數實部，繪制小波系數實部等值線圖。將預處理數據代入式（3）可以得到小波方差。通過上述計算可以得到水位和降雨量周期變化特征。3）小波相干分析。將預處理的水位和降雨量數據代入式（4）可以得到兩個時間序列在不同時間尺度下的相關性，分析導致水位呈現周期變化的外部影響因素。

3" 結果

對地下水動態監測數據進行連續小波變換可以得到小波系數方差圖（圖6）、小波系數等值線圖（圖7）和不同周期小波系數動態變化圖（圖8）。小波系數方差圖通過突出不同周期振幅強弱的差異性幫助分析時間序列周期變化特征，小波系數等值線圖能夠反映時間序列中所包含的各周期與振幅大小，不

同周期小波系數動態變化圖可以反映同一周期內振幅隨時間變化規律。

荻水監測井水位2005—2008年期間達到最大值，2010—2012年期間達到最小值（圖5）。由圖6可知該井水位存在年際周期變化，從圖7a可知井水位在1996—2015年間存在11 a周期變化，在2010—2015年間存在4～5 a周期變化。由圖8a結合圖5可知：荻水監測井11 a周期變化小波系數振幅最大且隨時間變化較為穩定，4～5 a周期小波系數振幅呈現由弱變強的變化趨勢，年際周期小波系數振幅最小；在2005—2008年期間，在4～5 a周期和11 a周期小波系數峰值共同影響下，該井水位達到最高值；在2010—2012年期間，由于4～5 a周期和11 a周期小波系數處于谷值，該井水位降至最低值。

莒縣地震辦監測井水位在2003—2004年期間達到最大值，在2010年之后呈現逐漸下降趨勢（圖5）。

圖6表明該井水位存在14 a強趨勢變化和11 a周期變化。由圖7b

可知該井水位在1996—2015年間存在14 a強趨勢變化，在2000—2015年間存在

11 a周期變化。由圖8b結合圖5可以看出：莒縣地震辦監測井

14 a強趨勢項振幅最大且隨時間變化較為穩定，11 a周期小波系數振幅較小且隨時間呈現強—弱—強的不穩定變化趨勢；在2003—2004年期間，受到14 a強趨勢項小波系數峰值影響，該井水位達到最高值，在2010年以后，由于14 a強趨勢項和11 a周期小波系數均處于向谷值變化趨勢，該井水位在2010年之后持續下降。

東莞中學監測井水位在1996—1997年期間達到最低值，在2001—2005年期間達到最高值（圖5）。由圖6結合圖7可知，該井水位在2000—2013年間存在14 a強趨勢變化，在1996—2015年間存在11 a較弱周期變化。由圖8c結合圖5可以看出：東莞中學監測井14 a強趨勢項小波系數振幅最大且隨時間變化較為穩定，11 a周期小波系數振幅較小且

隨時間變化呈現強—弱—強變化趨勢；在1996—1997年期間，受到14 a強趨勢項和11 a周期小波系數谷值影響，該井水位降至最低值；在2001—2005年期間，受14 a強趨勢項小波系數峰值影響，井水位達最高值，在這期間雖然受11 a周期小波系數谷值影響，但因其振幅較小，其對水位變化影響較小，因此井水位逐漸上升至歷史最高值。

對日照地區降雨量而言，結合小波方差圖（圖6）和小波系數等值線圖（圖7）可知，在1996—2015年間日照地區降雨量呈現年際周期、4～5 a周期和11 a周期變化，以及14 a強趨勢項變化特征。

4" 討論

4.1" 地下水位年際周期變化影響因素

通過計算不同監測井水位多年月均值可以發現，3口監測井水位均呈現夏季水位上升、冬季水位下降的年際動態變化趨勢。為進一步分析不同類型含水層地下水動態與降雨之間的關系，分別計算荻水監測井、莒縣地震辦監測井和東莞中學監測井水位與降雨量之間的小波相干函數。

圖9為 3口監測井水位和降雨量小波相干性分析結果，可以得出：對荻水監測井而言，1998—2015年間，在1 a時間尺度上，該井水位與降雨量呈現高度相關性（圖9a）；對

莒縣地震辦監測井而言，1997—2014年間在1 a時間尺度上，該井水位與降雨量呈現高度相關性（圖9b）；對東莞中學監測井而言，1998—2015年間在1 a時間尺度上以及2003—2012年間在4～5 a時間尺度上，該井水位與降雨量呈現高度相關性（圖9c）。因此，3口監測井水位年際動態變化特征均由降雨所致。

進一步分析小波相干結果圖（圖9），可以得出不同含水層中地下水位動態變化與降雨量時滯性關系。

圖9中半數以上箭頭指向右側，表明3口地下水監測井水位變化滯后于降雨量變化。再根據相位角確定滯后時間可知，荻水監測井、莒縣地震辦監測井、東莞中學監測井水位與降雨量之間平均相位角分別為45°、60°和30°；表明不同監測井水位變化對降雨響應時間存在差異，荻水監測井（第四系松散巖類孔隙水）水位變化滯后于降雨約1.5月，莒縣地震辦監測井（基巖裂隙水）水位變化滯后于降雨約2.0月，東莞中學監測井（碳酸鹽巖類裂隙巖溶水）水位變化滯后于降雨約1.0月。

4.2" 地下水位多年周期變化動態變化影響因素

小波系數等值線圖和小波方差結果表明，3口地下水監測井水位除呈現年際周期變化（荻水監測井）和年際動態變化（莒縣地震辦監測井和東莞中學

監測井）特征外，還呈現多年周期動態變化特征。通過計算得到3口地下水監測井水位與降雨量不同周期小波系數相關系數（圖10），用以分析井水位多年動態變化影響因素。

荻水監測井水位（圖10a，d）1996—2015年間呈現11 a周期變化，在2010—2015年間呈現4～5 a周期變化特征。計算水位與降雨量小波系數相關性可知，對于2010—2015年間的4～5 a周期變化而言，兩者相關系數為0.84，

具有較強相關性，

對于19962015年間的11 a周期變化而言，二者相關系數為0.82，

也具有較強相關性；說明該井水位4～5 a

周期和11 a周期變化主要受降雨量4～5 a

周期和11 a周期變化影響。

莒縣地震辦監測井水位（圖10b，e）在1996—2015年間呈現11 a周期變化和14 a強趨勢項變化。分別計算井水位11 a周期變化、14 a強趨勢項小波系數與對應的降雨量小波系數的相關性可知

，對11 a周期變化而言，兩者相關系數為0.51，有一定相關性

，

對于14 a強趨勢項而言，兩者相關系數為0.84，具有較強相關性；說明該井水位11 a周期變化

和14 a強趨勢項

均為降雨影響所致。

東莞中學監測井水位（圖10c，f）在1996—2015年間呈現11 a周期變化和14 a強趨勢項變化。分別計算井水位11 a周期變化、14 a強趨勢項小波系數與對應的降雨量小波系數的相關性可知，對于11 a周期變化而言，兩者相關系數為0.59，有一定相關性

，對于14 a強趨勢項而言，二者相關系數為0.83，具有較強相關性；說明該井水位11 a周期變化

和14 a強趨勢項均為降雨影響所致。

4.3" 不同類型含水層水位多年動態變化特征

前人［3536］研究表明，日照地區降雨周期特征主要受太平洋年代際濤動和南方濤動影響呈現年際周期變化和4～5 a周期變化，同時還受太陽黑子極值出現周期變化影響呈現11 a周期變化［37］。

荻水監測井連通第四系松散巖類孔隙含水層，為非承壓含水層，含水層埋藏深度較淺，頂板埋深為3.3 m，其上覆弱透水層為粉質黏土，厚度僅為1.1 m，降雨補給地下水速度較快，井水位受降雨的年際周期、4～5 a周期和11 a周期變化較為明顯；因此，第四系松散巖類孔隙含水層水位在不同時間段內呈現年際周期、4～5 a周期和11 a周期變化。

莒縣地震辦監測井連通基巖裂隙含水層，為承壓含水層，該井連通的兩個含水層埋深分別為40.0～45.0 m和115.0～118.0 m，其中40.0～45.0 m含水層巖性為花崗巖，115.0～118.0 m含水層巖性為石英巖，裂隙發育，上覆弱透水層厚度較大，降雨補給地下水速度較慢，含水層滲透系數較小；因此，基巖裂隙含水層水位未呈現年際周期變化，僅呈現11 a周期變化和14 a強趨勢項變化。

東莞中學監測井連通含水層為碳酸鹽巖類裂隙巖溶含水層，含水層巖性為砂頁巖，巖石風化嚴重，含水層埋深較淺，但由于裂隙巖溶含水層較為復雜且裂隙發育，導致含水層與深部含水層相連通，大氣降水補給深層含水層速度較慢；因此深部含水層地下水位呈現11 a周期變化和14 a強趨勢項變化，從而導致接受深部含水層越流補給的東莞中學監測井水位也呈現11 a周期變化和14 a強趨勢變化。

5" 結論

1） 日照地區不同類型含水層地下水位動態變化特征存在差異，荻水監測井（第四系松散巖類孔隙水）水位在不同時間段內呈現年際周期、4～5 a周期和11 a周期變化特征，莒縣地震辦監測井（基巖裂隙水）和東莞中學監測井（碳酸鹽巖類裂隙巖溶水）水位均呈現11 a周期變化和14 a強趨勢項變化特征。總體而言，該地區降雨量呈現年際周期、4～5 a周期、11 a周期變化以及14 a強趨勢項變化特征。

2） 小波相干分析結果表明，水位變化滯后于降雨量變化，降雨量周期變化是造成地下水位年際動態變化和多年周期變化的主要影響因素。

3）時滯性分析結果表明，荻水監測井（第四系松散巖類孔隙水）、莒縣地震辦監測井（基巖裂隙水）和東莞中學監測井（碳酸鹽巖類裂隙巖溶水）水位變化滯后于降雨的時間分別為約1.5月、約2.0月和約1.0月。不同含水層地下水對降雨量響應時間存在差異是造成其呈現不同周期變化特征的主要原因。

參考文獻（References）：

［1］" Zhang S C， Shi Z M， Wang G C， et al. Quantitative Assessment of the Mechanisms of Earthquake Induced Groundwater Level Change in the MP Well， Three Gorges Area［J］. Pure and Applied Geophysics， 2018， 175（1）： 24752484.

［2］" Wang C Y， Manga M， Shirzaei M， et al. Induced Seismicity in Oklahoma Affects Shallow Groundwater［J］. Seismological Research Letters， 2017， 88（1）： 956962.

［3］" 常廷改， 胡曉. 水庫誘發地震研究進展［J］. 水利學報， 2018， 49（1）： 11091122.

Chang Tinggai， Hu Xiao. Research Progress on Reservoir Induced Earthquake［J］. Journal of Hydraulic Engineering， 2018， 49（1）： 11091122.

［4］" 車用太， 魚金子， 劉成龍， 等. 水庫誘發地震地下水監測網建設的若干問題探討［J］. 地震學報， 2010， 10（1）： 203213.

Che Yongtai， Yu Jinzi， Liu Chenglong， et al. A Discussion on Some Problems in Groundwater Network Construction for Monitoring Reservoir-Induced Earthquakes［J］. Acta Seismological Sinica， 2010， 10（1）： 203213.

［5］" Qu S， Shi Z M， Wang G C， et al. Using Water-Level Fluctuations in Response to Earth-Tide and Barometric-Pressure Changes to Measure the in-Situ Hydrogeological Properties of an Overburden Aquifer in a Coalfield［J］. Hydrogeology Journal， 2020， 50（1）： 115.

［6］" 趙全升， 孔智涵， 胡舒婭，等.柴達木盆地馬海鹽湖地下鹵水地球物理探測及應用［J］.吉林大學學報（地球科學版）， 2023， 53（5）：15601572.

Zhao Quansheng， Kong Zhihan， Hu Shuya， et al. Geophysical Exploration and Application of Underground Brine of Mahai Salt Lake in Qaidam Basin［J］. Journal of Jilin University （Earth Science Edition）， 2023， 53（5）： 15601572.

［7］" 何冠儒， 史浙明. 地下水對氣壓和固體潮響應研究進展［J］. 地震研究， 2021， 44（1）： 541549.

He Guanru， Shi Zheming. Advance in the Groundwater Level Response to Barometric Pressure and Earth Tide［J］. Journal of Seismological Research， 2021， 44（1）： 541549.

［8］" 蘭雙雙， 遲寶明， 姜紀沂. 地下水位對近震和遠震異常響應的比較：以汶川地震和蘇門答臘地震為例［J］. 吉林大學學報（地球科學版）， 2011， 26（1）： 523532， 565.

Lan Shuangshuang， Chi Baoming， Jiang Jiyi. Comparison of Groundwater Level Response to Near Earthquake and Distant Earthquake： Taking Wenchuan Earthquake and Sumatra Earthquake for Example［J］. Journal of Jilin University （Earth Science Edition）， 2011， 26（1）： 523532， 565.

［9］" 董殿偉， 林沛， 晏嬰. 北京平原地下水水位監測網優化［J］. 水文地質工程地質， 2007， 34（1）： 1018.

Dong Dianwei， Lin Pei， Yan Ying. Optimum Design of Groundwater Level Monitoring Network of Beijing Plain［J］. Hydrogeology amp; Engineering Geology， 2007， 34（1）： 1018.

［10］" 侯金武. 精心實施國家地下水監測工程［J］. 水文地質工程地質， 2014， 41（5）： 13.

Hou Jinwu. The Implementation of National Groundwater Monitoring Project［J］. Hydrogeology amp; Engineering Geology， 2014， 41（5）： 13.

［11］" 李文鵬.“水文地質與水資源調查計劃”進展［J］. 水文地質工程地質， 2022， 49（1）：16.

Li Wenpeng. Achievements of Investigation Program on Hydrogeology and Water Resources of CGS［J］. Hydrogeology amp; Engineering Geology， 2022， 49（1）： 16.

［12］" 李振函， 張春榮， 朱偉. 日照市沿海地區海水入侵現狀與分析［J］. 水文地質工程地質， 2009， 36（1）： 419.

Li Zhenhan， Zhang Chunrong， Zhu Wei. The Present and Analysis of the Sea Water Intrusion in Coastland of Rizhao［J］. Hydrogeology amp; Engineering Geology， 2009， 36（1）： 419.

［13］" 高銘駿， 汪旭， 孟玲瓏， 等. 沿海發展區海水入侵與土壤鹽漬化空間關聯解析［J］. 環境化學， 2022， 41（1）： 1426.

Gao Mingjun， Wang Xu， Meng Linglong， et al. Spatial Correlation Analysis of Seawater Intrusion and Soil Salinization in Coastal Developed Area［J］. Environmental Chemistry， 2022， 41（1）： 1426.

［14］" 談恒文， 張子軍， 安豐芹. 日照市東港區地下水環境質量現狀評價［J］. 地球學報， 2005， 26（1）： 518.

Tan Hengwen， Zhang Zijun， An Fengqin. Underground Water Environmental Quality Evaluation of Donggang District， Rizhao City， Shandong Province［J］. Acta Geochemistry Sinica， 2005， 26（1）： 518.

［15］" 楊鵬， 袁杰， 秦鵬. 日照市地下水動態特征及演化規律［J］. 地質學報， 2019， 15（1）： 100110.

Yang Peng， Yuan Jie， Qin Peng. Dynamic Characteristic and Evolution Law of Groundwater in Rizhao City， Shandong［J］. Acta Geological Sinica， 2019， 15（1）： 100110.

［16］" Obergfell C， Bakker M， Zaadnoorduk W J， et al. Deriving Hydrogeological Parameters Through Time Series Analysis of Groundwater Head Fluctuations Around Well Fields［J］. Hydrogeology Journal， 2013， 21（1）： 987999.

［17］" Moritz B S， Bartz-Beielstein T. Time Series Missing Value Imputation［J］. R Journal， 2015， 9（1）： 822.

［18］" 楊忠平， 盧文喜， 李平. 時間序列模型在吉林西部地下水動態變化預測中的應用［J］. 水利學報， 2005， 36（1）： 14751479.

Yang Zhongping， Lu Wenxi， Li Ping. Application of Time Series Model to Predict Groundwater Regime［J］. Journal of Hydraulic Engineering， 2005， 36（1）： 14751479.

［19］" 于寧莉， 易東云， 涂先勤. 時間序列中自相關與偏相關函數分析［J］. 數學理論與應用， 2007， 27（1）： 5457.

Yu Ningli， Yi Dongyun， Tu Xianqin. Analyze Auto-Correlations and Partial-Correlations Function in Time Series［J］. Mathematical Theory and Applications， 2007， 27（1）： 5457.

［20］" Manga M. On the Timescales Characterizing Groundwater Discharge at Springs［J］. Journal of Hydrology， 1999， 219（1）： 5669.

［21］" Wu Z， Huang N E， Long S R， et al. On the Trend， Detrending， and Variability of Nonlinear and Nonstationary Time Series［J］. Proceedings of the National Academy of Sciences of the United States of America， 2007， 104（1）： 1488914894.

［22］" Siino M W， Scudero S J， Cannelli V， et al. Multiple Seasonality in Soil Radon Time Series［J］. Scientific Reports， 2019， 9： 8610.

［23］" 閆佰忠，蓋俊百，王昕洲，等.石家莊滹沱河山前沖洪積扇地下水位動態演變特征及影響機制［J］.吉林大學學報（地球科學版）， 2023， 53（6）：18801891.

Yan Baizhong， Gai Junbai， Wang Xinzhou， et al. Characteristics of Groundwater Level Dynamic Evolution and Influence Factors in Mountain-Front Alluvial Plain of Hutuo River in Shijiazhuang［J］. Journal of Jilin University （Earth Science Edition）， 2023， 53（6）： 18801891.

［24］" Oygur T， Unal G. Vector Wavelet Coherence for Multiple Time Series［J］. International Journal of Dynamics and Control， 2021， 9（1）： 403409.

［25］" Lee E， Kim S. Wavelet Analysis of Soil Moisture Measurements for Hillslope Hydrological Processes［J］. Journal of Hydrology， 2019， 575（1）： 8293.

［26］" 孔慶敏， 王廣才， 史浙明. 云南地區震前地下流體異常特征統計分析［J］. 地震學報， 2018， 40（1）： 632645， 689.

Kong Qingmin， Wang Guangcai， Shi Zheming. Statistical Analysis of Pre-Seismic Anomalous Characteristics of Subsurface Fluids in Yunnan Region［J］. Acta Seismological Sinica， 2018， 40（1）： 632645， 689.

［27］" 仵彥卿. 地下水與地質災害［J］. 地下空間， 1999， 20（1）： 303310，316.

Wu Yanqing. Groundwater and Geological Hazards［J］. Underground Space， 1999， 20（1）： 303310，316.

［28］" 平建華， 李升， 欽麗娟. 地下水動態預測模型的回顧與展望［J］. 水資源保護， 2006， 22（1）： 1115.

Ping Jianhua， Li Sheng， Qin Lijuan. Review and Prospect of Dynamic Prediction Model for Groundwater［J］. Water Resource Protection，2006， 22（1）： 1115.

［29］" 梁四海， 徐德偉， 萬力. 黃河源區基流量的變化規律及影響因素［J］. 地學前緣， 2008， 15（1）： 1019.

Liang Sihai， Xu Dewei， Wan Li. Periodic Regularity of the Base Flow in the Headwater Region of the Yellow River and Affecting Factors［J］. Earth Science Frontiers， 2008， 15（1）： 1019.

［30］" 祁曉凡， 李文鵬， 李海濤. 濟南巖溶泉域地下水位、降水、氣溫與大尺度氣象模式的遙相關［J］. 水文地質工程地質， 2015， 42（6）： 1828.

Qi Xiaofan， Li Wenpeng， Li Haitao. Teleconnections Between Groundwater Levels， Precipitation， Air Temperature of the Jinan Karst Springs Watershed and Large Scale Climatic Patterns［J］. Hydrogeology amp; Engineering Geology， 2015， 42（6）： 1828.

［31］" 鳳蔚， 祁曉凡， 李海濤. 雄安新區地下水水位與降水及北太平洋指數的小波分析［J］. 水文地質工程地質， 2017， 44（1）： 819.

Feng Wei， Qi Xiaofan， Li Haitao. Wavelet Analysis Between Groundwater Level Regimes and Precipitation， North Pacific Index in the Xiongan New Area［J］. Hydrogeology amp; Engineering Geology， 2017， 44（1）： 819.

［32］" Yan R， Woith H， Wang R J， et al. Decadal Radon Cycles in a Hot Spring［J］. Scientific Reports， 2017， 7（1）： 1829.

［33］" Zhang S C， Shi Z M， Wang G C， et al. Groundwater Radon Precursor Anomalies Identification by Decision Tree Method［J］. Applied Geochemistry， 2020， 16（1）： 104696.

［34］" Torrence C. A Practical Guide to Wavelet Analysis［J］. Bulletin of the American Meteorological Society， 1998， 79（1）： 6178.

［35］" 徐澤華， 韓美， 路廣， 等. 山東省降雨侵蝕力與氣候指數關系研究［J］. 生態與農村環境學報， 2018， 34（10）： 880889.

Xu Zehua， Han Mei， Lu Guang， et al. Relationship Between Rainfall Erosivity in Shandong Province and Climate Indices［J］. Journal of Ecology and Rural Environment，2018， 34（10）： 880889.

［36］" 王善華， 黃鎮. 赤道東太平洋海溫、南方濤動和我國東部降水的譜分析［J］. 氣象科學， 1987， 5（2）： 6673.

Wang Shanhua， Huang Zhen. The Spectral Analysis of the Sea Surface Temperature in the Eastern Equatorial Pacific the Southern Oscillation and the Rainfall in the Eastern China［J］. Scientia Meteorologica Sinica，1987， 5（2）： 6673.

［37］" 朱繼前， 韓美， 徐澤華，等. 淮河流域不同量級降雨時空分布特征及其影響因素［J］. 水土保持研究， 2019， 26（4）： 8795.

Zhu Jiqian， Han Mei， Xu Zehua， et al. Temporal-Spatial Distribution Characteristics and Factors of Different Magnitude Rainfall in Huaihe River Basin［J］. Research of Soil and Water Conservation，2019， 26（4）： 8795.