CN105891889B - A kind of method and device of gravity anomaly border enhancing - Google Patents
A kind of method and device of gravity anomaly border enhancing Download PDFInfo
- Publication number
- CN105891889B CN105891889B CN201610197300.2A CN201610197300A CN105891889B CN 105891889 B CN105891889 B CN 105891889B CN 201610197300 A CN201610197300 A CN 201610197300A CN 105891889 B CN105891889 B CN 105891889B
- Authority
- CN
- China
- Prior art keywords
- gravity anomaly
- mesh point
- value
- data
- enhancement
- 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
- 230000005484 gravity Effects 0.000 title claims abstract description 176
- 238000000034 method Methods 0.000 title claims abstract description 79
- 230000002708 enhancing effect Effects 0.000 title claims abstract description 38
- 238000004364 calculation method Methods 0.000 claims description 44
- 239000000284 extract Substances 0.000 claims description 15
- 238000000605 extraction Methods 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 abstract description 9
- 230000003321 amplification Effects 0.000 abstract description 6
- 238000003199 nucleic acid amplification method Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 15
- 230000002159 abnormal effect Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 230000002547 anomalous effect Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 238000011017 operating method Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000003708 edge detection Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000001965 increasing effect Effects 0.000 description 2
- 230000010365 information processing Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000004513 sizing Methods 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000009828 non-uniform distribution Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 229920001296 polysiloxane Polymers 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
- G06V10/443—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components by matching or filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Geophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Image Processing (AREA)
Abstract
A kind of method and device of gravity anomaly border enhancing provided by the embodiments of the present application.The described method includes：The gravity anomaly data in target work area are obtained, and gridding is carried out to the gravity anomaly data；Gravity anomaly boundary extremum information of first mesh point in preset direction is extracted using spatial domain operator one by one；Gravity anomaly boundary extremum information based on first mesh point, calculates the relevant enhancement value of the second mesh point.The gravity anomaly border Enhancement Method that the application provides carries out spatial domain processing to the gravity anomaly data, avoids amplification of the legacy frequencies domain algorithm to data noise, can obtain clear, continuous gravity anomaly boundary information.
Description
Technical field
The present invention relates to the technical field of geophysical exploration in oil and Solid Mineral exploration, more particularly to a kind of gravity
The method and device of exception boundary enhancing.
Background technology
In technical field of geophysical exploration, by gravimetric prospecting it will be seen that the distribution of geology field source, including field
The form and scale in source etc..Gravity anomaly is the combined reaction of ground lower density body uneven distribution, using in gravity anomaly data
Abnormal sideline property step band, it is abnormal distort, the information data such as abnormal partition characteristics can identify the abnormal geological border, than
Such as stratum breaking position.And the information related with geological boundry is usually expressed as high-frequency, low amplitude value, often be difficult directly from
The horizontal level of geological boundry is extracted in Bouguer map.It is generally necessary to first to the boundary information of gravity abnormal data into
The abnormal geological boundary information is extracted again after row enhancing.
The method of exception boundary enhancing of the prior art mainly has derivative class algorithm, and the derivative class algorithm mainly includes
Vertical Derivative method, aggregate level derivative method, total gradient modulus method, oblique derivative method, Xi Ta (Theta) figure method etc..Wherein, the aggregate level
Derivative method can use the method for carrying out evolution to the quadratic sum of transverse axis and plotted horizontal derivative to try to achieve.Derivative class algorithm
Common feature is amplified the radio-frequency component in signal in the combination using Vertical Derivative or derivative, makes in gravity anomaly map
Maximum or flex point be located approximately at the surface of geological objects boundary, subsequently exception can be extracted by the method for edge detection
Geological boundry position.
Derivative class algorithm is generally carried out using the method for frequency domain Fourier transform, and therefore, derivative class algorithm is often to original
The noise of data has amplification, and the directionality of derivative extreme value and continuity are poor, to the Boundary Recognition ability of weak anomaly
It is weaker, cause the extraction accuracy to the abnormal geological border relatively low.
The content of the invention
The purpose of the application is to provide a kind of method and device of gravity anomaly border enhancing, can reduce gravity anomaly
The noise of boundary information, improves the resolving power on gravity anomaly border.
The present invention provides a kind of method and device of gravity anomaly border enhancing, the method and device are specifically so
Realize：
A kind of method of gravity anomaly border enhancing, the described method includes：
The gravity anomaly data in target work area are obtained, and gridding is carried out to the gravity anomaly data；
Gravity anomaly boundary extremum information of first mesh point in preset direction is extracted using spatial domain operator one by one；
Gravity anomaly boundary extremum information based on first mesh point, calculates the relevant enhancement value of the second mesh point.
Optionally, in one embodiment of the invention, the gravity anomaly border pole based on first mesh point
Value information, calculates the relevant enhancement value of the second mesh point, including：
Gravity anomaly boundary extremum information based on first mesh point, is established centered on the second mesh point, size
For the enhancement value calculation window of the first preset value；
Adjacent data that extract eight directions adjacent with the enhancement value calculation window, that size is first preset value
Window；
The related coefficient of the enhancement value calculation window and the adjacent data window in adjacent eight directions is calculated respectively
Value；
Relevant enhancement value using eight evolution values of the product of eight correlation coefficient value as second mesh point.
Optionally, in one embodiment of the invention, after the progress gridding to the gravity data, also wrap
Include：
Noise in the gravity anomaly data is removed using low pass filtering method.
Optionally, in one embodiment of the invention, the method further includes：
Judge whether the relevant enhancement value of second mesh point is less than predetermined threshold value；
When judging result for when being, the relevant enhancing value information based on first mesh point, calculates second grid
The secondary relevant enhancement value of point.
Optionally, in one embodiment of the invention, the use spatial domain operator extracts the first mesh point and exists one by one
The gravity anomaly boundary extremum information of preset direction includes：
Spatial domain operator of the target work area on preset direction is obtained, the size of the spatial domain operator is pre- for second
If value；
Extract in the first mesh point centered on mesh point to be calculated, the extreme value calculation window that size is the second preset value；
Product using the extreme value calculation window and the spatial domain operator is different as the gravity of the mesh point to be calculated
Normal boundary extremum；
Other mesh points in first mesh point are traveled through, it is different in the gravity of preset direction to generate first mesh point
Normal boundary extremum information.
A kind of device of gravity anomaly border enhancing, described device include：
Gridding unit, net is carried out for obtaining the gravity anomaly data in target work area, and to the gravity anomaly data
Format；
Boundary extremum computing unit, for extracting gravity of first mesh point in preset direction one by one using spatial domain operator
Exception boundary extreme value information；
Enhancement value computing unit, for the gravity anomaly boundary extremum information based on first mesh point, calculates second
The relevant enhancement value of mesh point.
Optionally, in one embodiment of the invention, the enhancement value computing unit includes：
Enhancement value window establishes unit, for the gravity anomaly boundary extremum information based on first mesh point, establishes
Centered on the second mesh point, the enhancement value calculation window that size is the first preset value；
Adjacent window apertures extraction unit, for extract adjacent with the enhancement value calculation window eight directions, size be institute
State the adjacent data window of the first preset value；
Related coefficient computing unit, for calculating the enhancement value calculation window and the phase in adjacent eight directions respectively
The correlation coefficient value of adjacent data window；
Square root unit, for using eight evolution values of the product of eight correlation coefficient value as second net
The relevant enhancement value of lattice point.
Optionally, in one embodiment of the invention, described device further includes：
Denoising unit, for removing the noise in the gravity anomaly data using low pass filtering method.
Optionally, in one embodiment of the invention, described device further includes：
Whether enhancement value judging unit, the relevant enhancement value for judging second mesh point are less than predetermined threshold value；
Secondary enhancement unit, for when the enhancement value judging unit judging result for be when, based on first net
The relevant enhancing value information of lattice point, calculates the secondary relevant enhancement value of second mesh point.
Optionally, in one embodiment of the invention, the boundary extremum computing unit includes：
Operator acquiring unit, for obtaining spatial domain operator of the target work area on preset direction, the spatial domain
The size of operator is the second preset value；
Extreme value window extraction unit, for extracting in the first mesh point centered on mesh point to be calculated, size second
The extreme value calculation window of preset value；
Extreme value computing unit, based on the product of the extreme value calculation window and the spatial domain operator is treated as described in
Calculate the gravity anomaly boundary extremum of mesh point；
Computing unit is traveled through, for traveling through other mesh points in first mesh point, generates first mesh point
In the gravity anomaly boundary extremum information of preset direction.
The gravity anomaly border Enhancement Method and device that the application provides, can extract the gravity anomaly data after gridding
Gravity anomaly boundary extremum information on preset direction, the extreme value on the gravity anomaly border generally correspond to the border of model
Point.The gravity anomaly boundary extremum information is then based on, the relevant enhancement value of mesh point can be calculated, further obtain model
The enhancement value of boundary point.The gravity anomaly border Enhancement Method that the application provides carries out at spatial domain the gravity anomaly data
Reason, avoids amplification of the legacy frequencies domain algorithm to data noise, can obtain clear, continuous gravity anomaly border letter
Breath, the research to target work area stratum breaking data provide strong basic data.
Brief description of the drawings
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, below will be to embodiment or existing
There is attached drawing needed in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, in the premise of not making the creative labor property
Under, other attached drawings can also be obtained according to these attached drawings.
Fig. 1 is a kind of method flow schematic diagram of embodiment of gravity anomaly border provided by the invention Enhancement Method；
Fig. 2 is a kind of method flow diagram of embodiment of calculated gravity anomaly boundary extremum provided by the invention；
Fig. 3 be it is provided by the invention calculating mesh point be concerned with enhancement value a kind of embodiment method flow diagram；
Fig. 4 is the position relationship of enhancement value calculation window provided by the invention and the adjacent data window in adjacent eight directions
Figure；
Fig. 5 is the actual measurement Bouguer map of China northeast target work area A provided by the invention；
Fig. 6 is the Bouguer map after the progress upward continuation denoising provided by the invention to work area A；
Fig. 7 is after the gravity anomaly boundary information provided by the invention by work area A east northeasts, northwest (NW) both direction is merged
Fracture trend graph；
Fig. 8 is the gravity anomaly boundary information provided by the invention for being extracted and being strengthened to work area A using aggregate level derivative method
Figure；
Fig. 9 is a kind of modular structure schematic diagram of embodiment of gravity anomaly border provided by the invention intensifier；
Figure 10 is a kind of modular structure schematic diagram of embodiment of enhancement value computing unit provided by the invention；
Figure 11 is the modular structure schematic diagram of another embodiment of gravity anomaly border provided by the invention intensifier；
Figure 12 is a kind of modular structure schematic diagram of embodiment of boundary extremum computing unit provided by the invention.
Embodiment
It is in order to make those skilled in the art better understand the technical solutions in the application, real below in conjunction with the application
The attached drawing in example is applied, the technical solution in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example is merely a part but not all of the embodiments of the present application.It is common based on the embodiment in the application, this area
Technical staff's all other embodiments obtained without making creative work, should all belong to the application protection
Scope.
The method of gravity anomaly border described herein enhancing is described in detail below in conjunction with the accompanying drawings.Fig. 1 is
A kind of method flow schematic diagram of the embodiment for the gravity anomaly border Enhancement Method that the application provides.Although this application provides
Such as following embodiments or method operating procedure shown in the drawings, but based on routine or without performing creative labour in the method
In can include more or less operating procedures.In the step of necessary causality is not present in logicality, these steps
Rapid execution sequence is not limited to execution sequence provided by the embodiments of the present application.The gravity anomaly border of the method in practice
It can perform according to embodiment or method shown in the drawings order or perform parallel when either device performs during enhancing
(such as environment of parallel processor or multiple threads).
A kind of embodiment for the gravity anomaly border Enhancement Method that specific the application provides is as shown in Figure 1, the method
It can include：
S1：The gravity anomaly data in target work area are obtained, and gridding is carried out to the gravity anomaly data.
Gravimetric observation data are usually the gravity in the target work area that geological exploration personnel are observed by gravity measuring equipment
Observe data.When the gravimetric observation data of observation station and the normal gravity of the observation station are usually present deviation, it may be determined that
Gravity anomaly occurs for the observation station.Cause the gravity anomaly of observation station to generally there are many factors, for example, the observation station with
Geoid is there are certain height difference, and when observation station position is higher than the earth horizontal plane, gravity value is smaller less than normal；When observation station with
There are additional gravity value during surplus materials, can be also produced at observation station between geoid.
Gravity anomaly data in the embodiment of the present application can be the bouguer gravity anomaly data in the target work area.This reality
Apply in example, Bouguer correction can be carried out to the influence factor for causing gravity anomaly in the gravimetric observation data, if in Bouguer school
The exception being still had after just is known as bouguer gravity anomaly.Can by the gravimetric observation data after correction with corresponding observation station
Normal gravity difference as bouguer gravity anomaly, difference claims normal anomaly to be positive, and difference claims negative anomaly to be negative.To described heavy
After power observation data carry out Bouguer processing, some influences for causing the objective factor of gravity anomaly to be brought are eliminated, can be significantly
The huge data volume of gravimetric observation data is reduced, significantly improves the follow-up efficiency for calculating processing.
After the gravity anomaly data in target work area are obtained, gridding can be carried out to the gravity anomaly data.Tool
Body, Bouguer map can be divided into square regular grid, the square rule with certain sizing grid
The size and shape of grid can be depending on research purpose and working scale, and sizing grid is generally gravimetric network compartment
Away from several times to more than ten times.Using the average value of the observer nodes GRAVITY ANOMALIES included in square regular grid as it is described just
The typical value of square regular grid.Gridding is carried out to the gravity anomaly data, can be by the weight of spatially non-uniform distribution
Power abnormal data homogenizes, and on the one hand can suppress local noise, on the other hand can make up the numerical value of blank grid.
In the embodiment of the present application, the gravity anomaly data in target work area can be obtained, and to the gravity anomaly data into
The gravity anomaly data are carried out gridding, can suppress local noise to a certain extent by row gridding.
In the embodiment of the present application, after the progress gridding to the gravity data, it can also include：
Noise in the gravity anomaly data is removed using low pass filtering method.
The low pass filtering method can include moving average method, upward continuation method etc..Wherein, for the upward of gravity anomaly
Continuation can protrude the off-note of larger anomalous body, the larger anomalous body can include it is zonal or
The fairly large anomalous body in person deep, and the off-note of the anomalous body of scale is smaller is suppressed, the anomalous body of the scale is smaller
It can include local or shallow and small anomalous body.The upward continuation method can eliminate near-surface effect, and it is different to remove gravity
The influence of noise of regular data.
S2：Gravity anomaly boundary extremum information of first mesh point in preset direction is extracted using spatial domain operator one by one.
Usually, gravity anomaly has directionality, and the gravity anomaly direction of the same area is often consistent, different regions
Gravity anomaly direction generally there is larger difference.Such as south China in China, Lower Yangtze massif, the weight of northeast Song Liao plains regions
Power extremely in east northeast to distribution, based on the Inner Mongol, Gansu, the gravity anomaly in Xinjiang and loess plateau are arranged with East and West direction.This reality
Apply in example, can determine to carry out the preset direction of gravity anomaly border enhancing according to the provincial characteristics of objective area.Fig. 2 is this
Apply for a kind of method flow diagram of embodiment of the calculated gravity anomaly boundary extremum provided, as shown in Fig. 2, the method bag
Include：
S11：Obtain spatial domain operator of the target work area on preset direction, the size of the spatial domain operator is the
Two preset values.
Spatial domain operator in the present embodiment can include direction difference spatial domain operator, wherein, as shown in table 1-2, east
West to difference operator and the difference operator of North and South direction be Sobel (Sobel) operator, Sobel Operator can be to grid
Gravity anomaly data after change are weighted average filter, have certain noise inhibiting ability.As shown in table 3-4, the north
East is difference operator to difference operator and northwest (NW) direction difference operator, can equally be had to the gravity anomaly data after gridding
Noise suppressed acts on.Certainly, the spatial domain operator in the present embodiment is not limited only to direction difference spatial domain operator, can also wrap
Include Roberts operators, Canny operators, Prewitt operators, edge detection operator, Second-order Operator etc..Usually, spatial domain operator
Size be fixed, such as the size of Sobel Operator is 3 × 3 matrixes, can obtain the size of spatial domain operator as
Two preset values, for follow-up selection data window size.
It should be noted that according to the provincial characteristics in target work area, the general direction difference operator for only choosing preset direction
Extreme value processing is carried out to the gravity anomaly data after the gridding.For example, positioned at the target work area in China south China, can utilize
East northeast direction difference operator extracts the gravity anomaly boundary extremum information in the target work area.
S12：Extract in the first mesh point centered on mesh point to be calculated, the extreme value calculating window that size is the second preset value
Mouthful.
The first mesh point in the present embodiment can be the set of mesh point in the gravity anomaly data after gridding, described
Mesh point to be calculated belongs to the single mesh point in the first mesh point.Extraction is centered on the mesh point to be calculated, size is
The extreme value calculation window of second preset value.For example, extraction is with mesh point a to be calculatedi,jCentered on size be 3 × 3 extreme value meter
Window is calculated, the extreme value calculation window data are：
S13：Weight using the product of the extreme value calculation window and the spatial domain operator as the mesh point to be calculated
Power exception boundary extreme value.
In the embodiment of the present application, the data in the extreme value calculation window and the product of the spatial domain operator can be made
For the gravity anomaly boundary extremum of the mesh point to be calculated.For example, based on the direction difference spatial domain operator of above-mentioned table 1-4 with
And with mesh point a to be calculatedi,jCentered on extreme value calculation window, the mesh point a to be calculated can be calculated respectivelyi,jIn east
West, north and south, east northeast, the gravity anomaly boundary extremum in northwest (NW) direction.Specifically, for east-west direction, gravity anomaly boundary extremum
Ai,jFormula be：
Ai,j=ai-1,j+1-ai-1,j-1+2(ai,j+1-ai,j-1)+ai+1,j+1-ai+1,j-1；
For North and South direction, gravity anomaly boundary extremum Ai,jFormula be：
Ai,j=ai+1,j+1-ai-1,j+1+2(ai+1,j-ai-1,j)+ai+1,j-1-ai-1,j-1；
For east northeast direction, gravity anomaly boundary extremum Ai,jFormula be：
Ai,j=ai+1,j-ai,j+1+ai,j-1-ai-1,j；
For northwest (NW) direction, gravity anomaly boundary extremum Ai,jFormula be：
Ai,j=ai,j+1-ai-1,j+ai+1,j-ai,j-1。
S14：Other mesh points in first mesh point are traveled through, generate weight of first mesh point in preset direction
Power exception boundary extreme value information.
After the gravity anomaly boundary extremum of the mesh point to be calculated is calculated, it can calculate after the same method
Other mesh points in first mesh point, until the gravity anomaly border of all mesh points in the first mesh point is calculated
Extreme value, generates gravity anomaly boundary extremum information of first mesh point in preset direction.
It should be noted that in gravity anomaly data after the gridding, the mesh point in marginal position is general
Extreme value calculation window can not be generated, at this time, it may be necessary to carry out flared end to the gravity anomaly data., can be with the embodiment of the present application
The method that zero is decayed to using cosine carries out flared end to the gravity anomaly data after the gridding.If the gravity after gridding is different
Regular data is ZY (1:m,1:N), m is line number, and n is the points on every line, if flared end number is t, for example, as t=2, flared end
Data afterwards are ZK (1-t：M+t, 1-t：n+t).ZK (i, j)=ZY (i, j), (1≤i≤m, 1≤j≤n), then in ZK are set
Data can flared end as follows：
ZK (i, j)=ZY (i, 1) × cos ((π × (1-j))/2t), (1≤i≤m, 1-t≤j≤0)；
ZK (i, j)=ZY (i, n) × cos ((π × (j-n))/2t), (1≤i≤m, n+1≤j≤n+t)；
ZK (i, j)=ZY (1, j) × cos ((π × (1-i))/2t), (1-t≤i≤0,1-t≤j≤n+t)；
ZK (i, j)=ZY (m, j) × cos ((π × (i-m))/2t), (m+1≤i≤m+t, 1-t≤j≤n+t).
In the present embodiment, the gravity anomaly data after gridding can be converted into extreme value information, be concerned with for follow-up
Enhancing handle, can to avoid in the algorithm of legacy frequencies domain to the amplification of the gravity anomaly data noise after gridding.
S3：Gravity anomaly boundary extremum information based on first mesh point, calculates the relevant enhancing of the second mesh point
Value.
The second mesh point in the present embodiment is any one mesh point in first mesh point, based on the target
The gravity anomaly boundary extremum information of first mesh point in work area, can calculate the relevant enhancement value of the second mesh point.Fig. 3 is this
Application provide calculating mesh point be concerned with enhancement value a kind of embodiment method flow diagram, as shown in figure 3, the method bag
Include：
S101：Gravity anomaly boundary extremum information based on first mesh point, establish centered on the second mesh point,
Size is the enhancement value calculation window of the first preset value.
In the present embodiment, it can be established based on the gravity anomaly boundary extremum information of first mesh point with the second net
Centered on lattice point, size be the first preset value enhancement value calculation window.The first preset value size can be according to integral net
Lattice number is set, and in the present embodiment, first preset value can be 5 × 5.
S102：Phase that extract eight directions adjacent with the enhancement value calculation window, that size is first preset value
Adjacent data window.
In the present embodiment, the related coefficient of the enhancement value calculation window and the window in adjacent eight directions can be calculated
Value.Therefore, eight directions adjacent with the enhancement value calculation window, equal-sized adjacent data window can be extracted.Fig. 4
It is the location diagram of enhancement value calculation window and the adjacent data window in adjacent eight directions that the application provides, such as Fig. 4 institutes
Show, wherein, northwest adjacent data window is upwards then to one grid of left by the position of the enhancement value calculation window
Obtained the window's position.
S103：It is related to the adjacent data window in adjacent eight directions that the enhancement value calculation window is calculated respectively
Coefficient value.
In the present embodiment, the adjacent data of the enhancement value calculation window and adjacent eight directions can be calculated respectively
The correlation coefficient value of window.Wherein, the calculation formula of the correlation coefficient value can be：
Wherein, X can be enhancement value calculation window data, and Y can be the data in either one adjacent data windowThe average of window total data is respectively corresponded to, n is the first preset value, and k is 1~8.
For example, based on first mesh point in the gravity anomaly boundary extremum information data of preset direction, X and Y can be with
For data shown in following formula, according to the calculation formula of above-mentioned correlation coefficient value, the second mesh point X can be calculated13With eight phases
The correlation coefficient value of adjacent data window.
S104：Relevant increasing using eight evolution values of the product of eight correlation coefficient value as second mesh point
Intensity values.
Assuming that 8 correlation coefficient value being calculated are r1-8, the calculation formula of the relevant enhancement value of second mesh point
Can be：
Follow-up, can also relevant enhanced gravity anomaly data be carried out with secondary relevant enhancing, it is described secondary relevant
Enhancing can comprise the following steps：
SS1：Judge whether the relevant enhancement value of second mesh point is less than predetermined threshold value；
SS2：When judging result for when being, the relevant enhancing value information based on first mesh point, calculates described second
The secondary relevant enhancement value of mesh point.
When the relevant enhancement value for judging second mesh point is less than predetermined threshold value, first mesh point can be based on
Relevant enhancing value information, calculate the secondary relevant enhancing of second mesh point according to the method and step of above-mentioned S101-S104
Value.
Illustrate above-described embodiment method below by a specific application scenarios, Fig. 5 is China east provided by the invention
The actual measurement Bouguer map of northern target work area A, as shown in figure 5, the depth of color represents the height of GRAVITY ANOMALIES in figure,
Color is deeper, represents that GRAVITY ANOMALIES is higher, conversely, color is more shallow, represents that GRAVITY ANOMALIES is lower.Utilize upward continuation 2,000
The method of rice removes the noise of Bouguer map in Fig. 5, obtain shown in Fig. 6 upward continuation denoising is carried out to work area A after
Bouguer map.According to geological condition known to the A of work area, using east northeast, northwest (NW) direction difference operator extraction east northeast, north
The gravity anomaly boundary extremum information of western both direction.East northeast based on extraction, northwest (NW) both direction gravity anomaly boundary extremum
Information, carries out each mesh point in Fig. 6 relevant enhancing processing, and the first time for obtaining work area A is concerned with enhancing data.Then adopt
The enhancing data that are concerned with threshold value 0.985 to the first time carry out secondary relevant enhancing processing, obtain area's east northeast to and northwest (NW)
To gravity anomaly boundary information, Fig. 7 be will be after work area A east northeasts, the fusion of northwest (NW) both direction gravity anomaly boundary information it is disconnected
Trend graph is split, reflects the main body distribution situation of this area's rift fracture.Fig. 8 is utilization aggregate level derivative method provided by the invention
The gravity anomaly boundary information figure for extracting and strengthening to work area A, by Fig. 7 compared with Fig. 8, it can be found that the present invention is extracted and increased
Strong gravity anomaly boundary information becomes apparent from, continuously.In addition, embodiment method provided by the invention not only may be used as shown in Figure 7
With extract and strengthen work area A east northeasts to main faults, but also distortion to orthogonal direction (NW trending) and bad break fracture are anti-
What is reflected is clearer, is demonstrated by having stronger disposal ability than traditional class algorithm reciprocal.
The gravity anomaly border Enhancement Method that the application provides, can extract the gravity anomaly data after gridding default
Gravity anomaly boundary extremum information on direction, the extreme value on the gravity anomaly border generally correspond to the boundary point of model.So
The gravity anomaly boundary extremum information is based on afterwards, can be calculated the relevant enhancement value of mesh point, further be obtained model boundary
The enhancement value of point.The gravity anomaly border Enhancement Method that the application provides carries out spatial domain processing to the gravity anomaly data,
Amplification of the legacy frequencies domain algorithm to data noise is avoided, clear, continuous gravity anomaly boundary information can be obtained, it is right
The research of target work area stratum breaking data provides strong basic data.
Another aspect of the present invention also provides a kind of device of gravity anomaly border enhancing, and Fig. 9 is the gravity that the application provides
The modular structure schematic diagram of a kind of embodiment of exception boundary intensifier, as shown in figure 9, described device 90 can include：
Gridding unit 91, carries out for obtaining the gravity anomaly data in target work area, and to the gravity anomaly data
Gridding；
Boundary extremum computing unit 92, for extracting weight of first mesh point in preset direction one by one using spatial domain operator
Power exception boundary extreme value information；
Enhancement value computing unit 93, for the gravity anomaly boundary extremum information based on first mesh point, calculates the
The relevant enhancement value of two mesh points.
The gravity anomaly border intensifier that the application provides, can extract the gravity anomaly data after gridding default
Gravity anomaly boundary extremum information on direction, the extreme value on the gravity anomaly border generally correspond to the boundary point of model.So
The gravity anomaly boundary extremum information is based on afterwards, can be calculated the relevant enhancement value of mesh point, further be obtained model boundary
The enhancement value of point.The gravity anomaly border Enhancement Method that the application provides carries out spatial domain processing to the gravity anomaly data,
Amplification of the legacy frequencies domain algorithm to data noise is avoided, clear, continuous gravity anomaly boundary information can be obtained, it is right
The research of target work area stratum breaking data provides strong basic data.
Figure 10 is a kind of modular structure schematic diagram of the embodiment for the enhancement value computing unit that the application provides, such as Figure 10 institutes
Show, the enhancement value computing unit 93 can include：
Enhancement value window establishes unit 101, for the gravity anomaly boundary extremum information based on first mesh point, builds
Stand centered on the second mesh point, the enhancement value calculation window that size is the first preset value；
Adjacent window apertures extraction unit 102, for extract adjacent with the enhancement value calculation window eight directions, size be
The adjacent data window of first preset value；
Related coefficient computing unit 103, for calculating the enhancement value calculation window and adjacent eight directions respectively
Adjacent data window correlation coefficient value；
Square root unit 104, for using eight evolution values of the product of eight correlation coefficient value as described
The relevant enhancement value of two mesh points.
Optionally, in one embodiment of the invention, described device further includes：
Denoising unit, for removing the noise in the gravity anomaly data using low pass filtering method.
Figure 11 is the modular structure schematic diagram of another embodiment for the gravity anomaly border intensifier that the application provides,
As shown in figure 11, described device 110 can also include：
Whether enhancement value judging unit 111, the relevant enhancement value for judging second mesh point are less than predetermined threshold value；
Secondary enhancement unit 112, for when the enhancement value judging unit judging result for be when, based on described first
The relevant enhancing value information of mesh point, calculates the secondary relevant enhancement value of second mesh point.
Figure 12 be the application provide boundary extremum computing unit a kind of embodiment modular structure schematic diagram, such as Figure 12
Shown, the boundary extremum computing unit 92 can include：
Operator acquiring unit 121, for obtaining spatial domain operator of the target work area on preset direction, the space
The size of domain operator is the second preset value；
Extreme value window extraction unit 122, for extracting in the first mesh point centered on mesh point to be calculated, size the
The extreme value calculation window of two preset values；
Extreme value computing unit 123, for using the product of the extreme value calculation window and the spatial domain operator as described in
The gravity anomaly boundary extremum of mesh point to be calculated；
Computing unit 124 is traveled through, for traveling through other mesh points in first mesh point, generates first grid
Gravity anomaly boundary extremum information of the point in preset direction.
Each embodiment in this specification is described by the way of progressive, what each embodiment stressed be with
The difference of other embodiment, between each embodiment identical similar part mutually referring to.It is real especially for system
For applying example, since it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changes without departing from the application's
Spirit.
Although being mentioned in teachings herein, picture is cut, RGB channel color is changed, finely tunes (fine-tuning), convolution god
Pictorial information processing, the description of Artificial Neural Network Structures through network model Rotating fields, convolution, loss feedback etc., still, this
Application be not limited to must be complete standard or mentioned mode information processing, the feelings of Artificial Neural Network Structures
Condition.Foregoing description in the application involved by each embodiment is only the application in some embodiments in the application, some
Embodiment amended slightly can also carry out the scheme of each embodiment of above-mentioned the application on the basis of standard, model, method.
Certainly, other deformations without creativeness of the process method step described in the application the various embodiments described above are being met, still
It can realize identical application, details are not described herein.
Although this application provides the method operating procedure as described in embodiment or flow chart, based on conventional or noninvasive
The work for the property made can include more or less operating procedures.The step of being enumerated in embodiment order is only numerous steps
A kind of mode in execution sequence, does not represent unique execution sequence., can when device or client production in practice performs
With according to embodiment, either method shown in the drawings order is performed or parallel performed (such as at parallel processor or multithreading
The environment of reason).
The device or module that above-described embodiment illustrates, can specifically be realized by computer chip or entity, or by with
The product of certain function is realized.For convenience of description, it is divided into various modules during description apparatus above with function to describe respectively.
The function of each module can be realized in same or multiple softwares and/or hardware when implementing the application.It is of course also possible to
The module for realizing certain function is combined into realization by multiple submodule or subelement.
Method, apparatus described herein or module can be realized in a manner of computer readable program code controller by
Any appropriate mode is realized, for example, controller can take such as microprocessor or processor and storage can be by (micro-)
The computer-readable medium of computer readable program code (such as software or firmware) that processor performs, logic gate, switch, specially
With integrated circuit (Application Specific Integrated Circuit, ASIC), programmable logic controller (PLC) and embedding
Enter the form of microcontroller, the example of controller includes but not limited to following microcontroller：ARC 625D、Atmel AT91SAM、
Microchip PIC18F26K20 and Silicone Labs C8051F320, Memory Controller are also implemented as depositing
A part for the control logic of reservoir.It is also known in the art that except real in a manner of pure computer readable program code
, completely can be by the way that method and step be carried out programming in logic come controller with logic gate, switch, special beyond existing controller
The form of integrated circuit, programmable logic controller (PLC) and embedded microcontroller etc. realizes identical function.Therefore this controller
A kind of hardware component is considered, and hardware can also be considered as to the device for being used for realization various functions that its inside includes
Structure in component.Or even, the device for being used for realization various functions can be considered as the software either implementation method
Module can be the structure in hardware component again.
Part of module in herein described device can be in the general of computer executable instructions
Described in context, such as program module.Usually, program module includes performing particular task or realizes specific abstract data class
The routine of type, program, object, component, data structure, class etc..The application can also be put into practice in a distributed computing environment,
In these distributed computing environment, by performing task by communication network and connected remote processing devices.In distribution
In computing environment, program module can be located in the local and remote computer-readable storage medium including storage device.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required hardware.Based on such understanding, the technical solution of the application is substantially in other words
The part to contribute to the prior art can be embodied in the form of software product, can also pass through the implementation of Data Migration
During embody.The computer software product can be stored in storage medium, such as ROM/RAM, magnetic disc, CD, bag
Some instructions are included to use so that a computer equipment (can be that personal computer, mobile terminal, server, or network are set
It is standby etc.) perform method described in some parts of each embodiment of the application or embodiment.
Each embodiment in this specification is described by the way of progressive, the same or similar portion between each embodiment
Divide mutually referring to what each embodiment stressed is the difference with other embodiment.The whole of the application or
Person part can be used in numerous general or special purpose computing system environments or configuration.Such as：Personal computer, server calculate
Machine, handheld device or portable set, laptop device, mobile communication terminal, multicomputer system, based on microprocessor be
System, programmable electronic equipment, network PC, minicom, mainframe computer include the distribution of any of the above system or equipment
Formula computing environment etc..
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changes without departing from the application's
Spirit.
Claims (8)
- A kind of 1. method of gravity anomaly border enhancing, it is characterised in that the described method includes：The gravity anomaly data in target work area are obtained, and gridding is carried out to the gravity anomaly data；Gravity anomaly boundary extremum information of first mesh point in preset direction is extracted using spatial domain operator one by one；Gravity anomaly boundary extremum information based on first mesh point, is established centered on the second mesh point, size The enhancement value calculation window of one preset value；Adjacent data window that extract eight directions adjacent with the enhancement value calculation window, that size is first preset value Mouthful；The correlation coefficient value of the enhancement value calculation window and the adjacent data window in adjacent eight directions is calculated respectively；Relevant enhancement value using eight evolution values of the product of eight correlation coefficient value as second mesh point.
- 2. the method for a kind of gravity anomaly border enhancing according to claim 1, it is characterised in that described to described heavy After power data carries out gridding, further include：Noise in the gravity anomaly data is removed using low pass filtering method.
- 3. the method for a kind of gravity anomaly border enhancing according to claim 1, it is characterised in that the method is also wrapped Include：Judge whether the relevant enhancement value of second mesh point is less than predetermined threshold value；When judging result for when being, the relevant enhancing value information based on first mesh point, calculates second mesh point Secondary relevant enhancement value.
- 4. the method for a kind of gravity anomaly border enhancing according to claim 1, it is characterised in that described to use spatial domain Operator extracts the first mesh point one by one to be included in the gravity anomaly boundary extremum information of preset direction：Spatial domain operator of the target work area on preset direction is obtained, the size of the spatial domain operator is default for second Value；Extract in the first mesh point centered on mesh point to be calculated, the extreme value calculation window that size is the second preset value；Gravity anomaly side using the product of the extreme value calculation window and the spatial domain operator as the mesh point to be calculated Boundary's extreme value；Other mesh points in first mesh point are traveled through, generate first mesh point on the gravity anomaly side of preset direction Boundary's extreme value information.
- 5. a kind of device of gravity anomaly border enhancing, it is characterised in that described device includes：Gridding unit, gridding is carried out for obtaining the gravity anomaly data in target work area, and to the gravity anomaly data；Boundary extremum computing unit, for extracting gravity anomaly of first mesh point in preset direction one by one using spatial domain operator Boundary extremum information；Enhancement value computing unit, for the gravity anomaly boundary extremum information based on first mesh point, is established with the second net Centered on lattice point, size be the first preset value enhancement value calculation window；It is additionally operable to extraction and the enhancement value calculation window phase Adjacent eight directions, adjacent data window that size is first preset value；It is additionally operable to calculate the enhancement value calculating respectively The correlation coefficient value of window and the adjacent data window in adjacent eight directions；It is additionally operable to eight correlation coefficient value Relevant enhancement value of the eight evolution values of product as second mesh point.
- 6. the device of a kind of gravity anomaly border enhancing according to claim 5, it is characterised in that described device is also wrapped Include：Denoising unit, for removing the noise in the gravity anomaly data using low pass filtering method.
- 7. the device of a kind of gravity anomaly border enhancing according to claim 5, it is characterised in that described device is also wrapped Include：Whether enhancement value judging unit, the relevant enhancement value for judging second mesh point are less than predetermined threshold value；Secondary enhancement unit, for when the enhancement value judging unit judging result for be when, based on first mesh point Relevant enhancing value information, calculate the secondary relevant enhancement value of second mesh point.
- A kind of 8. device of gravity anomaly border enhancing according to claim 5, it is characterised in that the boundary extremum meter Calculating unit includes：Operator acquiring unit, for obtaining spatial domain operator of the target work area on preset direction, the spatial domain operator Size be the second preset value；Extreme value window extraction unit, for extracting in the first mesh point centered on mesh point to be calculated, size it is default for second The extreme value calculation window of value；Extreme value computing unit, for using the product of the extreme value calculation window and the spatial domain operator as the net to be calculated The gravity anomaly boundary extremum of lattice point；Computing unit is traveled through, for traveling through other mesh points in first mesh point, generates first mesh point pre- The gravity anomaly boundary extremum information of set direction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610197300.2A CN105891889B (en) | 2016-03-31 | 2016-03-31 | A kind of method and device of gravity anomaly border enhancing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610197300.2A CN105891889B (en) | 2016-03-31 | 2016-03-31 | A kind of method and device of gravity anomaly border enhancing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105891889A CN105891889A (en) | 2016-08-24 |
| CN105891889B true CN105891889B (en) | 2018-05-04 |
Family
ID=57011812
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610197300.2A Active CN105891889B (en) | 2016-03-31 | 2016-03-31 | A kind of method and device of gravity anomaly border enhancing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105891889B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109283589B (en) * | 2018-08-20 | 2019-12-06 | 桂林理工大学 | Method for acquiring horizontal component of gravity field |
| CN111007571B (en) * | 2019-11-28 | 2022-05-13 | 吉林大学 | Aeromagnetic data geologic body boundary identification method based on three-dimensional structure tensor |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102937725A (en) * | 2012-11-12 | 2013-02-20 | 中国科学院地质与地球物理研究所 | Potential field anomaly edge enhancement method based on transition zone and phase congruency |
| WO2013150464A1 (en) * | 2012-04-04 | 2013-10-10 | Geco Technology B.V. | Seismic data processing with frequency diverse de-aliasing filtering |
| CN103439755A (en) * | 2013-09-05 | 2013-12-11 | 杨洋 | Method for identifying geologic-body boundaries through use of potential-field high-order horizontal gradient model |
| CN103513288A (en) * | 2012-06-21 | 2014-01-15 | 中国石油天然气集团公司 | Compensation direction filtering method of two-dimensional grid data |
| CN104280784A (en) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | Method for recognizing small fracture through gravity |
| CN104965232A (en) * | 2015-06-04 | 2015-10-07 | 中国地质科学院矿产资源研究所 | Automatic extraction method of magnetic structure grillwork in low latitude region |
-
2016
- 2016-03-31 CN CN201610197300.2A patent/CN105891889B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013150464A1 (en) * | 2012-04-04 | 2013-10-10 | Geco Technology B.V. | Seismic data processing with frequency diverse de-aliasing filtering |
| CN103513288A (en) * | 2012-06-21 | 2014-01-15 | 中国石油天然气集团公司 | Compensation direction filtering method of two-dimensional grid data |
| CN102937725A (en) * | 2012-11-12 | 2013-02-20 | 中国科学院地质与地球物理研究所 | Potential field anomaly edge enhancement method based on transition zone and phase congruency |
| CN104280784A (en) * | 2013-07-12 | 2015-01-14 | 中国石油天然气集团公司 | Method for recognizing small fracture through gravity |
| CN103439755A (en) * | 2013-09-05 | 2013-12-11 | 杨洋 | Method for identifying geologic-body boundaries through use of potential-field high-order horizontal gradient model |
| CN104965232A (en) * | 2015-06-04 | 2015-10-07 | 中国地质科学院矿产资源研究所 | Automatic extraction method of magnetic structure grillwork in low latitude region |
Non-Patent Citations (1)
| Title |
|---|
| 重力异常图中的边界增强和提取技术;肖锋 等;《吉林大学学报(地球科学版)》;20110731;第41卷(第4期);第1198、1201最后一段及1202页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105891889A (en) | 2016-08-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| He et al. | An improved fractional-order differentiation model for image denoising | |
| CN103116875B (en) | Self-adaptation bilateral filtering image de-noising method | |
| Looney et al. | Multiscale image fusion using complex extensions of EMD | |
| CN104793253B (en) | Aviation electromagnetic data de-noising method based on mathematical morphology | |
| CN106441288A (en) | Adaptive wavelet denoising method for accelerometer | |
| JP2013235594A5 (en) | ||
| CN104933679B (en) | A kind of method and its correspondence system of enlarged drawing | |
| Song et al. | A spatially adaptive total variation regularization method for electrical resistance tomography | |
| CN105891889B (en) | A kind of method and device of gravity anomaly border enhancing | |
| CN103208104B (en) | A kind of image de-noising method based on nonlocal theory | |
| CN103065291A (en) | Image fusion method based on promoting wavelet transform and correlation of pixel regions | |
| CN109885903A (en) | A kind of ground nuclear magnetic resonance signal peaks noise remove method based on model | |
| CN109903254A (en) | Based on the improved bilateral filtering method of Poisson's kernel | |
| CN103455986B (en) | Random noise point detecting method based on fractional order differential gradient | |
| CN106682391B (en) | The screening of actual measurement stormy waves situation and theoretical spectrum approximating method based on P-M spectrum and JONSWAP spectrum | |
| CN109116423A (en) | A kind of diffraction multiple wave drawing method and device | |
| CN106814394B (en) | Station formula wave detector combines low frequency restoration method with conventional wave detector | |
| Li et al. | A novel medical image fusion approach based on nonsubsampled shearlet transform | |
| CN110335287A (en) | The extracting method and device of Architectural drawing data | |
| Kim et al. | Demosaicking using geometric duality and dilated directional differentiation | |
| Gupta | Enhanced edge detection technique for satellite images | |
| CN104156925A (en) | Processing method and system used for carrying out speckle removing and boundary enhancing on ultrasound image | |
| CN107527057B (en) | Wind speed and power abnormal data eliminating method and device | |
| CN105260249A (en) | Method for extracting calculation intensity features of spatial calculation domain | |
| He et al. | An automatic abrupt information extraction method based on singular value decomposition and higher-order statistics |
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 |