CN106383361B - A kind of method of speed data grid updating - Google Patents
A kind of method of speed data grid updating Download PDFInfo
- Publication number
- CN106383361B CN106383361B CN201610786593.8A CN201610786593A CN106383361B CN 106383361 B CN106383361 B CN 106383361B CN 201610786593 A CN201610786593 A CN 201610786593A CN 106383361 B CN106383361 B CN 106383361B
- Authority
- CN
- China
- Prior art keywords
- grid
- new
- track
- cross
- coordinate
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 238000005070 sampling Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 102100023943 Arylsulfatase L Human genes 0.000 description 1
- 101000975827 Homo sapiens Arylsulfatase L Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Image Processing (AREA)
Abstract
The present invention provides a kind of methods of speed data grid updating, it the described method comprises the following steps: new, the old two nested grids file of input, pass through the rule that the coordinate value using corresponding position in new and old two nested grid is constant, old mesh parameter is converted into actual coordinate, again coordinate value by new grid conversion at corresponding parameter, thus achieve the purpose that update grid.Method of the invention can be suitable for the seismic data process of earth's surface and subsurface picture complexity, have the characteristics that calculate step simplicity and high treating effect, have broad application prospects in current seismic data high Precision Processing.
Description
Technical field
The invention belongs to petroleum gas seismic prospectings to investigate field, more particularly, be related to one kind applied to petroleum
The method for the speed data grid updating that the seism processing of seismic exploration is explained.
Background technique
In general, many work in seismic data process are shot point coordinate, geophone station coordinate and root based on seismic channel
It is carried out according to processing grid defined in these coordinates, these grids are the bridges connected between field actual coordinate and acquisition data
The definition of beam, referred to as field layout.The field document number of each seismic channel is had recorded in the trace header of field seismic data
And Taoist monastic name, and the coordinate information of shot point and geophone station is then recorded in field class report.Geometry definition is exactly with field file
Number and record Taoist monastic name be index, assign the correct shot point coordinate of each seismic channel, geophone station coordinate, and according to the grid of definition
Situation calculates the center point coordinate and face element serial number of each face element, and these data are recorded on earthquake trace header or observation system
In database.
Since the imaging of seismic data later period is mainly related with speed, and speed is according to previously defined grid
Face element mode is expressed, therefore the correctness of speed grid is most important.And in some actual treatments, it often refers to same
The speed data of position needs to change its grid to be applied on different seismic datas, it is therefore desirable to pass through renewal speed
The grids of data is realized." speed data " it is that a series of record same spatial locations different reflection intervals or depth are corresponding
The data of velocity amplitude sampling point；" seismic data " it is that a series of record same spatial locations different reflection intervals or depth receive
The amplitudes sampling point information such as seismic reflection, refraction data.
Summary of the invention
Aiming at the problems existing in the prior art, it is an object of the invention to solve in deficiencies of the prior art
At least one of.
The purpose of the present invention is to provide one kind to be defined by speed grid control parameter, using based on the same coordinate system
Under, the constant thinking of the coordinate of face element ultimately forms speed data grid updating method.
To achieve the goals above, the present invention provides a kind of method of speed data grid updating, the method includes
Following steps: old speed data, the old corresponding old mesh parameter of speed data and new mesh parameter A, are inputted；B, according to old speed
Any record Taoist monastic name of the trace header of data and old mesh parameter calculate the corresponding face element coordinate of any record Taoist monastic name；C, it presses
According to the principle that face element coordinate under the same coordinate system is constant, joined according to any corresponding face element coordinate of record Taoist monastic name and new grid
Number calculates updated speed data trace header；D, repeating said steps B and C, until all records of the trace header of old speed data
Taoist monastic name has been handled.
In an exemplary embodiment of the method for speed data grid updating of the present invention, the old mesh parameter and new
Mesh parameter includes: main profile starting symbol, main profile termination number, the termination of main profile increment, cross-track starting symbol, cross-track
Number, cross-track increment, grid starting point X-coordinate, grid starting point Y coordinate, the projection of main profile grid spacing in the X direction away from
From, main profile grid spacing projector distance in the Y direction, cross-track grid spacing projector distance in the X direction and
The projector distance of cross-track grid spacing in the Y direction.
In an exemplary embodiment of the method for speed data grid updating of the present invention, the step B includes:
B1, the point number nx for calculating each main profile:
Nx=1+ (cross-track termination number-cross-track starting symbol)/cross-track increment.
The main profile number of lines steps before main profile where b2, calculating any record Taoist monastic name:
Steps=int ((any record Taoist monastic name-cross-track starting symbol)/nx).
Cross-track the number of lines stepx, stepx before cross-track where b3, calculating any record Taoist monastic name
=(any record Taoist monastic name-cross-track starting symbol)-steps*nx.
B4, any record Taoist monastic name corresponding face element coordinate (x, y) in old grid is calculated:
X-coordinate:
Projector distance+the stepx* of x=grid starting point X-coordinate+steps* main profile grid spacing in the X direction, which is got in touch with, to be surveyed
The projector distance of wire grid spacing in the X direction.
Y coordinate:
Projector distance+the stepx* of y=grid starting point Y coordinate+steps* main profile grid spacing in the Y direction, which is got in touch with, to be surveyed
The projector distance of wire grid spacing in the Y direction.
In an exemplary embodiment of the method for speed data grid updating of the present invention, the step C includes:
C1, updated main profile number is calculated according to any corresponding face element coordinate of record Taoist monastic name and new mesh parameter
With cross-track number.Specifically:
Dx=(the new grid starting point X-coordinate of current X-coordinate -).
Dy=(the new grid starting point Y coordinate of current Y coordinate -).
Dinv=1/ is (between the new grid cross-track grid of projector distance * of new grid main profile grid spacing in the X direction
Away from the new grid cross-track of projector distance * of the new grid main profile grid spacing of projector distance-in the Y direction in the Y direction
The projector distance of grid spacing in the X direction).
Del=(the new grid cross-track of projector distance-dy* of the new grid cross-track grid spacing of dx* in the Y direction
The projector distance * dinv of grid spacing in the X direction).
New main profile number=new grid main profile starting symbol+new grid main profile increment * del.
Del '=int ((the new grid main profile of projector distance-dx* of the new grid main profile grid spacing of dy* in the X direction
The projector distance of grid spacing in the Y direction) * dinv).
New cross-track number=new grid cross-track starting symbol+new grid cross-track increment * del '.
C2, the new main profile number being calculated according to step c1, new cross-track number calculate new note with new mesh parameter
Record Taoist monastic name.Specifically:
Steps=(the new grid main profile starting symbol of new main profile number -)/new grid main profile increment.
Stepx=(the new grid cross-track starting symbol of new cross-track number -)/new grid cross-track increment.
New record Taoist monastic name=steps* (1+ (new grid cross-track termination number-new grid cross-track starting symbol)/new net
Lattice cross-track increment)+stepx+1.
In an exemplary embodiment of the method for speed data grid updating of the present invention, after being updated in the step C
Speed data trace header include new main profile number, new cross-track number and new record Taoist monastic name.
In an exemplary embodiment of the method for speed data grid updating of the present invention, the record Taoist monastic name can be
No. CDP, No. CMP, No. CRP or No. CIP.
In an exemplary embodiment of the method for speed data grid updating of the present invention, the method can also include
After step D, new speed merging is carried out, the merging step includes: to reject to overlapping block, is carried out to blank block
Interpolation, to be merged into the speed data based on new grid.
Compared with prior art, the beneficial effect comprise that the ground of earth's surface and subsurface picture complexity can be suitable for
Data processing is shaken, has the characteristics that calculate step simplicity and high treating effect, have in current seismic data high Precision Processing
Wide application prospect.
Detailed description of the invention
Fig. 1 shows the flow chart of the speed data grid updating method of an exemplary embodiment of the present invention.
Fig. 2 shows old grid schematic diagrames.
Fig. 3 is after carrying out grid updating to Fig. 2 using the speed data grid updating method of exemplary embodiment of the present
New grid schematic diagram.
Description of symbols:
1- cross-track (XLINE), 2- main profile (INLINE), O point are old grid starting point, and O ' is new grid starting point.
Specific embodiment
It hereinafter, will be in conjunction with exemplary embodiment and attached drawing the speed data grid updating that the present invention will be described in detail
Method.
For ease of understanding, applicant first to this application involves basic concepts simply illustrated: exploration is on ground
Face lays the form of grid to realize, such as one piece of survey area, can use one of a series of adjacent small grid compositions
Big grid, small grid side length 50m*50m, naturally it is also possible to be 10m*20m, by carrying out these parameter designings before exploring.These are small
Grid is referred to as " face element ", corresponding CDP/CMP number (number in order), with the center position coordinates of face element it is equivalent this
Exploration information in a face element numbers recording with depth or time for corresponding depth direction with CDP/CMP in each face element
A series of velocity amplitudes of variation.CDP refers to that common depth point, CMP are common point.
The speed data of same location needs to change its grid to be applied on different seismic datas, it is therefore desirable to
It is realized by the grid of renewal speed data." same location " is that velocity space coordinate is identical, and corresponding is to use CMP (in altogether
Heart point) coordinate description, because the coordinate of speed is not first for recording each location point, but only recording the grid of definition
The coordinate of a start reference position, the coordinates of other positions further according to distance reference position the number of grids and mesh width, with
And the coordinate of start reference position converts, therefore corresponding originally in the speed of the same coordinate, with different grid computings,
Appear in different positions.
The present invention is defined using key parameter, using based under the same coordinate system, the constant thinking of the coordinate of face element, finally
Speed data grid updating method is formed, more efficient grid updating processing can be provided for seismic velocity data processing.Face element is
One rectangle or square in corresponding space level face, it is corresponding a series of from shot point excitation, by the reflection of a certain depth, to inspection
(same road seismic data necessarily corresponds to a big gun to the reflection point earthquake information in the identical face element region in wave point reception recording process
Point, as soon as a reflection point and receiving point), a series of face element interconnected forms the net of a covering certain area
Lattice, each face element correspond to a series of data set for the seismic channel that reflection points are all fallen in the face element, referred to as common midpoint gather (
That is CMP trace gather), corresponding CMP number, and the coordinate namely CMP of the corresponding face element center of each CMP number are sat
Mark.
Basic ideas of the invention are: new, the old two nested grids file of input passes through new, old two nested grids using input
Old mesh parameter is converted into actual coordinate (such as CDP coordinate) by the constant rule of the coordinate value of corresponding position in file, then
Coordinate value by new grid conversion at corresponding parameter, thus achieve the purpose that update grid.Mesh parameter is also known as speed net
Lattice control parameter.Wherein, grid file define 12 critical speed grid control parameters include: INLINE starting symbol,
INLINE termination number, INLINE increment (or being INLINE number spacing), XLINE starting symbol, XLINE starting symbol, XLINE increase
Measure (or being XLINE number spacing), grid starting point (or being grid origin) X-coordinate, grid starting point Y coordinate, INLINE net
Compartment away from the X direction projector distance, INLINE grid spacing projector distance in the Y direction, XLINE grid spacing is in X
Projector distance, the projector distance of XLINE grid spacing in the Y direction on direction.
The method of speed data grid updating according to an exemplary embodiment of the present invention the following steps are included:
S1, the old speed data of input, the old corresponding old mesh parameter of speed data and new mesh parameter.
S2, any No. CDP and old mesh parameter according to the trace header of old speed data, calculate any CDP correspondence
Face element coordinate.I.e. by No. CDP of the trace header of old speed data, according to old grid, calculates this No. CDP corresponding face element and sit
Mark, specifically includes the following steps:
S101, it is calculated by current corresponding No. CDP of old grid apart from grid origin (INLINE starting symbol, and XLINE starting
Number corresponding CDP point) interval INLINE line number and XLINE points, specifically:
(1) the point number nx of each INLINE line can be calculated by grid file:
Nx=1+ (XLINE terminates number-XLINE starting symbol)/XLINE increment.
(2) the INLINE number of lines steps before INLINE line where calculating current No. CDP:
Steps=int (the point number of (current No. CDP-XLINE starting symbol)/each INLINE line), int representative pair
Calculated result is rounded.
(3) the XLINE number of lines stepx before the XLINE line at current No. CDP place:
The INLINE line before INLINE line where stepx=(current No. CDP-XLINE starting symbol)-is No. CDP current
The point number of item number * each INLINE line.
S102, current No. CDP coordinate in old grid is calculated.
X-coordinate:
INLINE number of lines * INLINE grid before the INLINE line of x=grid starting point X-coordinate+current No. CDP place
XLINE number of lines * XLINE grid spacing before the projector distance+current No. CDP place XLINE line of spacing in the X direction
Projector distance in the X direction.
Y coordinate:
INLINE number of lines * INLINE grid before the INLINE line of y=grid starting point Y coordinate+current No. CDP place
XLINE number of lines * XLINE grid spacing before the projector distance+current No. CDP place XLINE line of spacing in the Y direction
Projector distance in the Y direction.
S3, the principle constant according to face element coordinate under the same coordinate system, according to the corresponding face element of any record Taoist monastic name
Coordinate and new mesh parameter calculate updated speed data trace header.
According to principle new, that old grid coordinate system is constant, i.e., identical coordinate is also identical in different grids.Pass through the seat
It marks with new grid and calculates other updated all speed data trace headers, as No. INLINE, No. XLINE, No. CDP, CDPX are sat
Mark and CDPY coordinate etc., specifically includes the following steps:
S301, according to changing coordinates and new No. INLINE and No. XLINE of new grid computing:
Dx=(the new grid starting point X-coordinate of current X-coordinate -).
Dy=(the new grid starting point Y coordinate of current Y coordinate -).
Dinv=1/ (the new grid XLINE grid spacing of the projector distance * of new grid INLINE grid spacing in the X direction
Between the new grid XLINE grid of projector distance * of the new grid INLINE grid spacing of projector distance-in the Y direction in the Y direction
Away from projector distance in the X direction).
Del=is (between the new grid XLINE grid of projector distance-dy* of the new grid XLINE grid spacing of dx* in the Y direction
Away from projector distance * dinv in the X direction).
Subline (new No. INLINE)=new grid INLINE starting symbol+new grid INLINE increment * del.
Del '=(int) (0.5+ (the new grid of projector distance-dx* of the new grid INLINE grid spacing of dy* in the X direction
The projector distance of INLINE grid spacing in the Y direction) * dinv).
Xline (new No. XLINE)=new grid XLINE starting symbol+new grid XLINE increment * del '.
S302, according to current new No. INLINE, No. XLINE No. CDP new with new grid computing.
Steps=(new No. INLINE-new grid INLINE starting symbol)/new grid INLINE increment.
Stepx=(new No. XLINE-new grid XLINE starting symbol)/new grid XLINE increment.
New No. CDP=steps* (1+ (new grid XLINE termination number-new grid XLINE starting symbol)/new grid XLINE
Increment)+stepx+1.
Refer to common-depth-point gather sequence according to CDP sequence in the above exemplary embodiments,.Before seismic data superposition
CDP is the trace gather in corresponding many same face elements, and the CDP of speed data only corresponds to one of speed data.But this hair
It is bright without being limited thereto, the trace gather of speed data, such as conventional CDP, CMP (the two meanings basic one can be generated according to the front
Cause) and CRP, CIP etc., grid updating is carried out, the sequence of their data numbers does not have essential distinction, only corresponds to trace gather
The influence for whether considering subsurface formations construction generated.That is method of the invention can also carry out CMP, CRP, CIP etc.
Processing.
Fig. 1 shows the flow chart of the speed data grid updating method of an exemplary embodiment of the present invention.Such as Fig. 1
Shown, in an exemplary embodiment of the present invention, the method for speed data grid updating can be used such as under type:
(1) N number of old speed data is inputted.Wherein, N is the positive integer more than or equal to 1.
(2) N number of grid file corresponding with speed data old in step (1) is inputted.Wherein, the numerical value of N and step (1)
The numerical value of middle N is equal.
(3) new speed grid file is inputted, old speed file is carried out according to the constant principle of coordinate according to CDP sequence
Grid trace header data calculate.
(4) judge whether all CDP are disposed, if "No", return step (3) continues, if "Yes", entrance is next
Step.That is, "No" expression has handled all CDP points of speed data not yet, continue to locate so still needing to return previous cycles
Reason；"Yes" expression has handled all CDP points of speed data, so continuing the processing of next link.
(5) currently after had been friends in the past speed file process, into carrying out new speed merging in next step.Currently had been friends in the past speed
After spending file process, into carrying out new speed merging in next step.
Old grid is some scrappy block of cells explored in the past, can have overlapping and blank area between block, newly
Grid is to cover these with regard to a big grid of grid, but after having updated respective speed data trace header in preceding step, is needed
Be merged into the speed data based on big new grid, so need to reject the block of overlapping, the block of blank into
Row interpolation.
(6) new speed trace header grid updating.
(7) new speed data is exported, calculating terminates.
The method of speed data grid updating of the invention is further described below in conjunction with specific example.
Example
The identical coordinate of the same coordinate system causes INLINE with XLINE reference numeral different due to the movement of grid, because
This, CDP number also changes.
The speed grid control parameter of old grid and new grid is as shown in following table, Fig. 2 and Fig. 3, wherein old grid origin is
O point in Fig. 1, new grid origin are the O ' point in Fig. 2.
| Speed grid control parameter | Old grid | New grid |
| INLINE starting symbol | 1 | 1 |
| INLINE termination number | 100 | 200 |
| INLINE increment | 1 | 1 |
| XLINE starting symbol | 101 | 101 |
| XLINE termination number | 2100 | 2100 |
| XLINE increment | 1 | 1 |
| Grid starting point X-coordinate | 0 | 100 |
| Grid starting point Y coordinate | 0 | 100 |
| The projector distance of INLINE grid spacing in the X direction | 0 | 0 |
| The projector distance of INLINE grid spacing in the Y direction | 10 | 10 |
| The projector distance of XLINE grid spacing in the X direction | 10 | 10 |
| The projector distance of XLINE grid spacing in the Y direction | 0 | 0 |
As shown in Fig. 2, handling one No. CDP=100151 in old grid:
Calculate point number nx=int (1+ (XLINE terminates number-XLINE starting symbol)/XLINE increasing of each INLINE line
Amount)=2000.
INLINE number of lines steps=(the current No. CDP-XLINE before INLINE line where calculating current No. CDP
Starting symbol)/each INLINE line point number=int ((100151-101)/2000)=50
Calculate it is described it is any record Taoist monastic name where cross-track before XLINE number of lines stepx=(current No. CDP-
XLINE starting symbol)-it is No. CDP current where INLINE line before INLINE number of lines * each INLINE line point number
=(100151-101) -50*2000=50.
I.e. as shown in Fig. 2, No. CDP=100051 in old grid, corresponding INLINE=51, XLINE=151, i.e. A point.
Calculate current No. CDP coordinate in old grid:
INLINE number of lines * INLINE grid before the INLINE line of x=grid starting point X-coordinate+current No. CDP place
XLINE number of lines * XLINE grid spacing before the projector distance+current No. CDP place XLINE line of spacing in the X direction
Projector distance=0+50*0+50*10=500 in the X direction.
INLINE number of lines * INLINE grid before the INLINE line of y=grid starting point Y coordinate+current No. CDP place
XLINE number of lines * XLINE grid spacing before the projector distance+current No. CDP place XLINE line of spacing in the Y direction
Projector distance=0+50*10+50*0=500 in the Y direction.
That is, current No. CDP (100151) coordinate in old grid are as follows: X=500 meters, Y=500 meters.
As shown in figure 3, according to changing coordinates and new No. INLINE and No. XLINE of new grid computing, wherein newly grid
Grid starting point is O ' (100,100):
Dx=(the new grid starting point X-coordinate of current X-coordinate -)=(500-100)=400.
Dy=(the new grid starting point Y coordinate of current Y coordinate -)=(500-100)=400.
Dinv=1/ (the new grid XLINE grid spacing of the projector distance * of new grid INLINE grid spacing in the X direction
Between the new grid XLINE grid of projector distance * of the new grid INLINE grid spacing of projector distance-in the Y direction in the Y direction
Away from projector distance in the X direction)=1/ (0*10-10*10)=- 0.01.
Del=is (between the new grid XLINE grid of projector distance-dy* of the new grid XLINE grid spacing of dx* in the Y direction
Away from projector distance * dinv in the X direction)=(400*0-400*10* (- 0.01))=40.
No. INLINE new=new grid INLINE starting symbol+new grid INLINE increment * del=1+1*40=41.
Del '=(int) (the new grid of projector distance-dx* of the new grid INLINE grid spacing of dy* in the X direction
The projector distance of INLINE grid spacing in the Y direction) * dinv)=int ((400*0-400*10) * (- 0.01))=40.
No. XLINE new=new grid XLINE starting symbol+new grid XLINE increment * del '=101+1*40=141.
According to currently new No. INLINE, No. XLINE No. CDP new with new grid computing:
Steps=(new No. INLINE-new grid INLINE starting symbol)/new grid INLINE increment=(41-1)/1=
40。
Stepx=(new No. XLINE-new grid XLINE starting symbol)/new grid XLINE increment=(141-101)/1=
40。
New No. CDP=steps* (1+ (new grid XLINE termination number-new grid XLINE starting symbol)/new grid XLINE
Increment)+stepx+1=40* (1+ (2100-101)/1)+40+1=80041.
That is, coordinate is X=500 meters in new grid, Y=500 meters of corresponding No. CDP=80041, INLINE=41, XLINE
=141, the A ' point in corresponding diagram 3.
It repeats the above steps, until all CDP are disposed.
Generally speaking, a kind of method of speed data grid updating provided by the invention is fixed by speed grid control parameter
Justice, using based under the same coordinate system, the constant thinking of the coordinate of face element ultimately forms speed data grid updating method.With it is existing
There is technology to compare, the invention has the following advantages that the many-one input pattern that (1) is flexible and convenient: the old grid and speed of input
File can be multipair, and be output under unified new mesh definition mode.(2) new data merges and interpolation automatically: for defeated
The old grid and speed file entered is after new grid conversion, and there are spatial sampling distances and time sampling apart from irregular existing
As automatic merging and interpolation form final unified speed file.Speed data grid updating method proposed by the present invention is special
Suitable for the seismic data process of earth's surface and subsurface picture complexity, have the characteristics that calculate step simplicity and high treating effect,
It has broad application prospects in seismic data high Precision Processing at present.
Although having been combined attached drawing above and exemplary embodiment describing the present invention, those of ordinary skill in the art
It will be apparent to the skilled artisan that can be carry out various modifications to above-described embodiment in the case where not departing from spirit and scope of the claims.
Claims (4)
1. a kind of method of speed data grid updating, which is characterized in that the described method comprises the following steps:
A, old speed data, the old corresponding old mesh parameter of speed data and new mesh parameter are inputted；
B, according to any record Taoist monastic name of the trace header of old speed data and old mesh parameter, it is corresponding to calculate any record Taoist monastic name
Face element coordinate；
C, the principle constant according to face element coordinate under the same coordinate system, according to the corresponding face element coordinate of any record Taoist monastic name and
New mesh parameter calculates updated speed data trace header；
D, repeating said steps B and C, until all record Taoist monastic names of the trace header of old speed data have been handled, wherein
The old mesh parameter and new mesh parameter include: main profile starting symbol, main profile termination number, main profile increment, contact
Survey line starting symbol, cross-track termination number, cross-track increment, grid starting point X-coordinate, grid starting point Y coordinate, main profile grid
Spacing projector distance in the X direction, main profile grid spacing projector distance in the Y direction, cross-track grid spacing are in X
The projector distance of projector distance and cross-track grid spacing in the Y direction on direction,
The step B includes:
B1, the point number nx for calculating each main profile,
Nx=1+ (cross-track termination number-cross-track starting symbol)/cross-track increment；
The main profile number of lines steps before main profile where b2, calculating any record Taoist monastic name,
Steps=int ((any record Taoist monastic name-cross-track starting symbol)/nx)；
The cross-track number of lines stepx before cross-track where b3, calculating any record Taoist monastic name,
Stepx=(any record Taoist monastic name-cross-track starting symbol)-steps*nx；
B4, any record Taoist monastic name corresponding face element coordinate (x, y) in old grid is calculated:
X-coordinate:
Projector distance+stepx* cross-track the net of x=grid starting point X-coordinate+steps* main profile grid spacing in the X direction
Compartment is away from projector distance in the X direction；
Y coordinate:
Projector distance+stepx* cross-track the net of y=grid starting point Y coordinate+steps* main profile grid spacing in the Y direction
Compartment away from projector distance in the Y direction,
The step C includes:
C1, updated main profile number and connection are calculated according to any corresponding face element coordinate of record Taoist monastic name and new mesh parameter
Network survey line number；
Dx=(the new grid starting point X-coordinate of current X-coordinate -)；
Dy=(the new grid starting point Y coordinate of current Y coordinate -)；
(the new grid cross-track grid spacing of the projector distance * of new grid main profile grid spacing in the X direction exists dinv=1/
The new grid cross-track grid of projector distance * of the new grid main profile grid spacing of projector distance-in the Y direction in Y-direction
The projector distance of spacing in the X direction)；
Del=(the new grid cross-track grid of projector distance-dy* of the new grid cross-track grid spacing of dx* in the Y direction
The projector distance * dinv of spacing in the X direction)；
New main profile number=new grid main profile starting symbol+new grid main profile increment * del；
Del '=int ((the new grid main profile grid of projector distance-dx* of the new grid main profile grid spacing of dy* in the X direction
The projector distance of spacing in the Y direction) * dinv)；
New cross-track number=new grid cross-track starting symbol+new grid cross-track increment * del '；
C2, the new main profile number being calculated according to step c1, new cross-track number calculate new recorded trace with new mesh parameter
Number；
Steps=(the new grid main profile starting symbol of new main profile number -)/new grid main profile increment；
Stepx=(the new grid cross-track starting symbol of new cross-track number -)/new grid cross-track increment；
New record Taoist monastic name=steps* (1+ (new grid cross-track termination number-new grid cross-track starting symbol)/new grid connection
Network survey line increment)+stepx+1.
2. the method for speed data grid updating according to claim 1, which is characterized in that after being updated in the step C
Speed data trace header include new main profile number, new cross-track number and new record Taoist monastic name.
3. the method for speed data grid updating according to claim 1, which is characterized in that the record Taoist monastic name is CDP
Number, No. CMP, No. CRP or No. CIP.
4. the method for speed data grid updating according to claim 1, which is characterized in that the method also includes in step
After rapid D, new speed merging is carried out, the merging step includes: to reject to overlapping block, is carried out to blank block slotting
Value, to be merged into the speed data based on new grid.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610786593.8A CN106383361B (en) | 2016-08-31 | 2016-08-31 | A kind of method of speed data grid updating |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610786593.8A CN106383361B (en) | 2016-08-31 | 2016-08-31 | A kind of method of speed data grid updating |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106383361A CN106383361A (en) | 2017-02-08 |
| CN106383361B true CN106383361B (en) | 2018-12-11 |
Family
ID=57939438
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610786593.8A Active CN106383361B (en) | 2016-08-31 | 2016-08-31 | A kind of method of speed data grid updating |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106383361B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110068863A (en) * | 2019-04-04 | 2019-07-30 | 中国科学院计算机网络信息中心 | A kind of determination method and device of common depth point CDP |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6201765B1 (en) * | 1995-11-27 | 2001-03-13 | Geco-Prakle (Uk) Limited | Method of monitoring quality of seismic data processing and method of processing vertical seismic profile data |
| CN102901985A (en) * | 2011-07-25 | 2013-01-30 | 中国石油化工股份有限公司 | Depth domain layer speed correcting method suitable for undulating surface |
| CN104570098A (en) * | 2015-01-28 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | Method and device for acquiring common-azimuth-angle seismic data |
| CN105158797A (en) * | 2015-10-16 | 2015-12-16 | 成都理工大学 | Staggered mesh wave equation forward modeling method based on actual seismic data |
| CN105549080A (en) * | 2016-01-20 | 2016-05-04 | 中国石油大学(华东) | Undulating ground surface waveform inversion method based on auxiliary coordinate system |
-
2016
- 2016-08-31 CN CN201610786593.8A patent/CN106383361B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6201765B1 (en) * | 1995-11-27 | 2001-03-13 | Geco-Prakle (Uk) Limited | Method of monitoring quality of seismic data processing and method of processing vertical seismic profile data |
| CN102901985A (en) * | 2011-07-25 | 2013-01-30 | 中国石油化工股份有限公司 | Depth domain layer speed correcting method suitable for undulating surface |
| CN104570098A (en) * | 2015-01-28 | 2015-04-29 | 中国石油集团东方地球物理勘探有限责任公司 | Method and device for acquiring common-azimuth-angle seismic data |
| CN105158797A (en) * | 2015-10-16 | 2015-12-16 | 成都理工大学 | Staggered mesh wave equation forward modeling method based on actual seismic data |
| CN105549080A (en) * | 2016-01-20 | 2016-05-04 | 中国石油大学(华东) | Undulating ground surface waveform inversion method based on auxiliary coordinate system |
Non-Patent Citations (1)
| Title |
|---|
| 三维地震网格参数提取方法研究与应用;张文书 等;《石油工业计算机应用》;20081231;第16卷(第3期);第6-7、24页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106383361A (en) | 2017-02-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106646598B (en) | A kind of FAST-AIC methods microseism picking up signal method | |
| CN102879820B (en) | Based on the three-dimensional table layer model construction method of triangle gridding | |
| CN105807316B (en) | Ground observation microseism velocity model corrections method based on amplitude superposition | |
| CN104036544A (en) | Building roof reconstruction method based on airborne LiDAR data | |
| CN105093319B (en) | Ground micro-seismic static correcting method based on 3D seismic data | |
| CN109738940B (en) | A method for locating acoustic emission/microseismic events in the presence of empty space | |
| CN102901985A (en) | Depth domain layer speed correcting method suitable for undulating surface | |
| CN106772577A (en) | Source inversion method based on microseism data and SPSA optimized algorithms | |
| CN106291725A (en) | A kind of method of fast inversion underground geologic bodies locus | |
| CN105954795A (en) | Grid successive dissection method used for microseismic positioning | |
| CN104360396B (en) | A kind of three kinds of preliminary wave Zoumaling tunnel methods of TTI medium between offshore well | |
| CN107678057A (en) | Determined in 3-D seismics or optimize wave detector or equipment investment amount and the method for arrangement | |
| CN113534259A (en) | Vibroseis efficient acquisition real-time prestack time migration imaging method | |
| CN106383361B (en) | A kind of method of speed data grid updating | |
| CN106054252B (en) | A kind of method and device of pre-stack time migration | |
| CN105717538B (en) | Undulating surface seismic data migration datum plane conversion method and device | |
| CN113960532B (en) | A microseismic location method based on secondary location calculation of hypothetical sources | |
| CN116660980A (en) | A Microseismic Location Method Based on Improved Equation | |
| CN105223630B (en) | Omnibearing observation systematic parameter Demonstration Method based on geological model | |
| CN104280766B (en) | A kind of direct offset method of utilization local data lineups slope | |
| RU2539097C2 (en) | Method of logging geomagnetic field variations on additional tie lines when performing magnetic survey on water bodies | |
| CN106291703B (en) | A kind of binary system processing method of seismic velocity data | |
| CN109188516A (en) | The micro-seismic event localization method of the domain Radon energy scan superposition | |
| CN105425290B (en) | A kind of method and device of pre-stack time migration | |
| CN103616720B (en) | A kind of refraction velocity acquisition methods and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan Applicant after: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180330 Address after: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei Applicant after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201110 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: BGP Inc., China National Petroleum Corp. Address before: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei Patentee before: BGP Inc., China National Petroleum Corp. |