信宜观测井位于粤西窿起带,该区地质构造复杂,总体格局为一组北东向的活动断层和一组北西向的活动断层互相切割,形成网格状构造,其中北东向的断层有合浦—玉林—梧州断裂、廉江—信宜断裂、吴川—四会断裂；北西向的断层有横县—博白—茂名断裂、高州—电白断裂、茶山—闸坡断裂。区域内历史上曾发生多次6级以上强震。

观测井孔所采取的基岩主要为花岗片麻岩变质后产生的混合岩,井孔总深101.90 m。以裂隙承压水为主,少量孔隙水。目前有水位、水温等观测项目,及辅助观测的气象三要素。

广东信宜井静水位测项是广东省内唯一观测资料地震预报效能评估A类井,现水位观测使用LN-3型数字化水位仪观测,就多年观测情况来看,信宜井具有较清晰的年变动态,观测曲线的固体潮效应显著,并且在2008年汶川地震, 2011年日本311地震后都出现了显著的震后效应。说明该观测井水位观测对区域应力场的变化反应较为灵敏。降雨是该井水位变化的主要影响因素。

图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（x₁）+ 0.27327x₂+(8.7727e-05)x₃+0.00070842x₄-0.0011766x₅

其中y:水位上高幅度(m);x₁:水位升高持续时间（d）;x₂:升高阶段最高水位(m);x₃: 上升开始前2d至结束前2d降雨量1(mm);x₄:升高期间最大单日降雨量(mm);x₅:上升开始前20天至开始前3d降雨量2(mm)。

对回归方程进行显著性检验p值为0.00514小于0.05 ,回归方程是显著的。图3为这12次水位升高的多元线性回归检验,最显著的一次高值异常是7号发生在苍梧地震前的2016年1月5日至2月2日,高值异常回落后约6个月发生了苍梧地震（震中距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

图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.018874x₁+ 0.00036502x₂+0.00076246x₃

y: 水位升高幅度(m)

x₁: 水位升高持续时间（天）

x₂: 升高时间段降雨量(m)

x₃: 最大单日降雨量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