抽水地面沉降模型研究的回顾与展望.pdf

Journal l一3249／P 京新出报刊(2005)Z166号／2005／13(增刊)CNI ofEngineeringGeology工程地质学报 抽水地面沉降模型研究的回顾与展望 王光亚 (江苏省地质调查研究院南京210018) 摘要本文综合了地面沉降机理及最近几十年来用于模拟抽取地下水所致地面沉降的数学模型，分析了一维垂直渗流一垂 直固结模型、三维水流模型一垂直固结的准三维模型(两步模型)和考虑孔隙水压力变化和土体变形的耦合效应的完全耦合 模型，论述了每种类型模型的主要特点、不足。指出三维渗流和三维变形的耦合模型是地面沉降模型进一步研究的方向和将 会砸对的困难。 关键词地面沉降模型耦合渗流变形生命旋回模型 l 引 言 2地面沉降机理 过量抽取地下水引起的地面沉降是一种发展过 由于地下水的过量开采所致的地面沉降，目前 程缓慢，持续时间长，影响范围广的不可逆转的地质 普遍采用、也广为接受的解释是土力学中的有效应 灾害类型，其对地质环境的破坏是渐近的，对经济发 力原理。地下水的开采过程实际上是在井点的某个 展、城市建设和人民生活所构成的负面影响是极其 特定深度段给承压含水层(也可以是微承压水或潜 巨大的。 水含水层也一样)施加一个负压，使含水层的该井段 1936年，墨西哥的Guevas发表了《墨西哥的地 与周围出现压力差，周围土层中的水在这种压力差 面沉降问题》一文后，地面沉降受到国际社会的广泛 的驱使下向井运移。由于含水层的渗透系数总是远 关注(阎世骏等，1996)。此前专门的国际地面沉降 大于其顶底板的弱透水层、甚至是不透水层(当渗透 学术会议已召开了6届，2005年10月将在中国的上系数极小时)，所以抽水(负压)首先在含水砂层中 海召开。 形成一个测压水头降落漏斗，继而在邻近的顶、底板 目前，中国的长江、黄河、珠江三角洲、松辽平 土层中，在漏斗区影响范围内出某种程度的流体 原、环渤海地区及东南沿海平原的大多数城市，地面 (水)的渗流，虽然这个“漏斗”在初始时刻只是微 沉降正在继续大面积的发生、发展之中(段永侯， 型的，但是，漏斗和渗流则是实实在在存在着的。 1998)。地面沉降范围在扩大，危害在不断加剧。 随着抽水过程的持续，如果地下水的补给量不 地面沉降预测计算模型可分为两大类：一类为 足以补偿开采量，即，出现了地下水的过量开采，或 多孔介质流体渗流模型，由二部分构成，即水流(孔 者说，如果地下水的补给速率小于抽水速率，那么漏 隙水压力)预测模型，如相关法、解析法和数值法等， 斗区的测压水头持续下降、范围持续扩大将不可避 及土力学模型，如含水层的线弹性计算模型、弱透砂 免。含水砂层中的水不断向漏斗中心渗流，并不断经 层(粘性土层)的最终沉降量模型、太沙基固结模型、 由抽水井被开发利用，相邻各粘性土层中的孔隙水 流变固结模型、比奥(Blot)固结理论、弹塑性固结模 将不断向含水砂层中释放，孔隙水压力也随之降低。 型、回归模型及半理论、半经验模型以及最优化计算 K．Terzaghi于最早发现并提出了有效应力原理 法等。此外，还有水位一时间直接构成的沉降预测模 和固结理论，认为饱和土体任一平面上受到的总应 型，如多元统计方法计算沉降与渗透固结模型和流 力由有效应力和孔隙水压力两部分构成，其间关系 变固结模型；另二类为生命旋回模型。即，直接由沉 满足：a二=O一7w)z， 降量和时间的关系式构成，是从地面沉降的整个发 3 o G zz zz+IT。7z’ 展过程来考虑沉降的，有泊松旋回模型、Verhulst生 物模型和灰色模型。