天津市地面沉降预测方法.pdfVIP

-186- 技术方法 文章编号：1007-1903(2006)Z1-0186-13 天津市地面沉降预测方法 陈文杰。崔晓东 (天津市地质调查研究院，天津，300191) 擒要：本文介绍了在天津市地面沉降研究中常用的预测方法。包括数值模拟方法和数理统计方法以及神经 网络预测方法。数值模型为MODFLOW有限差分模型与地层压缩模型(IDP)的耦合模型。统计模型主要有 回归分析模型以及灰色模型和VERHUI．ST模型等。神经网络模型主要介绍BP神经网络。文中给出了计算 实例。并对采用的各种方法的适用性进行综合评速。 关键词：地面沉降；MODFLOW；IDP；计模型；灰色系统；VER．HULST模型；BP神经网络模型 中图分类号：G241 文献标识码：A 天津市从1959年发现地面沉降，1986年开始治理并进行大范围的地面沉降监测工作， 目前已经积累了20多年的地下水水位和地下水开采量监测资料以及近20年的大地水准测 量资料。为了更好地控制地面沉降对地面沉降工作进行规划，地面沉降预测方法研究很有必 要。本文主要介绍天津市地面沉降预测预报研究工作中常用的预测方法。 常用的预测方法包括数值模型法以及统计方法。数值模型采用MODFLOW有限差分模 型与地层压缩模型(IDP)的耦合模型。计算区域为天津市区，根据天津市区的水文地质特征 将第四系含水岩组在垂向上概化为一个潜水一微承压水含水组、五个承压含水组。 统计模型主要有回归分析模型以及灰色模型和VERHULST模型等。其中回归模型包括 一元线性回归、一元非线性回归以及多元回归模型等。神经网络模型主要是BP网络(Back— Network)。 Propagation 1 数值模型 运用地下水运动与地层压缩的耦合模型预测地面沉降，模型由描述地下水三维运动的 MODFLOW与描述地层压缩的IDP耦合而成 1．1 MODFLOW地下水三维运动数学 模拟第四系因开采地下水引起的地面沉降，面积546平方公里，模拟深度550米。共划 分为一个潜水一微承压水含水组和五个承压含水组，各含水组之间为粘性土弱透水层。地下 水不仅有水平运动，也有垂向运动，砂层的弹性释水和粘性土的压密释水同时并存。 地下水拟三维运动方程为： ％砉(鲁)彩n寿＼O别h,I以+o’r-hi)KJ～m,=誓 式中：h一含水层水头，m；j1’一相邻含水层水头，米；“含水层渗透系数，m／d；肛源汇项，1／d；m’一 作者简介：陈文杰i天津市地质调查研究院副总工程师、信息中心主任，电话：022邮 箱：TJDDY@EYOU．COM；崔晓东，天津市地质调查研究院院长。 技术方法 一187一 弱透水层厚度，．m；k’一弱透水层垂向渗透系数，m／d；ss一含水层比弹性释水系数，1／m，j一 1，2，3，4，5。 1．2地层压缩模型 粘性土孔隙水压力： 他=器 粘性土变形量： AS=APxm，x△Z=Ah,r．xmo×衄 式中：△％一时段长度内粘性土孔隙水头变化，m；V_c-粘性土垂向传导力，m硼；危。一粘性土孔隙 水头，m；^，一含水层水头，m；一时段长度，d；驴水的溶重，g／m3；脚广粘性土垂向压缩系数， 段内粘性土垂向压缩量，mo． 1．3地面沉降模拟与预测 利用1981—1999年的地面水准测量资料和分层标监测资料对模拟结果进行了检验(图 1-6)，根据目前的开采现状，利用检验后的模型进行了沉降预测(图7-8)。 ∞M∞嘲■^POF SIMULATION鲫Bsl睫嗍∞¨∞3X呻 图1 1959-1983年模拟地面沉降等值线㈤ 图2 1959—1983年实测地面沉降等值线(m) -188- 技术方法 圈419591993年实舅地面沉悻等值线(n1) 圈3l’转～19蚣年模拟地面沉降暮值线fml 图3