CN105093300B - A kind of boundary recognition of geological body method and device - Google Patents
A kind of boundary recognition of geological body method and device Download PDFInfo
- Publication number
- CN105093300B CN105093300B CN201510445467.1A CN201510445467A CN105093300B CN 105093300 B CN105093300 B CN 105093300B CN 201510445467 A CN201510445467 A CN 201510445467A CN 105093300 B CN105093300 B CN 105093300B
- Authority
- CN
- China
- Prior art keywords
- value
- point
- mrow
- sampled point
- dip
- 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.)
- Active
Links
Abstract
The present invention provides a kind of boundary recognition of geological body method and device, and the boundary recognition of geological body method includes：Set a length L to be 2M+1 sliding window, and length L rank value is set for N, obtain 2M+1 sampled point, and Gauss function G (n) is calculated according to 2M+1 sampled point；The optimum angle of incidence pair at the local reflex interface where each sampled point is determined using dip scanning method；Respectively using each sampled point as analysis sampling point, and centered on analyzing sampling point, read the amplitude of 2M+1 point respectively from multiple directions according to the apparent dip of analysis sampling point, generate multiple seismic signal sequence Sk(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1；According to earthquake signal sequence SkAnd Gauss function G (n) calculates new seismic sequence G (n)Sk(n)；Calculate new seismic sequence GSK (n) generalized hilbert transform value GHk；According to generalized hilbert transform value GHkThe border detection value of generation analysis sampling point；A data volume is generated according to multiple border detection values, and horizon slice analysis is carried out to data volume, geological objects boundary is recognized.
Description
Technical field
The present invention is the interpretation technique on geophysical prospecting for oil geological data, specifically can on one kind
Recognize the boundary recognition of geological body method and device of small-scale geological objects boundary.
Background technology
At present, coherent technique is widely used in detection river course, fracture.Although coherent technique experienced three generations's development,
All there is average effect, it is difficult to highlight narrow river course and the fracture of small turn-off.Meanwhile, coherent algorithm is to noise-sensitive, in river course, fracture
Mutual cut place, seismic signal is mixed and disorderly, noise is strong, and coherent technique is difficult to accurately identify narrow river course and small fracture, can only identify
River course, the substantially development region of fracture.
The anti-noise ability of existing generalized hilbert transform bound test technology is strong, can highlight plastid side on a small scale
Boundary.In practice flow, generalized hilbert transform is carried out from main profile and cross-track both direction using window function
Detect different geological objects boundaries.However, often spatial distribution is complicated for the trend of actual geologic body (river course, fracture), only surveyed to main
Line and cross-track both direction, which carry out border detection, can not comprehensively detect geological objects boundary.
The content of the invention
The main purpose of the embodiment of the present invention is to provide a kind of boundary recognition of geological body method and device, complete more to like
Detect geological objects boundary to face.
To achieve these goals, the embodiment of the present invention provides a kind of boundary recognition of geological body method, described geologic body
Boundary Recognition method includes：Set a length L to be 2M+1 sliding window, and set length L rank value for N, obtain 2M+1
Sampled point, and Gauss function G (n) is calculated according to the 2M+1 sampled point, M is natural number；It is true using dip scanning method
The optimum angle of incidence pair at the local reflex interface where fixed each sampled point；It regard each sampled point as analysis sample respectively
Point, and centered on the analysis sampling point, 2M+1 are read respectively from multiple directions according to the apparent dip of the analysis sampling point
The amplitude of point, generates multiple seismic signal sequence Sk(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1；According to
The seismic signal sequence SkAnd Gauss function G (n) calculates new seismic sequence G (n)Sk(n)；Calculate the new seismic sequence
GSK (n) generalized hilbert transform value GHk；According to the generalized hilbert transform value GHkThe generation analysis sampling point
Border detection value；
A data volume is generated according to multiple described border detection values, and horizon slice analysis is carried out to the data volume,
Recognize geological objects boundary.
In one embodiment, above-mentioned utilization dip scanning method determines local reflex circle where each sampled point
The optimum angle of incidence pair in face, including：Setting one point, major and minor axis centered on sampled point (x, y) are respectively a, b horizontal window；
The seismologic parameter of the sampled point is obtained by the horizontal window；The seismologic parameter is substituted into sampled point coherent calculation public
Formula, obtains the coherent value of the sampled point；Line direction, the extreme value P_min, P_max, Q_ of road direction apparent formation dip are set respectively
Min, Q_max；Apparent dip P, Q of geologic body reflecting interface are determined by certain step delta P, Δ Q；It is fixed according to nyquist sampling
Rule obtains np*nqIndividual inclination angle to (P, Q) and corresponding coherent value c (t, P, Q), wherein, np=(P_max-P_min)/Δ P, nq=
(Q_max-Q_min)/ΔQ；With apparent dip to (PL,QM) centered on point choose c*d apparent dip (Pl, Qm) and its corresponding phase
Dry values c (t, Pl,Qm) a curved surface G (P, Q) is fitted to, wherein, Pl=PL± c* Δs P, Qm=QM±d*ΔP；The curved surface G
The corresponding inclination angle of maximum of (P, Q) is to the optimum angle of incidence pair of (P, Q) for the sampled point.
In one embodiment, above-mentioned Gauss function is：Wherein, α falling for standard deviation
Number.
In one embodiment, above-mentioned new seismic sequence GSK (n)=Sk(n) * G (n), k=1,2 ..., i.
In one embodiment, it is above-mentioned according to the generalized hilbert transform value GHkThe side of the generation analysis sampling point
Boundary's detected value, including：Calculate the generalized hilbert transform value GHkImaginary part absolute value IGHk；Judge multiple described absolute
Value IGHkIn maximum, and using the maximum be used as it is described analysis sampling point border detection value.
The embodiment of the present invention also provides a kind of boundary recognition of geological body device, described boundary recognition of geological body device bag
Include：Gauss function computing unit, for setting a length L to be 2M+1 sliding window, and sets length L rank value N, obtains
2M+1 sampled point, and Gauss function G (n) is calculated according to the 2M+1 sampled point, M is natural number；Optimum angle of incidence is to true
Order member, for the optimum angle of incidence pair using the local reflex interface where each sampled point of dip scanning device determination；
Seismic signal sequence generating unit, for respectively using each sampled point as analysis sampling point, and using it is described analysis sampling point as
Center, according to it is described analysis sampling point optimum angle of incidence to from multiple directions respectively read 2M+1 point amplitude, generation it is many
Individual seismic signal sequence Sk(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1；New seismic sequence computing unit, is used
According to the seismic signal sequence SkAnd Gauss function G (n) calculates new seismic sequence G (n)Sk(n)；Broad sense Hilbert becomes
Value computing unit is changed, for calculating the new seismic sequence GSK (n) generalized hilbert transform value GHk；Border detection value is given birth to
Into unit, for according to the generalized hilbert transform value GHkThe border detection value of the generation analysis sampling point；Geologic body side
Boundary's recognition unit, for generating a data volume according to multiple described border detection values, and carries out cutting along layer to the data volume
Piece is analyzed, and recognizes geological objects boundary.
In one embodiment, above-mentioned optimum angle of incidence to determining unit specifically for：Set one is with sampled point (x, y)
Central point, major and minor axis are respectively a, b horizontal window；The seismologic parameter of the sampled point is obtained by the horizontal window；Will
The seismologic parameter substitutes into sampled point coherent calculation formula, obtains the coherent value of the sampled point；Line direction, road side are set respectively
To the extreme value P_min, P_max, Q_min, Q_max of apparent formation dip；Determine that geologic body reflects boundary by certain step delta P, Δ Q
Apparent dip P, the Q in face；N is obtained according to Nyquist Sampling Theoremsp*nqIndividual inclination angle to (P, Q) and corresponding coherent value c (t, P, Q),
Wherein, np=(P_max-P_min)/Δ P, nq=(Q_max-Q_min)/Δ Q；With apparent dip to (PL,QM) centered on put choose
C*d apparent dip (Pl, Qm) and its corresponding coherent value c (t, Pl,Qm) a curved surface G (P, Q) is fitted to, wherein, Pl=PL±
C* Δs P, Qm=QM±d*ΔP；The corresponding inclination angle of maximum of the curved surface G (P, Q) is the optimal of the sampled point to (P, Q)
Inclination angle pair.
In one embodiment, above-mentioned Gauss function is：Wherein, α falling for standard deviation
Number.
In one embodiment, above-mentioned new seismic sequence GSK (n)=Sk(n) * G (n), k=1,2 ..., i.
In one embodiment, above-mentioned border detection value generation unit specifically for：The broad sense Hilbert is calculated to become
Change value GHkImaginary part absolute value IGHk；Judge multiple absolute value IGHkIn maximum, and using the maximum as
The border detection value of the analysis sampling point.
The beneficial effects of the present invention are with stronger noise suppression ability, can more accurately and comprehensively recognize variously
Geological objects boundary under matter environment.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, embodiment will be described below
In required for the accompanying drawing that uses be briefly described, it should be apparent that, drawings in the following description are only some of the present invention
Embodiment, for those of ordinary skill in the art, without having to pay creative labor, can also be according to these
Accompanying drawing obtains other accompanying drawings.
Fig. 1 is the flow chart of the boundary recognition of geological body method according to the embodiment of the present invention；
Fig. 2 is the window schematic diagram in the coherent calculation method according to the embodiment of the present invention；
Fig. 3 A and Fig. 3 B are the dip scanning schematic diagram according to the embodiment of the present invention；
Fig. 4 is the schematic diagram of the multi-direction reading amplitude process according to the embodiment of the present invention；
Fig. 5 is the seismic sequence figure in the direction 1 according to the embodiment of the present invention；
Fig. 6 is the analysis result for doing slice analysis along Silurian top surface on two data volumes according to the embodiment of the present invention
Schematic diagram；
Fig. 7 is to respond diagrammatic cross-section according to the seismic profile in the river course of the embodiment of the present invention；
Fig. 8 is the structural representation of the boundary recognition of geological body device according to the embodiment of the present invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made
Embodiment, belongs to the scope of protection of the invention.
The embodiment of the present invention provides a kind of boundary recognition of geological body method and device.The present invention is carried out below in conjunction with accompanying drawing
Describe in detail.
The embodiment of the present invention provides a kind of boundary recognition of geological body method, as shown in figure 1, the boundary recognition of geological body method
Mainly include the following steps that：
Step S101：Set a length L to be 2M+1 sliding window, and set length L rank value for N, obtain 2M+1
Sampled point, and Gauss function G (n) is calculated according to 2M+1 sampled point；
Step S102：The optimum angle of incidence at the local reflex interface where each sampled point is determined using dip scanning method
It is right；
Step S103：Respectively using each sampled point as analysis sampling point, and centered on analyzing sampling point, according to analysis sampling point
Optimum angle of incidence to reading the amplitude of 2M+1 point, the multiple seismic signal sequence S of generation respectively from multiple directionsk(n), k
=1,2 ..., i, n=0,1 ..., 2M, i=2M+1；
Step S104：According to earthquake signal sequence SkAnd Gauss function G (n) calculates new seismic sequence G (n)Sk(n)；
Step S105：Calculate new seismic sequence GSK (n) generalized hilbert transform value GHk；
Step S106：According to generalized hilbert transform value GHkThe border detection value of generation analysis sampling point；
Step S107：A data volume is generated according to multiple border detection values, and horizon slice analysis is carried out to data volume, is known
Other geological objects boundary.
By above-mentioned step S101~step S107, the optimum angle of incidence pair of sampled point and multiple sides are obtained by window function
Upward seismic sequence, and generalized hilbert transform is carried out to the seismic sequence, so as to obtain border detection value, recognize geology
Body border.Above-mentioned recognition methods has stronger noise suppression ability, can more accurately and comprehensively recognize under various geological environments
Geological objects boundary.
Below with reference to above steps, the boundary recognition of geological body method to the embodiment of the present invention is described in detail.
In above-mentioned step S101, set a length L to be 2M+1 sliding window, and length L rank value be set for N,
In embodiments of the present invention, the rank value N of the length of sliding window takes 1 or 2.2M+1 sampled point is obtained by the sliding window,
And Gauss function G (n) is calculated according to 2M+1 sampled point of acquisition,Wherein, M is natural number, n=
0,1 ..., 2M, α be standard deviation inverse, in embodiments of the present invention, α >=3.
Above-mentioned step S102, utilizes the local reflex interface where each sampled point of dip scanning method determination
Optimum angle of incidence pair.First, the coherent calculation method proposed according to Marfurt, one point centered on sampled point (x, y) of setting,
Major and minor axis is respectively a, b planar elliptical or the horizontal window of rectangle, as shown in Fig. 2 the plane window includes the earthquake of J roads
Road, analysis site coherent calculation formula：
Wherein, times (ms) or depth (m) of the τ for sampled point；P, q are respectively the local geology interface where sampled point
Along the apparent dip in x, y direction, unit is ms/m (time-domain) or m/m (Depth Domain)；U (t, x, y) be transverse and longitudinal coordinate be (x, y) when
Between be t earthquake number strong point amplitude；uH(t, x, y) is seismic channel u (t, x, y) Hilbert transform；K is that length includes w/
The vertical window of τ sampling point of Δ, the window size is 2w ms or m, and Δ τ is the sampling interval, and its unit is ms or m.
Explained in actual layer position in work, line, the extreme value P_min, P_max, Q_ of road direction apparent formation dip are set
Min, Q_max, P_min<p<P_max, Q_min<q<Q_max.Geologic body reflecting interface is searched by certain step delta P, Δ Q
Apparent dip P, Q, as shown in Fig. 3 A and Fig. 3 B, by Nyquist (Nyquist) Sampling Theorem：ΔP≤1/(2afmax), Δ Q≤1/
(2bfmax), wherein fmax is the highest frequency of seismic data.N can be obtainedp*nqIndividual inclination angle is to (P, Q) and its corresponding relevant
Value c (t, P, Q), wherein np=(P_max-P_min)/Δ P, nq=(Q_max-Q_min)/Δ Q.
When there is c (t, PL,QMDuring) >=c (t, P, Q), i.e. the maximum corresponding (P of coherent valueL,QM) it is sampled point location
Shake the corresponding apparent dip of reflecting surface.Due to taking apparent dip to (P, Q) by certain step-length (Δ P, Δ Q) discrete sampling, by step-length
In scanning process, even if c (t, PL,QM) >=c (t, P, Q), (PL,QM) also it is not optimal apparent dip pair.Therefore, with regarding
Inclination angle is to (PL,QM) centered on point choose c*d apparent dip (Pl, Qm), Pl=PL± c* Δs P, Qm=QM± d* Δ P, and its correspondingly
Coherent value c (t, Pl,Qm) a curved surface G (P, Q) is fitted to, G (P, Q) the corresponding inclination angle of maximum is exactly optimal to (P, Q)
Inclination angle pair.
Above-mentioned step S103：Respectively using each sampled point as analysis sampling point, and centered on analyzing sampling point, according to point
The optimum angle of incidence of sampling point is analysed to reading the amplitude of 2M+1 point, the multiple seismic signal sequence S of generation respectively from multiple directionsk
(n)。
Specifically, centered on sampled point, read the amplitude of 2M+1 sampled point respectively from nine directions, form 9
Seismic signal sequence.And 9 seismic signal sequences are designated as Sk(n), k=1,2 ..., 9, n=0,1 ..., 2M.Wherein, 9
Individual direction can specifically be obtained by direction as shown in Figure 4：Direction 1,2 is oversampled points and parallel with main profile and cross-track；
Direction 3 is oversampled points and parallel with depth direction；Direction 4,5 is respectively the same depth plane and and main profile of oversampled points
With direction of the cross-track in 45 degree of angles；Direction 6,7 for oversampled points and in the same cross-track plane, with depth
Direction or main profile direction are in the both direction of 45 degree of angles；Direction 8,9 is oversampled points and flat in same main profile
It is in the both direction of 45 degree of angles with depth direction or cross-track direction on face.
For example：On direction 1, exemplified by length of window is 9 sampled point, as shown in figure 5, intermediate point S1(4) it is analysis sample
Point, on direction 1, centered on analyzing sampling point, 4 sampled points is selected in both sides, and open from first point of opposite direction respectively
Begin, the amplitude of the sampled point in whole window is read, so that the seismic signal sequence S formed on direction 11(n), n=0,
1,....,8。
It should be noted that in embodiments of the present invention, be using read respectively from 9 directions the amplitude of sampled point as
Example, but the present invention is not limited thereto.
Seismic signal sequence S on above-mentioned multiple directions are generatedk(n) and after Gauss function G (n), by above-mentioned
Step S104, calculates new seismic sequence GSk(n).Specifically, the new seismic sequence GSK (n)=Sk(n) * G (n), k=1,
2 ..., 9, n=0,1 ..., 2M.
Further, step S105 is performed, new seismic sequence G is calculatedSK (n) generalized hilbert transform value GHk。GHk
For plural number, the absolute value IGH of its imaginary part is calculated respectively for 9 plural numbersk, k=1,2 ..., 9.Also, with reference to step S106,
Generate the border detection value of each sampled point.Specifically judging multiple absolute value IGHkIn maximum, and by the maximum
It is used as the border detection value of sampled point.
Above-mentioned steps S107, it is multi-faceted to each sampled point node-by-node algorithm on every master line, every contact side line
Generalized hilbert transform, ultimately forms a data volume, and horizon slice analysis, identification geologic body side are carried out to the data volume
Boundary, so as to clearly highlight the small-scale geological objects boundaries such as small fracture and narrow river course.
The boundary recognition of geological body method of the embodiment of the present invention is specifically described below in conjunction with an instantiation.
In this example, it is to be applied to boundary recognition of geological body method in the Siluric stratum in Harrar Ha Tang oil fields.Breathe out
The Siluric stratum for drawing Ha Tang oil fields is development Strike-slip faulted and meandering stream deposit system inside clastic deposited sediments stratum, stratum.
Seismic data cube is calculated using multi-faceted generalized hilbert transform Boundary Recognition technology and coherent technique, in two numbers
According to shown in (a) and (b) being on body along Silurian top surface in slice analysis, such as Fig. 6, the two ratio is adopted to same quantitative range
Shown with same colour code.The seismic response in meandering stream river course shows as arrow in seismic event drop-down, such as Fig. 7 on seismic profile
The signified position of head, seismic profile survey line position is located in Fig. 6 (a), (b) and is located in the signified survey line position of arrow, Fig. 7
Yellow seismic horizon at trough is Siluric stratum top surface seismic horizon.Seismic response shows as lineups drop-down at river course
(Fig. 6 (a)), the river course residing for E, F, G, H, I position shown in Fig. 7 is the river course of the succession sexual development of experience many phases, early stage
River course is broad, scale is big (in Fig. 6 (a) at ellipse), and later stage river channel evolution is the less branch channel of scale, and river course many phases are stacked
Development, seismic response is in mixed and disorderly strong transmitting, and noise is strong.Influenceed by horizontal, vertical average effect and noise, coherent technique is difficult to carve
The small-scale river location (at E, F, G, H, I position in Fig. 6 (b)) of sedimentation period development at the top of later stage Silurian is drawn,
And multi-faceted generalized hilbert transform bound test technology is not influenceed by early stage river course, small-scale river can be preferably recognized
Road (at E, F, G, H, I position in Fig. 6 (a)).
The embodiment of the present invention also provides a kind of boundary recognition of geological body device, as shown in figure 8, the boundary recognition of geological body is filled
Putting mainly includes：Gauss function computing unit 1, optimum angle of incidence are to determining unit 2, seismic signal sequence generating unit 3, newly
Sequence calculation sequence 4, generalized hilbert transform value computing unit 5, border detection value generation unit 6 and geological objects boundary is shaken to know
Other unit 7 etc..
Wherein, above-mentioned Gauss function computing unit 1 is used to set a length L to be 2M+1 sliding window, and sets
Length L rank value N, in embodiments of the present invention, the rank value N of the length of sliding window take 1 or 2.Obtained by the sliding window
2M+1 sampled point, and Gauss function G (n) is calculated according to 2M+1 sampled point,Wherein, M is certainly
So number, n=0,1 ..., 2M, α be standard deviation inverse, in embodiments of the present invention, α >=3.
Optimum angle of incidence is used to determine local reflex circle where each sampled point using dip scanning device to determining unit 2
The optimum angle of incidence pair in face.First, the coherent calculation method proposed according to Marfurt, sets one centered on sampled point (x, y)
Point, major and minor axis are respectively a, b planar elliptical or the horizontal window of rectangle, as shown in Fig. 2 the plane window is genuine comprising J
Shake road, analysis site coherent calculation formula：
Wherein, times (ms) or depth (m) of the τ for sampled point；P, q are respectively the local geology interface where sampled point
Along the apparent dip in x, y direction, unit is ms/m (time-domain) or m/m (Depth Domain)；U (t, x, y) be transverse and longitudinal coordinate be (x, y) when
Between be t earthquake number strong point amplitude；uH(t, x, y) is seismic channel u (t, x, y) Hilbert transform；K is that length includes w/
The vertical window of τ sampling point of Δ, the window size is 2w ms or m, and Δ τ is the sampling interval, and its unit is ms or m.
Explained in actual layer position in work, line, the extreme value P_min, P_max, Q_ of road direction apparent formation dip are set
Min, Q_max, P_min<p<P_max, Q_min<q<Q_max.Geologic body reflecting interface is searched by certain step delta P, Δ Q
Apparent dip P, Q, as shown in Fig. 3 A and Fig. 3 B, by Nyquist (Nyquist) Sampling Theorem：ΔP≤1/(2afmax), Δ Q≤1/
(2bfmax), wherein fmax is the highest frequency of seismic data.N can be obtainedp*nqIndividual inclination angle is to (P, Q) and its corresponding relevant
Value c (t, P, Q), wherein np=(P_max-P_min)/Δ P, nq=(Q_max-Q_min)/Δ Q.
When there is c (t, PL,QMDuring) >=c (t, P, Q), i.e. the maximum corresponding (P of coherent valueL,QM) it is sampled point location
Shake the corresponding apparent dip of reflecting surface.Due to taking apparent dip to (P, Q) by certain step-length (Δ P, Δ Q) discrete sampling, by step-length
In scanning process, even if c (t, PL,QM) >=c (t, P, Q), (PL,QM) also it is not optimal apparent dip pair.Therefore, with regarding
Inclination angle is to (PL,QM) centered on point choose c*d apparent dip (Pl, Qm), Pl=PL± c* Δs P, Qm=QM± d* Δ P, and its correspondingly
Coherent value c (t, Pl,Qm) a curved surface G (P, Q) is fitted to, G (P, Q) the corresponding inclination angle of maximum is exactly optimal to (P, Q)
Inclination angle pair.
Above-mentioned seismic signal sequence generating unit 3, for respectively using each sampled point as analysis sampling point, and with analyze
Centered on sampling point, according to analysis sampling point optimum angle of incidence to from multiple directions respectively read 2M+1 point amplitude, generate
Multiple seismic signal sequence Sk(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1.
Specifically, centered on sampled point, read the amplitude of 2M+1 sampled point respectively from nine directions, form 9
Seismic signal sequence.And 9 seismic signal sequences are designated as Sk(n), k=1,2 ..., 9, n=0,1 ..., 2M.Wherein, 9
Individual direction can specifically be obtained by direction as shown in Figure 4：Direction 1,2 is oversampled points and parallel with main profile and cross-track；
Direction 3 is oversampled points and parallel with depth direction；Direction 4,5 is respectively the same depth plane and and main profile of oversampled points
With direction of the cross-track in 45 degree of angles；Direction 6,7 for oversampled points and in the same cross-track plane, with depth
Direction or main profile direction are in the both direction of 45 degree of angles；Direction 8,9 is oversampled points and flat in same main profile
It is in the both direction of 45 degree of angles with depth direction or cross-track direction on face.
For example：On direction 1, exemplified by length of window is 9 sampled point, as shown in figure 5, intermediate point S1(4) it is analysis sample
Point, on direction 1, centered on analyzing sampling point, 4 sampled points is selected in both sides, and open from first point of opposite direction respectively
Begin, the amplitude of the sampled point in whole window is read, so that the seismic signal sequence S formed on direction 11(n), n=0,
1,....,8。
It should be noted that in embodiments of the present invention, be using read respectively from 9 directions the amplitude of sampled point as
Example, but the present invention is not limited thereto.
Seismic signal sequence S on above-mentioned multiple directions are generatedk(n) and after Gauss function G (n), by above-mentioned
New seismic sequence computing unit 4, calculate new seismic sequence GSk(n).Specifically, the new seismic sequence GSK (n)=Sk
(n) * G (n), k=1,2 ..., 9, n=0,1 ..., 2M.
Further, above-mentioned generalized hilbert transform value computing unit 5 is triggered, new seismic sequence G is calculatedSK's (n)
Generalized hilbert transform value GHk。GHkFor plural number, the absolute value IGH of its imaginary part is calculated respectively for 9 plural numbersk, k=1,
2 ..., 9.Also, with reference to above-mentioned border detection value generation unit 6, generate the border detection value of each sampled point.Specifically,
It is to judge multiple absolute value IGHkIn maximum, and using the maximum as sampled point border detection value.
Above-mentioned boundary recognition of geological body unit 7, for each sampled point on every master line, every contact side line
The multi-faceted generalized hilbert transform of node-by-node algorithm, ultimately forms a data volume, and carry out horizon slice to the data volume
Analysis, recognizes geological objects boundary, so as to clearly highlight the small-scale geological objects boundaries such as small fracture and narrow river course.
One of ordinary skill in the art will appreciate that realizing that all or part of step in above-described embodiment method can lead to
Cross program to instruct the hardware of correlation to complete, the program can be stored in a computer read/write memory medium, such as
ROM/RAM, magnetic disc, CD etc..
Particular embodiments described above, has been carried out further in detail to the purpose of the present invention, technical scheme and beneficial effect
Describe in detail it is bright, should be understood that the foregoing is only the present invention specific embodiment, the guarantor being not intended to limit the present invention
Scope is protected, within the spirit and principles of the invention, any modification, equivalent substitution and improvements done etc. should be included in this
Within the protection domain of invention.
Claims (10)
1. a kind of boundary recognition of geological body method, it is characterised in that described boundary recognition of geological body method includes：
Set a length L to be 2M+1 sliding window, and length L rank value is set for N, 2M+1 sampled point of acquisition, and according to
The 2M+1 sampled point calculates Gauss function G (n), and M is natural number；
Utilize the optimum angle of incidence pair at the local reflex interface where each sampled point of dip scanning method determination；
Respectively using each sampled point as analysis sampling point, and centered on the analysis sampling point, according to the analysis sampling point
Optimum angle of incidence to reading the amplitude of 2M+1 point, the multiple seismic signal sequence S of generation respectively from multiple directionsk(n), k
=1,2 ..., i, n=0,1 ..., 2M, i=2M+1；
According to the seismic signal sequence SkAnd Gauss function G (n) calculates new seismic sequence G (n)Sk(n)；
Calculate the new seismic sequence GSK (n) generalized hilbert transform value GHk；
According to the generalized hilbert transform value GHkThe border detection value of the generation analysis sampling point；
A data volume is generated according to multiple described border detection values, and horizon slice analysis, identification are carried out to the data volume
Geological objects boundary.
2. boundary recognition of geological body method according to claim 1, it is characterised in that determined using dip scanning method every
The optimum angle of incidence pair at the local reflex interface where the individual sampled point, including：
Setting one point, major and minor axis centered on sampled point (x, y) are respectively a, b horizontal window；
The seismologic parameter of the sampled point is obtained by the horizontal window；
The seismologic parameter is substituted into sampled point coherent calculation formula, the coherent value of the sampled point is obtained；
Line direction, the extreme value P_min, P_max, Q_min, Q_max of road direction apparent formation dip are set respectively；
Apparent dip P, Q of geologic body reflecting interface are determined by certain step delta P, Δ Q；
N is obtained according to nyquist sampling lawp*nqIndividual apparent dip to (P, Q) and corresponding coherent value c (t, P, Q), wherein, np
=(P_max-P_min)/Δ P, nq=(Q_max-Q_min)/Δ Q；
With apparent dip to (PL,QM) centered on point choose c*d apparent dip to (Pl, Qm) and its corresponding coherent value c (t, Pl,Qm)
A curved surface G (P, Q) is fitted to, wherein, Pl=PL± c* Δs P, Qm=QM±d*ΔQ；
The corresponding apparent dip of maximum of the curved surface G (P, Q) is to the optimal apparent dip pair of (P, Q) for the sampled point.
3. boundary recognition of geological body method according to claim 1, it is characterised in that described Gauss function is：
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>*</mo>
<mfrac>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mo>,</mo>
</mrow>
Wherein, α is the inverse of standard deviation.
4. boundary recognition of geological body method according to claim 1, it is characterised in that described new seismic sequence GSk(n)
=Sk(n) * G (n), k=1,2 ..., i.
5. boundary recognition of geological body method according to claim 1, it is characterised in that become according to the broad sense Hilbert
Change value GHkThe border detection value of the generation analysis sampling point, including：
Calculate the generalized hilbert transform value GHkImaginary part absolute value IGHk；
Judge multiple absolute value IGHkIn maximum, and using the maximum be used as it is described analysis sampling point border detection
Value.
6. a kind of boundary recognition of geological body device, it is characterised in that described boundary recognition of geological body device includes：
Gauss function computing unit, for setting a length L to be 2M+1 sliding window, and sets length L rank value N, obtains
2M+1 sampled point is taken, and Gauss function G (n) is calculated according to the 2M+1 sampled point, M is natural number；
Optimum angle of incidence is to determining unit, for utilizing local reflex circle where each sampled point of dip scanning device determination
The optimum angle of incidence pair in face；
Seismic signal sequence generating unit, for respectively using each sampled point as analysis sampling point, and with the analysis sample
Point centered on, according to it is described analysis sampling point optimum angle of incidence to from multiple directions respectively read 2M+1 point amplitude, give birth to
Into multiple seismic signal sequence Sk(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1；
New seismic sequence computing unit, for according to the seismic signal sequence SkAnd Gauss function G (n) calculates new earthquake (n)
Sequence GSk(n)；
Generalized hilbert transform value computing unit, for calculating the new seismic sequence GSK (n) generalized hilbert transform
Value GHk；
Border detection value generation unit, for according to the generalized hilbert transform value GHkThe border of the generation analysis sampling point
Detected value；
Boundary recognition of geological body unit, for generating a data volume according to multiple described border detection values, and to the data
Body carries out horizon slice analysis, recognizes geological objects boundary.
7. boundary recognition of geological body device according to claim 6, it is characterised in that the optimum angle of incidence is to determining unit
Specifically for：
Setting one point, major and minor axis centered on sampled point (x, y) are respectively a, b horizontal window；
The seismologic parameter of the sampled point is obtained by the horizontal window；
The seismologic parameter is substituted into sampled point coherent calculation formula, the coherent value of the sampled point is obtained；
Line direction, the extreme value P_min, P_max, Q_min, Q_max of road direction apparent formation dip are set respectively；
Apparent dip P, Q of geologic body reflecting interface are determined by certain step delta P, Δ Q；
N is obtained according to Nyquist Sampling Theoremsp*nqIndividual apparent dip to (P, Q) and corresponding coherent value c (t, P, Q), wherein, np
=(P_max-P_min)/Δ P, nq=(Q_max-Q_min)/Δ Q；
With apparent dip to (PL,QM) centered on point choose c*d apparent dip to (Pl, Qm) and its corresponding coherent value c (t, Pl,Qm)
A curved surface G (P, Q) is fitted to, wherein, Pl=PL± c* Δs P, Qm=QM±d*ΔQ；
The corresponding apparent dip of maximum of the curved surface G (P, Q) is to the optimal apparent dip pair of (P, Q) for the sampled point.
8. boundary recognition of geological body device according to claim 6, it is characterised in that described Gauss function is：
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>*</mo>
<mfrac>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mo>,</mo>
</mrow>
Wherein, α is the inverse of standard deviation.
9. boundary recognition of geological body device according to claim 6, it is characterised in that described new seismic sequence GSk(n)
=Sk(n) * G (n), k=1,2 ..., i.
10. boundary recognition of geological body device according to claim 6, it is characterised in that the border detection value generation is single
Member specifically for：
Calculate the generalized hilbert transform value GHkImaginary part absolute value IGHk；
Judge multiple absolute value IGHkIn maximum, and using the maximum be used as it is described analysis sampling point border detection
Value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510445467.1A CN105093300B (en) | 2015-07-27 | 2015-07-27 | A kind of boundary recognition of geological body method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510445467.1A CN105093300B (en) | 2015-07-27 | 2015-07-27 | A kind of boundary recognition of geological body method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105093300A CN105093300A (en) | 2015-11-25 |
| CN105093300B true CN105093300B (en) | 2017-10-17 |
Family
ID=54574154
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510445467.1A Active CN105093300B (en) | 2015-07-27 | 2015-07-27 | A kind of boundary recognition of geological body method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105093300B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107656309B (en) * | 2017-08-15 | 2019-09-06 | 中国石油天然气股份有限公司 | A kind of method and device of determining fault throw |
| CN108549102B (en) * | 2018-03-29 | 2020-03-17 | 西安交通大学 | Stratum structure curvature estimation method combining gradient structure tensor analysis and multi-window analysis |
| CN112505773B (en) * | 2020-12-15 | 2022-10-28 | 中石化石油工程技术服务有限公司 | Fracture inheritance quantitative calculation and judgment method |
| CN112749493A (en) * | 2021-01-25 | 2021-05-04 | 中国自然资源航空物探遥感中心 | Geologic body boundary detection method and system based on full magnetic gradient tensor eigenvalue |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6963804B2 (en) * | 2002-11-19 | 2005-11-08 | Saudi Arabian Oil Company | Seismic data processing method to enhance fault and channel display |
| CN104570078B (en) * | 2013-10-18 | 2017-10-27 | 中国石油化工股份有限公司 | A kind of solution cavity detection method of the similitude cross directional variations rate based on frequency domain inclination angle |
| CN103777241B (en) * | 2014-01-21 | 2016-06-08 | 成都理工大学 | Three dimensional seismic data rapid edge-detection method based on time domain generalized Hilbert transformation |
-
2015
- 2015-07-27 CN CN201510445467.1A patent/CN105093300B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105093300A (en) | 2015-11-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hall et al. | Fracture characterization at Valhall: Application of P-wave amplitude variation with offset and azimuth (AVOA) analysis to a 3D ocean-bottom data set | |
| CN103913776B (en) | The detection method of spherical boulder in shield-tunneling construction | |
| CN105425292B (en) | A kind of oil gas forecasting method and device | |
| CN105717542B (en) | The domain Wheeler based on 3-D seismics is sliced sand body engraving process | |
| CA2964893C (en) | Structure tensor constrained tomographic velocity analysis | |
| CN106772592B (en) | Diffracted wave focuses the analysis method and device of energy | |
| CN107817526B (en) | Prestack seismic gather segmented amplitude energy compensation method and system | |
| CN106556861B (en) | A kind of azimuthal AVO inversion method based on Omnibearing earthquake auto data | |
| CN105093300B (en) | A kind of boundary recognition of geological body method and device | |
| CN107728204A (en) | Based on the anisotropic crack prediction method of prestack compressional wave and system | |
| CN109188520A (en) | Thin reservoir thickness prediction method and device | |
| CN105259571A (en) | Stratum inclination angle detection method | |
| CN106772593B (en) | The imaging method and device of diffracted wave | |
| CN104375178B (en) | A kind of carbonate rock fractured cave method for predicting reservoir and device | |
| CN110320576A (en) | Method for forecasting advanced geology, device and storage medium | |
| CN105277976A (en) | Earthquake forward modelling method based on rock outcrop radar detection | |
| Guillemoteau et al. | Fast 3D multichannel deconvolution of electromagnetic induction loop-loop apparent conductivity data sets acquired at low induction numbers | |
| CN106772599B (en) | A kind of method and device calculating formation shear speed | |
| CN106842299B (en) | A method of the crack quantification prediction based on seismic properties | |
| Carpentier et al. | Semiautomated suppression of above-surface diffractions in GPR data | |
| CN108693560A (en) | A kind of scattering wave imaging method and system based on cross-correlation road | |
| US9791580B2 (en) | Methods and systems to separate wavefields using pressure wavefield data | |
| CN108089228B (en) | Interpretation data method and device for determining stratum rock characteristics | |
| CN112684498A (en) | Reservoir fracture prediction method and system based on wide-azimuth seismic data | |
| CN110261905A (en) | Complex value based on pitch angle control is concerned with microfault recognition methods |
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 |