CN107870365A - High frequency static correction value modification method and system - Google Patents

High frequency static correction value modification method and system Download PDF

Info

Publication number
CN107870365A
CN107870365A CN201610862929.4A CN201610862929A CN107870365A CN 107870365 A CN107870365 A CN 107870365A CN 201610862929 A CN201610862929 A CN 201610862929A CN 107870365 A CN107870365 A CN 107870365A
Authority
CN
China
Prior art keywords
high frequency
correction value
landform
static correction
frequency static
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.)
Granted
Application number
CN201610862929.4A
Other languages
Chinese (zh)
Other versions
CN107870365B (en
Inventor
林伯香
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201610862929.4A priority Critical patent/CN107870365B/en
Publication of CN107870365A publication Critical patent/CN107870365A/en
Application granted granted Critical
Publication of CN107870365B publication Critical patent/CN107870365B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking

Abstract

Disclose a kind of high frequency static correction value modification method and system.This method can include:Based on terrain elevation data, high frequency static correction value and smooth radius, landform altitude high fdrequency component is obtained;Based on high frequency static correction value and landform altitude high fdrequency component, similarity factor is obtained;Based on landform altitude high fdrequency component and similarity factor, landform correlated components are obtained;And based on high frequency static correction value and landform correlated components, obtain revised high frequency static correction value.The present invention obtains high-precision revised high frequency static correction value by extracting irrational static correction composition related to landform and being rejected.

Description

High frequency static correction value modification method and system
Technical field
The present invention relates to field of seismic exploration, is risen and fallen related point more particularly, to a kind of high frequency static correction value mesorelief Measure extracting method and system.
Background technology
In field of seismic exploration, static correction is an important step in seismic processing chain.Counted based on first arrival The high frequency static correction value of method calculates, be all so that first arrival time curve it is smoothened be target.It is but acute in some topography variations Strong mountain area, because first arrival ray path and reflection raypaths are widely different, first arrival smoothly with rational reflection line-ups when Away from not fully consistent between relation, these high frequency static correction value compositions related to hypsography are generally rational to acquisition Reflection line-ups time-distance equation is unfavorable, it should is differentiated as much as possible with rejecting.
Inventor's discovery, can be in a way using method smooth in certain limit is carried out to high frequency static correction value Extract this unfavorable high frequency static correction value composition and rejected, but this simple smoothing method be possible to those with Hypsography is unrelated, is probably the correct high frequency static correction value composition caused by the reasons such as near-surface model error also picks Remove.Therefore, it is necessary to develop a kind of high-precision high frequency static correction value modification method and system.
The information for being disclosed in background of invention part is merely intended to deepen the reason of the general background technology to the present invention Solution, and be not construed as recognizing or imply known to those skilled in the art existing of the information structure in any form Technology.
The content of the invention
The present invention proposes a kind of high frequency static correction value modification method and system, and it can be related to landform by extracting Irrational static correction composition is simultaneously rejected, and obtains high-precision revised high frequency static correction value.
According to an aspect of the invention, it is proposed that a kind of high frequency static correction value modification method.Methods described can include:Base In terrain elevation data, high frequency static correction value and smooth radius, landform altitude high fdrequency component is obtained;Based on the high frequency static correction Amount and the landform altitude high fdrequency component, obtain similarity factor;Based on the landform altitude high fdrequency component and the similarity factor, Obtain landform correlated components;And based on the high frequency static correction value and the landform correlated components, obtain revised high frequency Static correction value.
According to another aspect of the invention, it is proposed that a kind of high frequency static correction value update the system, the system can include: For based on terrain elevation data, high frequency static correction value and smooth radius, the unit of acquisition landform altitude high fdrequency component;For base In the high frequency static correction value and the landform altitude high fdrequency component, the unit of similarity factor is obtained;For based on the landform Elevation high fdrequency component and the similarity factor, obtain the unit of landform correlated components;And for based on the high frequency static correction Amount and the landform correlated components, obtain the unit of revised high frequency static correction value.
Methods and apparatus of the present invention has other characteristics and advantage, and these characteristics and advantage are attached from what is be incorporated herein It will be apparent in figure and subsequent embodiment, or by the accompanying drawing being incorporated herein and subsequent specific reality Apply in mode and stated in detail, these the drawings and specific embodiments are provided commonly for explaining the certain principles of the present invention.
Brief description of the drawings
Exemplary embodiment of the invention is described in more detail in conjunction with the accompanying drawings, it is of the invention above-mentioned and its Its purpose, feature and advantage will be apparent, wherein, in exemplary embodiment of the invention, identical reference number Typically represent same parts.
Fig. 1 shows the flow chart of the step of high frequency static correction value modification method according to the present invention.
Fig. 2 shows the schematic diagram of the receiving point landform altitude according to an embodiment of the invention.
Fig. 3 shows the schematic diagram of the high frequency static correction value according to an embodiment of the invention.
Fig. 4 is shown removes low frequency and the terrain elevation data and landform of superfrequency according to an embodiment of the invention The schematic diagram of elevation high fdrequency component.
Fig. 5 shows the landform correlated components and revised high frequency static correction value according to an embodiment of the invention Schematic diagram.
Fig. 6 shows the schematic diagram of the stacked section of the high frequency static correction value according to prior art.
Fig. 7 shows showing for the stacked section of the revised high frequency static correction value according to an embodiment of the invention It is intended to.
Fig. 8 is shown removes low frequency and the terrain elevation data and landform of superfrequency according to an embodiment of the invention The schematic diagram of elevation high fdrequency component.
Fig. 9 shows the landform correlated components and revised high frequency static correction value according to an embodiment of the invention Schematic diagram.
Figure 10 is shown removes low frequency and the terrain elevation data of superfrequency and ground according to an embodiment of the invention The schematic diagram of shape elevation high fdrequency component.
Figure 11 shows the landform correlated components and revised high frequency static correction according to an embodiment of the invention The schematic diagram of amount.
Embodiment
The present invention is more fully described below with reference to accompanying drawings.Although the side of being preferable to carry out of the present invention is shown in accompanying drawing Formula, however, it is to be appreciated that may be realized in various forms the present invention without should be limited by embodiments set forth herein.Phase Instead, there is provided these embodiments be in order that the present invention is more thorough and complete, and can be by the scope of the present invention intactly It is communicated to those skilled in the art.
Embodiment 1
Fig. 1 shows the flow chart of the step of high frequency static correction value modification method according to the present invention.
In this embodiment, can be included according to the high frequency static correction value modification method of the present invention:Step 101, it is based on Terrain elevation data, high frequency static correction value and smooth radius, obtain landform altitude high fdrequency component;Step 102, based on the quiet school of high frequency Positive quantity and landform altitude high fdrequency component, obtain similarity factor;Step 103, based on landform altitude high fdrequency component and similarity factor, obtain Obtain landform correlated components;And step 104, based on high frequency static correction value and landform correlated components, it is quiet to obtain revised high frequency Correcting value.
The embodiment obtains high accuracy by extracting irrational static correction composition related to landform and being rejected Revised high frequency static correction value.
The following detailed description of the specific steps of the high frequency static correction value modification method according to the present invention.
In one example, based on terrain elevation data, high frequency static correction value and smooth radius, landform altitude high frequency is obtained Component.
In one example, obtaining landform altitude high fdrequency component can include:Based on terrain elevation data and smooth radius, Obtain the terrain elevation data for removing low frequency;Based on the terrain elevation data for removing low frequency, special radio-frequency head is removed by smoothing processing Point, low frequency and the terrain elevation data of superfrequency are removed in acquisition;And based on high frequency static correction value and remove low frequency and the ground of superfrequency Shape altitude data, obtain landform altitude high fdrequency component.
Specifically, terrain elevation data, high frequency static correction value and smooth radius can be based on, wherein, smooth radius includes: First smooth radius, the second smooth radius and the 3rd smooth radius.Low frequency point is rejected from terrain elevation data by (1) formula Amount, obtain the terrain elevation data for removing low frequency:
Wherein, e1(x) terrain elevation data of low frequency, e are removed in expression0(x) landform altitude at the x measuring points of plan-position is represented, R represents the first smooth radius, and A (x) represents to fall into the set for the measuring point being less than or equal in the range of R relative to the distance of x points, w1 (η-x) is weight coefficient.
Can be to removing the terrain elevation data of low frequency, (2) formula of use is smoothed, and rejects the part of special high frequency, Low frequency and the terrain elevation data of superfrequency are removed in acquisition:
Wherein, e2(x) represent to remove low frequency and the terrain elevation data of superfrequency, RBThe second smooth radius is represented, B (x) is represented Fall into and be less than or equal to R relative to the distance of x pointsBIn the range of measuring point set, w2(η-x) is weight coefficient, RB[R/ can be taken 3,2R/3] value between, works as RBWhen=0, e is directly used2(x)=e1(x), i.e., without the smooth of this step, people in the art Member can carry out computing as the case may be.
Low frequency and superfrequency terrain elevation data ratio is gone into high-frequency components of statics to there is identical to put down with (3) formula handle The level of equal amplitude, obtains landform altitude high fdrequency component:
Wherein, ehf(x) landform altitude high fdrequency component, t are represented0(x) high frequency static correction value is represented.
In one example, based on high frequency static correction value and landform altitude high fdrequency component, similarity factor is obtained.
In one example, obtaining similarity factor can include:Based on high frequency static correction value and landform altitude high fdrequency component, Obtain preliminary similarity factor;And based on preliminary similarity factor, handled by weighting, obtain similarity factor.
Specifically, high frequency static correction value and landform altitude high fdrequency component can be based on, utilizes (4) formula, is obtained preliminary similar Coefficient:
Wherein, s1(x) preliminary similarity factor is represented.
Preliminary similarity factor is weighted with (5) formula again, high frequency static correction value and the high frequency division of landform altitude can be obtained The similarity factor of amount:
Wherein, scThreshold value is represented, can be with value between [0.5,0.7], to the preliminary similar system less than the threshold value Number is further weighted.Those skilled in the art can set threshold value as the case may be.
In one example, based on landform altitude high fdrequency component and similarity factor, landform correlated components are obtained.
In one example, obtaining landform correlated components can include:Based on landform altitude high fdrequency component and similarity factor, Obtain preliminary landform correlated components;And based on preliminary landform correlated components, by smoothing processing, obtain landform correlated components.
Specifically, landform altitude high fdrequency component and similarity factor can be based on, utilizes (6) formula, it is related to obtain preliminary landform Component:
t1(x)=ehf(x)s2(x) (6)
Wherein, t1(x) preliminary landform correlated components are represented.
Preliminary landform correlated components are carried out with (7) formula again smoothly, landform correlated components can be obtained:
Wherein, t2(x) landform correlated components, R are representedCThe 3rd smooth radius is represented, C (x) represents to fall into relative to x points Distance is less than or equal to RCIn the range of measuring point set, w3(η-x) is weight coefficient.RCThe value between [R/3,2R/3] can be taken, Work as RCWhen=0, t is directly used2(x)=t1(x), i.e., can be according to specific feelings without the smooth of this step, those skilled in the art Condition carries out value computing.
Wherein, for weight coefficient w1(η-x)、w2(η-x) and w3(η-x), can be unrelated with x relativeness with η One constant or with function relevant with x relativeness η.Work as w1(η-x)、w2(η-x) and w3(η-x) is with η and x Relativeness relevant function when, w1(η-x)、w2(η-x) and w3(η-x) is to increase with the distance between η and x and reduce Function, and w1(η-x)、w2(η-x) and w3(η-x) is all higher than being equal to 0.Those skilled in the art can set as the case may be Weight coefficient w1(η-x)、w2(η-x) and w3(η-x)。
In one example, based on high frequency static correction value and landform correlated components, the quiet school of revised high frequency can be obtained Positive quantity.
In one example, obtaining revised high frequency static correction value can include:High frequency static correction value is related to landform Component subtracts each other, and obtains revised high frequency static correction value.
Specifically, by (8) formula, high frequency static correction value and landform correlated components is subtracted each other, it is quiet to obtain revised high frequency Correcting value:
thf(x)=t0(x)-t2(x) (8)
Wherein, thf(x) revised high frequency static correction value is represented.
For ease of understanding the scheme of embodiment of the present invention and its effect, three concrete application examples given below.Ability Field technique personnel should be understood that the example only for the purposes of understanding the present invention, and its any detail is not intended in any way The limitation present invention.
Using example 1
Fig. 2 shows the schematic diagram of the receiving point landform altitude according to an embodiment of the invention.Fig. 3 shows root According to the schematic diagram of the high frequency static correction value of an embodiment of the invention.
In this example, R=600m, R are takenB=0, RC=300m, w1(η-x)、w2(η-x) and w3(η-x) three weight functions Use the constant 1 unrelated with x relativeness with η, threshold value sc=0.5.
Based on terrain elevation data, high frequency static correction value and smooth radius, wherein, smooth radius can include:First is flat Sliding radius, the second smooth radius and the 3rd smooth radius.Low frequency component is rejected from terrain elevation data by (1) formula, obtained Remove the terrain elevation data of low frequency:
Wherein, e1(x) terrain elevation data of low frequency, e are removed in expression0(x) landform altitude at the x measuring points of plan-position is represented, R represents the first smooth radius, and A (x) represents to fall into the set for the measuring point being less than or equal in the range of R relative to the distance of x points, w1 (η-x) is weight coefficient, w1(η-x)=1.
Terrain elevation data to removing low frequency, (2) formula of use are smoothed, and reject the part of special high frequency, are obtained Remove low frequency and the terrain elevation data of superfrequency:
Wherein, e2(x) represent to remove low frequency and the terrain elevation data of superfrequency, RBThe second smooth radius is represented, B (x) is represented Fall into and be less than or equal to R relative to the distance of x pointsBIn the range of measuring point set, w2(η-x) represents B (x) weight coefficient, w2 (η-x)=1.Due to RB=0, directly use e2(x)=e1(x).With (3) formula the terrain elevation data ratio for removing low frequency and superfrequency Example obtains landform altitude high fdrequency component into the level with high-frequency components of statics with identical average amplitude value:
Wherein, ehf(x) landform altitude high fdrequency component, t are represented0(x) high frequency static correction value is represented.
Fig. 4 is shown removes low frequency and the terrain elevation data and landform of superfrequency according to an embodiment of the invention The schematic diagram of elevation high fdrequency component, wherein, grey dotted line e2(x) represent to remove low frequency and the ground of superfrequency by what the calculating of (2) formula obtained Shape altitude data, due to RB=0, grey dotted line e2(x) it is actually that the landform altitude number for removing low frequency obtained is calculated by (1) formula According to dotted line ehf(x) it is that the landform altitude high fdrequency component obtained is calculated by (3) formula.
Based on high frequency static correction value and landform altitude high fdrequency component, (4) formula is utilized, obtains preliminary similarity factor:
Wherein, s1(x) preliminary similarity factor is represented.
Preliminary similarity factor is weighted with (5) formula again, obtains high frequency static correction value and landform altitude high fdrequency component Similarity factor:
Wherein, scRepresent threshold value, sc=0.5.
Specifically, landform altitude high fdrequency component and similarity factor can be based on, utilizes (6) formula, it is related to obtain preliminary landform Component:
t1(x)=ehf(x)s2(x) (6)
Wherein, t1(x) preliminary landform correlated components are represented.
Preliminary landform correlated components are carried out with (7) formula again smoothly, landform correlated components can be obtained:
Wherein, t2(x) landform correlated components, R are representedCThe 3rd smooth radius is represented, C (x) represents to fall into relative to x points Distance is less than or equal to RCIn the range of measuring point set, w3(η-x) is weight coefficient, w3(η-x)=1.
By (8) formula, high frequency static correction value and landform correlated components are subtracted each other, obtain revised high frequency static correction value:
thf(x)=t0(x)-t2(x) (8)
Wherein, thf(x) revised high frequency static correction value is represented.
Fig. 5 shows the landform correlated components and revised high frequency static correction value according to an embodiment of the invention Schematic diagram, wherein, dotted line t2(x) it is that the landform correlated components obtained are calculated by (7) formula, R during calculatingC=300m, and grey point Line thf(x) it is that the revised high frequency static correction value obtained is calculated by (8) formula.
Fig. 6 shows the schematic diagram of the stacked section of the high frequency static correction value according to prior art.Fig. 7 is shown according to this The schematic diagram of the stacked section of the revised high frequency static correction value of one embodiment of invention.From based on first arrival statistical method The irrational static correction composition related to landform is extracted in obtained high frequency static correction value and is rejected, improves and is applied to ground The precision of the first arrival statistics high frequency static correction value of seismic data processing.
Using example 2
It is identical with application example 1 using the receiving point landform altitude and high frequency static correction value of example 2.
In this example, R=600m, R are takenB=300m, RC=0, w1(η-x)、w2(η-x) and w3(η-x) three weight functions Use the constant 1 unrelated with x relativeness with η, threshold value sc=0.5.
Based on terrain elevation data, high frequency static correction value and smooth radius, wherein, smooth radius can include:First is flat Sliding radius, the second smooth radius and the 3rd smooth radius.Low frequency component is rejected from terrain elevation data by (1) formula, obtained Remove the terrain elevation data of low frequency:
Wherein, e1(x) terrain elevation data of low frequency, e are removed in expression0(x) landform altitude at the x measuring points of plan-position is represented, R represents the first smooth radius, and A (x) represents to fall into the set for the measuring point being less than or equal in the range of R relative to the distance of x points, w1 (η-x) is weight coefficient, w1(η-x)=1.
Terrain elevation data to removing low frequency, (2) formula of use are smoothed, and reject the part of special high frequency, are obtained Remove low frequency and the terrain elevation data of superfrequency:
Wherein, e2(x) represent to remove low frequency and the terrain elevation data of superfrequency, RBThe second smooth radius is represented, B (x) is represented Fall into and be less than or equal to R relative to the distance of x pointsBIn the range of measuring point set, w2(η-x) is weight coefficient, w2(η-x)=1. With (3) formula going the terrain elevation data ratio of low frequency and superfrequency there is identical mean amplitude of tide into high-frequency components of statics The level of value, obtain landform altitude high fdrequency component:
Wherein, ehf(x) landform altitude high fdrequency component, t are represented0(x) high frequency static correction value is represented.
Fig. 8 is shown removes low frequency and the terrain elevation data and landform of superfrequency according to an embodiment of the invention The schematic diagram of elevation high fdrequency component, wherein, grey dotted line e2(x) represent to remove low frequency and the ground of superfrequency by what the calculating of (2) formula obtained Shape altitude data, dotted line ehf(x) it is that the landform altitude high fdrequency component obtained is calculated by (3) formula.
Based on high frequency static correction value and landform altitude high fdrequency component, (4) formula is utilized, obtains preliminary similarity factor:
Wherein, s1(x) preliminary similarity factor is represented.
Preliminary similarity factor is weighted with (5) formula again, obtains high frequency static correction value and landform altitude high fdrequency component Similarity factor:
Wherein, scRepresent threshold value, sc=0.5.
Based on landform altitude high fdrequency component and similarity factor, (6) formula is utilized, obtains preliminary landform correlated components:
t1(x)=ehf(x)s2(x) (6)
Wherein, t1(x) preliminary landform correlated components are represented.
Preliminary landform correlated components are carried out with (7) formula again smoothly, to obtain landform correlated components:
Wherein, t2(x) landform correlated components, R are representedCThe 3rd smooth radius is represented, C (x) represents to fall into relative to x points Distance is less than or equal to RCIn the range of measuring point set, w3(η-x) is weight coefficient, w3(η-x)=1.Due to RC=0, directly use t2(x)=t1(x)。
By (8) formula, high frequency static correction value and landform correlated components are subtracted each other, obtain revised high frequency static correction value:
thf(x)=t0(x)-t2(x) (8)
Wherein, thf(x) revised high frequency static correction value is represented.
Fig. 9 shows the landform correlated components and revised high frequency static correction value according to an embodiment of the invention Schematic diagram, wherein, dotted line t2(x) be by (7) formula calculate obtain landform correlated components, due to calculate when RC=0, stain Line t2(x) the preliminary landform correlated components actually obtained by the calculating of (6) formula, and grey dotted line thf(x) it is to be calculated by (8) formula The revised high frequency static correction value obtained.
Using example 3
It is identical with application example 1 using the receiving point landform altitude and high frequency static correction value of example 3.
In this example, R=600m, R are takenB=300m, RC=0, threshold value sc=0.5, w1(η-x)、w2(η-x) and w3 (η-x) three weight coefficients are using the function relevant with η and x relativeness shown in (9) formula:
Wherein, | η-x | the distance of η points and x points in plane is represented,R1=R, R2=RB、R3=RC
Based on terrain elevation data, high frequency static correction value and smooth radius, wherein, smooth radius includes:First smooth half Footpath, the second smooth radius and the 3rd smooth radius.Low frequency component is rejected from terrain elevation data by (1) formula, acquisition is gone low The terrain elevation data of frequency:
Wherein, e1(x) terrain elevation data of low frequency, e are removed in expression0(x) landform altitude at the x measuring points of plan-position is represented, R represents the first smooth radius, and A (x) represents to fall into the set for the measuring point being less than or equal in the range of R relative to the distance of x points, w1 (η-x) is weight coefficient.
Terrain elevation data to removing low frequency, (2) formula of use are smoothed, and reject the part of special high frequency, are obtained Remove low frequency and the terrain elevation data of superfrequency:
Wherein, e2(x) represent to remove low frequency and the terrain elevation data of superfrequency, RBThe second smooth radius is represented, B (x) is represented Fall into and be less than or equal to R relative to the distance of x pointsBIn the range of measuring point set, w2(η-x) is weight coefficient.Gone with (3) formula handle Low frequency and superfrequency terrain elevation data ratio obtain into the level with high-frequency components of statics with identical average amplitude value Landform altitude high fdrequency component:
Wherein, ehf(x) landform altitude high fdrequency component, t are represented0(x) high frequency static correction value is represented.
Figure 10 is shown removes low frequency and the terrain elevation data of superfrequency and ground according to an embodiment of the invention The schematic diagram of shape elevation high fdrequency component, wherein, grey dotted line e2(x) represent to be calculated by (2) formula obtain go low frequency and superfrequency Terrain elevation data, dotted line ehf(x) it is that the landform altitude high fdrequency component obtained is calculated by (3) formula.
Based on high frequency static correction value and landform altitude high fdrequency component, (4) formula is utilized, obtains preliminary similarity factor:
Wherein, s1(x) preliminary similarity factor is represented.
Preliminary similarity factor is weighted with (5) formula again, obtains high frequency static correction value and landform altitude high fdrequency component Similarity factor:
Wherein, scRepresent threshold value, sc=0.5.
Based on landform altitude high fdrequency component and similarity factor, (6) formula is utilized, obtains preliminary landform correlated components:
t1(x)=ehf(x)s2(x) (6)
Wherein, t1(x) preliminary landform correlated components are represented.
Preliminary landform correlated components are carried out with (7) formula again smoothly, to obtain landform correlated components:
Wherein, t2(x) landform correlated components, R are representedCThe 3rd smooth radius is represented, C (x) represents to fall into relative to x points Distance is less than or equal to RCIn the range of measuring point set, w3(η-x) represents C (x) weight coefficient.Due to RC=0, directly use t2 (x)=t1(x)。
By (8) formula, high frequency static correction value and landform correlated components are subtracted each other, obtain revised high frequency static correction value:
thf(x)=t0(x)-t2(x) (8)
Wherein, thf(x) revised high frequency static correction value is represented.
Figure 11 shows the landform correlated components and revised high frequency static correction according to an embodiment of the invention The schematic diagram of amount, wherein, dotted line t2(x) represent by the landform correlated components of (7) formula calculating acquisition, and grey dotted line thf(x) table Show and the revised high frequency static correction value obtained is calculated by (8) formula.
In summary, this method is the process of refinement of the high frequency static correction value to being counted based on first arrival, from based on first arrival The irrational static correction composition related to landform is extracted in the high frequency static correction value that statistical method obtains and is rejected, is improved First arrival applied to seism processing counts the precision of high frequency static correction value.
It will be understood by those skilled in the art that the purpose of the description to embodiments of the present invention is only for exemplarily above Illustrate the beneficial effect of embodiments of the present invention, be not intended to embodiments of the present invention being limited to given any show Example.
Embodiment 2
According to the embodiment of the present invention, there is provided a kind of high frequency static correction value update the system, the system can include: For based on terrain elevation data, high frequency static correction value and smooth radius, the unit of acquisition landform altitude high fdrequency component;For base In high frequency static correction value and the landform altitude high fdrequency component, the unit of similarity factor is obtained;For based on landform altitude high frequency Component and similarity factor, obtain the unit of landform correlated components;And for based on high frequency static correction value and landform correlated components, Obtain the unit of revised high frequency static correction value.
The embodiment obtains high accuracy by extracting irrational static correction composition related to landform and being rejected Revised high frequency static correction value.
In one example, obtaining landform altitude high fdrequency component can include:Based on terrain elevation data and smooth radius, Obtain the terrain elevation data for removing low frequency;Based on the terrain elevation data for removing low frequency, special radio-frequency head is removed by smoothing processing Point, low frequency and the terrain elevation data of superfrequency are removed in acquisition;And based on high frequency static correction value and remove low frequency and the ground of superfrequency Shape altitude data, obtain landform altitude high fdrequency component.
In one example, obtaining similarity factor can include:Based on high frequency static correction value and landform altitude high fdrequency component, Obtain preliminary similarity factor;And based on preliminary similarity factor, handled by weighting, obtain similarity factor.
In one example, obtaining landform correlated components can include:Based on landform altitude high fdrequency component and similarity factor, Obtain preliminary landform correlated components;And based on preliminary landform correlated components, by smoothing processing, obtain landform correlated components.
In one example, obtaining revised high frequency static correction value can include:High frequency static correction value is related to landform Component subtracts each other, and obtains revised high frequency static correction value.
It will be understood by those skilled in the art that the purpose of the description to embodiments of the present invention is only for exemplarily above Illustrate the beneficial effect of embodiments of the present invention, be not intended to embodiments of the present invention being limited to given any show Example.
It is described above the embodiments of the present invention, described above is exemplary, and non-exclusive, and It is also not necessarily limited to disclosed each embodiment.It is right in the case of without departing from the scope and spirit of illustrated each embodiment Many modifications and changes will be apparent from for those skilled in the art.The choosing of term used herein Select, it is intended to best explain the principle, practical application or the improvement to the technology in market of each embodiment, or make this technology Other those of ordinary skill in field are understood that each embodiment disclosed herein.

Claims (10)

1. a kind of high frequency static correction value modification method, including:
Based on terrain elevation data, high frequency static correction value and smooth radius, landform altitude high fdrequency component is obtained;
Based on the high frequency static correction value and the landform altitude high fdrequency component, similarity factor is obtained;
Based on the landform altitude high fdrequency component and the similarity factor, landform correlated components are obtained;And
Based on the high frequency static correction value and the landform correlated components, revised high frequency static correction value is obtained.
2. high frequency static correction value modification method according to claim 1, wherein, obtain the landform altitude high fdrequency component bag Include:
Based on the terrain elevation data and the smooth radius, the terrain elevation data for removing low frequency is obtained;
Based on the terrain elevation data for removing low frequency, special HFS is removed by smoothing processing, low frequency and spy are removed in acquisition The terrain elevation data of high frequency;And
Low frequency and the terrain elevation data of superfrequency are gone with described based on the high frequency static correction value, it is high to obtain the landform altitude Frequency component.
3. high frequency static correction value modification method according to claim 1, wherein, obtaining the similarity factor includes:
Based on the high frequency static correction value and the landform altitude high fdrequency component, preliminary similarity factor is obtained;And
Based on the preliminary similarity factor, handled by weighting, obtain the similarity factor.
4. high frequency static correction value modification method according to claim 1, wherein, obtaining the landform correlated components includes:
Based on the landform altitude high fdrequency component and the similarity factor, preliminary landform correlated components are obtained;And
Based on the preliminary landform correlated components, by smoothing processing, the landform correlated components are obtained.
5. high frequency static correction value modification method according to claim 1, wherein, obtain the revised high frequency static correction Amount includes:The high frequency static correction value subtracts each other with the landform correlated components, obtains the revised high frequency static correction value.
6. a kind of high frequency static correction value update the system, including:
For based on terrain elevation data, high frequency static correction value and smooth radius, the unit of acquisition landform altitude high fdrequency component;
For based on the high frequency static correction value and the landform altitude high fdrequency component, obtaining the unit of similarity factor;
For based on the landform altitude high fdrequency component and the similarity factor, the unit of acquisition landform correlated components;And
For based on the high frequency static correction value and the landform correlated components, obtaining the list of revised high frequency static correction value Member.
7. high frequency static correction value update the system according to claim 6, wherein, obtain the landform altitude high fdrequency component bag Include:
Based on the terrain elevation data and the smooth radius, the terrain elevation data for removing low frequency is obtained;
Based on the terrain elevation data for removing low frequency, special HFS is removed by smoothing processing, low frequency and spy are removed in acquisition The terrain elevation data of high frequency;And
Low frequency and the terrain elevation data of superfrequency are gone with described based on the high frequency static correction value, it is high to obtain the landform altitude Frequency component.
8. high frequency static correction value update the system according to claim 6, wherein, obtaining the similarity factor includes:
Based on the high frequency static correction value and the landform altitude high fdrequency component, preliminary similarity factor is obtained;And
Based on the preliminary similarity factor, handled by weighting, obtain the similarity factor.
9. high frequency static correction value update the system according to claim 6, wherein, obtaining the landform correlated components includes:
Based on the landform altitude high fdrequency component and the similarity factor, preliminary landform correlated components are obtained;And
Based on the preliminary landform correlated components, by smoothing processing, the landform correlated components are obtained.
10. high frequency static correction value update the system according to claim 6, wherein, obtain the revised quiet school of high frequency Positive quantity includes:The high frequency static correction value subtracts each other with the landform correlated components, obtains the revised high frequency static correction value.
CN201610862929.4A 2016-09-28 2016-09-28 High frequency static correction value modification method and system Active CN107870365B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610862929.4A CN107870365B (en) 2016-09-28 2016-09-28 High frequency static correction value modification method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610862929.4A CN107870365B (en) 2016-09-28 2016-09-28 High frequency static correction value modification method and system

Publications (2)

Publication Number Publication Date
CN107870365A true CN107870365A (en) 2018-04-03
CN107870365B CN107870365B (en) 2019-04-12

Family

ID=61761650

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610862929.4A Active CN107870365B (en) 2016-09-28 2016-09-28 High frequency static correction value modification method and system

Country Status (1)

Country Link
CN (1) CN107870365B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110927797A (en) * 2018-09-20 2020-03-27 中国石油化工股份有限公司 Static correction reference surface calculation method and system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6424920B1 (en) * 1999-09-17 2002-07-23 Konstantin Sergeevich Osypov Differential delay-time refraction tomography
CN104090301A (en) * 2014-07-21 2014-10-08 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Three-dimensional high-frequency static correction value obtaining method
CN104330826A (en) * 2014-07-31 2015-02-04 中国石油化工股份有限公司 A method for removing various noises under the condition of complex surface
CN105093280A (en) * 2014-05-15 2015-11-25 中国石油化工股份有限公司 Method of decomposing low frequency and high frequency components of surface layer model influencing earthquake data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6424920B1 (en) * 1999-09-17 2002-07-23 Konstantin Sergeevich Osypov Differential delay-time refraction tomography
CN105093280A (en) * 2014-05-15 2015-11-25 中国石油化工股份有限公司 Method of decomposing low frequency and high frequency components of surface layer model influencing earthquake data
CN104090301A (en) * 2014-07-21 2014-10-08 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Three-dimensional high-frequency static correction value obtaining method
CN104330826A (en) * 2014-07-31 2015-02-04 中国石油化工股份有限公司 A method for removing various noises under the condition of complex surface

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110927797A (en) * 2018-09-20 2020-03-27 中国石油化工股份有限公司 Static correction reference surface calculation method and system

Also Published As

Publication number Publication date
CN107870365B (en) 2019-04-12

Similar Documents

Publication Publication Date Title
CN104091064B (en) PS-DInSAR ground surface deformation measurement parameter estimation method based on optimal solution space search method
CN104392206B (en) A kind of image processing method of readings of pointer type meters automatic identification
CN107203981B (en) A kind of image defogging method based on fog concentration feature
CN105277923A (en) Single channel radar signal sorting method
CN104964685B (en) A kind of decision method of mobile phone athletic posture
CN104545887B (en) The recognition methods of artifact ecg wave form and device
CN106225786B (en) A kind of adaptive pedestrian navigation system zero-speed section detecting method
CN110174056A (en) A kind of object volume measurement method, device and mobile terminal
CN102737370B (en) Method and device for detecting image foreground
CN103487788B (en) The fast automatic extracting method of a kind of train pulse signal
CN103837895B (en) Matching preliminary wave obtains the method for short-period static corrections
CN107066998A (en) A kind of pointer-type circular single instrument board real-time identification method of utilization mobile device
CN105225243B (en) One kind can antimierophonic method for detecting image edge
CN110455316A (en) A kind of adaptive zero-speed section detecting method
CN112361977B (en) Linear distance measuring method based on weight distribution
CN106296680A (en) A kind of multiple features fusion high-resolution remote sensing image dividing method based on region
CN106327536A (en) Neck line measuring method based on sectional point clouds
CN105956544A (en) Remote sensing image road intersection extraction method based on structural index characteristic
CN107870365B (en) High frequency static correction value modification method and system
CN108921170A (en) A kind of effective picture noise detection and denoising method and system
CN108491927A (en) A kind of data processing method and device based on neural network
CN107256421A (en) A kind of rice wheat seed rapid counting method
CN108426521A (en) A kind of quality determining method and device of component
CN105405134A (en) Camera calibration mark center coordinate extraction method
CN105093280A (en) Method of decomposing low frequency and high frequency components of surface layer model influencing earthquake data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant