CN105093280B - Surface-level model is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component - Google Patents
Surface-level model is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component Download PDFInfo
- Publication number
- CN105093280B CN105093280B CN201410207196.1A CN201410207196A CN105093280B CN 105093280 B CN105093280 B CN 105093280B CN 201410207196 A CN201410207196 A CN 201410207196A CN 105093280 B CN105093280 B CN 105093280B
- Authority
- CN
- China
- Prior art keywords
- mrow
- level model
- frequency component
- elevation
- low
- 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 provides a kind of surface-level model to the low frequency of earthquake data influence and the decomposition method of radio-frequency component, belong to seismic prospecting data data processing field.This method includes:(1) surface-level model data, smoothing range and work area replacement velocity parameter v are inputtedR;(2) the top circle elevation of surface-level model low-frequency component, i.e., smooth datum elevation are calculated;(3) the bottom circle elevation of surface-level model low-frequency component is calculated;(4) average speed of surface-level model low-frequency component is calculated;(5) radio-frequency component that surface-level model influences, i.e. time correcting value and elevation knots modification are calculated;(6) repeat step (2) completes the calculating of all measuring points to (5), then exports the time adjustment amount and elevation knots modification of the radio-frequency component that top circle's elevation, bottom circle elevation, average speed, the surface-level model of the low-frequency component influence.
Description
Technical field
The invention belongs to seismic prospecting data data processing field, and in particular to a kind of surface-level model is to earthquake data influence
Low frequency and radio-frequency component decomposition method.
Background technology
To be limited to by existing velocity modeling technology, relief surface skew is also needed to earthquake data application static correcting method,
Eliminate the HFS that surface-level model influences.The low frequency part that surface-level model influences is remained in data, by velocity modeling and
Skew solves.Correspondingly, after the HFS correction of surface-level model influence is carried out on geological data, surface-level model is only eliminated
The influence of horizontal fast changing portion, the low-frequency component that surface-level model influences also exist.The radio-frequency component that surface-level model influences includes
High frequency time correcting value and the part of elevation knots modification 2, the radio-frequency component correction that surface-level model influence is carried out on geological data include
Static shift to geological data and to shot point, receive the change of point height.Shot point and reception point height institute after change
Smooth reference plane is commonly referred to as in position.
The current surface-level model influence low frequency and the decomposition of radio-frequency component of carrying out has 2 routes.First, first calculate smooth benchmark
Face calculates high frequency time correcting value and elevation knots modification again, i.e., first earth's surface elevation is carried out smoothly obtaining a smooth reference plane,
Static correction value required for seismic data corrections to smooth reference plane is included into radio-frequency component, remaining is included into low-frequency component;Second,
First calculate high frequency time correcting value and calculate smooth reference plane and elevation knots modification again, i.e., first by certain principle from static correction value
High frequency time correcting value is decomposited, recycles the low-frequency component of time adjustment amount is counter to release smooth reference plane.This 2 route calculations
Time adjustment amount and elevation knots modification be all not only reflection to surface-level model transverse direction fast changing portion, contain to table
The stripping and filling of layer model weathering zone low-frequency component.
The content of the invention
It is an object of the invention to solve problem present in above-mentioned prior art, there is provided a kind of surface-level model is to earthquake number
According to the low frequency of influence and the decomposition method of radio-frequency component, according to actual surface-level model parameter and given smoothing range parameter, meter
The low-frequency component of surface-level model is calculated, and then the radio-frequency component for calculating surface-level model influence is that high frequency time correcting value and elevation change
Amount.The radio-frequency component represents during surface-level model influences laterally fast-changing part, available in being offset in relief surface to table
Layer model influence in laterally quick change and skew, modeling technique can not correct process part correction.To earthquake data application
After the high frequency time correcting value and elevation knots modification, influence of the surface-level model low-frequency component to geological data is not changed.
The present invention is achieved by the following technical solutions:
A kind of surface-level model to the low frequency of earthquake data influence and the decomposition method of radio-frequency component, including:
(1) surface-level model data, smoothing range and work area replacement velocity parameter v are inputtedR;
(2) the top circle elevation of surface-level model low-frequency component, i.e., smooth datum elevation are calculated;
(3) the bottom circle elevation of surface-level model low-frequency component is calculated;
(4) average speed of surface-level model low-frequency component is calculated;
(5) radio-frequency component that surface-level model influences, i.e. time correcting value and elevation knots modification are calculated;
(6) repeat step (2) completes the calculating of all measuring points to (5), then export the low-frequency component top circle elevation,
The time adjustment amount and elevation knots modification for the radio-frequency component that bottom circle elevation, average speed, surface-level model influence.
In the step (1), the surface-level model data include measuring point plane coordinates, push up boundary's elevation, bottom circle elevation, are averaged
Speed;The smoothing range is the smoothing range for calculating surface-level model low-frequency component;
For actual surface-level model by pushing up 3 boundary's elevation, bottom circle elevation and average speed parameter descriptions, top circle's elevation is measuring point height
Journey, average speed are the average speed for pushing up medium between boundary and bottom circle;Surface-level model low-frequency component is equally by pushing up boundary's elevation, bottom circle
The 3 parameter descriptions of elevation and average speed;
If the top circle elevation of actual surface-level model, bottom circle elevation, average speed are respectively e at the x measuring points of plan-positions(x)、
eb(x), v (x), top circle elevation, bottom circle elevation, the average speed of the surface-level model low-frequency component to be calculated use E respectivelys(x)、Eb
(x), V (x) is represented;
Given smoothing range is a two-dimentional window, and when seismic prospecting is a two-dimentional survey line, smoothing range deteriorates to
One-dimensional window,;
Relative to measuring point to be calculated, smoothing range includes measuring point to be calculated, and as the change of point position to be calculated is put down
Move.
What the step (2) was realized in:
The top circle elevation E of surface-level model low-frequency components(x) it is that actual surface-level model static correction value exists in given smoothing range
The top circle elevation of the surface-level model corresponding to weighted average in given smoothing range, i.e., in the smoothing range relative to x points
The weighted average of measuring point actual surface-level model top circle elevation:
A (x) is the set of the measuring point fallen into the smoothing range relative to x points in formula, and w (η-x) is weight coefficient.
What the step (3) was realized in:
The bottom circle elevation E of surface-level model low-frequency componentb(x) it is that actual surface-level model static correction value exists in given smoothing range
The bottom circle elevation of the surface-level model corresponding to weighted average in given smoothing range, i.e., in the smoothing range relative to x points
The weighted average of measuring point actual surface-level model bottom circle elevation:
A (x) is the set of the measuring point fallen into the smoothing range relative to x points in formula, and w (η-x) is weight coefficient.
What the step (4) was realized in:
The average speed V (x) of surface-level model low-frequency component is that actual surface-level model static correction value exists in given smoothing range
The speed of the surface-level model corresponding to weighted average in given smoothing range, i.e., with relative to being surveyed in the smoothing range of x points
The average speed of the actual surface-level model of point:
H (x)=e in formulas(x)-eb(x) be actual surface-level model thickness, A (x) is to fall into the smooth model relative to x points
The set of measuring point in enclosing, w (η-x) is weight coefficient.
The weight coefficient w (η-x) is a constant either with function relevant with x relativeness η;
When w (η-x) is with the η functions relevant with x relativeness, w (η-x) be with the distance between η and x increase and
The function of reduction, and w (η-x) is more than or equal to 0.
What the step (5) was realized in:
The radio-frequency component that surface-level model influences is the difference that actual surface-level model influences to influence with surface-level model low-frequency component
Embodiment in time adjustment amount and elevation knots modification;The time adjustment amount that surface-level model is influenceed in radio-frequency component is actual top layer
Model ing static correction amount and the difference of surface-level model low-frequency component static correction value:
Δ t (x)=- (es(x)-eb(x))/v(x)+(Es(x)-Eb(x))/V(x)+(Eb(x)-eb(x))/vR (6)
Elevation knots modification in the radio-frequency component that surface-level model influences is surface-level model low-frequency component crest level and actual table
The difference of layer model crest level:
Δ z (x)=Es(x)-es(x) (7)。
Compared with prior art, the beneficial effects of the invention are as follows:
The invention is applied to 3D real data, realizes the decomposition of the low frequency and radio-frequency component of surface-level model influence, obtains
The radio-frequency component that the low-frequency component and surface-level model of surface-level model influence.The radio-frequency component that surface-level model influences reflects top layer
Influence of the model transverse direction fast changing portion to geological data, carrying out this surface-level model to geological data influences radio-frequency component
Correction, does not change influence of the surface-level model low-frequency component to data.
The invention be relief surface migration imaging during more reasonably calculate with using the high frequency that surface-level model influences into
Divide and lay the foundation.
Brief description of the drawings
Fig. 1 is the actual surface-level model top in work area circle elevation (unit rice).
Fig. 2 is the actual surface-level model bottom in work area circle elevation (unit rice).
Fig. 3 is the actual surface-level model thickness in work area (unit rice).
Fig. 4 is the actual surface-level model average speed in work area (unit meter per second).
Fig. 5 is the top circle elevation (unit rice) for the work area surface-level model low-frequency component that embodiment 1 obtains.
Fig. 6 is the bottom circle elevation (unit rice) for the work area surface-level model low-frequency component that embodiment 1 obtains.
Fig. 7 is the thickness (unit rice) for the work area surface-level model low-frequency component that embodiment 1 obtains.
Fig. 8 is the average speed (unit meter per second) for the work area surface-level model low-frequency component that embodiment 1 obtains.
Fig. 9 is the time adjustment amount (unit millisecond) for the radio-frequency component that the work area surface-level model that embodiment 1 obtains influences.
Figure 10 is the elevation knots modification (unit rice) for the radio-frequency component that the work area surface-level model that embodiment 1 obtains influences.
Figure 11 is the top circle elevation (unit rice) for the work area surface-level model low-frequency component that embodiment 2 obtains.
Figure 12 is the bottom circle elevation (unit rice) for the work area surface-level model low-frequency component that embodiment 2 obtains.
Figure 13 is the thickness (unit rice) for the work area surface-level model low-frequency component that embodiment 2 obtains.
Figure 14 is the average speed (unit meter per second) for the work area surface-level model low-frequency component that embodiment 2 obtains.
Figure 15 is the time adjustment amount (unit millisecond) for the radio-frequency component that the work area surface-level model that embodiment 2 obtains influences.
Figure 16 is the elevation knots modification (unit rice) for the radio-frequency component that the work area surface-level model that embodiment 2 obtains influences.
Figure 17 is the top circle elevation (unit rice) for the work area surface-level model low-frequency component that embodiment 3 obtains.
Figure 18 is the bottom circle elevation (unit rice) for the work area surface-level model low-frequency component that embodiment 3 obtains.
Figure 19 is the thickness (unit rice) for the work area surface-level model low-frequency component that embodiment 3 obtains.
Figure 20 is the average speed (unit meter per second) for the work area surface-level model low-frequency component that embodiment 3 obtains.
Figure 21 is the time adjustment amount (unit millisecond) for the radio-frequency component that the work area surface-level model that embodiment 3 obtains influences.
Figure 22 is the elevation knots modification (unit rice) for the radio-frequency component that the work area surface-level model that embodiment 3 obtains influences.
Figure 23 is the step block diagram of the inventive method.
Embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings:
The radio-frequency component that surface-level model influences is the reflection of surface-level model transverse direction fast changing portion.Enter to geological data
After the radio-frequency component correction that row surface-level model influences, the influence of surface-level model transverse direction fast changing portion, top layer are only eliminated
Influence of the model low-frequency component to geological data also exists.The radio-frequency component that surface-level model influences includes high frequency time correcting value and height
The part of journey knots modification 2, the radio-frequency component that surface-level model influence is carried out on geological data correct the quiet school contained to geological data
Positive time shift and the change to shot point, reception point height.Shot point and reception point height position after change are commonly referred to as smoothly
Reference plane.
After given the smoothing range for calculating surface-level model low-frequency component, surface-level model low-frequency component is to give smooth model
Enclose the surface-level model corresponding to weighted average of the interior actual surface-level model static correction value in given smoothing range, surface-level model
The radio-frequency component of influence is that actual surface-level model is influenceed with the difference that surface-level model low-frequency component influences in time adjustment amount and height
Embodiment on journey knots modification.
For surface-level model by pushing up 3 boundary's elevation, bottom circle elevation and average speed parameter descriptions, top circle's elevation is measuring point elevation,
Average speed is the average speed for pushing up medium between boundary and bottom circle.Measuring point refers to shot point, receiving point.It is if real at the x measuring points of plan-position
Top circle elevation, bottom circle elevation, the average speed of border surface-level model are respectively es(x)、eb(x), v (x), then surface-level model thickness h
(x)=es(x)-eb(x), with elevation EdFor the measuring point static correction value of final reference plane
T (x)=- (es(x)-eb(x))/v(x)+(Ed-eb(x))/vR (1)
V in formulaRIt is replacement velocity.
Calculate the weighted average of all measuring point static correction values in the given smoothing range relative to x points(1) formula is substituted into, can be obtained
T (x)=- (Es(x)-Eb(x))/V(x)+(Ed-Eb(x))/vR (2)
Wherein Es(x)、Eb(x), V (x) represents the top circle elevation of surface-level model low-frequency component, bottom circle elevation and average respectively
Speed:
(3) formula is to fall into the test points set relatively in the given smoothing range of x points to the A (x) in (5) formula, and w (η-x) is
Weight coefficient.
The top circle elevation E of surface-level model low-frequency components(x) it is commonly known as smooth reference plane.Carried out to geological data
After the radio-frequency component correction that surface-level model influences, measuring point is corrected in the smooth reference plane.Surface-level model influence high frequency into
Time adjustment amount Δ t (x)=t (x)-T (x), elevation knots modification Δ z (x)=E divideds(x)-es(x)。
After given the smoothing range for calculating surface-level model low-frequency component, surface-level model low-frequency component is to give smooth model
Enclose the surface-level model corresponding to weighted average of the interior actual surface-level model static correction value in given smoothing range, surface-level model
The radio-frequency component of influence is that actual surface-level model is influenceed with the difference that surface-level model low-frequency component influences in time adjustment amount and height
Embodiment on journey knots modification.
As shown in figure 23, the inventive method is as follows:
(1) surface-level model data and smoothing range, work area replacement velocity parameter are inputted
Input includes measuring point plane coordinates, top circle's elevation, bottom circle elevation, the work area surface-level model and computational chart of average speed
The smoothing range of layer model low-frequency component;
For actual surface-level model by pushing up 3 boundary's elevation, bottom circle elevation and average speed parameter descriptions, top circle's elevation is measuring point height
Journey, average speed are the average speed for pushing up medium between boundary and bottom circle.Surface-level model low-frequency component is equally by pushing up boundary's elevation, bottom circle
The 3 parameter descriptions of elevation and average speed.
If the top circle elevation of actual surface-level model, bottom circle elevation, average speed are respectively e at the x measuring points of plan-positions(x)、
eb(x), v (x), top circle elevation, bottom circle elevation, the average speed of the surface-level model low-frequency component to be calculated use E respectivelys(x)、Eb
(x), V (x) is represented.
The given smoothing range for being used to calculate surface-level model low-frequency component is a two dimension (to be one-dimensional during two-dimentional survey line)
Window (two-dimentional window (circle, rectangle, ellipse etc., circle are given by length shaft length etc. by 2 length of sides, ellipse by radius, rectangle);
One-dimensional window is given by smooth length).
The size of smoothing range depends on purposes.Have when being used for those to surface relief and near-surface velocity cross directional variations
When having the relief surface migration processing of stronger processing power, relatively small smoothing range can be used, and works as and is used for those
Have to surface relief and near-surface velocity cross directional variations during the relief surface migration processing of weaker disposal ability it is necessary to using phase
To larger smoothing range.
Relative to measuring point to be calculated, smoothing range will include measuring point to be calculated, and with the change of point position to be calculated
Translation.
Input work area replacement velocity parameter vR。
(2) the top circle elevation of surface-level model low-frequency component, i.e., smooth datum elevation are calculated
The top circle elevation E of surface-level model low-frequency components(x) it is that actual surface-level model static correction value exists in given smoothing range
The top circle elevation of the surface-level model corresponding to weighted average in given smoothing range, is exactly specifically relative to the flat of x points
The weighted average of measuring point actual surface-level model top circle elevation in sliding scope:
A (x) is the set of the measuring point fallen into the smoothing range relative to x points in formula, and w (η-x) is weight coefficient.
(3) the bottom circle elevation of surface-level model low-frequency component is calculated
The bottom circle elevation E of surface-level model low-frequency componentb(x) it is that actual surface-level model static correction value exists in given smoothing range
The bottom circle elevation of the surface-level model corresponding to weighted average in given smoothing range, is exactly specifically relative to the flat of x points
The weighted average of measuring point actual surface-level model bottom circle elevation in sliding scope:
A (x) is the set of the measuring point fallen into the smoothing range relative to x points in formula, and w (η-x) is weight coefficient.
(4) average speed of surface-level model low-frequency component is calculated
The average speed V (x) of surface-level model low-frequency component is that actual surface-level model static correction value exists in given smoothing range
The speed of the surface-level model corresponding to weighted average in given smoothing range, is exactly specifically the smooth model relative to x points
Enclose the average speed of the actual surface-level model of interior measuring point:
H (x)=e in formulas(x)-eb(x) be actual surface-level model thickness, A (x) is to fall into the smooth model relative to x points
The set of measuring point in enclosing, w (η-x) is weight coefficient.
In above-mentioned (3), (4), (5) formula, weight coefficient w (η-x) can be a constant (such as 1) or with η and x
The relevant function of relativeness.When w (η-x) is with the η functions relevant with x relativeness, w (η-x) can be with η and x
The distance between increase and reduce function, and w (η-x) be more than or equal to 0.
Wherein, η and x is vector symbol, represents the point in plane.F (η), f (x) represent identical functional relation, simply
Argument character is different.Because (3), (4), (5) formula will calculate low frequency top layer data from original table layer data, so certainly
η is introduced outside variable symbol x again:X is used to represent point to be calculated, and η represents to participate in the point for calculating to be calculated result.
(5) radio-frequency component that surface-level model influences, i.e. time correcting value and elevation knots modification are calculated
The radio-frequency component that surface-level model influences is the difference that actual surface-level model influences to influence with surface-level model low-frequency component
Embodiment in time adjustment amount and elevation knots modification.The time adjustment amount that surface-level model is influenceed in radio-frequency component is actual top layer
Model ing static correction amount and the difference of surface-level model low-frequency component static correction value:
Δ t (x)=- (es(x)-eb(x))/v(x)+(Es(x)-Eb(x))/V(x)+(Eb(x)-eb(x))/vR (6)
(6) V in formulaRIt is the replacement velocity parameter for being used for static correction and calculating of input.
The elevation knots modification that surface-level model is influenceed in radio-frequency component is surface-level model low-frequency component crest level and actual top layer
The difference of model crest level:
Δ z (x)=Es(x)-es(x) (7)
(2) to (5) complete the calculating of all measuring points more than repeating.
(6) low-frequency component (top circle's elevation, bottom circle elevation, average speed) of work area measuring point surface-level model, top layer mould are exported
The radio-frequency component (time adjustment amount, elevation knots modification) that type influences.
In above-mentioned (1) formula into (7) formula, (1) formula is the fundamental formular of deriving static correction values, (2) formula to (5) formula be with
(1) based on formula, " weighted averages of all measuring point static correction values in the given smoothing range relative to x points " derived
As a result;(6) formula is the result that (1) formula subtracts (2) formula;(7) formula is that smooth datum elevation subtracts actual elevation.
The embodiment of the inventive method is as follows:
One three-dimensional work area surface-level model example.Fig. 1, Fig. 2 and Fig. 3 are the top circle elevation of actual surface-level model, bottom respectively
Boundary's elevation and average speed, Fig. 4 are the thickness (difference of Ji Ding circle elevations and bottom circle elevation) of actual surface-level model.
Embodiment 1
The smoothing range for calculating surface-level model low-frequency component is the circle of the radius R=1000m centered on measuring point to be calculated, is adopted
Weight coefficient is weighed with waiting.
(1) the top circle elevation of surface-level model low-frequency component is calculated
The top circle elevation of surface-level model low-frequency component is calculated using (3) formula, wherein A (x) is to fall into the radius centered on x points
The set of measuring point in R=1000m smoothing range, weight coefficient w (η-x)=1.
(2) the bottom circle elevation of surface-level model low-frequency component is calculated
The bottom circle elevation of surface-level model low-frequency component is calculated using (4) formula, wherein A (x) is to fall into the radius centered on x points
The set of measuring point in R=1000m smoothing range, weight coefficient w (η-x)=1.
(3) average speed of surface-level model low-frequency component is calculated
The average speed of surface-level model low-frequency component is calculated using (5) formula, wherein A (x) is to fall into the radius centered on x points
The set of measuring point in R=1000m smoothing range, weight coefficient w (η-x)=1.
(4) radio-frequency component that surface-level model influences is calculated
The time adjustment amount in the radio-frequency component of surface-level model influence is calculated using (6) formula, utilizes (7) formula to calculate top layer mould
Elevation knots modification in the radio-frequency component that type influences.
Fig. 5, Fig. 6, Fig. 7 and Fig. 8 be respectively the top circle elevation of the surface-level model low-frequency component being calculated, bottom circle elevation,
Average speed and thickness, Fig. 9 and Figure 10 be respectively the radio-frequency component that the surface-level model that is calculated influences time adjustment amount and
Elevation knots modification.
Embodiment 2
The smoothing range of surface-level model low-frequency component is the circle that radius is R=1000m centered on measuring point to be calculated, is used
A kind of weight coefficient with measuring point distance change to be calculated
In formula | η-x | it is the distance in plane between η points and x points, kc=1/(2R)、km=1/(8R)。
(1) the top circle elevation of surface-level model low-frequency component is calculated
The top circle elevation of surface-level model low-frequency component is calculated using (3) formula, wherein A (x) is to fall into the radius centered on x points
The set of measuring point in R=1000m smoothing range, weight coefficient w (η-x) are calculated by (8) formula.
(2) the bottom circle elevation of surface-level model low-frequency component is calculated
The bottom circle elevation of surface-level model low-frequency component is calculated using (4) formula, wherein A (x) is to fall into the radius centered on x points
The set of measuring point in R=1000m smoothing range, weight coefficient w (η-x) are calculated by (8) formula.
(3) average speed of surface-level model low-frequency component is calculated
The average speed of surface-level model low-frequency component is calculated using (5) formula, wherein A (x) is to fall into the radius centered on x points
The set of measuring point in R=1000m smoothing range, weight coefficient w (η-x) are calculated by (8) formula.
(4) radio-frequency component that surface-level model influences is calculated
The time adjustment amount in the radio-frequency component of surface-level model influence is calculated using (6) formula, utilizes (7) formula to calculate top layer mould
Elevation knots modification in the radio-frequency component that type influences.
Figure 11, Figure 12, Figure 13 and Figure 14 are the top circle elevation for the surface-level model low-frequency component being calculated, bottom circle height respectively
Journey, average speed and thickness, Figure 15 and Figure 16 are the time adjustment for the radio-frequency component that the surface-level model being calculated influences respectively
Amount and elevation knots modification.
Embodiment 3
The smoothing range of surface-level model low-frequency component is for centered on measuring point to be calculated, half length of side is respectively dx=
900m, dy=950m rectangle, the half of rectangle catercorner lengthUsing one kind with treating
Calculate the weight coefficient of measuring point distance change
In formula | η-x | it is the distance in plane between η points and x points, kc=1/(2R)、km=1/(8R)。
(1) the top circle elevation of surface-level model low-frequency component is calculated
The top circle elevation of surface-level model low-frequency component is calculated using (3) formula, wherein A (x) is to fall into one side of something centered on x points
The set of long measuring point respectively in dx=900m, dy=950m rectangle, weight coefficient w (η-x) are calculated by (9) formula.
(2) the bottom circle elevation of surface-level model low-frequency component is calculated
The bottom circle elevation of surface-level model low-frequency component is calculated using (4) formula, wherein A (x) is to fall into one side of something centered on x points
The set of long measuring point respectively in dx=900m, dy=950m rectangle, weight coefficient w (η-x) are calculated by (9) formula.
(3) average speed of surface-level model low-frequency component is calculated
The average speed of surface-level model low-frequency component is calculated using (5) formula, wherein A (x) is to fall into one side of something centered on x points
The set of long measuring point respectively in dx=900m, dy=950m rectangle, weight coefficient w (η-x) are calculated by (9) formula.
(4) radio-frequency component that surface-level model influences is calculated
The time adjustment amount in the radio-frequency component of surface-level model influence is calculated using (6) formula, utilizes (7) formula to calculate top layer mould
Elevation knots modification in the radio-frequency component that type influences.
Figure 17, Figure 18, Figure 19 and Figure 20 are the top circle elevation for the surface-level model low-frequency component being calculated, bottom circle height respectively
Journey, average speed and thickness, Figure 21 and Figure 22 are the time adjustment for the radio-frequency component that the surface-level model being calculated influences respectively
Amount and elevation knots modification.
Above-mentioned technical proposal is one embodiment of the present invention, for those skilled in the art, at this
On the basis of disclosure of the invention application process and principle, it is easy to make various types of improvement or deformation, be not limited solely to this
Invent the method described by above-mentioned embodiment, therefore previously described mode is simply preferable, and and without limitation
The meaning of property.
Claims (5)
1. a kind of surface-level model is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component, it is characterised in that:The side
Method includes:
(1) surface-level model data, smoothing range and work area replacement velocity parameter v are inputtedR;
(2) the top circle elevation of surface-level model low-frequency component, i.e., smooth datum elevation are calculated;
(3) the bottom circle elevation of surface-level model low-frequency component is calculated;
(4) average speed of surface-level model low-frequency component is calculated;
(5) radio-frequency component that surface-level model influences, i.e. time correcting value and elevation knots modification are calculated;
(6) repeat step (2) completes the calculating of all measuring points to (5), then exports top circle elevation, the bottom circle of the low-frequency component
The time adjustment amount and elevation knots modification for the radio-frequency component that elevation, average speed, surface-level model influence;
In the step (1), the surface-level model data include measuring point plane coordinates, top circle's elevation, bottom circle elevation, average speed
Degree;The smoothing range is the smoothing range for calculating surface-level model low-frequency component;
For actual surface-level model by pushing up 3 boundary's elevation, bottom circle elevation and average speed parameter descriptions, top circle's elevation is measuring point elevation,
Average speed is the average speed for pushing up medium between boundary and bottom circle;Surface-level model low-frequency component is equally by pushing up boundary's elevation, bottom circle height
The 3 parameter descriptions of journey and average speed;
If the top circle elevation of actual surface-level model, bottom circle elevation, average speed are respectively e at the x measuring points of plan-positions(x)、eb
(x), v (x), top circle elevation, bottom circle elevation, the average speed of the surface-level model low-frequency component to be calculated use E respectivelys(x)、Eb
(x), V (x) is represented;
Given smoothing range is a two-dimentional window, and when seismic prospecting is a two-dimentional survey line, smoothing range deteriorates to one
Tie up window;
Relative to measuring point to be calculated, smoothing range includes measuring point to be calculated, and is translated with the change of point position to be calculated;
What the step (5) was realized in:
The radio-frequency component that surface-level model influences be actual surface-level model influence with the difference that surface-level model low-frequency component influences when
Between embodiment on correcting value and elevation knots modification;The time adjustment amount that surface-level model is influenceed in radio-frequency component is actual surface-level model
The difference of static correction value and surface-level model low-frequency component static correction value:
Δ t (x)=- (es(x)-eb(x))/v(x)+(Es(x)-Eb(x))/V(x)+(Eb(x)-eb(x))/vR (6)
Elevation knots modification in the radio-frequency component that surface-level model influences is surface-level model low-frequency component crest level and actual top layer mould
The difference of type crest level:
Δ z (x)=Es(x)-es(x) (7)。
2. surface-level model according to claim 1 is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component, its
It is characterised by:What the step (2) was realized in:
The top circle elevation E of surface-level model low-frequency components(x) it is that actual surface-level model static correction value is giving in given smoothing range
The top circle elevation of the surface-level model corresponding to weighted average in smoothing range, i.e., relative to measuring point in the smoothing range of x points
The weighted average of actual surface-level model top circle elevation:
<mrow>
<msub>
<mi>E</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>&eta;</mi>
<mo>&Element;</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<msub>
<mi>e</mi>
<mi>s</mi>
</msub>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
<mi>w</mi>
<mo>(</mo>
<mrow>
<mi>&eta;</mi>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>/</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>&eta;</mi>
<mo>&Element;</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>-</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
A (x) is the set of the measuring point fallen into the smoothing range relative to x points in formula, and w (η-x) is weight coefficient.
3. surface-level model according to claim 1 is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component, its
It is characterised by:What the step (3) was realized in:
The bottom circle elevation E of surface-level model low-frequency componentb(x) it is that actual surface-level model static correction value is giving in given smoothing range
The bottom circle elevation of the surface-level model corresponding to weighted average in smoothing range, i.e., relative to measuring point in the smoothing range of x points
The weighted average of actual surface-level model bottom circle elevation:
<mrow>
<msub>
<mi>E</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>&Element;</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>e</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&eta;</mi>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>/</mo>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>&eta;</mi>
<mo>&Element;</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&eta;</mi>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
A (x) is the set of the measuring point fallen into the smoothing range relative to x points in formula, and w (η-x) is weight coefficient.
4. surface-level model according to claim 1 is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component, its
It is characterised by:What the step (4) was realized in:
The average speed V (x) of surface-level model low-frequency component is that actual surface-level model static correction value is giving in given smoothing range
The speed of the surface-level model corresponding to weighted average in smoothing range is that is, actual relative to measuring point in the smoothing range of x points
The average speed of surface-level model:
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>&eta;</mi>
<mo>&Element;</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<mi>h</mi>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
<mi>w</mi>
<mo>(</mo>
<mrow>
<mi>&eta;</mi>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>/</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>&eta;</mi>
<mo>&Element;</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>&eta;</mi>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>/</mo>
<mi>v</mi>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
H (x)=e in formulas(x)-eb(x) be actual surface-level model thickness, A (x) is fallen into the smoothing range relative to x points
Measuring point set, w (η-x) is weight coefficient.
5. the surface-level model according to any one of claim 2 to 4 is to the low frequency of earthquake data influence and point of radio-frequency component
Solution method, it is characterised in that:The weight coefficient w (η-x) is a constant either with letter relevant with x relativeness η
Number;
When w (η-x) is with the η functions relevant with x relativeness, w (η-x) is to increase with the distance between η and x and reduce
Function, and w (η-x) be more than or equal to 0.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410207196.1A CN105093280B (en) | 2014-05-15 | 2014-05-15 | Surface-level model is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410207196.1A CN105093280B (en) | 2014-05-15 | 2014-05-15 | Surface-level model is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105093280A CN105093280A (en) | 2015-11-25 |
CN105093280B true CN105093280B (en) | 2018-04-06 |
Family
ID=54574134
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410207196.1A Active CN105093280B (en) | 2014-05-15 | 2014-05-15 | Surface-level model is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105093280B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106932823B (en) * | 2015-12-31 | 2018-12-21 | 上海青凤致远地球物理地质勘探科技有限公司 | A kind of smooth datum level selection method |
CN107870365B (en) * | 2016-09-28 | 2019-04-12 | 中国石油化工股份有限公司 | High frequency static correction value modification method and system |
CN110967752B (en) * | 2018-09-29 | 2022-05-03 | 中国石油化工股份有限公司 | Static correction method for seismic data input by undulating surface migration |
CN112394406A (en) * | 2019-08-13 | 2021-02-23 | 中国石油化工股份有限公司 | Method for establishing simulated earth surface depth domain velocity model |
CN111827966B (en) * | 2020-03-25 | 2022-04-15 | 大庆油田有限责任公司 | Multi-well acoustic logging curve consistency processing method and device and storage medium |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6253157B1 (en) * | 1998-12-14 | 2001-06-26 | Exxonmobil Upstream Research Co. | Method for efficient manual inversion of seismic velocity information |
CN101072084A (en) * | 2006-05-11 | 2007-11-14 | 中国石油化工股份有限公司 | Wideband digital signal orthogonal beam frequency division method and its use |
WO2009077442A2 (en) * | 2007-12-14 | 2009-06-25 | Shell Internationale Research Maatschappij B.V. | Method of processing data obtained from seismic prospecting |
CN102841380B (en) * | 2011-06-22 | 2015-12-16 | 中国石油化工股份有限公司 | A kind of surface relief complex structure Earthquakes grouping of data is concerned with noise attenuation method |
-
2014
- 2014-05-15 CN CN201410207196.1A patent/CN105093280B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN105093280A (en) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105093280B (en) | Surface-level model is to the low frequency of earthquake data influence and the decomposition method of radio-frequency component | |
CN102879820B (en) | Based on the three-dimensional table layer model construction method of triangle gridding | |
CN102436671B (en) | Virtual viewpoint drawing method based on depth value non-linear transformation | |
CN104570102B (en) | Method for combining near-surface velocity model with middle-deep stratum velocity model | |
CN104330827B (en) | Surface model static correction processing method | |
CN103093223A (en) | Light spot image center quick locating method | |
CN103412333A (en) | Static correction base level determining method | |
CN104091065A (en) | Intermittent flow numerical simulation method for solving shallow water problem | |
CN104200453A (en) | Parallax image correcting method based on image segmentation and credibility | |
CN105242315B (en) | A kind of method and device of geological structure imaging | |
CN106776979A (en) | Vector Electronic Map increased quality automatic mode based on remote sensing | |
CN104392447B (en) | A kind of image matching method based on shade of gray | |
CN102338887A (en) | Irregular-size space-variant grid tomography imaging statics correction method | |
CN107870361A (en) | A kind of earthquake diving Wave chromatography imaging method, device and terminal device | |
CN108267743A (en) | Iteratively faster Underwater Navigation method based on fitting | |
CN105717538A (en) | Irregular surface seismic data migration datum plane transformation method and irregular surface seismic data migration datum plane transformation device | |
CN109671152B (en) | Line point cloud geometric characteristic estimation method based on Gaussian weighted discrete derivative | |
CN112016956B (en) | Ore grade estimation method and device based on BP neural network | |
Wirasaet et al. | Artificial boundary layers in discontinuous Galerkin solutions to shallow water equations in channels | |
CN106443829A (en) | Method and apparatus for constructing near-surface model | |
US10628552B2 (en) | Systems and methods for unfaulting point clouds | |
CN106446379B (en) | Any exploitation working face ground movement and deformation prediction method based on probability integration process | |
Zhu et al. | 3D unstructured spectral element method for frequency-domain airborne EM forward modeling based on Coulomb gauge | |
Jin et al. | Inversion study of vertical eddy viscosity coefficient based on an internal tidal model with the adjoint method | |
Stosic et al. | Computational approach to improving the efficiency of river discharge measurement |
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 |