CN106168678A - A kind of crack is propagated separation method and the system of shear wave - Google Patents
A kind of crack is propagated separation method and the system of shear wave Download PDFInfo
- Publication number
- CN106168678A CN106168678A CN201610562765.3A CN201610562765A CN106168678A CN 106168678 A CN106168678 A CN 106168678A CN 201610562765 A CN201610562765 A CN 201610562765A CN 106168678 A CN106168678 A CN 106168678A
- Authority
- CN
- China
- Prior art keywords
- component
- wave
- fast
- wave component
- slow
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
Abstract
The present invention provides separation method and the system propagating shear wave in a kind of crack, coordinate acquisition system is rotated according to set angle, if the anglec of rotation of coordinate acquisition system is the polarization azimuth of fast wave, coordinate acquisition system will overlap with natural system of coordinates.If the coordinate acquisition system anglec of rotation is β, postrotational radial component R (t) is slow wave component, and cross stream component T (t) is fast wave component.Now, fast, the energy ratio of slow wave component is the energy ratio measuring radial component R (t) and cross stream component T (t) obtained, and this ratio should tangent value azimuthal with fracture strike equal.This programme is with above-mentioned relation as theoretical basis, radial component R (t) obtained after analyzing rotating acquisition coordinate system and the energy ratio of cross stream component T (t) are the most equal with the tangent value of the anglec of rotation, when the two is equal, then the anglec of rotation can regard as the polarization azimuth of fast wave.The such scheme of the present invention, analytic operation process is simple, and agrees with mutually with existing theory and have the highest accuracy.
Description
Technical field
The present invention relates to seismic exploration technique field, be specifically related to a kind of crack is propagated the separation method of shear wave and is
System.
Background technology
When shear wave is propagated in the type anisotropic medium of crack, the polarization direction of incident shear wave changes, and is split into
Two mutually perpendicular independent shear waves in polarization direction.The two shear wave is propagated with different amplitudes and speed, is wherein parallel to split
The shear wave velocity of gap trend is very fast, and referred to as fast wave, by S1Represent;The shear wave velocity being perpendicular to fracture strike is relatively slow, referred to as slow wave,
By S2Represent.As in figure 2 it is shown, by fast wave S1With slow wave S2The crack azimuthal coordinates system that direction determines, referred to as natural system of coordinates, by
The coordinate system of single-point p-wave source and three-component (X, Y, Z) cymoscope composition is coordinate acquisition system.Generally coordinate acquisition system with from
So there is a unknown angle, i.e. fracture strike azimuth between coordinate system.Polarization azimuth and fracture strike due to fast wave
Azimuth is identical, as long as therefore obtaining the polarization azimuth of fast wave, separates Concerning With Fast-slow Waves, it becomes possible to Accurate Prediction fracture strike.Cause
This is precisely separating Concerning With Fast-slow Waves, is the key of existing fracture strike orientation judgement, and not yet has accuracy higher in prior art
Concerning With Fast-slow Waves separation method can be put in actual application.
Summary of the invention
The technical problem to be solved in the present invention is the Concerning With Fast-slow Waves separation method accuracy that in prior art, shear wave is propagated in crack
Low.
In order to solve above-mentioned technical problem, the present invention provides following technical scheme:
A kind of crack is propagated the separation method of shear wave, comprises the steps:
Rotating acquisition coordinate system, its anglec of rotation is counter clockwise direction angle [alpha];
On the basis of postrotational coordinate acquisition system, gather radial component R (t) and cross stream component T (t), the wherein t of shear wave
For the shear wave propagation time in crack;
Obtain energy ratio G of described radial component R (t) and cross stream component T (t) collected;
Judge whether anglec of rotation α meets tg (90 ° of-α)=G;
If being unsatisfactory for, after anglec of rotation α is adjusted, return the step of rotating acquisition coordinate system;
If meeting, obtain fast wave component S1(t) and slow wave component S2The expression formula of (t):
S1(t)=R (t) cos β-T (t) sin β;S2(t)=R (t) sin β+T (t) cos β;
Wherein, β=90 °-α, S1(t) and S2T () represents that fast wave component and slow wave component are when the propagation time is t pair respectively
Answer propagation distance;
According to fast wave component S1(t) and slow wave component S2T the expression formula of () judges fast wave component S1T whether () divide than slow wave
Amount S2T () first arrives ad-hoc location;
If then judging the β polarization azimuth as fast wave, returning rotating acquisition after otherwise anglec of rotation α being adjusted and sitting
The step of mark system.
Preferably, the separation method propagating shear wave in above-mentioned crack also comprises the steps:
According to fast wave component S1(t) and slow wave component S2T the expression formula of () obtains fast wave component and slow wave component propagation time
Difference.
Preferably, the separation method propagating shear wave in above-mentioned crack also comprises the steps:
Polaried orientation angle beta according to fast wave is to described fast wave component S1(t) and described slow wave component S2The expression formula of (t)
Amplitude is corrected processing.
Preferably, propagating the separation method of shear wave in above-mentioned crack, the polarization azimuth of fast wave is to described fast wave component
S1(t) and described slow wave component S2T the amplitude of the expression formula of () is corrected in the step processed, enter according to method for numerical simulation
Row correction process, obtains:
Fast wave component expression formula after correction: S01(t)=S1(t)/cosβ;
Slow wave weight expression after correction: S02(t)=S2(t)/sinβ。
Based on same inventive concept, the present invention also provides for the shear wave piece-rate system propagated in a kind of crack, including:
Rotary unit, rotating acquisition coordinate system, its anglec of rotation is counter clockwise direction angle [alpha];
Collecting unit, on the basis of postrotational coordinate acquisition system, gathers radial component R (t) and the cross stream component T of shear wave
T (), wherein t is the shear wave propagation time in crack;
Data capture unit, obtains energy ratio G of described radial component R (t) and cross stream component T (t) collected;
First judging unit, it is judged that whether anglec of rotation α meets tg (90 ° of-α)=G;
Processing unit, when the judged result of described judging unit is for being, obtains fast wave component S1(t) and slow wave component S2
The expression formula of (t): S1(t)=R (t) cos β-T (t) sin β;S2(t)=R (t) sin β+T (t) cos β;Wherein, β=90 °-α, S1
(t) and S2T () represents fast wave component and slow wave component corresponding propagation distance when the propagation time is t respectively;
Second judging unit, according to fast wave component S1(t) and slow wave component S2T the expression formula of () judges fast wave component S1(t)
Whether than slow wave component S2T () first arrives ad-hoc location;
Described rotary unit, after being adjusted anglec of rotation α when the judged result of described judging unit is no, rotates
Described coordinate acquisition system;
Described rotary unit, is no in the judged result of described first judging unit or the judged result of the second judging unit
Time, anglec of rotation α is adjusted rear rotating acquisition coordinate system;
Described processing unit, when the judged result of described second judging unit is for being, it is determined that β is the polarization side of fast wave
Parallactic angle, S1T ()=R (t) cos β-T (t) sin β is the fast wave expression formula after separating, S2T ()=R (t) sin β+T (t) cos β is for dividing
Slow wave expression formula after from.
Preferably, the shear wave piece-rate system propagated in above-mentioned crack, described data capture unit is additionally operable to according to fast wave
Component S1(t) and slow wave component S2T the expression formula of () obtains fast wave component and slow wave component propagation time difference.
Preferably, the shear wave piece-rate system propagated in above-mentioned crack, also include:
Correction unit, according to the polaried orientation angle beta of fast wave to described fast wave component S1(t) and described slow wave component S2(t)
The amplitude of expression formula is corrected processing.
Preferably, the shear wave piece-rate system propagated in above-mentioned crack, described correction unit, enters according to method for numerical simulation
Row correction process, the fast wave component expression formula after being corrected: S01(t)=S1(t)/cosβ;Slow wave component after correction is expressed
Formula: S02(t)=S2(t)/sinβ。
The technique scheme that the present invention provides, compared with prior art, at least has the advantages that
(1) crack of the present invention is propagated separation method and the system of shear wave, by coordinate acquisition system according to set angle
Degree rotates, if the anglec of rotation of coordinate acquisition system is the polarization azimuth (i.e. fracture strike azimuth) of fast wave, gathers
Coordinate system will overlap with natural system of coordinates.According to existing knowledge, can push away, if the coordinate acquisition system anglec of rotation is β, then
Postrotational radial component R (t) is slow wave component, and cross stream component T (t) is fast wave component.Now, fast, slow wave component
Energy ratio is the energy ratio measuring radial component R (t) and cross stream component T (t) obtained, and this ratio should be with crack
The tangent value of angle of strike complementary angle is equal.This programme, with above-mentioned relation as theoretical basis, obtains after analyzing rotating acquisition coordinate system
The energy ratio of radial component R (t) arrived and cross stream component T (t) is the most equal with the tangent value of anglec of rotation complementary angle, when the two
Time equal, then the anglec of rotation can regard as the polarization azimuth (i.e. fracture strike azimuth) of fast wave.The above-mentioned side of the present invention
Case, analytic operation process is simple, and agrees with mutually with existing theory and have the highest accuracy.
(2) crack of the present invention is propagated separation method and the system of shear wave, it is contemplated that different angles are likely to be of
Identical tangent value is therefore equal with the energy ratio of radial component and cross stream component at the tangent value judging anglec of rotation complementary angle
Time, according to the result obtained, the propagation relation of fast wave and slow wave is judged further.Due to fast wave spread speed ratio slowly
The spread speed of ripple is fast, therefore can get time, position relationship according to fast wave function expression and slow wave function expression,
Really when fast wave arrives ad-hoc location prior to slow wave, then it is believed that obtained result is accurately, otherwise result is inaccurate
True, the accuracy of separating resulting can be further ensured that by this step.
(3) crack of the present invention is propagated separation method and the system of shear wave, also includes fast wave component and slow wave
The amplitude of component is corrected the step processed, and is corrected fast wave and slow wave amplitude by method for numerical simulation, can obtain
Divide to incident along fissure direction the most accurately and polarization fast wave component with being perpendicular to incident and polarization the slow wave of fissure direction
Amount.
Accompanying drawing explanation
Fig. 1 is the flow chart of the separation method propagating shear wave in crack described in one embodiment of the invention;
Fig. 2 is angular relationship schematic diagram between coordinate acquisition system and Concerning With Fast-slow Waves;
Fig. 3 is crack discriminant curve;
Fig. 4 A, 4B are the azimuthal result Critical curve of fracture strike described in one embodiment of the invention;
Fig. 5 A-5D is that described in one embodiment of the invention, shear wave curve is transformed to R, T component, S successively by X, Y-component1、S2
Component, S01、S02The curve chart of component;
Fig. 6 invents the theory diagram of the separation method propagating shear wave in crack described in an embodiment;
Fig. 7 be the present invention one be coal geology model schematic described in embodiment;
Fig. 8 A is the wave field radial component record of Fig. 7 institute representation model;
Fig. 8 B is the wave field cross stream component record of Fig. 7 institute representation model;
Fig. 8 C is the wave field Z-direction component recording of Fig. 7 institute representation model;
Fig. 9 A, 9B are CCP Dao Ji azimuth superposition schematic diagram.
Reference therein is:
1-rotary unit, 2-collecting unit, 3-data capture unit, 4-the first judging unit, 5-processing unit, 6-second
Judging unit, 7-corrects unit.
Detailed description of the invention
Embodiment 1
The present embodiment provides the separation method propagating shear wave in a kind of crack, as it is shown in figure 1, comprise the steps:
S1: rotating acquisition coordinate system, its anglec of rotation is α, the described anglec of rotation in a counterclockwise direction on the basis of.This step
In, the direction of rotation for coordinate acquisition system does not has particular requirement, and the anglec of rotation simply finally given is in a counterclockwise direction
Obtain for standard, such as, can turn clockwise 240 °, finally give the anglec of rotation for rotation 120 °, the then anglec of rotation counterclockwise
Degree is 120 °.As another kind of metering method, it is believed that dextrorotation gyration is negative, and rotated counterclockwise by angle is just, example
As turned clockwise 30 °, it is believed that have rotated-30 ° counterclockwise.
S2: on the basis of postrotational coordinate acquisition system, gathers radial component R (t) and cross stream component T (t) of shear wave, its
Middle t is the shear wave propagation time in crack.In this step, available cymoscope of the prior art, can detect according to cymoscope
Obtain X (t) component and Y (t) component, according to X (t) component and Y (t) component of cymoscope, it is possible to obtain radial component R further
(t) and cross stream component T (t): R (t)=X (t) cos α-Y (t) sin α;T (t)=X (t) sin α+Y (t) cos α.
S3: obtain energy ratio G of described radial component R (t) and cross stream component T (t) collected.
S4: judge whether anglec of rotation α meets tg (90 ° of-α)=G;If being unsatisfactory for, enter step S8;If meeting, enter
Enter step S5.
S5: obtain fast wave component S1(t) and slow wave component S2The expression formula of (t): S1(t)=R (t) cos β-T (t) sin β;
S2(t)=R (t) sin β+T (t) cos β;Wherein, β=90 °-α, S1(t) and S2T () represents fast wave component and slow wave component respectively
The correspondence propagation distance when the propagation time is t.
S6: according to fast wave component S1(t) and slow wave component S2T the expression formula of () judges fast wave component S1T whether () be than slow wave
Component S2T () first arrives ad-hoc location;If then entering step S7, otherwise enter step S8.
S7: judge the β polarization azimuth as fast wave, S1T ()=R (t) cos β-T (t) sin β is that the fast wave after separating is expressed
Formula, S2T ()=R (t) sin β+T (t) cos β is the slow wave expression formula after separating.
S8: return step S1 after anglec of rotation α is adjusted.
Below in conjunction with the accompanying drawings the above results is derived in detail.In order to the Concerning With Fast-slow Waves after quantitative description shear wave splitting with point
Split the relation between front shear wave, illustrate on the basis of the coordinate acquisition system in Fig. 2.Wherein fracture strike and line direction
Angle is β, the present embodiment is ignored transmission potential loss, when shear wave vertical up time, it is assumed that its particle vibration direction and survey
Line direction is parallel, and the displacement at moment t is S, after entering upright creviced formation, divides the fast wave S1With slow wave S2Can be by public affairs
Formula (1) represents:
Obviously, X-component and the Y-component of physical record is:
From Fig. 1 it is seen that, if by 90 ° of-β of XOY coordinate system rotated counterclockwise by angle, postrotational components R is slow wave, point
Amount T is fast wave, component X after rotationRWith component YRExpression formula be:
Can verify that and obtain:
Assuming that the anglec of rotation is 90 ° of-α, after rotation, the expression formula of components R and component T is:
Convolution (1)-(5), obtain after the expression formula relevant with components R and component T being arranged:
It is hereby achieved that:
As α=β, if only considering fast, the energy of slow wave, calculating components R and energy ratio G of component T, then having:
Formula (8) illustrates, when coordinate acquisition system is rotated to natural system of coordinates, after rotation the energy of seismic wave than with rotation
The tangent value of the complementary angle of angle is equal, therefore, it can ask for the growth main formula position in crack, underground, i.e. fracture strike according to formula (8)
Azimuth.
The such scheme of the present embodiment, rotates coordinate acquisition system according to set angle, if the rotation of coordinate acquisition system
Gyration is the polarization azimuth (i.e. fracture strike azimuth) of fast wave, and coordinate acquisition system will overlap with natural system of coordinates.
According to existing knowledge, can push away, if the coordinate acquisition system anglec of rotation is β, then postrotational radial component R (t) is slowly
Wave component, cross stream component T (t) is fast wave component.Now, fast, the energy ratio of slow wave component is and measures the radial direction obtained
The energy ratio of components R (t) and cross stream component T (t), and this ratio should tangent value azimuthal with fracture strike equal.This
Scheme, with above-mentioned relation as theoretical basis, analyzes radial component R (t) and cross stream component T (t) obtained after rotating acquisition coordinate system
Energy ratio the most equal with the tangent value of the complementary angle of the anglec of rotation, when the two is equal, then the anglec of rotation can be regarded as
The polarization azimuth (i.e. fracture strike azimuth) of fast wave.The such scheme of the present embodiment, analytic operation process is simple, and with
Theoretical agreeing with mutually is had to have the highest accuracy.
During actual asking for, by a series of β value, components R and component T are changed, and corresponding each β value calculates phase
Energy ratio G answered and tangent value, formation G (β) and two curves of tg β, referred to as discriminant curve, as shown in Figure 3.Article two, song
The β value of line point of intersection is the azimuth of stratum fracture strike, namely the polarization azimuth of fast wave.
Time from figure it can also be seen that utilize said method to carry out crack orientation detection, in whole 0~180 ° of interval, have two
Individual solution 30 ° and 120 °, in order to differentiate the correctness solved, has carried out numerical simulation study.Fig. 4 A show fracture strike
When azimuth is 30 °, Concerning With Fast-slow Waves propagation curve, is learnt by relative analysis, fast wave S1Ratio slow wave S2First arrive ad-hoc location.Fig. 4 B
It show Concerning With Fast-slow Waves propagation curve when fracture strike azimuth is 120 °, learnt by relative analysis, fast wave S1Ratio slow wave S2After arrive
Reach ad-hoc location.Obviously the result rotating 30 ° is correct, and the result of 120 ° is incorrect.
It is further preferred that the said method of the present embodiment also comprises the steps:
S9: according to fast wave component S1(t) and slow wave component S2T the expression formula of () obtains fast wave component and slow wave component is propagated
Time difference.
As shown in Figure 4 A, once it is determined that the propagation curve of fast wave and slow wave, it is easy to obtain the propagation of fast wave and slow wave
Time difference, shown in Fig. 4, the propagation time difference of the two is 6s.
It is further preferred that the said method of the present embodiment also comprises the steps:
S10: according to the polaried orientation angle beta of fast wave to described fast wave component S1(t) and described slow wave component S2The expression of (t)
The amplitude of formula is corrected processing.It is corrected processing according to method for numerical simulation, obtains: the fast wave component after correction is expressed
Formula: S01(t)=S1(t)/cosβ;Slow wave weight expression after correction: S02(t)=S2(t)/sinβ。
The fast wave component utilizing two-dimentional two component circumgyrations to obtain is relevant with source receiver azimuth with the amplitude of slow wave component, and
Fast wave and vertical the crack incident and slow wave of polarization the amplitude with polarization incident along orientation, crack obtained is wanted only to depend on eventually
Unrelated with source receiver azimuth in source level.Use method for numerical simulation to carry out research to find according to above formula S1(t) and S2(t)
Component carries out amplitude recovery correction, can obtain and the fast wave S of polarization incident along fissure direction01T () and vertical fissure direction are incident
Slow wave S with polarization02(t).The conversion process of the component of the above shear wave of the present embodiment can be according to shown in Fig. 5 A to Fig. 5 D
Curve obtains more intuitively.
Embodiment 2
The shear wave piece-rate system propagated in a kind of crack of the present embodiment offer, as shown in Figure 6, including:
Rotary unit 1, rotating acquisition coordinate system, its anglec of rotation is counter clockwise direction angle [alpha];
Collecting unit 2, on the basis of postrotational coordinate acquisition system, gathers radial component R (t) and the cross stream component of shear wave
T (t), wherein t is the shear wave propagation time in crack;
Data capture unit 3, obtains energy ratio G of described radial component R (t) and cross stream component T (t) collected;
First judging unit 4, it is judged that whether anglec of rotation α meets tg (90 ° of-α)=R;
Processing unit 5, when the judged result of described judging unit 4 is for being, obtains fast wave component S1(t) and slow wave component
S2The expression formula of (t): S1(t)=R (t) cos β-T (t) sin β;S2(t)=R (t) sin β+T (t) cos β;Wherein, β=90 °-α,
S1(t) and S2T () represents fast wave component and slow wave component corresponding propagation distance when the propagation time is t respectively.
Second judging unit 6, according to fast wave component S1(t) and slow wave component S2T the expression formula of () judges fast wave component S1
T whether () be than slow wave component S2T () first arrives ad-hoc location;
Described rotary unit 1, in judged result or the judged result of the second judging unit 6 of described first judging unit 4
For time no, anglec of rotation α is adjusted rear rotating acquisition coordinate system;
Described processing unit 5, when the judged result of described second judging unit 6 is for being, it is determined that β is the polarization side of fast wave
Parallactic angle, S1T ()=R (t) cos β-T (t) sin β is the fast wave expression formula after separating, S2T ()=R (t) sin β+T (t) cos β is for dividing
Slow wave expression formula after from.The such scheme of the present embodiment is likely to be of identical tangent value in view of different angles, therefore exists
Judge when the tangent value of anglec of rotation complementary angle is equal with the energy ratio of radial component and cross stream component, further according to obtaining
The propagation relation of fast wave and slow wave is judged by result.Owing to the spread speed of fast wave is faster than the spread speed of slow wave, therefore
Can obtain the relation of the two propagation time, propagation distance according to fast wave expression formula and slow wave expression formula, certain fast wave is prior to slowly
When ripple arrives ad-hoc location, then it is believed that obtained result is accurately, otherwise result is inaccurate, by this step
The accuracy of separating resulting can be further ensured that.
Further, described data capture unit 3 is additionally operable to according to fast wave component S1(t) and slow wave component S2The expression of (t)
Formula obtains fast wave component and slow wave component propagation time difference.Once it is determined that the propagation curve of fast wave and slow wave, it is easy to obtain
Fast wave and the propagation time difference of slow wave.
As it is shown in figure 1, as preferred scheme, said system also includes: correction unit 7, according to the polarization side of fast wave
Parallactic angle β is to described fast wave component S1(t) and described slow wave component S2T the amplitude of the expression formula of () is corrected processing.According to numerical value
Analogy method is corrected processing, the fast wave component expression formula after being corrected: S01(t)=S1(t)/cosβ;Slow after correction
Wave component expression formula: S02(t)=S2(t)/sinβ。
The such scheme of the present embodiment, rotates coordinate acquisition system according to set angle, if the rotation of coordinate acquisition system
Gyration is the polarization azimuth (i.e. fracture strike azimuth) of fast wave, and coordinate acquisition system will overlap with natural system of coordinates.
According to existing knowledge, can push away, if the coordinate acquisition system anglec of rotation is β, then postrotational radial component R (t) is slowly
Wave component, cross stream component T (t) is fast wave component.Now, fast, the energy ratio of slow wave component is and measures the radial direction obtained
Components R (t) and the energy ratio of cross stream component T (t), and this ratio should be with the tangent value phase of fracture strike azimuth complementary angle
Deng.This programme is with above-mentioned relation as theoretical basis, and radial component R (t) obtained after analyzing rotating acquisition coordinate system is divided with horizontal
The energy ratio of amount T (t) is the most equal with the tangent value of anglec of rotation complementary angle, and when the two is equal, then the anglec of rotation can be recognized
It is set to the polarization azimuth (i.e. fracture strike azimuth) of fast wave.The such scheme of the present embodiment, analytic operation process is simple, and
Agree with mutually with existing theory and there is the highest accuracy.
Embodiment 3
The present embodiment provides the identification marking in a kind of cracks in coal seam orientation, splits to set up coal seam transformed wave wave field and coal seam
Variation relation between gap, carries out forward simulation to cracks in coal seam model.Coal geology model is as it is shown in fig. 7, all stratum are equal
For horizontal layer, cracks in coal seam is vertical crack parallel to each other, moves towards parallel with Y-axis, the parameter such as table 1 that each layer is corresponding.
Table 1 coal geology model physical parameter
Excite at O point, with O point as the center of circle, radius as r on the round survey line of (r=1400m in model) receive, crack side
Parallactic angle θ is started at X-axis (i.e. direct north), uses the wave field of reflectivity method computation model, and wavelet elects the Ricker wavelet of 60Hz as,
Obtain wave field record as shown in figures 8 a-8 c.
From the point of view of wave field record, as shown in figs. 8 a and 8b, when converted shear wave is through coal seam, walk with crack according to its incident orientation
To relativeness and there is different division situations: during when incident orientation vertical fracture strike (azimuth is 0 °), radially
In components R, shear wave whilst on tour is the longest, and on cross stream component T, amplitude is zero, only produces slow wave;When the parallel fracture strike of incident orientation
(i.e. azimuth is 90 °), in radial component R, shear wave whilst on tour is the shortest, and on cross stream component T, amplitude is zero, only produces fast wave;At it
Then can observe fast transverse wave and two kinds of ripples of slow shear-wave in its orientation, in the case of small incident, the arrival time difference of two kinds of ripples is with it
Become big through the lengthening in path, coal seam, when incident orientation and fracture strike are 45 ° of (or 225 °) angles, cross stream component T amplitude
Reaching maximum, incident orientation is every 90 °, and the phase place of cross stream component T occurs 180 ° of reversings.
As shown in Figure 8 C, when compressional wave is through coal seam, the Orientation differences of speed also occurs: when incident orientation is put down with fracture strike
During row, speed is fast, and amplitude is strong, and time vertical, speed is slow, and amplitude is weak.This azimuthal anisotropy feature is only on geophone offset farther out
Could substantially observe.
According to coal seam transformed wave wave field characteristics, it can be deduced that the recognition marks in cracks in coal seam orientation, meet radially the most simultaneously
The angle that components R shear wave whilst on tour is the shortest, cross stream component T amplitude is zero is exactly coal seam fast wave polarization azimuth.
In Three-dimendimal fusion earthquake data acquisition, there is different source receiver azimuth angles in common-conversion point gather Nei Ge road, because of
And different angles is had from fracture strike, result causes radial component and the amplitude of cross stream component that each road is corresponding and phase place all to produce
Raw corresponding change.For some common-conversion point gather (i.e. CCP road collection), can be formed different orientations radial component and
Cross stream component stack result, as it is shown in figure 9, radial component and cross stream component are respectively arranged with 36 roads, each road is a certain particular azimuth
In the range of the stack result in each road, the 1st road represents the superposition in each road, the 2nd road generation in the range of all source receiver azimuth angles are 0~10 °
Table all source receiver azimuths angle be 11~20 ° in the range of the superposition in each road, by that analogy.
Those skilled in the art are it should be appreciated that embodiments of the invention can be provided as method, system or computer program
Product.Therefore, the reality in terms of the present invention can use complete hardware embodiment, complete software implementation or combine software and hardware
Execute the form of example.And, the present invention can use at one or more computers wherein including computer usable program code
The upper computer program product implemented of usable storage medium (including but not limited to disk memory, CD-ROM, optical memory etc.)
The form of product.
The present invention is with reference to method, equipment (system) and the flow process of computer program according to embodiments of the present invention
Figure and/or block diagram describe.It should be understood that can the most first-class by computer program instructions flowchart and/or block diagram
Flow process in journey and/or square frame and flow chart and/or block diagram and/or the combination of square frame.These computer programs can be provided
Instruction arrives the processor of general purpose computer, special-purpose computer, Embedded Processor or other programmable data processing device to produce
A raw machine so that the instruction performed by the processor of computer or other programmable data processing device is produced for real
The device of the function specified in one flow process of flow chart or multiple flow process and/or one square frame of block diagram or multiple square frame now.
These computer program instructions may be alternatively stored in and computer or other programmable data processing device can be guided with spy
Determine in the computer-readable memory that mode works so that the instruction being stored in this computer-readable memory produces and includes referring to
Make the manufacture of device, this command device realize at one flow process of flow chart or multiple flow process and/or one square frame of block diagram or
The function specified in multiple square frames.
These computer program instructions also can be loaded in computer or other programmable data processing device so that at meter
Perform sequence of operations step on calculation machine or other programmable devices to produce computer implemented process, thus at computer or
The instruction performed on other programmable devices provides for realizing at one flow process of flow chart or multiple flow process and/or block diagram one
The step of the function specified in individual square frame or multiple square frame.
Although preferred embodiments of the present invention have been described, but those skilled in the art once know basic creation
Property concept, then can make other change and amendment to these embodiments.So, claims are intended to be construed to include excellent
Select embodiment and fall into all changes and the amendment of the scope of the invention.
Claims (8)
1. a crack is propagated the separation method of shear wave, it is characterised in that comprise the steps:
Rotating acquisition coordinate system, its anglec of rotation is counter clockwise direction angle [alpha];
On the basis of postrotational coordinate acquisition system, gathering radial component R (t) and cross stream component T (t) of shear wave, wherein t is horizontal
The ripple propagation time in crack;
Obtain energy ratio G of described radial component R (t) and cross stream component T (t) collected;
Judge whether anglec of rotation α meets tg (90 ° of-α)=G;
If being unsatisfactory for, after anglec of rotation α is adjusted, return the step of rotating acquisition coordinate system;
If meeting, obtain fast wave component S1(t) and slow wave component S2The expression formula of (t):
S1(t)=R (t) cos β-T (t) sin β;S2(t)=R (t) sin β+T (t) cos β;
Wherein, β=90 °-α, S1(t) and S2T () represents that fast wave component and slow wave component are corresponding when the propagation time is t respectively and passes
Broadcast distance;
According to fast wave component S1(t) and slow wave component S2T the expression formula of () judges fast wave component S1T whether () be than slow wave component S2
T () first arrives ad-hoc location;
If then judging the β polarization azimuth as fast wave, S1T ()=R (t) cos β-T (t) sin β is that the fast wave after separating is expressed
Formula, S2T ()=R (t) sin β+T (t) cos β is the slow wave expression formula after separating;Return after otherwise anglec of rotation α being adjusted
The step of rotating acquisition coordinate system.
Crack the most according to claim 1 is propagated the separation method of shear wave, it is characterised in that also comprise the steps:
According to fast wave component S1(t) and slow wave component S2T the expression formula of () obtains fast wave component and slow wave component propagation time difference.
Crack the most according to claim 1 is propagated the separation method of shear wave, it is characterised in that also comprise the steps:
Polaried orientation angle beta according to fast wave is to described fast wave component S1(t) and described slow wave component S2The amplitude of the expression formula of (t)
It is corrected processing.
Crack the most according to claim 3 is propagated the separation method of shear wave, it is characterised in that:
The polarization azimuth of fast wave is to described fast wave component S1(t) and described slow wave component S2T the amplitude of the expression formula of () carries out school
In the step just processed, it is corrected processing according to method for numerical simulation, obtains:
Fast wave component expression formula after correction: S01(t)=S1(t)/cosβ;
Slow wave weight expression after correction: S02(t)=S2(t)/sinβ。
5. the shear wave piece-rate system propagated in a crack, it is characterised in that including:
Rotary unit, rotating acquisition coordinate system, its anglec of rotation is counter clockwise direction angle [alpha];
Collecting unit, on the basis of postrotational coordinate acquisition system, gathers radial component R (t) and cross stream component T (t) of shear wave,
Wherein t is the shear wave propagation time in crack;
Data capture unit, obtains energy ratio G of described radial component R (t) and cross stream component T (t) collected;
First judging unit, it is judged that whether anglec of rotation α meets tg (90 ° of-α)=G;
Processing unit, when the judged result of described judging unit is for being, obtains fast wave component S1(t) and slow wave component S2(t)
Expression formula: S1(t)=R (t) cos β-T (t) sin β;S2(t)=R (t) sin β+T (t) cos β;Wherein, β=90 °-α, S1(t) and
S2T () represents fast wave component and slow wave component corresponding propagation distance when the propagation time is t respectively;
Second judging unit, according to fast wave component S1(t) and slow wave component S2T the expression formula of () judges fast wave component S1T whether ()
Than slow wave component S2T () first arrives ad-hoc location;
Described rotary unit, when the judged result of described first judging unit or the judged result of the second judging unit are no,
Anglec of rotation α is adjusted rear rotating acquisition coordinate system;
Described processing unit, when the judged result of described second judging unit is for being, it is determined that β is the polarization azimuth of fast wave, S1
T ()=R (t) cos β-T (t) sin β is the fast wave expression formula after separating, S2T ()=R (t) sin β+T (t) cos β is separation after
Slow wave expression formula.
The shear wave piece-rate system propagated in crack the most according to claim 5, it is characterised in that:
Described data capture unit is additionally operable to according to fast wave component S1(t) and slow wave component S2T the expression formula of () obtains fast wave component
With slow wave component propagation time difference.
The shear wave piece-rate system propagated in crack the most according to claim 5, it is characterised in that also include:
Correction unit, according to the polaried orientation angle beta of fast wave to described fast wave component S1(t) and described slow wave component S2The expression of (t)
The amplitude of formula is corrected processing.
The shear wave piece-rate system propagated in crack the most according to claim 7, it is characterised in that:
Described correction unit, is corrected processing according to method for numerical simulation, the fast wave component expression formula after being corrected: S01
(t)=S1(t)/cosβ;Slow wave weight expression after correction: S02(t)=S2(t)/sinβ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610562765.3A CN106168678B (en) | 2016-07-15 | 2016-07-15 | The separation method and system of shear wave are propagated in a kind of crack |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610562765.3A CN106168678B (en) | 2016-07-15 | 2016-07-15 | The separation method and system of shear wave are propagated in a kind of crack |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106168678A true CN106168678A (en) | 2016-11-30 |
CN106168678B CN106168678B (en) | 2018-08-17 |
Family
ID=58065944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610562765.3A Active CN106168678B (en) | 2016-07-15 | 2016-07-15 | The separation method and system of shear wave are propagated in a kind of crack |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106168678B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080062814A1 (en) * | 2006-09-12 | 2008-03-13 | Schlumberger Technology Corporation | Discriminating natural fracture- and stress-induced sonic anisotropy using a combination of image and sonic logs |
CN102096097A (en) * | 2009-12-14 | 2011-06-15 | 武侯区巅峰机电科技研发中心 | Method for testing fracture azimuth of fast and slow transverse wave splitting |
CN103076632A (en) * | 2012-09-20 | 2013-05-01 | 中国石油天然气集团公司 | Method and device for detecting developing degree of cracks in oil and gas reservoir |
CN104898162A (en) * | 2014-03-06 | 2015-09-09 | 中国石油化工股份有限公司 | Crack detection method in geological exploration |
CN105116448A (en) * | 2015-08-11 | 2015-12-02 | 中国石油天然气集团公司 | Conversion wave azimuthal anisotropy correction method and device thereof |
-
2016
- 2016-07-15 CN CN201610562765.3A patent/CN106168678B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080062814A1 (en) * | 2006-09-12 | 2008-03-13 | Schlumberger Technology Corporation | Discriminating natural fracture- and stress-induced sonic anisotropy using a combination of image and sonic logs |
CN102096097A (en) * | 2009-12-14 | 2011-06-15 | 武侯区巅峰机电科技研发中心 | Method for testing fracture azimuth of fast and slow transverse wave splitting |
CN103076632A (en) * | 2012-09-20 | 2013-05-01 | 中国石油天然气集团公司 | Method and device for detecting developing degree of cracks in oil and gas reservoir |
CN104898162A (en) * | 2014-03-06 | 2015-09-09 | 中国石油化工股份有限公司 | Crack detection method in geological exploration |
CN105116448A (en) * | 2015-08-11 | 2015-12-02 | 中国石油天然气集团公司 | Conversion wave azimuthal anisotropy correction method and device thereof |
Non-Patent Citations (2)
Title |
---|
何兵寿 等: ""煤体中直立裂隙的多波地震响应及预测"", 《地质学报》 * |
滕佃波 等: ""基于地震横波分裂理论的火成岩裂缝检测"", 《地球物理学进展》 * |
Also Published As
Publication number | Publication date |
---|---|
CN106168678B (en) | 2018-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105467440B (en) | A kind of omnidirectional's vector seismic data processing technique and device | |
CN102053261B (en) | Method for processing seismic data | |
CN102540250B (en) | Azimuth fidelity angle domain imaging-based fractured oil and gas reservoir seismic exploration method | |
CN1313841C (en) | Layer stripping converted reflected waveforms for dipping fractures | |
CN102053277A (en) | Method for detecting reservoir fissure development direction by utilizing seismic data | |
CN102879800B (en) | Method for detecting shear-wave splitting fracture | |
CN107678063B (en) | A kind of multi-component converted wave crack prediction method based on Rank correlation | |
CN104297784A (en) | Primary wave azimuthal anisotropy based fracture predicting method | |
CN107817523B (en) | The analysis method and device of diffracted wave migration velocity | |
CN104199103B (en) | Static correction method and static correction device for converted wave | |
CN105093281A (en) | Earthquake multi-wave modeling method under inverse framework | |
CN109212589A (en) | It is a kind of to cooperate with parallel earthquake-capturing observation system design method based on GPU/CPU | |
CN104570111A (en) | Common-posture gather azimuthal angle analyzing and correcting method and device | |
CN107132575A (en) | The method that fracture azimuth is predicted based on shear wave polarographic analysis | |
CN106468781A (en) | A kind of crack prediction method based on minimum entropy circumgyration | |
Dunn et al. | High-resolution earthquake relocation in the New Madrid seismic zone | |
CN103869366B (en) | A kind of method and device determining crack fracture strike | |
CN104199098B (en) | Method for predicting underground cracks by aid of polarization characteristics of reflected and converted shear waves | |
LU et al. | Imaging techniques of multi-component seismic data | |
CN105301638B (en) | A kind of method and apparatus in extraction base of weathering face | |
CN104898162B (en) | Crack detection method in geological exploration | |
US20190018156A1 (en) | Highly accurate focal mechanism for microseismic envents | |
Li et al. | SV-P extraction and imaging for far-offset vertical seismic profile data | |
CN102798888B (en) | Method for calculating velocity ratio of longitudinal wave to transverse wave by using non-zero wellhead distance data | |
CN106168678B (en) | The separation method and system of shear wave are propagated in a kind of crack |
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 |