En este documento, presentamos una solución aproximada de alto orden con precisión uniforme para ecuaciones integrales de Volterra no lineales en 3D. Este esquema numérico se construye basado en el método de interpolación Lagrangiana cúbica en bloque tridimensional. Al mismo tiempo, realizamos un análisis del error de truncamiento local del esquema numérico basado en el teorema de Taylor. A través del análisis teórico, llegamos a la conclusión de que el orden de convergencia óptimo de este esquema numérico de alto orden es 4. Finalmente, verificamos la efectividad y aplicabilidad del método a través de cuatro ejemplos numéricos.