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.