Multi-target impact point positioning method and device, computer equipment and storage medium

Abstract

The invention discloses a multi-target impact point positioning method and device, computer equipment and a storage medium, and relates to the technical field of seismic source positioning. The method comprises the following steps: determining a plurality of model point locations in a seismic source region according to a preset interval; obtaining vibration wave signals of vibration waves from each model point location to a plurality of detectors, and obtaining a vibration wave signal sequence corresponding to each detector under each model point location; obtaining an arrival time difference translation signal sequence corresponding to each detector under each model point location according to each vibration wave signal sequence; obtaining a plurality of brightness function values corresponding to each model point location according to each time difference translation signal sequence; and determining a plurality of target impact point coordinates according to the brightness function values. Therefore, the position of the impact point is found by using the shock wave signal to the time difference, and the problem of impact point positioning caused by double effects of disorder and aliasing of multi-target signals is solved.

Description

Technical Field

[0001] The present invention relates to the field of earthquake source positioning technology, and in particular to a multi-target impact point positioning method, device, computer equipment and storage medium. Background Art

[0002] The existing impact point positioning system uses a multi-station intersection positioning method to locate the impact point coordinates, but it is limited to single-target positioning of one impact point, and is powerless to locate multiple target impact points that occur in a short period of time. Summary of the invention

[0003] In view of this, the purpose of the present invention is to overcome the deficiencies in the prior art and to provide a multi-target impact point positioning method, device, computer equipment and storage medium for positioning multiple impact points by utilizing the detection time difference of a detector.

[0004] The present invention provides the following technical solutions:

[0005] In a first aspect, the present invention provides a multi-target impact point positioning method, comprising:

[0006] Determine multiple model points in the earthquake source area according to preset spacing;

[0007] Acquire the vibration wave signals from each of the model points to the multiple detectors, and obtain the vibration wave signal sequence corresponding to each of the detectors at each of the model points;

[0008] According to each of the vibration wave signal sequences, a time difference translation signal sequence corresponding to each of the detectors at each of the model points is obtained;

[0009] Obtaining a plurality of brightness function values ​​corresponding to each of the model points according to each of the time difference shift signal sequences;

[0010] A plurality of target impact point coordinates are determined according to the brightness function values.

[0011] In one embodiment, obtaining the time difference translation signal sequence corresponding to each detector at each model point according to each shock wave signal sequence includes:

[0012] For each of the model points, determining the relative arrival time difference of the vibration waves of each of the detectors according to the signal reception time corresponding to each of the vibration wave signal sequences;

[0013] Based on a preset time window length, each of the shock wave signal sequences is shifted according to the relative arrival time difference of the shock wave to obtain a plurality of arrival time difference shifted signal sequences.

[0014] In one embodiment, determining the relative arrival time difference of the vibration waves of each detector according to the signal reception time corresponding to each vibration wave signal sequence includes:

[0015] determining a reference object from each of said detectors;

[0016] The relative arrival time difference of the vibration wave of each detector relative to the reference object is determined according to the signal receiving time corresponding to each vibration wave signal sequence.

[0017] In one embodiment, determining the relative arrival time difference of the shock wave of each detector relative to the reference object according to the signal reception time corresponding to each shock wave signal sequence includes:

[0018] Normalizing the amplitude of each of the shock wave signal sequences to obtain a plurality of pre-processed signal sequences;

[0019] The relative arrival time difference of the vibration wave of each detector relative to the reference object is determined according to the signal receiving time corresponding to each pre-processed signal sequence.

[0020] In one embodiment, the determining of a plurality of target impact point coordinates according to each of the brightness function values ​​includes:

[0021] Determine the maximum value of the brightness function corresponding to each of the model points from the brightness function values;

[0022] Determining a plurality of objective function values ​​from the maximum values ​​of the brightness functions according to a preset source threshold;

[0023] The coordinates of each target impact point are determined according to each objective function value.

[0024] In one embodiment, the determining of a plurality of objective function values ​​from the maximum values ​​of the brightness functions according to a preset source threshold comprises:

[0025] Determine a brightness function peak group from each of the brightness function maximum values ​​according to the preset source threshold;

[0026] A plurality of real brightness peak values ​​are determined according to the brightness function peak value group, and the real brightness peak values ​​are used as the target function values.

[0027] In one embodiment, obtaining a plurality of brightness function values ​​corresponding to each of the model points according to each of the time difference shift signal sequences includes:

[0028] For each of the model points, based on multiple time points, the absolute value of the shock wave signal corresponding to each of the time points is obtained from each of the time difference translation signal sequences;

[0029] A plurality of brightness function values ​​corresponding to the model points are obtained according to the absolute values ​​of the shock wave signals.

[0030] In a second aspect, the present invention provides a multi-target impact point positioning device, comprising:

[0031] A determination module, used to determine a plurality of model points in the earthquake source area according to a preset spacing;

[0032] An acquisition module is used to acquire the vibration wave signals from each of the model points to a plurality of detectors, and obtain the vibration wave signal sequence corresponding to each of the detectors at each of the model points;

[0033] A translation module, used for obtaining a time difference translation signal sequence corresponding to each detector at each model point according to each vibration wave signal sequence;

[0034] A calculation module, used for obtaining a plurality of brightness function values ​​corresponding to each of the model points according to each of the time difference shift signal sequences;

[0035] A positioning module is used to determine the coordinates of multiple target impact points according to each of the brightness function values.

[0036] In a third aspect, the present invention proposes a computer device, comprising a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the multi-target impact point positioning method as described in the first aspect is implemented.

[0037] In a fourth aspect, the present invention proposes a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the multi-target impact point positioning method as described in the first aspect.

[0038] The multi-target impact point positioning method, device, computer equipment and storage medium disclosed in the present invention determine multiple model points according to a preset spacing in the epicenter area; obtain the shock wave signals of the shock wave from each of the model points to multiple detectors, and obtain the shock wave signal sequence corresponding to each of the detectors at each of the model points; obtain the time difference translation signal sequence corresponding to each of the detectors at each of the model points according to each of the shock wave signal sequences; obtain multiple brightness function values ​​corresponding to each of the model points according to each of the time difference translation signal sequences; and determine the coordinates of multiple target impact points according to each of the brightness function values. In this way, multiple model points are set in the epicenter area, and the original shock wave signal is translated using the time difference of the shock wave signal corresponding to each of the model points, so that the shock wave signal at the epicenter can be enhanced, and the impact point position corresponding to the enhanced shock wave signal is found using the brightness function, solving the impact point positioning problem caused by the dual effects of disorder and aliasing of multi-target signals. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] In order to more clearly illustrate the technical solution of the present invention, the following is a brief introduction to the drawings required for use in the embodiments. It should be understood that the following drawings only illustrate certain embodiments of the present invention and should not be regarded as limiting the scope of protection of the present invention. In each of the drawings, similar components are numbered similarly.

[0040] Figure 1 A schematic diagram of a process of locating the multi-target impact points proposed in this embodiment is shown;

[0041] Figure 2 A schematic diagram showing the model point layout proposed in this embodiment is shown;

[0042] Figure 3 A schematic diagram showing the impact point and the position of the detector proposed in this embodiment;

[0043] Figure 4 Another schematic diagram of the process of locating the multi-target impact points proposed in this embodiment is shown;

[0044] Figure 5 A schematic diagram of the shock wave signal sequence proposed in this embodiment is shown;

[0045] Figure 6 A schematic diagram of a time difference shift signal sequence proposed in this embodiment is shown;

[0046] Figure 7 A schematic diagram showing the shock wave signal sequence and the pre-processing signal sequence proposed in this embodiment is shown;

[0047] Figure 8 A schematic diagram showing multiple brightness function values ​​corresponding to a model point proposed in this embodiment;

[0048] Fig. 9 Another schematic diagram of the process of the multi-target impact point positioning method proposed in this embodiment is shown;

[0049] Fig.10 A schematic diagram showing a brightness function peak group proposed in this embodiment is shown;

[0050] Fig.11 A schematic diagram showing the multi-target positioning result proposed in this embodiment;

[0051] Fig.12 A schematic diagram of the structure of the multi-target impact point positioning device proposed in this embodiment is shown.

[0052] Description of the accompanying drawings:

[0053] 1200 - multi-target impact point positioning device; 1201 - determination module; 1202 - acquisition module; 1203 - translation module; 1204 - calculation module; 1205 - positioning module. DETAILED DESCRIPTION

[0054] The technical solutions in the embodiments of the present invention will be described clearly and completely below in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all the embodiments.

[0055] The components of the embodiments of the present invention generally described and shown in the drawings herein may be arranged and designed in a variety of different configurations. Therefore, the following detailed description of the embodiments of the present invention provided in the drawings is not intended to limit the scope of the claimed invention, but merely represents selected embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

[0056] Hereinafter, the terms "including", "having" and their cognates, which may be used in various embodiments of the present invention, are intended only to indicate specific features, numbers, steps, operations, elements, components or combinations of the foregoing items, and should not be understood as first excluding the existence of one or more other features, numbers, steps, operations, elements, components or combinations of the foregoing items or adding the possibility of one or more features, numbers, steps, operations, elements, components or combinations of the foregoing items.

[0057] Furthermore, the terms “first”, “second”, “third”, etc. are merely used for distinguishing descriptions and are not to be understood as indicating or implying relative importance.

[0058] Unless otherwise defined, all terms (including technical terms and scientific terms) used herein have the same meanings as those generally understood by those skilled in the art to which the various embodiments of the present invention belong. The terms (such as those defined in generally used dictionaries) will be interpreted as having the same meanings as the contextual meanings in the relevant technical field and will not be interpreted as having idealized meanings or overly formal meanings unless clearly defined in the various embodiments of the present invention.

[0059] Example 1

[0060] Before entering the countdown of the multi-target positioning test, the user runs the multi-target impact point positioning software; the data processing computer runs the multi-target impact point positioning software, starts initialization and self-test, and after the system is normal, imports the time difference model into the memory and waits for the target to be launched; according to the launch command, the user starts the signal acquisition function; the data processing computer sends the acquisition instruction to the wireless data transmission network communication interface; the wireless data transmission network transmits the acquisition instruction to each detector; each detector executes the signal acquisition instruction, and after filtering, digital / analog conversion, timing and other pre-processing of the signal, the data is packaged and sent to the network communication interface; after the wireless data transmission network receives the packaged data from the detector, it is transmitted to the data processing computer through the network; after receiving the packaged data, the data processing computer unpacks the data, and then displays the dynamic acquisition signal on the screen in real time, and the user can monitor the acquisition signal of each detector through the screen; according to the launch end command, after the user confirms that all shock wave signals have been received, the signal acquisition function is stopped, the positioning function is started, and the multi-target impact point positioning method is executed.

[0061] The embodiment of the present disclosure provides a multi-target impact point positioning method for positioning multiple impact points by utilizing the detection time difference of a detector.

[0062] See also Figure 1 The multi-target impact point positioning method includes steps S101 to S105, and each step is described in detail below.

[0063] Step S101, determining a plurality of model points in the earthquake source area according to a preset spacing.

[0064] In this embodiment, multiple model points are determined in the earthquake source area based on the preset spacing. For example, the 400m×400m earthquake source area is divided into 10201 (101×101) model points at a preset spacing of 0.5m. The model point layout diagram is as follows: Figure 2 shown.

[0065] For each pair of detectors, the 10201 model points have corresponding arrival time differences of seismic waves. This one-to-one correspondence can be regarded as an arrival time difference model with a preset spacing of 0.5 meters in the 400m×400m source area. The range of the preset spacing can be determined according to factors such as the wavelength of the seismic wave and the sensitivity of the detector.

[0066] Step S102, acquiring the shock wave signals from each of the model points to a plurality of detectors, and obtaining the shock wave signal sequence corresponding to each of the detectors at each of the model points.

[0067] In this embodiment, a plurality of detectors are arranged for detecting shock waves. The shock wave source in the source area generates shock waves. When each detector executes the signal acquisition instruction, each model point is used as the simulated propagation starting point of the shock wave, and each detector will record the shock wave signal transmitted from each model point. These shock wave signals include information such as the arrival time, waveform, and amplitude of the shock wave.

[0068] For each model point, each detector will receive and record the shock wave signals of different time scales transmitted from the model location, and form a shock wave signal sequence based on the shock wave signals. 10201 , and 8 detectors, then the model point A i The shock wave signal sequence under (i=1,2,……,10201) includes u1~u8, and the shock wave signal sequence u j (j=1, 2, ..., 8) is a sequence composed of the shock wave signals received by the jth detector.

[0069] See also Figure 3 The square in the figure is the earthquake source area, the circles inside the square represent the locations of the 10 impact points, and the diamonds outside the square represent the locations of the 8 detectors.

[0070] Step S103, obtaining a time difference translation signal sequence corresponding to each detector at each model point according to each vibration wave signal sequence.

[0071] In this embodiment, each model point of each pair of detectors has a corresponding vibration wave arrival time difference. Therefore, each vibration wave signal sequence can be translated according to the time information in the vibration wave signal sequence to obtain the corresponding arrival time difference translation signal sequence of each detector at each model point.

[0072] Since the shock wave signals of the detector at the source are superimposed and enhanced after being translated according to the time difference, the impact point can be located based on the translated signal sequence.

[0073] See also Figure 4 In a specific embodiment, step S103 includes steps S1031 to S1032, and each step is described in detail below.

[0074] Step S1031, for each of the model points, determining the relative arrival time difference of the vibration waves of each of the detectors according to the signal reception time corresponding to each of the vibration wave signal sequences.

[0075] In this embodiment, for each model point A i , the relative arrival time difference of the shock waves of each detector is determined according to the signal reception time corresponding to each shock wave signal sequence at the current model point.

[0076] In a specific embodiment, step S1031 includes: determining a reference object from each of the detectors; and determining a relative time difference of the vibration wave of each of the detectors relative to the reference object according to a signal reception time corresponding to each of the vibration wave signal sequences.

[0077] In this embodiment, a reference object is determined from each detector, and the relative arrival time difference of the shock wave of each other detector relative to the reference object is further determined according to the signal reception time corresponding to each shock wave signal sequence. Among them, the relative arrival time difference of the shock wave is obtained for each model point, and a unified reference object should be used.

[0078] Exemplarily, τ η1 and τ ηj They are respectively from model point A i The travel time of the shock wave η to detector 1 and detector j, with detector 1 as the reference object, the relative arrival time difference τ of the shock wave from detector j to detector 1 j1 is: j1 =τ ηj -τ η1 .

[0079] Step S1032, based on a preset time window length, each of the shock wave signal sequences is shifted according to the relative arrival time difference of the shock wave to obtain a plurality of arrival time difference shifted signal sequences.

[0080] In this embodiment, a signal time window is selected according to the measured waveform, and a preset time window length is set to m; further, each shock wave signal sequence is shifted according to the preset time window length m and the relative arrival time difference of the shock wave, and a plurality of corresponding arrival time difference shifted signal sequences are obtained.

[0081] For example, taking 8 detectors as an example, for model point A i , the time difference shift signal sequences are:

[0082] u n1 (T0+1),u n1 (T0+2),……,u n1 (T0+m);

[0083] u n2 (T0+1+τ 21 ),u n2 (T0+2+τ 21 ),……,u n2 (T0+m+τ 21 );

[0084] …

[0085] u n8 (T0+1+τ81 ),u n8 (T0+2+τ 81 ),……,u n8 (T0+m+τ 81 ) 。

[0086] Among them, T0 is the starting time of the preset time window, corresponding to the serial number of the starting time of the time window.

[0087] See also Figure 5 Due to the aliasing of different target signals, not every channel in the shock wave signal sequence can clearly distinguish 10 target signals. Taking detector 1 as the reference object, the remaining 7 detectors are shifted according to their respective arrival time differences. Figure 6 shown.

[0088] In a specific embodiment, the amplitude of each of the vibration wave signal sequences is normalized to obtain a plurality of pre-processed signal sequences; and the relative time difference of the vibration wave of each detector relative to the reference object is determined according to the signal reception time corresponding to each of the pre-processed signal sequences.

[0089] In this embodiment, since the distance from each impact point to each detector is different, the attenuation of the shock wave signal is also different, that is, the amplitude of the signal collected by the detector is also different.

[0090] For multi-target positioning, the influence caused by different distances must be eliminated to ensure that the signal from the same source has the same influence on the brightness function. Therefore, before the signal sequence is translated, the original shock wave signal sequence must be normalized. Figure 7 , (a) is the shock wave signal sequence, (b) is the preprocessed signal sequence.

[0091] Step S104, obtaining a plurality of brightness function values ​​corresponding to each of the model points according to each of the time difference shift signal sequences.

[0092] In this embodiment, a plurality of brightness function values ​​corresponding to each model point are calculated based on the signal data in each time difference shift signal sequence, and are used to analyze the positioning of a plurality of impact points.

[0093] In a specific embodiment, step S104 includes: for each of the model points, based on multiple time points, obtaining the absolute value of the shock wave signal corresponding to each of the time points from each of the time difference translation signal sequences; and obtaining multiple brightness function values ​​corresponding to the model point according to the absolute value of each of the shock wave signals.

[0094] In this embodiment, for each model point, the absolute value of the shock wave signal corresponding to different time points in the time difference shift signal sequence is obtained, and multiple brightness function values ​​corresponding to the model point are obtained according to the absolute value of each shock wave signal. Figure 8 shown.

[0095] For example, for each model point A i , taking 8 detectors as an example, the brightness function value calculation formula is: In the formula, N=8, i=(1,2,…,m).

[0096] Step S105, determining a plurality of target impact point coordinates according to the brightness function values.

[0097] In this embodiment, the point coordinates corresponding to the maximum value of the brightness function among all brightness function values ​​are the impact point coordinates of the positioning. When there are multiple impact points, it is necessary to screen all the maximum values ​​of the brightness function to determine the multiple target impact point coordinates.

[0098] See also Fig. 9 In a specific embodiment, step S105 includes steps S1051 to S1052, and each step is described in detail below.

[0099] Step S1051, determining the maximum value of the brightness function corresponding to each of the model points from the brightness function values.

[0100] In this embodiment, the maximum value of the brightness function corresponding to each model point is determined from all the brightness function values. There is only one maximum value of the brightness function corresponding to one model point.

[0101] Step S1052, determining a plurality of target function values ​​from the maximum values ​​of the brightness functions according to a preset source threshold.

[0102] In this embodiment, the maximum value of the brightness function may correspond to the real peak of the earthquake source, or it may be an incidental peak that is unrelated to the earthquake source. It is necessary to use a preset earthquake source threshold to determine multiple target function values ​​from each brightness function maximum value, and the target function value is the real brightness peak value.

[0103] In a specific embodiment, step S1052 includes: determining a plurality of brightness function peak groups from each of the brightness function maximum values ​​according to the preset source threshold; determining a plurality of real brightness peaks according to each of the brightness function peak groups, and using the real brightness peaks as the target function values.

[0104] In this embodiment, multiple brightness function peak groups are determined from the maximum values ​​of each brightness function according to a preset source threshold; each brightness function peak group is processed, the accompanying peaks and secondary peaks are screened out, and multiple real brightness peaks are determined, and the real brightness peaks are used as the target function values.

[0105] For example, see Fig.10 The horizontal dotted line is the preset source threshold line, and the empirical threshold is 0.85. Each brightness function peak group exceeding the preset source threshold corresponds to at least one impact point.

[0106] Step S1053, determining the coordinates of each target impact point according to each objective function value.

[0107] In this embodiment, the source scanning algorithm (SSA) is used to solve multiple objective function values ​​to obtain the corresponding target impact point coordinates. The multi-target positioning results are as follows: Fig.11 shown.

[0108] The multi-target impact point positioning method proposed in this embodiment is to determine multiple model points according to a preset spacing in the epicenter area; obtain the shock wave signals of the shock wave from each of the model points to multiple detectors, and obtain the shock wave signal sequence corresponding to each of the detectors at each of the model points; obtain the time difference translation signal sequence corresponding to each of the detectors at each of the model points according to each of the shock wave signal sequences; obtain multiple brightness function values ​​corresponding to each of the model points according to each of the time difference translation signal sequences; and determine the coordinates of multiple target impact points according to each of the brightness function values. In this way, multiple model points are set in the epicenter area, and the original shock wave signal is translated using the time difference of the shock wave signal corresponding to each of the model points, so that the shock wave signal at the epicenter can be enhanced, and the impact point position corresponding to the enhanced shock wave signal is found using the brightness function, solving the problem of impact point positioning caused by the dual effects of disorder and aliasing of multi-target signals.

[0109] Example 2

[0110] In addition, the present disclosure provides a multi-target impact point positioning device 1200, see Fig.12 , the device comprises:

[0111] A determination module 1201 is used to determine a plurality of model points in the earthquake source area according to a preset spacing;

[0112] An acquisition module 1202 is used to acquire the vibration wave signals from each of the model points to a plurality of detectors, and obtain a vibration wave signal sequence corresponding to each of the detectors at each of the model points;

[0113] The translation module 1203 is used to obtain the time difference translation signal sequence corresponding to each detector at each model point according to each vibration wave signal sequence;

[0114] A calculation module 1204 is used to obtain a plurality of brightness function values ​​corresponding to each of the model points according to each of the time difference shift signal sequences;

[0115] The positioning module 1205 is used to determine the coordinates of multiple target impact points according to each of the brightness function values.

[0116] Optionally, the translation module 1203 is also used to determine, for each of the model points, the relative arrival time difference of the shock wave of each detector according to the signal reception time corresponding to each of the shock wave signal sequences; based on a preset time window length, translate each of the shock wave signal sequences according to the relative arrival time difference of the shock wave to obtain a plurality of arrival time difference translation signal sequences.

[0117] Optionally, the translation module 1203 is further used to determine a reference object from each of the detectors; and determine the relative time difference of the shock wave of each detector relative to the reference object according to the signal reception time corresponding to each of the shock wave signal sequences.

[0118] Optionally, the translation module 1203 is further used to normalize the amplitude of each of the vibration wave signal sequences to obtain a plurality of pre-processed signal sequences; and determine the relative arrival time difference of the vibration wave of each detector relative to the reference object according to the signal reception time corresponding to each of the pre-processed signal sequences.

[0119] Optionally, the positioning module 1205 is also used to determine the maximum value of the brightness function corresponding to each of the model points from the brightness function values; determine multiple target function values ​​from the maximum values ​​of the brightness function according to a preset source threshold; and determine the coordinates of each target impact point according to each of the target function values.

[0120] Optionally, the positioning module 1205 is further used to determine a brightness function peak group from each of the brightness function maximum values ​​according to the preset source threshold; determine multiple real brightness peaks according to the brightness function peak group, and use the real brightness peaks as the target function value.

[0121] Optionally, the calculation module 1204 is also used to obtain, for each of the model points, based on multiple time points, the absolute value of the shock wave signal corresponding to each of the time points in the time difference translation signal sequence; and obtain multiple brightness function values ​​corresponding to the model point according to the absolute value of each of the shock wave signals.

[0122] The device provided in the embodiment of the present disclosure can execute the steps of the multi-target impact point positioning method provided in Example 1, which will not be described again to avoid repetition.

[0123] The multi-target impact point positioning device proposed in this embodiment determines multiple model points according to a preset spacing in the epicenter area; obtains the shock wave signals of the shock waves from each of the model points to multiple detectors, and obtains the shock wave signal sequence corresponding to each of the detectors at each of the model points; obtains the time difference translation signal sequence corresponding to each of the detectors at each of the model points according to each of the shock wave signal sequences; obtains multiple brightness function values ​​corresponding to each of the model points according to each of the time difference translation signal sequences; and determines the coordinates of multiple target impact points according to each of the brightness function values. In this way, multiple model points are set in the epicenter area, and the original shock wave signals are translated using the time difference of the shock wave signals corresponding to each of the model points, so that the shock wave signals at the epicenter can be enhanced, and the impact point position corresponding to the enhanced shock wave signals is found using the brightness function, solving the problem of impact point positioning caused by the dual effects of disorder and aliasing of multi-target signals.

[0124] Example 3

[0125] In addition, an embodiment of the present disclosure provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the multi-target impact point positioning method described in Example 1 is implemented.

[0126] The device provided in the embodiment of the present disclosure can execute the steps of the multi-target impact point positioning method provided in Example 1, which will not be described again to avoid repetition.

[0127] Example 4

[0128] The embodiment of the present disclosure provides a computer-readable storage medium storing a computer program. When the computer program is executed by a processor, the multi-target impact point positioning method described in Embodiment 1 is implemented.

[0129] In this embodiment, the computer-readable storage medium may be a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0130] The computer-readable storage medium provided in this embodiment can implement the multi-target impact point positioning method provided in Example 1, and will not be described again here to avoid repetition.

[0131] In all examples shown and described herein, any specific values ​​should be interpreted as merely exemplary and not as limiting, and thus other examples of the exemplary embodiments may have different values.

[0132] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, further definition and explanation thereof is not required in subsequent drawings.

[0133] The above-mentioned embodiments only express several implementation methods of the present invention, and the description thereof is relatively specific and detailed, but it cannot be understood as limiting the scope of the present invention. It should be pointed out that for ordinary technicians in this field, several modifications and improvements can be made without departing from the concept of the present invention, which all belong to the protection scope of the present invention.

Claims

1. A multi-target impact point positioning method, characterized in that: include: Determine multiple model points in the earthquake source area according to preset spacing; Acquire the vibration wave signals from each of the model points to the multiple detectors, and obtain the vibration wave signal sequence corresponding to each of the detectors at each of the model points; According to each of the vibration wave signal sequences, a time difference translation signal sequence corresponding to each of the detectors at each of the model points is obtained; Obtaining a plurality of brightness function values ​​corresponding to each of the model points according to each of the time difference shift signal sequences; A plurality of target impact point coordinates are determined according to the brightness function values.

2. The multi-target impact point positioning method according to claim 1, characterized in that: The step of obtaining the time difference translation signal sequence corresponding to each detector at each model point according to each vibration wave signal sequence includes: For each of the model points, determining the relative arrival time difference of the vibration waves of each of the detectors according to the signal reception time corresponding to each of the vibration wave signal sequences; Based on a preset time window length, each of the shock wave signal sequences is shifted according to the relative arrival time difference of the shock wave to obtain a plurality of arrival time difference shifted signal sequences.

3. The multi-target impact point positioning method according to claim 2, characterized in that: The step of determining the relative arrival time difference of the vibration waves of each detector according to the signal reception time corresponding to each vibration wave signal sequence comprises: determining a reference object from each of said detectors; The relative arrival time difference of the vibration wave of each detector relative to the reference object is determined according to the signal receiving time corresponding to each vibration wave signal sequence.

4. The multi-target impact point positioning method according to claim 3, characterized in that: The step of determining the relative arrival time difference of the vibration waves of each detector relative to the reference object according to the signal reception time corresponding to each vibration wave signal sequence comprises: Normalizing the amplitude of each of the shock wave signal sequences to obtain a plurality of pre-processed signal sequences; The relative arrival time difference of the vibration wave of each detector relative to the reference object is determined according to the signal receiving time corresponding to each pre-processed signal sequence.

5. The multi-target impact point positioning method according to claim 1, characterized in that: Determining a plurality of target impact point coordinates according to each of the brightness function values ​​comprises: Determine the maximum value of the brightness function corresponding to each of the model points from the brightness function values; Determining a plurality of objective function values ​​from the maximum values ​​of the brightness functions according to a preset source threshold; The coordinates of each target impact point are determined according to each objective function value.

6. The multi-target impact point positioning method according to claim 5, characterized in that: The step of determining a plurality of objective function values ​​from the maximum values ​​of the brightness functions according to a preset source threshold comprises: Determine a brightness function peak group from each of the brightness function maximum values ​​according to the preset source threshold; A plurality of real brightness peak values ​​are determined according to the brightness function peak value group, and the real brightness peak values ​​are used as the target function values.

7. The multi-target impact point positioning method according to claim 1, characterized in that: The step of obtaining a plurality of brightness function values ​​corresponding to each of the model points according to each of the time difference shift signal sequences comprises: For each of the model points, based on multiple time points, the absolute value of the shock wave signal corresponding to each of the time points is obtained from each of the time difference translation signal sequences; A plurality of brightness function values ​​corresponding to the model points are obtained according to the absolute values ​​of the shock wave signals.

8. A multi-target impact point positioning device, characterized in that: include: A determination module, used to determine a plurality of model points in the earthquake source area according to a preset spacing; An acquisition module is used to acquire the vibration wave signals from each of the model points to a plurality of detectors, and obtain the vibration wave signal sequence corresponding to each of the detectors at each of the model points; A translation module, used for obtaining a time difference translation signal sequence corresponding to each detector at each model point according to each vibration wave signal sequence; A calculation module, used for obtaining a plurality of brightness function values ​​corresponding to each of the model points according to each of the time difference shift signal sequences; A positioning module is used to determine the coordinates of multiple target impact points according to each of the brightness function values.

9. A computer device, characterized in that: It comprises a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the multi-target impact point positioning method as described in any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: It stores a computer program, which, when executed by a processor, implements the multi-target impact point positioning method as described in any one of claims 1 to 7.

Images