CN105425307B - A kind of method of 1D potential field anomalies curve structure 2D potential field anomaly sections - Google Patents
A kind of method of 1D potential field anomalies curve structure 2D potential field anomaly sections Download PDFInfo
- Publication number
- CN105425307B CN105425307B CN201510870525.5A CN201510870525A CN105425307B CN 105425307 B CN105425307 B CN 105425307B CN 201510870525 A CN201510870525 A CN 201510870525A CN 105425307 B CN105425307 B CN 105425307B
- Authority
- CN
- China
- Prior art keywords
- mrow
- potential field
- dimensional
- anomaly
- munderover
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 13
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 20
- 230000002159 abnormal effect Effects 0.000 claims abstract description 14
- 230000005484 gravity Effects 0.000 claims description 24
- 238000005070 sampling Methods 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims 2
- 238000000605 extraction Methods 0.000 abstract 1
- 239000000284 extract Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Complex Calculations (AREA)
- Peptides Or Proteins (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
A kind of method of 1D potential field anomalies curve structure 2D potential field anomaly sections, step are:The first step, it will containmIndividual sampled point, put away from for △xOne-dimensional potential field anomaly digital independent to one-dimension arrayfIn;Second step, it is rightfCarry out one-dimensional discrete Multiscale Wavelet Decomposition;3rd step, extraction is every respectivelyOneThe high fdrequency component of multi-scale wavelet details, and carry out one-dimensional discrete inverse wavelet transform;4th step, scale dimension is equivalent to intend depth dimension, and by the potential field anomaly data that multiple dimensioned high fdrequency component is set up by the order of small yardstick to large scale, forms 2-D dataF;5th step, to 2-D dataFGridding and interpolation processing are carried out, is then depicted as Two-dimensional Potential Field abnormal profile figure, just realizes one-dimensional potential field anomaly curve structure Two-dimensional Potential Field abnormal profile, the two dimensional cross-section of acquisition has that abnormal information displaying is directly perceived, is easy to the advantages of geologic interpretation.
Description
Technical field
The invention belongs to potential field data processing technical field, more particularly to a kind of 1D potential field anomalies curve structure 2D potential fields are different
The method of normal section.
Background technology
Potential field includes gravitational field, magnetic field and electric field.The exception of potential field is integrated from different depth, different scale field source body
Response and the result being superimposed jointly.Itself limited by potential field exploitation method theory, target seeker gyro collection in field is difficult to obtain deeply
(dimension) information is spent, causes the geologic interpretation of potential field anomaly not directly perceived enough, so as to increase the difficulty of potential field data geologic interpretation work
Degree, and reduce the accuracy of geologic interpretation.
The content of the invention
In order to overcome the above-mentioned deficiencies of the prior art, it is an object of the invention to propose a kind of 1D potential field anomalies curve structure
The method for building 2D potential field anomaly sections, this method have abnormal information displaying directly perceived, the advantages of being easy to geologic interpretation.
To achieve these goals, the technical solution adopted by the present invention is:A kind of 1D potential field anomalies curve builds 2D potential fields
The method of abnormal profile, comprises the following steps:
The first step, m sampled point will be contained, put away from the one-dimensional potential field anomaly digital independent for △ x into one-dimension array f;
Second step, one-dimensional discrete Multiscale Wavelet Decomposition direct transform, one-dimensional discrete multi-scale wavelet are carried out to one-dimension array f
Decomposing direct transform is:
In formula, ψ is referred to as wavelet (also referred to as morther wavelet);ψ*For ψ conjugate function;A represents coefficient of dilatation, reflects specific
The width (also referred to as yardstick) of basic function;B represents translation coefficient, specifies the position translated along x-axis, and
And make a0=2, b0=1, j ∈ Z, k ∈ Z;△ x for point away from;I=0,2 ..., m-1, it is sampling sequence number;
3rd step, extracts the high fdrequency component of each multi-scale wavelet details respectively, and carries out one-dimensional discrete multi-scale wavelet point
Inverse transformation is solved, one-dimensional discrete Multiscale Wavelet Decomposition contravariant is changed to:
In formula, ψ is referred to as wavelet (also referred to as morther wavelet);J ∈ Z, k ∈ Z;X=i △ x, sat for data sampling point
Mark, △ x for point away from;I=0,2 ..., m-1, it is sampling sequence number;
4th step, scale dimension is equivalent to intend depth dimension, and multiple dimensioned height is set up by by the order of small yardstick to large scale
The potential field anomaly data of frequency component, form 2-D data F;
5th step, gridding and interpolation processing are carried out to 2-D data F, are then depicted as Two-dimensional Potential Field abnormal profile
Figure, just realize one-dimensional two-dimentional (2D) potential field anomaly section of (1D) potential field anomaly curve structure.
The beneficial effects of the invention are as follows:
This method is of a relatively high based on potential field anomaly frequency caused by shallow-layer field source body, and potential field caused by deep layer field source body
Abnormal frequency is relatively low, and increases with depth, and the frequency of potential field anomaly is gradually being reduced to supposed premise, by it is one-dimensional from
Multiscale Wavelet Decomposition technology is dissipated, the potential field anomaly from different depth, different field sources is subjected to multi-resolution decomposition, then chi
Degree dimension is equivalent to intend depth dimension, and is combined potential field anomaly by the order by small yardstick to large scale, so that will be one-dimensional
(1D) potential field anomaly curve, which is expanded, constructs two-dimentional (2D) potential field anomaly section, and this two-dimentional (2D) section can be directly used for geology solution
Release;This method has abnormal information displaying directly perceived, the advantages of being easy to geologic interpretation.
Brief description of the drawings
Fig. 1 is actual one-dimensional (1D) bouguer gravity anomaly curve map.
Fig. 2 is the Scale Decomposition of one-dimensional discrete small echo 13 of the present invention and reconstructs the different scale high fdrequency component Bouguer gravity of acquisition
Abnormal curve figure.
Fig. 3 is two dimension (2D) Bouguer weight that the one-dimensional Scale Decomposition of (1D) bouguer gravity anomaly curve 13 structure of the present invention obtains
Power abnormal profile figure.
Embodiment
The present invention is described in more detail with reference to the accompanying drawings and examples.
Referring to Fig. 1,2,3, a kind of method of 1D potential field anomalies curve structure 2D potential field anomaly sections, comprise the following steps:
The first step, m sampled point will be contained, put away from the one-dimensional potential field anomaly digital independent for △ x into one-dimension array f,
Referring to Fig. 1;
Second step, one-dimensional discrete Multiscale Wavelet Decomposition direct transform, one-dimensional discrete multi-scale wavelet are carried out to one-dimension array f
Decomposing direct transform is:
In formula, ψ is referred to as wavelet (also referred to as morther wavelet);ψ*For ψ conjugate function;A represents coefficient of dilatation, reflects specific
The width (also referred to as yardstick) of basic function;B represents translation coefficient, specifies the position translated along x-axis, and
And make a0=2, b0=1, j ∈ Z, k ∈ Z;△ x for point away from;I=0,2 ..., m-1, it is sampling sequence number;
3rd step, extracts the high fdrequency component of each multi-scale wavelet details respectively, and carries out one-dimensional discrete multi-scale wavelet point
Inverse transformation is solved, one-dimensional discrete Multiscale Wavelet Decomposition contravariant is changed to:
In formula, ψ is referred to as wavelet (also referred to as morther wavelet);J ∈ Z, k ∈ Z;X=i △ x, sat for data sampling point
Mark, △ x for point away from;I=0,2 ..., m-1, it is sampling sequence number, referring to Fig. 2;
4th step, scale dimension is equivalent to intend depth dimension, and multiple dimensioned height is set up by by the order of small yardstick to large scale
The potential field anomaly data of frequency component, form 2-D data F;
5th step, gridding and interpolation processing are carried out to 2-D data F, are then depicted as Two-dimensional Potential Field abnormal profile
Figure, one-dimensional two-dimentional (2D) potential field anomaly section of (1D) potential field anomaly curve structure is just realized, referring to Fig. 3.
Embodiment
155 sampled points will be contained, put and implement step away from the one-dimensional bouguer gravity anomaly data instance explanation of actual measurement for 2km
Suddenly:
The first step, 155 sampled points will be contained, put away from the one-dimensional bouguer gravity anomaly digital independent for 2km a to dimension
In group f, Fig. 1;
Second step, the Scale Decomposition direct transform of one-dimensional discrete small echo 13, the yardstick of one-dimensional discrete small echo 13 are carried out to one-dimension array f
Decomposing direct transform is:
In formula, ψ is referred to as wavelet (also referred to as morther wavelet);ψ*For ψ conjugate function;A represents coefficient of dilatation, reflects specific
The width (also referred to as yardstick) of basic function;B represents translation coefficient, specifies the position translated along x-axis, and
And make a0=2, b0=1, j=0,1 ..., 12, k ∈ Z;△ x=2km, for point away from;I=0,1 ..., 154, it is sampling sequence number;
3rd step, extracts the high fdrequency component of each multi-scale wavelet details respectively, and carries out one-dimensional discrete multi-scale wavelet point
Inverse transformation is solved, one-dimensional discrete Multiscale Wavelet Decomposition contravariant is changed to:
In formula, ψ is referred to as wavelet (also referred to as morther wavelet);J=0,1 ..., 12, k ∈ Z;X=i △ x, are adopted for data
Sampling point coordinate;△ x=2km, for point away from;I=0,1 ..., 154, it is sampling sequence number, Fig. 2;
4th step, scale dimension is equivalent to intend depth dimension, and it is high by 13 yardsticks are set up by the order of small yardstick to large scale
The bouguer gravity anomaly data of frequency component, form 2-D data F;
5th step, gridding and interpolation processing are carried out to 2-D data F, are then depicted as two-dimentional bouguer gravity anomaly
Profile, just realize one-dimensional two-dimentional (2D) bouguer gravity anomaly section of (1D) bouguer gravity anomaly curve structure, Fig. 3.
Case effect explanation:
Fig. 1 is actual one-dimensional (1D) bouguer gravity anomaly curve, and wherein abscissa is sampled point sequence number, and ordinate is Bouguer
GRAVITY ANOMALIES △ g, the curve is more apparent along the off-note performance on line direction (transverse direction), but lacks the letter of Depth Domain
Which type of geologic body breath, can not clearly judge by, where cause these bouguer gravity anomalies.
Each curve represents the bouguer gravity anomaly curve of different scale in Fig. 2, and wherein abscissa is sampled point sequence number,
Ordinate is bouguer gravity anomaly value △ g.
Fig. 3 is to be obtained by one-dimensional (1D) bouguer gravity anomaly curve by one-dimensional discrete Multiscale Wavelet Decomposition and reconstruct
Two-dimentional (2D) bouguer gravity anomaly section, wherein abscissa are sampled point sequence number, and ordinate is yardstick (plan depth), and colour code represents
Bouguer gravity anomaly value △ g at two dimensional cross-section diverse location, due to the present invention by scale dimension be equivalent to intend depth dimension, so as to open up
Depth domain information is put on display, can effectively identify bouguer gravity anomaly in two dimension from the exceptional value size of the profile and abnormal profile
Distribution characteristics on section, and qualitatively judge out and cause the geologic body of bouguer gravity anomaly to be high density body or low-density.
Claims (2)
- A kind of 1. method of 1D potential field anomalies curve structure 2D potential field anomaly sections, it is characterised in that comprise the following steps:The first step, m sampled point will be contained, put away from the one-dimensional potential field anomaly digital independent for △ x into one-dimension array f;Second step, one-dimensional discrete Multiscale Wavelet Decomposition direct transform, one-dimensional discrete Multiscale Wavelet Decomposition are carried out to one-dimension array f Direct transform is:<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>W</mi> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>b</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <mo>|</mo> <mi>a</mi> <mo>|</mo> </mrow> </msqrt> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>m</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msup> <mi>&psi;</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>i</mi> <mo>&CenterDot;</mo> <mi>&Delta;</mi> <mi>x</mi> <mo>-</mo> <mi>b</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>&CenterDot;</mo> <mi>&Delta;</mi> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>a</mi> <mo>&NotEqual;</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced>In formula, ψ is referred to as wavelet;ψ*For ψ conjugate function;A represents coefficient of dilatation, reflects the width of specific basic function;B tables Show translation coefficient, specify the position translated along x-axis, andAnd make a0=2, b0=1, j ∈ Z, k ∈ Z;△x For point away from;I=0,1 ..., m-1, it is sampling sequence number;3rd step, the high fdrequency component of each multi-scale wavelet details is extracted respectively, and it is anti-to carry out one-dimensional discrete Multiscale Wavelet Decomposition Conversion, one-dimensional discrete Multiscale Wavelet Decomposition contravariant are changed to:<mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mo>-</mo> <mi>&infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&infin;</mi> </mrow> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mo>-</mo> <mi>&infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&infin;</mi> </mrow> </munderover> <mi>W</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>,</mo> <mi>k</mi> <mo>)</mo> </mrow> <mi>&psi;</mi> <mrow> <mo>(</mo> <msup> <mn>2</mn> <mrow> <mo>-</mo> <mi>j</mi> </mrow> </msup> <mi>x</mi> <mo>-</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow>In formula, ψ is referred to as wavelet;J ∈ Z, k ∈ Z;X=i △ x, it is data sampling point coordinates;△ x for point away from;I=0, 1 ..., m-1, it is sampling sequence number;4th step, scale dimension is equivalent to intend depth dimension, and multiple dimensioned high frequency division is set up by by the order of small yardstick to large scale The potential field anomaly data of amount, form 2-D data F;5th steps, gridding and interpolation processing are carried out to 2-D data F, are then depicted as Two-dimensional Potential Field abnormal profile figure, Just one-dimensional potential field anomaly curve structure Two-dimensional Potential Field abnormal profile is realized.
- 2. a kind of method of 1D potential field anomalies curve structure 2D potential field anomaly sections according to claim 1, its feature exist In comprising the following steps:The first step, 155 sampled points will be contained, put away from the one-dimensional bouguer gravity anomaly digital independent for 2km to one-dimension array f In;Second step, the Scale Decomposition direct transform of one-dimensional discrete small echo 13, the Scale Decomposition of one-dimensional discrete small echo 13 are carried out to one-dimension array f Direct transform is:<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>W</mi> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>b</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mrow> <mo>|</mo> <mi>a</mi> <mo>|</mo> </mrow> </msqrt> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>154</mn> </munderover> <msup> <mi>&psi;</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>i</mi> <mo>&CenterDot;</mo> <mi>&Delta;</mi> <mi>x</mi> <mo>-</mo> <mi>b</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>&CenterDot;</mo> <mi>&Delta;</mi> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mi>a</mi> <mo>&NotEqual;</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced>In formula, ψ is referred to as wavelet;ψ*For ψ conjugate function;A represents coefficient of dilatation, reflects the width of specific basic function;B tables Show translation coefficient, specify the position translated along x-axis, andAnd make a0=2, b0=1, j=0,1 ..., 12, k∈Z;△ x=2km, for point away from;I=0,1 ..., 154, it is sampling sequence number;3rd step, the high fdrequency component of each multi-scale wavelet details is extracted respectively, and it is anti-to carry out one-dimensional discrete Multiscale Wavelet Decomposition Conversion, one-dimensional discrete Multiscale Wavelet Decomposition contravariant are changed to:<mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>12</mn> </munderover> <munderover> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mo>-</mo> <mi>&infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&infin;</mi> </mrow> </munderover> <mi>W</mi> <mrow> <mo>(</mo> <mi>j</mi> <mo>,</mo> <mi>k</mi> <mo>)</mo> </mrow> <mi>&psi;</mi> <mrow> <mo>(</mo> <msup> <mn>2</mn> <mrow> <mo>-</mo> <mi>j</mi> </mrow> </msup> <mi>x</mi> <mo>-</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow>In formula, ψ is referred to as wavelet;J=0,1 ..., 12, k ∈ Z;X=i △ x, it is data sampling point coordinates;△ x=2km, For point away from;I=0,1 ..., 154, it is sampling sequence number;4th step, scale dimension is equivalent to intend depth dimension, and the high frequency division of 13 yardsticks is set up by by the order of small yardstick to large scale The bouguer gravity anomaly data of amount, form 2-D data F;5th step, gridding and interpolation processing are carried out to 2-D data F, are then depicted as two-dimentional bouguer gravity anomaly section Figure, just realize one-dimensional bouguer gravity anomaly curve and build two-dimentional bouguer gravity anomaly section.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510870525.5A CN105425307B (en) | 2015-12-01 | 2015-12-01 | A kind of method of 1D potential field anomalies curve structure 2D potential field anomaly sections |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510870525.5A CN105425307B (en) | 2015-12-01 | 2015-12-01 | A kind of method of 1D potential field anomalies curve structure 2D potential field anomaly sections |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105425307A CN105425307A (en) | 2016-03-23 |
CN105425307B true CN105425307B (en) | 2018-01-19 |
Family
ID=55503631
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510870525.5A Expired - Fee Related CN105425307B (en) | 2015-12-01 | 2015-12-01 | A kind of method of 1D potential field anomalies curve structure 2D potential field anomaly sections |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105425307B (en) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5673191A (en) * | 1995-04-10 | 1997-09-30 | Atlantic Richfield Company | Method and apparatus for identifying geological structures using wavelet analysis of potential fields |
US6691039B1 (en) * | 2002-08-30 | 2004-02-10 | John M. Robinson | Removal of noise from seismic data using improved radon transformations |
CN103245972B (en) * | 2012-02-07 | 2016-06-08 | 中国石油天然气集团公司 | A kind of method determining complex geological structure in two-dimensional space |
CN102809762B (en) * | 2012-08-13 | 2014-09-03 | 成都理工大学 | Reservoir imaging technique based on full-frequency-band seismic information mining |
-
2015
- 2015-12-01 CN CN201510870525.5A patent/CN105425307B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105425307A (en) | 2016-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110058317B (en) | Aviation transient electromagnetic data and aviation magnetotelluric data joint inversion method | |
Wang et al. | A wavelet-based baseline drift correction method for grounded electrical source airborne transient electromagnetic signals | |
CN103654789A (en) | Fast magnetic resonance parametric imaging method and system | |
CN104280784A (en) | Method for recognizing small fracture through gravity | |
CN105353428A (en) | Ground-air cooperative electromagnetic data correction method for magnetic field prolongation of ground reference region | |
CN101504659B (en) | Method for extracting geoscience spatial information based on generalized self-similarity principle | |
CN105372718B (en) | A kind of method that 2D potential fields build 3D potential fields | |
CN102944905B (en) | A kind of gravity-magnetic anomaly disposal route based on direction wavelet analysis | |
CN105629317A (en) | Magnetotelluric noise suppressing method based on intersite transfer function | |
Yang et al. | 3D inversion of total magnetic intensity data for time-domain EM at the Lalor massive sulphide deposit | |
CN104732076A (en) | Method for extracting energy trace characteristic of side channel | |
CN103675926A (en) | Conductivity-depth conversion method for aviation transient electromagnetic data | |
CN108008456B (en) | A method of delineation mesothermal gold deposits deep three-dimensional emphasis U metallogeny Favourable Target Areas | |
CN105425306B (en) | The method that 3D potential fields are built based on 2D discrete wavelets multi-resolution decomposition | |
Zhang et al. | 3D inversion of large-scale frequency-domain airborne electromagnetic data using unstructured local mesh | |
CN105425307B (en) | A kind of method of 1D potential field anomalies curve structure 2D potential field anomaly sections | |
CN104598672A (en) | Method for computing electromagnetic response error caused by emission source attitude change | |
CN106646638A (en) | Solid source transient electromagnetic three-dimensional tunnel advanced forecasting method | |
CN107356963A (en) | A kind of adaptive seismic signal coherence body property analysis method of data-driven | |
Wang et al. | A constrained scheme for high precision downward continuation of potential field data | |
CN114488327B (en) | Horizontal magnetic field and vertical magnetic field combined measurement method in well based on ground base point | |
CN112882124B (en) | Method for delineating beneficial parts of three-dimensional uranium ore formation | |
CN204731426U (en) | A kind of very low frequency (VLF) visits water instrument | |
Lawal et al. | Wavelet analysis of High resolution Aeromagnetic (HRAM) data over part of Chad Basin, Nigeria | |
CN115421212B (en) | Synchronous processing method and system for heavy magnetic map complement and denoising |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180119 Termination date: 20181201 |