Method for high-resolution determination of seismic event P-wave reverse azimuth and slowness

Abstract

The invention discloses a method for high-resolution determination of seismic event P-wave reverse azimuth and slowness. The working steps of the method are that: supposing that the array is composedof N stations, taking the nearest station from the center of the array as the reference station and setting the station as the origin; the vertical component seismic displacement waveform records of the N stations are read, and instrument response removing and de-mean processing is performed on the records and the same sampling rate is obtained through resampling; the waveform records of all the stations in the step 2 are aligned according to the transmitting array time and the time window is selected; and the selected seismic records are windowed and subjected to Fourier transform so as to obtain the record spectrum, i.e. the observation vector. The method has the beneficial effects that the high-resolution seismic event P-wave reverse azimuth and slowness parameters can be provided by the method in comparison with the conventional station array FK algorithm; and the P-wave reverse azimuth and slowness parameters of the seismic event can be inversed quickly by using the constraint ofthe sparse characteristic of the sources in space distribution.

Description

technical field [0001] The invention belongs to the field of P-wave processing of seismic events, and relates to a seismic array data processing method, in particular to a method for determining the inverse azimuth and slowness of P-waves of seismic events with high resolution. Background technique [0002] The P wave in the seismic event is the longitudinal wave, and the P wave is the wave that is emitted from the source and refracted and propagated all the way. It is the first wave in the teleseismic measurement. [0003] Usually, methods for determining the P-wave inverse azimuth and slowness parameters of seismic events include frequency-domain beamforming (FK) and time-domain beamforming. For frequency-domain beamforming (FK) or time-domain beamforming (beamforming) methods to determine the P-wave back azimuth and slowness parameters of seismic events, the resolution is low, which is not convenient for staff to analyze and view, and through FK or beamforming The method...

AI Technical Summary

Problems solved by technology

[0004] The purpose of the present invention is to solve the problem that the resolution of the P-wave inverse azimuth and slowness parameters of seismic events is determined by the method of frequency domain beamforming (FK) or time domain beamforming (beamforming), which is not convenient for staff to analyze and view , and the P-wave inverse azimuth and slowness parameters of multiple seismic events cannot be distinguished by FK or beamforming methods, and a method for determining the P-wave inverse azimuth and slowness of seismic events with high resolution is proposed

Images