阅读说明：本技术 一种含圆孔结构的声发射/微震事件定位方法、系统、终端及可读储存介质 (Acoustic emission/microseismic event positioning method, system, terminal and readable storage medium containing round hole structure ) 是由 尚雪义 刘彩云 陈勇 郑雨晴 王进华 于 2021-07-30 设计创作，主要内容包括：本发明公开了一种含圆孔结构的声发射/微震事件定位方法、系统、终端及可读存储介质。该方法利用三角函数推导出含圆孔平面内不同震源和传感器空间位置关系下的P波传播最短距离计算方法；再借助立体几何和三角函数等知识,推导三维空间含全通圆孔和未通圆孔结构震源到传感器的P波传播最短距离的计算方法,从而将该最短路径算法引入定位目标函数的计算中,最后基于拾取的P波初至数据、定位目标函数反演出含圆孔结构的声发射/微震事件位置,本发明适用于二维平面结构和三维空间结构,可准确计算含圆孔结构的声发射/微震信号最短传播路径的精确解,使用了数据质量更高的P波初至数据及网格搜索算法,更易获得全局最优解,定位精度高、稳定性好。(The invention discloses an acoustic emission/microseismic event positioning method, system, terminal and readable storage medium with a round hole structure. The method utilizes trigonometric functions to derive a P wave propagation shortest distance calculation method under different seismic source and sensor spatial position relations in a plane containing circular holes; and deducing a calculation method of the shortest distance of P wave propagation from a seismic source containing a full-through circular hole structure and a non-through circular hole structure to a sensor in a three-dimensional space by means of knowledge such as solid geometry, trigonometric function and the like, thereby introducing the shortest path algorithm into the calculation of a positioning objective function, and finally performing acoustic emission/microseismic event positions containing the circular hole structure based on the picked P wave first arrival data and the positioning objective function.)

1. A method for locating acoustic emission/microseismic events with a circular hole structure is characterized by comprising the following steps:

picking up P wave first-break data of the acoustic emission/microseismic event, wherein a plurality of acoustic emission/microseismic sensors are arranged;

constructing a positioning objective function based on a P wave travel time equation, wherein a P wave propagation shortest distance formula between an acoustic emission source/seismic source and a sensor is set according to different spatial position relations of the acoustic emission source/seismic source, the sensor and a circular hole, and the structure containing the circular hole is a planar structure, a three-dimensional non-through circular hole structure or a three-dimensional all-through circular hole structure;

and deducing the acoustic emission/microseismic event position of the structure containing the round hole based on the picked P wave first-motion data and the positioning objective function inversion.

2. The method of claim 1, wherein: when the structure containing the round hole is a three-dimensional full-through round hole structure or a round hole in a plane structure, determining the shortest P wave propagation distance between an acoustic emission source/seismic source and a sensor according to the following rules:

a. the straight line segment determined by the acoustic emission source/seismic source to the sensor does not intersect with the round hole, and the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is the straight line distance between the acoustic emission source/seismic source and the sensor;

b. a straight line segment determined by the acoustic emission source/seismic source to the sensor does not penetrate through the round hole, but the extension line of the straight line segment intersects with the round hole, and the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is the straight line distance between the acoustic emission source/seismic source and the sensor;

c. the straight line segment determined by the acoustic emission source/seismic source to the sensor passes through the round hole, and the shortest distance of P wave propagation between the acoustic emission source/seismic source and the sensor is as follows:

L′_min(i)＝min(le_i+rad′_i+ls_i)

of formula (II) to'_min(i)Minimum distance of P wave propagation from acoustic source/seismic source to ith sensor, le_iFor the straight-line distance, ls, of P-wave signal propagation from the acoustic source/seismic source to the tangent point G_iIs the straight-line segment distance, rad ', that the P-wave signal travels from the exit tangent point F to the ith sensor'_iThe short arc length between the tangent point G and the tangent point F is set;

in the three-dimensional structure, the tangent point G is a tangent point M₁(xe_i1，ye_i1，ze_i1) M 'tangent point'₁(xe_i2，ye_i2，ze_i2) One tangent point of (1), tangent point F is tangent point N₁(xs_i1，ys_i1，zs_i1) And tangent point N'₁(xs_i2，ys_i2，zs_i2) One tangent point of (a); tangent point M₁、M′₁、N₁、N′₁Four points of tangency, xe, with the circular hole when P-wave signal propagates in three-dimensional structure_i1，ye_i1，ze_i1,xe_i2，ye_i2，ze_i2,xs_i1，ys_i1，zs_i1,xs_i2，ys_i2，zs_i2All are tangent point position coordinates;

in a two-dimensional planar structure, when the plane coordinates are represented by x and z axes, the tangent point G is a tangent point M (xe)_i1，ze_i1) Tangent point M' (xe)_i2，ze_i2) Is one tangent point of (1), tangent point F is tangent point N (xs)_i1，zs_i1) Tangent point N' (xs)_i2，zs_i2) One tangent point of (a); the tangent points M, M 'and N, N' are four tangent points xe between the P-wave signal and the round hole when the P-wave signal propagates in the planar structure_i1，ze_i1,xe_i2，ze_i2,xs_i1，zs_i1,xs_i2，zs_i2All are tangent point position coordinates;

when the structure containing the round hole is a structure without passing through the round hole, the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is determined according to the following rules:

A. if the straight line section between the acoustic emission source/seismic source and the sensor does not intersect with the round hole, the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is the straight line distance between the acoustic emission source/seismic source and the sensor;

B. if the straight line segment between the acoustic emission source/seismic source and the sensor is intersected with the round hole: (1) if the acoustic emission source/seismic source and the sensor are positioned at the side of the circular hole according to c L'_minCalculating the shortest P wave propagation distance between an acoustic emission source/seismic source and a sensor by a formula; (2) if the acoustic emission source/seismic source and the sensor are both positioned on the solid side, the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is the linear distance between the acoustic emission source/seismic source and the sensor; (3) if the acoustic source/seismic source and the sensor are arranged on the solid side and the round hole side, the acoustic source/seismic source and the sensor do not participate in positioning calculation.

3. The method of claim 2, wherein: if the structure containing the round hole is a full-through round hole structure in a three-dimensional space, a tangent point G (xe)_i，ye_i，ze_i) Tangent point F (xs)_i，ys_i，zs_i) Short arc and long rad between'_iThe formula of (1) is as follows:

wherein, with tangent point G, F projection to three-dimensional containing round hole structure sample front, note that the projection point that tangent point G, F corresponds is: g' (xe)_i，C，ze_i)、F'(xs_i，C，zs_i) C is a constant, xe_i，ye_i，ze_i，xs_i，ys_i，zs_iIs a tangent point position coordinate;

taking the horizontal line on the left side of the circular hole as the center, respectively spreading the upper and lower cylindrical hole surfaces upwards and downwards, wherein l is the arc length between the projection points G 'and F', and is r theta₁，θ₁Is the included angle between the projection points G ', F' and the circle center of the circular hole, r is the radius of the circular hole, d_G′F′Representing the distance of the straight line segment between the projection points G ', F', existing：

If the structure containing the circular holes is a circular hole structure in a two-dimensional plane, the tangent point G (xe)_i，ze_i) Tangent point F (xs)_i，zs_i) Short arc and long rad between'_iThe formula of (1) is as follows:

rad′_i＝rθ，

wherein theta is an included angle between the tangent point G, F and the circle center of the circular hole which is less than 180 degrees, and d_qIndicating the straight line segment distance between tangent points G, F,

4. the method of claim 2, wherein: in a three-dimensional space, an acoustic emission source/seismic source, tangent points G ', F' and a sensor are positioned on the same plane x + my + nz + e which is 0, and m, n and e are plane equation coefficients and satisfy the following conditions:

wherein the coordinates of the acoustic emission source/seismic source are A (x)₀，y₀，z₀) Sensor coordinate is B (x)_i，y_i，z_i) I is the sensor number;

shortest propagation distance L'_min(i)The solving formula of (2) is as follows:

diff(le_i+rad′_i+ls_i)＝0

wherein diff is a partial derivative function.

5. The method of claim 2, wherein: if the coordinates of the acoustic source/seismic source are A (x)₀，y₀，z₀) Sensor coordinate is B (x)_i，y_i，z_i) I is the number of the sensor, the acoustic emission source/the seismic source and the sensor are projected to the front surface of the three-dimensional sample containing the round hole structure, and the projection coordinate of the acoustic emission source/the seismic source in a projection plane is A' (x)₀，C，z₀) Projection coordinate of the sensor is B' (x)_i，C，z_i) And C is a constant, namely the y-axis coordinate of the front face, the judgment of different spatial position relations of the acoustic emission source/seismic source, the sensor and the round hole is as follows:

if the center of circle O (x) of the circular hole in the projection plane_c，z_c) The distance d 'from the straight line a' x + b 'z + c' to 0 is more than or equal to r, and r is the radius of the circular hole, so that the straight line segment determined from the acoustic emission source/seismic source to the sensor in the three-dimensional space does not intersect with the circular hole; wherein, the straight line a 'x + B' z + c 'is 0, and a', B 'and c' are linear equation coefficients of a 'and B' in the projection plane;

if the distance d ' < r and the angle OA ' B ' and angle OB ' A ' are not equal to acute angles, a straight line segment determined by an acoustic emission source/seismic source-sensor in the three-dimensional space does not pass through the round hole, but the extension line of the straight line segment is intersected with the round hole;

if the distance d ' < r and the & lt OA ' B ' and & lt OB ' A ' are acute angles, a straight line segment determined by an acoustic emission source/seismic source to a sensor in the three-dimensional space passes through the round hole.

6. The method of claim 1, wherein: the construction process of the positioning target function based on the P wave travel time equation comprises the following steps: selecting a certain sensor as a reference, calculating the difference between a P wave travel time equation corresponding to other sensors and a P wave travel time equation corresponding to the reference sensor based on a double difference method, accumulating the difference values, and determining a positioning target by selecting the minimum value after accumulating the difference values;

wherein, if the first sensor is used as a reference, the formula of the positioning objective function is as follows:

where Minimize is the minimization function, n is the number of sensors, and the P-wave velocity is v_pThe first arrival time of the P wave picked up by the ith sensor is t_iThe first arrival time of the P wave picked up by the 1 st sensor is t₁，ΔT_i1P-wave first arrival system error, L ', of ith sensor relative to 1 st sensor'_min(i)Minimum distance of P wave propagation, L ', between acoustic emission source/seismic source to ith sensor'_min(1)The shortest distance of P wave propagation from an acoustic emission source/seismic source to the 1 st sensor, and the P wave first arrival system error corresponding to the ith sensor is T_i；

If the structure containing the circular hole is a three-dimensional structure without the circular hole, the formula of the positioning objective function is as follows:

7. the method of claim 1, wherein: the process of deducing the position of the acoustic emission source/seismic source containing the circular hole structure based on the picked P wave first-motion data and the positioning objective function inversion is to determine the position of the acoustic emission source/micro-seismic event containing the circular hole structure by adopting a grid search method, and the method specifically comprises the following steps:

firstly, obtaining an initial position of an acoustic emission/microseismic event by utilizing sparse grid search;

and then, determining a target area by taking the preliminary position as a center, and meshing the target area to search the acoustic emission/microseismic event position.

8. A system based on the method of any one of claims 1-7, characterized by: the method comprises the following steps:

the P wave first arrival data pickup unit is used for picking up P wave first arrival data of the acoustic emission/microseismic event, wherein a plurality of acoustic emission/microseismic sensors are arranged;

a positioning objective function establishing unit: constructing a positioning target function based on a P wave travel time equation;

the shortest path calculation unit is used for calculating the shortest P wave propagation distance between an acoustic emission source/seismic source with a round hole structure and a sensor, wherein the shortest P wave propagation distance formula between the acoustic emission source/seismic source and the sensor is set according to different spatial position relations between the acoustic emission source/seismic source, the sensor and the round hole, and the round hole structure is a planar structure, a three-dimensional non-through round hole structure or a three-dimensional full-through round hole structure;

and the acoustic emission/microseismic event positioning unit is used for deriving the acoustic emission/microseismic event position containing the circular hole structure based on the picked P wave first-motion data and the positioning objective function inversion.

9. A terminal, characterized by: comprising one or more processors and memory, the memory storing a computer program that is invoked by the processors to implement:

a method as claimed in any one of claims 1 to 7, wherein the method comprises the step of locating acoustic emission/microseismic events having a circular aperture structure.

10. A readable storage medium, characterized by: a computer program is stored, which is invoked by a processor to implement:

a method as claimed in any one of claims 1 to 7, wherein the method comprises the step of locating acoustic emission/microseismic events having a circular aperture structure.

Technical Field

The invention belongs to the field of acoustic emission/microseismic monitoring, and particularly relates to an acoustic emission/microseismic event positioning method, system, terminal and readable storage medium with a round hole structure.

Background

Acoustic emission/microseismic (AE/MS) source localization technology is widely used as an effective structural health monitoring method in the fields of rock mechanics, metal processing, and the like. The elastic waves generated by the rapid release of the structure during the rupture process propagate in different orientations and are captured by spatially arranged sensors. By means of P wave first-arrival data of waves recorded by the sensor, the time-space distribution characteristics (time and position) of acoustic emission/microseismic events can be obtained so as to deduce the mechanical state of the structure, and further effective prevention and control measures (the microseismic and acoustic emission are different in scale, but the positioning principle is similar) are taken. The seismic source positioning described above generally assumes that the structure has a uniform velocity and the wave travels in a straight path. However, for the hole structure such as a pipeline, a tunnel, a mine tunnel, etc., the acoustic emission/microseismic signal does not propagate along a straight line, and the straight line positioning method may have a large positioning error. The method has important significance on monitoring the structural health by realizing the accurate positioning of the acoustic emission/microseismic event of the three-dimensional structure containing the circular holes.

A travel-time-based ray tracing method is commonly used for AE/MS source positioning, and the method utilizes the difference between the observed time and the theoretical time of a P wave to establish a target function. Geiger (1910) establishes a classical linear positioning method, converts a time difference equation into a linear equation set, and then iteratively solves the position of a seismic source by using a least square method. Furthermore, Waldhauser and Ellsworth (2000) proposes a double-difference positioning method, which assumes that propagation paths of two similar seismic event excitation wave fields are similar, and effectively reduces the influence of structural abnormality on travel time on a common path when similar earthquakes propagate to a station; li et al (2016) gives four common localization objective functions based on a combination of time-difference, double-difference, L1 and L2 norms. The convergence of the positioning objective function is often closely related to the selected optimization inversion algorithm, and researchers introduce various iteration techniques to solve the optimal solution of the objective function, such as Newton's iteration, which has strong dependence on an initial point, is easy to enter a local optimal solution, and has poor positioning stability. Therefore, some more global methods are used to solve the above positioning objective function, such as simplex method, grid search method, particle swarm optimization algorithm, bayesian method, etc.

However, in complex structural geometries and complex materials, AE/MS events cannot be accurately located using the above-described location methods, typically assuming constant wave velocity and straight-line path propagation. To this end, Baxter (2007) developed a Delta T mapping technique (DTM) based on first-arrival time to overcome this drawback, which uses threshold crossings to determine the first-arrival time of a wave at the sensor, which can locate the AE source in different structural complexity materials, but which may generate false locations when the actual signal of the first threshold crossing begins to fall below a set threshold level. Therefore, Safaa Kh. Al-Juglali (2016) proposes an improved full-automatic DTM positioning method, a clustering algorithm is used for automatically identifying and selecting highly-relevant events of each grid point, meanwhile, a 'minimum difference' method is used for determining the position of a seismic source, and the positioning precision is improved. In addition, the donglong et al (2017) propose a multi-step positioning method (MLM) for unpredicted speed of heterogeneous complex propagation media and a comprehensive analytic solution without a preliminary speed measurement based on a positioning function of a time difference of arrival (TD) model, and further improve positioning accuracy. For the AE/MS seismic source location of a complex structure with a circular hole, researchers carry out a series of researches on the location of a cylinder without considering the thickness, two points on the cylinder are projected onto a circular ring, then the cylinder is unfolded to a plane, and the P wave propagation distance and time are easy to obtain according to the Pythagorean theorem. However, the structure containing the circular holes usually has a certain thickness, and the universality of the method is poor. Therefore, the grid with the round hole structure is discretized by the Dong Long et al (2020) and the Jiang et al (2017), and then the shortest path is found by respectively using the A search algorithm and the FMM algorithm, so that a better positioning effect is obtained on the round hole structure. However, the grid discretization of the process causes deviation of the obtained shortest path, the positioning accuracy is reduced, and the algorithm is complex. Therefore, it is necessary to provide a method for calculating the precise shortest distance between two points in the two-dimensional/three-dimensional circular hole-containing structure, so as to improve the positioning precision.

Disclosure of Invention

The invention aims to provide an acoustic emission/microseismic event seismic source positioning technology applied to a complex structure containing a round hole, in particular to an acoustic emission/microseismic event accurate positioning method, a system, a terminal and a readable storage medium containing a round hole structure, which utilize a trigonometric function to deduce a method for accurately calculating the shortest distance of P wave propagation in a plane containing the round hole under the spatial position relationship of different seismic sources and sensors; then deducing an accurate calculation method for the shortest distance from the seismic source to the sensor in a three-dimensional space containing a full-through round hole and a non-through round hole structure in a plane containing the round holes by means of knowledge such as solid geometry, trigonometric functions and the like; therefore, the method can accurately calculate the precise solution of the shortest propagation path of the acoustic emission/microseismic signal P wave with the circular hole structure, and is suitable for a plane structure and a three-dimensional space.

On one hand, the invention provides a method for accurately positioning an acoustic emission/microseismic event with a circular hole structure, which comprises the following steps:

picking up acoustic emission/microseismic P wave first arrival data, wherein a plurality of acoustic emission/microseismic sensors are arranged;

constructing a positioning objective function based on a P wave travel time equation, wherein a P wave propagation shortest distance formula between an acoustic emission source/seismic source and a sensor is set according to different spatial position relations of the acoustic emission source/seismic source, the sensor and a circular hole, and the structure containing the circular hole is a planar structure, a three-dimensional non-through circular hole structure or a three-dimensional all-through circular hole structure;

and deducing the position of an acoustic emission source/seismic source containing a circular hole structure based on the picked P wave first-motion data and the positioning target function inversion.

Optionally, when the structure containing the circular hole is a three-dimensional full-through circular hole structure or a circular hole in a planar structure, the shortest P-wave propagation distance between the acoustic emission source/seismic source and the sensor is determined according to the following rules:

a. the straight line segment determined by the acoustic emission source/seismic source to the sensor does not intersect with the round hole, and the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is the straight line distance between the acoustic emission source/seismic source and the sensor;

b. a straight line segment determined by the acoustic emission source/seismic source to the sensor does not penetrate through the round hole, but the extension line of the straight line segment intersects with the round hole, and the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is the straight line distance between the acoustic emission source/seismic source and the sensor;

c. the straight line segment determined by the acoustic emission source/seismic source to the sensor passes through the round hole, and the shortest distance of P wave propagation between the acoustic emission source/seismic source and the sensor is as follows:

L′_min＝min(le_i+rad'_i+ls_i)

of formula (II) to'_minMinimum distance of P wave propagation from acoustic source/seismic source to ith sensor, le_iFor the straight-line distance, ls, of P-wave signal propagation from the acoustic source/seismic source to the tangent point G_iIs the straight-line segment distance, rad ', that the P-wave signal travels from the exit tangent point F to the ith sensor'_iThe short arc length between the tangent point G and the tangent point F is set; in the three-dimensional structure, the tangent point G is a tangent point M₁(xe_i1，ye_i1，ze_i1) M 'tangent point'₁(xe_i2，ye_i2，ze_i2) One tangent point of (1), tangent point F is tangent point N₁(xs_i1，ys_i1，zs_i1) And tangent point N'₁(xs_i2，ys_i2，zs_i2) One tangent point of (a); tangent point M₁、M′₁、N₁、N′₁Four points of tangency, xe, with the circular hole when P-wave signal propagates in three-dimensional structure_i1，ye_i1，ze_i1,xe_i2，ye_i2，ze_i2,xs_i1，ys_i1，zs_i1,xs_i2，ys_i2，zs_i2All are tangent point position coordinates;

in a two-dimensional planar structure, when the plane coordinates are represented by x and z axes, the tangent point G is a tangent point M (xe)_i1，ze_i1) Tangent point M' (xe)_i2，ze_i2) Is one tangent point of (1), tangent point F is tangent point N (xs)_i1，zs_i1) Tangent point N' (xs)_i2，zs_i2) One tangent point of (a); the tangent points M, M 'and N, N' are four tangent points xe between the P-wave signal and the round hole when the P-wave signal propagates in the planar structure_i1，ze_i1,xe_i2，ze_i2,xs_i1，zs_i1,xs_i2，zs_i2All are tangent point position coordinates;

when the structure containing the round hole is a structure without passing through the round hole, the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is determined according to the following rules:

A. if the straight line section between the acoustic emission source/seismic source and the sensor does not intersect with the round hole, the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is the straight line distance between the acoustic emission source/seismic source and the sensor;

B. if the straight line segment between the acoustic emission source/seismic source and the sensor is intersected with the round hole: (1) if the acoustic emission source/seismic source and the sensor are both positioned at the circular hole side (the structure of the circular hole which is not communicated can be divided into a solid structure and a structure containing the full through hole, the circular hole side is the structure containing the full through hole, and the acoustic emission source/seismic source and the sensor can be judged to belong to the circular hole side or the solid side according to the coordinate range of the acoustic emission source/seismic source and the sensor in the circular hole direction), and the acoustic emission source/seismic source and the sensor are positioned according to c L'_minCalculating the shortest P wave propagation distance between an acoustic emission source/seismic source and a sensor by a formula; (2) if the acoustic emission source/seismic source and the sensor are both positioned on the solid side, the shortest P wave transmission distance between the acoustic emission source/seismic source and the sensor is the linear distance between the acoustic emission source/seismic source and the sensor; (3) if the acoustic emission source/seismic source and the sensor are arranged on the solid side and the round hole side respectively, the propagation condition of acoustic emission/micro-seismic signals is complex and the acoustic emission/micro-seismic signals do not participate in positioning calculation.

Alternatively, if the structure containing the circular hole is a circular hole structure in a three-dimensional space, the tangent point G (xe)_i，ye_i，ze_i) Tangent point F (xs)_i，ys_i，zs_i) Short arc and long rad between'_iThe formula of (1) is as follows:

wherein, with tangent point G, F projection to three-dimensional containing round hole structure sample front, note that the projection point that tangent point G, F corresponds is: g' (xe)_i，C，ze_i)、F'(xs_i，C，zs_i) And C is a constant,xe_i，ye_i，ze_i，xs_i，ys_i，zs_iis a tangent point position coordinate;

taking the horizontal line on the left side of the circular hole as the center, respectively spreading the upper and lower cylindrical hole surfaces upwards and downwards, wherein l is the arc length between the projection points G 'and F', and is r theta₁，θ₁Is the included angle between the projection points G ', F' and the circle center of the circular hole, r is the radius of the circular hole, d_G′F′Representing the straight-line segment distance between projection points G ', F', there exists:

if the structure containing the circular holes is a circular hole structure in a two-dimensional plane, the tangent point G (xe)_i，ze_i) Tangent point F (xs)_i，zs_i) Short arc and long rad between'_iThe formula of (1) is as follows:

rad'_i＝rθ，

wherein theta is an included angle between the tangent point G, F and the circle center of the circular hole which is less than 180 degrees, and d_qIndicating the straight line segment distance between tangent points G, F,

optionally, in a three-dimensional space, the acoustic source/seismic source, the tangent points G ', F' and the sensor are located on the same plane x + my + nz + d ═ 0, and m, n and e are plane equation coefficients, which satisfy:

wherein the coordinates of the acoustic emission source/seismic source are A (x)₀，y₀，z₀) Sensor coordinate is B (x)_i，y_i，z_i) I is the sensor number;

shortest propagation distance L'_minThe solving formula of (2) is as follows:

diff(le_i+rad'_i+ls_i)＝0

where dif is the partial derivative function.

Alternatively, if the coordinates of the acoustic source/seismic source are A (x)₀，y₀，z₀) Sensor coordinate is B (x)_i，y_i，z_i) I is the number of the sensor, the acoustic emission source/the seismic source and the sensor are projected to the front surface of the three-dimensional sample containing the round hole structure, and the projection coordinate of the acoustic emission source/the seismic source in a projection plane is A' (x)₀，C，z₀) Projection coordinate of the sensor is B' (x)_i，C，z_i) And C is a constant, namely the y-axis coordinate of the front face, the judgment of different spatial position relations of the acoustic emission source/seismic source, the sensor and the round hole is as follows:

if the center of circle O (x) of the circular hole in the projection plane_c，z_c) The distance d 'from the straight line a' x + b 'z + c' to 0 is more than or equal to r, and r is the radius of the circular hole, so that the straight line segment determined from the acoustic emission source/seismic source to the sensor in the three-dimensional space does not intersect with the circular hole; wherein, the straight line a 'x + B' z + c 'is 0, and a', B 'and c' are linear equation coefficients of a 'and B' in the projection plane;

if the distance d ' < r and the angle OA ' B ' and angle OB ' A ' are not equal to acute angles, a straight line segment determined by an acoustic emission source/seismic source-sensor in the three-dimensional space does not pass through the round hole, but the extension line of the straight line segment is intersected with the round hole;

if the distance d ' < r and the & lt OA ' B ' and & lt OB ' A ' are acute angles, a straight line segment determined by an acoustic emission source/seismic source to a sensor in the three-dimensional space passes through the round hole.

Optionally, the P-wave travel time equation-based positioning objective function is constructed by the following process: selecting a certain sensor as a reference, calculating the difference between a P wave travel time equation corresponding to other sensors and a P wave travel time equation corresponding to the reference sensor based on a double difference method, accumulating the difference values, and determining a positioning target by selecting the minimum value after accumulating the difference values;

wherein, if the first sensor is used as a reference, the formula of the positioning objective function is as follows:

where Minimize is the minimization function, n is the number of sensors, and the P-wave velocity is v_pThe first arrival time of the P wave picked up by the ith sensor is t_iThe first arrival time of the P wave picked up by the 1 st sensor is t₁，ΔT_i1P-wave first arrival system error, L ', of ith sensor relative to 1 st sensor'_min(i)Minimum distance of P wave propagation, L ', between acoustic emission source/seismic source to ith sensor'_min(1)The shortest distance of P wave propagation from an acoustic emission source/seismic source to the 1 st sensor, and the P wave first arrival system error corresponding to the ith sensor is T_i；

If the structure containing the circular hole is a three-dimensional structure without the circular hole, the formula of the positioning objective function is as follows:

it should be noted that the positioning objective functions corresponding to the three-dimensional full through hole and the planar circular hole-containing structure are also the same.

Optionally, the process of deriving the acoustic emission source/seismic source position containing the circular hole structure based on the picked P-wave first-arrival data and the positioning objective function inversion is to determine the acoustic emission source/microseismic event position containing the circular hole structure by using a grid search method, specifically as follows:

firstly, obtaining an initial position of an acoustic emission/microseismic event by utilizing sparse grid search;

and then, determining a target area by taking the preliminary position as a center, and meshing the target area to search the acoustic emission/microseismic event position.

In a second aspect, the present invention provides an acoustic emission/microseismic event localization system comprising a circular aperture structure, comprising:

the P wave first arrival data pickup unit is used for picking up acoustic emission/microseismic P wave first arrival data, wherein a plurality of acoustic emission/microseismic sensors are arranged;

a positioning objective function establishing unit: constructing a positioning target function based on a P wave travel time equation;

the shortest path calculation unit is used for calculating the shortest P wave propagation distance between an acoustic emission source/seismic source with a round hole structure and a sensor, wherein the shortest P wave propagation distance formula between the acoustic emission source/seismic source and the sensor is set according to different spatial position relations between the acoustic emission source/seismic source, the sensor and the round hole, and the round hole structure is a planar structure, a three-dimensional non-through round hole structure or a three-dimensional full-through round hole structure;

and the acoustic emission/microseismic event positioning unit is used for deriving the acoustic emission/microseismic event position containing the circular hole structure based on the picked P wave first-motion data and the positioning objective function inversion.

In a third aspect, the present invention provides a terminal comprising one or more processors and a memory, the memory storing a computer program for invocation by the processors to implement:

a method for positioning acoustic emission/microseismic events with a circular hole structure.

In a fourth aspect, the present invention provides a readable storage medium storing a computer program for invocation by a processor to implement:

a method for positioning acoustic emission/microseismic events with a circular hole structure.

Advantageous effects

The invention provides an acoustic emission/microseismic event positioning method with a round hole structure, which is applied to acoustic emission/microseismic event positioning of a complex structure with a round hole. The P wave propagation shortest distance formula between the acoustic emission source/seismic source and the sensor is set according to different spatial position relations between the acoustic emission source/seismic source, the sensor and the round hole, the type of the round hole-containing structure is a planar structure, a three-dimensional blind round hole structure or a three-dimensional full-through round hole structure, and particularly, the shortest distance accurate solution of the hole-containing structure in the plane containing the round hole and the three-dimensional space is derived by utilizing knowledge such as trigonometric function, solid geometry and the like, so that the method is suitable for the planar structure and the three-dimensional structure, and can effectively avoid positioning errors caused by the fact that acoustic emission/microseismic signals propagate along straight lines in the round hole structure.

In a further preferred scheme of the invention, a double-difference method is provided for inverting the acoustic emission/microseismic signal P wave first-motion system error to obtain corrected P wave first-motion data, so that the basic data obtained by positioning the target function is more accurate and effective.

In a further preferred scheme, the method adopts a grid search method to determine the position of the acoustic emission/microseismic event, starts to acquire a preliminary position through sparse grids, reduces the search range, and then utilizes dense grid search to obtain a final positioning result, so that the process is easy to obtain a global optimal solution, and has high positioning precision and good stability.

Drawings

FIG. 1 is a flow chart of a method of the present invention;

FIG. 2 is a flowchart of the calculation of the shortest propagation distance of P-wave signals with a circular hole structure;

FIG. 3 is a schematic diagram of the shortest path between two points of a planar circular hole-containing structure;

FIG. 4 is a schematic diagram of the shortest path between two points of a three-dimensional round hole structure;

FIG. 5 is a graph of the positioning results of test events for structures containing circular holes in a plane.

Detailed Description

The present invention will be further described with reference to the accompanying drawings 1-5.

The method idea of the invention is as follows: aiming at the problems that positioning errors possibly caused by linear path hypothesis when a seismic source with a circular hole complex structure is positioned and the obtained shortest path has deviation, positioning precision is reduced, the algorithm is complex and the like due to grid discretization in an A-x search algorithm and an FMM algorithm, the method for calculating the accurate shortest distance between two points in the circular hole structure is provided. The method utilizes knowledge such as trigonometric function, solid geometry and the like to deduce a calculation method of the shortest path from a plane and a three-dimensional seismic source with a circular hole structure to a sensor. Therefore, the invention searches the position of the acoustic emission source/seismic source based on the optimized shortest path calculation method, and improves the positioning precision. And finally, establishing a positioning target function after correcting the P wave first-motion system difference, and preferably searching the position of the sound emission source/seismic source by using a grid search algorithm with better global property, thereby further improving the positioning precision and having good stability.

The following will first set forth the reasoning process of the method for calculating the shortest path from an acoustic source/a seismic source to a sensor, which is provided by the present invention, and the seismic source is taken as an example.

As shown in FIG. 2, first, a trigonometric function is used to derive the shortest P-wave propagation distance L in the spatial position relationship between the seismic source and the sensor in the plane containing the circular hole_minThe accurate calculation method comprises the following steps:

i, calculating a linear equation from an acoustic emission/micro-seismic source to a sensor in a plane: let the coordinates of the in-plane acoustic emission/microseismic source be A (x)₀，z₀) Sensor coordinate is B (x)_i，z_i) And i is the sensor number. As shown in (1), (2), and (3) of fig. 3, if the five-pointed star represents a certain seismic source position and the triangle represents a certain sensor position, the equation of the straight line from the seismic source to the sensor i is:

z＝k(x-x₀)+z₀ (1)

wherein k is (z)_i-z₀)/(x_i-x₀). Further, the expression (1) is developed to obtain ax + bz + c as 0. Wherein, a is k, b is-1, c is-kx₀+z₀。

II, judging whether the straight line segment from the seismic source to the sensor is intersected with the inner circular hole of the plane:

calculating the center O (x) of the circular hole_c，z_c) Distance to the straight line ax + bz + c ═ 0:

a. if d is larger than or equal to r, and r is the radius of the circular hole, the circular hole is separated from or tangent to the straight line, as shown in the schematic diagram (1) in FIG. 3, the shortest transmission of the P wave signalDistance of broadcast L_minThe distance between the seismic source and the sensor is the straight-line distance:

b. if d is less than r, judging whether both the [ OAB ] and the [ OBA ] are acute angles by using a cosine law, wherein the method comprises the following steps:

in the formula, the angle OAB belongs to [0, pi ]]，∠OBA∈[0,π]. AB represents the distance of a straight line segment from the seismic source to the sensor, OA represents the distance of a straight line segment from the seismic source to the center of the circular hole, OB represents the distance of a straight line segment from the center of the circular hole to the sensor, and the distance can be determined according to the seismic source coordinate A (x)₀，z₀) Sensor coordinate B (x)_i，z_i) And center coordinates O (x)_c，z_c) And (4) calculating.

(1) If the [ OAB ] and the [ OBA ] are not both acute angles, the straight-line segment from the seismic source to the sensor does not pass through the circular hole, as shown in the schematic diagram of (2) in fig. 3, at this time, the shortest propagation path L_minCalculating the distance of a linear section of the seismic source and the sensor by using the formula (3);

(2) if both the [ OAB ] and the [ OBA ] are acute angles, the straight line segment from the seismic source to the sensor passes through the circular hole, which is a schematic diagram shown in (3) in fig. 3, at the moment, the tangent point of a P-wave signal propagation path and the circular hole needs to be calculated, and the tangent point of the P-wave signal starting from the seismic source and the circular hole meets the following formula:

two potential tangent point coordinates M (xe) of the P-wave signal from the seismic source to the round hole can be obtained by the formula (6)_i1，ze_i1) And M' (xe)_i2，ze_i2) (ii) a In the same way, the sensor reversely transmits to obtain two potential tangent points N and N', and the coordinates of the two potential tangent points are respectively marked as N (xs)_i1，zs_i1) And N' (zs)_i2，zs_i2) There are 4 combinations of tangents M, M ', N, N' with sensor B and source A, namely: AM + MN + NB, AM + MN '+ N' B, AM '+ M' N '+ N' B, and AM '+ M' N + NB. The short arc length of the corresponding combination on the circle is recorded as And

the distance of a P wave signal propagation straight line segment between the seismic source and the tangent point isAnd the P-wave signal travels a linear distance from the potential tangent point to the sensor ofWherein j (j ═ 1,2,3,4) corresponds to the above four combinations. Arc length between tangent pointsThe calculation formula is as follows:

wherein the content of the first and second substances,theta is an included angle between the two selected tangent points and the circle center of the circular hole of less than 180 degrees, and d_qIndicating the straight-line segment distance between the two selected tangent points,tangent point F (xs)_i，zs_i) Corresponding to N (xs)_i1，zs_i1) And N' (zs)_i2，zs_i2) Is a tangent point of (b), tangent point G (xe)_i，ze_i) Corresponding to M (xe)_i1，ze_i1) And M' (xe)_i2，ze_i2) One tangent point of (a);

then, in this case, the shortest propagation distance L of the P-wave signal from the seismic source to the sensor in the plane_minComprises the following steps:

the invention derives the shortest distance L 'from a seismic source to a sensor of an acoustic emission/microseismic P wave signal containing a circular hole structure in a three-dimensional space by virtue of knowledge such as solid geometry, trigonometric function and the like'_minThe precise calculation method of (1). The three-dimensional round hole containing structure is divided into a full-through round hole and a non-through round hole, wherein the non-through round hole is formed by drilling a part of depth of the round hole but not a through hole. As shown in figure 2, the sound emission/microseismic P-wave signal with the all-through circular hole structure has the shortest distance L 'from a seismic source to a sensor'_minThe precise calculation method comprises the following specific steps:

in this embodiment, the three-dimensional complex structure with all-through circular holes has a seismic source coordinate of A (x)₀，y₀，z₀) Sensor coordinate is B (x)_i，y_i，z_i) I is the sensor number, as shown in FIG. 4, the triangle represents a certain sensor position, the five-pointed star represents a certain seismic source position, and O is the center of the circle. Projecting the position of the seismic source and the position of the sensor to the front surface of the three-dimensional sample with the round hole structure (in the embodiment, the position is defined as a surface with the round hole), and judging whether a straight line section determined from the seismic source to the sensor intersects with the round hole according to the spatial position relation of the seismic source and the sensor in a projection plane (the condition in the projection plane can be shown in fig. 3).

In the projection plane: the projection coordinate of the seismic source is A' (x)₀，C，z₀) Projection coordinate of the sensor is B' (x)_i，C，z_i) C is constant (not involved in calculation), the source projection point of the projection plane is to the sensor projection pointThe equation of the straight line is:

z＝k′(x-x₀)+z₀ (9)

wherein k ═ z_i-z₀)/(x_i-x₀). Further, the formula (9) is expanded to obtain a 'x + b' z + c ═ 0. Wherein, a ' ═ k ', b ' ═ -1, c ' ═ -k ' x₀+z₀。

Calculating the center O (x) of the circle hole in the projection plane_c，z_c) Distance to the straight line a 'x + b' z + c ═ 0:

a. if in the projection plane: d' ≧ r, r is the radius of the circular hole, the straight line segment determined by the sensor from the seismic source of the P-wave signal in the three-dimensional space does not intersect with the circular hole, as shown in the schematic diagram (1) in FIG. 4. At this time, the P-wave propagates along a straight line in the three-dimensional space, and the shortest distance calculation method is as follows:

b. if in the projection plane: d ' < r, judging whether the & lt OA ' B ' and & lt OB ' A ' are both acute angles by using the cosine theorem, and referring to the formula (4) and the formula (5) by calculating the & lt OA ' B ' and & lt OB ' A '.

(1) If the angle OA 'B' and the angle OB 'A' are not both acute angles, the straight line segment from the three-dimensional space structure seismic source to the sensor does not pass through the round hole, as shown in the schematic diagram of (2) in FIG. 4, and the shortest propagation path L 'is formed at the moment'_minCalculating by using the formula (11);

(2) if < OA 'B' and < OB 'A' are both acute angles, the straight line segment from the seismic source to the sensor of the three-dimensional space structure passes through the circular hole, as shown in the schematic diagram of (3) in FIG. 4. The P-wave signal has two tangent points with the circular hole when propagating in the three-dimensional structure, and is set as G (xe) in a partial enlarged view at the cylindrical hole_i，ye_i，ze_i)、F(xs_i，ys_i，zs_i)。

Straight line distance le between seismic source and tangent point G_iAnd the linear distance ls that the P-wave signal travels from the tangent point F to the sensor_iThe calculation formula is as follows:

arc length rad between tangent points'_iAnd (3) calculating: the tangent point G, F was projected onto the front surface of the three-dimensional sample with a round hole structure, and the projected point was designated as G' (xe)_i，C，ze_i)、F'(xs_i，C，zs_i) Expanding the cylindrical hole, and projecting the arc length l between the points G 'and F', and the difference value | ye between the longitudinal coordinates of the tangent points_i-ys_iL the determined straight line segment GG ' and the arc length GF between the tangent points in the three-dimensional space form a right triangle, and then the rad ' is calculated according to the pythagorean theorem '_iThe following were used:

wherein l is the arc length between the projection points G 'and F' and the arc length of the derivation process of the shortest path of the structure containing the circular hole in the planeThe same applies to the calculation method, see formula (7), where,theta is the included angle between the projection points G 'and F' and the circle center of the circular hole, d_G′F′Representing the straight-line segment distance between the projection points G ', F',

in a three-dimensional space, if the path from a seismic source to a sensor position is shortest, the seismic source, the tangent points G ', F' and the sensor position should be on the same plane x + my + nz + d ═ 0, and m, n and e are plane equation coefficients, then:

according to the formula (14), the unknowns m, n, e, ye can be represented_i、ys_iExpressed as a function of some unknown number. Further, the shortest propagation distance L 'of the P-wave signal from the seismic source to the sensor in the three-dimensional space'_minCan be expressed as:

L′_min＝min(le_i+rad'_i+ls_i) (15)

calculating shortest propagation distance L'_minThis is equivalent to solving the solution corresponding to the above-mentioned unknown partial derivative being 0:

diff(le_i+rad'_i+ls_i)＝0 (16)

where dif is the partial derivative function.

It should be understood that for planar structures, equation (8) is equivalent to equation (15),and le_i、ls_iAre the same meaning and all refer to the straight line distance, ls, that the P-wave signal travels from the acoustic source/seismic source to the tangent point G_iThe straight-line distance that the P-wave signal travels from the exit tangent point F to the i-th sensor,and rad'_iBoth represent the arc length between tangent point G and tangent point F, the choice of which determines j. Therefore, the shortest propagation distance L 'is actually obtained from the equation (15) for the planar structure'_min。

Deducing the shortest path L' of acoustic emission/microseismic P wave signal containing blind circular hole_minThe specific analysis steps of the accurate calculation method are similar to the shortest propagation path of the P wave signal containing the round hole in the three-dimensional space.

a. If the shortest path between the seismic source and the sensor does not intersect with the cylindrical hole, the shortest path can be calculated by using a formula (11);

b. if the shortest path between the seismic source and the sensor intersects with the round hole: (1) if the seismic source position and the sensor position are positioned at the circular hole side, the calculation method of the shortest path containing the all-pass circular hole structure is the same as that of the shortest path containing the all-pass circular hole structure, and the following formulas (12) to (16) are included; (2) if the seismic source position and the sensor position are positioned on the solid side, the straight-line distance between the two points is the distance from the seismic source to the sensor of the P wave; (3) if the seismic source and the sensor are respectively positioned at the solid side and the round hole side, the acoustic emission/microseismic signal transmission condition is complex and does not participate in positioning calculation. In the subsequent grid searching process, if the seismic source and the sensor are respectively positioned at the solid side and the round hole side, the group of data does not participate in calculation.

Through the reasoning, a method for determining the shortest path is constructed, so that the technical problem provided by the invention can be solved by utilizing the shortest path and applying the existing method in the field for positioning. The following embodiments are examples of a positioning objective function and a grid search after removing P-wave first arrival system errors, and other possible embodiments are not limited to the implementation processes of the following embodiments.

The present invention will be further described with reference to the following examples.

Example 1:

as shown in fig. 1, the method for locating an acoustic emission/microseismic event with a circular hole structure provided in this embodiment includes the following steps:

step 1: and (3) an automatic and manual method is used for picking up P wave first arrival data of the acoustic emission/microseismic event.

In the embodiment, an AIC method is selected to pick up P-wave first-arrival data of acoustic emission/microseismic events, a Matlab is used for compiling a man-machine interaction program to display an automatic method picking-up effect, and finally P-wave first-arrival picking-up data with large picking-up errors are corrected manually.

Step 2: and constructing a positioning objective function for removing the P wave first arrival system error.

The invention provides a double-difference method inversion acoustic emission/microseismic signal P wave first-break system error based on an acoustic emission lead source/blasting event with a known position, and corrected P wave first-break data is obtained.

Let the acoustic emission/microseismic source origin time be t₀P wave velocity is v_pThe first arrival time of the picked P wave is t_iP wave first arrival system error is T_iThen the travel time equation is:

L_min＝v_p(t_i-t₀-T_i) (17)

subtracting the travel time equation corresponding to the ith sensor from the travel time equation corresponding to the 1 st sensor to obtain the P wave first arrival system error delta T of the ith sensor relative to the 1 st sensor_i1I.e. Delta T_i1＝T_i-T₁＝t_i-t₁-(L′_min(i)-L′_min(1))/v_p。

Based on the double difference method, the P wave is first arrived at the system error delta T before the acoustic emission positioning_i1Added to equation (17), the following localization objective function is obtained:

for three-dimensional non-through-hole holes, the mean value of the above objective functions is used as the objective function, i.e.n is the number of sensors involved in the positioning.

Wherein the shortest distance L of P wave propagation between the acoustic emission source/seismic source and the sensor_minCalculated according to the contents of the foregoing method.

It should be noted that, in this embodiment, the above formula is used as a positioning objective function, and in other feasible embodiments, other existing positioning objective functions may be selected, for example, the existing function with the smallest cumulative difference between the actual travel time and the theoretical travel time is used, and the shortest P-wave propagation distance L between the acoustic source/seismic source and the sensor is obtained by applying the improved method of the present invention_minThe selected positioning objective function is not limited.

And step 3: determining the position of an acoustic emission/microseismic event with a round hole structure by adopting a grid search method;

the mesh division is carried out on the structure containing the round holes, the denser the mesh is, the better the positioning effect is generally, but the calculation efficiency is reduced, and therefore, the multi-scale mesh search is adopted to determine the position of the acoustic emission/microseismic event: firstly, searching to obtain the approximate position of an acoustic emission/microseismic event by utilizing a sparse grid; and then determining a small area by taking the position as a center, dividing the small area into more dense grids again, and continuously searching to obtain the acoustic emission/microseismic event position with higher precision.

It should be noted that, when the grid search method is selected for searching in this embodiment, a grid central point may be selected as an assumed seismic source, so that a location objective function is used to search for a position of a acoustic emission/microseismic event, and a specific search process may refer to the prior art, which is not specifically limited in this respect; in other possible embodiments, the invention is not limited to this, i.e. the positioning objective function and the shortest distance L of P-wave propagation are set_minOn the basis, other existing searching methods can be selected for searching, and the method belongs to the technical means under the conception of the invention. Secondly, in this embodiment, step 1 is performed before step 2, and in other feasible embodiments, theoretical analysis may be performed first, that is, a positioning objective function is constructed first, and then data is acquired and actually acquired data is used for searching.

Example 2:

on the basis of the acoustic emission/microseismic event positioning method with the round hole structure, the invention provides an acoustic emission/microseismic event positioning system with the round hole structure, which comprises the following steps: the device comprises a shortest path calculation unit, a P wave first arrival data pickup unit, a positioning objective function establishing unit and an acoustic emission/microseismic event positioning unit.

Wherein the shortest path calculating unit: the method is used for calculating the shortest P wave propagation distance between an acoustic emission source/seismic source with a round hole structure and a sensor, wherein the shortest P wave propagation distance formula is set according to different spatial position relations between the acoustic emission source/seismic source, the sensor and the round hole, and the round hole structure is a planar structure, a three-dimensional non-through round hole structure or a three-dimensional all-through round hole structure. For the specific theoretical content and the calculation process, please refer to the content of the foregoing method.

P wave first arrival data pickup unit: p-wave first arrival data for picking up acoustic emission/microseismic events. In the embodiment, after P wave first arrival data of an acoustic emission/microseismic event is picked up by adopting an AIC method, a man-machine interaction program is compiled by using Matlab to display the picking effect of an automatic method, and finally, a unit for picking signals with larger errors and without P wave first arrival picking is corrected by using a manual mode;

a positioning objective function establishing unit: the positioning target function is used for establishing a positioning target function after the P wave first arrival system error is removed;

acoustic emission/microseismic event localization unit: and the acoustic emission/microseismic event position containing the round hole structure is deduced based on the picked P wave first-motion data and the positioning objective function inversion. In this embodiment, a grid search method is selected to determine the location of the acoustic emission/microseismic event with the circular hole structure.

It should be noted that, if the positioning objective function is the function described in embodiment 1, it preferably further includes a system error correction unit: and (3) inverting the error of the acoustic emission signal P wave first arrival system by using a double difference method to obtain corrected P wave first arrival data.

For the specific implementation process of each unit module, refer to the corresponding process of the foregoing method. It should be understood that, the specific implementation process of the above unit module refers to the method content, and the present invention is not described herein in detail, and the division of the above functional module unit is only a division of a logic function, and there may be another division manner in the actual implementation, for example, multiple units or components may be combined or may be integrated into another system, or some features may be omitted, or may not be executed. Meanwhile, the integrated unit can be realized in a hardware form, and can also be realized in a software functional unit form.

Example 3:

the present embodiment provides a terminal, which includes:

one or more processors;

a memory for storing one or more programs;

the one or more programs, when executed by the one or more processors, cause the one or more processors to implement:

step 1: picking up P wave first arrival data of the acoustic emission/microseismic event;

step 2: constructing a positioning objective function for removing P wave first arrival system errors;

and step 3: and deducing the acoustic emission/microseismic event position of the structure containing the round hole based on the picked P wave first-motion data and the positioning objective function inversion. In this embodiment, a grid search method is selected to determine the location of the acoustic emission/microseismic event with the circular hole structure.

The terminal further includes: and the communication interface is used for communicating with external equipment and carrying out data interactive transmission.

The memory may include high speed RAM memory, and may also include a non-volatile defibrillator, such as at least one disk memory.

If the memory, the processor and the communication interface are implemented independently, the memory, the processor and the communication interface may be connected to each other through a bus and perform communication with each other. The bus may be an industry standard architecture bus, a peripheral device interconnect bus, an extended industry standard architecture bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc.

Optionally, in a specific implementation, if the memory, the processor, and the communication interface are integrated on a chip, the memory, the processor, that is, the communication interface may complete communication with each other through the internal interface.

The specific implementation process of each step refers to the explanation of the foregoing method.

It should be understood that in the embodiments of the present invention, the Processor may be a Central Processing Unit (CPU), and the Processor may also be other general purpose processors, Digital Signal Processors (DSPs), Application Specific Integrated Circuits (ASICs), Field Programmable Gate Arrays (FPGAs) or other Programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, and the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The memory may include both read-only memory and random access memory, and provides instructions and data to the processor. The portion of memory may also include non-volatile random access memory. For example, the memory may also store device type information.

Example 4:

the present embodiments provide a readable storage medium storing a computer program for invocation by a processor to implement:

step 1: picking up P wave first arrival data of the acoustic emission/microseismic event;

step 2: constructing a positioning objective function for removing P wave first arrival system errors;

and step 3: and deducing the acoustic emission/microseismic event position of the structure containing the round hole based on the picked P wave first-motion data and the positioning objective function inversion. In this embodiment, a grid search method is selected to determine the location of the acoustic emission/microseismic event with the circular hole structure.

The specific implementation process of each step refers to the explanation of the foregoing method.

The readable storage medium is a computer readable storage medium, which may be an internal storage unit of the controller according to any of the foregoing embodiments, for example, a hard disk or a memory of the controller. The readable storage medium may also be an external storage device of the controller, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), and the like provided on the controller. Further, the readable storage medium may also include both an internal storage unit of the controller and an external storage device. The readable storage medium is used for storing the computer program and other programs and data required by the controller. The readable storage medium may also be used to temporarily store data that has been output or is to be output.

Based on such understanding, the technical solution of the present invention essentially or partially contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned readable storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.

Data validation

The effectiveness of the method is illustrated by taking an acoustic emission/microseismic theory event of a plane containing a circular hole structure as an example. Fig. 5 shows a cloud image of the positioning result of the test event of the planar structure with circular holes, wherein (a) and (b) in fig. 5 respectively correspond to the cloud images of the search time fitting error of the coarse grid and the fine grid. The length of a structure containing circular holes in a plane is set to be 30cm, the width of the structure is set to be 20cm, the radius r of the circular hole in the center of the structure is set to be 2.5cm, the position of an acoustic emission event is set to be (25, 5) cm, the coordinates of a sensor are shown in a table 1, and the P wave propagation speed is 0.6 cm/us. Firstly, a coarse grid is used for dividing a circular hole structure, time fitting errors of grid points obtained through searching are shown in fig. 5 (a), a pentagram is a test event position, and a rhombus is a positioning result. It can be seen that the closer to the source location, the smaller the temporal fitting error. Further, the time fitting error of the grid of 1cm around the rhombus position is searched as shown in fig. 5 (b), and the searched acoustic emission/microseismic event position is (25, 5) cm and is coincident with the test event position. In conclusion, the accurate positioning method for the acoustic emission/microseismic event with the circular hole structure provided by the invention has high positioning accuracy, and can obtain a stable positioning result by adopting a grid search method with better overall performance.

TABLE 1 sensor location coordinates

It should be emphasized that the examples described herein are illustrative and not restrictive, and thus the invention is not to be limited to the examples described herein, but rather to other embodiments that may be devised by those skilled in the art based on the teachings herein, and that various modifications, alterations, and substitutions are possible without departing from the spirit and scope of the present invention.