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

Abstract

The application 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 positions according to a preset interval in a seismic source area; obtaining vibration wave signals of vibration waves from each model point position to a plurality of detectors, to obtain vibration wave signal sequences corresponding to each detector under each model point position; obtaining time difference translation signal sequences corresponding to each detector under each model point position according to each vibration wave signal sequence; obtaining a plurality of brightness function values corresponding to each model point position according to each time difference translation signal sequence; and determining a plurality of target impact point coordinates according to each brightness function value. In this way, the time difference of the vibration wave signal is used to find the impact point position, and the problem of impact point positioning caused by the out-of-order and aliasing double effects of the multi-target signal is solved.

Description

Technical Field

[0001] This invention relates to the field of seismic source location technology, and in particular to a method, apparatus, computer equipment, and storage medium for locating the impact points of multiple targets. Background Technology

[0002] Existing bullet impact point positioning systems use multi-station intersection positioning methods to locate the bullet impact point coordinates, but they are only limited to single-target positioning of one bullet impact point and are powerless to locate multiple bullet 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 shortcomings of the prior art and provide a method, apparatus, computer equipment and storage medium for locating multiple projectile impact points by utilizing the detection time difference of the detector.

[0004] This invention provides the following technical solution:

[0005] In a first aspect, the present invention proposes a method for locating the impact point of multiple targets, comprising:

[0006] Within the seismic source area, multiple model points are determined according to preset spacing;

[0007] The vibration wave signals from each of the model points to multiple detectors are obtained, and the vibration wave signal sequence corresponding to each detector at each of the model points is obtained.

[0008] Based on the vibration wave signal sequences, the time difference translation signal sequences corresponding to the detectors at each of the model points are obtained;

[0009] Based on the time difference translation signal sequence, multiple brightness function values ​​corresponding to each model point are obtained;

[0010] The coordinates of multiple target impact points are determined based on the brightness function values ​​described above.

[0011] In one embodiment, obtaining the time difference translation signal sequence corresponding to each detector at each model point based on each of the vibration wave signal sequences includes:

[0012] For each of the model points, the relative arrival time difference of the vibration waves of each detector is determined according to the signal reception time corresponding to each vibration wave signal sequence;

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

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

[0015] Determine the reference object from each of the aforementioned detectors;

[0016] The relative time difference of the vibration waves of each detector relative to the reference object is determined based on the signal reception time corresponding to each vibration wave signal sequence.

[0017] In one embodiment, determining the relative time difference of the vibration waves of each detector relative to the reference object based on the signal reception time corresponding to each of the vibration wave signal sequences includes:

[0018] Amplitude normalization is performed on each of the aforementioned vibration wave signal sequences to obtain multiple preprocessed signal sequences;

[0019] The relative time difference of the vibration waves of each detector relative to the reference object is determined based on the signal reception time corresponding to each of the preprocessed signal sequences.

[0020] In one embodiment, determining the coordinates of multiple target impact points based on each of the brightness function values ​​includes:

[0021] Determine the maximum value of the brightness function corresponding to each of the aforementioned brightness function values;

[0022] Multiple target function values ​​are determined from the maximum values ​​of each of the aforementioned brightness functions based on a preset source threshold.

[0023] The coordinates of the impact points of each target are determined based on the objective function values.

[0024] In one embodiment, determining multiple target function values ​​from the maximum values ​​of each of the brightness functions based on a preset source threshold includes:

[0025] Based on the preset source threshold, determine the peak group of the brightness function from the maximum values ​​of each brightness function;

[0026] Multiple true brightness peaks are determined based on the peak group of the brightness function, and the true brightness peaks are used as the target function value.

[0027] In one embodiment, obtaining multiple brightness function values ​​corresponding to each model point based on each of the time difference translation signal sequences includes:

[0028] For each of the model points, based on multiple time points, the absolute value of the vibration wave signal corresponding to each time point is obtained from the time difference translation signal sequence.

[0029] Based on the absolute values ​​of each vibration wave signal, multiple brightness function values ​​corresponding to the model points are obtained.

[0030] Secondly, the present invention proposes a multi-target impact point positioning device, comprising:

[0031] The determination module is used to determine multiple model points within the seismic source area according to preset intervals.

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

[0033] The translation module is used to obtain the time difference translation signal sequence corresponding to each detector at each model point based on each of the vibration wave signal sequences;

[0034] The calculation module is used to obtain multiple brightness function values ​​corresponding to each model point based on the time difference translation signal sequence.

[0035] The positioning module is used to determine the coordinates of multiple target impact points based on the brightness function values.

[0036] Thirdly, the present invention proposes 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, it implements the multi-target impact point localization method as described in the first aspect.

[0037] Fourthly, the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the multi-target impact point localization method as described in the first aspect.

[0038] This invention discloses a method, apparatus, computer device, and storage medium for locating multiple target impact points. Within a seismic source region, multiple model points are determined at preset intervals. Seismic wave signals from each model point to multiple detectors are acquired, resulting in a seismic wave signal sequence corresponding to each detector at each model point. Based on these seismic wave signal sequences, a time difference translation signal sequence corresponding to each detector at each model point is obtained. Multiple brightness function values ​​corresponding to each model point are obtained based on these time difference translation signal sequences. Finally, the coordinates of multiple target impact points are determined based on these brightness function values. By setting multiple model points within the seismic source region and using the time difference of the seismic wave signals corresponding to each model point to translate the original seismic wave signals, the seismic wave signals at the seismic source are enhanced. Simultaneously, the brightness function is used to find the impact point position corresponding to the enhanced seismic wave signal, solving the problem of impact point location caused by the dual effects of disordered and aliased multi-target signals. Attached Figure Description

[0039] To more clearly illustrate the technical solution of the present invention, the accompanying drawings used in the embodiments will be briefly described below. It should be understood that the following drawings only show some embodiments of the present invention and should not be regarded as a limitation on the scope of protection of the present invention. In the various drawings, similar components are numbered similarly.

[0040] Figure 1 A flowchart of the multi-target impact point localization method proposed in this embodiment is shown;

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

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

[0043] Figure 4 Another flowchart of the multi-target impact point localization method proposed in this embodiment is shown;

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

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

[0046] Figure 7 A schematic diagram of the vibration wave signal sequence and preprocessed signal sequence proposed in this embodiment is shown;

[0047] Figure 8 This embodiment shows a schematic diagram of multiple brightness function values ​​corresponding to a model point.

[0048] Figure 9 This diagram illustrates another step in the multi-target impact point localization method proposed in this embodiment.

[0049] Figure 10 A schematic diagram of the peak group of the brightness function proposed in this embodiment is shown;

[0050] Figure 11 A schematic diagram of the multi-target localization results proposed in this embodiment is shown;

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

[0052] Explanation of reference numerals in the attached diagram:

[0053] 1200 - Multi-target impact point positioning device; 1201 - Determination module; 1202 - Acquisition module; 1203 - Translation module; 1204 - Calculation module; 1205 - Positioning module. Detailed Implementation

[0054] The technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments.

[0055] The components of the embodiments of the invention described and illustrated herein can typically be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely to illustrate selected embodiments of the invention. All other embodiments obtained by those skilled in the art based on the embodiments of the invention without inventive effort are within the scope of protection of the invention.

[0056] In the following, the terms “comprising,” “having,” and their cognates, which may be used in various embodiments of the invention, are intended only to indicate a particular feature, number, step, operation, element, component, or combination thereof, and should not be construed as excluding, firstly, the presence of one or more other features, numbers, steps, operations, elements, components, or combinations thereof, or adding the possibility of one or more features, numbers, steps, operations, elements, components, or combinations thereof.

[0057] Furthermore, the terms "first," "second," and "third" are used only to distinguish descriptions and should not be interpreted as indicating or implying relative importance.

[0058] Unless otherwise specified, all terms used herein (including technical and scientific terms) shall have the same meaning as commonly understood by one of ordinary skill in the art to which the various embodiments of the invention pertain. Terms (such as those defined in commonly used dictionaries) shall be interpreted as having the same meaning as in their contextual meaning in the relevant technical field and shall not be interpreted as having an idealized or overly formal meaning, unless clearly defined in the various embodiments of the invention.

[0059] Example 1

[0060] Before the countdown to the multi-target positioning test begins, the user runs the multi-target impact point positioning software. The data processing computer runs the software, initiates initialization and self-test, and imports the time difference model into memory after the system is running normally, awaiting target launch. Based on the launch command, the user activates the signal acquisition function. The data processing computer sends the acquisition command to the wireless data transmission network communication interface. The wireless data transmission network transmits the acquisition command to each detector. Each detector executes the signal acquisition command and performs preprocessing such as filtering, digital-to-analog conversion, and time synchronization on the signal before packaging and sending the data to the network communication interface. After receiving the packaged data from the detector, the wireless data transmission network transmits it to the data processing computer. The data processing computer receives the packaged data, unpacks it, and then displays the dynamically acquired signals in real time on the screen, allowing the user to monitor the acquired signals from each detector. Based on the launch termination command, the user confirms receipt of all vibration wave signals, stops the signal acquisition function, activates the positioning function, and executes the multi-target impact point positioning method.

[0061] This disclosure provides a method for locating multiple projectile impact points by utilizing the detection time difference of a detector.

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

[0063] Step S101: Determine multiple model points within the seismic source area according to a preset spacing.

[0064] In this embodiment, multiple model points are determined within the seismic source region based on preset spacing. For example, a 400m × 400m seismic source region is divided into 10201 (101 × 101) model points at preset spacings of 0.5m. A model point layout diagram is shown below. Figure 2 As shown.

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

[0066] Step S102: Obtain the vibration wave signals from each of the model points to the multiple detectors, and obtain the vibration wave signal sequence corresponding to each detector at each of the model points.

[0067] In this embodiment, multiple detectors are arranged to detect seismic waves. Seismic waves are generated by seismic sources within the seismic source region. When each detector executes a signal acquisition command, it uses each model point as the simulated propagation starting point of the seismic wave, and each detector records the seismic wave signal transmitted from each model point. These seismic wave signals include information such as the arrival time, waveform, and amplitude of the seismic wave.

[0068] For each model point, each detector receives and records seismic wave signals at different time scales transmitted for that model location, and assembles a seismic wave signal sequence based on these signals. For example, if there are model points A1 to A2... 10201 And 8 detectors, then model point A i The seismic wave signal sequence under (i = 1, 2, ..., 10201) includes u1 to u8, and the seismic wave signal sequence u j (j=1,2,……,8) is a sequence composed of vibration wave signals received by the j-th detector.

[0069] Please see Figure 3 In the diagram, the square represents the earthquake source area, the circles inside the square represent the 10 impact points, and the rhombuses outside the square represent the locations of the 8 detectors.

[0070] Step S103: Obtain the time difference translation signal sequence corresponding to each detector at each model point based on each vibration wave signal sequence.

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

[0072] Since the seismic wave signals from the detector at the earthquake 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] Please see Figure 4 In one specific embodiment, step S103 includes steps S1031 to S1032, and each step will be described in detail below.

[0074] Step S1031: For each of the model points, determine the relative arrival time difference of the vibration wave of each detector based on the signal reception time corresponding to each vibration wave signal sequence.

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

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

[0077] In this embodiment, a reference object is determined from among the various detectors, and the relative arrival time difference of the vibration waves of the other detectors relative to the reference object is further determined based on the signal reception time corresponding to each vibration wave signal sequence. Specifically, a single reference object should be used to obtain the relative arrival time difference of the vibration waves for each model point.

[0078] Exemplary, τ η1 and τ ηj From model point A respectively i The travel time of the vibration wave η to detectors 1 and j is calculated, with detector 1 as the reference object, and the relative travel time τ of the vibration wave received by detector j relative to detector 1. j1 For: τ j1 =τ ηj -τ η1 .

[0079] Step S1032: Based on the preset time window length, each vibration wave signal sequence is shifted according to the relative arrival time difference of the vibration wave to obtain multiple arrival time difference shifted signal sequences.

[0080] In this embodiment, a signal time window is selected based on the measured waveform, and a preset time window length of m is set; further, each vibration wave signal sequence is shifted according to the preset time window length m and the relative arrival time difference of the vibration wave, resulting in multiple arrival time difference shifted signal sequences.

[0081] As an example, taking 8 detectors as an example, for model point A i The individual time difference shift signal sequences are as follows:

[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] Where T0 is the start time of the preset time window, and the sequence number of the start time of the corresponding time window.

[0087] Please see Figure 5 Due to the aliasing of different target signals, not every channel in the seismic wave signal sequence can clearly distinguish 10 target signals. Taking detector 1 as the reference object, the other 7 detectors, according to their respective time difference shifted signal sequences, are as follows: Figure 6 As shown.

[0088] In one specific embodiment, the amplitude of each of the vibration wave signal sequences is normalized to obtain multiple preprocessed signal sequences; the relative time difference of the vibration waves of each detector relative to the reference object is determined according to the signal reception time corresponding to each of the preprocessed signal sequences.

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

[0090] For multi-target localization, the effects of distance differences must be eliminated to ensure that signals from the same seismic source have the same impact on the brightness function. Therefore, before translating the signal sequence, the original seismic wave signal sequence must undergo amplitude normalization. Please refer to [link to relevant documentation]. Figure 7 (a) is the seismic wave signal sequence, and (b) is the preprocessed signal sequence.

[0091] Step S104: Obtain multiple brightness function values ​​corresponding to each model point based on the time difference translation signal sequence.

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

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

[0094] In this embodiment, for each model point, the absolute values ​​of the seismic wave signals at different time points in the time difference translation signal sequence are obtained, and multiple brightness function values ​​corresponding to the model point are obtained based on the absolute values ​​of each seismic wave signal. The multiple brightness function values ​​corresponding to one model point are as follows: Figure 8 As shown.

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

[0096] Step S105: Determine the coordinates of multiple target impact points based on the brightness function values.

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

[0098] Please see Figure 9 In one specific embodiment, step S105 includes steps S1051 to S1052, and each step is described in detail below.

[0099] Step S1051: Determine the maximum value of the brightness function corresponding to each model point 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 brightness function values. Each model point has one and only one maximum value of the brightness function.

[0101] Step S1052: Determine multiple target function values ​​from the maximum values ​​of each brightness function according to the preset source threshold.

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

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

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

[0105] For example, please see Figure 10 The horizontal dashed line is the preset source threshold line. 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: Determine the coordinates of each target impact point based on each of the target function values.

[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 localization results are as follows: Figure 11 As shown.

[0108] The multi-target impact point localization method proposed in this embodiment determines multiple model points within the seismic source region according to a preset spacing; acquires seismic wave signals from each model point to multiple detectors, obtaining a seismic wave signal sequence corresponding to each detector at each model point; obtains a time difference translation signal sequence corresponding to each detector at each model point based on the seismic wave signal sequence; obtains multiple brightness function values ​​corresponding to each model point based on the time difference translation signal sequence; and determines the coordinates of multiple target impact points based on the brightness function values. In this way, by setting multiple model points within the seismic source region and using the time difference of the seismic wave signals corresponding to each model point to translate the original seismic wave signals, the seismic wave signals at the seismic source can be enhanced. Simultaneously, the brightness function is used to find the impact point position corresponding to the enhanced seismic wave signal, solving the problem of impact point localization caused by the dual effects of signal disorder and aliasing in multi-target systems.

[0109] Example 2

[0110] Furthermore, this disclosure provides a multi-target impact point positioning device 1200, please refer to [link to relevant documentation]. Figure 12 The device includes:

[0111] The module 1201 is used to determine multiple model points within the seismic source region according to a preset spacing.

[0112] The acquisition module 1202 is used to acquire the vibration wave signals from each of the model points to multiple detectors, and to obtain the vibration wave signal sequence corresponding to each detector 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 based on each of the vibration wave signal sequences;

[0114] Calculation module 1204 is used to obtain multiple brightness function values ​​corresponding to each model point based on each of the time difference translation signal sequences;

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

[0116] Optionally, the translation module 1203 is further configured to determine the relative arrival time difference of the vibration wave of each detector based on the signal reception time corresponding to each vibration wave signal sequence for each model point; and to translate each vibration wave signal sequence according to the relative arrival time difference based on a preset time window length to obtain multiple arrival time difference translation signal sequences.

[0117] Optionally, the translation module 1203 is further configured to determine a reference object from each of the detectors; and to determine the relative time difference of the vibration waves of each detector relative to the reference object based on the signal reception time corresponding to each vibration wave signal sequence.

[0118] Optionally, the translation module 1203 is further configured to perform amplitude normalization on each of the vibration wave signal sequences to obtain multiple preprocessed signal sequences; and to determine the relative time difference of the vibration waves of each detector relative to the reference object based on the signal reception time corresponding to each of the preprocessed signal sequences.

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

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

[0121] Optionally, the calculation module 1204 is further configured to, for each of the model points, obtain the absolute value of the vibration wave signal corresponding to each of the time points from the time difference translation signal sequence; and obtain multiple brightness function values ​​corresponding to the model points based on the absolute values ​​of the vibration wave signals.

[0122] The apparatus provided in this embodiment can execute the steps of the multi-target impact point positioning method provided in Embodiment 1. To avoid repetition, it will not be described again.

[0123] The multi-target impact point positioning device proposed in this embodiment determines multiple model points within the seismic source region according to a preset spacing; acquires seismic wave signals from each model point to multiple detectors, obtaining a seismic wave signal sequence corresponding to each detector at each model point; obtains a time difference translation signal sequence corresponding to each detector at each model point based on the seismic wave signal sequence; obtains multiple brightness function values ​​corresponding to each model point based on the time difference translation signal sequence; and determines the coordinates of multiple target impact points based on the brightness function values. In this way, by setting multiple model points within the seismic source region and using the time difference of the seismic wave signals corresponding to each model point to translate the original seismic wave signals, the seismic wave signals at the seismic source are enhanced. Simultaneously, the brightness function is used to find the impact point position corresponding to the enhanced seismic wave signal, solving the problem of impact point positioning caused by the dual effects of signal disorder and aliasing in multi-target systems.

[0124] Example 3

[0125] Furthermore, this disclosure provides a computer device including a memory and a processor. The memory stores a computer program, which, when executed by the processor, implements the multi-target impact point localization method described in Embodiment 1.

[0126] The device provided in this embodiment can execute the steps of the multi-target impact point positioning method provided in Embodiment 1. To avoid repetition, it will not be described again.

[0127] Example 4

[0128] This disclosure provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the multi-target impact point localization method described in Embodiment 1.

[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, etc.

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

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

[0132] It should be noted that similar labels and letters in the following figures indicate similar items. Therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures.

[0133] The above-described embodiments are merely illustrative of several implementations of the present invention, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention.

Claims

1. A multi-target impact point location method, characterized by, The method comprises the following steps: determining a plurality of model points in a source region according to a preset interval; obtaining vibration wave signals of vibration waves from each of the model points to a plurality of detectors, to obtain a vibration wave signal sequence corresponding to each of the detectors under each of the model points; obtaining a time difference translation signal sequence corresponding to each of the detectors under each of the model points according to each of the vibration wave signal sequences; obtaining a plurality of brightness function values corresponding to each of the model points according to each of the time difference translation signal sequences; determining a plurality of target impact point coordinates according to each of the brightness function values; the step of obtaining a plurality of brightness function values corresponding to each of the model points according to each of the time difference translation signal sequences comprises the following steps: for each of the model points, obtaining an absolute value of a vibration wave signal corresponding to each of a plurality of time points from each of the time difference translation signal sequences based on the time points; obtaining a plurality of brightness function values corresponding to the model point according to each of the vibration wave signal absolute values.

2. The multi-target impact point location method of claim 1, wherein, the step of obtaining a time difference translation signal sequence corresponding to each of the detectors under each of the model points according to each of the vibration wave signal sequences comprises the following steps: for each of the model points, determining a relative time difference of a vibration wave of each of the detectors according to a signal receiving time corresponding to each of the vibration wave signal sequences; based on a preset time window length, translating each of the vibration wave signal sequences according to the relative time difference of the vibration wave, to obtain a plurality of time difference translation signal sequences.

3. The multi-target impact point location method of claim 2, wherein, the step of determining a relative time difference of a vibration wave of each of the detectors according to a signal receiving time corresponding to each of the vibration wave signal sequences comprises the following steps: determining a reference object from each of the detectors; determining a relative time difference of a vibration wave of each of the detectors with respect to the reference object according to a signal receiving time corresponding to each of the vibration wave signal sequences.

4. The multi-target impact point location method of claim 3, wherein, the step of determining a relative time difference of a vibration wave of each of the detectors with respect to the reference object according to a signal receiving time corresponding to each of the vibration wave signal sequences comprises the following steps: amplitude normalization is performed on each of the vibration wave signal sequences to obtain a plurality of preprocessed signal sequences; determining a relative time difference of a vibration wave of each of the detectors with respect to the reference object according to a signal receiving time corresponding to each of the preprocessed signal sequences.

5. The multi-target impact point location method of claim 1, wherein, the step of determining a plurality of target impact point coordinates according to each of the brightness function values comprises the following steps: determining a maximum value of a brightness function corresponding to each of the model points from each of the brightness function values; determining a plurality of target function values from each of the maximum values of the brightness function according to a preset source threshold; determining a plurality of target impact point coordinates according to each of the target function values.

6. The multi-target impact point location method of claim 5, wherein, the step of determining a plurality of target function values from each of the maximum values of the brightness function according to a preset source threshold comprises the following steps: determining a peak value group of the brightness function from each of the maximum values of the brightness function according to the preset source threshold; determining a plurality of real brightness peak values from the peak value group of the brightness function, and taking the real brightness peak values as the target function values.

7. A multiple target impact point location device characterized by, The method comprises the following steps: a determining module is configured to determine a plurality of model points in a source region according to a preset interval; an obtaining module is configured to obtain vibration wave signals of vibration waves from each of the model points to a plurality of detectors, to obtain a vibration wave signal sequence corresponding to each of the detectors under each of the model points; a translation module configured to obtain a time-difference translation signal sequence corresponding to the detector at each model point according to each of the seismic wave signal sequences; a calculation module configured to obtain a plurality of brightness function values corresponding to each model point according to each of the time-difference translation signal sequences; a positioning module configured to determine a plurality of target impact point coordinates according to each of the brightness function values; the calculation module is further configured to, for each model point, obtain a seismic wave signal absolute value corresponding to each time point from each of the time-difference translation signal sequences based on a plurality of time points, and obtain a plurality of brightness function values corresponding to the model point according to each of the seismic wave signal absolute values.

8. A computer device, comprising: a memory and a processor, the memory storing a computer program, and the computer program being executed by the processor to implement the multi-target impact point positioning method according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, a memory storing a computer program, and the computer program being executed by a processor to implement the multi-target impact point positioning method according to any one of claims 1 to 6.

Images