CN104808244A - Ancient river course river bed erosion depth automatically recognizing method - Google Patents
Ancient river course river bed erosion depth automatically recognizing method Download PDFInfo
- Publication number
- CN104808244A CN104808244A CN201510173944.3A CN201510173944A CN104808244A CN 104808244 A CN104808244 A CN 104808244A CN 201510173944 A CN201510173944 A CN 201510173944A CN 104808244 A CN104808244 A CN 104808244A
- Authority
- CN
- China
- Prior art keywords
- seed points
- extrapolation
- seismic response
- response point
- main profile
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides an ancient river course river bed erosion depth automatically recognizing method. The method comprises setting river bed seed points; according to the coherence and the inclination of a seismic reflection stratum where the seed points are positioned, driving the seed points to automatically trace extrapolating seismic response points on a current main line, which are the same to the seed points in phase, in two directions with the crossline increasing and decreasing on the same main line; according the coherent mean value and the variance of the seismic reflection stratum where every extrapolating seismic response point is positioned, determining whether the tracing reaches river course boundaries, if so, stopping the tracing to obtain all the seed points on the current main line; according to the coherence and the inclination of all the seed points on the current main line, driving all the seed points on the current main line, automatically driving the seed points on all the main lines within the work area of the river course in the directions with the main lines increasing and decreasing. The ancient river course river bed erosion depth automatically recognizing method can help automatically recognize the erosion depth of ancient river beds, indirectly guide predication of underground river karst reservoir stratums and facilitates prediction of positions of efficient oil reservoirs.
Description
Technical field
The present invention relates to geophysical seismic reservoir prediction technique field, particularly a kind of paleo-streams River Bank Erosion degree of depth automatic identifying method.
Background technology
Underground river karst system is the Favorable Areas that efficient oil reservoir is grown.The growth of underground river karst system is controlled by earth's surface palaeodrainage pattern.The deep position of surface drainage system is the sluicing place of underground river karst system, for karstification provides hydrodynamic condition, is conducive to forming large-scale underground river karst Reservoir Body.
Although the carbonate karst reservoir growth of River Bank Erosion effect to both sides, riverbed and bottom position has building action, depth of erosion is excessively dark, easily destroys the Capped Ability of the cap rock of riverbed bottom position.Even if the karst reservoir of riverbed bottom part down is grown, be also unfavorable for Hydrocarbon Formation Reservoirs.Therefore, the identification of paleo-streams and depth of erosion thereof can instruct the underground river karst reservoir identifying and be conducive to Hydrocarbon Formation Reservoirs.
At present, the seismic properties such as relevant, amplitude, the frequency division of 3D seismic data are in conjunction with microtomy: dropping cut slice, horizon slice and strata slicing (Lei Ming, Wang Jiangong, Liu Caiyan etc.: the application of seismic slice technology in work area, Anda, geophysical prospecting for oil, 2010,45 (3): 418 ~ 422; Li Qingren, Chen Shoutian, Zhang Cai etc.: the paleo-streams sand body recognition methods of Hailaer Basin, Northeast China, geophysical prospecting for oil, 2006,41 (4): 402 ~ 4046; Han Xi, glad friend, Che Yanxin etc.: utilize seismic properties along layer analysis technique study sedimentary environment of fluvial facies, geophysical prospecting for oil, 2007,42 (1): 120 ~ 124) although river course plane or spatial can be portrayed, the depth of erosion in riverbed, paleo-streams can not be estimated.Large with country rock formation lithology difference bottom riverbed, carbonatite paleo-streams, wave impedance is the strong peak reflection of continuous print, and channel boundary place is mixed and disorderly weak reflection, and bottom whole riverbed, earthquake axis of reflection is continuous, can follow the trail of continuously along river course.But river course distribution is wide, tortuosity is high, branch's river channel is grown, and manual trace workload is very big.
Summary of the invention
Embodiments provide a kind of paleo-streams River Bank Erosion degree of depth automatic identifying method, in order to realize the automatic identification of the paleo-streams River Bank Erosion degree of depth, the method comprises:
A seismic response point is chosen as Seed Points at the maximum crest place of riverbed seismic section reflection wave;
Calculate the relevant of seismic reflection layer face, Seed Points place and inclination angle;
According to the relevant of seismic reflection layer face, Seed Points place and inclination angle, drive Seed Points to the both direction being positioned on same main profile, cross-track increases and reduces, on the current main profile of automatic tracing, identical with Seed Points phase place extrapolation seismic response point;
Calculate the relevant of seismic reflection layer face, each extrapolation seismic response point place and inclination angle, calculate each extrapolation seismic response point and 8 the seismic response points adjacent with this extrapolation seismic response point, the coherence average value of totally 9 seismic response points and variance;
According to coherence average value and variance, judge whether each extrapolation seismic response point has arrived channel boundary; When judging that extrapolation seismic response point does not arrive channel boundary, according to the relevant of seismic reflection layer face, each extrapolation seismic response point place and inclination angle, continuing with each extrapolation seismic response point is new Seed Points, to the both direction that cross-track increases and reduces, the extrapolation seismic response point that automatic tracing is identical with new Seed Points phase place; When judging that extrapolation seismic response point has arrived channel boundary place, stopping following the trail of, obtaining Seed Points all on current main profile;
According to the relevant of Seed Points all on current main profile and inclination angle, drive all Seed Points on current main profile, to the direction that main profile increases and reduces, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point, judge whether each extrapolation seismic response point has arrived channel boundary, when judging that extrapolation seismic response point has arrived channel boundary place, stop following the trail of, obtain all Seed Points on adjacent main survey line, until obtain the Seed Points within the scope of work area, river course on all main profiles;
According to the Seed Points on main profiles all within the scope of work area, river course, calculate the depth of erosion in riverbed, paleo-streams.
In one embodiment, according to the relevant of seismic reflection layer face, Seed Points place and inclination angle, drive Seed Points to being positioned on same main profile, both direction that cross-track increases and reduces, on the current main profile of automatic tracing, identical with Seed Points phase place extrapolation seismic response point, comprising:
Get in 6 seismic response points of cross-track augment direction, maximum coherent value and the inclination angle of correspondence thereof, to the direction being positioned on same main profile, cross-track increases, on the current main profile of automatic tracing, identical with Seed Points phase place extrapolation seismic response point, getting cross-track reduces in 6 seismic response points in direction, maximum coherent value and the inclination angle of correspondence thereof, to being positioned on same main profile, the cross-track direction of reducing, on the current main profile of automatic tracing, identical with Seed Points phase place extrapolation seismic response point.
In one embodiment, in being concerned with and inclination angle according to Seed Points all on current main profile, drive all Seed Points on current main profile, to the direction that main profile increases and reduces, automatic tracing is with on current main profile adjacent main survey line, in identical with Seed Points phase place extrapolation seismic response point, drive the direction that a Seed Points increases to main profile and reduces, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point comprises:
Get in 6 seismic response points of main profile augment direction, maximum coherent value and the inclination angle of correspondence thereof, to the direction being positioned on same cross-track, main profile increases, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point, getting main profile reduces in 6 seismic response points in direction, maximum coherent value and the inclination angle of correspondence thereof, to being positioned on same cross-track, the main profile direction of reducing, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point.
In one embodiment, according to the relevant of Seed Points all on current main profile and inclination angle, drive all Seed Points on current main profile, to the direction that main profile increases and reduces, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point, judge whether each extrapolation seismic response point has arrived channel boundary, when judging that extrapolation seismic response point has arrived channel boundary place, stop following the trail of, obtain all Seed Points on adjacent main survey line, until obtain the Seed Points within the scope of work area, river course on all main profiles, comprising:
According to the relevant of current main profile coboundary Seed Points and inclination angle, drive border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with border Seed Points phase place Seed Points;
According to the relevant of non-border Seed Points on current main profile and inclination angle, drive non-border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with non-border Seed Points phase place Seed Points.
In one embodiment, according to the relevant of current main profile coboundary Seed Points and inclination angle, drive border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with border Seed Points phase place Seed Points, comprising:
Judge that the extrapolation seismic response point of border Seed Points is positioned at the position on river course;
When judging that the extrapolation seismic response point of border Seed Points is positioned at inside, river course, to extrapolate, seismic response point is new Seed Points, to on same main profile, the cross-track identical extrapolation seismic response point of direction automatic tracing phase place that increases or reduce, until judge that extrapolation seismic response point has arrived channel boundary place, stop following the trail of, obtain the Seed Points on adjacent main survey line.
In one embodiment, according to the relevant of current main profile coboundary Seed Points and inclination angle, drive border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with border Seed Points phase place Seed Points, comprising:
When judging that the extrapolation seismic response point of border Seed Points is positioned at channel boundary or outside, determine that border Seed Points is invalid Seed Points, redefine germinative seed point, with germinative seed point for border Seed Points, to the extrapolation seismic response point that adjacent main survey line automatic tracing phase place is identical, obtain the Seed Points on adjacent main survey line; Germinative seed point is: adjacent with the invalid Seed Points finally determined and it is extrapolated, seismic response point is positioned at inside, river course.
In one embodiment, judge that extrapolation seismic response point has arrived channel boundary place, be specially:
When extrapolation seismic response point and 8 the seismic response points adjacent with extrapolation seismic response point, the coherence average value of totally 9 seismic response points is little, when variance yields is large, has proved channel boundary.
Embodiment of the present invention technical scheme is based on seismic event similarity principle, on the basis that riverbed Seed Points is set, Seed Points extrapolation is driven with stratigraphic dip and seismic waveshape similarity, adopt border, multiwindow rotation sweep identification riverbed simultaneously, riverbed boundary constraint extrapolation scope, finally realizes the depth of erosion prediction in the river course of phase more than many.Therefore, the paleo-streams River Bank Erosion degree of depth automatic identifying method provided by the embodiment of the present invention, can instruct the prediction of underground river karst reservoir indirectly.Meanwhile, the identification of the fossil river bed degree of depth is conducive to the position of predicting efficient oil reservoir.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a application's part, does not form limitation of the invention.In the accompanying drawings:
Fig. 1 is the schematic flow sheet of paleo-streams River Bank Erosion degree of depth automatic identifying method in the embodiment of the present invention;
Fig. 2 is the schematic diagram of the multiwindow rotation sweep identification channel boundary in the embodiment of the present invention centered by Seed Points;
Fig. 3 is the idiographic flow schematic diagram of paleo-streams River Bank Erosion degree of depth automatic identifying method in the embodiment of the present invention.
Fig. 4 is that the maximum crest place at riverbed seismic section reflection wave in the embodiment of the present invention arranges Seed Points schematic diagram;
Fig. 5 is based on riverbed Seed Points automatic tracing trajectory plane schematic diagram in the embodiment of the present invention;
Fig. 6 is the schematic diagram that the extrapolation seismic response point of border Seed Points in the embodiment of the present invention is positioned at inside, river course;
Fig. 7 is the schematic diagram that the extrapolation seismic response point of border Seed Points in the embodiment of the present invention is positioned at outside and border, river course.
Embodiment
For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with embodiment and accompanying drawing, the present invention is described in further details.At this, exemplary embodiment of the present invention and illustrating for explaining the present invention, but not as a limitation of the invention.
The embodiment of the present invention is based on seismic event similarity principle, on the basis that riverbed Seed Points is set, Seed Points extrapolation is driven with stratigraphic dip and seismic waveshape similarity, adopt border, multiwindow rotation sweep identification riverbed simultaneously, when occurring that average coherence value is little, during the large exception of variance yields, prove channel boundary, riverbed boundary constraint extrapolation scope, finally realizes the depth of erosion prediction in the river course of phase more than many.Be described in detail below.
Fig. 1 is the schematic flow sheet of paleo-streams River Bank Erosion degree of depth automatic identifying method in the embodiment of the present invention, and as shown in Figure 1, the method comprises the steps:
Step 101: choose a seismic response point at the maximum crest place of riverbed seismic section reflection wave as Seed Points;
Step 102: calculate the relevant of seismic reflection layer face, Seed Points place and inclination angle;
Step 103: according to the relevant of seismic reflection layer face, Seed Points place and inclination angle, drive Seed Points to being positioned on same main profile, both direction that cross-track increases and reduces, on the current main profile of automatic tracing, identical with Seed Points phase place extrapolation seismic response point;
Step 104: calculate the relevant of seismic reflection layer face, each extrapolation seismic response point place and inclination angle, calculate each extrapolation seismic response point and 8 the seismic response points adjacent with this extrapolation seismic response point, the coherence average value of totally 9 seismic response points and variance;
Step 105: according to coherence average value and variance, judges whether each extrapolation seismic response point has arrived channel boundary; When judging that extrapolation seismic response point does not arrive channel boundary, according to the relevant of seismic reflection layer face, each extrapolation seismic response point place and inclination angle, continuing with each extrapolation seismic response point is the both direction that new Seed Points increases to cross-track and reduces, the extrapolation seismic response point that automatic tracing is identical with new Seed Points phase place; When judging that extrapolation seismic response point has arrived channel boundary place, stopping following the trail of, obtaining Seed Points all on current main profile;
Step 106: according to the relevant of Seed Points all on current main profile and inclination angle, drive all Seed Points on current main profile, to the direction that main profile increases and reduces, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point, judge whether each extrapolation seismic response point has arrived channel boundary, when judging that extrapolation seismic response point has arrived channel boundary place, stop following the trail of, obtain all Seed Points on adjacent main survey line, until obtain the Seed Points within the scope of work area, river course on all main profiles;
Step 107: according to the Seed Points on main profiles all within the scope of work area, river course, calculates the depth of erosion in riverbed, paleo-streams.Particularly, after riverbed Seed Points has all been estimated, just can know the riverbed degree of depth.Meanwhile, we manual trace can go out the degree of depth of bunding, and the degree of depth deducting bunding by the degree of depth in riverbed is exactly the depth of erosion in riverbed.
In one embodiment, according to the relevant of seismic reflection layer face, Seed Points place and inclination angle, drive Seed Points to being positioned on same main profile, both direction that cross-track increases and reduces, on the current main profile of automatic tracing, identical with Seed Points phase place extrapolation seismic response point, comprising:
Get in 6 seismic response points of cross-track augment direction, maximum coherent value and the inclination angle of correspondence thereof, to the direction being positioned on same main profile, cross-track increases, on the current main profile of automatic tracing, identical with Seed Points phase place extrapolation seismic response point, getting cross-track reduces in 6 seismic response points in direction, maximum coherent value and the inclination angle of correspondence thereof, to being positioned on same main profile, the cross-track direction of reducing, on the current main profile of automatic tracing, identical with Seed Points phase place extrapolation seismic response point.
In one embodiment, in being concerned with and inclination angle according to Seed Points all on current main profile, drive all Seed Points on current main profile, to the direction that main profile increases and reduces, automatic tracing is with on current main profile adjacent main survey line, in identical with Seed Points phase place extrapolation seismic response point, drive the direction that a Seed Points increases to main profile and reduces, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point comprises:
Get in 6 seismic response points of main profile augment direction, maximum coherent value and the inclination angle of correspondence thereof, to the direction being positioned on same cross-track, main profile increases, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point, getting main profile reduces in 6 seismic response points in direction, maximum coherent value and the inclination angle of correspondence thereof, to being positioned on same cross-track, the main profile direction of reducing, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point.
In one embodiment, according to the relevant of Seed Points all on current main profile and inclination angle, drive all Seed Points on current main profile, to the direction that main profile increases and reduces, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point, judge whether each extrapolation seismic response point has arrived channel boundary, when judging that extrapolation seismic response point has arrived channel boundary place, stop following the trail of, obtain all Seed Points on adjacent main survey line, until obtain the Seed Points within the scope of work area, river course on all main profiles, comprising:
According to the relevant of current main profile coboundary Seed Points and inclination angle, drive border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with border Seed Points phase place Seed Points;
According to the relevant of non-border Seed Points on current main profile and inclination angle, drive non-border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with non-border Seed Points phase place Seed Points.
In one embodiment, according to the relevant of current main profile coboundary Seed Points and inclination angle, drive border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with border Seed Points phase place Seed Points, comprising:
Judge that the extrapolation seismic response point of border Seed Points is positioned at the position on river course;
When judging that the extrapolation seismic response point of border Seed Points is positioned at inside, river course, to extrapolate, seismic response point is new Seed Points, to on same main profile, the cross-track identical extrapolation seismic response point of direction automatic tracing phase place that increases or reduce, until judge that extrapolation seismic response point has arrived channel boundary place, stop following the trail of, obtain the Seed Points on adjacent main survey line.
In one embodiment, according to the relevant of current main profile coboundary Seed Points and inclination angle, drive border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with border Seed Points phase place Seed Points, comprising:
When judging that the extrapolation seismic response point of border Seed Points is positioned at channel boundary or outside, determine that border Seed Points is invalid Seed Points, redefine germinative seed point, with germinative seed point for border Seed Points, to the extrapolation seismic response point that adjacent main survey line automatic tracing phase place is identical, obtain the Seed Points on adjacent main survey line; Germinative seed point is: adjacent with the invalid Seed Points finally determined and it is extrapolated, seismic response point is positioned at inside, river course.
In one embodiment, judge that extrapolation seismic response point has arrived channel boundary place, be specially:
When extrapolation seismic response point and 8 the seismic response points adjacent with extrapolation seismic response point, the coherence average value of totally 9 seismic response points is little, when variance yields is large, has proved channel boundary.
Be described with example more below, so that understand how to implement the present invention.
To there are obvious earthquake radiation lineups in paleo-streams drop-down on seismic section, bottom riverbed with under to cover stratum reflecting interface strong.Paleo-streams provided by the invention River Bank Erosion degree of depth automatic identifying method, to Seed Points and each extrapolation seismic response point, 9 square-shaped planar windows are all adopted to calculate the relevant of the reflection aspect at Seed Points and extrapolation seismic response point place and inclination angle (as shown in Figure 2), each square-shaped planar window represents the geological state of a Seed Points or extrapolation seismic response point, a round dot in each square-shaped planar window corresponding diagram 5, Fig. 6 and Fig. 7.According to relevant and inclination angle, drive Seed Points and each extrapolation seismic response point on the survey line within the scope of work area, river course, the extrapolation seismic response point that automatic tracing phase place is identical, in the process of automatic tracing, calculate coherence average value and variance yields that the coherence average value at each extrapolation seismic response point place and variance yields are calculating 9 square-shaped planar windows, channel boundary is judged by coherence average value and variance yields, when occurring that coherence average value is little, when variance yields is large, proof has been extrapolated to border, riverbed, until obtain the Seed Points on all surveys line, realize whole district's riverbed seismic response lineups automatically to identify, finally realize the depth of erosion prediction in the river course of phase more than many.(implication of hereinafter described " extrapolation " is the identical extrapolation seismic response point of automatic tracing phase place) specific as follows:
Fig. 3 is the idiographic flow schematic diagram of paleo-streams River Bank Erosion degree of depth automatic identifying method in the embodiment of the present invention, and Fig. 4 is that the maximum crest place at riverbed seismic section reflection wave in the embodiment of the present invention arranges Seed Points schematic diagram; Fig. 5 is based on riverbed Seed Points automatic tracing trajectory plane schematic diagram in the embodiment of the present invention; As shown in Fig. 3, Fig. 4 and Fig. 5, choose a seismic response point at the maximum crest place of riverbed seismic section reflection wave as Seed Points (as shown in Figure 4), with Seed Points u (i, j, t
i,j) centered by, wherein, u is the seismic response point on boundary at the bottom of riverbed, and i is inline (main profile), and j is crossline (cross-track), t
i,jfor the time, adopt 9 relevant and inclination angles (as described in epimere) comprising the reflection aspect at the square-shaped planar window estimation Seed Points place of 9 sampled points.Relevant and apparent dip ck (Pk, Qk), the wherein k=1 of estimation 9 windows, 2 ..., 9.Drive Seed Points extrapolation by inclination angle, realize the automatic tracing of seismic event bottom riverbed.
As shown in Figure 2, Seed Points to the method for each line direction extrapolation (automatic tracing) is specifically: estimation window 1, 2, 4, 5, 7, the coherent value of 8, get wherein maximum coherent value and the apparent dip of correspondence thereof, to be positioned at same inline line crossline reduce driving riverbed, direction at the bottom of boundary's automatic formation trace, estimation window 2, 3, 5, 6, 8, 9 coherent value, get wherein maximum coherent value and the apparent dip of correspondence thereof, boundary's automatic formation trace at the bottom of the driving riverbed, direction that same inline line crossline increases, estimation window 1, 2, 3, 4, 5, 6 coherent value, get wherein maximum coherent value and the apparent dip of correspondence thereof, boundary's automatic formation trace at the bottom of the driving riverbed, direction that same crossline road inline increases, estimation window 4, 5, 6, 7, 8, 9 coherent value, get wherein maximum coherent value and the apparent dip of correspondence thereof, boundary's automatic formation trace at the bottom of the driving riverbed, direction that same crossline road inline reduces.When each extrapolation seismic response point is specifically to each line direction extrapolation, is also described above, does not repeat them here.
Simultaneously in conjunction with border, multiwindow rotation sweep identification riverbed, the variance and the mean value that calculate the coherent value of 9 windows are respectively
away from river course lateral boundaries, seismic amplitude, waveform similarity, coherence average value is large, and variance yields is little.At the seismic waveshape near lateral boundaries place, river course in mixed and disorderly reflection, its similarity degree is low, and coherence average value is little, and variance yields is large.Therefore, the variance of 9 window coherent value and the change identification channel boundary of mean value can be passed through.
Particularly, as shown in Fig. 3, Fig. 4 and Fig. 5, first, a sampled point is chosen at the maximum crest place of seismic reflection lineups middle part, riverbed as Seed Points (shown in Fig. 4 orbicular spot): u (i, j, t
i,j), wherein, u is the sampled point on boundary at the bottom of riverbed, and i is main profile inline, j is cross-track crossline, t
i,jfor riverbed is positioned at the time at drawing lines (i, j) place.In Fig. 5, black circle is that Seed Points is in plane projection.
Then, as shown in Figures 3 and 5, the both direction extrapolation (shown in the arrow of the transverse direction in Fig. 5 on Li line) main profile Li line increasing to cross-track crossline and reduces, whether be extrapolated to channel boundary place simultaneously if carrying out rotary window scanning recognition at each extrapolation seismic response point.In riverbed middle position, the average coherence mean value of 9 windows is large, and variance yields is little, and Seed Points and contiguous extrapolation seismic response point are all positioned on same riverbed, and average coherence value is large, and variance yields is little.In the process of both sides extrapolation, there is the exception that coherence average value is little and variance yields is large, show to be extrapolated to both sides, riverbed boundary u (i, j-m t
i, j-m), u (i, j+n, t
i, j+n), now, stop the extrapolation of both sides.
By u (i, j, the t of estimation
i,j) (j=j-m, j-m+1 ...., j+n-1, j+n) as new Seed Points, to same cross-track crossline road, the direction extrapolation (the longitudinal arrow being positioned at many cross-tracks such as cross-track Lj in Fig. 5 is shown) that main profile inline increases and reduces.With u (i, j, t
i,j) (j=j-m, j-m+1 ...., j+n-1, j+n) extrapolate in the process of the riverbed sampled point on adjacent threads, first from both sides, riverbed boundary u (i, j-m, t Li line
i, j-m), u (i, j+n, t
i, j+n) riverbed response point (seismic response point of namely extrapolating) u (i ± 1, j-m, t on the contiguous inline line that starts to extrapolate
i ± 1, j-m), u (i ± 1, j+n, t
i ± 1, j+n), rotation sweep is carried out to each extrapolation seismic response point and carries out riverbed Boundary Recognition.
As shown in Figure 6, if with frontier point u (i, j-m, t
i, j-m) be extrapolation seismic response point (being positioned at the top of black circle on cross-track Lj-m, longitudinal arrow indication open circles in Fig. 6) u (i+1, j-m, the t of Seed Points (black circle in Fig. 6)
i+1, j-m) be in inside, riverbed, then again with u (i+1, j-m, t
i+1, j-m) be Seed Points, to the direction extrapolation that same Li+1 line crossline reduces (in Fig. 6 on Li+1 line lateral arrows direction), be extrapolated to channel boundary always and stop.And non-frontier point u (i, j-m+k, the t on Li line
i, j-m+k), (k=1,2,3 ...) and (in such as Fig. 6 on main profile Li, be positioned at cross-track Lj-m+1, the soft dot on Lj-m+2, Lj-m+3), then according on Li line extrapolation obtain u (i, j-m+k, t
i, j-m+k) (k=1,2,3...) as Seed Points, the extrapolation seismic response point (round dot in Fig. 6 pointed by Lj-m+1, Lj-m+2 and Lj-m+3 line upward arrow) that obtains on Li+1 line of extrapolating.
As shown in Figure 7, if with frontier point u (i, j-m, t
i, j-m) be riverbed response point (namely extrapolate seismic response point) u (i+1, j-m, the t of Seed Points (in Fig. 7 black circle) extrapolation (in Fig. 7 Lj-m line upward arrow direction) adjacent line
i+1, j-m) be positioned at border, riverbed and outside, then look u (i, j-m, t
i, j-m) be invalid Seed Points, and with u (i, j-m+k, t
i, j-m+k), k=1,2,3 .., p-1, p, p+1... are extrapolation seismic response point u (i+1, j-m+k, the t that Seed Points extrapolation (in Fig. 7 Lj-m+1, Lj-m+2 and Lj-m+3 line upward arrow direction) is corresponding
i+1, j-m+k), k=1,2,3 .., p-1, p, p+1....Until as extrapolation seismic response u (i+1, j-m+p, t
i+1, j-m+p) be positioned at inside, river course, then by u (i, j-m+p, t
i, j-m+p) u (i+1, j-m+p, the t of its correspondence of extrapolating as effective Seed Points
i+1, j-m+p).Frontier point u (i, j-m, t is carried out by above-mentioned same method
i, j-m), u (i, j+n, t
i, j+n) successively extrapolate to adjacent line, finally realize whole district's riverbed seismic response lineups and automatically identify.
Technical solution of the present invention can identify fossil river bed depth of erosion automatically, by the identification of fossil river bed, indirectly instructs the prediction of underground river karst reservoir.Meanwhile, the identification of the fossil river bed degree of depth is conducive to the position of predicting efficient oil reservoir.The present invention proposes the paleo-streams River Bank Erosion degree of depth automatic identifying method based on seismic event similarity principle.The method, on the basis arranging riverbed Seed Points, drives Seed Points extrapolation with stratigraphic dip and seismic waveshape similarity, adopts border, multiwindow rotation sweep identification riverbed simultaneously, riverbed boundary constraint extrapolation scope.Finally realize the depth of erosion prediction in the river course of phase more than many.Identify that portraying of underground river karst reservoir can be instructed indirectly in the locus in deep riverbed.
In addition, in the application of Harrar Ha Tang, by paleo-streams provided by the invention River Bank Erosion degree of depth automatic identifying method and device, have identified 8 paleo-streamss and the riverbed degree of depth thereof.Based on this, identify multiple underground river karst system, underground river karst oil reservoir is efficient oil reservoir, and actual probing also confirms the validity of Forecasting Methodology.
Obviously, those skilled in the art should be understood that, each module of the above-mentioned embodiment of the present invention or each step can realize with general calculation element, they can concentrate on single calculation element, or be distributed on network that multiple calculation element forms, alternatively, they can realize with the executable program code of calculation element, thus, they can be stored and be performed by calculation element in the storage device, and in some cases, step shown or described by can performing with the order be different from herein, or they are made into each integrated circuit modules respectively, or the multiple module in them or step are made into single integrated circuit module to realize.Like this, the embodiment of the present invention is not restricted to any specific hardware and software combination.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (7)
1. a paleo-streams River Bank Erosion degree of depth automatic identifying method, is characterized in that, comprising:
A seismic response point is chosen as Seed Points at the maximum crest place of riverbed seismic section reflection wave;
Calculate the relevant of seismic reflection layer face, Seed Points place and inclination angle;
According to the relevant of seismic reflection layer face, Seed Points place and inclination angle, drive Seed Points to the both direction being positioned on same main profile, cross-track increases and reduces, on the current main profile of automatic tracing, identical with described Seed Points phase place extrapolation seismic response point;
Calculate the relevant of seismic reflection layer face, each extrapolation seismic response point place and inclination angle, calculate each extrapolation seismic response point and 8 the seismic response points adjacent with this extrapolation seismic response point, the coherence average value of totally 9 seismic response points and variance;
According to described coherence average value and variance, judge whether each extrapolation seismic response point has arrived channel boundary; When judging that extrapolation seismic response point does not arrive channel boundary, according to the relevant of seismic reflection layer face, each extrapolation seismic response point place and inclination angle, continuing with each extrapolation seismic response point is new Seed Points, to the both direction that cross-track increases and reduces, the extrapolation seismic response point that automatic tracing is identical with described new Seed Points phase place; When judging that extrapolation seismic response point has arrived channel boundary place, stopping following the trail of, obtaining Seed Points all on current main profile;
According to the relevant of Seed Points all on current main profile and inclination angle, drive all Seed Points on current main profile, to the direction that main profile increases and reduces, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point, judge whether each extrapolation seismic response point has arrived channel boundary, when judging that extrapolation seismic response point has arrived channel boundary place, stop following the trail of, obtain all Seed Points on adjacent main survey line, until obtain the Seed Points within the scope of work area, river course on all main profiles;
According to the Seed Points on main profiles all within the scope of work area, river course, calculate the depth of erosion in riverbed, paleo-streams.
2. the method for claim 1, it is characterized in that, described according to the relevant of seismic reflection layer face, Seed Points place and inclination angle, drive Seed Points to the both direction being positioned on same main profile, cross-track increases and reduces, on the current main profile of automatic tracing, identical with described Seed Points phase place extrapolation seismic response point, comprising:
Get in 6 seismic response points of cross-track augment direction, maximum coherent value and the inclination angle of correspondence thereof, to the direction being positioned on same main profile, cross-track increases, on the current main profile of automatic tracing, identical with described Seed Points phase place extrapolation seismic response point, getting cross-track reduces in 6 seismic response points in direction, maximum coherent value and the inclination angle of correspondence thereof, to being positioned on same main profile, the cross-track direction of reducing, on the current main profile of automatic tracing, identical with described Seed Points phase place extrapolation seismic response point.
3. the method for claim 1, it is characterized in that, described according to the relevant of Seed Points all on current main profile and inclination angle, drive all Seed Points on current main profile, to the direction that main profile increases and reduces, automatic tracing is with on current main profile adjacent main survey line, in identical with Seed Points phase place extrapolation seismic response point, drive the direction that a Seed Points increases to main profile and reduces, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point comprises:
Get in 6 seismic response points of main profile augment direction, maximum coherent value and the inclination angle of correspondence thereof, to the direction being positioned on same cross-track, main profile increases, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point, getting main profile reduces in 6 seismic response points in direction, maximum coherent value and the inclination angle of correspondence thereof, to being positioned on same cross-track, the main profile direction of reducing, automatic tracing with on current main profile adjacent main survey line, identical with Seed Points phase place extrapolation seismic response point.
4. the method for claim 1, it is characterized in that, described according to the relevant of Seed Points all on current main profile and inclination angle, drive all Seed Points on current main profile, to the direction that main profile increases and reduces, on automatic tracing and current main profile adjacent main survey line, the extrapolation seismic response point identical with Seed Points phase place, judge whether each extrapolation seismic response point has arrived channel boundary, when judging that extrapolation seismic response point has arrived channel boundary place, stop following the trail of, obtain all Seed Points on adjacent main survey line, until obtain the Seed Points within the scope of work area, river course on all main profiles, comprise:
According to the relevant of current main profile coboundary Seed Points and inclination angle, drive border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with described border Seed Points phase place Seed Points;
According to the relevant of non-border Seed Points on current main profile and inclination angle, drive non-border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with described non-border Seed Points phase place Seed Points.
5. method as claimed in claim 4, it is characterized in that, described according to the relevant of current main profile coboundary Seed Points and inclination angle, drive border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with described border Seed Points phase place Seed Points, comprising:
Judge that the extrapolation seismic response point of border Seed Points is positioned at the position on river course;
When judging that the extrapolation seismic response point of border Seed Points is positioned at inside, river course, be new Seed Points with described extrapolation seismic response point, to on same main profile, the cross-track identical extrapolation seismic response point of direction automatic tracing phase place that increases or reduce, until judge that extrapolation seismic response point has arrived channel boundary place, stop following the trail of, obtain the Seed Points on adjacent main survey line.
6. method as claimed in claim 5, it is characterized in that, described according to the relevant of current main profile coboundary Seed Points and inclination angle, drive border Seed Points to adjacent main survey line automatic tracing, obtain on adjacent main survey line, identical with described border Seed Points phase place Seed Points, also comprise:
When judging that the extrapolation seismic response point of border Seed Points is positioned at channel boundary or outside, determine that described border Seed Points is invalid Seed Points, redefine germinative seed point, with germinative seed point for border Seed Points, to the extrapolation seismic response point that adjacent main survey line automatic tracing phase place is identical, obtain the Seed Points on adjacent main survey line; Described germinative seed point is: adjacent with the invalid Seed Points finally determined and it is extrapolated, seismic response point is positioned at inside, river course.
7. the method as described in claim as arbitrary in claim 1 to 6, is characterized in that, described judgement extrapolation seismic response point has arrived channel boundary place, is specially:
When extrapolation seismic response point and 8 the seismic response points adjacent with extrapolation seismic response point, the coherence average value of totally 9 seismic response points is little, when variance yields is large, has proved channel boundary.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510173944.3A CN104808244B (en) | 2015-04-14 | 2015-04-14 | Ancient stream channel River Bank Erosion depth automatic identifying method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510173944.3A CN104808244B (en) | 2015-04-14 | 2015-04-14 | Ancient stream channel River Bank Erosion depth automatic identifying method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104808244A true CN104808244A (en) | 2015-07-29 |
| CN104808244B CN104808244B (en) | 2017-03-08 |
Family
ID=53693229
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510173944.3A Active CN104808244B (en) | 2015-04-14 | 2015-04-14 | Ancient stream channel River Bank Erosion depth automatic identifying method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104808244B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110133719A (en) * | 2019-05-23 | 2019-08-16 | 中国石油大学(华东) | Glutenite cascade correlation method based on 90 ° of phase shift data bodies, seed point tracking |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5056066A (en) * | 1990-06-25 | 1991-10-08 | Landmark Graphics Corporation | Method for attribute tracking in seismic data |
| US5251184A (en) * | 1991-07-09 | 1993-10-05 | Landmark Graphics Corporation | Method and apparatus for finding horizons in 3D seismic data |
| CN103592681A (en) * | 2013-09-16 | 2014-02-19 | 电子科技大学 | Signal classification based seismic image horizon tracking method |
| CN104181596A (en) * | 2014-08-27 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | Geologic horizon automatic tracking method and device |
-
2015
- 2015-04-14 CN CN201510173944.3A patent/CN104808244B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5056066A (en) * | 1990-06-25 | 1991-10-08 | Landmark Graphics Corporation | Method for attribute tracking in seismic data |
| US5251184A (en) * | 1991-07-09 | 1993-10-05 | Landmark Graphics Corporation | Method and apparatus for finding horizons in 3D seismic data |
| CN103592681A (en) * | 2013-09-16 | 2014-02-19 | 电子科技大学 | Signal classification based seismic image horizon tracking method |
| CN104181596A (en) * | 2014-08-27 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | Geologic horizon automatic tracking method and device |
Non-Patent Citations (2)
| Title |
|---|
| 朱桂娟等: "基于最佳时窗刻画技术的河道相储集砂体识别—以港西油田为例", 《勘探地球物理进展》 * |
| 龙隆等: "浅层河道砂体的地震准确识别与精细刻画", 《大庆石油地质与开发》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110133719A (en) * | 2019-05-23 | 2019-08-16 | 中国石油大学(华东) | Glutenite cascade correlation method based on 90 ° of phase shift data bodies, seed point tracking |
| CN110133719B (en) * | 2019-05-23 | 2020-08-04 | 中国石油大学(华东) | Glutenite period sub-division method based on 90-degree phase shift data volume and seed point tracking |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104808244B (en) | 2017-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Liu et al. | Sealing mechanisms in volcanic faulted reservoirs in Xujiaweizi extension, northern Songliao Basin, northeastern China | |
| Zweigel et al. | Reservoir geology of the Utsira Formation at the first industrial-scale underground CO2 storage site (Sleipner area, North Sea) | |
| Paull et al. | High-resolution bathymetry of the axial channels within Monterey and Soquel submarine canyons, offshore central California | |
| CN104698497B (en) | A kind of fine calibration of thin sand-mud interbed layer position | |
| Aarre et al. | Seismic detection of subtle faults and fractures | |
| CN110424955B (en) | Internal diving method for complex broken block | |
| CN105319585B (en) | A kind of method hidden using thin-layers interference amplitude recovery identification of hydrocarbon | |
| Young et al. | Exploration of an alluvial aquifer in Oman by time-domain electromagnetic sounding | |
| Liu et al. | Characterization of architectural elements of Ordovician fractured-cavernous carbonate reservoirs, Tahe Oilfield, China | |
| Eiken | Twenty years of monitoring CO2 injection | |
| Krois et al. | Miano field, Pakistan, A case history of model driven exploration | |
| CN109138837A (en) | Reservoir geology guidance method | |
| CN104049283A (en) | Method for identification and prediction of fore-land basin extremely-thick conglomerate body | |
| Rizzo et al. | Deep electrical resistivity tomography for the hydrogeological setting of Muro Lucano Mounts Aquifer (Basilicata, Southern Italy) | |
| Jenkins et al. | Fit for purpose monitoring-a progress report on the CO2CRC Otway Stage 3 project | |
| Jewell et al. | Phases of blowout initiation and stabilization on Padre Island revealed through ground-penetrating radar and remotely sensed imagery | |
| Ishido et al. | Changes in geophysical observables caused by CO2 injection into saline aquifers | |
| Hemmings-Sykes | The influence of faulting on hydrocarbon migration in the Kupe area, south Taranaki Basin, New Zealand | |
| Mabey | Geophysical studies in the intermontane basins in southern California | |
| DesRoches et al. | Influence of fracture anisotropy and lithological heterogeneity on wellfield response in a fluvial sandstone aquifer of the Carboniferous Moncton Subbasin, Canada | |
| CN104808244A (en) | Ancient river course river bed erosion depth automatically recognizing method | |
| Kolenković et al. | Regional capacity estimates for CO2 geological storage in deep saline aquifers–Upper Miocene sandstones in the SW part of the Pannonian basin | |
| CN110910267A (en) | Method and device for determining marine carbonate rock interlayer | |
| Li et al. | Characteristics of discontinuities in Heifangtai landslide area in Gansu, China | |
| Trofimov | Deep CMP seismic surveying along the Tatseis-2003 geotraverse across the Volga-Ural petroliferous province |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |