Abstract: The embodiments of the present application include acquiring a monitoring region and each observation point therein; partitioning the monitoring region into N layers of grids according to a seismic source positioning accuracy, wherein a side length of a grid cell of an i-th layer of grid is D/2i-1, i=1, . . . N, and D is an initial side length of the grid cell and not more than a double of a distance between the respective observation points; searching all nodes in a first layer of grid to acquire a node satisfying a preset condition therefrom; from i=2, determining and searching nodes satisfying a first preset requirement in the i-th layer of grid, to acquire a node satisfying the preset condition therefrom, until a search in an N-th layer of grid is completed, wherein a node satisfying the preset condition acquired in the N-th layer of grid is a seismic source point location.

Abstract: The embodiments of the present application include acquiring a monitoring region and each observation point therein; partitioning the monitoring region into N layers of grids according to a seismic source positioning accuracy, wherein a side length of a grid cell of an i-th layer of grid is D/2i-1, i=1, . . . N, and D is an initial side length of the grid cell and not more than a double of a distance between the respective observation points; searching all nodes in a first layer of grid to acquire a node satisfying a preset condition therefrom; from i=2, determining and searching nodes satisfying a first preset requirement in the i-th layer of grid, to acquire a node satisfying the preset condition therefrom, until a search in an N-th layer of grid is completed, wherein a node satisfying the preset condition acquired in the N-th layer of grid is a seismic source point location.

Abstract: Provided is a method for performing layer Q factor inversion by using an amplitude spectrum attribute of a down going wave of vertical seismic profile data in a geophysical exploration data processing technology. In the method, first an F-K (frequency-wave number) method is used to perform wave field separation on VSP original data, so as to obtain a down going wave; a down going wavelet and a monitoring wavelet are selected to undergone Fourier transform to obtain an amplitude spectrum, polynomial fitting is performed on the amplitude spectrum to obtain an equivalent Q, and a formula between the equivalent Q and a layer Q is used to perform inversion, so as to obtain the layer Q. The method has a strong capability of resisting random disturbance, and is capable of removing a difference of triggering wavelet. The algorithm is simple and can greatly save workload; moreover, the layer Q value obtained through inversion has a desirable stability and high precision.