A novel method for globally approximating the static voltage stability region boundaries (SVSRBs) of power systems is proposed by applying parametric polynomial approximation to the criterion equation which defines the SVSRB. Known as the Galerkin method, the implicit function portraying the SVSRB is described as a linear combination of basis polynomial functions, and the coefficients of the basis functions are obtained by projecting the SVSRB criterion onto each basis function. The approximation guarantees high precision globally in the whole domain of the parameter of interest rather than only in the neighborhood of a point, and the error can be controlled by the degree of polynomial basis functions. In the meantime, analytical expression of the left or right eigenvector of the system's Jacobian matrix corresponding to the zero eigenvalue is obtained in the form of polynomial, which provides valuable information for online voltage stability control or monitoring. Case studies in a 10-bus test system and IEEE 118-bus test system verifies the validity, accuracy, and flexibility of the proposed method.