CN102243319A - Detection method and apparatus for underground discontinuum - Google Patents

Abstract

The embodiment of the invention provides a detection method and an apparatus for underground discontinuum. The method comprises the following steps: obtaining earthquake data; selecting a size of a computation window according to a geological structure corresponding to the earthquake data; selecting a time range in which a time length is larger than the size of the computation window; obtaining a curvature of each to-be-detected data point in the earthquake data by utilizing a related algorithm, the computation window selected in advance, and the time range selected in advance; obtaining a general curvature of the whole earthquake data; and detecting an underground crack and an underground fault according to the general curvature of the whole earthquake data. According to the method, a calculating process of a general curvature is direct with small error and high efficiency, so that an identification capability for discontinuum can be improved substantially.

Description

A kind of detection method of underground discontinuum and device

Technical field

The present invention relates to the seismic exploration technique field, be specifically related to a kind of detection method and device of underground discontinuum, be particularly related to the new method of the multiple dimensioned seismic volume curvature of a kind of calculating in seismic data interpretation attribute, new multiple dimensioned earthquake curvature attribute can be applied in the oil exploration and development fields, discerns the fracture and the crack of different brackets.

Background technology

Earthquake curvature attribute and subterranean fracture etc. are in close relations.At present, Chang Yong earthquake curvature computing method have two kinds: a kind of curvature of face computing method that are based on layer bit data; Another kind is to calculate reflection inclination angle, stratum and orientation earlier, asks for the computing method of body curvature again by coordinate transform.

The curvature of face computing method based on layer bit data of prior art at first are described.In order to calculate curvature more specific on the curved surface, need the quadric surface of a part of match, in Fig. 1, eight points adjacent with it with current point (z5) come quadric surface of match on least squares sense.Adopt 3 * 3 grid methods, this local secondary curve surface definition is as follows:

z＝ax ²+by ²+cxy+dx+ey+f， (1)

Wherein

$a=\frac{1}{2}\frac{{\∂}^2}{\∂x}\frac{z}{2}=\frac{z1+z3+z4+z6+z7+z9}{6\mathrm{\Δ}{x}^2}-\frac{z2+z5+z8}{3\mathrm{\Δ}{x}^2},---\left(2\right)$

$b=\frac{1}{2}\frac{{\&PartialD;}^{2}z}{{\&PartialD;\mathrm{<figure-callout\; id="2"\; label="y"\; filenames="HSA00000123211900041.png"\; state="\{\{state\}\}">y</figure-callout>}}^2}=\frac{z1+z2+z3+z7+z8+z9}{6\mathrm{\&Delta;}{y}^2}-\frac{z4+z5+z6}{3\mathrm{\&Delta;}{\mathrm{<figure-callout\; id="2"\; label="y"\; filenames="HSA00000123211900041.png"\; state="\{\{state\}\}">y</figure-callout>}}^2},---\left(3\right)$

$c=\frac{1}{2}\frac{{\&PartialD;}^{2}z}{\&PartialD;x\&PartialD;y}=\frac{z3+z7-z1-z9}{4\mathrm{\&Delta;x\&Delta;y}},---\left(4\right)$

$d=\frac{1}{2}\frac{\&PartialD;z}{\&PartialD;x}=\frac{z3+z6+z9-z1-z4-z7}{6\mathrm{\&Delta;x}},---\left(5\right)$

$e=\frac{1}{2}\frac{\&PartialD;z}{\&PartialD;y}=\frac{\mathrm{<figure-callout\; id="1"\; label="z"\; filenames="HSA00000123211900031.png,HSA00000123211900041.png"\; state="\{\{state\}\}">z</figure-callout>}1+\mathrm{<figure-callout\; id="2"\; label="z"\; filenames="HSA00000123211900041.png"\; state="\{\{state\}\}">z</figure-callout>}2+z3-z7-z8-z9}{6\mathrm{\&Delta;y}},---\left(6\right)$

$f=\frac{2(z2+z4+z6+z8)-(z1+z3+z7+z9)+5z5}{9},---\left(7\right)$

Maximum positive curvature computing formula is

K _pos＝(a+b)+[(a-b) ²+c ²] ^1/2. (8)

Wherein, z1 is the value (as shown in Figure 1) of curved surface at the grid node place to z9, and Δ x and Δ y are the x direction of grid node and the interval on the y direction, by the coefficient a in the equation (1), and b, c, d, e, f can calculate the curvature of current sampling point.Also can adopt big grid node, for example 5 * 5, but this need separate bigger system of equations, needs more calculated amount.Grid cell is few more, and the calculating of curvature is fast more, and the local feature of the curvature that obtains is obvious more.Above-mentioned curvature of face computing method those skilled in the art can be with reference to following document: Robert A.Curvature attributes and their application to 3D interpreted horizons.First Break, 2001,19 (2): 85～100.

The computing method of body curvature of the prior art below are described.Al-Dossary and Marfurt (2006) propose computing method (the Al-Dossary S.and Marfurt K.J.3D volumetricmultispectral estimates of reflector curvature and rotation.Geophysics of body curvature, 2006,71 (5): 41～51), this method is calculated inclination angle, reflection horizon and bearing data earlier, obtains body curvature through rotation of coordinate with different differentiate formula again.This method specifically comprises:

Suppose to have the inclination data body:

$u(z,x,y)=p(z,x,y)\hat{X}+q(z,x,y)\hat{Y},---\left(9\right)$

Wherein p and q are the components at inclination angle,

With Be the vector of unit length in the Cartesian coordinate system, Al-Dossary and Marfurt are decomposed into two parts with it from mathematics

$\mathrm{div}\left(u\right)=\&PartialD;p/\&PartialD;x+\&PartialD;q/\&PartialD;y,---\left(10\right)$

With

$r_z=\&PartialD;p/\&PartialD;y-\&PartialD;q/\&PartialD;x,---\left(11\right)$

Use the reflection inclination data p and the q of input, at the x=y=0 place, the coefficient in the equation (1) becomes:

D _x(D _xz)＝2a＝D _xp， (12)

D _y(D _yz)＝2b＝D _yq， (13)

D _xq+D _yp＝2c， (14)

Wherein, operator D _xAnd D _yCan be equation (15), differentiate formula arbitrarily in (16) and (18).By selecting coordinate system, obtain the coefficient in the equation (1), calculate body curvature.

Adoptable differentiate formula is as follows:

$\frac{\mathrm{du}\left(x\right)}{\mathrm{dx}}=\frac{\frac{2}{3}u(x+h)-\frac{2}{3}u(x-h)-\frac{1}{12}u(x+2h)+\frac{1}{12}u(x-2h)}{h},---\left(15\right)$

$D23=\frac{1}{2\mathrm{\&mu;}}[F(-x)-F\left(x\right)],---\left(16\right)$

Wherein F (x) is:

$\frac{F\left(x\right)}{A}=\cos \mathrm{\&phi;exp}(-\frac{\mathrm{\&alpha;x}}{\mathrm{\&eta;\&mu;}})-\cos (\frac{\mathrm{\&gamma;x}}{\mathrm{\&eta;\&mu;}}+\mathrm{\&phi;})\exp (-\frac{\mathrm{\&beta;x}}{\mathrm{\&eta;\&mu;}})+\frac{\mathrm{\&eta;}}{\mathrm{\&alpha;}}C\left(\mathrm{\&phi;}\right)\sin \left(\frac{x}{\mathrm{\&mu;}}\right)\exp (-\frac{x}{\mathrm{\&mu;}}).---\left(17\right)$

$F(\&PartialD;u/\&PartialD;x)=-i{k}_{x}F\left[u\left(x\right)\right].---\left(18\right)$

The inventor finds that prior art has the following disadvantages at least in realizing process of the present invention:

For the curvature of face computing method of prior art one, need before carrying out curvature calculating, know layer position (curved surface) data earlier, the precision of this method is subjected to the influence of the picking errors in the interpretation process of layer position easily.

Body curvature computing method for prior art two, the body curvature method that Al-Dossary and Marfurt propose is too complicated: at first need to calculate the inclination data body, utilize the inclination data body to obtain coefficient in the equation (1) then by differentiate, after again coordinate system suitably being rotated, just obtain body curvature at last.But for yardstick decomposes is not a direct method.

To sum up, it is lower to adopt these two kinds of computing method to discern the precision of discontinuums such as subterranean fracture.

Summary of the invention

The purpose of the embodiment of the invention is, a kind of detection method and device of underground discontinuum is provided, to improve recognition capability and the detection efficiency to underground discontinuum.

On the one hand, the embodiment of the invention provides a kind of detection method of underground discontinuum, and described method comprises: obtain geological data; According to the tectonic structure of described geological data correspondence, choose the size of computation window; Access time, length was greater than the time range of described computation window; The computation window of utilize coherent algorithm, choosing in advance reaches the time range of choosing in advance, obtains the curvature of each testing data point in the described geological data; According to the curvature of each testing data point, obtain the body curvature of whole geological data; Body curvature according to described whole geological data detects underground crack and tomography.

On the other hand, the embodiment of the invention provides a kind of pick-up unit of underground discontinuum, and described device comprises: data capture unit is used to obtain geological data; The time window choose the unit, be used for tectonic structure according to described geological data correspondence, choose the size of computation window; Time range is chosen the unit, is used for access time length greater than the time range of described computation window; The curvature acquiring unit is used to computation window that utilizes coherent algorithm, chooses in advance and the time range of choosing in advance, obtains the curvature of each testing data point in the described geological data; Body curvature acquiring unit is used for the curvature according to each testing data point, obtains the body curvature of whole geological data; The discontinuum detecting unit is used for the body curvature according to described whole geological data, and underground crack and tomography are detected.

The technique scheme of the confession of the embodiment of the invention, having overcome prior art calculates the influence that the middle level bit data picked up and causes the precision to reduce in the curvature of face, and the deficiency of aspects such as body curvature computation process is loaded down with trivial details, complexity, make the computation process of body curvature more straightforward, error is little, precision is high, counting yield is high, thereby can greatly improve the recognition capability and the detection efficiency of discontinuum.

Description of drawings

In order to be illustrated more clearly in the embodiment of the invention or technical scheme of the prior art, to do one to the accompanying drawing of required use in embodiment or the description of the Prior Art below introduces simply, apparently, accompanying drawing in describing below only is some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain other accompanying drawing according to these accompanying drawings.

Fig. 1 is 3 * 3 grid method synoptic diagram of prior art;

Fig. 2 is the overall flow figure of the detection method of 1 one kinds of underground discontinuums of the embodiment of the invention;

Fig. 3 is a kind of particular flow sheet of the embodiment of the invention 1 according to method shown in Figure 2;

Fig. 3 A is the 3D seismic data space synoptic diagram of the embodiment of the invention 1;

Fig. 4 is the 2D section of the generated data of the embodiment of the invention 1;

Fig. 5 is the maximum positive curvature attribute of 2D section among Fig. 4 of the embodiment of the invention 1;

Fig. 6 is the process flow diagram of detection method of a kind of underground discontinuum of the embodiment of the invention 2;

Fig. 7 is a kind of particular flow sheet of the embodiment of the invention 2 according to method shown in Figure 6;

Fig. 8 is the 2D section of the actual geological data of the embodiment of the invention 2;

Fig. 9 is the coherence slice of the actual geological data of the embodiment of the invention 2;

Figure 10 is the maximum positive curvature section of full range band of the actual geological data of the embodiment of the invention 2;

Figure 11 is the maximum positive curvature section of low frequency of the actual geological data of the embodiment of the invention 2;

Figure 12 is the maximum positive curvature section of high frequency of the actual geological data of the embodiment of the invention 2;

Figure 13 is the functional block diagram of the pick-up unit of 3 one kinds of underground discontinuums of the embodiment of the invention;

Figure 13 A is the detailed functional block diagram of the embodiment of the invention 3 mean curvature acquiring units 140;

Figure 13 B is the detailed functional block diagram of the search subelement 141 of the embodiment of the invention 3;

Figure 14 is the concrete function block diagram of the pick-up unit of the embodiment of the invention 3 another kind of underground discontinuums;

Figure 14 A is the concrete function block diagram of the embodiment of the invention 3 multiband data capture units 240;

Figure 14 B is the concrete function block diagram of the embodiment of the invention 3 multiple dimensioned body curvature computing units 250.

Embodiment

For the purpose, technical scheme and the advantage that make the embodiment of the invention clearer, below in conjunction with the accompanying drawing in the embodiment of the invention, technical scheme in the embodiment of the invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, rather than whole embodiment.Based on the embodiment among the present invention, those of ordinary skills belong to the scope of protection of the invention not making the every other embodiment that is obtained under the creative work prerequisite.

The method of the embodiment of the invention is based on following principle: suppose that subsurface reflective boundary is continuous, therefore the seismic reflection lineups of same reflecting interface are continuous, similar.In seismic data volume, the reflection line-ups at interface (events) forms a continuous curved surface, and the similarity of lineups can be calculated by coherent algorithm.In one embodiment, for the current sampling point that will calculate, at first adopt coherent algorithm to ask for a local curved surface, utilize this local curved surface to calculate the curvature of this sampled point then, according to the curvature of each data point, obtain the body curvature of whole geological data then.In another embodiment, geological data is decomposed, and, obtain the body curvature on the different scale in conjunction with the body curvature computing method among the last embodiment by application window Fourier Fourier conversion (WFT).

Embodiment 1:

The embodiment of the invention 1 provides a kind of detection method of underground discontinuum.Fig. 2 is the process flow diagram of detection method of a kind of underground discontinuum of the embodiment of the invention 1.As shown in Figure 2, this method comprises the steps:

S201, obtain geological data;

S202, according to the tectonic structure of described geological data correspondence, choose the size of computation window;

Particularly, the detailed process of S202 can comprise: when the tectonic structure on the earthquake section is complicated, choose little computation window; Perhaps, when the tectonic structure on the earthquake section is simple, window when choosing computation.The basis for estimation of tectonic structure complexity is: degree and factors such as quantity, fold morphology and formation contact according to the mature fault on the seismic section are judged.For example can choose computation window based on following principle: when the Geological Structural Forms on the seismic section of described geological data only comprises parallel laminar (parallelstratified), window when choosing computation; Perhaps, Geological Structural Forms on the seismic section of described geological data is except that comprising parallel laminar, when also comprising anticline (anticline), at least a in tiltedly (syncline), tomography (fault), unconformability (unconformity) or fold (fold), choose little computation window.For example, little computation window can be the time window less than 5 sampled points, and the time window that comprises 5 or above sampled point window when being computation, but the embodiment of the invention is not restricted to this.

S203, access time, length was greater than the default time range of described computation window;

S204, utilize coherent algorithm, the computation window chosen in advance and the time range chosen in advance, obtain the curvature of each testing data point in the described geological data;

Particularly, when obtaining the curvature of a testing data point, the process of step S204 specifically can comprise: for each the testing data point in the geological data, road with described testing data point place is the center, on each phase neighboring trace, in described time range, search for, seek on described each phase neighboring trace based on coherent algorithm and put maximally related data point with described testing data according to described computation window; According to putting maximally related a plurality of data point and local curved surface of described testing data point match with described testing data on a plurality of phase neighboring traces; Calculate the curvature of described testing data point according to described local curved surface.

Particularly, when on seeking a phase neighboring trace, putting maximally related data point with described testing data, road with described testing data point place is the center, on described phase neighboring trace, in described time range, search for, seek on the described phase neighboring trace process of putting maximally related data point with described testing data based on coherent algorithm and specifically comprise the steps: according to described computation window

On described phase neighboring trace, in described time range, search for according to computation window, obtain a plurality of coherent value based on coherent algorithm, from described a plurality of coherent value, select maximum coherent value, and the data point of the coherent value correspondence of described maximum is defined as the relevant data points of described testing data point on described phase neighboring trace.

When obtaining a coherent value based on coherent algorithm, detailed process comprises the steps:

Choose first data acquisition in the computation window on the road at described testing data point place; In the described time range of described phase neighboring trace, choose upper second data acquisition in computation window of described phase neighboring trace; Described first data acquisition and described second data acquisition are concerned with calculating to obtain a coherent value based on coherent algorithm; After executing once relevant calculating on the described phase neighboring trace, the computation window on the described phase neighboring trace is moved along predetermined direction in described time range.

All data in the described time range have been traveled through repeating the computation window of aforesaid operations on described phase neighboring trace on the described phase neighboring trace: will obtain a plurality of coherent value after repeating aforesaid operations;

From described a plurality of coherent value, select maximum coherent value, and the data point of the coherent value correspondence of described maximum is defined as the relevant data points of described testing data point on described phase neighboring trace.

Alternatively, described predetermined direction is the direction that increases along the quantity of data point, perhaps the direction that reduces along the quantity of data point.When described predetermined direction is during along direction that the quantity of data point increases, the starting point of described second computation window is arranged at the lower limit place of described default time range; Perhaps, when described predetermined direction is during along direction that the quantity of data point reduces, the starting point of described second computation window is positioned at the upper vault of described default time range.

Particularly, according to putting in the process of a maximally related a plurality of data point and a local curved surface of described testing data point match on a plurality of phase neighboring traces with described testing data, if adopt " 3 * 3 " gridding method, this local curved surface is to put the local curved surface that maximally related 8 data points constitute by described testing data point and with described testing data, certainly, the embodiment of the invention is not restricted to this, as adopting " 5 * 5 " gridding method, then this local curved surface is to put the local curved surface that maximally related 24 data points constitute by described testing data point and with described testing data.Adopt the grid node of other scale can realize equally.In the process of the curvature of calculating described testing data point according to described local curved surface, the method for calculating the curvature of testing data point can see also formula (1) to formula (8).

S205, according to the curvature of each data point, obtain the body curvature of whole geological data;

Particularly, the detailed process of S205 can comprise: according to the curvature of each the testing data point on the road at described testing data point place, obtain the curvature in this road; According to the curvature in each road, obtain the body curvature of whole geological data.Described body curvature comprises maximum positive curvature.The embodiment of the invention is not restricted to this, and the detailed process of S205 also can comprise: calculate the curvature of current data point, like this each point of the whole data volume of cycle calculations is just entirely shaken the body curvature of data.

S206, according to the body curvature of described whole geological data, underground crack and tomography are detected.

Particularly, the place that body curvature absolute value is big is corresponding to underground crack and tomography.Can detect or discern underground crack and tomography according to this principle.

Particularly, in the flow process shown in Figure 2 based on coherent algorithm can be first generation coherent body computational algorithm, second generation coherent body computational algorithm or third generation coherent body computational algorithm.Below further specify the coherent algorithm of using in the embodiment of the invention method shown in Figure 2.

The embodiment of the invention has adopted coherent algorithm to seek local curved surface in the algorithm that calculates curvature.Relevant is the mathematical measure of two signal similar degree, has three kinds of methods to calculate coherent body at present.Preferably adopt third generation coherent body technique C3 method in the embodiment of the invention.C3 coherent algorithm those skilled in the art can list of references: Gersztenkorn A.and Marfurt K.J.Eigenstructure-basedcoherence computations as an aid to 3D structural and stratigraphic mapping.Geophysics, 1999,64 (5): 1446～1479.

The detailed description of C3 coherent algorithm below is provided.In the 3D geological data, the window data body that analyze is a little 3D cube, comprises the J road, and per pass has N sampled point.The amplitude of each sampled point in this cube is arranged in a matrix D in order, and wherein the index of sampled point is N, and the index of Taoist monastic name is J:

Its covariance matrix is E=D ^TD.The mark of covariance matrix E is:

$\mathrm{Tr}\left(E\right)=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{e}_{\mathrm{jj}}=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\mathrm{\&lambda;}}_j---\left(20\right)$

λ j (j=1 ..., J) be the eigenwert of covariance matrix E.Relevant being expressed as:

$C=\frac{{\mathrm{\&lambda;}}_1}{\mathrm{Tr}\left(E\right)},---\left(21\right)$

Wherein, λ ₁It is the eigenvalue of maximum of covariance matrix E.

Method shown in Figure 2 below is described in further detail by way of example.

Fig. 3 is a kind of particular flow sheet of the embodiment of the invention 1 according to method shown in Figure 2.As shown in Figure 3, this method comprises the steps:

S301, at the architectonic complexity on the seismic section, choose suitable computation window size; For complex structure, choose little computation window, for simple structure, window when choosing computation.

Particularly, computation window is meant that to need selected time range, its purposes be the sampling number certificate in the window in order to obtain in order to calculate.

S302, for the current sampling point that calculates, be the center of computation window; Road with the current sampling point place is the center, chooses the default time range on the phase neighboring trace, adopts coherent algorithm to calculate on the phase neighboring trace and the maximally related point of current sampling point.

The implication in road is meant the geological data that writes down on the single acceptance point.This road is meant the current road that calculates.Because the road of geological data distributes as net shape on the plane, the phase neighboring trace is the road with the value difference 1 of the main profile in this road and/or cross-track.

S303, on the phase neighboring trace repeating step S302, obtain current 8 maximally related points on the phase neighboring trace; Utilize these 9 points to come fit equation (1), obtain each coefficient (for example a, b, c, d, e, f) of equation (1), calculate the curvature value of current point then according to each coefficient of obtaining and equation (8) based on " 3 * 3 " gridding method.

S304, on whole road repeating step S302 and step S303, obtain the result of calculation in this road.

The process of S305, repeating step S304 is obtained the result of calculation in each road, to obtain the curvature result of whole data volume.

Below describe among above-mentioned steps S302 and the S303 search procedure in detail based on coherent algorithm:

Fig. 3 A is the 3D seismic data space synoptic diagram of the embodiment of the invention 1.Please, suppose that 3D seismic data has L bar survey line, the M road is arranged on the every survey line, have N data point on the per pass, form a 3D seismic data space representation and be in conjunction with consulting Fig. 3 A: L * M * N, wherein L, M and N are positive integer; Suppose the current location tables of data point in described 3D seismic data space to be calculated be shown (i, j, k), the location tables of road, current data point place (hereinafter to be referred as current road) in described 3D seismic data space be shown (i, j), a plurality of phase neighboring traces of road, current data point place correspondence are expressed as (i-1 respectively in the position in 3D seismic data space, j-1), (i-1, j), (i-1, j+1), (i, j-1), (i, j+1), (i+1, j-1), (i+1, j), (i+1, j+1), wherein, i ∈ [1, L], j ∈ [1, M], k ∈ [1, N].

The lengths table of supposing the computation window chosen is shown l, and then described current data point k is set to the center of described computation window l, and chooses described current road (i, the data of the computation window l that to be positioned at described current data point k on j) be the center.

Preset Time scope on each phase neighboring trace is set at [k-r, k+r], and chooses computation window on each phase neighboring trace, its lengths table is shown l, wherein, and r 〉=l/2;

Respectively at the described Preset Time scope [k-r of each phase neighboring trace, k+r] in, choose the phase neighboring trace upper in computation window l data and described current road (i, that chooses on j) is positioned at computation window l data, based on the coherent algorithm calculating that is concerned with, obtain a coherent value.Wherein, the starting point of the computation window l of described phase neighboring trace is from k-r, after the once relevant calculating of every execution, make computation window l at Preset Time scope [k-r, k+r] interior data point of starting point increase, and in whole Preset Time scope [k-r, k+r] the repeatedly relevant calculation process of interior circulation execution, i.e. continuous window l during mobile computing on each phase neighboring trace, and choose on the phase neighboring trace on data in the computation window l after moving and the current road calculating that is concerned with of the data in the computation window l, all interior data points of described whole Preset Time scope [k-r, k+r] have been traveled through until described computation window l.

On each phase neighboring trace, carry out repeatedly relevant calculating obtaining a plurality of coherent value, and choose coherent value corresponding sampling points maximum in described a plurality of coherent value, as on the described phase neighboring trace with the maximally related data point of described current data point.

Particularly, if find maximum coherence value on the phase neighboring trace, then with the center of the computation window of this maximum coherence value correspondence as maximally related data point.

Further, the proof procedure that the embodiment of the invention 1 is carried out body curvature Calculation Method based on coherent algorithm below is described.

Fig. 4 is the 2D section of the generated data of the embodiment of the invention 1.Generated data is the data that theoretical formula calculates, and geological data in general sense is meant the data that construction collects.As shown in Figure 4, in this synthetic curved surface geological data, curved surface left-half shape is the arched roof of a projection, the right half part shape be one recessed hollow, calculate the body curvature (maximum positive curvature) of these earthquake data with the coherent algorithm of the embodiment of the invention.Two dimensional cross-section such as Fig. 4 along arched roof summit and hollow low spot, the maximum positive curvature attribute of this two dimensional cross-section as shown in Figure 5, Fig. 5 is the maximum positive curvature attribute of 2D section among Fig. 4 of the embodiment of the invention 1, as shown in Figure 5, the value of maximum positive curvature attribute has been portrayed the morphological feature of this two dimensional cross-section exactly: on the occasion of the expression protuberance, negative value is represented depression, and the curvature attribute just/extreme point of negative value is corresponding consistent with the arched roof summit of structure and hollow low spot, has verified embodiment of the invention algorithm validity thus.

In embodiment of the invention Fig. 2 and the method shown in Figure 3, use coherent algorithm can effectively calculate the body curvature of geological data, the algorithm explicit physical meaning, principle is simple, and travelling speed is fast, can effectively improve recognition capability and detection efficiency to underground discontinuum.

Embodiment 2:

The embodiment of the invention 2 provides the detection method of another kind of underground discontinuum.Before carrying out detection, calculate multiple dimensioned earthquake curvature attribute earlier, this process is mainly by application window Fourier conversion, geological data is decomposed into different frequency bands, then the geological data behind the frequency division is calculated its corresponding body curvature respectively, to obtain the body curvature of different scale.In the process of calculating body curvature, also be based on coherent algorithm.

Fig. 6 is the process flow diagram of detection method of a kind of underground discontinuum of the embodiment of the invention 2.As shown in Figure 6, this method comprises the steps:

S601, obtain geological data;

S602, described geological data is carried out spectrum analysis, determine the frequency band range of described geological data;

S603, according to the frequency band range of described geological data, application window Fourier Fourier conversion is decomposed into geological data in the different sub-band scopes with described geological data;

Particularly, by data are carried out spectrum analysis, obtain the whole frequency band scope of data at S602.In the process of S603, can be as required, with the whole frequency band scope division is at least two different sub-band scopes (promptly littler frequency band range), thereby be selected first input parameter of window Fourier conversion, be about to first input parameter that described at least two different sub-band scopes are chosen to be window Fourier conversion; Then the described geological data that obtains is chosen to be second input parameter of window Fourier conversion; Then, according to described first input parameter and described second input parameter, obtain to correspond respectively to the geological data of at least two different sub-band scopes based on window Fourier conversion.After carrying out said process, the output result of window Fourier conversion is the geological data corresponding to the different frequency bands parameter, and whole thus geological data is broken down into the geological data in the different frequency band ranges.

For example, one group of frequency band parameter comprises the high frequency value and the low frequency value of frequency band, again for example, by given many group frequency band parameters, geological data can be decomposed into geological data in the different frequency bands by window Fourier conversion.

S604, calculate the body curvature of the geological data correspondence of each frequency band respectively based on coherent algorithm, to obtain the body curvature of different scale.

Particularly, the detailed process of S604 can comprise: carry out curvature respectively based on the geological data of coherent algorithm after to frequency division and calculate.Computation process specifically can comprise: adopt coherent algorithm to obtain 9 points (as shown in Figure 1) on the local curved surface, calculate the curvature of current point, this process of iteration can obtain multiple dimensioned curvature attribute.And please those skilled in the art's combination consult the method shown in the embodiment of the invention 1.

Fig. 7 is a kind of particular flow sheet of the embodiment of the invention 2 according to method shown in Figure 6, and as shown in Figure 7, this method comprises the steps:

S701, obtain geological data; The same S601 of this step repeats no more.

S702, described geological data is carried out spectrum analysis, to determine the frequency range of described geological data; The same S602 of this step repeats no more.

S703, judge whether to carry out frequency division, when judged result enters step S704 for need carry out frequency division the time, otherwise, enter step S706.

S704, according to the frequency range of described geological data, application window Fourier Fourier conversion is decomposed into different frequency bands with described geological data; The same S603 of this step repeats no more, and enters S705 after the S704.

S705, calculate the body curvature of the geological data correspondence of each frequency band respectively based on coherent algorithm, to obtain the body curvature of different scale.The same S604 of this step repeats no more.

Particularly, the size of yardstick image study object.The geological data of the corresponding high frequency band of the curvature attribute of small scale, vice versa.

S706, calculate the body curvature of the geological data correspondence of full range band based on coherent algorithm, to obtain full range belt body curvature.Particularly, the method for institute's foundation is also based on coherent algorithm in the computation process of full range belt body curvature.

Alternatively, S703 not necessarily also can both calculate multiple dimensioned body curvature according to concrete research needs, also calculated full range belt body curvature.The purpose of calculating full range belt body curvature is to compare with the curvature of other yardsticks, can provide more fully information for seismic interpretation.

Below describe further by a concrete example and real data to be carried out multiple dimensioned body curvature is calculated and the processing procedure of full range belt body curvature calculating.

At first, calculate maximum positive curvature (body curvature) attribute of this earthquake data full range band.

Fig. 8 is a two-dimension earthquake longitudinal profile of the embodiment of the invention 2, and the straight line of locating in 1.8 seconds has marked the time location (1.8s) of the time slice of being got (Fig. 9 to Figure 12), longitudinal axis express time among Fig. 8, and transverse axis is represented Taoist monastic name.Section from Fig. 8 as can be seen, located to have passed through in 1.8 seconds two shatter belts (as shown in elliptical section among Fig. 8 is divided) and tomography (as square frame among Fig. 8 partly shown in).As everyone knows, coherence properties can be discerned the discontinuum in the geological data well, as tomography and crack etc.Fig. 9 is the horizontal coherence slice of the actual geological data of the embodiment of the invention 2, and as shown in Figure 9, coherence properties section has clearly reflected the information of two shatter belts and tomography.Figure 10 is the maximum positive curvature section of full range band of the actual geological data of the embodiment of the invention 2, and as shown in figure 10, maximum positive curvature attribute section has also clearly reflected the information of two shatter belts and tomography.This explanation curvature attribute is the same with coherence properties, can be used for discerning discontinuums such as tomography and crack, and good correlativity is arranged between curvature attribute and fold, tomography, anticline and the shatter belt.

Then, geological data being carried out multiple dimensioned body curvature calculates.

Because the frequency band of this geological data is not wide, only its (geological data) is decomposed into two data volumes of high and low frequency with window Fourier conversion, these two data volumes are calculated maximum positive curvature (body curvature) attribute respectively.The low-frequency band of geological data is geology result's a background, and the high frequency band of geological data has reflected little architectonic feature.Therefore the curvature (as Figure 11 and shown in Figure 12) of different scale has reflected the different characteristic of geologic structure.Wherein, Figure 11 is the maximum positive curvature section of low frequency of the actual geological data of the embodiment of the invention 2; Figure 12 is the maximum positive curvature section of high frequency of the actual geological data of the embodiment of the invention 2.

At last, the curvature of the curvature of the curvature of comparative analysis full range band data, low-frequency band data, high frequency band data.

Particularly, compare with the curvature (as shown in figure 10) of full range band data, the curvature of low-frequency data (as shown in figure 11) has comprised less information, only portrayed architectonic background characteristics, for example, rectangular area among contrast Figure 10 and Figure 11, information is few in the curvature of Figure 11 medium and low frequency data as can be seen, has only portrayed architectonic background characteristics; The curvature of high-frequency data (as shown in figure 12) can demonstrate more little architectonic feature, therefore can reflect little geologic feature better, for example, the rectangular area among contrast Figure 10 and Figure 12, the curvature of Figure 12 medium-high frequency data can show more little geologic body feature as can be seen.Confirmed that thus the multi-dimension curvature attribute can disclose different geologic features, can improve the resolution of explaining, disclosed the more information in the data.

The curvature of face computing method based on layer bit data of prior art, the information (curved surface data) of necessary known layer bit data before calculating, promptly use the layer bit data that provides as input, the precision of this method is subjected to the influence of the picking errors in the interpretation process of layer position easily.And the method for the embodiment of the invention is for being input as whole geological data, by finding a local curved surface based on coherent algorithm in seismic data volume, and utilizes the data of this local curved surface to carry out curvature and calculates, and precision is higher, error is littler.And the embodiment of the invention is towards the multiple dimensioned earthquake curvature attribute acquisition methods based on coherent algorithm of seismic data interpretation, and the use coherent algorithm can effectively calculate the body curvature of geological data, the algorithm explicit physical meaning, and principle is simple, and travelling speed is fast.And application window Fourier conversion can obtain multiple dimensioned body curvature attribute, and this new multi-dimension curvature attribute can greatly improve the recognition capability of geologic bodies such as crack.

Embodiment 3:

The embodiment of the invention also provides a kind of pick-up unit of underground discontinuum.

Figure 13 is the functional block diagram of pick-up unit of a kind of underground discontinuum of the embodiment of the invention 3.As Figure 13 institute method, this device comprises:

Data capture unit 110 is used to obtain geological data;

The time window choose unit 120, be used for tectonic structure according to described geological data correspondence, choose the size of computation window;

Particularly, window selected cell 120 when described specifically can be used for when the Geological Structural Forms on the seismic section of described geological data only comprises parallel laminar window when choosing computation; Perhaps, the Geological Structural Forms on the seismic section of described geological data when also comprising anticline, at least a in oblique, tomography, unconformability or fold, is chosen little computation window except that comprising parallel laminar.

Time range is chosen unit 130, is used for access time length greater than the time range of described computation window;

Curvature acquiring unit 140 is used to computation window that utilizes coherent algorithm, chooses in advance and the time range of choosing in advance, obtains the curvature of each testing data point in the described geological data;

Body curvature acquiring unit 150 is used for the curvature according to each testing data point, obtains the body curvature of whole geological data;

Particularly, described body curvature acquiring unit 150 specifically can be used for the curvature according to each the testing data point on the road at described testing data point place, obtains the curvature in this road; And, obtain the body curvature of whole geological data according to the curvature in each road.

Discontinuum detecting unit 160 is used for the body curvature according to described whole geological data, and underground crack and tomography are detected.

Figure 13 A is the detailed functional block diagram of the curvature acquiring unit 140 of the embodiment of the invention 3.As shown in FIG. 13A, curvature acquiring unit 140 comprises: search subelement 141, surface fitting subelement 142 and curvature computation subunit 143;

Described search subelement 141, the road that is used for described testing data point place is the center, on each phase neighboring trace, in described time range, search for, seek on described each phase neighboring trace based on coherent algorithm and put maximally related data point with described testing data according to described computation window;

Described surface fitting subelement 142 is used for according to putting maximally related a plurality of data point and local curved surface of described testing data point match with described testing data on a plurality of phase neighboring traces;

Described curvature computation subunit 143 is used for calculating according to described local curved surface the curvature of described testing data point.

Alternatively, described local curved surface comprises: put the local curved surface that maximally related 8,15 or 24 data points constitute by described testing data point and with described testing data.Described coherent algorithm comprises the first generation, the second generation or third generation coherent body computational algorithm.Described body curvature comprises maximum positive curvature.

Particularly, described search subelement 141, specifically can be used on described phase neighboring trace, in described time range, searching for according to computation window, obtain a plurality of coherent value based on coherent algorithm, from described a plurality of coherent value, select maximum coherent value, and the data point of the coherent value correspondence of described maximum is defined as the relevant data points of described testing data point on described phase neighboring trace.

Figure 13 B is the detailed functional block diagram of the search subelement 141 of the embodiment of the invention 3.Shown in Figure 13 B, search subelement 141 comprises: first data acquisition is chosen module 1411, second data acquisition is chosen module 1412, relevant computing module 1413 and computation window mobile module 1414;

Described first data acquisition is chosen module 1411, is used to choose first data acquisition in the computation window on the road at described testing data point place;

Described second data acquisition is chosen module 1412, is used in the described time range of described phase neighboring trace, chooses upper second data acquisition in computation window of described phase neighboring trace;

Described relevant computing module 1413 is used for described first data acquisition and described second data acquisition are concerned with calculating to obtain a coherent value based on coherent algorithm;

Described computation window mobile module 1414 is used for after executing once relevant calculating on the described phase neighboring trace computation window on the described phase neighboring trace being moved along predetermined direction in described time range.

Alternatively, described predetermined direction is the direction that increases along the quantity of data point, perhaps the direction that reduces along the quantity of data point.When described predetermined direction is during along direction that the quantity of data point increases, the starting point of described second computation window is arranged at the lower limit place of described default time range; Perhaps, when described predetermined direction is during along direction that the quantity of data point reduces, the starting point of described second computation window is positioned at the upper vault of described default time range.

Particularly, when on seeking a phase neighboring trace, putting maximally related data point with described testing data, first data acquisition is chosen module 1411, second data acquisition is chosen module 1412, relevant computing module 1413 and computation window mobile module 1414 and will be repeated multi-pass operations all data in computation window has traveled through described time range, will obtain a plurality of coherent value after above-mentioned functions unit 1411-1414 repeats repeatedly aforesaid operations;

Then, described search subelement 141 is selected maximum coherent value from described a plurality of coherent value, and the data point of the coherent value correspondence of described maximum is defined as the relevant data points of described testing data point on described phase neighboring trace.

Figure 14 is the concrete function block diagram of pick-up unit of the another kind of underground discontinuum of the embodiment of the invention 3.As shown in figure 14, this device 20 comprises:

Data capture unit 210 is used to obtain geological data; Its function is with data capture unit 110.

Spectral analysis unit 220 is used for described geological data is carried out spectrum analysis, determines the frequency band range of described geological data;

Multiband data capture unit 240 is used for the frequency band range according to described geological data, and application window Fourier Fourier conversion is decomposed into geological data in the different sub-band scopes with described geological data; Multiple dimensioned body curvature computing unit 250 is used for calculating respectively based on coherent algorithm the body curvature of the geological data correspondence of each sub-band scope, to obtain the body curvature of different scale.

Discontinuum detecting unit 270 is used for according to the body curvature that obtains different scale underground crack and tomography being detected.

Figure 14 A is the concrete function block diagram of the multiband data capture unit 240 of the embodiment of the invention 3.Shown in Figure 14 A, described multiband data capture unit 240 can comprise that band decomposition subelement 241, first input parameter are chosen subelement 242, second input parameter is chosen subelement 243 and window Fourier varitron unit 244;

Described band decomposition subelement 241 can be used for the frequency band range of described geological data is divided at least two different sub-band scopes; Described first input parameter is chosen subelement 242, can be used for described at least two different sub-band scopes are chosen to be first input parameter of window Fourier conversion; Described second input parameter is chosen subelement 243, can be used for the described geological data that obtains is chosen to be second input parameter of window Fourier conversion; Described window Fourier varitron unit 244 can be used for according to described first input parameter and described second input parameter, obtains to correspond respectively to the geological data of at least two different sub-band scopes based on window Fourier conversion.

Alternatively, pick-up unit 20 can further include:

Frequency division judging unit 230 is connected with band decomposition unit 240 with spectral analysis unit 220, is used to judge whether geological data to be decomposed into different frequency ranges.And,

Full range belt body curvature computing unit 260, be connected with discontinuous detecting unit 270 with frequency division judging unit 230, be used for when the judged result of frequency division judging unit 230 when not needing frequency division, calculate the body curvature of the geological data correspondence of full range band based on coherent algorithm, to obtain full range belt body curvature.

Above-mentioned discontinuum detecting unit 270 can also be used for according to the result of calculation of multiple dimensioned body curvature computing unit 250 and full range belt body curvature computing unit 260 underground discontinuum being detected.

Particularly, the function of full range belt body curvature computing unit 260 can by time window selected cell 120, time range shown in Figure 13 choose unit 130, curvature acquiring unit 140 and body curvature acquiring unit 150 common intact with.

Figure 14 B is the concrete function block diagram of multiple dimensioned body curvature computing unit 250 among the embodiment of the invention 3 Figure 14.As shown in Figure 14B, multiple dimensioned body curvature computing unit 250 can further comprise the first yardstick body curvature computing unit 251, the second yardstick body curvature computing unit 252 and the 3rd yardstick body curvature computing unit 253, is respectively applied for the first body curvature, the second body curvature and the trisome curvature of the data correspondence of calculating first frequency band, second frequency band and the 3rd frequency band.The function of the wherein above-mentioned first yardstick body curvature computing unit 251, the second yardstick body curvature computing unit 252 and the 3rd yardstick body curvature computing unit 253 all can adopt the time window among Figure 13 choose unit 120, time range choose unit 130, curvature acquiring unit 140, and body curvature acquiring unit 150 realize, accordingly, the data of above-mentioned each functional unit 120-150 processing are the data from first frequency band, second frequency band or the 3rd frequency band.

The detailed functions of each functional unit shown in embodiment of the invention Figure 13-Figure 14 B, and the device principle of work in aforesaid method embodiment, describe, so do not repeat them here.

The described pick-up unit of the embodiment of the invention, the use coherent algorithm can effectively calculate the body curvature of geological data, the algorithm explicit physical meaning, principle is simple, and travelling speed is fast.And application window Fourier conversion can obtain multiple dimensioned body curvature attribute, and this new multi-dimension curvature attribute can greatly improve the recognition capability of geologic bodies such as crack.

One of ordinary skill in the art will appreciate that all or part of flow process that realizes in the foregoing description method, be to instruct relevant hardware to finish by computer program, described program can be stored in the computer read/write memory medium, this program can comprise the flow process as the embodiment of above-mentioned each side method when carrying out.Wherein, described storage medium can be magnetic disc, CD, read-only storage memory body (Read-OnlyMemory, ROM) or at random store memory body (Random Access Memory, RAM) etc.

Above embodiment only in order to the technical scheme of the explanation embodiment of the invention, is not intended to limit; Although the embodiment of the invention is had been described in detail with reference to previous embodiment, those of ordinary skill in the art is to be understood that: it still can be made amendment to the technical scheme that aforementioned each embodiment put down in writing, and perhaps part technical characterictic wherein is equal to replacement; And these modifications or replacement do not make the essence of appropriate technical solution break away from the spirit and scope of each embodiment technical scheme of the embodiment of the invention.

Claims (18)

1. the detection method of a underground discontinuum is characterized in that, described method comprises:

Obtain geological data;

According to the tectonic structure of described geological data correspondence, choose the size of computation window;

Access time, length was greater than the time range of described computation window;

The computation window of utilize coherent algorithm, choosing in advance reaches the time range of choosing in advance, obtains the curvature of each testing data point in the described geological data;

According to the curvature of each testing data point, obtain the body curvature of whole geological data;

Body curvature according to described whole geological data detects underground crack and tomography.

2. method according to claim 1 is characterized in that, when obtaining the curvature of a testing data point, comprising:

Road with described testing data point place is the center, searches in described time range according to described computation window on each phase neighboring trace, seeks on described each phase neighboring trace based on coherent algorithm and puts maximally related data point with described testing data;

According to putting maximally related a plurality of data point and local curved surface of described testing data point match with described testing data on a plurality of phase neighboring traces;

Calculate the curvature of described testing data point according to described local curved surface.

3. method according to claim 2 is characterized in that, when putting maximally related data point with described testing data on seeking a phase neighboring trace, specifically comprises:

On described phase neighboring trace, in described time range, search for according to computation window, obtain a plurality of coherent value based on coherent algorithm, from described a plurality of coherent value, select maximum coherent value, and the data point of the coherent value correspondence of described maximum is defined as the relevant data points of described testing data point on described phase neighboring trace.

4. method according to claim 3 is characterized in that, when obtaining a coherent value based on coherent algorithm, specifically comprises:

Choose first data acquisition in the computation window on the road at described testing data point place;

In the described time range of described phase neighboring trace, choose upper second data acquisition in computation window of described phase neighboring trace;

Described first data acquisition and described second data acquisition are concerned with calculating to obtain a coherent value based on coherent algorithm;

After executing once relevant calculating on the described phase neighboring trace, the computation window on the described phase neighboring trace is moved along predetermined direction in described time range.

5. method according to claim 4 is characterized in that, described predetermined direction is the direction that increases along the quantity of data point, perhaps the direction that reduces along the quantity of data point.

6. method according to claim 5 is characterized in that,

When described predetermined direction is during along direction that the quantity of data point increases, the starting point of described second computation window is arranged at the lower limit place of described default time range; Perhaps,

When described predetermined direction is during along direction that the quantity of data point reduces, the starting point of described second computation window is positioned at the upper vault of described default time range.

7. method according to claim 1 is characterized in that, described tectonic structure according to described geological data correspondence is chosen the size of computation window; Comprise:

When the Geological Structural Forms on the seismic section of described geological data only comprises parallel laminar, window when choosing computation; Perhaps,

Geological Structural Forms on the seismic section of described geological data when also comprising anticline, at least a in oblique, tomography, unconformability or fold, is chosen little computation window except that comprising parallel laminar.

8. method according to claim 1 is characterized in that, described obtain geological data after, described method also comprises:

Described geological data is carried out spectrum analysis, determine the frequency band range of described geological data;

According to the frequency band range of described geological data, application window Fourier Fourier conversion is decomposed into geological data in the different sub-band scopes with described geological data;

Calculate the body curvature of the geological data correspondence of each sub-band scope respectively based on coherent algorithm, to obtain the body curvature of different scale.

9. method according to claim 8 is characterized in that, according to the frequency band range of described geological data, application window Fourier Fourier conversion is decomposed into geological data in the different sub-band scopes with described geological data; Comprise:

The frequency band range of described geological data is divided at least two different sub-band scopes;

Described at least two different sub-band scopes are chosen to be first input parameter of window Fourier conversion;

The described geological data that obtains is chosen to be second input parameter of window Fourier conversion;

According to described first input parameter and described second input parameter, obtain to correspond respectively to the geological data of at least two different sub-band scopes based on window Fourier conversion.

10. according to each described method among the claim 1-9, it is characterized in that described curvature according to each data point obtains the body curvature of whole geological data; Comprise:

According to the curvature of each the testing data point on the road at described testing data point place, obtain the curvature in this road; According to the curvature in each road, obtain the body curvature of whole geological data.

11. the pick-up unit of a underground discontinuum is characterized in that, described device comprises:

Data capture unit is used to obtain geological data;

The time window choose the unit, be used for tectonic structure according to described geological data correspondence, choose the size of computation window;

Time range is chosen the unit, is used for access time length greater than the time range of described computation window;

The curvature acquiring unit is used to computation window that utilizes coherent algorithm, chooses in advance and the time range of choosing in advance, obtains the curvature of each testing data point in the described geological data;

Body curvature acquiring unit is used for the curvature according to each testing data point, obtains the body curvature of whole geological data;

The discontinuum detecting unit is used for the body curvature according to described whole geological data, and underground crack and tomography are detected.

12. device according to claim 11 is characterized in that, described curvature acquiring unit comprises search subelement, surface fitting subelement and curvature computation subunit;

Described search subelement, the road that is used for described testing data point place is the center, on each phase neighboring trace, in described time range, search for, seek on described each phase neighboring trace based on coherent algorithm and put maximally related data point with described testing data according to described computation window;

Described surface fitting subelement is used for according to putting maximally related a plurality of data point and local curved surface of described testing data point match with described testing data on a plurality of phase neighboring traces;

Described curvature computation subunit is used for calculating according to described local curved surface the curvature of described testing data point.

13. device according to claim 12 is characterized in that,

Described search subelement, specifically be used on described phase neighboring trace, in described time range, searching for according to computation window, obtain a plurality of coherent value based on coherent algorithm, from described a plurality of coherent value, select maximum coherent value, and the data point of the coherent value correspondence of described maximum is defined as the relevant data points of described testing data point on described phase neighboring trace.

14. device according to claim 13 is characterized in that, described search subelement comprises: first data acquisition is chosen module, second data acquisition is chosen module, relevant computing module and computation window mobile module;

Described first data acquisition is chosen module, is used to choose first data acquisition in the computation window on the road at described testing data point place;

Described second data acquisition is chosen module, is used in the described time range of described phase neighboring trace, chooses upper second data acquisition in computation window of described phase neighboring trace;

Described relevant computing module is used for described first data acquisition and described second data acquisition are concerned with calculating to obtain a coherent value based on coherent algorithm;

Described computation window mobile module is used for after executing once relevant calculating on the described phase neighboring trace computation window on the described phase neighboring trace being moved along predetermined direction in described time range.

15. device according to claim 11 is characterized in that,

Window is chosen the unit when described, specifically is used for when the Geological Structural Forms on the seismic section of described geological data only comprises parallel laminar window when choosing computation; Perhaps, the Geological Structural Forms on the seismic section of described geological data when also comprising anticline, at least a in oblique, tomography, unconformability or fold, is chosen little computation window except that comprising parallel laminar.

16. device according to claim 11 is characterized in that, described device also comprises:

Spectral analysis unit is used for described geological data is carried out spectrum analysis, determines the frequency band range of described geological data;

The multiband data capture unit is used for the frequency band range according to described geological data, and application window Fourier Fourier conversion is decomposed into geological data in the different sub-band scopes with described geological data;

Multiple dimensioned body curvature computing unit is used for calculating respectively based on coherent algorithm the body curvature of the geological data correspondence of each sub-band scope, to obtain the body curvature of different scale.

17. device according to claim 16 is characterized in that,

Described multiband data capture unit comprises that band decomposition subelement, first input parameter are chosen subelement, second input parameter is chosen subelement and window Fourier varitron unit;

Described band decomposition subelement is used for the frequency band range of described geological data is divided at least two different sub-band scopes;

Described first input parameter is chosen subelement, is used for described at least two different sub-band scopes are chosen to be first input parameter of window Fourier conversion;

Described second input parameter is chosen subelement, is used for the described geological data that obtains is chosen to be second input parameter of window Fourier conversion;

Described window Fourier varitron unit is used for according to described first input parameter and described second input parameter, obtains to correspond respectively to the geological data of at least two different sub-band scopes based on window Fourier conversion.

18. according to each described device among the claim 11-17, it is characterized in that,

Described body curvature acquiring unit specifically is used for the curvature according to each the testing data point on the road at described testing data point place, obtains the curvature in this road; And, obtain the body curvature of whole geological data according to the curvature in each road.

Images