基于最小二乘的地面沉降预测评估方法研究.docxVIP

基于最小二乘的地面沉降预测评估方法研究 1 关于地面沉降量的建模 由于大量地下水的开采，导致地面沉降的风险之一是建筑工地风险风险评估的重要内容之一，尤其是建筑工地附近的城市。目前,从国内实际工程建设的地质灾害危险性评估报告和有关文献来看,对于那些缺乏实际地面沉降观测资料的场地区域,地面沉降量的估算方法和模型较为庞杂,尚无一套简明有效的通用手段。究其原因,一方面是由于地下水位下降所引起的地面沉降机制较为复杂,不同地区、不同环境的具体情况差异较大,难以用几种简约通用的方法一概而论;另一方面,尽管有一些理论适用范围广、假设和模型均较为严密的方法,但在实际应用中,往往难以有效地获取其所需的各种资料和众多参数,仍然无法达到预期的效果。因此,关于地下水水位下降所引起的地面沉降预测,一些可行的作法是针对具体的区域特点和场地类型,建立最适合本区域的简约模型。 本文以北京东南部亦庄轻轨工程建设场地为例,探讨如何根据区域内数量有限的水位下降和地面沉降观测资料,确定具有代表性的区域骨架成分弹性储水因子和非弹性储水因子,进而采用有效应力法对该场地2006-2015年的地面沉降进行评估预测。 2 线路场地现状条件 拟建的北京亦庄轻轨是北京市城市轨道交通线网规划中连接北京市中心城和亦庄新城的轻轨线路,跨越丰台、朝阳、大兴和通县4个辖区和亦庄经济开发区,全长约23.8 km,其中高架线为13.90km,地下线为8.75 km,U型槽为0.65 km,地面线为0.50 km。共设车站14座,其中地下车站6座,高架车站8座。另外,起点站宋家庄设停车场一处,占地约5×104m2;终点站亦庄设车辆段一处,占地约21×104m2(图1)。 从工程地质上看,拟建场地位居永定河冲洪积平原中下部,沿线穿越古漯水河故道及古河道之间的台地。场地地形随地貌单元的变化略有起伏,自然地面标高范围在23.0~40.0 m之间,整体趋势为自西北向东南逐渐降低。覆盖层主要为第四纪松散冲洪积物,厚度在70~150 m左右,由黏性土、砂类土和卵石层交互沉积,具有较明显的沉积韵律。第四系孔隙水主要为台地潜水、层间潜水和承圧水,第四系下伏基岩主要为页岩、泥灰岩和砂岩等。 北京亦庄轻轨工程属重要建设项目,资料分析及野外勘察表明,场地区域存在或潜在的主要地质灾害包括活动断裂、砂土液化、地面沉降等。其中,由于地下水超采导致的地面沉降现象在北京市已非常突出,而拟建的线路场地又地处北京八里庄-大郊亭沉降中心南侧,个别区段的历史累积沉降量甚至高达200 mm。因此,地面沉降的定量的评估和合理的预测,成为本工程场地地质灾害评价的最重要内容之一。 3 含水层砂性及岩性砂性土复合运动学模型的建立 本文将首先利用本工程场地周围1999-2005年期间已知的水位下降量和对应的地面沉降量,反演获取有效应力法地面沉降模型中适合于本区域的两个关键参数:骨架成分弹性储水因子Sske和非弹性储水因子Sskv。在此基础上,利用该模型对北京亦庄轻轨工程建设场地2006-2015年的地面沉降进行预测评估。 有效应力法假设土体变形仅发生在垂向,而水平向的应变可以忽略不计。基于这种假设和有效应力原理,孔隙水压力有任何变化,有效应力都将跟着等值反方向变化,则有效应力变化为 式中:Δσ′为有效应力的变化值;ρw为水的密度;g为重力常数;Δh为含水层水位(头)的变化。 假如垂向有效应力的增加导致垂向压缩量呈线性增长,并考虑到对于砂性土,多为弹性变形;而对于黏性土,在有效应力大于其前期固结应力条件下,为非弹性变形,则在理论上可用下列公式计算地面沉降量: 式中:Δbe为含水层弹性压缩量(m);Δbv为含水层非弹性压缩量(m);Δh为水头变化值(m);eb为含水层砂性土层的累计厚度(m);vb为含水层黏性土层的累计厚度(m)。 于是,总的地面沉降量ΔB为 通常情况下,通过建设场地的地层概化模型,可分别获得场地不同位置砂性土层和黏性土层的累计厚度eb和vb,若可获得适合本区域的砂性土骨架成分弹性储水因子Sske和黏性土骨架成分非弹性储水因子Sskv,则可根据水头变化值Δh计算出相应地面位置的总沉降量ΔB。为有效获得适合本区域的Sske和Sskv,将基于建设场地所在区域的已有钻孔资料和实际地面沉降观测结果,采用最小二乘拟合方法,具体方法如下。 假设建设场所在区域已有n个钻孔,其中第i个钻孔的砂性土厚度和黏性土厚度分别为bei和bvi,在水头下降量为Δhi的情况下,实际观测的地面总沉降量为ΔBi,则可建立以下联立方程: 式中:iv为第i个钻孔处利用拟合反演结果Sske和Sskv计算获得的地表总沉降ΔB′与实际观测值ΔB之间的残差。所谓最小二乘拟合反演,就是确定出这样一套Sske和Sskv,使上述方程组中所有残差的加权平方和最小,即 式中:ip为第i个钻