CN105093280A - Method of decomposing low frequency and high frequency components of surface layer model influencing earthquake data - Google Patents

Abstract

The invention provides a method of decomposing the low frequency and high frequency components of a surface layer model influencing earthquake data, belonging to the seismic prospecting data processing field. The method comprises: (1) inputting surface layer model data, a smooth scope and a work area replacement velocity parameter vR; (2) calculating the top boundary elevation of surface layer model low frequency components, i.e., smooth reference plane elevation; (3) calculating the bottom boundary elevation of the surface layer model low frequency components; (4) calculating the average speed of the surface layer model low frequency components; (5) calculating surface layer model influence high frequency components, i.e., a time correcting value and an elevation changing value; and (6) repeating steps (2) to (5) to complete the calculation of all measuring points, and outputting the top boundary elevation, bottom boundary elevation, and average speed of the low frequency components, and the time correcting value and elevation changing value of the surface layer model influence high frequency components.

Description

The low frequency of surface-level model to earthquake data influence and the decomposition method of radio-frequency component

Technical field

The invention belongs to seismic prospecting data data processing field, be specifically related to the low frequency of a kind of surface-level model to earthquake data influence and the decomposition method of radio-frequency component.

Background technology

By the limitation of existing velocity modeling technology, relief surface skew also needs earthquake market demand static correcting method, eliminates the HFS of surface-level model impact.The low frequency part of surface-level model impact still retains in the data, is solved by velocity modeling and skew.Correspondingly, after HFS geological data being carried out to surface-level model impact corrects, only eliminate the impact of the horizontal fast changing portion of surface-level model, the low-frequency component of surface-level model impact also exists.The radio-frequency component of surface-level model impact comprises high frequency time correcting value and elevation knots modification 2 part, and radio-frequency component correction geological data being carried out to surface-level model impact contains the static shift of geological data and the change to shot point, acceptance point elevation.Shot point after change and acceptance point elevation position are often called as level and smooth reference field.

The current decomposition carried out surface-level model and affect low frequency and radio-frequency component has 2 routes.One is first calculate level and smooth reference field to calculate high frequency time correcting value and elevation knots modification again, namely first a level and smooth reference field is obtained to earth's surface elevation is smoothing, static correction value required for seismic data corrections to level and smooth reference field is included into radio-frequency component, and all the other are included into low-frequency component; Two is first calculate high frequency time correcting value to calculate level and smooth reference field and elevation knots modification again, namely first from static correction value, decomposites high frequency time correcting value by certain principle, and the low-frequency component of recycling time adjustment amount is counter releases level and smooth reference field.The time adjustment amount of these 2 route calculation and elevation knots modification are not only all the reflection to the horizontal fast changing portion of surface-level model, contain the stripping to surface-level model weathering zone low-frequency component and filling.

Summary of the invention

The object of the invention is to solve the difficult problem existed in above-mentioned prior art, the low frequency of a kind of surface-level model to earthquake data influence and the decomposition method of radio-frequency component are provided, according to actual surface-level model parameter and given smoothing range parameter, the low-frequency component of reckoner layer model, and then the radio-frequency component of reckoner layer model impact and high frequency time correcting value and elevation knots modification.This radio-frequency component to represent in surface-level model impact laterally fast-changing part, to be used in relief surface skew horizontal Rapid Variable Design in surface-level model impact and skew, modeling technique cannot the correction of correct processing section.After this high frequency time correcting value of earthquake market demand and elevation knots modification, do not change the impact of surface-level model low-frequency component on geological data.

The present invention is achieved by the following technical solutions:

The low frequency of surface-level model to earthquake data influence and a decomposition method for radio-frequency component, comprising:

(1) input table layer model data, smoothing range and work area replacement velocity parameter v _r;

(2) the top circle elevation of reckoner layer model low-frequency component, i.e. level and smooth datum elevation;

(3) end circle elevation of reckoner layer model low-frequency component;

(4) average velocity of reckoner layer model low-frequency component;

(5) radio-frequency component of reckoner layer model impact, i.e. time correcting value and elevation knots modification;

(6) repeat step (2) and complete the calculating of all measuring points to (5), then export the top circle elevation of described low-frequency component, the time adjustment amount of radio-frequency component that end circle elevation, average velocity, surface-level model affect and elevation knots modification.

In described step (1), described surface-level model data comprise measuring point planimetric coordinates, top circle elevation, end circle elevation, average velocity; Described smoothing range is the smoothing range of reckoner layer model low-frequency component;

Actual surface-level model describes by pushing up boundary's elevation, end circle elevation and average velocity 3 parameters, top circle's elevation and measuring point elevation, the average velocity of medium between average velocity Shi Ding circle and end circle; Surface-level model low-frequency component is described by top circle elevation, end circle elevation and average velocity 3 parameters equally;

If the top circle elevation of the actual surface-level model in x measuring point place, planimetric position, end circle elevation, average velocity are respectively e _s(x), e _bx (), v (x), the top circle elevation of the surface-level model low-frequency component that calculate, end circle elevation, average velocity use E respectively _s(x), E _bx (), V (x) represent;

Given smoothing range is a two-dimentional window, and when seismic prospecting is a two-dimentional survey line, smoothing range deteriorates to one dimension window;

Relative to measuring point to be calculated, smoothing range comprises measuring point to be calculated, and along with the change translation of point position to be calculated.

Described step (2) is achieved in that

The top circle elevation E of surface-level model low-frequency component _sx () is the top circle elevation of surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, namely relative to the weighted mean value of measuring point actual surface-level model top circle elevation in the smoothing range of x point:

$E_s\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(e_s\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}w(\mathrm{\η}-x)---\left(3\right)$

In formula, A (x) falls into the set relative to the measuring point in the smoothing range of x point, and w (η-x) is weighting coefficient.

Described step (3) is achieved in that

The end circle elevation E of surface-level model low-frequency component _bx () is the end circle elevation of surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, namely relative to the weighted mean value of boundary's elevation at the bottom of the actual surface-level model of measuring point in the smoothing range of x point:

$E_b\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(e_b\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}w(\mathrm{\η}-x)---\left(4\right)$

In formula, A (x) falls into the set relative to the measuring point in the smoothing range of x point, and w (η-x) is weighting coefficient.

Described step (4) is achieved in that

Average velocity V (x) of surface-level model low-frequency component is the speed of surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, namely with the average velocity of the actual surface-level model of measuring point in the smoothing range relative to x point:

$V\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(h\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}(\left(h\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/v\left(\mathrm{\η}\right))---\left(5\right)$

H (x) in formula=e _s(x)-e _bx () is the thickness of actual surface-level model, A (x) falls into the set relative to the measuring point in the smoothing range of x point, and w (η-x) is weighting coefficient.

The function that described weighting coefficient w (η-x) is a constant or follows the relativeness of η and x relevant;

When w (η-x) is the function of following the relativeness of η and x relevant, w (η-x) is the function increasing with the distance between η and x and reduce, and w (η-x) is more than or equal to 0.

Described step (5) is achieved in that

The radio-frequency component of surface-level model impact is the embodiment of difference in time adjustment amount and elevation knots modification that the impact of actual surface-level model affects with surface-level model low-frequency component; The surface-level model time adjustment amount affected in radio-frequency component is the difference of actual surface-level model static correction value and surface-level model low-frequency component static correction value:

Δt(x)＝-(e _s(x)-e _b(x))/v(x)+(E _s(x)-E _b(x))/V(x)+(E _b(x)-e _b(x))/v _R(6)

Elevation knots modification in the radio-frequency component of surface-level model impact is the difference of surface-level model low-frequency component crest level and actual surface-level model crest level:

Δz(x)＝E _s(x)-e _s(x)(7)。

Compared with prior art, the invention has the beneficial effects as follows:

This invention is applied to 3D real data, achieves the low frequency of surface-level model impact and the decomposition of radio-frequency component, obtains the low-frequency component of surface-level model and the radio-frequency component of surface-level model impact.The radio-frequency component of surface-level model impact reflects the horizontal fast changing portion of surface-level model to the impact of geological data, carries out the correction that this surface-level model affects radio-frequency component, do not change the impact of surface-level model low-frequency component on data to geological data.

This invention is that the radio-frequency component of more reasonably calculating and the impact of application table layer model in relief surface migration imaging process lays the foundation.

Accompanying drawing explanation

Fig. 1 is actual surface-level model top, work area circle elevation (unit rice).

Fig. 2 is boundary's elevation (unit rice) at the bottom of the actual surface-level model in work area.

Fig. 3 is the actual surface-level model thickness in work area (unit rice).

Fig. 4 is the actual surface-level model average velocity in work area (unit meter per second).

Fig. 5 is the top circle elevation (unit rice) of the work area surface-level model low-frequency component that embodiment 1 obtains.

Fig. 6 is the end circle elevation (unit rice) of the work area surface-level model low-frequency component that embodiment 1 obtains.

Fig. 7 is the thickness (unit rice) of the work area surface-level model low-frequency component that embodiment 1 obtains.

Fig. 8 is the average velocity (unit meter per second) of the work area surface-level model low-frequency component that embodiment 1 obtains.

Fig. 9 is the time adjustment amount (unit millisecond) of the radio-frequency component of the work area surface-level model impact that embodiment 1 obtains.

Figure 10 is the elevation knots modification (unit rice) of the radio-frequency component of the work area surface-level model impact that embodiment 1 obtains.

Figure 11 is the top circle elevation (unit rice) of the work area surface-level model low-frequency component that embodiment 2 obtains.

Figure 12 is the end circle elevation (unit rice) of the work area surface-level model low-frequency component that embodiment 2 obtains.

Figure 13 is the thickness (unit rice) of the work area surface-level model low-frequency component that embodiment 2 obtains.

Figure 14 is the average velocity (unit meter per second) of the work area surface-level model low-frequency component that embodiment 2 obtains.

Figure 15 is the time adjustment amount (unit millisecond) of the radio-frequency component of the work area surface-level model impact that embodiment 2 obtains.

Figure 16 is the elevation knots modification (unit rice) of the radio-frequency component of the work area surface-level model impact that embodiment 2 obtains.

Figure 17 is the top circle elevation (unit rice) of the work area surface-level model low-frequency component that embodiment 3 obtains.

Figure 18 is the end circle elevation (unit rice) of the work area surface-level model low-frequency component that embodiment 3 obtains.

Figure 19 is the thickness (unit rice) of the work area surface-level model low-frequency component that embodiment 3 obtains.

Figure 20 is the average velocity (unit meter per second) of the work area surface-level model low-frequency component that embodiment 3 obtains.

Figure 21 is the time adjustment amount (unit millisecond) of the radio-frequency component of the work area surface-level model impact that embodiment 3 obtains.

Figure 22 is the elevation knots modification (unit rice) of the radio-frequency component of the work area surface-level model impact that embodiment 3 obtains.

Figure 23 is the step block diagram of the inventive method.

Embodiment

Below in conjunction with accompanying drawing, the present invention is described in further detail:

The radio-frequency component of surface-level model impact is the reflection of the horizontal fast changing portion of surface-level model.After radio-frequency component geological data being carried out to surface-level model impact corrects, only eliminate the impact of the horizontal fast changing portion of surface-level model, the impact of surface-level model low-frequency component on geological data also exists.The radio-frequency component of surface-level model impact comprises high frequency time correcting value and elevation knots modification 2 part, and radio-frequency component correction geological data being carried out to surface-level model impact contains the static shift of geological data and the change to shot point, acceptance point elevation.Shot point after change and acceptance point elevation position are often called as level and smooth reference field.

After the given smoothing range of reckoner layer model low-frequency component, surface-level model low-frequency component is surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, and the radio-frequency component of surface-level model impact is the embodiment of difference in time adjustment amount and elevation knots modification that the impact of actual surface-level model affects with surface-level model low-frequency component.

Surface-level model describes by pushing up boundary's elevation, end circle elevation and average velocity 3 parameters, top circle's elevation and measuring point elevation, the average velocity of medium between average velocity Shi Ding circle and end circle.Measuring point refers to shot point, acceptance point.If the top circle elevation of the actual surface-level model in x measuring point place, planimetric position, end circle elevation, average velocity are respectively e _s(x), e _b(x), v (x), so surface-level model thickness h (x)=e _s(x)-e _bx (), with elevation E _dfor this measuring point static correction value of final reference field

t(x)＝-(e _s(x)-e _b(x))/v(x)+(E _d-e _b(x))/v _R(1)

V in formula _rit is replacement velocity.

Calculate the weighted mean value of all measuring point static correction values in the given smoothing range relative to x point $T\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(t\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}w(\mathrm{\η}-x),$ (1) formula is substituted into, can obtain

T(x)＝-(E _s(x)-E _b(x))/V(x)+(E _d-E _b(x))/v _R(2)

Wherein E _s(x), E _bx (), V (x) represent the top circle elevation of surface-level model low-frequency component, end circle elevation and average velocity respectively:

$E_s\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(e_s\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}w(\mathrm{\η}-x)---\left(3\right)$

$E_b\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(e_b\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}w(\mathrm{\η}-x)---\left(4\right)$

$V\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(h\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}(\left(h\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/v\left(\mathrm{\η}\right))---\left(5\right)$

(3) formula is fall into the test points set in the given smoothing range of relative x point to the A (x) in (5) formula, and w (η-x) is weighting coefficient.

The top circle elevation E of surface-level model low-frequency component _sx () is commonly called level and smooth reference field.After radio-frequency component geological data being carried out to surface-level model impact corrects, measuring point is corrected on this level and smooth reference field.Time adjustment amount Δ t (x)=t (x)-T (x), elevation knots modification Δ z (the x)=E of the radio-frequency component of surface-level model impact _s(x)-e _s(x).

After the given smoothing range of reckoner layer model low-frequency component, surface-level model low-frequency component is surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, and the radio-frequency component of surface-level model impact is the embodiment of difference in time adjustment amount and elevation knots modification that the impact of actual surface-level model affects with surface-level model low-frequency component.

As shown in figure 23, the inventive method is as follows:

(1) input table layer model data and smoothing range, work area replacement velocity parameter

Input comprises measuring point planimetric coordinates, top circle elevation, end circle elevation, the work area surface-level model of average velocity and the smoothing range of reckoner layer model low-frequency component;

Actual surface-level model describes by pushing up boundary's elevation, end circle elevation and average velocity 3 parameters, top circle's elevation and measuring point elevation, the average velocity of medium between average velocity Shi Ding circle and end circle.Surface-level model low-frequency component is described by top circle elevation, end circle elevation and average velocity 3 parameters equally.

If the top circle elevation of the actual surface-level model in x measuring point place, planimetric position, end circle elevation, average velocity are respectively e _s(x), e _bx (), v (x), the top circle elevation of the surface-level model low-frequency component that calculate, end circle elevation, average velocity use E respectively _s(x), E _bx (), V (x) represent.

The given smoothing range for reckoner layer model low-frequency component is a two dimension (being one dimension during two-dimentional survey line) window (two-dimentional window (circle, rectangle, ellipse etc., circle by radius, rectangle by 2 length of sides, oval given by major and minor axis length etc.); One dimension window is given by smooth length).

The size of smoothing range depends on purposes.When being used to those, to surface relief and near-surface velocity horizontal change, there is the relief surface migration processing of stronger processing power, relatively little smoothing range can be used, and when being used to those, to surface relief and near-surface velocity horizontal change, there is the relief surface migration processing of more weak processing power, relatively large smoothing range will be used.

Relative to measuring point to be calculated, smoothing range will comprise measuring point to be calculated, and along with the change translation of point position to be calculated.

Input work area replacement velocity parameter v _r.

(2) the top circle elevation of reckoner layer model low-frequency component, i.e. level and smooth datum elevation

The top circle elevation E of surface-level model low-frequency component _sx () is the top circle elevation of surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, be exactly the weighted mean value of measuring point actual surface-level model top circle elevation in the smoothing range relative to x point particularly:

$E_s\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(e_s\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}w(\mathrm{\η}-x)---\left(3\right)$

In formula, A (x) falls into the set relative to the measuring point in the smoothing range of x point, and w (η-x) is weighting coefficient.

(3) end circle elevation of reckoner layer model low-frequency component

The end circle elevation E of surface-level model low-frequency component _bx () is the end circle elevation of surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, be exactly the weighted mean value of boundary's elevation at the bottom of the actual surface-level model of measuring point in the smoothing range relative to x point particularly:

$E_b\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(e_b\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}w(\mathrm{\η}-x)---\left(4\right)$

In formula, A (x) falls into the set relative to the measuring point in the smoothing range of x point, and w (η-x) is weighting coefficient.

(4) average velocity of reckoner layer model low-frequency component

Average velocity V (x) of surface-level model low-frequency component is the speed of surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, is exactly the average velocity of the actual surface-level model of measuring point in the smoothing range relative to x point particularly:

$V\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(h\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}(\left(h\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/v\left(\mathrm{\η}\right))---\left(5\right)$

H (x) in formula=e _s(x)-e _bx () is the thickness of actual surface-level model, A (x) falls into the set relative to the measuring point in the smoothing range of x point, and w (η-x) is weighting coefficient.

In above-mentioned (3), (4), (5) formula, the function that weighting coefficient w (η-x) can be a constant (such as 1) or follow the relativeness of η and x relevant.When w (η-x) is the function of following the relativeness of η and x relevant, w (η-x) can be the function increasing with the distance between η and x and reduce, and w (η-x) is more than or equal to 0.

Wherein, η and x is vector symbol, represents the point in plane.F (η), f (x) represent identical funtcional relationship, and just argument character is different.Because (3), (4), (5) formula will calculate low frequency top layer data from original table layer data, so introduce again η: x for representing point to be calculated outside argument character x, η represents the point participating in calculating to be calculated some result.

(5) radio-frequency component of reckoner layer model impact, i.e. time correcting value and elevation knots modification

The radio-frequency component of surface-level model impact is the embodiment of difference in time adjustment amount and elevation knots modification that the impact of actual surface-level model affects with surface-level model low-frequency component.The surface-level model time adjustment amount affected in radio-frequency component is the difference of actual surface-level model static correction value and surface-level model low-frequency component static correction value:

Δt(x)＝-(e _s(x)-e _b(x))/v(x)+(E _s(x)-E _b(x))/V(x)+(E _b(x)-e _b(x))/v _R(6)

(6) V in formula _rit is the replacement velocity parameter calculated for static correction of input.

The surface-level model elevation knots modification affected in radio-frequency component is the difference of surface-level model low-frequency component crest level and actual surface-level model crest level:

Δz(x)＝E _s(x)-e _s(x)(7)

Repeat above (2) complete all measuring points calculating to (5).

(6) low-frequency component (top circle elevation, end circle elevation, average velocity) of work area measuring point surface-level model is exported, the radio-frequency component (time adjustment amount, elevation knots modification) of surface-level model impact.

In above-mentioned (1) formula in (7) formula, (1) formula is the fundamental formular of deriving static correction values, (2) formula is based on (1) formula to (5) formula, the result of the weighted mean value of all measuring point static correction values " in the given smoothing range relative to x point " derived; (6) formula is the result that (1) formula subtracts (2) formula; (7) formula is that level and smooth datum elevation subtracts actual elevation.

The embodiment of the inventive method is as follows:

A 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, end circle elevation and average velocity respectively, and Fig. 4 is the thickness (difference of Ji Ding circle elevation and end circle elevation) of actual surface-level model.

Embodiment 1

The smoothing range of reckoner layer model low-frequency component is the circle of radius R=1000m centered by measuring point to be calculated, adopts weighting coefficient such as the power of grade.

(1) the top circle elevation of reckoner layer model low-frequency component

Adopt the top circle elevation of (3) formula reckoner layer model low-frequency component, wherein A (x) be fall into radius R=1000m centered by x point smoothing range in the set of measuring point, weighting coefficient w (η-x)=1.

(2) end circle elevation of reckoner layer model low-frequency component

Adopt the end circle elevation of (4) formula reckoner layer model low-frequency component, wherein A (x) be fall into radius R=1000m centered by x point smoothing range in the set of measuring point, weighting coefficient w (η-x)=1.

(3) average velocity of reckoner layer model low-frequency component

Adopt the average velocity of (5) formula reckoner layer model low-frequency component, wherein A (x) be fall into radius R=1000m centered by x point smoothing range in the set of measuring point, weighting coefficient w (η-x)=1.

(4) radio-frequency component of reckoner layer model impact

Adopt the time adjustment amount in the radio-frequency component of (6) formula reckoner layer model impact, utilize the elevation knots modification in the radio-frequency component of (7) formula reckoner layer model impact.

Fig. 5, Fig. 6, Fig. 7 and Fig. 8 are the top circle elevation of the surface-level model low-frequency component calculated, end circle elevation, average velocity and thickness respectively, Fig. 9 and Figure 10 is time adjustment amount and the elevation knots modification of the radio-frequency component of the surface-level model impact calculated respectively.

Embodiment 2

The smoothing range of surface-level model low-frequency component is the circle that radius is R=1000m centered by measuring point to be calculated, adopts a kind of weight coefficient with changing with measuring point to be calculated distance

$w(\mathrm{\η}-x)=\frac{\sin \left(2\mathrm{\π}{k}_c\right|\mathrm{\η}-x\left|\right)}{\mathrm{\π}|\mathrm{\η}-x|}\cos \left(2\mathrm{\π}{k}_m\right|\mathrm{\η}-x\left|\right)---\left(8\right)$

In formula | η-x| is the distance in plane between η point and x point, k _c=1/ (2R), k _m=1/ (8R).

(1) the top circle elevation of reckoner layer model low-frequency component

Adopt the top circle elevation of (3) formula reckoner layer model low-frequency component, wherein A (x) be fall into radius R=1000m centered by x point smoothing range in the set of measuring point, weighting coefficient w (η-x) is calculated by (8) formula.

(2) end circle elevation of reckoner layer model low-frequency component

Adopt the end circle elevation of (4) formula reckoner layer model low-frequency component, wherein A (x) be fall into radius R=1000m centered by x point smoothing range in the set of measuring point, weighting coefficient w (η-x) is calculated by (8) formula.

(3) average velocity of reckoner layer model low-frequency component

Adopt the average velocity of (5) formula reckoner layer model low-frequency component, wherein A (x) be fall into radius R=1000m centered by x point smoothing range in the set of measuring point, weighting coefficient w (η-x) is calculated by (8) formula.

(4) radio-frequency component of reckoner layer model impact

Adopt the time adjustment amount in the radio-frequency component of (6) formula reckoner layer model impact, utilize the elevation knots modification in the radio-frequency component of (7) formula reckoner layer model impact.

Figure 11, Figure 12, Figure 13 and Figure 14 are the top circle elevation of the surface-level model low-frequency component calculated, end circle elevation, average velocity and thickness respectively, Figure 15 and Figure 16 is time adjustment amount and the elevation knots modification of the radio-frequency component of the surface-level model impact calculated respectively.

Embodiment 3

The smoothing range of surface-level model low-frequency component is that half length of side is respectively the rectangle of dx=900m, dy=950m for centered by measuring point to be calculated, the half of rectangle catercorner length adopt a kind of weight coefficient with changing with measuring point to be calculated distance

$w(\mathrm{\η}-x)=\frac{\sin \left(2\mathrm{\π}{k}_c\right|\mathrm{\η}-x\left|\right)}{\mathrm{\π}|\mathrm{\η}-x|}\cos \left(2\mathrm{\π}{k}_m\right|\mathrm{\η}-x\left|\right)---\left(9\right)$

In formula | η-x| is the distance in plane between η point and x point, k _c=1/ (2R), k _m=1/ (8R).

(1) the top circle elevation of reckoner layer model low-frequency component

Adopt the top circle elevation of (3) formula reckoner layer model low-frequency component, wherein A (x) is the set falling into measuring point in rectangle that half length of side centered by x point is respectively dx=900m, dy=950m, and weighting coefficient w (η-x) is calculated by (9) formula.

(2) end circle elevation of reckoner layer model low-frequency component

Adopt the end circle elevation of (4) formula reckoner layer model low-frequency component, wherein A (x) is the set falling into measuring point in rectangle that half length of side centered by x point is respectively dx=900m, dy=950m, and weighting coefficient w (η-x) is calculated by (9) formula.

(3) average velocity of reckoner layer model low-frequency component

Adopt the average velocity of (5) formula reckoner layer model low-frequency component, wherein A (x) is the set falling into measuring point in rectangle that half length of side centered by x point is respectively dx=900m, dy=950m, and weighting coefficient w (η-x) is calculated by (9) formula.

(4) radio-frequency component of reckoner layer model impact

Adopt the time adjustment amount in the radio-frequency component of (6) formula reckoner layer model impact, utilize the elevation knots modification in the radio-frequency component of (7) formula reckoner layer model impact.

Figure 17, Figure 18, Figure 19 and Figure 20 are the top circle elevation of the surface-level model low-frequency component calculated, end circle elevation, average velocity and thickness respectively, Figure 21 and Figure 22 is time adjustment amount and the elevation knots modification of the radio-frequency component of the surface-level model impact calculated respectively.

Technique scheme is one embodiment of the present invention, for those skilled in the art, on the basis that the invention discloses application process and principle, be easy to make various types of improvement or distortion, and the method be not limited only to described by the above-mentioned embodiment of the present invention, therefore previously described mode is just preferred, and does not have restrictive meaning.

Claims (7)

1. the low frequency of surface-level model to earthquake data influence and a decomposition method for radio-frequency component, is characterized in that: described method comprises:

(1) input table layer model data, smoothing range and work area replacement velocity parameter v _r;

(2) the top circle elevation of reckoner layer model low-frequency component, i.e. level and smooth datum elevation;

(3) end circle elevation of reckoner layer model low-frequency component;

(4) average velocity of reckoner layer model low-frequency component;

(5) radio-frequency component of reckoner layer model impact, i.e. time correcting value and elevation knots modification;

(6) repeat step (2) and complete the calculating of all measuring points to (5), then export the top circle elevation of described low-frequency component, the time adjustment amount of radio-frequency component that end circle elevation, average velocity, surface-level model affect and elevation knots modification.

2. the low frequency of surface-level model according to claim 1 to earthquake data influence and the decomposition method of radio-frequency component, it is characterized in that: in described step (1), described surface-level model data comprise measuring point planimetric coordinates, top circle elevation, end circle elevation, average velocity; Described smoothing range is the smoothing range of reckoner layer model low-frequency component;

Actual surface-level model describes by pushing up boundary's elevation, end circle elevation and average velocity 3 parameters, top circle's elevation and measuring point elevation, the average velocity of medium between average velocity Shi Ding circle and end circle; Surface-level model low-frequency component is described by top circle elevation, end circle elevation and average velocity 3 parameters equally;

If the top circle elevation of the actual surface-level model in x measuring point place, planimetric position, end circle elevation, average velocity are respectively e _s(x), e _bx (), v (x), the top circle elevation of the surface-level model low-frequency component that calculate, end circle elevation, average velocity use E respectively _s(x), E _bx (), V (x) represent;

Given smoothing range is a two-dimentional window, and when seismic prospecting is a two-dimentional survey line, smoothing range deteriorates to one dimension window;

Relative to measuring point to be calculated, smoothing range comprises measuring point to be calculated, and along with the change translation of point position to be calculated.

3. the low frequency of surface-level model according to claim 2 to earthquake data influence and the decomposition method of radio-frequency component, is characterized in that: described step (2) is achieved in that

The top circle elevation E of surface-level model low-frequency component _sx () is the top circle elevation of surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, namely relative to the weighted mean value of measuring point actual surface-level model top circle elevation in the smoothing range of x point:

$E_s\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(e_s\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}w(\mathrm{\η}-x)---\left(3\right)$

In formula, A (x) falls into the set relative to the measuring point in the smoothing range of x point, and w (η-x) is weighting coefficient.

4. the low frequency of surface-level model according to claim 2 to earthquake data influence and the decomposition method of radio-frequency component, is characterized in that: described step (3) is achieved in that

The end circle elevation E of surface-level model low-frequency component _bx () is the end circle elevation of surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, namely relative to the weighted mean value of boundary's elevation at the bottom of the actual surface-level model of measuring point in the smoothing range of x point:

$E_b\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(e_b\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}w(\mathrm{\η}-x)---\left(4\right)$

In formula, A (x) falls into the set relative to the measuring point in the smoothing range of x point, and w (η-x) is weighting coefficient.

5. the low frequency of surface-level model according to claim 2 to earthquake data influence and the decomposition method of radio-frequency component, is characterized in that: described step (4) is achieved in that

Average velocity V (x) of surface-level model low-frequency component is the speed of surface-level model weighted mean value corresponding to of actual surface-level model static correction value in given smoothing range in given smoothing range, namely relative to the average velocity of the actual surface-level model of measuring point in the smoothing range of x point:

$V\left(x\right)=\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}\left(h\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/\underset{\mathrm{\η}\∈A\left(x\right)}{\mathrm{\Σ}}(\left(h\left(\mathrm{\η}\right)w(\mathrm{\η}-x)\right)/v\left(\mathrm{\η}\right))---\left(5\right)$

H (x) in formula=e _s(x)-e _bx () is the thickness of actual surface-level model, A (x) falls into the set relative to the measuring point in the smoothing range of x point, and w (η-x) is weighting coefficient.

6. according to the arbitrary described low frequency of surface-level model to earthquake data influence of claim 3 to 5 and the decomposition method of radio-frequency component, it is characterized in that: the function that described weighting coefficient w (η-x) is a constant or follows the relativeness of η and x relevant;

When w (η-x) is the function of following the relativeness of η and x relevant, w (η-x) is the function increasing with the distance between η and x and reduce, and w (η-x) is more than or equal to 0.

7. the low frequency of surface-level model according to claim 2 to earthquake data influence and the decomposition method of radio-frequency component, is characterized in that: described step (5) is achieved in that

The radio-frequency component of surface-level model impact is the embodiment of difference in time adjustment amount and elevation knots modification that the impact of actual surface-level model affects with surface-level model low-frequency component; The surface-level model time adjustment amount affected in radio-frequency component is the difference of actual surface-level model static correction value and surface-level model low-frequency component static correction value:

Δt(x)＝-(e _s(x)-e _b(x))/v(x)+(E _s(x)-E _b(x))/V(x)+(E _b(x)-e _b(x))/vR(6)

Elevation knots modification in the radio-frequency component of surface-level model impact is the difference of surface-level model low-frequency component crest level and actual surface-level model crest level:

Δz(x)＝E _s(x)-e _s(x)(7)。