As one of the key status of gas in shale reservoir, adsorption gas accounts for considerable percentage of total gas amount. Due to the complexity and nanostructure of shale gas reservoir, it is very challenging to represent adsorption gas through traditional methods. However, the integration of the fractal theory and molecular dynamics (MD) simulation may provide a new perspective of understanding such nanostructure and the micro-phenomenon happening in it. The key purpose of this paper is to investigate the adsorption phenomenon in shale kerogen. By using MD simulation and grand canonical Monte Carlo (GCMC) algorithm, the adsorption of methane in 2, 5 and 10nm slit-like pores is simulated for different temperature and pressure status. According to the results, the average gas density in smaller pores is higher than that in bigger pores, and multilayer adsorption presents on some areas of pore surfaces. Then, the simulation results are analyzed using the multilayer fractal adsorption model. The analysis indicates that the number of adsorption layer increases with pressure increase: four-layer adsorption presents in 10nm pores while three-layer adsorption shows up in 2nm and 5nm pores due to pore volume limit. Fractal dimension of pore wall surface generated in this study is in the range of 2.31–2.63. Moreover, high temperature could decrease the adsorption behavior in reservoir condition.
