Molecular dynamics simulation of the intercalation behaviors of methane hydrate in montmorillonite

The formation and mechanism of CH4 hydrate intercalated in montmorillonite are investigated by molecular dynamics (MD) simulation. The formation process of CH4 hydrate in montmorillonite with 1 ~ 8 H2O layers is observed. In the montmorillonite, the “surface H2O” constructs the network by hydrogen bonds with the surface Si-O ring of clay, forming the surface cage. The “interlayer H2O” constructs the network by hydrogen bonds, forming the interlayer cage. CH4 molecules and their surrounding H2O molecules form clathrate hydrates. The cation of montmorillonite has a steric effect on constructing the network and destroying the balance of hydrogen bonds between the H2O molecules, distorting the cage of hydrate in clay. Therefore, the cages are irregular, which is unlike the ideal CH4 clathrate hydrates cage. The pore size of montmorillonite is another impact factor to the hydrate formation. It is quite easier to form CH4 hydrate nucleation in montmorillonite with large pore size than in montmorillonite with small pore. The MD work provides the constructive information to the investigation of the reservoir formation for natural gas hydrate (NGH) in sediments.