Para estimar con precisión el Retraso Troposférico Zenith (ZTD) para el posicionamiento de alta precisión del Sistema Global de Navegación por Satélite (GNSS), este estudio propone un método de modelado de ZTD basado en el modelo UNB3m y la interpolación espacial Kriging, en el cual los parámetros óptimos de interpolación espacial se determinan en función de los errores correspondientes a diferentes combinaciones de los parámetros de interpolación, y la distribución espacial de las estaciones de modelado GNSS se determina por los errores de interpolación de las estaciones GNSS seleccionadas aleatoriamente en varias ocasiones. Para verificar la precisión y fiabilidad del modelo propuesto, se utilizan como entradas las estimaciones de ZTD de 132,685 épocas con una resolución temporal de 1 h o 2 h durante 28 años, desde 1997 hasta 2025, de la red global de estaciones de seguimiento GNSS en operación continua; los resultados de ZTD en cualquier posición y el momento de observación correspondiente se pueden obtener con el modelo propuesto. Los resultados experimentales muestran que el error del modelo es inferior a 30 mm en más del 85% de las épocas de observación, los resultados de estimación de ZTD se ven menos afectados por la posición horizontal y la altura de las estaciones GNSS que los modelos tradicionales, y el error de interpolación de ZTD se mejora en 10-40 mm en comparación con los modelos GPT3 y UNB3m en las cuatro estaciones de verificación GNSS. Por lo tanto, esta tecnología puede proporcionar resultados de estimación de ZTD para el monitoreo de deformaciones híbridas de frecuencia simple y dual, así como datos densos de ZTD para la inversión de Vapor de Agua Precipitable (PWV). Dado que el método propuesto tiene las ventajas de una implementación simple, alta precisión, alta fiabilidad y facilidad de promoción, se espera que se aplique completamente en otras aplicaciones de posicionamiento de alta precisión.