CN106249295A - A kind of borehole microseismic P, S ripple associating method for rapidly positioning and system - Google Patents

Abstract

The present invention relates to a kind of borehole microseismic P, S ripple associating method for rapidly positioning and system, described method includes: a, obtain same event microseism data, same event microseism data is carried out resolution of vectors, it is thus achieved that the vertical azimuth of two three-component geophones；B, according to snell law, set up the functional relationship between vertical depth H and horizontal range L with three-component geophone as starting point；C, foundation P&S depth of convolution degree object function relevant with focal point position, pass through scan depths so that object function is minimum, calculates focal point relative to three-component geophone vertical depth and horizontal range, it is achieved P, S ripple associating quickly location.The present invention is just with two vertical azimuths of cymoscope, according to snell law, sets up focal point position relevant P&S depth of convolution degree object function, passes through scan depths so that object function is minimum, can be achieved with P, S ripple associating quickly location.

Description

A kind of borehole microseismic P, S ripple associating method for rapidly positioning and system

Technical field

The present invention relates to borehole microseismic localization process field, particularly relate to a kind of borehole microseismic P, S ripple associating method for rapidly positioning and system.

Background technology

Microseism Fracturing Monitoring technology is one of key technology in non-traditional dense sandstone gas, shale gas reservoir reservoir oil gas field development, the seismic source information oriented according to inverting, crack attribute (principal stress trend, fracture width, density etc.) can be obtained, it is used for evaluating fracturing effect, analyze crack induction rule, optimize well spacing etc..Therefore, in microseism signal processing, final purpose is seismic source location, also known as the most crucial technology of microseism signal processing.

Borehole microseismic monitoring is one of microseism observed pattern, and feature is that down-hole three-component receives, and microseism data is higher, and the focus number received is abundanter with type.At present, borehole microseismic location technology mainly has two kinds of thinkings:

One is just to drill based on P ripple, S ripple event whilst on tour, represent algorithm and have network searching method, simulated annealing, geiger method etc., advantage is easily to realize, and shortcoming is owing to the first arrival phase signal weak P of causing ripple, S ripple event whilst on tour are difficult to accurately pick up, and affects positioning result；

The second location thinking is based on wave equation convolution, and representing algorithm has interferometric method, reverse-time migration method, passive source imaging method, and advantage is to need not pickup event first arrival, and shortcoming is to data signal to noise ratio, rate pattern requirement height, calculates cost high.

Summary of the invention

The technical problem to be solved is for the deficiencies in the prior art, it is provided that a kind of borehole microseismic P, S ripple associating method for rapidly positioning and system.

The technical scheme is that a kind of borehole microseismic P, S ripple associating method for rapidly positioning, comprise the steps:

A, obtain same event microseism data, same event microseism data are carried out resolution of vectors, it is thus achieved that the vertical azimuth of two three-component geophones；

B, according to snell law, set up the functional relationship between vertical depth H and horizontal range L that three-component geophone is starting point；

C, foundation P&S depth of convolution degree object function relevant with focal point position, pass through scan depths so that object function is minimum, calculates focal point relative to three-component geophone vertical depth and horizontal range, it is achieved P, S ripple associating quickly location.

The invention has the beneficial effects as follows: the invention provides a kind of simple borehole microseismic event recognition method based on resolution of vectors, just with two vertical azimuths of cymoscope, according to snell law, set up focal point position relevant P&S depth of convolution degree object function, pass through scan depths, make object function minimum, can be achieved with P, S ripple associating quickly location.

Another technical scheme that the present invention solves above-mentioned technical problem is as follows: a kind of borehole microseismic P, S ripple associating quick positioning system, sets up module and result computing module including data acquisition module, resolution of vectors module, function；

Described data acquisition module, it is used for obtaining same event microseism data；

Described resolution of vectors module, it is for carrying out resolution of vectors to same event microseism data, it is thus achieved that the vertical azimuth of two three-component geophones；

Described function sets up module, and it, for according to snell law, sets up the functional relationship between vertical depth H and horizontal range L that three-component geophone is starting point；

Described result computing module, it, for setting up P&S depth of convolution degree object function relevant with focal point position, passes through scan depths so that object function is minimum, calculate focal point relative to three-component geophone vertical depth and horizontal range, it is achieved P, S ripple associating quickly location.

Accompanying drawing explanation

Fig. 1-1 is that the present invention realizes borehole microseismic P, S ripple rapid positioning operation flow chart；

Fig. 1-2 is that the present invention realizes borehole microseismic P, S ripple rapid positioning operation system block diagram；

Fig. 2 is borehole microseismic Vertical Observation schematic diagram: multiple focus, two cymoscopes ▲；

Fig. 3-1 is the X-component in borehole microseismic model three component seismic data；

Fig. 3-2 is the Y-component in borehole microseismic model three component seismic data；

Fig. 3-3 is the Z component in borehole microseismic model three component seismic data；

Fig. 4-1 is the P component in borehole microseismic model data resolution of vectors；

Fig. 4-2 is the S component in borehole microseismic model data resolution of vectors；

Fig. 5 is the vertical azimuth list obtained after all focal points resolution of vectors of noiselessness " X " type；

Fig. 6 is p-and s-wave velocity model list；

Fig. 7 is that noiselessness " X " type all focal points P, S ripple quick positioning result schematic diagram of associating is shown by the present invention；

Fig. 8 is noiselessness " X " type model points P, the S ripple quick position error rectangular histogram of associating (abscissa is focus sequence number, and vertical coordinate is inverting value and actual value error)；

Fig. 9-1 is that borehole microseismic model Z component data are increased random noise；

Fig. 9-2 is that borehole microseismic model X-component data are increased random noise；

Fig. 9-3 is that borehole microseismic model Y-component data are increased random noise；

Figure 10-1 is containing the P component in noise model resolution of vectors；

Figure 10-2 is containing the S component in noise model resolution of vectors；

Figure 11 is the vertical azimuth list obtained after the model containing noise has focal point resolution of vectors；

Figure 12 is that the model containing noise is had focal point P, S ripple quick positioning result schematic diagram of associating to show by the present invention；

Figure 13 is to combine quick position error rectangular histogram containing noise model points P, S ripple.

Detailed description of the invention

Being described principle and the feature of the present invention below in conjunction with accompanying drawing, example is served only for explaining the present invention, is not intended to limit the scope of the present invention.

Below in conjunction with the accompanying drawings the present invention is described in further detail:

Borehole microseismic three-component pretreatment includes: polarization analysis, rotation, orientation Concordance, resolution of vectors, denoising etc..Wherein, resolution of vectors can regard " with joint efforts and component " relation as, and by three-component vector signal, " with joint efforts " becomes along direction of wave travel P ripple or vertical direction of wave travel S ripple scalar information, produces vertical azimuth simultaneously.The present invention utilizes based on the vertical azimuth of resolution of vectors just, it is intended to utilizes different cymoscope P ripple, S ripple propagation rays path geometry common factor mode, realizes seismic source location.

Borehole microseismic three component signal, there is the features such as pressure break hypocentral location is unknown, focus seismic phase type is unknown, P/S ripple event mixes, its resolution of vectors process is exactly that Vector Message on three components " merges " one scalar of one-tenth, and numerous P ripples, S ripple event is separated.Meanwhile, this process can produce the vertical azimuth relevant with direction of wave travel, and the present invention utilizes this feature just, in conjunction with snell law, along the different reverse directions of cymoscope propagation path, set up the target equation relevant with focus, rapid solving, it is achieved P, S ripple co-located.

Assume initially that, focal point excites P ripple and S ripple signal, before reaching inspection well three-component geophone, still belong to scalar wave, after three-component X, Y, Z cymoscope receives signal, become three direction vector waves perpendicular to each other, and resolution of vectors can regard its inverse process as, i.e. made a concerted effort R by X, Y-component computational geometry and tangentially remain T, calculate Z component again and R geometry is made a concerted effort P and tangential S, producing two attitudes simultaneously: horizontal azimuth, vertical azimuth.Wherein, in physical significance, vertical azimuth represents the angle of cut between direction of wave travel and cymoscope Z component vertically downward.

As Figure 1-1, a kind of borehole microseismic P, S ripple associating method for rapidly positioning, comprise the steps:

Same event microseism data are carried out resolution of vectors by the first step: obtain same event microseism data, calculate two three-component geophone P ripple Vertical Square parallactic angles θ₁、θ₂With a cymoscope vertical azimuthal angle beta of S ripple₂；

Specifically, first borehole microseismic two cymoscopes three-component X, Y, Z are comprised P, S ripple event data and carry out resolution of vectors (wherein, two cymoscope three-components receive combines seismic source information from same P&S ripple, show as in data two groups of three-components have received represent P ripple, the impulse waveform data of S ripple event), obtain two groups of vertical azimuth angle theta of P ripple along direction of wave travel P ripple, vertical direction of wave travel S ripple, and two cymoscopes₁、θ₂With a cymoscope vertical azimuthal angle beta of S ripple₂；

Wherein, described resolution of vectors uses hodograph-histogram method to realize, specific as follows:

Window when arranging one, contains P ripple event, to three-component geophone X, Y-component amplitude (x_i,y_i) make polarization analysis；Definition instantaneous energy E_iWith transient bearing φ_i:

$E_i=x_i^2+y_i^2---\left(1\right)$

tgφ_i=y_i/x_i (2)

First, according to above formula, calculate each sampling point transient bearing φ_iWith instantaneous energy E_i, make instantaneous energy E simultaneously_iTo transient bearing φ_iRectangular histogram；Then, (x is drawn_i,y_i) amplitude coordinate system line graph, i.e. hodograph figure, it is carried out linear fit, estimates inclination angle scope；Finally, with reference to inclination angle scope, at E_i、φ_iRectangular histogram is found the position of instantaneous energy maximum, the most corresponding angle φ, is three-component geophone relative to direction of wave travel horizontal azimuth；

Utilize three-component geophone relative to direction of wave travel horizontal azimuth φ, X, Y-component carried out rotation processing, it is thus achieved that horizontal radial R and tangential T:

R_i=x_i cos(φ)+y_i sin(φ)

(3)

T_i=-x_i sin(φ)+y_i cos(φ)

Wherein, R_i、T_iIt is respectively horizontal radial, tangential component instantaneous amplitude after rotating；

Equally, to cymoscope Z component, horizontal radial R component amplitude (Z_i,R_i) make polarization analysis, draw (Z_i,R_i) amplitude line graph, i.e. hodograph figure, it is carried out linear fit, estimates inclination angle scope；With reference to inclination angle scope, rectangular histogram is found the position of instantaneous energy maximum, obtain the cymoscope of its correspondence relative to the vertical azimuth angle theta of direction of wave travel, i.e. complete resolution of vectors and calculate vertical azimuth；

Inspection well any two three-component geophone P ripple is carried out resolution of vectors, it is thus achieved that corresponding vertical azimuth is θ₁、θ₂, meanwhile, one of them cymoscope S ripple is carried out resolution of vectors, it is thus achieved that vertical azimuthal angle beta₂。

Second step: by known formation medium velocity model, according to snell law, from three-component geophone, along propagation path reverse direction, set up the functional relationship between vertical depth H and horizontal range L with three-component geophone as starting point, to there being two P wave functions L_P1=f_P1(H)、L_P2=f_P2(H) with S wave function L_S2=f_S2(H)；

Specifically, being implemented as follows of functional relationship is set up:

Assume that between two cymoscopes, vertical dimension is h₀, provide second any depth H in cymoscope underground of distance, according to acoustic logging, be divided into from shallow to deep n layer p-and s-wave velocity layer [Vp₁,Vp₂,......,Vp_N-1,Vp_N]、[Vs₁,Vs₂,......,Vs_N-1,Vs_N], and depth H is also divided into corresponding N shell [H₁,H₂,......,H_N-1,H_N]；

According to snell law, by two vertical azimuth angle theta of P ripple₁、θ₂With a vertical azimuthal angle beta of S ripple₂, obtain the horizontal range relevant with depth H:

$L_{P1}=f_{P1}\left(H\right)=(H_1+h_0)\·tan\left({\mathrm{\θ}}_1\right)+\underset{i=2}{\overset{N}{\Σ}}{H}_i\·\frac{{\mathrm{Vp}}_i/{\mathrm{Vp}}_1\·sin\left({\mathrm{\θ}}_1\right)}{\sqrt{1-{({\mathrm{Vp}}_i/{\mathrm{Vp}}_1\·\sin ({\mathrm{\θ}}_1\left)\right)}^2}}---\left(4\right)$

$L_{P2}=f_{P2}\left(H\right)=\underset{i=1}{\overset{N}{\Σ}}{H}_i\·\frac{{\mathrm{Vp}}_i/{\mathrm{Vp}}_1\·sin\left({\mathrm{\θ}}_2\right)}{\sqrt{1-{({\mathrm{Vp}}_i/{\mathrm{Vp}}_1\·sin({\mathrm{\θ}}_2\left)\right)}^2}}---\left(5\right)$

$L_{S2}=f_{S2}\left(H\right)=\underset{i=1}{\overset{N}{\Σ}}{H}_i\·\frac{{\mathrm{Vs}}_i/{\mathrm{Vs}}_1\·sin\left({\mathrm{\β}}_2\right)}{\sqrt{1-{({\mathrm{Vs}}_i/{\mathrm{Vs}}_1\·sin({\mathrm{\β}}_2\left)\right)}^2}}---\left(6\right)$

Wherein, L_P1、L_P2It is respectively the P-wave distance cymoscope radial level distance that shallow-layer, any depth H of deep layer three-component geophone are corresponding, L_S2For the S ripple propagation distance cymoscope radial level distance that any depth H of deep layer cymoscope is corresponding.

3rd step: the object function OPJ=relevant with hypocentral location is set | | f_P1(H)-f_P2(H)||+||f_P1(H)-f_S2(H) | |, when OPJ → 0, quickly calculate optimal solution depth H₀With respective horizontal distance L₀=(f_P1(H₀)+f_P2(H₀)+f_S2(H₀))/3, i.e. required is that focal point is relative to three-component geophone vertical depth and horizontal range.

Specifically, focal point being implemented as relative to three-component geophone vertical depth and horizontal range is calculated:

Set up the object function relevant with hypocentral location, quickly calculate optimal solution depth H₀With respective horizontal distance L₀=(f_P1(H₀)+f_P2(H₀)+f_S2(H₀))/3, it is achieved P, S ripple co-located, object function is as follows:

OPJ (H)=(L_P1-L_P2)+(L_P1-L_S2)

(7)

=(f_P1(H)-f_P2(H))+(f_P1(H)-f_S2(H))

When OPJ → 0, representing that a cymoscope P ripple horizontal range intersects with another cymoscope P ripple, S ripple horizontal range simultaneously, this position of intersecting point is hypocentral location.

Wherein, optimal solution depth H is quickly calculated₀With respective horizontal distance L₀Use Newton's dichotomy, be implemented as follows:

Rapidly find out two depth H₁、H₂So that OPJ (H₁)·OPJ(H₂)≤0；

Choose intermediate value H₃=(H₁+H₂)/2, compare OPJ (H₁)、OPJ(H₂)、OPJ(H₃), if OPJ is (H₁)·OPJ(H₃)≤0, then make new degree of depth combination H₁=H₁、H₂=H₃If, OPJ (H₂)·OPJ(H₃)≤0, then make new degree of depth combination H₁=H₂、H₂=H₃；

Recalculate intermediate value object function, then compare, such iteration, until finally meeting convergence error, obtain optimal solution H₀=(H₁+H₂)/2, substitute into object function, calculate optimal level distance L₀=(f_P1(H₀)+f_P2(H₀)+f_S2(H₀))/3, it is achieved P, S ripple associating quickly location.

As shown in Figure 1-2, a kind of borehole microseismic P, S ripple associating quick positioning system, set up module and result computing module including data acquisition module, resolution of vectors module, function；

Described data acquisition module, it is used for obtaining same event microseism data；

Described resolution of vectors module, it is for carrying out resolution of vectors to same event microseism data, it is thus achieved that the vertical azimuth of two three-component geophones；

Described function sets up module, and it, for according to snell law, sets up the functional relationship between vertical depth H and horizontal range L that three-component geophone is starting point；

Described result computing module, it, for setting up P&S depth of convolution degree object function relevant with focal point position, passes through scan depths so that object function is minimum, calculate focal point relative to three-component geophone vertical depth and horizontal range, it is achieved P, S ripple associating quickly location.

P, S ripple face of the present invention quickly speed localization process effect is verified below by borehole microseismic model data.

As in figure 2 it is shown, be made up of two three-component geophones and one group of " X " type focal point.

Fig. 3-1 to Fig. 3-3 is borehole microseismic three component seismic data at the center focal point of simulation, comprises one group of P, S ripple information, excites P ripple, S ripple the most simultaneously.

Fig. 4-1 and Fig. 4-2 is that X, Y, Z three-component in Fig. 3-1 to Fig. 3-3 is carried out resolution of vectors result.It can be seen that in figure, P component mainly comprises P ripple information, S component mainly comprises S ripple information, it is achieved that propagate and vertical transmission direction scalar P ripple, S ripple to along ripple from original vector.All focal points are carried out resolution of vectors, in addition to obtaining respective P ripple, S ripple information, obtains one group of vertical azimuth, as shown in Figure 5 simultaneously.

Then, by known speed model (can essentially obtain from acoustic logging) (participation Fig. 6), arranging deeper cymoscope is initial point, according to formula (4), (5), (6), the calculating degree of depth is from the beginning of H=0, and corresponding cymoscope P ripple, S ripple horizontal range, further according to formula (7), calculating three relevant H function intersections of complex curve, the corresponding degree of depth, horizontal range are required focal point position.Fig. 7 is for calculate all focal point position views by the inventive method, and Fig. 8 is then corresponding rectangular histogram by mistake, it can be seen that the biggest number focal point position error, within 1 meter, illustrates the feasibility of seismic source location of the present invention.

Finally, increase random noise to model, if Fig. 9-1 is to shown in Fig. 9-3, analyze noise to seismic source location result sensitivity of the present invention.Equally, according to flow chart 1, resolution of vectors (such as Figure 10-1 and Figure 10-2) is carried out to containing noise data, it is thus achieved that vertical azimuth (such as Figure 11).Further according to vertical azimuth, carry out P, S ripple co-located, be finally inversed by all focal point positions, as shown in figure 12.And Figure 13 is Figure 12 inversion result and actual position histogram of error, it can be seen that maximum error is, but most of error is all within 3 meters, illustrates that the present invention exists certain noise immunity to noise, has a certain degree of practical value.

In sum, the present invention proposes " with joint efforts " with " separation " it is assumed that complete borehole microseismic three-component resolution of vectors, it is thus achieved that key, represent P ripple, the vertical azimuth of S direction of wave travel and the cymoscope Z component angle of cut.Utilizing this azimuth, in conjunction with snell law, set up the target equation relevant with hypocentral location, by Newton's dichotomy, fast inversion goes out focal point position, it is achieved borehole microseismic P&S ripple co-located processes.Model data is tested, and illustrates that localization process of the present invention is simple, easy, and simultaneously to noise-sensitive analysis, also the explanation present invention has to a certain degree noise immunity, processes offer technical support real-time for borehole microseismic field condition.

The foregoing is only presently preferred embodiments of the present invention, not in order to limit the present invention, all within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. made, should be included within the scope of the present invention.

Claims (9)

1. borehole microseismic P, S ripple associating method for rapidly positioning, it is characterised in that include as follows Step:

A, obtain same event microseism data, same event microseism data carried out resolution of vectors, Obtain the vertical azimuth of two three-component geophones；

B, according to snell law, set up vertical depth H and horizontal range that three-component geophone is starting point Functional relationship between L；

C, foundation P&S depth of convolution degree object function relevant with focal point position, pass through scan depths so that Object function is minimum, calculates focal point relative to three-component geophone vertical depth and horizontal range, it is achieved P, S ripple associating quickly location.

A kind of borehole microseismic P, S ripple associating method for rapidly positioning, It is characterized in that, being implemented as of described step a:

Obtain same event microseism data, same event microseism data are carried out resolution of vectors, calculate Go out two three-component geophone P ripple Vertical Square parallactic angles θ₁、θ₂With a cymoscope vertical azimuthal angle beta of S ripple₂。

A kind of borehole microseismic P, S ripple associating method for rapidly positioning, It is characterized in that, being implemented as of step b:

By known formation medium velocity model, according to snell law, from three-component geophone, edge Propagation path reverse direction, sets up the vertical depth H with three-component geophone as starting point and horizontal range L Between functional relationship, to there being two P wave functions L_P1=f_P1(H)、L_P2=f_P2(H) with a S ripple Function L_S2=f_S2(H)。

A kind of borehole microseismic P, S ripple associating method for rapidly positioning, It is characterized in that, being implemented as of described step c:

The object function relevant with hypocentral location is set OPJ=‖ f_P1(H)-f_P2(H)‖+‖f_P1(H)-f_S2(H) ‖, when OPJ → 0, quickly calculates optimal solution deep Degree H₀With respective horizontal distance L₀=(f_P1(H₀)+f_P2(H₀)+f_S2(H₀₎)/3, i.e. required is that focal point is relative Three-component geophone vertical depth and horizontal range.

A kind of borehole microseismic P, S ripple associating method for rapidly positioning, It is characterized in that, described resolution of vectors uses hodograph-histogram method to realize, specific as follows:

Window when arranging one, contains P ripple event, to three-component geophone X, Y-component amplitude (x_i,y_i) make polarization analysis；Definition instantaneous energy E_iWith transient bearing φ_i:

$E_i=x_i^2+y_i^2---\left(1\right)$

tgφ_i=y_i/x_i (2)

First, according to above formula, calculate each sampling point transient bearing φ_iWith instantaneous energy E_i, make simultaneously Instantaneous energy E_iTo transient bearing φ_iRectangular histogram；Then, (x is drawn_i,y_i) amplitude coordinate system line graph, I.e. hodograph figure, carries out linear fit to it, estimates inclination angle scope；Finally, with reference to inclination angle scope, At E_i、φ_iThe position of searching instantaneous energy maximum in rectangular histogram, the most corresponding angle φ, i.e. For three-component geophone relative to direction of wave travel horizontal azimuth；

Utilize three-component geophone relative to direction of wave travel horizontal azimuth φ, X, Y-component are revolved Turn and process, it is thus achieved that horizontal radial R and tangential T:

R_i=x_icos(φ)+y_isin(φ) (3)

T_i=-x_i sin(φ)+y_i cos(φ)

Wherein, R_i、T_iIt is respectively horizontal radial, tangential component instantaneous amplitude after rotating；

Equally, to cymoscope Z component, horizontal radial R component amplitude (Z_i,R_i) make polarization analysis, draw Go out (Z_i,R_i) amplitude line graph, i.e. hodograph figure, it is carried out linear fit, estimates inclination angle model Enclose；With reference to inclination angle scope, rectangular histogram is found the position of instantaneous energy maximum, obtains its correspondence Cymoscope, relative to the vertical azimuth angle theta of direction of wave travel, i.e. completes resolution of vectors and calculates vertical azimuth；

Inspection well any two three-component geophone P ripple is carried out resolution of vectors, it is thus achieved that corresponding vertical orientation Angle is θ₁、θ₂, meanwhile, one of them cymoscope S ripple is carried out resolution of vectors, it is thus achieved that vertical azimuth β₂。

A kind of borehole microseismic P, S ripple associating method for rapidly positioning, It is characterized in that, set up being implemented as follows of functional relationship:

Assume that between two cymoscopes, vertical dimension is h₀, provide second cymoscope underground of distance the deepest Degree H, according to acoustic logging, is divided into from shallow to deep n layer p-and s-wave velocity layer [Vp₁,Vp₂,......,Vp_N-1,Vp_N]、 [Vs₁,Vs₂,......,Vs_N-1,Vs_N], and depth H is also divided into corresponding N shell [H₁,H₂,......,H_N-1,H_N]；

According to snell law, by two vertical azimuth angle theta of P ripple₁、θ₂With a vertical azimuthal angle beta of S ripple₂, Obtain the horizontal range relevant with depth H:

$L_{P1}=f_{P1}\left(H\right)=(H_1+h_0)\·\tan \left({\mathrm{\θ}}_1\right)+\underset{i=2}{\overset{N}{\Σ}}{H}_i\·\frac{{\mathrm{Vp}}_i/{\mathrm{Vp}}_1\·\sin \left({\mathrm{\θ}}_1\right)}{\sqrt{1-{({\mathrm{Vp}}_i/{\mathrm{Vp}}_1\·\sin \left({\mathrm{\θ}}_1\right))}^2}}---\left(4\right)$

$L_{P2}=f_{P2}\left(H\right)=\underset{i=1}{\overset{N}{\Σ}}{H}_i\·\frac{{\mathrm{Vp}}_i/{\mathrm{Vp}}_1\·sin\left({\mathrm{\θ}}_2\right)}{\sqrt{1-{({\mathrm{Vp}}_i/{\mathrm{Vp}}_1\·sin({\mathrm{\θ}}_2\left)\right)}^2}}---\left(5\right)$

$L_{S2}=f_{S2}\left(H\right)=\underset{i=1}{\overset{N}{\Σ}}{H}_{\mathrm{i}}\·\frac{{\mathrm{Vs}}_{\mathrm{i}}\·/{\mathrm{Vs}}_1\·sin\left({\mathrm{\β}}_2\right)}{\sqrt{1-{({\mathrm{Vs}}_{\mathrm{i}}\·/{\mathrm{Vs}}_1\·sin({\mathrm{\β}}_2\left)\right)}^2}}---\left(6\right)$

Wherein, L_P1、L_P2It is respectively the P ripple that shallow-layer, any depth H of deep layer three-component geophone are corresponding Propagation distance cymoscope radial level distance, L_S2The S ripple corresponding for any depth H of deep layer cymoscope passes Broadcast distance cymoscope radial level distance.

A kind of borehole microseismic P, S ripple associating method for rapidly positioning, It is characterized in that, calculate focal point implementing relative to three-component geophone vertical depth and horizontal range For:

Set up the object function relevant with hypocentral location, quickly calculate optimal solution depth H₀With respective horizontal Distance L₀=(f_P1(H₀)+f_P2(H₀)+f_S2(H₀))/3, it is achieved P, S ripple co-located, object function is such as Under:

OPJ (H)=(L_P1-L_P2)+(L_P1-L_S2) (7)

=(f_P1(H)-f_P2(H))+(f_P1(H)-f_S2(H))

When OPJ → 0, represent cymoscope P ripple horizontal range simultaneously with another cymoscope P ripple, S ripple horizontal range intersects, and this position of intersecting point is hypocentral location.

A kind of borehole microseismic P, S ripple associating method for rapidly positioning, It is characterized in that, quickly calculate optimal solution depth H₀With respective horizontal distance L₀Use Newton's dichotomy, It is implemented as follows:

Rapidly find out two depth H₁、H₂So that OPJ (H₁)·OPJ(H₂)≤0；

Choose intermediate value H₃=(H₁+H₂)/2, compare OPJ (H₁)、OPJ(H₂)、OPJ(H₃), if OPJ(H₁)·OPJ(H₃)≤0, then make new degree of depth combination H₁=H₁、H₂=H₃If, OPJ(H₂)·OPJ(H₃)≤0, then make new degree of depth combination H₁=H₂、H₂=H₃；

Recalculating intermediate value object function, then compare, such iteration, until finally meeting convergence Error, obtains optimal solution H₀=(H₁+H₂)/2, substitute into object function, calculate optimal level distance L₀=(f_P1(H₀)+f_P2(H₀)+f_S2(H₀))/3, it is achieved P, S ripple associating quickly location.

9. borehole microseismic P, S ripple associating quick positioning system, it is characterised in that include data Acquisition module, resolution of vectors module, function set up module and result computing module；

Described data acquisition module, it is used for obtaining same event microseism data；

Described resolution of vectors module, it is for carrying out resolution of vectors to same event microseism data, it is thus achieved that The vertical azimuth of two three-component geophones；

Described function sets up module, and it is for according to snell law, and setting up three-component geophone is starting point Vertical depth H and horizontal range L between functional relationship；

Described result computing module, it is for setting up P&S ripple depth targets letter relevant with focal point position Number, passes through scan depths so that object function is minimum, calculates focal point and hangs down relative to three-component geophone The straight degree of depth and horizontal range, it is achieved P, S ripple associating quickly location.