In recent years, the incidence rate of pulmonary embolism (PE) has increased dramatically. Currently, the correct diagnosis rate of PE in China is relatively low, and the diagnosis error rate and missed diagnosis rate were as high as about 80%. The most standard method of PE detection is pulmonary artery digital subtraction angiography (DSA), but pulmonary artery DSA is an invasive examination, and patients can have certain risks and discomfort. Noninvasive monitoring of PE remains challenging in cardiovascular medicine. We attempt to study the light propagation in human thoracic tissues and explore the possibility of near-infrared spectroscopy (NIRS) in noninvasive detection of PE. In this study, by utilizing the Monte Carlo simulation method for voxelized media and the Visible Chinese Human dataset, we quantified and visualized the photon migration in human thoracic region. The influence of the development (three levels) of PE on the light migration was observed. Results showed that around 4.6% light fluence was absorbed by the pulmonary tissue. The maximum signal sensitivity distribution reached 0.073% at the 2.8- to 3.1-cm light source-detector separation. The normalized light intensity was significantly different among different PE levels and formed a linear relationship ( , ). The study found that photons could reach the pulmonary artery tissue, the light intensity was linearly related to the degrees of embolism, PE could be quantitatively diagnosed by NIRS. Meanwhile, the optimized distance in between the light source and detector, 2.8 to 3.1cm, was recommended to be used in future potential noninvasive optical diagnosis of PE.