Theoretical description of Doppler-broadened noise-immune cavity-enhanced optical heterodyne molecular spectroscopy under optically saturated conditions