多孔介质中地面沉降与地下水头下降的关系.docxVIP

多孔介质中地面沉降与地下水头下降的关系 1 储水层变形特性分析 在jacob假设下，对地面沉降机研究，尤其是对储水层垂直压缩量与地下水头沉降深度之间的关系，则储水层的弹性储水率是必要的参数。储水层变形后，如果参数发生变化，则对进一步的推挤有重要影响。但仅由该参数的传统或一般定义式无法回答这个问题。本文在Jacob假定下通过推导储水率与多孔弹性常参数的关系,将发现储水率的骨架弹性贡献部分为与时间无关的常数,这将使地面沉降机理得以用几个主要控制要素构成的简明关系式来表达,特别是导致地面沉降与地下水位之间具有线性关系。 2 储水率的定义 储水率是指孔隙水头上升(下降)一个单位时从单位体积的多孔介质中储存(释放)的体积水量。其定义式最早在半个世纪前由Jacob(1940,1950)在含水层只存在垂向弹性形变、上覆荷载保持不变及含水层骨架颗粒不可压缩的假定下通过推导以水头为因变量的渗流连续性方程而给出 其中ρw为流体密度,g为重力加速度,β为流体体积压缩系数,α为多孔介质的体积压缩系数,n为孔隙率。上述的假定后来被总称为Jacob假定,其核心是含水层只发生垂向弹性形变。β和α分别被定义作 及 其中p为孔隙流体压力,式(3)成立的前提是介质颗粒本身不可压缩。Irmay(1968)注意到定义式(1)在推导过程中由于假定控制单元形变而导致的错误,从而给出了它的正确定义式 该式表明含水层的弹性储水率由骨架弹性储水率(1-n)αρwg与孔隙水弹性储水率nβρwg两部分构成。在高压缩性的多孔介质体系中,与α相比,水的弹性储水率部分常常可以忽略;但在高渗透性单元中必须考虑水的压缩性,因为它与骨架的压缩性处于同样量级水平。 储水率作为含水层的重要特性参数之一,主要被水文地质学家和石油工程师们所使用,且包含在几乎所有地下水流模型所采用的水流方程中。水流方程的基本形式如下 其中q为Darcy比流速或孔隙水相对于多孔介质颗粒的相对比流速,·表示散度,H为水头,t为时间。 3 含水层弹性储水率的机制 设压应力取正,压缩应变取正,根据广义Hooke定律,假定初始时刻含水层系统不受应力作用,那么均质多孔介质中弹性应力应变的本构关系可以表达为(Verruijt,1969) 式中为有效主应力;σii为总主应力;p为流体压力;εii为应变分量;σ≡εxx+εyy+εzz为骨架体应变;G=E/2(1+v)为剪切模量,E其中为杨氏模量,v为干土泊松比;λ为拉梅弹性常数,其与E和v的关系为λ=Ev/(1+v)(1-2v)。应用Jacob的水平方向零应变(εxx+εyy=0)及垂直方向总荷载的零变化(σzz=σ=+p=C)假定,由方程(6)中的第三式得到流体压力与体应变(或垂向应变)的微分关系为 在忽略骨架颗粒的压缩性且假定骨架颗粒与水的均质性的前提下,冉兴龙(2003)将Helm(1987)给出的以孔隙流速和颗粒流速表达的三维固结多孔介质的质量守恒方程改写成如下的形式 值得特别指出的是,当忽略水的压缩性时,式(8)右端项为零,这时式(8)实际上就是三维固结时的Terzaghi-Rendulic固结方程(Rendulic,1935)。将式(7)、(8)结合得 孔隙水压力p可用水头H与空间位置高度z表达为 其中g为重力加速度,z为位置高度。将式(10)代入式(9)得 对比方程(11)与方程(5)得到 实际上,该式括弧中的第二部分可以看作一个新的多孔弹性常参数cm,即令 这样一来,式(12)简化为 这就是采用多孔弹性介质的常参数建立起来的储水率的定义式。 比较式(14)与式(4)得 由于α定义为多孔介质(包括骨架与孔隙)的垂向体积压缩系数或垂向压缩系数,所以cm是多孔介质骨架的垂向体积压缩系数或骨架垂向压缩系数,式(13)及式(15)均显示cm恒大于零。 在Jacob假定下,由于cm为与时间无关的常数,加之水的重率ρwg随温度压力的变化也很小,所以含水层弹性储水率的骨架弹性贡献部分(cmρwg)为与时间无关的常数;而弹性储水率的孔隙水弹性贡献部分(nβρwg)却是随着含水层形变而变化的。 4 地面沉降与地下水头沉降之间的线性关系 4.1 垂向变形与狱内地面沉降 忽略颗粒本身的压缩性时,不管是否考虑水的弹性形变,冉兴龙(2003)在Jacob假定条件下通过将储水层形变的质量连续性方程(Helm.1987)与孔隙水流连续性方程(式(5))的耦合推导出利用体应变ε(压缩应变取正)与孔隙水头H表出的储水层形变关系方程 式中为储水层的骨架弹性储水率,其定义由储水层弹性储水率SS的骨架弹性部分构成,即 注意到式(15)并代入,使式(17)变成 据上节,正是含水层的骨架弹性储水率,为与时间无关的常数。式(16)表明:任一点处单位体积储水层在单位时间内由于水头上升而由介质骨架膨胀弹性储存的体积水