大坝及其周围地质体中渗流-温度-应力耦合的研究综述.docVIP

大坝及其周围地质体中渗流-温度-应力耦合的研究综述 周慧敏 梅明荣 (工程力学系 河海大学 南京 210098) 摘要:分析了大坝周围地质岩体中渗流场、温度场与应力场的相互作用，对三场耦合重要参数及控制方程进行了系统的分类与总结,并对工程使用进行归纳展望。 关键词：大坝地质岩体、渗流场、温度场、应力场、耦合分析 1 引 言 据统计，30 % ~40%的水电工程大坝失事是由地下水流渗透压力、温度应力引起岩体变形，而岩体变形又使其渗透性能发生变化，即渗流场、温度场和应力场（THM）相互作用是溃坝的主要原因。考虑单一力的作用解决问题效果往往并不理想，各种作用力相互交叉和渗透十分突出。因此，在设计时为了保证大坝长期安全运行有必要进行大坝-岩体系统THM 全耦合分析研究。近年来，THM研究方兴未艾，在核废料储存库防污染问题上已经扩展到THMC。经过大量阅读和深入思考，作者对工程中THM耦合分析的研究成果进行系统的评述与小结。 2 渗流场、温度场和应力场之间的相互作用 渗流场－应力场耦合：坝体-岩体系统中渗流场的变化引起渗流作用力包括渗流体积力和渗透压的改变,其中的应力状态改变使渗透系数改变进而改变渗流场。有两种可能：（1）贯通裂隙之中水压力的传递过程和扩容劈裂作用,诱使裂纹不断地扩展,对岩石强度和破裂模式的影响十分显著，裂隙的扩展导致孔隙水压力和梯度分布的非均匀性,重新改变应力分布形式,相互贯通的裂隙作为新的渗流通道，将使渗透系数增大；（2）裂隙岩体受法向、侧向应力作用时，岩体变得密实，导致裂隙岩体渗透系数降低。 -岩体系统中渗流场改变，将使渗透水流参与坝体-岩体系统中的热量传递与交换，从而影响温度场的分布。另一方面，坝体-岩体系统中温度场的改变，形成的温度梯度引起水的运动改变坝体-岩体系统渗透系数。而岩石的弹性模量及强度都是温度的函数,因此温度的改变也将导致坝体-岩体系统应力场的变化。 , 将在坝体-,比热、导温系数则随温度的升高而增加，裂隙岩体温度升高，裂隙封闭，开通裂隙数密度减少，导水性降低，温度场对应力场的影响还体现在温度的变化将导致孔隙流体压力较大的变化；当应力场发生变化时,存在岩体内部力学能量转换等。 (包括热传导和热对流)时，在合适的地质环境下渗流场、温度场与应力场之间会产生耦合作用。Darcy定律： （1） 其中：V为渗流场中任一点的渗流速度，K为渗透系数，J为水力梯度。 假定岩体中热传导符合Fourier定律： （2） 其中T—温度，J—热流密度，K—为导热率。 能量守恒方程：根据热力学第一定律，单位体积中单位时间内，由外界传入系统的能量与内部热源产生的能量之和等于物质内能的增量与对外界做功之和。 3.1 渗流场-应力场耦合 双场耦合方程多采用机理分析法。此两场耦合方程已经较为成熟，控制方程格式基本如下[1]： （3） （1）渗流连续性方程，根据质量守恒定理，假设介质变形-应力为线弹性关系，将连续方程、状态方程、Darcy定律结合起来：水压梯度引起的流动和渗流源之和恒等于流体在介质内的累积率。所谓的流体在介质内的累积率分别由密度变化与孔隙度变化引起。孔隙变化有三个直接原因：体积应变；液体孔隙压力引起的颗粒大小改变；有效应力引起的颗粒大小改变。其中后两个影响相对较小可以忽略。 （2）弹性力学应力与应变协调方程，渐变性材料应力引起的应变用位移表示，包括弹、塑性应变所引起的应力。其建立的基础需要形成合适的能够描述受力表面变化的硬化规律及能够描述裂隙材料本构关系的流动规律。 裂隙岩体渗透系数是最重要的水力参数，裂隙岩体渗流渗透系数的非均匀性非常突出，导致其各向异性非常明显，岩体裂隙隙宽、单一裂隙体表比等影响因素难以确定。叶源新[2]提出渗透系数k是孔隙变化量的函数，参照己有的试验成果，设定渗透系数与孔隙率成指数关系： （4） （5） 式中：a为耦合参数，表征应力、应变对渗透系数的影响程度，由试验确定；p为孔隙水压力；为孔隙变化量；为体积应变；Q、H、为Biot常数。 3.2 渗流场-温度场耦合 其耦合方程也较为成熟。方程意义清晰。控制方程格式基本如下[3、4]： (6) （1）渗流连续性方程与渗流场-温度场耦合相似，不同之处在于荷载由温度等效。水头梯度与温度梯度引起的水流通量之和，再计算进水源，恒等于流体在介质内的累积率。 2）热量平衡方程：渗流夹带的热量与坝体-岩体本身热传导作用的热量之和为单位时间内流入单位体积净热量，再算进热源，恒等于单位时间内岩体温度升高所吸收的热