The energy diffusion process in a few two-dimensional Fermi-Pasta-Ulam-type lattices is numerically simulated via the equilibrium local energy spatiotemporal correlation. Just as the nonlinear fluctuating hydrodynamic theory suggested, the diffusion propagator consists of a bell-shaped central heat mode and a sound mode extending with a constant speed. The profiles of the heat and sound modes satisfy the scaling properties from a random-walk-with-velocity-fluctuation process very well. An effective phonon approach is proposed, which expects the frequencies of renormalized phonons as well as the sound speed with quite good accuracy. Since many existing analytical and numerical studies indicate that heat conduction in such two-dimensional momentum-conserving lattices is divergent and the thermal conductivity κ increases logarithmically with lattice length, it is expected that the mean-square displacement of energy diffusion grows as tlnt. Discrepancies, however, are noticeably observed.