Three-Dimensional Magnetic Inversion Based on an Adaptive Quadtree Data Compression

: Three-dimensional magnetic inversion allows the distribution of magnetic parameters to be obtained, and it is an important tool for geological exploration and interpretation. However, because of the redundancy of the data obtained from large-scale investigations or high-density sampling, it is very computationally intensive to use these data for iterative inversion calculations. In this paper, we propose a method for compressing magnetic data by using an adaptive quadtree decomposition method, which divides the two-dimensional data region into four quadrants and progressively subdivides them by recursion until the data in each quadrant meets the regional consistency criterion. The method allows for dense sampling at the abnormal boundaries with large amplitude changes and sparse sampling at regions with small amplitude changes, and achieves the best approximation to the original data with the least amount of data, thus retaining more anomalous information while achieving the purpose of data compression. In addition, assigning values to the data in the quadrants using the averaging method is essentially equivalent to average ﬁltering, which reduces the noise of the magnetic data. Testing the synthetic model and applying the method to mineral exploration a prove that it can e ﬀ ectively compress the magnetic data and greatly improve the computational e ﬃ ciency.


Introduction
The magnetic survey is widely used in geological mapping, geological structure surveys, mineral resources exploration, sedimentary basin basement detection and other fields [1]. Three-dimensional magnetic inversion refers to the computation of the spatial distribution and properties of underground sources from magnetic anomalies of an observation surface. According to the different inversion algorithms, magnetic inversion can be summarized into two categories [2], one is to obtain the equivalent magnetism through inversion calculations to outline the distribution of magnetic anomaly field sources, called the probability tomography [3][4][5]. This is a semi-quantitative inversion interpretation method that is simple, stable, and fast to compute. It can be used to locate the geometric position of the anomaly and solve complex geological problems. Another type of inversion method is to directly calculate the magnetic distribution in the imaging space through linear or nonlinear algorithms to obtain information about the shape, volume, magnetic parameter distribution and properties of anomalies. This is called the physical property inversion method, and it can quantitatively interpret the geometric and physical properties of the target geologic map and provide a strong basis for subsequent geological the size of the kernel function matrix. A quadtree is a tree-like data structure, with four sub-blocks on each node. Quadtree decomposition is the process of dividing the data into four quadrants and then recursively dividing them until each quadrant meets the precision requirement. The concept of the image quadtree was proposed by Klinger and Dyer in the early 1970s [41]. The quadtree is an important method in image processing. It has the advantages of fast decomposition, low time complexity, and flexibility. It is widely used in computer graphics, image segmentation, image fusion, remote sensing image processing, geographic information system, geological modeling, and other fields [42][43][44].
This paper first reviews the development of 3D inversion and then analyzes the solutions proposed by different scholars for large-scale data. Finally, an adaptive quadtree decomposition technology is proposed, which can compress the data while preserving the anomalous information of the observation data, thereby greatly improving the inversion efficiency. The effectiveness of the method is demonstrated by experiments on synthetic models and measured data.

Inversion Method
The inversion is the calculation of the physical or geometric parameters of a subsurface structure or an anomaly from anomalous data on the observation surface. For the linear inversion problem, the subterranean space discretizes into units that satisfy the following Equation (1): In the magnetic inversion, d denotes the vector of observed magnetic data, d = (d 1 , d 2 , . . . , d N ) T , where N is the number of data. m is the vector of model parameters, m = (m 1 , m 2 , . . . , m M ) T , where M is the number of dissected grids, and A denotes the sensitivity matrix of N × M dimensions, expressed as A i,j (i = 1,..., N); j = 1,..., M), which represents the response of the j-th unit at the i-th observation point.
Theoretically, it can be solved for the linear group of Equations (1) using Equation (2).
The solution for A in Equation (2) only applies if A is a square matrix, and no longer applies if A is an overdetermined or underdetermined matrix. In geophysics, the amount of data collected is often much smaller than the parameters of the required solution, so the solution to Equation (1) is an underdetermined problem with multiple solutions. The Tikhonov regularization method, which has become a classic strategy, effectively solves this problem by minimizing the objective function, [45]. The objective function is composed of fitting items, model constraints, and regularization parameters, which is expressed as: Of these, φ d represents the data fitting item, φ m represents the model fitting item, and µ represents the regularization parameter, which is used to control the weight of the model fitting item and the data fitting item. When µ is large, the model fitting term has a greater impact on the equation, and when µ is small, the data fitting term has a greater impact on the equation. The data fitting term is expressed as follows: where 2 2 denotes the square of L2 norm, and W d is the data weighting matrix, usually composed of the covariance of the data noise. If the data noise is less, it can be taken as 1. The difference between the existing different forms of constraint inversion methods lies in the way of constructing model constraints, which requires the selection of appropriate constraints in relation to the geological problem. A general form of the model fitting term is as follows: where m 0 is the reference model and W m is the constraint weight term of the model, including the depth weight constraint W z [13] and the minimum support constraint W e [12]. Thus, Equation (5) can be written in the form of Equation (6).
W z is to counteract the influence of skin effect so that the inversion results are not concentrated at the surface. z is the central buried depth of the model element, and z 0 is the value related to the model element and height. β is a positive number used to control the size of the weighting function. W e is introduced to make the inversion parameters more focused. In Equation (8), m represents the value of any grid, and e is called the focusing factor, which is a constant with a small value. The smaller e is, the more focused the model parameters are. If m is very small, the value of w e tends to 0, and if m is very large, the value of w e tends to 1. When minimizing the objective function, it makes the value of w e smaller, resulting in more m tending to 0 and thus obtaining a sparse solution.
In this paper, the regularized conjugate gradient (RCG) method was used to solve the objective function [16]. Since the magnetism in the crust tends to be within a specific range, the addition of a physical upper and lower bound constraint can ensure that the values of the model parameters are within a reasonable range during iterative calculations [46], as shown in Equation (9), where m min and m max denote the upper and lower bounds of the magnetic parameter values, which can be obtained from geological data or logging data, respectively.
if m j < m min , then m j = m min if m j > m max , then m j = m max (9) It is worth noting that the amount of observational data is often very large, resulting in a large sensitivity matrix A. An oversized A places high demands on computer storage when solving the minimum of the objective function. In addition, the sensitivity matrix A increases the computational load of the solution process, resulting in a very slow efficiency [36]. In this paper, we solve this problem by compressing the amount of data.

Adaptive Quadtree Data Compression Algorithm
Adaptive quadtree decomposition is an image segmentation algorithm based on uniformity detection. The segmentation scale of the data area with a large gradient change and uneven value is very small, so more data can be obtained, whereas the segmentation scale of the data region with a small gradient change and uniform value is coarse, so less data is obtained. Compared with the regular sampling method with equal intervals down-sampling, this method retains more detail and achieves the best approximation to the original data with the least amount of data. Figure 1 shows the decomposition diagram of the three-level quadtree model. The quadtree decomposition algorithm is used to subdivide the two-dimensional space into four regions by recursion. The first level decomposition divides the region R into four sub-blocks R 1 , R 2 , R 3 and R 4 . In the second decomposition, the region consistency of each sub-block is judged in turn. If R 2 and R 3 meet the regional consistency standard, they are no longer decomposed while if R 1 and R 2 do not meet the regional consistency criteria, they continue to be decomposed into four sub-blocks, namely, R 1-1 , R 1-2 , R 1-3 , R 1-4 and R 4-1 , R 4-2 , R 4-3 , R 4-4 . meet the regional consistency standard, they are no longer decomposed while if R1 and R2 do not meet the regional consistency criteria, they continue to be decomposed into four sub-blocks, namely, R1-1, R1-2, R1-3, R1-4 and R4-1, R4-2, R4-3, R4-4. The uniformity of the data in the sub-blocks is called regional consistency. In this paper, the regional consistency standard is measured by variance, which is used to measure the magnitude and uniformity of fluctuations in the data. The larger the variance, the greater the difference in the data. Assuming that the amount of data in the segmented sub-block is n × n, the formula for calculating the variance is shown on the left side of Equation (9).
where ( , ) f i j represents the value of the point ( , ) i j , 1 i n ≤ ≤ ,1 j n ≤ ≤ , and f represents the average value of the data in the sub-block. T is the variance threshold that is used to judge whether the sub-block is uniform. The variance threshold is not fixed, and the selection of threshold affects the error between the compressed data and the original data and the data compression ratio. When the variance threshold is large, it means that more data will be discarded, and the data error will increase and the compression ratio will decrease. By drawing a profile of error and threshold and a profile of the compression ratio and threshold, the intersection of the two sections can be selected as the best threshold. Figure 2 shows a flowchart of the adaptive quadtree decomposition for data compression. First, given the segmentation size and variance threshold, when the sub-block size is larger than the minimum size and the variance Ti is greater than the given variance threshold T, the segmentation is continued until all the sub-blocks meet the conditions. In this way, the data from different sub-blocks of various sizes are obtained, and the average value and the coordinates of the center point are calculated, and are then stored as a compressed matrix. The method of using the average value of all data points of a sub-block to replace all of the points in the whole sub-block is called the mean method. It is essentially a low-pass filter that filters out the high-frequency part of the data for de-noising. The uniformity of the data in the sub-blocks is called regional consistency. In this paper, the regional consistency standard is measured by variance, which is used to measure the magnitude and uniformity of fluctuations in the data. The larger the variance, the greater the difference in the data. Assuming that the amount of data in the segmented sub-block is n × n, the formula for calculating the variance is shown on the left side of Equation (9).
where f (i, j) represents the value of the point (i, j), 1 ≤ i ≤ n,1 ≤ j ≤ n, and f represents the average value of the data in the sub-block. T is the variance threshold that is used to judge whether the sub-block is uniform. The variance threshold is not fixed, and the selection of threshold affects the error between the compressed data and the original data and the data compression ratio. When the variance threshold is large, it means that more data will be discarded, and the data error will increase and the compression ratio will decrease. By drawing a profile of error and threshold and a profile of the compression ratio and threshold, the intersection of the two sections can be selected as the best threshold. Figure 2 shows a flowchart of the adaptive quadtree decomposition for data compression. First, given the segmentation size and variance threshold, when the sub-block size is larger than the minimum size and the variance T i is greater than the given variance threshold T, the segmentation is continued until all the sub-blocks meet the conditions. In this way, the data from different sub-blocks of various sizes are obtained, and the average value and the coordinates of the center point are calculated, and are then stored as a compressed matrix. The method of using the average value of all data points of a sub-block to replace all of the points in the whole sub-block is called the mean method. It is essentially a low-pass filter that filters out the high-frequency part of the data for de-noising.

Synthetic Model Test
To illustrate the effectiveness and feasibility of data compression based on the adaptive quadtree decomposition algorithm in 3D magnetic inversion, two synthetic models are used in this paper.
Model 1: As shown in Figure 3a, the synthetic model is a simple cuboid with a length, width, and height of 400 m, 400 m, and 200 m, respectively, the buried depth of the top surface is 200 m, the magnetization is 5 A/m, the magnetic inclination is 60°, and the magnetic declination is 10°. Then, the forward calculation is carried out for the model. The observation height is set at 80 m, the survey network is 128 × 128, the spacing in x−axis and y−axis directions is 10 m, and the magnetic anomaly map with 7% Gaussian noise is added, as shown in Figure 3b.

Synthetic Model Test
To illustrate the effectiveness and feasibility of data compression based on the adaptive quadtree decomposition algorithm in 3D magnetic inversion, two synthetic models are used in this paper.
Model 1: As shown in Figure 3a, the synthetic model is a simple cuboid with a length, width, and height of 400 m, 400 m, and 200 m, respectively, the buried depth of the top surface is 200 m, the magnetization is 5 A/m, the magnetic inclination is 60 • , and the magnetic declination is 10 • . Then, the forward calculation is carried out for the model. The observation height is set at 80 m, the survey network is 128 × 128, the spacing in x-axis and y-axis directions is 10 m, and the magnetic anomaly map with 7% Gaussian noise is added, as shown in Figure 3b.

Synthetic Model Test
To illustrate the effectiveness and feasibility of data compression based on the adaptive quadtree decomposition algorithm in 3D magnetic inversion, two synthetic models are used in this paper.
Model 1: As shown in Figure 3a, the synthetic model is a simple cuboid with a length, width, and height of 400 m, 400 m, and 200 m, respectively, the buried depth of the top surface is 200 m, the magnetization is 5 A/m, the magnetic inclination is 60°, and the magnetic declination is 10°. Then, the forward calculation is carried out for the model. The observation height is set at 80 m, the survey network is 128 × 128, the spacing in x−axis and y−axis directions is 10 m, and the magnetic anomaly map with 7% Gaussian noise is added, as shown in Figure 3b.   The magnetic data for synthetic model 1 is compressed by the adaptive quadtree data compress algorithm. The parameter settings are as follows: the minimum grid is 4, the maximum grid is 16, and the variance threshold is 2000. Figure 4a,b show the cut image and the extracted data scatter map after data compression, respectively. The amount of data is reduced from 16384 points to 277 points. The compressed data is only 1.69% of the original data, and the error is 0.040.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 14 The magnetic data for synthetic model 1 is compressed by the adaptive quadtree data compress algorithm. The parameter settings are as follows: the minimum grid is 4, the maximum grid is 16, and the variance threshold is 2000. Figure 4a,b show the cut image and the extracted data scatter map after data compression, respectively. The amount of data is reduced from 16384 points to 277 points. The compressed data is only 1.69% of the original data, and the error is 0.040.  Figure 5b. The magnetic data for synthetic model 2 is compressed by adaptive quadtree decomposition algorithm. The parameter settings are as follows: the minimum grid is 4, the maximum grid is 16, and the variance threshold is 1600. Figure  6a shows the cut image using quadtree decomposition algorithm, and Figure 6b shows the compression algorithm extracts data scatter diagram, the data volume is reduced from 16384 points to 505 points, thus the compressed data volume is only 3.08% of the original data, and the error is 0.067.   Figure 5b. The magnetic data for synthetic model 2 is compressed by adaptive quadtree decomposition algorithm. The parameter settings are as follows: the minimum grid is 4, the maximum grid is 16, and the variance threshold is 1600. Figure 6a shows the cut image using quadtree decomposition algorithm, and Figure 6b shows the compression algorithm extracts data scatter diagram, the data volume is reduced from 16384 points to 505 points, thus the compressed data volume is only 3.08% of the original data, and the error is 0.067.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 14 The magnetic data for synthetic model 1 is compressed by the adaptive quadtree data compress algorithm. The parameter settings are as follows: the minimum grid is 4, the maximum grid is 16, and the variance threshold is 2000. Figure 4a,b show the cut image and the extracted data scatter map after data compression, respectively. The amount of data is reduced from 16384 points to 277 points. The compressed data is only 1.69% of the original data, and the error is 0.040.  Figure 5b. The magnetic data for synthetic model 2 is compressed by adaptive quadtree decomposition algorithm. The parameter settings are as follows: the minimum grid is 4, the maximum grid is 16, and the variance threshold is 1600. Figure  6a shows the cut image using quadtree decomposition algorithm, and Figure 6b shows the compression algorithm extracts data scatter diagram, the data volume is reduced from 16384 points to 505 points, thus the compressed data volume is only 3.08% of the original data, and the error is 0.067.  As can be seen from Figures 4b and 6b, the data points at the anomaly boundary are intensively sampled, which fully preserves the detailed features while only a small number of data points are reserved at the places where the anomaly changes slowly, thus greatly reducing the magnitude of the data. The compressed data are reconstructed as shown in Figures 4c and 6c. The results show that the compressed data retains the abnormal details and effectively eliminates the noise compared with the original data.
Taking the inversion calculation of model 2 as an example, we can compare the inversion using the compressed data and the inversion using the original data. The underground grid is divided into 64 × 64 × 32, and the grid spacing is 20 m × 20 m × 20 m. The horizontal and vertical slices of the inversion results are shown in Figure 7. When comparing the inversion results of the compressed data with the original data, it can be seen that the resolution of the compressed data is better than that of the original data. The results show that although some information is lost after data compression, the detailed features of the data are still retained. Since the step length and direction are calculated from all the data in the inversion iteration, when the inefficient information is reduced, the effective information is given more weight, which improves the efficiency of the iteration and makes the resolution of the anomaly boundary higher than the other locations, makes the inversion results more convergent, and can also show smaller anomalies. As can be seen from Figures 4b and 6b, the data points at the anomaly boundary are intensively sampled, which fully preserves the detailed features while only a small number of data points are reserved at the places where the anomaly changes slowly, thus greatly reducing the magnitude of the data. The compressed data are reconstructed as shown in Figures 4c and 6c. The results show that the compressed data retains the abnormal details and effectively eliminates the noise compared with the original data.
Taking the inversion calculation of model 2 as an example, we can compare the inversion using the compressed data and the inversion using the original data. The underground grid is divided into 64 × 64 × 32, and the grid spacing is 20 m × 20 m × 20 m. The horizontal and vertical slices of the inversion results are shown in Figure 7. When comparing the inversion results of the compressed data with the original data, it can be seen that the resolution of the compressed data is better than that of the original data. The results show that although some information is lost after data compression, the detailed features of the data are still retained. Since the step length and direction are calculated from all the data in the inversion iteration, when the inefficient information is reduced, the effective information is given more weight, which improves the efficiency of the iteration and makes the resolution of the anomaly boundary higher than the other locations, makes the inversion results more convergent, and can also show smaller anomalies. Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 14

Application in Mineral Exploration
The adaptive quadtree compression algorithm was applied to the magnetic data set collected for the Ma'anshan magnetite deposit. The magnetite deposit is located in the Ningwu mining area in the middle and lower reaches of the Yangtze River. The magmatic activity in this area was a strong, mainly volcanic eruption and magmatic intrusion in the Yanshanian period. The main subvolcanic rock body in the survey area is diorite porphyrite, which is the mineralization mother rock of iron and sulfur deposits in the mining area. Diorite porphyrite with its medium magnetism and magnetite with high magnetism have an obvious magnetic difference from the surrounding rocks, which are the basis of magnetic exploration. There are many sulfur and iron ore mining points in this area (see Figure 8). To find ore in the periphery and deep part of the old mine, the Unmanned Aerial Vehicle (UAV) aeromagnetic survey system was used to investigate this area. The total survey line length is 111.9 km, the average flight height is 84 m, the line distance is 100 m, and the average point distance is 0.86 m. A total of 145495 data points were collected. After data preprocessing and IGRF International Geomagnetic Reference Field (IGRF) correction, the minimum curvature gridding was carried out with a grid size of 25 m. Finally, the contour map of the total magnetic anomaly obtained is shown in Figure 8a, and the amplitude of anomaly is −1800 nT-3400 nT. The total magnetic field data is the superposition effect of different magnetic geological bodies. To separate the local field from the background field, the upward continuation method was used. Upward continuation was used to calculate the abnormal value of a certain height above the

Application in Mineral Exploration
The adaptive quadtree compression algorithm was applied to the magnetic data set collected for the Ma'anshan magnetite deposit. The magnetite deposit is located in the Ningwu mining area in the middle and lower reaches of the Yangtze River. The magmatic activity in this area was a strong, mainly volcanic eruption and magmatic intrusion in the Yanshanian period. The main subvolcanic rock body in the survey area is diorite porphyrite, which is the mineralization mother rock of iron and sulfur deposits in the mining area. Diorite porphyrite with its medium magnetism and magnetite with high magnetism have an obvious magnetic difference from the surrounding rocks, which are the basis of magnetic exploration. There are many sulfur and iron ore mining points in this area (see Figure 8). To find ore in the periphery and deep part of the old mine, the Unmanned Aerial Vehicle (UAV) aeromagnetic survey system was used to investigate this area. The total survey line length is 111.9 km, the average flight height is 84 m, the line distance is 100 m, and the average point distance is 0.86 m. A total of 145495 data points were collected. After data preprocessing and IGRF International Geomagnetic Reference Field (IGRF) correction, the minimum curvature gridding was carried out with a grid size of 25 m. Finally, the contour map of the total magnetic anomaly obtained is shown in Figure 8a, and the amplitude of anomaly is −1800 nT-3400 nT. The total magnetic field data is the superposition effect of different magnetic geological bodies. To separate the local field from the background field, the upward continuation method was used. Upward continuation was used to calculate the abnormal value of a certain height above the observation plane according to the measured anomaly [46,47]. Since the potential field value is inversely proportional to the square of the distance from the field source to the observation point, the field value of the field source with different depths has different changes in speed when the same height is extended. When the upward continuation is performed, the field source with large depth is highlighted because of its slow attenuation speed, thus the deep anomaly can be highlighted. In the upward continuation of the data, when the continuation height is greater than 300 m, the field change is very small, which can be approximately regarded as the regional magnetic field. As shown in Figure 8b, the residual gravity anomaly can be separated by subtracting the upward continuation anomaly from the total magnetic anomaly. There are two main North-East (N-E) trending high magnetic anomaly belts, in which there are many elliptical and equiaxed local anomalies with large amplitude variations. Some of the local anomalies correspond to the known iron deposits in the area.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 14 observation plane according to the measured anomaly [46,47]. Since the potential field value is inversely proportional to the square of the distance from the field source to the observation point, the field value of the field source with different depths has different changes in speed when the same height is extended. When the upward continuation is performed, the field source with large depth is highlighted because of its slow attenuation speed, thus the deep anomaly can be highlighted. In the upward continuation of the data, when the continuation height is greater than 300 m, the field change is very small, which can be approximately regarded as the regional magnetic field. As shown in Figure 8b, the residual gravity anomaly can be separated by subtracting the upward continuation anomaly from the total magnetic anomaly. There are two main North-East (N-E) trending high magnetic anomaly belts, in which there are many elliptical and equiaxed local anomalies with large amplitude variations. Some of the local anomalies correspond to the known iron deposits in the area. To restore the distribution of magnetic field sources, the underground was divided into 64 × 64 × 20 grids, each with a length of 50 m. Firstly, the adaptive quadtree decomposition method was used to compress the magnetic data, and the variance threshold was set to 7800 (see Figure 9). The location map of the sampling points is shown in Figure 10. The black cross in the graph represents the position of the sampling data, with 1534 sampling points, accounting for 1.05% of the original data. The compressed data were used for 3D magnetic inversion, and the inversion results are shown in Figure 11.  To restore the distribution of magnetic field sources, the underground was divided into 64 × 64 × 20 grids, each with a length of 50 m. Firstly, the adaptive quadtree decomposition method was used to compress the magnetic data, and the variance threshold was set to 7800 (see Figure 9). The location map of the sampling points is shown in Figure 10. The black cross in the graph represents the position of the sampling data, with 1534 sampling points, accounting for 1.05% of the original data. The compressed data were used for 3D magnetic inversion, and the inversion results are shown in Figure 11.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 14 observation plane according to the measured anomaly [46,47]. Since the potential field value is inversely proportional to the square of the distance from the field source to the observation point, the field value of the field source with different depths has different changes in speed when the same height is extended. When the upward continuation is performed, the field source with large depth is highlighted because of its slow attenuation speed, thus the deep anomaly can be highlighted. In the upward continuation of the data, when the continuation height is greater than 300 m, the field change is very small, which can be approximately regarded as the regional magnetic field. As shown in Figure 8b, the residual gravity anomaly can be separated by subtracting the upward continuation anomaly from the total magnetic anomaly. There are two main North-East (N-E) trending high magnetic anomaly belts, in which there are many elliptical and equiaxed local anomalies with large amplitude variations. Some of the local anomalies correspond to the known iron deposits in the area. To restore the distribution of magnetic field sources, the underground was divided into 64 × 64 × 20 grids, each with a length of 50 m. Firstly, the adaptive quadtree decomposition method was used to compress the magnetic data, and the variance threshold was set to 7800 (see Figure 9). The location map of the sampling points is shown in Figure 10. The black cross in the graph represents the position of the sampling data, with 1534 sampling points, accounting for 1.05% of the original data. The compressed data were used for 3D magnetic inversion, and the inversion results are shown in Figure 11.  A total of 10 medium and high magnetization anomalies were obtained by 3D inversion, and were labeled M1-M10. The spatial location map is shown in Figure 11a. Figure 11b shows the inversion results at different depths. The distribution of anomalies corresponds to the surface magnetic anomalies, and their buried depth is between 100 m and 500 m. Figure 11c-f show the location characteristics of ore bodies at different sections. Combining existing geological data and inversion results, the M1 anomaly was located in the Xiang Shan (XS) pyrite area, which is a mined-out area above 100 m underground. The inversion results show that there are magnetite ore bodies under the mined-out area, inclined to the southeast with a deepest depth of 500 m. The M2 anomalous body is located in the south of the Nan Shan (NS) iron mining pit. A large number of diorite porphyrite, magnetite orebody, and pyrite ore body are exposed on the surface. It is speculated that the magnetite (mineralized) body exists underground at a buried depth of 200-400 m. Diorite porphyry and magnetite ore bodies are exposed on the surface where the M3 anomaly is located. Ao Shan (AS) iron ore is located near the southeast side. It is speculated that the anomaly is a magnetized ore body. M4 anomaly body is located on the north side of the AS iron ore pit. The exposed stratum is the Lower Cretaceous, and the residual magnetization is large. It is presumed to be a magnetized mineralized body. The buried depth is about 100 m, the thickness is about 400 m, and the width is more than 300 m. The M5 anomalous body is located in the Luo Bo Shan (LBS) mining area. The magnetite at the top of the location has been mined. There is magnetite in the deep rock body at a depth of more than 300 m, and the thickness is less than 100 m. The M8 anomaly shows obvious high magnetization, and it is speculated that there are magnetite ore bodies in the lower part. The anomalous bodies (M6, M7, M9, M10) in the northeast and southeast of the survey area have low magnetization and are presumed to be diorite porphyrite. Therefore, there are five magnetite bodies, M1, M2, M3, M4, and M8 in this area, which are large in scale and they have potential for development.

Conclusions
In this paper, an adaptive quadtree decomposition method is proposed to compress magnetic data, which is decomposed recursively until the data in each quadrant meets the regional consistency standard. The model tests show that the method can be used for fine sampling where the gradient of magnetic anomaly data is large, and for sparse sampling at the place with a gentle anomaly change, to realize the effective compression of magnetic data and greatly improve the calculation efficiency of 3D inversion. In addition, the mean method can suppress noise and improve the resolution of inversion results.
The method was applied to the magnetic data for the Ma'anshan mining area, and five possible areas of magnetite orebodies were delineated in the surrounding and deep part of the existing mining area, which provides some significant guidance for the next prospecting work.
The 3D magnetic inversion algorithm is an important means of data interpretation. However, the calculation efficiency is very low. The adaptive quadtree decomposition calculation is only used to compress the data space. When the amount of data is large, the underground subdivision grid is very large, which still leads to the calculation speed being very slow. The calculation efficiency can be further improved by combining wavelet compression technology, adaptive mesh generation A total of 10 medium and high magnetization anomalies were obtained by 3D inversion, and were labeled M1-M10. The spatial location map is shown in Figure 11a. Figure 11b shows the inversion results at different depths. The distribution of anomalies corresponds to the surface magnetic anomalies, and their buried depth is between 100 m and 500 m. Figure 11c-f show the location characteristics of ore bodies at different sections. Combining existing geological data and inversion results, the M1 anomaly was located in the Xiang Shan (XS) pyrite area, which is a mined-out area above 100 m underground. The inversion results show that there are magnetite ore bodies under the mined-out area, inclined to the southeast with a deepest depth of 500 m. The M2 anomalous body is located in the south of the Nan Shan (NS) iron mining pit. A large number of diorite porphyrite, magnetite orebody, and pyrite ore body are exposed on the surface. It is speculated that the magnetite (mineralized) body exists underground at a buried depth of 200-400 m. Diorite porphyry and magnetite ore bodies are exposed on the surface where the M3 anomaly is located. Ao Shan (AS) iron ore is located near the southeast side. It is speculated that the anomaly is a magnetized ore body. M4 anomaly body is located on the north side of the AS iron ore pit. The exposed stratum is the Lower Cretaceous, and the residual magnetization is large. It is presumed to be a magnetized mineralized body. The buried depth is about 100 m, the thickness is about 400 m, and the width is more than 300 m. The M5 anomalous body is located in the Luo Bo Shan (LBS) mining area. The magnetite at the top of the location has been mined. There is magnetite in the deep rock body at a depth of more than 300 m, and the thickness is less than 100 m. The M8 anomaly shows obvious high magnetization, and it is speculated that there are magnetite ore bodies in the lower part. The anomalous bodies (M6, M7, M9, M10) in the northeast and southeast of the survey area have low magnetization and are presumed to be diorite porphyrite. Therefore, there are five magnetite bodies, M1, M2, M3, M4, and M8 in this area, which are large in scale and they have potential for development.

Conclusions
In this paper, an adaptive quadtree decomposition method is proposed to compress magnetic data, which is decomposed recursively until the data in each quadrant meets the regional consistency standard. The model tests show that the method can be used for fine sampling where the gradient of magnetic anomaly data is large, and for sparse sampling at the place with a gentle anomaly change, to realize the effective compression of magnetic data and greatly improve the calculation efficiency of 3D inversion. In addition, the mean method can suppress noise and improve the resolution of inversion results.
The method was applied to the magnetic data for the Ma'anshan mining area, and five possible areas of magnetite orebodies were delineated in the surrounding and deep part of the existing mining area, which provides some significant guidance for the next prospecting work.
The 3D magnetic inversion algorithm is an important means of data interpretation. However, the calculation efficiency is very low. The adaptive quadtree decomposition calculation is only used to compress the data space. When the amount of data is large, the underground subdivision grid is very large, which still leads to the calculation speed being very slow. The calculation efficiency can be further improved by combining wavelet compression technology, adaptive mesh generation technology, parallel computing, and so on.