CN112904277B - Tunnel surrounding rock breaking point positioning method based on improved gray wolf algorithm - Google Patents

Tunnel surrounding rock breaking point positioning method based on improved gray wolf algorithm Download PDF

Info

Publication number
CN112904277B
CN112904277B CN202110096316.5A CN202110096316A CN112904277B CN 112904277 B CN112904277 B CN 112904277B CN 202110096316 A CN202110096316 A CN 202110096316A CN 112904277 B CN112904277 B CN 112904277B
Authority
CN
China
Prior art keywords
space
point
coordinates
acoustic emission
fitness function
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
CN202110096316.5A
Other languages
Chinese (zh)
Other versions
CN112904277A (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.)
Chongqing City Construction Investment Group Co ltd
China Merchants Chongqing Communications Research and Design Institute Co Ltd
Third Engineering Co Ltd of China Railway Construction Bridge Engineering Bureau Group Co Ltd
Original Assignee
Chongqing City Construction Investment Group Co ltd
China Merchants Chongqing Communications Research and Design Institute Co Ltd
Third Engineering Co Ltd of China Railway Construction Bridge Engineering Bureau Group Co Ltd
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 Chongqing City Construction Investment Group Co ltd, China Merchants Chongqing Communications Research and Design Institute Co Ltd, Third Engineering Co Ltd of China Railway Construction Bridge Engineering Bureau Group Co Ltd filed Critical Chongqing City Construction Investment Group Co ltd
Priority to CN202110096316.5A priority Critical patent/CN112904277B/en
Publication of CN112904277A publication Critical patent/CN112904277A/en
Application granted granted Critical
Publication of CN112904277B publication Critical patent/CN112904277B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/18Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves

Abstract

The invention discloses a tunnel surrounding rock breaking point positioning method based on an improved wolf algorithm, aiming at the defects of lower efficiency and poorer precision of the traditional acoustic emission positioning algorithm, the invention provides a method for searching the space coordinates of the surrounding rock breaking point by taking the sum of residual absolute values of travel time equations as an adaptability function of the wolf algorithm on the basis of acoustic emission monitoring data, and simultaneously, in order to overcome the problem of poor local searching capability of the wolf algorithm, a population memory elimination mechanism is introduced to improve the wolf algorithm, and a cosine function decreasing form is adopted to improve convergence factors. The method and the system have the advantages that the calculation efficiency, accuracy and stability of the positioning algorithm are better improved, the occurrence position of cracking and breaking of surrounding rock in the tunnel excavation process is displayed in an intuitive and image mode, and a rapid and accurate real-time monitoring and early warning means is provided for tunnel construction.

Description

Tunnel surrounding rock breaking point positioning method based on improved gray wolf algorithm
Technical Field
The invention relates to the technical field of tunnel safety monitoring, in particular to a tunnel surrounding rock breaking point positioning method based on an improved gray wolf algorithm.
Background
In the tunnel excavation process, serious personal and property losses can be caused by instability damage of surrounding rocks, so that the method has important engineering significance for monitoring the cracking of the surrounding rocks in real time in the construction process, and can provide basis for disaster prevention and control, thereby avoiding the occurrence of disasters.
The acoustic emission monitoring technology not only can detect the generation of microcracks in the rock in real time, but also can realize more accurate positioning, and is widely applied to geotechnical engineering geological disaster prevention and control. However, the conventional acoustic emission positioning algorithm has respective limitations, low efficiency and poor precision, and needs to be improved.
With the rapid development of computer technology, more and more intelligent algorithms are developed successively, and the intelligent algorithms are widely focused by students once coming out due to the advantages of strong searching capability, rapid convergence, simple structure, easy realization and the like, and are successfully applied to the fields of parameter optimization, image processing and the like, but are still freshly applied to the field of acoustic emission positioning.
In conclusion, the traditional acoustic emission positioning method is improved by using the intelligent algorithm, so that the acoustic emission positioning of the tunnel surrounding rock fracture can be realized more quickly and accurately, and more reliable monitoring data is provided for the tunnel engineering disaster prevention and control work.
Disclosure of Invention
Aiming at the defects existing in the prior art, the invention provides a tunnel surrounding rock breaking point positioning method based on an improved gray wolf algorithm, the position of a surrounding rock breaking point is positioned by adopting an intelligent algorithm through acoustic emission monitoring data, and the efficiency, accuracy and stability of positioning are improved.
The specific technical scheme is as follows:
in a first aspect, a method for positioning a breaking point of a surrounding rock of a tunnel based on an improved wolf algorithm is provided, which comprises the following steps:
collecting P-wave signals generated when surrounding rock of a tunnel is broken through a plurality of acoustic emission probes which are uniformly distributed, and recording the receiving time of each acoustic emission probe for receiving the P-wave signals;
and positioning the breaking point coordinates corresponding to the breaking points through a gray wolf algorithm based on the space coordinates of each acoustic emission probe in the positioning space and the receiving time of the received P-wave signals.
With reference to the first aspect, in a first implementation manner of the first aspect, the locating, by a gray wolf algorithm, the rupture point coordinates corresponding to the rupture point includes:
randomly selecting a plurality of space points in a positioning space, and determining the coordinates of each space point;
calculating the fitness function value of each space point through the coordinates of the space point, the space coordinates corresponding to each acoustic emission probe and the receiving time;
selecting m optimal space points from all the space points according to the fitness function value;
updating the fitness function value of each space point based on the coordinates of the optimal space point;
judging whether the ending condition is met or not, if not, re-selecting the fitness function value of each space point of the m optimal space points, iterating in the way until the ending condition is met, and selecting the coordinate corresponding to the minimum fitness function value as the coordinate of the breaking point.
With reference to the first implementation manner of the first aspect, in a second implementation manner of the first aspect, the calculating an fitness function value of each spatial point includes:
calculating travel time equation residual errors of the space points relative to each acoustic emission probe according to the coordinates of the space points, the space coordinates corresponding to the acoustic emission probes and the receiving time;
and determining the fitness function value of the space point according to the travel time equation residual error of the space point relative to each acoustic emission probe.
With reference to the first implementation manner of the first aspect, in a third implementation manner of the first aspect, the following method is adopted to select m optimal spatial points:
sequencing all the fitness function values obtained by calculation according to the sequence from small to large;
and selecting the space points corresponding to the m fitness function values as optimal space points according to the arrangement sequence.
With reference to the first implementation manner of the first aspect, in a fourth implementation manner of the first aspect, the fitness function value of the spatial point is updated by adopting the following method:
calculating the position distance between the space point and each optimal space point and the convergence factor corresponding to the current iteration step number;
determining candidate position coordinates of the space points according to the convergence factor and the position distance;
updating coordinates of the space points according to the candidate position coordinates;
and calculating to obtain the fitness function value of the updated space point through the updated space point coordinates, the space coordinates corresponding to each acoustic emission probe and the receiving time.
With reference to the fourth implementation manner of the first aspect, in a fifth implementation manner of the first aspect, based on a preset total iteration step number and a current iteration step number, determining a convergence factor corresponding to the current iteration step number by using a cosine function.
With reference to the first aspect, in a sixth implementation manner of the first aspect, after updating the fitness function value of each spatial point, the method further includes: and performing space point replacement through a population memory elimination mechanism.
With reference to the first aspect, in a seventh implementation manner of the first aspect, the method further includes:
recording the breaking point coordinates of each determined breaking point;
and drawing a three-dimensional image according to the recorded data and the space coordinates of all the acoustic emission probes.
In a second aspect, a tunnel surrounding rock breaking point positioning method based on an improved wolf algorithm is provided, which comprises the following steps:
a plurality of acoustic emission probes are uniformly arranged on the inner wall of the surrounding rock of the tunnel;
determining the space coordinates of each acoustic emission probe in the positioning space;
recording the moment when each acoustic emission probe receives a P-wave signal generated by surrounding rock fracture;
and positioning the breaking point coordinates corresponding to the breaking points through a gray wolf algorithm based on the space coordinates of each acoustic emission probe in the positioning space and the receiving time of the received P-wave signals.
With reference to the second aspect, in a first implementation manner of the second aspect, the determining spatial coordinates of each acoustic emission probe in the positioning space includes:
constructing a three-dimensional rectangular coordinate system corresponding to the positioning space by taking the center of the tunnel face as an origin;
and determining the space coordinates of the acoustic emission probes according to the relative positions of each acoustic emission probe and the center of the face.
The beneficial effects are that:
(1) The gray wolf algorithm in the intelligent algorithm is utilized, acoustic emission monitoring data are combined, and the sum of residual absolute values of travel time equations is used as an adaptability function of the gray wolf algorithm to search space coordinates of surrounding rock fracture points, so that the calculation efficiency, accuracy and stability of the traditional positioning algorithm are improved.
(2) The linear decreasing form of the convergence factor in the traditional wolf algorithm is changed into the decreasing form of the cosine function, so that the wolf algorithm has higher global searching capability in the early stage and higher local searching capability in the later stage, and is more in line with the actual searching process.
(3) In order to solve the problem of poor local searching capability of the wolf algorithm, a population memory elimination mechanism is introduced to improve the wolf algorithm, the optimal individuals of the population are always reserved in the iteration process, and the worst individuals are eliminated, so that the local searching capability is improved, and the searching efficiency is also accelerated.
(4) According to the optimal space point position coordinates obtained by searching the improved gray wolf algorithm, a three-dimensional acoustic emission positioning schematic diagram of the surrounding rock fracture of the tunnel is drawn in real time, the occurrence position of the surrounding rock fracture in the tunnel excavation process can be shown for engineering construction personnel intuitively and vividly, and the effect of real-time monitoring and early warning is achieved.
Drawings
In order to more clearly illustrate the embodiments of the present invention, the drawings that are required to be used in the embodiments will be briefly described. Throughout the drawings, the elements or portions are not necessarily drawn to actual scale.
FIG. 1 is a flow chart of a method for locating a breaking point according to an embodiment of the present invention;
FIG. 2 is a flow chart of the gray wolf algorithm shown in FIG. 1 for determining the coordinates of a rupture point;
FIG. 3 is a flow chart of updating fitness function values of other spatial points shown in FIG. 2;
FIG. 4 is a flow chart of determining a rupture point according to the gray wolf algorithm provided by an embodiment of the present invention;
FIG. 5 is a flow chart of a method for locating a rupture point according to an embodiment of the present invention;
fig. 6 is a flowchart of a method for locating a breaking point according to an embodiment of the present invention.
Detailed Description
Embodiments of the technical scheme of the present invention will be described in detail below with reference to the accompanying drawings. The following examples are only for more clearly illustrating the technical aspects of the present invention, and thus are merely examples, and are not intended to limit the scope of the present invention.
The method for locating the breaking point of the surrounding rock of the tunnel based on the improved gray wolf algorithm as shown in fig. 1 comprises the following steps:
collecting P-wave signals generated when surrounding rock of a tunnel is broken through a plurality of acoustic emission probes which are uniformly distributed, and recording the moment when each acoustic emission probe receives the P-wave signals;
and step two, positioning the coordinates of the breaking points corresponding to the breaking points through a gray wolf algorithm based on the space coordinates of each acoustic emission probe in the positioning space and the receiving time of the received P-wave signals.
Specifically, when a local crack of each place of the surrounding rock of the tunnel is broken, the acoustic emission probes receive a P-wave signal in real time, the space coordinates of each acoustic emission probe and the time when the P-wave signal is received are output to an external computer, and the computer can position the breaking point coordinates of the breaking points by adopting the gray wolf algorithm in the existing intelligent algorithm according to the space coordinates of each acoustic emission probe and the time when the P-wave signal is received, so that the calculation efficiency, accuracy and stability of the positioning algorithm are improved.
In this embodiment, preferably, as shown in fig. 2, the locating, by a wolf algorithm, the fracture point coordinates corresponding to the fracture point includes:
step 1, randomly selecting a plurality of space points in a positioning space, and determining the coordinates of each space point;
step 2, calculating the fitness function value of each space point through the coordinates of the space point, the space coordinates corresponding to each acoustic emission probe and the receiving time;
step 3, selecting m optimal space points from all space points according to the fitness function value;
step 4, updating the fitness function value of each space point based on the coordinates of the optimal space point;
and 5, judging whether the ending condition is met, if not, re-selecting the fitness function value of each space point of the m optimal space points, iterating in the way until the ending condition is met, and selecting the coordinate corresponding to the minimum fitness function value as the coordinate of the breaking point.
In particular, the method comprises the steps of,
firstly, the computer can randomly initialize three-dimensional position coordinates of 100 space points, namely, initialize a gray wolf population, in the monitoring range of each acoustic emission probe, namely, in a positioning space, by using a rand function, wherein the coordinates of each space point are x j ,y j ,z j ,j=1,2,3…,100;
And secondly, calculating the fitness function value of each space point, namely the fitness of each individual wolf by the computer through the corresponding coordinates, the space coordinates corresponding to each acoustic emission probe and the receiving time.
In this embodiment, in order to improve the accuracy and stability of positioning calculation, the travel time equation residual error of the space point relative to each acoustic emission probe is calculated according to the coordinates of the space point, the space coordinates corresponding to the acoustic emission probes, and the receiving time, and the fitness function value of the space point is determined according to the travel time equation residual error of the space point relative to each acoustic emission probe, where the specific calculation formula is as follows:
Figure BDA0002914371100000061
wherein X is i ,Y i ,Z i For the spatial position coordinates of the acoustic emission probe, i=1, 2,3 …, N, t i For the receiving moment of the P wave signal received by the acoustic emission probe, vp is the P wave speed, t 1 The receiving moment of the acoustic emission probe which receives the P-wave signal first.
Then, the computer selects 3 (m=3) space points with the smallest fitness function value as the optimal space points by comparing the fitness function values of the space points, namely 3 head wolves alpha, beta and delta.
In this embodiment, in order to accelerate the computing efficiency of the computer, the following method is adopted to select the optimal spatial point:
sequencing all the fitness function values obtained by calculation according to the sequence from small to large;
and selecting the space points corresponding to the m fitness function values as optimal space points according to the arrangement sequence.
Specifically, the computer uses the sort function to sequence all fitness function values from small to large, and then selects the spatial point corresponding to the fitness function value of the first 3 bits from the selected spatial point as the optimal spatial point according to the sequence.
Then, as shown in fig. 3, the computer updates the fitness function value of the spatial point by the following method:
step 4-1, calculating the position distance between the space point and each optimal space point and the convergence factor corresponding to the current iteration step number;
step 4-2, determining candidate position coordinates of the space points according to the convergence factor and the position distance;
step 4-3, updating coordinates of the space points according to the candidate position coordinates;
and 4-4, calculating to obtain the fitness function value of the updated space point through the updated space point coordinates, the space coordinates corresponding to each acoustic emission probe and the receiving time.
Specifically, first, the distances between each other spatial point, i.e., candidate wolf, and 3 optimal spatial points α, β, δ are calculated as follows:
Figure BDA0002914371100000071
Figure BDA0002914371100000072
Figure BDA0002914371100000073
r is a random number within the range of 0, 1.
Then according to the distance from the space point to 3 optimal space points alpha, beta and delta, 3 candidate position coordinates of each candidate wolf are determined according to the following calculation formula:
Figure BDA0002914371100000074
Figure BDA0002914371100000075
Figure BDA0002914371100000076
a is a convergence factor decreasing from 2 to 0, in this embodiment, in order to enable the positioning method to have higher global searching capability in the early stage and higher local searching capability in the later stage, so as to be more suitable for an actual searching process, in this embodiment, a cosine function is used to determine the convergence factor corresponding to the current iteration step number, and a specific calculation formula is as follows:
Figure BDA0002914371100000077
where t is the current iteration step number and maxt is the total iteration step number.
After determining 3 candidate position coordinates corresponding to the space point, calculating coordinates after the space point update by adopting the following calculation formula:
Figure BDA0002914371100000078
Figure BDA0002914371100000079
Figure BDA0002914371100000081
and finally, calculating by using the coordinates updated by the space points, the space coordinates corresponding to each acoustic emission probe and the receiving time by using the calculation formula to obtain the fitness function value updated by the space points.
And finally, the computer judges whether the optimal space point, namely the fitness function value of the head wolf is smaller than a preset error or the current iteration number is larger than a preset total iteration number, if not, the step 3 is returned, and if so, the coordinate corresponding to the minimum fitness function value is recorded as the coordinate of the breaking point.
The second embodiment is substantially the same as the first embodiment, and the main difference is that: as shown in fig. 4, in this embodiment, preferably, after updating the fitness function value of each spatial point, the method further includes:
and 6, replacing space points through a population memory elimination mechanism.
Specifically, after the fitness function values of all the space points are updated, the computer ranks the space points by using the sort function, finds 3 space points with the largest fitness function values after updating, and replaces the space points with the largest fitness function values by using the 3 optimal space points before updating, namely the 3 top wolves. Therefore, optimal individuals of the population are always reserved in the iteration process, and the worst individuals are eliminated, so that the local searching capability is improved, and the searching efficiency is also improved.
In the third embodiment, as shown in fig. 5, in this embodiment, it is preferable that the method further includes:
and 7, recording the breaking point coordinates of the breaking points determined each time, and drawing a three-dimensional image according to the recorded data and the space coordinates of all the acoustic emission probes.
Specifically, after each time the computer determines a breaking point coordinate, a three-dimensional acoustic emission positioning schematic diagram of the breaking of the surrounding rock of the tunnel can be drawn in real time according to the breaking point coordinate and the space coordinates of all acoustic emission probes, so that the occurrence position of the breaking of the surrounding rock in the tunnel excavation process can be shown for engineering construction personnel more intuitively and vividly, and the effect of real-time monitoring and early warning is achieved.
The method for locating the breaking point of the surrounding rock of the tunnel based on the improved gray wolf algorithm as shown in fig. 6 comprises the following steps:
s1, uniformly arranging a plurality of acoustic emission probes on the inner wall of surrounding rock of a tunnel;
s2, determining the space coordinates of each acoustic emission probe in a positioning space;
s3, recording the time when each acoustic emission probe receives a P-wave signal generated by surrounding rock fracture;
and S4, positioning the coordinates of the breaking points corresponding to the breaking points through a gray wolf algorithm based on the space coordinates of each acoustic emission probe and the receiving time of the received P-wave signals.
In particular, the method comprises the steps of,
firstly, a plurality of acoustic emission probes can be uniformly arranged on the inner wall of a tunnel manually;
secondly, constructing a three-dimensional rectangular coordinate system of a positioning space by taking the center of the tunnel face as an origin, establishing a tunnel three-dimensional model, and determining the space coordinate of each acoustic emission probe in the positioning space according to the actual arrangement position of the acoustic emission probe;
then, recording the moment when each acoustic emission probe receives a P-wave signal generated by surrounding rock fracture;
finally, based on the space coordinates of each acoustic emission probe and the receiving time of the received P-wave signal, the corresponding breaking point coordinates of the breaking point are positioned by adopting a gray-wolf algorithm, and the specific flow of positioning the breaking point coordinates by adopting the gray-wolf algorithm is the same as that of the positioning method, so that the detailed description is omitted.
The above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention, and are intended to be included within the scope of the appended claims and description.

Claims (9)

1. The tunnel surrounding rock breaking point positioning method based on the improved gray wolf algorithm is characterized by comprising the following steps of:
collecting P-wave signals generated when surrounding rock of a tunnel is broken through a plurality of acoustic emission probes which are uniformly distributed, and recording the receiving time of each acoustic emission probe for receiving the P-wave signals;
based on the space coordinates of each acoustic emission probe in the positioning space and the receiving time of the received P-wave signal, positioning the coordinates of the breaking point corresponding to the breaking point through a gray wolf algorithm;
the locating the corresponding breaking point coordinates of the breaking point through the gray wolf algorithm comprises the following steps:
randomly selecting a plurality of space points in a positioning space, and determining the coordinates of each space point;
calculating the fitness function value of each space point through the coordinates of the space point, the space coordinates corresponding to each acoustic emission probe and the receiving time, wherein the specific calculation formula is as follows:
Figure FDA0004201193300000011
wherein X is i ,Y i ,Z i For spatial position coordinates of the acoustic emission probe, i=1, 2, 3..n, t i For the receiving moment of the P wave signal received by the acoustic emission probe, vp is the P wave speed, t 1 The receiving moment of the acoustic emission probe for receiving the P-wave signal is the first time;
selecting m optimal space points from all the space points according to the fitness function value;
updating the fitness function value of each space point based on the coordinates of the optimal space point;
judging whether the ending condition is met or not, if not, re-selecting the fitness function value of each space point of the m optimal space points, iterating in the way until the ending condition is met, and selecting the coordinate corresponding to the minimum fitness function value as the coordinate of the breaking point.
2. The tunnel surrounding rock fracture point positioning method based on the improved wolf algorithm of claim 1, wherein the calculating the fitness function value of each spatial point comprises:
calculating travel time equation residual errors of the space points relative to each acoustic emission probe according to the coordinates of the space points, the space coordinates corresponding to the acoustic emission probes and the receiving time;
and determining the fitness function value of the space point according to the travel time equation residual error of the space point relative to each acoustic emission probe.
3. The tunnel surrounding rock fracture point positioning method based on the improved wolf algorithm as claimed in claim 1, wherein the following method is adopted to select m optimal space points:
sequencing all the fitness function values obtained by calculation according to the sequence from small to large;
and selecting the space points corresponding to the m fitness function values as optimal space points according to the arrangement sequence.
4. The tunnel surrounding rock fracture point positioning method based on the improved wolf algorithm as claimed in claim 1, wherein the fitness function value of the space point is updated by adopting the following method:
calculating the position distance between the space point and each optimal space point and the convergence factor corresponding to the current iteration step number;
determining candidate position coordinates of the space points according to the convergence factor and the position distance;
updating coordinates of the space points according to the candidate position coordinates;
and calculating to obtain the fitness function value of the updated space point through the updated space point coordinates, the space coordinates corresponding to each acoustic emission probe and the receiving time.
5. The tunnel surrounding rock breaking point positioning method based on the improved wolf algorithm of claim 4, wherein the convergence factor corresponding to the current iteration step number is determined by adopting a cosine function based on the preset total iteration step number and the current iteration step number.
6. The tunnel surrounding rock breaking point positioning method based on the improved wolf algorithm according to claim 1, further comprising, after updating the fitness function value of each spatial point: and performing space point replacement through a population memory elimination mechanism.
7. The tunnel surrounding rock breaking point positioning method based on the improved wolf algorithm of claim 1, further comprising:
recording the breaking point coordinates of each determined breaking point;
and drawing a three-dimensional image according to the recorded data and the space coordinates of all the acoustic emission probes.
8. The tunnel surrounding rock breaking point positioning method based on the improved gray wolf algorithm is characterized by comprising the following steps of:
a plurality of acoustic emission probes are uniformly arranged on the inner wall of the surrounding rock of the tunnel;
determining the space coordinates of each acoustic emission probe in the positioning space;
recording the moment when each acoustic emission probe receives a P-wave signal generated by surrounding rock fracture;
based on the space coordinates of each acoustic emission probe and the receiving time of the received P-wave signal, locating the coordinates of the breaking point corresponding to the breaking point through a gray wolf algorithm;
the positioning of the rupture point coordinates corresponding to the rupture point by the gray wolf algorithm comprises:
randomly selecting a plurality of space points in a positioning space, and determining the coordinates of each space point;
calculating the fitness function value of each space point through the coordinates of the space point, the space coordinates corresponding to each acoustic emission probe and the receiving time;
selecting m optimal space points from all the space points according to the fitness function value;
and updating the fitness function value of each space point based on the coordinates of the optimal space point, wherein the specific calculation formula is as follows:
Figure FDA0004201193300000031
wherein X is i ,Y i ,Z i For spatial position coordinates of the acoustic emission probe, i=1, 2, 3..n, t i For the receiving moment of the P wave signal received by the acoustic emission probe, vp is the P wave speed, t 1 For the first sound of receiving P-wave signalsThe receiving moment of the transmitting probe;
judging whether the ending condition is met or not, if not, re-selecting the fitness function value of each space point of the m optimal space points, iterating in the way until the ending condition is met, and selecting the coordinate corresponding to the minimum fitness function value as the coordinate of the breaking point.
9. The method for locating a tunnel surrounding rock fracture point based on the improved wolf algorithm of claim 8, wherein the determining the spatial coordinates of each acoustic emission probe in the locating space comprises:
constructing a three-dimensional rectangular coordinate system corresponding to the positioning space by taking the center of the tunnel face as an origin;
determining the relative position of each acoustic emission probe and the center of the face;
and determining the space coordinates of each acoustic emission probe according to the corresponding relative positions.
CN202110096316.5A 2021-01-25 2021-01-25 Tunnel surrounding rock breaking point positioning method based on improved gray wolf algorithm Active CN112904277B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110096316.5A CN112904277B (en) 2021-01-25 2021-01-25 Tunnel surrounding rock breaking point positioning method based on improved gray wolf algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110096316.5A CN112904277B (en) 2021-01-25 2021-01-25 Tunnel surrounding rock breaking point positioning method based on improved gray wolf algorithm

Publications (2)

Publication Number Publication Date
CN112904277A CN112904277A (en) 2021-06-04
CN112904277B true CN112904277B (en) 2023-06-02

Family

ID=76118943

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110096316.5A Active CN112904277B (en) 2021-01-25 2021-01-25 Tunnel surrounding rock breaking point positioning method based on improved gray wolf algorithm

Country Status (1)

Country Link
CN (1) CN112904277B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113917287A (en) * 2021-11-22 2022-01-11 国家电网有限公司 Substation bus joint discharge heating fault monitoring system

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2580768A (en) * 1947-08-14 1952-01-01 Ibm Data look-up apparatus for computing or other machines
US4858462A (en) * 1989-01-20 1989-08-22 The Babcock & Wilcox Company Acoustic emission leak source location
JP2004219075A (en) * 2003-01-09 2004-08-05 Civil Engineering Research Institute Of Hokkaido Ae measurement body, function testing method for ae measurement body, ae measuring method, and active natural ground stabilization evaluating method
US8693286B1 (en) * 2009-10-16 2014-04-08 Snap-On Incorporated Position measurement for collision repair systems
CN108717201A (en) * 2018-06-20 2018-10-30 成都理工大学 A kind of tunnel surrounding microquake sources localization method
CN110345392A (en) * 2019-07-16 2019-10-18 辽宁石油化工大学 The localization method and device of oil pipeline leak source
CN110986747A (en) * 2019-12-20 2020-04-10 桂林电子科技大学 Landslide displacement combined prediction method and system

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2567307A (en) * 1946-06-12 1951-09-11 Western Union Telegraph Co System and apparatus involving optical scanning
US20050159403A1 (en) * 2003-04-22 2005-07-21 Pharmacia Corporation Compositions of a cyclooxygenase-2 selective inhibitor and a calcium modulating agent for the treatment of central nervous system damage
WO2010105136A2 (en) * 2009-03-13 2010-09-16 Cornell University Spoofing detection for civilian gnss signals
CN109283047B (en) * 2018-11-29 2023-10-20 四川大学 Rock mass damage monitoring system and evaluation method in deep engineering environment
CN109738519B (en) * 2019-01-04 2021-08-17 国网四川省电力公司广安供电公司 Denoising method for ultrasonic detection of lead of high-voltage bushing of transformer
CN110167138B (en) * 2019-05-23 2021-01-01 西安电子科技大学 Station distribution optimization method of passive time difference positioning system based on improved wolf optimization algorithm
CN110944342B (en) * 2019-10-24 2023-03-10 江西理工大学 Wireless sensor network deployment optimization method, device, system and storage medium
CN110913404B (en) * 2019-11-11 2022-05-31 沈阳理工大学 UWSNs node positioning method based on node movement prediction
CN110909947B (en) * 2019-11-29 2022-08-09 华中科技大学 Rectangular part layout method and equipment based on wolf algorithm
CN111882106A (en) * 2020-06-16 2020-11-03 江苏大学 Short-term power load prediction method based on comprehensive factors and CEEMD-IGWO-GRNN

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2580768A (en) * 1947-08-14 1952-01-01 Ibm Data look-up apparatus for computing or other machines
US4858462A (en) * 1989-01-20 1989-08-22 The Babcock & Wilcox Company Acoustic emission leak source location
JP2004219075A (en) * 2003-01-09 2004-08-05 Civil Engineering Research Institute Of Hokkaido Ae measurement body, function testing method for ae measurement body, ae measuring method, and active natural ground stabilization evaluating method
US8693286B1 (en) * 2009-10-16 2014-04-08 Snap-On Incorporated Position measurement for collision repair systems
CN108717201A (en) * 2018-06-20 2018-10-30 成都理工大学 A kind of tunnel surrounding microquake sources localization method
CN110345392A (en) * 2019-07-16 2019-10-18 辽宁石油化工大学 The localization method and device of oil pipeline leak source
CN110986747A (en) * 2019-12-20 2020-04-10 桂林电子科技大学 Landslide displacement combined prediction method and system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于GWO-ELM的逆变器开路故障诊断;姚芳;姜涛;刘明宇;董超群;郑帅;;电源学报(第01期) *

Also Published As

Publication number Publication date
CN112904277A (en) 2021-06-04

Similar Documents

Publication Publication Date Title
CN110609321B (en) Micro seismic source positioning method based on speed model database
CN110119994B (en) Metal ore blasting pile displacement extraction and prediction method supported by GIS
CN107340516B (en) Combined logic fast track starting method based on Doppler velocity
CN113553763A (en) Tunnel surrounding rock rapid grading method and system
CN112904277B (en) Tunnel surrounding rock breaking point positioning method based on improved gray wolf algorithm
CN108896962A (en) Iteration localization method based on sound position fingerprint
CN104374828A (en) Ultrasonic tomography imaging method of detection on hidden defect
CN104749630B (en) Method for constructing microseism monitoring velocity model
CN110006994B (en) Nondestructive testing method for multiple defects in built building structure
CN108574290B (en) Method, device, terminal and readable storage medium for positioning oscillation source of forced oscillation
CN108008353B (en) Method for ensuring anchor position stability by utilizing anchor point mutual ranging
CN112285650B (en) Method, system and storage medium for positioning unknown wave velocity sound emission source in presence of abnormal TDOA
Li et al. Automatic arrival-time picking of P-and S-waves of microseismic events based on object detection and CNN
CN112415571B (en) Microseism positioning method, storage medium and device
CN111767546B (en) Deep learning-based input structure inference method and device
CN111221036B (en) Target area seismic source positioning method and system containing unknown cavity
CN110441815B (en) Simulated annealing Rayleigh wave inversion method based on differential evolution and block coordinate descent
CN114063150A (en) Seismic source-station velocity model selection method based on ML-KNN algorithm
CN110243939B (en) Dam defect detection system based on mechanical wave positioning and working method thereof
CN111189926A (en) Method and system for identifying structure hole position based on global search
CN111965729A (en) Real-time monitoring method, system and device for vibroseis combination center
CN116593295B (en) Method and device for improving acoustic emission positioning precision by utilizing rock anisotropic wave velocity
CN116299374B (en) Sonar imaging underwater automatic calibration positioning method and system based on machine vision
CN116561563B (en) Slope displacement prediction method and related device based on residual prediction model
Huo et al. A study of impact point detecting method based on seismic signal

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