We investigate the spectral energy distributions (SEDs) of Massive Young Stellar Objects (MYSOs) using the various infrared observational data including the Infrared Space Observatory (ISO) data. We model the dust envelopes around the stars using a radiative transfer model for spherically symmetric geometry. Comparing the model results with the observed SEDs of the two MYSOs (AFGL 4176 and AFGL 2591), we derive the relevant dust shell parameters including the dust opacity, the dust density distribution, and dust temperature distribution. We find that the spherical model can produce the SEDs roughly similar to the observations. We expect that the results would be helpful for making more realistic non-spherical dust envelope models for MYSOs.