李利波 張 立 毛慧玲 羅睿潔

?

滇南地區(qū)降雨量與水位變化的關系——以澄江井、建水井為例1

李利波 張 立 毛慧玲 羅睿潔

（云南省地震局，昆明 650224）

2015年1月9日，云南地區(qū)出現(xiàn)大幅度水位同步轉折上升現(xiàn)象，而地下流體群體異常的聚集區(qū)可能是未來強震的孕育區(qū)，及時排除環(huán)境等干擾因素，對確定地震危險區(qū)很有幫助。本文選取澄江井和建水井的水位觀測數(shù)據(jù)，分析降雨量與水位變化的關系，認為滇南地區(qū)的水位年變化與雨季降雨量之間有較好的線性關系。本文對澄江井和建水井2012—2014年的水位異常信息進行提取，對其映震情況進行了總結，研究表明滇南地區(qū)降雨量對水位引起的異常，主要表現(xiàn)在變化幅度和變化時間等方面。此次出現(xiàn)的水位同步轉折上升現(xiàn)象多為降雨干擾影響。

降雨量 水位變化量 澄江井 建水井 滇南地區(qū)

引言

云南地區(qū)自2015年1月9日開始出現(xiàn)了一批短期同步地下流體宏微觀異常變化現(xiàn)象，如滇南地區(qū)的建水、澄江、瀘西水位，江川、曲江水溫，滇西地區(qū)洱源水位，滇西南地區(qū)的孟連水位、鎮(zhèn)源水溫、普洱離子等（圖1）。而云南地區(qū)2015年1月8、9日普降大雨、雪，該群體異常的出現(xiàn)，是降雨影響還是屬于前兆異常，引起許多學者的關注。

地下水位易受降雨、氣壓、固體潮、地應力等多種因素影響（張昭棟等，1989；Zhang等，1993；魏煥等，2003；車用太等，2006）。前人對氣壓、固體潮、地應力等因素對地下水位變化的影響進行了許多深入的研究，并嘗試利用相關數(shù)學模型方法，半定量或定量地分離和排除其引起的水位變化（魚金子等，1994；車用太等，2006；付虹等，2007，2014；楊竹轉等；2008；張立等，2009，2011；陳亮等，2012；張體移等，2012；趙棟，2013；孫小龍等，2013），而地下水位變化與降雨量的關系較為復雜，受多種因素影響，表現(xiàn)為水位變化與降雨同步、滯后或基本無變化，難以較好地定量分離降雨量對水位的影響。巨幅流體異常的出現(xiàn)是判定未來強震孕育區(qū)的依據(jù)之一，本文試圖分析降雨量與水位變化的關系，排除干擾因素，確認前兆異常，跟蹤分析強震孕育區(qū)。

1 數(shù)據(jù)選取和井孔概況

本文以位于小江斷裂帶、石屏-建水斷裂帶、紅河斷裂帶等交匯部位的澄江井和建水井為例（圖2），選取自2012年以來的年、月、日水位和降雨量數(shù)據(jù)，分析了滇南地區(qū)的降雨量與水位變化的相關關系。由于澄江井的日降雨量沒有觀測數(shù)據(jù)，以附近江川井（距離約46km）的數(shù)據(jù)代替。下面介紹3個井孔的基本情況。

澄江井（一中）位于澄江盆地西北部，小江斷裂帶南段西支的次級斷裂上，坐標為102°53′E，24°40′N，井深310m?；鶐r巖性為侏羅系下統(tǒng)馮家河組灰色、深灰色灰?guī)r和硅質條帶灰?guī)r，地下水為基巖裂隙水，水位主要受降雨量影響。

江川井（漁村）位于小江斷裂帶南段西支的次級斷裂上，坐標為102°45′E，24°20′N，井深180.42m。基巖巖性為新近系灰色泥巖夾細砂巖層。水位受降雨量、北側的河流、地下水等因素影響，變化較大。本文只分析降雨量對其水位的影響。

建水井（黃龍寺）位于石屏-建水斷裂的東端南緣，坐標為102°45′E，23°39′N，井深180m?；鶐r巖性為泥盆系中統(tǒng)灰色、深灰色泥巖夾砂礫巖層，地下水為基巖裂隙水，徑流途徑長，補給來源較遠，水位基本不受地下水影響。

2 資料分析

從數(shù)據(jù)中分別提取水位變化量較大時對應的降雨量，以及降雨量較大時對應的水位變化量，剔除空值和水位儀器校正值等人為因素造成的影響，統(tǒng)計分析降雨量與水位變化量間的關系，確認前兆異?，F(xiàn)象。

澄江井和建水井所處地區(qū)均屬亞熱帶大陸性氣候，氣候溫和濕潤。據(jù)2000—2014年的數(shù)據(jù)統(tǒng)計，建水井和澄江井的年降雨量分別為548—1163mm和693—1068mm。各年間的降雨量差異較大，降雨量≥100mm的月份主要集中在5—9月。

2.1 降雨量與水位年變化的關系

降雨對澄江井水位的影響主要表現(xiàn)為2012年在雨季為峰谷型（升—降—升），2013和2014年隨著雨季到來顯著上升，雨季過后下降（圖3）；水位年變幅度逐漸變大，由0.55m變大為1.93m。建水井2012年為谷峰型（降—升—降），2013年水位變化較?。ㄏ陆?.05m），2014年隨著雨季到來顯著上升，雨季過后下降（圖4）；水位年變幅度在2014年達到0.40m。

提取水位年變化量較大時對應的降雨量，對其關系進行一元回歸分析（圖5），得到澄江井的水位年變化與雨季降雨量的相關系數(shù)為0.9941；建水井的水位年變化與雨季降雨量的相關系數(shù)為0.9391。反映了滇南地區(qū)的水位年變化受雨季降雨影響較大，且相對穩(wěn)定。

2.2 降雨量與水位月變化的關系

降雨對水位月動態(tài)影響多集中表現(xiàn)在雨季（5—9月），水位隨著季節(jié)性降雨過程而起伏，差異較大（圖6和圖7）。從2012—2014年的統(tǒng)計數(shù)據(jù)可以看出，江川井有74次日降雨量≥10mm，最多一次為2014年6月6日，降水量71.8mm，而同期澄江井有35次日靜水位變化≥5cm；建水井日降雨量≥10mm有32次，而同期日靜水位變化≥5cm只有7次?？梢姲l(fā)生較大降雨不一定對水位變化有影響。

前人研究認為，氣壓、固體潮等因素對地下水位的影響一般在5cm左右（Zhang等，1993；魏煥等，2003；張昭棟等，1993；趙棟等，2013），因此本文選取水位月變化量≥5cm時對應的降雨量數(shù)據(jù)進行統(tǒng)計。對降雨量及水位月變化量進行一元回歸分析，得到澄江井的水位月變化與月降雨量的相關系數(shù)為0.1366（圖8（a）），建水井的水位月變化與月降雨量的相關系數(shù)為0.1277（圖8（b））。這反映了滇南地區(qū)的月降雨量對水位月變化受有一定的影響，但相關性不高。

2.3 日降雨量與水位變化的關系

水位日變化與降雨量的關系較為復雜。一般來說，降雨量較大時，水位會明顯上升，但有時候反而下降，表現(xiàn)為異常現(xiàn)象（圖9和圖10）。提取水位日變化量較大時對應的降雨量，對其關系進行一元回歸分析（圖11），得到澄江井的水位日變化與江川井日降雨量的相關系數(shù)為0.002；建水井的水位日變化與日降雨量的相關系數(shù)0.0018。這反映了滇南地區(qū)的日降雨量與水位日變化之間難以確立定量關系，隨機分布較大。但澄江井的數(shù)據(jù)多集中在降雨量5—30mm，水位日變化5—20cm，占統(tǒng)計數(shù)據(jù)的77.1%。

此外水位變化表現(xiàn)為與降雨過程同步、滯后或者基本無變化等，具有多變性，這可能與降雨時間和空間的不均勻性，以及方式的多變性有關。建水井在發(fā)生較大的降雨時，水位變化一般滯后1—2小時（圖12）；澄江井沒有整點觀測的降雨量數(shù)據(jù)，不作說明。

水位日變化與降雨量的關系復雜，還與觀測井的水文地質條件有關。澄江井和建水井基巖主要是灰?guī)r，內(nèi)部可能發(fā)育溶蝕現(xiàn)象，使其滲透性發(fā)生改變。此外周圍小江次級斷裂和石屏-建水斷裂帶、與降雨補給區(qū)的距離等因素可能影響水位的變化。

3 前兆異常信息提取

與降雨量相關的水位變化前兆異常主要包括水位變化幅度和變化時間等方面。剔除觀測井的水文地質背景、降雨時空分布、外界環(huán)境干擾等影響因素后（車用太等，2011），可作為前兆異常跟蹤。

3.1 水位變化幅度的異常

在雨季（5—9月），2014年澄江井的水位年變幅度為1.93m，偏離2012年和2013年同期幅度的103.3%和102.8%；建水井的水位年變幅度為0.40m，偏離2012年和2013年同期幅度的66.9%和170.2%，破壞正常年變規(guī)律，可提取為前兆異常信息。小江斷裂帶北段發(fā)生的2014年8月3日魯?shù)閟6.5級地震和紅河斷裂帶以西發(fā)生的2014年12月6日景谷s6.6級地震前，水位變化幅度較大（圖3及圖4）。

3.2 水位變化時間的異常

一般來說，澄江井和建水井的水位在雨季上升，旱季下降。2012年澄江井在7月呈下降狀態(tài)，建水井在5—7月呈下降狀態(tài)；而2013年建水井在6—12月呈趨勢下降狀態(tài)，破壞正常雨季上升規(guī)律，可提取為前兆異常信息，對應小江斷裂帶北段發(fā)生的2012年9月7日彝良s5.6、5.7級地震和紅河斷裂帶北段發(fā)生的2013年8月31日中甸s5.9級地震（圖3及圖4）。

2012—2014年澄江井和建水井的水位前兆異常特征見表1。2014年景谷s6.6級地震的水位異常距離發(fā)震時間稍長，大于3個月。其它地震，如彝良s5.6、5.7級地震、中甸s5.9級地震和魯?shù)?.5級地震的水位異常時間在3個月內(nèi)，屬于短期異常，這對確定地震危險區(qū)很有幫助。此外，水位發(fā)生較大變化時，若沒有降雨補給或降雨量很小，或早于降雨開始時間等，可作為前兆異常跟蹤，但要注意井區(qū)與降雨補給區(qū)的距離因素。

3.3 2015年1月云南地區(qū)水位同步轉折上升現(xiàn)象的分析

云南地區(qū)2015年1月8、9日普降大雨、雪，洱源、建水、澄江、瀘西、孟連地區(qū)這兩天的降雨量較大，分別為20.6mm、59.8mm、61mm、44.2mm、117.3cm。洱源、建水和孟連井水位在1月9日開始轉折上升，澄江井水位在1月9日先下降、10日才開始上升，孟連井水位在1月10日才顯著上升，反映了水位變化受降雨影響的多變性。洱源和建水井水位在1月10日開始下降并逐漸恢復到正常水平，澄江和瀘西井水位在1月12日轉折后未恢復（圖1）。孟連地區(qū)在1月17—20日又有一次降雨過程，降雨量達55.1mm，井水位在15日小幅轉折后，19日再次上升，24日開始下降并逐漸恢復到正常水平（圖13），反映了該群體水位同步轉折上升現(xiàn)象是由于降雨影響所致。

表1 水位前兆異常特征

Table 1 Characteristics of the precursory anomalies of water level

4 結論與討論

本文通過對澄江井和建水井的降雨量和水位觀測數(shù)據(jù)進行統(tǒng)計分析，得到以下認識：

（1）滇南地區(qū)的降雨量與水位變化的關系。水位年變化與雨季降雨量之間的線性關系較好，影響相對穩(wěn)定；水位月變化與月降雨量之間有一定的線性關系，但相關性一般；水位日變化與日降雨量之間則難以確定定量關系，隨機性較大。

（2）滇南地區(qū)的降雨量對水位變化的影響較為復雜，表現(xiàn)出與降雨過程同步、滯后或者基本無變化等多變性。這可能與降雨時間和空間的不均勻性、降雨方式的多變性、觀測井的水文地質條件、區(qū)域斷裂帶活動影響以及降雨補給區(qū)的距離等因素有關。

（3）對澄江井和建水井2012—2014年的水位異常信息進行提取，并對其映震情況進行了總結，認為滇南地區(qū)與降雨量相關的水位變化前兆異常，主要包括水位變化幅度和變化時間等方面。

（4）分析認為2015年1月云南地區(qū)水位同步轉折上升現(xiàn)象，是受降雨干擾影響。因此及時排除降雨等相關干擾因素，對確定地震危險區(qū)很有幫助。

此外，觀測井內(nèi)部水文地質條件的變化，井區(qū)降雨量數(shù)據(jù)的不準確以及水位觀測資料不連續(xù)等，可能對結果有影響，需要對其進行進一步分析研究。

車用太，魚金子，張大維，1993．降雨對深井水位動態(tài)的影響．地震，（4）：8—15．

車用太，魚金子，2006．地震地下流體學．北京：氣象出版社．

車用太，魚金子，劉成龍等，2011．判別地下水異常的干擾性與前兆性的原則及其應用實例．地震學報，33（6）：800—808．

陳亮，丁富雄，張蕾等，2012．建水地下流體監(jiān)測井水位驟降異常分析．地震研究，35（4）：471—476．

付虹，鄔成棟，劉強等，2007．印尼巨大地震引起的云南水位異常記錄及其意義．地震地質，29（4）：873—882．

付虹，鄔成棟，趙小艷等，2014．云南開遠井水位異常分析．地震學報，36（3）：292—298．

孫小龍，劉耀煒，晏銳，2013．云南姚安井2009年10月后水位下降的成因分析．地震學報，35（3）：410—420．

魏煥，張昭棟，耿杰等，2003．井水位氣壓加卸載響應比．西北地震學報，25（1）：82—85．

魚金子，車用太，張大維等，1994．降雨對湯坑群井水位動態(tài)影響的剖析．地震研究，17（3）：273—280．

楊竹轉，鄧志輝，劉春國等，2008．中國大陸井水位與水溫動態(tài)對汶川S8.0地震的同震響應特征分析．地震地質，30（4）：895—905．

趙棟，易立新，王廣才等，2013．地下水位中地震前兆信息提取方法研究．地震工程學報，35（2）：334—341．

張昭棟，鄭金涵，馮初剛，1989．井水位的固體潮效應和氣壓效應與含水層參數(shù)間的定量關系．西北地震學報，11（3）：47—52．

張昭棟，耿杰，高玉斌等，1993．井水位降雨影響的定量改正．地震學報，15（2）：202—207．

張立，趙洪聲，劉耀煒等，2009．云南會澤井水位與水溫相關關系及其變異的地震預測意義．地震研究，32（3）：228—230．

張立，趙洪聲，蘇有錦等，2011．云南強震的水位異常指標研究．高原地震，23（4）：10—12．

張體移，吳富煥，畢青等，2012．通海高大水位異常與地震分析．云南大學學報（自然科學版），34（S2）：86—92．

Zhang Z. D., Zheng J. H., 1993. Response of water level in deep well to earth tide and barometric pressure. Journal of Seismological Research, 16(4): 431—439.

Relationship between Rainfall and Water Level in Deep Wells: A Case Study of Chenjiang Well and Jianshui Well in South of Yunnan Province

Li Libo, Zhang li, Mao Huiling and Luo Ruijie

（Earthquake Administration of Yunnan Province, Kunming 650224, China）

Water levels were raised synchronously and sharply by a large margin in the Yunnan area since January 9, 2015. The related area with systematic anomaly of underground fluid may be the signal for strong earthquake preparation. However, such anomaly could also be caused by non-earthquake related factors. In this paper, we selected the observation data of Chenjiang well and Jianshui well, to study the relationship between rainfall and the variable quantity of water level. We found that there exists a clear linear relationship between rainfall and the variation of water level in south of the Yunnan area. Then, we extracted the abnormal information of water level of Chenjiang well and Jianshui well from 2012 to 2014, and analyzed the earthquake reflecting features. Finally, we concluded that abnormal changes of water level are mainly caused by rainfall variation.

Rainfall amount；The variable quantity of water level；Chenjiang well；Jianshui well；South of Yunnan area

1基金項目 云南省人民政府十項重點工程“云南大震短臨跟蹤工作”和中國地震局監(jiān)測預報司“2015年地下流體指標建設”資助

2015-09-25

李利波，男，生于1984年。碩士，工程師。主要從事地下流體分析和地震預測研究。E-mail：licags@126.com

李利波，張立，毛慧玲，羅睿潔，2017．滇南地區(qū)降雨量與水位變化的關系——以澄江井、建水井為例．震災防御技術，12（1）：230—239．doi：10.11899/zzfy20170124