CN110058299A - Earthquake positioning method and device and terminal equipment - Google Patents

Abstract

The embodiment of the present invention is applicable to the field of seismic technology, and discloses a seismic positioning method, device, terminal equipment and storage medium, wherein the method includes: calculating the surface wave empirical Green function between pairs of stations through the cross-correlation of seismic background noise; calculating The first group travel time of the surface wave in the empirical Green function of the surface wave between the first type of station and the second type of station, and the second group travel time of the seismic surface wave is calculated; according to the first group travel time, the earthquake occurrence area is constructed by interpolation The group travel time table of surface waves to the stations in the target area; the waveform of the surface wave empirical Green function between the pairs of stations is corrected to the seismic waveform of the virtual hypocenter, and the depth travel time correction scale is constructed; by constructing the objective function, the hypocenter is determined Location. The embodiments of the present invention can enhance the uniqueness of the positioning solution, reduce the dependence of the positioning accuracy on the velocity model, the given hypocenter depth and the hypocenter mechanism, and improve the positioning efficiency, positioning accuracy and positioning accuracy.

Description

Seismic positioning method, device and terminal equipment

technical field

The present invention belongs to the technical field of earthquakes, and in particular relates to an earthquake positioning method, device, terminal device and computer-readable storage medium.

Background technique

Earthquake location is one of the most classic and basic problems in the field of seismology. Accurate earthquake location results are not only of great significance to the study of seismic activity structure, the internal structure of the earth, the rupture process and geometric structure of the source, etc., but also to the disaster relief after the earthquake. Work is also crucial.

The traditional earthquake location method calculates the travel time or waveform of a given velocity model, and finds the spatial point that best matches the observed travel time or waveform as the earthquake occurrence location. This method is highly dependent on the accuracy of the velocity model. In order to reduce the serious influence of the velocity model error on the positioning results, the noise surface wave Green function extracted between pairs of stations can be used to constrain the earthquake location, which usually includes the following steps: selecting stations near the earthquake epicenter, Correlation, calculate the surface wave empirical Green function between the near station and the far station far away from the earthquake location, the surface wave empirical Green function includes the Rayleigh wave empirical Green function EGF _Ray and the Love wave empirical Green function EGF _Lov ; 3D velocity model, synthesizing the surface wave Green functions between near and far stations, including GF _Ray and GF _Lov , and given the focal depth and focal mechanism at the same time, calculate the synthetic theoretical seismogram; and then calculate the Rayleigh wave separately through waveform cross-correlation and Love wave traveltime corrections; finally, combined with the seismic records actually recorded by the far station and the synthesized theoretical seismogram, the traveltime corrections of Rayleigh waves and Love waves are introduced, and the two traveltime corrections are reversed by least squares iteration. to determine the epicenter of the earthquake.

It can be seen that the existing method of using surface wave seismic positioning is carried out on the premise that the general position of the earthquake to be located is known in advance, and the prior information is required to determine which nearby station to use for traveltime correction. When multiple nearby stations are close to the earthquake distance, multiple earthquake epicenter positions will be obtained by using different stations, and the solution is not unique; and when there is no suitable station near the earthquake location or the distance between the stations is relatively far, the travel time correction amount at this time There is a large deviation from the actual required correction amount. Especially when the source is deep, the deviation between the travel time correction between the near and far stations on the surface and the actual correction required by the deep source is large, and the positioning is inaccurate. In addition, waveform cross-correlation calculation needs to be performed, because the calculation amount of waveform cross-correlation calculation is large, which seriously affects the calculation efficiency. To sum up, the location solutions of the existing seismic location methods may not be unique, and depend on the given focal depth and focal mechanism, and the location accuracy, precision, and location efficiency are low.

SUMMARY OF THE INVENTION

In view of this, embodiments of the present invention provide a seismic positioning method, device, terminal device, and computer-readable storage medium, so as to solve the problem that the localized seismic source of the existing seismic positioning method may not be unique, and the positioning accuracy depends on a given The accuracy of the focal depth and focal mechanism, the positioning accuracy, precision and efficiency are all low, and the real focal mechanism of the earthquake can also be inverted through waveform fitting.

A first aspect of the embodiments of the present invention provides an earthquake positioning method, including:

Through the cross-correlation of the seismic background noise recorded by the stations, the surface wave empirical Green function between the station pairs is calculated; wherein the station pairs include the stations belonging to the first type of stations and the stations belonging to the second type of stations. , the first type of station is a station located in a seismic zone, the second type of station is a station located in a target area with a distance from the seismic zone greater than a preset distance threshold, and the A dense array is set up in the earthquake-generating area;

Calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station, and calculate the seismic surface according to the surface wave data recorded by the second type of station the travel time of the second group of waves;

According to the first group travel time, construct the group travel time table of the surface waves from the seismic region to the station in the target area by interpolation;

Correcting the waveform of the surface wave empirical Green function between the pair of stations to the seismic waveform of the double-couple type virtual hypocenter, and constructing a depth traveltime correction scale according to the group travel time of the surface wave of the virtual hypocenter;

An objective function is constructed according to the first group travel time, the second group travel time, the depth travel time correction scale, and a preset correction time amount, and the epicenter position is determined according to the objective function.

Optionally, after the determining the location of the epicenter, the method further includes:

According to the seismic waveform of the virtual epicenter, the focal depth and the focal mechanism are determined.

Optionally, constructing a depth traveltime correction scale according to the group traveltime of the surface wave of the virtual source, including:

calculating the travel time of the third group of surface waves of the virtual hypocenter;

calculating the difference between the travel time of the first group and the travel time of the third group;

According to the difference value, the focal depth is traversed, and the depth traveltime correction scale is constructed.

Optionally, the determining the location of the epicenter according to the objective function includes:

Based on the objective function Φ(X,Y,h)=min||T _EQ -T _EGFs +T _Dep +T ₀ || ₂ , the traveltime residual of the surface wave group is minimized, and the source is determined by a grid search algorithm where, T _EGFS is the first group travel time, T _EQ is the second group travel time, T _Dep is the depth travel time correction scale, T ₀ is the preset correction time amount, the preset The amount of correction time is the amount of time that needs to be corrected based on the source time function and the existence of the seismic moment.

Optionally, calculating the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station, and according to the recorded time of the second type of station. Surface wave data calculates the second group traveltime of seismic surface waves, including:

Calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station by performing narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave;

By performing narrowband frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave, the second group travel time of the seismic surface wave is calculated according to the surface wave data recorded by the second type of station.

A second aspect of the embodiments of the present invention provides a seismic positioning device, including:

A calculation module for calculating the surface wave empirical Green function between pairs of stations through cross-correlation of seismic background noises recorded by stations; wherein the pairs of stations include stations belonging to the first type of stations and stations belonging to the second type The station of the station, the first type of station is the station set in the seismic zone, and the second type of station is set in the target area whose distance from the seismic zone is greater than the preset distance threshold. a station, wherein the earthquake-generating area is provided with a dense array;

The travel time calculation module is used to calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station, and according to the recorded time of the second type of station Surface wave data to calculate the second group travel time of seismic surface waves;

a travel time table calculation module, configured to construct a group travel time table of surface waves from the seismic region to the station in the target area by interpolation according to the first group travel time;

A building module is used to correct the waveform of the surface wave empirical Green function between the pair of stations to the seismic waveform of the double-couple type virtual hypocenter, and build a depth traveltime correction scale according to the group travel time of the surface wave of the virtual hypocenter ;

The first determination module is configured to construct an objective function according to the first group travel time, the second group travel time, the depth travel time correction scale and a preset correction time amount, and determine the epicenter position according to the objective function.

Optionally, also include:

The second determining module is configured to determine the focal depth and the focal mechanism according to the seismic waveform of the virtual hypocenter.

Optionally, the building blocks include:

a computing unit for computing the third group travel time of the surface wave of the virtual hypocenter;

a difference calculation unit, configured to calculate the difference between the first group travel time and the third group travel time;

A construction unit, configured to traverse the hypocenter depth according to the difference value, and construct the depth traveltime correction scale.

A third aspect of the embodiments of the present invention provides a terminal device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, which is implemented when the processor executes the computer program The steps of the method according to any one of the above first aspects.

A fourth aspect of the embodiments of the present invention provides a computer-readable storage medium, where the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, implements the method according to any one of the foregoing first aspects A step of.

The beneficial effects that the embodiment of the present invention has compared with the prior art are:

The embodiment of the present invention provides more data for positioning by setting up dense arrays in the earthquake-generating area, enhances the uniqueness of understanding, and improves positioning accuracy and positioning accuracy; the deviation caused by the velocity model is corrected by using the seismic background noise, And by correcting the waveform of the surface wave empirical Green function as a virtual source, the dependence of the positioning accuracy on the given focal depth and focal mechanism is reduced, and the accuracy and positioning accuracy are improved. The amount of calculation improves the positioning efficiency.

Description of drawings

In order to illustrate the technical solutions in the embodiments of the present invention more clearly, the following briefly introduces the accompanying drawings that need to be used in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only for the present invention. In some embodiments, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without any creative effort.

1 is a schematic flowchart of a seismic positioning method according to an embodiment of the present invention;

FIG. 2 is a schematic diagram of station distribution according to an embodiment of the present invention;

3 is a schematic diagram of an earthquake positioning effect provided by an embodiment of the present invention;

FIG. 4 is a schematic block diagram of another flowchart of an earthquake positioning method provided by an embodiment of the present invention;

5 is a schematic block diagram of the structure of a seismic positioning device provided by an embodiment of the present invention;

FIG. 6 is a schematic diagram of a terminal device according to an embodiment of the present invention.

Detailed ways

In the following description, for the purpose of illustration rather than limitation, specific details such as specific system structures and technologies are set forth in order to provide a thorough understanding of the embodiments of the present invention. However, it will be apparent to those skilled in the art that the present invention may be practiced in other embodiments without these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.

In order to illustrate the technical solutions of the present invention, the following specific embodiments are used for description.

Please refer to FIG. 1 , which is a schematic flowchart of an earthquake positioning method according to an embodiment of the present invention. The method may include the following steps:

Step S101: Calculate the surface wave empirical Green's function between the pair of stations through the cross-correlation of the seismic background noise recorded by the stations; wherein, the pair of stations includes the stations belonging to the first type of stations and the stations belonging to the second type of stations. The first type of station is the station set in the earthquake-generating area, the second type of station is the station set in the target area whose distance from the earthquake-generating area is greater than the preset distance threshold, and the earthquake-generating area is set with dense stations array.

It can be understood that the seismic background noise refers to the seismic coda signal and weak seismic signal recorded at the station, which is distinguished from the seismic signal with a larger magnitude. Because the signal is very weak and has the same amplitude as the noise, it is called background. noise.

It should be noted that the above-mentioned earthquake-generating areas refer to areas with relatively active earthquakes, which are relatively close to the epicenter. The first type of stations located in the seismic zone can be called near stations. The first-type stations are set up in the earthquake-generating area, and a plurality of the first-type stations are in a dense array. It can be seen that the uniqueness of the positioning solution is enhanced by using the dense array in the seismic area to obtain more positioning data.

The second type of station can be called a remote station, which is a station far away from the seismic zone. The remote station is set in the target area, and the distance between the target area and the seismic zone is greater than the preset distance threshold. The preset distance The threshold can be set according to the actual application scenario.

The above station pairs include near stations and far stations, and the surface wave empirical Green's Function (EGFs) between the station pairs is calculated, that is, the surface wave empirical Green's function EGFs between the near and far stations is calculated. The calculation process can be calculated by the following formula (1).

The left-hand term of formula (1) is is the derivative of the background noise cross-correlation to time t, and the right-hand term is is the surface wave EGFs between station A and station B, where C _AB is the cross-correlation function between station A and station B.

In order to better introduce the station distribution situation, the following description will be given with reference to the schematic diagram of the station distribution provided by the embodiment of the present invention shown in FIG. 2 .

As shown in Fig. 2, in the plane composed of the X axis and the Y axis, "○" represents the center position, "□" represents the far stage, and "△" represents the near stage. Near the center is the earthquake-generating area, and there are multiple near-stations, that is, a dense array, and at a position far from the center, there are multiple far-stations scattered in a staggered manner.

It can be understood that the above-mentioned FIG. 2 is only an exemplary distribution diagram, and the station distribution may also be in other forms, which are not limited herein.

Step S102: Calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station, and calculate the first group of surface waves according to the surface wave data recorded by the second type of station. When the two groups leave.

Understandably, earthquake travel time refers to the time it takes for seismic waves to travel from the epicenter to the observation point. The group travel time of a surface wave refers to the travel time of the envelope or energy of the surface wave.

In specific applications, the corresponding surface wave group travel time can be calculated by performing narrow-band frequency domain Gaussian filtering and Hilbert transform on the surface wave waveform.

In an embodiment, the group travel time can be obtained by performing narrow-band frequency domain Gaussian filtering and Hilbert transform on the waveform, so optionally, in this step, the above-mentioned calculation of the plane between the first type of station and the second type of station. The specific process of calculating the traveltime of the first group of surface waves in the empirical Green function of the wave and calculating the traveltime of the second group of seismic surface waves according to the surface wave data recorded by the second type of station can be as follows: Domain high-pass Gaussian filtering and Hilbert transform to calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station; through the narrowband frequency domain high-pass Gaussian filtering of the surface wave waveform and Hilbert transform to calculate the travel time of the second group of seismic surface waves based on the surface wave data recorded by the second type of stations.

Among them, the narrow-band frequency-domain Gaussian filtering formula is formula (2) to formula (4), and the Hilbert transform is the following formula 3.

in,

w is the original waveform data before transformation, is the Fourier transform of w, H is the Gaussian filter factor, Hilbert-transformed data for the filtered waveform.

Step S103 , according to the first group travel time, construct a group travel time table of the surface waves between the earthquake generating area and the station in the target area by interpolation.

It is understood that the travel time table refers to the timetable of the propagation of seismic waves at different epicentral distances. Although a dense array has been set up in the earthquake-generating area, dividing the earthquake-generating area into finer grids for interpolation and constructing a travel time table can further improve the positioning accuracy and positioning accuracy.

Step S104 , correcting the waveform of the surface wave empirical Green function between the pair of stations to the seismic waveform of the double-couple type virtual hypocenter, and constructing a depth traveltime correction scale according to the group traveltime of the surface wave of the virtual hypocenter.

Specifically, the waveform of the surface wave empirical Green function between pairs of stations, that is, the surface wave EGFs waveform obtained from the seismic background noise, is corrected to the seismic waveform of a double-couple type virtual hypocenter with a certain depth. Then, calculate the difference between the surface wave group travel time of the virtual hypocenter and the EGFs surface wave travel time, that is, the first group travel time, and then traverse the focal depth to obtain the surface wave travel time correction scale that needs to be corrected due to the real focal depth. .

In the specific application, the frequency domain correction formula from the EGFs waveform to the virtual hypocenter waveform is as follows:

in, and are the frequency domain representation of the corrected three-component waveform, r ₁ , r ₂ and l ₁ are the eigenfunction values of the Rayleigh wave and Love wave, respectively, M is the focal mechanism solution of the earthquake, is the empirical Green's function of surface waves corresponding to different components.

Through the above formulas (6) to (8), the surface wave EGFs obtained by the seismic background noise can be corrected to the seismic waveform of the virtual hypocenter.

Step S105 , construct an objective function according to the first group travel time, the second group travel time, the depth travel time correction scale and the preset correction time amount, and determine the epicenter position according to the objective function.

Specifically, by traversing the depth of the hypocenter, an objective function is constructed to minimize the traveltime residual of the surface wave group, and then the location of the hypocenter can be determined by a grid search algorithm.

In an embodiment, optionally, the above-mentioned specific process of determining the location of the epicenter according to the objective function may include: based on the objective function Φ(X, Y, h)=min||T _EQ -T _EGFs +T _Dep +T ₀ || ₂ , minimize the surface wave group traveltime residual, and determine the hypocenter location by grid search algorithm; among them, T _EGFS is the first group travel time, T _EQ is the second group travel time, T _Dep is the depth travel time correction scale, T ₀ is the preset correction time amount, that is, T ₀ is the time amount that needs to be corrected due to the existence of the source time function and the earthquake occurrence time.

In order to better introduce the positioning result, the following description will be given with reference to the schematic diagram of the earthquake positioning effect provided by the embodiment of the present invention shown in FIG. 3 . As shown in Figure 3, the closer to the epicenter, the darker the color, and the deepest point is the true epicenter. The dots in the epicenter are the epicenter positions determined by this method, and it can be seen that it completely coincides with the real epicenter. The schematic diagram of the positioning effect is an effect diagram obtained by using the earthquake positioning method provided by the embodiment of the present invention, and it can be seen that the epicenter position is accurately located.

In this embodiment, by setting up a dense array in the earthquake-generating area, more data is provided for positioning, the uniqueness of understanding is enhanced, and the positioning accuracy and positioning accuracy are improved; the deviation caused by the velocity model is corrected by using the background noise of the earthquake , and by correcting the waveform of the surface wave empirical Green function as a virtual hypocenter, the dependence of the positioning accuracy on the velocity model, the given hypocenter depth and hypocenter mechanism is reduced, and the accuracy and positioning accuracy are improved; no waveform cross-correlation is required. The calculation reduces the amount of calculation and improves the positioning efficiency.

Embodiment 2

Please refer to FIG. 4 , which is another schematic flowchart of the seismic positioning method provided by the embodiment of the present invention. The method may include the following steps:

Step S401: Calculate the surface wave empirical Green's function between pairs of stations through the cross-correlation of the seismic background noise recorded by the stations.

Step S402: Calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station, and calculate the first group of surface waves according to the surface wave data recorded by the second type of station. When the two groups leave.

Step S403 , according to the first group travel time, construct a group travel time table of the surface waves between the stations from the seismic region to the target area by interpolation.

It should be noted that, the above steps S301 to S303 are the same as the steps S201 to S203 of the above-mentioned first embodiment. For related introduction, please refer to the corresponding content above, which will not be repeated here.

Step S404, correcting the waveform of the surface wave empirical Green function between the pair of stations to the seismic waveform of the double-couple type virtual hypocenter.

Step S405: Calculate the traveltime of the third group of the surface wave group of the virtual hypocenter, calculate the difference between the traveltime of the first group and the traveltime of the third group, traverse the depth of the source according to the difference, and construct a depth traveltime correction scale.

Step S406: Construct an objective function according to the first group travel time, the second group travel time, the depth travel time correction scale and the preset correction time amount, and determine the epicenter position according to the objective function.

Step S407: Determine the focal depth and the focal mechanism according to the seismic waveform of the virtual hypocenter.

Specifically, according to the seismic waveform of the virtual hypocenter, the true depth of the hypocenter and the mechanism of the hypocenter can be simultaneously determined through traversal algorithm and waveform fitting.

In this embodiment, by setting up a dense array in the earthquake-generating area, more data is provided for positioning, the uniqueness of understanding is enhanced, and the positioning accuracy and positioning accuracy are improved; the deviation caused by the velocity model is corrected by using the background noise of the earthquake , and by correcting the waveform of the surface wave empirical Green function as a virtual hypocenter, the dependence of the positioning accuracy on the velocity model, the given hypocenter depth and hypocenter mechanism is reduced, and the accuracy and positioning accuracy are improved; no waveform cross-correlation is required. The calculation reduces the amount of calculation and improves the positioning efficiency.

It should be understood that the size of the sequence numbers of the steps in the above embodiments does not mean the sequence of execution, and the execution sequence of each process should be determined by its functions and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.

Embodiment 3

Please refer to FIG. 5 , which is a schematic structural block diagram of a seismic positioning apparatus provided by an embodiment of the present invention, and the apparatus may include:

The calculation module 51 is used to calculate the surface wave empirical Green function between the station pairs through the cross-correlation of the seismic background noise recorded by the stations; wherein, the station pairs include the stations belonging to the first type of stations and the stations belonging to the second type. The station of the station, the first type of station is the station set in the seismic zone, the second type of station is the station set in the target area whose distance from the seismic zone is greater than the preset distance threshold, the seismic zone is set There is a dense array;

The travel time calculation module 52 is used to calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station, and calculate the earthquake according to the surface wave data recorded by the second type of station The second group travel time of the surface wave;

The travel time table calculation module 53 is used for constructing the group travel time table of the surface waves from the seismic region to the station in the target area by interpolation according to the first group travel time;

The building module 54 is used to correct the waveform of the surface wave empirical Green function between the station pairs to the seismic waveform of the double-couple type virtual hypocenter, and build a depth traveltime correction scale according to the group travel time of the surface wave of the virtual hypocenter;

The first determination module 55 is configured to construct an objective function according to the first group travel time, the second group travel time, the depth travel time correction scale and the preset correction time amount, and determine the epicenter position according to the objective function.

In a possible implementation, the above device may further include:

The second determining module is used for determining the focal depth and the focal mechanism according to the seismic waveform of the virtual hypocenter.

In one possible implementation, the above building blocks may include:

a calculation unit for calculating the travel time of the third group of surface waves of the virtual hypocenter;

A difference calculation unit for calculating the difference between the first group of travel times and the third group of travel times;

The construction unit is used to traverse the focal depth according to the difference value, and construct the depth traveltime correction scale.

In a possible implementation, the above-mentioned first determining module may include:

A determination unit for minimizing the traveltime residual of the surface wave group based on the objective function Φ(X,Y,h)=min||T _EQ -T _EGFs +T _Dep +T ₀ || ₂ , determined by a grid search algorithm The location of the hypocenter; among them, T _EGFS is the first group travel time, T _EQ is the second group travel time, T _Dep is the depth travel time correction scale, and T ₀ is the preset correction time amount.

In a possible implementation, the above travel time calculation module may include:

The first travel time calculation unit is used to calculate the first time of the surface wave in the empirical Green function of the surface wave between the first type of station and the second type of station by performing narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the surface wave waveform. when a group goes

The second traveltime calculation unit calculates the second group traveltime of seismic surface waves according to the surface wave data recorded by the second type of station by performing narrowband frequency domain high-pass Gaussian filtering and Hilbert transform on the surface wave waveform.

In this embodiment, by setting up a dense array in the earthquake-generating area, more data is provided for positioning, the uniqueness of understanding is enhanced, and the positioning accuracy and positioning accuracy are improved; the deviation caused by the velocity model is corrected by using the background noise of the earthquake , and by correcting the waveform of the surface wave empirical Green function as a virtual hypocenter, the dependence of the positioning accuracy on the velocity model, the given hypocenter depth and hypocenter mechanism is reduced, and the accuracy and positioning accuracy are improved; no waveform cross-correlation is required. The calculation reduces the amount of calculation and improves the positioning efficiency.

Embodiment 4

FIG. 6 is a schematic diagram of a terminal device provided by an embodiment of the present invention. As shown in FIG. 6 , the terminal device 6 in this embodiment includes: a processor 60 , a memory 61 , and a computer program 62 stored in the memory 61 and running on the processor 60 . When the processor 60 executes the computer program 62 , the steps in each of the above-mentioned embodiments of the seismic positioning method are implemented, for example, steps S101 to S105 shown in FIG. 1 . Alternatively, when the processor 60 executes the computer program 62, the functions of the modules or units in the above-mentioned apparatus embodiments, for example, the functions of the modules 51 to 55 shown in FIG. 5 are implemented.

Exemplarily, the computer program 62 can be divided into one or more modules or units, and the one or more modules or units are stored in the memory 61 and executed by the processor 60 to complete the this invention. The one or more modules or units may be a series of computer program instruction segments capable of performing specific functions, and the instruction segments are used to describe the execution process of the computer program 62 in the terminal device 6 . For example, the computer program 62 can be divided into a calculation module, a travel time calculation module, a travel time table calculation module, a construction module and a first determination module, and the specific functions of each module are as follows:

A calculation module for calculating the surface wave empirical Green function between pairs of stations through cross-correlation of seismic background noises recorded by stations; wherein the pairs of stations include stations belonging to the first type of stations and stations belonging to the second type The station of the station, the first type of station is the station set in the seismic zone, and the second type of station is set in the target area whose distance from the seismic zone is greater than the preset distance threshold. a station, a dense array is arranged in the earthquake-generating area; a travel time calculation module is used to calculate the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station travel time, and calculate the travel time of the second group of seismic surface waves according to the surface wave data recorded by the second type of station; travel time table calculation module is used for constructing the earthquake occurrence area by interpolation according to the travel time of the first group The group travel time table of the surface waves to the stations in the target area; the building module is used to correct the waveform of the surface wave empirical Green function between the station pairs to the seismic waveform of the double-couple type virtual source, and according to the According to the group travel time of the surface wave of the virtual hypocenter, a depth travel time correction scale is constructed; the first determination module is used for correcting the depth travel time according to the first group travel time, the second group travel time, the depth travel time correction scale and the preset correction. The amount of time constructs an objective function, and based on the objective function, the location of the hypocenter is determined.

The terminal device 6 may be a computing device such as a desktop computer, a notebook, a palmtop computer, and a cloud server. The terminal device may include, but is not limited to, the processor 60 and the memory 61 . Those skilled in the art can understand that FIG. 6 is only an example of the terminal device 6, and does not constitute a limitation on the terminal device 6, and may include more or less components than the one shown, or combine some components, or different components For example, the terminal device may further include an input and output device, a network access device, a bus, and the like.

The so-called processor 60 may be a central processing unit (Central Processing Unit, CPU), and may also be other general-purpose processors, digital signal processors (Digital Signal Processors, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), Off-the-shelf programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, and the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.

The memory 61 may be an internal storage unit of the terminal device 6 , such as a hard disk or a memory of the terminal device 6 . The memory 61 may also be an external storage device of the terminal device 6, such as a plug-in hard disk equipped on the terminal device 6, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, flash card (Flash Card) and so on. Further, the memory 61 may also include both an internal storage unit of the terminal device 6 and an external storage device. The memory 61 is used to store the computer program and other programs and data required by the terminal device. The memory 61 can also be used to temporarily store data that has been output or will be output.

Those skilled in the art can clearly understand that, for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example. Module completion, that is, dividing the internal structure of the device into different functional units or modules to complete all or part of the functions described above. Each functional unit and module in the embodiment may be integrated in one processing unit, or each unit may exist physically alone, or two or more units may be integrated in one unit, and the above-mentioned integrated units may adopt hardware. It can also be realized in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing from each other, and are not used to limit the protection scope of the present application. For the specific working processes of the units and modules in the above-mentioned system, reference may be made to the corresponding processes in the foregoing method embodiments, which will not be repeated here.

In the foregoing embodiments, the description of each embodiment has its own emphasis. For parts that are not described or described in detail in a certain embodiment, reference may be made to the relevant descriptions of other embodiments.

Those of ordinary skill in the art can realize that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Skilled artisans may implement the described functionality using different methods for each particular application, but such implementations should not be considered beyond the scope of the present invention.

In the embodiments provided by the present invention, it should be understood that the disclosed apparatus, terminal device, and method may be implemented in other manners. For example, the above-described embodiments of the apparatus and terminal equipment are only illustrative. For example, the division of the modules or units is only a logical function division. In actual implementation, there may be other division methods, such as multiple units. Or components may be combined or may be integrated into another system, or some features may be omitted, or not implemented. On the other hand, the shown or discussed mutual coupling or direct coupling or communication connection may be through some interfaces, indirect coupling or communication connection of devices or units, and may be in electrical, mechanical or other forms.

The units described as separate components may or may not be physically separated, and components displayed as units may or may not be physical units, that is, may be located in one place, or may be distributed to multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution in this embodiment.

In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically alone, or two or more units may be integrated into one unit. The above-mentioned integrated units may be implemented in the form of hardware, or may be implemented in the form of software functional units.

The integrated modules or units, if implemented in the form of software functional units and sold or used as independent products, may be stored in a computer-readable storage medium. Based on this understanding, the present invention can implement all or part of the processes in the methods of the above embodiments, and can also be completed by instructing relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer When the program is executed by the processor, the steps of the foregoing method embodiments can be implemented. Wherein, the computer program includes computer program code, and the computer program code may be in the form of source code, object code, executable file or some intermediate form, and the like. The computer-readable medium may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer memory, a read-only memory (ROM, Read-Only Memory) , Random Access Memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium, etc. It should be noted that the content contained in the computer-readable media may be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction, for example, in some jurisdictions, according to legislation and patent practice, the computer-readable media Electric carrier signals and telecommunication signals are not included.

The above-mentioned embodiments are only used to illustrate the technical solutions of the present invention, but not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: it is still possible to implement the foregoing implementations. The technical solutions described in the examples are modified, or some technical features thereof are equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should be included in the within the protection scope of the present invention.

Claims (10)

1. an earthquake location method, is characterized in that, comprises: Through the cross-correlation of the seismic background noise recorded by the stations, the surface wave empirical Green function between the station pairs is calculated; wherein the station pairs include the stations belonging to the first type of stations and the stations belonging to the second type of stations. , the first type of station is a station located in a seismic zone, the second type of station is a station located in a target area with a distance from the seismic zone greater than a preset distance threshold, and the A dense array is set up in the earthquake-generating area; Calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station, and calculate the seismic surface according to the surface wave data recorded by the second type of station the travel time of the second group of waves; According to the first group travel time, construct the group travel time table of the surface waves from the seismic region to the station in the target area by interpolation; Correcting the waveform of the surface wave empirical Green function between the pair of stations to the seismic waveform of the double-couple type virtual hypocenter, and constructing a depth traveltime correction scale according to the group travel time of the surface wave of the virtual hypocenter; An objective function is constructed according to the first group travel time, the second group travel time, the depth travel time correction scale, and a preset correction time amount, and the epicenter position is determined according to the objective function. 2. The earthquake location method according to claim 1, characterized in that, after said determining the location of the epicenter, further comprising: According to the seismic waveform of the virtual epicenter, the focal depth and the focal mechanism are determined. 3. The earthquake location method according to claim 1, wherein the construction of a depth traveltime correction scale according to the group traveltime of the surface wave of the virtual hypocenter, comprises: calculating the travel time of the third group of surface waves of the virtual hypocenter; calculating the difference between the travel time of the first group and the travel time of the third group; According to the difference value, the focal depth is traversed, and the depth traveltime correction scale is constructed. 4. The earthquake location method according to any one of claims 1 to 3, wherein the determining the epicenter position according to the objective function comprises: Based on the objective function Φ(X,Y,h)=min||T _EQ -T _EGFs +T _Dep +T ₀ || ₂ , the traveltime residual of the surface wave group is minimized, and the source is determined by a grid search algorithm where, T _EGFS is the first group travel time, T _EQ is the second group travel time, T _Dep is the depth travel time correction scale, T ₀ is the preset correction time amount, the preset The amount of correction time is the amount of time that needs to be corrected based on the source time function and the existence of the seismic moment. 5 . The seismic positioning method according to claim 4 , wherein the calculating the first group travel time of surface waves in the surface wave empirical Green function between the first type of station and the second type of station. 6 . , and calculate the second group traveltime of seismic surface waves according to the surface wave data recorded by the second type of station, including: Calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station by performing narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave; By performing narrowband frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave, the second group travel time of the seismic surface wave is calculated according to the surface wave data recorded by the second type of station. 6. A seismic positioning device, characterized in that, comprising: A calculation module for calculating the surface wave empirical Green function between pairs of stations through cross-correlation of seismic background noises recorded by stations; wherein the pairs of stations include stations belonging to the first type of stations and stations belonging to the second type The station of the station, the first type of station is the station set in the seismic zone, and the second type of station is set in the target area whose distance from the seismic zone is greater than the preset distance threshold. a station, wherein the earthquake-generating area is provided with a dense array; The travel time calculation module is used to calculate the travel time of the first group of surface waves in the surface wave empirical Green function between the first type of station and the second type of station, and according to the recorded time of the second type of station Surface wave data to calculate the second group travel time of seismic surface waves; a travel time table calculation module, configured to construct a group travel time table of surface waves from the seismic region to the station in the target area by interpolation according to the first group travel time; A building module is used to correct the waveform of the surface wave empirical Green function between the pair of stations to the seismic waveform of the double-couple type virtual hypocenter, and build a depth traveltime correction scale according to the group travel time of the surface wave of the virtual hypocenter ; The first determination module is configured to construct an objective function according to the first group travel time, the second group travel time, the depth travel time correction scale and a preset correction time amount, and determine the epicenter position according to the objective function. 7. The seismic positioning device according to claim 6, further comprising: The second determining module is configured to determine the focal depth and the focal mechanism according to the seismic waveform of the virtual hypocenter. 8. The seismic positioning device of claim 6, wherein the building module comprises: a computing unit for computing the third group travel time of the surface wave of the virtual hypocenter; a difference calculation unit, configured to calculate the difference between the first group travel time and the third group travel time; A construction unit, configured to traverse the hypocenter depth according to the difference value, and construct the depth traveltime correction scale. 9. A terminal device, comprising a memory, a processor, and a computer program stored in the memory and running on the processor, the processor implementing the computer program as claimed in the claims Steps of any one of 1 to 5 of the method. 10. A computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 5 are implemented .