The content of the invention
For the deficiencies in the prior art, an object of the present invention is solve present in above-mentioned prior art
One or more problems.For example, an object of the present invention is in the preferable side of pip in a kind of seismic reflection tomography is provided
Method, by the guidance at subsurface structure inclination angle, preferably go out and meet architectonic continuity pip.
To achieve these goals, it is of the invention to provide pip method for optimizing in a kind of seismic reflection tomography.It is described
Pip method for optimizing comprises the following steps:Form dip section;Point of the pickup with residual move out time on offset gather;Will
The spot projection with residual move out time picked up is on the dip section, obtaining chromatographing pip;Along dip section direction,
Choose successively and deviate the minimum chromatography pip of upper chromatography pip inclination direction as preferable chromatography pip.
Pip method for optimizing in seismic reflection tomography according to an exemplary embodiment of the present invention, it is described to choose preferable layer
Analysis pip step can include:Along dip section direction, determine a chromatography pip and adjacent thereto treat preferably
Chromatography pip line and horizontal line between angle, with it is described on the adjacent multiple chromatography pips of a pip
In, the minimum chromatography pip of the angle is selected as preferable chromatography pip.
Pip method for optimizing in seismic reflection tomography according to an exemplary embodiment of the present invention, a upper pip
The distance between multiple chromatography pips adjacent with a upper pip can be no more than 100~200 meters.
Pip method for optimizing in seismic reflection tomography according to an exemplary embodiment of the present invention, by obtaining subsurface reflective
The structure dip of point, forms the dip section.
Pip method for optimizing in seismic reflection tomography according to an exemplary embodiment of the present invention, it can use based on relevant
Property scan method, that the scan method based on complex seismic trace or the method based on gradient-structure tensor obtain the underground is anti-
The structure dip of exit point.
Another aspect of the present invention provides a kind of method of pre-stack depth migration velocity modeling.The pre-stack depth migration speed
The method of degree modeling includes:During pre-stack depth migration velocity modeling, using anti-in seismic reflection tomography as described above
The preferred pip of exit point method for optimizing.
The method of pre-stack depth migration velocity modeling according to an exemplary embodiment of the present, may comprise steps of:
(A1) pre-stack depth migration, is carried out to pre-stack depth migration rate pattern, output offset is away from trace gather;(A2) offset distance road, is judged
Whether collection, which reaches default, is evened up degree;(A3) if, offset gather be not reaching to it is default even up degree, carry out it is described preferably
The step of chromatographing pip;(A4) indirect ray tracking, is carried out on preferable pip;(A5) reflection tomography equation, is established
And solve;(A6), according to chromatography non trivial solution renewal pre-stack depth migration rate pattern;(A7), repeating said steps (A1) are extremely
(A6), until offset gather, which reaches default, evens up degree.Wherein it is possible to according to pre-stack time migration velocity model building prestack
Depth migration initial velocity model.
Compared with prior art, the present invention accurately efficiently preferably can go out to chromatograph pip, and control is anti-on the whole
The distribution of exit point, enhance the distribution rationality of chromatography pip.
Embodiment
Hereinafter, reflected in the seismic reflection tomography that the present invention will be described in detail with reference to accompanying drawing and exemplary embodiment
Point method for optimizing.
One aspect of the present invention provides a kind of method for optimizing of pip in seismic reflection tomography, and the method for optimizing passes through
Subsurface structure inclination angle, which preferably goes out, meets architectonic continuity pip.
Pip method for optimizing may include following steps in seismic reflection tomography according to an exemplary embodiment of the present invention:
A, the structure dip of subsurface reflection point is obtained by related algorithm (such as dip scanning method), inclination angle is formed and cuts open
Face.
Wherein, related algorithm may include dip scanning method, and conventional dip scanning method may include:Based on coherence's
Scan method, the scan method based on complex seismic trace or the method based on gradient-structure tensor, it is preferred that be based on relevant
The scan method of property.Scan method based on complex seismic trace and the method based on gradient-structure tensor are due to its scanning accuracy
Limitation is not widely used.Used simultaneously by majority with its higher computational accuracy based on the scan method of similitude
Development.Wherein, earliest is to propose the non-normalized cross-correlation of two dimension by Picou and Utzmann (1962) to scan, the calculating
Method is to carry out dip scanning along two-dimension earthquake survey line, and three-dimension curved surface is abstracted into two dimension carries out inclination angle estimation, resulting
Inclination data is in ribbon, and not accurate enough;Finn (1986) proposes three-dimensional dip and azimuthal definition, and right in two dimension
It is put into practice;Marfurt (1998,2000) proposes the three-dimensional dip scanning method based on similitude, with coherence properties
In C2 algorithms it is consistent, Finn and Backus (1986) algorithm refer on a real three dimensions;Due to the calculation
Method is to use centre scan window, a smoothing effect be present to the construction of analysis site both sides when causing to calculate, in order to improve this
Point, Marfurt (2006) proposes the mode of more window scannings, as a kind of sane inclination angle and azimuth method of estimation, the party
Method can be to the angle of heel of tomography two change scanning it is relatively accurate.
B, point of the pickup with residual move out time (RMO) on offset gather.Wherein, that picks up has residual move out time
The positional information of point has application follow-up.
C, the point (also known as pickup point) with residual move out time picked up is projected on dip section.
In this step, the pickup point with residual move out time on offset gather is projected on dip section, obtained
Chromatograph pip.
D, along dip section direction, choose deviate the minimum chromatography reflection of upper chromatography pip inclination direction successively
Point is as preferable chromatography pip.
Along dip section direction screening pickup point (chromatography pip).Point of proximity (i.e. layer can be judged according to inclination direction
Analyse pip) choice, specifically, the principle for screening cut-off be deviates a upper pip inclination direction minimum, that is, exist
Adjacent with last layer analysis pip treats preferably to chromatograph in pip, chooses the inclination direction for deviateing upper chromatography pip
Minimum chromatography pip is as preferable chromatography pip.Here, the selection of neighboring reflection point can determine according to distance,
For example, it is not more than 100~200 meters of point with last layer analysis pip distance, but the invention is not restricted to this, can also be as needed
Adjacent layer position pip is determined using other manner.
Wherein, dip section is that have value on each position, i.e., each chromatography pip on dip section
There is structure dip value.Inclination direction refers to the direction of the inclination value along the position with residual move out time point, as along the inclination angle
Direction, for example certain point inclination angle is 20 °, then inclination direction is to initially form 20 ° of direction with horizontal plane from certain described point.
In the present embodiment, choosing preferable chromatography pip step includes:Along dip section direction, one is determined
Angle between chromatography pip and the adjacent thereto line and horizontal line for treating preferable chromatography pip, with it is described on one
In the adjacent multiple chromatography pips of individual pip, the minimum chromatography pip of the angle is selected to be reflected as preferable chromatography
Point.
As shown in figure 1, point p is selected chromatography pip, carried out in point a, b, c, d, e adjacent with point p preferably,
Wherein, point p structure dip θ is 30 °, point a and point p angle δap(i.e. angle between point p and point a line and horizontal line)
For 51 °, point b and point p angle δbpFor 38 °, point c and point p angle δcpFor 32 °, point d and point p angle δdpFor 18 °, point e
With point p angle δepFor -23 °.Wherein, δcpDifference between θ is minimum, thus, point c is to be selected along structure dip direction
The point gone out, remaining point a, b, d and e point are the point given up.
It should be noted that the present invention preferably pip is by choosing angle closest to the angle of structure dip, if such as
The angle of two and above angle is identical, and is all then chosen and the closest layer of pip closest to the angle of structure dip
Analyse pip.
Another aspect of the present invention provides a kind of pre-stack depth migration velocity modeling method, and this method is in pre-stack depth migration
During velocity modeling, suitable pip is determined using the method for preferred pip in seismic reflection tomography as described above.
In one exemplary embodiment, can be included according to the prestack depth velocity shifts modeling method of the present invention following
Step:
(A1) pre-stack depth migration, is carried out to pre-stack depth migration rate pattern, output offset is away from trace gather.
(A2), judge whether offset gather reaches default and even up degree.
(A3) if, offset gather be not reaching to it is default even up degree, using in seismic reflection tomography as described above
It is preferred that the preferred pip of the method for pip, it is determined that rational pip.
(A4) indirect ray tracking, is carried out on preferable pip.Specifically, chased after using preferable pip as ray
Track point, using last round of pre-stack time migration rate pattern as ray tracing model, carry out indirect ray tracking;Wherein, the first round
Pre-stack time migration rate pattern is pre-stack depth migration initial velocity model according to pre-stack time migration velocity model building.
(A5), establish reflection tomography equation and solve.
(A6), according to chromatography non trivial solution renewal pre-stack depth migration rate pattern.
(A7), repeating said steps (A1) to (A6), until offset gather, which reaches default, evens up degree.It is default to even up journey
Degree can be actually needed the reasonable of determination according to and even up degree value, for example, offset gather is all evened up.
In summary, pip method for optimizing in seismic reflection tomography of the invention, preferably chromatographed by structure dip anti-
The distribution of exit point, on the whole control pip, the distribution rationality of chromatography pip is enhanced, and in pre-stack depth migration
Preferable application effect is obtained in modeling.
Although above by describing the present invention with reference to exemplary embodiment, those skilled in the art should be clear
Chu, in the case where not departing from the spirit and scope that claim is limited, the exemplary embodiment of the present invention can be carried out each
Kind modifications and changes.