CN104199107B - Depth prediction approach and system before brill based on the many wave datum of vertical seismic - Google Patents
Depth prediction approach and system before brill based on the many wave datum of vertical seismic Download PDFInfo
- Publication number
- CN104199107B CN104199107B CN201410383847.2A CN201410383847A CN104199107B CN 104199107 B CN104199107 B CN 104199107B CN 201410383847 A CN201410383847 A CN 201410383847A CN 104199107 B CN104199107 B CN 104199107B
- Authority
- CN
- China
- Prior art keywords
- time
- wave
- depth
- function
- lineups
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides depth prediction approach and system before a kind of brill based on the many wave datum of vertical seismic, and described method includes: gather the many wave datum of vertical seismic;The described many wave datum of vertical seismic are processed, the many wave datum of vertical seismic after being processed;Obtained well surface seismic section;Position and the wave group feature of target zone is identified well surface seismic section from described crossing;Direct wave corresponding to described target zone, converted wave and reflection line-ups are identified in position according to described target zone and the many wave datum of wave group feature vertical seismic after described process;Time depth function is set up according to described direct wave, converted wave and reflection line-ups;Predetermined depth and predicted time is determined according to described time depth function.Solve the problem that existing Forecasting Methodology precision is inadequate, in time, accurately for depth data before the brill of drilling engineering offer formation at target locations.
Description
Technical field
The present invention is about technical field of geophysical exploration, especially with regard to ripple application treatment technologies many in vertical seismic data,
It is concretely depth prediction approach and system before a kind of brill based on the many wave datum of vertical seismic.
Background technology
Utilizing vertical seismic data to carry out boring front depth prediction is a traditional geophysical exploration technology, it almost along with
The birth of vertical seismic data observation method and occur.Even to this day, vertical seismic data are utilized to carry out boring front depth prediction still
It it is so the important earth physical means obtaining drill bit front target zone information.
In recent years, vertical seismic data technique obtains quick development, but it bores front depth prediction approach but almost without appointing
What is progressive.In prior art, the method that only Zhu Guangming in 1986 is previously mentioned in " vertical seismic method " book has
Representativeness, its method particularly includes: even up on NMO section vertically shaking school, by first arrival and target zone reflection line-ups
With fitting a straight line and downwardly extend, the degree of depth corresponding to two straight line joints be exactly prediction the degree of depth.
The method being previously mentioned in " vertical seismic method " book is widely used in field of geophysical exploration for many years, solves
Depth prediction problem before a lot of actual brills, but the method there is also following defect.
(1), the operation of the method itself need to rely on perusal to determine straight-line extension direction, it exists the biggest the most true
Qualitative, operating error can be produced;
(2), the method only uplink and downlink compressional wave to target zone carries out curve tracing, and can only obtain a joint, especially
It is in the bigger area of stratigraphic dip, and descending compressional wave and up compressional wave angle in vertical seismic data section are the least, cross
The point degree of depth is more difficult to determine, it was predicted that precision receives severely restricts.
Therefore, how to propose a kind of brand-new vertical seismic data many ripples depth prediction approach, before it can improve brill, the degree of depth is pre-
Survey precision, provide reliable data support to be this area technical barrier urgently to be resolved hurrily to drilling engineering.
Summary of the invention
For the above-mentioned technological deficiency overcoming vertical seismic data many depth of convolutions degree Predicting Technique of the prior art to exist, the present invention
Provide depth prediction approach and system before a kind of brill based on the many wave datum of vertical seismic, by the many wave datum of vertical seismic
Following the trail of the multiple lineups relevant to target zone in wave field, the most multiple reflection line-ups set up multiple time depth letter simultaneously
Number, and then obtain predetermined depth and predicted time, decrease the occasionality predicted the outcome, precision of prediction is significantly increased.
It is an object of the invention to provide depth prediction approach before a kind of brill based on the many wave datum of vertical seismic, including:
Gather the many wave datum of vertical seismic;Processing the described many wave datum of vertical seismic, the vertical seismic after being processed is many
Wave datum;Obtained well surface seismic section;From described cross identify well surface seismic section target zone position and
Wave group feature;Position according to described target zone and the many wave datum of wave group feature vertical seismic after described process identify
Go out direct wave corresponding to described target zone, converted wave and reflection line-ups;According to described direct wave, converted wave and
Reflection line-ups sets up time depth function;When determining predetermined depth and prediction according to described time depth function
Between.
It is an object of the invention to provide depth prediction system before a kind of brill based on the many wave datum of vertical seismic, including:
Seismic data acquisition device, is used for gathering the many wave datum of vertical seismic;Seismic data process device, for described vertical
Seismic Multi-wave data process, the many wave datum of vertical seismic after being processed;Seismic profile acquisition device, is used for obtaining
Cross well surface seismic section;Target zone identification device, for identifying target zone from described crossing well surface seismic section
Position and wave group feature;Reflection line-ups identification device, is used for the position according to described target zone and wave group feature from institute
State and the many wave datum of the vertical seismic after process identify direct wave corresponding to described target zone, converted wave and reflected wave in phase
Axle;Time deep function set up device, for set up the time according to described direct wave, converted wave and reflection line-ups deep
Degree function;Device is deeply determined, for determining predetermined depth and predicted time according to described time depth function during prediction.
The beneficial effects of the present invention is, it is provided that depth prediction approach and be before a kind of brill based on the many wave datum of vertical seismic
System, by following the trail of the multiple lineups relevant to target zone, the most multiple reflections in vertical seismic many wave datum wave field simultaneously
Lineups set up multiple time depth function, and then obtain predetermined depth and predicted time, decrease predict the outcome accidental
Property, precision of prediction is significantly increased, improves depth prediction precision before brill, provide reliable data support to drilling engineering, solve
The problem that certainly existing Forecasting Methodology precision is inadequate.
For the above and other objects, features and advantages of the present invention can be become apparent, preferred embodiment cited below particularly, and
Coordinate institute's accompanying drawings, be described in detail below.
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing skill
In art description, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only the present invention
Some embodiments, for those of ordinary skill in the art, on the premise of not paying creative work, it is also possible to root
Other accompanying drawing is obtained according to these accompanying drawings.
The flow process of the front depth prediction approach of a kind of based on the many wave datum of vertical seismic brill that Fig. 1 provides for the embodiment of the present invention
Figure;
Fig. 2 is the particular flow sheet of step S102 in Fig. 1;
Fig. 3 is the particular flow sheet of step S106 in Fig. 1;
Fig. 4 is the particular flow sheet of step S301 in Fig. 3;
Fig. 5 is the particular flow sheet of step S302 in Fig. 3;
Fig. 6 is the particular flow sheet of step S107 in Fig. 1;
The structural frames of the front depth prediction system of a kind of based on the many wave datum of vertical seismic brill that Fig. 7 provides for the embodiment of the present invention
Figure;
A kind of based on the many wave datum of vertical seismic the earthquake bored in front depth prediction system that Fig. 8 provides for the embodiment of the present invention
The concrete structure block diagram of data processing equipment 200;
Before a kind of based on the many wave datum of vertical seismic the brill that Fig. 9 provides for the embodiment of the present invention in depth prediction system time deep
Function sets up the concrete structure block diagram of device 600;
Before a kind of based on the many wave datum of vertical seismic the brill that Figure 10 provides for the embodiment of the present invention in depth prediction system time
The concrete structure block diagram of deep sequence pickup model 601;
Before a kind of based on the many wave datum of vertical seismic the brill that Figure 11 provides for the embodiment of the present invention in depth prediction system time
Deep function sets up the concrete structure block diagram of module 602;
Pre-in depth prediction system before a kind of based on the many wave datum of vertical seismic the brill that Figure 12 provides for the embodiment of the present invention
The concrete structure block diagram of device 700 is deeply determined during survey;
Figure 13 is depth prediction principle schematic before the vertical seismic many ripples of many wave datum bore;
Figure 14 is many ripples type schematic diagram of vertical seismic multimode type data;
Figure 15 is that the Tarim Basin vertical seismic data VSP many ripples of well logging bore VSP section after front depth prediction example polarized rotation
Schematic diagram;
Figure 16 is the partial enlarged drawing that before the VSP many ripples of well logging in Tarim Basin bore, depth prediction example predicts the A point that crosses;
Figure 17 is the partial enlarged drawing that before the VSP many ripples of well logging in Tarim Basin bore, depth prediction example predicts the B point that crosses;
Figure 18 is depth prediction example schematic before the VSP many ripples of well logging in Caspian Sea region, shore bore.
Detailed description of the invention
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clearly and completely
Describe, it is clear that described embodiment is only a part of embodiment of the present invention rather than whole embodiments.Based on this
Embodiment in invention, the every other reality that those of ordinary skill in the art are obtained under not making creative work premise
Execute example, broadly fall into the scope of protection of the invention.
It is an object of the invention to the wave field information effectively utilizing VSP abundant, precognition bores the front target zone degree of depth, solves existing prediction
The problem that method precision is inadequate, in time, accurately for depth data before the brill of drilling engineering offer formation at target locations.
Fig. 1 is the particular flow sheet of the front depth prediction approach of a kind of based on the many wave datum of vertical seismic brill that the present invention proposes,
As shown in Figure 1, described method includes:
S101: gather the many wave datum of vertical seismic.
In the particular embodiment, available existing instrument collects the many wave datum of vertical seismic.What the present invention was gathered hangs down
Straight Seismic Multi-wave data are broadly divided into two groups, and first group is descending compressional wave (P ripple), up compressional wave (PP ripple), up conversion
Shear wave (PPs ripple), its common ground is that to arrive at the down going wave of target zone be compressional wave;Second group is descending shear wave (S ripple or descending
Converted shear wave Ps ripple), the up reflection wave (PsPs ripple) of shear wave, the up conversion compressional wave (PsP ripple) of shear wave, it is altogether
It is that to arrive at the down going wave of target zone be shear wave with point.
S102: the described many wave datum of vertical seismic are processed, the many wave datum of vertical seismic after being processed, Fig. 2
Particular flow sheet for step S102.
S103: obtained well surface seismic section.In the particular embodiment, available existing instrument collected well ground
Face seismic profile.
S104: identify position and the wave group feature of target zone from described crossing well surface seismic section.
S105: identify according to the position of described target zone and the many wave datum of wave group feature vertical seismic after described process
Go out the reflection line-ups that described target zone is corresponding.In a particular embodiment, can be according to seismic data interpretation principle vertically
Seismic Multi-wave data wave field identifies and corresponding to this target zone, respectively organizes lineups.Direct wave that the present invention is identified, turn
Change ripple and reflection line-ups includes that descending through compressional wave lineups, up reflected P-wave lineups, up converted shear wave are same
Phase axle, descending through shear wave or downward converted P-SV wave lineups, the up reflection wave lineups of shear wave, up turn of shear wave
Change compressional wave lineups.
S106: set up time depth function according to described reflection line-ups.Fig. 3 is the particular flow sheet of step S106.
S107: determine predetermined depth and predicted time according to described time depth function.Fig. 6 is step S107
Particular flow sheet.
Fig. 2 is the particular flow sheet of step S102, and as shown in Figure 2, this step specifically includes:
S401: the described many wave datum of vertical seismic are carried out dephasing process, to ensure that each group of wave field all can follow the trail of peak
Be worth to reliable time depth relation curve, the step for rely on conventional seismic processing software to realize.
S402: the many wave datum of vertical seismic after processing dephasing carry out noise attentuation process;
S403: the many wave datum of vertical seismic after processing noise attentuation carry out energy compensating process.
S404: the many wave datum of vertical seismic after processing energy compensating carry out spin orientation process.
S405: the many wave datum of vertical seismic after processing spin orientation carry out deconvolution process, after being processed vertically
Shake many wave datum.
In a particular embodiment, available existing method and software carry out noise attentuation, energy compensating, spin orientation,
Deconvolution processes, and obtains the various wave fields many wave datum of VSP the most clearly in section.Each processing links described can not change
Phase place, it is necessary to assure data are always zero phase.
Fig. 3 is the particular flow sheet of step S106, from the figure 3, it may be seen that this step specifically includes:
S301: the crest of direct wave, converted wave and reflection line-ups described in pickup, obtains time depth relation sequence
Row, Fig. 4 is the particular flow sheet of step S301, and as shown in Figure 4, this step specifically includes:
S401: pick up the crest of descending through compressional wave lineups, obtain very first time depth relationship sequence;
S402: pick up the crest of up reflected P-wave lineups, obtains the second time depth relational sequence;
S403: pick up the crest of up converted shear wave lineups, obtain the 3rd time depth relational sequence;
S404: pick up descending through shear wave lineups or the crest of downward converted P-SV wave lineups, obtain the 4th time depth
Relational sequence, described descending through shear wave or downward converted P-SV wave producing method are different, but communication effect in the earth formation
Identical with section performance, therefore merge into same feature at this and be analyzed;
S405: the crest of the up reflection wave lineups of pickup shear wave, obtains the 5th time depth relational sequence;
S406: the crest of the up conversion compressional wave lineups of pickup shear wave, obtains the 6th time depth relational sequence.
There is difference in various lineups developmental state in real data, the time depth relational sequence of its correspondence not necessarily may be used
All to obtain.Following steps are all to illustrate under the obtainable supposed premise of total data, when individual wave fields cannot be distinguished
Not, when causing time depth relational sequence to obtain, the computational accuracy of the present invention can be affected on a small quantity, but it realizes
Method is identical.
From the figure 3, it may be seen that step S106 also includes:
S302: set up time depth function according to described time depth relational sequence.Fig. 5 is the tool of step S302
Body flow chart, as shown in Figure 5, this step specifically includes:
S501: set up the time depth letter of descending through compressional wave lineups according to described very first time depth relationship sequence
Number, referred to as very first time depth function, in a particular embodiment, available discrete data quadratic polynomial matching is former
Reason obtains the time depth function curve of P ripple lineups.Described discrete data quadratic polynomial approximating method is ripe skill
Art, available various digital processing softwares realize, and very first time depth function is shown below:
Y=apX2+bpX+cp
Wherein: Y is reflex time, X is the corresponding degree of depth, ap、bp、cpTime depth curve matching for P ripple
Multinomial coefficient, these coefficients can be obtained by the software of any support data fitting operation.
S502: set up the time depth function of up reflected P-wave lineups according to the second time depth relational sequence, claims
It it is the second time depth function.
In a particular embodiment, available discrete data quadratic polynomial fitting theory obtains the time of PP ripple lineups
Depth function curve.Described discrete data quadratic polynomial approximating method is mature technology, available various digital processings
Software realizes, and the second time depth function is shown below:
Y=appX2+bppX+cpp
Wherein: Y is reflex time, X is the corresponding degree of depth, app、bpp、cppTime depth curve for PP ripple is intended
The multinomial coefficient closed, these coefficients can be obtained by the software of any support data fitting operation.
S503: set up the time depth function of up converted shear wave lineups according to the 3rd time depth relational sequence, claims
Being the 3rd time depth function, in a particular embodiment, available discrete data quadratic polynomial fitting theory obtains
The time depth function curve of PPs ripple lineups.Described discrete data quadratic polynomial approximating method is mature technology,
Available various digital processing softwares realize, and the 3rd time depth function is shown below:
Y=appsX2+bppsX+cpps
Wherein: Y is reflex time, X is the corresponding degree of depth, apps、bpps、cppsTime depth curve for PPs ripple
The multinomial coefficient of matching, these coefficients can be obtained by the software of any support data fitting operation.
S504: set up the time depth function of descending shear wave lineups according to the 4th time depth relational sequence, is referred to as
4th time depth function, the most descending through shear wave of time depth function of descending shear wave lineups or downward converted P-SV wave
The time depth function of lineups, in a particular embodiment, available discrete data quadratic polynomial fitting theory obtains
To S ripple or the time depth function curve of Ps ripple lineups.Described discrete data quadratic polynomial approximating method is for becoming
Cooking technique, available various digital processing softwares realize, and the 4th time depth function is shown below:
Y=asX2+bsX+csOr Y=apsX2+bpsX+cps
Wherein: Y is reflex time, X is the corresponding degree of depth, as、bs、cs、aps、bps、cpsFor S ripple or Ps ripple
The multinomial coefficient of time depth curve matching, these coefficients can be obtained by the software of any support data fitting operation
Arrive.
S505: set up the time depth letter of the up reflection wave lineups of shear wave according to the 5th time depth relational sequence
Number, the referred to as the 5th time depth function, in a particular embodiment, available discrete data quadratic polynomial matching is former
Reason obtains the time depth function curve of PsPs ripple lineups.Described discrete data quadratic polynomial approximating method is ripe
Technology, available various digital processing softwares realize, and the 5th time depth function is shown below:
Y=apspsX2+bpspsX+cpsps
Wherein: Y is reflex time, X is the corresponding degree of depth, apsps、bpsps、cpspsTime depth for PsPs ripple is bent
The multinomial coefficient of line matching, these coefficients can be obtained by the software of any support data fitting operation.
S506: set up the time depth letter of the up conversion compressional wave lineups of shear wave according to the 6th time depth relational sequence
Number, the referred to as the 6th time depth function.In a particular embodiment, available discrete data quadratic polynomial matching is former
Reason obtains the time depth function curve of PsP ripple lineups.Described discrete data quadratic polynomial approximating method is ripe
Technology, available various digital processing softwares realize, and the 6th time depth function is shown below:
Y=apspX2+bpspX+cpsp
Wherein: Y is reflex time, X is the corresponding degree of depth, apsp、bpsp、cpspTime depth curve for PsP ripple
The multinomial coefficient of matching, these coefficients can be obtained by the software of any support data fitting operation.
Fig. 6 is the particular flow sheet of step S107 in Fig. 1, it will be appreciated from fig. 6 that this step specifically includes:
S601: according to very first time depth function and the second time depth function obtain the first depth prediction result,
One time predicted the outcome, and the most in a particular embodiment, solved very first time depth function and the second time depth
The binary quadratic equation group of function composition, i.e. can get the first depth prediction result X1, the very first time predicts the outcome Y1。
S602: according to very first time depth function and the 3rd time depth function obtain the second depth prediction result,
Two time prediction results, the most in a particular embodiment, solve very first time depth function and the 3rd time depth
The binary quadratic equation group of function composition, i.e. can get the second depth prediction result X2, the second time prediction result Y2;
S603: according to the second time depth function and the 3rd time depth function obtain the 3rd depth prediction result,
Three-time estimates result, the most in a particular embodiment, solves the second time depth function and the 3rd time depth
The binary quadratic equation group of function composition, i.e. can get the 3rd depth prediction result X3, three-time estimates result Y3;
S604: according to the 4th time depth function and the 5th time depth function obtain the 4th depth prediction result,
Four time prediction results, the most in a particular embodiment, solve the 4th time depth function and the 5th time depth
The binary quadratic equation group of function composition, i.e. can get the 4th depth prediction result X4, the 4th time prediction result Y4;
S605: according to the 4th time depth function and the 6th time depth function obtain the 5th depth prediction result,
Five time prediction results, the most in a particular embodiment, solve the 4th time depth function and the 6th time depth
The binary quadratic equation group of function composition, i.e. can get the 5th depth prediction result X5, the 5th time prediction result Y5;
S606: according to the 5th time depth function and the 6th time depth function obtain the 6th depth prediction result,
Six time prediction results, the most in a particular embodiment, solve the 5th time depth function and the 6th time depth
The binary quadratic equation group of function composition, i.e. can get the 6th depth prediction result X6, the 6th time prediction result Y6;
The negative result not having physical significance has all been given up during above-mentioned solving equations.
S607: to the first described depth prediction result, the second depth prediction result, the 3rd depth prediction result, the 4th deep
Degree predicts the outcome, the 5th depth prediction result, the 6th depth prediction result are averaged value computing, obtain predetermined depth Xfinal,
As follows:The wherein number of the n depth prediction result by being tried to achieve.
S608: predict the outcome the described very first time, the second time prediction result and three-time estimates result are put down
Mean operation, obtains compressional wave predicted time Yfinal, as follows:Wherein n by tried to achieve first to
The number of three-time estimates result.Pre-to the 4th described time prediction result, the 5th time prediction result, the 6th time
Survey result to be averaged value computing, obtain shear wave predicted time Y 'final, as follows:Wherein m is institute
The number of the 4th to the 6th time prediction result tried to achieve.
As it has been described above, be depth prediction approach before a kind of based on the many wave datum of vertical seismic the brill that the present invention provides, pass through
Following the trail of the multiple lineups relevant to target zone in vertical seismic many wave datum wave field, the most multiple reflection line-ups are built simultaneously
Vertical multiple time depth functions, and then obtain predetermined depth and predicted time, decrease the occasionality predicted the outcome, significantly
Increase precision of prediction.
The structural frames of the front depth prediction system of a kind of based on the many wave datum of vertical seismic brill that Fig. 7 provides for the embodiment of the present invention
Figure, as shown in Figure 7, described system includes:
Seismic data acquisition device 100, is used for gathering the many wave datum of vertical seismic.
In the particular embodiment, available existing instrument collects the many wave datum of vertical seismic.What the present invention was gathered hangs down
Straight Seismic Multi-wave data are broadly divided into two groups, and first group is descending compressional wave (P ripple), up compressional wave (PP ripple), up conversion
Shear wave (PPs ripple), its common ground is that to arrive at the down going wave of target zone be compressional wave;Second group is descending shear wave (S ripple or descending
Converted shear wave Ps ripple), the up reflection wave (PsPs ripple) of shear wave, the up conversion compressional wave (PsP ripple) of shear wave, it is altogether
It is that to arrive at the down going wave of target zone be shear wave with point.
Seismic data process device 200, for the described many wave datum of vertical seismic are processed, hanging down after being processed
Straight Seismic Multi-wave data, Fig. 8 is the concrete structure block diagram of seismic data process device 200.
Seismic profile acquisition device 300, is used for obtaining well surface seismic section.In the particular embodiment, available existing
Instrument is had to collect well surface seismic section.
Target zone identification device 400, for identifying position and the ripple of target zone from described crossing well surface seismic section
Stack features.
Reflection line-ups identification device 500, is used for the position according to described target zone and wave group feature after described process
The many wave datum of vertical seismic identify the reflection line-ups that described target zone is corresponding.In a particular embodiment, can basis
Seismic data interpretation principle identifies in vertical seismic many wave datum wave field respectively organizes lineups corresponding to this target zone.This
Bright identified direct wave, converted wave and reflection line-ups include descending through compressional wave lineups, up longitudinal wave reflection
Lineups, up converted shear wave lineups, descending through shear wave lineups, the up reflection wave reflection line-ups of shear wave,
The up conversion compressional wave lineups of shear wave.
Time deep function set up device 600, for setting up time depth function according to described reflection line-ups.When Fig. 9 is
Deep function sets up the concrete structure block diagram of device 600.
Device 700 is deeply determined, for determining predetermined depth and predicted time according to described time depth function during prediction.
Figure 12 is the concrete structure block diagram deeply determining device 700 during prediction.
Fig. 8 is the concrete structure block diagram of seismic data process device 200, as shown in Figure 8, seismic data process device 200
Specifically include:
Dephasing processing module 201, for carrying out dephasing process to the described many wave datum of vertical seismic, to ensure
Each group wave field all can be followed the trail of peak value and obtain reliable time depth relation curve, the step for rely on conventional seismic processing
Software can realize.
Noise attentuation processing module 202, the many wave datum of vertical seismic after processing dephasing are carried out at noise attentuation
Reason;
Energy compensating processing module 203, the many wave datum of vertical seismic after processing noise attentuation are carried out at energy compensating
Reason.
Spin orientation processing module 204, the many wave datum of vertical seismic after processing energy compensating are carried out at spin orientation
Reason.
Deconvolution processing module 205, the many wave datum of vertical seismic after processing spin orientation carry out deconvolution process,
The many wave datum of vertical seismic after being processed.
In a particular embodiment, available existing method and software carry out noise attentuation, energy compensating, spin orientation,
Deconvolution processes, and obtains the various wave fields many wave datum of VSP the most clearly in section.Each processing links described can not change
Phase place, it is necessary to assure data are always zero phase.
Fig. 9 when being deep function set up the concrete structure block diagram of device 600, as shown in Figure 9, time deep function set up device 600
Specifically include:
Time deep sequence pickup model 601, for picking up the crest of described direct wave, converted wave and reflection line-ups,
Obtaining time depth relational sequence, Figure 10 is the concrete structure block diagram of deep sequence pickup model 601 when being, as shown in Figure 10,
This step specifically includes:
Deep sequence pickup unit 6011 when first, for picking up the crest of descending through compressional wave lineups, obtains the very first time
Depth relationship sequence;
Deep sequence pickup unit 6012 when second, for picking up the crest of up reflected P-wave lineups, obtained for the second time
Depth relationship sequence;
Deep sequence pickup unit 6013 when the 3rd, for picking up the crest of up converted shear wave lineups, obtained for the 3rd time
Depth relationship sequence;
Deep sequence pickup unit 6014 when the 4th, for picking up descending through shear wave lineups or downward converted P-SV wave lineups
Crest, obtains the 4th time depth relational sequence;
Deep sequence pickup unit 6015 when the 5th, for picking up the crest of the up reflection wave lineups of shear wave, obtains the 5th
Time depth relational sequence;
Deep sequence pickup unit 6016 when the 6th, for picking up the crest of the up conversion compressional wave lineups of shear wave, obtains the 6th
Time depth relational sequence.
There is difference in various lineups developmental state in real data, the time depth relational sequence of its correspondence not necessarily may be used
All to obtain.Following steps are all to illustrate under the obtainable supposed premise of total data, when individual wave fields cannot be distinguished
Not, when causing time depth relational sequence to obtain, the computational accuracy of the present invention can be affected on a small quantity, but it realizes
Method is identical.
As shown in Figure 9, time, deep function is set up device 600 and is also included:
Time deep function set up module 602, for setting up time depth function according to described time depth relational sequence.
Figure 11 when being deep function set up the concrete structure block diagram of module 602, as shown in Figure 11, time deep function set up module 602 and have
Body includes:
When first, function sets up unit 6021, descending through vertical for setting up according to described very first time depth relationship sequence
The time depth function of ripple lineups, referred to as very first time depth function, in a particular embodiment, may utilize from
Dissipate data quadratic polynomial fitting theory and obtain the time depth function curve of P ripple lineups.Described discrete data secondary
Polynomial fitting method is mature technology, and available various digital processing softwares realize, very first time depth function such as following formula
Shown in:
Y=apX2+bpX+cp
Wherein: Y is reflex time, X is the corresponding degree of depth, ap、bp、cpTime depth curve matching for P ripple
Multinomial coefficient, these coefficients can be obtained by the software of any support data fitting operation.
When second, function sets up unit 6022, for setting up up reflected P-wave homophase according to the second time depth relational sequence
The time depth function of axle, the referred to as second time depth function.
In a particular embodiment, available discrete data quadratic polynomial fitting theory obtains the time of PP ripple lineups
Depth function curve.Described discrete data quadratic polynomial approximating method is mature technology, available various digital processings
Software realizes, and the second time depth function is shown below:
Y=appX2+bppX+cpp
Wherein: Y is reflex time, X is the corresponding degree of depth, app、bpp、cppTime depth curve for PP ripple is intended
The multinomial coefficient closed, these coefficients can be obtained by the software of any support data fitting operation.
When the 3rd, function sets up unit 6023, for setting up up converted shear wave homophase according to the 3rd time depth relational sequence
The time depth function of axle, the referred to as the 3rd time depth function, in a particular embodiment, available discrete data
Quadratic polynomial fitting theory obtains the time depth function curve of PPs ripple lineups.Described discrete data secondary is multinomial
Formula approximating method is mature technology, and available various digital processing softwares realize, and the 3rd time depth function is shown below:
Y=appsX2+bppsX+cpps
Wherein: Y is reflex time, X is the corresponding degree of depth, apps、bpps、cppsTime depth curve for PPs ripple
The multinomial coefficient of matching, these coefficients can be obtained by the software of any support data fitting operation.
When the 4th, function sets up unit 6024, for setting up descending shear wave lineups according to the 4th time depth relational sequence
Time depth function, the referred to as the 4th time depth function, the time depth function of descending shear wave lineups is the most descending directly
Reach the time depth function of shear wave or downward converted P-SV wave lineups, in a particular embodiment, available discrete data
Quadratic polynomial fitting theory obtains S ripple or the time depth function curve of Ps ripple lineups.Described discrete data two
Order polynomial approximating method is mature technology, and available various digital processing softwares realize, and the 4th time depth function is as follows
Shown in formula:
Y=asX2+bsX+csOr Y=apsX2+bpsX+cps
Wherein: Y is reflex time, X is the corresponding degree of depth, as、bs、cs、aps、bps、cpsFor S ripple or Ps ripple
The multinomial coefficient of time depth curve matching, these coefficients can be obtained by the software of any support data fitting operation
Arrive.
When the 5th, function sets up unit 6025, for setting up the up reflection horizontal stroke of shear wave according to the 5th time depth relational sequence
The time depth function of ripple lineups, the referred to as the 5th time depth function, in a particular embodiment, may utilize from
Dissipate data quadratic polynomial fitting theory and obtain the time depth function curve of PsPs ripple lineups.Described discrete data two
Order polynomial approximating method is mature technology, and available various digital processing softwares realize, and the 5th time depth function is as follows
Shown in formula:
Y=apspsX2+bpspsX+cpsps
Wherein: Y is reflex time, X is the corresponding degree of depth, apsps、bpsps、cpspsTime depth for PsPs ripple is bent
The multinomial coefficient of line matching, these coefficients can be obtained by the software of any support data fitting operation.
When the 6th, function sets up unit 6026, indulges for setting up the up conversion of shear wave according to the 6th time depth relational sequence
The time depth function of ripple lineups, the referred to as the 6th time depth function.In a particular embodiment, may utilize from
Dissipate data quadratic polynomial fitting theory and obtain the time depth function curve of PsP ripple lineups.Described discrete data two
Order polynomial approximating method is mature technology, and available various digital processing softwares realize, and the 6th time depth function is as follows
Shown in formula:
Y=apspX2+bpspX+cpsp
Wherein: Y is reflex time, X is the corresponding degree of depth, apsp、bpsp、cpspTime depth curve for PsP ripple
The multinomial coefficient of matching, these coefficients can be obtained by the software of any support data fitting operation.
Figure 12 is deeply to determine the concrete structure block diagram of device 700, as shown in Figure 12 during prediction, it was predicted that time deeply determine device 700
Specifically include:
Module 701 is deeply determined, for obtaining according to very first time depth function and the second time depth function when first
First depth prediction result, the very first time predict the outcome, and the most in a particular embodiment, solve very first time degree of depth letter
Number and the binary quadratic equation group of the second time depth function composition, i.e. can get the first depth prediction result X1, first
Time prediction result Y1。
Module 702 is deeply determined, for obtaining according to very first time depth function and the 3rd time depth function when second
Second depth prediction result, the second time prediction result, the most in a particular embodiment, solve very first time degree of depth letter
Number and the binary quadratic equation group of the 3rd time depth function composition, i.e. can get the second depth prediction result X2, second
Time prediction result Y2;
Module 703 is deeply determined, for obtaining according to the second time depth function and the 3rd time depth function when the 3rd
3rd depth prediction result, three-time estimates result, the most in a particular embodiment, solve the second time depth letter
Number and the binary quadratic equation group of the 3rd time depth function composition, i.e. can get the 3rd depth prediction result X3, the 3rd
Time prediction result Y3;
Module 704 is deeply determined, for obtaining according to the 4th time depth function and the 5th time depth function when the 4th
4th depth prediction result, the 4th time prediction result, the most in a particular embodiment, solve the 4th time depth letter
Number and the binary quadratic equation group of the 5th time depth function composition, i.e. can get the 4th depth prediction result X4, the 4th
Time prediction result Y4;
Module 705 is deeply determined, for obtaining according to the 4th time depth function and the 6th time depth function when the 5th
5th depth prediction result, the 5th time prediction result, the most in a particular embodiment, solve the 4th time depth letter
Number and the binary quadratic equation group of the 6th time depth function composition, i.e. can get the 5th depth prediction result X5, the 5th
Time prediction result Y5;
Module 706 is deeply determined, for obtaining according to the 5th time depth function and the 6th time depth function when the 6th
6th depth prediction result, the 6th time prediction result, the most in a particular embodiment, solve the 5th time depth letter
Number and the binary quadratic equation group of the 6th time depth function composition, i.e. can get the 6th depth prediction result X6, the 6th
Time prediction result Y6;
The negative result not having physical significance has all been given up during above-mentioned solving equations.
Predetermined depth determines module 707, for the first described depth prediction result, the second depth prediction result, the 3rd
Depth prediction result, the 4th depth prediction result, the 5th depth prediction result, the 6th depth prediction result be averaged value fortune
Calculate, obtain predetermined depth Xfinal, as follows:Wherein n is the depth prediction result tried to achieve
Number.
Predicted time determines module 708, for predicting the outcome the described very first time, the second time prediction result and the
Three-time estimates result is averaged value computing, obtains compressional wave predicted time Yfinal, as follows:Its
Middle n by tried to achieve first to the number of three-time estimates result.To the 4th described time prediction result, the 5th time
Predict the outcome, the 6th time prediction result is averaged value computing, obtains shear wave predicted time Y 'final, as follows:The wherein number of the m the 4th to the 6th time prediction result by being tried to achieve.
As it has been described above, be depth prediction system before a kind of based on the many wave datum of vertical seismic the brill that the present invention provides, pass through
Following the trail of the multiple lineups relevant to target zone in vertical seismic many wave datum wave field, the most multiple reflection line-ups are built simultaneously
Vertical multiple time depth functions, and then obtain predetermined depth and predicted time, decrease the occasionality predicted the outcome, significantly
Increase precision of prediction.
Below in conjunction with specific embodiment, technical scheme is discussed in detail.
In order to improve the front depth prediction precision of brill, providing reliable data support to drilling engineering, the present invention proposes a kind of complete
New VSP many ripples depth prediction approach.Follow the trail of the multiple lineups relevant to target zone, then in VSP wave field simultaneously
Utilize many group traveltime-distance equation to try to achieve a series of prediction data, and then statistics obtains predetermined depth.The application of the method is permissible
Reduce the occasionality of result of calculation, precision of prediction is significantly increased.Figure 13 is that before the vertical seismic many ripples of many wave datum bore, the degree of depth is pre-
Surveying principle schematic, Figure 14 is many ripples type schematic diagram of vertical seismic multimode type data, and Figure 15 is that Tarim Basin is vertical
The geological data VSP many ripples of well logging bore VSP generalized section after front depth prediction example polarized rotation, and Figure 16 is Tarim Basin
The VSP many ripples of well logging bore front depth prediction example and predict the partial enlarged drawing of the A point that crosses, and Figure 17 is Tarim Basin VSP well logging
Many ripples bore front depth prediction example and predict the partial enlarged drawing of the B point that crosses, and Figure 18 is Caspian Sea region, the shore VSP many ripples of well logging
Depth prediction example schematic before boring.
Being the split section of two wave fields shown in figure, top half is the longitudinal wave field after wave field separation, and the latter half is ripple
Converted wave wave field after the separation of field, two wave fields are evened up according to descending compressional wave and down conversion preliminary wave and splice respectively, therefore spell
The place of connecing is descending through compressional wave position and down conversion ripple first arrival position, and have again in top half wave field is up vertical simultaneously
Up converted shear wave in ripple and the latter half wave field, totally four kinds of ripples, crossed by four ripples, has been calculated down-hole accurately
The position of two target zones.
Figure 15 to Figure 17 show five kinds of wave fields and crosses the example of predetermined depth, and the prediction of relative Prediction distance is mistaken for 0.82%;
Figure 18 show four kinds of ripples and crosses the example of predetermined depth, and the prediction of relative Prediction distance is mistaken for 1.31%;Two examples all take
Good prediction effect.
In sum, depth prediction approach and system before a kind of based on the many wave datum of vertical seismic the brill that the present invention proposes, logical
Cross in vertical seismic many wave datum wave field, follow the trail of the multiple lineups relevant to target zone, the most multiple reflection line-ups simultaneously
Set up multiple time depth function, and then obtain predetermined depth and predicted time, decrease the occasionality predicted the outcome, greatly
Width increases precision of prediction, improves depth prediction precision before brill, provides reliable data support to drilling engineering, solves existing
The problem that Forecasting Methodology precision is inadequate.
The beneficial effects of the present invention is:
1), the present invention is applied to multiple wave field, is mutually authenticated with tracing process in identifying, can effectively reduce axis of reflection knowledge
Not, repeatedly wave interference etc. brings the error that target zone is assert;
2), for the bigger area of stratigraphic dip, the up-downgoing compressional wave that existing Forecasting Methodology is utilized homophase in VSP section
Axle clamp angle is less, it is impossible to obtain reliable depth prediction result;Meanwhile, in the VSP seismic profile in this kind of area, converted wave is but
Grow very much, hence with many ripples Forecasting Methodology in the present invention, large dip angle strata depth prediction problem can be solved;
3), the present invention through statistical average obtain finally predict the outcome, decrease the occasionality of calculating, can be significantly increased pre-
Survey precision.
One of ordinary skill in the art will appreciate that all or part of flow process realizing in above-described embodiment method, can be by meter
Calculation machine program instructs relevant hardware and completes, and described program can be stored in general computer read/write memory medium,
This program is upon execution, it may include such as the flow process of the embodiment of above-mentioned each method.Wherein, described storage medium can be magnetic disc,
CD, read-only store-memory body (Read-Only Memory, ROM) or random store-memory body (Random Access
Memory, RAM) etc..
Those skilled in the art are it will also be appreciated that the various functions that the embodiment of the present invention is listed are to be come by hardware or software
Realize depending on specifically applying the design requirement with whole system.Those skilled in the art can specifically should for every kind
With, it is possible to use the function described in the realization of various methods, but this realization is understood not to protect beyond the embodiment of the present invention
Scope.
The present invention applies specific embodiment principle and the embodiment of the present invention are set forth, saying of above example
Bright method and the core concept thereof being only intended to help to understand the present invention;Simultaneously for one of ordinary skill in the art, depend on
According to the thought of the present invention, the most all will change, in sum, in this specification
Hold and should not be construed as limitation of the present invention.
Claims (14)
1. a depth prediction approach before brill based on the many wave datum of vertical seismic, is characterized in that, described method specifically includes:
Gather the many wave datum of vertical seismic;
The described many wave datum of vertical seismic are processed, the many wave datum of vertical seismic after being processed;
Obtained well surface seismic section;
Position and the wave group feature of target zone is identified well surface seismic section from described crossing;
Institute is identified in position according to described target zone and the many wave datum of wave group feature vertical seismic after described process
State direct wave corresponding to target zone, converted wave and reflection line-ups;
Time depth function is set up according to described direct wave, converted wave and reflection line-ups;
Predetermined depth and predicted time is determined according to described time depth function.
Method the most according to claim 1, is characterized in that, processes the described many wave datum of vertical seismic,
The many wave datum of vertical seismic after being processed specifically include:
The described many wave datum of vertical seismic are carried out dephasing process;
The many wave datum of vertical seismic after processing dephasing carry out noise attentuation process;
The many wave datum of vertical seismic after processing noise attentuation carry out energy compensating process;
The many wave datum of vertical seismic after processing energy compensating carry out spin orientation process;
The many wave datum of vertical seismic after processing spin orientation carry out deconvolution process, the many ripples of vertical seismic after being processed
Data.
Method the most according to claim 1, is characterized in that, the many wave datum of described vertical seismic include descending compressional wave,
The up conversion compressional wave of up compressional wave, up converted shear wave, descending shear wave, the up reflection wave of shear wave and shear wave;Institute
Direct wave, converted wave and the reflection line-ups stated include descending through compressional wave lineups, up reflected P-wave lineups,
Up converted shear wave lineups, up reflection wave lineups, up conversion compressional wave lineups and descending through shear wave homophase
Axle or downward converted P-SV wave lineups.
Method the most according to claim 3, is characterized in that, according to described direct wave, converted wave and echo
Lineups are set up time depth function and are specifically included:
The crest of direct wave, converted wave and reflection line-ups described in pickup, obtains time depth relational sequence;
Time depth function is set up according to described time depth relational sequence.
Method the most according to claim 4, is characterized in that, direct wave, converted wave and the echo described in pickup
The crest of lineups, obtains time depth relational sequence and specifically includes:
Pick up the crest of descending through compressional wave lineups, obtain very first time depth relationship sequence;
Pick up the crest of up reflected P-wave lineups, obtain the second time depth relational sequence;
Pick up the crest of up converted shear wave lineups, obtain the 3rd time depth relational sequence;
Pick up descending through shear wave lineups or the crest of downward converted P-SV wave lineups, obtain the 4th time depth relation sequence
Row;
The crest of the up reflection wave lineups of pickup shear wave, obtains the 5th time depth relational sequence;
The crest of the up conversion compressional wave lineups of pickup shear wave, obtains the 6th time depth relational sequence.
Method the most according to claim 5, is characterized in that, when setting up according to described time depth relational sequence
Between depth function specifically include:
Set up the time depth function of descending through compressional wave lineups according to described very first time depth relationship sequence, claim
For very first time depth function;
The time depth function of up reflected P-wave lineups is set up according to the second time depth relational sequence, referred to as second
Time depth function;
The time depth function of up converted shear wave lineups is set up according to the 3rd time depth relational sequence, the referred to as the 3rd
Time depth function;
The time depth function of descending shear wave lineups, the referred to as the 4th time is set up according to the 4th time depth relational sequence
Depth function;
Set up the time depth function of the up reflection wave lineups of shear wave according to the 5th time depth relational sequence, claim
It it is the 5th time depth function;
Set up the time depth function of the up conversion compressional wave lineups of shear wave according to the 6th time depth relational sequence, claim
It it is the 6th time depth function.
Method the most according to claim 6, is characterized in that, determines pre-depth measurement according to described time depth function
Degree and predicted time specifically include:
The first depth prediction result, the very first time is obtained according to very first time depth function and the second time depth function
Predict the outcome;
The second depth prediction result, the second time is obtained according to very first time depth function and the 3rd time depth function
Predict the outcome;
The 3rd depth prediction result, the 3rd time is obtained according to the second time depth function and the 3rd time depth function
Predict the outcome;
The 4th depth prediction result, the 4th time is obtained according to the 4th time depth function and the 5th time depth function
Predict the outcome;
The 5th depth prediction result, the 5th time is obtained according to the 4th time depth function and the 6th time depth function
Predict the outcome;
The 6th depth prediction result, the 6th time is obtained according to the 5th time depth function and the 6th time depth function
Predict the outcome;
To the first described depth prediction result, the second depth prediction result, the 3rd depth prediction result, the 4th depth prediction
Result, the 5th depth prediction result and the 6th depth prediction result are averaged value computing, obtain predetermined depth;
Predict the outcome the described very first time, the second time prediction result and three-time estimates result are averaged value fortune
Calculate, obtain compressional wave predicted time, pre-to described the 4th time prediction result, the 5th time prediction result and the 6th time
Survey result to be averaged value computing, obtain shear wave predicted time.
8. a depth prediction system before brill based on the many wave datum of vertical seismic, is characterized in that, described system specifically includes:
Seismic data acquisition device, is used for gathering the many wave datum of vertical seismic;
Seismic data process device, for processing, after being processed vertically the described many wave datum of vertical seismic
Shake many wave datum;
Seismic profile acquisition device, is used for obtaining well surface seismic section;
Target zone identification device, special for the position and ripple group identifying target zone from described crossing well surface seismic section
Levy;
Reflection line-ups identification device, for according to vertical after described process of the position of described target zone and wave group feature
Seismic Multi-wave data identify direct wave corresponding to described target zone, converted wave and reflection line-ups;
Time deep function set up device, for set up the time according to described direct wave, converted wave and reflection line-ups deep
Degree function;
Device is deeply determined, for determining predetermined depth and predicted time according to described time depth function during prediction.
System the most according to claim 8, is characterized in that, described seismic data process device specifically includes:
Dephasing processing module, for carrying out dephasing process to the described many wave datum of vertical seismic;
Noise attentuation processing module, the many wave datum of vertical seismic after processing dephasing carry out noise attentuation process;
Energy compensating processing module, the many wave datum of vertical seismic after processing noise attentuation carry out energy compensating process;
Spin orientation processing module, the many wave datum of vertical seismic after processing energy compensating carry out spin orientation process;
Deconvolution processing module, the many wave datum of vertical seismic after processing spin orientation carry out deconvolution process, obtain
The many wave datum of vertical seismic after process.
System the most according to claim 8, is characterized in that, the many wave datum of described vertical seismic include descending compressional wave,
The up conversion compressional wave of up compressional wave, up converted shear wave, descending shear wave, the up reflection wave of shear wave and shear wave;Institute
Direct wave, converted wave and the reflection line-ups stated include descending through compressional wave lineups, up reflected P-wave lineups,
Up converted shear wave lineups, up reflection wave lineups, up conversion compressional wave lineups and descending through shear wave homophase
Axle or downward converted P-SV wave lineups.
11. systems according to claim 10, is characterized in that, described time deep function is set up device and specifically included:
Time deep sequence pickup model, for picking up the crest of described direct wave, converted wave and reflection line-ups, obtain
Time depth relational sequence;
Time deep function set up module, for setting up time depth function according to described time depth relational sequence.
12. systems according to claim 11, is characterized in that, described time deep sequence pickup model specifically includes:
Deep sequence pickup unit when first, for picking up the crest of descending through compressional wave lineups, obtains the very first time degree of depth
Relational sequence;
Deep sequence pickup unit when second, for picking up the crest of up reflected P-wave lineups, obtains the second time depth
Relational sequence;
Deep sequence pickup unit when the 3rd, for picking up the crest of up converted shear wave lineups, obtains the 3rd time depth
Relational sequence;
Deep sequence pickup unit when the 4th, for picking up descending through shear wave lineups or the ripple of downward converted P-SV wave lineups
Peak, obtains the 4th time depth relational sequence;
Deep sequence pickup unit when the 5th, for picking up the crest of the up reflection wave lineups of shear wave, obtained for the 5th time
Depth relationship sequence;
Deep sequence pickup unit when the 6th, for picking up the crest of the up conversion compressional wave lineups of shear wave, obtained for the 6th time
Depth relationship sequence.
13. systems according to claim 12, is characterized in that, described time deep function is set up module and specifically included:
When first, function sets up unit, same for setting up descending through compressional wave according to described very first time depth relationship sequence
The time depth function of phase axle, referred to as very first time depth function;
When second, function sets up unit, for setting up up reflected P-wave lineups according to the second time depth relational sequence
Time depth function, the referred to as second time depth function;
When the 3rd, function sets up unit, for setting up up converted shear wave lineups according to the 3rd time depth relational sequence
Time depth function, the referred to as the 3rd time depth function;
When the 4th, function sets up unit, for setting up the time of descending shear wave lineups according to the 4th time depth relational sequence
Depth function, the referred to as the 4th time depth function;
When the 5th, function sets up unit, same for setting up the up reflection wave of shear wave according to the 5th time depth relational sequence
The time depth function of phase axle, the referred to as the 5th time depth function;
When the 6th, function sets up unit, same for setting up the up conversion compressional wave of shear wave according to the 6th time depth relational sequence
The time depth function of phase axle, the referred to as the 6th time depth function.
14. systems according to claim 13, is characterized in that, deeply determine that device specifically includes during described prediction:
Module is deeply determined, for obtaining first according to very first time depth function and the second time depth function when first
Depth prediction result, the very first time predict the outcome;
Module is deeply determined, for obtaining second according to very first time depth function and the 3rd time depth function when second
Depth prediction result, the second time prediction result;
Module is deeply determined, for obtaining the 3rd according to the second time depth function and the 3rd time depth function when the 3rd
Depth prediction result, three-time estimates result;
Module is deeply determined, for obtaining the 4th according to the 4th time depth function and the 5th time depth function when the 4th
Depth prediction result, the 4th time prediction result;
Module is deeply determined, for obtaining the 5th according to the 4th time depth function and the 6th time depth function when the 5th
Depth prediction result, the 5th time prediction result;
Module is deeply determined, for obtaining the 6th according to the 5th time depth function and the 6th time depth function when the 6th
Depth prediction result, the 6th time prediction result;
Predetermined depth determines module, for the first described depth prediction result, the second depth prediction result, the 3rd degree of depth
Predict the outcome, the 4th depth prediction result, the 5th depth prediction result, the 6th depth prediction result are averaged value computing,
Obtain predetermined depth;
Predicted time determines module, for the described very first time is predicted the outcome, the second time prediction result, the 3rd time
Predict the outcome the value computing that is averaged, and obtains compressional wave predicted time, pre-to the 4th described time prediction result, the 5th time
Survey result, the 6th time prediction result is averaged value computing, obtains shear wave predicted time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410383847.2A CN104199107B (en) | 2014-08-06 | 2014-08-06 | Depth prediction approach and system before brill based on the many wave datum of vertical seismic |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410383847.2A CN104199107B (en) | 2014-08-06 | 2014-08-06 | Depth prediction approach and system before brill based on the many wave datum of vertical seismic |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104199107A CN104199107A (en) | 2014-12-10 |
CN104199107B true CN104199107B (en) | 2016-08-31 |
Family
ID=52084419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410383847.2A Active CN104199107B (en) | 2014-08-06 | 2014-08-06 | Depth prediction approach and system before brill based on the many wave datum of vertical seismic |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104199107B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107765309B (en) * | 2017-10-09 | 2019-03-29 | 中国石油天然气集团有限公司 | Depth prediction approach before a kind of brill based on vertical seismic data |
CN108107485B (en) * | 2017-11-17 | 2019-11-08 | 中国石油天然气集团公司 | The appraisal procedure and device of seismic data multiple wave pollution level |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0871904B1 (en) * | 1995-11-27 | 2003-05-28 | WesternGeco Limited | Method of monitoring quality of seismic data processing and method of processing vertical seismic profile data |
GB2359364B (en) * | 1999-11-11 | 2003-07-02 | Bp Amoco Corp | Vector recomposition of seismic 3D converted-wave data |
CN102798888A (en) * | 2011-05-27 | 2012-11-28 | 中国石油天然气集团公司 | Method for calculating velocity ratio of longitudinal wave to transverse wave by using non-zero wellhead distance data |
CN103809216A (en) * | 2012-11-12 | 2014-05-21 | 中国石油天然气集团公司 | Resistivity data and earthquake data combined speed field construction method |
-
2014
- 2014-08-06 CN CN201410383847.2A patent/CN104199107B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0871904B1 (en) * | 1995-11-27 | 2003-05-28 | WesternGeco Limited | Method of monitoring quality of seismic data processing and method of processing vertical seismic profile data |
GB2359364B (en) * | 1999-11-11 | 2003-07-02 | Bp Amoco Corp | Vector recomposition of seismic 3D converted-wave data |
CN102798888A (en) * | 2011-05-27 | 2012-11-28 | 中国石油天然气集团公司 | Method for calculating velocity ratio of longitudinal wave to transverse wave by using non-zero wellhead distance data |
CN103809216A (en) * | 2012-11-12 | 2014-05-21 | 中国石油天然气集团公司 | Resistivity data and earthquake data combined speed field construction method |
Non-Patent Citations (3)
Title |
---|
利用零井源距VSP资料自动预测钻头前方目的层深度;王成礼;《石油地球物理勘探》;19980630;第23卷(第3期);355-357 * |
广义非线性反演及其应用;胡光岷,等;《石油物探》;20000630;第39卷(第2期);44-51 * |
零偏VSP叠前资料中的低频阴影检测;张固澜,等;《石油地球物理勘探》;20110630;第46卷(第3期);438-441 * |
Also Published As
Publication number | Publication date |
---|---|
CN104199107A (en) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104375188B (en) | Seismic wave transmission attenuation compensation method and device | |
CN104181596B (en) | Geologic horizon automatic tracking method and device | |
CN107765298B (en) | A kind of method and device of determining near-surface velocity model | |
CN108254780A (en) | A kind of microseism positioning and anisotropic velocity structure tomographic imaging method | |
CA2641867A1 (en) | Methods and apparatus for predicting the hydrocarbon production of a well location | |
CN106990438A (en) | Complicated earth surface region micro logging point position and well depth determine method and device | |
CN103675907A (en) | AVO inversion hydrocarbon detection method based on petrographic constraints | |
CN104793247A (en) | Near-surface speed modeling method and device | |
CN104297784A (en) | Primary wave azimuthal anisotropy based fracture predicting method | |
CN107942404B (en) | A kind of method and device in determining crack | |
CN105301639A (en) | Speed field inversion method and device based on VSP double-weight travel time tomography | |
CN106405644B (en) | Crack determines method and apparatus | |
CN111983683B (en) | Prediction method and system for lake-facies limestone reservoir under low-well condition | |
CN108333627A (en) | Igneous rock area is broken recognition methods and the device of the true and false | |
WO2022011104A1 (en) | Generating a model for seismic velocities in a subsurface region using inversion with lateral variations | |
CN105259581A (en) | Seismic data time-depth conversion method | |
CN104199107B (en) | Depth prediction approach and system before brill based on the many wave datum of vertical seismic | |
CN105866839B (en) | Static correction method and static correction device on basis of common-attitude trace gathers | |
US20210208295A1 (en) | Iterative migration velocity optimization for a vsp survey using semblance | |
CN109444960A (en) | Determination method for fault tectonic | |
CN114861515A (en) | Method, device, equipment and medium for calculating layer speed data volume | |
CN109188522B (en) | Velocity field construction method and device | |
CN109991663B (en) | Work area seismic velocity body correction method and device | |
CN103076629B (en) | Oil-gas exploration method and device based on longitudinal and horizontal seismic wave velocity ratio | |
CN106526668B (en) | Original waveform extraction and imaging method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |