Crystallization is widely used for the isolation and purification of active pharmaceutical ingredients (APIs) in the final manufacturing stage. Bazedoxifene (BZD) acetate, classified as a third-generation selective estrogen receptor modulator, is employed in the treatment of postmenopausal osteoporosis. This work reports the equilibrium solubility data of BZD acetate in fifteen different mono-solvents, including methanol, ethanol, n-propanol, i-propanol, n-butanol, ethyl formate, methyl acetate, ethyl acetate, acetone, 2-butanone, cyclohexanone, tetrahydrofuran (THF), 2-methyltetrahydrofurane (2-MeTHF), 1,4-dioxane and acetonitrile (ACN), within the temperature range “T = 283.15–323.15 K” under p = 0.1 MPa. The solubilities of BZD acetate in various solvents increase with increasing temperature. Among the alcoholic solvents, the solubility order is as follows: methanol > ethanol > n-propanol > n-butanol > i-propanol. In esters solvent, the solubility of BZD acetate decreases in the following order: methyl acetate > ethyl acetate > ethyl formate. For ketones solvent, the solubility decreases in the following order: cyclohexanone > acetone > 2-butanone. The maximum value of BZD acetate solubility was 0.011 mol·mol−1 in tetrahydrofuran at T = 323.15 K, but the least data was found in ethyl formate at 283.15 K. Four thermodynamic models, i.e., the Yaws model, modified Apelblat model, λh model, and Wilson model were applied to fit the solubility data of BZD acetate in different solvents. Among them, the Wilson model represents the most accurate fitting outcomes. Finally, the thermodynamic functions of mixing (including enthalpy, entropy, and Gibbs energy) were calculated, and the results indicate a mixing process that is both spontaneous and driven by entropy. The solubility data, correlated models, and derived thermodynamic functions present essential thermodynamic fundamentals for the industrial production of BZD acetate crystallization, facilitating efficient separation and purification.