The invention discloses a seismic
observation data automatic calibration method based on median filtering, and the method comprises the steps: firstly analyzing an original
signal, determining a time window of a seismic
signal, extracting a
signal segment in the time window of the seismic signal, carrying out the band-pass filtering of the signal segment according to the type of the original signal, and obtaining an initial band-pass filtering signal;
processing the initial band-pass filtering signal by using median filtering windows with different lengths, calculating indexes such as a signal-to-
noise ratio and an effective signal
retention rate, and screening an optimal
initial window through second-order difference and sliding standard deviation; and
processing the initial band-pass filtering signal by using an optimal signal-to-
noise ratio window to obtain an initial median filtering signal, calculating a point-by-
point correlation coefficient of the initial band-pass filtering signal and a residual error, carrying out regional dynamic window adjustment according to the
correlation coefficient, carrying out median filtering
processing again, and finally carrying out band-pass filtering processing again to obtain a final calibration signal. According to the design, the selected length of the window is not fixed, signals of different sampling rates and observation frequency bands can be adapted, the universality is high, and the
noise reduction effect and effective signal retention are considered at the same time.