Method and system for determining azimuth of marsquake/moonquake event

1 . A method for determining an azimuth of a marsquake or moonquake event, comprising:

collecting original waveform data of a marsquake or moonquake by a plurality of seismometers in a single seismic station, wherein when the single seismic station is situated on Mars, the original waveform data collected by the plurality of the seismometers in the single seismic station is original waveform data of a marsquake; and when the single seismic station is situated on Moon, the original waveform data collected by the plurality of the seismometers in the single seismic station is original waveform data of a moonquake;

acquiring, by a computer, the original waveform data of the marsquake or moonquake from the plurality of the seismometers in the single seismic station;

filtering, by the computer, the original waveform data to obtain filtered waveform data;

intercepting, by the computer, filtered waveform data of a set time length to obtain first arrival waveform data of a P wave, wherein a starting time of the set time length is an arrival time of the P wave;

performing, by the computer, data analysis on the first arrival waveform data of the P wave, and determining the first arrival waveform data of the P wave of an east (E) component obtained after the data analysis as X-column data and the first arrival waveform data of the P wave of a north (N) component obtained after the data analysis as Y-column data; and

performing, by the computer, optimal fitting on the X-column data and the Y-column data using a direct least square ellipse fitting algorithm to obtain a first arrival elliptical locus of the P wave obtained after the optimal fitting, and determining an azimuth of the marsquake or moonquake relative to the single seismic station according to the first arrival elliptical locus of the P wave obtained after the optimal fitting.

2 . The method for determining an azimuth of a marsquake or moonquake event according to claim 1 , wherein the filtering the original waveform data to obtain filtered waveform data specifically comprises:

rotating the original waveform data to a ZNE orthogonal coordinate system to obtain ZNE coordinate system data; and

filtering the ZNE coordinate system data using a filter to obtain the filtered waveform data.

3 . The method for determining an azimuth of a marsquake or moonquake event according to claim 1 , wherein the determining an azimuth of the marsquake or moonquake relative to the single seismic station according to the first arrival elliptical locus of the P wave obtained after the optimal fitting specifically comprises:

determining a quadrant into which a major axis angle of the first arrival elliptical locus of the P wave obtained after the optimal fitting falls according to a positive or negative polarity of a maximum value of absolute values in the X-column data and a positive or negative polarity of a maximum value of absolute values in the Y-column data;

determining a major axis orientation of the first arrival elliptical locus of the P wave obtained after the optimal fitting according to the quadrant into which the major axis angle falls; and

converting the major axis orientation of the first arrival elliptical locus of the P wave obtained after the optimal fitting into a geographic azimuth, and determining the geographic azimuth as the azimuth of the marsquake or moonquake relative to the single seismic station.

4 . A system for determining an azimuth of a marsquake or moonquake event, comprising: a plurality of seismometers in a single seismic station and a computer,

wherein the plurality of the seismometers are configured to collect original waveform data of a marsquake or moonquake; when the single seismic station is situated on Mars, the original waveform data collected by the plurality of the seismometers in the single seismic station is original waveform data of a marsquake; and when the single seismic station is situated on Moon, the original waveform data collected by the plurality of the seismometers in the single seismic station is original waveform data of a moonquake

the computer is configured to acquire the original waveform data of the marsquake or moonquake from the plurality of the seismometers for executing operations comprising:

filtering the original waveform data to obtain filtered waveform data;

intercepting filtered waveform data of a set time length to obtain first arrival waveform data of a P wave, wherein a starting time of the set time length is an arrival time of the P wave;

performing data analysis on the first arrival waveform data of the P wave, and determining the first arrival waveform data of the P wave of an E component obtained after the data analysis as X-column data and the first arrival waveform data of the P wave of an N component obtained after the data analysis as Y-column data; and

performing optimal fitting on the X-column data and the Y-column data using a direct least square ellipse fitting algorithm to obtain a first arrival elliptical locus of the P wave obtained after the optimal fitting, and determining an azimuth of the marsquake or moonquake relative to the single seismic station according to the first arrival elliptical locus of the P wave obtained after the optimal fitting.

5 . The system for determining an azimuth of a marsquake or moonquake event according to claim 4 , wherein the filtering the original waveform data to obtain filtered waveform data specifically comprises:

rotating the original waveform data to a ZNE orthogonal coordinate system to obtain ZNE coordinate system data; and

filtering the ZNE coordinate system data using a filter to obtain filtered waveform data.

6 . The system for determining an azimuth of a marsquake or moonquake event according to claim 4 , wherein for determining an azimuth of the marsquake or moonquake relative to the single seismic station according to the first arrival elliptical locus of the P wave obtained after the optimal fitting comprises:

determining a quadrant into which a major axis angle of the first arrival elliptical locus of the P wave obtained after the optimal fitting falls according to a positive or negative polarity of a maximum value of absolute values in the X-column data and a positive or negative polarity of a maximum value of absolute values in the Y-column data;

determining a major axis orientation of the first arrival elliptical locus of the P wave obtained after the optimal fitting according to the quadrant into which the major axis angle falls; and

converting the major axis orientation of the first arrival elliptical locus of the P wave obtained after the optimal fitting into a geographic azimuth, and determine the geographic azimuth as the azimuth of the marsquake or moonquake relative to the single seismic station.