CN104932021A  Constrained tomography speed modeling method based on reverse ray tracing  Google Patents
Constrained tomography speed modeling method based on reverse ray tracing Download PDFInfo
 Publication number
 CN104932021A CN104932021A CN201410100436.8A CN201410100436A CN104932021A CN 104932021 A CN104932021 A CN 104932021A CN 201410100436 A CN201410100436 A CN 201410100436A CN 104932021 A CN104932021 A CN 104932021A
 Authority
 CN
 China
 Prior art keywords
 ray tracing
 velocity
 angle
 modeling
 common image
 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
 238000003325 tomography Methods 0.000 title abstract description 22
 238000000034 method Methods 0.000 claims abstract description 19
 230000005012 migration Effects 0.000 claims abstract description 19
 239000011159 matrix material Substances 0.000 claims abstract description 16
 230000035945 sensitivity Effects 0.000 claims abstract description 16
 238000004458 analytical method Methods 0.000 claims abstract description 12
 238000004587 chromatography analysis Methods 0.000 claims description 27
 238000003384 imaging method Methods 0.000 claims description 16
 239000000284 extract Substances 0.000 claims description 9
 230000000875 corresponding Effects 0.000 claims description 6
 230000003459 antidromic Effects 0.000 claims description 5
 238000005553 drilling Methods 0.000 claims description 5
 238000005286 illumination Methods 0.000 claims description 5
 238000000605 extraction Methods 0.000 claims description 3
 238000006243 chemical reaction Methods 0.000 abstract description 4
 238000010586 diagram Methods 0.000 description 18
 238000005516 engineering process Methods 0.000 description 3
 238000001514 detection method Methods 0.000 description 2
 230000005284 excitation Effects 0.000 description 2
 230000004048 modification Effects 0.000 description 2
 238000006011 modification reaction Methods 0.000 description 2
 230000000644 propagated Effects 0.000 description 2
 229910052704 radon Inorganic materials 0.000 description 2
 SYUHGPGVQRZVTBUHFFFAOYSAN radon(0) Chemical group data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nMzAwcHgnIGhlaWdodD0nMzAwcHgnIHZpZXdCb3g9JzAgMCAzMDAgMzAwJz4KPCEtLSBFTkQgT0YgSEVBREVSIC0tPgo8cmVjdCBzdHlsZT0nb3BhY2l0eToxLjA7ZmlsbDojRkZGRkZGO3N0cm9rZTpub25lJyB3aWR0aD0nMzAwJyBoZWlnaHQ9JzMwMCcgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMTM4JyB5PScxNzAnIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6NDBweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID5SPC90ZXh0Pgo8dGV4dCB4PScxNjUuNicgeT0nMTcwJyBjbGFzcz0nYXRvbS0wJyBzdHlsZT0nZm9udC1zaXplOjQwcHg7Zm9udC1zdHlsZTpub3JtYWw7Zm9udC13ZWlnaHQ6bm9ybWFsO2ZpbGwtb3BhY2l0eToxO3N0cm9rZTpub25lO2ZvbnQtZmFtaWx5OnNhbnMtc2VyaWY7dGV4dC1hbmNob3I6c3RhcnQ7ZmlsbDojM0I0MTQzJyA+bjwvdGV4dD4KPC9zdmc+Cg== data:image/svg+xml;base64,PD94bWwgdmVyc2lvbj0nMS4wJyBlbmNvZGluZz0naXNvLTg4NTktMSc/Pgo8c3ZnIHZlcnNpb249JzEuMScgYmFzZVByb2ZpbGU9J2Z1bGwnCiAgICAgICAgICAgICAgeG1sbnM9J2h0dHA6Ly93d3cudzMub3JnLzIwMDAvc3ZnJwogICAgICAgICAgICAgICAgICAgICAgeG1sbnM6cmRraXQ9J2h0dHA6Ly93d3cucmRraXQub3JnL3htbCcKICAgICAgICAgICAgICAgICAgICAgIHhtbG5zOnhsaW5rPSdodHRwOi8vd3d3LnczLm9yZy8xOTk5L3hsaW5rJwogICAgICAgICAgICAgICAgICB4bWw6c3BhY2U9J3ByZXNlcnZlJwp3aWR0aD0nODVweCcgaGVpZ2h0PSc4NXB4JyB2aWV3Qm94PScwIDAgODUgODUnPgo8IS0tIEVORCBPRiBIRUFERVIgLS0+CjxyZWN0IHN0eWxlPSdvcGFjaXR5OjEuMDtmaWxsOiNGRkZGRkY7c3Ryb2tlOm5vbmUnIHdpZHRoPSc4NScgaGVpZ2h0PSc4NScgeD0nMCcgeT0nMCc+IDwvcmVjdD4KPHRleHQgeD0nMzUuMDQ1NScgeT0nNTMuNTkwOScgY2xhc3M9J2F0b20tMCcgc3R5bGU9J2ZvbnQtc2l6ZToyM3B4O2ZvbnQtc3R5bGU6bm9ybWFsO2ZvbnQtd2VpZ2h0Om5vcm1hbDtmaWxsLW9wYWNpdHk6MTtzdHJva2U6bm9uZTtmb250LWZhbWlseTpzYW5zLXNlcmlmO3RleHQtYW5jaG9yOnN0YXJ0O2ZpbGw6IzNCNDE0MycgPlI8L3RleHQ+Cjx0ZXh0IHg9JzUxLjA0MDknIHk9JzUzLjU5MDknIGNsYXNzPSdhdG9tLTAnIHN0eWxlPSdmb250LXNpemU6MjNweDtmb250LXN0eWxlOm5vcm1hbDtmb250LXdlaWdodDpub3JtYWw7ZmlsbC1vcGFjaXR5OjE7c3Ryb2tlOm5vbmU7Zm9udC1mYW1pbHk6c2Fucy1zZXJpZjt0ZXh0LWFuY2hvcjpzdGFydDtmaWxsOiMzQjQxNDMnID5uPC90ZXh0Pgo8L3N2Zz4K [Rn] SYUHGPGVQRZVTBUHFFFAOYSAN 0.000 description 2
 210000003229 CMP Anatomy 0.000 description 1
 230000002159 abnormal effect Effects 0.000 description 1
 239000003990 capacitor Substances 0.000 description 1
 238000002591 computed tomography Methods 0.000 description 1
 238000010276 construction Methods 0.000 description 1
 230000001419 dependent Effects 0.000 description 1
 238000005259 measurement Methods 0.000 description 1
 238000011084 recovery Methods 0.000 description 1
 238000004088 simulation Methods 0.000 description 1
Abstract
The invention provides a constrained tomography speed modeling method based on reverse ray tracing. The constrained tomography speed modeling method based on reverse ray tracing comprises the steps of: extracting angle domain common image gathers; automatically fitting the residual curvature of the extracted angle gathers, and calculating a travel time residual by means of a conversion relationship between the travel time residual and a depth residual; by means of a reverse ray tracing method, recording a ray arc length, and calculating a sensitivity matrix; calculating a slow degree disturbance and updating a speed model; and stacking the front depth migration and the extracted angle domain common image gathers, iterating the steps until all controlled layers are updated. The constrained tomography speed modeling method based on reverse ray tracing realizes high migration speed analysis efficiency and high precision in waveform inversion. On condition that a waveform inversion technique cannot be applied in practice, tomography speed inversion has certain advantages and theoretical advantages as a speed updating method.
Description
Technical field
The present invention relates to field of geophysical exploration, particularly relate to a kind of constraint chromatography velocity modeling method based on Inverse ray tracing.
Background technology
Tomography is applied to medical science the earliest, and early 1980s is introduced into applied geophysics field, and how main research utilizes seismic event, electromagnetic wave or other data to earth internal imaging, has entered practical stage since the nineties in 20th century.
The theoretical foundation of tomography is Radon conversion, i.e. u (ρ, θ)=∫
_{l}f (x, y) ds.From the angle of interior of articles image reconstruction, the image of an object slice is two space variable x, and the function of y, is called image function, is designated as f (x, y); With incident wave " irradiation " object of different directions, the wave field information observed is at least the function of incident wave direction θ and observation station position ρ Two Variables, is called projection function, is designated as u (ρ, θ).Within 1917, Radon proves, the projection function u (ρ, θ) of known all incidence angle θs, can Recovery image function f (x, y) uniquely.From the meaning of chromatography, the character such as such as slowness exception, decay of model can add up by the seismic signal (at least can be approximate like this in high frequency limit) propagated along reflection paths, when multiple tracks reflection paths is propagated through this model from many directions, the information being enough to reconstruct this model just can be provided.
For geophysics medium velocity chromatography conversion method, projection function is generally (or disturbance quantity when walking) when walking, slowness (or slowness the is abnormal) distribution of picture function normally underground medium.Therefore, Tomography Velocity inverting mainly solves a system of linear equations: L Δ s=Δ t
Wherein, L is the matrix of m × n dimension, element l
_{ij}the physical meaning of representative is the ray length of ith ray in a jth grid, providing, being commonly referred to sensitivity kernel function (or sensitivity matrix) by just drilling; Δ s is a n × 1 dimension group, and the physical meaning of representative is the slowness value (slowness changing value) in each grid; Δ t is a m × 1 dimension group, and the physical meaning of representative is the whilst on tour of ray (traveltime difference); What dimension m represented is the sum of ray, and n is the lattice number of rate pattern.
Speed is the important parameter of migration imaging and geologic interpretation.In recent years due to the widespread use of prestack depth domain imaging, the Method and Technology setting up rate pattern from geological data is developed rapidly and is applied.Different from the acquiring method of sensitivity kernel function according to the wave field characteristics utilized, chromatography conversion method can be respectively the tomographic inversion modeling based on ray, the limited frequency chromatography based on wave theory.In ray tomography, whilst on tour is obtained by focusdetection field test, then follows the trail of with the reflection of given background velocity field the whilst on tour obtained and compares.According to Fermat principle, the relation between whilst on tour disturbance and velocity disturbance is linear, carries out the velocity field that back projection just obtains renewal along reflection paths to whilst on tour disturbance.Wave theory chromatography postpones to realize by scattered wave field being replaced whilst on tour.Wave field is also obtained by focusdetection field test, and the wave field obtained with given background velocity field forward simulation contrasts, then by the velocity field wave field disturbance back projection of measurement being obtained in the medium of propagation path renewal.From this view point, both implementation procedures are unified.Difference is, wave equation be utilize the approximate or Rytov of Born approximate carry out linearizing.Ray tomography method implementation procedure is relatively simple, and precision is relatively high, uses comparatively extensive in application.
Because ray tomography principle is simple, the data volume of utilization is little (only using whilst on tour information), and the nonlinear degree of inverse problem is relatively weak, forward modelling and inverting solve relatively stable, and counting yield is high, more weak to the dependence of initial model, therefore, ray tomography is most widely used.But, ray tomography is by the restriction of its theory, can only through the finite part of object, high to the susceptibility of limited perspective, usual requirement change of velocity of wave in the scope of a wavelength is very little, the yardstick of heterogeneous body with seismic event appearance than much larger, and requiring that dielectric model is level and smooth, these all limit its range of application.For the relative merits of ray tomography, we have invented a kind of constraint chromatography velocity modeling method based on Inverse ray tracing newly, solve above technical matters.
Summary of the invention
The object of this invention is to provide a kind of take into account migration velocity analysis efficiency and the constraint chromatography velocity modeling method based on Inverse ray tracing of precision of waveform inversion.
Object of the present invention realizes by following technical measures: based on the constraint chromatography velocity modeling method of Inverse ray tracing, should comprise: step 1, and extracted angle domain common image gathers based on the constraint chromatography velocity modeling method of Inverse ray tracing; Step 2, automatic Fitting pickup angle gathers residue curvature, utilizes the transformational relation of whilst on tour residual sum degree of depth residual error to ask for whilst on tour residual error; Step 3, utilizes Inverse ray tracing method, and recording ray arc length, asks for sensitivity matrix; Step 4, asks for slowness disturbance by LSQR method, renewal speed model; And step 5, prestack depth migration and extraction angle domain common image gathers, iteration carries out abovementioned steps, until all keycourses complete renewal.
Object of the present invention also realizes by following technical measures:
In step 1, the rate pattern obtained by general velocity analysis, as initial velocity model, uses the skew of relief surface Gaussian beam carry out migration imaging and extract angle domain common image gathers.
In step 2, on the basis of conventional automatic pick method, utilize the relational expression of ADCIGs road collection residue curvature and angle, residue curvature corresponding to each angle is picked up in matching.
In step 3, migrated section picks up interface, layer position, set up ray tracing and just drilling required reference velocity model, then utilize antidromic illumination raytracing scheme to ask for sensitivity matrix.
In step 3, every bar directions of rays of ray tracing is corresponding with the shooting angle of angle domain common image gathers, carries out Inverse ray tracing using the imaging point at interface as the starting point of ray, makes it to match with observation data by adjustment shooting angle.
In step 4, by whilst on tour residual error, sensitivity matrix, set up the inversion equation group adding regularization factors; Utilize LSQR method solving equation group and renewal speed model.
In step 4, utilize Boundary Using Sonic Logging control refutation process and retrain, amendment inversion equation, with log data sound wave curve constraint velocity inversion result, refutation process iteration is carried out, and finally obtains the tomographic inversion method for solving under borehole restraint.
In steps of 5, the rate pattern that renewal obtains is for prestack depth migration and extract angle domain common image gathers, whether even up according to angle domain common image gathers and accuracy requirement to speed, judge whether also to need iteration again, if also do not evened up, then repeat step 1 to step 5, until meet accuracy requirement, otherwise carry out the velocity inversion of lower one deck, successively analyze from top to bottom, until all layer analysis are complete.
The constraint chromatography velocity modeling method based on Inverse ray tracing in the present invention, belong to the Tomography Velocity Inverse modeling method in seismic data process, the method is a kind of velocity inversion method between wave equation migration velocity analysis and waveform inversion, it has taken into account the efficiency of migration velocity analysis and the precision of waveform inversion, when waveform inversion technology also cannot practical, it is advanced that Tomography Velocity inverting has certain advantage and theory as speed update method.Under study for action, the present invention adopts typical exploratory area model, will carry out precision analysis and compare, with the accuracy of verification method and applicability Tomography Velocity inversion method with waveform inversion method.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of a specific embodiment of the constraint chromatography velocity modeling method based on Inverse ray tracing of the present invention;
Fig. 2 is the schematic diagram of sandconglomerate body model initial velocity field in one embodiment of the invention;
Fig. 3 is the schematic diagram of sandconglomerate body model initial offset result in one embodiment of the invention;
Fig. 4 is the schematic diagram of initial angle territory common imaging gather in one embodiment of the invention;
Fig. 5 is antidromic illumination raypath figure in one embodiment of the invention;
Fig. 6 is the schematic diagram of the velocity field in one embodiment of the invention after the modeling of sandconglomerate body model modification;
Fig. 7 is the schematic diagram of migration result after one embodiment of the invention medium velocity upgrades;
Fig. 8 is the schematic diagram of the angle domain common image gathers after one embodiment of the invention medium velocity upgrades;
Fig. 9 is the schematic diagram of the correct velocity field of sandconglomerate body model in one embodiment of the invention;
Figure 10 is rate curve comparison diagram before and after one embodiment of the invention medium velocity upgrades.
Embodiment
For making above and other object of the present invention, feature and advantage can become apparent, cited below particularly go out preferred embodiment, and coordinate institute's accompanying drawings, be described in detail below.
As shown in Figure 1, Fig. 1 is the process flow diagram of the constraint chromatography velocity modeling method based on Inverse ray tracing of the present invention.
In step 101, extract angle domain common image gathers.The rate pattern obtained by general velocity analysis, as initial velocity model (as Fig. 2), uses the skew of relief surface Gaussian beam carry out migration imaging (as Fig. 3) and extract angle domain common image gathers (as Fig. 4).Fig. 2 is the schematic diagram of sandconglomerate body model initial velocity field in one embodiment of the invention; Obtain the stacking velocity field of time domain based on CMP road collection, then utilize DIX formula to be converted to the interval velocity field of Depth Domain, set up initial velocity field.Fig. 3 is the schematic diagram of sandconglomerate body model initial offset result in one embodiment of the invention; Because initial velocity is inaccurate, the construction location of relative complex is inaccurate, imaging precision is poor, there is the phenomenons such as obvious diffraction interference.Fig. 4 is the schematic diagram of initial angle territory common imaging gather in one embodiment of the invention; When speed is inaccurate, the rally of imaging road shows certain curvature, upwards or curved downwards.Angle gathers presents and significantly upwarps trend as shown in the figure, illustrates that initial velocity is less than normal.Flow process enters into step 102.
In step 102, automatic Fitting pickup angle gathers residue curvature, utilizes the transformational relation of whilst on tour residual sum degree of depth residual error to ask for whilst on tour residual error.The angle domain common image gathers pickup residue depth difference extracted, asks for whilst on tour residual error in the relation of degree of depth surplus and whilst on tour residual error.
The present invention is on the basis of conventional automatic pick method, and utilize the relational expression of ADCIGs road collection residue curvature and angle, residue curvature corresponding to each angle is picked up in matching.Different from the automatic pick method of routine, the method by the impact of capacitor gap, can not obtain the residue curvature of accurate ADCIGs road collection.
The calculated amount calculating whilst on tour residual error due to traditional chromatography velocity inversion is large, efficiency and precision comparison low, and the whilst on tour residual error in the residue curvature information of angle domain common image gathers and Tomography Velocity inversion equation group has and directly contacts, the present invention is by the whilst on tour residual error in the residue curvature information calculating computed tomography velocity inversion system of equations of angle domain common image gathers.Flow process enters into step 103.
In step 103, utilize Inverse ray tracing method, recording ray arc length, asks for sensitivity matrix.Migrated section picks up interface, layer position, sets up ray tracing and just drilling required reference velocity model, then utilize antidromic illumination raytracing scheme to ask for sensitivity matrix.Based on the Tomography Velocity inverting of angle domain common image gathers, the acquisition of its traveltime difference is no longer dependent on ray tracing and is just drilling and provide whilst on tour information, but only need obtain sensitivity matrix.Every bar directions of rays of ray tracing must be corresponding with the shooting angle of angle domain common image gathers, carries out Inverse ray tracing using the imaging point at interface as the starting point of ray, makes it to match with observation data by adjustment shooting angle.Fig. 5 is antidromic illumination raypath figure in one embodiment of the invention, threw model arrive ground acceptance point from imaging point excitation ray.All imaging point excitation ray records obtain ray arc length and are recorded in sensitivity matrix.Flow process enters into step 104.
In step 104, ask for slowness disturbance by LSQR method, renewal speed model.By whilst on tour residual error, sensitivity matrix, set up the inversion equation group adding regularization factors; Utilize LSQR method solving equation group and renewal speed model.Consider the Very Illconditioned of chromatography equation, the present invention adopts LSQR method to solve chromatography system of equations.In order to improve the stability of tomographic inversion, add the regularization of first order derivative type.Meanwhile, Boundary Using Sonic Logging is utilized to control refutation process and retrain, amendment inversion equation, with log data sound wave curve constraint velocity inversion result.Refutation process iteration is carried out, and finally obtains the tomographic inversion method for solving under borehole restraint.Flow process enters into step 105.
In step 105, prestack depth migration and extraction angle domain common image gathers; Iteration carries out said process, until all keycourses complete renewal.Fig. 6 is the schematic diagram of the velocity field in one embodiment of the invention after the modeling of sandconglomerate body model modification; As can be seen from the figure the layer position after upgrading is more clear, and it is more accurate that structure is portrayed.The rate pattern that renewal obtains is for prestack depth migration (as Fig. 7) and extract angle domain common image gathers.Whether even up according to angle domain common image gathers and accuracy requirement to speed, judge whether also to need again iteration (as Fig. 8).If also do not evened up, then repeat flow process above, until meet accuracy requirement; Otherwise carry out the velocity inversion of lower one deck.Successively analyze from top to bottom, until all layer analysis are complete.Fig. 7 is the schematic diagram of migration result after one embodiment of the invention medium velocity upgrades; Visible to utilize Tomography Velocity field to carry out offseting the result obtained more clear, and underground structure is recovered well.Fig. 8 is the schematic diagram of the angle domain common image gathers after one embodiment of the invention medium velocity upgrades; The lineups of the angle domain common image gathers in figure are all well evened up, and velocity analysis is complete, and the precision of inversion speed field is higher.Flow process terminates.
Fig. 9 is the schematic diagram of the correct velocity field of sandconglomerate body model in one embodiment of the invention; Figure 10 is rate curve comparison diagram before and after one embodiment of the invention medium velocity upgrades.Can find that Tomography Velocity and true velocity are substantially identical by contrast, velocity error, in the scope that systematic error allows, proves validity and the accuracy of the method.
The constraint chromatography velocity modeling method based on Inverse ray tracing in the present invention, based on Inverse ray tracing method, theoretical based on tomographic inversion, finally setting up a migration velocity field is prestack depth migration service.The method is the ray tomography velocity modeling method based on Angle Domain Common Image Gather, based on Inverse ray tracing, ask for sensitivity matrix, extracts whilst on tour residual error from the residue curvature of imaging road collection.Simultaneously in order to improve the precision of inverting, when inversion equation is set up, logging speed is added regularization matrix, to reduce the multisolution of inverting.
Claims (8)
1. based on the constraint chromatography velocity modeling method of Inverse ray tracing, it is characterized in that, should comprise based on the constraint chromatography velocity modeling method of Inverse ray tracing:
Step 1, extracts angle domain common image gathers;
Step 2, automatic Fitting pickup angle gathers residue curvature, utilizes the transformational relation of whilst on tour residual sum degree of depth residual error to ask for whilst on tour residual error;
Step 3, utilizes Inverse ray tracing method, and recording ray arc length, asks for sensitivity matrix;
Step 4, asks for slowness disturbance by LSQR method, renewal speed model; And
Step 5, prestack depth migration and extraction angle domain common image gathers, iteration carries out abovementioned steps, until all keycourses complete renewal.
2. the constraint chromatography velocity modeling method based on Inverse ray tracing according to claim 1, it is characterized in that, in step 1, the rate pattern obtained by general velocity analysis, as initial velocity model, uses the skew of relief surface Gaussian beam carry out migration imaging and extract angle domain common image gathers.
3. the constraint chromatography velocity modeling method based on Inverse ray tracing according to claim 1, it is characterized in that, in step 2, on the basis of conventional automatic pick method, utilize the relational expression of ADCIGs road collection residue curvature and angle, residue curvature corresponding to each angle is picked up in matching.
4. the constraint chromatography velocity modeling method based on Inverse ray tracing according to claim 1, it is characterized in that, in step 3, migrated section picks up interface, layer position, set up ray tracing and just drilling required reference velocity model, then utilize antidromic illumination raytracing scheme to ask for sensitivity matrix.
5. the constraint chromatography velocity modeling method based on Inverse ray tracing according to claim 4, it is characterized in that, in step 3, every bar directions of rays of ray tracing is corresponding with the shooting angle of angle domain common image gathers, carry out Inverse ray tracing using the imaging point at interface as the starting point of ray, make it to match with observation data by adjustment shooting angle.
6. the constraint chromatography velocity modeling method based on Inverse ray tracing according to claim 1, is characterized in that, in step 4, by whilst on tour residual error, sensitivity matrix, sets up the inversion equation group adding regularization factors; Utilize LSQR method solving equation group and renewal speed model.
7. the constraint chromatography velocity modeling method based on Inverse ray tracing according to claim 6, it is characterized in that, in step 4, Boundary Using Sonic Logging is utilized to control refutation process and retrain, amendment inversion equation, with log data sound wave curve constraint velocity inversion result, refutation process iteration is carried out, and finally obtains the tomographic inversion method for solving under borehole restraint.
8. the constraint chromatography velocity modeling method based on Inverse ray tracing according to claim 1, it is characterized in that, in steps of 5, the rate pattern that renewal obtains is for prestack depth migration and extract angle domain common image gathers, whether even up according to angle domain common image gathers and accuracy requirement to speed, judge whether also to need iteration again, if also do not evened up, then repeat step 1 to step 5, until meet accuracy requirement, otherwise carry out the velocity inversion of lower one deck, successively analyze from top to bottom, until all layer analysis are complete.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201410100436.8A CN104932021A (en)  20140318  20140318  Constrained tomography speed modeling method based on reverse ray tracing 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201410100436.8A CN104932021A (en)  20140318  20140318  Constrained tomography speed modeling method based on reverse ray tracing 
Publications (1)
Publication Number  Publication Date 

CN104932021A true CN104932021A (en)  20150923 
Family
ID=54119271
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201410100436.8A Pending CN104932021A (en)  20140318  20140318  Constrained tomography speed modeling method based on reverse ray tracing 
Country Status (1)
Country  Link 

CN (1)  CN104932021A (en) 
Cited By (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN106646600A (en) *  20161228  20170510  中国石油化工股份有限公司  Method for shifting seismic depth by utilizing transverse isotropic medium with vertical symmetric axis 
CN107589446A (en) *  20160706  20180116  中国石油化工股份有限公司  The tomography velocity modeling method of wave path is calculated using Gaussian beam 
CN109143366A (en) *  20170627  20190104  中国石油化工股份有限公司  Near surface first arrival tomographic statics method and computer readable storage medium 
CN110967745A (en) *  20180929  20200407  中国石油化工股份有限公司  Depth domain velocity modeling method for igneous rock 
CN111025397A (en) *  20200106  20200417  中国石油化工股份有限公司  Method for obtaining depth domain velocity model by combining seismic data reflected wave and scattered wave 
Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

US20020049540A1 (en) *  20000807  20020425  Dimitri Bevc  Velocity analysis using angledomain common image gathers 
US20090257308A1 (en) *  20080411  20091015  Dimitri Bevc  Migration velocity analysis methods 
US20120075954A1 (en) *  20100924  20120329  CGGVeritas Services (U.S.) Inc.  Device and Method for Calculating 3D Angle Gathers from Reverse Time Migration 
CN102841376A (en) *  20120906  20121226  中国石油大学(华东)  Retrieval method for chromatography speed based on undulating surface 
CN102841375A (en) *  20120906  20121226  中国石油大学(华东)  Method for tomography velocity inversion based on angle domain common imaging gathers under complicated condition 
CN102841379A (en) *  20120906  20121226  中国石油大学(华东)  Method for analyzing prestack time migration and speed based on common scatter point channel set 

2014
 20140318 CN CN201410100436.8A patent/CN104932021A/en active Pending
Patent Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

US20020049540A1 (en) *  20000807  20020425  Dimitri Bevc  Velocity analysis using angledomain common image gathers 
US20090257308A1 (en) *  20080411  20091015  Dimitri Bevc  Migration velocity analysis methods 
US20120075954A1 (en) *  20100924  20120329  CGGVeritas Services (U.S.) Inc.  Device and Method for Calculating 3D Angle Gathers from Reverse Time Migration 
CN102841376A (en) *  20120906  20121226  中国石油大学(华东)  Retrieval method for chromatography speed based on undulating surface 
CN102841375A (en) *  20120906  20121226  中国石油大学(华东)  Method for tomography velocity inversion based on angle domain common imaging gathers under complicated condition 
CN102841379A (en) *  20120906  20121226  中国石油大学(华东)  Method for analyzing prestack time migration and speed based on common scatter point channel set 
NonPatent Citations (2)
Title 

秦宁 等: ""基于角道集的井约束层析速度反演"", 《石油地球物理勘探》 * 
秦宁 等: ""自动拾取的成像空间域走时层析速度反演"", 《石油地球物理勘探》 * 
Cited By (6)
Publication number  Priority date  Publication date  Assignee  Title 

CN107589446A (en) *  20160706  20180116  中国石油化工股份有限公司  The tomography velocity modeling method of wave path is calculated using Gaussian beam 
CN107589446B (en) *  20160706  20190312  中国石油化工股份有限公司  The tomography velocity modeling method of wave path is calculated using Gaussian beam 
CN106646600A (en) *  20161228  20170510  中国石油化工股份有限公司  Method for shifting seismic depth by utilizing transverse isotropic medium with vertical symmetric axis 
CN109143366A (en) *  20170627  20190104  中国石油化工股份有限公司  Near surface first arrival tomographic statics method and computer readable storage medium 
CN110967745A (en) *  20180929  20200407  中国石油化工股份有限公司  Depth domain velocity modeling method for igneous rock 
CN111025397A (en) *  20200106  20200417  中国石油化工股份有限公司  Method for obtaining depth domain velocity model by combining seismic data reflected wave and scattered wave 
Similar Documents
Publication  Publication Date  Title 

CN104330828B (en)  Dessert reservoir forecasting method and forecasting device  
CN104932021A (en)  Constrained tomography speed modeling method based on reverse ray tracing  
CN102937721B (en)  Limited frequency tomography method for utilizing preliminary wave travel time  
CN102841375A (en)  Method for tomography velocity inversion based on angle domain common imaging gathers under complicated condition  
CN103293552B (en)  A kind of inversion method of Prestack seismic data and system  
CN102841379B (en)  Method for analyzing prestack time migration and speed based on common scatter point channel set  
CN102841376A (en)  Retrieval method for chromatography speed based on undulating surface  
CN104536043A (en)  Depth domain overall velocity model combination method and device  
CN104459782B (en)  Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit  
CN105388518A (en)  Centroid frequency and spectral ratio integrated borehole seismic quality factor inversion method  
CN104570125A (en)  Method for improving imaging velocity model precision by utilizing well data  
CN106597533B (en)  A kind of Depth Domain velocity modeling method for mountain front seism processing  
CN102565856B (en)  Nearsurface noise suppression method on basis of wave equation forward modeling  
CN105301636B (en)  The method for building up and device of rate pattern  
CN103913768A (en)  Method and device for modeling superficial layer in earth surface based on seismic wave data  
CN102116869A (en)  Highprecision prestack domain least square migration seismic imaging technology  
CN102901985A (en)  Depth domain layer speed correcting method suitable for undulating surface  
CN103576200A (en)  Low signaltonoise ratio zone shallow wave impedance interface static correction method  
Parsons et al.  Crustal structure of the Colorado Plateau, Arizona: Application of new long‐offset seismic data analysis techniques  
CN101201409B (en)  Method for revising earthquake data phase  
CN104730572A (en)  Diffracted wave imaging method and device based on L0 seminorm  
CN105093281A (en)  Earthquake multiwave modeling method under inverse framework  
CN102877828A (en)  CT (Computed Tomography) imaging method of threedimensional multiwell combined well land  
CN103852789B (en)  Nonlinear chromatography method and device for seismic data  
CN102565852B (en)  Angle domain prestack offset data processing method aiming to detect oilgasbearing property of reservoir 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
C06  Publication  
SE01  Entry into force of request for substantive examination  
C10  Entry into substantive examination  
WD01  Invention patent application deemed withdrawn after publication 
Application publication date: 20150923 

WD01  Invention patent application deemed withdrawn after publication 