论文部分内容阅读
摘要:采用抽水、放水实验方法,研究地倾斜的变化量与地下水抽取量之间的关系、蓄水池池水变化量与地倾斜变化量之问的关系。通过对抽水实验数据分析研究可知:在放水速率(V放水<9m3/h)、总放水量(T<82m3)有限的前提下,弥渡地震台观测山洞正上方30m处的蓄水池水量变化不会引起地倾斜观测资料的明显变化;而地下水抽取是造成地倾斜异常变化的主要原因,且地倾斜异常变化在时间上与抽水作业(起止时间)保持高度同步性;抽水活动持续时间与地倾斜变化量是成非线性关系:随抽水时间延长,地倾斜变化速率开始较快,随后减小。但在开始抽水2小时时段内表现为近似线性关系,即NS测向YNS=-0.04553ms/min,EW测向VEW=-0.07971ms/min。通过频域分析可知:弥渡地倾斜观测受抽水干扰的异常信号的频率在0.0018-0.003Hz之间。由此,笔者设计了带阻滤波处理软件对抽水干扰信号频段进行处理,降低了抽水作业对弥渡地倾斜观测的影响。
关键词:抽水实验;地倾斜;抽水干扰;带阻滤波处理软件
中图分类号:P315.6
文献标识码:A
文章编号:1000-0666(2015)03-0416-10
0 引言
地壳在地球内力和外力作用下不断运动,由于地球的非均匀性,这种运动必然会在某些地区引起较大的地壳形变,在一定条件下导致地壳的突然破裂(地震发生)。地震与地壳形变关系最为直接,故而受到国内外地震预报探索者的普遍重视。地壳形变是一个复杂的系统,形变量的大小不仅与地壳内部的运动特征有关,而且与许多外界因素,如气温、气压、降雨量、蒸发量、地下水位、地表温度以及地下不同深度的地温等有关。囚抽水产生的地下水水位的动态变化会造成所在区域局部地表形变。作为观测精度很高(10-8-10 -10)的定点形变观测在采集形变前兆观测数据时,经常会记录到由于台站附近水井抽水而引起的异常现象。深井抽水会产生一定的地倾斜变化,其变化的幅度和持续时间由抽水点至观测仪器的方位、距离、抽水量以及地下水降深等因素决定。地震前兆观测的精度直接影响了地震前兆的分析,从而影响了地震分析预报的准确度。而抽水干扰就是影响地震前兆观测精度较大的因素之一。
倾斜观测仪主要包括水管倾斜仪、石英摆倾斜仪和钻孔垂直摆倾斜仪等。目前这些倾斜仪器的观测精度可达10-9,且连续性和稳定性较好。弥渡地震台地倾斜观测仪器为FSQ型水管倾斜仪,该观测于2013年12月在模拟观测基础上升级改造成DSQ型数字化观测。从2004年起,在弥渡地震台观测山洞旁的弥渡一中校园内挖了一口深层抽水井,并在山洞的顶上方修建了一个蓄水池,每天把井水抽到蓄水池供给弥渡一中使用。该学校的抽水、蓄水作业导致水管倾斜观测数据大幅变化,造成了观测质量下降以及地震前兆异常信号识别困难的问题。
弥渡地震台这样的观测环境无意中营造了一个天然的地倾斜受地下卸载与地面增压的小实验场,这样的实验场在全国同类观测中,并不多见。本文利用现有的实验场地,深入了解困扰弥渡地震台形变观测受抽水下扰的问题,了解其影响机制。进一步研究地下水抽取、蓄水池池水变化量与地倾斜变化之间的关系,为类似抽水、蓄水使观测场地增压、卸压,从而导致地倾斜观测量变化的相关研究提供参考。
1 台站观测环境概况
弥渡地震台位于弥渡冲击倾斜盆地东南部边沿近山地带,处在金沙江一红河断裂系上,位于红河断裂与程海断裂及楚雄一建水断裂的交汇点处(图1)。台站海拔1680m,地貌多为侵蚀低山地貌。
形变观测山洞是1985年在防空洞基础上改造而成的地震专用观测山洞。其中,各仪器室和引道总长288m,宽3m,高3m。北洞口至东西仪器室西端的进深为20m,西洞口至北东、南西仪器室的南西端进深为65m。平均覆盖厚度约为40m。山洞岩性以轻度风化玄武岩为主,较破碎。洞内日温差小于0.01℃,年温差小于0.4℃,相对湿度约为89%。经过多年观测使用,该山洞对天气变化的抗干扰能力较好。
2 抽水、放水实验
2.1 实验相关资料
2.1.1 抽水井及其附近水文地质资料
根据《弥渡县第一完全中学供水井竣工报告》的钻探揭露:本区范围内岩石种类单一,厚度人,地层主要由第四系浮土(主要为粉土、粉质粘土)及玄武岩层组成。其中粉土、粉质粘土层为相对隔水层,含水层主要为破碎的玄武岩层。玄武岩为中等风化,上部岩芯较完整,多为长柱状,含水性相对较差。中下部岩芯较为破碎,裂隙发育,地下水水蚀明显。其中64.1-72.0m为层问断层破碎带,可见黑色断层泥及滑动面。含水较为丰富,地下水类型为裂隙性地下水(图2)。
2.1.2 蓄水池设计资料
根据蓄水池设计图(图3)可知,蓄水池形状近似于长方体(19.5mx5.5m×1.62m)。距池底1.15m、1.41m处的池壁上分别安装有水位传感器a和b,当自动控制阀检测到蓄水池水位下降至a传感器位置时,水泵开始抽水作业;当水位超过b传感器位置时,水泵停止抽水作业。通过计算可知,在抽水实验模式下,蓄水池对山洞所产生的最火理论荷载压力变化在174m3水所产生的压力范围内;在自动抽水模式下,蓄水池对观测山洞的荷载压力范围是同定的,即蓄水池所产生的压力变化幅度始终保持在约28m3水所产生的压力范围内。
2.1.3 实验系统设计
为了摸清弥渡地倾斜观测受弥渡一中抽水作业的详细干扰情况,研究地下水抽取量与地倾斜变化量,蓄水池池水变化量与地倾斜变化量之间存在的定量关系。结合现有场地环境,通过与水井拥有方弥渡一巾沟通、协调后决定:在弥渡一中水管管网系统的不同位置布设水表,对水泵抽水量、蓄水池蓄水量、水池放水量进行实时监控。 2.2 实验过程及数据
实验场地布置如图4所示,与本次实验有关的水管线路共有3条,一条为抽水井一蓄水池之间的供水管道,另两条为蓄水池供水管道。通过在水管线路上安装水表(A、B、C、D)、闸阀(a、b、c、d)来对抽水作业、蓄水作业进行定量监控。与弥渡一中进行协调,在不影响弥渡一中日常用水需求的前提下,我们将实验日期选在学生离校的寒假、暑假期间进行。抽水实验时,通过分时段读取各水表读数,并结合蓄水池设计参数,可以计算出每一次手动抽水实验时的实验时段、实验时长、水泵抽水量、池进水量、池放水量、池水变化量、抽水速率、放水速率等参数(表1)。
3 实验数据定性分析
3.1 蓄水池池水变化量与地倾斜观测之间的关系
根据表2蓄水池放水实验数据统计情况,并结合弥渡台水管倾斜仪在蓄水池放水实验时段实测的观测数据(图5、6)进行仔细分析,在放水速率、总放水量有限的前提下,蓄水池放水时段前后数据连续性较好,并末出现明显数据异常情况。为此,通过初步研究可表明:在弥渡台形变观测山洞上的蓄水池蓄水、放水不会对弥渡台地倾斜观测造成明显干扰。
3.2 地下水抽取与地倾斜变化之间的关系
表3为本次抽水实验所获得的详细实验数据,将表3中的抽水实验数据与抽水实验时段弥渡台Dso型地倾斜仪实测数据(图7-9)进行对比分析,可发现如下规律:抽水作业可导致单位时间内的弥渡地倾斜变化量相对抽水前发生明显改变,在开始抽水和停止抽水两时刻均同时出现拐点;抽水所引起的地倾斜变化在时间上与抽水作业(起止时间)保持高度同步性。实验结果与李恩建等研究结果一致。抽水活动引起的地倾斜变化,还具有随抽水时间延长先加速、后减速的变化过程。上述变化过程可以用王庆良等提出的含水层水平变形机理作出较为合理的解释:即在抽水开始阶段,下伏断层尚处于抽水引起的地层径向拉张应变区域内,断层面上的正应力减小;随着抽水过程加长和内部径向挤压区的不断向外发展,下伏断层将逐渐处于抽水引起的径向挤压区域内,断层面上的正应力又将增大。
4 抽水对地倾斜影响的定量分析
在对实验数据定性分析之后,确定了弥渡台地倾斜异常变化是由深井抽水引起,且抽水时间与弥渡台地倾斜异常变化时间高度同步(图9),利用该实验成果,笔者以2014年1-8月弥渡台DSQ型水管倾斜仪观测数据随机选取的38个抽水时段原始数据为研究对象进行统计分析,以便为下一步定量研究抽水量与地倾斜变化量关系提供更加充分、客观、科学的数据支撑。详细统计情况如表4所示。
4.1 抽水干扰异常曲线的时域定量分析
使用Matlab数据处理软件对弥渡台DSQ型地倾斜仪原始观测数据进行处理:将表4中38个抽水时段的原始观测数据曲线汇总在一起,可得到弥渡台2014年1-8月巾随机的38个抽水时段地倾斜实测曲线图(图10);将38个抽水时段的原始观测数据进行平均计算,可得到在开始抽水2h内弥渡台平均地倾斜变化量曲线图(图lla);最后,将抽水时段平均地倾斜变化量进行分钟值差分处理,可得到平均地倾斜变化率统计图(图llb)。
从图8可看出弥渡台整个抽水时段的地倾斜变化量成非线性变化关系。但在抽水2h内,其因抽水而导致的地倾斜变化量几乎成线性关系(图lla),其地倾斜变化率基本上为一个同定值(图llb)。通过计算可得到弥渡地倾斜观测受抽水下扰引起的地倾斜变化率为:NS测向VNS=-0.04553ms/min,EW测向VFW=-0.07971ms/min.
4.2 抽水干扰异常曲线的频域分析
结合弥渡一巾日用水量(平均90m3/d)统计情况,可初步计算出弥渡一中在日常内动抽水模式下,平均抽水时长一般在2.5-4.5h范围内(表3的实际抽水统计情况可验证),而日抽水次数与弥渡一中日用水量、用水高峰期分布有关,据表3统计可知弥渡一巾抽水外日抽水次数(即日抽水下扰时段数)约为2.5次。假设数据恢复时段时长与抽水时段时长相等的情况下,弥渡台受抽水干扰的地倾斜异常信号的周期应是平均抽水时长的两倍,即5-9h,转换成频率可知弥渡台受抽水下扰的地倾斜信号频率范围在0.0018-0.003Hz之间,故受抽水干扰的地倾斜异常信号频率范围处于地倾斜同体潮的1/3日潮汐附近,通过实测曲线频谱与同体潮曲线频谱对比分析可得到验证(图12)。其中,弥渡台地倾斜同体潮理论值数据是采用EIS2000专业软件计算所得。为此,给我们提供了一种从原始信号中去除抽水下扰信号的方法。笔者通过设计带阻滤波处理软件把抽水干扰信号频段滤除,带阻滤波处理软件的设计原理如下:把受抽水下扰的地倾斜观测信号通过快速Fourier变换对其进行频率分析,然后在频率域中将抽水干扰信号的频率成分(0.0018-0.003Hz)的振幅置零,接着运用Fourier逆变换(IFFT)到时间域从而达到滤波的效果。滤波效果如图13所示。详细的matlab滤波程序如下:加载弥渡台地倾斜观测数据间序列及频率序列
要滤去频率的上限和下限
yy=zeros(size(y));%设置与y相同元素的数组
form=0:N-1;%将频率落在该频率范围及其人于Nyquist频率的波滤去
if(m/(N*dt)>fl&m/(N*dt) yy(m+l)=0;%置在此频率范围内的振动振幅为零
else
yy(m+1)=y(m+1);%其余频率范围的振动振幅不变;%运用IFFT变换同时间域并绘制滤波后的数据曲线
5 结论
经过对抽水实验数据的分析研究可知:在放水速率(V放水<9m3/h)、总放水量(T<82m3)有限的前提下,弥渡台观测山洞正上方30m处的蓄水池池水变化不会对弥渡地倾斜观测造成明显干扰;而抽水作业是造成地倾斜异常变化的主要原因,且在时间上与抽水作业保持高度同步性。此外,抽水活动引起的地倾斜变化是成非线性关系的,具体表现为随抽水时间延长先加速、后减速的变化过程。但在开始抽水2h内表现为近似线性关系,即NS测向vNS=-0.04553ms/min,EW测向VEW=-0.07971ms/min。通过频域分析可知:弥渡地倾斜观测受抽水下扰的异常信号的频率在0.0018-0.003Hz之间。由此,笔者设计了带阻滤波处理软件对抽水干扰信号频段进行处理,大人降低了抽水作业对弥渡地倾斜观测的影响,其效果较为明显。
关键词:抽水实验;地倾斜;抽水干扰;带阻滤波处理软件
中图分类号:P315.6
文献标识码:A
文章编号:1000-0666(2015)03-0416-10
0 引言
地壳在地球内力和外力作用下不断运动,由于地球的非均匀性,这种运动必然会在某些地区引起较大的地壳形变,在一定条件下导致地壳的突然破裂(地震发生)。地震与地壳形变关系最为直接,故而受到国内外地震预报探索者的普遍重视。地壳形变是一个复杂的系统,形变量的大小不仅与地壳内部的运动特征有关,而且与许多外界因素,如气温、气压、降雨量、蒸发量、地下水位、地表温度以及地下不同深度的地温等有关。囚抽水产生的地下水水位的动态变化会造成所在区域局部地表形变。作为观测精度很高(10-8-10 -10)的定点形变观测在采集形变前兆观测数据时,经常会记录到由于台站附近水井抽水而引起的异常现象。深井抽水会产生一定的地倾斜变化,其变化的幅度和持续时间由抽水点至观测仪器的方位、距离、抽水量以及地下水降深等因素决定。地震前兆观测的精度直接影响了地震前兆的分析,从而影响了地震分析预报的准确度。而抽水干扰就是影响地震前兆观测精度较大的因素之一。
倾斜观测仪主要包括水管倾斜仪、石英摆倾斜仪和钻孔垂直摆倾斜仪等。目前这些倾斜仪器的观测精度可达10-9,且连续性和稳定性较好。弥渡地震台地倾斜观测仪器为FSQ型水管倾斜仪,该观测于2013年12月在模拟观测基础上升级改造成DSQ型数字化观测。从2004年起,在弥渡地震台观测山洞旁的弥渡一中校园内挖了一口深层抽水井,并在山洞的顶上方修建了一个蓄水池,每天把井水抽到蓄水池供给弥渡一中使用。该学校的抽水、蓄水作业导致水管倾斜观测数据大幅变化,造成了观测质量下降以及地震前兆异常信号识别困难的问题。
弥渡地震台这样的观测环境无意中营造了一个天然的地倾斜受地下卸载与地面增压的小实验场,这样的实验场在全国同类观测中,并不多见。本文利用现有的实验场地,深入了解困扰弥渡地震台形变观测受抽水下扰的问题,了解其影响机制。进一步研究地下水抽取、蓄水池池水变化量与地倾斜变化之间的关系,为类似抽水、蓄水使观测场地增压、卸压,从而导致地倾斜观测量变化的相关研究提供参考。
1 台站观测环境概况
弥渡地震台位于弥渡冲击倾斜盆地东南部边沿近山地带,处在金沙江一红河断裂系上,位于红河断裂与程海断裂及楚雄一建水断裂的交汇点处(图1)。台站海拔1680m,地貌多为侵蚀低山地貌。
形变观测山洞是1985年在防空洞基础上改造而成的地震专用观测山洞。其中,各仪器室和引道总长288m,宽3m,高3m。北洞口至东西仪器室西端的进深为20m,西洞口至北东、南西仪器室的南西端进深为65m。平均覆盖厚度约为40m。山洞岩性以轻度风化玄武岩为主,较破碎。洞内日温差小于0.01℃,年温差小于0.4℃,相对湿度约为89%。经过多年观测使用,该山洞对天气变化的抗干扰能力较好。
2 抽水、放水实验
2.1 实验相关资料
2.1.1 抽水井及其附近水文地质资料
根据《弥渡县第一完全中学供水井竣工报告》的钻探揭露:本区范围内岩石种类单一,厚度人,地层主要由第四系浮土(主要为粉土、粉质粘土)及玄武岩层组成。其中粉土、粉质粘土层为相对隔水层,含水层主要为破碎的玄武岩层。玄武岩为中等风化,上部岩芯较完整,多为长柱状,含水性相对较差。中下部岩芯较为破碎,裂隙发育,地下水水蚀明显。其中64.1-72.0m为层问断层破碎带,可见黑色断层泥及滑动面。含水较为丰富,地下水类型为裂隙性地下水(图2)。
2.1.2 蓄水池设计资料
根据蓄水池设计图(图3)可知,蓄水池形状近似于长方体(19.5mx5.5m×1.62m)。距池底1.15m、1.41m处的池壁上分别安装有水位传感器a和b,当自动控制阀检测到蓄水池水位下降至a传感器位置时,水泵开始抽水作业;当水位超过b传感器位置时,水泵停止抽水作业。通过计算可知,在抽水实验模式下,蓄水池对山洞所产生的最火理论荷载压力变化在174m3水所产生的压力范围内;在自动抽水模式下,蓄水池对观测山洞的荷载压力范围是同定的,即蓄水池所产生的压力变化幅度始终保持在约28m3水所产生的压力范围内。
2.1.3 实验系统设计
为了摸清弥渡地倾斜观测受弥渡一中抽水作业的详细干扰情况,研究地下水抽取量与地倾斜变化量,蓄水池池水变化量与地倾斜变化量之间存在的定量关系。结合现有场地环境,通过与水井拥有方弥渡一巾沟通、协调后决定:在弥渡一中水管管网系统的不同位置布设水表,对水泵抽水量、蓄水池蓄水量、水池放水量进行实时监控。 2.2 实验过程及数据
实验场地布置如图4所示,与本次实验有关的水管线路共有3条,一条为抽水井一蓄水池之间的供水管道,另两条为蓄水池供水管道。通过在水管线路上安装水表(A、B、C、D)、闸阀(a、b、c、d)来对抽水作业、蓄水作业进行定量监控。与弥渡一中进行协调,在不影响弥渡一中日常用水需求的前提下,我们将实验日期选在学生离校的寒假、暑假期间进行。抽水实验时,通过分时段读取各水表读数,并结合蓄水池设计参数,可以计算出每一次手动抽水实验时的实验时段、实验时长、水泵抽水量、池进水量、池放水量、池水变化量、抽水速率、放水速率等参数(表1)。
3 实验数据定性分析
3.1 蓄水池池水变化量与地倾斜观测之间的关系
根据表2蓄水池放水实验数据统计情况,并结合弥渡台水管倾斜仪在蓄水池放水实验时段实测的观测数据(图5、6)进行仔细分析,在放水速率、总放水量有限的前提下,蓄水池放水时段前后数据连续性较好,并末出现明显数据异常情况。为此,通过初步研究可表明:在弥渡台形变观测山洞上的蓄水池蓄水、放水不会对弥渡台地倾斜观测造成明显干扰。
3.2 地下水抽取与地倾斜变化之间的关系
表3为本次抽水实验所获得的详细实验数据,将表3中的抽水实验数据与抽水实验时段弥渡台Dso型地倾斜仪实测数据(图7-9)进行对比分析,可发现如下规律:抽水作业可导致单位时间内的弥渡地倾斜变化量相对抽水前发生明显改变,在开始抽水和停止抽水两时刻均同时出现拐点;抽水所引起的地倾斜变化在时间上与抽水作业(起止时间)保持高度同步性。实验结果与李恩建等研究结果一致。抽水活动引起的地倾斜变化,还具有随抽水时间延长先加速、后减速的变化过程。上述变化过程可以用王庆良等提出的含水层水平变形机理作出较为合理的解释:即在抽水开始阶段,下伏断层尚处于抽水引起的地层径向拉张应变区域内,断层面上的正应力减小;随着抽水过程加长和内部径向挤压区的不断向外发展,下伏断层将逐渐处于抽水引起的径向挤压区域内,断层面上的正应力又将增大。
4 抽水对地倾斜影响的定量分析
在对实验数据定性分析之后,确定了弥渡台地倾斜异常变化是由深井抽水引起,且抽水时间与弥渡台地倾斜异常变化时间高度同步(图9),利用该实验成果,笔者以2014年1-8月弥渡台DSQ型水管倾斜仪观测数据随机选取的38个抽水时段原始数据为研究对象进行统计分析,以便为下一步定量研究抽水量与地倾斜变化量关系提供更加充分、客观、科学的数据支撑。详细统计情况如表4所示。
4.1 抽水干扰异常曲线的时域定量分析
使用Matlab数据处理软件对弥渡台DSQ型地倾斜仪原始观测数据进行处理:将表4中38个抽水时段的原始观测数据曲线汇总在一起,可得到弥渡台2014年1-8月巾随机的38个抽水时段地倾斜实测曲线图(图10);将38个抽水时段的原始观测数据进行平均计算,可得到在开始抽水2h内弥渡台平均地倾斜变化量曲线图(图lla);最后,将抽水时段平均地倾斜变化量进行分钟值差分处理,可得到平均地倾斜变化率统计图(图llb)。
从图8可看出弥渡台整个抽水时段的地倾斜变化量成非线性变化关系。但在抽水2h内,其因抽水而导致的地倾斜变化量几乎成线性关系(图lla),其地倾斜变化率基本上为一个同定值(图llb)。通过计算可得到弥渡地倾斜观测受抽水下扰引起的地倾斜变化率为:NS测向VNS=-0.04553ms/min,EW测向VFW=-0.07971ms/min.
4.2 抽水干扰异常曲线的频域分析
结合弥渡一巾日用水量(平均90m3/d)统计情况,可初步计算出弥渡一中在日常内动抽水模式下,平均抽水时长一般在2.5-4.5h范围内(表3的实际抽水统计情况可验证),而日抽水次数与弥渡一中日用水量、用水高峰期分布有关,据表3统计可知弥渡一巾抽水外日抽水次数(即日抽水下扰时段数)约为2.5次。假设数据恢复时段时长与抽水时段时长相等的情况下,弥渡台受抽水干扰的地倾斜异常信号的周期应是平均抽水时长的两倍,即5-9h,转换成频率可知弥渡台受抽水下扰的地倾斜信号频率范围在0.0018-0.003Hz之间,故受抽水干扰的地倾斜异常信号频率范围处于地倾斜同体潮的1/3日潮汐附近,通过实测曲线频谱与同体潮曲线频谱对比分析可得到验证(图12)。其中,弥渡台地倾斜同体潮理论值数据是采用EIS2000专业软件计算所得。为此,给我们提供了一种从原始信号中去除抽水下扰信号的方法。笔者通过设计带阻滤波处理软件把抽水干扰信号频段滤除,带阻滤波处理软件的设计原理如下:把受抽水下扰的地倾斜观测信号通过快速Fourier变换对其进行频率分析,然后在频率域中将抽水干扰信号的频率成分(0.0018-0.003Hz)的振幅置零,接着运用Fourier逆变换(IFFT)到时间域从而达到滤波的效果。滤波效果如图13所示。详细的matlab滤波程序如下:加载弥渡台地倾斜观测数据间序列及频率序列
要滤去频率的上限和下限
yy=zeros(size(y));%设置与y相同元素的数组
form=0:N-1;%将频率落在该频率范围及其人于Nyquist频率的波滤去
if(m/(N*dt)>fl&m/(N*dt)
else
yy(m+1)=y(m+1);%其余频率范围的振动振幅不变;%运用IFFT变换同时间域并绘制滤波后的数据曲线
5 结论
经过对抽水实验数据的分析研究可知:在放水速率(V放水<9m3/h)、总放水量(T<82m3)有限的前提下,弥渡台观测山洞正上方30m处的蓄水池池水变化不会对弥渡地倾斜观测造成明显干扰;而抽水作业是造成地倾斜异常变化的主要原因,且在时间上与抽水作业保持高度同步性。此外,抽水活动引起的地倾斜变化是成非线性关系的,具体表现为随抽水时间延长先加速、后减速的变化过程。但在开始抽水2h内表现为近似线性关系,即NS测向vNS=-0.04553ms/min,EW测向VEW=-0.07971ms/min。通过频域分析可知:弥渡地倾斜观测受抽水下扰的异常信号的频率在0.0018-0.003Hz之间。由此,笔者设计了带阻滤波处理软件对抽水干扰信号频段进行处理,大人降低了抽水作业对弥渡地倾斜观测的影响,其效果较为明显。