CN104280772A - Recognition method for microseism phase in well - Google Patents
Recognition method for microseism phase in well Download PDFInfo
- Publication number
- CN104280772A CN104280772A CN201310292019.3A CN201310292019A CN104280772A CN 104280772 A CN104280772 A CN 104280772A CN 201310292019 A CN201310292019 A CN 201310292019A CN 104280772 A CN104280772 A CN 104280772A
- Authority
- CN
- China
- Prior art keywords
- component
- perforation
- data
- wave
- ripple
- 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.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a recognition method for a microseism phase in a well and belongs to the field of three-component signal processing of a microseism in the well. The recognition method includes the steps that (1), according to known perforation three-component data of the microseism in the well, rotation treatment is performed on a component X and a component Y of the perforation data to obtain a rotated new component Rx, a rotated new component Ry and a perforation orientation angle theta, relative to the perforation position of the fractured well, of a component X of a monitoring well three-component detector; (2), the perforation orientation angle theta is used for performing orientation uniformity correction treatment on a component X and a component Y of fracturing data of the microseism in the well to obtain a fracturing new component Cx and a fracturing new component Cy uniform in polarity; (3), Cx is selected as a reference channel K; (4), the microseism phase type is recognized through judging positive and negative polarities of amplitudes of the reference trace K and a component Z at the same time, then a P wave and an S wave in an effective event of the component Z are separated, and finally P wave type data and S wave type data are output.
Description
Technical field
The invention belongs to borehole microseismic three component signal process field, be specifically related to a kind of borehole microseismic phase identification method.
Background technology
The seismic phase of seismic signal is separated along with the development of seismic prospecting polarization technology, multi-components wave detector widespread use and progressively developing, and is carry out seismic phase identification with relation between seismic event spatial polarization characteristic and direction of wave travel.Lot of experiments and actual data application analysis are all carried out both at home and abroad in this respect.
Abroad to earthquake signal polarization Quality Research comparatively early, famous seismologist B.B. Ge Licen first Application three-component seismic data in 1909, determines the position angle of direction of wave travel horizontal projection according to the particle movement direction of first component direct P ripple.Nineteen fifty-two, adopt husband by Γ .A. Clarence Gamble and apply wave polarization feature, propose orientation--Phase contrast, utilize observation station vibration parameters to follow the trail of the phase place of seismic event.B. И in 1964. Bamda row husband has further developed seismic location observation method, proposes particle displacement direction calculating, elliptic polarization vibration parameters analytical approach.Afterwards between 1969 to 1979, the people such as И .B. Bo Meilancai, А .H. Mo Zirenke utilize multi-component exploration, identify P ripple, PS ripple, the information of penetrating the corresponding polarization of shear wave SH ripple and other kinetic characters.Phase late 1980s, earthquake wave polarization is used to the effective sensitivity increasing instrument, works out polarized filtering method, carries out filtering according to significant wave, interference wave by polarization characteristic; Since the nineties in 20th century, due to crosshole seismic, VSP exploration and the widespread use of multi-components wave detector, complex wave field more can be obtained, developed and utilize polarization--position versus method realizes seismic phase and is effectively separated.
The domestic research be separated for seismic phase is started late, mainly to the earth earthquake, petroleum prospecting multi-component data carries out pie slice, Radon transform, polarization filterings etc. realize wave field in length and breadth and are separated, as Liu Xiqiang etc. (1998) proposes the three-component seismologic record phase identification method proposing based on single seismographic station based on weak shock phase identification method and (2000) of wavelet package transforms, the three-component seismic phase that Li Huiting etc. (2000) propose identifies automatically, the application of Wang Juan etc. (2004) neural network in seismic phase identification, Li Shanyou etc. (2006) propose the automatic recognition of seismic phase based on amplitude and instantaneous frequency, the near earthquake phase identification method based on " the acceleration array " that Liu's equality (2011) proposes, Song Weiqi etc. (2008) microearthquake data frequency field is relevant-time domain the polarization filtering method identification of single seismic phase micro-seismic event and inverting etc. that propose with (2011).
For microseismic, it be a kind of by observation, analyze in activity in production the small seismic events produced and monitor the impact of activity in production, effect and underground state, in the oil-gas field development stage, the fracturing yield increasing of reservoir can be realized, optimize well spacing, improve recovery ratio.In addition, in hydraulically created fracture process, the micro-seismic event seismic phase type produced is random, P ripple and S ripple may be produced simultaneously, also single P ripple or S ripple may only be produced, be subject to formation factor, the impact of well surrounding enviroment and the interference of all kinds of noise, the microearthquake data that physical record is arrived is more complicated simultaneously.Wherein, seismic phase identification is one of based process of micro-seismic event inverting, how to identify from the microseismograms of complexity and isolate effective micro-seismic event, determine the seismic phase type of these effective micro-seismic event simultaneously, most important to the work of microearthquake inverting from now on.
Summary of the invention
The object of the invention is to solve the difficult problem existed in above-mentioned prior art, a kind of borehole microseismic phase identification method is provided, library track is obtained according to perforation data orientation consistency treatment, carry out simple polarity judgement again, just can realize automatically identifying separating P wave from S wave seismic phase categorical data from the multiple event of pressure break.
The present invention is achieved through the following technical solutions borehole microseismic P ripple, the identification of S ripple, and sub-argument goes out corresponding data:
A kind of borehole microseismic phase identification method, comprises the following steps:
(1), by known borehole microseismic perforation three component seismic data, correlation hole data X component and perforation data Y-component carry out the perforation orientation angle θ rotating the rear new components R x of process acquisition rotation and rotate rear newly components R y and the relative fractured well perforating site of monitor well three-component seismometer X component;
(2), utilize the perforation orientation angle θ that step (1) obtains, borehole microseismic pressure break data X, Y two component are carried out to orientation Concordance process and obtain polarity consistent pressure break new component Cx and the new component Cy of pressure break;
(3), from the C that step (2) obtains
x, C
yin, choose C
xas reference road K.
(4), by judging that the positive-negative polarity of library track K and the same time amplitude of borehole microseismic Z component identifies seismic phase type, then isolating the P ripple in Z component validity event and S ripple, finally exporting P ripple categorical data and S ripple categorical data.
The perforation orientation angle θ of the relative fractured well perforating site of monitor well three-component seismometer X component in described step (1) adopts horizontal component hodograph-histogram to ask for out.
Described step (2) is decomposed by coordinate geometry and is realized, specific as follows:
Input azimuth angle theta, re-starts vertical coordinate system geometric angle decompose pressure break data X, Y two component:
C
X＝x
i cos(θ)+y
i sin(θ)
C
Y＝-x
i sin(θ)+y
i cos(θ)
Wherein, x
i, y
ibe respectively X, Y-component amplitude, C
x, C
yfor two component amplitude values new after geometry value decomposition, and C
xdirection is that perforation is to monitor well wave detector horizontal projection direction.
Described by judging that the positive-negative polarity of library track K and borehole microseismic Z component same time amplitude identifies that seismic phase type is specific as follows in described step (4):
If current wave field is upward traveling wave, as the amplitude A of library track K
kwith the amplitude A of borehole microseismic Z component
zbe multiplied when being greater than 0, i.e. A
ka
zduring > 0, positive-negative polarity performance is consistent, and now seismic phase is judged as P ripple, works as A
ka
zduring < 0, positive-negative polarity performance is contrary, and now seismic phase is judged as S ripple;
If current wave field is down going wave, judge that conclusion is then contrary with upward traveling wave.
Compared with prior art, the invention has the beneficial effects as follows: the present invention can identify borehole microseismic pressure break event seismic phase type exactly, and its operating process is simple and practical, accurately locating for follow-up specific aim P ripple or S ripple source location provides authentic data basis.
Accompanying drawing explanation
Fig. 1 is that the present invention realizes borehole microseismic seismic phase identifying operation process flow diagram;
Fig. 2 sets up borehole microseismic model views system schematic;
Fig. 3-1 is the perforation X component in borehole microseismic model fractured well water filling perforation three-component microseism section;
Fig. 3-2 is the perforation Y-component in borehole microseismic model fractured well water filling perforation three-component microseism section;
Fig. 3-3 is the perforation Z components in borehole microseismic model fractured well water filling perforation three-component microseism section;
Fig. 4-1 is that borehole microseismic model produces the X component in the three-component microseism section of pressure break event one correspondence by perforation;
Fig. 4-2 is that borehole microseismic model produces the Y-component in the three-component microseism section of pressure break event one correspondence by perforation;
Fig. 4-3 is that borehole microseismic model produces the Z component in the three-component microseism section of pressure break event one correspondence by perforation;
Fig. 5-1 is that borehole microseismic model produces the X component in the three-component microseism section of pressure break event two correspondence by perforation;
Fig. 5-2 is that borehole microseismic model produces the Y-component in the three-component microseism section of pressure break event two correspondence by perforation;
Fig. 5-3 is that borehole microseismic model produces the Z component in the three-component microseism section of pressure break event two correspondence by perforation;
Fig. 6-1 is the new components R x that in borehole microseismic model perforation data, level two component obtains after rotating process;
Fig. 6-2 is new components R y that in borehole microseismic model perforation data, level two component obtains after rotating process;
Fig. 7 is borehole microseismic model perforation data: the relative geophone orientation angle of perforation obtained after rotating;
Fig. 8-1 is borehole microseismic model pressure break event one: new component Cx after perforation orientation Concordance;
Fig. 8-2 is model pressure break events one: new component Cy after perforation orientation Concordance;
Fig. 9-1 is borehole microseismic model pressure break event one: identifying seismic phase in Z component is P wave datum;
Fig. 9-2 is borehole microseismic model pressure break events one: identifying seismic phase in Z component is S wave datum;
Figure 10-1 is borehole microseismic model pressure break event two: the new component Cx obtained after perforation orientation Concordance;
Figure 10-2 is borehole microseismic model pressure break events two: the new component Cy obtained after perforation orientation Concordance;
Figure 11-1 is borehole microseismic model pressure break event two: identifying seismic phase in Z component is P wave datum;
Figure 11-2 is borehole microseismic model pressure break events two: identifying seismic phase in Z component is S wave datum;
Figure 12-1 is borehole microseismic real data: perforation X component;
Figure 12-2 is borehole microseismic real data: perforation Y-component;
Figure 12-3 is borehole microseismic real data: perforation Z component;
Figure 13-1 is borehole microseismic real data: the actual pressure break X component of pressure break three-component microseism section produced by perforation;
Figure 13-2 is borehole microseismic real data: the actual pressure break Y-component of pressure break three-component microseism section produced by perforation;
Figure 13-3 is borehole microseismic real data: the actual pressure break Z component of pressure break three-component microseism section produced by perforation;
Figure 14-1 is borehole microseismic real data: the new components R x obtained after perforation X, Y-component rotate;
Figure 14-2 is borehole microseismic real data: the new components R y obtained after perforation X, Y-component rotate;
Figure 15 is borehole microseismic real data: the relative geophone orientation angle of perforation;
Figure 16-1 is borehole microseismic real data: obtain new component Cx after pressure break data perforation orientation Concordance;
Figure 16-2 is borehole microseismic real data: obtain new component Cy after pressure break data perforation orientation Concordance;
Figure 17-1 is borehole microseismic real data: identify pressure break data Z component microearthquake seismic phase P wave datum;
Figure 17-2 is borehole microseismic real data: identify pressure break data Z component microearthquake seismic phase S wave datum;
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail:
The invention provides one and judge borehole microseismic phase identification method based on perforation orientation angle consistance library track polarity, as shown in Figure 1, realize the method and comprise following four steps:
(1), by known borehole microseismic perforation three component seismic data, level two component (i.e. perforation data X component and perforation data Y-component) is carried out rotating process acquisition and rotates rear new components R x new components R y and monitor well three-component seismometer X component relative fractured well perforating site relative bearing θ (perforation orientation angle namely in Fig. 1) rear with rotation;
(2), utilize the perforation orientation angle θ that step (1) obtains, orientation Concordance process is carried out to borehole microseismic pressure break data X, Y two component, make wave detector X, Y-component direction points to the surface level projecting direction of detector position from unordered becoming at random by perforating site, and obtain new level two component that polarity is consistent, i.e. pressure break new component Cx and the new component Cy of pressure break;
(3), from the C obtained
x, C
yin, choose C
xas reference road K.
(4), by judging that the positive-negative polarity of library track K and the same time amplitude of borehole microseismic Z component identifies seismic phase type, after decomposing according to P ripple, S ripple vertical coordinate system, polarity rule judges that current event seismic phase type is P ripple or S ripple, finally realize borehole microseismic all event seismic phase type identification, then export Z component microearthquake shake to shake as S ripple categorical data, so that follow-up borehole microseismic P ripple or the process of S ripple specific aim as P ripple categorical data and Z component microearthquake.
Described by judging that the positive-negative polarity of library track K and borehole microseismic Z component same time amplitude identifies that seismic phase type is specific as follows in described step (4):
If current wave field is upward traveling wave, when two component amplitude values are multiplied A
ka
zduring > 0, positive-negative polarity performance is consistent, and now seismic phase is judged as compressional wave, works as A
ka
zduring < 0, positive-negative polarity performance is contrary, and now seismic phase is judged as shear wave; And if current wave field is down going wave, correspondence judges that conclusion is then just contrary with upward traveling wave.
Monitor well three-component seismometer X component relative fractured well perforating site relative bearing θ in described step (1) adopts that existing algorithm---horizontal component hodograph-histogram is asked for out, specific as follows:
If two component X, Y i-th sampling point amplitudes are (x
i, y
i), then this instantaneous energy E
iwith transient bearing θ
ibe expressed as:
tgθ
i＝y
i/x
i
First, according to above formula, calculate each sampling point transient bearing θ
iwith instantaneous energy E
i, make instantaneous energy E simultaneously
ito transient bearing θ
ihistogram; Then, (x is drawn
i, y
i) amplitude coordinate system line graph (also known as hodograph), according to hodograph linear relationship, estimation inclination angle scope; Finally, with reference to inclination angle scope, at E
i, θ
ifind the position of instantaneous energy peak-peak in histogram, now corresponding with it angle, be required azimuth angle theta.
Described step (2) is decomposed by coordinate geometry and is realized, specific as follows:
Input azimuth angle theta, re-starts vertical coordinate system geometric angle decompose pressure break data X, Y two component:
C
X＝x
i cos(θ)+y
i sin(θ)
C
Y＝-x
i sin(θ)+y
i cos(θ)
Wherein, x
i, y
ibe respectively X, Y-component amplitude, C
x, C
yfor two component amplitude values new after geometry value decomposition, and C
xdirection is that perforation is to monitor well wave detector horizontal projection direction.
Key of the present invention judges that library track K and borehole microseismic Z component same time amplitude positive-negative polarity are to identify seismic phase type, and after decomposing using P ripple, S ripple vertical coordinate system, polarity rule is as seismic phase type basis for estimation:
First, set up vertical coordinate system, be made up of vertical downward direction M and horizontal direction N.
Then, suppose that direction of wave travel and the coordinate M angle of cut are
wave amplitude value size is A, if seismic phase is P ripple, because direction of vibration is consistent with the direction of propagation, then seismic event decomposes M-N coordinate components and is:
If seismic phase is S ripple, because direction of vibration is vertical with the direction of propagation, then now seismic event decomposes M-N coordinate components and is:
Finally, according to recording geometry in real well, if the microearthquake wave field that wave detector receives is upward traveling wave, then the respective coordinates angle of cut
scope is 0 ~ 90 degree, now two component amplitude P after P Wave Decomposition
m, P
npositive-negative polarity is identical, and two component amplitude S after S Wave Decomposition
m, S
npositive-negative polarity is contrary; If the microearthquake wave field received is down going wave, then the respective coordinates angle of cut
scope is 90 ~ 180 degree, now P
mwith P
npositive-negative polarity is contrary, and S
mwith S
npositive-negative polarity is identical.
From contrasting above, P ripple, S ripple component polarity after vertical coordinate system M-N decomposes exists regular, we can using this rule as blind source microearthquake seismic phase basis for estimation, if that is: microearthquake wave field is upward traveling wave, and vertical coordinate two component M, N amplitude A
ma
n> 0, positive-negative polarity is identical, judges that this seismic phase is P ripple type, works as A
ma
nduring < 0, positive-negative polarity is contrary, judges that this seismic phase is S ripple; If microearthquake wave field is down going wave, work as A
ma
nduring > 0, judge that this seismic phase is S ripple, work as A
ma
nduring < 0, then now seismic phase is P ripple, and conclusion and upward traveling wave are just contrary.
The microearthquake three component signal that monitor well well seismometer receives, there is the features such as pressure break source location is unknown, focus seismic phase type is unknown, wave detector horizontal direction is pointed at random, and in succession there is multiple event in the even short period, the present invention is directed to borehole microseismic seismic phase identification difficult point, propose one and judge seismic phase method of identification based on perforation orientation angle consistance library track polarity, judge by reference to road and pressure break Z component amplitude positive-negative polarity, can automatically identify seismic phase type is P ripple or S ripple.
Lower mask body introduces the present invention.First, by known borehole microseismic perforation three component seismic data, utilize hodograph--histogram method, carries out rotation process to level two component, obtains perforation to wave detector direction and the wave detector X component direction angle of cut, i.e. perforation relative bearing θ.
If X, Y two component i-th sampling point amplitude is (x
i, y
i), then this instantaneous energy E
iwith transient bearing θ
ibe expressed as:
tgθ
i＝y
i/x
i (2)
According to formula (1), (2), calculate each sampling point transient bearing θ respectively
iwith instantaneous energy E
i, make all instantaneous energy E simultaneously
ito transient bearing θ
ihistogram; Meanwhile, each sampling point (x is utilized
i, y
i) amplitude, with x
ifor horizontal ordinate, y
ifor ordinate, make all sampling point line graphs (namely representing the hodograph in particle movement direction), according to component direct P ripple hodograph, there is linear relationship feature, estimation position angle value probable ranges; Within the scope of the valuation of position angle, from E
i, θ
ithe position of instantaneous energy peak-peak is found, now corresponding with it θ in histogram
iangle, is required perforation relative bearing θ.
Then, utilize azimuth angle theta, orientation Concordance process is carried out to borehole microseismic pressure break data level X, Y two component, make all wave detector X, Y-component direction becomes the relative geophone orientation of perforation from unordered at random, and according to coordinate angle decomposition, obtain the two new component C that polarity is consistent
x, C
y:
C
X＝x
i cos(θ)+y
i sin(θ) (3)
C
Y＝-x
i sin(θ)+y
i cos(θ)
Secondly, from two new component C
x, C
yin, choose C
xas reference component K.
Finally, according to borehole microseismic Z component TRAVEL TIME, with wave detector change in depth feature, (if whilst on tour diminishes with depth shallower, this wave field is upward traveling wave; If whilst on tour becomes large with depth shallower, this wave field is then down going wave), first determine that wave field is upward traveling wave or down going wave; Then, compare Z component and reference component K sampling point amplitude of same time positive-negative polarity, so that " for microearthquake upward traveling wave, if amplitude positive-negative polarity is identical between vertical coupled two components, judging that this seismic phase is P ripple type, if polarity is contrary, is then S ripple; For microearthquake down going wave, if amplitude positive-negative polarity is contrary between vertical coupled two components, judging that this seismic phase is P ripple type, if polarity is identical, is then S ripple " judgment criterion identifies now seismic phase type.
The present invention judges that seismic phase type is that two component polarity exist the inherent law of " identical " or " on the contrary " because P ripple, S ripple are after vertical coupled coordinate decomposes according to key:
First, set up vertical coordinate system M-N, wherein, direction M is vertically downward, and direction N is horizontal direction.Then, suppose that direction of wave travel and the coordinate M angle of cut are
wave amplitude value size is A, if seismic phase is P ripple, because direction of vibration is consistent with the direction of propagation, then seismic event decomposes M-N coordinate components and is:
If seismic phase is S ripple, because direction of vibration is vertical with the direction of propagation, then now seismic event decomposes M-N coordinate components and is:
Finally, because Z component in recording geometry in real well is always downward along well track, direction of wave travel and the coordinate M angle of cut are
span is divided into two kinds of situations: work as the angle of cut
scope is 0 ~ 90 degree, and the microearthquake wave field that now wave detector receives is upward traveling wave, judges, two component amplitude P after P Wave Decomposition according to formula (4), (5)
m, P
npolarity is identical, but two component amplitude S after S Wave Decomposition
mwith S
npolarity is contrary; If the microearthquake wave field received is down going wave, the angle of cut
scope is 90 ~ 180 degree, now two component P after P Wave Decomposition
mwith P
npositive-negative polarity is contrary, and two component S after S Wave Decomposition
mwith S
npositive-negative polarity is identical.
Because borehole microseismic wave field can be divided into single upward traveling wave, down going wave or both combinations, by perforation orientation Concordance, obtain the library track K mutually vertical coupled with microearthquake Z component, according to P ripple, the vertical coupled two component positive-negative polarity rules of S ripple, just can realize blind source microearthquake seismic phase type and accurately identify.
From borehole microseismic theoretical model data and actual data application, seismic phase recognition effect of the present invention is described respectively below.
Carry out borehole microseismic model data seismic phase discriminance analysis, the theoretical accuracy of main checking the inventive method.First, set up borehole microseismic monitoring system, as shown in Figure 2, design two pressure break points of a perforation and generation, and be all the two focus of P, S ripple microearthquake, monitor well places 16 grades of three-component seismometers and horizontal component direction has randomness.The borehole microseismic perforation three component seismic data whilst on tour section of Fig. 3-1 to Fig. 3-3 for simulating, Fig. 4-1 to Fig. 4-3, Fig. 5-1 are two pressure break three-component whilst on tour sections of corresponding generation to Fig. 5-3, difference is that Fig. 4-1 to Fig. 4-3 represents upward traveling wave wave field, and Fig. 5-1 to Fig. 5-3 represents down going wave wave field.
Then, operation steps according to Fig. 1, takes hodograph-histogram method to rotate process (as shown in Fig. 6-1 to Fig. 6-2) to model perforated horizontal component.After rotating, concentration of energy is to R
xand R
yenergy is almost nil, shows to obtain the relative geophone orientation angle (as shown in Figure 7 form) of perforation exactly.
Secondly, utilize perforation orientation angle, respectively orientation Concordance process (as shown in Fig. 8-1 to Fig. 8-2, Figure 10-1 to Figure 10-2) is carried out to upward traveling wave, down going wave two pressure break event horizontal components, choose C simultaneously
xcomponent is as reference road K.
Finally, relatively amplitude positive-negative polarity between pressure break data Z component and library track K, according to judgment criterion of the present invention, identify which seismic phase meets P ripple vertical coupled coordinate decomposition polarity rule, which seismic phase meets the vertical coupled coordinate of S ripple and decomposes polarity rule, automatic identification also isolates corresponding microseism data, and as shown in Fig. 9-1 to Fig. 9-2, Figure 11-1 is respectively upward traveling wave, down going wave different seismic phase recognition result comparison diagram with Figure 11-2.Can find out, per pass P ripple, S ripple seismic source information are clearly identified separation, and its result is consistent with two Source Model design theory, describes the accuracy of seismic phase recognizer of the present invention.
In addition, carry out borehole microseismic real data data seismic phase discriminance analysis, mainly verify that the present invention can to the automatic recognition capability of multiple micro-seismic event real data seismic phase.Figure 12-1 to Figure 12-3 is borehole microseismic real data perforation three-component section, and Figure 13-1 to Figure 13-3 is real data pressure break three-component section, similarly, according to Fig. 1, carries out operational processes.First, utilize hodograph-histogram method, obtain perforation and rotate rear two new components and counterparty's parallactic angle (as Figure 14-1 to Figure 14-2, Figure 15); Then, perforation orientation Concordance (as shown in Figure 16-1 to Figure 16-2) is carried out, by C to pressure break data horizontal component
xcomponent is defined as library track K; Secondly, amplitudes egales contrast is carried out to library track K and pressure break data Z component; Finally, according to polarity judgment criterion of the present invention, identify all pressure break event seismic phase types in pressure break data Z component, as shown in Figure 17-1 to Figure 17-2.Can find out, by phase identification method of the present invention, the most P ripple of pressure break data Z component, S wave energy are enough accurately identified and separate, and demonstrate the present invention and have comparatively strong, applicability.
Technique scheme is one embodiment of the present invention, for those skilled in the art, on the basis that the invention discloses application process and principle, be easy to make various types of improvement or distortion, and the method be not limited only to described by the above-mentioned embodiment of the present invention, therefore previously described mode is just preferred, and does not have restrictive meaning.
Claims (4)
1. a borehole microseismic phase identification method, is characterized in that: said method comprising the steps of:
(1), by known borehole microseismic perforation three component seismic data, correlation hole data X component and perforation data Y-component carry out the perforation orientation angle θ rotating the rear new components R x of process acquisition rotation and rotate rear newly components R y and the relative fractured well perforating site of monitor well three-component seismometer X component;
(2), utilize the perforation orientation angle θ that step (1) obtains, borehole microseismic pressure break data X, Y two component are carried out to orientation Concordance process and obtain polarity consistent pressure break new component Cx and the new component Cy of pressure break;
(3), from the C that step (2) obtains
x, C
yin, choose C
xas reference road K.
(4), by judging that the positive-negative polarity of library track K and the same time amplitude of borehole microseismic Z component identifies seismic phase type, then isolating the P ripple in Z component validity event and S ripple, finally exporting P ripple categorical data and S ripple categorical data.
2. borehole microseismic phase identification method according to claim 1, is characterized in that: the perforation orientation angle θ of the relative fractured well perforating site of monitor well three-component seismometer X component in described step (1) adopts horizontal component hodograph-histogram to ask for out.
3. borehole microseismic phase identification method according to claim 2, is characterized in that: described step (2) is decomposed by coordinate geometry and realized, specific as follows:
Input azimuth angle theta, re-starts vertical coordinate system geometric angle decompose pressure break data X, Y two component:
C
X＝x
i cos(θ)+y
i sin(θ)
C
Y＝-x
i sin(θ)+y
i cos(θ)
Wherein, x
i, y
ibe respectively X, Y-component amplitude, C
x, C
yfor two component amplitude values new after geometry value decomposition, and C
xdirection is that perforation is to monitor well wave detector horizontal projection direction.
4. borehole microseismic phase identification method according to claim 3, is characterized in that: described by judging that the positive-negative polarity of library track K and borehole microseismic Z component same time amplitude identifies that seismic phase type is specific as follows in described step (4):
If current wave field is upward traveling wave, as the amplitude A of library track K
kwith the amplitude A of borehole microseismic Z component
zbe multiplied when being greater than 0, i.e. A
ka
zduring > 0, positive-negative polarity performance is consistent, and now seismic phase is judged as P ripple, works as A
ka
zduring < 0, positive-negative polarity performance is contrary, and now seismic phase is judged as S ripple;
If current wave field is down going wave, judge that conclusion is then contrary with upward traveling wave.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310292019.3A CN104280772A (en) | 2013-07-12 | 2013-07-12 | Recognition method for microseism phase in well |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310292019.3A CN104280772A (en) | 2013-07-12 | 2013-07-12 | Recognition method for microseism phase in well |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104280772A true CN104280772A (en) | 2015-01-14 |
Family
ID=52255834
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310292019.3A Pending CN104280772A (en) | 2013-07-12 | 2013-07-12 | Recognition method for microseism phase in well |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104280772A (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016115763A1 (en) * | 2015-01-21 | 2016-07-28 | 中国矿业大学(北京) | Common-posture gather azimuthal angle analysis and correction method and device |
| CN106154332A (en) * | 2015-05-13 | 2016-11-23 | 中国石油化工股份有限公司 | A kind of borehole microseismic ripple event first arrival recognition methods in length and breadth |
| CN106154321A (en) * | 2016-07-20 | 2016-11-23 | 中国石油天然气集团公司 | The detection method of perforation signal and device |
| CN106168675A (en) * | 2015-05-18 | 2016-11-30 | 中国石油化工股份有限公司 | Borehole microseismic P/S ripple event recognition method and device |
| CN107843924A (en) * | 2016-09-21 | 2018-03-27 | 中国石油化工股份有限公司 | Utilize the seismic source location of moving constraint at the beginning of P ripples and focal mechanism joint inversion method |
| CN107918157A (en) * | 2016-10-08 | 2018-04-17 | 中国石油化工股份有限公司 | Three-component P ripples first motion focal mechanism inversion method and device |
| CN110531412A (en) * | 2019-09-27 | 2019-12-03 | 中国石油大学(北京) | A method of calculating borehole microseismic event relative bearing |
| CN112558147A (en) * | 2019-09-25 | 2021-03-26 | 中国石油化工股份有限公司 | Polarization analysis method and system for micro seismic data in well |
| CN114185081A (en) * | 2021-11-11 | 2022-03-15 | 北京科技大学 | Method for monitoring hydraulic fracturing range of ground straight well |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006106856A (en) * | 2004-09-30 | 2006-04-20 | Japan Petroleum Exploration Co Ltd | Three-dimensional image generation method and program |
| CN102841373A (en) * | 2012-08-23 | 2012-12-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Microseism positioning method based on azimuth angle constraint |
-
2013
- 2013-07-12 CN CN201310292019.3A patent/CN104280772A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006106856A (en) * | 2004-09-30 | 2006-04-20 | Japan Petroleum Exploration Co Ltd | Three-dimensional image generation method and program |
| CN102841373A (en) * | 2012-08-23 | 2012-12-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Microseism positioning method based on azimuth angle constraint |
Non-Patent Citations (6)
| Title |
|---|
| 何伟兵: "基于多分量的井间地震纵横波分离方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 刘希强，等: "用于三分向记录震相识别的小波变换方法", 《地震学报》 * |
| 宋维琪，等: "微地震资料自适应滤波方法研究", 《石油物探》 * |
| 朱卫星，等: "微地震信号的偏振—位置对比法震相分离技术", 《石油地球物理勘探》 * |
| 朱卫星，等: "自适应极化滤波在微地震信号处理中的应用", 《勘探地球物理进展》 * |
| 朱卫星: "微地震信号的震相分离", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2016115763A1 (en) * | 2015-01-21 | 2016-07-28 | 中国矿业大学(北京) | Common-posture gather azimuthal angle analysis and correction method and device |
| CN106154332A (en) * | 2015-05-13 | 2016-11-23 | 中国石油化工股份有限公司 | A kind of borehole microseismic ripple event first arrival recognition methods in length and breadth |
| CN106154332B (en) * | 2015-05-13 | 2018-08-31 | 中国石油化工股份有限公司 | A kind of borehole microseismic longitudinal and shear wave event first arrival recognition methods |
| CN106168675B (en) * | 2015-05-18 | 2018-05-08 | 中国石油化工股份有限公司 | Borehole microseismic P/S ripples event recognition method and device |
| CN106168675A (en) * | 2015-05-18 | 2016-11-30 | 中国石油化工股份有限公司 | Borehole microseismic P/S ripple event recognition method and device |
| CN106154321A (en) * | 2016-07-20 | 2016-11-23 | 中国石油天然气集团公司 | The detection method of perforation signal and device |
| CN107843924A (en) * | 2016-09-21 | 2018-03-27 | 中国石油化工股份有限公司 | Utilize the seismic source location of moving constraint at the beginning of P ripples and focal mechanism joint inversion method |
| CN107918157A (en) * | 2016-10-08 | 2018-04-17 | 中国石油化工股份有限公司 | Three-component P ripples first motion focal mechanism inversion method and device |
| CN107918157B (en) * | 2016-10-08 | 2019-07-23 | 中国石油化工股份有限公司 | Three-component P wave first motion focal mechanism inversion method and device |
| CN112558147A (en) * | 2019-09-25 | 2021-03-26 | 中国石油化工股份有限公司 | Polarization analysis method and system for micro seismic data in well |
| CN112558147B (en) * | 2019-09-25 | 2023-10-10 | 中国石油化工股份有限公司 | Polarization analysis method and system for microseism data in well |
| CN110531412A (en) * | 2019-09-27 | 2019-12-03 | 中国石油大学(北京) | A method of calculating borehole microseismic event relative bearing |
| CN114185081A (en) * | 2021-11-11 | 2022-03-15 | 北京科技大学 | Method for monitoring hydraulic fracturing range of ground straight well |
| CN114185081B (en) * | 2021-11-11 | 2023-10-27 | 北京科技大学 | Ground vertical well hydraulic fracturing range monitoring method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104280772A (en) | Recognition method for microseism phase in well | |
| CN106353792B (en) | Method suitable for positioning micro-seismic source of hydraulic fracturing | |
| CN104216008B (en) | Downhole fracturing microseismic event identification method | |
| CN109143374B (en) | Method and system for imaging scattering body around well | |
| CN110529087B (en) | Method and device for evaluating hydraulic fracturing effect of stratum | |
| CN104280775B (en) | Microseism monitoring and positioning method based on full-waveform vector offset superposition | |
| CN107490808B (en) | A kind of method for building up of high reliability seismic prospecting observation system | |
| CN102879817B (en) | Control method for acquiring underground crack information on basis of ground earthquake data | |
| CN102053261A (en) | Method for processing seismic data | |
| CN105785440A (en) | Dispersion curve extracting method for mine channel wave double-component seismic signal | |
| CN103758511A (en) | Method and device for identifying hidden reservoir through underground reverse time migration imaging | |
| CN105093301B (en) | The generation method and device of common imaging point angle of reflection angle gathers | |
| CN1837859A (en) | Three-dimensional seismic data processing quality monitoring technology | |
| CN110824563A (en) | Reservoir lithology prediction method based on Xgboost algorithm | |
| CN105093319A (en) | Ground micro-seismic static correction method based on three-dimensional seismic data | |
| CN102721979B (en) | Seismic data-based thin layer automatic interpretation and thickness prediction method and device | |
| CN104142518A (en) | Method for analyzing pre-stack time migration response of seismological observation system | |
| CN104316965A (en) | Prediction method and system for fissure azimuth and intensity | |
| CN106646609A (en) | Multi-scan microseism multi-parameter joint quick inversion method | |
| CN112305591B (en) | Tunnel advanced geological prediction method and computer readable storage medium | |
| CN107678064A (en) | A kind of interval transit time real time extracting method | |
| CN111045077A (en) | Full waveform inversion method of land seismic data | |
| CN105093314B (en) | A kind of method for determining microseism focus | |
| CN104360384B (en) | Micro-seismic event localization method based on wave energy automatically scanning in length and breadth and device | |
| CN103576197A (en) | Method for extracting converted wave angle channel set |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150114 |