随机介质固热耦合数学模型与岩石热破裂数值实验.pdf

里主三堡垫查盔堂堡主堂焦丝壅—— 摘 要 众所周知，在高温环境下，岩石物理力学特性是一个极重要的问题。在石油 的三次开采中，采用火烧油层的方法，降低油的粘度，同时诱发岩石破裂，达到 提高储层渗透性的目的；高放射性核废料的深层处置，在相当长的时间，美国、 法国～直致力于花岗岩岩层储存核废料，担心的就是岩石热破裂后，使核废料储 存区域与地下水相通，最终危及人类生存环境，同时核贮库的围岩是吸收放射性 核素并阻碍其迁移的最关键部分，其在高温作用下的长期力学性质及力学行为的 演变规律对贮存库的选址、设计及长期安全性预测具有决定性的影响，是需要研 究的主要问题之一；高温岩体地热资源开发中，从300。C以上的高温岩体中，提 取地热，由于温度降低，而产生热破裂；液体的地下储藏以及对地壳演变过程的 研究和模拟中，所有这些，都追切需要考虑工程围岩在高温下的物理力学性质和 力学行为的长期演变规律。所以，研究在高温条件下岩石的破裂以及物理力学性 质对高放射性核废料的深层处置、高温岩体地热资源开发等深部开采问题有着重 大的实际工程意义。 目前岩石破裂数值实验研究仅限于在外力载荷下对岩石进行变形、破裂和物 理力学性质的研究，在高温下对岩石的物理力学性质的研究又主要依赖于现场观 测和实验室物理实验。现场观测对工程而言是非常必要的，但由于这种方法受到 现场条件、人力、物力和财力的限制，在较大工程中难以充分发挥其作用。到目 前为止，还没有对岩石的热破裂及其在温度作用下岩石的物理力学性质的数值实 验研究。由于岩石材料的非均匀性、非连续性、以及几何结构的复杂性，现有的 解析方法尚缺少有效的手段对此过程进行研究，很难对岩石的热破裂与物理力学 性质的变化规律做相对准确的描述。因此，既考虑岩石材料的随机非均匀性，又 考虑了固体变形场与温度场的耦合作用和数值实验节省经费、直观演示等特点， 建立随机介质固热耦合数学模型以及有限元数值解法，采用数值实验方法来研究 岩石热破裂问题具有更重要的理论意义和现实意义。为此，本文从岩石介质的非 均匀性特征入手，在引入岩石非均匀性细观统计描述和物理学愈渗理论的基础 上，主要的研究内容如下： 1．考虑岩石颗粒的物理力学性质的随机非均匀性，从理论上建立了随机非均 质热弹塑性力学模型和随机介质固热耦合数学模型，并推导随机介质模型的有限 元数值解法。这些方程与方法，奠定了岩石热破裂研究的基础。 2．对岩石随机介质的热弹性力学模型的平面轴对称问题和球对称问题分别 进行了解析分析，提出了岩石物理力学参数作为随机变量在解析分析中的处理方 法，并给出在随机指数分布和韦泊分布下的解析解；对随机介质固热耦合数学模 摘要 型的平面轴对称问题也求出了解析解，包括位移、变形和热应力的精确解。 3．在二维随机非均质热弹塑性力学模型和四种随机分布(均匀分布、指数分 布、正态分布和Weibull分布)下，在高温条件下，对热膨胀系数变化引起的岩 石破裂门槛值温度的变化做了详细研究，揭示了岩石热破裂规律及门槛值温度随 分布参数i／1的变化规律，得到了一致的结论：岩石热破裂门槛值温度随分布参数 m的增加而升高，呈幂函数或直线关系；对岩石的力学性质进行了数值实验研究： 包括门槛值温度随着弹性模量和泊松比的变化规律，采用正态分布和Weibull分 布数值实验的结果与相应的实验结果相吻合；并再现了岩石热破裂的变化过程， 表明文中所采用的数值实验的方法是可行和可靠的。 4．在平面随机介质固热耦合数学模型和三种随机分布(指数分布、正态分布 和Weibull分布)下，对岩石的热物理特性(包括温度分布随热传导系数和比热 容的变化规律)、岩石的渗透率和热破裂细观机理进行了数值实验的研究。通过 数值实验详细讨论和分析了由于温度和岩石的非均匀性的作用下岩石试件中的 裂纹的产生、形成、扩展、汇合以及贯通过程的变化规律，从随机介质固热耦合 方面揭示了岩石热破裂的本质。这为深部开采提供了可靠的理论依据，同时对岩 土工程也具有～定的指导意义。 5．采用随机介质固热耦台数学模型，在正态分布、指数分布和Weibull分布 下，对高温岩体地热开发人工储留层二次破裂进行了数值模拟，揭示了岩体温度 场和应力场的特征和变化规律，同时直观地显示出人工储留层二次破裂随分布参 数111和时间变化的全过程，并表明热膨胀系数的概率分布形式对岩体热破裂规律 有重要影响