多元线性回归方法在地下水水位预测中的应用.docxVIP

多元线性回归方法在地下水水位预测中的应用陈志宏(北京市地质工程勘察院, 100037)摘要　在对已有长期开采资料的水源地进行资源评价时, 充分利用长系列观测资料, 可采用多元线性 回归方法对地下水水位的变化规律进行研究并对其变化趋势进行预测。 本文重点对相关因素的确定和多 元线性回归方法的应用进行了研究。关键词　多元线性回归, 地下水水位, 预测中图分类号:P 33319　文献标识码: A 　文章编号:100721903 ( 1999)0320020207对于已经长期开采的水源地而言,能否充分利用其开采量、 地下水水位动态、 降水量等长系列实际观测资料, 是决定资源评价及水位预报结果是否可靠的主要因素。北 京 市 第 三 水 厂 自 1959 年 开 采 至 今, 已 开 采 了 40 年。 该 厂 原 设 计 供 水 能 力 29万 m 3 /d, 在供水高峰期 (5～ 8 月) 已不能满足供水需要。 由于水厂位于严重超采区, 客观 条件已不允许增加地下水开采量, 从而提出了在年内不增加总开采量的条件下, 用水高峰期增加 10 万 m 3 /d 的开采量的调峰供水方案。 为此, 我们采用多元线性回归方法对地下水 水位的变化规律、 影响因素进行了研究, 并对调峰供水方案实施后、 第三水厂水源地地区第四系地下水水位可能发生的变化进行预测。1　水源地地下水动态特征水源地区地下水属于潜水动态区, 其主要补给源来自于侧向补给和大气降水补给。111　年内变化特征1) 最高水位滞后时间长根据区内长期观测孔多年情况观测资料分析, 由于地下水开采和大气降水补给减少的 影响, 水源地地区地下水水位的最低值一般出现在 6 月下旬, 之后, 由于大气降水的补给 作用, 地下水水位开始回升。 但与一般潜水动态不同的是, 年内最高水位不是出现在降水 集中的 7～ 8 月之后, 而是出现在 12 月底。 这是由于作为主要补给源的侧向流入量径流途 径较长, 使得水位回升出现滞后现象。 以区内长期观测孔 1994 年观测为例, 见表 1。2) 水位变幅小由于地下水以侧向流入补给为主, 使得年内水位变幅不大, 一般在 2～ 5m , 年变幅由 西北向东南, 呈现由大到小的变化规律。收稿日期: 1999206230112　年际变化特征地下水多年变化受人工开采、 大气降水、 永定河放水及回灌综合影响。表 1　1994 年水位峰值出现时间统计表孔号孔　位最高水位 (m )日　期最低水位 (m )日　期年变幅 (m )104农科院31159- - 12 3027191- - 06冉35151- - 12 3031111- - 06 204140100团校29133- - 12 2526137- - 06 302196通过四季青乡东冉村 98 号观测孔 (1959～ 1997 年) 动态曲线看出, 第四系地下水位变化可分三个时段; 1959～ 1971 年为地下水相对平稳时段; 1972～ 1988 年为地下水急剧下降 时段; 1989～ 1997 年为地下水相对回升时段 (图 1)。图 1　98 号长观孔水位动态观测曲线1) 相对平稳时段: 1959～ 1971 年, 第三水厂水源地地区开采主要以农业及生活用水为 主。三厂开采量 13 年累计 56167 万 m 3 , 平均每年开采 432015 万 m 3。但该时段年均降水量 为 63713mm , 接近建国以来平均降水量 638122mm 。 由于降水入渗的充分补给, 开采量较 小, 地下水位下降幅度不大, 累计水位下降 1193m , 平均每年下降 0116m 。2) 急剧下降时段: 1972～ 1988 年, 由于北京市第三水厂地区工业自备井、农灌井开采 量逐渐增加, 三厂开采量 17 年累计 165265133 万 m 3。平均每年开采 9721 万 m 3 , 年均开采 量 是 六 十 年 代 的 2125 倍。 时 段 年 平 均 降 水 量 为 57811mm , 比 多 年 平 均 降 水 量 少 了6011mm , 年均降水量减少 9142% 。由于开采量急剧增大和降水量的减少, 该时段地下水位 共下降了 14169m , 平均每年下降 0186m , 下降幅度是六十年代的 5176 倍。 造成水源地地区第四系地下水严重超采。3) 相对回升时段: 1989～ 1997 年, 三厂开采量 9 年累计万 m 3 , 平均每年开 采 9722164 万 m 3。 由于 1994～ 1996 年降水量平均达 69515mm ,