For the models of stress corrosion and pressure solution used in the 2D finite element method(FEM) code of thermo-hydro-mechanical coupling analysis for dual-porosity medium developed by the authors

a correction factor for water saturation is introduced. Aiming at a hypothetical model for geological disposal of nuclear waste located at a water-bearing stratum

three computation conditions with different initial fracture saturations were designed(Sw20=1.0、0.8、0.2)

then the corresponding two-dimensional numerical simulation for the coupled thermo-hydro-mechanical processes were carried out; and the states of temperatures

rates and magnitudes of aperture closure

pore and fracture pressures

flow velocities and stresses in the rock mass were investigated. The results show：With the change of initial fracture saturation from high value to low value

the closure rates caused by stress corrosion and pressure solution vary from fast to slow

and the times of fracture aperture tending to its residual value and contact-area ratio of asperity tending to its nominal area from their initial values also increase

as well as the drop in stress intensity factors slows down; the changes and distributions of fracture and pore pressures in near-field as well as the states of flow vectors have obvious distinctions; the magnitudes and distributions of stresses within the rock mass in three calculation cases make little difference.