PurposeNatural gas hydrate (NGH) has been regarded as one of the most important resources due to NGH's large amounts of reserve. However, NGH development still faces many technical challenges, such as low production rate and reservoir instability resulting from NGH decomposition. Therefore, developing a fully coupled THMC model for simulating the hydrate decomposition and studying its mechanical behavior is very important and necessary. The purpose of this article is to develop and solve a multi-phase, strong nonlinearity and large-scale fully coupled thermal-hydro-mechanical–chemical (THMC) model for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH decomposition.Design/methodology/approachIn this paper, a multi-phase, strong nonlinearity and large-scale fully coupled THMC model is developed for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH dissociation. The fully coupled THMC model is solved by using a fully implicit finite element method, in which the gas pressure, water pressure, temperature and displacement are taken as basic unknown variables. The proposed model is validated against with the experimental data, showing high accuracy and reliability.FindingsA multi-phase, strong nonlinearity and large-scale fully coupled THMC model is developed for simulating the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and rock deformation during NGH decomposition. The proposed model is validated against with the experimental data, showing high accuracy and reliability.Research limitations/implicationsSome assumptions are made to make the model tractable, including (1) the composition gas of hydrate is pure methane; (2) the gas-liquid multi-phase flow in the pore obeys Darcy's law; (3) hydrate occurs on the surface of soil particles, both of them form the composite consolidation material; (4) the small-strain assumption is applied to composite solid materials, which are treated as skeletons and cannot be moved; (5) momentum change caused by phase change is not considered.Practical implications NGH has been regarded as one of the most important resources due to its large amounts of reserve. However, NGH development still faces many technical challenges, such as low production rate and reservoir instability resulting from NGH decomposition. Most of the existing studies decouple the process with solid deformation and seepage behavior, but the accuracy of the numerical results will be sacrificed to certain extent. Therefore, it is very important and necessary to develop a fully coupled THMC model for simulating the hydrate decomposition and studying its mechanical behavior.Social implications NGH, widely distributed in shallow seabed or permanent frozen region, has the characteristics of high energy density and high combustion efficiency (Yan et al., 2020). A total of around 7.5 × 1,018 m3 has been proved to exist around the world and 1 m3 of NGH can release about 160–180 m3 of natural gas (Kvenvolden and Lorenson) under normal conditions. Safely and sustainably extracting NGH commercially can effectively relieve global energy pressure and contribute to achieving carbon reduction goals.Originality/valueThe novelty of the present work lies in mainly two aspects. First, a fully coupled THMC model is developed for studying the multi-physics processes involving solid-liquid-gas flow, heat transfer, NGH phase change and solid deformation during NGH dissociation. Second, the numerical solution is obtained by using a fully implicit finite element method (FEM) and is validated against experimental data.