Deep Mineral Exploration of the Jinchuan Cu–Ni Sulfide Deposit Based on Aeromagnetic, Gravity, and CSAMT Methods

: The exploration of deep mineral resources is an important prerequisite for meeting the continuous demand of resources. The geophysical method is one of the most effective means of exploring the deep mineral resources with a large depth and a high resolution. Based on the study of the geological background, petrophysical properties, and aeromagnetic anomaly characteristics of the Jinchuan Cu–Ni sulfide deposit, which is famous throughout the world, this paper uses the widely used gravity, aeromagnetic, and CSAMT (controlled source audio ‐ frequency magnetotellurics) methods with a complementary resolution to reveal the favorable prospecting position. In order to obtain better inversion results, the SL0 norm tight support focusing regularization inversion method is introduced to process the section gravity and aeromagnetic data of the mining area. By combining the results with CSAMT, it is found that the medium ‐ low resistivity, high density, and the high magnetic anomaly areas near the structural belt can nicely correspond with the known ore ‐ bearing rock masses in the mining area. At the same time, according to the geophysical exploration model and geological and physical property data, four favorable ore ‐ forming prospect areas are delineated in the deep part of the known mining area. the rock and ore bodies. The intersecting parts of faults in different directions form irregular columnar ore bodies.


Introduction
Mineral resources are the material basis of human science and technology progress and social and economic development [1]. In recent years, with the continuous and stable growth of China's economy, the contradiction between the supply and demand of mineral resources has become increasingly prominent, among which all kinds of metal mineral resources are in short supply. The effective way to alleviate the shortage of resource supply is to carry out deep prospecting and strengthen the resource reserve. The Jinchuan Cu-Ni sulfide deposit is one of the three largest Ni-Cu-(PGE) deposits in the world. It is of great significance to carry out deep prospecting in this area to ensure the supply of copper and nickel resources. Moreover, the available geological data show that the deep parts of the known mining area and the surrounding area have a good prospecting potential, mainly based on the following: according to the metallogenic model, structural characteristics, and the spatial location of ore-bearing magma emplacement, there may be more ore bodies in the deep part of the mining area; the lower end of the main rock mass in the first, second, and third mining areas is not completely revealed; the new ore body is indeed found in the deep of the II mining area [2]. An independent ore body has been found in the surrounding rock at the bottom of No.24 ore body, with an increase of more than 600,000 tons of nickel metal. It is believed that there may be a relatively rich copper-nickel sulfide ore bodies in the surrounding rocks in the footwall direction of the western rock bodies of Jinchuan, and there may be a large number of sulfide residues in the deep magma chamber corresponding to the eastern rock body [3]. For the surrounding area of the mining area, the previous data show that the rock mass of the mining area shows the characteristics of echelon arrangement in space, and there are new ore-bearing rock masses near the main ore body, such as No. 58 ore body of the third mining area [4].
The deep prospecting methods mainly include geological, geophysical, geochemical, and drilling methods. The geological prospecting method relies on the observation and analysis of surface outcrop or drilling core to infer the underground geological conditions, which has limitations for deep prospecting, especially for the concealed characteristics. Geochemical exploration is an important technical support condition for deep prospecting, and drilling engineering is the realization condition. However, the depth and scope that these two methods can reach is also limited. Geophysical exploration is the basic means to obtain the information of concealed parts beyond other prospecting methods, which has a large detection depth, a high resolution, and various means [5][6][7][8]. It can carry out multi-scale detection in the target area and provide rich information for deep ore prospecting [9,10]. Kheyrollahi et al. [11] discovered and predicted the distribution pattern of porphyry copper deposits in the tertiary magmatic belt by the upward extension and boundary enhancement of magnetic anomalies. Xiao and Wang [12] used Bouguer gravity and aeromagnetic data to further understand the geological and mineral resources near the porphyry copper molybdenum polymetallic mineralization in the Tianshan area of China. Hu et al. [13] explored potential iron and polymetallic lead-zinc-copper deposits in the Longmen area by the CSAMT method and found high-grade lead-zinc-silver-titanium ore through drilling based on inversion results. Guo et al. [14] applied the CSAMT method to the exploration of the Jianshui lead-zinc mine and drew the underground resistivity distribution map through data processing and inversion. According to the CSAMT results, the location of the ore body is inferred, and the results are verified by drilling. The lead-zinc ore body is 373.70-407.35 m in the well. Shah et al. [15] comprehensively used aeromagnetic induced polarization, magnetotelluric and borehole geological alteration, magnetic susceptibility, and density data to explore the copper-gold molybdenum Pebble porphyry deposit, and achieved good prospecting results. In order to accurately detect the underground structure of complex deposits and solve the problems of uniqueness and inconsistency in the single parameter inversion model, Zhang and Li [16] proposed a two-dimensional gravity gradient and a magnetotelluric joint inversion method based on data space and normalized cross gradient constraints. Melo et al. [17] proposed a geological characterization method that can identify copper deposits based on geophysical inversion. This method can use geophysical data and sparse geological information to evaluate the target quickly, especially for the first stage of deep target or concealed target exploration. The success of this method is verified by the inversion of magnetic data and the direct current resistivity data of Cristalino iron oxide copper gold deposit in northern Brazil. Lee et al. [18] obtained the resistivity model consistent with the regional geology through the 3D joint inversion of magnetotelluric and Z-axis tipper electromagnetic data, which not only shows the mineralization belt interpreted for the Morrison porphyry Cu-Au-Mo deposit but is also conducive to the exploration of the disseminated sulfide of other porphyry deposits.
Some researchers also used geophysical methods to carry out prospecting work in the Jinchuan copper-nickel mining area and its surrounding areas and obtained some knowledge or achievements. Through the joint interpretation of gravity and magnetic data, it is considered that the M-15 anomaly is caused by ultrabasic rocks with a buried depth of more than 1200 m, which has a positive effect on the indication of deep Cu-Ni deposits [19]. According to the comprehensive prospecting model of geology, geophysics, and geochemistry, Wen and Luo [20] carried out prospecting and prediction work in the deep and edges of the Jinchuan copper-nickel mining area and found five potential target areas. In 2006, Fu and Li [21] established a comprehensive geological geophysical prospecting model based on the characteristics of the geophysical geochemical field of rock masses and ore deposits in different mining areas. On the basis of a systematic analysis of metallogenic geological conditions and comprehensive geophysical, geochemical, and remote sensing information, Gao [2] believed that the joint area of the first and second mining areas, the joint area of the I and II mining areas, the joint area of No.1 and No.2 ore bodies in the II mining area, and the overlap areas of geophysical and geochemical anomalies in the III mining area are important locations for ore body tracing.
On the basis of systematically summarizing the geological background of mineralization, this paper uses the aeromagnetic and gravity methods with low exploration cost and high efficiency, and the CSAMT method with a large exploration depth and a high vertical resolution to indicate the favorable metallogenic locations in the deep of the Jinchuan Cu-Ni deposit and its surrounding area. In order to overcome the weakness of deep signal in deep gravity and magnetic exploration, the processing method of potential field data is studied. Additionally, the focus inversion method based on the SL0 norm with a good convergence effect is used to process the aeromagnetic or gravity data of the profile, and, based on these results and the CSAMT inversion results, geological and physical properties data and geophysical profile data are interpreted. The results show that the range of the abnormal bodies corresponds well with the ore-bearing bodies of the known mining areas, and four favorable metallogenic targets are delineated according to the results.

Geological Background of the Survey Area
The Jinchuan Cu-Ni sulfide deposit is a part of the Longshoushan metallogenic belt, which is located in Longshoushan terrane in the southwest of the Alxa block of the North China Craton (NCC) [22,23]. The NCC is one of the three major Precambrian blocks in China, which formed from an amalgamation of micro blocks [24][25][26]. The Alxa block is located in the westernmost part of the NCC and is in fault contact with the Tarim craton to the west, bounded by the North Qilian orogenic belt in the south and the Central Asia orogenic belt in the north (Figure 1a) [27,28].The Longshoushan terrane is a long narrow northwest trending uplift. It is 195 km long and 30-35 km wide, which is controlled by deep faults on both sides of the north and the south (F1 and F2). Its north side is adjacent to the Chaoshui basin, and its south side is separated from the Qilian orogenic belt (Figure 1b). In the Longshoushan terrance, the main outcropping strata are the Paleoproterozoic Longshoushan group, the late Mesoproterozoic Dunzigou group [29], and the Neoproterozoic-Cambrian Hanmushan group [30,31]. The Longshoushan group is the oldest metamorphic basement in the Longshoushan terrane [32]. Strong metamorphism and deformation in multiple periods [33] make the strata fragmented and it is difficult to judge the original stratigraphic sequence [34]. It mainly consists of schlieren and homogenic migmatites, marble, biotite-plagioclase gneiss, granulites, quartz schist, amphibolite, and pyroclastic rock [27,35]. The Dunzigou group is the earliest sedimentary overlying strata in the Longshoushan terrane, which is in angular unconformity contact with the lower Longshoushan group. In the geological history, the Longshoushan terrane has experienced multiple structural changes. The present NW trending faults and folds are mainly Caledonian and later structures. Faults in the Longshoushan terrane are mainly NW-and NE-trending structures, cutting the metamorphic formation [23]. Magmatic rocks are widely developed in the region, and magmatism occurred from Paleoproterozoic to Neoproterozoic, mainly in the Paleoproterozoic and the Paleozoic [36][37][38][39]. With regard to the formation age of Jinchuan intrusion, different researchers have adopted different methods to obtain a large number of isotopic age data, which can be roughly divided into two ranges: 1400-1600 Ma and 800-1000 Ma, representing Mesoproterozoic and Neoproterozoic, respectively [40,41].
The ore bearing ultrabasic rock bodies unconformity intrudes into the Baijiazuizi formation of pre-Great Wall system, which is in direct contact with gneiss, marble and banded migmatite, in the form of wall. It is about 6500 m long, 20-527 m wide on the surface (Figure 2a) and has a southwestern downward extension of more than 1000 m from the ground surface ( Figure 2b) [28]. The overall strike is N50° W, inclined to the SW, with a dip angle of 50°-80°.The rock mass area is about 1.34 km 2 . The ultramafic intrusion is divided into four sections by F8, F16-1, and F23andisnumbered as III, I, II, and IV from west to east (Figure 2a), which are corresponding to the four mining areas, respectively. The three main ore bodies with proved reserves of great economic value are respectively hosted in No.1 and No.2 ore bodies of Segment II and No.24 ore body of Segment I, and No.58 ore body is hosted in an independent rock body in the southwest side of Segment III (Figure 2b). The main strata in the mining area are the Baijiazuizi formation and the Quaternary system of the Longshoushan group of the pre-Great Wall system. The Quaternary sediments are mainly distributed in the east and west ends of the Jinchuan intrusion and the north of F1. Baijiazuizi formation is the direct surrounding rock of the Jinchuan deposit, which is distributed in a NW-SE direction, consistent with the regional structural direction. Baijiazuizi formation underwent multiple magmatic intrusions and multiple metamorphisms, forming a series of rocks mainly composed of migmatite, gneiss, and marble ( Figure 2a) [42].Among them, gneiss with stable chemical properties and poor water permeability is a good barrier layer, which enables the ore-forming materials to fully crystallize and differentiate in ultrabasic magma. Marble is active in chemical properties. It is favorable for the formation of contact metasomatic mineralization [43]. Therefore, Baijiazuizi formation is an important prospecting indicator.
The mining area has experienced structural activities many times. The structures of different periods and different directions superposed each other, making the mining area fold and fracture developed. The axial near the EW fold group includes the anticline where the deposit is located and a large syncline in the south of the mining area. The axial near the NE fold group is nearly vertical to the NW direction main structural lines of the mining area, among which the NE direction fold group across the ultramafic rock mass is the most significant, which plays an important role in the shape change and mineralization re-enrichment of the ore body, and the rich ore body is obviously thickened at the turning part of the fold. The NW trending faults are the most developed, followed by the NE and nearly EW trending faults. As one of the most important ore-controlling factors, faults not only control the emplacement of an ore bearing rock but also control the re-enrichment of mineralization and the spatial position of ore body. The NW ore-controlling faults are related to the spatial distribution of ore bodies. The NE or near EW ore controlling faults mainly cut the rock and ore bodies. The intersecting parts of faults in different directions can form irregular columnar ore bodies.
The ore-forming materials in the mining area mainly come from ultramafic magma. Ultrabasic rock is the ore-forming parent rock and the surrounding rock of the main ore body. The relationship between ultrabasic rock and mineralization is mainly reflected in the spatial change of rock mass and the relationship between lithofacies and mineralization. Only the spatial relationship between the two is introduced here. The shape of ore-bearing rock mass is irregular, and the development of the ore body is closely related to the floor. Generally, the concave part of the floor is favorable for the accumulation of ore-forming materials, and the ore body is thick. The occurrence of ore-bearing rock mass controls the occurrence of the stratoid ore body, the ore body is the same as the rock mass, and the strike is NW. The thickness of the ore-bearing rock mass is related to the thickness of the rich ore body. At the bottom of the thick rock mass, the rich ore body is also thick, which can be seen in the first and second mining areas.

Aeromagnetic, Gravity, and CSAMT Surveys
Based on the characteristics of density, magnetism, and resistivity (Table 1), the rocks and ores in the mining area can be roughly divided into three categories. The first category is copper-nickel ore, showing the characteristics of high density, high magnetism, and low resistivity; the second category is ultrabasic rock, showing the characteristics of high density, strong magnetism, and medium resistivity; the third category is the rock surrounding the ultrabasic rock, with the characteristics of low density, weak magnetism, and high resistivity. The differences of these physical properties provide a precondition for geophysical exploration work such as gravity, magnetic, and electrical methods in the study area [21]. 16  Convenient and efficient aeromagnetic exploration has been carried out in the main mining area. The aeromagnetic work uses the power glider as the carrier, and the measuring instrument is the helium optical-pumping magnetometer with a sensitivity of 0.001 nT. The average flight height is 93 m, the measurement scale is 1:10,000, the distance between survey lines is 100 m, and the distance between points is 2.7-3 m. The maximum dynamic noise level of the survey line in the survey area is 0.052 nT, most of which is less than 0.04 nT, and the average value is 0.027 nT, meeting the measurement requirements. Reduction to the pole can eliminate the asymmetry of the magnetic anomaly position caused by the declination and inclination of the magnetization field. After reducing magnetic anomaly to the pole, the anomaly information is more abundant, including the anomalies of different properties, scales, and depths. It is the basic data for anomaly interpretation. The induced magnetization of ultrabasic rocks and ores with high magnetic susceptibility in the Jinchuan Cu-Ni mining area is obviously greater than the residual magnetization [19], so reduction to the pole can be carried out, and the result is shown in Figure 3. The negative aeromagnetic anomaly of the mining area is located in the northeast, the isoline is relatively disordered, and the minimum negative anomaly is less than-140 nT. The positive aeromagnetic anomaly is mainly located in the southwest and central part, showing a significant northwest distribution. The known mining areas III, I, II, and IV (magenta curve range in Figure 3) are all in the high positive aeromagnetic anomaly area. The high positive anomaly in the III mining area is nearly circular, with a diameter of about 750 m, an area of about 0.4 km 2 , and a maximum anomaly intensity of more than 350 nT. The high positive anomaly in the I mining area extends northwestward in a belt, with a length of about 1400 m, a width of about 600 m, and an area of about 0.75 km 2 . The maximum anomaly value is located in the southeast end of the mining area, and the maximum anomaly value is more than 600 nT. The two ends of the high normal abnormal morphology in the II mining area have obvious distortion, but generally it is a strip extending northwestward, with a length of about 3000 m, a width of about 900 m in the west section, a width of about 650 m in the east section, an area of about 2.23 km 2 , and the maximum abnormal value is more than 1500 nT. In addition, the cascade zones on both sides of the I and II mining areas are relatively steep, showing the characteristics of steepness in the northeast and slowness in the southwest, suggesting that the abnormal body is steeply inclined to the southwest. The corresponding high positive anomaly of the IV mining area has a irregular ellipse shape, with long axis in east-west direction, about 1000 m long, 850 m wide, and an area of about 0.75 km 2 . Combined with the geological map and the characteristics of physical parameters of rocks and ores, the high positive aeromagnetic anomaly in the mining area is mainly caused by the ore bearing ultrabasic rocks, so the high positive aeromagnetic anomaly is an important indicator of ultrabasic rocks. In order to understand the scope of deep ore-bearing rock masses in the mining area, we use the boundary enhancement method presented by Zhang et al. [44] to detect the boundaries of aeromagnetic anomaly in the mining area. It can be seen from the results that the scope of the boundaries of the underground rock masses is determined (Figure 4), especially the edge positions of the deep rock masses of the four mining areas completely covered by the Quaternary are delineated, which provides the exploration scope for the prospecting of the deep Cu-Ni deposit. In order to better understand the characteristics of the deep rock anomalies, we first introduce the dualtree complex wavelet into the multi-scale anomaly separation of aeromagnetic anomalies. The dualtree complex wavelet not only has the advantages of wavelet transformation but also the characteristics of approximate translation invariance, more directional selectivity, and limited data redundancy. The results show that with the increase of decomposition scale, the detail information in the shallow part decreases gradually ( Figure 5). The range of the high normal anomaly in the second mining area in the southwest side is gradually expanding, indicating that with the increase of burial depth, the range of ore bearing ultrabasic rock masses is gradually expanding, which shows that the deep mining area has a good prospecting potential.  On this basis, we use the gravity and CSAMT data to acquire the locations of the deep potential ore-bearing rock masses in the II and IV mining areas, and the locations of the survey lines are shown in Figure 6.The profiles of gravity and CSAMT data are designed according to the characteristics of magnetic anomalies and the structure of the mining area. The direction of profile Lmg1~3 is 38.38° and that of Lmg4 is 0°. The parameters of CSAMT measurement of the four profiles are determined according to the proposed exploration depth. The minimum receiving and transmitting distance is 12 km, and the maximum is about 16 km. The power supply electrodes are arranged parallel to the survey line, the electrodes distance is 2 km, and the azimuth error is less than 3°. The sampling frequency is 1-9600 Hz, and there are 41 sampling frequency points. The station distance is 50 m. The quality inspection of the CSAMT measurement adopts the method of data observation on the inspection point again. The data quality evaluation is to calculate the mean square relative error of the resistivity of the inspection point. The mean square relative error of the resistivity of the single point in this work is 1.3-4.9%, which meets the specification and design requirements of less than or equal to 5%.The scale of gravity profile (Lmg-1 and 4) work is 1:5000, and the distance between profile points is 20 m. The instrument used is a CG-5 high-precision gravimeter with a reading resolution of 1 μgal. Before field operation, in addition to various checks and adjustments, static and dynamic tests are carried out to ensure the good performance of the gravimeter in use. The total mean square error of the Bouguer gravity anomaly is ±0.079 × 10 5 m/s 2 , which meets the design requirements and has reliable quality.

Inversion Methods
In this part, the section gravity, corresponding aeromagnetic data and CSAMT data are inverted to understand the distribution of different geological bodies in the depth of the second and fourth mining areas, so as to predict the favorable prospecting prospect combining with geological data. In this paper, the constraint inversion method based on the SL0 norm tight support is used to process the section gravity and aeromagnetic data.
The following is the principle of tightly supported focused inversion based on the SL0 norm. If m is model space and d is data space, the relationship between the two is F, and the forward process is expressed as: The inversion is expressed as: Potential field data inversion is an underdetermined problem. In order to reduce multiple solutions, the Tikhonov regularization method is commonly used. The inversion process can be written as: Among them,   m  is the two norm of the difference between the observed data and the theoretical forward data,  is the regularization parameter,   m s is the stabilizer, which represents the model objective function based on the prior information constraint. In this paper, the minimum compactly supported functional is used as the stabilizer, which can make the inversion result have a better focusing effect [45].
The integral equation is used to express the minimum support functional stabilizers as follows: Among them,｛a, b｝represents the internal product of a and b, which is the focusing factor, and e is the focus factor, which is related to the focus effect of the inversion output. m apr is a model based on prior information.
To simplify the above formula, variable weight functional ) (m e  is introduced, which is expressed as follows: Equation (4) The objective function given by Equation (7) is similar to the traditional objective function form. The difference is that a variable weight matrix needs to be introduced into the model parameters of Equation (7). This paper uses the conjugate gradient method to solve the problem of parameter functional minimization given by Equation (7). In Equation (5), can be regarded as a regularization parameter , and since it also has a focusing effect, it can be called a regularizationfocusing factor. As this parameter becomes smaller, the corresponding stabilizer can minimize the non-zero deviation of the model parameters from the prior information.
The smooth L0 algorithm (SL0 algorithm) comes from the sparse signal recovery theory and is used to solve the problem of how to accurately solve m in the inverse problem. It uses a suitable smooth continuous function to approximate the discontinuous L0 norm and minimizes it by using a minimization algorithm on the smooth function, thereby obtaining the minimum L0 norm and obtaining a sparse solution. In this paper, a Gaussian function with an expected value of 0 is selected to approximate the smooth function of the L0 norm. The continuous function be expressed as: Among them,  represents the approximate degree of continuous and discontinuous L0 norm. Then there are: (11) Or approximately: Define a new function: Then: The above formula shows that is true when  is small, and,when 0   , this approximate relationship tends to be equal. Therefore, in order to find the solution with the smallest L0 norm, we can take a small value of , and make ( ) F m  the maximum. For small values of  , F  is highly uneven and contains many local maxima, so it is difficult to maximize it. For a large value of  , F  is smooth and contains fewer local maxima, so it is easier to maximize it. In order to have the largest F  for any value of  , this paper uses a decreasing sequence of  to maximize F  . For each  with a large front value, the initial value of the maximization algorithm of F  is the maximum value of the corresponding F  . When  gradually decreases, the initial value of F  corresponding to each  starts from the maximum value close to the actual F  . Therefore, the SL0 algorithm does not fall into the local maximum problem and can find the actual maximum value of F  for a small value of  , and give the solution of the smallest L0 norm. Compared with a tightlysupported focused inversion, SL0 norm-constrained tightly-supported focused inversion continuously adjusts W e based on a priori information in the form of a weighting function, making the inversion results more accessible to actual physical parameter models. Therefore, the objective function of the inversion method based on SL0 norm tight support focus can be expressed as follows: For the CSAMT data, this paper uses the widely used conventional SCS2D software for inversion. This program is an active audio magnetotelluric data processing program developed on the basis of magnetotelluric data processing. Its development level is relatively mature. The suitability selection of SCS2D software inversion parameters is an important part of data processing. The correct selection of inversion parameters will directly affect the accuracy of subsequent data interpretation. The initial background model selected in this paper is a moving average model.

Inversion Results and Interpretation
In this section, the section aeromagnetic and gravity data of the II and IV mining areas are inversed based on the SL0 tight support focus inversion method. At the same time, combined with the CSAMT inversion results and geological and rocks' physical properties, the corresponding structures of the survey lines are inferred and interpreted, and the favorable positions of deep mineralization are delineated.

Survey Results of the II Mining Area
The II mining area is located in the southeast of the F16-1 fault and the northwest of the No.56 exploration line, with the largest copper-nickel ore body developed in the Jinchuan Cu-Ni deposit. The ore bearing strata are mainly pre-Sinian Baijiazuizi formations. Faults are developed, mainly including three groups of faults in the NW, NE, and nearly EW directions. The ultrabasic rock body is the ore-bearing parent rock. Under the control of the project, the rock body is in the shape of a rock wall, trending to the northwest, inclining to the southwest, with a length of more than 3000 m. the horizontal thickness is shown as thin at both ends and thick in the middle, up to 1550 m. West of line 26 of the mining area, the occurrence of the ore body is relatively steep, plate like and lens like, and is in the form of a completely intrusive contact or mixed gradual intrusive contact with the early lithofacies. The ore body to the east of line 26 is in lenticular or stratoid shape. For the deep prospecting of the mining area, the previous study shows that the deep part of No.2 ore body in the II mining area has the possibility of a branch compound and a pinch-out reappearance of ore body. The main reason for this is that the extension of No.2 ore body below the 1000 m level is not revealed, and the geological sections of line [28][29][30] show that the ore body below the 1100 m level has not been pinched out [2]. Based on the multiple geophysical data, this paper investigates the deep of No.1 ore body so as to find out whether there is the possibility of new ore body in the deep.
The three survey lines of Lmg-1-3 arranged in No.1 ore body of the II mining area coincide with the No.8, 12, and 14 exploration lines, respectively, which are close to each other and are arranged in parallel along the southeast direction. It is of great significance for indicating the deep resistivity, density, and magnetic variations in the profiles. Gravity and CSAMT explorations have been carried out along the Lmg-1 line, respectively. Figure 7 shows the inversion results of gravity and corresponding aeromagnetic data by using the SL0 method. The density and magnetism of the media under the line are obviously different. There are two obvious high-density abnormal areas in the profile, which extend to the deep of the southwest part. The high-magnetism abnormal areas correspond to the high-density abnormal areas and show similar changing characteristics towards the deep. The inversion result of CSAMT shows that the high resistivity areas are distributed on both sides of the profile, the medium-low resistivity areas are mainly located in the middle, tend to the southwest, the dip angle changes from steep to slow, and the extension is large, which is consistent with the high density and high magnetic areas in Figure7 (Figure8). CSAMT exploration was also carried out along the Lmg-2 and Lmg-3 lines. The resistivity inversion results and the corresponding aeromagnetic anomaly inversion results show similar resistivity and magnetism distribution characteristics to Lmg-1 (Figure9).   Combined with the geological map, the petrophysical properties of the study area and the results of the known exploration profiles, it is obvious that the high resistivity, low density, and low magnetism areas on both sides of the profiles are caused by the Baijiazuizi formation, the shallow medium-low resistivity, high density, and high magnetism areas are caused by the ore-bearing ultrabasic rock masses, and the transition zone on both sides of the middle abnormal area are the fault zones where ultrabasic rocks intrude into Baijiazuizi formation, as shown by the red dotted line in Figure 10. At the same time, we can see that the single geophysical method has limitations. For example, the resolution of CSAMT in the deep of the profile is insufficient, and it cannot clearly indicate the location of the deep abnormal target body. Therefore, we have roughly determined the target locations based on the high-density, high-magnetism, medium-low-resistivity geophysical exploration model, and the metallogenic law that is easy to form ore at a low-lying structure place, and we have inferred that the favorable metallogenic location of the line Lmg-1 is about 400-900 m and the burial depth is about 1100-1500 m, the favorable metallogenic area of the line Lmg-2 is about 800-1300 m and the burial depth is about 1200-1700 m, and the favorable metallogenic area of the line Lmg-3 is about 500-1000 m and the burial depth is about 1000-1500 m. The approximate location of the target body is shown by the black dotted line in Figure 10, and three drilling verification locations are designed.

Survey Results of the IV Mining Area
The known ultramafic rock mass in the fourth mining area is 1160 m long and completely covered by the Quaternary. The thickness of the cover layer is 60-140 m. The Baijiazuizi formation and ore-bearing ultrabasic rock mass are developed below the cover layer. The ore body is dominated by low grade ore, and No.1 ore body is the main ore body. It is produced in the concave section of the bottom of the rock body and is lenticular. The upper and lower parts are small and the middle is large. The strike of the Lmg-4 exploration line in the mining area is just south, passing through No.1 ore body, close to the No.10 exploration line. Gravity and CSAMT surveys were carried out along this line, their lengths are slightly different. Figure 11 is the inversion results of gravity and aeromagnetic anomalies based on the SL0 algorithm. The corresponding positions in the middle and lower parts of the profile have obvious large high-density and high-magnetic anomaly areas, and they all have the characteristics of extending to the upper left. The inversion result of CSAMT shows that the resistivity on the right side of the profile is relatively high, the resistivity in the middle is relatively low, and that their contact zone changes from steep to slow towards the deep (Figure12). Compared with the inversion results in Figure 11, it was found that the high-density and high-magnetic area is consistent with the middle medium-low resistivity area. The inversion results are interpreted in combination with a geological map, the results of a nearby geological exploration line, and the physical parameters of the rocks. It can be seen that the shallow low-density, low-magnetic and low-resistivity areas of the profile are mainly caused by the Quaternary, while the high-resistivity area on the right side is the reflection of Baijiazuizi formation, the high-density, high-magnetic, and low-resistivity area in the upper part of the contact zone corresponds to the known ore-bearing ultrabasic rocks in the four mining areas. At the same time, it is inferred that the transition zone is a fault zone (red dotted line on the right side of Figure 12), which is the channel for ultrabasic magma to intrude into Baijiazuizi formation. Based on the geophysical prospecting model and metallogenic law, we speculate that the deep high density, high magnetism, and medium-low resistivity area in a large range near the channel is the favorable target area for prospecting, and the location is roughly within the survey line 2600-3200 m and its buried depth is 1200-1700 m. At the same time, the drilling hole location is designed, as shown in Figure 12.

Discussion and Conclusions
The geophysical method is an important supporting technology for the exploration of deep metal mineral resources. In this paper, the deep exploration of the Jinchuan Cu-Ni sulfide deposit is carried out based on the gravity, aeromagnetic, and CSAMT methods with a complementary resolution. Among them, to overcome the ill-posed problem of inversion and reduce the multiplicity of solutions, the focus inversion method based on the SL0 norm is introduced to the inversion of gravity and aeromagnetic data, and the widely used SCS2D inversion software is used for the inversion of CSAMT data, both of which have achieved good inversion results. The medium-low resistivity, high density, and high magnetic areas shown in the inversion results can correspond to the known ore-bearing rocks in the shallow part well. In addition, four favorable target areas are delineated in the deep part of the mining area.
However, these methods are indirect and have some limitations. Because the ultrabasic rock is the parent rock and the surrounding rock of the ore body, the density and magnetism of the two are similar, and the ore body is often located in the middle or lower part of the ultrabasic rock body, so the geophysical signal produced by the ore body is easy to be covered by the ultrabasic rock body, meaning that it is difficult to distinguish the ultrabasic rock body and the ore body. Therefore, more geological data are needed to increase the reliability of target area prediction. In addition, the geophysical data in this paper are limited. It is expected to carry out seismic exploration in the study area and to obtain the velocity structure in the deep part of the mining area, so as to further supplement the supporting evidence of favorable prospective areas. At the same time, it is expected to carry out drilling work at the predicted favorable prospective areas, so as to further verify the results of geophysical deep exploration in this paper.
In general, it is of great significance to study the deep exploration of the Cu-Ni deposit in Jinchuan based on multiple geophysical methods. CSAMT, gravity, and aeromagnetic data not only indicate the known ore locations of the Ⅱ and Ⅳ mining areas but also indicate the favorable ore locations in the deep. At the same time, this achievement provides a good reference for further deep exploration of the mining area and its surrounding areas, as well as a good demonstration of the feasibility and effectiveness of the geophysical methods used to detect deep metal mines.