Molecular dynamic simulation was used to investigate kerogen extraction.
Extraction mechanisms of kerogen moieties with supercritical CO2 are elucidated.
Electrostatic interaction played a significant role in kerogen–shale interaction.
Effects of pressure and temperature on kerogen extraction were interpreted.
The simulation results are consistent with experimental results.