CN106249295A - A combined rapid positioning method and system for microseismic P and S waves in a borehole - Google Patents

A combined rapid positioning method and system for microseismic P and S waves in a borehole Download PDF

Info

Publication number
CN106249295A
CN106249295A CN201510330765.6A CN201510330765A CN106249295A CN 106249295 A CN106249295 A CN 106249295A CN 201510330765 A CN201510330765 A CN 201510330765A CN 106249295 A CN106249295 A CN 106249295A
Authority
CN
China
Prior art keywords
wave
depth
vertical
component
microseismic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510330765.6A
Other languages
Chinese (zh)
Other versions
CN106249295B (en
Inventor
余波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201510330765.6A priority Critical patent/CN106249295B/en
Publication of CN106249295A publication Critical patent/CN106249295A/en
Application granted granted Critical
Publication of CN106249295B publication Critical patent/CN106249295B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

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 θ1、θ2With a cymoscope vertical azimuthal angle beta of S ripple2
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 cymoscopes1、θ2With a cymoscope vertical azimuthal angle beta of S ripple2
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 (xi,yi) make polarization analysis;Definition instantaneous energy EiWith transient bearing φi:
E i = x i 2 + y i 2 - - - ( 1 )
tgφi=yi/xi (2)
First, according to above formula, calculate each sampling point transient bearing φiWith instantaneous energy Ei, make instantaneous energy E simultaneouslyiTo transient bearing φiRectangular histogram;Then, (x is drawni,yi) 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 Ei、φ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:
Ri=xi cos(φ)+yi sin(φ)
(3)
Ti=-xi sin(φ)+yi cos(φ)
Wherein, Ri、TiIt is respectively horizontal radial, tangential component instantaneous amplitude after rotating;
Equally, to cymoscope Z component, horizontal radial R component amplitude (Zi,Ri) make polarization analysis, draw (Zi,Ri) 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 θ1、θ2, meanwhile, one of them cymoscope S ripple is carried out resolution of vectors, it is thus achieved that vertical azimuthal angle beta2
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 LP1=fP1(H)、LP2=fP2(H) with S wave function LS2=fS2(H);
Specifically, being implemented as follows of functional relationship is set up:
Assume that between two cymoscopes, vertical dimension is h0, 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 [Vp1,Vp2,......,VpN-1,VpN]、[Vs1,Vs2,......,VsN-1,VsN], and depth H is also divided into corresponding N shell [H1,H2,......,HN-1,HN];
According to snell law, by two vertical azimuth angle theta of P ripple1、θ2With a vertical azimuthal angle beta of S ripple2, obtain the horizontal range relevant with depth H:
L P 1 = f P 1 ( H ) = ( H 1 + h 0 ) · t a n ( θ 1 ) + Σ i = 2 N H i · Vp i / Vp 1 · s i n ( θ 1 ) 1 - ( Vp i / Vp 1 · sin ( θ 1 ) ) 2 - - - ( 4 )
L P 2 = f P 2 ( H ) = Σ i = 1 N H i · Vp i / Vp 1 · s i n ( θ 2 ) 1 - ( Vp i / Vp 1 · s i n ( θ 2 ) ) 2 - - - ( 5 )
L S 2 = f S 2 ( H ) = Σ i = 1 N H i · Vs i / Vs 1 · s i n ( β 2 ) 1 - ( Vs i / Vs 1 · s i n ( β 2 ) ) 2 - - - ( 6 )
Wherein, LP1、LP2It is respectively the P-wave distance cymoscope radial level distance that shallow-layer, any depth H of deep layer three-component geophone are corresponding, LS2For 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 | | fP1(H)-fP2(H)||+||fP1(H)-fS2(H) | |, when OPJ → 0, quickly calculate optimal solution depth H0With respective horizontal distance L0=(fP1(H0)+fP2(H0)+fS2(H0))/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 H0With respective horizontal distance L0=(fP1(H0)+fP2(H0)+fS2(H0))/3, it is achieved P, S ripple co-located, object function is as follows:
OPJ (H)=(LP1-LP2)+(LP1-LS2)
(7)
=(fP1(H)-fP2(H))+(fP1(H)-fS2(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 calculated0With respective horizontal distance L0Use Newton's dichotomy, be implemented as follows:
Rapidly find out two depth H1、H2So that OPJ (H1)·OPJ(H2)≤0;
Choose intermediate value H3=(H1+H2)/2, compare OPJ (H1)、OPJ(H2)、OPJ(H3), if OPJ is (H1)·OPJ(H3)≤0, then make new degree of depth combination H1=H1、H2=H3If, OPJ (H2)·OPJ(H3)≤0, then make new degree of depth combination H1=H2、H2=H3
Recalculate intermediate value object function, then compare, such iteration, until finally meeting convergence error, obtain optimal solution H0=(H1+H2)/2, substitute into object function, calculate optimal level distance L0=(fP1(H0)+fP2(H0)+fS2(H0))/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.一种井中微地震P、S波联合快速定位方法,其特征在于,包括如下步骤:1. a microseismic P, S wave joint fast location method in a well, is characterized in that, comprises the steps: a、获取同一事件微地震数据,对同一事件微地震数据进行矢量分解,获得两个三分量检波器的垂向方位角;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、根据snell定律,建立三分量检波器为起点的垂向深度H与水平距离L之间的函数关系;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、建立与震源点位置有关P&S波深度目标函数,通过扫描深度,使得目标函数最小,计算出震源点相对三分量检波器垂直深度与水平距离,实现P、S波联合快速定位。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.根据权利要求1所述一种井中微地震P、S波联合快速定位方法,其特征在于,所述步骤a的具体实现为: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: 获取同一事件微地震数据,对同一事件微地震数据进行矢量分解,计算出两个三分量检波器P波垂直方位角θ1、θ2和一个检波器S波垂向方位角β2The microseismic data of the same event are obtained, and the microseismic data of the same event are vector decomposed to calculate the vertical azimuths θ 1 and θ 2 of two three-component geophones for P waves and the vertical azimuth β 2 of one geophone for S waves. 3.根据权利要求2所述一种井中微地震P、S波联合快速定位方法,其特征在于,步骤b的具体实现为: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: 由已知地层介质速度模型,根据snell定律,从三分量检波器出发,沿波传播路径逆方向,建立以三分量检波器为起点的垂向深度H与水平距离L之间的函数关系,对应有两个P波函数LP1=fP1(H)、LP2=fP2(H)与一个S波函数LS2=fS2(H)。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.根据权利要求3所述一种井中微地震P、S波联合快速定位方法,其特征在于,所述步骤c的具体实现为: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: 设置与震源位置有关的目标函数OPJ=‖fP1(H)-fP2(H)‖+‖fP1(H)-fS2(H)‖,当OPJ→0时,快速计算出最优解深度H0与相应水平距离L0=(fP1(H0)+fP2(H0)+fS2(H0))/3,即所求为震源点相对三分量检波器垂直深度与水平距离。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 0 and corresponding horizontal distance L 0 = (f P1 (H 0 )+f P2 (H 0 )+f S2 (H 0 ) )/3, that is, the vertical depth and horizontal distance. 5.根据权利要求2所述一种井中微地震P、S波联合快速定位方法,其特征在于,所述矢量分解采用矢端曲线-直方图法实现,具体如下: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: 设置一个时窗,包含了P波事件,对三分量检波器X、Y分量振幅值(xi,yi)作偏振分析;定义瞬时能量Ei和瞬时方位φiSet 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 : EE. ii == xx ii 22 ++ ythe y ii 22 -- -- -- (( 11 )) tgφi=yi/xi (2)tgφ i =y i /x i (2) 首先,根据上式,计算出每个样点瞬时方位φi与瞬时能量Ei,同时作出瞬时能量Ei对瞬时方位φi的直方图;然后,画出(xi,yi)振幅值坐标系连线图,即矢端曲线图,对其进行线性拟合,估算出倾角范围;最后,参考倾角范围,在Ei、φi直方图中寻找瞬时能量最大值的位置,此时与其对应的角度φ,即为三分量检波器相对波传播方向水平方位角;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; 利用三分量检波器相对波传播方向水平方位角φ,对X、Y分量进行旋转处理,获得水平径向R与切向T: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: Ri=xicos(φ)+yisin(φ) (3)R i =x i cos(φ)+y i sin(φ) (3) Ti=-xi sin(φ)+yi cos(φ)T i =-x i sin(φ)+y i cos(φ) 其中,Ri、Ti分别为旋转后水平径向、切向分量瞬时振幅;Among them, R i and T i are the instantaneous amplitudes of horizontal radial and tangential components respectively after rotation; 同样,对检波器Z分量、水平径向R分量振幅值(Zi,Ri)作偏振分析,画出(Zi,Ri)振幅值连线图,即矢端曲线图,对其进行线性拟合,估算出倾角范围;参考倾角范围,在直方图中寻找瞬时能量最大值的位置,求出其对应的检波器相对波传播方向垂向方位角θ,即完成矢量分解与计算垂向方位角;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; 对观测井任意两个三分量检波器P波进行矢量分解,获得相应垂向方位角为θ1、θ2,同时,对其中一个检波器S波进行矢量分解,获得垂向方位角β2Perform vector decomposition on the P waves of any two three-component geophones in the observation well to obtain the corresponding vertical azimuths θ 1 and θ 2 , and at the same time, perform vector decomposition on the S wave of one of the geophones to obtain the vertical azimuth β 2 . 6.根据权利要求3所述一种井中微地震P、S波联合快速定位方法,其特征在于,建立函数关系的具体实现如下: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: 假设两个检波器之间垂直距离为h0,给出距离第二个检波器地下任意深度H,根据声波测井,分成从浅到深N层纵横波速度层[Vp1,Vp2,......,VpN-1,VpN]、[Vs1,Vs2,......,VsN-1,VsN],且深度H亦分为相应N层[H1,H2,......,HN-1,HN];Assuming that the vertical distance between the two geophones is h 0 , 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 1 , Vp 2 ,. .....,Vp N-1 ,Vp N ], [Vs 1 ,Vs 2 ,......,Vs N-1 ,Vs N ], and the depth H is also divided into corresponding N layers [H 1 ,H 2 ,...,H N-1 ,H N ]; 根据snell定律,由两个P波垂向方位角θ1、θ2与一个S波垂向方位角β2,得到与深度H有关的水平距离:According to Snell's law, the horizontal distance related to the depth H is obtained from two P-wave vertical azimuths θ 1 and θ 2 and one S-wave vertical azimuth β 2 : LL PP 11 == ff PP 11 (( Hh )) == (( Hh 11 ++ hh 00 )) ·· tanthe tan (( θθ 11 )) ++ ΣΣ ii == 22 NN Hh ii ·· VpVp ii // VpVp 11 ·· sinsin (( θθ 11 )) 11 -- (( VpVp ii // VpVp 11 ·· sinsin (( θθ 11 )) )) 22 -- -- -- (( 44 )) LL PP 22 == ff PP 22 (( Hh )) == ΣΣ ii == 11 NN Hh ii ·&Center Dot; VpVp ii // VpVp 11 ·· sthe s ii nno (( θθ 22 )) 11 -- (( VpVp ii // VpVp 11 ·&Center Dot; sthe s ii nno (( θθ 22 )) )) 22 -- -- -- (( 55 )) LL SS 22 == ff SS 22 (( Hh )) == ΣΣ ii == 11 NN Hh ii ·&Center Dot; Vsvs. ii ·&Center Dot; // Vsvs. 11 ·&Center Dot; sthe s ii nno (( ββ 22 )) 11 -- (( Vsvs. ii ·&Center Dot; // Vsvs. 11 ·&Center Dot; sthe s ii nno (( ββ 22 )) )) 22 -- -- -- (( 66 )) 其中,LP1、LP2分别为浅层、深层三分量检波器任意深度H对应的P波传播距离检波器径向水平距离,LS2为深层检波器任意深度H对应的S波传播距离检波器径向水平距离。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.根据权利要求4所述一种井中微地震P、S波联合快速定位方法,其特征在于,计算震源点相对三分量检波器垂直深度与水平距离的具体实现为: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: 建立与震源位置有关的目标函数,快速计算出最优解深度H0与相应水平距离L0=(fP1(H0)+fP2(H0)+fS2(H0))/3,实现P、S波联合定位,目标函数如下:Establish the objective function related to the source position, quickly calculate the optimal solution depth H 0 and the corresponding horizontal distance L 0 =(f P1 (H 0 )+f P2 (H 0 )+f S2 (H 0 ))/3, To achieve joint positioning of P and S waves, the objective function is as follows: OPJ(H)=(LP1-LP2)+(LP1-LS2) (7)OPJ(H)=(L P1 -L P2 )+(L P1 -L S2 ) (7) =(fP1(H)-fP2(H))+(fP1(H)-fS2(H))=(f P1 (H)-f P2 (H))+(f P1 (H)-f S2 (H)) 当OPJ→0时,表示一个检波器P波水平距离同时与另一个检波器P波、S波水平距离相交,该交点位置即为震源位置。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.根据权利要求7所述一种井中微地震P、S波联合快速定位方法,其特征在于,快速计算出最优解深度H0与相应水平距离L0采用牛顿二分法,具体实现如下: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 0 and corresponding horizontal distance L 0 adopt Newton's dichotomy method, concrete realization is as follows: 快速找出两个深度H1、H2,使得OPJ(H1)·OPJ(H2)≤0;Quickly find two depths H 1 and H 2 so that OPJ(H 1 )·OPJ(H 2 )≤0; 选取中间值H3=(H1+H2)/2,比较OPJ(H1)、OPJ(H2)、OPJ(H3),如果OPJ(H1)·OPJ(H3)≤0,则令新的深度组合H1=H1、H2=H3,如果OPJ(H2)·OPJ(H3)≤0,则令新的深度组合H1=H2、H2=H3Select the middle value H 3 =(H 1 +H 2 )/2, compare OPJ(H 1 ), OPJ(H 2 ), OPJ(H 3 ), if OPJ(H 1 )·OPJ(H 3 )≤0, Then set the new depth combination H 1 =H 1 , H 2 =H 3 , if OPJ(H 2 )·OPJ(H 3 )≤0, then set the new depth combination H 1 =H 2 , H 2 =H 3 ; 重新计算中间值目标函数,再进行比较,如此迭代,直至最后符合收敛误差,得到最优解H0=(H1+H2)/2,代入目标函数,计算出最优水平距离L0=(fP1(H0)+fP2(H0)+fS2(H0))/3,实现P、S波联合快速定位。Recalculate the intermediate value objective function, and then compare, and iterate until the convergence error is finally met, and the optimal solution H 0 =(H 1 +H 2 )/2 is obtained, which is substituted into the objective function to calculate the optimal horizontal distance L 0 = (f P1 (H 0 )+f P2 (H 0 )+f S2 (H 0 ))/3, realizing joint fast positioning of P and S waves. 9.一种井中微地震P、S波联合快速定位系统,其特征在于,包括数据采集模块、矢量分解模块、函数建立模块和结果计算模块;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; 所述函数建立模块,其用于根据snell定律,建立三分量检波器为起点的垂向深度H与水平距离L之间的函数关系;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; 所述结果计算模块,其用于建立与震源点位置有关P&S波深度目标函数,通过扫描深度,使得目标函数最小,计算出震源点相对三分量检波器垂直深度与水平距离,实现P、S波联合快速定位。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.
CN201510330765.6A 2015-06-15 2015-06-15 A kind of borehole microseismic P, S wave joint method for rapidly positioning and system Active CN106249295B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510330765.6A CN106249295B (en) 2015-06-15 2015-06-15 A kind of borehole microseismic P, S wave joint method for rapidly positioning and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510330765.6A CN106249295B (en) 2015-06-15 2015-06-15 A kind of borehole microseismic P, S wave joint method for rapidly positioning and system

Publications (2)

Publication Number Publication Date
CN106249295A true CN106249295A (en) 2016-12-21
CN106249295B CN106249295B (en) 2018-08-31

Family

ID=57626379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510330765.6A Active CN106249295B (en) 2015-06-15 2015-06-15 A kind of borehole microseismic P, S wave joint method for rapidly positioning and system

Country Status (1)

Country Link
CN (1) CN106249295B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772591A (en) * 2017-04-05 2017-05-31 吉林大学 A kind of combined positioning-method suitable for improving microseism reliability of positioning
CN109655901A (en) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 A kind of adaptive angle of polarization calculation method of frequency domain and system
CN110146924A (en) * 2019-07-03 2019-08-20 中国地质大学(北京) Position and azimuth inversion method of submarine seismograph based on first arrival polarization azimuth of water wave
CN110716230A (en) * 2018-07-13 2020-01-21 中国石油化工股份有限公司 Well-ground combined micro-seismic positioning method
CN112068204A (en) * 2019-06-10 2020-12-11 中国石油化工股份有限公司 Remote borehole microseism monitoring and positioning method and computer storage medium
CN112068205A (en) * 2019-06-10 2020-12-11 中国石油化工股份有限公司 Microseism event rapid positioning method for full-coverage well-ground combined monitoring
CN114170750A (en) * 2021-12-15 2022-03-11 姚珍栋 An intrusion alarm system and method based on P-S wave detection
CN114428353A (en) * 2020-10-10 2022-05-03 中国石油化工股份有限公司 Microseismic compression and shear wave positioning method, device, electronic equipment and medium for inclined wells
CN114428357A (en) * 2020-10-10 2022-05-03 中国石油化工股份有限公司 Vertical well microseism longitudinal and transverse wave positioning method and device, electronic equipment and medium
CN116027395A (en) * 2021-10-25 2023-04-28 中国石油化工股份有限公司 Multi-component seismic data R-T rotation method and device, electronic equipment and medium
CN117412279A (en) * 2023-12-14 2024-01-16 江西北斗云智慧科技有限公司 Beidou emergency call method, system, computer and storage medium

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106772591B (en) * 2017-04-05 2018-08-14 吉林大学 A kind of joint positioning method being suitable for improving microseism reliability of positioning
CN106772591A (en) * 2017-04-05 2017-05-31 吉林大学 A kind of combined positioning-method suitable for improving microseism reliability of positioning
CN109655901B (en) * 2017-10-11 2020-08-25 中国石油化工股份有限公司 Frequency domain self-adaptive polarization angle calculation method and system
CN109655901A (en) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 A kind of adaptive angle of polarization calculation method of frequency domain and system
CN110716230B (en) * 2018-07-13 2021-08-24 中国石油化工股份有限公司 Well-ground combined micro-seismic positioning method
CN110716230A (en) * 2018-07-13 2020-01-21 中国石油化工股份有限公司 Well-ground combined micro-seismic positioning method
CN112068205A (en) * 2019-06-10 2020-12-11 中国石油化工股份有限公司 Microseism event rapid positioning method for full-coverage well-ground combined monitoring
CN112068204A (en) * 2019-06-10 2020-12-11 中国石油化工股份有限公司 Remote borehole microseism monitoring and positioning method and computer storage medium
CN110146924B (en) * 2019-07-03 2020-05-26 中国地质大学(北京) Submarine seismograph position and orientation inversion method based on water wave first arrival polarization orientation
CN110146924A (en) * 2019-07-03 2019-08-20 中国地质大学(北京) Position and azimuth inversion method of submarine seismograph based on first arrival polarization azimuth of water wave
CN114428353A (en) * 2020-10-10 2022-05-03 中国石油化工股份有限公司 Microseismic compression and shear wave positioning method, device, electronic equipment and medium for inclined wells
CN114428357A (en) * 2020-10-10 2022-05-03 中国石油化工股份有限公司 Vertical well microseism longitudinal and transverse wave positioning method and device, electronic equipment and medium
CN114428357B (en) * 2020-10-10 2025-04-18 中国石油化工股份有限公司 Vertical well microseismic longitudinal and transverse wave positioning method, device, electronic equipment and medium
CN116027395A (en) * 2021-10-25 2023-04-28 中国石油化工股份有限公司 Multi-component seismic data R-T rotation method and device, electronic equipment and medium
CN114170750A (en) * 2021-12-15 2022-03-11 姚珍栋 An intrusion alarm system and method based on P-S wave detection
CN117412279A (en) * 2023-12-14 2024-01-16 江西北斗云智慧科技有限公司 Beidou emergency call method, system, computer and storage medium

Also Published As

Publication number Publication date
CN106249295B (en) 2018-08-31

Similar Documents

Publication Publication Date Title
CN106249295A (en) A combined rapid positioning method and system for microseismic P and S waves in a borehole
CN100549731C (en) Secondary positioning method of detection point
CN106772577A (en) Source inversion method based on microseism data and SPSA optimized algorithms
CN107290722B (en) The localization method and device of microquake sources
CN106094029A (en) The method utilizing offset distance vector sheet geological data Predicating Reservoir Fractures
CN105510880A (en) Microseism focus positioning method based on double-difference method
CN104730574B (en) The method for building near surface structure model
CN105607119B (en) Near-surface model construction method and static correction value acquiring method
CN111257941B (en) A combined submarine seismograph azimuth automatic identification device and method
CN103149588B (en) Method and system for calculating VTI anisotropic parameters by using well seismic calibration
CN104166159B (en) The fracture morphology disposal route of four-dimensional micro-seismic monitoring and system
CN106468781B (en) A kind of crack prediction method based on minimum entropy rotary process
CN106468782A (en) A kind of based on ceiling capacity than method crack prediction method
CN105738952A (en) Horizontal well region reservoir rock facies modeling method
Poliannikov et al. Interferometric hydrofracture microseism localization using neighboring fracture
CN103954996B (en) Method and device for determining fracture strike of stratum fracture based on travel time method
CN109212594B (en) Combined positioning method for longitudinal waves and transverse waves of anisotropic medium
Zhao et al. Borehole azimuthal acoustic imaging using 3D spatial scanning: Application in acoustic detection of nearby wells
CN107765306B (en) A kind of VSP initial velocity modeling method and device
CN104297790A (en) Quantitative evaluation method for offset distribution in fluctuating submarine earthquake observation system
CN103558637A (en) Far detection method based on three-component sensor
CN107340537A (en) A kind of method of P-SV converted waves prestack reverse-time depth migration
CN102914797A (en) Method and device for obtaining formation anisotropy coefficient
CN109085642B (en) Anisotropic medium microseism event positioning method
CN113960532A (en) Microseism positioning method based on secondary positioning calculation of imaginary source

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant