CN104360384A - Microseism event positioning method and device based on automatic scanning of longitudinal and transverse wave energy - Google Patents

Abstract

The invention relates to a microseism event positioning method and device based on automatic scanning of longitudinal and transverse wave energy. The method includes the steps that microseism event recognition is carried out on underground monitoring data; a microseism event travel-time table of a searched area is built; sliding scanning is carried out on a recognized microseism event according to the microseism event travel-time table, and accordingly diffraction stacked energy is obtained; sliding scanning is carried out on the polarization direction of the recognized microseism event according to the microseism event travel-time table, so information of the polarization direction is obtained; an objective function of microseism event direction consistence and scanned and stacked energy maximization is solved according to the diffraction stacked energy and the information of the polarization direction, and accordingly an optimal solution of the objective function is obtained. The optimal solution represents the position where direction stability and consistence are best and longitudinal and transverse wave scanned energy is strongest in the direction and the position is where the microseism event happens.

Description

Based on micro-seismic event localization method and the device of wave energy autoscan in length and breadth

Technical field

The present invention relates to pressure break micro-seismic monitoring data technique field, particularly a kind of micro-seismic event localization method based on wave energy autoscan in length and breadth and device.

Background technology

Pressure break micro-seismic monitoring signal is subject to the impact of rock burst energy size and background noise, and the Signal-to-Noise received differs greatly, and pickup micro-seismic event time of arrival and data process in real time affect greatly.

Pressure break micro-seismic monitoring signal duration is long, data volume is huge, the micro-seismic event that signal to noise ratio (S/N ratio) is high accurately can be picked up its time of arrival by automatic Picking algorithm, but the micro-seismic event of signal to noise ratio (S/N ratio) medium and on the lower side can not reach good precision by automatic Picking, there is larger impact to the direction calculating of micro-seismic event and positioning result.On-the-spot microseismic signal transacting needs the microearthquake not up to standard to first break pickup to carry out man-machine interactively correction to reach industrial treatment standard, largely be delayed the real-time of microearthquake signal processing results, adjust the ageing of pressing crack construction scheme in real time by micro-seismic monitoring result like this and receive impact.

Summary of the invention

For solving the problem of prior art, the present invention proposes a kind of micro-seismic event localization method based on wave energy autoscan in length and breadth and device, by ripple automatic energy scan method in length and breadth, micro-seismic event direction of earthquake source and position are calculated automatically, improves the ageing of micro-seismic monitoring data process.

For achieving the above object, the invention provides a kind of micro-seismic event localization method based on wave energy autoscan in length and breadth, the method comprises:

Micro-seismic event identification is carried out to underground monitoring data;

Set up the micro-seismic event whilst on tour table of region of search;

According to micro-seismic event whilst on tour table, slip scan is carried out to the micro-seismic event identified, obtain diffraction stack energy;

According to micro-seismic event whilst on tour table, slip scan is carried out to the polaried orientation of the micro-seismic event identified, obtain polaried orientation information;

According to described diffraction stack energy and described polaried orientation information, maximize objective function to micro-seismic event orientation consistance and scanning stack power to solve, obtain the optimum solution of objective function, this optimum solution is expressed as: the most stable and consistent in orientation and the position that ripple scanning energy is the strongest in length and breadth in this orientation, and this position is the position of micro-seismic event.

Preferably, the expression formula of described micro-seismic event orientation consistance and scanning stack power maximization objective function is:

f(Azimuth,E)＝f1(Azimuth)→min&&f2(E)→max；

Wherein, $f1\left(\mathrm{Azimuth}\right)=\mathrm{abs}(\mathrm{Azimuth}(x_i,y_j,z_k,t_p)-\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}\mathrm{Azimuth}(x_i,y_j,z_k,t_p)/N);$ F2 (E)=E (x _i, y _j, z _k, t); X above _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _i∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; Azimuth is orientation function, and E is energy function, and t is compressional wave or shear wave time, t _pit is the compressional wave time; Min represents that f1 (Azimuth) function gets minimum value, max represents that f2 (E) function gets maximal value, and min & & f2 (E) represents that need meet function f 1 (Azimuth) gets minimum value sum functions f2 (E) and get maximal value simultaneously; → be mathematic sign, represent infinite approach.

Preferably, describedly according to micro-seismic event whilst on tour table, the step that the micro-seismic event identified carries out slip scan to be comprised:

In region of search, three-dimensional each coordinate axis of sliding carries out diffraction energy superposition; Wherein, the computing formula of diffraction energy superposition is:

$E(x_i,y_j,z_k)=\underset{t=t(x_i,y_j,z_k)-\mathrm{\τ}}{\overset{t=t(x_i,y_j,z_k)+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t(x_i,y_j,x_k)\right)$

Wherein, x _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _j∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; A is amplitude, and t is the function of time, and τ is time window, and E is energy function;

According to micro-seismic event whilst on tour table sliding stack wave energy in length and breadth on signal time axle; Wherein, the expression formula superposing compressional wave energy is: $E\left(p\right)=\underset{{t=t}_p-\mathrm{\τ}}{\overset{{t=t}_p+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t_p\right);$ The expression formula of superposition shear wave energy is: $E\left(s\right)=\underset{{t=t}_s-\mathrm{\τ}}{\overset{{t=t}_s+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t_s\right);$ N is the number of wave detector, and E is energy function, and A is amplitude, t _pthe compressional wave time, t _sbe the shear wave time, t is the function of time.

Preferably, describedly according to micro-seismic event whilst on tour table, the step that the polaried orientation of the micro-seismic event identified carries out slip scan to be comprised:

According to micro-seismic event whilst on tour table, signal time axle slides and calculates the polarization angle of each time point micro-seismic event accordingly.

Preferably, described micro-seismic event whilst on tour table is set up based on correction rate model.

For achieving the above object, present invention also offers a kind of micro-seismic event locating device based on wave energy autoscan in length and breadth, this device comprises:

Recognition unit, for carrying out micro-seismic event identification to underground monitoring data;

Unit set up by whilst on tour table, for setting up the micro-seismic event whilst on tour table of region of search;

First slip scan unit, for carrying out slip scan according to micro-seismic event whilst on tour table to the micro-seismic event identified, obtains diffraction stack energy;

Second slip scan unit, for carrying out slip scan according to micro-seismic event whilst on tour table to the polaried orientation of the micro-seismic event identified, obtains polaried orientation information;

Positioning unit, for according to described diffraction stack energy and described polaried orientation information, maximize objective function to micro-seismic event orientation consistance and scanning stack power to solve, obtain the optimum solution of objective function, this optimum solution is expressed as: the most stable and consistent in orientation and the position that ripple scanning energy is the strongest in length and breadth in this orientation, and this position is the position of micro-seismic event.

Preferably, the micro-seismic event orientation consistance that described positioning unit uses with the expression formula of scanning stack power maximization objective function is:

f(Azimuth,E)＝f1(Azimuth)→min&&f2(E)→max；

Wherein, $f1\left(\mathrm{Azimuth}\right)=\mathrm{abs}(\mathrm{Azimuth}(x_i,y_j,z_k,t_p)-\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}\mathrm{Azimuth}(x_i,y_j,z_k,t_p)/N);$ F2 (E)=E (x _i, y _j, z _k, t); X above _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _i∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; Azimuth is orientation function, and E is energy function, and t is compressional wave or shear wave time, t _pit is the compressional wave time; Min represents that f1 (Azimuth) function gets minimum value, max represents that f2 (E) function gets maximal value, and min & & f2 (E) represents that need meet function f 1 (Azimuth) gets minimum value sum functions f2 (E) and get maximal value simultaneously; → be mathematic sign, represent infinite approach.

Preferably, described first slip scan unit comprises:

Three-dimensional coordinate slip scan module, in region of search, three-dimensional each coordinate axis of sliding carries out diffraction energy superposition; Wherein, the computing formula of diffraction energy superposition is:

$E(x_i,y_j,z_k)=\underset{t=t(x_i,y_j,z_k)-\mathrm{\τ}}{\overset{t=t(x_i,y_j,z_k)+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t(x_i,y_j,x_k)\right)$

Wherein, x _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _j∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; A is amplitude, and t is the function of time, and τ is time window, and E is energy function;

Signal time axle slip scan module, for according to micro-seismic event whilst on tour table on signal time axle sliding stack in length and breadth wave energy wherein, the expression formula of superposition compressional wave energy is: $E\left(p\right)=\underset{{t=t}_p-\mathrm{\τ}}{\overset{{t=t}_p+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t_p\right);$ The expression formula of superposition shear wave energy is: $E\left(s\right)=\underset{{t=t}_s-\mathrm{\τ}}{\overset{{t=t}_s+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t_s\right);$ N is the number of wave detector, and E is energy function, and A is amplitude, t _pthe compressional wave time, t _sbe the shear wave time, t is the function of time.

Preferably, described second slip scan unit specifically for:

According to micro-seismic event whilst on tour table, signal time axle slides and calculates the polarization angle of each time point micro-seismic event accordingly.

Preferably, described whilst on tour table is set up unit and is set up micro-seismic event whilst on tour table based on correction rate model.

Technique scheme has following beneficial effect: the technical program is in pressure break micro-seismic monitoring data process location, utilizes the mode of ripple scanning in length and breadth to carry out inverting location to the orientation of micro-seismic event and position.Examples prove the technical program can carry out generation orientation and the process of diffraction energy autoscan to the microearthquake signal identified, finally by the locus of orientation consistance and energy maximization objective function determination micro-seismic event, process does not need artificial adjustment, significantly improves the ageing of process.

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

Fig. 1 is a kind of micro-seismic event localization method process flow diagram based on wave energy autoscan in length and breadth that the present invention proposes;

Fig. 2 is a kind of micro-seismic event locating device block diagram based on wave energy autoscan in length and breadth that the present invention proposes;

Fig. 3 is microearthquake signal ripple signal schematic representation in length and breadth in the present embodiment;

Fig. 4 is the micro-seismic event RMS amplitude figure of the present embodiment;

Fig. 5 is the micro-seismic event amplitude ratio schematic diagram of the present embodiment;

Fig. 6 is the present embodiment wave diffraction energy scan result schematic diagram in length and breadth;

Fig. 7 is the present embodiment polarization scans result schematic diagram;

Fig. 8 is the micro-seismic event locus schematic diagram that the present embodiment is determined.

Embodiment

Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.

The principle of work of technical scheme of the present invention is: in pressure break micro-seismic monitoring data localization process, utilizes the mode of ripple scanning in length and breadth to carry out inverting location to the orientation of micro-seismic event and position.First be that event recognition is carried out to raw readings, set up the micro-seismic event whilst on tour table of region of search, carry out polaried orientation according to compressional wave whilst on tour to scan, carry out diffraction stack energy scan by ripple whilst on tour in length and breadth again, finally jointly determine the locus of micro-seismic event according to stable scan position consistance and energy maximization objective function.

Based on above-mentioned principle of work, the present invention proposes a kind of micro-seismic event localization method based on wave energy autoscan in length and breadth, as shown in Figure 1.The method comprises:

Step 101): micro-seismic event identification is carried out to underground monitoring data;

Step 102): the micro-seismic event whilst on tour table setting up region of search;

Step 103): according to micro-seismic event whilst on tour table, slip scan is carried out to the micro-seismic event identified, obtain diffraction stack energy;

Step 103) described in slip scan divide two stages: be first region of search slide three-dimensional each coordinate axis carry out diffraction energy superposition; Moreover, according to ripple whilst on tour table in length and breadth at the upper sliding stack wave energy in length and breadth of signal time axle (t direction).

Step 104): according to micro-seismic event whilst on tour table, slip scan is carried out to the polaried orientation of the micro-seismic event identified, obtain polaried orientation information;

Step 104) slip scan be calculate the polarization angle of each time point upper slip of signal time axle (t direction) based on micro-seismic event whilst on tour table;

Step 105): according to described diffraction stack energy and described polaried orientation information, maximize objective function to micro-seismic event orientation consistance and scanning stack power to solve, obtain the optimum solution of objective function, this optimum solution is expressed as: the most stable and consistent in orientation and the position that ripple scanning energy is the strongest in length and breadth in this orientation, this position is the position of micro-seismic event, completes the location of micro-seismic event.

Step 105) described in optimum solution refer to that namely position that in the orientation determined in geophone orientation consistance, diffraction stack energy is maximum is the locus (determining the X-coordinate of micro-seismic event, Y-coordinate and Z coordinate) of this micro-seismic event.

As shown in Figure 2, be a kind of micro-seismic event locating device block diagram based on wave energy autoscan in length and breadth that the present invention proposes.This device comprises:

Recognition unit 201, for carrying out micro-seismic event identification to underground monitoring data;

Unit 202 set up by whilst on tour table, for setting up the micro-seismic event whilst on tour table of region of search;

First slip scan unit 203, for carrying out slip scan according to micro-seismic event whilst on tour table to the micro-seismic event identified, obtains diffraction stack energy;

Second slip scan unit 204, for carrying out slip scan according to micro-seismic event whilst on tour table to the polaried orientation of the micro-seismic event identified, obtains polaried orientation information;

Positioning unit 205, for according to described diffraction stack energy and described polaried orientation information, maximize objective function to micro-seismic event orientation consistance and scanning stack power to solve, obtain the optimum solution of objective function, this optimum solution is expressed as: the most stable and consistent in orientation and the position that ripple scanning energy is the strongest in length and breadth in this orientation, and this position is the position of micro-seismic event.

Preferably, the micro-seismic event orientation consistance that described positioning unit 205 uses with the expression formula of scanning stack power maximization objective function is:

The expression formula that described micro-seismic event orientation consistance and scanning stack power maximize objective function is:

f(Azimuth,E)＝f1(Azimuth)→min&&f2(E)→max；

Wherein, $f1\left(\mathrm{Azimuth}\right)=\mathrm{abs}(\mathrm{Azimuth}(x_i,y_j,z_k,t_p)-\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}\mathrm{Azimuth}(x_i,y_j,z_k,t_p)/N);$ F2 (E)=E (x _i, y _j, z _k, t); X above _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _i∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; Azimuth is orientation function, and E is energy function, and t is compressional wave or shear wave time, t _pit is the compressional wave time; Min represents that f1 (Azimuth) function gets minimum value, max represents that f2 (E) function gets maximal value, and min & & f2 (E) represents that need meet function f 1 (Azimuth) gets minimum value sum functions f2 (E) and get maximal value simultaneously; → be mathematic sign, represent infinite approach.

Preferably, described first slip scan unit 203 comprises:

Three-dimensional coordinate slip scan module, in region of search, three-dimensional each coordinate axis of sliding carries out diffraction energy superposition; Wherein, the computing formula of diffraction energy superposition is:

$E(x_i,y_j,z_k)=\underset{t=t(x_i,y_j,z_k)-\mathrm{\τ}}{\overset{t=t(x_i,y_j,z_k)+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t(x_i,y_j,x_k)\right)$

Wherein, x _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _j∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; A is amplitude, and t is the function of time, and τ is time window, and E is energy function;

Signal time axle slip scan module, for according to micro-seismic event whilst on tour table sliding stack wave energy in length and breadth on signal time axle; Wherein, the expression formula superposing compressional wave energy is: the expression formula of superposition shear wave energy is: n is the number of wave detector, and E is energy function, and A is amplitude, t _pthe compressional wave time, t _sbe the shear wave time, t is the function of time.

Preferably, described second slip scan unit 204 specifically for:

According to micro-seismic event whilst on tour table, signal time axle slides and calculates the polarization angle of each time point micro-seismic event accordingly.Adopt the method for eigenwert-proper vector to calculate, the method has open source information to inquire about, and does not relate to new computing formula, and this step repeatedly adopts this algorithm to calculate in sliding process.

Preferably, described whilst on tour table is set up unit 202 and is set up micro-seismic event whilst on tour table based on correction rate model.

Following examples are described with reference to the accompanying drawings the technical program, and the present embodiment concrete steps are as follows:

1) microearthquake signals collecting continues to carry out, in magnanimity well, carry out microearthquake Signal analysis in Monitoring Data, and extracting micro-seismic event, as shown in Figure 3, is the ripple signal schematic representation in length and breadth of microearthquake signal in the present embodiment.Be extracted micro-seismic event signal segment in magnanimity micro-seismic monitoring data, microearthquake signal in length and breadth ripple signal is clearly.As shown in Figure 3, p wave interval velocity comparatively S wave velocity is fast, and P-wave comparatively S-wave first arrives wave detector, and S-wave is lower slightly compared with the frequency of P-wave.

2) by checking signal, fine correction is carried out to the rate pattern of target area, set up the micro-seismic event ripple whilst on tour table in length and breadth of scanning area based on correction rate model;

3) according to the RMS amplitude of microearthquake magnitude determinations signal, as shown in Figure 4, be the micro-seismic event RMS amplitude figure of the present embodiment; The synthesis of three-component (XYZ component) scalar is carried out to microearthquake signal, outstanding amplitude susceptibility, thus the amplitude information of outstanding microearthquake signal; Calculated amplitude energy Ratios on root mean square basis as shown in Figure 5, is the micro-seismic event amplitude ratio schematic diagram of the present embodiment.RMS amplitude basis is asked for per pass amplitude ratio, highlights the energy of micro-seismic event time of arrival.By micro-seismic event ripple whilst on tour table slip scan micro-seismic event diffraction stack energy in length and breadth on amplitude energy ratio; As shown in Figure 6, be the present embodiment wave diffraction energy scan result schematic diagram in length and breadth.According to micro-seismic event in length and breadth the in length and breadth wave diffraction energy of ripple whilst on tour table to three-dimensional search region carry out convergence superposition, Fig. 6 is the continuous multiple side view at depth direction.

4) according to the polaried orientation of micro-seismic event ripple whilst on tour table slip scan micro-seismic event in length and breadth; As shown in Figure 7, be the present embodiment polarization scans result schematic diagram.

5) in scanning area, according to step 3) the diffraction stack energy that obtains and step 4) the polaried orientation information that obtains, maximize objective function to micro-seismic event orientation consistance and scanning stack power to solve, optimum solution represents the most stable and consistent in orientation and the position that ripple scanning energy is the strongest in length and breadth in the direction in which, namely completes the space orientation of micro-seismic event.As shown in Figure 8, be micro-seismic event locus schematic diagram that the present embodiment is determined.Microearthquake located space position coordinates is: (595.2,5.3,3465.1).

Can be found by the present embodiment: the technical program can carry out autoscan location to the orientation of micro-seismic event and position automatically, and save artificial regulation time, on-the-spot real-time is better.

Above-described embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only the specific embodiment of the present invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (10)

1., based on a micro-seismic event localization method for wave energy autoscan in length and breadth, it is characterized in that, the method comprises:

Micro-seismic event identification is carried out to underground monitoring data;

Set up the micro-seismic event whilst on tour table of region of search;

According to micro-seismic event whilst on tour table, slip scan is carried out to the micro-seismic event identified, obtain diffraction stack energy;

According to micro-seismic event whilst on tour table, slip scan is carried out to the polaried orientation of the micro-seismic event identified, obtain polaried orientation information;

According to described diffraction stack energy and described polaried orientation information, maximize objective function to micro-seismic event orientation consistance and scanning stack power to solve, obtain the optimum solution of objective function, this optimum solution is expressed as: the most stable and consistent in orientation and the position that ripple scanning energy is the strongest in length and breadth in this orientation, and this position is the position of micro-seismic event.

2. the method for claim 1, is characterized in that, the expression formula that described micro-seismic event orientation consistance and scanning stack power maximize objective function is:

f(Azimuth,E)＝f1(Azimuth)→min&&f2(E)→max；

Wherein, $f1\left(\mathrm{Azimuth}\right)=\mathrm{abs}(\mathrm{Azimuth}(x_i,y_j,z_k,t_p)-\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}\mathrm{Azimuth}(x_i,y_j,z_k,t_p)/N);$ F2 (E)=E (x _i, y _j, z _k, t); X above _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _i∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; Azimuth is orientation function, and E is energy function, and t is compressional wave or shear wave time, t _pit is the compressional wave time; Min represents that f1 (Azimuth) function gets minimum value, max represents that f2 (E) function gets maximal value, and min & & f2 (E) represents that need meet function f 1 (Azimuth) gets minimum value sum functions f2 (E) and get maximal value simultaneously; → be mathematic sign, represent infinite approach.

3. method as claimed in claim 1 or 2, is characterized in that, describedly comprises the step that the micro-seismic event identified carries out slip scan according to micro-seismic event whilst on tour table:

In region of search, three-dimensional each coordinate axis of sliding carries out diffraction energy superposition; Wherein, the computing formula of diffraction energy superposition is:

$E(x_i,y_j,z_k)=\underset{t=t(x_i,y_j,z_k)-\mathrm{\τ}}{\overset{t=t(x_i,y_j,z_k)+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t(x_i,y_j,z_k)\right)$

Wherein, x _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _j∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; A is amplitude, and t is the function of time, and τ is time window, and E is energy function;

According to micro-seismic event whilst on tour table sliding stack wave energy in length and breadth on signal time axle; Wherein, the expression formula superposing compressional wave energy is: $E\left(p\right)=\underset{t=t_p-\mathrm{\τ}}{\overset{t=t_p+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t_p\right);$ The expression formula of superposition shear wave energy is: $E\left(s\right)=\underset{t=t_s-\mathrm{\τ}}{\overset{t=t_s+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t_s\right);$ N is the number of wave detector, and E is energy function, and A is amplitude, t _pthe compressional wave time, t _sbe the shear wave time, t is the function of time.

4. method as claimed in claim 1 or 2, is characterized in that, describedly comprises the step that the polaried orientation of the micro-seismic event identified carries out slip scan according to micro-seismic event whilst on tour table:

According to micro-seismic event whilst on tour table, signal time axle slides and calculates the polarization angle of each time point micro-seismic event accordingly.

5. method as claimed in claim 1 or 2, it is characterized in that, described micro-seismic event whilst on tour table is set up based on correction rate model.

6., based on a micro-seismic event locating device for wave energy autoscan in length and breadth, it is characterized in that, this device comprises:

Recognition unit, for carrying out micro-seismic event identification to underground monitoring data;

Unit set up by whilst on tour table, for setting up the micro-seismic event whilst on tour table of region of search;

First slip scan unit, for carrying out slip scan according to micro-seismic event whilst on tour table to the micro-seismic event identified, obtains diffraction stack energy;

Second slip scan unit, for carrying out slip scan according to micro-seismic event whilst on tour table to the polaried orientation of the micro-seismic event identified, obtains polaried orientation information;

Positioning unit, for according to described diffraction stack energy and described polaried orientation information, maximize objective function to micro-seismic event orientation consistance and scanning stack power to solve, obtain the optimum solution of objective function, this optimum solution is expressed as: the most stable and consistent in orientation and the position that ripple scanning energy is the strongest in length and breadth in this orientation, and this position is the position of micro-seismic event.

7. device as claimed in claim 6, is characterized in that, the expression formula that the micro-seismic event orientation consistance that described positioning unit uses and scanning stack power maximize objective function is:

f(Azimuth,E)＝f1(Azimuth)→min&&f2(E)→max；

Wherein, $f1\left(\mathrm{Azimuth}\right)=\mathrm{abs}(\mathrm{Azimuth}(x_i,y_j,z_k,t_p)-\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}\mathrm{Azimuth}(x_i,y_j,z_k,t_p)/N);$ F2 (E)=E (x _i, y _j, z _k, t); X above _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _i∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; Azimuth is orientation function, and E is energy function, and t is compressional wave or shear wave time, t _pit is the compressional wave time; Min represents that f1 (Azimuth) function gets minimum value, max represents that f2 (E) function gets maximal value, and min & & f2 (E) represents that need meet function f 1 (Azimuth) gets minimum value sum functions f2 (E) and get maximal value simultaneously; → be mathematic sign, represent infinite approach.

8. device as claimed in claims 6 or 7, it is characterized in that, described first slip scan unit comprises:

Three-dimensional coordinate slip scan module, in region of search, three-dimensional each coordinate axis of sliding carries out diffraction energy superposition; Wherein, the computing formula of diffraction energy superposition is:

$E(x_i,y_j,z_k)=\underset{t=t(x_i,y_j,z_k)-\mathrm{\τ}}{\overset{t=t(x_i,y_j,z_k)+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t(x_i,y_j,z_k)\right)$

Wherein, x _irepresent the x-axis coordinate figure of cartesian coordinate system, y _jrepresent the y-axis coordinate figure of cartesian coordinate system, z _krepresent the z-axis coordinate figure of cartesian coordinate system; x _i∈ [x _min, x _max], y _j∈ [y _min, y _max], z _k∈ [z _min, z _max], N is the number of wave detector, x _minrepresent the minimum value of x coordinate axis in region of search, x _maxrepresent the maximal value of x coordinate axis in region of search; y _minrepresent the minimum value of y coordinate axis in region of search, y _maxrepresent the maximal value of y coordinate axis in region of search; z _minrepresent the minimum value of z coordinate axle in region of search, z _maxrepresent the maximal value of z coordinate axle in region of search; A is amplitude, and t is the function of time, and τ is time window, and E is energy function;

Signal time axle slip scan module, for according to micro-seismic event whilst on tour table on signal time axle sliding stack in length and breadth wave energy wherein, the expression formula of superposition compressional wave energy is: $E\left(p\right)=\underset{t=t_p-\mathrm{\τ}}{\overset{t=t_p+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t_p\right);$ The expression formula of superposition shear wave energy is: $E\left(s\right)=\underset{t=t_s-\mathrm{\τ}}{\overset{t=t_s+\mathrm{\τ}}{\mathrm{\Σ}}}\underset{n=1}{\overset{n=N}{\mathrm{\Σ}}}{A}^2\left(t_s\right);$ N is the number of wave detector, and E is energy function, and A is amplitude, t _pthe compressional wave time, t _sbe the shear wave time, t is the function of time.

9. device as claimed in claims 6 or 7, is characterized in that, described second slip scan unit specifically for:

According to micro-seismic event whilst on tour table, signal time axle slides and calculates the polarization angle of each time point micro-seismic event accordingly.

10. device as claimed in claims 6 or 7, it is characterized in that, described whilst on tour table is set up unit and is set up micro-seismic event whilst on tour table based on correction rate model.