CN102338887A - Irregular-size space-variant grid tomography imaging statics correction method - Google Patents

Irregular-size space-variant grid tomography imaging statics correction method Download PDF

Info

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
Application number
CN2010102314876A
Other languages
Chinese (zh)
Other versions
CN102338887B (en
Inventor
李培明
冯泽元
闫智慧
马青坡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN 201010231487 priority Critical patent/CN102338887B/en
Publication of CN102338887A publication Critical patent/CN102338887A/en
Application granted granted Critical
Publication of CN102338887B publication Critical patent/CN102338887B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Irregular size space-variant grid tomography static correcting method
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.
CN 201010231487 2010-07-14 2010-07-14 Irregular-size space-variant grid tomography imaging statics correction method Active CN102338887B (en)

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 (6)

* Cited by examiner, † Cited by third party
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

Citations (6)

* Cited by examiner, † Cited by third party
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
CN1797030A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for plotting structural diagram of contour line of complex normal fault
CN1797028A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for eliminating imaging error of undulating the earth's surface
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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
CN1797030A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for plotting structural diagram of contour line of complex normal fault
CN1797028A (en) * 2004-12-29 2006-07-05 中国石油天然气集团公司 Method for eliminating imaging error of undulating the earth's surface
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)

* Cited by examiner, † Cited by third party
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 (10)

* Cited by examiner, † Cited by third party
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

Also Published As

Publication number Publication date
CN102338887B (en) 2013-09-04

Similar Documents

Publication Publication Date Title
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
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) Cross-well 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 three-dimensional 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 three-dimensional multi-well 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 three-dimensional 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 mid-deep 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) Multi-information 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