In the complex magnetic anomalies under undulating terrain, multiple source superposition and significant noise conditions, the “noise avoidance inversion” depends on three experiments. The displacement experiment accurately solves the prediction problem of the form of the target magnetic source, finds the matching relationship with the maximum thickness of the target magnetic source, and the noise change experiment demonstrates the reliability of the noise avoidance inversion. It can be seen that “noise avoidance inversion” creates a “chaotic image back to true” inversion route for complex magnetic anomalies. Noise avoidance inversion of multi-source superimposed complex magnetic anomalies in Erhai Lake is an example.
In the topic of complex magnetic anomaly, it is difficult to predict the magnetic source morphology of the target anomalies, and then comprehensively explain the abnormalities in the inversion field of complex magnetic anomalies. This is difficult to mention, should be the wish of many scholars.
Studying the target anomalies in avoidance noise inversion in complex abnormalities seems an inconceivable puzzle. However, when a “recognition function” can identify the abnormal magnetic source form of the target on the spot, the bottleneck problem of solving the problem will make substantial progress. This is the original idea of “noise avoidance reaction”.
Variable model is the model basis of noise avoidance inversion; elevation tracking is the method support of noise avoidance inversion; remutation experiment is the experimental basis of noise avoidance inversion discovery, “recognition function” and noise avoidance.
Theoretically, variable models can be equivalent to any magnet model of independent morphology and produce anomalies with unique features on arbitrary undulating topography. Importantly, the following formula can easily give three models of the thickening of the magnet, the book, the unchanged tip body, the lens body, and the plate body. Under the rolling terrain conditions, the magnetic field expression of lens shape, plate body, pointed extinguishing body and other forms can be written as:
Δ T x y = 1 n ∑ i = 1 n 2 b i J S sin β 1 − cos 2 I 0 cos 2 A F i ( h x y x , y ± L i , l i ) (1)
F i ( h x y , x , y ± L i , l i ) = f ( h x y , x , y + L i , l i ) − f ( h x y , x , y − L i , l i )
f ( h x y , x , y + L ) = L z c 2 + L z 2 ( x − l cos β ) cos ξ − ( h x y + l sin β ) cos ξ ( x − l cos β ) 2 + ( h x y + l sin β ) 2 − L z ( x sin ξ − h x y cos ξ ) ( D 2 + L 2 ) D 2 − 1 B 2 + L z 2 ( P + 1 c 2 + L z 2 − P D 2 + L z 2 ) ( ( cos ξ cos β − cos ξ − cos 2 ξ + K 2 ) L z − 2 k sin ( i s − β ) B ) + 2 K cos ( i s − β ) ( 1 c 2 + L z 2 − 1 D 2 + L z 2 )
where: L z = Y + L ; L f = Y − L ; P = h x y sin β − x cos β ; B = x sin β + h x y cos β ; D 2 = x 2 + h x y 2 ; c 2 = ( x − l cos β ) 2 + ( h x y + l sin β ) 2 ; K = cos i s ctg A ; ξ = 2 i s − β ; i s = arctg ( tg I 0 csc A ) ; h x y = h + H ( x , y ) ; b i = b 1 , b 2 , b 3 , ⋯ , b n ; L i = L 1 , L 2 , L 3 , ⋯ , L n ; l i = l 1 , l 2 , l 3 , ⋯ , l n ; H ( x , y ) for height data with zero origin; J s , effective magnetization within the profile; ξ , magnetization feature angle; L, half length of the direction of the magnetic source; l, Magnetic source depth elongation; b, half thickness of magnetic source; β , magnetic source inclination; I 0 , geomagnetic inclination; i s , effective magnetization inclination within the profile; A, magnetic source azimuth.
Using the above variable model formula, different sizes of book plates are selected, parallel superimposed, constructed in the profile: the “pointed body with thick top and narrow bottom”; narrow “lens body” up to the bottom and thick in the middle; up and down equally thick “thick plate body”. The three types of models are the basis of forming the “model three methods” to identify magnet morphology.
“Model three methods” refers to the combination of methods using the exclusive interpretation method of three models to identify the same target abnormality in complex abnormalities.
Elevation tracking is to use the theoretical abnormalities with the terrain influence to explain the field abnormalities affected by the terrain, therefore, for the terrain impact, no need to handle, but the use, is also resolved.
Elevation tracking, “tracking” refers to the magnetic source with different buried depths in the profile. It can automatically use the elevation of the surface projection on the magnet as zero, and take the buried depth scale as the “relative elevation” explaining the corresponding magnetic source anomalies, which can virtually resolve the impact of terrain on the anomalous inversion process.
Elevation tracking formula for resolving the terrain influence for noise avoidance and inversion:
H m ( i ) = H ( i ) − H ( I N T N 2 + J J ) (2)
In type, H m ( i ) , elevation tracking array; H ( i ) , section elevation array; H ( I N T N 2 + J J ) , exception origin elevation; J J , number of points where the origin deviates from the point in the profile.
For the three replacement experimental sections, plate abnormal interpretation method, lens abnormal interpretation method (model three methods) were used, and three groups and nine interpretation results were obtained.
Using the integration difference of magnetic moment and outlier absolute value of nine interpretation results in the front and after three groups, through alignment, screening and combination, and finally found an “identification function” that can correctly distinguish the abnormal magnetic source after the target anomalous magnetic source morphology is replaced φ .
φ = Δ T 00 − Δ T S Y Δ T 00 + Δ T D J M S M B 3 (3)
Δ T 00 , integrvalues of complex outliers; Δ T S Y , the remaining outliers are absolute integral value; Δ T D J , superimposed outliers absolute values; M S , the model method explains the magnetic moments in the target anomaly results; M B , the bookkeeping method explains the magnetic moment results of the target exception.
The Recognition Parameter consists of two information: curve ratio and magnetic moment ratio.
Curve ratio refers to the three types of curves obtained by the model three-way method, where the superposition anomalies involved include the superimposed anomalies prior to the interpretation of the target anomalies. The remaining anomalies refer to the difference between the complex outlier absolute integral value and the superimposed outlier absolute integral value
Magnetic moment ratio refers to the magnetic moment ratio obtained by model 3 method and accounting interpretation method explaining target anomalies respectively.
The recognition function interpretation of the structure is the reverse difference between the curve and the magnetic moment of the three model methods to explain the target anomaly. Only the recognition function of this combination is chosen to directly discriminate the anomalous magnetic source morphology of the target.
The results in
| Replacement experiments | recognition methods | Identify parameters | Identify results |
|---|---|---|---|
| Abnormal thick plate body replacement “D” abnormal | Abnormal thick plate interpretation | 23.85 | Abnormal thickness plate |
| Lens body abnormal interpretation | 22.98 | ||
| Interpretation method of cone | 22.10 | ||
| Lens body abnormal replacement “D” abnormal | Abnormal thick plate interpretation | 6.79 | Lens body abnormal |
| Lens body abnormal interpretation | 17.03 | ||
| Interpretation method of cone | 9.52 | ||
| Conical abnormalities replacement “D” abnormal | Abnormal thick plate interpretation | 22.84 | Conical abnormalities |
| Lens body abnormal interpretation | 20.41 | ||
| Interpretation method of cone | 25.5 |
Factors affecting the three types of model interpretation methods, including topographic influence, superposition effects of anomalies, and high-frequency interference. When these effects explain the substitution anomaly by using the model three method, they are the invariant state of instantaneous “fixed frame”, which is called “curing noise”, and the test of the model three method does not mask the characteristic information of the magnetic source. It is therefore possible to be revealed and distinguished by the recognition function. It is reasonable to argue that the replacement experiment demonstrates the “curing noise”.
At this point, the replacement experiment finally found the way to determine the attribution of the magnetic source category by relying on the model three-method plus recognition function to determine the target anomalies in the measured abnormality.
After the “model three-method-recognition function” solves the target anomalous magnetic source morphology, the interpretation process of complex magnetic anomalies can be completed objectively, and we give close to the actual interpretation results. This is the theoretical basis for noise avoidance inversion to solve the bottleneck problem.
Model three-method-recognition function tests the “D” anomaly in W01, whose field source morphology is pointed body, and the deep source anomaly of G” is lens.
The displacement experiment gives an effective method to identify the morphology of target abnormal magnetic sources in complex abnormalities: “model three methods-parameter identification”. It can be seen that the “replacement experiment” is the experimental theoretical basis for the noise avoidance inversion.
Under the premise of the target abnormal morphology, the complex abnormal inversion process is suspended. At that time, the thickness can be adjusted to observe the changes of the remaining abnormal, and the thickness of the remaining abnormal minimum point is, the maximum thickness of the target abnormal magnetic source. The principle is that the remaining anomaly is only minimized when the magnetic source morphology of the target anomaly is consistent with the thickness.
The so-called noise that explains the target anomalies, includes many factors. During topographic effects, the effects of adjacent abnormal errors, the effects of remaining superimposed anomalies in residual abnormalities that have not been explained, and the high-frequency interference in the presence of shallow layers. These effects are completely different in different stages of the interpretation of complex anomaly profiles, but are “fixed”, and transiently invariant “cured noise”.
W01 explained section “D” abnormal morphology, identified as pointed body, the maximum thickness of the top is 92 m by thickness (see
| Try the thickness m | Residual outliers nT | Identify parameters | Determine thickness results m | Relative error % |
|---|---|---|---|---|
| 20 | 7204.05 | 20.8726 | 92 | 2.17 |
| 40 | 7192.11 | 20.9911 | ||
| 80 | 7142.50 | 21.4822 | ||
| 90 | 7124.26 | 21.6626 | ||
| 92 | 7079.21 | 21.9133 | ||
| 94 | 7108.05 | 21.9111 | ||
| 100 | 7102.51 | 21.8814 | ||
| 110 | 7104.79 | 21.8806 |
The model experiment found an identification function that directly distinguishes the magnetic source morphology attribution of target anomalies in complex anomalies, and gave effective methods to determine the magnetic source thickness.
Model trimethod, used to test target exception in W01 profile “D”, results show that the magnetic source form in the W01 section is the pointed body, through the thickness of the maximum thickness of the pointed body is 92 meters. In the expanded application of deep source anomaly, the results of deep source anomaly “The magnetic source form of G” is the lens body. Thickness is 75 m.
Interference noise of different orders is integrated into complex anomalies to test whether it affects the determination of the magnetic source thickness to measure the varying noise. Whether in noise avoidance inversion, noise can still be considered as a constant “cured noise” without affecting the inversion result.
The comparison
During the process of magnetic source thickness determination, the noise interference is large or small and both transient and does not affect the thickness determination result. Doudoubled noise is incorporated into complex anomalies and its thickness determination does not change. See
It can be seen that the “displacement experiment” utilizes a systematic experiment.
A theoretical and methodological basis for the noise avoidance inversion is provided.
| Try the thickness m | Residual outliers nT | Recognition function | Determine thickness results m | Relative error % |
|---|---|---|---|---|
| 20 | 10941 | 0.18161 | 92 | 2.17 |
| 40 | 10988 | 0.19672 | ||
| 80 | 10930 | 0.27534 | ||
| 90 | 10924 | 0.32121 | ||
| 92 | 10921 | 0.34851 | ||
| 94 | 11006 | 0.34534 | ||
| 100 | 10922 | 0.34113 | ||
| 110 | 10937 | 0.21438 |
As shown from
Importantly, the thickness determination process also demonstrates the complex multi-element noise and does not mask the weak magnetic source thickness information; it demonstrates that the noise avoidance inversion avoids the noise effects and gives reliable inversion results.
Noise avoidance inversion was completed by three stages: “bookboard fitting”, “morphological inversion” and “geological interpretation”.
It was found that independent morphological magnet anomalies with certain buried depth are not significantly different from the profile morphology of the book board anomaly, therefore, the book board anomalies are unified to simulate complex magnetic anomalies.
Under the premise of known latitude and longitude and determination of magnetic inclination angle, the abnormal fitting process can be easily completed according to the magnetic field decomposition, superposition and fitting principle, origin position, magnetic moment size, magnetic source depth, magnet inclination angle, length and delay change (see
Its task: to identify the abnormal magnetic source form of the target through the “identification parameters”, and on the basis of determining the thickness, supplemented by necessary adjustment, and strive to keep the form and thickness of the target abnormal magnetic source close to the actual state, to provide a reliable geophysical foundation for geological analysis.
The optimized interpretation shows that in the profile, the “D” anomalous magnetic source is the pointed body, the upper end is 92 m thick, and the magnetic source of the “G” anomaly is the lens body, with an intermediate thickness of 70 m.
Magnet data explaining the model anomalies, except for the thickness data, came from the bookkeeping fitting phase and were slightly modified if necessary (
Combined with the existing geological data, the morphological inversion results are used to draw the “geological interpretation profile of noise avoidance inversion”. Verification recommendations are also made (
Erhai magnetic anomaly is located at the foot of the Cangshan racetrack, undulating terrain. This anomaly is superimposed by multiple magnetic sources with significant high-frequency interference, and is a typical complex magnetic anomaly (
In the same interpretation section, the abnormalities of different directions have different effective magnetization directions in the section, so they need to be handled differently.
The profiles with two abnormal reverse complexes from north to east show unique platform features. Therefore, in the process of book board fitting, its abnormal direction selection, we must make objective different from other target abnormalities. Through the heterogeneous treatment, the effective magnetization inclination in the profile and the effective magnetization inclination of the adjacent magnetic source. Ectropic treatment results: The abnormal platform features were well fitted (
In the original superposition anomalies of W01, the “D”, “E” anomalies are superimposed as similar to monomer anomalies etc. The absolute value of the strong negative value of the exception “D” is greater than the positive, releasing the “E” exception it overrides. To determine its authenticity, the corresponding “underwater peaks” (relative peaks in the “negative region” are easily found in the negative region below the “anomalous” peaks (
When the anomaly explains that “a, B, C, D, E” reaches “F”, the “residual anomaly” presents a slow positive and negative low frequency feature, which is clearly important information for the presence of a deep source “G” (
The magnetic measurement data noise avoidance inversion results combined with the good prospecting prospect.
The complex morphology should be related to structure and hydrothermal activity after sub-magmatic phase.
The “D” magnet, located at the deep extension bend of the rock body, is a larger so-called pointed fire extinguisher, which is suspected to be an “iron cap”. Coincidentally, in the depths of the rock, there is a deep magnet “G”. Between the two, although there are no magnets, this magnetic-free space, is indeed the most advantageous space to explore medium-heat polymetallic sediments.
The innovative way of “chaos inversion” comes from the establishment and demonstration of the new concept of “curing noise”; relies on the identification function of the target magnetic source morphology and thickness in complex magnetic anomalies; and finally develops the three-step combination of “noise avoidance inversion” through the model three method.
The “three methods” of the variable model is the method of noise avoidance inversion; the “identification parameter” obtained by the replacement experiment is the experimental theoretical basis of “noise avoidance inversion”.
Hidden by many shallow magnet anomalies, buried nearly kilometer deep source abnormality can be found; successfully decompose complex anomalies with overlapping conditions and different directions, and complete the full parameter inversion respectively, which fully shows that noise avoidance inversion has important practical and generalized value.
Although the practice of noise avoidance inversion is insufficient, the advanced idea of simplifying complicated is worth continuing to deepen research.
The authors declare no conflicts of interest regarding the publication of this paper.
Cai, J.X. and Yan, Y.P. (2021) Complex Magnetic Anomalous Chaotic Inversion. International Journal of Geosciences, 12, 1071-1083. https://doi.org/10.4236/ijg.2021.1212057