CN105093300A - Geologic body boundary identification method and apparatus - Google Patents

Abstract

The invention provides a geologic body boundary identification method and an apparatus. The geologic body boundary identification method comprises the following steps of setting a slide window whose length L is 2M+1, setting an order value of the length L to be N, acquiring 2M+1 sampling points and calculating a Gauss window function G(n) according to the 2M+1 sampling points; using an inclination angle scanning method to determine an optimum inclination angle pair of a local reflection interface where each sampling point is located; taking each sampling point as an analysis sampling point respectively and taking the analysis sampling point as a center, reading amplitude values of the 2M+1 points from multiple directions according to an viewing dip angle of the analysis sampling point respectively and generating a plurality of earthquake signal sequences Sk(n), wherein the K is equal to 1, 2,..., I; the n is equal to 0, 1,..., 2M; and the i is equal to 2M+1; according to the earthquake signal sequences Sk(n) and the Gauss window function G(n), calculating a new earthquake sequence GSk(n); calculating a generalized Hilbert transform value GHk of the new earthquake sequence GSk(n); according to the generalized Hilbert transform value GHk, generating a boundary detection value of each analysis sampling point; according to the plurality of boundary detection values, generating a data body, carrying out across-layer slice analysis on the data body and identifying a geologic body boundary.

Description

A kind of boundary recognition of geological body method and device

Technical field

The invention relates to the interpretation technique of geophysical prospecting for oil geological data, particularly, is about a kind of boundary recognition of geological body method and the device that can identify geological objects boundary on a small scale.

Background technology

At present, coherent technique is widely used in detecting river course, fracture.Although coherent technique experienced by three generations's development, all there is average effect, be difficult to highlight narrow river course and the fracture of little turn-off.Meanwhile, coherent algorithm is to noise-sensitive, and in river course, mutual cut place of rupturing, seismic signal is mixed and disorderly, noise is strong, and coherent technique is difficult to accurately identify narrow river course and little fracture, can only identify the roughly growth region of river course, fracture.

The anti-noise ability of existing generalized hilbert transform bound test technology is strong, can highlight geological objects boundary on a small scale.In practice flow process, adopt window function to carry out generalized hilbert transform from main profile and cross-track both direction and detect different geological objects boundary.But often space distribution is complicated for the trend of plastid (river course, fracture) practically, and only carrying out Boundary Detection to main profile and cross-track both direction can not detect geological objects boundary all sidedly.

Summary of the invention

The fundamental purpose of the embodiment of the present invention is to provide a kind of boundary recognition of geological body method and device, more to like to detect geological objects boundary all sidedly.

To achieve these goals, the embodiment of the present invention provides a kind of boundary recognition of geological body method, described boundary recognition of geological body method comprises: arrange the moving window that a length L is 2M+1, and the rank value of length L is set for N, obtain 2M+1 sampled point, and calculating Gauss function G (n) according to described 2M+1 sampled point, M is natural number; Dip scanning method is utilized to determine the optimum angle of incidence pair at the local reflex interface at each described sampled point place; Respectively using sampled point described in each as analysis sampling point, and centered by described analysis sampling point, the apparent dip according to described analysis sampling point reads 2M+1 the amplitude put from multiple directions respectively, generates multiple seismic signal sequence S _k(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1; According to described seismic signal sequence S _kn () and Gauss function G (n) calculate new seismic sequence G _sk (n); Calculate described new seismic sequence G _sthe generalized hilbert transform value GH of k (n) _k; According to described generalized hilbert transform value GH _kgenerate the Boundary Detection value of described analysis sampling point;

Generate a data volume according to multiple described Boundary Detection value, and horizon slice analysis is carried out to described data volume, identify geological objects boundary.

In one embodiment, the above-mentioned dip scanning method that utilizes determines the optimum angle of incidence pair at the local reflex interface at each described sampled point place, comprising: arrange a point centered by sampled point (x, y), horizontal window that major and minor axis is respectively a, b; The seismologic parameter of described sampled point is obtained by described horizontal window; Described seismologic parameter is substituted into sampled point coherent calculation formula, obtain the coherent value of described sampled point; The extreme value P_min of line direction, direction, road apparent formation dip is set respectively, P_max, Q_min, Q_max; By apparent dip P, Q of certain step delta P, Δ Q plastid reflecting interface definitely; N is obtained according to nyquist sampling law _p* n _qindividual inclination angle to the coherent value c (t, P, Q) of (P, Q) and correspondence, wherein, n _p=(P_max-P_min)/Δ P, n _q=(Q_max-Q_min)/Δ Q; With apparent dip to (P _l, Q _m) centered by point choose c*d apparent dip (P _l, Q _m) and coherent value c (t, the P of correspondence _l, Q _m) fit to a curved surface G (P, Q), wherein, P _l=P _l± c* Δ P, Q _m=Q _m± d* Δ P; Inclination angle corresponding to the maximal value of described curved surface G (P, Q) is to (P, the Q) optimum angle of incidence pair for described sampled point.

In one embodiment, above-mentioned Gauss function is: wherein, α is the inverse of standard deviation.

In one embodiment, above-mentioned new seismic sequence G _sk (n)=S _k(n) * G (n), k=1,2 ..., i.

In one embodiment, above-mentioned according to described generalized hilbert transform value GH _kgenerate the Boundary Detection value of described analysis sampling point, comprising: calculate described generalized hilbert transform value GH _kthe absolute value IGH of imaginary part _k; Judge multiple described absolute value IGH _kin maximal value, and using the Boundary Detection value of described maximal value as described analysis sampling point.

The embodiment of the present invention also provides a kind of boundary recognition of geological body device, described boundary recognition of geological body device comprises: Gauss function computing unit, for arranging the moving window that a length L is 2M+1, and the rank value N of length L is set, obtain 2M+1 sampled point, and calculating Gauss function G (n) according to described 2M+1 sampled point, M is natural number; Optimum angle of incidence to determining unit, for the optimum angle of incidence pair utilizing dip scanning device to determine the local reflex interface at each described sampled point place; Seismic signal sequence generating unit, for respectively using sampled point described in each as analysis sampling point, and centered by described analysis sampling point, according to the optimum angle of incidence of described analysis sampling point to reading 2M+1 the amplitude put respectively from multiple directions, generate multiple seismic signal sequence S _k(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1; New seismic sequence computing unit, for according to described seismic signal sequence S _kn () and Gauss function G (n) calculate new seismic sequence G _sk (n); Generalized hilbert transform value computing unit, for calculating described new seismic sequence G _sthe generalized hilbert transform value GH of k (n) _k; Boundary Detection value generation unit, for according to described generalized hilbert transform value GH _kgenerate the Boundary Detection value of described analysis sampling point; Boundary recognition of geological body unit, for generating a data volume according to multiple described Boundary Detection value, and carrying out horizon slice analysis to described data volume, identifying geological objects boundary.

In one embodiment, above-mentioned optimum angle of incidence to determining unit specifically for: a point centered by sampled point (x, y) is set, horizontal window that major and minor axis is respectively a, b; The seismologic parameter of described sampled point is obtained by described horizontal window; Described seismologic parameter is substituted into sampled point coherent calculation formula, obtain the coherent value of described sampled point; The extreme value P_min of line direction, direction, road apparent formation dip is set respectively, P_max, Q_min, Q_max; By apparent dip P, Q of certain step delta P, Δ Q plastid reflecting interface definitely; N is obtained according to Nyquist Sampling Theorem _p* n _qindividual inclination angle to the coherent value c (t, P, Q) of (P, Q) and correspondence, wherein, n _p=(P_max-P_min)/Δ P, n _q=(Q_max-Q_min)/Δ Q; With apparent dip to (P _l, Q _m) centered by point choose c*d apparent dip (P _l, Q _m) and coherent value c (t, the P of correspondence _l, Q _m) fit to a curved surface G (P, Q), wherein, P _l=P _l± c* Δ P, Q _m=Q _m± d* Δ P; Inclination angle corresponding to the maximal value of described curved surface G (P, Q) is to (P, the Q) optimum angle of incidence pair for described sampled point.

In one embodiment, above-mentioned Gauss function is: wherein, α is the inverse of standard deviation.

In one embodiment, above-mentioned new seismic sequence G _sk (n)=S _k(n) * G (n), k=1,2 ..., i.

In one embodiment, above-mentioned Boundary Detection value generation unit specifically for: calculate described generalized hilbert transform value GH _kthe absolute value IGH of imaginary part _k; Judge multiple described absolute value IGH _kin maximal value, and using the Boundary Detection value of described maximal value as described analysis sampling point.

Beneficial effect of the present invention is, has stronger noise suppression ability, can identify the geological objects boundary under various geologic media more accurately, all sidedly.

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, below the accompanying drawing used required in describing embodiment is briefly described, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

Fig. 1 is the process flow diagram of the boundary recognition of geological body method according to the embodiment of the present invention;

Fig. 2 is according to the window schematic diagram in the coherent calculation method of the embodiment of the present invention;

Fig. 3 A and Fig. 3 B is 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 schematic diagram doing slice analysis on two data volumes along Silurian end face according to the embodiment of the present invention;

Fig. 7 is the seismic section response diagrammatic cross-section according to 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, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong 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.Below in conjunction with accompanying drawing, the present invention is described in detail.

The embodiment of the present invention provides a kind of boundary recognition of geological body method, and as shown in Figure 1, this boundary recognition of geological body method mainly comprises the following steps:

Step S101: the moving window that a length L is 2M+1 is set, and the rank value of length L is set for N, obtain 2M+1 sampled point, and calculate Gauss function G (n) according to 2M+1 sampled point;

Step S102: utilize dip scanning method to determine the optimum angle of incidence pair at the local reflex interface at each sampled point place;

Step S103: respectively using each sampled point as analysis sampling point, and to analyze centered by sampling point, according to analyzing the optimum angle of incidence of sampling point to reading 2M+1 the amplitude put respectively from multiple directions, generate multiple seismic signal sequence S _k(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1;

Step S104: according to earthquake burst S _kn () and Gauss function G (n) calculate new seismic sequence G _sk (n);

Step S105: calculate new seismic sequence G _sthe generalized hilbert transform value GH of k (n) _k;

Step S106: according to generalized hilbert transform value GH _kgenerate the Boundary Detection value analyzing sampling point;

Step S107: generate a data volume according to multiple Boundary Detection value, and carry out horizon slice analysis to data volume, identifies geological objects boundary.

By above-mentioned step S101 ~ step S107, the optimum angle of incidence being obtained sampled point by window function to and multiple directions on seismic sequence, and generalized hilbert transform is carried out to this seismic sequence, thus obtains Boundary Detection value, identify geological objects boundary.Above-mentioned recognition methods has stronger noise suppression ability, can identify the geological objects boundary under various geologic media more accurately, all sidedly.

Below with reference to above steps, the boundary recognition of geological body method of the embodiment of the present invention is described in detail.

In above-mentioned step S101, arrange the moving window that a length L is 2M+1, and arrange the rank value of length L for N, in embodiments of the present invention, the rank value N of the length of moving window gets 1 or 2.Obtain 2M+1 sampled point by this moving window, and calculate Gauss function G (n) according to 2M+1 the sampled point obtained, wherein, M is natural number, n=0,1 ..., 2M, α are the inverse of standard deviation, in embodiments of the present invention, and α>=3.

Above-mentioned step S102, utilizes dip scanning method to determine the optimum angle of incidence pair at the local reflex interface at each described sampled point place.First, according to the coherent calculation method that Marfurt proposes, one is arranged with sampled point (x, y) centered by, point, major and minor axis are respectively the planar elliptical of a, b or the horizontal window of rectangle, as shown in Figure 2, this plane window comprises J road seismic trace, analysis site coherent calculation formula:

$\mathrm{\σ}(\mathrm{\τ},p,q)=\frac{\underset{k=-K}{\overset{+K}{\Σ}}\left\{{\[\underset{j=1}{\overset{J}{\Σ}}u(\mathrm{\τ}+k\mathrm{\Δ}t-{\mathrm{px}}_j-{\mathrm{qy}}_j,x_j,y_j)\]}^2+{\[\underset{j=1}{\overset{J}{\Σ}}{u}^H(\mathrm{\τ}+k\mathrm{\Δ}t-{\mathrm{px}}_j-{\mathrm{qy}}_j,x_j,y_j)\]}^2\right\}}{J\underset{k=-K}{\overset{+K}{\Σ}}\underset{j=1}{\overset{J}{\Σ}}\{{\[u(\mathrm{\τ}+k\mathrm{\Δ}t-{\mathrm{px}}_j-{\mathrm{qy}}_j,x_j,y_j)\]}^2+{\[u^H(\mathrm{\τ}+k\mathrm{\Δ}t-{\mathrm{px}}_j-{\mathrm{qy}}_j,x_j,y_j)\]}^2\}},$

Wherein, τ is time (ms) or the degree of depth (m) of sampled point; P, q are respectively the apparent dip of local geological interface along x, y direction at sampled point place, and unit is ms/m (time domain) or m/m (Depth Domain); U (t, x, y) is the earthquake number strong point amplitude of t for transverse and longitudinal coordinate is (x, y) time; u ^h(t, x, y) is the Hilbert transform of seismic trace u (t, x, y); K is the vertical window that length comprises w/ Δ τ sampling point, and this window size is 2wms or m, and Δ τ is sampling interval, and its unit is ms or m.

In the layer position interpretation work of reality, the extreme value P_min of line, direction, road apparent formation dip is set, P_max, Q_min, Q_max, P_min<p<P_max, Q_min<q<Q_max.Apparent dip P, Q of plastid reflecting interface is raked out, as shown in Fig. 3 A and Fig. 3 B, by Nyquist (Nyquist) Sampling Theorem: Δ P≤1/ (2af by certain step delta P, Δ Q _max), Δ Q≤1/ (2bf _max), wherein fmax is the highest frequency of seismic data.N can be obtained _p* n _qindividual inclination angle is to the coherent value c (t, P, Q) of (P, Q) and correspondence thereof, wherein n _p=(P_max-P_min)/Δ P, n _q=(Q_max-Q_min)/Δ Q.

When there is c (t, P _l, Q _mtime)>=c (t, P, Q), (P that namely maximum coherent value is corresponding _l, Q _m) be apparent dip corresponding to sampled point place earthquake reflecting surface.Owing to getting apparent dip to (P, Q), by step scan process, even if c is (t, P by certain step-length (Δ P, Δ Q) discrete sampling _l, Q _m)>=c (t, P, Q), (P _l, Q _m) be also not best apparent dip pair.Therefore, with apparent dip to (P _l, Q _m) centered by point choose c*d apparent dip (P _l, Q _m), P _l=P _l± c* Δ P, Q _m=Q _m± d* Δ P, and coherent value c (t, the P of correspondence _l, Q _m) fitting to a curved surface G (P, Q), the inclination angle that the maximal value of G (P, Q) is corresponding is exactly best inclination angle pair to (P, Q).

Above-mentioned step S103: respectively using each sampled point as analysis sampling point, and to analyze centered by sampling point, according to analyzing the optimum angle of incidence of sampling point to reading 2M+1 the amplitude put respectively from multiple directions, generate multiple seismic signal sequence S _k(n).

Particularly, centered by sampled point, read the amplitude of 2M+1 sampled point respectively from nine directions, form 9 seismic signal sequences.And 9 seismic signal sequences are designated as S _k(n), k=1,2 ..., 9, n=0,1 ..., 2M.Wherein, 9 directions specifically can obtain by direction as shown in Figure 4: direction 1,2 is oversampled points and parallel with cross-track with main profile; Direction 3 is oversampled points and parallel with depth direction; Direction 4,5 is respectively the same depth plane of oversampled points and is the direction of 45 degree of angles with main profile and cross-track; That direction 6,7 is oversampled points and be in same cross-track plane, be the both direction of 45 degree of angles with depth direction or main profile direction; That direction 8,9 is oversampled points and be in same main profile plane, be the both direction of 45 degree of angles with depth direction or cross-track direction.

Such as: on direction 1, length of window is 9 sampled points is example, as shown in Figure 5, and intermediate point S ₁(4) for analyzing sampling point, on direction 1, to analyze centered by sampling point, select 4 sampled points respectively in both sides, and from opposite direction first point, read the amplitude of the sampled point in whole window, thus form the seismic signal sequence S on direction 1 ₁(n), n=0,1 ...., 8.

It should be noted that, in embodiments of the present invention, be the amplitude to read sampled point from 9 directions respectively, but the present invention is not limited to this.

Generating the seismic signal sequence S in above-mentioned multiple directions _kafter (n) and Gauss function G (n), by above-mentioned steps S104, calculate new seismic sequence G _sk (n).Particularly, this new seismic sequence G _sk (n)=S _k(n) * G (n), k=1,2 ..., 9, n=0,1 ..., 2M.

Further, perform step S105, calculate new seismic sequence G _sthe generalized hilbert transform value GH of k (n) _k.GH _kfor plural number, 9 plural numbers are calculated respectively to the absolute value IGH of its imaginary part _k, k=1,2 ..., 9.Further, integrating step S106, generates the Boundary Detection value of each sampled point.Particularly, be judge multiple absolute value IGH _kin maximal value, and using the Boundary Detection value of this maximal value as sampled point.

Above-mentioned steps S107, the generalized hilbert transform multi-faceted to each sampled point node-by-node algorithm on every bar master line, every bar contact side line, final formation data volume, and horizon slice analysis is carried out to this data volume, identify geological objects boundary, thus the small-scale geological objects boundaries such as little fracture and narrow river course can be highlighted clearly.

Be specifically described below in conjunction with the boundary recognition of geological body method of an instantiation to the embodiment of the present invention.

In this example, be that boundary recognition of geological body method is applied in the Siluric stratum in Harrar Ha Tang oil field.The Siluric stratum in Harrar Ha Tang oil field is clastic deposited sediments stratum, and Strike-slip faulted and meandering stream deposit system are grown in inside, stratum.Multi-faceted generalized hilbert transform Boundary Recognition technology and coherent technique is adopted to calculate seismic data volume, two data volumes do slice analysis along Silurian end face, as shown in (a) and (b) in Fig. 6, the two ratio, to same quantitative range, adopts same colour code to show.It is drop-down that the seismic response in meandering river river course shows as seismic event on seismic section, as the position of arrow indication in Fig. 7, seismic section survey line position is arranged in the survey line position of Fig. 6 (a), (b) arrow indication, and the yellow seismic horizon being positioned at trough place in Fig. 7 is Siluric stratum end face seismic horizon.River course place seismic response shows as lineups drop-down (Fig. 6 (a)), the river course residing for E, F, G, H, I position shown in Fig. 7 is the river course that the inheritance experiencing many phases is grown, early stage river course is broad, scale large (in Fig. 6 (a) oval place), later stage river channel evolution is the branch channel that scale is less, the stacked growth of river course many phases, seismic response is in mixed and disorderly strong transmitting, and noise is strong.By the impact of horizontal stroke, vertical average effect and noise, coherent technique is difficult to the small-scale river location (E, F, G, H, I position in Fig. 6 (b)) portraying the growth in period of later stage Silurian top sedimentation, and multi-faceted generalized hilbert transform bound test technology is not by the impact in early stage river course, can identify preferably river course (E, F, G, H, I position in Fig. 6 (a)) on a small scale.

The embodiment of the present invention also provides a kind of boundary recognition of geological body device, as shown in Figure 8, this boundary recognition of geological body device mainly comprises: Gauss function computing unit 1, optimum angle of incidence are to determining unit 2, seismic signal sequence generating unit 3, new seismic sequence computing unit 4, generalized hilbert transform value computing unit 5, Boundary Detection value generation unit 6 and boundary recognition of geological body unit 7 etc.

Wherein, above-mentioned Gauss function computing unit 1 for arranging the moving window that a length L is 2M+1, and arranges the rank value N of length L, and in embodiments of the present invention, the rank value N of the length of moving window gets 1 or 2.Obtain 2M+1 sampled point by this moving window, and calculate Gauss function G (n) according to 2M+1 sampled point, wherein, M is natural number, n=0,1 ..., 2M, α are the inverse of standard deviation, in embodiments of the present invention, and α>=3.

Optimum angle of incidence is to the optimum angle of incidence pair of determining unit 2 for utilizing dip scanning device to determine the local reflex interface at each sampled point place.First, according to the coherent calculation method that Marfurt proposes, one is arranged with sampled point (x, y) centered by, point, major and minor axis are respectively the planar elliptical of a, b or the horizontal window of rectangle, as shown in Figure 2, this plane window comprises J road seismic trace, analysis site coherent calculation formula:

$\mathrm{\&sigma;}(\mathrm{\&tau;},p,q)=\frac{\underset{k=-K}{\overset{+K}{\&Sigma;}}\left\{{\&lsqb;\underset{j=1}{\overset{J}{\&Sigma;}}u(\mathrm{\&tau;}+k\mathrm{\&Delta;}t-{\mathrm{px}}_j-{\mathrm{qy}}_j,x_j,y_j)\&rsqb;}^2+{\&lsqb;\underset{j=1}{\overset{J}{\&Sigma;}}{u}^H(\mathrm{\&tau;}+k\mathrm{\&Delta;}t-{\mathrm{px}}_j-{\mathrm{qy}}_j,x_j,y_j)\&rsqb;}^2\right\}}{J\underset{k=-K}{\overset{+K}{\&Sigma;}}\underset{j=1}{\overset{J}{\&Sigma;}}\{{\&lsqb;u(\mathrm{\&tau;}+k\mathrm{\&Delta;}t-{\mathrm{px}}_j-{\mathrm{qy}}_j,x_j,y_j)\&rsqb;}^2+{\&lsqb;u^H(\mathrm{\&tau;}+k\mathrm{\&Delta;}t-{\mathrm{px}}_j-{\mathrm{qy}}_j,x_j,y_j)\&rsqb;}^2\}},$

Wherein, τ is time (ms) or the degree of depth (m) of sampled point; P, q are respectively the apparent dip of local geological interface along x, y direction at sampled point place, and unit is ms/m (time domain) or m/m (Depth Domain); U (t, x, y) is the earthquake number strong point amplitude of t for transverse and longitudinal coordinate is (x, y) time; u ^h(t, x, y) is the Hilbert transform of seismic trace u (t, x, y); K is the vertical window that length comprises w/ Δ τ sampling point, and this window size is 2wms or m, and Δ τ is sampling interval, and its unit is ms or m.

In the layer position interpretation work of reality, the extreme value P_min of line, direction, road apparent formation dip is set, P_max, Q_min, Q_max, P_min<p<P_max, Q_min<q<Q_max.Apparent dip P, Q of plastid reflecting interface is raked out, as shown in Fig. 3 A and Fig. 3 B, by Nyquist (Nyquist) Sampling Theorem: Δ P≤1/ (2af by certain step delta P, Δ Q _max), Δ Q≤1/ (2bf _max), wherein fmax is the highest frequency of seismic data.N can be obtained _p* n _qindividual inclination angle is to the coherent value c (t, P, Q) of (P, Q) and correspondence thereof, wherein n _p=(P_max-P_min)/Δ P, n _q=(Q_max-Q_min)/Δ Q.

When there is c (t, P _l, Q _mtime)>=c (t, P, Q), (P that namely maximum coherent value is corresponding _l, Q _m) be apparent dip corresponding to sampled point place earthquake reflecting surface.Owing to getting apparent dip to (P, Q), by step scan process, even if c is (t, P by certain step-length (Δ P, Δ Q) discrete sampling _l, Q _m)>=c (t, P, Q), (P _l, Q _m) be also not best apparent dip pair.Therefore, with apparent dip to (P _l, Q _m) centered by point choose c*d apparent dip (P _l, Q _m), P _l=P _l± c* Δ P, Q _m=Q _m± d* Δ P, and coherent value c (t, the P of correspondence _l, Q _m) fitting to a curved surface G (P, Q), the inclination angle that the maximal value of G (P, Q) is corresponding is exactly best inclination angle pair to (P, Q).

Above-mentioned seismic signal sequence generating unit 3, for respectively using each sampled point as analysis sampling point, and to analyze centered by sampling point, according to analyzing the optimum angle of incidence of sampling point to the amplitude reading 2M+1 point from multiple directions respectively, generate multiple seismic signal sequence S _k(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1.

Particularly, centered by sampled point, read the amplitude of 2M+1 sampled point respectively from nine directions, form 9 seismic signal sequences.And 9 seismic signal sequences are designated as S _k(n), k=1,2 ..., 9, n=0,1 ..., 2M.Wherein, 9 directions specifically can obtain by direction as shown in Figure 4: direction 1,2 is oversampled points and parallel with cross-track with main profile; Direction 3 is oversampled points and parallel with depth direction; Direction 4,5 is respectively the same depth plane of oversampled points and is the direction of 45 degree of angles with main profile and cross-track; That direction 6,7 is oversampled points and be in same cross-track plane, be the both direction of 45 degree of angles with depth direction or main profile direction; That direction 8,9 is oversampled points and be in same main profile plane, be the both direction of 45 degree of angles with depth direction or cross-track direction.

Such as: on direction 1, length of window is 9 sampled points is example, as shown in Figure 5, and intermediate point S ₁(4) for analyzing sampling point, on direction 1, to analyze centered by sampling point, select 4 sampled points respectively in both sides, and from opposite direction first point, read the amplitude of the sampled point in whole window, thus form the seismic signal sequence S on direction 1 ₁(n), n=0,1 ...., 8.

It should be noted that, in embodiments of the present invention, be the amplitude to read sampled point from 9 directions respectively, but the present invention is not limited to this.

Generating the seismic signal sequence S in above-mentioned multiple directions _kafter (n) and Gauss function G (n), by above-mentioned new seismic sequence computing unit 4, calculate new seismic sequence G _sk (n).Particularly, this new seismic sequence G _sk (n)=S _k(n) * G (n), k=1,2 ..., 9, n=0,1 ..., 2M.

Further, trigger above-mentioned generalized hilbert transform value computing unit 5, calculate new seismic sequence G _sthe generalized hilbert transform value GH of k (n) _k.GH _kfor plural number, 9 plural numbers are calculated respectively to the absolute value IGH of its imaginary part _k, k=1,2 ..., 9.Further, in conjunction with above-mentioned Boundary Detection value generation unit 6, the Boundary Detection value of each sampled point is generated.Particularly, be judge multiple absolute value IGH _kin maximal value, and using the Boundary Detection value of this maximal value as sampled point.

Above-mentioned boundary recognition of geological body unit 7, for the generalized hilbert transform multi-faceted to each sampled point node-by-node algorithm on every bar master line, every bar contact side line, final formation data volume, and horizon slice analysis is carried out to this data volume, identify geological objects boundary, thus the small-scale geological objects boundaries such as little fracture and narrow river course can be highlighted clearly.

One of ordinary skill in the art will appreciate that the hardware that all or part of step realized in above-described embodiment method can carry out instruction relevant by program has come, this program can be stored in a computer read/write memory medium, such as ROM/RAM, magnetic disc, CD etc.

Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (10)

1. a boundary recognition of geological body method, is characterized in that, described boundary recognition of geological body method comprises:

Arrange the moving window that a length L is 2M+1, and arrange the rank value of length L for N, obtain 2M+1 sampled point, and calculate Gauss function G (n) according to described 2M+1 sampled point, M is natural number;

Dip scanning method is utilized to determine the optimum angle of incidence pair at the local reflex interface at each described sampled point place;

Respectively using sampled point described in each as analysis sampling point, and centered by described analysis sampling point, according to the optimum angle of incidence of described analysis sampling point to reading 2M+1 the amplitude put respectively from multiple directions, generate multiple seismic signal sequence S _k(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1;

According to described seismic signal sequence S _kn () and Gauss function G (n) calculate new seismic sequence G _sk (n);

Calculate described new seismic sequence G _sthe generalized hilbert transform value GH of k (n) _k;

According to described generalized hilbert transform value GH _kgenerate the Boundary Detection value of described analysis sampling point;

Generate a data volume according to multiple described Boundary Detection value, and horizon slice analysis is carried out to described data volume, identify geological objects boundary.

2. boundary recognition of geological body method according to claim 1, is characterized in that, utilizes dip scanning method to determine the optimum angle of incidence pair at the local reflex interface at each described sampled point place, comprising:

A point centered by sampled point (x, y) is set, horizontal window that major and minor axis is respectively a, b;

The seismologic parameter of described sampled point is obtained by described horizontal window;

Described seismologic parameter is substituted into sampled point coherent calculation formula, obtain the coherent value of described sampled point;

The extreme value P_min of line direction, direction, road apparent formation dip is set respectively, P_max, Q_min, Q_max;

By apparent dip P, Q of certain step delta P, Δ Q plastid reflecting interface definitely;

N is obtained according to nyquist sampling law _p* n _qindividual inclination angle to the coherent value c (t, P, Q) of (P, Q) and correspondence, wherein, n _p=(P_max-P_min)/Δ P, n _q=(Q_max-Q_min)/Δ Q;

With apparent dip to (P _l, Q _m) centered by point choose c*d apparent dip (P _l, Q _m) and coherent value c (t, the P of correspondence _l, Q _m) fit to a curved surface G (P, Q), wherein, P _l=P _l± c* Δ P, Q _m=Q _m± d* Δ P;

Inclination angle corresponding to the maximal value of described curved surface G (P, Q) is to (P, the Q) optimum angle of incidence pair for described sampled point.

3. boundary recognition of geological body method according to claim 1, is characterized in that, described Gauss function is:

Wherein, α is the inverse of standard deviation.

4. boundary recognition of geological body method according to claim 1, is characterized in that, described new seismic sequence G _sk (n)=S _k(n) * G (n), k=1,2 ..., i.

5. boundary recognition of geological body method according to claim 1, is characterized in that, according to described generalized hilbert transform value GH _kgenerate the Boundary Detection value of described analysis sampling point, comprising:

Calculate described generalized hilbert transform value GH _kthe absolute value IGH of imaginary part _k;

Judge multiple described absolute value IGH _kin maximal value, and using the Boundary Detection value of described maximal value as described analysis sampling point.

6. a boundary recognition of geological body device, is characterized in that, described boundary recognition of geological body device comprises:

Gauss function computing unit, for arranging the moving window that a length L is 2M+1, and arranges the rank value N of length L, obtains 2M+1 sampled point, and calculates Gauss function G (n) according to described 2M+1 sampled point, and M is natural number;

Optimum angle of incidence to determining unit, for the optimum angle of incidence pair utilizing dip scanning device to determine the local reflex interface at each described sampled point place;

Seismic signal sequence generating unit, for respectively using sampled point described in each as analysis sampling point, and centered by described analysis sampling point, according to the optimum angle of incidence of described analysis sampling point to reading 2M+1 the amplitude put respectively from multiple directions, generate multiple seismic signal sequence S _k(n), k=1,2 ..., i, n=0,1 ..., 2M, i=2M+1;

New seismic sequence computing unit, for according to described seismic signal sequence S _kn () and Gauss function G (n) calculate new seismic sequence G _sk (n);

Generalized hilbert transform value computing unit, for calculating described new seismic sequence G _sthe generalized hilbert transform value GH of k (n) _k;

Boundary Detection value generation unit, for according to described generalized hilbert transform value GH _kgenerate the Boundary Detection value of described analysis sampling point;

Boundary recognition of geological body unit, for generating a data volume according to multiple described Boundary Detection value, and carrying out horizon slice analysis to described data volume, identifying geological objects boundary.

7. boundary recognition of geological body device according to claim 6, is characterized in that, described optimum angle of incidence to determining unit specifically for:

A point centered by sampled point (x, y) is set, horizontal window that major and minor axis is respectively a, b;

The seismologic parameter of described sampled point is obtained by described horizontal window;

Described seismologic parameter is substituted into sampled point coherent calculation formula, obtain the coherent value of described sampled point;

The extreme value P_min of line direction, direction, road apparent formation dip is set respectively, P_max, Q_min, Q_max;

By apparent dip P, Q of certain step delta P, Δ Q plastid reflecting interface definitely;

N is obtained according to Nyquist Sampling Theorem _p* n _qindividual inclination angle to the coherent value c (t, P, Q) of (P, Q) and correspondence, wherein, n _p=(P_max-P_min)/Δ P, n _q=(Q_max-Q_min)/Δ Q;

With apparent dip to (P _l, Q _m) centered by point choose c*d apparent dip (P _l, Q _m) and coherent value c (t, the P of correspondence _l, Q _m) fit to a curved surface G (P, Q), wherein, P _l=P _l± c* Δ P, Q _m=Q _m± d* Δ P;

Inclination angle corresponding to the maximal value of described curved surface G (P, Q) is to (P, the Q) optimum angle of incidence pair for described sampled point.

8. boundary recognition of geological body device according to claim 6, is characterized in that, described Gauss function is:

Wherein, α is the inverse of standard deviation.

9. boundary recognition of geological body device according to claim 6, is characterized in that, described new seismic sequence G _sk (n)=S _k(n) * G (n), k=1,2 ..., i.

10. boundary recognition of geological body device according to claim 6, is characterized in that, described Boundary Detection value generation unit specifically for:

Calculate described generalized hilbert transform value GH _kthe absolute value IGH of imaginary part _k;

Judge multiple described absolute value IGH _kin maximal value, and using the Boundary Detection value of described maximal value as described analysis sampling point.