基于温度场-渗流场-应力场耦合的双重介质流固热耦合分析.docxVIP

基于温度场-渗流场-应力场耦合的双重介质流固热耦合分析 1 裂隙岩体耦合 对岩体（尤其是断裂岩体）的温度场耦合研究具有重要的理论意义，并具有广泛的实践背景。高温干热岩体的开发、核废料贮存库的安全评估、环境工程中污染源的迁移及深部开采等都迫切需要研究岩体的温度场–渗流场–应力场三场耦合作用。从20世纪80年代以来,对该问题的研究取得了重要进展。J.Noorisha等提出了饱和裂隙岩体的THM耦合基本方程式,N.Barton和S.C.Bandis针对工程岩体对三场耦合进行初步探讨;L.Jing等给出了相对系统的三场耦合模型,但对模型简化和实用方面的研究还不够深入。近十几年来,提出了许多耦合模型,如单纯的热传导模型,水流动模型、水热耦合模型、固流热耦合模型;考虑岩体介质性态,有均质模型、裂隙介质模型等。相比之下较全面的模型是T.Kohl等提出的固流热耦合数学模型,T.A.Sprecker等提出并完善了裂隙岩体的温度场、渗流场与应力场耦合的有限元解法。当前,随着能源储库、核废料贮存库及干热岩体地热开发等大型岩土工程的出现,三场耦合研究成果已逐步服务于工程,用于评价工程的长期稳定性、安全性[9~12]。国内学者李宁等建立了一些耦合模型,但大多是基于多孔连续介质或离散介质假设。本文采用双重介质(拟连续介质+裂隙介质)提出研究裂隙岩体温度场–渗流场–应力场耦合数学模型,开发了三场耦合数值模拟软件,通过与渗流场–应力场耦合,应力场–温度场耦合进行数值比较可知,进行双重介质温度场–渗流场–应力场耦合分析是十分必要的。 2 裂隙岩体膜渗流规律模型 裂隙岩体是岩土工程中广泛遇到的一类复杂岩体,是由随机分布的裂隙和被裂隙切割的岩块组成的不连续介质。裂隙岩体温度场–渗流场–应力场耦合作用是一个相对复杂的研究问题,主要表现在裂隙岩体赋存地质环境的各个组成部分,即渗流场、应力场与温度场都随时间、空间而发生变化;与此同时,以上各个组成部分之间的耦合作用处于一种复杂动态变化过程中。传统的离散介质假设(忽略岩块孔隙系统的透水性)和连续介质假设(将裂隙中的渗流平均到岩体的渗流中去)往往很难真实地反映裂隙岩体的温度场—渗流场–应力场耦合作用,而双重介质(拟连续介质+裂隙介质)三场耦合模型在一定程度上能真实地模拟裂隙岩体渗流、传热过程,对大型裂隙(断层、节理、人工裂缝等)可以查明,按裂隙介质对待;对被这些大中型裂隙切割而成的含众多低序次的小裂隙介质则不可能查明,认为是不确定的,采用拟连续介质模型模拟。然后根据两类介质接触处水头、温度、位移相等来建立裂隙岩体的三场耦合模型。 本文建立的三场耦合模型基于如下假设: (1)岩体是由拟连续介质、裂隙介质组成的结构体,拟连续介质为均质同性的弹塑性体,裂缝介质服从节理单元模型。 (2)渗流规律服从非线性达西定律。 (3)不考虑水的汽化,岩体被单相水饱和。 (4)热质在固相、液相介质中传递方式以传导和强迫对流为主,各介质的比热容及热传导系数不随温度和压力而变化。 (5)介质的变形以孔隙变形为主。 2.1 裂隙介质拟连续介质及kfi 式中:e为体积变形;βT,βP分别为流体热膨胀系数和压缩系数。 对于拟连续介质: 对于裂隙介质: 式(3)~(7)中,kij,kfi分别拟连续介质和裂隙介质渗透系数;Θe为有效体积应力;b为裂隙宽度;v(T)为运动黏滞系数;φ为孔隙率;qi为比流量;p为水压;si为切向坐标。 2.2 耦合的能量平衡方程 根据热力学第一定律:外界传入系统的能量与内部热源产生的能量之和等于物质内能的增量与对外做功之和,建立分别考虑岩体(固相)与温度场、液相与温度场耦合的能量平衡方程,该方程考虑了液体热对流(强迫对流为主)和固相热传导而导致的能量迁移。 对于拟连续介质,有 裂缝中主要介质为水,能量的迁移是通过裂隙水流的热传导、热对流(强迫对流为主)及岩体裂隙表面和裂隙水流之间的热传导的双过程实现的,即 式(8),(9)中:λr,λw分别为岩体和水的热传导系数;ρr,ρw分别为岩体和水的密度;Cpr,Cpw分别为岩体和水的热容系数;Tr,Tw分别为岩体和水的温度;Trb为岩体裂缝边缘的温度。 2.3 温度场–渗流场—静力平衡控制方程 对于拟连续介质,静力平衡控制方程为 对于裂隙介质,静力平衡控制方程为 式(10)~(12)中:λ,μ均为拉梅常数;uj,ji,ui,jj均为岩体位移;F为体力;β为岩体热膨胀系数;σn′,σs′分别为裂缝法向和切向有效应力;kn,ks分别为裂缝法向与切向刚度;εn,εs分别为裂缝法向与切向变形;α为等效孔隙压系数。 式(10)~(12)这3组控制方程,再辅以初始和边界条件,就构成了基于双重介质假设的裂隙岩体温度场–渗流场–应力场耦合数学模型。 这一模型具有如下特点:流体渗流方程中考虑了岩体的变形和