将Taron等提出的颗粒聚集体的压力溶解模型引入笔者所研制的孔隙介质热-水-应力耦合有限元程序中，并使用摩尔-库仑准则，针对一个假设的实验室尺度且位于饱和石英颗粒聚集岩体中的高放废物地质处置库模型，拟定弹性分析和弹塑性分析两种计算工况，进行4 a处置时段的数值模拟，考察了岩体中的温度、颗粒界面水膜及孔隙中的溶质浓度、迁移和沉淀质量、孔隙率及渗透系数、孔隙水压力、地下水流速和应力及塑性区的变化、分布情况。结果主要显示：弹塑性分析中由于应力调整和增大了分子扩散系数，使得塑性区的颗粒介质的溶解、迁移和沉淀有明显的变化，并对渗流场（孔隙水压力及流速）和应力场产生显著的影响。但两种工况弹性区中的颗粒介质的溶解、迁移和沉淀差别较小。

This paper introduces the model of pressure solution for granular aggregates established by Taron et al. into the finite element code for analysis of thermo-hydro-mechanical (T-H-M) coupling in porous media developed by the first author. Using the Mohr-Coulomb yield criterion, a hypothetical disposal model for nuclear waste located in a saturated quartz aggregate rock mass with a laboratory scale is simulated. Two computation cases, the elastic analysis and the elastoplastic analysis are designed. Then the corresponding numerical simulation for a disposal period of 4 years is carried out. The states of temperatures, solute concentrations in the intergranular fluid film and at the pore space, migration and precipitation masses, porosities and permeabilities, pore pressures, flow velocities and stresses and plastic zones in the rock mass are investigated. It is shown that, because of the stress adjustment and the increased molecular diffusivity in the elastoplastic analysis, there are obvious changes of solution, migration and precipitation of aggregate medium in the plastic zones, and the seepage field (including pressures and flow velocities of pore water) and the stress field are markedly influenced. But the differences of solution, migration and precipitation of aggregate medium in the elastic zones for two cases are not significant.