CN111880220B - Seismic source positioning method, device, equipment and storage medium - Google Patents

Seismic source positioning method, device, equipment and storage medium Download PDF

Info

Publication number
CN111880220B
CN111880220B CN202010927795.6A CN202010927795A CN111880220B CN 111880220 B CN111880220 B CN 111880220B CN 202010927795 A CN202010927795 A CN 202010927795A CN 111880220 B CN111880220 B CN 111880220B
Authority
CN
China
Prior art keywords
positioning
wave
calculating
wave velocity
positioning error
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.)
Active
Application number
CN202010927795.6A
Other languages
Chinese (zh)
Other versions
CN111880220A (en
Inventor
李涛
陈炳瑞
王睿
朱新豪
王旭
王庆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan Tibet Railway Co ltd
Wuhan Institute of Rock and Soil Mechanics of CAS
China State Railway Group Co Ltd
Original Assignee
Sichuan Tibet Railway Co ltd
Wuhan Institute of Rock and Soil Mechanics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sichuan Tibet Railway Co ltd, Wuhan Institute of Rock and Soil Mechanics of CAS filed Critical Sichuan Tibet Railway Co ltd
Priority to CN202010927795.6A priority Critical patent/CN111880220B/en
Publication of CN111880220A publication Critical patent/CN111880220A/en
Application granted granted Critical
Publication of CN111880220B publication Critical patent/CN111880220B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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. analysis, for interpretation, for correction
    • 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 application belongs to earthquake monitoring, and relates to a seismic source positioning method, which comprises the following steps: acquiring time-of-arrival information of a seismic source and at least one microseismic signal; calculating at least one positioning error corresponding to the positioning function through at least one positioning function and at least one arrival time information of the microseismic signals; repeatedly calculating a plurality of corresponding wave velocities through each positioning error to obtain at least one wave velocity set corresponding to the positioning error; calculating the range of each wave velocity set, and taking the wave velocity set with the minimum range as a target set; if the target set is an empty set, selecting a corresponding positioning function with the minimum positioning error as a target function; if the target set is not an empty set, taking a CM (canonical method) as a target function; the source location is calculated by an objective function. The positioning precision under the complex geological condition is improved.

Description

Seismic source positioning method, device, equipment and storage medium
Technical Field
The invention relates to the field of seismic monitoring, in particular to a method, a device, equipment and a storage medium for seismic source positioning.
Background
The micro seismic source positioning is an important component of micro seismic monitoring and disaster early warning, and the selection of a micro seismic source positioning target function is one of key problems of micro seismic source positioning. Experts and scholars at home and abroad do much research work aiming at the objective function, and the main points are as follows: knowing a wave velocity model, and solving a classical method of the origin time and the origin position; a combination method of solving the microseismic source position, the origin time and the wave velocity model together, and the like. The classical method is most widely applied in the seismic field and mining engineering, and inaccurate speed model feeding is the greatest defect of the method; the combination method well solves the problem of inaccurate velocity model, improves the positioning precision of the micro seismic source to a greater extent, but the mutual correlation of the parameters of the position of the micro seismic source, the seismic time and the medium velocity brings about the problem of unstable positioning result. In summary, although a great deal of research on objective functions has been conducted before, in practical engineering, there are advantages and disadvantages of each objective function, and how to select an objective function for positioning, what is different in positioning accuracy of different objective functions, and a systematic method is still needed for research.
Therefore, aiming at the defects, the invention provides the adaptive identification algorithm of the micro seismic source positioning objective function, which is very necessary for improving the positioning precision of the micro seismic source under the complex geological conditions and further improving the timeliness and accuracy of the early warning of geological disasters such as rock burst, mine earthquake, collapse and the like.
Disclosure of Invention
The method and the device can adapt to complex geological conditions, and improve the positioning precision under the complex geological conditions, thereby improving the timeliness and accuracy of the early warning of geological disasters such as rock burst, mine earthquake, collapse and the like.
In a first aspect, the present application provides a method of seismic source localization, comprising:
acquiring time-of-arrival information of a seismic source and at least one microseismic signal;
calculating at least one positioning error corresponding to the positioning function through at least one positioning function and at least one arrival time information of the microseismic signals;
repeatedly calculating a plurality of corresponding wave velocities through each positioning error to obtain at least one wave velocity set corresponding to the positioning error;
calculating the range of each wave velocity set, and taking the wave velocity set with the minimum range as a target set;
if the target set is an empty set, selecting the corresponding positioning function with the minimum positioning error as a target function;
if the target set is not an empty set, taking a CM (classical method) as a target function;
the source location is calculated by the objective function.
In some possible designs, after the obtaining arrival time information for the seismic source and the at least one microseismic signal, the method further comprises:
judging the uniqueness of the positioning results of the joint method JM, the improved joint method IJM and the time method TM according to the arrival time information division of the microseismic signals and the positions of the sensors, and taking the CM as a target function if the positioning results of the positioning methods are not unique.
In some possible designs, the calculating at least one positioning error corresponding to the positioning function by at least one positioning function and arrival time information of the at least one microseismic signal includes:
calculating a first positioning error w by JMJ；
Calculating a second positioning error w by IJMI；
Calculating a third positioning error w by TMT；
Calculating the fourth positioning error w by LM methodL。
In some possible designs, the repeatedly calculating a plurality of wave velocities corresponding to each of the positioning errors to obtain at least one wave velocity set corresponding to the positioning error includes:
by Δ w ═ xn-wiI calculate the wave velocity error, where wiIs wJ、wI、wTAnd wLAny one of (1), xnPositioning error of CM method;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset lower limit value of the wave speed, and obtaining at least one lower limit value of the wave speed through mode search;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset upper limit value of the wave speed, and obtaining an upper limit value corresponding to at least one lower limit value of the wave speed through the mode search;
and obtaining the at least one wave velocity set by combining the at least one wave velocity lower limit value and the upper limit value corresponding to the at least one wave velocity lower limit value.
In some possible designs, the calculating a range of each of the wave velocity sets, with the wave velocity set with the smallest range as a target set, includes:
obtaining the wJCorresponding first set of wave velocities (v)Jmin,vJmax) W toICorresponding first set of wave velocities (v)Imin,vImax) W toTCorresponding first set of wave velocities (v)Tmin,vTmax) And said wLCorresponding first set of wave velocities (v)Lmin,vLmax)；
Through min (v)Jmax-vJmin，vImax-vImin，vTmax-vTmin，vLmax-vLmin) Calculating the minimum value of the wave speed range in the at least one wave speed set;
and taking the wave velocity set corresponding to the minimum value of the wave velocity range as a target set.
In some possible designs, the calculating at least one positioning error corresponding to the positioning function from at least one positioning function and at least one arrival time information of the microseismic signals comprises:
by passing
Figure BDA0002669052020000031
Calculating at least one positioning error corresponding to said positioning function, wherein (x)1,y1,z1) Is the source coordinate, (x)0,y0,z0) The microseismic signal coordinates.
In some possible designs, before the acquiring time-of-arrival information for the seismic source and the at least one microseismic signal, the method further comprises:
and acquiring the coordinates of the seismic source through a particle swarm algorithm.
In a second aspect, the present application provides a seismic source positioning apparatus having the functionality of implementing a method corresponding to the seismic source positioning platform provided in the first aspect above. The functions can be realized by hardware, and the functions can also be realized by executing corresponding software by hardware. The hardware or software includes one or more modules corresponding to the above functions, which may be software and/or hardware.
The device for seismic source positioning comprises:
the input and output module is used for acquiring time-of-arrival information of the seismic source and the at least one microseismic signal;
a processing module for calculating at least one positioning error corresponding to the positioning function by at least one positioning function and at least one arrival time information of the microseismic signals; repeatedly calculating a plurality of corresponding wave velocities through each positioning error to obtain at least one wave velocity set corresponding to the positioning error; calculating the range of each wave velocity set, and taking the wave velocity set with the minimum range as a target set; if the target set is an empty set, selecting the corresponding positioning function with the minimum positioning error as a target function; if the target set is not an empty set, taking a CM (classical method) as a target function; the source location is calculated by the objective function.
In some possible designs, the processing module is further to:
judging the uniqueness of the positioning results of the joint method JM, the improved joint method IJM and the time method TM according to the arrival time information division of the microseismic signals and the positions of the sensors, and taking the CM as a target function if the positioning results of the positioning methods are not unique.
In some possible designs, the processing module is further to:
calculating a first positioning error w by JMJ；
Calculating a second positioning error w by IJMI；
Calculating a third positioning error w by TMT；
Calculating the fourth positioning error w by LM methodL。
In some possible designs, the processing module is further to:
by Δ w ═ xn-wiI calculate the wave velocity error, where wiIs wJ、wI、wTAnd wLAny one of (1), xnPositioning error of CM method;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset lower limit value of the wave speed, and obtaining at least one lower limit value of the wave speed through mode search;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset upper limit value of the wave speed, and obtaining an upper limit value corresponding to at least one lower limit value of the wave speed through the mode search;
and obtaining the at least one wave velocity set by combining the at least one wave velocity lower limit value and the upper limit value corresponding to the at least one wave velocity lower limit value.
In some possible designs, the processing module is further to:
obtaining the wJCorresponding first set of wave velocities (v)Jmin,vJmax) W toICorresponding first set of wave velocities (v)Imin,vImax) W toTCorresponding first set of wave velocities (v)Tmin,vTmax) And said wLCorresponding first set of wave velocities (v)Lmin,vLmax)；
Through min (v)Jmax-vJmin，vImax-vImin，vTmax-vTmin，vLmax-vLmin) Calculating the minimum value of the wave speed range in the at least one wave speed set;
and taking the wave velocity set corresponding to the minimum value of the wave velocity range as a target set.
In some possible designs, the processing module is further to:
by passing
Figure BDA0002669052020000041
Calculating at least one positioning error corresponding to said positioning function, wherein (x)1,y1,z1) Is the source coordinate, (x)0,y0,z0) The microseismic signal coordinates.
In some possible designs, the processing module is further to:
and acquiring the coordinates of the seismic source through a particle swarm algorithm.
In yet another aspect, the present application provides a seismic source locating apparatus, which includes at least one connected processor, a memory, and an input/output unit, wherein the memory is used for storing program codes, and the processor is used for calling the program codes in the memory to execute the method of the above aspects.
Yet another aspect of the present application provides a computer storage medium comprising instructions which, when executed on a computer, cause the computer to perform the method of the above-described aspects.
Compared with the prior art, the invention creatively provides an adaptive identification algorithm of the micro seismic source positioning target function aiming at the condition that how the micro seismic source positioning target function is selected to achieve better positioning effect, and the algorithm determines each partition positioning target function by comparing the positioning accuracy of a JM method, an IJM method, a TM method, an LM method of unknown wave velocity with that of a CM method of known wave velocity and judging whether the CM method has a 'wave velocity range'. The method can adapt to complex geological conditions, and simultaneously improves the positioning precision under complex geological conditions, thereby improving the timeliness and the accuracy of the early warning of geological disasters such as rock burst, mine earthquake, collapse and the like.
Drawings
FIG. 1-1 is a schematic flow chart of a method for seismic source location in an embodiment of the present application;
FIGS. 1-2 are schematic monitoring planes and zones for a method of seismic source location in an embodiment of the present application;
FIGS. 1-3 are tables of mine blasting test event information for a method of seismic source location in an embodiment of the present application;
FIGS. 1-4 are flowcharts of adaptive identification algorithms for micro-seismic source location objective functions for a method of seismic source location in an embodiment of the present application;
FIGS. 1-5 are comparison plots of positioning results of various positioning algorithms of a method of seismic source positioning in an embodiment of the present application;
FIGS. 1-6 illustrate the convergence of solutions for methods of the method of seismic source localization in embodiments of the present application;
FIGS. 1-7 are schematic monitoring planes and zones for a method of seismic source location in an embodiment of the present application;
FIGS. 1-8 are plots of positioning results versus wave velocity for methods of a method of seismic source positioning in an embodiment of the present application;
FIGS. 1-9 illustrate sensor coordinates and time of arrival of a shot event for a method of seismic source location in an embodiment of the present application;
FIG. 2 is a schematic diagram of an apparatus for seismic source location in an embodiment of the present application;
fig. 3 is a schematic structural diagram of a computer device in an embodiment of the present application.
The implementation, functional features and advantages of the objectives of the present application will be further explained with reference to the accompanying drawings.
Detailed Description
It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application. The terms "first," "second," and the like in the description and in the claims of the present application and in the above-described drawings are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It will be appreciated that the data so used may be interchanged under appropriate circumstances such that the embodiments described herein may be practiced otherwise than as specifically illustrated or described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or modules is not necessarily limited to those steps or modules explicitly listed, but may include other steps or modules not explicitly listed or inherent to such process, method, article, or apparatus, and such that a division of modules presented in this application is merely a logical division that may be implemented in an actual application in a different manner, such that multiple modules may be combined or integrated into another system, or some features may be omitted, or may not be implemented.
Referring to fig. 1-1, a method for seismic source localization provided by the present application is illustrated, the method comprising:
101. time-of-arrival information of the seismic source and the at least one microseismic signal is obtained.
In this embodiment, first, let i equal to 1, input the ith burst test signal of j partitions, the wave velocity range is (3000,9000), the wave velocity unit is m/s, and generate shockThe time range t e (t)min-0.5,tmin) Wherein t isminAnd setting the coordinate range according to the monitoring range of the j subarea, wherein the unit of the minimum arrival time of the trigger sensor is s. Dividing the whole mine into two monitoring areas, i.e. Q, according to the area and sensor array which are mainly considered to be monitored1Region and Q2The schematic views of zones, monitoring planes and zones are shown in FIGS. 1-2. At Q1Region and Q2The test zones are subjected to a burst test, test point 1 and test point 2, and the test point information is shown in table 2, wherein the test point 1 is (1876.68, 8582.30, -198) in coordinates and the test point 2 is (2073.15, 8579.18, -198) in coordinates. It can be found from table 2 that the number of sensors triggered by 2 blasting test events is 12 and 10 respectively, which are both greater than the requirement of 4 sensors. By dividing the sensor arrival times of the 2 times blasting test events, as shown in fig. 1-6, it can be found that the sensor arrival times of the 12 times blasting test events exceed 4 different arrival values, and by combining the sensor coordinates, it can be inferred that no multi-solution condition exists in the 2 times blasting test events.
102. Calculating at least one positioning error corresponding to the positioning function from at least one positioning function and at least one arrival time information of the microseismic signals.
In this embodiment, positioning is performed by using a Joint Method (JM), an Improved Joint Method (IJM), a Time Method (TM), and a positioning method (LM), and the positioning errors of the algorithms are wJ，wI，wT，wL。
103. And repeatedly calculating a plurality of corresponding wave velocities through each positioning error to obtain at least one wave velocity set corresponding to the positioning error.
In this embodiment, the positioning error of the CM method is calculated to reach wJThe wave velocity of the time is searched from v 3000m/s by using a pattern search method, the target function of the search is an equation (18), and the corresponding wave velocity when the equation (18) is 0 is vimin. The actual seismic source of the blasting is (x'0,y'0,z'0) The position of the micro seismic source is (x)0,y0,z0)。
Figure BDA0002669052020000071
Then, a pattern search method is used to start a search from v 9000m/s, and when equation (18) is 0, the corresponding wave velocity is defined as vimax(ii) a Obtaining the wave velocity range (v) required by the CM method that the positioning error is smaller than that of the JM methodimin,vimax) Put the wave velocity range into the wave velocity set VJIn (1). Similarly, the wave velocity range required by the CM method positioning error smaller than the IJM method positioning error is obtained through calculation, and the wave velocity range is put into the wave velocity set VI(ii) a Calculating to obtain a wave velocity range required by the CM method positioning error smaller than the TM method positioning error, and putting the wave velocity range into a wave velocity set VT(ii) a Calculating to obtain a wave velocity range required by the CM method positioning error smaller than the LM method positioning error, and putting the wave velocity range into a wave velocity set VLIn (1).
104. And calculating the range of each wave velocity set, and taking the wave velocity set with the minimum range as a target set.
In this embodiment, the wave velocity sets V are separately filteredJ、VI、VT、VLTaking intersection to obtain the wave velocity range (v) required by JM methodJmin,vJmax) IJM method required wave velocity range (v)Imin,vImax) The wave velocity range (v) required for the TM methodTmin,vTmax) Obtaining the wave velocity range (v) required by the LM methodLmin,vLmax). Contrast wave velocity range (v)Jmin,vJmax)、(vImin,vImax)、(vTmin,vTmax) And (v)Lmin,vLmax) Taking the minimum set as (v)min,vmax)。
105-1, if the target set is an empty set, selecting the corresponding positioning function with the minimum positioning error as the target function.
In the present embodiment, the range of the velocity of waves (v)min,vmax) When the method exists, the j partition is positioned by adopting a CM method, and the optimal uniform wave velocity required by the CM method is calculated according to mode search.
105-2, if the target set is not an empty set, taking a classical method CM as a target function;
in this example, w is comparedJ，wI，wT，wLTo obtain the minimum positioning error wminSelecting wminAnd carrying out positioning by a corresponding positioning algorithm.
106. The source location is calculated by the objective function.
The specific process is shown in fig. 1-4, after 9 times of blasting, the 9 times of blasting events of the Q1 region and the Q2 region are located by adopting the CM method, the JM method, the IJM method, the TM method and the LM method, and the spatial location comparison graph is shown in fig. 1-5. Comparing with FIGS. 1-5, the average positioning errors of the method of the present invention, CM method, JM method, IJM method, TM method and LM method are 14.62m, 20.47m, 103.23m, 103.30m, 19.04m and 208.31m, respectively; compared with a CM method, a JM method, an IJM method, a TM method and an LM method, the method of the invention is respectively improved by 25.58%, 85.84%, 85.85%, 23.21% and 92.98%, and therefore, the method of the invention has better positioning effect by adopting a self-adaptive identification algorithm of a micro seismic source positioning target function.
Compared with the prior art, the invention creatively provides an adaptive identification algorithm of the micro seismic source positioning target function aiming at the condition that how the micro seismic source positioning target function is selected to achieve better positioning effect, and the algorithm determines each partition positioning target function by comparing the positioning accuracy of a JM method, an IJM method, a TM method, an LM method of unknown wave velocity with that of a CM method of known wave velocity and judging whether the CM method has a 'wave velocity range'. The method can adapt to complex geological conditions, and simultaneously improves the positioning precision under complex geological conditions, thereby improving the timeliness and the accuracy of the early warning of geological disasters such as rock burst, mine earthquake, collapse and the like.
In some embodiments, after the obtaining the arrival time information of the seismic source and the at least one microseismic signal, the method further comprises:
judging the uniqueness of the positioning results of the JM method, the IJM method and the TM method according to the arrival time information division of the microseismic signals and the positions of the sensors, and if the positioning results of the positioning methods are not unique, taking the CM as a target function.
In the above embodiments, as shown in fig. 1-6 and 1-7, the uniqueness of the positioning solutions of the JM method, IJM method, TM method and LM method is determined according to the position relationship between the sensor array and the seismic source, as shown in table 1, and if there are multiple solutions in the JM method, IJM method, TM method and LM method, the CM method is used for positioning.
In some embodiments, said calculating at least one positioning error corresponding to said positioning function from at least one positioning function and at least one arrival time relationship of said microseismic signals comprises:
calculating a first positioning error w by JMJ；
Calculating a second positioning error w by IJMI；
Calculating a third positioning error w by TMT；
Calculating the fourth positioning error w by LM methodL。
In the above embodiment, the objective function of the JM method is:
Figure BDA0002669052020000081
wherein the content of the first and second substances,
Figure BDA0002669052020000082
v is the propagation velocity of the P wave, t0The occurrence time n of the seismic source is the number of the sensors, (x)0,y0,z0) As the source coordinates, each monitor sensor coordinates are (x)i,yi,zi)(i＝1,2,...,n)，Ti(i ═ 1, 2.. times.n) is the i-th sensor, ti(i ═ 1, 2., n) are observed times for each sensor.
The objective function of the IJM method is:
Figure BDA0002669052020000091
wherein, the theoretical arrival time difference of the sensor i and the sensor j (i ≠ j) is as follows:
Figure BDA0002669052020000092
the observed time differences for sensor i and sensor j (i ≠ j) are:
Δtij＝ti-tj (5)
li( i 1, 2.. n.) is the distance of each sensor to the seismic source.
The objective function of the TM method is:
Figure BDA0002669052020000093
wherein the content of the first and second substances,
Δfij＝li(tj-t0)-lj(ti-t0) (7)
the objective function of the LM method is:
Figure BDA0002669052020000094
wherein:
Δfijmk＝(li-lj)(tm-tk)-(ti-tj)(lm-lk) (9)
the objective function of the CM method is:
Figure BDA0002669052020000095
adopting JM method, IJM method, TM method and LM method to carry out positioning, and obtaining the positioning error of each algorithm as wJ，wI，wT，wL. Judging the convergence of the positioning results of the JM method, the IJM method, the TM method, the LM method and the CM method according to the time-of-arrival division condition of the sensor, and triggering the convergence of the time-of-arrival division condition and the solution of the sensor by each algorithm, such asFIGS. 1-5, schematic plan views of sensors in relation to the location of a seismic source, such as FIGS. 1-7.
In some embodiments, said repeatedly calculating a plurality of wave velocities corresponding to each of said positioning errors to obtain at least one wave velocity set corresponding to said positioning error comprises:
by Δ w ═ xn-wiI calculate the wave velocity error, where wiIs wJ、wI、wTAnd wLAny one of (1), xnPositioning error of CM method;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset lower limit value of the wave speed, and obtaining at least one lower limit value of the wave speed through mode search;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset upper limit value of the wave speed, and obtaining an upper limit value corresponding to at least one lower limit value of the wave speed through the mode search;
and obtaining the at least one wave velocity set by combining the at least one wave velocity lower limit value and the upper limit value corresponding to the at least one wave velocity lower limit value.
In the above embodiment, as shown in fig. 1 to 8, the positioning error of the CM method is calculated to reach wJThe wave velocity of the time is searched from v 3000m/s by using a pattern search method, the target function of the search is an equation (18), and the corresponding wave velocity when the equation (18) is 0 is vimin. The actual seismic source of the blasting is (x'0,y'0,z'0) The position of the micro seismic source is (x)0,y0,z0)。
Figure BDA0002669052020000101
Then, a pattern search method is used to start a search from v 9000m/s, and when equation (18) is 0, the corresponding wave velocity is defined as vimax(ii) a Obtaining the wave velocity range (v) required by the CM method that the positioning error is smaller than that of the JM methodimin,vimax) Put the wave velocity range into the wave velocity set VJIn (1).
In some embodiments, the calculating a range of each of the wave velocity sets, and taking the wave velocity set with the smallest range as a target set, includes:
obtaining the wJCorresponding first set of wave velocities (v)Jmin,vJmax) W toICorresponding first set of wave velocities (v)Imin,vImax) W toTCorresponding first set of wave velocities (v)Tmin,vTmax) And said wLCorresponding first set of wave velocities (v)Lmin,vLmax)；
Through min (v)Jmax-vJmin，vImax-vImin，vTmax-vTmin，vLmax-vLmin) Calculating the minimum value of the wave speed range in the at least one wave speed set;
and taking the wave velocity set corresponding to the minimum value of the wave velocity range as a target set.
In the above embodiment, as shown in FIG. 1-2, 4 shots were performed in the Q1 region, 5 shots were performed in the Q2 region, the CM method with a wave velocity of 5544m/s was used in the Q1 region, and the TM method was used in the Q2 region, and the positioning was performed based on the particle swarm optimization. Particle swarm algorithm parameter setting: learning factor c1＝c22.05, inertial weight w 0.5, population size N pop4000, flight number N g5000, adaptation value condition e 1.0 × 10-50The coordinate range is set as x e (1500,2200), y e (7500,9000), z e (-500,100) in m, v e (3000,9000) in m/s, t e (t)min-0.5,tmin) Wherein t isminThe minimum arrival time of the trigger event is in the unit of s, and the positioning result is shown in fig. 1-5.
In some embodiments, said calculating at least one positioning error corresponding to said positioning function from at least one positioning function and at least one arrival time information of said microseismic signals comprises:
by passing
Figure BDA0002669052020000111
Calculating at least one positioning error corresponding to said positioning function, wherein (x)1,y1,z1) Is composed ofThe source coordinates, (x)0,y0,z0) The microseismic signal coordinates.
In the above embodiment, as shown in fig. 1-2, the localization objective functions in the Q1 region and the Q2 region can be obtained by analyzing 2 blasting test events by using the adaptive identification algorithm of the microseismic source localization objective function, and the results are: the Q1 area should be positioned by CM method, and the optimal uniform wave velocity should be 5544 m/s; the Q2 region should be located by selecting TM method, and the relationship between the location results and wave velocity of CM method, JM method, IJM method, TM method and LM method is shown in FIGS. 1-6.
In some embodiments, prior to the acquiring time-of-arrival information for the seismic source and the at least one microseismic signal, the method further comprises:
and acquiring the coordinates of the seismic source through a particle swarm algorithm.
In the above embodiment, the blasting time is as shown in fig. 1 to 9, and after the microseismic event is triggered, the target function of each partition is adopted to perform positioning based on the particle swarm optimization. Particle swarm algorithm parameter setting: learning factor c1＝c22.05, inertial weight w 0.5, population size N pop4000, flight number N g5000, adaptation value condition e 1.0 × 10-50The search range setting is determined according to the monitoring area, v is equal to (3000,9000) and is in the unit of m/s and t is equal to (t)min-0.5,tmin) Wherein t isminThe minimum arrival time of the trigger event is in units of s.
Fig. 2 shows a schematic diagram of a seismic source locating apparatus 20, which can be used for seismic source location. The apparatus for source localization in an embodiment of the present application is capable of implementing the steps of the method for source localization corresponding to the embodiments corresponding to fig. 1-1 described above. The functions performed by the seismic source locating apparatus 20 may be implemented by hardware, or may be implemented by hardware executing corresponding software. The hardware or software includes one or more modules corresponding to the above functions, which may be software and/or hardware. The device for positioning the seismic source may include an input/output module 201 and a processing module 202, and the processing module 202 and the input/output module 201 may refer to operations executed in the embodiment corresponding to fig. 1-1, which are not described herein again. The input-output module 201 may be used to control input, output, and acquisition operations of the input-output module 201.
In some embodiments, the input-output module 201 may be configured to obtain arrival time information of the seismic source and the at least one microseismic signal;
the processing module 202 may be configured to calculate at least one positioning error corresponding to at least one positioning function from the at least one positioning function and the at least one arrival time information of the microseismic signals; repeatedly calculating a plurality of corresponding wave velocities through each positioning error to obtain at least one wave velocity set corresponding to the positioning error; calculating the range of each wave velocity set, and taking the wave velocity set with the minimum range as a target set; if the target set is an empty set, selecting the corresponding positioning function with the minimum positioning error as a target function; if the target set is not an empty set, taking a CM (classical method) as a target function; the source location is calculated by the objective function.
In some embodiments, the processing module 202 is further configured to:
judging the uniqueness of the positioning results of the joint method JM, the improved joint method IJM and the time method TM according to the arrival time information division of the microseismic signals and the positions of the sensors, and taking the CM as a target function if the positioning results of the positioning methods are not unique.
In some embodiments, the processing module 202 is further configured to:
calculating a first positioning error w by JMJ；
Calculating a second positioning error w by IJMI；
Calculating a third positioning error w by TMT；
Calculating the fourth positioning error w by LM methodL。
In some embodiments, the processing module 202 is further configured to:
by Δ w ═ xn-wiI calculate the wave velocity error, where wiIs wJ、wI、wTAnd wLIn (1)One term, xnPositioning error of CM method;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset lower limit value of the wave speed, and obtaining at least one lower limit value of the wave speed through mode search;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset upper limit value of the wave speed, and obtaining an upper limit value corresponding to at least one lower limit value of the wave speed through the mode search;
and obtaining the at least one wave velocity set by combining the at least one wave velocity lower limit value and the upper limit value corresponding to the at least one wave velocity lower limit value.
In some embodiments, the processing module 202 is further configured to:
obtaining the wJCorresponding first set of wave velocities (v)Jmin,vJmax) W toICorresponding first set of wave velocities (v)Imin,vImax) W toTCorresponding first set of wave velocities (v)Tmin,vTmax) And said wLCorresponding first set of wave velocities (v)Lmin,vLmax)；
Through min (v)Jmax-vJmin，vImax-vImin，vTmax-vTmin，vLmax-vLmin) Calculating the minimum value of the wave speed range in the at least one wave speed set;
and taking the wave velocity set corresponding to the minimum value of the wave velocity range as a target set.
In some embodiments, the processing module 202 is further configured to:
by passing
Figure BDA0002669052020000131
Calculating at least one positioning error corresponding to said positioning function, wherein (x)1,y1,z1) Is the source coordinate, (x)0,y0,z0) The microseismic signal coordinates.
In some embodiments, the processing module 202 is further configured to:
and acquiring the coordinates of the seismic source through a particle swarm algorithm.
The creating apparatus in the embodiment of the present application is described above from the perspective of the modular functional entity, and the following describes a computer device from the perspective of hardware, as shown in fig. 3, which includes: a processor, a memory, an input-output unit (which may also be a transceiver, not identified in fig. 3), and a computer program stored in the memory and executable on the processor. For example, the computer program may be a program corresponding to the method of seismic source localization in the embodiment corresponding to fig. 1-1. For example, when the computer device implements the functions of the apparatus for source localization 20 as shown in fig. 2, the processor, when executing the computer program, implements the steps of the method for source localization performed by the apparatus for source localization 20 in the embodiment corresponding to fig. 2 described above. Alternatively, the processor, when executing the computer program, implements the functions of the modules in the apparatus 20 for seismic source localization according to the embodiment corresponding to fig. 2. Also for example, the computer program may be a program corresponding to the method of seismic source localization in the embodiment corresponding to fig. 1-1.
The Processor 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. The general purpose processor may be a microprocessor or the processor may be any conventional processor or the like which is the control center for the computer device and which connects the various parts of the overall computer device using various interfaces and lines.
The memory may be used to store the computer programs and/or modules, and the processor may implement various functions of the computer device by running or executing the computer programs and/or modules stored in the memory and invoking data stored in the memory. The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, video data, etc.) created according to the use of the cellular phone, etc. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
The input-output unit may also be replaced by a receiver and a transmitter, which may be the same or different physical entities. When they are the same physical entity, they may be collectively referred to as an input-output unit. The input and output may be a transceiver.
The memory may be integrated in the processor or may be provided separately from the processor.
Through the above description of the embodiments, those skilled in the art will clearly understand that the method of the above embodiments can be implemented by software plus a necessary general hardware platform, and certainly can also be implemented by hardware, but in many cases, the former is a better implementation manner. Based on such understanding, the technical solutions of the present application may be embodied in the form of a software product, which is stored in a storage medium (e.g., ROM/RAM), and includes several instructions for enabling a terminal (e.g., a mobile phone, a computer, a server, or a network device) to execute the method according to the embodiments of the present application.
The embodiments of the present application have been described above with reference to the drawings, but the present application is not limited to the above-mentioned embodiments, which are only illustrative and not restrictive, and those skilled in the art can make many changes and modifications without departing from the spirit and scope of the present application and the protection scope of the claims, and all changes and modifications that come within the meaning and range of equivalency of the claims are to be embraced within their scope.

Claims (5)

1. A method for seismic source localization, the method comprising:
acquiring time-of-arrival information of a seismic source and at least one microseismic signal;
calculating at least one positioning error corresponding to the positioning function through at least one positioning function and at least one arrival time information of the microseismic signals;
repeatedly calculating a plurality of corresponding wave velocities through each positioning error to obtain at least one wave velocity set corresponding to the positioning error;
calculating the range of each wave velocity set, and taking the wave velocity set with the minimum range as a target set;
if the target set is an empty set, selecting the corresponding positioning function with the minimum positioning error as a target function;
if the target set is not an empty set, taking a CM (classical method) as a target function;
calculating a seismic source position through the objective function;
after the time of arrival information of the seismic source and the at least one microseismic signal is obtained, the method further comprises:
judging the uniqueness of the positioning results of the combination method JM, the improved combination method UM and the time method TM according to the arrival time information division of the microseismic signals and the positions of the sensors, and taking the CM as a target function if the positioning results of the positioning methods are not unique;
calculating at least one positioning error corresponding to the positioning function by at least one positioning function and arrival time information of the at least one microseismic signal, comprising:
calculating a first positioning error w by JMJ；
Calculating a second positioning error w by IJMI；
Calculating a third positioning error w by TMT；
By passingLM method for calculating fourth positioning error wL；
Repeatedly calculating a plurality of corresponding wave velocities through each positioning error to obtain at least one wave velocity set corresponding to the positioning error, wherein the wave velocity set comprises:
by Δ w ═ xn-wiI calculate the wave velocity error, where wiIs wJ、wI、wTAnd wLAny one of (1), xnPositioning error of CM method;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset lower limit value of the wave speed, and obtaining at least one lower limit value of the wave speed through mode search;
when the delta w is equal to 0, enabling the wave speed to be equal to a preset upper limit value of the wave speed, and obtaining an upper limit value corresponding to at least one lower limit value of the wave speed through the mode search;
obtaining the at least one wave velocity set by combining the at least one wave velocity lower limit value and an upper limit value corresponding to the at least one wave velocity lower limit value;
the calculating the range of each wave velocity set, and taking the wave velocity set with the minimum range as a target set, includes:
obtaining the wJCorresponding first set of wave velocities (v)Jmin，vJmax) W toICorresponding first set of wave velocities (v)Imin，vImax) W toTCorresponding first set of wave velocities (v)Tmin，vTmax) And said wLCorresponding first set of wave velocities (v)Lmin，vLmax)；
Through min (v)Jmax-vJmin，vImax-vImin，vTmax-vTmin，vLmax-vLmin) Calculating the minimum value of the wave speed range in the at least one wave speed set;
and taking the wave velocity set corresponding to the minimum value of the wave velocity range as a target set.
2. The method of claim 1, wherein said calculating at least one positioning error corresponding to said positioning function from at least one positioning function and at least one arrival time information of said microseismic signals comprises:
by passing
Figure FDA0003226437020000021
Calculating at least one positioning error corresponding to said positioning function, wherein (x)1，y1，z1) Is the source coordinate, (x)0，y0，z0) The microseismic signal coordinates.
3. The method of any of claims 1-2, wherein prior to the obtaining time-of-arrival information for the seismic source and the at least one microseismic signal, the method further comprises:
and acquiring the coordinates of the seismic source through a particle swarm algorithm.
4. A computer device, characterized in that the computer device comprises:
at least one processor, a memory, and an input-output unit;
wherein the memory is configured to store program code and the processor is configured to invoke the program code stored in the memory to perform the method of any of claims 1-2.
5. A computer storage medium comprising instructions which, when executed on a computer, cause the computer to perform the method of any one of claims 1-2.
CN202010927795.6A 2020-09-07 2020-09-07 Seismic source positioning method, device, equipment and storage medium Active CN111880220B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010927795.6A CN111880220B (en) 2020-09-07 2020-09-07 Seismic source positioning method, device, equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010927795.6A CN111880220B (en) 2020-09-07 2020-09-07 Seismic source positioning method, device, equipment and storage medium

Publications (2)

Publication Number Publication Date
CN111880220A CN111880220A (en) 2020-11-03
CN111880220B true CN111880220B (en) 2022-02-22

Family

ID=73199552

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010927795.6A Active CN111880220B (en) 2020-09-07 2020-09-07 Seismic source positioning method, device, equipment and storage medium

Country Status (1)

Country Link
CN (1) CN111880220B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112731523B (en) * 2020-12-22 2022-09-16 北京环境特性研究所 Seismic source positioning method, computer equipment and computer readable storage medium
CN112946735B (en) * 2021-02-23 2023-03-31 石家庄铁道大学 Rockfall impact positioning method and device based on micro-seismic monitoring system
CN116243379B (en) * 2023-02-28 2023-11-10 中国矿业大学 Strong mineral earthquake prediction method based on earthquake focus mechanism and positioning error calibration

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8902707B2 (en) * 2007-04-09 2014-12-02 Baker Hughes Incorporated Analysis of uncertainty of hypocenter location using the combination of a VSP and a subsurface array
CN106680867B (en) * 2016-11-17 2018-07-24 大连理工大学 A kind of pinpoint dynamic parameter method of microseismic event
CN106990435B (en) * 2017-06-07 2018-11-23 中煤科工集团西安研究院有限公司 It is a kind of to weaken the microseism localization method and device for relying on first break pickup precision
CN107884822B (en) * 2017-11-13 2019-09-27 北京矿冶研究总院 Method for improving positioning precision of mining micro-seismic source
CN108717201B (en) * 2018-06-20 2019-10-25 成都理工大学 A kind of tunnel surrounding microquake sources localization method
CN109738940B (en) * 2019-01-18 2021-01-29 东北大学 Acoustic emission/microseismic event positioning method under condition of existing empty zone
CN110609321B (en) * 2019-09-24 2020-11-13 中国科学院武汉岩土力学研究所 Micro seismic source positioning method based on speed model database

Also Published As

Publication number Publication date
CN111880220A (en) 2020-11-03

Similar Documents

Publication Publication Date Title
CN111880220B (en) Seismic source positioning method, device, equipment and storage medium
Raykar et al. Position calibration of microphones and loudspeakers in distributed computing platforms
Ajdler et al. Acoustic source localization in distributed sensor networks
CN107968987B (en) RSSI weighted centroid positioning method based on fixed integral combined with environmental parameters
CN112245923A (en) Collision detection method and device in game scene
CN115048824B (en) Collision detection method and device and computer readable medium
Zhou et al. A new acoustic emission source location method using tri-variate kernel density estimator
Patel et al. Modeling and analysis of a manufacturing system with deadlocks to generate the reachability tree using petri net system
Machado et al. A generic approach to build plant models for DES verification purposes
CN109410198A (en) A kind of timing motion detection method, device and equipment
Singh et al. Node localization in wireless sensor networks using the M5P tree and SMOreg algorithms
Lee et al. Uncertainty self-management with perception net based geometric data fusion
US10331803B2 (en) Computer aided modeling
US11079468B2 (en) Detection of acoustic events
CN112584985A (en) Simulated signature key for robot simulation
CN113552536B (en) Acoustic emission/microseismic event positioning method, system, terminal and readable storage medium containing round hole structure
US10406439B1 (en) Method and system for facilitating generation of a virtual world including one or more objects and one or more agents based on a user intention
JP2010152863A (en) System and method for extracting boundary node, and program
Wendeberg et al. Polynomial time approximation algorithms for localization based on unknown signals
CN111679298A (en) Integrity monitoring method and integrity monitoring device of navigation system and electronic equipment
Yairi et al. Fast and simple topological map construction based on cooccurrence frequency of landmark observation
Andresen Target Impact Location Determination: An Analytic Method
US20240005049A1 (en) Systems and methods for generation of exploded views of computer-aided design models
US11089037B2 (en) Method for automatically retrieving workflow and method for detecting intrusion in distributed system
KR101606283B1 (en) Appratus for automatic generation of expected results for testing component-based software system

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
TR01 Transfer of patent right

Effective date of registration: 20230914

Address after: 430071 small Hongshan 2, fruit lake street, Wuchang District, Wuhan, Hubei.

Patentee after: INSTITUTE OF ROCK AND SOIL MECHANICS, CHINESE ACADEMY OF SCIENCES

Patentee after: Sichuan Tibet Railway Co.,Ltd.

Patentee after: China National Railway Group Co.,Ltd.

Address before: 430071 small Hongshan 2, fruit lake street, Wuchang District, Wuhan, Hubei.

Patentee before: INSTITUTE OF ROCK AND SOIL MECHANICS, CHINESE ACADEMY OF SCIENCES

Patentee before: Sichuan Tibet Railway Co.,Ltd.

TR01 Transfer of patent right