Half-space

Objetivo

With a simple half-space model we explore electric field, magnetic flux, and conduction currents in the domain due to the vertical magnetic dipole (VMD) source. By analyzing those fields with variable transmitting frequency, Tx height, conductivity, and susceptibility values, we understand how the inductive source exicte the earth. In addition, we discuss how the measured magnetic flux data behave as a function of frequency.

Oscillating magnetic flux is generated by injected sinusoidal current in the loop, and will generate induced electric field everywhere due to Faraday’s law:

(341)\[\nabla \times \mathbf{E} = -\imath \omega \mathbf{B}\]

Because conductivity of the air is zero (\(\sigma_{air} \approx 0\)), conduction currents in the earth is zero (\(\mathbf{J}_f=\sigma \mathbf{E}\)). From Ampere’s law:

(342)\[\nabla \times \mathbf{B} = \mu_0\mathbf{J} = \mu_0(\mathbf{J}_e + \mathbf{J}_f + \mathbf{J}_m)\]

where \(\mathbf{J}_m=\nabla \times \mathbf{M}\) is magnetic current density. Biot-Savart Law shown in Biot-Savart allows us to compute magnetic flux density from current density. For the free-space case, conduction and magnetic currents are both zero hence only the current source term \(\mathbf{J}_e\) will be taken account hence by applying Biot-Savart law to this we can obtain primary magnetic flux \(\mathbf{B}^{p}\). Similarly, by only taking into account conduction current, \(\mathbf{J}_f\), secondary magnetic flux \(\mathbf{B}^{s}\) can be evaluated. So if we define Biot-Savart operator \(G_{BS}[\cdot]\) which takes currents, then primary, secondary, and total magnetic fields can be defined as

(343)\[\mathbf{B}^{p} = G_{BS}[\mathbf{J}_e]\]

(344)\[\mathbf{B}^{s} = G_{BS}[\mathbf{J}_f + \mathbf{J}_m]\]

(345)\[\mathbf{B} =\mathbf{B}^{p}+\mathbf{B}^{s} =G_{BS}[\mathbf{J}_e+\mathbf{J}_f+\mathbf{J}_m]\]

Separating total magentic fields as primary and secondary will be useful for our following analyses to see impact of the eddy currents induced in the earth.

Figura 89 A loop source in a half-space earth. \(\sigma_{air}\) and \(\sigma_{earth}\) conductivity of the air and earth, respectively. \(\mu_{0}\) and \(\mu_{earth}\) susceptibility of the air and earth, respectively. \(z\) and offset correspondingly indicate Tx height (m) and Tx-Rx offset (m).

Using the setup shown in Figura 89, we perform a simulation, compute electric, magnetic flux, and current density. Here \(z\) = 30 m, \(\sigma_{earth}=0.01\) S/m, \(\mu_{earth}=\mu_0\), and \(f=10^3\) Hz. Figura 90 shows imaginary part of electric fielc and current density in \(y\)-direction.

Rotating electric field induced by oscillating magnetic field exist everywhere and show great amplitude at the near surface. Current density at an air region is zero, and this will make impact to the magnetic flux density. Figura 91 shows imaginary part of \(\mathbf{B}\). Imaginary part of magnetic flux is localized at the near surface, and mostly propagating downward.

Note

Imaginary part of \(\mathbf{B}\) only includes secondary field i.e. it is soely from the earth. This will be directly related to real part of \(\mathbf{E}\) from Eq. (341).

Figura 90 Real part of \(E_y\) a) and \(J_y\) b) in \(x\text{-}z\) plane. Postive and negative values indicate vector directing into the page and out of the page.

Figura 91 Imaginary part of \(\mathbf{B}\) in \(x\text{-}z\) plane

Different from the imaginary part of \(\mathbf{B}\), real part includes both primary and secondary magnetic flux. Figura 92 a and b show the real part of total and secondary magnetic flux density as a vector. Primary flux is dominant in \(\mathbf{B}\) hence signficantly localized close to Tx location, while secondary flux induced in the conductive earth are much more spreaded. Direction of the total (upward) and secondary flux (downward) is opposite in general. In practice, we cannot measure magnetic fields everywhere, but

Figura 92 Real part of total a) and secondary magnetic flux density b).

Pergunta

Compare imaginary and real part of \(\mathbf{B}_s\) Figura 91 and Figura 92 b. What is the main difference? Which component is more sensitive to the near surface?

Considering airborne EM survey or small-loop EM survey often

Figura 93 Meausured real component of \(B_z\) data as a function of frequency at Rx when Tx-Rx offset is 10 m. Total a) and secondary b) \(B_z\).

You can use below link to play with the app that we used to run a given simulation.

Using the app introduced, now we are going to run suite of simulations to see impact of four main factors to EM fields due to inductive source excitation:

• Frequency

• Tx height

• Conductivity

• Susceptibility

Frequency

With the same setup, now we vary frequency from 10 ⁶ to 10 ³. We may expect that the lower the frequency the deeper information. Figura 94 shows current density in \(y\)-direction at four different frequencies: 10 ⁵, 10 ⁴, 10 ³, and 10 ². As the frequency decreases, distribution of currents are more spreaded indicating the lower frequency has deeper depth information. Amplitude of the current increases as frequency increases.

Figura 94 Real part of \(E_y\) in \(x\text{-}z\) plane at 10 ⁵ Hz a), 10 ⁴ Hz a), 10 ³ Hz a), and 10 ² Hz d).

Tx height

Previously Tx height (\(z\)) was 30 m similar to airborne EM survey. Then what happens if we lower the \(z\) to 0 m? This can be considered as a small-loop EM survey such as EM-31 on the ground. Here, we use different \(\sigma\) = 0.1 S/m to emphasize difference in the results of two simulations, and fix frequency as 10 KHz. Figura 95 a and b respecively show the the imaginary part of \(B_z\) at \(z\) = 0 and \(z\) = 30 m. When the source is located right at the surface, magnetic flux is much more localized close to the near surface. However, when Tx height is increased to 30 m, magnetic flux has much distributed to both deeper depth and greater horizontal extent. This will indicate different resolution power in a) small-loop and b) airborne EM systems. First a small-loop EM survey will have much higher resolution to small near surface structures such as re-bar and UXO compared to an airborne EM survey. In contrast, an airborne EM survey will cover much greater volume of earth both in horizontal and vertical directions (i.e. greater depth of investigation). Similar difference can be recognized in the observed data. Figura 96 shows corresponding \(B_z\) data at a receiver location marked as red in Figura 95. Here the offset between Tx and Rx is 10 m. Changes in measured response are much more significant when \(z\) = 0 m compared to \(z\) = 30 m. Especially, when \(z\) = 0 m imaginary part of \(B_z\) shows sign reversal near 40 KHz , whereas when \(z\) = 30 m no sign reversal are occurred.

Figura 95 Imaginary part of \(B_z\) in \(x \text{-} z\) plane. At a) z = 0 m and b) z = 30 m.

Figura 96 Meausured imaginary component of \(B_z\) data as a function of frequency at Rx when Tx-Rx offset is 10 m.

Conductivity

Often a main goal of inductive source EM survey is imaging conductivity structure of the earth and for this, understanding how EM fields changes due to conductivity will be crucial. Here we assumed half-space model, hence we consider changes in a half-space conductivity value. Figura 97 a and b show real part of \(B_z\) component when \(\sigma_{earth}\) is 0.1 S/m and 0.01 S/m, respectively. When \(\sigma_{earth}\) is more conductive, EM fields decays faster, and hence Figura 97 a shows much localized distribution compared to Figura 97 b. This also indicates, conductive ground has much limited depth of investigation compared to resistive ground. Figura 98 show shows corresponding \(B_z\) data for both conductivity values.

Figura 97 Real part of \(B_z\)-direction in \(x\text{-}z\) plane. \(\sigma_{earth}\) = 0.1 S/m a) and 0.01 S/m b).

Figura 98 Meausured real component of \(B_z\) data as a function of frequency at Rx when Tx-Rx offset is 10 m. \(\sigma_{earth}\) = 0.1 S/m a) and 0.01 S/m b).

Susceptibility

We consider effect of magnetic currents shown in Eq. (345). Upon the definition of primary magnetic flux \(\mathbf{B}^p\), effects of susceptibility is not contained in \(\mathbf{B}^p\), but included in secondary flux \(\mathbf{B}^s\). Similar to the physics of magnetic survey, even static magnetic flux (zero frequency) will excite susceptible medium in the direction of primary magnetic flux indicating this secondary magnetic flux will be in the same direction of primary magnetic flux. Figura 99 a and b show real part of magnetic flux when \(\chi\)chi`=0, respectively. Difference between these two cases are drastic especially at near surface. And at location where Rx is located direction of the magnetic flux is opposite. Figura 100 shows measured \(B_z\) component for those two cases. When frequency is low, real part of \(B_z\) when \(\chi\) = 0.1 has positive sign, whereas the sign changes to negative when frequency gets higher. Imaginary part of \(B_z\) for two cases do not show significant difference indicating effects of susceptibility is minor in this component.

Pergunta

Why susceptibility effects is significant at low frequency, real part of \(B_z\) data?

Figura 99 Real part of \(B_z\)-direction in \(x\text{-}z\) plane. \(\chi_{earth}\) = 0.1 a) and 0 b). Frequency is 1000 Hz, and \(\sigma_{earth}\) is fixed to 0.01 S/m.

Figura 100 Meausured real component of \(B_z\) data as a function of frequency at Rx when Tx-Rx offset is 10 m. \(\chi_{earth}\) = 0.1 a) and 0 b). \(\sigma{earth}\) is fixed to 0.01 S/m.