Open access peer-reviewed chapter

Wide-band Rock and Ore Samples Complex Permittivity Measurement

Written By

Sixin Liu, Junjun Wu, Lili Zhang and Hang Dong

Submitted: October 20th, 2010 Reviewed: March 16th, 2011 Published: July 5th, 2011

DOI: 10.5772/17304

Chapter metrics overview

2,881 Chapter Downloads

View Full Metrics

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εr, its capacitance Cis εrtimes of vacuum capacitance C0. If we measure the capacitance before and after the capacitor is filled with the medium, we can obtain the permittivity. This method is mainly fitted for low-frequency range between 20MHz-200MHz.

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:

  1. 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;

  2. 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;

  3. altered K-feldspar granite, high-grade ore, low-grade ore from the Nianzigou Molybdenum mine, Chifeng, Inner Mongolian, China;

  4. albite rhyolite porphyry, breccia porphyry, quartz albitophyre, albitophyre copper, malachite copper oxide ore from the Qunji copper mine, Xinjiang, China;

  5. vesicular amygdaloidal andesite, massive diorite, andesitic copper ore from the Musi copper mine, Xinjiang, China;

  6. 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.

Figure 1.

Measurement flow chart

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. “ε1”and “μ1” are relative permittivity and relative magnetic conductivity of the medium filled in coaxial probe. “εs” and “μs” are permittivity and relative magnetic conductivity of the sample measured.

Figure 2.

Open-ended coaxial probe schematic diagram

It is assumed that only the TEMmode wave is traveling in coaxial probe and the time-harmonic factorejωtis ignored. The complex electric field and complex magnetic field of forward wave and backward wave in coaxial probe can be expressed respectively as,


where,r,φand zare the coordinate variable of cylindrical coordinate system in the equations.k1=ωμ0ε0μ1ε1,η1=μ0μ1/(ε0ε1).ε0, μ0are the permittivity and magnetic conductivity in the vacuum. ωis the angular frequency. Γis the reflection coefficient. A is the electric field amplitude of the forward wave on probe terminal surface.

The electromagnetic fieldErs,Hφsin the material measured can be expressed as the integral sum of all plane waves (including high-order mode) in spectral domain.


In the equations,γ=kc2-ks2, andRe(γ)0,ks=ωμ0ε0μsεs,Y(kc)=jωε0εs/γ, Γb(kc)=-exp(-2γd), J1(x)is first-order first-kind Bessel function, B(kc)is the field amplitude expressed in the spectral domain. kcis the continuous eigenvalue.

As the boundary conditions at the z=0plane are satisfied, the tangent components of the field are equal; following equations can be derived,


Multiply J1(kcr)rto two sides of equation (5), and then take an integral with the form0dr, we get,


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 δfunction:


Then the equation (8) can be modified:


SinceddxJ0(x)=-J1(x), the right of the equation (7) can be calculated:


Then the equation (7) can be rearranged:


Letkc=kc, the equation (13) becomes,


The equation (14) is substituted to equation (6) and takes an integral withabdrwith the equation (6):


Eliminating A and rearranging, (15) becomes,


The result is:

J0(x)is a Zero-order Bessel function of the first kind. It becomes the following:

Because the other parameters are all known. The equation (18) is a nonlinear equation about independent variable εsand variableΓ. The question of solving permittivity becomes a process of solving this nonlinear equation.

Figure 3.

Schematic diagram of rock sample holder

Figure 4.

VNA and the rock sample holder

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.

  1. 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.

  2. 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.

If Γmis the reflection coefficient obtained through measuring and Γais the practical reflection coefficient of probe terminal, Γmcan be expressed as (Blackham & Pollard, 1997):

where, edis the limited directivity error; eris frequency response error; esis equivalent source matching error. The reflection coefficient of the material Γacan be calculated through equation (17) or (18). Through measuring the reflection coefficient of three kinds of materialsΓm, the three equations abouted,er,escan be obtained. There are three variables and three equations, the error coefficientsed,er,escan be obtained.

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,Γa=-1; when the calibration material is air, the reflection coefficient of every frequency can be calculated through equation (17) or (18), because the permittivity of air is 1. According to the same theory, the reflection coefficient of de-ionized water can be calculated. Here, the reflection coefficient of water is obtained through the Cole-Cole formula.


where, εsis a direct current permittivity. εis an optical frequency permittivity. ω0is a Debye relaxation angle frequency. αis a Cole-Cole factor.

By substituting the reflection coefficient of air filmΓair_aandΓair_m, the reflection coefficient of de-ionized water Γwater_aand Γwater_mand the reflection coefficient of short-circuit Γshort_a=-1and Γshort_minto the equation (19) separately. We get,


where,C=Γwater_m-AA-B, A=Γair_m,B=Γshort_m.

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.

Figure 5.

Reflection coefficient before and after calibration

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

αis a weighting coefficient in this equation, ΓmandΓcare measured and the calculated reflection coefficients. The real part and the imaginary part should be treated equally to avoid that the large part dominates over the small part too much. This is the typical optimization problem. Here, εcan be thought as a complex-single variable. But the most mathematical software optimization tool can not process complex variable optimization question. So the complex permittivity is divided to real part and imaginary part. The variable x is a vector array, where,x1=Re(ε),x2=Im(ε). The selection of weighting coefficient is based on the numerous tests.

We solve the optimization process using the simplex method. The value of f(ε)after the optimization for every frequency is displayed in Fig. 6 for the material PTFE. It can be seen that the precision is very well. When the optimization stops, the objective function of minimum point satisfy the error requirement.

Figure 6.

The value of optimization objective function

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.

Figure 7.

Permittivity of PTFE sample

Figure 8.

Real part relative error of the permittivity of PTFE sample

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.

Figure 9.

Permittivity of the air

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.

Figure 10.

Permittivity of methanol

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.

Figure 11.

Photographs of the rocks and ores samples from 6 metal mines

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.

RocksRock or Ore namesFig. no.Measured permittivitySample number
3hybrid diorite11(1c)5-1016
4altered hornblende pyroxenite11(1d)5-177
5pyroxene peridotite11(1f)10-2010
6low-grade ore11(1g)9-239
7medium-grade ore11(1f)20-7010
8high-grade ore11(1h)5-9512

Table 1.

Rocks and ores from the Changren nickel-copper mine, Jilin, China

Figure 12.

Permittivity of marbles. (a) 8 Marble’s samples permittivities; (a) average of mable samples’ permittivities

Figure 13.

Averaged relative permittivities of rocks and ores from the Changren nickel-cooper mine

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.

RocksRock or Ore namesFig. nopermittivitySamples number
1tuffaceous fine-grained sandstone11(2a)5-75
2tuffaceous breccia11(2b)5-65
3dacitoid crystal tuff11(2c)5.5-813
5crystal tuff11(2e)5-7.511
6high-grade ore11(2f)10-7011
7medium-grade ore11(2g)10-125
8low-grade ore11(2h)5-108

Table 2.

Messages of the Huanghuagou lead-zinc mine, Chifeng, Inner Mongolian, China

Figure 14.

Averaging permittivities of ores and rocks from the Huanghuagou lead-zinc mine, Chifeng, Inner Mongolian, China

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.

RocksRock or Ore namesFig. nopermittivitySamples number
1altered K-feldspar granite11(3a)4.5-7.523 samples
2high-grade ore11(3b)5-157 (No: 02, 05, 07, 08, 09, 10, 11)
3low-grade ore11(3c)4-104 (No: 01, 03, 04, 06)

Table 3.

Messages of rocks and ores from the Nianzigou molybdenum mine, Chifeng, Inner Mongolian, China

Figure 15.

Averaging permittivities of the rocks and ores from the Nianzigou molybdenum mine, Chifeng, Inner Mongolian, China

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.

RocksRock or Ore namesFig. noPermittivitySamples number
1albite rhyolite porphyry (core)11( 4a)5-5.58
2breccia porphyry11(4b)5-5.514
3quartz albitophyre11(4c)5-7.511
4albitophyre ore11(4d)5-1016(No:01-09,11-17)
5malachite oxide ore11(4e)5-5.514(No.:01-14)

Table 4.

Messages of rocks and ores from the Qunji Copper mine, Xinjiang, China

Figure 16.

Average permittivities of rocks and ores from the Qunji Copper mine, Xinjiang, China

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.

Figure 17.

Averaging permittivities of rocks and ores from the Musi copper mine, Xinjiang, China

RocksRock or Ore namesFig. noPermittivitySamples number
1vesicular amygdaloidal andesite11(5a)5.5-12.519 samples
2massive diorite11(5b)7.5-1112 samples
3andesitic copper ore11(5c_1; 5c_2)5-1024 samples

Table 5.

Messages of rocks and ores from the Musi copper mine, Xinjiang, China

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.

RocksRock or Ore namesFig. noPermittivitySamples number
1glutenite11( 6a)8-146
2copper ore11(6b)8-4011
3lead-zinc ore11(6c)15-454

Table 6.

Messages of rocks and ores from the Zengnan copper mine, Xinjiang, China

Figure 18.

Averaging permittivities of rocks and ores from the Zengnan copper mine, Xinjiang, China


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.


  1. 1. BlackhamD. V.PollardR. D.19971952 An improved technique for permittivity measurements using a coaxial probe.IEEE Transactions on Instrumentation and Measurement,46510931099,0018-9456
  2. 2. Coutanceau-MonteilN.JacquinC.1993Improvements of the coaxial line technique for measuring complex dielectric permittivity of centimetric samples in the 20 to 1000 MHz range: application to sedimentary rocks.Log Analyst,3452133,1529-9074
  3. 3. DanielsD.2004Ground Penetrating Radar, 2nd Edition, The Institution of Engineering and Technology,0-86341-360-9United Kingdom.
  4. 4. FanS.StaebellkK.MisraD.19901952 Static analysis of an open-ended coaxial line terminated by layered media.IEEE Transactions on Instrumentation and Measurement,392435437,0018-9456
  5. 5. HoshinaS.KanaiY.MiyakawaM.2001A numerical study on the measurement region of a coaxial probe used for complex permittivity measurement.IEEE Trans. Magn,37533113314,0018-9464
  6. 6. JolH. M.2009Ground Penetrating Radar Theory and Applications,ELSEVIER,0-44453-348-6the Netherlands.
  7. 7. JordanB. P. .SheppardR. J.SzarnowskiS.1978The dielectric properties of formamide, ethanediol and methanol.J. Phys. D: Appl. Phys,115695701,0022-3727
  8. 8. KraszewskiA. W.1980Microwave aquametry- A review.Journal of Microwave Power,154209220,0022-2739
  9. 9. LiC. L.ChenK. M.1995Determination of electromagnetic properties of materials using flanged open-ended coaxial probe-full-wave analysis.IEEE Transactions on Instrumentation and Measurement,4411927,0018-9456
  10. 10. MisraD. K.1987A quasi-static analysis of open-ended coaxial lines.IEEE Transactions on Microwave Theory and Techniques,3510925928,0018-9480
  11. 11. NelsonS. O.BartleyP. G.1998Open-ended coaxial line permittivity measurements on pulverized materials.IEEE Transactions on Instrumentation and Measurement,471133137,0018-9456
  12. 12. NicolsonA. M.RossG.1970Measurement of intrinsic properties of Materials by time domain techniques.IEEE Trans on Instrument & Measurement,IM-19,4377382,0018-9456
  13. 13. NyshadhamA.SibbaldclC. L.StuchlyS. S.1992Permittivity measurements using open-ended sensors and reference liquid calibration--an uncertainty analysis.IEEE Transactions on Microwave Theory and Techniques,402305314,0018-9480
  14. 14. RobertsS. R.HippelA. Von1946A new method for measuring dielectric constant and loss in the range of centimeter wave.Journal of Applied Physics,177610616,0021-8979
  15. 15. ShenL. C.1985A laboratory technique for measuring dielectric properties of core samples at ultrahigh frequencies.Society of Petroleum Engineers Journal,254502514,0197-7520
  16. 16. ShiX. D.ShenL. C.1989Dielectric properties of reservoir rocks at very-high frequencies.ACTA GEOPHYSICA SINICA,32199110,0001-5733
  17. 17. ShenL. C.MarouniH.ZhangY. X.ShiX. D.1987Analysis of the parallel disk sample holder for dielectric permittivity measurement,IEEE transaction on geoscience and remote sensing,GE-25,5534540,0196-2892
  18. 18. StuchlyM. A.StuchlyS. S.1980Coaxial line reflection methods for measuring dielectric properties of biological substances at radio and microwave frequencie-a review.IEEE Trans. Instrum. Meas.,293176183,0018-9456
  19. 19. VenkateshM. S.RaghavanG. S. V.2005An overview of dielectric properties measuring techniques.CANADIAN BIOSYSTEMS ENGINEERING,4771530,1492-9058
  20. 20. WeirW. B.1974Automatic Measurement of Complex Dielectric Constant and Permeability at Microwave Frequency.Proc. IEEE,6213336,0018-9219
  21. 21. WeiW.SridharW.1991Radiation-corrected open-ended coax line technique for dielectric measurements of liquids up to 20 GHz.IEEE Transactions on Microwave Theory and Techniques,393526531,0018-9480
  22. 22. WuM. Z.YaoX.ZhangL. Y.2000An improved coaxial probe technique for measuring microwave permittivity of thin dielectric materials.Measurement Science and Technology,111116171622,0957-0233
  23. 23. ZhengH.SmithC. E.1991Permittivity measurement using a short open-ended coaxial line probe.IEEE Microw. Guided Wave Lett,111337339,1051-8207

Written By

Sixin Liu, Junjun Wu, Lili Zhang and Hang Dong

Submitted: October 20th, 2010 Reviewed: March 16th, 2011 Published: July 5th, 2011