Abstract: Because of the good reflection of the ionosphere to the VLF wave, the heating model suitable for the high frequency band is not fully applicable in the VLF band. In this paper, the energy distribution law of the plane wave in the horizontal ionosphere layered structure is analyzed by the recursive algorithm of the full wave reflection coefficient and the wave mode amplitude. Combined with the ionosphere chemical reaction model, the energy distribution of the VLF wave in the ionosphere is analyzed. Compared with the methods used in other literatures, the full wave method takes into account the superposition of reflection and transmission waves of each layer. The simulation results show that the energy flow absorption and electronic temperature disturbance amplitude obtained by the full wave method are larger. The ionosphere disturbance under different incidence angle, azimuth angle and geomagnetic angle are compared, and it is found that changing the incident angle of the radio wave will change the reflection and absorption height of the radio wave in the ionosphere; while changing the azimuth angle and the geomagnetic inclination will only affect the absorption amplitude of the energy flow in the ionosphere, and will not affect the reflection height of the radio wave.