On the Theory of Methane Hydrate Decomposition in a One-Dimensional Model in Porous Sediments: Numerical Study

The purpose of this paper is to present a one-dimensional model that simulates the thermo-physical processes for methane hydrate decomposition in porous media. The mathematical model consists of equations for the conservation of energy, gas, and liquid as well as the thermodynamic equilibrium equation for temperature and pressure (P−T) in the hydrate stability region. The developed model is solved numerically by using the implicit finite difference technique on the grid system, which correctly describes the appearance of phase, latency, and boundary conditions. The Newton–Raphson method was employed to solve a system of nonlinear algebraic equations after defining and preparing the Jacobean matrix. Additionally, the proposed model describes the decomposition of methane hydrate by thermal catalysis of the components that make up the medium through multiple phases in porous media. In addition, the effect of thermodynamic processes during the hydrate decomposition on the pore saturation rate with hydrates a7nd water during different time periods was studied in a one-dimensional model. Finally, in a one-dimensional model over various time intervals, t=1, 10, 50 s, the pressure and temperature distributions during the decomposition of methane hydrates are introduced and investigated. The obtained results include more accurate solutions and are consistent with previous models based on the analysis of simulations and system stability.