The experimental devices were established using a vapor-liquid phase equilibrium sampling method. The solubility of hydrogen in the model components of coal liquid, such as hexadecane, methylbenzene, 1,2,3,4-tetrahydronaphthalene, naphthalene, quinoline, and 1-naphthol, was measured. A mathematical model based on Pierotti method and Henry law is proposed. The expression of the model is ln x (H2) = -a/T + bT + c ln T + d ln P-H2 + e, in which related parameters in the model are obtained by regression of known solubility data for different solvents. Using the model, the literature data of hydrogen solubility were regressed, and related parameters and mathematical expressions were obtained. The average absolute deviation between regression values and experimental results is 1.4-5%.