渗透固结系数变化对地面沉降计算影响.pdfVIP

渗透固结系数的变化对地面沉降计算的影响 张炳峰，王国体 合肥工业大学，安徽合肥 (230009) 摘 要：通过对经典 Biot 渗流力学进一步的考察，讨论在土体饱和状态下，渗透系数随孔 隙率变化而变化对固结系数的影响，从而考虑固结系数变化对固结度，沉降计算的影响。在 已有的地面沉降计算方程中增加一个反映渗流系数和孔隙变化率关系的耦合方程。对阜阳市 的实测沉降资料进行相应的分析计算和对比，其结果表明这个方程能很好的反映地面沉降过 程。 关键词：渗流与应力耦合，地面沉降，固结系数 1. 引言 对地下水日益增加的开采，尤其是在松散冲积层、湖泊或浅海沉积地带，其后果之一 是导致地表下降或沉陷，即地面沉降。由于地下水开采而导致较大的地面沉降，在高度发达 地区是比较常见的。对地面沉降状况及补救办法的事例研究，对于将要面临同样问题的发展 中地区是有帮助的。地下水开采引起的地下水压力变化将导致渗流介质体及其相关地质体受 力状态的变化，进而导致渗流介质体变形(这方面最为典型的是地下水开采诱发地面沉降)； 渗流介质体的变形过程中，其物性参数(孔隙度n 、压缩系数ε、渗透系数K 等)不再是经典 理论中的常量，而是与力学状态有关的变量；如此，渗流问题的研究就必须既要考虑介质体 的应力应变，又要考虑渗流特性的变异；这类问题，显然不是经典渗流理论所能解决的，然 而Biot 三维固结理论能够较好的解决这一问题。 2. 地面沉降过程中渗流与应力耦合模型 2.1 经典Biot 渗流力学耦合方程分析 1941 年，美国物理学和应用数学家Maurice Anthony Biot ，将Terzaghi 有效应力理论推 广到了三维空间上的固结问题研究中，并在1942、1954 和1956 年，将该问题的研究进行了 [1] 深入和完善，由此建立的Biot 孔弹性理论 ，为流固耦合问题的研究奠定了理论基础。 Biot 孔弹性理论模型主要有四个基础方程组成：平衡方程(1 式) 、几何方程(2 式) 、本构 方程(3 式和4 式)和渗流方程(5 式) ； 在三维空间上，基本方程可写为以下模型。 ∂σ 平衡方程： ij +ρX j 0 （i,j=1,2,3 ） （1） ∂x ij 几何方程： ε (u =+u ) / 2 ，ε ε +ε +ε （2 ） ij i , j j ,i v 11 22 33 ′ =− + 本构方程： σ σ αp δε 2Gε （3 ） ij ij ij v ij ∆n p / Q =−αε p / H =−σ / 3H （4 ） v ii 渗流方程： 2 1 ∂p ∂εt （5 ） K ij ∇ p =− α Q ∂t ∂t 3λ+2G 2(1−µ)G E K