基于PR方程以及Vander Waals混合规则，对多元有机混合工质的气液相平衡计算方法进行了研究，以二元混合有机工质R410A、三元混合有机工质R401A、四元混合有机工质R405A为示例工质，进行了等温泡点计算，考察了两种不同的计算相互作用系数方法对计算结果的影响，并将计算结果与美国标准局的NIST REFPROP 7.0的计算结果进行比较。结果表明，两种方法对于压力及气相组分的计算的平均相对误差均小于10%，满足工程计算需要，第一种计算相互作用系数的方法有着更高的精度，其计算压力的平均相对误差最大为3.14%；在温度较低以及气相摩尔组分较小时，拥有较大的误差，其平均误差满足精度要求。

Based on the PR equation and the Vander Waals mixing rule, the vapor-liquid equilibrium calculation method of multi-component organic mixed working fluid was studied. Taking binary mixed organic containing R410A, ternary mixed organic containing R401A and quaternary mixed organic containing R405A as example working fluid, the isothermal bubble point calculation was carried out, and the influence of two different methods of calculating the interaction coefficient on the calculation results was investigated. And compare the calculation results with the calculation results of NIST REFPROP 7.0 of the US Bureau of Standards. The results show that the average relative error of the two methods for the calculation of pressure and gas phase components is less than 10%, which meets the needs of Engineering calculation. The first method for calculating the interaction coefficient has higher precision, and the average relative error of the calculated pressure is 3.14%; When the temperature is low and the gas phase molar fraction is small, there is a large error, and the average error satisfies the accuracy requirement.

云南省教育厅科学研究(2019J0025)；国家自然科学基金(51736005)