DEM simulations of methane hydrate exploitation by thermal recovery and depressurization methods

• Published in: Computers and geotechnics Vol. 80; pp. 410 - 426
• Main Authors: Jiang, Mingjing, Fu, Chang, Cui, Liang, Shen, Zhifu, Zhu, Fangyuan
• Format: Journal Article
• Language: English
• Published: New York Elsevier Ltd 01-12-2016; Elsevier BV
• Subjects: Axial strain; Axial stress; Bonding; Cementing; Compression; Compressive strength; Computer simulation; Confining; Contraction; Depressurization; Discrete element method; Dissociation; Distinct element method; Distribution; Exploitation; Hydrates; Hydrostatic pressure; Ice; Ice formation; Ice loads; Macro- and micro mechanical responses; Mechanical properties; Methane; Methane hydrate; Pore pressure; Pore water; Pore water pressure; Pressure reduction; Recovery; Sand; Sediment load; Sediments; Simulation; Soils; Studies; Temperature; Temperature effects; Thermal recovery; Two dimensional models; Void ratio; Water pressure; Methane hydrate; Macro- and micro mechanical responses; Thermal recovery; Depressurization; Distinct element method
• ISSN: 0266-352X; 1873-7633
• DOI: 10.1016/j.compgeo.2016.05.011

Methane hydrate (MH, also called fiery ice) exists in forms of pore filling, cementing and load-bearing skeleton in the methane hydrate bearing sediment (MHBS) and affects its mechanical behavior greatly. To study the changes of macro-scale and micro-scale mechanical behaviors of MHBS during exploitation by thermal recovery and depressurization methods, a novel 2D thermo-hydro-mechanical bonded contact model was proposed and implemented into a platform of distinct element method (DEM), PFC2D. MHBS samples were first biaxially compressed to different deviator stress levels to model different in-situ stress conditions. With the deviator stress maintained at constant, the temperature was then raised to simulate the thermal recovery process or the pore water pressure (i.e. confining pressure for MH bond) was decreased to simulate the depressurization process. DEM simulation results showed that: during exploitation, the axial strain increased with the increase of temperature (in the thermal recovery method) or decrease of pore water pressure (in the depressurization method); sample collapsed during MH dissociation if the deviator stress applied was larger than the compression strength of a pure host sand sample; sample experienced volume contraction but its void ratio was slightly larger than the pure host sand sample at the same axial strain throughout the test. By comparison with the laboratory test results, the new model was validated to be capable of reproducing the exploitation process by thermal recovery and depressurization methods. In addition, some micro-scale parameters, such as contact distribution, bond distribution, and averaged pure rotation rate, were also analyzed to investigate their relationships with the macroscopic responses.