In order to explore the effect of the attachment site of the benzene ring in the backbone of the surfactant on its diffusion characteristics on the surface of anthracite, the molecular dynamics simulation method was used, and the four isomers (m-C16, m = 2,4,6,8; m represents the attachment site of the benzene ring in the backbone) of sodium hexadecyl benzene sulfonate (SHS) were selected. Binary models of surfactant/anthracite, surfactant/graphene modified by oxygen-containing functional groups, and a ternary model of water/surfactant/anthracite were constructed. By analyzing a series of properties such as interaction energy, contact surface area, relative concentration distribution, radial distribution function, hydrophobic tail chain order parameter, etc., it is concluded that the adsorption strength of 4-C16 on the surface of anthracite is the highest; the reason is that 4-C16 has the highest degree of aggregation near the oxygen-containing functional groups on the surface of anthracite. Further investigations find that 4-C16 can be densely covered on the ketone group, and the longer branch chain of 4-C16 has the highest degree of order in the Z-axis direction.