The solubility of 2-ethylanthraquinone in binary mixtures of oligooxymethylene dimethyl ethers with different number of CH2O groups of n = 2, 3, and 4 (PODE3 + PODE4, PODE2 + PODE4, and PODE2 + PODE3) were measured by a synthetic method with a laser monitoring observation technique at the temperature range from 293.15 to 343.15 K at 101.325 kPa. The solubility data were correlated by the modified Apelblat model, lambda h equation, Wilson model, and the nonrandom two liquid (NRTL) model. The correlation results demonstrated that the NRTL model was the proper model with higher accuracy. Therefore, NRTL model was used to represent the temperature dependence of the molar fraction solubility. Furthermore, the partial molar excess Gibbs free energy, the partial molar excess enthalpy, and the partial molar excess entropy for 2-ethylanthraquinone dissolution were calculated using the activity coefficients (NRTL model).