CN105974488A - De-noising method and device for time shifting gravity increment - Google Patents
De-noising method and device for time shifting gravity increment Download PDFInfo
- Publication number
- CN105974488A CN105974488A CN201610517205.6A CN201610517205A CN105974488A CN 105974488 A CN105974488 A CN 105974488A CN 201610517205 A CN201610517205 A CN 201610517205A CN 105974488 A CN105974488 A CN 105974488A
- Authority
- CN
- China
- Prior art keywords
- gravity
- measuring point
- denoising
- increment
- time shift
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
Abstract
The invention provides a de-noising method and device for a time shifting gravity increment. The method comprises the step of weighting the gravity variances of all work area measured points measured multiple times within a first time period and the gravity variances and time shifting gravity increments measured multiple times within a second time period to obtain a de-noised time shifting gravity increment of the current measured point. By adopting the weighting method, the problem that the relation between the data obtained by observing time shifting microgravity multiple times and the data obtained by single observation and multiple readings is not considered during de-noising in the prior art is effectively solved, and the observation error of the time shifting gravity increment, namely observation noise, is removed, so that the goodness of fit between the finally obtained data and the practical geologic status is higher.
Description
Technical field
The present invention relates to oil exploration technology field, particularly to denoising method and the device of a kind of time shift gravity increment.
Background technology
During production of hydrocarbons, the oil gas material of underground can redistribute, thus causes the gravitational field generation office of underground
The change in portion.The measurement carried out for the change of Subsurface Gravity field is referred to as time shift microgravimetry, to earth Microgravity
The data that the localized variation produced is carried out obtained by repetition measurement are referred to as time shift microgravity data.Therefore, time shift is carried out
Need during microgravimetry same work area is repeatedly observed, identical measuring point is repeatedly observed repeatedly reading, from
And all there is time shift microgravity increment in the measuring point in work area.
At present, microgravity exploration data are carried out denoising, i.e. the indivedual kick points in survey data are put down
The method that sliding filtering is used mainly has: the denoising method of upward continuation and the denoising side of low-pass filtering
Method.
But above two method does not consider that carrying out time shift microgravity repeatedly observes the data and single observation obtained
The repeatedly contact between the data obtained by reading, thus it is not particularly suited for the denoising feelings of time shift microgravity increment
Condition.
For the problems referred to above, effective solution is the most not yet proposed.
Summary of the invention
Embodiments provide denoising method and the device of a kind of time shift gravity increment, to solve prior art is adopted
When time shift gravimetric data being carried out denoising by conventional denoising method, it is impossible to consider that time shift microgravity is repeatedly seen
The data recorded and single observe repeatedly the problem contacted between the data obtained by reading.
Embodiments provide the denoising method of a kind of time shift gravimetric data, including: each measuring point in acquisition work area
Time shift gravity increment, wherein, described time shift gravity increment be measuring point in first time period the gravity of repetitive measurement put down
Average and within the second time period difference between the gravity meansigma methods of repetitive measurement；Obtain each measuring point described first
The gravity variance of repetitive measurement in the time period and the gravity variance of repetitive measurement within the second time period；To described each
In measuring point, each measuring point uses following operation to obtain in each measuring point described time shift gravity increment after each measuring point denoising:
Select current measuring point as central point, will be limited as border using default denoising interval centered by described central point
In the range of all measuring points in first time period the gravity variance of repetitive measurement, within the second time period repetitive measurement
Gravity variance and time shift gravity increment be weighted process, obtain the described current measuring point after denoising time shift gravity increase
Amount.
In one embodiment, in obtaining each measuring point described after time shift gravity increment after each measuring point denoising,
Also include: the time shift gravity increment after measuring point denoising each in each measuring point described is increased as the time shift gravity of non-denoising
Amount, repeats following denoising and operates until reaching default number of repetition: choosing each measuring point in each measuring point described
Select current measuring point as central point, will be limited as border using default denoising interval centered by described central point
In the range of all measuring points in first time period the gravity variance of repetitive measurement, within the second time period repetitive measurement
The time shift gravity increment of gravity variance and non-denoising is weighted processing, and obtains the time shift of the described current measuring point after denoising
Gravity increment.
In one embodiment, described default number of repetition is not less than 10 times.
In one embodiment, the time shift gravity increasing of the described current measuring point after denoising it is calculated according to below equation
Amount:
When A < (D/2) and B < time (D/2), M=N；
When (D/2) < A < (3D/4) and (D/2) < B < time (3D/4), M=N × a+P × b；
When (D/2)<A<(3D/4) and B>(3D/4), (D/2)<when B<(3D/4) and A>(3D/4) meets alternative one,
M=N × c+P × d；
As A > (3D/4) and B > (3D/4) time, M=N × e+P × f；
Wherein, A represents using described default denoising interval as all measuring points in the limited range of border first
The gravity variance of repetitive measurement in time period, B represents the gravity variance of repetitive measurement within the second time period, and D represents
Interval as all measuring points repetitive measurement in first time period in the limited range of border using described default denoising
Gravity variance and within the second time period the maximum in the gravity variance of repetitive measurement, M represents the institute after denoising
Stating the time shift gravity increment of current measuring point, N represents described time shift gravity increment, and P represents interval with described default denoising
As the meansigma methods of the time shift gravity increment of all measuring points in the limited range of border, a, b, c, d, e, f
All represent the weight coefficient between 0% to 100%.
In one embodiment, the time shift gravity increment of each measuring point in acquisition work area, including: to moving weight time described
Power increment carries out gridding process；Obtain described each measuring point repetitive measurement in first time period gravity variance and
The gravity variance of repetitive measurement within the second time period, including: each measuring point described is repeatedly surveyed in first time period
Amount gravity variance and within the second time period the gravity variance of repetitive measurement carry out gridding process.
In one embodiment, described gridding processing method includes: Kriging regression method.
In one embodiment, described default denoising interval includes: the matrix of n × n, and wherein, described n represents institute
State the measuring point number in default denoising interval.
The embodiment of the present invention additionally provides the denoising device of a kind of time shift gravity increment, including: meansigma methods acquisition module,
The time shift gravity increment of each measuring point in obtaining work area, wherein, described time shift gravity increment is that measuring point is when first
Between in section repetitive measurement gravity meansigma methods and within the second time period difference between the gravity meansigma methods of repetitive measurement；
Variance acquisition module, for obtaining the gravity variance of described each measuring point repetitive measurement in first time period and the
The gravity variance of repetitive measurement in two time periods；Weighting processing module, for adopting each measuring point in each measuring point described
Time shift gravity increment after each measuring point denoising is obtained in each measuring point described: select current measuring point conduct with following operation
Central point, will be interval as all surveys in the limited range of border using default denoising centered by described central point
Point is the gravity variance of repetitive measurement, the gravity variance of repetitive measurement and time shift within the second time period in first time period
Gravity increment is weighted processing, and obtains the time shift gravity increment of the described current measuring point after denoising.
In one embodiment, the denoising device of described time shift gravity increment also includes: replicated blocks, for obtaining
In each measuring point described after time shift gravity increment after each measuring point denoising, each measuring point in each measuring point described is gone
Each measuring point in each measuring point described, as the time shift gravity increment of non-denoising, is repeated to hold by the time shift gravity increment after making an uproar
Row following denoising operation is until reaching default number of repetition: single replicated blocks, be used for selecting current measuring point as in
Heart point, will be interval as all measuring points in the limited range of border using default denoising centered by described central point
The gravity variance of repetitive measurement, the gravity variance of repetitive measurement and non-denoising within the second time period in first time period
Time shift gravity increment be weighted process, obtain the time shift gravity increment of the described current measuring point after denoising.
In one embodiment, described default number of repetition is not less than 10 times.
In the above-described embodiments, during time shift gravity increment is carried out denoising, by all measuring points first
The gravity variance of repetitive measurement in time period, within the second time period the gravity variance of repetitive measurement and time shift gravity increment
Being weighted processing, obtain the time shift gravity increment of the current measuring point after denoising, the method i.e. using weighting to process is effective
Solve time prior art carries out denoising, do not account for carrying out time shift microgravity repeatedly observe the data that obtain with
Single observes repeatedly the problem of the contact between the data obtained by reading, eliminates the observation of time shift gravity increment by mistake
Difference, i.e. observation noise so that the data finally given are higher with the actual geological condition goodness of fit.
Accompanying drawing explanation
Accompanying drawing described herein is used for providing a further understanding of the present invention, constitutes the part of the application, not
Constitute limitation of the invention.In the accompanying drawings:
Fig. 1 is the flow chart of the denoising method of a kind of time shift gravity increment of the embodiment of the present invention；
Fig. 2 is a kind of structured flowchart of the denoising device of a kind of time shift gravity increment of the embodiment of the present invention.
Detailed description of the invention
For making the object, technical solutions and advantages of the present invention clearer, below in conjunction with embodiment and accompanying drawing,
The present invention is described in further details.Here, the exemplary embodiment of the present invention and explanation thereof are used for explaining this
Bright, but not as a limitation of the invention.
For using conventional denoising method that time shift gravimetric data is carried out denoising time institute produced problem, send out
A person of good sense proposes the denoising method of a kind of time shift gravity increment, with solve use upward continuation denoising method and
When the denoising method of low-pass filtering carries out denoising, consideration does not carries out time shift microgravity and repeatedly observes and obtaining
Data and single observe repeatedly the problem contacted between the data obtained by reading.Specifically, this one time shift weight
The denoising method of power increment can be as it is shown in figure 1, comprise the following steps:
Step 101: the time shift gravity increment of each measuring point in acquisition work area, wherein, described time shift gravity increment is for surveying
The point gravity meansigma methods of repetitive measurement in first time period and the gravity meansigma methods of repetitive measurement within the second time period
Between difference；
Step 102: obtain the gravity variance of described each measuring point repetitive measurement in first time period and when second
Between the gravity variance of repetitive measurement in section；
Wherein, relevant with gravity in the present embodiment data all refer in particular to microgravity.Microgravimetry refers specifically to, micro-
The measurement technology and method of weak gravity anomaly.Microgravity is also referred to as microgal level gravity, and milligal is gravimetric conventional list
Position, microgal is the one thousandth of milligal.
When carrying out gravity measurement during reality is applied, need to use multiple stage instrument that same measuring point is carried out repeatedly
Measure.Wherein, gravity meansigma methods refers to the meansigma methods of multiple stage instrument repeatedly gravity measurement result.Similarly, gravity
Variance refers to the variance yields of multiple stage instrument repeatedly gravity measurement result.
Wherein, above-mentioned first time period and the second time period are two independent time periods, do not have therebetween overlap.
Time shift gravity increment, the first gravity variance and the second gravity side to each measuring point that above-mentioned observation obtains respectively
Difference carries out gridding process.In the present embodiment, can use that above-mentioned observation obtains by Kriging regression method respectively is each
The time shift gravity increment of individual measuring point, the first gravity variance and the second gravity variance carry out gridding process, such that it is able to
Obtain time shift gravity increment grid file, the first gravity variance grid file and the second gravity variance grid file.
Step 103: each measuring point in each measuring point described uses following operation obtain each survey in each measuring point described
Point denoising after time shift gravity increment: select current measuring point as central point, will be centered by described central point, with in advance
If denoising is interval as all measuring points gravity side of repetitive measurement in first time period in the limited range of border
Gravity variance poor, repetitive measurement within the second time period and time shift gravity increment are weighted processing, after obtaining denoising
The time shift gravity increment of current measuring point.
Further, in obtaining each measuring point after time shift gravity increment after each measuring point denoising, can be by each
In measuring point, time shift gravity increment after each measuring point denoising is as the time shift gravity increment of non-denoising, to every in each measuring point
Individual measuring point repeat following denoising operate until reaching default number of repetition: select current measuring point as central point,
Will centered by this central point, using default denoising interval as all measuring points in the limited range of border when first
Between the gravity variance of repetitive measurement in section, within the second time period the gravity variance of repetitive measurement and the time shift weight of non-denoising
Power increment is weighted processing, and obtains the time shift gravity increment of the current measuring point after denoising.
Wherein, above-mentioned default denoising interval can include but not limited to: the matrix of n × n, and wherein, described n represents
Preset the measuring point number in denoising interval；Meanwhile, above-mentioned default number of repetition can be set as not less than 10 times.
In the present embodiment, the time shift gravity increment of the current measuring point after denoising can be calculated according to below equation:
When A < (D/2) and B < time (D/2), M=N；
When (D/2) < A < (3D/4) and (D/2) < B < time (3D/4), M=N × a+P × b；
When (D/2)<A<(3D/4) and B>(3D/4), (D/2)<when B<(3D/4) and A>(3D/4) meets alternative one,
M=N × c+P × d；
As A > (3D/4) and B > (3D/4) time, M=N × e+P × f；
Wherein, A represented using default denoising interval as all measuring points in the limited range of border in the very first time
The gravity variance of repetitive measurement in section, B represents the gravity variance of repetitive measurement within the second time period, and D represents with in advance
If denoising is interval as all measuring points gravity side of repetitive measurement in first time period in the limited range of border
Differing from and the maximum in the gravity variance of repetitive measurement within the second time period, M represents the current measuring point after denoising
Time shift gravity increment, N represents time shift gravity increment, and P represents interval as border limited range using default denoising
The meansigma methods of the time shift gravity increment of interior all measuring points, a, b, c, d, e, f all represent between 0% to 100%
Weight coefficient.
Such as: a could be arranged to 10%, now, accordingly, b could be arranged to 90%；C could be arranged to 21%,
Now, accordingly, d could be arranged to 79%；E could be arranged to 41%, and now, accordingly, f could be arranged to
59%.It should be noted that the above-mentioned value to a, b, c, d, e, f is intended merely to better illustrate the application,
Can also take other numerical value, this is not construed as limiting by the application.But, a, b, c, d, e, f need below satisfied
Require: a+b=100%, c+d=100%, e+f=100%.
The denoising process of the present embodiment is described below as a example by a certain measuring point is carried out time shift gravity increment denoising, can
To include: first, determine the current measuring point carrying out denoising the number of repetition determining denoising；Again with this survey
Centered by Dian, using default denoising interval as all measuring points in the limited range of border in first time period repeatedly
Measure gravity variance, within the second time period, the gravity variance of repetitive measurement and the time shift gravity increment of non-denoising are carried out
Weighting processes, and obtains the time shift gravity increment of the current measuring point after denoising；Then, the current survey after the denoising that will obtain
The time shift gravity increment of point, as the time shift gravity increment of non-denoising, repeats the process of above-mentioned denoising, until it reaches
The number of repetition of the above-mentioned denoising determined；Finally, the time shift gravity increment of the current measuring point after denoising is obtained.Right
Time shift gravity increment grid file, the first gravity variance grid file and the second gravity variance grid file measured
In all measuring points be performed both by above-mentioned denoising process, the time shift gravity increment grid file after denoising can be obtained.
Below with time shift gravity increment grid file, first gravity variance grid file and the second gravity side in somewhere
The denoising method of a kind of time shift gravity increment that the present embodiment is provided as a example by carrying out denoising by difference grid file is entered
Row explanation, it is important to note, however, that this specific embodiment is merely to be better described the present invention, it is right to be not intended that
The improper restriction of the present invention.
In the present embodiment, the matrix of 5 × 5, i.e. n=5 can be selected, as carrying out time shift gravity increment denoising fortune
The interval calculated.It is of course also possible to select the matrix of 3 × 3 as carrying out the interval of time shift gravity increment denoising computing,
This is not construed as limiting by the present invention.
Further, the number of repetition arranging time shift gravity increment denoising is 50 times.In the present embodiment, choosing
Select time shift gravity increment grid literary composition obtained after 50 numbers of repetition as denoising are because carrying out 50 denoisings
The denoising effect of part is optimal.Similarly, 60 denoisings when carrying out time shift gravity denoising, can also be carried out,
This is not construed as limiting by the present invention.
Correspondingly, when carrying out denoising, to selected to treat that denoising measuring point is interval, can carry out according to following rule
Weighting processes:
When A < 4 microgal and B < during 4 microgal, M=N；
When 4 microgal < A < 6 microgal and 4 microgal < B < during 6 microgal, M=N × 80%+P × 20%；
When 4 microgal<A<6 microgal and B>6 microgal or 4 microgal<B<6 microgal and A>6 microgal, the two
When condition meets alternative one, M=N × 60%+P × 40%；
As A > 6 microgal and during B > 6 microgal, M=N × 30%+P × 70%；
Wherein, A represents using default denoising interval as all measuring points in the limited range of border in first time period
The gravity variance of interior repetitive measurement, B represents the gravity variance of repetitive measurement within the second time period, and D represents to preset
Denoising is interval as the gravity variance of repetitive measurement in first time period of all measuring points in the limited range of border
And within the second time period the maximum in the gravity variance of repetitive measurement, M represents the current measuring point after denoising
Time shift gravity increment, N represents time shift gravity increment, and P represents interval as border limited range using default denoising
The meansigma methods of the time shift gravity increment of interior all measuring points.
Time shift gravity increment grid file, the first gravity variance grid file and the second gravity variance net to somewhere
During lattice file carries out denoising.It is possible, firstly, to determine the current measuring point carrying out denoising；Again with this survey
Centered by Dian, matrix using 5 × 5 interval as all measuring points in the limited range of border in first time period
The gravity variance of repetitive measurement, within the second time period the gravity variance of repetitive measurement and the time shift gravity increment of non-denoising
It is weighted processing, obtains the time shift gravity increment of the current measuring point after denoising；Then, working as after the denoising that will obtain
The time shift gravity increment of front measuring point, as the time shift gravity increment of non-denoising, repeats the process of above-mentioned denoising, until
After number of repetition reaches 50 times, the time shift gravity increment of the current measuring point after denoising can be obtained.Somewhere is surveyed
In fixed time shift gravity increment grid file, the first gravity variance grid file and the second gravity variance grid file
All measuring points are performed both by above-mentioned denoising process, can obtain the time shift gravity increment grid file after denoising.
Based on same inventive concept, the embodiment of the present invention additionally provides the denoising device of a kind of time shift gravity increment, as
Described in the following examples.Owing to the denoising device of time shift gravity increment solves principle and the time shift gravity increment of problem
Denoising method is similar, and therefore the enforcement of the denoising device of time shift gravity increment may refer to the denoising side of time shift gravity increment
The enforcement of method, repeats no more in place of repetition.Used below, term " unit " or " module " can realize
The software of predetermined function and/or the combination of hardware.Although the device described by following example preferably comes real with software
Existing, but hardware, or the realization of the combination of software and hardware also may and be contemplated.Fig. 2 is that the present invention is real
Execute a kind of structured flowchart of the denoising device of the time shift gravity increment of example, as in figure 2 it is shown, include: meansigma methods obtains mould
Block 201, variance acquisition module 202 and weighting processing module 203, illustrate this structure below.
Meansigma methods acquisition module 201, the time shift gravity increment of each measuring point in obtaining work area, wherein, time described
Move gravity increment to be the gravity meansigma methods of measuring point repetitive measurement in first time period and repeatedly survey within the second time period
Difference between the gravity meansigma methods of amount；
Variance acquisition module 202, for obtaining described each measuring point gravity variance of repetitive measurement in first time period
And within the second time period the gravity variance of repetitive measurement；
Weighting processing module 203, for each measuring point in each measuring point described is used following operation obtain described each
Time shift gravity increment after each measuring point denoising in measuring point: select current measuring point as central point, will be with described central point
Centered by, repeatedly survey in first time period using default denoising interval as all measuring points in the limited range of border
The gravity variance of amount, within the second time period, the gravity variance of repetitive measurement and time shift gravity increment are weighted processing,
Obtain the time shift gravity increment of the described current measuring point after denoising.
In one embodiment, described device can also include: replicated blocks, in obtaining each measuring point described
After time shift gravity increment after each measuring point denoising, by the time shift gravity after measuring point denoising each in each measuring point described
Increment, as the time shift gravity increment of non-denoising, repeats following denoising operation to each measuring point in each measuring point described
Until reaching default number of repetition: select current measuring point as central point, will be centered by described central point, with in advance
If denoising is interval as all measuring points gravity side of repetitive measurement in first time period in the limited range of border
Differ from, within the second time period, the gravity variance of repetitive measurement and the time shift gravity increment of non-denoising are weighted processing,
The time shift gravity increment of the described current measuring point after denoising.
In one embodiment, described default number of repetition is not less than 10 times.
In one embodiment, work as described in after described weighting processing module can be calculated denoising according to below equation
The time shift gravity increment of front measuring point:
When A < (D/2) and B < time (D/2), M=N；
When (D/2) < A < (3D/4) and (D/2) < B < time (3D/4), M=N × a+P × b；
When (D/2)<A<(3D/4) and B>(3D/4), (D/2)<when B<(3D/4) and A>(3D/4) meets alternative one,
M=N × c+P × d；
As A > (3D/4) and B > (3D/4) time, M=N × e+P × f；
Wherein, A represents using described default denoising interval as all measuring points in the limited range of border first
The gravity variance of repetitive measurement in time period, B represents the gravity variance of repetitive measurement within the second time period, and D represents
Interval as all measuring points repetitive measurement in first time period in the limited range of border using described default denoising
Gravity variance and within the second time period the maximum in the gravity variance of repetitive measurement, M represents the institute after denoising
Stating the time shift gravity increment of current measuring point, N represents described time shift gravity increment, and P represents interval with described default denoising
As the meansigma methods of the time shift gravity increment of all measuring points in the limited range of border, a, b, c, d, e, f
All represent the weight coefficient between 0% to 100%.
In one embodiment, described meansigma methods acquisition module is for carrying out at gridding moving gravity increment time described
Reason；Described variance acquisition module for the gravity variance of each measuring point described repetitive measurement in first time period and
Within the second time period, the gravity variance of repetitive measurement carries out gridding process.
In one embodiment, described gridding processing method includes: Kriging regression method.
In one embodiment, described default denoising interval includes: the matrix of n × n, and wherein, described n represents institute
State the measuring point number in default denoising interval.
As can be seen from the above description, the embodiment of the present invention achieves following technique effect: to time shift gravity
During increment carries out denoising, by all measuring points in first time period the gravity variance of repetitive measurement,
In two time periods, gravity variance and the time shift gravity increment of repetitive measurement are weighted processing, and obtain the current survey after denoising
The time shift gravity increment of point, the method i.e. using weighting to process efficiently solves and carries out in prior art not having during denoising
Have and consider to carry out time shift microgravity and repeatedly observe the data obtained and single is observed repeatedly between the data obtained by reading
The problem of contact, eliminate the observation error of time shift gravity increment, i.e. observation noise so that the data finally given
Higher with the actual geological condition goodness of fit.
Obviously, those skilled in the art should be understood that each module of the above-mentioned embodiment of the present invention or each step are permissible
Realizing with general calculating device, they can concentrate on single calculating device, or is distributed in multiple calculating
On the network that device is formed, alternatively, they can realize with calculating the executable program code of device, thus,
Can be stored in storing in device and be performed by calculating device, and in some cases, can be to be different from
The step shown or described by order execution herein, or they are fabricated to respectively each integrated circuit modules, or
Multiple modules in them or step are fabricated to single integrated circuit module and realize by person.So, the embodiment of the present invention
It is not restricted to any specific hardware and software combine.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for the skill of this area
For art personnel, the embodiment of the present invention can have various modifications and variations.All within the spirit and principles in the present invention,
Any modification, equivalent substitution and improvement etc. made, should be included within the scope of the present invention.
Claims (10)
1. the denoising method of a time shift gravity increment, it is characterised in that including:
The time shift gravity increment of each measuring point in acquisition work area, wherein, described time shift gravity increment is that measuring point is when first
Between in section repetitive measurement gravity meansigma methods and within the second time period difference between the gravity meansigma methods of repetitive measurement；
Obtain the gravity variance of described each measuring point repetitive measurement in first time period and many within the second time period
The gravity variance of secondary measurement；
In each measuring point described in the following operation of each measuring point employing in each measuring point described is obtained after each measuring point denoising
Time shift gravity increment: select current measuring point as central point, will be centered by described central point, with default denoising district
Between as all measuring points in the limited range of border in first time period the gravity variance of repetitive measurement, second
In time period, gravity variance and the time shift gravity increment of repetitive measurement are weighted processing, and obtain the current measuring point after denoising
Time shift gravity increment.
2. the method for claim 1, it is characterised in that each measuring point goes in obtaining each measuring point described
After time shift gravity increment after making an uproar, described method also includes:
Using the time shift gravity increment after measuring point denoising each in each measuring point described as the time shift gravity increment of non-denoising,
Each measuring point in each measuring point described repeats following denoising operate until reaching default number of repetition: select to work as
Front measuring point, will be interval as border limited range using default denoising centered by described central point as central point
In all measuring points in first time period the gravity variance of repetitive measurement, within the second time period the gravity of repetitive measurement
The time shift gravity increment of variance and non-denoising is weighted processing, and the time shift gravity obtaining the current measuring point after denoising increases
Amount.
3. method as claimed in claim 2, it is characterised in that described default number of repetition is not less than 10 times.
4. the method for claim 1, it is characterised in that be calculated the institute after denoising according to below equation
State the time shift gravity increment of current measuring point:
When A < (D/2) and B < time (D/2), M=N；
When (D/2) < A < (3D/4) and (D/2) < B < time (3D/4), M=N × a+P × b；
When (D/2)<A<(3D/4) and B>(3D/4), (D/2)<when B<(3D/4) and A>(3D/4) meets alternative one,
M=N × c+P × d；
As A > (3D/4) and B > (3D/4) time, M=N × e+P × f；
Wherein, A represents using described default denoising interval as all measuring points in the limited range of border first
The gravity variance of repetitive measurement in time period, B represents the gravity variance of repetitive measurement within the second time period, and D represents
Interval as all measuring points repetitive measurement in first time period in the limited range of border using described default denoising
Gravity variance and within the second time period the maximum in the gravity variance of repetitive measurement, M represents the institute after denoising
Stating the time shift gravity increment of current measuring point, N represents described time shift gravity increment, and P represents interval with described default denoising
As the meansigma methods of the time shift gravity increment of all measuring points in the limited range of border, a, b, c, d, e, f
All represent the weight coefficient between 0% to 100%.
5. the method for claim 1, it is characterised in that in acquisition work area, the time shift gravity of each measuring point increases
Amount, including: carry out gridding process to moving gravity increment time described；
Obtain the gravity variance of described each measuring point repetitive measurement in first time period and many within the second time period
The gravity variance of secondary measurement, including: to the gravity variance of each measuring point described repetitive measurement in first time period and
Within the second time period, the gravity variance of repetitive measurement carries out gridding process.
6. method as claimed in claim 5, it is characterised in that described gridding processing method includes: Ke Lijin
Interpolation method.
7. the method for claim 1, it is characterised in that described default denoising interval includes: the square of n × n
Battle array, wherein, described n represents the measuring point number in described default denoising interval.
8. the denoising device of a time shift gravity increment, it is characterised in that including:
Meansigma methods acquisition module, the time shift gravity increment of each measuring point in obtaining work area, wherein, described time shift weight
Power increment be measuring point repetitive measurement in first time period gravity meansigma methods and within the second time period repetitive measurement
Difference between gravity meansigma methods；
Variance acquisition module, for obtain described each measuring point repetitive measurement in first time period gravity variance and
The gravity variance of repetitive measurement within the second time period；
Weighting processing module, for using following operation to obtain each measuring point described to each measuring point in each measuring point described
In time shift gravity increment after each measuring point denoising: select current measuring point as central point, in will being with described central point
The heart, interval as all measuring points repetitive measurement in first time period in the limited range of border using default denoising
Gravity variance, within the second time period, the gravity variance of repetitive measurement and time shift gravity increment are weighted processing, and obtain
The time shift gravity increment of the described current measuring point after denoising.
9. device as claimed in claim 8, it is characterised in that described device also includes:
Replicated blocks, after the time shift gravity increment after each measuring point denoising in obtaining each measuring point described, will
In each measuring point described, time shift gravity increment after each measuring point denoising is as the time shift gravity increment of non-denoising, to described
In each measuring point, each measuring point repeats following denoising operation until reaching default number of repetition: select current measuring point
As central point, will be interval as the institute in the limited range of border using default denoising centered by described central point
Have measuring point in first time period the gravity variance of repetitive measurement, within the second time period repetitive measurement gravity variance and
The time shift gravity increment of non-denoising is weighted processing, and obtains the time shift gravity increment of the described current measuring point after denoising.
10. device as claimed in claim 9, it is characterised in that described default number of repetition is not less than 10 times.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610517205.6A CN105974488B (en) | 2016-07-04 | 2016-07-04 | A kind of denoising method and device of time shift gravity increment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610517205.6A CN105974488B (en) | 2016-07-04 | 2016-07-04 | A kind of denoising method and device of time shift gravity increment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105974488A true CN105974488A (en) | 2016-09-28 |
| CN105974488B CN105974488B (en) | 2018-02-02 |
Family
ID=56954726
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610517205.6A Active CN105974488B (en) | 2016-07-04 | 2016-07-04 | A kind of denoising method and device of time shift gravity increment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105974488B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115393547A (en) * | 2022-08-15 | 2022-11-25 | 中国地质大学（北京） | Omnidirectional filtering method and system for lunar satellite gravity anomaly data |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101769743A (en) * | 2010-01-04 | 2010-07-07 | 北京航空航天大学 | Distributed filtering device for MIMU and GPS combined navigation system |
| CN102236109A (en) * | 2010-05-06 | 2011-11-09 | 中国石油天然气集团公司 | Method for correcting interference of variance coefficient in gravity and magnetic interference zone |
| US20120082001A1 (en) * | 2010-10-01 | 2012-04-05 | Welker Kenneth E | Monitoring the Quality of Particle Motion Data During a Seismic Acquisition |
| US20140316708A1 (en) * | 2013-04-19 | 2014-10-23 | The Board Of Trustees Of The Leland Stanford Junior University | Oriented Wireless Structural Health and Seismic Monitoring |
| CN104570140A (en) * | 2013-10-29 | 2015-04-29 | 中国石油天然气集团公司 | High-precision time shifting microgravity acquiring method |
-
2016
- 2016-07-04 CN CN201610517205.6A patent/CN105974488B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101769743A (en) * | 2010-01-04 | 2010-07-07 | 北京航空航天大学 | Distributed filtering device for MIMU and GPS combined navigation system |
| CN102236109A (en) * | 2010-05-06 | 2011-11-09 | 中国石油天然气集团公司 | Method for correcting interference of variance coefficient in gravity and magnetic interference zone |
| US20120082001A1 (en) * | 2010-10-01 | 2012-04-05 | Welker Kenneth E | Monitoring the Quality of Particle Motion Data During a Seismic Acquisition |
| US20140316708A1 (en) * | 2013-04-19 | 2014-10-23 | The Board Of Trustees Of The Leland Stanford Junior University | Oriented Wireless Structural Health and Seismic Monitoring |
| CN104570140A (en) * | 2013-10-29 | 2015-04-29 | 中国石油天然气集团公司 | High-precision time shifting microgravity acquiring method |
Non-Patent Citations (2)
| Title |
|---|
| YUNXIANG LIU EL: ""Preliminary experiment of new time-lapse microgravity monitoring technique in T-gas reservoir"", 《SEG DENVER 2014 ANNUAL MEETING》 * |
| 刘云祥: ""重磁弱异常处理研究与应用"", 《勘探地球物理进展》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115393547A (en) * | 2022-08-15 | 2022-11-25 | 中国地质大学（北京） | Omnidirectional filtering method and system for lunar satellite gravity anomaly data |
| CN115393547B (en) * | 2022-08-15 | 2023-03-10 | 中国地质大学（北京） | Omnidirectional filtering method and system for lunar satellite gravity anomaly data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105974488B (en) | 2018-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wellmann et al. | Uncertainties have a meaning: Information entropy as a quality measure for 3-D geological models | |
| Shahzad et al. | TecDEM: A MATLAB based toolbox for tectonic geomorphology, Part 1: Drainage network preprocessing and stream profile analysis | |
| Deutsch | A sequential indicator simulation program for categorical variables with point and block data: BlockSIS | |
| EP3489725B1 (en) | Method for estimating geological architecture of a geologic volume | |
| Olalotiti-Lawal et al. | A multiobjective Markov chain Monte Carlo approach for history matching and uncertainty quantification | |
| Page et al. | The UCERF3 grand inversion: Solving for the long‐term rate of ruptures in a fault system | |
| Englund | Spatial simulation: environmental applications | |
| Wang et al. | Comparison of probabilistic seismic‐hazard models using averaging‐based factorization | |
| Julio et al. | Sampling the uncertainty associated with segmented normal fault interpretation using a stochastic downscaling method | |
| Gomez-Hernandez et al. | Stochastic characterization of gridblock permeabilities | |
| Console et al. | A physics-based earthquake simulator and its application to seismic hazard assessment in Calabria (Southern Italy) region | |
| Eidsvik et al. | Simulation–regression approximations for value of information analysis of geophysical data | |
| EP1606727A2 (en) | Method for quickly forming a stochastic method representing the distribution of a physical variable in a heterogeneous environment by appropriate selection of geostatic realizations | |
| Kiær et al. | Calendar time interpolation of amplitude maps from 4D seismic data | |
| AU2017219462B8 (en) | Method of calculating radiogenic heat production | |
| CN105974488A (en) | De-noising method and device for time shifting gravity increment | |
| Smith | Earthquake hazard and risk assessment in New Zealand by Monte Carlo methods | |
| AU2013303068A1 (en) | Methods and corresponding software module for quantifying risks or likelihoods of hydrocarbons being present in a geological basin or region | |
| Schueth et al. | The Probable Datum Method (PDM): a technique for estimating the age of origination or extinction of nannoplankton | |
| Artman et al. | Wave-equation microseismic imaging and event selection in the image domain | |
| Alpak et al. | Stochastic history matching of a deepwater turbidite reservoir | |
| Agi et al. | Impact of geological interpretation on reservoir 3D static model: Workflow, methodology approach and delivery process | |
| Macé et al. | Improved Resources and Reserves Estimation with the Combined Use of Scenario-Based Geological Models and Flow-Based Connectivity | |
| Echeverria et al. | A robust scheme for spatio-temporal inverse modeling of oil reservoirs | |
| MX2014015955A (en) | Truncation diagram determination for a pluri-gaussian estimation. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |