CN110209989A - A kind of anisotropy singularity exponents calculation method based on spatial weighting technology - Google Patents
A kind of anisotropy singularity exponents calculation method based on spatial weighting technology Download PDFInfo
- Publication number
- CN110209989A CN110209989A CN201910512555.7A CN201910512555A CN110209989A CN 110209989 A CN110209989 A CN 110209989A CN 201910512555 A CN201910512555 A CN 201910512555A CN 110209989 A CN110209989 A CN 110209989A
- Authority
- CN
- China
- Prior art keywords
- anisotropy
- spatial weighting
- singularity
- space
- singularity exponents
- 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.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 38
- 238000005516 engineering process Methods 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 25
- 239000002184 metal Substances 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 2
- 230000001089 mineralizing effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The anisotropy singularity exponents calculation method based on spatial weighting technology that the invention discloses a kind of, is related to modeling technique field；Method includes the following steps: S10, definition space weighted factor；S20 sets the rectangular calculation window of multiple side length sizes for position to be calculated a certain on geochemical map, calculates the anisotropy singularity space size constrained based on spatial weighting；S30 calculates the average content of spatial weighting geochemical elements corresponding to each rectangular calculation window；S40, the linear relationship using least square method, on log-log graph between fitting scatterplot；S50 finds out the slope of fit linear relationship, and thus extrapolates the space anisotropic singularity exponents value of position to be calculated；S60 slides rectangular calculation window to the next position, repeats step S20-S50, until space anisotropic singularity exponents value is calculated at by all positions.The present invention has abandoned the defect that SOM method has ignored space anisotropic when calculating singularity exponents.
Description
Technical field
The present invention relates to modeling technique field, specifically a kind of anisotropy singularity exponents based on spatial weighting technology
Calculation method.
Background technique
Angle based on GEOLOGICAL APPLICATION, this concept of singularity are initially to be proposed by Cheng (2007), and it is defined as
The geological phenomenon or geological process occurred within a narrow time or space with incorrect energy release or substance accumulation
A kind of property.Singularity identifies the part of (such as magmatic hydrotherm mineralizing process) Geological Variable in non-linear geology process exception
It has great significance in the characterization of singular structure feature.
Based on multi-fractal Theory, it is quantitative that the prior art proposes a kind of singular index charting technology based on sliding window
The size of non-linear geological process or the unusual attribute of event is had estimated, referred to as common local singular exponent calculation method
(OSM, OrdinarySingularityMapping).By taking geochemistry data as an example, OSM method calculates the tool of singularity exponents
Body process is as follows:
(1) for a certain position to be calculated, a series of side length sizes are set as rectangular calculation window；
(2) average content of geochemical elements corresponding to each window is calculated；
(3) least square method is utilized, the linear relationship being fitted between scatterplot on log-log graph；
(4) slope of fit linear relationship is found out, and thus extrapolates the singularity exponents value of position to be calculated；
(5) the step of rectangular calculation window arrives the next position, repeats (2)-(5) is slided, up to will be unusual at all positions
Until value is calculated.
Since common local singular exponent calculation method proposes, it is widely used in enhancing and identifying geochemistry
Exception information.But certain limitation that this method is deposited, i.e., its have ignored space anisotropic.
In recent years, space anisotropic is increasingly paid attention in anomalous identification including the importance in Singularity Analysis.Tool
The geochemical anomaly for showing mine meaning is often to have coupling in temporal-spatial evolution with control mine geologic elements, because they are
The enrichment (mine) of element provides the conditions such as necessary space, heat, fluid.The a large amount of geology fact verified control mine
Geologic elements are upper with apparent anisotropic character respectively in space.Therefore, space anisotropic is conducive to identification with bright
The geochemical anomaly information of true geological Significance combines Singularity Analysis with space anisotropic analysis, it is ensured that institute
The abnormal exception in geological Significance than being drawn a circle to approve merely with one of property of delineation is more reasonable.
Summary of the invention
The anisotropy singularity exponents calculation method based on spatial weighting technology that the purpose of the present invention is to provide a kind of,
To solve the above problems.
To achieve the above object, the invention provides the following technical scheme:
A kind of anisotropy singularity exponents calculation method based on spatial weighting technology, comprising the following steps:
S10, definition space weighted factor；
S20 sets the rectangular calculation window of multiple side length sizes for position to be calculated a certain on geochemical map, meter
Calculate the anisotropy singularity space size constrained based on spatial weighting；
S30 calculates the average content of spatial weighting geochemical elements corresponding to each rectangular calculation window；
S40, the linear relationship using least square method, on log-log graph between fitting scatterplot；
S50 finds out the slope of fit linear relationship, and thus extrapolates the space anisotropic singularity of position to be calculated
Index value；
S60 slides rectangular calculation window to the next position, repeats step S20-S50, until space at all positions is each
Until anisotropy singularity exponents value is calculated.
Based on the above technical solution, the present invention also provides following optional technical solutions:
In further embodiment: the spatial weighting factor is w [dBx(ε)](0≤w[dBx(ε)]≤1)。
In further embodiment: the side length size of multiple rectangular calculation windows is εi(εi-1< εi< εi+1, I=1,
2 ...), anisotropy singularity space size calculation formula is
In further embodiment: the calculation formula of the average content of the spatial weighting geochemical elements isIts
Middle μw[Bx(ε)] and ρw[Bx(ε)] it respectively represents included in BxThe total metal quality of spatial weighting and space in (ε) spatial dimension add
Weight average density.
In further embodiment: describedWherein E represent Europe it is several in
Obtain the dimension in space.
In further embodiment: describedIts
Middle E represents the dimension of Euclidean space.
In further embodiment: the slope of the fit linear relationship is ∈ (x) -2.
Compared to the prior art, beneficial effects of the present invention are as follows:
The present invention has abandoned the defect that SOM method has ignored space anisotropic when calculating singularity exponents, utilizes sky
Between weighting technique constrain the shapes and sizes of sliding window, the sliding window of arbitrary shape size can be simulated, to measure surprise
The anisotropic feature of anisotropic spatial model and intensity, for local singular value calculating provide a kind of general computation model and
Method.
Specific embodiment
Each embodiment cited by the present invention is only to illustrate the present invention, is not used to limit the scope of the present invention.To this
Any modification apparent easy to know or change are without departure from spirit and scope of the invention made by invention.
The shortcomings that in order to overcome OSM method to ignore anisotropy when identifying geochemical anomaly, it is new that we have proposed SWSM
Method, and carried out case verification.SWSM method mutually ties OSM and spatial weighting technology (Spatially Weighting, SW)
Close, for calculating space anisotropic singular index, can with comprehensive characterization and measure non-linear mineralizing process it is intrinsic unusual
Property and anisotropy the two critical natures.
Spatial weighting technology (Spatially Weighting, SW) is a kind of highly useful spacial analytical method, its energy
Be enough in quantitatively characterizing have anisotropy spatial model geologic elements, to avoid when calculating certain geological informations extract with
It using incoherent sample or variable, the uncorrelated variables such as uses, SW technology can be used, utilize the space based on geologic property
Weight factor is weighted sample contribution, to calculate related coefficient or covariance matrix.Therefore, SW technology be realize sample by
The spatial data analysis of the spatial characters constraint such as anisotropy of geologic feature control provides a kind of useful tool.
Specifically, SWSM method of the invention the following steps are included:
S10, definition space weighted factor w [dBx(ε)](0≤w[dBx(ε)]≤1)。
S20 sets the rectangular calculation window of multiple side length sizes for position to be calculated a certain on geochemical map, meter
Calculate the anisotropy singularity space size constrained based on spatial weighting.
Wherein rectangular calculation window side length size is εi(εi-1< εi< εi+1, i=1,2 ...), anisotropy singularity is empty
Between size calculation formula be
S30 calculates the average content of spatial weighting geochemical elements corresponding to each rectangular calculation window.
The calculation formula of the average content of spatial weighting geochemical elements is
Wherein μw[Bx(ε)] and ρw[Bx(ε)] it respectively represents included in BxThe total metallic of spatial weighting in (ε) spatial dimension
Amount and Space weighted mean density.
It is describedWherein E represents the dimension of Euclidean space.
It is describedWherein E represents Euclid's sky
Between dimension.
S40, the linear relationship using least square method, on log-log graph between fitting scatterplot.
S50 finds out the slope of fit linear relationship, and thus extrapolates the space anisotropic singularity of position to be calculated
Index value.
The slope of the fit linear relationship is ∈ (x) -2.
S60 slides rectangular calculation window to the next position, repeats step S20-S50, until space at all positions is each
Until anisotropy singularity exponents value is calculated.
It is obvious to a person skilled in the art that invention is not limited to the details of the above exemplary embodiments, Er Qie
In the case where without departing substantially from spirit or essential attributes of the invention, the present invention can be realized in other specific forms.Therefore, no matter
From the point of view of which point, the present embodiments are to be considered as illustrative and not restrictive, and the scope of the present invention is by appended power
Benefit requires rather than above description limits, it is intended that all by what is fallen within the meaning and scope of the equivalent elements of the claims
Variation is included within the present invention.
In addition, it should be understood that although this specification is described in terms of embodiments, but not each embodiment is only wrapped
Containing an independent technical solution, this description of the specification is merely for the sake of clarity, and those skilled in the art should
It considers the specification as a whole, the technical solutions in the various embodiments may also be suitably combined, forms those skilled in the art
The other embodiments being understood that.
Claims (7)
1. a kind of anisotropy singularity exponents calculation method based on spatial weighting technology, which is characterized in that including following step
It is rapid:
S10, definition space weighted factor；
S20 sets the rectangular calculation window of multiple side length sizes for position to be calculated a certain on geochemical map, calculates base
In the anisotropy singularity space size of spatial weighting constraint；
S30 calculates the average content of spatial weighting geochemical elements corresponding to each rectangular calculation window；
S40, the linear relationship using least square method, on log-log graph between fitting scatterplot；
S50 finds out the slope of fit linear relationship, and thus extrapolates the space anisotropic singularity exponents of position to be calculated
Value；
S60 slides rectangular calculation window to the next position, repeats step S20-S50, until by space at all positions respectively to different
Until property singularity exponents value is calculated.
2. the anisotropy singularity exponents calculation method according to claim 1 based on spatial weighting technology, feature
It is, the spatial weighting factor is w [dBx(ε)](0≤w[dBx(ε)]≤1)。
3. the anisotropy singularity exponents calculation method according to claim 2 based on spatial weighting technology, feature
It is, the side length size of multiple rectangular calculation windows is εi(εi-1< εi< ε+1, i=1,2 ...), anisotropy singularity
Space size calculation formula is
4. the anisotropy singularity exponents calculation method according to claim 3 based on spatial weighting technology, feature
It is, the calculation formula of the average content of the spatial weighting geochemical elements is
Wherein μw[Bx(ε)] and ρw[Bx(ε)] it respectively represents included in BxThe total metal quality of spatial weighting in (ε) spatial dimension and
Space weighted mean density.
5. the anisotropy singularity exponents calculation method according to claim 4 based on spatial weighting technology, feature
It is, it is describedWherein E represents the dimension of Euclidean space.
6. the anisotropy singularity exponents calculation method according to claim 5 based on spatial weighting technology, feature
It is, it is describedWherein E represents Euclid's sky
Between dimension.
7. the anisotropy singularity exponents calculation method according to claim 6 based on spatial weighting technology, feature
It is, the slope of the fit linear relationship is ∈ (x) -2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910512555.7A CN110209989A (en) | 2019-06-13 | 2019-06-13 | A kind of anisotropy singularity exponents calculation method based on spatial weighting technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910512555.7A CN110209989A (en) | 2019-06-13 | 2019-06-13 | A kind of anisotropy singularity exponents calculation method based on spatial weighting technology |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110209989A true CN110209989A (en) | 2019-09-06 |
Family
ID=67792702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910512555.7A Pending CN110209989A (en) | 2019-06-13 | 2019-06-13 | A kind of anisotropy singularity exponents calculation method based on spatial weighting technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110209989A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113065091A (en) * | 2021-04-12 | 2021-07-02 | 中国地质科学院地质力学研究所 | Method and device for analyzing anisotropic distribution rule of geological information and electronic equipment |
| CN116756674A (en) * | 2023-06-28 | 2023-09-15 | 中国地质科学院地质力学研究所 | Multi-scale anisotropic anomaly identification method, device, equipment and medium |
| CN117891883A (en) * | 2024-03-14 | 2024-04-16 | 山东观和集团有限公司 | Mineral exploration data optimal storage method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110137565A1 (en) * | 2009-12-08 | 2011-06-09 | Chevron U.S.A. Inc. | System and method for lacunarity analysis |
| WO2015119569A1 (en) * | 2014-02-10 | 2015-08-13 | Cr Development Ab | Method for quantifying isotropic diffusion and/or anisotropic diffusion in a sample |
| CN107480460A (en) * | 2017-08-25 | 2017-12-15 | 中国地质科学院矿产资源研究所 | Multi-element geochemical anomaly identification method for large-area complex background area |
-
2019
- 2019-06-13 CN CN201910512555.7A patent/CN110209989A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110137565A1 (en) * | 2009-12-08 | 2011-06-09 | Chevron U.S.A. Inc. | System and method for lacunarity analysis |
| WO2015119569A1 (en) * | 2014-02-10 | 2015-08-13 | Cr Development Ab | Method for quantifying isotropic diffusion and/or anisotropic diffusion in a sample |
| CN107480460A (en) * | 2017-08-25 | 2017-12-15 | 中国地质科学院矿产资源研究所 | Multi-element geochemical anomaly identification method for large-area complex background area |
Non-Patent Citations (1)
| Title |
|---|
| FAN XIAO等: ""A spatially weighted singularity mapping method applied to identify epithermal Ag and Pb-Zn polymetallic mineralization associated geochemical anomaly in Northwest Zhejiang, China"", 《JOURNAL OF GEOCHEMICAL EXPLORATION》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113065091A (en) * | 2021-04-12 | 2021-07-02 | 中国地质科学院地质力学研究所 | Method and device for analyzing anisotropic distribution rule of geological information and electronic equipment |
| CN113065091B (en) * | 2021-04-12 | 2021-11-05 | 中国地质科学院地质力学研究所 | Method and device for analyzing anisotropic distribution rule of geological information and electronic equipment |
| CN116756674A (en) * | 2023-06-28 | 2023-09-15 | 中国地质科学院地质力学研究所 | Multi-scale anisotropic anomaly identification method, device, equipment and medium |
| CN116756674B (en) * | 2023-06-28 | 2024-02-27 | 中国地质科学院地质力学研究所 | Multi-scale anisotropic anomaly identification method, device, equipment and medium |
| CN117891883A (en) * | 2024-03-14 | 2024-04-16 | 山东观和集团有限公司 | Mineral exploration data optimal storage method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110209989A (en) | A kind of anisotropy singularity exponents calculation method based on spatial weighting technology | |
| CN102930533B (en) | A kind of semi-supervised Hyperspectral imaging dimension reduction method based on improving K-mean cluster | |
| Illian et al. | A toolbox for fitting complex spatial point process models using integrated nested Laplace approximation (INLA) | |
| CN108446711A (en) | A kind of Software Defects Predict Methods based on transfer learning | |
| CN113821986B (en) | Vortex sea surface signal and underwater key parameter inversion method based on random forest | |
| CN110119556B (en) | Space-time evolution analysis method for regional water source conservation function | |
| CN109978872B (en) | White matter microstructure characteristic screening system and method based on white matter fiber tracts | |
| CN107728231A (en) | One kind prediction nuclear magnetic resonance log T2 T2The method of distribution | |
| Yang et al. | Analyses of extreme climate events over China based on CMIP5 historical and future simulations | |
| CN104268581A (en) | Remote sensing sub-pixel map-making method based on integrated pixel level and sub-pixel level spatial correlation characteristics | |
| CN103336305A (en) | Method for dividing petrophysical facies of tight sand reservoir based on grey theory | |
| CN115082439B (en) | Tropical cyclone strength determining method, medium and equipment fused with satellite cloud picture space-time information | |
| CN103902798B (en) | Data preprocessing method | |
| CN102938069A (en) | Pure and mixed pixel automatic classification method based on information entropy | |
| Maskey et al. | Temporal downscaling rainfall and streamflow records through a deterministic fractal geometric approach | |
| CN103268423B (en) | A kind of method that geography phenomenon multi-point simulation spatial scale selects | |
| Lewis et al. | Validation of a hydrodynamic-ecosystem model simulation with time-series data collected in the western English Channel | |
| CN102609721B (en) | Remote sensing image clustering method | |
| Vandeskog et al. | Quantile based modeling of diurnal temperature range with the five‐parameter lambda distribution | |
| CN108090653B (en) | Reservoir type identification method and device for reservoir | |
| CN104573339B (en) | The geologic parameter of shale gas reservoir determines method and apparatus | |
| CN105260944A (en) | Method for calculating statistical line loss based on LSSVM (Least Square Support Vector Machine) algorithm and association rule mining | |
| Jin et al. | A trend analysis of seasonal average temperatures over 40 years in South Korea using Mann-Kendall test and sen's slope | |
| Ranasinghe et al. | Does a more sophisticated storm erosion model improve probabilistic erosion estimates? | |
| CN104809336B (en) | A kind of regional feature sampling approach for considering spatial coherence |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190906 |
|
| RJ01 | Rejection of invention patent application after publication |