Molecular Simulation of Adsorption in Deep Marine Shale Gas Reservoirs

Abstract

Deep marine shale gas reservoirs are extremely rich in the Sichuan basin in China. However, due to the in situ conditions with high temperature and high pressure (HTHP), in particular reservoir pressure being usually much higher than the test pressure, it is difficult to accurately clarify the adsorption behavior, as seepage theory plays an important role in shale gas reserves evaluation. Therefore, three kinds of sorbent, including illite, quartz and kerogen, and two simulation methods, containing the grand canonical ensemble Monte Carlo method and molecular dynamics method, are synthetically used to determine the methane adsorption behavior under HTHP. The results show that both absolute adsorption and excess adsorption decrease with the increase of temperature. When the pressure increases, the absolute adsorption increases quickly and then slowly, and the excess adsorption first increases and then decreases. The superposition of wall potential energy is strongest in a circular hole, second in a square hole, and weakest in a narrow slit. The effect of pore size increases with the decrease of the pore diameter. Under HTHP, multi-layer adsorption can occur in shale, but the timing and number of layers are related to the sorbent type.