CN102338887A - Irregular-size space-variant grid tomography imaging statics correction method - Google Patents
Irregular-size space-variant 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
- model
- horizontal direction
- node
- vertical direction
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses an irregular-size space-variant 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 space-variant 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 near-surface 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 rule-based 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 space-variant.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 stair-stepping.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 space-variant 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 2-4 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+(B-Ai)*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 set-point 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 space-variant 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 unequal-interval 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 set-point 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 large-scale 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 space-variant 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 2-4 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+(B-Ai)*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 set-point 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) | 2010-07-14 | 2010-07-14 | Irregular-size space-variant grid tomography imaging statics correction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010231487 CN102338887B (en) | 2010-07-14 | 2010-07-14 | Irregular-size space-variant grid tomography imaging statics correction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102338887A true CN102338887A (en) | 2012-02-01 |
| CN102338887B CN102338887B (en) | 2013-09-04 |
Family
ID=45514714
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201010231487 Active CN102338887B (en) | 2010-07-14 | 2010-07-14 | Irregular-size space-variant grid tomography imaging statics correction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102338887B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102937721A (en) * | 2012-11-07 | 2013-02-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Limited frequency tomography method for utilizing preliminary wave travel time |
| CN103616722A (en) * | 2013-11-28 | 2014-03-05 | 中国石油天然气股份有限公司 | First-motion wave traveling time picking-up method and device |
| CN104181593A (en) * | 2014-08-28 | 2014-12-03 | 中国石油天然气集团公司 | Three-dimensional ray-free tracing contraflexure wave tomographic imaging method and device |
| CN106353810A (en) * | 2016-08-12 | 2017-01-25 | 中国石油天然气股份有限公司 | Geological lithological difference recognizing method and apparatus |
| CN104459774B (en) * | 2014-11-05 | 2017-04-05 | 中国石油天然气股份有限公司 | A kind of geologic lithology difference identification method and system |
| CN107783185A (en) * | 2017-09-14 | 2018-03-09 | 中国石油天然气股份有限公司 | A kind of processing method and processing device of tomographic statics |
| CN112036008A (en) * | 2020-07-28 | 2020-12-04 | 中国石油天然气集团有限公司 | Method and system for recording inversion speed by using hybrid model |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5394325A (en) * | 1993-04-07 | 1995-02-28 | Exxon Production Research Company | Robust, efficient three-dimensional finite-difference traveltime calculations |
| US20040073370A1 (en) * | 2002-08-20 | 2004-04-15 | Shivaji Dasgupta | Use of drill bit energy for tomographic modeling of near surface layers |
| CN1797028A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for eliminating imaging error of undulating the earth's surface |
| CN1797030A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for plotting structural diagram of contour line of complex normal fault |
| CN101021568A (en) * | 2007-02-07 | 2007-08-22 | 匡斌 | Three-dimensional integral prestack depth migration method |
| CN101086534A (en) * | 2006-06-07 | 2007-12-12 | 中国石油集团东方地球物理勘探有限责任公司 | Demodulator probe secondary localization method |
-
2010
- 2010-07-14 CN CN 201010231487 patent/CN102338887B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5394325A (en) * | 1993-04-07 | 1995-02-28 | Exxon Production Research Company | Robust, efficient three-dimensional finite-difference traveltime calculations |
| US20040073370A1 (en) * | 2002-08-20 | 2004-04-15 | Shivaji Dasgupta | Use of drill bit energy for tomographic modeling of near surface layers |
| CN1797028A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for eliminating imaging error of undulating the earth's surface |
| CN1797030A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for plotting structural diagram of contour line of complex normal fault |
| CN101086534A (en) * | 2006-06-07 | 2007-12-12 | 中国石油集团东方地球物理勘探有限责任公司 | Demodulator probe secondary localization method |
| CN101021568A (en) * | 2007-02-07 | 2007-08-22 | 匡斌 | Three-dimensional integral prestack depth migration method |
Non-Patent Citations (4)
| Title |
|---|
| 于岚: "层析反演静校正技术应用研究", 《中国石油大学工程硕士学位论文》, 31 December 2008 (2008-12-31) * |
| 何光明等: "非线性层析静校正在川西地区资料处理中的应用", 《石油物探》, vol. 45, no. 01, 28 February 2006 (2006-02-28), pages 88 - 92 * |
| 冯泽元等: "利用层析反演技术解决山地复杂区静校正问题", 《石油物探》, vol. 44, no. 03, 30 June 2005 (2005-06-30), pages 284 - 287 * |
| 赵妍: "层析成像静校正方法研究", 《中国石油大学硕士学位论文》, 31 May 2009 (2009-05-31) * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102937721A (en) * | 2012-11-07 | 2013-02-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Limited frequency tomography method for utilizing preliminary wave travel time |
| CN102937721B (en) * | 2012-11-07 | 2015-07-08 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Limited frequency tomography method for utilizing preliminary wave travel time |
| CN103616722A (en) * | 2013-11-28 | 2014-03-05 | 中国石油天然气股份有限公司 | First-motion wave traveling time picking-up method and device |
| CN103616722B (en) * | 2013-11-28 | 2016-04-06 | 中国石油天然气股份有限公司 | A kind of travel-times of seismic first breaks pick-up method and device |
| CN104181593A (en) * | 2014-08-28 | 2014-12-03 | 中国石油天然气集团公司 | Three-dimensional ray-free tracing contraflexure wave tomographic imaging method and device |
| CN104459774B (en) * | 2014-11-05 | 2017-04-05 | 中国石油天然气股份有限公司 | A kind of geologic lithology difference identification method and system |
| CN106353810A (en) * | 2016-08-12 | 2017-01-25 | 中国石油天然气股份有限公司 | Geological lithological difference recognizing method and apparatus |
| CN106353810B (en) * | 2016-08-12 | 2018-10-16 | 中国石油天然气股份有限公司 | A kind of geologic lithology difference identification method and device |
| CN107783185A (en) * | 2017-09-14 | 2018-03-09 | 中国石油天然气股份有限公司 | A kind of processing method and processing device of tomographic statics |
| CN107783185B (en) * | 2017-09-14 | 2019-05-07 | 中国石油天然气股份有限公司 | A kind of processing method and processing device of tomographic statics |
| CN112036008A (en) * | 2020-07-28 | 2020-12-04 | 中国石油天然气集团有限公司 | Method and system for recording inversion speed by using hybrid model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102338887B (en) | 2013-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106597533B (en) | A kind of Depth Domain velocity modeling method for mountain front seism processing | |
| CN102338887B (en) | Irregular-size space-variant grid tomography imaging statics correction method | |
| CN104459782B (en) | Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit | |
| CN104133245B (en) | The static correcting method and system of a kind of seismic data | |
| CN108064348A (en) | Chromatography conversion method when a kind of earthquake based on two spots ray tracing is walked | |
| Waheed et al. | A fast sweeping algorithm for accurate solution of the tilted transversely isotropic eikonal equation using factorization | |
| CN108005646B (en) | Stratum anisotropic resistivity extraction method based on electromagnetic wave logging while drilling data | |
| CN108710148B (en) | The steady phase prestack depth migration method in three-dimensional dip domain and device | |
| CN107843922A (en) | One kind is based on seismic first break and the united chromatography imaging method of Travel time | |
| CN106353793A (en) | Cross-well seismic tomography inversion method on basis of travel time incremental bilinear interpolation ray tracing | |
| CN113552625B (en) | Multi-scale full waveform inversion method for conventional land-domain seismic data | |
| Zhou | Multiscale deformable-layer tomography | |
| CN104360396B (en) | A kind of three kinds of preliminary wave Zoumaling tunnel methods of TTI medium between offshore well | |
| BR112015000879B1 (en) | System and method for modeling migration speed | |
| CN108303736B (en) | Ray tracing forward method for shortest path of anisotropic TI medium | |
| CN111221035A (en) | Seismic reflection wave slope and gravity anomaly data joint inversion method | |
| CN102877828A (en) | CT (Computed Tomography) imaging method of three-dimensional multi-well combined well land | |
| CN109884700B (en) | Multi-information fusion seismic velocity modeling method | |
| CN105137479A (en) | Method and device for calculating number of coverage times of surface elements | |
| CN102866422A (en) | Depth domain geological entity model generation method | |
| CN104155694A (en) | Residual static correction method of reflection converted shear wave common geophone stack section | |
| CN103217715B (en) | Multiple dimensioned regular grid Static Correction of Tomographic Inversion method | |
| CN105204064A (en) | Mixed domain Fourier finite difference migration method based on coefficient optimization | |
| CN1292263C (en) | Ray traction in earthquake prospection | |
| CN106896408A (en) | A kind of angle domain prestack time migration method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |