廣東信宜井水位受降雨影響的回歸分析

劉 錦,陳大慶, 趙微微,鐘天任,楊 雪

（1.廣東省地震局,廣州 510070；2.大慶油田有限責任公司第六采油廠第四油礦地面工程組,黑龍江 大慶 163114；3.中山大學,廣州 510275 )

0 引言

地下水位作為流體學科重要的物理觀測量,是一個包含大氣降水、氣壓、固體潮、地應力等多種影響因素的復合參數[1-2]。前人對降雨、氣壓、固體潮、地應力等因素對水位影響進行了研究[3-4],而降雨是地下水位變化的主要影響因素,由于水位受降雨影響的形成機理較為復雜,定量排除降雨干擾是一個難題,一些學者對水位的降雨排除做過一些研究[5-6],嘗試利用地下水位同一年的峰值與谷值之差（即年變化幅度）和降雨量、地下水位谷值變化逐年差值與年降雨量、水位年變化量與雨季降雨量進行一元回歸相關分析,得到一些有益結果。本文在前人研究的基礎上根據水位長時間斷續升高和短時間單調升高兩種情況統計了降雨參數與水位升高幅度的關系,分別建立了兩種情況下降雨影響水位升高幅度的多元回歸方程,從而為信宜井水位數據的日常分析提供了一種基于回歸統計的定量排除降雨干擾的方法[7-8]。

1 信宜井觀測條件及區域地質概況

1.1 區域地質背景

信宜觀測井位于粵西窿起帶,該區地質構造復雜,總體格局為一組北東向的活動斷層和一組北西向的活動斷層互相切割,形成網格狀構造,其中北東向的斷層有合浦—玉林—梧州斷裂、廉江—信宜斷裂、吳川—四會斷裂；北西向的斷層有橫縣—博白—茂名斷裂、高州—電白斷裂、茶山—閘坡斷裂。區域內歷史上曾發生多次6級以上強震。

1.2 觀測條件

觀測井孔所采取的基巖主要為花崗片麻巖變質后產生的混合巖,井孔總深101.90 m。以裂隙承壓水為主,少量孔隙水。目前有水位、水溫等觀測項目,及輔助觀測的氣象三要素。

廣東信宜井靜水位測項是廣東省內唯一觀測資料地震預報效能評估A類井,現水位觀測使用LN-3型數字化水位儀觀測,就多年觀測情況來看,信宜井具有較清晰的年變動態,觀測曲線的固體潮效應顯著,并且在2008年汶川地震, 2011年日本311地震后都出現了顯著的震后效應。說明該觀測井水位觀測對區域應力場的變化反應較為靈敏。降雨是該井水位變化的主要影響因素。

2 回歸分析

水位觀測受降雨影響的形成機理較為復雜,不僅與降雨的分布、大小、速率等氣象因素有關,而且與降雨滲入補給區的地形地貌、植被及包氣帶的厚度、透水性、含水量等多種區域水文地質因素有關（流體百科）。實際觀測中取得上述參數較為困難,特別是區域水文地質因素。本文利用回歸分析統計了降雨量、降雨持續時間等氣象因素與水位變化的關系,建立了回歸方程。但對于區域水文地質因素無法及時取得其數據,且包氣帶的透水性、含水量可能受到前期降雨的影響而隨時間變化,因此文中分兩種情況來分析,一種是水位長時間斷續升高,另外一種是水位短時間單調升高。前者反映了雨季長時間降雨對水位的影響,后者反映了單次持續降雨對水位的影響[9-12]。

2.1 水位長時間斷續升高期間降雨影響的回歸分析

圖1中水位觀測曲線在雨季大多數出現了斷續的升高,表1中選取了2014—2017年7月信宜3號井水位的12次長時間上升區間（部分時段有小幅的回落）,選擇了升高期間降雨量（降雨量1）、降雨持續時間,最大單日降雨量3個降雨的氣象因素,包氣帶的透水性、含水量等水文地質因素可能和前期的降雨量及水位的高度有關,因此選取升高前期降雨量（降雨量2）、升高階段最高水位兩個因素來替代以上區域水文地質因素,用以上5個自變量對水位升高幅度因變量進行多元回歸分析。

圖1 信宜井水位與日降雨量Fig.1 Water level and daily rainfall in Xinyi well

表1 信宜井2014—2017年12個顯著上升期間及水位上升幅度檢驗Table 1 12 significant rising periods and water level rise tests in Xinyi well from 2014 to 2017

由于降雨對水位的影響在時間上有滯后,因此需要考慮對于升高期間降雨量（降雨量1）和升高前期降雨量（降雨量2）兩個自變量在時間段上的選擇,圖2是計算兩個降雨量和水位升高幅度的相關系數,當水位升高期間（降雨量1）選擇上升開始前2 d至結束前2 d相關系數最大為0.8279（圖2）,開始前3 d正相關系數下降的非?？?說明信宜井水位對降雨的滯后在2 d。升高前期降雨量（降雨量2）選擇上升開始前20 d至開始前3 d負相關系數最大為-0.5412。但從前期降雨開始時間選擇從前26 d到前19 d負相關系數變化不大。因此文中對于升高期間降雨量（降雨量1）選擇為上升開始前2 d至結束前2 d；升高前期降雨量（降雨量2）選擇上升開始前20 d至開始前3 d。

圖2 2個降雨量不同的選擇時間段與水位升高幅度的相關系數Fig.2 The correlation coefficient between rainfall and the rise of water level at two periods with different rainfall

表2為以上6個參數之間的相關系數。從中可以看到上升持續時間、升高階段最高水位、升高期間降雨量（降雨量1）、最大單日降雨量4個參數與水位升高幅度是正相關,其中上升持續時間、升高期間降雨量（降雨量1）、最大單日降雨量與水位升高幅度的相關系數幾乎達到或者超過了0.8。而升高前期降雨量（降雨量2）與水位升高幅度是負相關。

表2 水位升高幅度與持續時間、升高階段最高水位、水位升高區間降雨量（降雨量1）、最大單日降雨量、升高前期降雨量（降雨量2）的相關系數Table 2 The correlation coefficient between the increasing amplitude of water level and the duration,the highest water level in the rising stage, the rainfall in the rising interval (rainfall 1) ,the maximum daily rainfall and the rainfall in the early rising stage (rainfall 2)

y=1.3464+0.002464（x1）+ 0.27327x2+(8.7727e-05)x3+0.00070842x4-0.0011766x5

其中y：水位上高幅度(m);x1：水位升高持續時間（d）;x2：升高階段最高水位 (m);x3： 上升開始前2 d至結束前2 d降雨量1(mm);x4：升高期間最大單日降雨量(mm);x5：上升開始前20天至開始前3 d降雨量2(mm)。

對回歸方程進行顯著性檢驗p值為0.00514小于0.05 ,回歸方程是顯著的。圖3為這12次水位升高的多元線性回歸檢驗,最顯著的一次高值異常是7號發生在蒼梧地震前的2016年1月5日至2月2日,高值異?；芈浜蠹s6個月發生了蒼梧地震（震中距182 km）。

圖3 信宜井水位升高幅度與持續時間、升高階段最高水位、水位升高區間（降雨量1）、最大單日降雨量、升高前期（降雨量2）降雨量的多元回歸檢驗Fig.3 Multiple regression test on the range and duration of water level rise,the highest water level in the rising stage,the range of water level rise (rainfall 1) ,the maximum daily rainfall,and the rainfall in the early rising stage (rainfall 2) in Xinyi well

2.2 水位短時間單調升高降雨的回歸分析

圖4為2014—2017年131次水位連續上升3 d（包含3 d）以上,水位升高持續時間、水位升高期間降雨量（降雨量1）、最大單日降雨量3個參數和水位升高幅度之間的關系,以上3個參數和水位升高幅度有正相關的關系,相關系數分別為0.78、0.70 ,和0.68。降雨量1 選擇水位開始升高前一天到升高結束當天相關系數最大。由于水位短時間單調升高前期的降雨量（降雨量2）時間上多數會和前一次升高期的降雨量（降雨量1）重合,因此降雨量2對水位升高幅度的負相關系數并不好,另外升高階段最高水位相關系數也不好,因此在對水位短時間單調升高的回歸分析中舍棄了降雨量2和最高水位這兩個自變量。

圖4 信宜井水位短時單調升高幅度與降雨量、持續時間、最大單日降雨量的關系Fig.4 Relationship between short-term rise of water level and rainfall, duration and maximum daily rainfall in Xinyi well

表3 水位短時升高幅度與持續時間、水位升高區間降雨量（降雨量1）、最大單日降雨量的相關系數Table 3 The correlation coefficient of short-time rising range and duration, interval rainfall (rainfall 1) and maximum daily rainfall

以上參數建立的多元回歸方程為

y=-0.0066511+0.018874x1+ 0.00036502x2+0.00076246x3

y: 水位升高幅度(m)

x1: 水位升高持續時間（天）

x2: 升高時間段降雨量(m)

x3: 最大單日降雨量1(mm)

對回歸方程進行顯著性檢驗p值為2.44e-35小于0.05,回歸方程是顯著的。圖5為時間序列上這131次水位上升幅度與根據回歸方程擬合的上升幅度,圖6為兩者的差值,從中可以看到2016年7月31日蒼梧地震前2015年開始水位上升的幅度明顯偏高,而震后升高的幅度減小。

圖5 信宜水位上升幅度與回歸擬合得到的上升幅度Fig.5 The extent of water level rise in Xinyi and its regression fitting

圖6 信宜水位上升幅度與回歸擬合上升幅度差值Fig.6 Difference between water level rise and regression fitting rise in Xinyi

3 結論和討論

我國的流體前兆觀測大多受到降雨影響,因此在日常分析中如何排除降雨的影響一直是困擾前兆流體觀測的一個難題。由于地下水位變化與降雨量的關系較為復雜,不僅與降雨的分布、大小、速率等氣象因素有關,而且與降雨滲入補給區的地形地貌、植被及包氣帶的厚度、透水性、含水量等多種區域水文地質因素有關,這些因素的獲取都較為困難。本文在前人研究的基礎上根據水位長時間斷續升高和短時間單調升高兩種情況統計了降雨參數與水位升高幅度的關系,分別建立了兩種情況下降雨影響水位升高幅度的多元回歸方程,從而為信宜井水位數據的日常分析提供了一種基于回歸統計的定量排除降雨干擾的方法。

前一種情況水位長時間斷續升高反映了雨季長時間降雨對水位的影響,回歸分析中選擇了高期間降雨量（降雨量1）、降雨持續時間,最大單日降雨量3個降雨的氣象因素和升高前期降雨量（降雨量2）、升高階段最高水位兩個和水文地質因素的參量對水位升高幅度因變量進行多元回歸分析。而后一種情況水位短時間單調升高反映了單次持續降雨對水位的影響,在回歸分析中沒有選取升高前期降雨量（降雨量2）、升高階段最高水位兩個參量是由于在時間上多數會和前一次升高期的降雨量（降雨量1）重合,因此降雨量2和最高水位對水位升高幅度的負相關系數不好。

兩種回歸分析方法都顯示了在蒼梧地震前從2015年開始信宜井水位升高的幅度就偏高,震后偏低。這可能和地震前后區域應力場的調整有關。