The research studies on the adsorption of surfactants on graphene help us to know how to use surfactants to exfoliate graphene from graphite or functionalize the graphene surface. Among them, molecular dynamics (MD) simulation has been widely used to investigate the adsorption of organic molecules and surfactants on graphene. In particular, coarse-grained (CG) MD simulation greatly improves the computational efficiency by simplifying the complexity of the studied systems, allowing us to explore the structure and dynamics of complex systems on larger spatial scales and longer time scales. However, an accurate prediction of the adsorption of surfactants on graphene is required by optimizing the interaction between surfactants and graphene, which is often overlooked by some CG models. In this work, we found that an accurate prediction of the adsorption enthalpies of organic molecules on graphene can be achieved by optimizing the interactions between organic molecules and benzene. Meanwhile, we simulated the adsorption of a surfactant on single-layer and double-layer graphene nanosheets, respectively. Our results revealed that increasing the temperature would favor the interactions between hydrophilic groups of surfactants. In addition, we discovered that the surfactant prefers to be adsorbed on the inner surfaces of double-layer graphene compared with the outer surfaces, and this is owing to the dehydration in the middle of double-layer graphene, which is beneficial to the hydrophilic interactions between surfactant molecules inside the double-layer graphene.