针对高放射核废深地质处置库近场环境，建立分布热源作用下饱和裂隙岩体三维水流-传热过程中位移和应力的一种半解析计算方法：采用Goodier热弹性位移势和Laplace变换计算由温度梯度产生的温梯位移和应力；考虑单一裂隙的情况，利用经典弹性力学的Boussinesq解和Cerruti解计算为满足边界条件的约束位移和应力，与温梯位移和应力叠加，可得总体热位移和应力；把裂隙面离散为矩形单元集合，采用极坐标系下的解析法计算包含奇点的单元积分，采用数值法计算与分布热源有关和不含奇点的单元积分。与基于裂隙面法向一维热传导假设的一种解析解对比，结果表明，半解析法与解析法的计算结果基本一致，但由于半解析法考虑岩石的三维热传导，因温度时空分布和演变的不同而导致不同的温梯应力。针对一个假想单裂隙岩体三维水流-传热过程，计算温梯位移和应力、约束位移和应力、总体位移和应力；结果表明，裂隙水流-传热可能对位移和应力的分布和演变有显著影响，距离分布热源较近的岩石因升温膨胀受到约束而出现压应力，而距离分布热源较远的岩石则可能因协调收缩受到约束而出现拉应力。

With regard to the near-field environment of deep geological repositories of high-level radioactive wastes, this study formulated a semi-analytical method for calculating the displacement and stress in the 3D processes of water flow and heat transfer in saturated sparsely fractured rocks. Goodier’s thermo-elastic displacement potential and Laplace transform are employed to calculate the temperature-gradient induced displacement and stress. Considering the case of a rock mass containing one single fracture, Boussinesq’s solution and Cerruti’s solution in the classic theory of elasticity are utilized to obtain the constraint-induced displacement and stress for complying with the boundary conditions, which supplement the temperature-gradient induced displacement and stress to arrive at the total thermal displacement and stress. The fracture is discretized into a set of rectangular elements, on which the integrals involving singularities are evaluated using an analytical method with polar coordinates, whereas numerical integration is employed for non-singular integrals and for the integral concerning the distributed heat source. Comparison with an analytical solution, which assumes one-dimensional heat conduction normal to the fracture walls, indicates that the semi-analytical results are essentially agreeable with minor differences only in the region where the effect of 3D heat conduction tends to be more significant. For a hypothetical 3D water flow and heat transfer process in a rock containing one single fracture, the distributions of temperature gradient induced, the constraint displacement induced, and the total displacements and stresses are calculated and examined, which reveal, among other things, that water flow and heat transfer may exert significant influences on the spatiotemporal variations of the displacements and stresses; and that the compressive stresses may occur in the region that is close to the distributed heat source due to constrained thermal expansion, while tensile stresses may exist in the region that is relatively further away from the distributed heat source due to constrained compatible contraction.