CN105988134A - Multi-wave synthesis record calibration method - Google Patents
Multi-wave synthesis record calibration method Download PDFInfo
- Publication number
- CN105988134A CN105988134A CN201510071248.1A CN201510071248A CN105988134A CN 105988134 A CN105988134 A CN 105988134A CN 201510071248 A CN201510071248 A CN 201510071248A CN 105988134 A CN105988134 A CN 105988134A
- Authority
- CN
- China
- Prior art keywords
- wave
- ray
- converted
- converted wave
- reflection coefficient
- 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.)
- Pending
Links
Abstract
The invention provides a multi-wave synthesis record calibration method. The multi-wave synthesis recording calibration method comprises the following steps of step1, carrying out multi-component ray tracking earthquake wave field simulation and analyzing a multiple-wave field response characteristic and a longitudinal wave and transverse wave difference; step2, carrying out transverse wave speed estimation on a well without all-wave line logging information and providing transverse wave curve data for converted wave synthesis record making; step3, carrying out converted wave reflection coefficient calculation, multi-wave earthquake wavelet extraction and convolution; and step4, carrying out multi-wave synthesis record integration calibration. By using the multi-wave synthesis recording calibration method, indication meanings of different types of multi-wave synthesis records to an oil gas are clear, and the method plays an important role in multi-wave multi-component earthquake data integration interpretation and has a certain guidance effect on pertinency of a multi-wave earthquake acquisition block and oil gas detection.
Description
Technical field
The present invention relates to oil field development technical field, especially relate to a kind of many ripples composite seismogram calibration method.
Background technology
Multi-wave seismic survey is a kind of new method of seismic prospecting growing up the 1980s, and it is as three
The appearance of component wave detector and occur.Theoretical research in early days to be laid particular emphasis on and set up basic process framework.With seabed
Cable is succeeded in developing, going deep into of transformed wave seismic prospecting technique study, and converted wave seismic technology has had new progress, promotes
The upsurge of external marine multi-component exploration, particularly in the North Sea and the Gulf of Mexico.And by land, static correction and face ripple are done
Problem of disturbing is more serious than sea much, while it is true, people are still devoted to the research of multicomponent seismic survey.Many points of many ripples
Amount seismic wave field meets space (three-dimensional) elastic wave propagation rule, the seismic prospecting that its all-wave field elastic wave excites and receives
Can reflect that underground medium is elastic and formation lithology information, with multiwave multicomponent earthquake source investigation, these are concealed comprehensively
Exploration targets, is a kind of good method.Wherein, multi-component converted wave seismic technology both had the longitudinal wave exploration degree of depth big,
Data collection is relatively easy and the feature of small investment, can reflect that again the shear wave velocity of underground medium changes.Multi-component converted wave
This feature of earthquake, makes the Direct Recognition of lithology exploration and oil gas be possibly realized.Simultaneously because multicomponent data acquisition,
On the premise of recording two horizontal component geological datas, it is possible to use the Shear Waves Splitting time difference reflection that shear wave splitting produces is split
The principal direction that seam is grown and growth density so that the exploration and development of crack crack type oil-gas reservoir is possibly realized.Nowadays multi-components
Converted wave seismic technology and the anisotropy theory and method research closely coupled with this technology have become domestic and international earthquake and have surveyed
One of the study hotspot in spy field.
In recent years, under oil-gas exploration industry demand promotion, multi-wave seismic Comprehensive interpretation development in recent years is very fast, this
Planting seismic technology uses p-wave source to excite, and receives at ground three-component geophone, obtains compressional wave and converted wave data simultaneously.
The feature utilizing multi-component seismic data informative can improve imaging precision, carry out anisotropy and crack identification,
Extract more physical property characteristic parameter (such as wave-amplitude difference, P-S wave velocity ratio, all kinds of elastic parameters etc. in length and breadth), thus
Detailed predicting reservoir lithology, even Direct Recognition oil-gas reservoir.Compared with longitudinal wave earthquake data, due to multi component signal
Vector characteristic, anisotropy and converted wave propagation path asymmetry, reflect the reflection line-ups at same subsurface geology interface
Time Inconsistency, therefore ripe conventional compressional wave is processed and interpretation technique cannot be applied directly to multiwave multicomponent earthquake data
, as many difficult problems such as converted wave statics, skew, the faults comparison demarcation of many ripples still become multiwave multicomponent earthquake skill
The bottleneck of art development.
It is the important step that in seismic data interpretation, many ripples data is combined with geology that faults comparison is demarcated, and is multi-wave seismic
Most critical, most basic work in data interpretation and application.Making converted wave Prof. Du Yucang record is that faults comparison is demarcated
, there is recurrence relation between longitudinal wave reflection coefficient and wave impedance in key, therefore, relies on convolution model, and its Forward Modeling and Inversion is asked
Topic is all easier to solve.And for converted wave, due to pip not at (being positioned at same level) focus and connecing
It on midpoint between sink, and is typically to receive in the range of cannon offset, thus Converted wave reflection coefficient is unlike compressional wave
Meeting recurrence relation like that, this results in difficulty and strengthens, and currently without very complete converted wave synthetic seismic record mark
Determine method.At present, continuing to increase with difficulties in exploration, the research of new technology and application become inevitable, therefore study simultaneously
The converted wave synthetic seismic record calibration technique forming complete set very must for the explanation application of multi-wave seismic data
Want.We have invented a kind of new many ripples composite seismogram calibration method for this, solve above technical problem.
Content of the invention
It is an object of the invention to provide a kind of by multi-wave seismic wave-field simulation, shear wave curve estimation, converted wave reflection being
Researching and analysing of the number key technique such as calculating, converted wave synthetic seismic record making, defines many ripples synthesis of the whole series
Record makes and comparative analysis method.
The purpose of the present invention can be achieved by the following technical measures: many ripples composite seismogram calibration method, and converted wave reflection is
Number Sequence and the convolution of wavelet, this many ripples composite seismogram calibration method includes: step 1, carries out multi-components ray tracing earthquake
Wave-field simulation, analysis Multiwave field response characteristic and in length and breadth ripple otherness;Step 2, does not has the well of full wave train log data to enter
Row shear wave velocity is estimated, provides shear wave profile for converted wave synthetic record;Step 3, carrying out converted wave reflection is
Number calculates, the extraction of multi-wave seismic wavelet and convolution;Step 4, carries out many ripples composite traces composite calibration.
The purpose of the present invention realizes also by following technical measures:
In step 1, by ray casting to the travel routes calculating ray or a certain geophone station of wave field arrival and travelling
When, thus these are many to single interface, double interface, multiple solutions, wedge shape lithologic body, the double interfaces with tomography, gassiness sand water sand
Individual model carries out compressional wave, the simulation of converted wave wave field characteristics, thus researchs and analyses Multiwave field response characteristic and ripple difference in length and breadth
Property.
In step 1, first mesh generation is carried out to communication space, sets up the network of a two dimension or three-dimensional so that
Each lattice has different velocity amplitude;Then from known focus, ray or ripple are calculated with ray equation
Field reaches travel routes and the whilst on tour of a certain geophone station.
In step 1, during the travel routes calculating ray or a certain geophone station of wave field arrival and whilst on tour, from mould
A certain focus in type net boundary is lighted, a selected initial transmissions angle, makes ray advance along this direction, when running into net
During the border of compartment, calculate the angle of emergence through ray according to Snell law, then make ray advance along new direction, directly
To the border that ray reaches geophone station place, this terminal and constraints geophone station position are made one and compare,
If trial fire exists error, after adjusting the angle of departure according to error condition, repeat the above steps, until the terminal of ray and inspection
Till error between the position of wave point meets the requirements prescribed, once meet above-mentioned requirements, can be according to ray path meter
Calculate whilst on tour and ray length.
In step 2, study area is had to the well of actual measurement shear wave, utilize multiple empirical equation to carry out estimation contrast, preferably go out
Empirical equation the most close, other do not have the well of shear wave data to be applied to study area;Shear wave is not surveyed for study area but
It is to have converted wave seismic data, carried out with actual converted wave profile by calculating the converted wave composite traces of different Forecasting Methodology
Contrast, window comparative analysis when being related to by target interval, preferably go out optimal shear wave Forecasting Methodology and carry out the horizontal stroke of remaining well
Ripple is predicted.
In step 3, Ostrander gas sand mo(u)ld type is set up, near by analyzing the different reflectance factor in AB, BC interface in model
Like the computational accuracy of formula, specify each equation approximate transform wave reflection coefficient pluses and minuses, preferably Aki-richards formula,
Carry out simplified formula, complete the calculating of Converted wave reflection coefficient.
In step 4, the thinking of equivalent layer is used to ask for incident compressional angle, by Converted wave reflection coefficient computing formula, meter
Calculate the Converted wave reflection coefficient sequence road collection of different geophone offset x, then these Converted wave reflection coefficient sequence road collection are added
Summation, obtains the Converted wave reflection coefficient sequence identical with actual superposition road, then with selected seismic wavelet convolution,
Obtain converted wave artificial synthesized E-selectin.
In step 4, use the thinking of equivalent layer to ask for P ripple incidence angle θ, i.e. use average speed Va and geophone offset X to represent
The incidence angle θ of any j-th reflecting interface can be written as
Wherein, Vp is velocity of longitudinal wave, t0When being longitudinal wave propagation outward journey.Here θ angle, corresponding j-th reflecting interface,
I.e. θ angle in Converted wave reflection coefficient approximate formula, utilizes Converted wave reflection coefficient approximate expression, calculates at geophone offset x solid
Converted wave reflection coefficient sequence under conditions of determining.
Many ripples composite seismogram calibration method in the present invention, develops and utilizes ray-trace modeling method to multiwave multicomponent earthquake
Wave field is analyzed research, the compressional wave under clearly different forward models, converted wave wave field response feature;The proposition of novelty
The strategy of shear wave estimation and method in the case of two kinds, by calculate the converted wave composite traces of different Forecasting Methodology with actual
Converted wave sectional plane contrasts, and preferably goes out optimal shear wave Forecasting Methodology;At high-precision shear wave curve estimation technique study
On the basis of, by setting up Ostrander gas sand mo(u)ld type, use Ostrander gas sand mo(u)ld type to numerous reflectance factor calculating sides
Method has carried out modeling computation, the precision and stability of each approximate formula in comparative analysis model, researchs and analyses each approximation side
After journey Converted wave reflection coefficient pluses and minuses, preferably approximate formula, and write corresponding calculation procedure, completing converted wave reflection is
The calculating of number;The thinking taking equivalent layer asks for incidence angle, calculates the Converted wave reflection coefficient sequence road of different geophone offset
Collection, and be added summation to these roads collection, obtains the Converted wave reflection coefficient sequence identical with actual superposition road, then with select
Seismic wavelet convolution, complete the making of converted wave synthetic seismogram.The method converted wave wave field response feature and routine
Compressional wave is different;Reflectance factor is unlike the reflectance factor of conventional compressional wave and there is recurrence relation between wave impedance.Can not only
Enough carry out converted wave composite traces Fine calibration, and reflected on composite traces by research compressional wave and converted wave further
The difference of feature, in conjunction with Multiwave field response characteristic, specify that the indicative significance to oil gas for the dissimilar many ripples composite traces.
Multiwave multicomponent earthquake Information integration is explained have important function by the research of the method, gathers the pin of block to multi-wave seismic
Certain directive function is also had to property and oil and gas detection.
Brief description
Fig. 1 is the flow chart of a specific embodiment of many ripples composite seismogram calibration method of the present invention;
Fig. 2 be the present invention a specific embodiment in the thick wedge shape lithologic body rate pattern figure of multiple solutions;
Fig. 3 be the present invention a specific embodiment in the thick wedge shape lithologic body longitudinal wave reflection figure of multiple solutions;
Fig. 4 be the present invention a specific embodiment in the thick wedge shape lithologic body converted wave reflectogram of multiple solutions;
Fig. 5 be the present invention a specific embodiment in Ostrander gas-bearing sandstone model figure;
Fig. 6 be the present invention a specific embodiment in Ostrander gas sand AB critical transition wave reflection coefficient figure;
Fig. 7 be the present invention a specific embodiment in Ostrander gas sand BC critical transition wave reflection coefficient figure;
Fig. 8 be the present invention a specific embodiment in cultivate 71-inspection 41 well compressional wave composite seismogram calibration figures;
Fig. 9 be the present invention a specific embodiment in cultivate 71-inspection 41 well converted wave composite seismogram calibration figures.
Detailed description of the invention
For making the above and other objects, features and advantages of the present invention to become apparent, cited below particularly go out preferred embodiment,
And coordinate institute's accompanying drawings, it is described in detail below.
As it is shown in figure 1, the flow chart of a specific embodiment of many ripples composite seismogram calibration method that Fig. 1 is the present invention.
In step 101, carry out multi-components ray tracing Simulation of Seismic Wave.By ray casting to calculating ray or wave field
Reach travel routes and the whilst on tour of a certain geophone station, thus to single interface, double interface, multiple solutions, wedge shape lithologic body (as
Shown in Fig. 2), the double interfaces with tomography, multiple model such as gassiness sand water sand carry out compressional wave, the simulation of converted wave wave field characteristics,
Thus research and analyse Multiwave field response characteristic and ripple otherness in length and breadth.
Carry out multi-wave and multi-component Processing and Interpretation Technology research, it is necessary first to its reflected wave field feature, reflectance mechanism is carried out deeply
Enter research.In order to Study of Seismic ripple is at the propagation law of underground, first mesh generation is carried out to communication space, set up one
Two-dimentional or three-dimensional network so that each lattice has different velocity amplitude.Then from known focus,
Calculate ray with ray equation or wave field reaches travel routes and the whilst on tour of a certain geophone station.The core of its computational methods
It is to iterate, Approach by inchmeal.It is exactly specifically to light from the borderline a certain focus of model meshes, selected one
Initial transmissions angle, makes ray advance along this direction, when running into the border between grid, just calculates according to Snell law and passes through
The angle of emergence of ray, then makes ray advance along new direction.So on, until ray reaches geophone station place
Border till.At this moment we need to compare this terminal and constraints geophone station position work one, in general,
Trial fire for the first time always has error.After adjusting the angle of departure according to error condition, repeat the above steps, until ray
Till error between terminal and the position of geophone station meets the requirement of our regulation.Once meet above-mentioned requirements, Ji Kegen
Calculate whilst on tour and ray length according to ray path.
By this analogue technique, the propagation time of gained seismic wave is relatively more accurate, intuitive display, is suitable for calculating the trip of ripple
During row, PP ripple and PS ripple can be simulated respectively, be particularly advantageous in judging large-scale geological structure.It is easy to our clear and definite ripple in length and breadth
Wave field reflectance signature.The such as lithologic body rate pattern of multiple solutions wedge shape shown in Fig. 2, can be divided by above wave-field simulation method
Do not obtain longitudinal wave reflection figure (as shown in Figure 3) and converted wave reflectogram (as shown in Figure 4), can be obvious by this two figures
Finding out that compressional wave and converted wave have notable difference on reflectance signature, converted wave is stronger than longitudinal wave reflection energy, and reflectance signature
By vertical magnification, thus wedge shape lithologic body is easier to identified, and stratum is easier to follow the trail of.Specify that wave field is anti-in length and breadth
After penetrating feature, flow process enters into step 102.
In step 102, carry out shear wave velocity estimation.The accurately acquisition of shear wave logging data is to many ripples composite seismogram calibration very
Crucial.Take Different Strategies and method for different situations, complete the estimation of high-precision shear wave velocity, close for converted wave
Record making is become to provide shear wave profile accurately.
The spread speed size of seismic wave be by the skeleton of rock, hole, fluid and the relevant environment such as temperature, pressure because of
The common comprehensive restriction of element.P-and s-wave velocity has important application at aspects such as reservoir lithology, physical property and fluid identification, horizontal
Accurately obtaining of ripple well-log information is very crucial to many ripples composite seismogram calibration, and the shear wave velocity being obtained by well logging at present is main
It is to extract from full wave train log.But the well typically doing full wave train log is fewer, even doing full wave train log, also
Only measure a certain depth segment, and when we do multi-wave seismic composite seismogram calibration, generally be directed to full well section, therefore at shear wave
The well logging interval of speed disappearance must carry out the estimation of shear wave velocity curve.
Domestic and international at present shear wave Forecasting Methodology has integrated two big classes, and one is empirical equation, and two is based on petrophysics
Method.The strategy proposing the estimation of shear wave in the case of two kinds of novelty and method in research process: one is for grinding
Studying carefully district has the well of actual measurement shear wave, utilizes multiple empirical equation to carry out estimation contrast, preferably goes out empirical equation the most close, should
For study area, other do not have the well of shear wave data;Two is study area not to be surveyed to shear wave still to have converted wave earthquake to provide
Material, is contrasted with actual converted wave profile by calculating the converted wave composite traces of different Forecasting Methodology, passes through destination layer
Window comparative analysis during Duan Xiangguan, preferably goes out optimal shear wave Forecasting Methodology and carries out the shear wave prediction of remaining well.Flow process enters
To step 103.
In step 103, carry out Converted wave reflection coefficient calculating.Set up Ostrander gas sand mo(u)ld type, as it is shown in figure 5, design
Model for folder gas sand in the middle of shale up and down, relative to the shale of A and C layer, the velocity of longitudinal wave of B layer gas sand subtracts
Having lacked 20%, density about decreases 10%.By analyzing the calculating of the different reflectance factor approximate formula in AB, BC interface in model
Precision, specifies each equation approximate transform wave reflection coefficient pluses and minuses, final preferred Aki-richards formula, has carried out public affairs
Formula simplifies and has write corresponding calculation procedure, completes the calculating of Converted wave reflection coefficient.
Due on converted wave the pip not midpoint between focus and receiving point, so the reflectance factor of converted wave is unlike indulging
Like that and there is recurrence relation between wave impedance in the reflectance factor of ripple.It is Zoeppritz to Converted wave reflection coefficient estimation basis
Equation, that reflects plane wave incidence to solid dielectric interface produced reflected P-wave (PP ripple), reflection wave (PS
Ripple), transmitted P-wave, transmitted shear wave, and give bit shift equation group.But, the analytic solutions form of equation group is very
Complexity, needs to go to describe the relation between them with a matrix, is difficult to the impact on amplitude for the Direct Analysis medium parameter,
Expressed physical significance is not directly perceived yet.Numerous scholars derive a series of approximate equation based on it, describe from all angles
The reflectance factor of PS ripple.The present invention establishes Ostrander gas sand mo(u)ld type and carries out model examination to multiple reflectance factors calculating method
Calculate, analyze different approximate equation respectively in the precision at AB interface and BC interface reflection coefficients.Fig. 6 is gas sand AB interface
Converted wave reflection coefficient figure, three kind approximate formulas have relatively less than 70 ° of result of calculations in incidence angle as we can see from the figure
Good uniformity, the strongest at 40 ° of reflected energies, cross 70 ° of Zoeppritz and calculate reflectance factor nearly zero with Aki formula,
And BP method reflectance factor polarity becomes negative, BP method is suitable for the incident back wave approximate calculation of high angle.Fig. 7 is gas sand
BC critical transition wave reflection coefficient figure, as we can see from the figure in the ranges of incidence angles less than 35 °, each equation calculates knot
Fruit has preferable uniformity, and more than about 45 °, Aki has bigger error with BP formula and Zoeppritz formula, exceedes
60 ° of reflectance factor polarity invert.Comprehensively think in the ranges of incidence angles less than 40 °, the AB that each equation calculates and
BC interface P-SV wave reflection coefficient, all has higher precision.
By the analysis to three kinds of approximate formula precision, research has been carried out approximate expression preferably, simplify, and write mutually accrued
Calculation program, completes independent Converted wave reflection coefficient computational methods.Flow process enters into step 104.
In step 104, carry out many ripples composite traces composite calibration.The thinking using equivalent layer asks for incident compressional angle, passes through
The Converted wave reflection coefficient computing formula that upper step research obtains, calculates the Converted wave reflection coefficient sequence road of different geophone offset x
Then these Converted wave reflection coefficient sequence road collection are added summation, obtain the converted wave identical with actual superposition road anti-by collection
Penetrate coefficient sequence.Again with selected seismic wavelet convolution, i.e. can get converted wave artificial synthesized E-selectin.Such as Fig. 8, figure
9 respectively cultivate 71-inspection 41 well compressional waves and converted wave composite seismogram calibration figure, can be seen that converted wave synthesis note from two figure contrasts
Record polarity inverts, and wavelet dominant frequency reduces, and the time has substantially stretching.
Make the key that converted wave Prof. Du Yucang record is converted wave and compressional wave joint interpretation, for converted wave, due to
Pip is not on (being positioned at same level) midpoint between focus and receiving point, and is typically at cannon offset model
Enclosing interior reception, thus Converted wave reflection coefficient meeting recurrence relation unlike compressional wave, this results in converted wave Forward Modeling and Inversion
Difficulty strengthen.Therefore calculate converted wave Prof. Du Yucang record when incident compressional angle to be solved ask for problem.For altogether
For every one that CCP road is concentrated, the converted wave receiving is under conditions of geophone offset is fixed, from each interface echo
Superposition, this has certain difference with along the Converted wave reflection coefficient of directions of rays.
The thinking using equivalent layer asks for P ripple incidence angle θ, i.e. uses root mean sequare velocity VrmsAny jth representing with geophone offset X
The incidence angle θ of individual reflecting interface can be written as
But between well head to log data starting point, there is no interval transit time data, carrying out synthetic seismogram calculating
In, this part data all use an average speed to replace, and so when calculating θ by formula (1), well head is to log data
Speed between Duan Qidian can only use average speed, and remainder uses root mean sequare velocity, but effect is not necessarily preferable.
Therefore, reality can use root mean sequare velocity V in average speed Va alternate form (2)rmsCalculate incident compressional angle
Here θ angle, the θ angle in corresponding j-th reflecting interface, i.e. Converted wave reflection coefficient approximate formula.So, profit
The Converted wave reflection coefficient approximate expression obtaining with research institute before, it is possible to calculate conversion under conditions of geophone offset x fixes
Wave reflection coefficient sequence.For a GongCCPDian road collection, the Converted wave reflection coefficient sequence of different geophone offset x can be calculated
Then these Converted wave reflection coefficient sequence road collection are added summation, obtain the converted wave identical with actual superposition road by road collection
Reflection coefficient sequence.Again with selected seismic wavelet convolution, i.e. can get converted wave artificial synthesized E-selectin.Flow process is tied
Bundle.
Many ripples composite seismogram calibration method in the present invention, by multi-wave seismic wave-field simulation, shear wave curve estimation, turn
Change researching and analysing of the key techniques such as the calculating of wave reflection coefficient, converted wave synthetic seismic record making, ultimately form
A whole set of many ripples synthetic record and comparative analysis method, establish good for the application of Multi-component seismic exploration integrated interpretation
Good basis, provides power-assisted for the development of multiwave multicomponent earthquake technology.
Claims (8)
1. more than ripple composite seismogram calibration method, it is characterised in that the convolution of Converted wave reflection coefficient sequence and wavelet, this many ripples close
Record scaling method is become to include:
Step 1, carries out multi-components ray tracing Simulation of Seismic Wave, analysis Multiwave field response characteristic and in length and breadth ripple otherness;
Step 2, does not has the well of full wave train log data to carry out shear wave velocity estimation, provides shear wave bent for converted wave synthetic record
Line data;
Step 3, carries out Converted wave reflection coefficient calculating, the extraction of multi-wave seismic wavelet and convolution;
Step 4, carries out many ripples composite traces composite calibration.
2. many ripples composite seismogram calibration method according to claim 1, it is characterised in that in step 1, chased after by ray
Track method reaches travel routes and the whilst on tour of a certain geophone station to calculating ray or wave field, thus to single interface, double interface,
Multiple solutions, wedge shape lithologic body, the double interfaces with tomography, gassiness this multiple models of sand water sand carry out compressional wave, converted wave wave field
Characteristic simulation, thus research and analyse Multiwave field response characteristic and ripple otherness in length and breadth.
3. many ripples composite seismogram calibration method according to claim 1, it is characterised in that in step 1, first to propagation
Space carries out mesh generation, sets up the network of a two dimension or three-dimensional so that each lattice has different
Velocity amplitude;Then from known focus, calculate ray with ray equation or wave field reaches the travelling road of a certain geophone station
Footpath and whilst on tour.
4. many ripples composite seismogram calibration method according to claim 3, it is characterised in that in step 1, is calculating ray
Or during wave field reaches the travel routes of a certain geophone station and whilst on tour, from the borderline a certain focal point of model meshes
Rise, a selected initial transmissions angle, make ray advance along this direction, when running into the border between grid, fixed according to Snell
Rule calculates the angle of emergence through ray, then makes ray advance along new direction, until ray reaches the limit at geophone station place
Till boundary, this terminal and constraints geophone station position are made one and compares, if trial fire exists error, according to error
After situation adjusts the angle of departure, repeat the above steps, until the error between the terminal of ray and the position of geophone station meets
Till the requirement of regulation, once meet above-mentioned requirements, whilst on tour and ray length can be calculated according to ray path.
5. many ripples composite seismogram calibration method according to claim 1, it is characterised in that in step 2, for study area
There is the well of actual measurement shear wave, utilize multiple empirical equation to carry out estimation contrast, preferably go out empirical equation the most close, be applied to
Study area other there is no the well of shear wave data;Study area is not surveyed to shear wave but have converted wave seismic data, pass through
The converted wave composite traces calculating different Forecasting Methodology is contrasted with actual converted wave profile, when being related to by target interval
Window comparative analysis, preferably goes out optimal shear wave Forecasting Methodology and carries out the shear wave prediction of remaining well.
6. many ripples composite seismogram calibration method according to claim 1, it is characterised in that in step 3, set up Ostrander
Gas sand mo(u)ld type, by analyzing the computational accuracy of the different reflectance factor approximate formula in AB, BC interface in model, specifies each equation
Approximate transform wave reflection coefficient pluses and minuses, preferably Aki-richards formula, carry out simplified formula, complete converted wave reflection
The calculating of coefficient.
7. many ripples composite seismogram calibration method according to claim 1, it is characterised in that in step 4, uses equivalent layer
Thinking ask for incident compressional angle, by Converted wave reflection coefficient computing formula, the converted wave calculating different geophone offset x is anti-
Penetrate coefficient sequence road collection, then summation is added to these Converted wave reflection coefficient sequence road collection, obtain superposing phase with reality
With Converted wave reflection coefficient sequence, then with selected seismic wavelet convolution, i.e. can get converted wave Prof. Du Yucang earthquake note
Record.
8. many ripples composite seismogram calibration method according to claim 7, it is characterised in that in step 4, uses equivalent layer
Thinking ask for P ripple incidence angle θ, i.e. use the incidence of any j-th reflecting interface that average speed Va and geophone offset x represents
Angle θ can be written as
Wherein, Vp is velocity of longitudinal wave, t0When being longitudinal wave propagation outward journey, θ angle here, corresponding j-th reflecting interface,
I.e. θ angle in Converted wave reflection coefficient approximate formula, utilizes Converted wave reflection coefficient approximate expression, calculates at geophone offset x solid
Converted wave reflection coefficient sequence under conditions of determining.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510071248.1A CN105988134A (en) | 2015-02-11 | 2015-02-11 | Multi-wave synthesis record calibration method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510071248.1A CN105988134A (en) | 2015-02-11 | 2015-02-11 | Multi-wave synthesis record calibration method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105988134A true CN105988134A (en) | 2016-10-05 |
Family
ID=57041645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510071248.1A Pending CN105988134A (en) | 2015-02-11 | 2015-02-11 | Multi-wave synthesis record calibration method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105988134A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109143335A (en) * | 2018-07-27 | 2019-01-04 | 中国地质调查局沈阳地质调查中心 | A kind of production method of synthetic seismogram, system, medium and equipment |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101013161A (en) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | Seismic exploration position calibration method based on prestack wave field simulation |
CN101329407A (en) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | Method for quick switching wave direct simulation to determine formation lithology and lithofacies change |
EP2261691A2 (en) * | 2002-08-30 | 2010-12-15 | Robinson, John M. | Removal of noise from seismic data using improved radon transformations |
CN104237946A (en) * | 2014-09-19 | 2014-12-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Single-layer reflected P-wave and reflection converted shear wave amplitude matching method based on well control |
-
2015
- 2015-02-11 CN CN201510071248.1A patent/CN105988134A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2261691A2 (en) * | 2002-08-30 | 2010-12-15 | Robinson, John M. | Removal of noise from seismic data using improved radon transformations |
CN101013161A (en) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | Seismic exploration position calibration method based on prestack wave field simulation |
CN101329407A (en) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | Method for quick switching wave direct simulation to determine formation lithology and lithofacies change |
CN104237946A (en) * | 2014-09-19 | 2014-12-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Single-layer reflected P-wave and reflection converted shear wave amplitude matching method based on well control |
Non-Patent Citations (2)
Title |
---|
石建新: "多分量地震处理解释关键方法应用研究", 《中国博士学位论文全文数据库•基础科学辑》 * |
马劲风 等: "转换波人工合成地震记录与速度反演方法研究", 《石油地球物理勘探》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109143335A (en) * | 2018-07-27 | 2019-01-04 | 中国地质调查局沈阳地质调查中心 | A kind of production method of synthetic seismogram, system, medium and equipment |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Veeken | Seismic stratigraphy and depositional facies models | |
US6269310B1 (en) | System for eliminating headwaves in a tomographic process | |
US6826484B2 (en) | 3D prestack time migration method | |
CN102692645B (en) | Method for performing joint inversion on P-wave and S-wave velocity ratio of reservoir by utilizing P-wave and converted wave data | |
CN108415075B (en) | A kind of method for predicting reservoir under condition of less wells | |
CN101329407B (en) | Method for quick switching wave direct simulation to determine formation lithology and lithofacies change | |
CN112883564B (en) | Water body temperature prediction method and prediction system based on random forest | |
CN103645503B (en) | A kind of three-dimensional time territory illumination analysis and vibration amplitude compensation method | |
US20130223187A1 (en) | Geological Structure Contour Modeling and Imaging | |
CN101329405A (en) | Simple method of multi-parameter seismic inversion | |
CN102508293A (en) | Pre-stack inversion thin layer oil/gas-bearing possibility identifying method | |
CN103513277B (en) | A kind of earthquake stratum fracture crack density retrieval method and system | |
CN109557582B (en) | A kind of two dimension multi-component seismic data offset imaging method and system | |
CN104570079A (en) | Time matching method of longitudinal wave and converted shear wave seismic data | |
CN106199704B (en) | A kind of Three-dimendimal fusion submarine cable seismic data velocity modeling method | |
CN104237937A (en) | Pre-stack seismic inversion method and system thereof | |
CN110007340A (en) | Salt dome speed density estimation method based on the direct envelope inverting of angle domain | |
CN104570085B (en) | A kind of wave ray parameter field joint inversion method in length and breadth | |
Al-Khazraji et al. | Uncertainty analysis to assess depth conversion accuracy: A case study of Subba Oilfield, Southern Iraq | |
CN104977609A (en) | Prestack longitudinal wave and transverse wave combined inversion method based on rapid simulated annealing | |
Luo et al. | Registration-free multicomponent joint AVA inversion using optimal transport | |
CN103472482A (en) | Multi-wave seismic data time domain matching method and system based on gene ordering system | |
Velásquez et al. | Depth-conversion techniques and challenges in complex sub-Andean provinces | |
CN105988134A (en) | Multi-wave synthesis record calibration method | |
Heilmann | CRS-stack-based seismic reflection imaging for land data in time and depth domains |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161005 |
|
WD01 | Invention patent application deemed withdrawn after publication |