Se ha informado sobre radiación electromagnética anormal asociada con la actividad sísmica en una amplia gama de frecuencias, pero su energía principal se concentra en las bandas de frecuencia superbaja (SLF) y extremadamente baja (ELF). Para estimar el efecto de la fuente de radiación sísmica, se modeló un radiador en una fundación de suelo como un dipolo eléctrico horizontal (HED), y se estudiaron las características de propagación de los campos electromagnéticos en la cavidad Tierra-ionosfera. Las expresiones de los campos electromagnéticos se pudieron obtener de acuerdo con el teorema de reciprocidad. Por lo tanto, se propuso un nuevo algoritmo llamado algoritmo integral numérico, que es adecuado tanto para las bandas SLF como ELF. El nuevo algoritmo se comparó con el algoritmo de aproximación asintótica cuando el punto de recepción no estaba cerca de la fuente de campo y el antípoda. Se encontró que los dos algoritmos estaban en excelente acuerdo, confirmando la validez del nuevo algoritmo para las bandas SLF y ELF.