The invention discloses a resonance inclusion and a coupling prediction earthquake method thereof; the method comprises the following steps: obtaining earth electricity resonance inclusion information, and determining the magnitudes and epicenter distances according to the earth electric resonance inclusion information; obtaining terrestrial magnetism resonance inclusion information including thetime of blast chimney terrestrial magnetism fast and slow waves reaching the monitoring stations; determining the magnitudes and epicenter distances according to the terrestrial magnetism resonance inclusion information; using each monitoring center as the center, using corresponding epicenter distances as the radius so as to determine a circle or an arc, and using the intersection point formed byall circles or arcs to serve as the primary specific earthquake epicenter position; screening to determine each earthquake unique epicenter and a strong earthquake sequence. The method uses the earthquake inclusion resonance wave periods to determine the magnitude, uses the resonance inclusion growth to determine the origin time, uses quantification mathematics to calculate and accurately determine the earthquake magnitude, the epicenter and the origin time, thus accurately predicting earthquakes. The method uses blast chimney high frequency resonance waves to determine the pre-earthquake early warning time.