CN105572732A - Fracture-developed zone detection method based on gradual increase of attribute change rate - Google Patents

Abstract

The invention provides a fracture-developed zone detection method based on gradual increase of attribute change rate, and belongs to the field of oilfield exploration and development. According to the method, variable-lattice gather relevant correction is conducted on the seismic trace of a post-stack center point, variable-scale high-order azimuth amplitude change rate calculation is conducted on the seismic traces of adjacent points, fractures of certain development intensity and different scales are accumulated to form an effective fracture-developed zone, and eventually, the steps are generalized to the whole work area to get a comprehensive prediction result of the development of different-scale fractures characterized by multi-order gradual increase of maximum-azimuth high-order attribute change rate for the whole seismic data.

Description

The Fractured Zone detection method that a kind of attribute change rate is successively promoted

Technical field

The invention belongs to In Oil Field Exploration And Development field, be specifically related to the Fractured Zone detection method that a kind of attribute change rate is successively promoted to.

Background technology

The rate of change of 3-D seismics amplitude attribute, as a kind of means of non-homogeneous layer description, has become the important technical of carbonate rock area fractured-vuggy reservoir intensive development band and the prediction of karst reservoir growth favo(u)rable target.But mathematically continuous print rate of change calculate due to the discrete sampling adopted of having in being sampled by geological observation net point replace, the step-length that different sampled datas is different from employing, make the calculating of derivative relevant to horizontal step-length, the calculating in the past adopted can not contain the entirety of different scale component, and make the amplitude change rate district of delineation in a lot of effectively zone can not become band distribution in flakes, can not the result of calculation of accurate description different scale, and make that the delineation boundary of the effective development belt in crack is fuzzy, validity is indefinite.

Technology based on attribute (or amplitude) rate of change is all launch on the basis of earthquake poststack seismic data, but traditional method is in adjacent several directions usually, selects one to apply a formula, carries out calculating between the seismic trace of fixing grid section, this is normally incomplete for different scale, different mutative scale FRACTURE PREDICTION of growing rank, result is also limitation, only represents extensive crack and can not obtain the entirety in the crack of different scale from small to large.1) or can only carry out calculating between the seismic trace in adjacent direction the defect that prior art exists is:, and this is normally inaccurate for having the FRACTURE PREDICTION of an orientation; 2) or a kind of high-order calculating formula (as 4 rank) can only be selected, a kind of array (as 5x5 point), a row filter operator calculating fracture development can only be used, the crack of different scale can not be predicted, especially the crack of small scale, also will need the growth roughly yardstick predicting crack.

In sum, the calculating of attribute change rate all launches on poststack seismic data, original continuous print calculates by the restriction of limited observation station by discretize, traditional method usually can only by manually choosing and a kind of fixing difference observing step-length relevant or higher-order wave equation, calculating is carried out between the adjacent seismic trace of fixing dot matrix, mathematically there is suppression in the crack of this treatment technology to small scale (high wave number), do not detect the relative change in different scale especially small scale crack, therefore the achievement detected does not embody the multiple dimensioned property in crack, can only testing result intermittently in a lot of effectively zone, and the entirety that can not detect the different cracks becoming band in flakes grows region, the overall distribution region of plane space fracture development can not accurately be described.

Summary of the invention

The object of the invention is to solve the difficult problem existed in above-mentioned prior art, the Fractured Zone detection method that a kind of attribute change rate is successively promoted is provided, on the basis comprising the information such as azimuth information, high-order rate of change of the poststack seismic data of inheriting tradition, introduce different scale, uprise time calculating of the various level in rank, and make this information stacked, excavate secondary fissure seam growth information not at the same level from big to small, the delineation of space effective range is provided, predicting the outcome of more accurate comprehensively fracture development effetive zone can be obtained.

The present invention is achieved by the following technical solutions:

The Fractured Zone detection method that a kind of attribute change rate is successively promoted, poststack central point seismic trace is launched to the road collection correlation-corrected of height battle array, the seismic trace of consecutive point carried out to the effective fracture development zone of different scale cracking initiation of the high-order orientation amplitude change rate calculating of mutative scale, accumulative certain growth strength, and be finally generalized to whole work area, obtain the maximum orientation high-order attribute change rate attribute of whole geological data, the multistage different scale fracture development integrated forecasting result be successively promoted to.

Said method comprising the steps of:

(1) poststack 3-d seismic data set is obtained;

(2) kth rank calculating required road collection square formation is selected to reset with coordinate: each sampling point of each seismic trace in the poststack 3-d seismic data set obtain step (1), get point centered by it, select adjacent in length and breadth to (k+1) * (k+1) neighboring track, build kJie road collection square formation, and this central point is set to initial point (0,0), then the coordinate of its adjacent spots is reset;

(3) the computing center road time difference relevant to neighboring track: computing center road with adjacent every together with cross-correlation, time shift amount when getting its maximum correlation carry out between TEC time error correction, eliminate the impact that stratum rises and falls;

(4) kth rank filter operator is utilized to calculate different azimuth attribute change rate;

(5) kth rank different azimuth attribute change rate maximal value is extracted;

(6) kth rank different azimuth attribute change rate mean value is extracted;

(7) each rank different azimuth attribute maximal value and mean value is extracted: repeat step (2)-(6) at each time sampling point, can k=2 be obtained, 4,6 ... N each rank different azimuth attribute maximal value and mean value; Wherein, N is any maximum integer that user selectes, and can select exponent number high arbitrarily, but after arriving certain rank, basic just constant, when specifically implementing, better by N=30 effect;

(8) extract be greater than the accumulative multistage different azimuth attribute mean value of preset value with (like this part of hydrocarbon content low yield or poor efficiency is not just added up, in case accumulative higher, but the part that hydrocarbon content is little is introduced into): to k=2 on each time sampling point, 4, 6, the each rank of N obtain the data acquisition of different azimuth attribute maximal value and mean value by step (7), get the maximal value in maximum value data set again, and each mean value being greater than certain given preset value in data acquisition of averaging carries out accumulative summation, namely be successively promoted, to obtain on this time sampling point of this road multistage different azimuth attribute mean value accumulative and,

(9) different scale FRACTURE PREDICTION achievement is exported: to each the time sampling point in every one in three-dimensional data, repetition step (2)-(8) complete multistage different azimuth attribute mean value and are successively promoted to extraction, then the multistage different azimuth attribute mean value obtained added up and export, be different scale fracture development integrated forecasting result, the result of calculation that the multistage different azimuth attribute change rate Fractured Zone be namely successively promoted to detects.

Described step (1) is achieved in that

After carrying out conventional processing by a geological data for distance prestack collection, obtain poststack imaging seismic data, form the poststack 3-d seismic data set in longitudinal direction, transverse direction, the degree of depth three directions.

The coordinate to its adjacent spots in described step (2) carries out replacement and is achieved in that

By with described initial point in length and breadth direction relative unit step-length measure the coordinate of adjacent spots reset, the initial value of k is 2; Coordinate is reset to the efficiency that can improve and calculate addressing.

Described step (4) is achieved in that

Each seismic trace of consecutive point after the correlation-corrected that the seismic trace of this central point and step (3) are obtained, utilize the kth rank filtering factor coefficient of first order derivative, calculate the kth rank attribute change rate in length and breadth, on oblique 4 different directions, obtain the kth high-order attribute change rate of the seismic trace of this central point on this time sampling point on four direction.

Described step (5) is achieved in that

Get the maximal value in the kth rank attribute change rate of the seismic trace on the four direction that obtains in step (4), be the maximum orientation characteristics rate of change in the kth rank of central point seismic trace on this time sampling point.

Described step (6) is achieved in that

Get the mean value in the kth rank attribute change rate of the seismic trace on the four direction that obtains in step (4), be the kth rank orientation characteristics rate of change of central point seismic trace on this time sampling point;

The inventive method has the following advantages:

(1) calculating of attribute change rate, considers in length and breadth, the calculating of oblique 4 different azimuth, and along with the change of yardstick dot matrix, can not neglect fracture orientation characteristic, length variations characteristic, and the space adapting to fracture development changes arbitrarily;

(2) calculating of attribute change rate, consider the multistage calculating (as from 2 rank to 20 rank) of first order derivative, obtain the automatic computing not relying on fixed size, simultaneously high-order again and again goes forward one by one, on mathematical principle, the detection of further outstanding small scale fracture development, improves the accuracy of detection of small scale Fractured Zone;

(3) calculating of attribute change rate, consider the comprehensive of different order, 1. get orientation maximal value the strongest outstanding one group of fracture development intensity of different order, the main body size of the yardstick reflection fracture development that maximal value is corresponding, and its orientation represents the main direction of this advantage crack group; 2. the successively promotion cumulative amount of the azimuth average of different order is got, the bulk strength of outstanding different scale fracture development, or the reflection of crack set, provide higher validity for the overall zone of delineation fracture development.

Compared with prior art, the invention has the beneficial effects as follows:

1) inherit direction characteristic, different azimuth fracture development state can be predicted;

2) can predict that the growth strength in maximum a kind of yardstick crack in multiple yardstick, scale and main body grow orientation;

3) introducing of many arrays and many operators the successively promotion to effective fracture, without the need to priori dimensional information, can carry out unitary set to the crack that difference grows scale and growth strength, the delineation the most effective development area of reservoir, evaluates " dessert " of the most favourable oil gas drilling.On the basis of traditional poststack seismic data, excavate more orientation, dimensional information as far as possible, automatically scale selection is completed and orientation is comprehensive, it is made to be embodied in multistage attribute change rate attribute, compare the attribute change rate body that traditional post-stack seismic data obtains, more accurate results of fracture prediction can be obtained.

Accompanying drawing explanation

Fig. 1, in 4 rank attribute change rate computation processes, 5x5 dot matrix selects and central point, consecutive point seismic trace position reset process, and by length and breadth, oblique 4 directions get 5 property calculation rate of change.For the calculating of each k rank attribute change rate provides similar guidance;

Fig. 2, the step block diagram of the inventive method;

Fig. 3-1, certain orientation, rank, oil field district 2 peak swing rate of change attributed graph;

Fig. 3-2, certain orientation, rank, oil field district 4 peak swing rate of change attributed graph;

Fig. 3-3, certain orientation, rank, oil field district 12 peak swing rate of change attributed graph;

Fig. 3-4, certain orientation, rank, oil field district 120 peak swing rate of change attributed graph;

Fig. 4 is that in the embodiment of the present invention, certain multistage orientation characteristics rate of change in rank, oil field district 2 ~ 20 is successively promoted to attributed graph.

Embodiment

Below in conjunction with accompanying drawing, the present invention is described in further detail:

The present invention is directed to effective fracture layer description problem in slit formation oilfield prospecting developing, devise and a kind ofly utilize successively being promoted of seismic properties rate of change to carry out the method for fracture detection development belt, can be applicable to the delineation of the effective expanded range of slit formation hydrocarbon-bearing pool in In Oil Field Exploration And Development.

The present invention is by adopting the multistage expression formula of high precision mutative scale, multi-faceted, the attribute change rate computing method be successively promoted to, can excavate and add up in superposition poststack seismic data to Orientation differences, dimensional variation, accumulative change relevant with relevant information that is fracture development, the validity information that the methods such as acquisition and existing traditional coherent algorithm, orientation amplitude change rate are higher, the extension making fracture development effetive zone is further clear and definite, provides better effect for the delineation of oil reservoir scope and reserves calculate.

(1) know-why

1. attribute change rate computing method

Usual attribute change rate is defined as:

$\mathrm{VAR}(x,y,t)=\sqrt{{\left(\frac{\mathrm{dA}(x,y,t)}{\mathrm{dx}}\right)}^2+{\left(\frac{\mathrm{dA}(x,y,t)}{\mathrm{dy}}\right)}^2}---\left(1\right)$

Wherein A (x, y, t) certain property value for a bit locating in 3-D seismics body, VAR (x, y, t) represents the Vector Mode along the variable quantity on x, y direction on this aspect, represents the intensity of this attribute change.Because seismological observation is in plane, longitudinal discrete sampling, can only be undertaken calculating and adjustment by different difference griddings and Orthogonal Composite in actual applications.By the Computing Principle of mathematics, differentiate difference gridding can be selected to comprise: difference coefficient before and after front difference coefficient, rear difference coefficient, center, the orthogonal and various ways such as difference coefficient before and after intersecting, also can consider the average of two or more difference coefficients in each direction, improve effective accuracy; Need in addition to carry out longitudinal adjustment to data, improve noise suppression and the impact of stabilizing longitudinal sampled point.Usual employing computing center's point both sides are attribute change rate average in certain hour window up and down, successively slides and reflect the result of calculation of this point.

2. association attributes rate of change computing method

When stratum fluctuations, in high-precision calculating, need the obliquity information introducing any point, what reflect to make calculating is the change of reservoir but not the change of stratal configuration.Be incorporated herein and calculate relevant time shift between road thus by completely along in length and breadth to the attribute change rate of plane computations, consider the space-variant vector of attribute along exhibition to three dimensions, the rate of change of its mould and attribute is:

$\mathrm{VAR}(x,y,t,{\mathrm{\Δt}}_x,{\mathrm{\Δt}}_y)=\sqrt{{\left(\frac{\mathrm{dA}(x,y,t,{\mathrm{\Δt}}_x)}{\mathrm{dx}}\right)}^2+{\left(\frac{\mathrm{dA}(x,y,t,{\mathrm{\Δt}}_y)}{\mathrm{dy}}\right)}^2}---\left(2\right)$

Wherein Δ tx, Δ ty represent that difference coefficient calculates the relevant time shift amount of wiggle trace of point-to-point transmission, get central point road and are calculated another road waveform correlation computations and obtain the time shift amount of maximal value.Road high to correlativity in computation process just carries out part time shift process, represents that base value reflection truly has the impact of fluctuations; General related coefficient then will not carry out the process in time shift, thus consider some mixed and disorderly mutually under the change of part attribute.Such improvement makes the calculating of attribute change rate be more suitable for stratum reflection the crack area of certain fluctuations.

3. 4 rank attribute change rate computing method

The mathematically calculating of a simple one-time continuous derivative, due to discrete relation, can adopt the difference scheme of different grid to calculate.In order to improve computational accuracy and reflect the impact of yardstick, as an illustration of the high-order approximation of first order derivative, to derive its computing formula with 4 rank attribute change rate computing method in high order.

(wherein h is neighboring track sampling interval for function u (x), u (x ± h) and u (x ± 2h), i.e. CDP interval) for any one direction being put centered by u (x) adjacent 2 attributes of both sides, launch the expression formula that can obtain first order derivative by Taylor.Have by Taylor expansion:

$u(x+h)=u\left(x\right)+h\frac{d}{\mathrm{dx}}u\left(x\right)+\frac{h^2}{2!}\frac{d^2}{{\mathrm{dx}}^2}u\left(x\right)+\frac{h^3}{3!}\frac{d^3}{{\mathrm{dx}}^3}u\left(x\right)+\frac{h^4}{4!}\frac{d^4}{{\mathrm{dx}}^4}u\left(x\right)+O\left(h^5\right)$

$u(x-h)=u\left(x\right)-h\frac{d}{\mathrm{dx}}u\left(x\right)+\frac{h^2}{2!}\frac{d^2}{{\mathrm{dx}}^2}u\left(x\right)-\frac{h^3}{3!}\frac{d^3}{{\mathrm{dx}}^3}u\left(x\right)+\frac{h^4}{4!}\frac{d^4}{{\mathrm{dx}}^4}u\left(x\right)+O\left(h^5\right)$

$u(x-2h)=u\left(x\right)-2h\frac{d}{\mathrm{dx}}u\left(x\right)+\frac{{\left(2h\right)}^2}{2!}\frac{d^2}{{\mathrm{dx}}^2}u\left(x\right)-\frac{{\left(2h\right)}^3}{3!}\frac{d^3}{{\mathrm{dx}}^3}u\left(x\right)+\frac{{\left(2h\right)}^4}{4!}\frac{d^4}{{\mathrm{dx}}^4}u\left(x\right)+O\left(h^5\right)$

Above-mentioned four expansions are associated with:

$\frac{1}{6}[u(x+2h)-u(x-2h)]-\frac{4}{3}[u(x+h)-u(x-h)]=2h\frac{d}{\mathrm{dx}}u\left(x\right)$

$\frac{d}{\mathrm{dx}}u\left(x\right)=\frac{\frac{1}{6}(u(x-2h)-\frac{4}{3}(u(x-h)+\frac{4}{3}(u(x+h)-\frac{1}{6}(u(x+2h)}{2h}---\left(4\right)$

4 rank expression formula filter operators of first order derivative are [1/6-4/304/3-1/6]=[0.167-1.3301.33-0.167].

4. multistage attribute change rate computing method

Copy 4 derivation modes of above-mentioned (3) formula, press from 2,4,6 respectively ..., successively recursion item by item, just can obtain k rank expression formula (k=2,4,6 of first order derivative, 8 ...) under the filter operator (see table 1) of different order.Due to the odd symmetry of difference scheme, each road, initial point both sides coefficient value is identical and symbol contrary, and as 4 rank [0.167-1.3301.33-0.167] operator, available [01.33-0.167] identifies R-branch coefficient, and P-branch coefficient gets its contrary sign.So just obtain all expression formulas of multistage high-precision attribute change rate, this provides calculating operator for developing mutative scale attribute change rate, and table 1 is first order derivative each rank myopia filter operator tables.

Table 1

(2) technology performing step

The present invention is directed to the crack elimination development belt entirety faced in actual In Oil Field Exploration And Development and iris wipe problem, the method adopted is carried out for actual poststack seismic data and is had orientation, multistage attribute change rate calculates, multistage cumulative review descriptive study, poststack central point seismic trace is launched to the road collection correlation-corrected of height battle array, the high-order orientation amplitude change rate seismic trace of consecutive point being carried out to mutative scale calculates, the effective fracture development zone of different scale cracking initiation of accumulative certain growth strength, and this step is generalized to whole work area the most at last, obtain the maximum orientation high-order attribute change rate attribute of whole geological data, the multistage different scale fracture development integrated forecasting result be successively promoted, and then provide means to the delineation of the effective fracture growth zone in this work area.

Technology realization flow as shown in Figure 2, choose collection when comprising the calculating of kth rank to carry out dip move-out, calculate different azimuth attribute change rate, get orientation maximal value and azimuth average, k is circulated, and obtain and multistagely accumulative be successively promoted to summing value, complete the calculating of a sampling point, after time and space completes whole data volume, form calculating achievement.Concrete steps are as follows:

1. earthquake three-dimensional poststack data volume

By the geological data that point gathers apart from prestack, the collection pre-service of conventional road, relative amplitude preserved processing and imaging are carried out to it, the poststack imaging seismic data volume of reflection subsurface formations reflectance signature can be obtained, form the 3-d seismic data set in longitudinal direction, transverse direction, the degree of depth three directions; Usually this data volume covers the floor area of-thousands of square kilometres up to a hundred, is made up of hundreds of thousands-millions of seismic trace;

2. different azimuth 4 rank attribute change rate calculates and realizes (4 rank illustration)

Sampling point involved by operator is calculated by k=4 rank attribute change rate, the calculating square formation (as shown in Figure 1) that to be a grid scope be (5 point × 5 point), (0 is set to by computing center's point attribute u coordinate in figure, 0), adjacent 24 attribute coordinates are followed successively by u (-2 ,-2), u (-2,-1), u (-2,0), u (-2,1), u (-2,2); U (-1 ,-2), u (-1 ,-1), u (-1,0), u (-1,1), u (-1,2); U (0 ,-2), u (0 ,-1), u (0,1), u (0,2); U (1 ,-2), u (1 ,-1), u (1,0), u (1,1), u (1,2); U (2 ,-2), u (2 ,-1), u (2,0), u (2,1), u (2,2); 4 rank expression formula filter operators of first order derivative are [01.33-0.167], also be [0.167-1.3301.33-0.167], respectively by 4 directions: west-east orientation (w-e), northwest-southeast calculate the attribute change rate of different azimuth to (n-s), east northeast-Nan Xi to (ne-sw) to (wn-es), North-south:

$u_{w-e}^4=0.167u(-\mathrm{2,0})-1.33u(-\mathrm{1,0})+0*u\left(\mathrm{0,0}\right)+1.33u\left(\mathrm{1,0}\right)-0.167u\left(\mathrm{2,0}\right);$

$u_{\mathrm{wn}-\mathrm{es}}^4=0.167u(-2,-2)-1.33u(-1,-1)+0*u\left(\mathrm{0,0}\right)+1.33u(1,1)-0.167u\left(\mathrm{2,2}\right);$

$u_{n-s}^4=0.167u(0,-2)-1.33u(0,-1)+0*u\left(\mathrm{0,0}\right)+1.33u(0,1)-0.167u\left(\mathrm{0,2}\right);---\left(5\right)$

$u_{\mathrm{ne}-\mathrm{sw}}^4=0.167u(2,-2)-1.33u(1,-1)+0*u\left(\mathrm{0,0}\right)+1.33u(-1,1)-0.167u(-2,2);$

Form different azimuth 4 rank attribute change rate result of calculation.

3. the multistage attribute change rate of different azimuth calculates and realizes

4 rank calculate and use 5x5 dot matrix, and 6 rank calculate and use 7x7 dot matrix ..., k rank use (k+1) x (k+1) dot matrix, Respectively by length and breadth, oblique 4 crossing directions, by 2 rank corresponding shown in table 1,4 rank, 6 rank ... k rank filter operator Continuous plus, will utilize the change order orientation characteristics rate of change of different range adjacent seismic trace computing center's point seismic trace under acquisition different scale; They represent respectively: different azimuth (4 orthogonal orientation), different order (k=2,4,6 ...) data acquisition

4. different scale and different azimuth multistage maximum attribute rate of change computing method

To the data set of above-mentioned formation, calculate in the achievement obtained at each often, choose 4 orientation, 2 ~ k rank (as 2 ~ 20 rank) maximal value is as by calculation level attribute change rate end value, calculated on whole 3-D data volume to step 3 by step 1, the result of calculation of the multistage maximum attribute rate of change in different scale orientation can be obtained, record obtains dot matrix corresponding to maximal value simultaneously, orientation, the most advantage fracture development intensity of this point can be obtained, grow yardstick (scale), grow the three-dimensional space data body of main direction, analyze the main body scale of different oil reservoir positions fracture development.

5. different scale and the multistage attribute change rate of different azimuth are successively promoted to computing method

In the computation process of kth rank attribute change rate, 4 orientation result of calculations are got to it is average, obtain the mean value of kth rank attribute change rate, represent the average growth strength in crack under a kind of fracture scale if to all different scale k there are all interested words of fracture development of certain growth strength (fracture development intensity is greater than certain preset value E), by being greater than the accumulative summation of preset value different scale fracture development intensity, the results of fracture prediction of full size just can be obtained.

${\stackrel{\‾}{u}}^k=0.25(u_{w-e}^k+u_{\mathrm{wn}-\mathrm{es}}^k+u_{n-s}^k+u_{\mathrm{ne}-\mathrm{sw}}^k)---\left(6\right)$

$U=\underset{k}{\mathrm{\Σ}}{\stackrel{\‾}{u}}^k$ (when time)

6. disposal of stability in attribute change rate

The various data types that computation attribute can select seismic processing to obtain, effective attribute comprises: seismic amplitude, frequency, phase place, the speed that seismic inversion obtains, elastic parameter and lithology, reservoir and other variable quantity, for suppressing noise and making extremely there is certain continuity before calculating, various process (three winks, trace integral differential+filtering, level and smooth) can be added to suppress the local jump of data and greatly indivedual value.

Fig. 3-1 to Fig. 3-4 is order orientation, system in Tahe Oilfield district 2,4,12,120 peak swing rate of change attributed graph, 4 width different order orientation amplitude change rate band unanimously indicates this district's solution fissure exactly and grows favourable distributive province, careful observation can find the raising along with order, and the small scale seam hole abnormal belt near large scale seam hole is further obvious.Tu4Wei Gai district poststack multistage (2 ~ 12 rank) orientation amplitude change rate is successively promoted to composite diagram, wherein every 4, rank orientation average amplitude change rate must be greater than preset value 1E8, composite diagram could be entered by accumulative superposition, successively be promoted by 2 ~ 12 rank and define continuous block shape extremely, indicate the abnormal block on obvious border, reflect corrosion Dong Feng district, and in little wire, dendritic exception is then and orientation relevant portion, be generally fracture developing zone, these wire as can be seen from Fig., dendritic abnormal population distribution is corrosion Dong Feng district neighbour, for effective development area in the secondary one-level crack with seam hole association.

The present invention is directed to effective fracture development belt delineation problem in slit formation oilfield prospecting developing, design and a kind ofly utilize successively being promoted of seismic properties rate of change to carry out the method for fracture detection development belt, adopt method for actual poststack seismic data carry out have pendage correct, growth direction, orientation, the multistage operator of different scale, successively a kind of calculating being successively promoted to attribute change rate of accumulative effective attribute rate of change, and this step is generalized to whole work area the most at last, the fracture development that the attribute change rate obtaining whole geological data is successively promoted predicts the outcome, can be applicable to the delineation of the effective expanded range of slit formation hydrocarbon-bearing pool in In Oil Field Exploration And Development.

Technique scheme is one embodiment of the present invention, for those skilled in the art, on the basis that the invention discloses application process and principle, be easy to make various types of improvement or distortion, and the method be not limited only to described by the above-mentioned embodiment of the present invention, therefore previously described mode is just preferred, and does not have restrictive meaning.

Claims (7)

1. the Fractured Zone detection method be successively promoted of an attribute change rate, it is characterized in that: described method is launched the road collection correlation-corrected of height battle array to poststack central point seismic trace, the seismic trace of consecutive point carried out to the effective fracture development zone of different scale cracking initiation of the high-order orientation amplitude change rate calculating of mutative scale, accumulative certain growth strength, and be finally generalized to whole work area, obtain the maximum orientation high-order attribute change rate attribute of whole geological data, the multistage different scale fracture development integrated forecasting result be successively promoted to.

2. the Fractured Zone detection method be successively promoted of attribute change rate according to claim 1, is characterized in that: said method comprising the steps of:

(1) poststack 3-d seismic data set is obtained;

(2) kth rank calculating required road collection square formation is selected to reset with coordinate: each sampling point of each seismic trace in the poststack 3-d seismic data set obtain step (1), get point centered by it, select adjacent in length and breadth to (k+1) * (k+1) neighboring track, build kJie road collection square formation, and this central point is set to initial point (0,0), then the coordinate of its adjacent spots is reset;

(3) the computing center road time difference relevant to neighboring track: computing center road with adjacent every together with cross-correlation, time shift amount when getting its maximum correlation carry out between TEC time error correction, eliminate the impact that stratum rises and falls;

(4) kth rank filter operator is utilized to calculate different azimuth attribute change rate;

(5) kth rank different azimuth attribute change rate maximal value is extracted;

(6) kth rank different azimuth attribute change rate mean value is extracted;

(7) each rank different azimuth attribute maximal value and mean value is extracted: repeat step (2)-(6) at each time sampling point, can k=2 be obtained, 4,6 ... N each rank different azimuth attribute maximal value and mean value; Wherein, N is any maximum integer that user selectes;

(8) extract be greater than the accumulative multistage different azimuth attribute mean value of preset value with: to k=2 on each time sampling point, 4,6, the each rank of N obtain the data acquisition of different azimuth attribute maximal value and mean value by step (7), get the maximal value in maximum value data set again, and each mean value being greater than certain given preset value in data acquisition of averaging carries out accumulative summation, namely be successively promoted, to obtain on this time sampling point of this road multistage different azimuth attribute mean value accumulative and;

(9) different scale FRACTURE PREDICTION achievement is exported: to each the time sampling point in every one in three-dimensional data, repetition step (2)-(8) complete multistage different azimuth attribute mean value and are successively promoted to extraction, then the multistage different azimuth attribute mean value obtained added up and export, being different scale fracture development integrated forecasting result.

3. the Fractured Zone detection method be successively promoted of attribute change rate according to claim 2, is characterized in that: described step (1) is achieved in that

After carrying out conventional processing by a geological data for distance prestack collection, obtain poststack imaging seismic data, form the poststack 3-d seismic data set in longitudinal direction, transverse direction, the degree of depth three directions.

4. the Fractured Zone detection method be successively promoted of attribute change rate according to claim 2, is characterized in that: the coordinate to its adjacent spots in described step (2) carries out replacement and is achieved in that

By with described initial point in length and breadth direction relative unit step-length measure the coordinate of adjacent spots reset, the initial value of k is 2; Coordinate is reset to the efficiency that can improve and calculate addressing.

5. the Fractured Zone detection method be successively promoted of attribute change rate according to claim 4, is characterized in that: described step (4) is achieved in that

Each seismic trace of consecutive point after the correlation-corrected that the seismic trace of this central point and step (3) are obtained, utilize the kth rank filtering factor coefficient of first order derivative, calculate the kth rank attribute change rate in length and breadth, on oblique 4 different directions, obtain the kth high-order attribute change rate of the seismic trace of this central point on this time sampling point on four direction.

6. the Fractured Zone detection method be successively promoted of attribute change rate according to claim 5, is characterized in that: described step (5) is achieved in that

Get the maximal value in the kth rank attribute change rate of the seismic trace on the four direction that obtains in step (4), be the maximum orientation characteristics rate of change in the kth rank of central point seismic trace on this time sampling point.

7. the Fractured Zone detection method be successively promoted of attribute change rate according to claim 6, is characterized in that: described step (6) is achieved in that

Get the mean value in the kth rank attribute change rate of the seismic trace on the four direction that obtains in step (4), be the kth rank orientation characteristics rate of change of central point seismic trace on this time sampling point.