基于格林函数法的裂隙岩体渗流 - 传热耦合计算与参数敏感度研究.docxVIP

基于格林函数法的裂隙岩体渗流-传热耦合计算与参数敏感度研究

一、引言

1.1研究背景与意义

随着全球能源需求的增长以及对环境保护的日益重视，高效开发和利用新能源成为当务之急。在这一背景下，裂隙岩体渗流-传热耦合问题在多个关键领域展现出极其重要的研究价值。

在核废料处置领域，安全处置核废料是核能可持续发展的关键环节。国际上普遍采用地下处置库的方式，利用深部岩石洞室永久储存核废料。然而，核废料中含有大量毒性大、发热量大的核素，其衰变产生的热量会导致岩体温度升高，进而影响岩体中裂隙水的流动形态和运动规律。由于岩体节理和裂隙的渗透率远高于完整岩石，裂隙水的运动可能导致核素迁移，危害人类生存环境。因此，深入研究核废物处置库裂隙岩体内渗流-传热耦合作用，对于准确分析评价处置库的阻滞行为，保障核废料处置的安全性至关重要。

在地热开发领域，地热能作为一种清洁、可再生能源，开发利用地热能对于减少温室气体排放和实现能源结构转型具有重要意义。裂隙岩体是地热储层的重要组成部分，其渗流-传热特性直接影响着地热能的提取效率和开发成本。通过研究裂隙岩体渗流-传热耦合机制，能够更好地理解地热储层的物理过程，为优化地热开发方案、提高地热资源利用率提供理论支持，有助于推动地热能的大规模商业化开发。

格林函数法作为一种强大的数学工具，在解决复杂的偏微分方程问题中具有独特优势。在裂隙岩体渗流-传热耦合研究中，运用格林函数法可以将复杂的物理问题转化为数学上可求解的形式，通过建立半解析法计算流程，能够准确地描述裂隙岩体渗流-传热耦合的真实情况，便于考虑如边界热源等与实际工程相关的复杂因素，为工程设计和分析提供精确的理论依据。

参数敏感度分析则是研究不同参数对裂隙岩体渗流-传热耦合计算结果影响程度的重要方法。通过参数敏感度分析，可以确定影响计算结果的关键参数，明确各参数的影响规律和敏感程度。这不仅有助于深入理解裂隙岩体渗流-传热耦合的物理机制，还能为实际工程中的参数测量和控制提供指导，在工程设计阶段能够根据参数的敏感度合理选择参数范围，优化工程方案，降低工程风险和成本，提高工程的可靠性和稳定性。

1.2国内外研究现状

在裂隙岩体渗流-传热耦合理论方面，国内外学者开展了大量研究。国外学者Parts早在1966年就用解析方法探讨了在多孔介质中自由对流的规律，并讨论了水平流对自由热对流的影响，为后续研究奠定了理论基础。Harlan于1973年提出的土体冻结过程中的热质迁移数学模型，首次将渗流场与温度场耦合起来分析，推动了耦合理论的发展。Hodgkinson在1980年建立了核废料围岩水热对流的数学模型，采用级数方法探讨其运动规律，进一步深化了对特定工程背景下耦合问题的认识。国内方面，宋晓晨、徐卫亚对裂隙岩体进行等效化处理，基于混合物理论推导出裂隙岩体等效连续介质温度场、渗流场和应力场三场耦合的全耦合数学模型及其控制方程，丰富了国内在该领域的理论体系。张树光基于传热学和渗流理论建立了深埋巷道围岩的热扩散数学模型，并通过数值求解分析了风流和渗流耦合作用下围岩的温度场分布特征。

在格林函数法应用方面，国外研究起步较早，已经将其应用于多种复杂物理问题的求解。在裂隙岩体渗流-传热耦合研究中，也逐渐开始采用格林函数法来建立半解析模型，以提高计算精度和解决复杂边界条件问题。国内学者也在积极探索格林函数法在该领域的应用，通过改进和优化算法，使其更适用于我国实际工程中的裂隙岩体渗流-传热问题，如针对饱和裂隙岩体中高放核废料地下处置库裂隙岩体渗流-传热耦合问题，采用格林函数法建立计算流程，为工程规划建设提供服务。

关于参数敏感度分析，国内外学者针对裂隙岩体渗流-传热耦合问题，对不同参数进行了敏感度研究。例如，路威、项彦勇以高放废物处置库设计概念模型为背景，研究发现裂隙宽度、裂隙水流流速是影响裂隙水温度的高敏感度因素，而岩石的热导率、比热和密度则为低敏感度因素。通过这些研究，明确了各参数对渗流-传热过程的影响程度，为工程实践中的参数控制和优化提供了依据。

1.3研究内容与方法

本文的研究内容主要涵盖以下几个方面：首先，深入研究裂隙岩体渗流-传热耦合的基本理论，明确渗流和传热过程中的物理机制以及它们之间的相互作用关系。其次，基于核废料处置库或地热开发等实际工程的概念模型，建立简化的裂隙岩体渗流-传热耦合模型，合理考虑岩体的几何特征、物理参数以及边界条件等因素。然后，运用格林函数法建立半解析法计算流程，编写相应的计算机程序，实现对裂隙岩体渗流-传热耦合问题的数值求解。最后，通过具体算例，计算和分析单裂隙岩体渗流-传热耦合的基本特征，开展参数敏感度分析，确定影响渗流-传热计算结果的关键参数及其