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

楊鵬，張壽川，郭本力，等.山東日照地區(qū)不同類型含水層地下水位多年動態(tài)變化特征及其影響因素.吉林大學(xué)學(xué)報（地球科學(xué)版）， 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.

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

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

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

關(guān)鍵詞：地下水動態(tài); 連續(xù)小波變換; 小波相干; 時滯性；日照地區(qū)

doi：10.13278/j.cnki.jjuese.20220322

中圖分類號：P641.8

文獻標志碼：A

收稿日期：20221128

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

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

基金項目：日照市城市地質(zhì)調(diào)查項目（SDGP371100202102000475）；中國地質(zhì)調(diào)查局項目（DD202216772）；中央財政基本科研業(yè)務(wù)費資助項目（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" 引言

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

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

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

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

1 "區(qū)域概況

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

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

1.1" 水文地質(zhì)概況

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

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

荻水監(jiān)測井，位于日照市嵐山區(qū)荻水村，井深6.1 m，含水層類型為第四系松散巖類孔隙水，鉆孔柱狀圖如圖2a所示：0～1.1 m為耕植土；1.1～2.2 m為粉質(zhì)黏土，粉土與黏土體積比約為3∶7；2.2～6.1 m為強風化花崗巖。莒縣地震辦監(jiān)測井，位于日照市莒縣城陽街道，井深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為花崗巖，含水量較少。東莞中學(xué)監(jiān)測井，位于日照市莒縣東莞鎮(zhèn)，井深5.9 m，含水層為碳酸鹽巖類裂隙巖溶水，鉆孔柱狀圖如圖2c所示，0～0.8 m為耕植土，0.8～5.9 m為強風化砂頁巖。

1.2" 區(qū)域地下水動態(tài)

日照地區(qū)地下水動態(tài)監(jiān)測工作始于1993年，積

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

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

2" 研究方法

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

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

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

圖中彩色曲線為動態(tài)監(jiān)測示例數(shù)據(jù)。

2.1" 連續(xù)小波變換

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

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

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

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

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

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

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

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

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

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

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

2.2" 小波相干

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

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

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

2.3" 水位數(shù)據(jù)收集及處理

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

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

3" 結(jié)果

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

同周期小波系數(shù)動態(tài)變化圖可以反映同一周期內(nèi)振幅隨時間變化規(guī)律。

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

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

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

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

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

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

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

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

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

4" 討論

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

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

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

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

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

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

4.2" 地下水位多年周期變化動態(tài)變化影響因素

小波系數(shù)等值線圖和小波方差結(jié)果表明，3口地下水監(jiān)測井水位除呈現(xiàn)年際周期變化（荻水監(jiān)測井）和年際動態(tài)變化（莒縣地震辦監(jiān)測井和東莞中學(xué)

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

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

具有較強相關(guān)性，

對于19962015年間的11 a周期變化而言，二者相關(guān)系數(shù)為0.82，

也具有較強相關(guān)性；說明該井水位4～5 a

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

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

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

，對11 a周期變化而言，兩者相關(guān)系數(shù)為0.51，有一定相關(guān)性

，

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

和14 a強趨勢項

均為降雨影響所致。

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

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

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

4.3" 不同類型含水層水位多年動態(tài)變化特征

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

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

莒縣地震辦監(jiān)測井連通基巖裂隙含水層，為承壓含水層，該井連通的兩個含水層埋深分別為40.0～45.0 m和115.0～118.0 m，其中40.0～45.0 m含水層巖性為花崗巖，115.0～118.0 m含水層巖性為石英巖，裂隙發(fā)育，上覆弱透水層厚度較大，降雨補給地下水速度較慢，含水層滲透系數(shù)較?。灰虼耍鶐r裂隙含水層水位未呈現(xiàn)年際周期變化，僅呈現(xiàn)11 a周期變化和14 a強趨勢項變化。

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

5" 結(jié)論

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

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

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

參考文獻（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］" 常廷改， 胡曉. 水庫誘發(fā)地震研究進展［J］. 水利學(xué)報， 2018， 49（1）： 11091122.

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

［4］" 車用太， 魚金子， 劉成龍， 等. 水庫誘發(fā)地震地下水監(jiān)測網(wǎng)建設(shè)的若干問題探討［J］. 地震學(xué)報， 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］" 趙全升， 孔智涵， 胡舒婭，等.柴達木盆地馬海鹽湖地下鹵水地球物理探測及應(yīng)用［J］.吉林大學(xué)學(xué)報（地球科學(xué)版）， 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］" 何冠儒， 史浙明. 地下水對氣壓和固體潮響應(yīng)研究進展［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］" 蘭雙雙， 遲寶明， 姜紀沂. 地下水位對近震和遠震異常響應(yīng)的比較：以汶川地震和蘇門答臘地震為例［J］. 吉林大學(xué)學(xué)報（地球科學(xué)版）， 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］" 董殿偉， 林沛， 晏嬰. 北京平原地下水水位監(jiān)測網(wǎng)優(yōu)化［J］. 水文地質(zhì)工程地質(zhì)， 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］" 侯金武. 精心實施國家地下水監(jiān)測工程［J］. 水文地質(zhì)工程地質(zhì)， 2014， 41（5）： 13.

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

［11］" 李文鵬.“水文地質(zhì)與水資源調(diào)查計劃”進展［J］. 水文地質(zhì)工程地質(zhì)， 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］" 李振函， 張春榮， 朱偉. 日照市沿海地區(qū)海水入侵現(xiàn)狀與分析［J］. 水文地質(zhì)工程地質(zhì)， 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］" 高銘駿， 汪旭， 孟玲瓏， 等. 沿海發(fā)展區(qū)海水入侵與土壤鹽漬化空間關(guān)聯(lián)解析［J］. 環(huán)境化學(xué)， 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］" 談恒文， 張子軍， 安豐芹. 日照市東港區(qū)地下水環(huán)境質(zhì)量現(xiàn)狀評價［J］. 地球?qū)W報， 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］" 楊鵬， 袁杰， 秦鵬. 日照市地下水動態(tài)特征及演化規(guī)律［J］. 地質(zhì)學(xué)報， 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］" 楊忠平， 盧文喜， 李平. 時間序列模型在吉林西部地下水動態(tài)變化預(yù)測中的應(yīng)用［J］. 水利學(xué)報， 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］" 于寧莉， 易東云， 涂先勤. 時間序列中自相關(guān)與偏相關(guān)函數(shù)分析［J］. 數(shù)學(xué)理論與應(yīng)用， 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］" 閆佰忠，蓋俊百，王昕洲，等.石家莊滹沱河山前沖洪積扇地下水位動態(tài)演變特征及影響機制［J］.吉林大學(xué)學(xué)報（地球科學(xué)版）， 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］" 孔慶敏， 王廣才， 史浙明. 云南地區(qū)震前地下流體異常特征統(tǒng)計分析［J］. 地震學(xué)報， 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］" 仵彥卿. 地下水與地質(zhì)災(zāi)害［J］. 地下空間， 1999， 20（1）： 303310，316.

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

［28］" 平建華， 李升， 欽麗娟. 地下水動態(tài)預(yù)測模型的回顧與展望［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］" 梁四海， 徐德偉， 萬力. 黃河源區(qū)基流量的變化規(guī)律及影響因素［J］. 地學(xué)前緣， 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］" 祁曉凡， 李文鵬， 李海濤. 濟南巖溶泉域地下水位、降水、氣溫與大尺度氣象模式的遙相關(guān)［J］. 水文地質(zhì)工程地質(zhì)， 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］" 鳳蔚， 祁曉凡， 李海濤. 雄安新區(qū)地下水水位與降水及北太平洋指數(shù)的小波分析［J］. 水文地質(zhì)工程地質(zhì)， 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］" 徐澤華， 韓美， 路廣， 等. 山東省降雨侵蝕力與氣候指數(shù)關(guān)系研究［J］. 生態(tài)與農(nóng)村環(huán)境學(xué)報， 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］" 王善華， 黃鎮(zhèn). 赤道東太平洋海溫、南方濤動和我國東部降水的譜分析［J］. 氣象科學(xué)， 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.