CN106249295A - A combined rapid positioning method and system for microseismic P and S waves in a borehole - 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. a microseismic P, S wave joint fast location method in a well, is characterized in that, comprises the steps: a. Acquire the microseismic data of the same event, perform vector decomposition on the microseismic data of the same event, and obtain the vertical azimuth angles of the two three-component geophones; b. According to Snell's law, establish the functional relationship between the vertical depth H and the horizontal distance L with the three-component detector as the starting point; c. Establish the P&S wave depth objective function related to the source point position. By scanning the depth, the objective function is minimized, and the vertical depth and horizontal distance between the source point and the three-component geophone are calculated to realize the joint rapid positioning of P and S waves. 2. according to claim 1, a kind of in-well microseismic P, S wave joint rapid positioning method is characterized in that, the specific realization of described step a is: The microseismic data of the same event are obtained, and the microseismic data of the same event are vector decomposed to calculate the vertical azimuths θ ₁ and θ ₂ of two three-component geophones for P waves and the vertical azimuth β ₂ of one geophone for S waves. 3. according to claim 2, a kind of in-well microseismic P, S wave joint rapid location method is characterized in that, the concrete realization of step b is: Based on the known velocity model of the formation medium, according to Snell’s law, starting from the three-component geophone, along the reverse direction of the wave propagation path, the functional relationship between the vertical depth H and the horizontal distance L starting from the three-component geophone is established, corresponding to There are two P wave functions L _P1 =f _P1 (H), L _P2 =f _P2 (H) and one S wave function L _S2 =f _S2 (H). 4. according to claim 3 described a kind of in-well microseismic P, S wave joint rapid location method, it is characterized in that, the concrete realization of described step c is: Set the objective function OPJ=‖f _P1 (H)-f _P2 (H)‖+‖f _P1 (H)-f _S2 (H)‖ related to the source position, and when OPJ→0, quickly calculate the optimal solution Depth H ₀ and corresponding horizontal distance L ₀ = (f _P1 (H ₀ )+f _P2 (H ₀ )+f _S2 (H _{0 )} )/3, that is, the vertical depth and horizontal distance. 5. according to claim 2 described a kind of in-well microseismic P, S wave joint rapid positioning method, it is characterized in that, described vector decomposition adopts end curve-histogram method to realize, specifically as follows: Set a time window, including the P-wave event, and perform polarization analysis on the amplitude values ( _xi , y _i ) of the X and Y components of the three-component detector; define the instantaneous energy E _i and the instantaneous orientation φ _i : ${\mathrm{<span\; class="notranslate">\; E.</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; x</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}^{\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; the\; y</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}^{\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; )</span>$ tgφ _i =y _i /x _i (2) First, according to the above formula, calculate the instantaneous orientation φ _i and instantaneous energy E _i of each sample point, and make a histogram of instantaneous energy E _i against instantaneous orientation φ _i ; then, draw the ( _xi , y _i ) amplitude value The line diagram of the coordinate system, that is, the arrow-end curve diagram, is linearly fitted to estimate the range of inclination; finally, refer to the range of inclination to find the position of the maximum instantaneous energy in the E _i and φ _i histograms. The corresponding angle φ is the horizontal azimuth angle of the three-component detector relative to the direction of wave propagation; Using the horizontal azimuth angle φ of the three-component geophone relative to the direction of wave propagation, the X and Y components are rotated to obtain the horizontal radial R and tangential T: R _i ＝x _i cos(φ)+y _i sin(φ) (3) T _i ＝-x _i sin(φ)+y _i cos(φ) Among them, R _i and T _i are the instantaneous amplitudes of horizontal radial and tangential components respectively after rotation; Similarly, conduct polarization analysis on the amplitude values (Z _i , R _i ) of the Z component and the horizontal radial R component of the detector, and draw the line diagram of the (Z _i , R _i ) amplitude value, that is, the sagittal graph, and carry out Linear fitting, estimate the range of inclination; refer to the range of inclination, find the position of the maximum instantaneous energy in the histogram, and find the corresponding vertical azimuth θ of the geophone relative to the direction of wave propagation, that is, complete the vector decomposition and calculate the vertical Azimuth; Perform vector decomposition on the P waves of any two three-component geophones in the observation well to obtain the corresponding vertical azimuths θ ₁ and θ ₂ , and at the same time, perform vector decomposition on the S wave of one of the geophones to obtain the vertical azimuth β ₂ . 6. according to claim 3 a kind of in-well microseismic P, S wave joint rapid positioning method, it is characterized in that, the concrete realization of setting up functional relationship is as follows: Assuming that the vertical distance between the two geophones is h ₀ , given the distance H from the second geophone to the ground, according to the acoustic logging, it is divided into N layers of compressional and shear wave velocity layers from shallow to deep [Vp ₁ , Vp ₂ ,. .....,Vp _N-1 ,Vp _N ], [Vs ₁ ,Vs ₂ ,......,Vs _N-1 ,Vs _N ], and the depth H is also divided into corresponding N layers [H ₁ ,H ₂ ,...,H _N-1 ,H _N ]; According to Snell's law, the horizontal distance related to the depth H is obtained from two P-wave vertical azimuths θ ₁ and θ ₂ and one S-wave vertical azimuth β ₂ : ${\mathrm{<span\; class="notranslate">\; L</span>}}_{\mathrm{<span\; class="notranslate">\; P</span>}\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; f</span>}}_{\mathrm{<span\; class="notranslate">\; P</span>}\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; h</span>}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span><span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; 0</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; \&CenterDot;</span>\mathrm{<span\; class="notranslate">\; the\; tan</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 2</span>}}{\overset{\mathrm{<span\; class="notranslate">\; N</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; \&CenterDot;</span>\frac{{\mathrm{<span\; class="notranslate">\; Vp</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; /</span>{\mathrm{<span\; class="notranslate">\; Vp</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; \&CenterDot;</span>\mathrm{<span\; class="notranslate">\; sin</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; )</span>}{\sqrt{\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; -</span>{<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; Vp</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; /</span>{\mathrm{<span\; class="notranslate">\; Vp</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; \&CenterDot;</span>\mathrm{<span\; class="notranslate">\; sin</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; )</span>}^{\mathrm{<span\; class="notranslate">\; 2</span>}}}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 4</span>}<span\; class="notranslate">\; )</span>$ ${\mathrm{<span\; class="notranslate">\; L</span>}}_{\mathrm{<span\; class="notranslate">\; P</span>}\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; f</span>}}_{\mathrm{<span\; class="notranslate">\; P</span>}\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; h</span>}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; N</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; \&Center\; Dot;</span>\frac{{\mathrm{<span\; class="notranslate">\; Vp</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; /</span>{\mathrm{<span\; class="notranslate">\; Vp</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; \&CenterDot;</span>\mathrm{<span\; class="notranslate">\; the\; s</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; )</span>}{\sqrt{\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; -</span>{<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; Vp</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; /</span>{\mathrm{<span\; class="notranslate">\; Vp</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; \&Center\; Dot;</span>\mathrm{<span\; class="notranslate">\; the\; s</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&theta;</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; )</span>}^{\mathrm{<span\; class="notranslate">\; 2</span>}}}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 5</span>}<span\; class="notranslate">\; )</span>$ ${\mathrm{<span\; class="notranslate">\; L</span>}}_{\mathrm{<span\; class="notranslate">\; S</span>}\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; =</span>{\mathrm{<span\; class="notranslate">\; f</span>}}_{\mathrm{<span\; class="notranslate">\; S</span>}\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; h</span>}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; N</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; h</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; \&Center\; Dot;</span>\frac{{\mathrm{<span\; class="notranslate">\; vs.</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; \&Center\; Dot;</span><span\; class="notranslate">\; /</span>{\mathrm{<span\; class="notranslate">\; vs.</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; \&Center\; Dot;</span>\mathrm{<span\; class="notranslate">\; the\; s</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&beta;</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; )</span>}{\sqrt{\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; -</span>{<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; vs.</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; \&Center\; Dot;</span><span\; class="notranslate">\; /</span>{\mathrm{<span\; class="notranslate">\; vs.</span>}}_{\mathrm{<span\; class="notranslate">\; 1</span>}}<span\; class="notranslate">\; \&Center\; Dot;</span>\mathrm{<span\; class="notranslate">\; the\; s</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; \&beta;</span>}}_{\mathrm{<span\; class="notranslate">\; 2</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; )</span>}^{\mathrm{<span\; class="notranslate">\; 2</span>}}}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 6</span>}<span\; class="notranslate">\; )</span>$ Among them, L _P1 and L _P2 are the radial and horizontal distances of the P wave propagation distance corresponding to the arbitrary depth H of the shallow and deep three-component geophones respectively, and L _S2 is the S wave propagation distance corresponding to the arbitrary depth H of the deep geophone radial horizontal distance. 7. according to claim 4 described a kind of microseismic P in the well, the S wave joint rapid location method, it is characterized in that, the specific realization of calculating source point relative three-component geophone vertical depth and horizontal distance is: Establish the objective function related to the source position, quickly calculate the optimal solution depth H ₀ and the corresponding horizontal distance L ₀ =(f _P1 (H ₀ )+f _P2 (H ₀ )+f _S2 (H ₀ ))/3, To achieve joint positioning of P and S waves, the objective 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, it means that the P-wave horizontal distance of one geophone intersects with the P-wave and S-wave horizontal distances of another geophone at the same time, and the position of the intersection point is the source position. 8. according to claim 7 described a kind of microseismic P in the well, the joint rapid positioning method of S wave, it is characterized in that, calculate optimal solution depth H fast ₀ and corresponding horizontal distance L ₀ adopt Newton's dichotomy method, concrete realization is as follows: Quickly find two depths H ₁ and H ₂ so that OPJ(H ₁ )·OPJ(H ₂ )≤0; Select the middle value H ₃ =(H ₁ +H ₂ )/2, compare OPJ(H ₁ ), OPJ(H ₂ ), OPJ(H ₃ ), if OPJ(H ₁ )·OPJ(H ₃ )≤0, Then set the new depth combination H ₁ =H ₁ , H ₂ =H ₃ , if OPJ(H ₂ )·OPJ(H ₃ )≤0, then set the new depth combination H ₁ =H ₂ , H ₂ =H ₃ ; Recalculate the intermediate value objective function, and then compare, and iterate until the convergence error is finally met, and the optimal solution H ₀ =(H ₁ +H ₂ )/2 is obtained, which is substituted into the objective function to calculate the optimal horizontal distance L ₀ = (f _P1 (H ₀ )+f _P2 (H ₀ )+f _S2 (H ₀ ))/3, realizing joint fast positioning of P and S waves. 9. A microseismic P and S wave joint rapid positioning system in a well, characterized in that it includes a data acquisition module, a vector decomposition module, a function establishment module and a result calculation module; The data acquisition module is used to obtain microseismic data of the same event; The vector decomposition module is used for vector decomposition of the microseismic data of the same event to obtain the vertical azimuths of two three-component geophones; Described function establishment module, it is used for according to Snell's law, establishes the functional relation between the vertical depth H of starting point and the horizontal distance L of three-component detector; The result calculation module is used to establish the P&S wave depth objective function related to the position of the seismic source point. By scanning the depth, the objective function is minimized, and the vertical depth and horizontal distance of the seismic source point relative to the three-component geophone are calculated to realize P and S waves. Joint rapid positioning.