CN110058299A - Earthquake positioning method and device and terminal equipment - Google Patents
Earthquake positioning method and device and terminal equipment Download PDFInfo
- Publication number
- CN110058299A CN110058299A CN201811071892.9A CN201811071892A CN110058299A CN 110058299 A CN110058299 A CN 110058299A CN 201811071892 A CN201811071892 A CN 201811071892A CN 110058299 A CN110058299 A CN 110058299A
- Authority
- CN
- China
- Prior art keywords
- travel time
- station
- seismic
- surface wave
- group travel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000012937 correction Methods 0.000 claims abstract description 45
- 238000004613 tight binding model Methods 0.000 claims abstract description 29
- 230000007246 mechanism Effects 0.000 claims abstract description 17
- 230000006870 function Effects 0.000 claims description 55
- 238000004364 calculation method Methods 0.000 claims description 28
- 238000004590 computer program Methods 0.000 claims description 19
- 238000001914 filtration Methods 0.000 claims description 10
- 238000010276 construction Methods 0.000 claims description 6
- 238000010845 search algorithm Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/65—Source localisation, e.g. faults, hypocenters or reservoirs
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Emergency Management (AREA)
- Business, Economics & Management (AREA)
- Acoustics & Sound (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the invention is suitable for the technical field of earthquakes and discloses an earthquake positioning method, an earthquake positioning device, terminal equipment and a storage medium, wherein the method comprises the following steps: calculating a surface wave empirical Green function between station pairs through seismic background noise cross-correlation; calculating a first group travel time of a surface wave in a surface wave empirical Green's function between the first type station and the second type station, and calculating a second group travel time of a seismic surface wave; according to the first group travel time, a group travel time table of surface waves from the earthquake-generating area to the station of the target area is constructed through interpolation; correcting the waveform of the surface wave empirical Green function between the station pairs into the seismic waveform of a virtual seismic source, and constructing a depth travel time correction scale; the source location is determined by constructing an objective function. The embodiment of the invention can enhance the uniqueness of the positioning solution, reduce the dependence of the positioning accuracy on the velocity model, the given seismic source depth and the seismic source mechanism, and improve the positioning efficiency, the positioning accuracy and the positioning precision.
Description
Technical Field
The invention belongs to the technical field of earthquakes, and particularly relates to an earthquake positioning method, an earthquake positioning device, terminal equipment and a computer readable storage medium.
Background
Seismic positioning is one of the most classical and basic problems in the field of seismology, and an accurate seismic positioning result has important significance for researching a seismic activity structure, an internal structure of the earth, a seismic source fracture process, a geometrical structure and the like and is also important for disaster relief work after an earthquake.
The traditional seismic positioning method searches a space point which is most matched with the observed travel time or waveform as an earthquake occurrence position by calculating the travel time or waveform of a given velocity model, and has strong dependence on the accuracy of the velocity model. In order to reduce the serious influence of the velocity model error on the positioning result, the seismic position may be constrained by using a noise surface wave green function extracted between the station pairs, which generally includes the following steps: selecting stations near the earthquake epicenter position, and calculating surface wave empirical Green functions between a near station and a far station far away from the earthquake position through earthquake background noise cross-correlation, wherein the surface wave empirical Green functions comprise Rayleigh wave empirical Green functions EGFRayAnd Love wave empirical Green function EGFLov(ii) a Then, an approximate one-dimensional velocity model is given, and a surface wave Green function between a near platform and a far platform is synthesized, wherein the surface wave Green function comprises GFRayAnd GFLovAnd at the same time, the depth of the seismic source is givenAnd a seismic source mechanism for calculating and synthesizing a theoretical seismic map; respectively calculating the travel time correction quantity of Rayleigh waves and Love waves through waveform cross correlation; and finally, introducing the travel time correction quantities of Rayleigh waves and Love waves by combining the seismic records actually recorded by the remote station and the synthesized theoretical seismic image, and determining the seismic epicenter position by performing least square iterative inversion on the two travel time correction quantities.
It can be seen that the existing method for positioning by using surface wave earthquake is carried out on the premise that the general position of the earthquake to be positioned is known in advance, the prior information is needed to determine which station nearby is used for time-lapse correction, and when a plurality of stations nearby are in the same distance with the earthquake, a plurality of earthquake epicenter positions can be obtained by using different stations, and the solution is not unique; and when no proper station exists near the earthquake occurrence position or the station is far away, the deviation between the travel time correction amount and the actually required correction amount is large. Especially, when the seismic source is deep, the deviation between the travel time correction amount between the near table and the far table on the earth surface and the actually required correction amount of the deep seismic source is large, and the positioning is inaccurate. In addition, waveform cross-correlation calculation is required, and the calculation efficiency is seriously influenced due to large operation amount of the waveform cross-correlation calculation. In summary, the existing seismic positioning method may not have a unique positioning solution, and depends on the given seismic source depth and seismic source mechanism, so that the positioning accuracy, precision and positioning efficiency are low.
Disclosure of Invention
In view of this, embodiments of the present invention provide a seismic positioning method, an apparatus, a terminal device, and a computer-readable storage medium, so as to solve the problems that a seismic source located by an existing seismic positioning method may not be unique, positioning accuracy depends on a given seismic source depth and accuracy of a seismic source mechanism, and positioning accuracy, and efficiency are low, and a real seismic source mechanism can also be inverted through waveform fitting.
A first aspect of an embodiment of the present invention provides a seismic positioning method, including:
calculating a surface wave empirical Green function between station pairs through the seismic background noise cross-correlation recorded by the stations; the station pair comprises a station belonging to a first type station and a station belonging to a second type station, the first type station is the station arranged in an earthquake region, the second type station is the station arranged in a target region with the distance to the earthquake region larger than a preset distance threshold value, and the earthquake region is provided with a dense array;
calculating first group travel time of surface waves in a surface wave empirical Green's function between the first type station and the second type station, and calculating second group travel time of seismic surface waves according to surface wave data recorded by the second type station;
according to the first group travel time, constructing a group travel time table of surface waves from the earthquake-generating area to the station of the target area through interpolation;
correcting the waveform of the surface wave empirical Green function between the station pairs into a seismic waveform of a double-couple type virtual seismic source, and constructing a depth travel time correction scale according to the group travel time of the surface waves of the virtual seismic source;
and constructing an objective function according to the first group travel time, the second group travel time, the depth travel time correction table and preset correction time, and determining the position of the seismic source according to the objective function.
Optionally, after the determining the position of the seismic source, further comprising:
and determining the depth of the seismic source and the mechanism of the seismic source according to the seismic waveform of the virtual seismic source.
Optionally, the constructing a depth travel time correction table according to the group travel time of the surface wave of the virtual seismic source includes:
calculating a third group travel time of the surface wave of the virtual seismic source;
calculating a difference between the first group travel time and the third group travel time;
and traversing the depth of the seismic source according to the difference value, and constructing the depth travel time correction table.
Optionally, said determining a source location according to said objective function comprises:
based on the objective function phi (X, Y, h) ═ min | | | TEQ-TEGFs+TDep+T0||2Minimizing the travel time residual of the surface wave group, and determining the seismic source position through a grid search algorithm; wherein, TEGFSIs the first group travel time, TEQIs the second group travel time, TDepFor the depth travel time correction scale, T0The preset correction time amount is the amount of time that needs to be corrected based on the time function of the seismic source and the existence of the seismic moment.
Optionally, the calculating a first group travel time of a surface wave in a surface wave empirical green's function between the first type station and the second type station, and calculating a second group travel time of a seismic surface wave according to surface wave data recorded by the second type station includes:
calculating first group travel time of the surface wave in the empirical green function of the surface wave between the first type station and the second type station by carrying out narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave;
and calculating a second group travel time of the seismic surface wave according to the surface wave data recorded by the second type station by carrying out narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave.
A second aspect of an embodiment of the present invention provides an earthquake locating apparatus, including:
the calculation module is used for calculating the surface wave empirical Green function between the station pairs through the seismic background noise cross-correlation recorded by the station; the station pair comprises a station belonging to a first type station and a station belonging to a second type station, the first type station is the station arranged in an earthquake region, the second type station is the station arranged in a target region with the distance to the earthquake region larger than a preset distance threshold value, and the earthquake region is provided with a dense array;
the travel time calculation module is used for calculating first group travel time of surface waves in a surface wave empirical Green function between the first type station and the second type station and calculating second group travel time of seismic surface waves according to surface wave data recorded by the second type station;
the travel time table calculation module is used for constructing a group travel time table of surface waves from the earthquake-generating area to the station of the target area through interpolation according to the first group travel time;
the construction module is used for correcting the waveform of the surface wave empirical Green function between the station pairs into a seismic waveform of a double couple type virtual seismic source and constructing a depth travel time correction scale according to the group travel time of the surface wave of the virtual seismic source;
and the first determining module is used for constructing an objective function according to the first group travel time, the second group travel time, the depth travel time correction table and preset correction time, and determining the position of the seismic source according to the objective function.
Optionally, the method further comprises:
and the second determination module is used for determining the depth of the seismic source and the mechanism of the seismic source according to the seismic waveform of the virtual seismic source.
Optionally, the building module comprises:
the computing unit is used for computing the third group travel time of the surface wave of the virtual seismic source;
a difference calculation unit for calculating a difference between the first group travel time and the third group travel time;
and the construction unit is used for traversing the depth of the seismic source according to the difference value and constructing the depth travel time correction scale.
A third aspect of 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, where the processor implements the steps of the method according to any one of the above first aspects when executing the computer program.
A fourth aspect of embodiments of the present invention provides a computer-readable storage medium, in which a computer program is stored, which, when executed by a processor, performs the steps of the method according to any one of the first aspect.
Compared with the prior art, the embodiment of the invention has the following beneficial effects:
according to the embodiment of the invention, the dense array is arranged in the earthquake generating area, so that more data are provided for positioning, the uniqueness of understanding is enhanced, and the positioning accuracy and the positioning precision are improved; the deviation brought by the velocity model is corrected by using the seismic background noise, and the waveform of the surface wave empirical Green's function is corrected into a virtual seismic source, so that the dependence of the positioning accuracy on the given seismic source depth and seismic source mechanism is reduced, and the accuracy and the positioning accuracy are improved; the waveform cross-correlation calculation is not needed, so that the calculation amount is reduced, and the positioning efficiency is improved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the embodiments or the prior art descriptions will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without inventive exercise.
Fig. 1 is a schematic flow chart of a seismic positioning method according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a station distribution according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of seismic positioning effects provided by embodiments of the present invention;
FIG. 4 is another schematic block flow diagram of a seismic positioning method provided by an embodiment of the invention;
FIG. 5 is a block diagram illustrating a seismic locating device according to 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 Description
In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular system structures, techniques, etc. in order to provide a thorough understanding of the embodiments of the invention. It will be apparent, however, to one skilled in the art that the present invention may be practiced in other embodiments that depart from 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 explain the technical means of the present invention, the following description will be given by way of specific examples.
Referring to fig. 1, a schematic flow chart of a seismic positioning method according to an embodiment of the present invention is shown, where the method includes the following steps:
s101, calculating a surface wave empirical Green function between station pairs through seismic background noise cross correlation recorded by the stations; wherein, the station is to including the station that belongs to first kind station and the station that belongs to second kind station, and first kind station is for setting up in the station of sending out the earthquake district, and second kind station is for setting up in the station that is greater than the target area's of preset distance threshold value apart from sending out the earthquake district, sends out the earthquake district and is provided with intensive platform array.
It is understood that the seismic background noise refers to the seismic wake wave signal and the weak seismic signal recorded on the station, which are distinguished from the seismic signal with a larger magnitude, and are called background noise because the signal is weak and has an amplitude equivalent to that of the noise.
The earthquake region is an area where an earthquake is active, and is close to the epicenter position. The first type of station located within the origin area may be referred to as a near station. A first type station is arranged in the earthquake area, and the first type stations are in a dense array. It can be seen that by using a dense array of seismic zones to obtain more positioning data, the uniqueness of the positioning solution is enhanced.
The second type of station may be referred to as a far station, which is a station far from the earthquake region, the far station is disposed in the target region, and the distance from the target region to the earthquake region is greater than a preset distance threshold, which may be set according to an actual application scenario.
The station pair includes a near station and a far station, and surface wave Empirical Green's Functions (EGFs) between the station pair are calculated, that is, surface wave Empirical Green's Functions (EGFs) between the near station and the far station are calculated. Wherein, the calculation process can be calculated by the following formula (1).
The left term of the formula (1) isThe derivative of the background noise cross-correlation with time t, the right termSurface waves EGFs between station A and station B, where CABIs a cross-correlation function between station a and station B.
To better describe the station distribution, the following description will be made with reference to the station distribution diagram provided in the embodiment of the present invention shown in fig. 2.
As shown in fig. 2, in the plane formed by the X axis and the Y axis, "○" indicates the center position, "□" indicates the far stage, "△" indicates the near stage, a plurality of near stages, i.e., a dense stage matrix, are collectively provided in the vicinity of the center position as an oscillation region, and a plurality of far stages are scattered in a staggered manner at positions farther from the center position.
It is to be understood that fig. 2 is an exemplary distribution diagram, and that the distribution of stations may be in other forms, and is not limited thereto.
Step S102, calculating first group travel time of surface waves in the empirical Green' S function of the surface waves between the first type station and the second type station, and calculating second group travel time of seismic surface waves according to surface wave data recorded by the second type station.
It is understood that seismic travel time refers to the time required for seismic waves to reach an observation point from a seismic source. And the group travel time of the surface wave refers to the wave envelope of the surface wave or the propagation time of energy.
In specific application, the corresponding surface wave group travel time can be calculated by carrying out narrow-band frequency domain Gaussian filtering and Hilbert transform on the surface wave waveform.
In an embodiment, the group travel time may be obtained by performing narrow-band frequency domain gaussian filtering and hilbert transform on the waveform, so optionally, this step, that is, calculating the first group travel time of the surface wave in the empirical green's function of the surface wave between the first type station and the second type station, and calculating the second group travel time of the seismic surface wave according to the surface wave data recorded by the second type station may specifically be: calculating first group travel time of the surface wave in the empirical Green's function of the surface wave between the first type station and the second type station by carrying out narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave; and calculating a second group travel time of the seismic surface wave according to the surface wave data recorded by the second type station by carrying out narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave.
The narrow-band frequency domain gaussian filter formula is formula (2) to formula (4), and the hilbert transform is as shown in formula 3 below.
Wherein,
w is the original waveform data prior to transformation,is a Fourier transform of w, H is a Gaussian filter factor,and performing Hilbert transform on the filtered waveform.
And S103, constructing a group travel time table of surface waves between the station from the vibration generating area to the target area through interpolation according to the first group travel time.
It is understood that the chronology refers to the schedule of seismic wave propagation over different epicenter distances. Although the earthquake area is provided with the dense array, the earthquake area is divided into finer grids for interpolation, and a travel time table is constructed, so that the positioning precision and the positioning accuracy can be further improved.
And S104, correcting the waveform of the surface wave empirical Green function between the station pairs into the seismic waveform of the double-couple type virtual seismic source, and constructing a depth travel time correction scale according to the group travel time of the surface waves of the virtual seismic source.
Specifically, the waveform of the empirical green's function of the surface wave between the station pairs, that is, the waveform of the surface wave EGFs obtained by the seismic background noise is corrected to the seismic waveform of a virtual source of a Double-couple (Double-couple) type having a certain depth. And then, calculating the difference value between the group travel time of the surface wave of the virtual seismic source and the EGFs surface wave group travel time, namely the first group travel time, traversing the seismic source depth, and obtaining a surface wave group travel time correction table required to be corrected and caused by the real seismic source depth.
In a specific application, the frequency domain correction formula from the EGFs waveform to the virtual seismic source waveform is as follows:
wherein,andrespectively, the frequency domain representation of the corrected three-component waveform, r1，r2And l1Characteristic function values of Rayleigh waves and Love waves respectively, M is a seismic source mechanism solution of the earthquake,empirical green's functions for the surface waves corresponding to the different components.
Through the formulas (6) to (8), the surface wave EGFs wave obtained through the seismic background noise can be corrected into the seismic waveform of the virtual seismic source.
And S105, constructing an objective function according to the first group travel time, the second group travel time and the depth travel time calibration table and preset calibration time, and determining the position of the seismic source according to the objective function.
Specifically, an objective function is constructed by traversing the depth of the seismic source, the travel time residual of the surface wave group is minimized, and then the position of the seismic source can be determined through a grid search algorithm.
In one embodiment, optionally, the specific process of determining the position of the seismic source according to the objective function may include: based on the objective function phi (X, Y, h) ═ min | | | TEQ-TEGFs+TDep+T0||2Minimizing the travel time residual of the surface wave group, and determining the position of a seismic source through a grid search algorithm; wherein, TEGFSIs the first group, TEQIs the second group of travel time, TDepFor depth travel-time correction scales, T0For a preset correction time, i.e. T0The amount of time that needs to be corrected due to the source time function and the presence of the origin time.
In order to better describe the positioning result, the following description will be made with reference to the schematic diagram of the seismic positioning effect provided by the embodiment of the invention shown in fig. 3. As shown in fig. 3, the closer the epicenter position is, the darker the color is, and the deepest position is the true epicenter position. The dots in the origin region are the epicenter positions determined by the method, and can be seen to be completely overlapped with the real epicenter positions. The positioning effect schematic diagram is an effect diagram obtained by using the earthquake positioning method provided by the embodiment of the invention, and can be seen that the epicenter position is accurately positioned.
In the embodiment, the dense array is arranged in the earthquake-generating area, so that more data are provided for positioning, the uniqueness of understanding is enhanced, and the positioning accuracy and the positioning precision are improved; the deviation brought by the velocity model is corrected by using the seismic background noise, and the waveform of the surface wave empirical Green's function is corrected into a virtual seismic source, so that the dependence of the positioning accuracy on the velocity model, the given seismic source depth and a seismic source mechanism is reduced, and the accuracy and the positioning accuracy are improved; the waveform cross-correlation calculation is not needed, so that the calculation amount is reduced, and the positioning efficiency is improved.
Example two
Referring to fig. 4, another schematic flow chart of a seismic positioning method according to an embodiment of the present invention is shown, where the method may include the following steps:
step S401, through the seismic background noise cross-correlation recorded by the station, the empirical Green function of the surface wave between the station pairs is calculated.
Step S402, calculating first group travel time of surface waves in the empirical Green' S function of the surface waves between the first type station and the second type station, and calculating second group travel time of seismic surface waves according to surface wave data recorded by the second type station.
And S403, constructing a group travel time table of surface waves between the stations from the origin area to the target area through interpolation according to the first group travel time.
It should be noted that the steps S301 to S303 are the same as the steps S201 to S203 of the first embodiment, and for related description, reference is made to the above corresponding contents, which are not repeated herein.
And S404, correcting the waveform of the surface wave empirical Green function between the station pairs into the seismic waveform of the double-couple type virtual seismic source.
And S405, calculating a third group travel time of the surface wave group of the virtual seismic source, calculating a difference value between the first group travel time and the third group travel time, traversing the seismic source depth according to the difference value, and constructing a depth travel time correction table.
Step S406, constructing an objective function according to the first group travel time, the second group travel time, the depth travel time correction table and the preset correction time, and determining the position of the seismic source according to the objective function.
Step S407, according to the seismic waveform of the virtual seismic source, the seismic source depth and the seismic source mechanism are determined.
Specifically, the real depth of the seismic source and the mechanism of the seismic source can be determined simultaneously by a traversal algorithm and waveform fitting according to the seismic waveform of the virtual seismic source.
In the embodiment, the dense array is arranged in the earthquake-generating area, so that more data are provided for positioning, the uniqueness of understanding is enhanced, and the positioning accuracy and the positioning precision are improved; the deviation brought by the velocity model is corrected by using the seismic background noise, and the waveform of the surface wave empirical Green's function is corrected into a virtual seismic source, so that the dependence of the positioning accuracy on the velocity model, the given seismic source depth and a seismic source mechanism is reduced, and the accuracy and the positioning accuracy are improved; the waveform cross-correlation calculation is not needed, so that the calculation amount is reduced, and the positioning efficiency is improved.
It should be understood that, the sequence numbers of the steps in the foregoing embodiments do not imply an execution sequence, and the execution sequence of each process should be determined by its function and inherent logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.
EXAMPLE III
Referring to fig. 5, a schematic block diagram of a seismic positioning apparatus according to an embodiment of the present invention is provided, where the apparatus may include:
a calculating module 51, configured to calculate a surface wave empirical green function between pairs of stations through seismic background noise cross-correlation recorded by the stations; the station pair comprises a station belonging to a first type station and a station belonging to a second type station, the first type station is a station arranged in an earthquake generating area, the second type station is a station arranged in a target area, the distance between the target area and the earthquake generating area is greater than a preset distance threshold value, and the earthquake generating area is provided with a dense array;
the travel time calculation module 52 is configured to calculate a first group travel time of a surface wave in a surface wave empirical green's function between the first type station and the second type station, and calculate a second group travel time of a seismic surface wave according to surface wave data recorded by the second type station;
a travel time table calculation module 53, configured to construct a group travel time table of a surface wave of a station from the origin area to the target area by interpolation according to the first group travel time;
a building module 54, configured to correct a waveform of a surface wave empirical green's function between the station pairs to a seismic waveform of a double couple type virtual seismic source, and build a depth travel time correction scale according to a group travel time of a surface wave of the virtual seismic source;
and the first determining 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 calibration table and the preset calibration time, and determine the seismic source position according to the objective function.
In one possible implementation, the apparatus may further include:
and the second determination module is used for determining the depth of the seismic source and the mechanism of the seismic source according to the seismic waveform of the virtual seismic source.
In a possible implementation, the building module may include:
the computing unit is used for computing the third group travel time of the surface wave of the virtual seismic source;
a difference calculation unit for calculating a difference between the first group travel time and the third group travel time;
and the construction unit is used for traversing the depth of the seismic source according to the difference value and constructing a depth travel time correction scale.
In a possible implementation, the first determining module may include:
a determination unit for determining a target function phi (X, Y, h) ═ min | | | TEQ-TEGFs+TDep+T0||2Minimizing the travel time residual of the surface wave group, and determining the position of a seismic source through a grid search algorithm; wherein, TEGFSIs the first group, TEQIs the second group of travel time, TDepFor correcting depth travel timeWatch, T0Is a preset amount of correction time.
In a possible implementation, the travel time calculation module may include:
the first travel time calculation unit is used for calculating first group travel time of the surface wave in the empirical green's function of the surface wave between the first type station and the second type station by carrying out narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave;
and the second travel time calculating unit is used for calculating the second group travel time of the seismic surface wave according to the surface wave data recorded by the second type station by carrying out narrow-band frequency domain high-pass Gaussian filtering and Hilbert conversion on the waveform of the surface wave.
In the embodiment, the dense array is arranged in the earthquake-generating area, so that more data are provided for positioning, the uniqueness of understanding is enhanced, and the positioning accuracy and the positioning precision are improved; the deviation brought by the velocity model is corrected by using the seismic background noise, and the waveform of the surface wave empirical Green's function is corrected into a virtual seismic source, so that the dependence of the positioning accuracy on the velocity model, the given seismic source depth and a seismic source mechanism is reduced, and the accuracy and the positioning accuracy are improved; the waveform cross-correlation calculation is not needed, so that the calculation amount is reduced, and the positioning efficiency is improved.
Example four
Fig. 6 is a schematic diagram of a terminal device according to an embodiment of the present invention. As shown in fig. 6, the terminal device 6 of this embodiment includes: a processor 60, a memory 61 and a computer program 62 stored in said memory 61 and executable on said processor 60. The processor 60, when executing the computer program 62, implements the steps in the various seismic positioning method embodiments described above, such as steps S101-S105 shown in fig. 1. Alternatively, the processor 60, when executing the computer program 62, implements the functions of the modules or units in the above-described device embodiments, such as the functions of the modules 51 to 55 shown in fig. 5.
Illustratively, the computer program 62 may be divided into one or more modules or units, which are stored in the memory 61 and executed by the processor 60 to accomplish the present invention. The one or more modules or units may be a series of computer program instruction segments capable of performing specific functions, which are used to describe the execution of the computer program 62 in the terminal device 6. For example, the computer program 62 may be divided into a computing module, a travel time table computing module, a building module, and a first determining module, and each module has the following specific functions:
the calculation module is used for calculating the surface wave empirical Green function between the station pairs through the seismic background noise cross-correlation recorded by the station; the station pair comprises a station belonging to a first type station and a station belonging to a second type station, the first type station is the station arranged in an earthquake region, the second type station is the station arranged in a target region with the distance to the earthquake region larger than a preset distance threshold value, and the earthquake region is provided with a dense array; the travel time calculation module is used for calculating first group travel time of surface waves in a surface wave empirical Green function between the first type station and the second type station and calculating second group travel time of seismic surface waves according to surface wave data recorded by the second type station; the travel time table calculation module is used for constructing a group travel time table of surface waves from the earthquake-generating area to the station of the target area through interpolation according to the first group travel time; the construction module is used for correcting the waveform of the surface wave empirical Green function between the station pairs into a seismic waveform of a double couple type virtual seismic source and constructing a depth travel time correction scale according to the group travel time of the surface wave of the virtual seismic source; and the first determining module is used for constructing an objective function according to the first group travel time, the second group travel time, the depth travel time correction table and preset correction time, and determining the position of the seismic source according to the objective function.
The terminal device 6 may be a desktop computer, a notebook, a palm computer, a cloud server, or other computing devices. The terminal device may include, but is not limited to, a processor 60, a memory 61. Those skilled in the art will appreciate that fig. 6 is merely an example of a terminal device 6 and does not constitute a limitation of terminal device 6 and may include more or less components than those shown, or some components in combination, or different components, for example, the terminal device may also include input output devices, network access devices, buses, etc.
The Processor 60 may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. 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, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), and the like, which are provided on the terminal device 6. Further, the memory 61 may also include both an internal storage unit and an external storage device of the terminal device 6. The memory 61 is used for storing the computer program and other programs and data required by the terminal device. The memory 61 may also be used to temporarily store data that has been output or is to be output.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-mentioned division of the functional units and modules is illustrated, and in practical applications, the above-mentioned function distribution may be performed by different functional units and modules according to needs, that is, the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-mentioned functions. Each functional unit and module in the embodiments may be integrated in one processing unit, or each unit may exist alone physically, or two or more units are integrated in one unit, and the integrated unit may be implemented in a form of hardware, or in a form of software functional unit. In addition, specific names of the functional units and modules are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present application. The specific working processes of the units and modules in the system may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the above embodiments, the descriptions of the respective embodiments have respective emphasis, and reference may be made to the related descriptions of other embodiments for parts that are not described or illustrated in a certain embodiment.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
In the embodiments provided in the present invention, it should be understood that the disclosed apparatus, terminal device and method may be implemented in other ways. For example, the above-described embodiments of the apparatus and the terminal device are merely illustrative, and for example, the division of the module or the unit is only one logical function division, and there may be another division in actual implementation, for example, a plurality of units or components may be combined or may be integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated modules or units, if implemented in the form of software functional units and sold or used as separate products, may be stored in a computer readable storage medium. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method embodiments may be implemented. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution medium, and the like. It should be noted that the computer readable medium may contain content that is subject to appropriate increase or decrease as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, computer readable media does not include electrical carrier signals and telecommunications signals as is required by legislation and patent practice.
The above-mentioned embodiments are only used for illustrating the technical solutions of the present invention, and not for limiting the same; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications and substitutions do not substantially depart from the spirit and scope of the embodiments of the present invention, and are intended to be included within the scope of the present invention.
Claims (10)
1. A seismic positioning method, comprising:
calculating a surface wave empirical Green function between station pairs through the seismic background noise cross-correlation recorded by the stations; the station pair comprises a station belonging to a first type station and a station belonging to a second type station, the first type station is the station arranged in an earthquake region, the second type station is the station arranged in a target region with the distance to the earthquake region larger than a preset distance threshold value, and the earthquake region is provided with a dense array;
calculating first group travel time of surface waves in a surface wave empirical Green's function between the first type station and the second type station, and calculating second group travel time of seismic surface waves according to surface wave data recorded by the second type station;
according to the first group travel time, constructing a group travel time table of surface waves from the earthquake-generating area to the station of the target area through interpolation;
correcting the waveform of the surface wave empirical Green function between the station pairs into a seismic waveform of a double-couple type virtual seismic source, and constructing a depth travel time correction scale according to the group travel time of the surface waves of the virtual seismic source;
and constructing an objective function according to the first group travel time, the second group travel time, the depth travel time correction table and preset correction time, and determining the position of the seismic source according to the objective function.
2. The seismic positioning method of claim 1, further comprising, after said determining the source location:
and determining the depth of the seismic source and the mechanism of the seismic source according to the seismic waveform of the virtual seismic source.
3. The seismic positioning method of claim 1, wherein constructing a depth travel time correction scale from the group travel times of the surface waves of the virtual seismic source comprises:
calculating a third group travel time of the surface wave of the virtual seismic source;
calculating a difference between the first group travel time and the third group travel time;
and traversing the depth of the seismic source according to the difference value, and constructing the depth travel time correction table.
4. A seismic positioning method according to any of claims 1 to 3, wherein said determining a source location from said objective function comprises:
based on the objective function phi (X, Y, h) ═ min | | | TEQ-TEGFs+TDep+T0||2Minimizing the travel time residual of the surface wave group, and determining the seismic source position through a grid search algorithm; wherein, TEGFSIs the first group travel time, TEQIs the second group travel time, TDepFor the depth travel time correction scale, T0The preset correction time amount is the amount of time that needs to be corrected based on the time function of the seismic source and the existence of the seismic moment.
5. The seismic positioning method of claim 4, wherein said calculating a first group travel time of a surface wave in a surface wave empirical green's function between said first type of station and said second type of station, and calculating a second group travel time of a seismic surface wave from surface wave data recorded by said second type of station comprises:
calculating first group travel time of the surface wave in the empirical green function of the surface wave between the first type station and the second type station by carrying out narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave;
and calculating a second group travel time of the seismic surface wave according to the surface wave data recorded by the second type station by carrying out narrow-band frequency domain high-pass Gaussian filtering and Hilbert transform on the waveform of the surface wave.
6. A seismic locating device, comprising:
the calculation module is used for calculating the surface wave empirical Green function between the station pairs through the seismic background noise cross-correlation recorded by the station; the station pair comprises a station belonging to a first type station and a station belonging to a second type station, the first type station is the station arranged in an earthquake region, the second type station is the station arranged in a target region with the distance to the earthquake region larger than a preset distance threshold value, and the earthquake region is provided with a dense array;
the travel time calculation module is used for calculating first group travel time of surface waves in a surface wave empirical Green function between the first type station and the second type station and calculating second group travel time of seismic surface waves according to surface wave data recorded by the second type station;
the travel time table calculation module is used for constructing a group travel time table of surface waves from the earthquake-generating area to the station of the target area through interpolation according to the first group travel time;
the construction module is used for correcting the waveform of the surface wave empirical Green function between the station pairs into a seismic waveform of a double couple type virtual seismic source and constructing a depth travel time correction scale according to the group travel time of the surface wave of the virtual seismic source;
and the first determining module is used for constructing an objective function according to the first group travel time, the second group travel time, the depth travel time correction table and preset correction time, and determining the position of the seismic source according to the objective function.
7. The seismic positioning apparatus of claim 6, further comprising:
and the second determination module is used for determining the depth of the seismic source and the mechanism of the seismic source according to the seismic waveform of the virtual seismic source.
8. The seismic positioning apparatus of claim 6, wherein the building module comprises:
the computing unit is used for computing the third group travel time of the surface wave of the virtual seismic source;
a difference calculation unit for calculating a difference between the first group travel time and the third group travel time;
and the construction unit is used for traversing the depth of the seismic source according to the difference value and constructing the depth travel time correction scale.
9. A terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of the method according to any one of claims 1 to 5 when executing the computer program.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program which, when being executed by a processor, carries out the steps of the method according to any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811071892.9A CN110058299B (en) | 2018-09-14 | 2018-09-14 | Earthquake positioning method and device and terminal equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811071892.9A CN110058299B (en) | 2018-09-14 | 2018-09-14 | Earthquake positioning method and device and terminal equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110058299A true CN110058299A (en) | 2019-07-26 |
| CN110058299B CN110058299B (en) | 2020-09-11 |
Family
ID=67315792
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811071892.9A Active CN110058299B (en) | 2018-09-14 | 2018-09-14 | Earthquake positioning method and device and terminal equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110058299B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110988987A (en) * | 2019-11-19 | 2020-04-10 | 南方科技大学 | Seismic source positioning method under single-well observation system and related device and equipment |
| CN111257928A (en) * | 2020-02-12 | 2020-06-09 | 中国科学技术大学 | Seismic source positioning method and positioning system |
| CN112255687A (en) * | 2020-10-26 | 2021-01-22 | 中煤科工集团西安研究院有限公司 | Seismic source function reconstruction method and device of seismic shearer during mining of stope face |
| CN112861721A (en) * | 2021-02-09 | 2021-05-28 | 南方科技大学 | Method and device for automatically extracting background noise dispersion curve |
| CN113221395A (en) * | 2021-03-16 | 2021-08-06 | 禁核试北京国家数据中心 | Construction method and application of seismic travel time parameterized grid model based on layered medium |
| CN113238280A (en) * | 2021-06-24 | 2021-08-10 | 成都理工大学 | Green function-based earthquake monitoring method |
| CN113359186A (en) * | 2021-06-01 | 2021-09-07 | 西北核技术研究所 | Observation signal amplitude measuring method based on natural seismic source radiation intensity correction |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000241558A (en) * | 1999-02-23 | 2000-09-08 | Sekisui Chem Co Ltd | Method for ground vibration measurement and method for ground structure estimation |
| CN102183786A (en) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | Fidelity amplitude gaussian beam pre-stack depth migration method under double complex conditions |
| CN102395902A (en) * | 2009-04-16 | 2012-03-28 | 兰德马克图形公司 | Seismic imaging systems and methods employing a fast target-oriented illumination calculation |
| CN104678435A (en) * | 2014-10-27 | 2015-06-03 | 李欣欣 | Method for extracting Rayleigh surface wave frequency dispersion curve |
| WO2015134379A1 (en) * | 2014-03-03 | 2015-09-11 | Westerngeco Llc | Modeling and filtering coherent noise in seismic surveys |
| CN106597536A (en) * | 2016-12-07 | 2017-04-26 | 中国石油天然气股份有限公司 | Processing method and device for seismic imaging |
| EP3232234A1 (en) * | 2016-04-13 | 2017-10-18 | CGG Services SAS | Method and apparatus for estimating surface wave coda using time-reversal experiments |
-
2018
- 2018-09-14 CN CN201811071892.9A patent/CN110058299B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000241558A (en) * | 1999-02-23 | 2000-09-08 | Sekisui Chem Co Ltd | Method for ground vibration measurement and method for ground structure estimation |
| CN102395902A (en) * | 2009-04-16 | 2012-03-28 | 兰德马克图形公司 | Seismic imaging systems and methods employing a fast target-oriented illumination calculation |
| CN102183786A (en) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | Fidelity amplitude gaussian beam pre-stack depth migration method under double complex conditions |
| WO2015134379A1 (en) * | 2014-03-03 | 2015-09-11 | Westerngeco Llc | Modeling and filtering coherent noise in seismic surveys |
| CN104678435A (en) * | 2014-10-27 | 2015-06-03 | 李欣欣 | Method for extracting Rayleigh surface wave frequency dispersion curve |
| EP3232234A1 (en) * | 2016-04-13 | 2017-10-18 | CGG Services SAS | Method and apparatus for estimating surface wave coda using time-reversal experiments |
| CN106597536A (en) * | 2016-12-07 | 2017-04-26 | 中国石油天然气股份有限公司 | Processing method and device for seismic imaging |
Non-Patent Citations (2)
| Title |
|---|
| JIE ZHANG 等: ""搜索引擎方法在实时地震监测中的应用"", 《国际地震动态》 * |
| 王椿镛 等: ""南北构造带岩石圈结构与地震的研究"", 《地球物理学报》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110988987A (en) * | 2019-11-19 | 2020-04-10 | 南方科技大学 | Seismic source positioning method under single-well observation system and related device and equipment |
| CN111257928B (en) * | 2020-02-12 | 2021-08-13 | 中国科学技术大学 | Seismic source positioning method and positioning system |
| CN111257928A (en) * | 2020-02-12 | 2020-06-09 | 中国科学技术大学 | Seismic source positioning method and positioning system |
| CN112255687A (en) * | 2020-10-26 | 2021-01-22 | 中煤科工集团西安研究院有限公司 | Seismic source function reconstruction method and device of seismic shearer during mining of stope face |
| CN112255687B (en) * | 2020-10-26 | 2024-01-19 | 中煤科工集团西安研究院有限公司 | Method and device for reconstructing seismic source function of stope face along with mining earthquake coal mining machine |
| CN112861721A (en) * | 2021-02-09 | 2021-05-28 | 南方科技大学 | Method and device for automatically extracting background noise dispersion curve |
| CN112861721B (en) * | 2021-02-09 | 2024-05-07 | 南方科技大学 | Method and device for automatically extracting background noise dispersion curve |
| CN113221395B (en) * | 2021-03-16 | 2023-10-20 | 禁核试北京国家数据中心 | Method for constructing seismic travel time parameterized grid model based on layered medium and application |
| CN113221395A (en) * | 2021-03-16 | 2021-08-06 | 禁核试北京国家数据中心 | Construction method and application of seismic travel time parameterized grid model based on layered medium |
| CN113359186A (en) * | 2021-06-01 | 2021-09-07 | 西北核技术研究所 | Observation signal amplitude measuring method based on natural seismic source radiation intensity correction |
| CN113359186B (en) * | 2021-06-01 | 2023-07-21 | 西北核技术研究所 | Observed signal amplitude measurement method based on natural seismic source radiation intensity correction |
| CN113238280A (en) * | 2021-06-24 | 2021-08-10 | 成都理工大学 | Green function-based earthquake monitoring method |
| CN113238280B (en) * | 2021-06-24 | 2023-02-24 | 成都理工大学 | Green function-based earthquake monitoring method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110058299B (en) | 2020-09-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110058299B (en) | Earthquake positioning method and device and terminal equipment | |
| CN111164462B (en) | Artificial source surface wave exploration method, surface wave exploration device and terminal equipment | |
| CN109451532B (en) | Method and device for checking position of base station | |
| CN112685928B (en) | Noise prediction method and system based on three-phase reactor sound source model | |
| CN112288254B (en) | Construction quality detection method, terminal and storage medium | |
| US20110010098A1 (en) | Method of pre-stack two-dimension-like transformation of three-dimensional seismic record | |
| CN109828302B (en) | Seismic source positioning method and device based on multiple vibration sensors | |
| WO2019242045A9 (en) | Method for calculating virtual source two-dimensional wavefront construction seismic wave travel time | |
| CN105759311A (en) | Near-real time earthquake source position positioning method | |
| CN103686995A (en) | Non-line-of-sight environment positioning method and device | |
| CN111460071A (en) | Deflection method, device and equipment of high-precision map and readable storage medium | |
| CN112393730A (en) | Magnetic beacon positioning method and system based on magnetic induction intensity and characteristic vector | |
| CN111624658B (en) | Depth domain imaging simulation method and system | |
| CN116415652A (en) | Data generation method and device, readable storage medium and terminal equipment | |
| CN117095002B (en) | Hub defect detection method and device and storage medium | |
| CN107870361B (en) | Seismic inflection wave tomography method and device and terminal equipment | |
| CN106597535A (en) | Method of improving elastic wave reverse time migration offset computation rate and space resolution | |
| CN101840001B (en) | Acquiring method and device of geological structure three-dimensional imaging data | |
| CN106054252A (en) | Pre-stack time migration method and device | |
| CN108052754B (en) | Analysis method of mechanical response of tunnel lining and terminal equipment | |
| CN114201875A (en) | Method for determining multi-sound-source noise equivalent model of transformer, terminal and storage medium | |
| CN105572741A (en) | Method for calculating 3D high-frequency static correction value | |
| CN116165706A (en) | Underground structure imaging method based on background noise surface wave double-wave beam focusing | |
| CN115239520A (en) | Method and device for eliminating errors of model component, electronic device and storage medium | |
| CN115270579A (en) | Second-order acoustic wave equation finite difference numerical simulation parameter selection method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |