The Parametric Instabilities embedded Propagation and Evolution of RF Spectrum program has been developed to integrate the propagation and absorption of lower hybrid wave (LHW) around the scrape-off layer (SOL) of magnetically confined plasmas with parametric decay instabilities (PDIs). Based on an improved nonlocal convective amplification model of a quasi-mode PDI and the ray-tracing method, the program can be applied to self-consistently analyse the energy transfer and spectral evolution of LHW during its propagation. For the first time, we have successfully simulated the density limit and spectral broadening of LHW around the SOL across different experimental scenarios.