The invention provides a green function-SPWVD time-frequency analysis method of ENFEMF data. The green function-SPWVD time-frequency analysis method comprises the steps of: firstly, performing smooth compression pretreatment on original ENPEMF data; then adopting a green function based on a Poisson equation for processing to obtain ENFEMF data processed by green function; adopting a smooth pseudo WVD algorithm for frequency sweeping processing; and finally carrying out time-frequency analysis summary to obtain a time-frequency joint distribution map of the original data. The green function-SPWVD time-frequency analysis method fully utilizes the superposition idea of the green function, reduces noise and reduces cross terms generated during WVD frequency sweeping, is novel in structure and simple in step, can improve processing efficiency when aiming at data analysis of large number of multifarious earth natural pulse electromagnetic fields, and has good application prospect.