CN104808244B - Ancient stream channel River Bank Erosion depth automatic identifying method - Google Patents
Ancient stream channel River Bank Erosion depth automatic identifying method Download PDFInfo
- Publication number
- CN104808244B CN104808244B CN201510173944.3A CN201510173944A CN104808244B CN 104808244 B CN104808244 B CN 104808244B CN 201510173944 A CN201510173944 A CN 201510173944A CN 104808244 B CN104808244 B CN 104808244B
- Authority
- CN
- China
- Prior art keywords
- point
- main profile
- seed point
- extrapolation
- seismic response
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The invention provides a kind of ancient stream channel River Bank Erosion depth automatic identifying method, the method includes：Riverbed seed point is set；According to the relevant and inclination angle in seed point place seismic reflection layer face, seed point is driven to the both direction being located on same main profile, cross-track increases and reduces, on the current main profile of automatic tracing and seed point phase place identical extrapolation seismic response point；According to coherence average value and the variance in each extrapolation seismic response point place seismic reflection layer face, judge whether, to channel boundary, when having determined at channel boundary, to stop following the trail of, obtain all of seed point on current main profile；According to the relevant and inclination angle of all seed points on current main profile, all seed points on current main profile are driven, the direction for increasing to main profile and reducing, the seed point in the range of automatic tracing river course work area on all main profiles.Program energy automatic identification fossil river bed depth of erosion, instructs the prediction of underground river karst reservoir indirectly, beneficial to the efficient oil reservoir position of prediction.
Description
Technical field
The present invention relates to geophysical seismic reservoir prediction technique field, more particularly to a kind of ancient stream channel River Bank Erosion depth
Automatic identifying method.
Background technology
Underground river karst system is the Favorable Areas of efficient oil reservoir development.The development of underground river karst system is controlled by earth's surface Gu water
System.The deep position of surface drainage is the sluicing place of underground river karst system, is that karstification provides hydrodynamic condition, is conducive to shape
Become large-scale underground river karst Reservoir Body.
Although River Bank Erosion acts on the carbonate karst reservoir development to riverbed both sides and bottom position building action,
But depth of erosion is too deep, the Capped Ability of the cap rock of riverbed bottom position is easily destroyed.Even if the karst of riverbed bottom part down
Reservoir development, is also unfavorable for Hydrocarbon Formation Reservoirs.Therefore, the identification of ancient stream channel and its depth of erosion can instruct identification to be conducive to oil gas
Become the underground river karst reservoir of Tibetan.
At present, the relevant, seismic properties such as amplitude, frequency dividing of 3D seismic data combine microtomy：Dropping cut slice, along layer
Section and strata slicing (Lei Ming, Wang Jiangong, Liu Caiyan etc.：Seismic slice technology is in application, the petroleum geophysics in Anda work area
Exploration, 2010,45 (3)：418～422；Li Qingren, Chen Shoutian, Zhang Cai etc.：The ancient stream channel sand body recognition methods of Hailaer Basin, Northeast China,
Geophysical prospecting for oil, 2006,41 (4)：402～4046；Han Xi, glad friend, car prolong letter etc.：Divided along layer using seismic properties
Analysis technique study sedimentary environment of fluvial facies, geophysical prospecting for oil, 2007,42 (1)：120～124) although river can be portrayed
Road plane or spatial, but the depth of erosion in ancient stream channel riverbed can not be estimated.Carbonate rock ancient stream channel riverbed bottom with enclose
Rock formation lithology difference is big, and wave impedance, in continuous strong wave peak reflection, is in weak reflection in a jumble at channel boundary, whole riverbed bottom
Seismic reflection axle is continuous, can continuously follow the trail of along river course.But river course distribution is wide, tortuosity is high, branch's river channel development, manually
Follow the trail of workload very big.
Content of the invention
A kind of ancient stream channel River Bank Erosion depth automatic identifying method is embodiments provided, in order to realize ancient stream channel river
The automatic identification of bed depth of erosion, the method include：
A seismic response point is chosen at the maximum crest of riverbed seismic profile back wave as seed point；
Calculate the relevant and inclination angle in seed point place seismic reflection layer face；
According to the relevant and inclination angle in seed point place seismic reflection layer face, drive seed point to being located on same main profile,
Cross-track increase and reduce both direction, on the current main profile of automatic tracing and seed point phase place identical extrapolation earthquake
Response point；
The relevant and inclination angle in each extrapolation seismic response point place seismic reflection layer face is calculated, is rung with calculating each extrapolation
Should put and the 8 seismic response points adjacent with the extrapolation seismic response point, the coherence average value of totally 9 seismic response points and side
Difference；
According to coherence average value and variance, judge whether each extrapolation seismic response point has arrived channel boundary；Outer when judging
When pushing away seismic response point and not arriving channel boundary, according to the relevant of each extrapolation seismic response point place seismic reflection layer face and incline
Angle, it is new seed point to continue with each extrapolation seismic response point, and the both direction for increasing to cross-track and reducing is chased after automatically
Track and new seed point phase place identical extrapolation seismic response point；When judging that extrapolation seismic response point has been arrived at channel boundary,
Stop following the trail of, obtain all of seed point on current main profile；
According to the relevant and inclination angle of all seed points on current main profile, all seed points on current main profile are driven,
The direction for increasing to main profile and reducing, on automatic tracing main profile adjacent with current main profile and seed point phase place identical
Extrapolation seismic response point, judges whether each extrapolation seismic response point has arrived channel boundary, when judgement extrapolation seismic response point is arrived
When at channel boundary, stop following the trail of, all seed points on adjacent main profile are obtained, until obtaining institute in the range of river course work area
The seed point having on main profile；
According to the seed point on all main profiles in the range of river course work area, the depth of erosion in ancient stream channel riverbed is calculated.
In one embodiment, according to the relevant and inclination angle in seed point place seismic reflection layer face, seed point is driven to position
On same main profile, the cross-track both direction that increases and reduce, on the current main profile of automatic tracing and seed point phase place
Identical extrapolation seismic response point, including：
Take in 6 seismic response points of cross-track augment direction, maximum coherent value and its corresponding inclination angle, to being located at
On same main profile, cross-track increase direction, on the current main profile of automatic tracing and seed point phase place identical extrapolation ground
Shake response point, takes in 6 seismic response points that cross-track reduces direction, maximum coherent value and its corresponding inclination angle, Xiang Wei
On same main profile, the direction of cross-track minimizing, with the extrapolation of seed point phase place identical on the current main profile of automatic tracing
Seismic response point.
In one embodiment, at the relevant and inclination angle according to all seed points on current main profile, currently main survey is driven
All seed points on line, to main profile increase and reduce direction, on automatic tracing main profile adjacent with current main profile and
In seed point phase place identical extrapolation seismic response point, the direction that one seed point of driving increases to main profile and reduces, automatically
Following the trail of on main profile adjacent with current main profile is included with seed point phase place identical extrapolation seismic response point：
Take in 6 seismic response points of main profile augment direction, maximum coherent value and its corresponding inclination angle, to be located at
On one cross-track, the direction that increases of main profile, on automatic tracing main profile adjacent with current main profile and seed point phase place phase
With extrapolation seismic response point, take in 6 seismic response points that main profile reduces direction, maximum coherent value and its corresponding incline
Angle, to be located at same cross-track on, main profile reduce direction, on automatic tracing main profile adjacent with current main profile and
Seed point phase place identical extrapolation seismic response point.
In one embodiment, according to the relevant and inclination angle of all seed points on current main profile, current main profile is driven
On all seed points, to main profile increase and reduce direction, on automatic tracing main profile adjacent with current main profile and plant
Son point phase place identical extrapolation seismic response point, judges whether each extrapolation seismic response point has arrived channel boundary, outer when judging
When pushing away seismic response point and having arrived at channel boundary, stop following the trail of, all seed points on adjacent main profile are obtained, until obtaining river
Seed point in the range of road work area on all main profiles, including：
According to the relevant and inclination angle of current main profile coboundary seed point, drive border seed point automatic to adjacent main profile
Follow the trail of, obtain on adjacent main profile and border seed point phase place identical seed point；
According to the relevant and inclination angle of non-border seed point on current main profile, non-border seed point is driven to adjacent main profile
Automatic tracing, obtains on adjacent main profile and non-border seed point phase place identical seed point.
In one embodiment, according to the relevant and inclination angle of current main profile coboundary seed point, border seed point is driven
To adjacent main profile automatic tracing, on adjacent main profile and border seed point phase place identical seed point is obtained, including：
Judge that the extrapolation seismic response point of border seed point is located at the position on river course；
When the extrapolation seismic response point for judging border seed point is located at river course inside, to extrapolate, seismic response point is new
Seed point, on same main profile, the cross-track direction automatic tracing phase place identical extrapolation seismic response that increases or reduce
Point, until judgement extrapolation seismic response point has been arrived at channel boundary, stops following the trail of, obtains the seed point on adjacent main profile.
In one embodiment, according to the relevant and inclination angle of current main profile coboundary seed point, border seed point is driven
To adjacent main profile automatic tracing, on adjacent main profile and border seed point phase place identical seed point is obtained, including：
When the extrapolation seismic response point for judging border seed point is located at channel boundary or outside, determine that border seed point is
Invalid seed point, redefines germinative seed point, with germinative seed point as border seed point, to adjacent main profile automatic tracing phase
Position identical extrapolation seismic response point, obtains the seed point on adjacent main profile；Germinative seed point is：Invalid with last determination
Seed point is adjacent and its extrapolation seismic response point is located inside river course.
In one embodiment, judge that extrapolation seismic response point has been arrived at channel boundary, specially：
When extrapolation seismic response point and with the extrapolation adjacent 8 seismic response points of seismic response point, totally 9 seismic response points
Coherence average value little, when variance yields is big, it was demonstrated that arrived channel boundary.
Embodiment of the present invention technical scheme is based on seismic wave similarity principle, on the basis of riverbed seed point is set, with
Stratigraphic dip and seismic waveform similitude drive seed point extrapolation, while recognizing riverbed border, river using multiwindow rotation sweep
Bedside bound constrained extrapolation scope, finally realizes the depth of erosion prediction in a plurality of river course of many phases.Therefore, carried by the embodiment of the present invention
For ancient stream channel River Bank Erosion depth automatic identifying method, the prediction of underground river karst reservoir can be instructed indirectly.Meanwhile, fossil river bed
The identification of depth is conducive to predicting the position of efficient oil reservoir.
Description of the drawings
Accompanying drawing described herein is used for providing a further understanding of the present invention, constitutes the part of the application, not
Constitute limitation of the invention.In the accompanying drawings：
Fig. 1 is the schematic flow sheet of River Bank Erosion depth automatic identifying method in ancient stream channel 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 on seed point；
Fig. 3 is the idiographic flow schematic diagram of River Bank Erosion depth automatic identifying method in ancient stream channel in the embodiment of the present invention.
Fig. 4 is to arrange seed point schematic diagram at the maximum crest of riverbed seismic profile back wave in the embodiment of the present invention；
Fig. 5 is to be based on riverbed seed point automatic tracing trajectory plane schematic diagram in the embodiment of the present invention；
Fig. 6 is that the extrapolation seismic response point of border seed point in the embodiment of the present invention is located at the schematic diagram inside river course；
Fig. 7 is the signal that the extrapolation seismic response point of border seed point in the embodiment of the present invention is located at the outside and border in river course
Figure.
Specific embodiment
For making the object, technical solutions and advantages of the present invention become more apparent, with reference to embodiment and accompanying drawing, right
The present invention is described in further details.Here, the exemplary embodiment of the present invention and its illustrate for explaining the present invention, but simultaneously
Not as a limitation of the invention.
The embodiment of the present invention is based on seismic wave similarity principle, on the basis of riverbed seed point is arranged, with stratigraphic dip
Seed point extrapolation is driven with seismic waveform similitude, while riverbed border is recognized using multiwindow rotation sweep, average when occurring
Coherent value is little, during the big exception of variance yields, it was demonstrated that arrived channel boundary, and riverbed boundary constraint is extrapolated scope, finally realizes a plurality of
The depth of erosion prediction in river course of many phases.It is specifically described below.
Fig. 1 is the schematic flow sheet of River Bank Erosion depth automatic identifying method in ancient stream channel in the embodiment of the present invention, as Fig. 1 institute
Show, the method comprises the steps：
Step 101：A seismic response point is chosen at the maximum crest of riverbed seismic profile back wave as seed point；
Step 102：Calculate the relevant and inclination angle in seed point place seismic reflection layer face；
Step 103：According to the relevant and inclination angle in seed point place seismic reflection layer face, seed point is driven to positioned at same master
On survey line, cross-track increase and reduce both direction, on the current main profile of automatic tracing and outside seed point phase place identical
Push away seismic response point；
Step 104：The relevant and inclination angle in each extrapolation seismic response point place seismic reflection layer face is calculated, is calculated outside each
Seismic response point and the 8 seismic response points adjacent with the extrapolation seismic response point are pushed away, the coherence average of totally 9 seismic response points
Value and variance；
Step 105：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 phase in each extrapolation seismic response point place seismic reflection layer face
Dry and inclination angle, it is both direction that new seed point increases to cross-track and reduces to continue with each extrapolation seismic response point,
Automatic tracing and new seed point phase place identical extrapolation seismic response point；When judgement extrapolation seismic response point has arrived channel boundary
During place, stop following the trail of, obtain all of seed point on current main profile；
Step 106：According to the relevant and inclination angle of all seed points on current main profile, drive all on current main profile
Seed point, the direction for increasing to main profile and reducing, on automatic tracing main profile adjacent with current main profile and seed point phase place
Identical extrapolation seismic response point, judges whether each extrapolation seismic response point has arrived channel boundary, when judgement rings with extrapolating
When should put at channel boundary, stop following the trail of, all seed points on adjacent main profile are obtained, until obtaining river course work area model
Enclose the seed point on interior all main profiles；
Step 107：According to the seed point on all main profiles in the range of river course work area, the erosion depth in ancient stream channel riverbed is calculated
Degree.Specifically, after riverbed seed point has all been estimated, it becomes possible to know riverbed depth.Meanwhile, we can manual trace go out
The depth of bunding, the depth for deducting bunding by the depth in riverbed are exactly the depth of erosion in riverbed.
In one embodiment, according to the relevant and inclination angle in seed point place seismic reflection layer face, seed point is driven to position
On same main profile, the cross-track both direction that increases and reduce, on the current main profile of automatic tracing and seed point phase place
Identical extrapolation seismic response point, including：
Take in 6 seismic response points of cross-track augment direction, maximum coherent value and its corresponding inclination angle, to being located at
On same main profile, cross-track increase direction, on the current main profile of automatic tracing and seed point phase place identical extrapolation ground
Shake response point, takes in 6 seismic response points that cross-track reduces direction, maximum coherent value and its corresponding inclination angle, Xiang Wei
On same main profile, the direction of cross-track minimizing, with the extrapolation of seed point phase place identical on the current main profile of automatic tracing
Seismic response point.
In one embodiment, at the relevant and inclination angle according to all seed points on current main profile, currently main survey is driven
All seed points on line, to main profile increase and reduce direction, on automatic tracing main profile adjacent with current main profile and
In seed point phase place identical extrapolation seismic response point, the direction that one seed point of driving increases to main profile and reduces, automatically
Following the trail of on main profile adjacent with current main profile is included with seed point phase place identical extrapolation seismic response point：
Take in 6 seismic response points of main profile augment direction, maximum coherent value and its corresponding inclination angle, to be located at
On one cross-track, the direction that increases of main profile, on automatic tracing main profile adjacent with current main profile and seed point phase place phase
With extrapolation seismic response point, take in 6 seismic response points that main profile reduces direction, maximum coherent value and its corresponding incline
Angle, to be located at same cross-track on, main profile reduce direction, on automatic tracing main profile adjacent with current main profile and
Seed point phase place identical extrapolation seismic response point.
In one embodiment, according to the relevant and inclination angle of all seed points on current main profile, current main profile is driven
On all seed points, to main profile increase and reduce direction, on automatic tracing main profile adjacent with current main profile and plant
Son point phase place identical extrapolation seismic response point, judges whether each extrapolation seismic response point has arrived channel boundary, outer when judging
When pushing away seismic response point and having arrived at channel boundary, stop following the trail of, all seed points on adjacent main profile are obtained, until obtaining river
Seed point in the range of road work area on all main profiles, including：
According to the relevant and inclination angle of current main profile coboundary seed point, drive border seed point automatic to adjacent main profile
Follow the trail of, obtain on adjacent main profile and border seed point phase place identical seed point；
According to the relevant and inclination angle of non-border seed point on current main profile, non-border seed point is driven to adjacent main profile
Automatic tracing, obtains on adjacent main profile and non-border seed point phase place identical seed point.
In one embodiment, according to the relevant and inclination angle of current main profile coboundary seed point, border seed point is driven
To adjacent main profile automatic tracing, on adjacent main profile and border seed point phase place identical seed point is obtained, including：
Judge that the extrapolation seismic response point of border seed point is located at the position on river course；
When the extrapolation seismic response point for judging border seed point is located at river course inside, to extrapolate, seismic response point is new
Seed point, on same main profile, the cross-track direction automatic tracing phase place identical extrapolation seismic response that increases or reduce
Point, until judgement extrapolation seismic response point has been arrived at channel boundary, stops following the trail of, obtains the seed point on adjacent main profile.
In one embodiment, according to the relevant and inclination angle of current main profile coboundary seed point, border seed point is driven
To adjacent main profile automatic tracing, on adjacent main profile and border seed point phase place identical seed point is obtained, including：
When the extrapolation seismic response point for judging border seed point is located at channel boundary or outside, determine that border seed point is
Invalid seed point, redefines germinative seed point, with germinative seed point as border seed point, to adjacent main profile automatic tracing phase
Position identical extrapolation seismic response point, obtains the seed point on adjacent main profile；Germinative seed point is：Invalid with last determination
Seed point is adjacent and its extrapolation seismic response point is located inside river course.
In one embodiment, judge that extrapolation seismic response point has been arrived at channel boundary, specially：
When extrapolation seismic response point and with the extrapolation adjacent 8 seismic response points of seismic response point, totally 9 seismic response points
Coherence average value little, when variance yields is big, it was demonstrated that arrived channel boundary.
Illustrated with example again below, in order to understand how to implement the present invention.
Ancient stream channel exist on seismic profile obvious earthquake radiation lineups drop-down, riverbed bottom with cover stratum reflection
Interface is strong.The ancient stream channel River Bank Erosion depth automatic identifying method that the present invention is provided, to seed point and each extrapolation seismic response
Point, all calculated using 9 square-shaped planar windows the relevant of the reflection aspect that seed point and extrapolation seismic response point are located and
Inclination angle (as shown in Figure 2), each square-shaped planar window represent the geological state of a seed point or extrapolation seismic response point, per
Individual square-shaped planar window corresponds to a round dot in Fig. 5, Fig. 6 and Fig. 7.According to relevant and inclination angle, drive seed point and each
, on the survey line in the range of river course work area, automatic tracing phase place identical extrapolation seismic response point, certainly for extrapolation seismic response point
During dynamic tracking, calculate the coherence average value at each extrapolation seismic response point and variance yields is and calculates 9 squares
The coherence average value of plane window and variance yields, judge channel boundary by coherence average value and variance yields, relevant flat when occurring
Average is little, when variance yields is big, it was demonstrated that be extrapolated to riverbed border, until obtaining the seed point on all surveys line, realizes whole district river
Bed seismic response lineups automatic identification, finally realizes the depth of erosion prediction in a plurality of river course of many phases.(hereinafter institute specific as follows
The implication of " extrapolation " for stating is automatic tracing phase place identical extrapolation seismic response point)：
Fig. 3 is the idiographic flow schematic diagram of River Bank Erosion depth automatic identifying method in ancient stream channel in the embodiment of the present invention, Fig. 4
It is, in the embodiment of the present invention, seed point schematic diagram is set at the maximum crest of riverbed seismic profile back wave；Fig. 5 is the present invention
Riverbed seed point automatic tracing trajectory plane schematic diagram is based in embodiment；As shown in Fig. 3, Fig. 4 and Fig. 5, the earthquake in riverbed is cutd open
A seismic response point is chosen at the maximum crest of face back wave as seed point (as shown in Figure 4), with seed point u (i, j,
ti,j) centered on, wherein, u is the seismic response point on riverbed bottom circle, and it is that crossline (survey by contact that i is inline (main profile), j
Line), ti,jThe reflection aspect being located for time, the square-shaped planar windows estimation seed points for including 9 sampled points using 9
Relevant and inclination angle (as described in epimere).Relevant and apparent dip ck (Pk, Qk), the wherein k=1,2 ..., 9 of 9 windows of estimation.Logical
Cross inclination angle and seed point extrapolation is driven, realize the automatic tracing of riverbed bottom seismic event.
As shown in Fig. 2 method of the seed point specifically to each line direction extrapolation (automatic tracing) is：Estimation window 1,2,4,
5th, 7,8 coherent value, takes maximum of which coherent value and its corresponding apparent dip, to positioned at same inline line crossline
The direction driving riverbed bottom circle automatic formation trace of minimizing, 2,3,5,6,8,9 coherent value of estimation window, take maximum of which and be concerned with
Value and its corresponding apparent dip, drive riverbed bottom circle automatic formation trace to the direction that same inline line crossline increases,
1,2,3,4,5,6 coherent value of estimation window, takes maximum of which coherent value and its corresponding apparent dip, to same crossline
The direction driving riverbed bottom circle automatic formation trace that road inline increases, 4,5,6,7,8,9 coherent value of estimation window, take wherein most
Big coherent value and its corresponding apparent dip, drive riverbed bottom interlayer position to the direction that same crossline road inline is reduced
Automatic tracing.Each extrapolation seismic response point specifically to each line direction extrapolate when, as also mentioned above, will not be described here.
Riverbed border is recognized in combination with multiwindow rotation sweep, calculate variance and the mean value of the coherent value of 9 windows
RespectivelyAway from river course lateral boundaries, seismic amplitude, waveform similarity, it is concerned with flat
Average is big, and variance yields is little.Reflect in mixed and disorderly in the seismic waveform at river course lateral boundaries, its similarity degree is low, coherence average
Value is little, and variance yields is big.The change that therefore, it can the variance by 9 window coherent value and mean value recognizes channel boundary.
Specifically, as shown in Fig. 3, Fig. 4 and Fig. 5, first, in the maximum crest of riverbed seismic reflection lineups middle part
Place chooses a sampled point as seed point (shown in Fig. 4 orbicular spot)：u(i,j,ti,j), wherein, u is the sampling on riverbed bottom circle
Point, it is cross-track crossline, t that i is main profile inline, ji,jThe time at drawing lines (i, j) place is located at for riverbed.In Fig. 5
Black circle is seed point in plane projection.
Then, as shown in Figures 3 and 5, two sides for increasing to cross-track crossline on main profile Li line and reducing
(shown in the horizontal arrow in Fig. 5 on Li line) is pushed out, while rotary window scanning is carried out in each extrapolation seismic response point
Identify whether to be extrapolated at channel boundary.In riverbed middle position, the average coherence mean value of 9 windows is big, variance yields
Little, seed point and neighbouring extrapolation seismic response point are all located on same riverbed, and average coherence value is big, and variance yields is little.In both sides
Occur the big exception of the little and variance yields of coherence average value during extrapolation, show to be extrapolated to riverbed both sides boundary u
(i,j-m ti,j-m)、u(i,j+n,ti,j+n), now, stop the extrapolation of both sides.
U (i, j, t by estimationi,j) (j=j-m, j-m+1 ...., j+n-1, j+n) as new seed point, to same
Cross-track crossline road, main profile inline increases and the direction extrapolation of minimizing (is located at cross-track Lj etc. many in Fig. 5
Shown in longitudinal arrow of bar cross-track).With u (i, j, ti,j) (j=j-m, j-m+1 ...., j+n-1, j+n) extrapolation adjacent threads
On riverbed sampled point during, first from riverbed both sides boundary u (i, j-m, t on Li linei,j-m)、u(i,j+n,ti,j+n)
Start to extrapolate adjacent to riverbed response point on inline line (seismic response point of extrapolating) u (i ± 1, j-m, ti±1,j-m)、u(i±1,j
+n,ti±1,j+n), carrying out rotation sweep to each extrapolation seismic response point carries out riverbed Boundary Recognition.
If as shown in fig. 6, with boundary point u (i, j-m, ti,j-m) for seed point (black circle in Fig. 6) extrapolation ground
Shake response point (positioned at the top of black circle, longitudinal arrow indication open circles on cross-track Lj-m in Fig. 6) u (i+1, j-m,
ti+1,j-m) inside riverbed, then again with u (i+1, j-m, ti+1,j-m) it is seed point, reduce to same Li+1 line crossline
Direction extrapolation (lateral arrows direction on Li+1 line in Fig. 6), be extrapolated to always channel boundary termination.And the non-border on Li line
Point u (i, j-m+k, ti,,j-m+k), (k=1,2,3 ...) (in such as Fig. 6 on main profile Li, be located at cross-track Lj-m+1,
Soft dot on Lj-m+2, Lj-m+3), then u (i, j-m+k, t is obtained according to extrapolating on Li linei,,j-m+k) (k=1,
2,3...) as seed point, extrapolation obtains extrapolation seismic response point on Li+1 line (Lj-m+1, Lj-m+2 and Lj-m+3 in Fig. 6
Round dot pointed by line upward arrow).
If as shown in fig. 7, with boundary point u (i, j-m, ti,j-m) extrapolate (in Fig. 7 for seed point (black circle in Fig. 7)
Lj-m line upward arrow direction) adjacent line riverbed response point (extrapolate seismic response point) u (i+1, j-m, ti+1,j-m) it is located at riverbed
Border and outside, then regard u (i, j-m, ti,j-m) it is invalid seed point, and with u (i, j-m+k, ti,j-m+k), k=1,2,3 .., p-
1, p, p+1... are seed point extrapolation (Lj-m+1, Lj-m+2 and Lj-m+3 line upward arrow direction in Fig. 7) corresponding extrapolation earthquake
Response point u (i+1, j-m+k, ti+1,j-m+k), k=1,2,3 .., p-1, p, p+1....Until when extrapolation seismic response u (i+1,
j-m+p,ti+1,j-m+p) be located inside river course, then by u (i, j-m+p, ti,j-m+p) extrapolate its corresponding u as effective seed point
(i+1,j-m+p,ti+1,j-m+p).Boundary point u (i, j-m, t are carried out in same way as described abovei,j-m)、u(i,j+n,ti,j+n) to neighbour
Line is gradually extrapolated, and finally realizes whole district's riverbed seismic response lineups automatic identification.
Technical solution of the present invention can automatic identification fossil river bed depth of erosion, by the identification of fossil river bed, indirectly instruct dark
The prediction of river karst reservoir.Meanwhile, the identification of Furukawa bed depth is conducive to predicting the position of efficient oil reservoir.The present invention proposes base
Ancient stream channel River Bank Erosion depth automatic identifying method in seismic wave similarity principle.The method is arranging the base of riverbed seed point
On plinth, seed point extrapolation is driven with stratigraphic dip and seismic waveform similitude, while riverbed is recognized using multiwindow rotation sweep
Border, riverbed boundary constraint extrapolation scope.The depth of erosion prediction in a plurality of many phase river course is finally realized.Recognize the sky in deep riverbed
Between position can instruct indirectly portraying for underground river karst reservoir.
In addition, in the application of Harrar Ha Tang, the ancient stream channel River Bank Erosion depth automatic identification side that provided by the present invention
Method and device, have identified 8 ancient stream channels and its riverbed depth.Based on this, multiple underground river karst systems, underground river are identified
Karst oil reservoir is efficient oil reservoir, and actual probing also confirms that 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 permissible
Realized with general computing device, they can be concentrated on single computing device, or be distributed in multiple computing devices
On the network for being constituted, alternatively, they can be realized with the executable program code of computing device, it is thus possible to by it
Store and executed by computing device in the storage device, and in some cases, can be held with the order being different from herein
The shown or described step of row, or they are fabricated to each integrated circuit modules respectively, or will be multiple in them
Module or step are fabricated to single integrated circuit module to realize.So, the embodiment of the present invention is not restricted to any specific hard
Part and software are combined.
The preferred embodiments of the present invention are the foregoing is only, the present invention is not limited to, for the skill of this area
For art personnel, the embodiment of the present invention can have various modifications and variations.All within the spirit and principles in the present invention, made
Any modification, equivalent substitution and improvement etc., should be included within the scope of the present invention.
Claims (7)
1. a kind of ancient stream channel River Bank Erosion depth automatic identifying method, it is characterised in that include：
A seismic response point is chosen at the maximum crest of riverbed seismic profile back wave as seed point u (i, j, ti,j), its
In, i is main profile, and j is cross-track, ti,jFor the time；
Calculate the relevant and inclination angle in seed point place seismic reflection layer face；
According to the relevant and inclination angle in seed point place seismic reflection layer face, seed point is driven on same main profile, contact
Survey line j increase and reduce both direction, on the current main profile of automatic tracing and the seed point phase place identical extrapolation earthquake
Response point；
The relevant and inclination angle in each extrapolation seismic response point place seismic reflection layer face is calculated, calculates each extrapolation seismic response point
And the 8 seismic response points adjacent with the extrapolation seismic response point, the coherence average value of totally 9 seismic response points and variance；
According to the coherence average value and variance, judge whether each extrapolation seismic response point has arrived channel boundary；Outer when judging
When pushing away seismic response point and not arriving channel boundary, according to the relevant of each extrapolation seismic response point place seismic reflection layer face and incline
Angle, it is new seed point to continue with each extrapolation seismic response point, the both direction for increasing to cross-track j and reducing, automatically
Follow the trail of and the new seed point phase place identical extrapolation seismic response point；When judgement extrapolation seismic response point has arrived channel boundary
During place, stop following the trail of, obtain all of seed point on current main profile；
According to the relevant and inclination angle of all seed points on current main profile, all seed points on current main profile are driven, to master
The direction that survey line i increases and reduces, is extrapolated with seed point phase place identical on automatic tracing main profile adjacent with current main profile
Seismic response point, judges whether each extrapolation seismic response point has arrived channel boundary, when judgement extrapolation seismic response point has arrived river
During road boundary, stop following the trail of, all seed points on adjacent main profile are obtained, until obtaining all masters in the range of river course work area
Seed point on survey line；
According to the seed point on all main profiles in the range of river course work area, the depth of erosion in ancient stream channel riverbed is calculated.
2. the method for claim 1, it is characterised in that described according to the relevant of seed point place seismic reflection layer face and
Inclination angle, drives seed point to the both direction being located on same main profile, cross-track j increases and reduces, and automatic tracing is current
On main profile and the seed point phase place identical extrapolation seismic response point, including：
Take in 6 seismic response points of cross-track j augment direction, maximum coherent value and its corresponding inclination angle, to be located at
On one main profile, cross-track j increase direction, on the current main profile of automatic tracing and outside the seed point phase place identical
Seismic response point is pushed away, is taken in 6 seismic response points that cross-track j reduces direction, maximum coherent value and its corresponding inclination angle,
To the direction being located on same main profile, cross-track j is reduced, on the current main profile of automatic tracing and the seed point phase place
Identical extrapolation seismic response point.
3. the method for claim 1, it is characterised in that all seed points is relevant on the current main profile of the basis
And inclination angle, drive all seed points on current main profile, the direction for increasing to main profile i and reducing, automatic tracing with current
On the adjacent main profile of main profile and in seed point phase place identical extrapolation seismic response point, a seed point is driven to main profile i
Increase and reduce direction, on automatic tracing main profile adjacent with current main profile and seed point phase place identical extrapolation earthquake
Response point includes：
Take in 6 seismic response points of main profile i augment direction, maximum coherent value and its corresponding inclination angle, same to being located at
On cross-track, the direction that increases of main profile i, on automatic tracing main profile adjacent with current main profile and seed point phase place phase
With extrapolation seismic response point, take in 6 seismic response points that main profile i reduces direction, maximum coherent value and its corresponding
Inclination angle, to the direction being located on same cross-track, main profile i is reduced, on automatic tracing main profile adjacent with current main profile,
With seed point phase place identical extrapolation seismic response point.
4. the method for claim 1, it is characterised in that on the current main profile of the basis all seed points relevant and
Inclination angle, drives all seed points on current main profile, the direction for increasing to main profile i and reducing, and automatic tracing is main with current
On the adjacent main profile of survey line and seed point phase place identical extrapolation seismic response point, whether each extrapolation seismic response point is judged
Channel boundary is arrived, when judging that extrapolation seismic response point has been arrived at channel boundary, has been stopped following the trail of, obtains on adjacent main profile
All seed points, until the seed point in the range of river course work area on all main profiles is obtained, including：
According to the relevant and inclination angle of current main profile coboundary seed point, border seed point is driven to be chased after to adjacent main profile automatically
Track, obtains on adjacent main profile and described border seed point phase place identical seed point；
According to the relevant and inclination angle of non-border seed point on current main profile, drive non-border seed point automatic to adjacent main profile
Follow the trail of, obtain on adjacent main profile and described non-border seed point phase place identical seed point.
5. method as claimed in claim 4, it is characterised in that the current main profile coboundary seed point of the basis relevant and
Inclination angle, drives border seed point to adjacent main profile automatic tracing, obtain on adjacent main profile with the border seed point phase
Position identical seed point, including：
Judge that the extrapolation seismic response point of border seed point is located at the position on river course；
When the extrapolation seismic response point for judging border seed point is located at river course inside, it is new with the extrapolation seismic response point
Seed point, on same main profile, the cross-track j direction automatic tracing phase place identical extrapolation seismic response that increases or reduce
Point, until judgement extrapolation seismic response point has been arrived at channel boundary, stops following the trail of, obtains the seed point on adjacent main profile.
6. method as claimed in claim 5, it is characterised in that the current main profile coboundary seed point of the basis relevant and
Inclination angle, drives border seed point to adjacent main profile automatic tracing, obtain on adjacent main profile with the border seed point phase
Position identical seed point, also includes：
When the extrapolation seismic response point for judging border seed point is located at channel boundary or outside, determine that the border seed point is
Invalid seed point, redefines germinative seed point, with germinative seed point as border seed point, to adjacent main profile automatic tracing phase
Position identical extrapolation seismic response point, obtains the seed point on adjacent main profile；The germinative seed point is：With last determination
Invalid seed point is adjacent and its extrapolation seismic response point is located inside river course.
7. the method as described in claim 1 to 6 any claim, it is characterised in that the judgement is extrapolated seismic response point
Arrive at channel boundary, specially：
When extrapolation seismic response point and with the extrapolation adjacent 8 seismic response points of seismic response point, the phase of totally 9 seismic response points
Dry mean value is little, when variance yields is big, it was demonstrated that arrived 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 CN104808244A (en) | 2015-07-29 |
| CN104808244B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
Family Cites Families (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 |
| CN103592681B (en) * | 2013-09-16 | 2016-05-04 | 电子科技大学 | A kind of seismic image tracing of horizons method based on signal classification |
| CN104181596B (en) * | 2014-08-27 | 2017-01-11 | 中国石油集团东方地球物理勘探有限责任公司 | Geologic horizon automatic tracking method and device |
-
2015
- 2015-04-14 CN CN201510173944.3A patent/CN104808244B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104808244A (en) | 2015-07-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Riedel et al. | Seismic imaging of a fractured gas hydrate system in the Krishna–Godavari Basin offshore India | |
| Paull et al. | High-resolution bathymetry of the axial channels within Monterey and Soquel submarine canyons, offshore central California | |
| Stuevold et al. | Polygonal faults in the Ormen Lange field, Møre basin, offshore mid Norway | |
| Mackay et al. | Cold‐based debris‐covered glaciers: evaluating their potential as climate archives through studies of ground‐penetrating radar and surface morphology | |
| CN105319585B (en) | A kind of method hidden using thin-layers interference amplitude recovery identification of hydrocarbon | |
| Tasianas et al. | High-resolution 3D seismic study of pockmarks and shallow fluid flow systems at the Snøhvit hydrocarbon field in the SW Barents Sea | |
| BR102016003517A2 (en) | amplitude versus angle analysis for quantitative interpretation | |
| Arnold et al. | Modeled subglacial water flow routing supports localized intrusive heating as a possible cause of basal melting of Mars' south polar ice cap | |
| Liu et al. | Characterization of architectural elements of Ordovician fractured-cavernous carbonate reservoirs, Tahe Oilfield, China | |
| Glas et al. | Hydrogeology of an alpine talus aquifer: Cordillera Blanca, Peru | |
| Rizzo et al. | Deep electrical resistivity tomography for the hydrogeological setting of Muro Lucano Mounts Aquifer (Basilicata, Southern Italy) | |
| Ahmad et al. | High-resolution seismic and resistivity profiling of a buried Quaternary subglacial valley: Northern Alberta, Canada | |
| Kunath et al. | Episodic venting of a submarine gas seep on geological time scales: Formosa Ridge, northern South China Sea | |
| Gómez-Ortiz et al. | The internal structure of modern barchan dunes of the Ebro River Delta (Spain) from ground penetrating radar | |
| CN104808244B (en) | Ancient stream channel River Bank Erosion depth automatic identifying method | |
| Chen et al. | Occurrence of submarine canyons, sediment waves and mass movements along the northern continental slope of the South China Sea | |
| Hemmings-Sykes | The influence of faulting on hydrocarbon migration in the Kupe area, south Taranaki Basin, New Zealand | |
| McBride et al. | A structural study of thermal tufas using ground-penetrating radar | |
| 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 | |
| Kaminskas et al. | New evidence of an early Pridoli barrier reef in the southern part of the Baltic Silurian basin based on three-dimensional seismic survey, Lithuania | |
| CN111856579B (en) | Method for comprehensively identifying basin edge sediment sector in deep water area of oil-gas-containing basin in south China sea | |
| Martinez et al. | Utility of ground-penetrating radar in near-surface, high-resolution imaging of Lansing-Kansas City (Pennsylvanian) limestone reservoir analogs | |
| Tibbetts et al. | Seismic-refraction measurements in Jackson hole, Wyoming | |
| Pinghua et al. | Comprehensive Dissection of Meandering-River Sedimentary System in Qinhuangdao 32-6 Oilfield, Bohai Bay Basin | |
| Sollie | Controls on hydrocarbon column-heights in the north-eastern North Sea |
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 |