Traditional methods for estimating ocean heat flux often introduce large uncertainties due to complex parameterizations. To tackle this issue, we developed a novel framework based on maximum entropy production (MEP) theory. By incorporating heat storage effects and refining the Bowen ratio, we enhanced the MEP method's accuracy. This research derives a new long-term global ocean latent heat flux dataset that offers high accuracy, enhancing our understanding of ocean energy dynamics.