The valve-side bushing of a converter transformer is a critical component in ultra-high-voltage direct current (UHVDC) systems, making its monitoring through digital twin technology highly significant. However, the complex structure and spatio-temporal nonlinearity of the bushing result in a large computational demand for its digital twin model, which requires an effective order reduction algorithm. This paper proposes a Spatio-temporal Non-uniformity Proper Orthogonal Decomposition (SN-POD) algorithm considering the inhomogeneity of space and time consumption to meet the reduced-order computational requirements of UHV valve-side bushings. This proposed method reduces the calculation time to 10% of the full-order simulation model while controlling the error range of the key research area less than 0.1%. The test results show that this method has good robustness, calculation speed, and accuracy. This research can significantly enhance the computational efficiency of digital twin modeling for valve-side bushings and provide a technical foundation for constructing digital twin models for UHV valve-side bushings.