CN102338887A  Irregularsize spacevariant grid tomography imaging statics correction method  Google Patents
Irregularsize spacevariant grid tomography imaging statics correction method Download PDFInfo
 Publication number
 CN102338887A CN102338887A CN2010102314876A CN201010231487A CN102338887A CN 102338887 A CN102338887 A CN 102338887A CN 2010102314876 A CN2010102314876 A CN 2010102314876A CN 201010231487 A CN201010231487 A CN 201010231487A CN 102338887 A CN102338887 A CN 102338887A
 Authority
 CN
 China
 Prior art keywords
 grid
 horizontal
 node
 vertical
 speed
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Granted
Links
 230000003068 static Effects 0.000 title claims abstract description 17
 238000003325 tomography Methods 0.000 title claims abstract description 6
 238000003384 imaging method Methods 0.000 title claims abstract description 5
 239000011295 pitches Substances 0.000 claims description 7
 230000001788 irregular Effects 0.000 claims description 5
 230000015572 biosynthetic process Effects 0.000 claims description 3
 230000000875 corresponding Effects 0.000 claims description 3
 238000005755 formation reactions Methods 0.000 claims description 3
 280000398338 Seismic companies 0.000 claims 1
 230000000694 effects Effects 0.000 abstract description 4
 238000004587 chromatography analysis Methods 0.000 description 6
 238000005516 engineering process Methods 0.000 description 2
 238000003874 inverse correlation nuclear magnetic resonance spectroscopy Methods 0.000 description 2
 238000005553 drilling Methods 0.000 description 1
 238000005070 sampling Methods 0.000 description 1
Abstract
The invention discloses an irregularsize spacevariant grid tomography imaging statics correction method in geophysical exploration. The method is characterized by: arranging different grid lengths along a horizontal direction during performing mesh generation; arranging novel nodes on four sides of the grid; performing equal layout on a same distance in a horizontal direction or a vertical direction; or carrying out the equal layout of the nodes on the each grid side in the horizontal direction or the vertical direction; initializing a model speed and performing ray path forward modeling; calculating a difference value between actual picked first break time and ray travelling time and the speed of an inversion model; completing imaging statics correction till inversion speed field is stable. By using the method of the invention, a near surface characteristic of an actual stratum can be better described; propagation rules of a first break wave on the near surface can be better reflected; forward and inversion precision of the speed model can be raised and a statics correction effect can be raised too.
Description
Technical field
The present invention relates to geophysical exploration method, is a kind of irregular size spacevariant grid tomography static correcting method.
Background technology
Static corrections processing mainly is to solve seismic event near surface medium variable effect hourage reflection wave stacking image problem in the seismic prospecting.Primary wave reflection near surface stratigraphic fluctuation utilizes primary wave to ask for the major technology means that static correction value becomes static correction.The nearsurface model that utilizes the method for chromatography FORWARD AND INVERSE PROBLEMS can set up relatively accurately, can response speed change in length and breadth, chromatography method has obtained using widely and having obtained good effect in many complex area aborning in recent years.
Be applied to all rulebased and fixed measure grid of conventional chromatography method in the static correction at present, and suppose that the speed in the grid is constant, size of mesh opening can not spacevariant.This condition restriction the authenticity and the rationality of tomographic inversion model.Especially rise and fall or stratum when having pinching when the face of land and subsurface interface have, the equidimension regular grid can't the accurate description geologic body or the boundary position of body of velocity, and the geologic body that it is described or the border of body of velocity can only be stairstepping.Like the subdivision result of the conventional chromatography method of the grid among Fig. 1 to model, continuous black line is the geologic body or the body of velocity interface of realistic model, the subdivision result the interface location place with have than mistake really.
Summary of the invention
The object of the invention is to provide a kind of can accurately describe the interfacial configuration of geologic body or body of velocity, the time error of ray tracing, improves the irregular size spacevariant grid tomography static correcting method of the precision of chromatography static correction inversion speed field.
The present invention realizes through following technological means:
1) acquiring seismic data carries out the mesh generation of horizontal direction and vertical direction after the processing, the horizontal direction grid node is horizontal direction meshes number+1, and described horizontal direction is provided with different grid length;
Described horizontal direction is provided with grid length and equals 1 times of track pitch in the violent position grid length of velocity variations, equals 24 times of track pitch in the inviolent position grid length of velocity variations, and grid length can be greater than 5 times of track pitches on the model both sides.
The mesh generation of described vertical direction calculates according to following formula:
Zij＝Ai+(BAi)*j/(N+1)，
In the formula: Ai physical points elevation makes up ground surface, and i is a horizontal direction grid node index,
J is a vertical direction grid node index,
Meshes number N,
Grid node j＜=N+1,
The bottom boundary B of model.
If known formation interface elevation, the j of vertical direction grid elevation Zij is a setpoint that is not more than the vertical direction meshes number, overlaps with this interface.
Near ground surface, in vertical direction, it is littler that the user can be provided with mesh spacing, and in the model deep, the user can be provided with mesh spacing more greatly.
2) on four limits of a grid, lay new node;
Lay new node and be in the horizontal direction or vertical direction, the distance between the adjacent node equates to lay, or:
In the horizontal direction or vertical direction, the node number on each grid edge is to equate to lay.
3) with the model velocity initialization;
Described model velocity initialization is with the minimum speed and the gradient factor, calculates the speed of each grid,
Or with the maximal rate of correspondence at the bottom of the corresponding minimum speed of ground surface and the model, interpolation model midrange speed.
4) raypath is just drilled;
It is to search apart from the shot point ray of all nodes recently that coexists in the grid that described raypath is just being drilled, and is new focal point with these nodes again, searches and coexist other node raypath of a grid of these new focal points; Search new focal point to the raypath of short whilst on tour between the shot point;
5) calculate actual difference of picking up first break time and ray hourage, the speed of inverse model;
6) iteration 4) and 5), stable up to the velocity field of inverting, accomplish the imaging static correction.
The present invention adopts irregular size spacevariant grid that model is dispersed can describe the near surface characteristic of actual formation better, can reflect the propagation law of primary wave near surface better, has improved rate pattern FORWARD AND INVERSE PROBLEMS precision and static correction effect.
Description of drawings
Fig. 1 is conventional grid model subdivision result;
Use the regular grid subdivision, in the face of land and model layer at the interface, subdivision result and actual has than big difference.
Fig. 2 is the chromatography subdivision result of the present invention to a model.
Embodiment
Specify instance below in conjunction with accompanying drawing.
The present invention adopts following steps to realize:
1, model meshes subdivision:
Behind the data sampling, the coordinate range of all acceptance points and shot point is confirmed the length of model, carries out horizontal direction unequalinterval subdivision; As shown in Figure 2,100 meters of track pitches, on the model left side, 500 meters of horizontal direction mesh spacings are at model middle part, 100 meters of mesh spacings.In vertical direction, confirm a ground surface with the elevation of acceptance point and shot point, the elevation that guarantees all acceptance points and shot point is all under this ground surface.Provide the bottom boundary elevation of model and the meshes number of vertical direction, linear interpolation goes out the height value of four angle points of intermediate mesh.
If known certain bed boundary elevation; Vertical direction grid elevation Zij value can be set, and (i is a horizontal direction grid node index; J is a setpoint that is not more than vertical direction grid node number); Overlap with this interface, adjust simultaneously among other Zij, guarantee from small to large in vertical direction Zij value.
Near ground surface, in vertical direction, it is littler that the user can be provided with mesh spacing, and in the model deep, the user can be provided with mesh spacing more greatly, and is as shown in Figure 2.
2, node is laid:
On four limits of a grid, lay new node.The present invention provides two kinds of nodes to lay scheme.A) equidistantly lay: promptly in the horizontal direction or the vertical direction of each grid, between adjacent two nodes apart from the time equate; B) etc. number is laid: in each grid, in the horizontal direction or vertical direction, the number of node equates.In these two kinds of schemes, all the minimum node spacing can be set, as 0.25 meter, purpose is to reduce the node number, improves arithmetic speed.These nodes can only be passed in the path of ray when just drilling.
3, model velocity initialization:
Provide the minimum speed on the model face of land, provide the maximal rate at the bottom of the gradient factor or the model again, calculate the speed in each grid in the middle of the model; Or with the maximal rate of correspondence at the bottom of the corresponding minimum speed of ground surface and the model, interpolation model midrange speed.
4, the shortest raypath search:
Search is new focal point with these nodes apart from the shot point ray of all nodes recently that coexists in the grid again, searches and coexist other node raypath of a grid of these new focal points; Be as the criterion with minimum time, search for them, write down this raypath ray length and hourage in each grid to the raypath of short whilst on tour between the shot point;
5, model inversion:
Calculate the difference of actual ray hourage of picking up first break time and search; Make up largescale system of equations again, the speed of utilizing this difference and the ray travel distance in grid to come inverse model; Upgrade model velocity.
6, iteration 4) and 5), stable up to the velocity field of inverting, handle the static correction result that obtains forming images.
Claims (5)
1. irregular size spacevariant grid tomography static correcting method, characteristic are to realize through following technological means:
1) acquiring seismic data carries out the mesh generation of horizontal direction and vertical direction after the processing, the horizontal direction grid node is horizontal direction meshes number+1, and described horizontal direction is provided with different grid length;
2) on four limits of a grid, lay new node;
Lay new node and be in the horizontal direction or vertical direction, the distance between the adjacent node equates to lay, or:
In the horizontal direction or vertical direction, the node number on each grid edge is to equate to lay;
3) with the model velocity initialization;
4) raypath is just drilled;
5) calculate actual difference of picking up first break time and ray hourage, the speed of inverse model;
6) iteration 4) and 5), stable up to the velocity field of inverting, accomplish the imaging static correction.
2. according to the method for claim 1; The described horizontal direction of step 1) is provided with grid length and equals 1 times of track pitch in the violent position grid length of velocity variations; Equal 24 times of track pitch in the inviolent position grid length of velocity variations, grid length can be greater than 5 times of track pitches on the model both sides.
3. according to the method for claim 1, the mesh generation of the described vertical direction of step 1) calculates according to following formula:
Zij＝Ai+(BAi)*j/(N+1)，
In the formula: Ai physical points elevation makes up ground surface, and i is a horizontal direction grid node index,
J is the vertical direction grid node,
Meshes number N,
Grid node j＜=N+1, if known formation interface elevation, j is a setpoint that is not more than the vertical direction meshes number, overlaps with this interface;
The bottom boundary B of model.
4. according to the method for claim 1, the described model velocity initialization of step 3) is with the minimum speed and the gradient factor, calculates the speed of each grid;
Or, with the maximal rate of correspondence at the bottom of the corresponding minimum speed of ground surface and the model, interpolation model midrange speed.
5. according to the method for claim 1; It is to search apart from the shot point ray of all nodes recently that coexists in the grid that the described raypath of step 4) is just being drilled; Be new focal point with these nodes again, search and coexist other node raypath of a grid of these new focal points; Search new focal point to the raypath of short whilst on tour between the shot point.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN 201010231487 CN102338887B (en)  20100714  20100714  Irregularsize spacevariant grid tomography imaging statics correction method 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN 201010231487 CN102338887B (en)  20100714  20100714  Irregularsize spacevariant grid tomography imaging statics correction method 
Publications (2)
Publication Number  Publication Date 

CN102338887A true CN102338887A (en)  20120201 
CN102338887B CN102338887B (en)  20130904 
Family
ID=45514714
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN 201010231487 Active CN102338887B (en)  20100714  20100714  Irregularsize spacevariant grid tomography imaging statics correction method 
Country Status (1)
Country  Link 

CN (1)  CN102338887B (en) 
Cited By (6)
Publication number  Priority date  Publication date  Assignee  Title 

CN102937721A (en) *  20121107  20130220  中国石油集团川庆钻探工程有限公司地球物理勘探公司  Limited frequency tomography method for utilizing preliminary wave travel time 
CN103616722A (en) *  20131128  20140305  中国石油天然气股份有限公司  Firstmotion wave traveling time pickingup method and device 
CN104181593A (en) *  20140828  20141203  中国石油天然气集团公司  Threedimensional rayfree tracing contraflexure wave tomographic imaging method and device 
CN106353810A (en) *  20160812  20170125  中国石油天然气股份有限公司  Geological lithological difference recognizing method and apparatus 
CN104459774B (en) *  20141105  20170405  中国石油天然气股份有限公司  A kind of geologic lithology difference identification method and system 
CN107783185A (en) *  20170914  20180309  中国石油天然气股份有限公司  A kind of processing method and processing device of tomographic statics 
Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

US5394325A (en) *  19930407  19950228  Exxon Production Research Company  Robust, efficient threedimensional finitedifference traveltime calculations 
US20040073370A1 (en) *  20020820  20040415  Shivaji Dasgupta  Use of drill bit energy for tomographic modeling of near surface layers 
CN1797030A (en) *  20041229  20060705  中国石油天然气集团公司  Method for plotting structural diagram of contour line of complex normal fault 
CN1797028A (en) *  20041229  20060705  中国石油天然气集团公司  Method for eliminating imaging error of undulating the earth's surface 
CN101021568A (en) *  20070207  20070822  匡斌  Threedimensional integral prestack depth migration method 
CN101086534A (en) *  20060607  20071212  中国石油集团东方地球物理勘探有限责任公司  Demodulator probe secondary localization method 

2010
 20100714 CN CN 201010231487 patent/CN102338887B/en active Active
Patent Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

US5394325A (en) *  19930407  19950228  Exxon Production Research Company  Robust, efficient threedimensional finitedifference traveltime calculations 
US20040073370A1 (en) *  20020820  20040415  Shivaji Dasgupta  Use of drill bit energy for tomographic modeling of near surface layers 
CN1797030A (en) *  20041229  20060705  中国石油天然气集团公司  Method for plotting structural diagram of contour line of complex normal fault 
CN1797028A (en) *  20041229  20060705  中国石油天然气集团公司  Method for eliminating imaging error of undulating the earth's surface 
CN101086534A (en) *  20060607  20071212  中国石油集团东方地球物理勘探有限责任公司  Demodulator probe secondary localization method 
CN101021568A (en) *  20070207  20070822  匡斌  Threedimensional integral prestack depth migration method 
NonPatent Citations (4)
Title 

于岚: "层析反演静校正技术应用研究", 《中国石油大学工程硕士学位论文》, 31 December 2008 (20081231) * 
何光明等: "非线性层析静校正在川西地区资料处理中的应用", 《石油物探》, vol. 45, no. 01, 28 February 2006 (20060228), pages 88  92 * 
冯泽元等: "利用层析反演技术解决山地复杂区静校正问题", 《石油物探》, vol. 44, no. 03, 30 June 2005 (20050630), pages 284  287 * 
赵妍: "层析成像静校正方法研究", 《中国石油大学硕士学位论文》, 31 May 2009 (20090531) * 
Cited By (10)
Publication number  Priority date  Publication date  Assignee  Title 

CN102937721A (en) *  20121107  20130220  中国石油集团川庆钻探工程有限公司地球物理勘探公司  Limited frequency tomography method for utilizing preliminary wave travel time 
CN102937721B (en) *  20121107  20150708  中国石油集团川庆钻探工程有限公司地球物理勘探公司  Limited frequency tomography method for utilizing preliminary wave travel time 
CN103616722A (en) *  20131128  20140305  中国石油天然气股份有限公司  Firstmotion wave traveling time pickingup method and device 
CN103616722B (en) *  20131128  20160406  中国石油天然气股份有限公司  A kind of traveltimes of seismic first breaks pickup method and device 
CN104181593A (en) *  20140828  20141203  中国石油天然气集团公司  Threedimensional rayfree tracing contraflexure wave tomographic imaging method and device 
CN104459774B (en) *  20141105  20170405  中国石油天然气股份有限公司  A kind of geologic lithology difference identification method and system 
CN106353810A (en) *  20160812  20170125  中国石油天然气股份有限公司  Geological lithological difference recognizing method and apparatus 
CN106353810B (en) *  20160812  20181016  中国石油天然气股份有限公司  A kind of geologic lithology difference identification method and device 
CN107783185A (en) *  20170914  20180309  中国石油天然气股份有限公司  A kind of processing method and processing device of tomographic statics 
CN107783185B (en) *  20170914  20190507  中国石油天然气股份有限公司  A kind of processing method and processing device of tomographic statics 
Also Published As
Publication number  Publication date 

CN102338887B (en)  20130904 
Similar Documents
Publication  Publication Date  Title 

CN102338887B (en)  Irregularsize spacevariant grid tomography imaging statics correction method  
CN104459782B (en)  Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit  
CN102914791B (en)  Kirchhoff prestack time migration method for processing seismic data of undulating surface  
CN104133245B (en)  The static correcting method and system of a kind of seismic data  
CN106353793A (en)  Crosswell seismic tomography inversion method on basis of travel time incremental bilinear interpolation ray tracing  
CN106597533B (en)  A kind of Depth Domain velocity modeling method for mountain front seism processing  
Waheed et al.  A fast sweeping algorithm for accurate solution of the tilted transversely isotropic eikonal equation using factorization  
CN104360396B (en)  A kind of three kinds of preliminary wave Zoumaling tunnel methods of TTI medium between offshore well  
CN108710148B (en)  The steady phase prestack depth migration method in threedimensional dip domain and device  
CN102866422B (en)  A kind of depth domain geological entity model generation method  
CN105301636B (en)  The method for building up and device of rate pattern  
CN107843922A (en)  One kind is based on seismic first break and the united chromatography imaging method of Travel time  
CN102877828A (en)  CT (Computed Tomography) imaging method of threedimensional multiwell combined well land  
CN105974479A (en)  Chromatographic 2D/3D anisotropic depth domain velocity modeling method for GPU (Graphics Processing Unit) spatial grids  
CN108680964A (en)  A kind of normalization weight magnetoelectricity shake joint inversion method based on structural constraint  
CN1292263C (en)  Ray traction in earthquake prospection  
CN108064348A (en)  Chromatography conversion method when a kind of earthquake based on two spots ray tracing is walked  
CN103217715B (en)  Multiple dimensioned regular grid Static Correction of Tomographic Inversion method  
CN108303736B (en)  Ray tracing forward method for shortest path of anisotropic TI medium  
CN102841374A (en)  Pseudo threedimensional fast microseism forward modeling method based on scanning surface forward modeling  
CN106896408A (en)  A kind of angle domain prestack time migration method  
CN109655890A (en)  A kind of shallow middeep strata joint chromatography inversion speed modeling method of Depth Domain and system  
CN104155694A (en)  Residual static correction method of reflection converted shear wave common geophone stack section  
CN109884700B (en)  Multiinformation fusion seismic velocity modeling method  
CN102466818B (en)  A kind of well earthquake data that utilizes is to the method for anisotropic medium imaging 
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 