频谱分析法在张掖城区地下水水位预报中的应用.pdfVIP

第 47卷第 2期 甘 肃 水 利 水 电技 术 Vo1．47．No．2 2011年 2月 GamuW aterConservancyandHydropowerTechnology Feb．，2011 · 水文水资源 · 频谱分析法在张掖城区地下水水位预报中的应用 曹永红 ，田辽西 (甘肃省地矿局 水文地质工程地质勘察院，甘肃 张掖 734000) 摘要：通过对张掖城区以南2km 的潜水观测孔 25a(1985—2009年)地下水年平均水位动态分析，拟合地下水水位 趋势方程，利用傅里叶级数将消除趋势项后的水位序列展开，分析谱参数之间的函数关系，建立地下水位预报模型， 并对模型进行显著性检验 ，可用于张掖城 区地下水水位预测。 关键词 ：地下水；预报模 型 ；张掖城 区 中图分类号：P641．74 文献标识码：B 文章编号：2095—0144(2011)02—003—03 1 引言 2 预报方法 张掖古称甘州 ，地处我国西北的河西走廊 中部 ，以其水 2．1 数据预处理 土资源丰富而闻名于世 ，素有 “塞上江南——金张掖”的美 选取张掖市城南梁家墩镇 田家小庙 2-3潜水观测孔历 称。历史上，甘州城区曾经是 半“城芦苇半城庙”的景观，地下 年动态观测资料 ，通过对 1985—2009年该孔 25a的年平均 高水位的现象一直持续到 20世纪70年代末 ，8O年代中期 地下水位观测资料和对应的年份序列 (1985年为第一年 ，之 以来，随着水利化程度的不断提高和水资源利用、分布格局 后依次类推 )进行曲线拟合 ，发现 3次曲线能很好的反映该 的变化，导致地下水总补给量不断减少和地下水开采量的持 孔水位变动趋势，利用最小二乘法法建立资料系列长期变化 续增加，张掖盆地地下水位处于区域性持续下降过程 ，城区 的趋势方程 ，求得各年份序列对应的趋势值，实测水位减去 一 带地下水位年均下降幅度为0．2～0．3m，到2003年 8月已 趋势值得到对应波动值见表 I_3、图2。 累计下降4～6m，著名的 “甘泉”和北大池(现甘泉公园)相继 经计算确定非线性 回归方程为 ： 干涸…。2003年 11月以来 ，甘州区大部分地段地下水位停止 ￡)=0．0007t3-0．0257t2+0．2079t+1480．7 (1) 下降开始出现上升趋势。特别是 2005年 9月以来，张掖市甘 式 中：F(t)——年平均水位 的趋势函数 ； 州城区及外围地下水位突然大面积上升 (图 1)，至2009年 r 年份，其中t=l为 1985年，t=5为 1989年，之后依 已累计上升 7-9rn。突出的表现在城区的地下室、地下仓库 次类推 。 渗水 ，积水深度达到 1．O～1．5m，部分平房出现裂缝或倒塌， 式(1)的判定系数R：=0．8875，t=l，2，3… 。可见该回归 居民住宅遭受严重威胁 ；农村乌江一三闸一带部分耕地被淹 方程相关度较高，能够很好地反应地下水位的变动趋势，故 没，直接经济损失 1500万元 ，已经形成环境地质灾害 ]。因 确定趋势项回归方程为 t)。 此，以系列地下水动态监测资料为基础 采用科学合理的方 2．2 傅里叶级数转换 法 ，建立地下水位时间序列数学模型 ，预测预报未来地下水 对地下水位变幅可假设其服从波函数 ： 位的变化趋势 ，对于张掖市制定合理的防灾减灾措施 ，实现 经济的可持续发展均有一定的现实意义。 (f)=