CN104618033B - A kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method - Google Patents
A kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method Download PDFInfo
- Publication number
- CN104618033B CN104618033B CN201510010552.5A CN201510010552A CN104618033B CN 104618033 B CN104618033 B CN 104618033B CN 201510010552 A CN201510010552 A CN 201510010552A CN 104618033 B CN104618033 B CN 104618033B
- Authority
- CN
- China
- Prior art keywords
- sequence
- gravitational cue
- self
- layer
- filter
- 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
- Noise Elimination (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method, after collection related data information, first multilayer decomposition is carried out to original gravitational cue, then the relevance degree of every one deck component and reference noise is calculated, determine that each layer component carries out the filter weights of self-adaptive harmonics detection according to relevance degree, finally complete the self-adaptive harmonics detection of each layer component, and use the filter result of every layer of component to complete the reconstruct of gravitational cue.The gravitational cue noise suppressed processing method that the present invention proposes overcomes filtering parameter in traditional filtering method can not the problem of self-adaptative adjustment; The present invention can retain effective gravitational cue, realizes high-precision gravitational cue and extracts.
Description
Technical field
The present invention relates to terrestrial gravitation Detection Techniques field, specifically a kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method, this invention effectively can suppress the noise in gravity measurement signal, extracts useful gravitational cue.
Background technology
Terrestrial gravitation data are grand strategy resources of country, have important effect at numerous areas such as national defense construction, space technology, resource detection and geophysics science.Gravitational cue process utilizes ocean or airborne gravitormeter to carry out gravitational cue collection on the spot exactly, and extracts useful gravitational cue by effective data processing method, for above-mentioned field provides gravimetric data support.In order to the measurement noises in gravitational cue effectively can be eliminated, the methods such as multiplex FIR or the IIR low pass filter of conventional process, Kalman filter.Although these methods achieve certain filter effect, but there is the problems such as the such as Parameter uncertainties such as cut-off frequency, exponent number, these uncertain factors bring adverse effect to the noise suppressed of filter, can not provide effective technical guarantee for gravity detection and data processing.
Summary of the invention
Goal of the invention: in order to overcome the deficiency existed in traditional filtering technique, the invention provides a kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method, can effectively eliminate gravity measurement noise and extract effective gravity signal.
Technical scheme: for achieving the above object, the technical solution used in the present invention is:
A kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method, comprises the steps:
(1) utilize ocean or airborne gravitormeter to measure according to effective survey line, obtain original gravitational cue sequence g (t)=[g (1) g (2) ... g (N)] of t, N is sequence length; Meanwhile, a high-precision accelerometer installed by gravimeter stabilized platform, measure interference and the noise of vertical direction, obtain noise reference signal sequence n (t) that length is N;
(2) initialization represents that the parameter k of number of repetition be the span of 1, k be 1 to K, K is arbitrary value in 3 ~ 10;
(3) constructing a pair length is N, and amplitude is identical, the Gaussian sequence x of carrier phase shift 180 °
(k)(t) and-x
(k)t (), then by x
(k)(t) and-x
(k)t () is superimposed with g (t) respectively obtains two new signal sequences
with
namely
Wherein, Gaussian sequence x
(k)(t) and-x
(k)t the amplitude of () selects the standard deviation of 0.2 ~ 0.5 times of original gravitational cue sequence g (t);
(4) to new signal sequence
with
carry out multilayer decomposition respectively, obtain two groups of decomposition result, often organize results set and be designated as respectively
with
wherein
M is Decomposition order;
(5) k adds 1, if k is less than or equal to K, then returns step (3), otherwise enters step (6);
(6) set of computations
with
the mean value of middle respective components, obtains new set G (t)={ g
1(t) g
2(t) ... g
m(t) }, wherein:
(7) each layer component g obtained in calculation procedure (6)
ithe relevance degree of noise reference signal sequence n (t) obtained in (t) and step (1), the relatedness computation formula of i-th layer of component is as follows:
Wherein || g
i(t) ||
2represent sequence g
ithe 2-norm of (t), || n (t) ||
2represent the 2-norm of sequence n (t);
(8) calculate each layer component g obtained in step (6) according to the relevance degree of step (7)
it filter weights R that () processes further
i, filter weights calculates as follows by formula:
(9) by R
ias the filter weights of self-adaptive harmonics detection process, to each layer component g obtained in step (6)
it () carries out self-adaptive harmonics detection process, obtain the filtering output sequence of each layer
(10) the filtering output sequence of each layer obtained by step (9)
carry out the reconstruct of gravitational cue sequence, the effective gravity burst after the measurement noises that can be eliminated
that is:
Multilayer decomposition in described step (4) specifically comprises the steps:
(41) new signal sequence is asked for
all extreme points, determine all maximum points and the minimum point of this sequence;
(42) all maximum points, minimum point are coupled together with cubic spline curve respectively, obtain upper and lower envelope
with
sequence between upper and lower envelope
(43) the average line of upper and lower envelope is calculated
and obtain
(44) h is judged
1t whether () meet following two conditions:
(441) signal zero crossing number is equal with limit number or be more or less the same in 1;
(442) average of the envelope formed respectively by maximum point and minimum point is 0;
If do not meet (441) and (442) two conditions, then use
substitute new signal sequence
repeat step (41) ~ (44), until meet (441) and (442) two conditions, namely obtain first decomposed component, be expressed as
(45) by sequence
repeat step (41) ~ (44) as new signal sequence, obtain each layer decomposed component
until m layer has decomposed, gathered
(46) to new signal sequence
carry out multilayer decomposition according to step (41) ~ (45), obtain multilayer and decompose set
Self-adaptive harmonics detection process in described step (9) specifically comprises the steps:
(91) designated length is the shape filtering structural element b of L, according to the feature of gravimeter signal, can choose sinusoidal configuration element, such as b={0 7.7071 10 7.7071 0}, now L=5;
(92) to g
it () employing form open-close filter and form are closed-Kai filter and are carried out data filtering, wherein:
Form open-close filter is defined as:
OC(g
i(t))=g
i(t)ο(R
ib)·(R
ib),
Form is closed-Kai filter and is defined as:
CO(g
i(t))=g
i(t)·(R
ib)ο(R
ib),
Open-close operation definition is:
Close-opening operation is defined as:
Operator Θ is defined as:
Operator
be defined as:
(93) to open-close filter results with close-Kai filter results and average, obtain self-adaptive harmonics detection and export
(94) (92) ~ (93) are repeated, until g
i(t), i=1,2 ..., the whole filtering of m is complete, obtains the filtering output sequence of each layer
i=1,2 ..., m.
Beneficial effect: multilayer self-adaptive harmonics detection gravitational cue noise suppressing method provided by the invention, has following advantage relative to prior art: 1, overcoming filtering parameter in traditional filtering method can not the problem of self-adaptative adjustment; 2, effective gravitational cue can be retained, realize high-precision gravitational cue and extract.
Accompanying drawing explanation
Fig. 1 is the realization flow figure of the inventive method;
Fig. 2 is the implementation result figure of the inventive method.
Embodiment
Below in conjunction with drawings and Examples, technical solutions according to the invention are further elaborated.
A kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method, comprises the steps:
(1) utilize ocean or airborne gravitormeter to measure according to effective survey line, obtain original gravitational cue sequence g (t)=[g (1) g (2) ... g (N)] of t, N is sequence length; Meanwhile, a high-precision accelerometer installed by gravimeter stabilized platform, measure interference and the noise of vertical direction, obtain noise reference signal sequence n (t) that length is N;
(2) initialization represents that the parameter k of number of repetition be the span of 1, k be 1 to K, K is arbitrary value in 3 ~ 10;
(3) constructing a pair length is N, and amplitude is identical, the Gaussian sequence x of carrier phase shift 180 °
(k)(t) and-x
(k)t (), then by x
(k)(t) and-x
(k)t () is superimposed with g (t) respectively obtains two new signal sequences
with
namely
Wherein, Gaussian sequence x
(k)(t) and-x
(k)t the amplitude of () selects the standard deviation of 0.2 ~ 0.5 times of original gravitational cue sequence g (t);
(4) to new signal sequence
with
carry out multilayer decomposition respectively, obtain two groups of decomposition result, often organize results set and be designated as respectively
with
wherein
M is Decomposition order;
Multilayer decomposition specifically comprises the steps:
(41) new signal sequence g is asked for
1 (k)t all extreme points of (), determine all maximum points and the minimum point of this sequence;
(42) all maximum points, minimum point are coupled together with cubic spline curve respectively, obtain upper and lower envelope
with
sequence between upper and lower envelope
(43) the average line of upper and lower envelope is calculated
and obtain
(44) h is judged
1t whether () meet following two conditions:
(441) signal zero crossing number is equal with limit number or be more or less the same in 1;
(442) average of the envelope formed respectively by maximum point and minimum point is 0;
If do not meet (441) and (442) two conditions, then use
substitute new signal sequence
repeat step (41) ~ (44), until meet (441) and (442) two conditions, namely obtain first decomposed component, be expressed as
(45) by sequence
repeat step (41) ~ (44) as new signal sequence, obtain each layer decomposed component
until m layer has decomposed, gathered
(46) to new signal sequence
carry out multilayer decomposition according to step (41) ~ (45), obtain multilayer and decompose set
(5) k adds 1, if k is less than or equal to K, then returns step (3), otherwise enters step (6);
(6) set of computations
with
the mean value of middle respective components, obtains new set G (t)={ g
1(t) g
2(t) ... g
m(t) }, wherein:
(7) each layer component g obtained in calculation procedure (6)
ithe relevance degree of noise reference signal sequence n (t) obtained in (t) and step (1), the relatedness computation formula of i-th layer of component is as follows:
Wherein || g
i(t) ||
2represent sequence g
ithe 2-norm of (t), || n (t) ||
2represent the 2-norm of sequence n (t);
(8) calculate each layer component g obtained in step (6) according to the relevance degree of step (7)
it filter weights R that () processes further
i, filter weights calculates as follows by formula:
(9) by R
ias the filter weights of self-adaptive harmonics detection process, to each layer component g obtained in step (6)
it () carries out self-adaptive harmonics detection process, obtain the filtering output sequence of each layer
Self-adaptive harmonics detection process specifically comprises the steps:
(91) designated length is the shape filtering structural element b of L, according to the feature of gravimeter signal, can choose sinusoidal configuration element, such as b={07.7071107.70710}, now L=5;
(92) to g
it () employing form open-close filter and form are closed-Kai filter and are carried out data filtering, wherein:
Form open-close filter is defined as:
OC(g
i(t))=g
i(t)ο(R
ib)·(R
ib),
Form is closed-Kai filter and is defined as:
CO(g
i(t))=g
i(t)·(R
ib)ο(R
ib),
Open-close operation definition is:
Close-opening operation is defined as:
Operator Θ is defined as:
Operator
be defined as:
(93) to open-close filter results with close-Kai filter results and average, obtain self-adaptive harmonics detection and export
(94) (92) ~ (93) are repeated, until g
i(t), i=1,2 ..., the whole filtering of m is complete, obtains the filtering output sequence of each layer
i=1,2 ..., m;
(10) the filtering output sequence of each layer obtained by step (9)
carry out the reconstruct of gravitational cue sequence, the effective gravity burst after the measurement noises that can be eliminated
that is:
Performance evaluation
The multilayer self-adaptive harmonics detection gravitational cue noise suppressed processing method that the present invention proposes overcomes filtering parameter in traditional filtering method can not the problem of self-adaptative adjustment, each noise like can be suppressed preferably, retain effective gravitational cue, realize high-precision gravitational cue and extract.
Accompanying drawing 2 is gravitational cue noise suppressing method implementation result figure.In figure, solid line represents gravity measurement signal before treatment, and dotted line represents the output gravitational cue after noise suppressing method of the present invention.As can be seen from accompanying drawing 2, high-frequency noise is effectively suppressed, and gravitational cue and variation tendency thereof are extracted efficiently, and reaches the object of noise suppressed.
The above is only main embodiment of the present invention; be noted that for those skilled in the art; under the premise without departing from the principles of the invention, can also make some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.
Claims (1)
1. a multilayer self-adaptive harmonics detection gravitational cue noise suppressing method, is characterized in that: comprise the steps:
(1) utilize ocean or airborne gravitormeter to measure according to effective survey line, obtain original gravitational cue sequence g (t)=[g (1) g (2) of t ... g (N)], N is sequence length; Meanwhile, a high-precision accelerometer installed by gravimeter stabilized platform, measure interference and the noise of vertical direction, obtain noise reference signal sequence n (t) that length is N;
(2) initialization represents that the parameter k of number of repetition be the span of 1, k be 1 to K, K is arbitrary value in 3 ~ 10;
(3) constructing a pair length is N, and amplitude is identical, the Gaussian sequence x of carrier phase shift 180 °
(k)(t) and-x
(k)t (), then by x
(k)(t) and-x
(k)t () is superimposed with g (t) respectively obtains two new signal sequences
with
namely
Wherein, Gaussian sequence x
(k)(t) and-x
(k)t the amplitude of () selects the standard deviation of 0.2 ~ 0.5 times of original gravitational cue sequence g (t);
(4) to new signal sequence
with
carry out multilayer decomposition respectively, obtain two groups of decomposition result, often organize results set and be designated as respectively
with
wherein
M is Decomposition order; Multilayer decomposition specifically comprises the steps:
(41) new signal sequence is asked for
all extreme points, determine all maximum points and the minimum point of this sequence;
(42) all maximum points, minimum point are coupled together with cubic spline curve respectively, obtain upper and lower envelope
with
sequence between upper and lower envelope
(43) the average line of upper and lower envelope is calculated
and obtain
(44) h is judged
1t whether () meet following two conditions:
(441) signal zero crossing number is equal with limit number or be more or less the same in 1;
(442) average of the envelope formed respectively by maximum point and minimum point is 0;
If do not meet (441) and (442) two conditions, then use
substitute new signal sequence
repeat step (41) ~ (44), until meet (441) and (442) two conditions, namely obtain first decomposed component, be expressed as
(45) by sequence
repeat step (41) ~ (44) as new signal sequence, obtain each layer decomposed component
until m layer has decomposed, gathered
(46) to new signal sequence
carry out multilayer decomposition according to step (41) ~ (45), obtain multilayer and decompose set
(5) k adds 1, if k is less than or equal to K, then returns step (3), otherwise enters step (6);
(6) set of computations
with
the mean value of middle respective components, obtains new set G (t)={ g
1(t) g
2(t) ... g
m(t) }, wherein:
(7) each layer component g obtained in calculation procedure (6)
ithe relevance degree of noise reference signal sequence n (t) obtained in (t) and step (1), the relatedness computation formula of i-th layer of component is as follows:
Wherein || g
i(t) ||
2represent sequence g
ithe 2-norm of (t), || n (t) ||
2represent the 2-norm of sequence n (t);
(8) calculate each layer component g obtained in step (6) according to the relevance degree of step (7)
it filter weights R that () processes further
i, filter weights calculates as follows by formula:
(9) by R
ias the filter weights of self-adaptive harmonics detection process, to each layer component g obtained in step (6)
it () carries out self-adaptive harmonics detection process, obtain the filtering output sequence of each layer
self-adaptive harmonics detection process specifically comprises the steps:
(91) designated length is the shape filtering structural element b of L;
(92) to g
it () employing form open-close filter and form are closed-Kai filter and are carried out data filtering, wherein:
Form open-close filter is defined as:
OC(g
i(t))=g
i(t)ο(R
ib)·(R
ib),
Form is closed-Kai filter and is defined as:
CO(g
i(t))=g
i(t)·(R
ib)ο(R
ib),
Open-close operation definition is:
g
i(t)ο(R
ib)·(R
ib)=g
i(t)Θ(R
ib)⊕(R
ib),
Close-opening operation is defined as:
g
i(t)·(R
ib)ο(R
ib)=g
i(t)⊕(R
ib)Θ(R
ib),
Operator Θ is defined as:
Operator ⊕ is defined as:
(93) to open-close filter results with close-Kai filter results and average, obtain self-adaptive harmonics detection and export
(94) (92) ~ (93) are repeated, until g
i(t), i=1,2 ..., the whole filtering of m is complete, obtains the filtering output sequence of each layer
i=1,2 ..., m
(10) the filtering output sequence of each layer obtained by step (9)
carry out the reconstruct of gravitational cue sequence, the effective gravity burst after the measurement noises that can be eliminated
that is:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510010552.5A CN104618033B (en) | 2015-01-08 | 2015-01-08 | A kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510010552.5A CN104618033B (en) | 2015-01-08 | 2015-01-08 | A kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104618033A CN104618033A (en) | 2015-05-13 |
CN104618033B true CN104618033B (en) | 2015-09-09 |
Family
ID=53152340
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510010552.5A Active CN104618033B (en) | 2015-01-08 | 2015-01-08 | A kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104618033B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105823492B (en) * | 2016-03-18 | 2018-08-21 | 北京卫星环境工程研究所 | Weak target signal extracting method in a kind of interference of ocean current |
CN108761565B (en) * | 2018-03-30 | 2019-11-12 | 中国国土资源航空物探遥感中心 | Method for evaluating quality, device, equipment and the storage medium of airborne gravity measurement |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4538251A (en) * | 1983-03-04 | 1985-08-27 | Digicon, Inc. | Marine seismic streamer cable for providing selectable detector array connections |
US4882713A (en) * | 1988-09-06 | 1989-11-21 | Exxon Production Research Company | Method for noise suppression in the stacking of seismic traces |
CN102338890A (en) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | Circular window band-pass amplitude-preserving filtering data processing method in geophysical exploration |
CN102338884A (en) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | Elliptic window direction band-pass amplitude-preserving filtering data processing method in geophysical prospecting |
CN102368094A (en) * | 2010-10-22 | 2012-03-07 | 中国石油天然气股份有限公司 | Rectangular window band-pass amplitude-preserving filtering data processing method in geophysical exploration |
WO2013177340A1 (en) * | 2012-05-22 | 2013-11-28 | Nxt Energy Solutions, Inc. | Gravity transducer and application to hydrocarbon exploration |
CN103472493A (en) * | 2013-09-06 | 2013-12-25 | 北京航天控制仪器研究所 | Weak signal detecting circuit used for gravity gradient measurement |
-
2015
- 2015-01-08 CN CN201510010552.5A patent/CN104618033B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4538251A (en) * | 1983-03-04 | 1985-08-27 | Digicon, Inc. | Marine seismic streamer cable for providing selectable detector array connections |
US4882713A (en) * | 1988-09-06 | 1989-11-21 | Exxon Production Research Company | Method for noise suppression in the stacking of seismic traces |
CN102338890A (en) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | Circular window band-pass amplitude-preserving filtering data processing method in geophysical exploration |
CN102338884A (en) * | 2010-10-22 | 2012-02-01 | 中国石油天然气股份有限公司 | Elliptic window direction band-pass amplitude-preserving filtering data processing method in geophysical prospecting |
CN102368094A (en) * | 2010-10-22 | 2012-03-07 | 中国石油天然气股份有限公司 | Rectangular window band-pass amplitude-preserving filtering data processing method in geophysical exploration |
WO2013177340A1 (en) * | 2012-05-22 | 2013-11-28 | Nxt Energy Solutions, Inc. | Gravity transducer and application to hydrocarbon exploration |
CN103472493A (en) * | 2013-09-06 | 2013-12-25 | 北京航天控制仪器研究所 | Weak signal detecting circuit used for gravity gradient measurement |
Non-Patent Citations (2)
Title |
---|
自适应Kalman和零相移滤波算法在重力信号处理中的对比;罗骋等人;《中国惯性技术学报》;20110601;全文 * |
自适应形态小波滤波在重力仪信号处理中的应用;赵立业等人;《2010年江苏省仪器仪表学会学术年会论文集》;20100101;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104618033A (en) | 2015-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102749648B (en) | Frequency division matched filtering method for improving offshore seismic data resolution by utilizing seismic focuses of different depths | |
RU2500002C2 (en) | Method of collecting data through three-dimensional regular electromagnetic array of small integration cells | |
CN104793253B (en) | Aviation electromagnetic data de-noising method based on mathematical morphology | |
CN106054244B (en) | The LPF of window multiple dimensioned full waveform inversion method when blocking | |
Luo et al. | Reduction to the pole at the geomagnetic equator | |
CN103995288B (en) | Gauss beam prestack depth migration method and device | |
CN106786561B (en) | A kind of Low-frequency Oscillation Modal Parameters discrimination method based on adaptive Kalman filter | |
CN102879817B (en) | Control method for acquiring underground crack information on basis of ground earthquake data | |
CN103630933A (en) | Nonlinear optimization based time-space domain staggered grid finite difference method and device | |
CN104122588A (en) | Spectral decomposition based post-stack seismic data resolution ratio increasing method | |
CN108387933A (en) | Method, device and system for determining formation quality factor | |
CN102854533A (en) | Wave field separation principle based denoising method for increasing signal to noise ratio of seismic data | |
Houtani et al. | Generation of a spatially periodic directional wave field in a rectangular wave basin based on higher-order spectral simulation | |
CN105676291B (en) | Multiple wave matching attenuation method based on optimized phase axis tracking | |
CN103674001A (en) | Fiber gyroscope denoising method based on enhanced adaptive time-frequency peak value filtration | |
CN104597488B (en) | Optimum design method of finite difference template of non-equiangular long-grid wave equation | |
CN104316958A (en) | Coherent processing method for identifying different scales of formation fractures | |
CN104614769A (en) | Beam-forming filtering method for suppressing seismic surface waves | |
CN104749628A (en) | Absorbing boundary reflection method based on dispersal viscosity wave equation | |
CN116520419B (en) | Hot fluid crack channel identification method | |
CN104618033B (en) | A kind of multilayer self-adaptive harmonics detection gravitational cue noise suppressing method | |
CN104330826A (en) | A method for removing various noises under the condition of complex surface | |
CN104297792B (en) | The phased inversion method of water channel reservoir is stacked on a kind of fan | |
CN108279437A (en) | Variable density ACOUSTIC WAVE EQUATION time higher order accuracy staggering mesh finite-difference method | |
CN106950600A (en) | A kind of minimizing technology of near surface scattering surface ripple |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |