CN105301639B - The method and its device of double weighting chromatography inversion speeds when being travelled based on VSP - Google Patents
The method and its device of double weighting chromatography inversion speeds when being travelled based on VSP Download PDFInfo
- Publication number
- CN105301639B CN105301639B CN201510685710.7A CN201510685710A CN105301639B CN 105301639 B CN105301639 B CN 105301639B CN 201510685710 A CN201510685710 A CN 201510685710A CN 105301639 B CN105301639 B CN 105301639B
- Authority
- CN
- China
- Prior art keywords
- mrow
- region
- grid
- vsp
- msub
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses it is a kind of travelled based on VSP when double weighting chromatography inversion speeds method and devices, be related to technical field of geophysical exploration, it can improve the precision and reliability of tomographic inversion result.This method includes：The initial velocity model that VSP geological datas obtain region is gathered in region；Each position in region sets up the credibility model in region according to its distance away from well and its depth away from earth's surface；Observed parameter, the initial velocity model in region and credibility model during based on collection VSP geological datas obtain grid slowness disturbed value in region；The new rate pattern for obtaining region is modified to the speed of grid in the initial velocity model in region by grid slowness disturbed value in region；When obtaining first arrival theory based on new rate pattern and walking, whether difference when judging when the first arrival theory is walked and actually walking meets predetermined threshold values, if meeting, output chromatography inversion speed.This method effectively increases the precision and reliability of tomographic inversion result.
Description
Technical field
The present invention relates to technical field of geophysical exploration, double weighting chromatographies are anti-during more particularly to a kind of travelling based on VSP
Drill the method and its device of velocity field.
Background technology
So far, China's western region pre-stack depth migration still not yet obtains remarkable effect, its essential reason be by
In the velocity field for being difficult to determine underground medium.Under this kind of area, velocity field is difficult to determination and there are factors, including：
In China's western region, arid massif low velocity layer thickness is big, reaches 500-1200 meters；The common side such as little refraction and micro logging
Method can only the shallower part in inverting area velocity field, the velocity field in the regional deep of inverting is unable to completely；Top layer is asked using deep-well is beaten
The method cost of velocity field is excessive, and the velocity field so obtained after operation is insufficient to accurately；It can be run into common chromatography method
After ray goes down, return less than ground, can not form the situation of diving Wave.Many factors are finally resulted in and can not met above
It is required that velocity field.
Due to there are a variety of unfavorable factors, conventional method can not all prove effective, and this area researcher begins attempt to new
Method carrys out the velocity field that inverting obtains underground medium, including based on vertical seismic profiling (VSP) VSP (Vertical Seismic
Profiling) method of first arrival time chromatography carrys out the velocity field of inverting underground medium.Tomography is a kind of using in detection mesh
The mark signal that observes of surface asks for the inversion method of target internal information, earliest applied to medical domain, then expands to
Other necks.In seismic exploration problem, tomography is exactly the travelling using the seismic wave that recorded in earth's surface or downhole receiving
When, the information such as amplitude and waveform, the distribution of the parameter such as the speed to rebuild underground medium.At present, in field of seismic exploration application
Widest chromatography imaging method is ray travel time tomography, that is, utilizes seimic travel time and the ray road of its propagation
Inverting underground medium speed is carried out in footpath.Conventional travel time tomography method has Inverse Projection, algebraic reconstruction technique, Joint iteration weight
Build method, singular value decomposition method, least square QR decomposition methods and least-squared conjugate gradient method etc..Inventors herein have recognized that on
These algorithms are stated all in accordance with the ray length through unit grid to distribute walking time error, walking time error is to fire to point out from ground
Hair looks for ray to reach VSP undergrounds point, according to the speed for iterating velocity field, when calculating its travelling and walking, then with real VSP
The error time that direct wave arrival time tries to achieve, may finally accurately be obtained although distributing walking time error according to ray length
The underground medium speed of underground distant place is obtained, but is distributed according to ray length by the underground well that walking time error can cause acquisition on the contrary
Underground medium speed it is more and more inaccurate, in this way, the inaccuracy of the underground medium speed by the well of underground can not meet demand,
And then cause the precision and reliability of the result of medium velocity progress tomographic inversion under base area not enough.
The content of the invention
Add in order to double when overcoming in the drawbacks described above of prior art, the embodiment of the present invention there is provided a kind of travelling based on VSP
The method and device of tomographic inversion velocity field is weighed, it can improve the precision and reliability of tomographic inversion result.
The concrete technical scheme of the embodiment of the present invention is：
A kind of methods of weighting chromatography inversion speeds double when being travelled based on VSP, it includes：
Vertical seismic profiling (VSP) VSP geological datas are gathered in region so as to obtain the first arrival time of VSP earthquake records, according to
The first arrival time of the VSP earthquake records obtains the initial velocity model in the region；
The region is set up in each position in the region according to its distance away from well and its depth away from earth's surface
Credibility model；
Observed parameter, the initial velocity model in the region and the confidence level mould during based on collection VSP geological datas
Type obtains grid slowness disturbed value in the region；
The speed of grid in the initial velocity model in the region is carried out by grid slowness disturbed value in the region
Amendment obtains the new rate pattern in the region；
When obtaining first arrival theory based on the new rate pattern and walking, judge when the first arrival theory is walked with when actually walking it
Whether difference meets predetermined threshold values, if meeting, output chromatography inversion speed.
A kind of devices of weighting chromatography inversion speeds double when being travelled based on VSP, it includes：
Initial velocity model sets up module, and it sets up institute according to the vertical seismic profiling (VSP) VSP geological datas gathered in region
State the initial velocity model in region；
Credibility model sets up module, its according to distance of each position in the region away from well and its away from earth's surface
Depth sets up the credibility model in the region；
Grid slowness disturbed value computing module, it is according to the initial velocity model, credibility model and collection in the region
Observed parameter during VSP geological datas obtains grid slowness disturbed value in the region；
Tomographic inversion velocity field module, it passes through initial velocity of the grid slowness disturbed value in the region to the region
The speed of grid is modified the new rate pattern for obtaining the region in model, and then calculates when obtaining first arrival theory and walking,
As defined in difference satisfaction when walking when the first arrival theory is walked and actually during threshold values, output chromatography inversion speed, when described first
Difference when when being walked to theory and actually walking be unsatisfactory for as defined in threshold values when, then based on observed parameter, the new speed in the region
Model and the credibility model retrieve grid slowness disturbed value in the region, and then are iterated, until described first
Difference when when being walked to theory and actually walking meet as defined in threshold values.
The method and device of double weighting chromatography inversion speeds when being travelled based on VSP in the embodiment of the present invention, it is being walked
When mistiming difference is with calculating, in addition to the ray length of unit grid, also introduce credibility model, range site grid is penetrated
Line length and confidence level correct grid slowness disturbed value, control the distribution of walking time error, so as to improve tomographic inversion result
Reliability and precision, improve topographic imaging by velocity result.
Brief description of the drawings
Accompanying drawing described here is only used for task of explanation, and is not intended in any way limit model disclosed by the invention
Enclose.In addition, shape and proportional sizes of each part in figure etc. are only schematical, the understanding of the present invention is used to help, and
It is not the shape and proportional sizes for specifically limiting each part of the invention.Those skilled in the art under the teachings of the present invention, can
To select various possible shapes and proportional sizes to implement the present invention as the case may be.
Flow charts of the Fig. 1 for the present invention in embodiment.
Fig. 2 is the geological model and VSP inspection wells of forward simulation use in the embodiment of the present invention.
Fig. 3 is forward simulation observation system schematic diagram in the embodiment of the present invention.
Fig. 4 is the initial velocity model set up in the embodiment of the present invention.
Fig. 5 is the credibility model set up in the embodiment of the present invention.
The velocity field that Fig. 6 obtains for tomographic inversion in the embodiment of the present invention.
Fig. 7 is the structural representation of device in the embodiment of the present invention.
Embodiment
With reference to the description of accompanying drawing and the specific embodiment of the invention, the details of the present invention can be clearly understood.But
It is, the embodiment of invention described herein to be only used for explaining the purpose of the present invention, and can not understands in any way
Into being limitation of the present invention.Under the teachings of the present invention, technical staff is contemplated that any possible change based on the present invention
Shape, these are regarded as belonging to the scope of the present invention.
Flow charts of the Fig. 1 for the present invention in embodiment, as shown in figure 1, double weighting chromatographies are anti-during a kind of travelling based on VSP
The method for drilling velocity field, it comprises the following steps：
S101：When first arrival of the vertical seismic profiling (VSP) VSP geological datas so as to obtain VSP earthquake records is gathered in region
Between, the initial velocity model in region is obtained according to the first arrival time of VSP earthquake records, it comprises the following steps：
S201：Deep-well vertical seismic profiling (VSP) VSP earthquake numbers are gathered by the method for ray-tracing forward modeling in region
According to, utilize first arrival time pickup software pick up VSP earthquake records first arrival time.Fig. 2 is forward simulation in the embodiment of the present invention
The geological model and VSP inspection wells of use, as shown in Fig. 2 Fig. 2 middle and upper parts are derrick, vertical black line is inspection well in Fig. 2.
3 to 5 mouthfuls of wells are looked in region, different depth will be arrived under multi-level downhole wave detector, GPS time services are disposed close to the borehole wall, and in well head
Seismic detector；Then using the method for putting small big gun, 0.5kg to 1kg dose can be used, according to the complicated journey of regional tectonics
Degree, can be blown out with GPS time service blasters on mountain；Fig. 3 is forward simulation observation system signal in the embodiment of the present invention
Figure, as shown in figure 3, big gun falls away from 25 meters away from for 100 meters, the scale blown out is the roads of 101 big gun * 215.It is complete using GPS time service seismic detectors
Weather receives VSP earthquake records.When the VSP earthquake records received according to GPS time services seismic detector pick up the first arrival of zero-offset VSP
Between, detailed process is as follows, picks up software using first arrival time, pickup zero-offset VSP is per first arrival time TR (I), I=together
1 ... N, wherein, I is Seismic Traces sequence number, and N is the total road number of VSP earthquake records.
S202：Initial construction model is set up using geological structure information known to region.
The geological structure information according to known to region sets up initial construction model, it is known that geological structure information at least include
One of the depth of subsurface reflective boundary, the form of reflecting interface, position of fault and fault strike.
S203：First arrival time inverting based on VSP earthquake records obtains the interval velocity in well week, by interval velocity to initial structure
The progress of modeling type handles the initial model speed for obtaining grid.
First arrival time based on pickup zero-offset VSP carries out the interval velocity that inverting obtains well week, specifically in zero well
Source is away under the conditions of, the time difference that the first arrival for obtaining different depth according to the first arrival time of zero-offset VSP is reached, then calculates each depth
Interval velocity on degree, has also just obtained the interval velocity in well week.Then initial construction model is handled by interval velocity
The initial model speed of grid, Fig. 4 is the initial velocity model set up in the embodiment of the present invention.For well in initial construction model
Model meshes near footpath, are filled using the interval velocity of zero well constraint inversion of the well, are somebody's turn to do for remote in initial construction model
The model meshes of well, according to the distance of the grid to neighbouring two mouthfuls of wells, to the layer speed of zero well constraint inversion of two mouthfuls of wells near the grid
Degree carries out distance weighted calculating, obtains the initial model speed of the grid.In the mistake of the above-mentioned initial model speed for obtaining grid
Cheng Zhong, the data for relying primarily on VSP wells are filled processing to obtain the first of grid from the centre of initial construction model to both sides
Beginning model velocity.
S204：Initial model speed based on grid carries out layered encryption and Mesh Smoothing processing obtains the original net in region
Lattice rate pattern.Initial model speed to grid carries out layered encryption and Mesh Smoothing processing obtains initial mesh rate pattern
V (ix, jz), wherein ix represent the grid sequence number of rate pattern X-direction, and ix=1 ... NX, NX represent the net of model in the X direction
Lattice number, jz represents the grid sequence number of rate pattern Z-direction, and jz=1 ... NZ, NZ represent the grid number of model in z-direction.
S102：Each position in region according to its distance away from well and its depth away from earth's surface set up region can
Credit model.
Fig. 5 is the credibility model set up in the embodiment of the present invention, as shown in figure 5, in region by well and one in earth's surface
The individual clear position of top layer speed gives in a confidence value, region other positions with the distance away from well and apart from earth's surface
Depth is incremented by, and confidence value is gradually reduced, so that the credibility model N (x, z) of forming region, wherein x are rate pattern X side
To grid sequence number, x=1 ... NX, NX is the grid number of model in the X direction, and z is the grid sequence number of rate pattern Z-direction, z
=1 ... NZ, NZ are the grid number of model in z-direction.
S103：Observed parameter, the initial velocity model in region and credibility model during based on collection VSP geological datas are obtained
Grid slowness disturbed value in region, it specifically includes following steps：
S301：The theory of the VSP synthetic seismograms together with is obtained according to observed parameter and the initial velocity model in region
First arrival time, the propagation path of every one of ray and ray through unit grid in length.
Following parameter is obtained by the method for ray tracing according to the initial velocity model and observed parameter in region：Per together
The theoretical first arrival time Ts (i) of VSP synthetic seismograms, wherein i=1 ... N, i is Seismic Traces sequence number, and N remembers for VSP earthquakes
Lu Zong roads number；Per ray propagation path (Rx (j), Ry (j)) together, wherein, j=1 ... M, M are the grid that current ray is passed through
Sum, Rx (j) and Ry (j) are respectively the X-direction for j-th of grid that ray is passed through and the grid sequence number of Z-direction；Ray through
Unit grid in length Lij, wherein, i be ray sequence number, j be model unit grid sequence number.Observed parameter is specially
Involved by step S201 Fig. 3：Big gun falls away from 25 meters away from for 100 meters, and the scale blown out is the roads of 101 big gun * 215.
S302：Theoretical first arrival time based on every one of VSP synthetic seismograms, the propagation path of every one of ray, penetrate
Line through unit grid in length and credibility model obtain grid slowness disturbed value in region.
Theoretical first arrival time based on every one of VSP synthetic seismograms, the propagation path of every one of ray, ray are being worn
Length and credibility model in the unit grid crossed are obtained in region in grid slowness disturbed value, are specially for there is J net
The region of lattice, during the situation of I bar ray informations, the calculation formula of grid slowness disturbed value is as follows in its region：
Wherein, LijRepresent ray length of i-th ray on j-th of grid, NijRepresent i-th ray in j-th of net
Confidence level on lattice, SLiRepresent the total length of i-th ray, SNiTotal confidence level of i-th ray is represented, η and α represent control
The relaxation factor of convergence rate and stability processed, I represents total number of rays, and J represents total grid number.
Because the azimuthal distribution that walking time error distributes not only relevant with ray length ray also and in unit grid has
Close, simple ray length loses control of walking time error distribution, so walking time error distribution is carried out in order to more accurate, in region
Interior grid slowness disturbed value is introduced into the confidence level of unit grid in credibility model when calculating.
In above-mentioned zone in the calculation formula of grid slowness disturbed value ray order be it is ascending according to well constraint inversion,
The order of receiving point depth from shallow to deep is ranked up.η's in the calculation formula of the grid slowness disturbed value of above-mentioned model takes
Value scope can be (0,2), and α span can be (0,2).
S104：The speed of grid in the initial velocity model in region is modified by the grid slowness disturbed value of model
Obtain the new rate pattern in region.
The speed of corresponding grid in the initial velocity model in region is repaiied using grid slowness disturbed value in region
The new rate pattern in region is just being obtained, speed correction formula is specially：
Wherein, ix represents the grid sequence number of the initial velocity model X-direction in region, ix=1 ..., NX, NX represents region
Grid number in initial velocity model X-direction, jz represents the grid sequence number of the initial velocity model Z-direction in region, jz=1 ...,
NZ, NZ represent the grid number in the initial velocity model Z-direction in region, Δ SjRepresent grid slowness disturbed value in region.
In the speed modification of grid in initial model, the reasonability of amended speed can be judged, it is right
Speed correction beyond given speed scope is defined.During confidence level revised at the same time it can also speed to grid
Row judge, if amended confidence level is too low, for example through current grid ray very little, can be generally considered as be less than 5
When, then the speed of current grid is not corrected.
S105：The velocity field in the new rate pattern in region is smoothed by linear filtering method.
When new rate pattern uses Kp × Lp smooth grid, then the formula in linear filtering method is specially：
Wherein, m, n represent the abscissa and ordinate at grid midpoint respectively, and V (m, n) represents that grid midpoint is (m, n)
Grid Velocity.The step can be with the exceptional value in release rate, to ensure the reliability of ray tracing result.
S106：When obtaining first arrival theory based on new rate pattern and walking, difference when judging when first arrival theory is walked and actually walking
Whether defined threshold values is met, if meeting, output chromatography inversion speed.Based on observed parameter, region if being unsatisfactory for
New rate pattern and credibility model retrieves the grid slowness disturbed value of model, is thus iterated, until first famous dictum
By difference when when walking and actually walking meet as defined in threshold values.
The velocity field that Fig. 6 obtains for tomographic inversion in the embodiment of the present invention, as shown in fig. 6, using new rate pattern just
Drill and calculate when obtaining the theories of VSP first arrivals and walking, judge when first arrival theory is walked with the VSP geological datas that gather when actually walking
Whether the difference of average value meets predetermined threshold values, if meeting predetermined threshold values, output chromatography inversion speed；If being unsatisfactory for
Predetermined threshold values, then repeat step S103 to S105 is to be iterated, until VSP geological datas when first arrival theory is walked with collection
In the difference of average value when actually walking meet predetermined threshold values, then export chromatography inversion speed.
In the embodiment of the present invention according to VSP preliminary wave radiation profiles it is extremely uneven the characteristics of, in tomographic inversion calculating process
In, except correcting grid slowness according to the ray length in unit grid, also introduce credibility model and come to grid matrix
Different confidence levels are set, and it is slow to correct grid using the ray length and double weighted formulas of the confidence level of grid in grid
Spend disturbed value.Ascending according to well constraint inversion in grid slowness amendment, ray is carried out the order of depth from shallow to deep one by one,
The reliability and precision of tomographic inversion result can so be improved.Two parameters of ray length and confidence level are utilized by this method
To control the distribution of walking time error, the reliability and precision of tomographic inversion result can be effectively improved, and then improve speed chromatography
Imaging results.
Fig. 7 is the structural representation of device in the embodiment of the present invention, as shown in fig. 7, double weightings during a kind of travelling based on VSP
The device of tomographic inversion velocity field, it includes：
Initial velocity model sets up module, and it sets up area according to the vertical seismic profiling (VSP) VSP geological datas gathered in region
The initial velocity model in domain；
Credibility model sets up module, and it is according to distance of each position in region away from well and its depth away from earth's surface
Set up the credibility model in region；
The grid slowness disturbed value computing module of model, it is according to the initial velocity model in region, credibility model and adopts
Observed parameter during collection VSP geological datas obtains grid slowness disturbed value in region；
Tomographic inversion velocity field module, it is by grid slowness disturbed value in region to net in the initial velocity model in region
The speed of lattice is modified the new rate pattern for obtaining region, and then calculates when obtaining first arrival theory and walking, and is walked originally to theory
When and difference when actually walking meet as defined in threshold values when, output chromatography inversion speed, when being walked originally to theory with when actually walking
Difference be unsatisfactory for as defined in threshold values when, then retrieved based on observed parameter, the new rate pattern in region and credibility model
Grid slowness disturbed value in region, and then be iterated, until difference when first arrival theory is walked and when actually walking meets defined valve
Value.
System, device, module or unit that above-described embodiment is illustrated, can specifically be realized by computer chip or entity,
Or realized by the product with certain function.
For convenience of description, it is divided into various units during description apparatus above with function to describe respectively.Certainly, this is being implemented
The function of each unit can be realized in same or multiple softwares and/or hardware during application.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.It is real especially for system
Apply for example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Remote processing devices connected by communication network perform task.In a distributed computing environment, program module can be with
Positioned at including in the local and remote computer-readable storage medium including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's
Spirit.
Claims (12)
1. a kind of methods of weighting chromatography inversion speeds double when being travelled based on VSP, it is characterised in that it includes：
Vertical seismic profiling (VSP) VSP geological datas are gathered in region so as to obtain the first arrival time of VSP earthquake records, according to described
The first arrival time of VSP earthquake records obtains the initial velocity model in the region；
Each position in the region according to its distance away from well and its depth away from earth's surface set up the region can
Credit model；
Observed parameter, the initial velocity model in the region and described during based on collection vertical seismic profiling (VSP) VSP geological datas
Credibility model obtains grid slowness disturbed value in the region；
The speed of grid in the initial velocity model in the region is modified by grid slowness disturbed value in the region
Obtain the new rate pattern in the region；
When obtaining first arrival theory based on the new rate pattern and walking, difference when judging when the first arrival theory is walked and actually walking is
It is no to meet predetermined threshold values, if meeting, output chromatography inversion speed.
2. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 1, it is characterised in that
When obtaining first arrival theory based on the new rate pattern and walking, whether difference when judging when the first arrival theory is walked and actually walking is full
The predetermined threshold values of foot, if it is not satisfied, then based on observed parameter, the new rate pattern in the region and the credibility model weight
The new grid slowness disturbed value obtained in region, and then loop iteration is carried out, until when the first arrival theory is walked and when actually walking
Difference meet predetermined threshold values.
3. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 1, it is characterised in that
Observed parameter, the initial velocity model in the region and institute when step is based on collection vertical seismic profiling (VSP) VSP geological datas
Credibility model is stated to obtain in the region in grid slowness disturbed value, specifically including following steps：
When obtaining the theoretical first arrival of VSP synthetic seismograms together with according to observed parameter and the initial velocity model in the region
Between, the propagation path of every one of ray and ray through unit grid in length；
Theoretical first arrival time based on every one of VSP synthetic seismograms, the propagation path of every one of ray, ray through
Length and the credibility model in unit grid obtain grid slowness disturbed value in the region.
4. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 3, it is characterised in that
Theoretical first arrival time of the step based on every one of VSP synthetic seismograms, the propagation path of every one of ray, ray through
Unit grid in length and the credibility model obtain in the grid slowness disturbed value of model, be specially for there is J
The region of grid, during the situation of I bar ray informations, the calculation formula of grid slowness disturbed value is as follows in its described region：
<mrow>
<msub>
<mi>&Delta;S</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>&eta;</mi>
<msub>
<mi>&lambda;</mi>
<mi>j</mi>
</msub>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mfrac>
<mrow>
<msub>
<mi>dT</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>SL</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>SN</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>SLN</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msubsup>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
<mi>&alpha;</mi>
</msubsup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mi>J</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>SL</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>SN</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>SLN</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>j</mi>
</munderover>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>SL</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>SN</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>...</mo>
<mi>I</mi>
<mo>)</mo>
</mrow>
</mrow>
1
Wherein, LijRepresent ray length of i-th ray on j-th of grid, NijRepresent i-th ray on j-th of grid
Confidence level, SLiRepresent the total length of i-th ray, SNiTotal confidence level of i-th ray is represented, η represents control convergence speed
Degree and the relaxation factor of stability, α represent the relaxation factor of control convergence speed and stability, and I represents total number of rays, and J is represented
Total grid number, error when dTi is represented when i-th ray theory forward modeling is walked and actually walked.
5. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 4, it is characterised in that
Ray order in the region in the calculation formula of grid slowness disturbed value is ascending according to its well constraint inversion, and receiving point is deep
The order of degree from shallow to deep is ranked up.
6. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 4, it is characterised in that
The span of η in the calculation formula of the grid slowness disturbed value of the model is (0,2), and α span is (0,2).
7. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 1, it is characterised in that：
The speed of grid is repaiied in initial velocity model in step by grid slowness disturbed value in the region to the region
In the new rate pattern for just obtaining the region, specially using grid slowness disturbed value in the region to the region
The speed of corresponding grid is modified the new rate pattern for obtaining the region, speed correction formula in initial velocity model
Specially：
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>x</mi>
<mo>,</mo>
<mi>j</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mi>x</mi>
<mo>,</mo>
<mi>j</mi>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>&Delta;S</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
</mrow>
Wherein, ix represents the grid sequence number of the initial velocity model X-direction in the region, ix=1 ..., NX, NX represents the area
Grid number in the initial velocity model X-direction in domain, jz represents the grid sequence number of the initial velocity model Z-direction in the region,
Jz=1 ..., NZ, NZ represents the grid number in the initial velocity model Z-direction in the region, Δ SjRepresent the region Intranet
Lattice slowness disturbed value.
8. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 1, it is characterised in that
Vertical seismic profiling (VSP) VSP geological datas are gathered in region in step so as to obtain the first arrival time of VSP earthquake records, according to
The first arrival time of VSP earthquake records is obtained in the initial velocity model in the region, is specifically included：
Vertical seismic profiling (VSP) VSP geological datas are gathered by the method for ray-tracing forward modeling in region, during using first arrival
Between pickup software pick up VSP earthquake records first arrival time；
Initial construction model is set up using geological structure information known to the region；
First arrival time inverting based on VSP earthquake records obtains the interval velocity in well week, and initial construction model is entered by interval velocity
Row processing obtains the initial model speed of grid；
Initial model speed based on grid carries out layered encryption and Mesh Smoothing processing obtains the initial velocity mould in the region
Type.
9. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 1, it is characterised in that
Before step obtains when first arrival theory is walked based on new rate pattern, methods described also includes：By linear filtering method to institute
The velocity field stated in the new rate pattern in region is smoothed.
10. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 8, it is characterised in that
When the new rate pattern uses Kp × Lp smooth grid, then the formula in linear filtering method is specially：
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>K</mi>
<mi>p</mi>
<mo>&times;</mo>
<mi>L</mi>
<mi>p</mi>
</mrow>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>K</mi>
<mi>p</mi>
</mrow>
</munderover>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>L</mi>
<mi>p</mi>
</mrow>
</munderover>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>+</mo>
<mi>i</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
</mrow>
Wherein, m, n represent the abscissa and ordinate at grid midpoint respectively, and V (m, n) represents the grid that grid midpoint is (m, n)
Speed, i represents smooth grid abscissa sequence number, and j represents smooth grid ordinate sequence number.
11. the methods of weighting chromatography inversion speeds double when being travelled based on VSP according to claim 1, it is characterised in that：
The region is set up according to its distance away from well and its depth away from earth's surface in each position of step in the region
In credibility model, it is specially：Given one of a clear position of top layer speed can by well and in earth's surface in the region
Other positions being incremented by with the distance away from well and the depth away from earth's surface in certainty value, the region, confidence value is gradually reduced,
So as to form the credibility model N (x, z) in the region, wherein x represents the grid sequence number of credibility model X-direction, and z is represented can
The grid sequence number of credit model Z-direction.
12. a kind of devices of weighting chromatography inversion speeds double when being travelled based on VSP, it is characterised in that it includes：
Initial velocity model sets up module, and it sets up the area according to the vertical seismic profiling (VSP) VSP geological datas gathered in region
The initial velocity model in domain；
Credibility model sets up module, and it is according to distance of each position in the region away from well and its depth away from earth's surface
Set up the credibility model in the region；
Grid slowness disturbed value computing module, it is vertical according to the initial velocity model in the region, credibility model and collection
Observed parameter during seismic profile VSP geological datas obtains grid slowness disturbed value in the region；
Tomographic inversion velocity field module, it passes through initial velocity model of the grid slowness disturbed value in the region to the region
The speed of middle grid is modified the new rate pattern for obtaining the region, and then calculates when obtaining first arrival theory and walking, and works as institute
When difference when stating when first arrival theory is walked and actually walking meets predetermined threshold values, output chromatography inversion speed, when the just famous dictum
When being unsatisfactory for predetermined threshold values by difference when when walking and actually walking, then the new rate pattern based on observed parameter, the region
Grid slowness disturbed value in the region is retrieved with the credibility model, and then is iterated, until the just famous dictum
Predetermined threshold values is met by difference when when walking and actually walking.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510685710.7A CN105301639B (en) | 2015-10-21 | 2015-10-21 | The method and its device of double weighting chromatography inversion speeds when being travelled based on VSP |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510685710.7A CN105301639B (en) | 2015-10-21 | 2015-10-21 | The method and its device of double weighting chromatography inversion speeds when being travelled based on VSP |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105301639A CN105301639A (en) | 2016-02-03 |
| CN105301639B true CN105301639B (en) | 2017-09-01 |
Family
ID=55199115
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510685710.7A Active CN105301639B (en) | 2015-10-21 | 2015-10-21 | The method and its device of double weighting chromatography inversion speeds when being travelled based on VSP |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105301639B (en) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107942372A (en) * | 2016-10-12 | 2018-04-20 | 中国石油化工股份有限公司 | Between well and well combine seismic CT imaging method and device |
| CN109655884B (en) * | 2017-10-10 | 2021-01-12 | 中国石油化工股份有限公司 | Method and system for removing speed pickup abnormal points |
| CN109884710B (en) * | 2019-03-20 | 2021-02-26 | 中国石油化工股份有限公司 | Micro-logging tomography method aiming at excitation well depth design |
| CN110058306B (en) * | 2019-04-11 | 2020-10-09 | 中国石油天然气股份有限公司 | Method, apparatus and computer storage medium for determining three-dimensional velocity volume |
| CN110162981B (en) * | 2019-04-18 | 2020-10-02 | 阿里巴巴集团控股有限公司 | Data processing method and device |
| US10846413B2 (en) | 2019-04-18 | 2020-11-24 | Advanced New Technologies Co., Ltd. | Data processing method and device |
| CN112305595B (en) * | 2019-07-24 | 2024-05-17 | 中国石油化工股份有限公司 | Method for analyzing geologic body structure based on refraction wave and storage medium |
| CN112415580B (en) * | 2019-08-21 | 2024-01-23 | 中国石油化工股份有限公司 | Method for eliminating abrupt interface of velocity model and prestack depth migration processing method |
| CN113589375B (en) * | 2020-04-30 | 2023-06-30 | 中国石油化工股份有限公司 | VSP layer speed inversion method based on calculation during constraint travel of inclined layer |
| CN112036008B (en) * | 2020-07-28 | 2024-06-25 | 中国石油天然气集团有限公司 | Method and system for recording inversion speed by using hybrid model |
| CN114814949B (en) * | 2021-01-21 | 2023-09-01 | 中国石油化工股份有限公司 | Shallow reverse VSP first arrival chromatography and stratum prediction method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101630016A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | Method for improving imaging quality of vertical seismic profile |
| GB2479990A (en) * | 2010-04-30 | 2011-11-02 | Schlumberger Holdings | Multicomponent seismic inversion of VSP data |
| CN103217715A (en) * | 2012-01-18 | 2013-07-24 | 中国石油天然气集团公司 | Multi-scale regular grid tomography inversion statics correction method |
| CN104459782A (en) * | 2014-12-10 | 2015-03-25 | 中国石油天然气股份有限公司 | Seismic velocity modeling method and modeling unit using thin layer chromatography inversion |
| CN104570106A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Near-surface tomographic velocity analysis method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2868167B1 (en) * | 2004-03-23 | 2006-05-19 | Inst Francais Du Petrole | METHOD FOR IMAGING HIGHLY SLOPED GEOLOGICAL INTERFACES IN SUBTERRANEAN FORMATION, PROVIDING FOR PRISMATIC REFLEXIONS |
| US7751279B2 (en) * | 2006-05-03 | 2010-07-06 | Baker Hughes Incorporated | Sub-salt reflection tomography and imaging by walkaway VSP survey |
-
2015
- 2015-10-21 CN CN201510685710.7A patent/CN105301639B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101630016A (en) * | 2008-07-16 | 2010-01-20 | 中国石油天然气集团公司 | Method for improving imaging quality of vertical seismic profile |
| GB2479990A (en) * | 2010-04-30 | 2011-11-02 | Schlumberger Holdings | Multicomponent seismic inversion of VSP data |
| CN103217715A (en) * | 2012-01-18 | 2013-07-24 | 中国石油天然气集团公司 | Multi-scale regular grid tomography inversion statics correction method |
| CN104570106A (en) * | 2013-10-29 | 2015-04-29 | 中国石油化工股份有限公司 | Near-surface tomographic velocity analysis method |
| CN104459782A (en) * | 2014-12-10 | 2015-03-25 | 中国石油天然气股份有限公司 | Seismic velocity modeling method and modeling unit using thin layer chromatography inversion |
Non-Patent Citations (3)
| Title |
|---|
| Variable damping constraint tomography and its application in VSP Data;Huang Guang-Nan et al.;《APPLIED GEOPHYSICS》;20120630;第9卷(第2期);第177-185页 * |
| 井间地震层析成像的应用技术;黄真萍等;《CT和三维成像学术年会论文集》;20041031;第106-109页 * |
| 变网格间距速度层析成像方法;黄光南等;《石油地球物理勘探》;20130630;第48卷(第3期);第379-389页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105301639A (en) | 2016-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105301639B (en) | The method and its device of double weighting chromatography inversion speeds when being travelled based on VSP | |
| CN104133245B (en) | The static correcting method and system of a kind of seismic data | |
| CN108064348A (en) | Seismic travel time tomography inversion method based on two-point ray tracing | |
| CN106772577B (en) | Source inversion method based on microseism data and SPSA optimization algorithm | |
| CN105386756B (en) | A method of brittle formation porosity is calculated using dependent variable | |
| CN102879820B (en) | Based on the three-dimensional table layer model construction method of triangle gridding | |
| CN106814391A (en) | Ground micro-seismic state event location method based on Fresnel zone tomographic inversion | |
| CN104459782B (en) | Seismic velocity modeling method and modeling unit using thin layer chromatography inversion | |
| CN103869363B (en) | Microseism localization method and device | |
| CN102901985B (en) | A kind of Depth Domain interval velocity modification method being applicable to relief surface | |
| CN103869368A (en) | Cannon first-arrival comprehensive modeling static correction method without surface layer survey data constraint | |
| CN105093319B (en) | Ground micro-seismic static correcting method based on 3D seismic data | |
| CN103088803B (en) | Method for identifying collapse factors caused by karst collapse | |
| CN107703540B (en) | A kind of microseism positioning and chromatography imaging method | |
| Küperkoch et al. | Long‐term monitoring of induced seismicity at the Insheim geothermal site, Germany | |
| Jolie et al. | The development of a 3D structural–geological model as part of the geothermal exploration strategy—a case study from the Brady’s geothermal system, Nevada, USA | |
| CN104360396B (en) | A kind of three kinds of preliminary wave Zoumaling tunnel methods of TTI medium between offshore well | |
| CN109884710A (en) | For the micro logging chromatography imaging method of excitation well depth design | |
| CN105242328A (en) | Paleo-thermal lithosphere thickness determination method and device | |
| CN104237946B (en) | Single-layer reflected P-wave and reflection converted shear wave amplitude matching method based on well control | |
| CN106569259A (en) | Regularized tomographic velocity inversion method and device based on structural inclination angle | |
| CN105717538B (en) | Undulating surface seismic data migration datum plane conversion method and device | |
| CN106353799A (en) | Inversion method of united chromatography speed of longitudinal and cross waves | |
| CN107817524A (en) | The method and apparatus of three-dimensional seismic tomography | |
| CN116956420A (en) | Tunnel structure integrated analysis method based on three-dimensional geological model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |