In wind turbine wake models, Gaussian models depend on multidimensional integration to ascertain the distribution of wake velocity deficits. These integrations, which often involve complex boundary conditions, significantly enhance the complexity of mathematical computations. Due to the difficulty of obtaining analytical solutions, numerical integration methods such as Monte Carlo or other numerical integration techniques are commonly employed. This study presents a three-dimensional wake model (3DJW) for horizontal axis wind turbines, utilizing the Weibull function to simplify wake deficit characterization instead of traditional Gaussian distribution methods. The 3DJW model considers wind shear effects and mass conservation laws to enhance predictions of vertical wake velocities. By integrating incoming wind conditions and turbine parameters, the model efficiently computes downstream wake velocities, improving computational efficiency. To enhance predictions in the ultra-far wake region, an improved three-dimensional Weibull wake model is proposed using the exponential fitting method. Validation through wind tunnel experiments and wind farm data demonstrates the model's accuracy in predicting wake deficits at the hub height, with relative errors in horizontal and vertical profiles mostly within 5% and 3%, respectively. The proposed model enables accurate and rapid calculation of wake velocities at any spatial location downstream, facilitating enhanced energy utilization and reduced costs.