浅谈西青区地面沉降的原因及治理措施.docxVIP

浅谈西青区地面沉降的原因及治理措施 天津市西青区位于市中心西南部。每年开采地下水约2000.10m，年平均沉降量为50m，中心城区年平均沉降量为25mm。中心城区南部（西青区交界处）年平均沉降量为30.50mm（西青区年平均沉降量）。结果表明，西青区的地面沉降对中心城区的地面沉降有明显影响。西青区2009年度平均沉降量43mm,2010年度平均沉降量53mm,其年沉降速率远远大于同期的天津市中心城区。因此本文以西青区郑庄子分层标监测数据和本区地下水开采量以及水位动态数据为分析基础,将观测孔水位动态曲线与模型计算曲线拟合对比,分析总结沉降量与地下水开采关系,为今后西青区控制地面沉降和调整地下水开采布局提供理论支撑。 1 西青区地下水开发模型建设 1.1 沉降速率分析 2004年7 月西青区郑庄子分层标建成,2006、2007、2009年沉降速率分别为64.43mm·a－１、57.06mm·a－１、47.27 mm·a－１,此处沉降速率较西青区全区平均沉降量要大。本次为研究地下水开采对地面沉降的影响,选取以西青郑庄子分层标为中心的数值模拟范围(图1),模拟区范围10km×10km×566m(长×宽×高);其模拟范围内沉降有明显的分带性,模拟范围内东北部沉降在30~50mm·a－１之间,西南部沉降大于50mm·a－１;与本区水位等值线比较后,初步判断该模拟范围地面沉降主要受本区地下水开采影响。 1.2 地下水流场及地下水均衡项 本区处于平原区,通过收集大量的地质资料,深入分析模拟区的含水层结构、地下水流场和流动特征、地下水均衡项,并建立水文地质概念模型,为建立数值模型提供依据。 (1) 改性组底界埋深 第一承压含水组底界埋深90~130 m,含水层普遍分布总厚10~20 m。导水系数20~61.87m３·d－１。第二承压含水组底界埋深170~185m,导水系数68.44~203.9 m３·d－１。第三承压含水组底界埋深300~323 m,含水层总厚30~55 m。导水系数28.13~127.9 m３·d－１。第四承压含水组底界埋深410~435 m,含水层总厚10~70 m。导水系数28.44~208.6 m３·d－１。第五承压含水组底界埋深566m,含水层总厚30~80m。导水系数64.36~215.6m３·d－１。 (2) 计算边界交换量 侧向边界:模拟区地势平坦但不是完整的水文地质单元,边界的流入流出量难以计算,因此模拟区侧向边界均设置为通用水头边界。通过给入每个应力期边界的水流参数、模型外的水位,模型自行计算流入或流出水量,从而调整边界上的水位,来计算边界的侧向交换量。 垂向边界:整个研究区在垂向上分为5个含水层组,根据含水层富水特征和开采条件来看,第二至第四含水层组研究意义较大,这几个含水层组在上部主要接受上一含水层组的越流补给,第五含水层组下部边界因下层没有分层标监测数据,不考虑沉降,模型中处理成隔水边界。第一含水层组接受大气降水入渗补给,蒸发和越流是其主要排泄方式,开采量小。从整个系统来看,垂向上只接受大气降水补给,蒸发、开采是其主要的排泄方式,系统内部各含水层组之间在垂向上通过越流产生部分交换量。 (3) 地下水系统的非均质性 从空间上看,模拟区地下水流整体上以水平运动为主、垂向运动为辅,地下水系统符合质量守恒定律和能量守恒定律,在常温常压下地下水运动符合达西定律,考虑含水层之间的流量交换,地下水运动可以概化为空间三维流。地下水系统的垂向运动是由层间水头差异引起的,地下水系统的输入输出随时间、空间变化,故地下水为非稳定流。参数随空间变化,体现了系统的非均质性,所以含水介质概化为非均质各向同性介质。综上所述,模拟区可概化成非均质、水平方向各向同性、垂向上存在变异、空间三维结构、非稳定地下水流系统,即地下水系统的概念模型。 (4) 地下水补气量参数 水文参数和水文地质参数为计算评价地下水的重要数据,也是影响评价结果的主要因素。用于地下水流模型的水文地质参数主要有两类,一类是用于计算各种地下水补排量的参数和经验系数,如大气降水入渗系数、灌溉入渗系数等,本次研究中这类参数主要影响降水入渗补给和蒸发量。另一类是含水层的水文地质参数,主要包括潜水含水层的渗透系数、给水度。此外,由于此次模型结合地面沉降模块,必须考虑弹性储水系数和非弹性储水系数。弹性储水系数主要参考经验得到。非弹性储水系数参考固结实验的计算结果,由经验和周边地区研究成果参考取得。 2 西青区地下水价值模型的建立 2.1 调查厚度的确定 计算区内层的划分按含水组划分5个层,含水组的厚度即为各层厚度,各层剖分采用相同的形式,主要是便于计算垂向各层之间的越流量。根据模拟区范围,将模拟区平面上分为50行,50 列,网格大小为200m×200m,总共2 500个单元格。