皖23井水位的强震映震效应

祝怡红 周振贵 蒋世清 胡江辉

摘 要：地下水作为一种极为活跃的物质，参与了地震的孕育、发生和调整的全过程。而且在地震的前兆显示中起着重要的作用。本研究介绍了皖23井的基本情况，立足于皖23井水位仪对于强震映震现象，根据建立多元线性回归模型，得出了皖23井的强震映震效应与强震震级，震中距的关系，给出了回归方程，并通过检验得出方程显著。该方程表明皖23井对应强震的映震效应跟强震的震级成正比，跟井震距的对数成反比，此种定量分析也有助于我们以后研究皖23井的水位变化及地震异常情况。

关键词：井水位 映震效应 多元线性回归 皖23井

中图分类号：P3 文献标识码：A 文章编号：1674-098X（2015）07（b）-0011-02

Abstract：As an extremely active substances， the groundwater participates the whole process of the gestation， occurrence and adjustment in the earthquake . Also， water level of the well plays an important role in the earthquake precursory . Based on WAN 23 well water level for the strong earthquake effection， this paper introduces the basic situation of WAN 23 well， and set up the multiple linear regression model.Then we can found the relationship between the strong earthquake effection and the magnitude，and establish a regression equation to pass the test. This equation shows that the strong earthquake effection of WAN 23 well water level is proportional to the earthquake magnitude，and inversely proportional to logarithmic epicenter distance.Thus quantitative analysis also help us study the change of WAN 23 well water level and seismic anomalies.

Key Words：well water level；earthquake effection；multiple linear regression；Wan23 Well

利用地下水进行地震预报和研究，有比较明确的物理基础。地下水作为一种极为活跃的物质，参与了地震的孕育、发生和调整的全过程，故而地下水观测是一种重要的前兆观测手段。地震波可在固体介质中传播，当它传播到含水层时，由于含水层弹性形变，可引发孔隙压力的波动并表现为含水层中水位的震荡变化。而强震来临之后，地下活动的变化，也可能引发空隙压力的波动，从而引起水位的变化，对于井水位的强震应震效应分析研究有着重要的意义。安庆皖23井属于安徽省十五网络项目地下流体的改造工程，水位观测比较稳定，数据可靠。本研究通过皖23井自2008年汶川地震以来5次7级以上强震水位同震水震波最大振幅、震级、井震距的关系建立了多元线性回归模型。

1 皖23井概况

皖23井始建于1987年，位于安庆市东北郊十公里处的北东向头坡断裂上。该井位于大龙山与长江之间的冲积平原，西北边有大龙山，东边为石塘湖水域，地下水补给区主要来自大龙山，见图1。井孔地理坐标：北纬30.58度，东经：117.07度。该井位于空军水厂附近农田中，其东边距石塘湖约100米，湖水位对井水位有荷载效应，呈正相关，区内有肥中、肥西—韩摆渡断裂、郯庐断裂、宿松—枞阳断裂带通过。井深650.69米，井口标高8.7米，水温19℃。密封套管0米～293米，以下为裸孔，大都为石英正长斑岩。主要含水层位于350米～369.5米和399.5米～401米处，另外470.9米～550米之间有四段弱含水层。地下水属裂隙承压水，地表水对井水位影响不大。

2 皖23井水位强震映震效应

2.1 皖23井水位同震变幅

皖23井数字化观测虽然才三四年时间，但是仪器运行较为平稳，连续率和完整率都较高，数据资料也较为真实可靠。

在2007年6月1日至2011年9月30日的水位资料中，对五次强震同震现象进行分析（地震资料来自中国地震台网），统计结果见表1。

例如，我们绘制皖23井对于日本9.0级地震的同震水位，得到图2。

2.2 水位应震响应特征分析

计算方法：多元线性回归模型的建立

一般称

为高斯—马尔柯夫线性模型（K元线性回归模型），并简记为

称为回归平面方程。

用试验值（样本值）对未知参数和作点估计和假设检验，从而建立y与之间的数量关系，利用最小二乘法求的估计量（作离差平方和）：

选择使Q达到最小。

解得估计值。

得到的代入回归平面方程得：

称为经验回归平面方程，称为经验回归系数。

假设

当H0成立时，

如果F>F1-α（k，n-k-1），则拒绝H0，认为y与x1，…，xk之间显著地有线性关系；否则就接受H0，认为y与x1，…，xk之间线性关系不显著。

我们根据皖23井所体现的同震效应来分析强震震级（Ms）及井震距（D）对于水位振幅（ΔH）影响的关系。假设

代入表1数据，应用上述回归分析计算方法得：

相关系数r2=0.9669，接近于1，说明回归方程显著；

F>F1-α（k，n-k-1）时拒绝H0，F=29.1878，说明回归方程越显著；

与F对应的概率p=0.0002<α（α为显著性水平，设为0.05），拒绝H0，所以回归模型成立。

从残差图可以看出，除第二个数据（青海玉树地震）外，其余数据的残差离零点均较近，且残差的置信区间均包含零点，这说明回归模型

能较好的符合原始数据，而第二个数据可视为异常点。

由以上分析得出，强震震级越大，引发的水位震荡最大振幅越大，井震距越大，引发的水位震荡最大振幅越小。

3 结论

（1）地震波是弹性波，它在固体介质中可传播，当这种波传播到含水层时，由于含水层弹性变形，也可引发孔隙压力的波动并表现为含水层中水位的振荡变化，皖23井对应强震的映震效应跟强震的震级成正比，跟井震距的对数成反比，这种效应可以根据建立多元线性回归模型得出关系式：

。

（2）根据对近几年的水位数据分析得出，皖23井水位的同震能力较弱，由于样本的局限性，本研究对于映震效应的分析有待于进一步完善。

参考文献

[1]晏锐.影响井水位变化的几种因素研究[D].北京：中国地震局地震预测研究所，2008：59-63.

[2]马建英，刘学领，杜宗和.天津潘庄井水位异常与地震关系的分析[J].华北地震科学，2003，25（2）：298-306.

[3]尹宝军，马丽等.汶川8.0级地震及其强余震引起的唐山井水位同震响应特征分析[J].地震学报，2009，31（2）：195-204.