CN107861154A - Fault strike recognition methods and device - Google Patents
Fault strike recognition methods and device Download PDFInfo
- Publication number
- CN107861154A CN107861154A CN201711089298.8A CN201711089298A CN107861154A CN 107861154 A CN107861154 A CN 107861154A CN 201711089298 A CN201711089298 A CN 201711089298A CN 107861154 A CN107861154 A CN 107861154A
- Authority
- CN
- China
- Prior art keywords
- mrow
- edge diffraction
- msub
- diffraction coefficient
- azimuth
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000003325 tomography Methods 0.000 claims abstract description 29
- 230000008859 change Effects 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000010304 firing Methods 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 15
- 238000005516 engineering process Methods 0.000 abstract description 10
- 230000008569 process Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000019637 foraging behavior Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a kind of fault strike recognition methods and device, it is related to the technical field of geological prospecting, this method includes:The edge diffraction coefficient in cylindrical coordinate is obtained, edge diffraction coefficient is edge diffraction coefficient caused by tomography edge;Edge diffraction coefficient in cylindrical coordinate is changed, be converted to the edge diffraction coefficient in spherical coordinate system, edge diffraction coefficient in spherical coordinate system includes multiple angle informations, wherein, edge diffraction coefficient in cylindrical coordinate includes outgoing angle information, and multiple angle informations include outgoing angle information and incidence angle information;Coordinate based on shot point in spherical coordinate system calculates the azimuth of shot point, and/or, the coordinate based on geophone station in spherical coordinate system calculates the azimuth of geophone station;Fault strike is determined based on the changing rule between azimuth and edge diffraction coefficient, when the application is alleviated by the progress fault strike judgement of existing sectional analysis technology, the poor technical problem of accuracy of analysis.
Description
Technical field
The present invention relates to the technical field of geological prospecting, more particularly, to a kind of fault strike recognition methods and device.
Background technology
When seismic wave is propagated in the earth formation, because the geologic structure of complexity and the suddenly change of lithology can produce diffraction
Ripple, therefore, the important information of a large amount of heterogeneous geologic bodies is often carried in diffracted wave.In seismic prospecting, due to using
Reflective wave method carrys out identifying purpose layer, ignores the little structure geologic body information of diffracted wave carrying, for identification tomography, pinching, vertical
The levels of precision of the little structure geologic body such as contact layer has certain influence.For a long time, people surround the accurate knowledge of tomographic systems
A large amount of effort are not made that, are proposed and are applied many recognition methods, and are successfully extracted perhaps from 3D data volume
The discontinuous attribute volume of more prominent fault informations, such as attribute of coherent data volume technology and ant body tracer technique.Wherein, coherence properties
The it is proposed of analytical technology make it that the noncoherent boundary analytical technology towards lithology has brought up to a brand-new level, but it is not
All influenceed with degree by analytical parameters, and the influence come with stratum remnants in response to belt that can not fundamentally abate the noise.
Ant body algorithm is a kind of chess game optimization algorithm by simulating ant colony foraging behavior developed in recent years, it is necessary to very
Strong computation capability.Above method is that the coherence of basis signal differentiates the presence of tomography mostly, does not consider side
Edge diffraction coefficient is with azimuthal variation.Therefore, during by existing analytical technology identifying fault strike, identification often occurs
Inaccurate phenomenon.In view of the above-mentioned problems, do not propose effective solution also.
The content of the invention
In view of this, it is an object of the invention to provide a kind of fault strike recognition methods and device, passed through with alleviating
When existing sectional analysis technology carries out fault strike identification, the poor technical problem of recognition accuracy.
In a first aspect, the embodiments of the invention provide a kind of fault strike recognition methods, including:Obtain in cylindrical coordinate
Edge diffraction coefficient, the edge diffraction coefficient are edge diffraction coefficient caused by tomography edge;To in the cylindrical coordinate
Edge diffraction coefficient is changed, and is converted to the edge diffraction coefficient in spherical coordinate system, the edge diffraction in the spherical coordinate system
Coefficient includes multiple angle informations, wherein, the edge diffraction coefficient in the cylindrical coordinate includes outgoing angle information, described more
Individual angle information includes outgoing angle information and incidence angle information;Coordinate based on shot point in the spherical coordinate system calculates the big gun
The azimuth of point, and/or, the coordinate based on geophone station in the spherical coordinate system calculates the azimuth of the geophone station;It is based on
Changing rule between the azimuth and the edge diffraction coefficient determines fault strike.
Further, the edge diffraction coefficient in the cylindrical coordinate is changed, is converted to the side in spherical coordinate system
Edge diffraction coefficient includes:The first coordinate information and the second coordinate information are obtained, wherein, first coordinate information is the shot point
Coordinate information in the spherical coordinate system, second coordinate information are coordinate of the geophone station in the spherical coordinate system
Information, the coordinate information include:The angle of emergence of diffracted wave, the incidence angle of diffracted wave, the shot point or the geophone station
The distance between origin of coordinates;Based on first coordinate information and second coordinate information in the cylindrical coordinate
Edge diffraction coefficient is changed, and is converted to the edge diffraction coefficient in spherical coordinate system.
Further, based on first coordinate information and second coordinate information to the edge in the cylindrical coordinate
Diffraction coefficient is changed, and the edge diffraction coefficient be converted in spherical coordinate system includes:According to formula in the cylindrical coordinate
Edge diffraction coefficient changed, be converted to the edge diffraction coefficient in spherical coordinate system, wherein, the formula is described as follows:Wherein, RL=rS+rR, (rS,θS,φS) it is described
First coordinate information, (rR,θR,φR) it is second coordinate information, RSRRepresent the shot point with the geophone interval from RIR
Represent the distance of image source point and geophone station, RLPassed through for the ray road of diffracted wave.
Further, the azimuth of the coordinate calculating shot point based on shot point in the spherical coordinate system includes:Pass through
FormulaThe azimuth of the shot point is calculated, wherein,For the orientation of the shot point
Angle, rSFor the distance between shot point and the origin of coordinates, θSFor the incidence angle of the diffracted wave, φSFor the angle of emergence of diffracted wave, its
In, θSAnd φSAngular range be 0 to 180 °,Angular range be 0 to 360 °.
Further, the azimuth of the coordinate calculating geophone station based on geophone station in the spherical coordinate system includes:
Pass through formulaThe azimuth of the geophone station is calculated, wherein,For the detection
The azimuth of point, rRFor the distance between the geophone station and the origin of coordinates, θRFor the incidence angle of the diffracted wave, φRFor diffraction
The angle of emergence of ripple, wherein, θRAnd φRAngular range be 0 to 180 °,Angular range be 0 to 360 °.
Further, fault strike bag is determined based on the changing rule between the azimuth and the edge diffraction coefficient
Include:The corresponding relation between the azimuth and the edge diffraction coefficient is determined, wherein, the corresponding relation is:When positioned at
During the surface of the tomography, if incidence angle is φSFor 90 °, then the edge diffraction coefficient VT(t) maximum is reached;Hanging down
Directly on the fault strike and section by the edge diffraction coefficient of the tomography, edge can occur directly over Diffraction Point
The phenomenon of diffraction coefficient polarity inversion;Work as output angle θR=90 °, incidence angle φSAt=90 °, parallel to the institute of the fault strike
State edge diffraction coefficient VT(t) maximum in three dimensions is reached;Graph of a relation is built based on the corresponding relation, wherein, institute
Stating graph of a relation includes changing rule between diffracted wave and the azimuth;Determine that the tomography is walked based on the graph of a relation
To.
Further, it is determined that the corresponding relation between the azimuth and the edge diffraction coefficient includes:Determine first
Corresponding relation, wherein, first corresponding relation is corresponding between the azimuth of the shot point and the edge diffraction coefficient
Relation;And/or the second corresponding relation is determined, wherein, azimuth and the side of second corresponding relation for the geophone station
Corresponding relation between edge diffraction coefficient.
Further, included based on corresponding relation structure graph of a relation:Based on first corresponding relation and described
Two corresponding relations build the graph of a relation.
Second aspect, the embodiments of the invention provide a kind of fault strike identification device, including:Acquiring unit, for obtaining
The edge diffraction coefficient in cylindrical coordinate is taken, the edge diffraction coefficient is edge diffraction coefficient caused by tomography edge;Conversion
Unit, for being changed to the edge diffraction coefficient in the cylindrical coordinate, be converted to the edge diffraction system in spherical coordinate system
Number, the edge diffraction coefficient in the spherical coordinate system include multiple angle informations, wherein, edge in the cylindrical coordinate around
Penetrate coefficient and include outgoing angle information, the multiple angle information includes outgoing angle information and incidence angle information;Computing unit, it is used for
Coordinate based on shot point in the spherical coordinate system calculates the azimuth of the shot point, and/or, sat based on geophone station in the ball
Coordinate in mark system calculates the azimuth of the geophone station;Determining unit, for based on the azimuth and the edge diffraction
Changing rule between coefficient determines fault strike.
The third aspect, the embodiments of the invention provide a kind of meter for the non-volatile program code that can perform with processor
Calculation machine computer-readable recording medium, described program code make computing device method described above.
In embodiments of the present invention, the edge diffraction coefficient in cylindrical coordinate is obtained first, then, in cylindrical coordinate
Edge diffraction coefficient is changed, and obtains the edge diffraction coefficient in spherical coordinate system, next, based on shot point in spherical coordinate system
Coordinate calculate the azimuth of the shot point, and/or, the coordinate based on geophone station in spherical coordinate system calculates the geophone station
Azimuth;Finally, fault strike is determined based on the changing rule between azimuth and edge diffraction coefficient.
In embodiments of the present invention, due in the edge diffraction coefficient value that is calculated in conventional cylindrical coordinate only comprising outgoing
Angle, and be transformed into spherical coordinate system, angle information can be increased to two angle informations of the angle of emergence and incidence angle by the angle of emergence, and
And respective party parallactic angle can be tried to achieve by the coordinate of shot point in spherical coordinate system and geophone station, and then by gained edge diffraction coefficient and orientation
Angle is mapped, and gained edge diffraction coefficient is transformed into three-dimensional perspective domain, thus, it becomes possible to taken by using diffracted wave
The information of band accurately analyzes fault strike, so alleviate by existing sectional analysis technology carry out sectional analysis when, analysis
The poor technical problem of the degree of accuracy.
Other features and advantages of the present invention will illustrate in the following description, also, partly become from specification
Obtain it is clear that or being understood by implementing the present invention.The purpose of the present invention and other advantages are in specification, claims
And specifically noted structure is realized and obtained in accompanying drawing.
To enable the above objects, features and advantages of the present invention to become apparent, preferred embodiment cited below particularly, and coordinate
Appended accompanying drawing, is described in detail below.
Brief description of the drawings
, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical scheme of the prior art
The required accompanying drawing used is briefly described in embodiment or description of the prior art, it should be apparent that, in describing below
Accompanying drawing is some embodiments of the present invention, for those of ordinary skill in the art, before creative work is not paid
Put, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is a kind of flow chart of fault strike recognition methods according to embodiments of the present invention;
Fig. 2 is the flow chart of step S104 in a kind of fault strike recognition methods according to embodiments of the present invention;
Fig. 3 is the flow chart of step S108 in a kind of fault strike recognition methods according to embodiments of the present invention;
Fig. 4 is a kind of schematic diagram of fault strike identification device according to embodiments of the present invention.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with accompanying drawing to the present invention
Technical scheme be clearly and completely described, it is clear that described embodiment is part of the embodiment of the present invention, rather than
Whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art are not making creative work premise
Lower obtained every other embodiment, belongs to the scope of protection of the invention.
Embodiment one:
According to embodiments of the present invention, there is provided a kind of embodiment of fault strike recognition methods is, it is necessary to illustrate, attached
The step of flow of figure illustrates can perform in the computer system of such as one group computer executable instructions, though also,
So logical order is shown in flow charts, but in some cases, can be with different from shown by order execution herein
Or the step of description.
Fig. 1 is a kind of flow chart of fault strike recognition methods according to embodiments of the present invention, as shown in figure 1, this method
Comprise the following steps:
Step S102, obtains the edge diffraction coefficient in cylindrical coordinate, and the edge diffraction coefficient produces for tomography edge
Edge diffraction coefficient;
Step S104, the edge diffraction coefficient in the cylindrical coordinate is changed, is converted to the side in spherical coordinate system
Edge diffraction coefficient, the edge diffraction coefficient in the spherical coordinate system include multiple angle informations;
In embodiments of the present invention, due in the edge diffraction coefficient value that is calculated in conventional cylindrical coordinate only comprising outgoing
Angle, and be transformed into spherical coordinate system, angle information can be increased to two angle informations of the angle of emergence and incidence angle by the angle of emergence, and
And respective party parallactic angle can be tried to achieve by the coordinate of shot point in spherical coordinate system and geophone station, and then by gained edge diffraction coefficient and orientation
Angle is mapped, and gained edge diffraction coefficient is transformed into three-dimensional perspective domain, thus, it becomes possible to taken by using diffracted wave
The information of band accurately analyzes fault strike, so alleviate by existing sectional analysis technology carry out fault strike analysis when,
The poor technical problem of recognition accuracy.
Step S106, the coordinate based on shot point in the spherical coordinate system calculate the azimuth of the shot point, and/or, base
The azimuth of the geophone station is calculated in coordinate of the geophone station in the spherical coordinate system;
Step S108, fault strike is determined based on the changing rule between the azimuth and the edge diffraction coefficient.
It should be noted that the edge diffraction coefficient described in step S108 refer to edge in above-mentioned cylindrical coordinate around
The edge diffraction coefficient penetrated in coefficient and spherical coordinate system.
In three-dimensional homogeneous medium space, if in edge diffraction caused by the edge in tomography or crack, then with
Azimuthal difference, the feature shown in the section at different directions angle have very big difference.For example, in diffraction edge, around
Polarity inversion, and the maximum in Diffraction Point tip edge diffraction coefficient value reaches three dimensions can occur for ejected wave.Work as section
In when diffracted wave occurs polarity inversion and diffraction amplitude and reaches maximum, the angle and diffraction edge-perpendicular, therefore can be whereby
Tomography and crack edge trend are judged, for identifying that little structure geologic body has significant application value in seismic prospecting.Therefore, grind
Study carefully edge diffraction coefficient and fault strike relation in omnibearing angel domain, it is small using diffracted wave identification in being explored for actual seismic
The tool that moves towards of tectonic geology body is of great significance.
The embodiment of the present invention is described in detail below in conjunction with embodiment.
In embodiments of the present invention, edge diffraction ripple caused by tomography edge when single shot point excites is obtained first;Then,
Calculate the value V of edge diffraction coefficient in cylindrical coordinateT(t) the edge diffraction coefficient, namely in cylindrical coordinate.
Specifically, edge diffraction coefficient VT(t) expression formula in cylindrical coordinate is equation below:
Wherein,And in above-mentioned formula, t be around
When walking of ejected wave, c are constant velocities, and shot point coordinate in cylindrical coordinate is (ρ0,φ0,z0), geophone station coordinate in cylindrical coordinate
For (ρ, φ, z), R is the distance of focal point and geophone station, and R' is that (image source point position shot point is on horizontal plane minute surface for image source point
Symmetrical point) with the distance of geophone station, R1For diffracted wave ray path (ray path is i.e. by shot point, by Diffraction Point,
Return to the shortest path of geophone station).
In embodiments of the present invention, R, R' and R1Expression formula can be to be described as follows:
In embodiments of the present invention, after the edge diffraction coefficient in acquiring cylindrical coordinate, it is possible to which post is sat
Edge diffraction coefficient in mark system is changed, and is converted to the edge diffraction coefficient in spherical coordinate system.
In an optional embodiment, as shown in Fig. 2 being changed to first edge diffraction coefficient, ball is converted to
Second edge diffraction coefficient in coordinate system comprises the following steps:
Step S1041, the first coordinate information and the second coordinate information are obtained, wherein, first coordinate information is described
Coordinate information of the shot point in the spherical coordinate system, second coordinate information are the geophone station in the spherical coordinate system
Coordinate information, the coordinate information include:The angle of emergence of diffracted wave, the incidence angle of diffracted wave, the shot point or the inspection
The distance between wave point and the origin of coordinates;
Step S1042, based on first coordinate information and second coordinate information to the side in the cylindrical coordinate
Edge diffraction coefficient is changed, and is converted to the edge diffraction coefficient in spherical coordinate system.
In embodiments of the present invention, the first coordinate information is coordinate information of the shot point in spherical coordinate system, and the second coordinate is believed
Cease the coordinate information for being geophone station in spherical coordinate system.First coordinate information can be expressed as:(rS,θS,φS), the second coordinate letter
Breath can be expressed as:(rR,θR,φR)。
First coordinate information (rS,θS,φS) in three variables be followed successively by:Between the shot point and the origin of coordinates away from
From the angle of emergence of diffracted wave, the incidence angle of diffracted wave, diffracted wave now is the diffracted wave detected in shot point;
Second coordinate information (rR,θR,φR) in three variables be followed successively by:Between the geophone station and the origin of coordinates away from
From the angle of emergence of diffracted wave, the incidence angle of diffracted wave, diffracted wave now is the diffracted wave detected in geophone station.
After the first coordinate information and the second coordinate information got, it is possible to based on the first coordinate information and second
Coordinate information is changed to the edge diffraction coefficient in cylindrical coordinate, obtains the edge diffraction coefficient in spherical coordinate system.
By foregoing description, due to only including an angle information in cylindrical coordinate, therefore by cylindrical coordinate around
Ejected wave is transformed under spherical coordinate system, and now, spherical coordinate system interrupting layer diffracted wave includes two angle informations (angle of emergence and incidences
Angle).The azimuth of shot point or geophone station is calculated again, now, it becomes possible to which edge diffraction coefficient is transformed into three-dimensional perspective
In domain, this is significant for analysis fault strike and azimuthal relation.
In another optional embodiment, based on first coordinate information and second coordinate information to described
Edge diffraction coefficient in cylindrical coordinate is changed, and is converted to the process description of edge diffraction coefficient in spherical coordinate system such as
Under:
The edge diffraction coefficient in cylindrical coordinate is changed according to formula, is converted to the edge diffraction in spherical coordinate system
Coefficient, wherein, expression formula of the edge diffraction coefficient in spherical coordinate system is equation below:
Wherein, RL=rS+rR, (rS,θS,φS) it is described
First coordinate information, (rR,θR,φR) it is second coordinate information, RSRRepresent the shot point with the geophone interval from RIR
Represent the distance of image source point and geophone station, RLPassed through for the ray road of diffracted wave, t is diffracted wave when walking, and c is constant velocity.
It should be noted that the edge diffraction coefficient V being as located under spherical coordinate system described by the formulaT(t)。
In embodiments of the present invention, after the edge diffraction coefficient in cylindrical coordinate is converted under spherical coordinate system, just
The azimuth of the shot point can be calculated based on coordinate of the shot point in the spherical coordinate system, and/or, seat is being asked based on geophone station
Coordinate in mark system calculates the azimuth of geophone station.
In another optional embodiment, the coordinate based on shot point in the spherical coordinate system calculates the shot point
Azimuthal detailed process is described as follows:
Pass through formulaThe azimuth of the shot point is calculated, wherein,For institute
State the azimuth of shot point, rSFor the distance between shot point and the origin of coordinates, θSFor the incidence angle of the diffracted wave, φSFor diffracted wave
The angle of emergence, wherein, θSAnd φSAngular range be 0 to 180 °,Angular range be 0 to 360 °.
In embodiments of the present invention, when calculating the azimuth of shot point, the main coordinate by shot point under spherical coordinate system
Information (rS,θS,φS) (that is, first coordinate information) calculate the azimuth of shot point.Specifically, formula can be passed throughTo calculate the azimuth of shot pointIt is up due to only being received in seismic prospecting
Ripple, therefore the angle of emergence and incidence angle θS、φSAngular range be 0 to 180 °, and due to arrangement be stereo observing system, therefore shot point
Project to azimuth coverage in the planes be 0 to 360 °.
By foregoing description, in embodiments of the present invention, then due to only including an angle information in cylindrical coordinate,
Therefore the diffracted wave in cylindrical coordinate is transformed under spherical coordinate system, now, spherical coordinate system interrupting layer edge diffraction coefficient includes two
Individual angle information (angle of emergence and incidence angle).The azimuth of shot point or geophone station is calculated again, now, it becomes possible to by edge
Diffraction coefficient is transformed into three-dimensional perspective domain, and this is significant for analysis fault strike and azimuthal relation.
In another optional embodiment, the coordinate based on geophone station in the spherical coordinate system calculates the detection
Azimuthal detailed process of point is described as follows:
Pass through formulaThe azimuth of the geophone station is calculated, wherein,For
The azimuth of the geophone station, rRFor the distance between the geophone station and the origin of coordinates, θRFor the incidence angle of the diffracted wave,
φRFor the angle of emergence of diffracted wave, wherein, θRAnd φRAngular range be 0 to 180 °,Angular range be 0 to 360 °.
In embodiments of the present invention, when calculating the azimuth of geophone station, mainly by geophone station under spherical coordinate system
Coordinate information (rR,θR,φR) (that is, second coordinate information) calculate the azimuth of geophone station.Specifically, formula can be passed throughTo calculate the azimuth of geophone stationIt is up due to only being received in seismic prospecting
Ripple, therefore the angle of emergence and incidence angle θR、φRAngular range be 0 to 180 °, and due to arrangement be stereo observing system, therefore shot point
Project to azimuth coverage in the planes be 0 to 360 °.
By foregoing description, in embodiments of the present invention, then by the edge diffraction coefficient in cylindrical coordinate change to
After under spherical coordinate system, it becomes possible to obtain two angle informations, i.e. the angle of emergence and incidence angle.Shot point or inspection are calculated again
The azimuth of wave point, now, it becomes possible to edge diffraction coefficient is transformed into three-dimensional perspective domain, this for analysis fault strike with
Azimuthal relation is significant.
In another optional embodiment, as shown in figure 3, based on the azimuth and the edge diffraction coefficient it
Between changing rule determine that fault strike comprises the following steps:
Step S1081, the one-to-one relationship between the azimuth and the edge diffraction coefficient is determined,
Wherein, the corresponding relation is:When positioned at the surface of the tomography, if incidence angle is φSFor 90 °, then
The edge diffraction coefficient VT(t) maximum is reached;Perpendicular to the fault strike and passing through the edge diffraction of the tomography
On the section of coefficient, the phenomenon of edge diffraction coefficient polarity inversion can occur directly over Diffraction Point;Work as output angle θR=90 °,
Incidence angle φSAt=90 °, parallel to the edge diffraction coefficient V of the fault strikeT(t) maximum in three dimensions is reached
Value;
Step S1082, graph of a relation is built based on the corresponding relation, wherein, the graph of a relation includes edge diffraction system
Changing rule between several and described azimuth;
Step S1083, the fault strike is determined based on the graph of a relation.
In embodiments of the present invention, after the azimuth of shot point or geophone station is obtained, it is possible to determine azimuth and
Corresponding relation between edge diffraction coefficient, and based on relation structure graph of a relation, so that it is determined that fault strike.
Wherein it is determined that the corresponding relation between the azimuth and the edge diffraction coefficient includes:
Determine the first corresponding relation, wherein, first corresponding relation for the shot point azimuth and the edge around
The corresponding relation penetrated between coefficient;And/or
The second corresponding relation is determined, wherein, azimuth and the edge of second corresponding relation for the geophone station
Corresponding relation between diffraction coefficient.
That is, in embodiments of the present invention, it is determined that after the first corresponding relation and the second corresponding relation, it is possible to
The graph of a relation is built based on first corresponding relation and second corresponding relation.
For example, building the first graph of a relation based on the first corresponding relation, then, the second relation is built based on the second corresponding relation
Figure.It is then possible to by the first graph of a relation to determine fault strike, and fault strike determined by the second graph of a relation.Also
The first graph of a relation and the second graph of a relation can be combined to integrate determination fault strike.
Specifically, in embodiments of the present invention, edge diffraction coefficient and incidence angle, the changing rule of the angle of emergence are analyzed, i.e.,:
(parallel to fault strike) when positioned at the surface of the tomography, incidence angle φS=90 °, edge diffraction coefficient VT(t) reach
To maximum;On the section perpendicular to the fault strike and by the edge diffraction coefficient of the tomography, in Diffraction Point just
The phenomenon of edge diffraction coefficient polarity inversion can occur for top;Work as output angle θR=90 °, incidence angle φSAt=90 °, parallel to disconnected
The edge diffraction coefficient profile values V of layer trendT(t) maximum in three dimensions is reached.
In embodiments of the present invention, when building graph of a relation, set along x-axis positive direction as 0 ° of azimuth, fault strike and Y
Axle is parallel.According to gained azimuthWith edge diffraction coefficient VT(t) graph of a relation is drawn, analysis edge diffraction coefficient is with azimuth
Changing rule:During azimuth progressively increases to 360 ° by 0 °, the changing pattern of edge diffraction coefficient be by strong (0 °) →
Weak → strong (90 °) → weak → strong (180 °) → weak → strong (270 °) → weak → strong (360 °);Parallel to fault strike edge around
Penetrate coefficient section and reach maximum (azimuth is 90 ° or 270 °).
In embodiments of the present invention, after fault strike is obtained, fault strike can also be verified.Authenticated
Cheng Zhong, it can carry out the mode of three-dimensional all-directional angle domain skew by taking out angle gathers in angle domain, it is deep to obtain actual tomography
Degree, the relation between the angle of emergence, azimuthAnd analysis obtains tomography edge diffraction coefficient and the angle of emergence, orientation
The relation at angle.Between the edge diffraction coefficient being calculated in conjunction with the algorithm provided using the embodiment of the present invention and azimuth
Changing rule come judge fault strike predict whether it is correct.
To sum up, the present invention is directed to diffracted wave dynamic characteristic in seismic wave field and proposes a kind of fault strike recognition methods.
The present invention calculates edge diffraction coefficient caused by tomography under cylindrical coordinate, due to post on the basis of diffracted wave dynamic characteristic
An angle information only is included in coordinate, therefore diffracted wave in cylindrical coordinate is transformed under spherical coordinate system.Spherical coordinate system interrupting layer
Edge diffraction coefficient includes two angle informations, can be calculated according to the coordinate of the two angle informations and shot point (geophone station) corresponding
Azimuth, using the azimuth relation corresponding with edge diffraction coefficient value of each shot point (geophone station) in stereo observing system,
Graph of a relation of the edge diffraction coefficient with azimuthal variation is drawn, reaches maximum and side so as to which basis works as edge diffraction amplitude in section
The angle recognition fault strike of polarity inversion occurs for edge diffraction coefficient.
Embodiment two:
The embodiment of the present invention additionally provides a kind of fault strike identification device, and the fault strike identification device is mainly used in holding
The fault strike recognition methods that row the above of the embodiment of the present invention is provided, is walked to tomography provided in an embodiment of the present invention below
Specific introduction is done to identification device.
Fig. 4 is a kind of schematic diagram of fault strike identification device according to embodiments of the present invention, as shown in figure 4, the tomography
Moving towards identification device mainly includes:Acquiring unit 10, converting unit 20, computing unit 30 and determining unit 40, wherein:
Acquiring unit 10, for obtaining the edge diffraction coefficient in cylindrical coordinate, the edge diffraction coefficient is tomography side
Edge diffraction coefficient caused by edge;
Converting unit 20, for being changed to the edge diffraction coefficient in the cylindrical coordinate, be converted to spherical coordinate system
In edge diffraction coefficient, the edge diffraction coefficient in the spherical coordinate system includes multiple angle informations, wherein, the post is sat
Edge diffraction coefficient in mark system includes outgoing angle information, and the multiple angle information includes outgoing angle information and incidence angle letter
Breath;
Computing unit 30, for the azimuth of the coordinate calculating shot point based on shot point in the spherical coordinate system, and/
Or, the coordinate based on geophone station in the spherical coordinate system calculates the azimuth of the geophone station;
Determining unit 40, for determining tomography based on the changing rule between the azimuth and the edge diffraction coefficient
Trend.
In embodiments of the present invention, due in the edge diffraction coefficient value that is calculated in conventional cylindrical coordinate only comprising outgoing
Angle, and be transformed into spherical coordinate system, angle information can be increased to two angle informations of the angle of emergence and incidence angle by the angle of emergence, and
And respective party parallactic angle can be tried to achieve by the coordinate of shot point in spherical coordinate system and geophone station, and then by gained edge diffraction coefficient and orientation
Angle is mapped, and gained edge diffraction coefficient is transformed into three-dimensional perspective domain, thus, it becomes possible to by using edge diffraction institute
The information of carrying accurately analyzes fault strike, and then alleviates and carry out fault strike analysis by existing sectional analysis technology
When, the poor technical problem of recognition accuracy.
Alternatively, converting unit 20 includes:Acquisition module, for obtaining the first coordinate information and the second coordinate information, its
In, first coordinate information is coordinate information of the shot point in the spherical coordinate system, and second coordinate information is institute
Coordinate information of the geophone station in the spherical coordinate system is stated, the coordinate information includes:The angle of emergence of diffracted wave, diffracted wave
Incidence angle, the shot point or the distance between the geophone station and the origin of coordinates;Modular converter, for being sat based on described first
Mark information and second coordinate information are changed to the edge diffraction coefficient in the cylindrical coordinate, are converted to spherical coordinate system
In edge diffraction coefficient.
Alternatively, modular converter is used for:The edge diffraction coefficient in the cylindrical coordinate is changed according to formula, turned
The edge diffraction coefficient being changed in spherical coordinate system, wherein, the formula is described as follows:Wherein, RL=rS+rR, (rS,θS,φS) it is described
First coordinate information, (rR,θR,φR) it is second coordinate information, RSRRepresent the shot point with the geophone interval from RIR
Represent the distance of image source point and geophone station, RLPassed through for the ray road of diffracted wave.
Alternatively, computing unit 30 is used for:Pass through formulaCalculate the big gun
The azimuth of point, wherein,For the azimuth of the shot point, rSFor the distance between shot point and the origin of coordinates, θSFor it is described around
The incidence angle of ejected wave, φSFor the angle of emergence of diffracted wave, wherein, θSAnd φSAngular range be 0 to 180 °,Angular range
For 0 to 360 °.
Alternatively, computing unit 30 is used for:Pass through formulaCalculate the inspection
The azimuth of wave point, wherein,For the azimuth of the geophone station, rRFor the distance between the geophone station and the origin of coordinates,
θRFor the incidence angle of the diffracted wave, φRFor the angle of emergence of diffracted wave, wherein, θRAnd φRAngular range be 0 to 180 °,
Angular range be 0 to 360 °.
Optionally it is determined that unit 40 is used for:The corresponding relation between the azimuth and the edge diffraction coefficient is determined,
Wherein, the corresponding relation is:When positioned at the surface of the tomography, if incidence angle is φSFor 90 °, then the edge
Diffraction coefficient VT(t) maximum is reached;Perpendicular to the fault strike and the edge diffraction coefficient that passes through the tomography cuts open
On face, the phenomenon of edge diffraction coefficient polarity inversion can occur directly over Diffraction Point;Work as output angle θR=90 °, incidence angle φS
At=90 °, parallel to the edge diffraction coefficient V of the fault strikeT(t) maximum in three dimensions is reached;It is based on
The corresponding relation builds graph of a relation, wherein, the graph of a relation includes the change between edge diffraction coefficient and the azimuth
Law;The fault strike is determined based on the graph of a relation.
Optionally it is determined that unit 40 is additionally operable to:The first corresponding relation is determined, wherein, first corresponding relation is described
Corresponding relation between the azimuth of shot point and the edge diffraction coefficient;And/or the second corresponding relation is determined, wherein, it is described
Corresponding relation of second corresponding relation between the azimuth of the geophone station and the edge diffraction coefficient.
Optionally it is determined that unit 40 is additionally operable to:Based on first corresponding relation and second corresponding relation structure institute
State graph of a relation.
The device that the embodiment of the present invention is provided, its realization principle and caused technique effect and preceding method embodiment phase
Together, to briefly describe, device embodiment part does not refer to part, refers to corresponding contents in preceding method embodiment.
In addition, in the description of the embodiment of the present invention, unless otherwise clearly defined and limited, term " installation ", " phase
Even ", " connection " should be interpreted broadly, for example, it may be being fixedly connected or being detachably connected, or be integrally connected;Can
To be mechanical connection or electrical connection;Can be joined directly together, can also be indirectly connected by intermediary, Ke Yishi
The connection of two element internals.For the ordinary skill in the art, with concrete condition above-mentioned term can be understood at this
Concrete meaning in invention.
In the description of the invention, it is necessary to explanation, term " " center ", " on ", " under ", "left", "right", " vertical ",
The orientation or position relationship of the instruction such as " level ", " interior ", " outer " be based on orientation shown in the drawings or position relationship, merely to
Be easy to the description present invention and simplify description, rather than instruction or imply signified device or element must have specific orientation,
With specific azimuth configuration and operation, therefore it is not considered as limiting the invention.In addition, term " first ", " second ",
" the 3rd " is only used for describing purpose, and it is not intended that instruction or hint relative importance.
A kind of fault strike recognition methods and the computer program product of device that the embodiment of the present invention is provided, including deposit
The computer-readable recording medium of the executable non-volatile program code of processor, the instruction that described program code includes are stored up
Available for the method performed described in previous methods embodiment, specific implementation can be found in embodiment of the method, will not be repeated here.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description,
The specific work process of device and unit, the corresponding process in preceding method embodiment is may be referred to, will not be repeated here.
In several embodiments provided herein, it should be understood that disclosed systems, devices and methods, can be with
Realize by another way.Device embodiment described above is only schematical, for example, the division of the unit,
Only a kind of division of logic function, can there is other dividing mode when actually realizing, in another example, multiple units or component can
To combine or be desirably integrated into another system, or some features can be ignored, or not perform.It is another, it is shown or beg for
The mutual coupling of opinion or direct-coupling or communication connection can be by some communication interfaces, device or unit it is indirect
Coupling or communication connection, can be electrical, mechanical or other forms.
The unit illustrated as separating component can be or may not be physically separate, show as unit
The part shown can be or may not be physical location, you can with positioned at a place, or can also be distributed to multiple
On NE.Some or all of unit therein can be selected to realize the mesh of this embodiment scheme according to the actual needs
's.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing unit, can also
That unit is individually physically present, can also two or more units it is integrated in a unit.
If the function is realized in the form of SFU software functional unit and is used as independent production marketing or in use, can be with
It is stored in the executable non-volatile computer read/write memory medium of a processor.Based on such understanding, the present invention
The part that is substantially contributed in other words to prior art of technical scheme or the part of the technical scheme can be with software
The form of product is embodied, and the computer software product is stored in a storage medium, including some instructions are causing
One computer equipment (can be personal computer, server, or network equipment etc.) performs each embodiment institute of the present invention
State all or part of step of method.And foregoing storage medium includes:USB flash disk, mobile hard disk, read-only storage (ROM, Read-
Only Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD etc. are various can be with
The medium of store program codes.
Finally it should be noted that:Embodiment described above, it is only the embodiment of the present invention, to illustrate the present invention
Technical scheme, rather than its limitations, protection scope of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, it will be understood by those within the art that:Any one skilled in the art
The invention discloses technical scope in, it can still modify to the technical scheme described in previous embodiment or can be light
Change is readily conceivable that, or equivalent substitution is carried out to which part technical characteristic;And these modifications, change or replacement, do not make
The essence of appropriate technical solution departs from the spirit and scope of technical scheme of the embodiment of the present invention, should all cover the protection in the present invention
Within the scope of.Therefore, protection scope of the present invention described should be defined by scope of the claims.
Claims (10)
- A kind of 1. fault strike recognition methods, it is characterised in that including:The edge diffraction coefficient in cylindrical coordinate is obtained, the edge diffraction coefficient is edge diffraction system caused by tomography edge Number;Edge diffraction coefficient in the cylindrical coordinate is changed, is converted to the edge diffraction coefficient in spherical coordinate system, institute The edge diffraction coefficient stated in spherical coordinate system includes multiple angle informations, wherein, the edge diffraction system in the cylindrical coordinate Number includes outgoing angle information, and the multiple angle information includes outgoing angle information and incidence angle information;Coordinate based on shot point in the spherical coordinate system calculates the azimuth of the shot point, and/or, based on geophone station described Coordinate in spherical coordinate system calculates the azimuth of the geophone station;Fault strike is determined based on the changing rule between the azimuth and the edge diffraction coefficient.
- 2. according to the method for claim 1, it is characterised in that the edge diffraction coefficient in the cylindrical coordinate is turned Change, the edge diffraction coefficient be converted in spherical coordinate system includes:The first coordinate information and the second coordinate information are obtained, wherein, first coordinate information is that the shot point is sat in the ball Coordinate information in mark system, second coordinate information is coordinate information of the geophone station in the spherical coordinate system, described Coordinate information includes:The angle of emergence of diffracted wave, the incidence angle of diffracted wave, the shot point or the geophone station and the origin of coordinates The distance between;The edge diffraction coefficient in the cylindrical coordinate is carried out based on first coordinate information and second coordinate information Conversion, is converted to the edge diffraction coefficient in spherical coordinate system.
- 3. according to the method for claim 2, it is characterised in that believed based on first coordinate information and second coordinate Breath is changed to the edge diffraction coefficient in the cylindrical coordinate, and the edge diffraction coefficient be converted in spherical coordinate system includes:The edge diffraction coefficient in the cylindrical coordinate is changed according to formula, is converted to the edge diffraction in spherical coordinate system Coefficient, wherein, the formula is described as follows:<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>V</mi> <mi>T</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>H</mi> <mrow> <mo>(</mo> <mi>cos</mi> <mo>&lsqb;</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>S</mi> </msub> <mo>-</mo> <msub> <mi>&phi;</mi> <mi>R</mi> </msub> </mrow> <mn>2</mn> </mfrac> <mo>&rsqb;</mo> <mo>)</mo> </mrow> <mfrac> <mrow> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mfrac> <msub> <mi>R</mi> <mrow> <mi>S</mi> <mi>R</mi> </mrow> </msub> <mi>C</mi> </mfrac> <mo>)</mo> </mrow> </mrow> <msub> <mi>R</mi> <mrow> <mi>S</mi> <mi>R</mi> </mrow> </msub> </mfrac> <mo>&PlusMinus;</mo> <mi>H</mi> <mrow> <mo>(</mo> <mi>cos</mi> <mo>&lsqb;</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>S</mi> </msub> <mo>+</mo> <msub> <mi>&phi;</mi> <mi>R</mi> </msub> </mrow> <mn>2</mn> </mfrac> <mo>&rsqb;</mo> <mo>)</mo> </mrow> <mfrac> <mrow> <mi>&delta;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mfrac> <msub> <mi>R</mi> <mrow> <mi>I</mi> <mi>R</mi> </mrow> </msub> <mi>C</mi> </mfrac> <mo>)</mo> </mrow> </mrow> <msub> <mi>R</mi> <mrow> <mi>I</mi> <mi>R</mi> </mrow> </msub> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mfrac> <mrow> <mi>H</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mfrac> <msub> <mi>R</mi> <mi>L</mi> </msub> <mi>C</mi> </mfrac> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&pi;</mi> <msqrt> <mrow> <msup> <mi>t</mi> <mn>2</mn> </msup> <mo>-</mo> <mfrac> <mrow> <msup> <msub> <mi>R</mi> <mi>L</mi> </msub> <mn>2</mn> </msup> </mrow> <msup> <mi>C</mi> <mn>2</mn> </msup> </mfrac> </mrow> </msqrt> </mrow> </mfrac> <mo>{</mo> <mfrac> <mrow> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>cos</mi> <mo>&lsqb;</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>S</mi> </msub> <mo>-</mo> <msub> <mi>&phi;</mi> <mi>R</mi> </msub> </mrow> <mn>2</mn> </mfrac> <mo>&rsqb;</mo> <mo>)</mo> </mrow> <msqrt> <mrow> <msubsup> <mi>R</mi> <mi>L</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>S</mi> <mi>R</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> <mrow> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mi>t</mi> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>I</mi> <mi>R</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>&PlusMinus;</mo> <mfrac> <mrow> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>cos</mi> <mo>&lsqb;</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>S</mi> </msub> <mo>+</mo> <msub> <mi>&phi;</mi> <mi>R</mi> </msub> </mrow> <mn>2</mn> </mfrac> <mo>&rsqb;</mo> <mo>)</mo> </mrow> <msqrt> <mrow> <msubsup> <mi>R</mi> <mi>L</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>I</mi> <mi>R</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> <mrow> <msup> <mi>c</mi> <mn>2</mn> </msup> <msup> <mi>t</mi> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>I</mi> <mi>R</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> <mo>,</mo> </mrow>Wherein, RL=rS+rR, (rS,θS,φS) it is described First coordinate information, (rR,θR,φR) it is second coordinate information, RSRRepresent the shot point with the geophone interval from RIR Represent the distance of image source point and geophone station, RLPassed through for the ray road of diffracted wave.
- 4. according to the method for claim 1, it is characterised in that institute is calculated based on coordinate of the shot point in the spherical coordinate system Stating the azimuth of shot point includes:Pass through formulaThe azimuth of the shot point is calculated, wherein,For the big gun The azimuth of point, rSFor the distance between shot point and the origin of coordinates, θSFor the incidence angle of the diffracted wave, φSFor going out for diffracted wave Firing angle, wherein, θSAnd φSAngular range be 0 to 180 °,Angular range be 0 to 360 °.
- 5. the method according to claim 1 or 4, it is characterised in that based on coordinate of the geophone station in the spherical coordinate system Calculating the azimuth of the geophone station includes:Pass through formulaTo calculate the azimuth of the geophone stationWherein,For The azimuth of the geophone station, rRFor the distance between the geophone station and the origin of coordinates, θRFor the incidence angle of the diffracted wave, φRFor the angle of emergence of diffracted wave, wherein, θRAnd φRAngular range be 0 to 180 °,Angular range be 0 to 360 °.
- 6. according to the method for claim 1, it is characterised in that based between the azimuth and the edge diffraction coefficient Changing rule determine that fault strike includes:The corresponding relation between the azimuth and the edge diffraction coefficient is determined,Wherein, the corresponding relation is:When positioned at the surface of the tomography, if incidence angle is φSFor 90 °, then the side Edge diffraction coefficient VT(t) maximum is reached;Perpendicular to the fault strike and the edge diffraction coefficient that passes through the tomography On section, the phenomenon of edge diffraction coefficient polarity inversion can occur directly over Diffraction Point;Work as output angle θR=90 °, incidence angle φSAt=90 °, parallel to the edge diffraction coefficient V of the fault strikeT(t) maximum in three dimensions is reached;Graph of a relation is built based on the corresponding relation, wherein, the graph of a relation is included between diffracted wave and the azimuth Changing rule;The fault strike is determined based on the graph of a relation.
- 7. according to the method for claim 6, it is characterised in that determine between the azimuth and the edge diffraction coefficient Corresponding relation include:The first corresponding relation is determined, wherein, azimuth and the edge diffraction system of first corresponding relation for the shot point Corresponding relation between number;And/orThe second corresponding relation is determined, wherein, azimuth and the edge diffraction of second corresponding relation for the geophone station Corresponding relation between coefficient.
- 8. according to the method for claim 7, it is characterised in that included based on corresponding relation structure graph of a relation:The graph of a relation is built based on first corresponding relation and second corresponding relation.
- A kind of 9. fault strike identification device, it is characterised in that including:Acquiring unit, for obtaining the edge diffraction coefficient in cylindrical coordinate, the edge diffraction coefficient produces for tomography edge Edge diffraction coefficient;Converting unit, for being changed to the edge diffraction coefficient in the cylindrical coordinate, be converted to the side in spherical coordinate system Edge diffraction coefficient, the edge diffraction coefficient in the spherical coordinate system include multiple angle informations, wherein, in the cylindrical coordinate Edge diffraction coefficient include outgoing angle information, the multiple angle information includes outgoing angle information and incidence angle information;Computing unit, the azimuth of the shot point is calculated for the coordinate based on shot point in the spherical coordinate system, and/or, base The azimuth of the geophone station is calculated in coordinate of the geophone station in the spherical coordinate system;Determining unit, for determining fault strike based on the changing rule between the azimuth and the edge diffraction coefficient.
- 10. a kind of computer-readable medium for the non-volatile program code that can perform with processor, it is characterised in that described Program code makes any described method in the computing device the claims 1-8.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711089298.8A CN107861154B (en) | 2017-11-08 | 2017-11-08 | Fault strike recognition methods and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711089298.8A CN107861154B (en) | 2017-11-08 | 2017-11-08 | Fault strike recognition methods and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107861154A true CN107861154A (en) | 2018-03-30 |
CN107861154B CN107861154B (en) | 2018-11-20 |
Family
ID=61701346
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711089298.8A Active CN107861154B (en) | 2017-11-08 | 2017-11-08 | Fault strike recognition methods and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107861154B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108845358A (en) * | 2018-09-26 | 2018-11-20 | 中国矿业大学(北京) | Tomography and the recognition methods of structural anomaly body and device |
CN112255679A (en) * | 2020-10-26 | 2021-01-22 | 中国石油天然气集团有限公司 | Seismic data diffraction depth migration processing method and device |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050232077A1 (en) * | 2003-08-19 | 2005-10-20 | Naum Marmalyevskyy | Method, system and apparatus for interpreting seismic data using duplex waves |
CN104237940A (en) * | 2014-09-29 | 2014-12-24 | 中国石油天然气股份有限公司 | Diffracted wave imaging method and device based on dynamic characteristics |
CN105259577A (en) * | 2015-10-13 | 2016-01-20 | 中国石油天然气集团公司 | Method and device for determining angle information of formation boundary |
CN106772583A (en) * | 2017-01-10 | 2017-05-31 | 中国科学院地质与地球物理研究所 | A kind of earthquake diffracted wave separation method and device |
-
2017
- 2017-11-08 CN CN201711089298.8A patent/CN107861154B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050232077A1 (en) * | 2003-08-19 | 2005-10-20 | Naum Marmalyevskyy | Method, system and apparatus for interpreting seismic data using duplex waves |
CN104237940A (en) * | 2014-09-29 | 2014-12-24 | 中国石油天然气股份有限公司 | Diffracted wave imaging method and device based on dynamic characteristics |
CN105259577A (en) * | 2015-10-13 | 2016-01-20 | 中国石油天然气集团公司 | Method and device for determining angle information of formation boundary |
CN106772583A (en) * | 2017-01-10 | 2017-05-31 | 中国科学院地质与地球物理研究所 | A kind of earthquake diffracted wave separation method and device |
Non-Patent Citations (2)
Title |
---|
孙建国: "《高频渐进散射理论及其在地球物理场数值模拟与反演成像中的应用》", 《吉林大学学报(地球科学版)》 * |
赵惊涛等: "《地震勘探中的边缘绕射波及其动力学识别方法》", 《地球物理学进展》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108845358A (en) * | 2018-09-26 | 2018-11-20 | 中国矿业大学(北京) | Tomography and the recognition methods of structural anomaly body and device |
CN112255679A (en) * | 2020-10-26 | 2021-01-22 | 中国石油天然气集团有限公司 | Seismic data diffraction depth migration processing method and device |
CN112255679B (en) * | 2020-10-26 | 2023-09-26 | 中国石油天然气集团有限公司 | Method and device for processing diffraction depth deviation of seismic data |
Also Published As
Publication number | Publication date |
---|---|
CN107861154B (en) | 2018-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106094029B (en) | Utilize the method for offset distance vector piece geological data Predicating Reservoir Fractures | |
CN105467440B (en) | Omnidirectional vector seismic data processing method and device | |
CN104237940B (en) | Diffracted wave imaging method and device based on dynamic characteristics | |
CN107817523B (en) | The analysis method and device of diffracted wave migration velocity | |
CN107783185B (en) | Chromatography static correction processing method and device | |
CN105093292A (en) | Data processing method and device for seismic imaging | |
CN102937721A (en) | Limited frequency tomography method for utilizing preliminary wave travel time | |
CN103995288A (en) | Gauss beam prestack depth migration method and device | |
CN105425299B (en) | Method and apparatus for determining formation fracture distribution | |
MX2015000034A (en) | Methods and systems for optimizing generation of seismic images. | |
CN101923175A (en) | Method for directly generating angle gathers by using wave-equation migration | |
CN103364833A (en) | High-precision dip estimation method | |
CN107656308B (en) | A kind of common scattering point pre-stack time migration imaging method based on time depth scanning | |
CN107015275B (en) | Karst collapse col umn detection method and device | |
CN104199088B (en) | Incident angle gather extraction method and system | |
CN107861154A (en) | Fault strike recognition methods and device | |
CN106199704A (en) | A kind of Three-dimendimal fusion submarine cable seismic data velocity modeling method | |
CN105353406B (en) | A kind of method and apparatus for generating angle gathers | |
CN105093300B (en) | Geologic body boundary identification method and device | |
CN104459787A (en) | Speed analysis method through seismic record of vertical receiving array | |
CN108693560A (en) | A kind of scattering wave imaging method and system based on cross-correlation road | |
CN106908838A (en) | The method for building target area stratum inelastic attenuation quality factor three-dimensional model | |
US9791580B2 (en) | Methods and systems to separate wavefields using pressure wavefield data | |
CN107422375B (en) | The determination method and apparatus of the CRP gather of subsea node | |
US20150293244A1 (en) | Estimation of a far field signature in a second direction from a far field signature in a first direction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |