Three-phase vapor-liquid-aqueous (VLA) equilibria can frequently appear in hydrocarbon reservoirs. The presence of the aqueous phase, which can commonly exist as a pure water phase, can dramatically increase the probability of encountering convergence issues in the three-phase equilibrium calculation algorithms. The primary reason is that there is a lack of good initial estimates of equilibrium ratios in both stability tests and flash calculations. In this research, we improve the three-phase VLA equilibrium calculation algorithm for water/hydrocarbon mixtures by developing a new scheme for initializing equilibrium ratios as well as a new procedure for three-phase equilibrium calculations. By applying this new initialization scheme, one can securely identify the aqueous-like stationary point on the tangent plane distance function. In our algorithm, convergence problems, which can frequently appear in the flash calculations of three-phase VLA equilibrium calculations, are avoided. We demonstrate the robustness and efficiency of the newly developed algorithm by showing a number of example calculations.