Rocks and ores from the Changren nickel-copper mine, Jilin, China
1. Introduction
Ground penetrating radar (GPR) is based on high-frequency electromagnetic wave propagation and its detecting targets are below the ground surface (Daniels, 2004; Jol, 2009). Velocity and attenuation are two important factors describing the electromagnetic wave in the media composed of rocks or soils. Velocity is inverse proportional to the root of the permittivity while other parameters are fixed generally. In order to understand the performance of GPR, permittivity testing and analysis are critical. In addition, during the metal ore exploration by borehole radar which is an operating mode of GPR, the permittivity difference between ore body and surrounding rock is the foundation for exploration. As the sampling site, geological environment, and seasons are changing, the permittivity are different even for the same rock. Therefore, permittivity measurement is very important.
Currently, measurement methods are basically indirect methods which are based on transmission line theory, characteristic impedance, and propagation constant. These variables have intrinsic relationship with permittivity which can be inverted from the measured data by certain calculation procedure. At the RF frequency band, common measurement methods include short-circuited wave-guide measurement, coaxial line transmission/reflection method, open-ended coaxial probe, resonant cavity method, free-space transmission technique, parallel-plate capacitance method, etc.
Roberts and von Hippel (1946) developed the short-circuited wave-guide measurement, sample is inserted at the end of the wave-guide or coaxial line, the standing wave is formed as the incident wave and the reflected wave coexist in the wave-guide. The sanding wave ratios (SWR’s) were required to measure in the case with and without sample. Permittivity can be determined by the change in the widths of nodes, sample length, and the waveguide dimension.
Resonant cavity method is a perturbation technique, which is frequently used for measuring permittivity because of its simplicity, accuracy, and high temperature capability (Venkatesh & Raghavan, 2005). This technique is based on the resonant frequency shift, and the change in absorption characteristics due to the insertion of sample material. The measurement is made by placing a sample completely in the center of a waveguide. The size of the cavity is designed for special frequency. Calibration is required before measurement, once the calibration is finished, the measurement is fast. The samples preparation is little difficult.
Free-space transmission technique is carried out by placing the sample between the transmitting and the receiving antennas (Kraszewski, 1980). The electrical parameters can be calculated by measuring the attenuation and the phase shift in the media. Free-space measurement can measure permittivity in a wide frequency band. It is a no-destructive and contact-less method. The frequency can be as high as 30GHz. The usual assumption for this technique is that a uniform plane wave incident upon the flat surface of a homogeneous material, and the planar sample has infinite extent laterally (Venkatesh & Raghavan, 2005).
Parallel disk capacitor technique (Shi & Shen, 1989; Shen, Marouni, Zhang, & shi, 1987) is based on the theory that the capacitance is dependant on the medium permittivity filling the capacitor. As a capacitor is filled with medium having permittivity of
Coaxial reflection/transmission method was introduced by Nicolson, Rose, and Wire, in 1970s, therefore, it is also called NRW method (Nicolson &Ross, 1970; Weir, 1974). The material under test (MUT) is inserted into coaxial line (waveguide). As the electromagnetic wave is traveling in the line, it meets the MUT, a part of the wave is transmitted, and the other is reflected. Attenuation and phase shift occur at the same time. Vector network analyzer (VNA) can measure the reflection and transmission coefficients at a wide frequency range. The permittivity can be inverted by electromagnetic theory from these measured data. This technique was first used for rock measurement (Shen, 1985), and the water-bearing, oil-bearing, or gas-bearing rocks can be distinguished from the measured data. Coutanceau-Monteil and Jacquin improved this technique in 1993. The rock samples are coated with a low melting point alloy and this improvement enables measurements to be made on larger samples and avoids air-gap correction at the solid-cell interface. However, the sample preparation is still difficult.
Open-ended coaxial probe technique (Zhen &Smith, 1991; Blackham&Pollard, 1997; Hoshina, Kanai, Miyakawa, 2001) is a technique of which the open-end of a coaxial is touched on the MUT, and the other end is connected to a VNA. The measured reflection coefficient is used to invert the permittivity. A flange is connected to the open-end to improve the accuracy. The circular waveguide probe, the rectangular waveguide probe, and the coaxial probe are often used. The coaxial probe can be used at frequency-domain, time-domain, point frequency, and it is suitable for kinds of media with different electric characteristic (low to high permittivity, lossy or lossless, magnetic lossy materials). This method is no-destructive, no-intrusive, and also characterized with wide frequency band, simple and open structure. It is fitted for on-line, in-vivo, and in-situ testing, it is receiving the attention from many scientists and engineers.
Open-ended probe electromagnetic radiation dates back to 1950s or earlier. Open-ended measurement technique was improved rapidly in recent 20 years (Stuchly M.A. & Stutch, S.S, 1980). It was first introduced and applied for biological tissue measurement at RF frequency, and was used for dielectric material testing at RF band shortly, and the lumped parameter model is used at the time for data processing. Later, as the testing frequency expanded to X and Ku band, quisi-static analysis model (Misra, 1987; Fans et al., 1952; Nyshadham et al., 1992; Wei, & Sridhar, 1991; Blackham & Pollard, 1997; Nelson& Bartley, 1998) and full-wave analysis model (Li & Chen, 1995) appeared, and the testing material also expanded to high permittivity from low permittivity. But the sample thickness was infinite. Until the early stage of 1990s, spectral-domain quisi-static analysis model appeared, and it was suitable for multi-layered material with finite thickness. Then, many researchers analyzed the high modes in the coaxial line, and reached the spectral-domain full-wave analysis model which is a complete mathematic model.
Several rules must be considered before we choose the techniques for measurement, such as frequency range preferred, sample preparation easiness, the equipments accessible. We choose open-ended coaxial probe technique backed by metal plate because the model is easy to be analyzed, the samples making is relatively easy, and the VNA is accessible for us.
We used open-ended coaxial probe to measure 342 samples within 31 categories from 6 mine sites, including:
granite, marble, hybrid diorite, altered hornblende pyroxenite, pyroxene peridotite, low-grade ore, medium-grade ore, high-grade ore from the Changren nickel-cupper mine, Jilin province, China;
tuffaceous fine-grained sandstone, tuffaceous breccia, dacitoid crystal tuff, dacite, crystal tuff, low-grade ore, medium-grade ore, high-grade ore, and pyrite from the Huanghuagou lead-zinc mine, Chifeng, Inner Mongolian, China;
altered K-feldspar granite, high-grade ore, low-grade ore from the Nianzigou Molybdenum mine, Chifeng, Inner Mongolian, China;
albite rhyolite porphyry, breccia porphyry, quartz albitophyre, albitophyre copper, malachite copper oxide ore from the Qunji copper mine, Xinjiang, China;
vesicular amygdaloidal andesite, massive diorite, andesitic copper ore from the Musi copper mine, Xinjiang, China;
glutenite, copper ore, lead-zinc ore from the Zengnan copper mine, Xinjiang, China.
It is found that high-grade ore can be distinguished from surrounding rocks by permittivity, low-grade ore shows similar permittivity as other veins. We first measured and analyzed so many samples including different rocks and ores at wide frequency range according to current references world widely.
2. Open-ended coaxial technique
Firstly, the theory of open-coaxial probe technique is introduced. Then, the single-layered short-circuit spectral domain quasi-static model is established. Finally, the permittivity is calculated through the optimized inversion procedure.
2.1. The open-coaxial probe technique
The process of open-ended coaxial probe measuring medium includes three steps.
The first step is experimental testing (hardware operation). Reflection coefficient, which contains electromagnetic parameters information of the material, is measured by an open-coaxial probe with a flange-plane connected to reflection equipment which may be vector reflectometry (vector network analyzer or six ports) or scalar reflectometry. The vector network analyzer is used here (Agilent E5071B). The second step is the theoretical modeling. The physics model for the measurement is built based on electromagnetic knowledge. Namely, the relationship between the reflection coefficient of open-coaxial probe and the electromagnetic parameters of the sample measured is usually complicated and nonlinear one. And this modeling process is also called forward process. The third step is about the rebuilding electromagnetic parameters. In this step, the electromagnetic parameters corresponding with reflection coefficient of the material measured are obtained. This process is often called inverse process. Because the equations established are often nonlinear, we often use numerical solution.
Both theoretical modeling and numerical inversion are important steps in electromagnetic parameters measurement technology. Their accuracy decides the reliability of the measured results. For open-coaxial probe technique, modeling process is a mathematics analysis process, while the inversion is a process of appropriate optimization iteration.
2.2. Spectral-domain quasi-static reflection coefficient model
At present, popular analysis models used are full wave model and quasi-static model. The former is a theory model solved according to electromagnetic field. But the latter ignores the influences of high-order mode to reflection coefficient near the aperture of a coaxial probe. In theory, full wave model can accurately measure dielectric properties of materials. In fact, some theories and test indicate that the model error introduced by ignoring high-order mode is equal to the measuring error induced by the hardware and the operation. So, we can use quasi-static model (Wu et al., 2000), if the measurement error cannot be ignored and the rapid data processing is wanted.
We use spectral-domain quasi-static method here. The flow chart for the measurement is as shown in Fig. 1. Firstly, the reflection coefficient of a MUT is measured by a VNA. Then the measured reflection coefficient is calibrated (or corrected) to achieve real reflection coefficient. Finally the permittivity is calculated through inversion from the reflection coefficient.
We used a single-layered short-circuit model here. Fig. 2 shows a schematic diagram of open-ended coaxial probe. “a” is the radius of inner conductor of coaxial probe. “b” is the radius of outer conductor. “d” is the thickness of sample measured. “
It is assumed that only the
where,
The electromagnetic field
In the equations,
As the boundary conditions at the
Multiply
The integral order of the left of equation (7) is rearranged as,
According to the orthogonality of Bessel functions, we get
According to the sifting properties of
Then the equation (8) can be modified:
Since
Then the equation (7) can be rearranged:
Let
The equation (14) is substituted to equation (6) and takes an integral with
Eliminating A and rearranging, (15) becomes,
The result is:
Because the other parameters are all known. The equation (18) is a nonlinear equation about independent variable
To avoid the influence of air-gap to the testing result, a rock sample holder is make as shown in fig. 3.
The practical testing equipment including a VNA is shown in Fig. 4.
In equation (17) and (18), because the Bessel function has oscillating property, the main difficulty focuses on the Bessel function with integral variable. Obviously these nonlinear equations have no analytical solution. So we uses numerical solution here. A big number (12500) is used for the positive infinite of the upper limit during the numerical intergal. The value of the big number is determined by different testing of many conditions.
2.3. Calibration
The calibration in this measurement includes two steps, one is the transmission line calibration, and the other is the probe calibration.
2.3.1. Transmission line calibration
The VNA is a exact device and is connected to coaxial probe through a coaxial cable. According to the operation requirement of the network analyzer, the coaxial line is calibrated using calibrating kits. The detailed operation process as follows.
VNA parameters setting. The frequency range is between 1MHz-1GHz in this test according to our request. Power can be selected by our need, for example, 0dBm. Large power is believed to sense large sample volume. We choose 1000 sampling points here.
Calibration process. We use single port calibration here because we only measure S11. The calibration kits which including the device SHORT, LOAD, and OPEN are used to calibrate the VNA. After this process, the reference plane for VNA is at the end of the coaxial cable. However, because the reflection surface is on the flange surface, not the end of the cable, further calibration is still needed.
2.3.2. Probe calibration
Probe calibration is an indirect method. We use short-circuit, the air, and the de-ionized water to calibrate the probe.
Ifwhere,
Short-circuit, and air are ideal calibration materials. The third material must have known permittivity. The de-ionized water is selected as the third calibration material here. When it is of short-circuit,
where,
By substituting the reflection coefficient of air film
where,
We can also get,
It can be concluded based on the equation (19) that:
Equation (24) determines the second step calibration. Fig. 7 shows the comparison among results before and after calibration for PTFE and de-ionized water, separately.
It can be seen that the real part of PTFE measured can be calibrated to around but different from 1.
2.4. Inversion calculation and error evaluation
If the permittivity of a material measured is known, the interface reflection coefficient (or admittance) can be calculated. This process is a forward one. The reverse process can be solved numerically. The following equation can be obtained from equation (18),
The solution is:
Because it is a complex calculation, the objective function is defined as
We solve the optimization process using the simplex method. The value of
We testify this technique using a standard material PTFE, air, and methanol.
We first test this technique with PTFE whose thickness is 10.50mm in this paper and has permittivity of 2.1-j0.0004 (Li & Chen, 1995) in microwave band. Because the imaginary part can not be measured exactly for lowly lossy medium (Wu et al., 2000) by this technique, we ignore the analysis for the imaginary part. The inverted permittivity is displayed in Fig. 7. The real part relative error at every frequency is displayed in Fg. 8.
We noticed that the arisen relative error is within 5% basically. The average relative error is 1.2749%. One of the many reasons leading to the error is the air gap between the flange and the sample. The main reasons of producing air gap are that the upper surface and down surface are not parallel and clean enough, and the upper surface and the down surface do not touch enough with coaxial probe flange-plane and short-circuit board, although we already tried our best.
The permittivity calculated by the air film is displayed in Fig. 9.
The relative error is 0.7692%. Because the air is a kind of calibration material, the permittivity of air calculated should be theoretical value 1.The relative error is below 0.8%. It proves the validity of inversion process.
The measured permittivity for methanol is displayed in Fig. 10.
The measured permittivity for methanol is compared with the theoritical values which is calculated by the debye equation or cole-cole equation (Jordan et al., 1978) as shown in Fig. 10. The measured data is accetable except that they have clear difference with the theotitical ones at high frequency range. The reducement of this error could be the future topic.
3. Measured results and analysis
342 rocks and ores sample within 31 categories from 6 mines are measured and analyzed in this part by using open-coaxial probe technique. The photos for these rocks and ores samples are shown in Fig. 11.
3.1. Samples from the Changren nickel-copper mine, Jilin, China
Table 1 shows the messages of rocks and ores from the Changren nickel-copper mine, Jilin, China.
Fig.12 shows marbles permittivities as an example, the solid and the dashed lines denote the real parts and the imagery parts. We find the values are diverse for the same rock. We think this kind of diversity is due to the fact of that the probe senses a small range and the samples are in-homogeneous. Therefore, we use the averaging value of these data to represent this sample, because the averaging could reflect the total characteristic.
Fig. 13 shows the average permittivities of all rocks and ores from the Changren nickel-cooper mine, China. We find high grade ore and medium grade ore have highest values, then the values range from high to low are the pyroxene peridotite, low grade ore, light alterative bornblende pyroxenite, marble, hybrid diorite, granitization granite.
Rocks | Rock or Ore names | Fig. no. | Measured permittivity | Sample number |
1 | granite | 11(1a) | 5-7.5 | 25 |
2 | marble | 11(1b) | 5-10 | 8 |
3 | hybrid diorite | 11(1c) | 5-10 | 16 |
4 | altered hornblende pyroxenite | 11(1d) | 5-17 | 7 |
5 | pyroxene peridotite | 11(1f) | 10-20 | 10 |
ores: | ||||
6 | low-grade ore | 11(1g) | 9-23 | 9 |
7 | medium-grade ore | 11(1f) | 20-70 | 10 |
8 | high-grade ore | 11(1h) | 5-95 | 12 |
Actually, the pyroxene peridotite, light alterative bornblende pyroxenite are basic rocks and ultra-basic rock which were ore carrier. When ore’s grade is low, the permittivity represents the carrier rock’s property. These basic rocks and ultra-basic rock come from tectonic emplacement. The granitized granite is the host rock which has distinguished lower values. These measured data show optimistic aspect for borehole radar detection for metal ore-body.
3.2. The samples from the Huanghuagou lead-zinc mine, Chifeng, Inner Mongolian, China
The table 2 shows the message of rocks and ores from the Huanghuagou Lead-Zinc mine Chifeng, China. Ores and rocks ranked by permittivity from high to low are high-grade ore, pyrite, medium-grade ore, dacitoid crystal tuff, low-grade ore, crystal tuff, tuffaceous breccia, tuffaceous sandstone, and dacite. The high-grade ore, pyrite, and the medium-grade ore are distinguishable from each other and the others.
Rocks | Rock or Ore names | Fig. no | permittivity | Samples number |
1 | tuffaceous fine-grained sandstone | 11(2a) | 5-7 | 5 |
2 | tuffaceous breccia | 11(2b) | 5-6 | 5 |
3 | dacitoid crystal tuff | 11(2c) | 5.5-8 | 13 |
4 | dacite | 11(2d) | 5.5-6 | 10 |
5 | crystal tuff | 11(2e) | 5-7.5 | 11 |
Ores | ||||
6 | high-grade ore | 11(2f) | 10-70 | 11 |
7 | medium-grade ore | 11(2g) | 10-12 | 5 |
8 | low-grade ore | 11(2h) | 5-10 | 8 |
9 | pyrite | 11(2i) | 20-40 | 4 |
3.3. Samples from the Nianzigou molybdenum mine, Chifeng, Inner Mongolian, China
The table 3 shows the messages of rocks and ores from the Nianzigou molybdenum mine, Chifeng, Inner Mongolian, China. Ores and rocks ranked by permittivity from high to low are high-grade ore, low-grade ore, and altered K-feldspar granite. The high-gride ore is distinguishable from other two, and the low-grade ore shows the nearly same permittivity as altered K-feldspar grinate.
Rocks | Rock or Ore names | Fig. no | permittivity | Samples number |
1 | altered K-feldspar granite | 11(3a) | 4.5-7.5 | 23 samples |
Ores | ||||
2 | high-grade ore | 11(3b) | 5-15 | 7 (No: 02, 05, 07, 08, 09, 10, 11) |
3 | low-grade ore | 11(3c) | 4-10 | 4 (No: 01, 03, 04, 06) |
3.4. Samples from the Qunji copper mine, Xinjiang, China
The table 4 shows the messages of rocks and ores from the Qunji Copper mine, Xinjiang, China. Ores and rocks ranked by permittivity from high to low are albitophyre ore, quartz albitophyre, breccia porphyry, malachite copper oxide ore, and albite rhyolite porphyry. The albitophyre ore is clearly distinguishable from the others in the real part. Other rocks and ore are ambitious in permittivity.
Rocks | Rock or Ore names | Fig. no | Permittivity | Samples number |
1 | albite rhyolite porphyry (core) | 11( 4a) | 5-5.5 | 8 |
2 | breccia porphyry | 11(4b) | 5-5.5 | 14 |
3 | quartz albitophyre | 11(4c) | 5-7.5 | 11 |
Ores | ||||
4 | albitophyre ore | 11(4d) | 5-10 | 16(No:01-09,11-17) |
5 | malachite oxide ore | 11(4e) | 5-5.5 | 14(No.:01-14) |
3.5. Samples from the Musi copper mine, Xinjiang, China
The table 5 shows the messages of rocks and ores from the Musi copper mine, Xinjiang, China. Ores and rocks ranked by permittivity from high to low are vesicular amygdaloidal andesite, massive diorite, and andesitic copper ore. The andesitic copper ore is distinguishable from the others and shows low permittivity characteristic which is opposite to other mines.
Rocks | Rock or Ore names | Fig. no | Permittivity | Samples number |
1 | vesicular amygdaloidal andesite | 11(5a) | 5.5-12.5 | 19 samples |
2 | massive diorite | 11(5b) | 7.5-11 | 12 samples |
Ores | Total:24 | |||
3 | andesitic copper ore | 11(5c_1; 5c_2) | 5-10 | 24 samples |
3.6. Samples from the Zengnan copper mine, Xinjiang, China
The table 6 shows the messages of rocks and ores from the Zengnan copper mine, Xinjiang, China. Ores and rocks ranked by permittivity from high to low are lead-zinc ore, copper ore, and glutenite. Three of them can be distinguished from each other.
Rocks | Rock or Ore names | Fig. no | Permittivity | Samples number |
1 | glutenite | 11( 6a) | 8-14 | 6 |
Ores | ||||
2 | copper ore | 11(6b) | 8-40 | 11 |
3 | lead-zinc ore | 11(6c) | 15-45 | 4 |
4. Conclusion
Open-ended coaxial technique can measure the permittivity in wide frequency range quickly. The sample machining is relatively simple, and only the smooth surfaces of the sample sheets are required. Because the sensing range of the probe concentrates mainly at the center of the probe and the samples measured are no so homogeneous, we use averaging value from several samples of a rock or ore to reduce the random effect due to their in-homogeneity. It is shown that permittivity of metal ore is higher than other rocks, and high-grade ore is distinguishable from surrounding rocks. These measurements provide insights into the wide-frequency permittivity of metal ores and rocks, and also provide basis for electromagnetic exploration by borehole radar.
There are still couple of problems with the current research. The sizes of the flange, the aperture of the probe, sheet sample thickness, are not optimized yet. The sensing area for the current probe is small for the inhomogeneous rocks and ores. These are all future works for us.
References
- 1.
Blackham D. V. Pollard R. D. 1997 1952 An improved technique for permittivity measurements using a coaxial probe . ,46 5 1093 1099 ,0018-9456 - 2.
Coutanceau-Monteil N. Jacquin C. 1993 Improvements of the coaxial line technique for measuring complex dielectric permittivity of centimetric samples in the 20 to 1000 MHz range: application to sedimentary rocks. ,34 5 21 33 ,1529-9074 - 3.
Daniels D. 2004 Ground Penetrating Radar, 2nd Edition , The Institution of Engineering and Technology,0-86341-360-9 United Kingdom. - 4.
Fan S. Staebellk K. Misra D. 1990 1952 Static analysis of an open-ended coaxial line terminated by layered media . ,39 2 435 437 ,0018-9456 - 5.
Hoshina S. Kanai Y. Miyakawa M. 2001 A numerical study on the measurement region of a coaxial probe used for complex permittivity measurement. ,37 5 3311 3314 ,0018-9464 - 6.
Jol H. M. 2009 Ground Penetrating Radar Theory and Applications, ELSEVIER,0-44453-348-6 the Netherlands. - 7.
Jordan B. P. . Sheppard R. J. Szarnowski S. 1978 The dielectric properties of formamide, ethanediol and methanol . ,11 5 695 701 ,0022-3727 - 8.
Kraszewski A. W. 1980 Microwave aquametry- A review . ,15 4 209 220 ,0022-2739 - 9.
Li C. L. Chen K. M. 1995 Determination of electromagnetic properties of materials using flanged open-ended coaxial probe-full-wave analysis . ,44 1 19 27 ,0018-9456 - 10.
Misra D. K. 1987 A quasi-static analysis of open-ended coaxial lines. ,35 10 925 928 ,0018-9480 - 11.
Nelson S. O. Bartley P. G. 1998 Open-ended coaxial line permittivity measurements on pulverized materials . ,47 1 133 137 ,0018-9456 - 12.
Nicolson A. M. Ross G. 1970 Measurement of intrinsic properties of Materials by time domain techniques. ,IM-19 ,4 377 382 ,0018-9456 - 13.
Nyshadham A. Sibbaldcl C. L. Stuchly S. S. 1992 Permittivity measurements using open-ended sensors and reference liquid calibration--an uncertainty analysis . ,40 2 305 314 ,0018-9480 - 14.
Roberts S. R. Hippel A. Von 1946 A new method for measuring dielectric constant and loss in the range of centimeter wave. ,17 7 610 616 ,0021-8979 - 15.
Shen L. C. 1985 A laboratory technique for measuring dielectric properties of core samples at ultrahigh frequencies . ,25 4 502 514 ,0197-7520 - 16.
Shi X. D. Shen L. C. 1989 Dielectric properties of reservoir rocks at very-high frequencies. ,32 1 99 110 ,0001-5733 - 17.
Shen L. C. Marouni H. Zhang Y. X. Shi X. D. 1987 Analysis of the parallel disk sample holder for dielectric permittivity measurement , ,GE-25 ,5 534 540 ,0196-2892 - 18.
Stuchly M. A. Stuchly S. S. 1980 Coaxial line reflection methods for measuring dielectric properties of biological substances at radio and microwave frequencie-a review. .,29 3 176 183 ,0018-9456 - 19.
Venkatesh M. S. Raghavan G. S. V. 2005 An overview of dielectric properties measuring techniques. ,47 7 15 30 ,1492-9058 - 20.
Weir W. B. 1974 Automatic Measurement of Complex Dielectric Constant and Permeability at Microwave Frequency. ,62 1 33 36 ,0018-9219 - 21.
Wei W. Sridhar W. 1991 Radiation-corrected open-ended coax line technique for dielectric measurements of liquids up to 20 GHz . ,39 3 526 531 ,0018-9480 - 22.
Wu M. Z. Yao X. Zhang L. Y. 2000 An improved coaxial probe technique for measuring microwave permittivity of thin dielectric materials. ,11 11 1617 1622 ,0957-0233 - 23.
Zheng H. Smith C. E. 1991 Permittivity measurement using a short open-ended coaxial line probe. ,1 11 337 339 ,1051-8207