地铁车站井点降水引起地面沉降计算

井点降水引起地面沉降计算分析摘要：井点降水所引起的地下水位的变化，使得土体中的有效应力发生变化，也即是附加应力的变化是基坑周围土体发生沉降的根本原因，本文以稳定渗流基本理论为基础，同时结合分层综合法提出了一种计算降水引起周围土体沉降方法，并且相应的给出了其计算步骤。关键词基坑降水稳定渗流沉降1前言随着社会的发展，科技越来越复杂，我们人类的建筑也越来越复杂，各种各样的建筑朝着地下发展，其中就涉及地下岩土体问题，所以研究地下岩土体的结构与性质和地下水的变化用着重要的意义，它关系着我们所修建的建筑物的安全，其中地下水中的井点降水对地面的沉降又是重要的研究方向，近年来，由于基坑降水引起的事故屡见不鲜，已经引起了科研、设计、施工人员的高度关注，对于降水引起的地面沉降的研究，一些学者已经做了很多的研究。其中，按照积水建筑物的延伸方向，可分为垂直集水建筑物和水平集水建筑物。前者包括井和竖井等，统称为井，地下水想垂直集水建筑物的运动简称为“井流”，本文以渗流理论作为基础，研究井点降水中地下水位的变化，接着利用分层总和法计算地面沉降的一些问题。井点降水引起的土层中水位的下降，其中涉及到水在土层中的流动情况，包括渗透理论等等(1)确定降水以后水位线的问题；(2)计算由于水位线变化引起含水层本身的变形问题。本文将依次对这两个问题进行讨论。基坑降水对地面沉降的影响的实质就是需要解决两个主要问题：2计算理论2.1渗流的基本理论水受固体边界的约束，只能在空隙中流动。由于固体边界的几何形状十分复杂，使得空隙中地下水的运动要素（例如流速矢量）的分布变化无常，若从这个微观水平上研究地下水的运动规律，实际上是不可能的，也是没有必要的。人们研究地下水的运动规律，必须从宏观水平上来考察，为此设计一个假想的流场。这个流场首先不能将水流约束在空隙之中，否则不仅涉及复杂的固体边界表面的刻画，而且水流在空间上是不连续的，使得一切基于连续函数的微积分手段都不能利用。因此，我们必须引入一个假象的水流代替真实的地下水流，这种假想水流是：充满整个多孔介质的连续体；而且这种假想水流的阻力与实际水流在空隙中所受的阻力相同；它的任意一点水头H和水流速度v等要素与实际水流在该点周围一个小范围内的平均值相等。这种假想水流便是宏观水平的地下水流，我们称之为“渗流”，它所占据的空间称之为“渗流场“。在这种假象的基础上，法国水力工程师亨利.达西在装有均质沙质的圆柱形桶中做了大量大实验，最后他的出来了一个理论，叫做达西定律，也即为渗流基本定律，其形式为：Q=KA=若将达西定律用于二维的或者是三维的地下水运动，则水力坡度不是常量，沿流向可以变大也可以变小，它应该用微分方式来表示，即在直角坐标系中也可表示为2.2裘布衣井流2.2.1裘布依稳定潜水井流在下列假定的条件下建立的：均质、各项同性、隔水地板水平的圆柱形潜水含水层，外侧面保持定水头，中心一口完整抽水井在此称之为圆岛模型，没有垂向入渗补给和蒸发排泄，且渗流服从线性定律的稳定流动。在上述条件下的潜水井中进行定流量抽水，经过一段时间后，渗流将会趋向稳定；潜水面由原来的水平状态变成漏斗状，在此称它为水位降落漏斗。依渗流连续性原理，这时各渗流断面的流量都相同，并等于抽水井的流量。裘布依稳定井流裘布依稳定潜水井流从平面上看：流线沿径向指想井轴，等水位线是同心圆，这种流动是径向流动。由于靠近孔处得水力坡度大，远离孔底的水力坡度小，所以等水位线在井孔附近密集，往外变疏。从剖面上看：最底部一根线是水平的直线，最上（潜水面处）的一根流线（也也称浸润曲线）是曲率最大的凸行曲线；中间的流线则过度。由上至下，从曲率最大的曲线逐渐变为水平的直线。剖面上的等水头线是也是一系列的弯曲程度不等的曲线，外围的等势线趋向铅锤的直线。从空间讲，等水头面围绕井轴旋转的一系列曲面，这些曲面的方程预先是难以想到的，为使问题简化起见，引入裘布衣假定，把剖面上的等水头线视为铅垂线，即忽略流速的垂向分量，从而把三维井流问题简化为二维流动来解决，这时的渗流断面被视为圆柱面。对流网有了基本的认识之后，可以建立有关的方程。我们以隔水底板为基准面，因此潜水面处的水头值等于渗流厚度h。从上面分析，这种情况下地下水流属于轴对称问题，采用极坐标更为方便，在此，取井轴为h轴，向上为正；沿隔水底板取r轴，向外为正。根据达西定律和裘布衣假定，任意渗流断面的流量由于h随r的增大而增大，因而，>0,如上分析，将渗流断面视为圆柱面，所以A=则求积分，将r由至R，由至。依无入渗补给、蒸发排泄以及稳定流的条件，各断面间的流量相等，则得上式是裘布衣稳定潜水井流的用水量方程，若引进水中水位降深，则上式可写成假如积分上、下限改为：r由至r；h由至h，则可得到降落漏斗曲线方程，即该式表明，降落漏斗曲线取决于内外边界的水位和，与Q和K无关.2.2.2裘布依稳定承压井流此部分与裘布依潜水井流基本相同，只是含水层是等厚的承压含水层。在这种条件下抽水剖面上的流线是相互平行的直线等水头线是铅垂线，等水头面则是真正的圆柱面，在这种情况下，渗流的水力坡度是相同的，其流动方向沿着r轴向。根据流网分析，渗流断面为圆柱面，即代入达西定律，得分离变量积分取积分极限：r为至R，H由至H0，得这就是裘布依稳定完整承压井流的涌水量方程若积分上下限改为：r由至r；H由至H，则可得到降落漏斗曲线方程，即2.3无限含水层中单个定流量井流含水层是有限的，或者说含水层都有它的外边界；如果含水层（水平方向）是如此之大，以致外边界对于含水层研究区段的水头分布没有明显的影响，则可称它为无限含水层；由此可见，有限含水层和无限含水层是相对的。对于压力传导系数小的含水层进行短时间的抽水的情况，可视为无限的含水层；但对于压力传导系数大的含水层做长时间的抽水情况，则可能要看成有限的含水层。有些抽水试验可视为无限的含水层，但对于长期开采地下水，计算时则可能要考虑位有限含水层。无限含水层中单个定流量井流方程基于下列假定条件：eq\o\ac(○,1)含水层是均质的、各项同性、等厚且水平分布，水和含水层均假定为弹性体；eq\o\ac(○,2)无垂向补给、排泄，即W=0；eq\o\ac(○,3)渗流满足达西定律；eq\o\ac(○,4)完整井，假定流量沿井壁均匀进水；eq\o\ac(○,5)水头下降引起地下水从储量中的释放是瞬时完成的；eq\o\ac(○,6)抽水前水头面是水平的；eq\o\ac(○,7)井径无限小且定流量抽水；eq\o\ac(○,8)含水层侧向无限延伸。根据条件eq\o\ac(○,1)和eq\o\ac(○,5)，可以应用轴对称流基本微分方程或，eq\o\ac(○,6)是初始条件；eq\o\ac(○,7)和eq\o\ac(○,8)分别是内、外边界条件。因此，该定解问题可写为，（0≤r<，t>0）（0≤r<）（t>0）（常量）（t>0）该定解可用积分变换法、分离变量法博尔兹门变换法求解，它的解为其中地下水动力学上习惯上记称为泰斯井流的井函数。于是可得计算承压井流的三条基本方程为和对于潜水完整井流问题，显然比承压井流要复杂的多。这里要注意一个问题，即承压井流的厚度M是不变的，而潜水井流的厚度h是变化的。如果忽视其他方面的差异，那么只要寻求潜水井流的h与承压井流的M之间的对应关系（水位h与水头H以及重力给水度的对应关系是显然的）之后，便可以将承压井流的解用于潜水井流。如果潜水井流满足前述承压井流的8个假定条件，其中第eq\o\ac(○,1)条改为“含水层是均质、各项同性、等厚且含水层底板水平”，再加第eq\o\ac(○,9)个条件，降深值远远小于潜水含水层厚度，流动满足于裘布衣假定，则潜水井流与承压井流可以对应起来。满足上述9个假定条件的潜水完整流，其流动微分方程为，即【=2\*ROMANII】（0≤r<,t>0）其中初始条件和边界条件为t>0（常量）t>0我们通过变换将承压井流问题【=1\*ROMANI】改变为潜水井流的问题【=2\*ROMANII】的形式，就可以获得参数间的对应关系。为此，在承压流的基本微分方程式两端同乘以M，并引入承压流势函数的定义则微分方程变为【=3\*ROMANIII】边界条件分别为t>0（常量）t>0由此可见，承压井流引入势函数的定义之后，承压完整井的定解问题【=3\*ROMANIII】与潜水完整井的定解问题的形式完全一样，于是其解的形式也应相同，只是要注意两者变量间的关系：潜水井流承压井流含水层厚度：给水度：势函数：由承压完整井流的解改写为其中利用上述三个变换关系，直接可以求得潜水完整井流定流量抽水的方程若用水位降深表示，则依则有公式对比承压井流方程可见，潜水井流的平均厚度值可按下式近似计算换句话说，我们可以引用抽水时潜水层的平均厚度的定义将潜水井流方程与承压井流方程对应起来或相互转化。在欧美一些文献上往往采用雅克布修正降深来转化上述两类井流方程关系，即记称为修正降深，为此，含水层的厚度保持不变而降深做出相应的修改。为此对应定流量抽水的承压完整井流三个基本方程，定流量的抽水潜水完整井流的三个基本微分方程为和上面所述的承压完整井流与潜水完整井流的三个基本方程均称为泰斯公式.2.4关于影响半径一些文章从，认为s=0处的r就是影响半径R，即或从而得出实际上，我们并不需要去计算影响半径R,因为它并不是含水层的参数。假如一定要去讨论降深为零处的径距r（r=R）,那么泰斯公式的回答就是：在开始抽水的那一瞬间，含水层处的水头处处都有下降；含水层的范围就是影响范围。当然，这已被理想化了，如果真实一些，这影响半径R为式中：为声音在含水层种的传播速度。3算例3.1工程概况XX站西南风道道位于XXX街下，其其出口在XXX对面，XXX总公司和和XX大楼之之间夹道中中风道分南南北遭和斜斜风道，斜斜风道走向向为偏东65。，在南北北风道南端端施工竖井井。双层结结构，开挖挖最大断面面13.11m12..6m。风道由由西南向东东北3%。下坡，南南端开挖标标高22.220m，挖深20.990m；北端开开挖标高21.770m，挖深21.440m。本本工程处于于XX冲积扇扇脊部的中中下部，围围岩分类为为j-1类。根据据XX站西半半部勘察剖剖面资料，揭揭示地层情情况如下：：土层名称埋深厚度人工填土3.5m轻亚粘土13-16m园砾16.5-177.5m亚粘土19-21m中砂22-24m厚厚度2-3mm卵石29-30m轻亚粘土3.2计算降水水曲线及确确定沉降地地层根据文献降水计计算结果，水水位变化涉涉及到三个个地层即==1\\*GBB3①园砾层、==2\\*GBB3②亚粘土层层、=33\*GB3③中砂层。3.3降水水一年期的的地面沉降降计算（1）园砾层对于园园砾层属砂砂层，可不不考虑时间间考虑因素素，取其层层厚H=2..5m。根据降降水计算结结果，水位位降深均在在园砾层底底以下，故故水位变化化值均取为为=0.556m，水的容容重=10，压缩模模量=455MPa。代入式，可可得此层降降深值=0.3311mmm（2）亚亚粘土层考虑时时间因素，其其层厚为H=3m，同样水水位变化值值均取为==3m，水的容容重=100KN/mm3压缩模量=777.7MPPa。渗透系系数k=1..73m//d，计算固固结度时按按三角形附附加应力分分布双面排排水，取H=1..5m，代入可可得此层最最终沉降量量为=111.6888mm。由式可可得任意时时间的沉降降量。（3）中中砂层不考虑时时间因素，其其层厚H=3m，水位变变化值随位位置变化为为一变量，水水的容重==10KNN/m3，压缩模量=500MPa。代入，即即可得中砂砂层沉降值值3.4结果汇汇总由降水引起的地地面沉降计计算值，可可得任一点点任意时刻刻的由降水水引起的沉沉降值。其其中仅列出出沿风道中中心线各点点由降水引引起的地面面沉降计算算结果如表表1所示。表1沿沿风道中心心线各点由由降水引起起的地面沉沉降计算结结果表45.6603X/mY/m初沉降/mm终沉降/mm001.99513.6890102.17