GPSFlow/Hydrate: A New Numerical Simulator for Modeling Subsurface Multicomponent and Multiphase Flow Behavior of Hydrate-Bearing Geologic Systems

Abstract

Numerical simulation has played a crucial role in modeling the behavior of natural gas hydrate (NGH). However, the existing numerical simulators worldwide have exhibited limitations in functionality, convergence, and computational efficiency. In this study, we present a novel numerical simulator, GPSFlow/Hydrate, for modeling the behavior of hydrate-bearing geologic systems and for addressing the limitations in the existing simulators. It is capable of simulating multiphase and multicomponent flow in hydrate-bearing subsurface reservoirs under ambient conditions. The simulator incorporates multiple mass components, various phases, as well as heat transfer, and sand is treated as an independent non-Newtonian flow and modeled as a Bingham fluid. The CH₄ or binary/ternary gas hydrate dissociation or formation, phase changes, and corresponding thermal effects are fully accounted for, as well as various hydrate formation and dissociation mechanisms, such as depressurization, thermal stimulation, and sand flow behavior. In terms of computation, the simulator utilizes a domain decomposition technology to achieve hybrid parallel computing through the use of distributed memory and shared memory. The verification of the GPSFlow/Hydrate simulator are evaluated through two 1D simulation cases, a sand flow simulation case, and five 3D gas production cases. A comparison of the 1D cases with various numerical simulators demonstrated the reliability of GPSFlow/Hydrate, while its application in modeling the sand flow further highlighted its capability to address the challenges of gas hydrate exploitation and its potential for broader practical use. Several successful 3D gas hydrate reservoir simulation cases, based on parameters from the Shenhu region of the South China Sea, revealed the correlation of initial hydrate saturation and reservoir condition with hydrate decomposition and gas production performance. Furthermore, multithread parallel computing achieved a 2–4-fold increase in efficiency over single-thread approaches, ensuring accurate solutions for complex physical processes and large-scale grids. Overall, the development of GPSFlow/Hydrate constitutes a significant scientific contribution to understanding gas hydrate formation and decomposition mechanisms, as well as to advancing multicomponent flow migration modeling and gas hydrate resource development.