2D FEM analysis for coupled thermo-hydro-mechanical-migratory processes in near field of hypothetical nuclear waste repository

Abstract

In order to consider the influence of temperature and stress fields on the migration of radioactive nuclide with underground water movement, an elastoplastic model and a 2D FEM code for analysis of coupled thermo-hydro-mechanical (THM) processes in saturated and unsaturated porous media were extended and improved through introducing the percolation and migration equation, so that the code can be used for solving the temperature field, flow field, stress field and nuclide concentration field simultaneously. The states of temperatures, pore pressures and nuclide concentrations in the near field of a hypothetical nuclear waste repository were investigated. The influence of the half life of the radioactive nuclide on the temporal change of nuclide concentration was analyzed considering the thermo-hydro-mechanical-migratory coupling. The results show that, at the boundary of the vitrified waste, the concentration of radioactive nuclide with a half life of 10 a falls after a period of rising, with the maximum value of 0.182 mol/m³ and the minimum value of 0.181 mol/m³ at the end of computation. For a half life of 1 000 a, the concentration of radioactive nuclide always increases with the increase of the time during the computation period; and the maximum value is 1.686 mol/m³ at the end of the computation. Therefore, under the condition of THM coupling, the concentration of radioactive nuclide with a shorter half life will decrease more quickly with water flow; but for the radioactive nuclide with a longer half life, its concentration will keep at a higher level for a longer time in the migration process.

Keywords

radioactive nuclide / concentration / thermo-hydro-mechanical-migratory coupling / 2D FEM analysis

References

[1] ChijimatsuM., NguyenT. S., JingL., JongeJ., KohlmeierM., MillardA., RejebA., RutqvistJ., SouleyM., SugitaY.. Numerical study of the THM effects on the near-field safty of a hopothetical nuclear waste repository—BMT1 of the DECOVALEX III Project. Part 1: Conceptualization and characterization of the problems and summary of results [J]. International Journal of Rock Mechanics and Mining Sciences, 2005, 42(5/6): 720-730

[2] MillardA., RejebA., ChijimatsuM., JingL., JongeJ., KohlmeierM., NguyenT. S., RutqvistJ., SouleyM., SugitaY.. Numerical study of the THM effects on the near-field safty of a hopothetical nuclear waste repository—BMT1 of the DECOVALEX III Project. Part 2: Effects of THM coupling in continuous and homogeneous rocks [J]. International Journal of Rock Mechanics and Mining Sciences, 2005, 42(5/6): 731-744

[3] RutqvistJ., ChijimatsuM., JingL., MillardA., NguyenT. S., RejebA., SugitaY., TsangC. F.. Numerical study of the THM effects on the near-field safety of a hopothetical nuclear waste repository—BMT1 of the DECOVALEX III Project. Part 3: Effects of THM coupling in sparsely fractured rocks [J]. International Journal of Rock Mechanics and Mining Sciences, 2005, 42(5/6): 745-755

[4] LiangB., LiuL., XueQ., SunW.-ji.. Study on numerical simulation of nuclides leakage for groundwater pollution [J]. Journal of System Simulation, 2007, 19(2): 261-311

[5] ZhaoY.-l., CaoP., ZhaoY.-s., Linh., WangY.-xian.. Dual media model for thermo-hydro-mechanical coupling and 3D numerical simulation [J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(S2): 4024-4031

[6] LiuY.-r., YangQ., HuangY.-s., LiX.-qin.. Parallel numerical analysis of coupled fluid flow and stress based on dual porosity media model [J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(4): 705-711

[7] FengZ.-c., ZhaoY.-s., LüZ.-xing.. Study on percolation law of 2D porous and fractured double-medium [J]. Acta Physica Sinica, 2007, 56(5): 2796-2801

[8] OlivellaS., GensA.. Double structure THM analyses of a heating test in a fractured tuff incorporating intrinsic permeability variations [J]. International Journal of Rock Mechanics and Mining Sciences, 2005, 42(5/6): 667-679

[9] RutqvistJ., BarrD., DattaR., GensA., MillardA., OlivellaS., TsangC. F., TsangY.. Coupled thermalhydrological-mechanical analyses of the Yucca Mountain Drift Scale Test—Comparison of field measurements to predictions of four different numerical models [J]. International Journal of Rock Mechanics and Mining Sciences, 2005, 42(5/6): 680-697

[10] SonnenthalE., ItoA., SpycherN., YuiM., AppsJ., SugitaY., ConradM., KawakamiS.. Approaches to modeling coupled thermal, hydrological, and chemical processes in the drift scale heater test at Yucca Mountain [J]. International Journal of Rock Mechanics and Mining Sciences, 2005, 42(5/6): 698-719

[11] LiuJ., ShengJ., PolakA., ElsworthD., YasuharaH., GraderA.. A fully-coupled hydrological-mechanical-chemical model for fracture sealing and preferential opening [J]. International Journal of Rock Mechanics and Mining Sciences, 2006, 43(1): 23-36

[12] RutqvistJ., FreifeldB., MinK. B., ElsworthD., TsangY.. Analysis of thermally induced changes in fractured rock permeability during 8 years of heating and cooling at the Yucca Mountain Drift Scale Test [J]. International Journal of Rock Mechanics and Mining Sciences, 2008, 45(8): 1373-1389

[13] ZhuW. C., BruhnsO. T.. Simulating excavation damaged zone around a circular opening under hydromechanical conditions [J]. International Journal of Rock Mechanics and Mining Sciences, 2008, 45(5): 815-830

[14] IvarsD. M.. Water inflow into excavations in fractured rock—A three-dimensional hydro-mechanical numerical study [J]. International Journal of Rock Mechanics and Mining Sciences, 2006, 43(5): 705-725

[15] MinK. B., RutqvistJ., TsangC. F., JingL.. Thermally induced mechanical and permeability changes around a nuclear waste reposity — A far-field study based on equivalent properties determined by a discrete approach [J]. International Journal of Rock Mechanics and Mining Sciences, 2005, 42(5/6): 765-780

[16] YuanS. C., HarrisonJ. P.. Development of a hydro-mechanical local degradation approach and its application to modelling fluid flow during progressive fracturing of heterogeneous rocks [J]. International Journal of Rock Mechanics and Mining Sciences, 2005, 42(7/8): 961-984

[17] LeiwsR. W., SchreflerB. A.The finite element method in the deformation and consolidation of porous media [M], 1987, New York, John Wiley

[18] CHIJIMATSU M, FUJITA T, KURIKAMI H, ITO A, SUGITA Y. Subtask BMT1-B: Implication of THM coupling on the near-field of a nuclear waste repository in a homogeneous rock mass [R]. DECOVALES III-Task3-Bench Mark Test 1(BMT1)-2002: 1–43.

[19] ZhangY.-jun.. Coupled thermo-hydro-mechanical model and 2D FEM analysis for FEBEX in situ test [J]. Chinese Journal of Geotechnical Engineering, 2007, 29(3): 313-318

[20] ZhangY.-jun.. Coupled thermo-hydro-mechanical model and relative numerical analysis for near field of conceptual nuclear waste repository [J]. Rock and Soil Mechanics, 2007, 28(1): 17-22

[21] KurikamiH.Study on modeling for coupled thermo-hydromechanical phenomena in rock mass [D], 2003, Kyoto, Kyoto University

[22] NISHIGAKI M. Density dependent transport analysis saturated-unsaturated porous media—3. Dimensional Eulerian Lagrangian method [R]. Okayama University, 2001.