After lengthy diagenesis and tectonic movement, a rock mass inevitably develops many pores and micro-fissures. A numerical simulation method was employed to study the thermal response characteristics of the rock mass under temperature seepage coupling by treating it as a pore-fissure medium and considering its anisotropic properties. Based on the mixed finite volume method (FVM), a numerical scheme of the governing equation for the temperature seepage coupling of the pore-fissure medium is derived, with the program solution module independently written in C++. On this basis, a numerical test model of the fissured mudstone is established to analyze the distribution of the rock mass temperature field under various thermal conductivities and the influence of fissure permeability on the seepage field. The mixed FVM results revealed that the temperature and water pressure distributions near the fissure were closely related to the directionality of thermal conductivity in the rock mass, as well as the thermal conductivity and permeability coefficient, respectively, of the fissure itself. Comparison with results from the finite element software ABAQUS demonstrated significant advantages of the proposed method when solving temperature and seepage problems in discontinuous geological bodies containing hiatuses, mutations, and fissures.

经过漫长的成岩作用和构造运动，岩体不可避免地会发育出许多孔隙和微裂缝。采用数值模拟方法研究岩体在温度渗流耦合作用下的热响应特性，将岩体视为孔裂介质并考虑其各向异性特性。基于混合有限体积法(FVM)，推导了孔裂介质温度渗流耦合控制方程的数值格式，程序求解模块采用C++独立编写。在此基础上，建立裂隙泥岩数值试验模型，分析不同热导率下岩体温度场分布及裂隙渗透率对渗流场的影响。混合FVM结果表明，裂隙附近的温度和水压分布与岩体中热导率的方向性以及裂隙本身的热导率和渗透系数密切相关。与有限元软件 ABAQUS 的结果比较表明，在解决包含裂隙、突变和裂缝的不连续地质体的温度和渗流问题时，该方法具有显着优势。