饱和岩体热流固耦合数学模型的研究.pdfVIP

饱和岩体热流固耦合数学模型研究 1 2 刘善利 ，赵坚 1. 长江南京航道工程局，南京(210011) 2. 河海大学水利水电工程学院，南京 (210098) E-mail: liulanxun123245@ 摘 要：在高压注水油采作业、冻土冷融及核废料库围岩等工程中，都可能会由于地质环境 条件（如内部温度、压力等）产生比较剧烈改变而导致岩体介质密度发生改变。基于对此类 岩体热流固耦合问题的认识和研究，本文把液固介质的密度作为随温度、压力改变而变化的 参数，把粘性耗散能热源项和温度梯度引起的流量变化流源项反映在热流固耦合作用过程 中，推导建立了考虑介质密度变化的饱和岩体热流固耦合数学模型。 关键词：饱和岩体；耦合；热流固；数学模型 0. 引言 岩体热流固（Thermo-Hydro-Mechanical，简写为THM ）三场耦合研究已深入到地热资 源开发利用、核废料安全处置、非等温煤层瓦斯渗流、冻土及油气田开采等领域中。建立反 映温度场、渗流场和应力场各场之间相互作用的耦合数学模型是研究岩体热流固耦合问题的 基础，也是重点。国内外很多学者根据不同研究对象和目的，先后建立了形式各异的热流固 耦合数学模型，如非等温瓦斯渗流数学模型、三维裂隙网络块裂介质模型、连续介质参数模 型、饱和小变形线热弹性模型、污染物运移模型、经典流固耦合方程组加上温度场经验方程 建立的模型、基于含油气盆地的多场耦合模型、冻土模型、饱和热弹性多孔介质瞬态全耦合 模型、核废料地下处置多相耦合模型等。在此类模型中，一般把液固介质的密度项作为常数 来考虑，即不考虑或忽略介质密度随温度、压力变化影响因素，也忽略源汇项的作用。研究 和实验[1-2]表明，这样处理与实际情况有出入，特殊情况下，分析结果会有较大误差。因此， 研究在温差和压力变化比较剧烈情况下的介质密度项和源项改变影响是必要的。 本文根据高压注水油采作业、冻土冷融及核废料库围岩等工程中，由于地质环境条件（如 内部温度、压力等）产生比较剧烈改变，可能导致岩体介质密度发生改变的情况，以孔祥言 [3] [4] 提出的“饱和小变形线热弹性模型”模型 为基础，运用混合分析法 ，推导建立固液介质密 度随温度、压力变化和源项等影响因素的饱和岩体温度场、应力场和流体渗流场耦合的数学 模型，提出了模型应用的环境条件，并在模型中加入反映粘性耗散能的热源项和温度梯度引 起的流量变化流源项。 1. 基本思路 建立热流固三场耦合数学模型的基础是构造反映物质特性的物性方程。物性方程包括状 态方程和本构方程。状态方程是联系物质(如地下水、岩体骨架等)和各种热力学参数之间关 系的方程，也称平衡态热力学特性方程，它反映物质特性参数随压力和温度变化的关系；而 本构方程是描述依赖于特定物质内部结构固有反应的方程，它反映物质分子热运动特性或内 在的激发-反应关系，即非平衡态热力学特性方程。 1.1 状态方程 [3] 孔祥言 等基于饱和多孔材料小变形情形的线性热弹性理论建立了热流固耦合数学模 型，模型考虑了流体和固体密度以及孔隙度随压力和温度的变化关系和液体粘度随温度变化 - 1 - 等因素。本文在建立地下水、固体介质状态方程中，采用孔祥言的理论方法，考虑介质密度 随压力和温度的函数关系。即 1 ρ ： ○地下水密度 f ρ ρ exp[c (p −p ) −β (T −T )] ρ +ρ c Δp −ρ β ΔT (1) f f 0 f 0 f 0