from geoscilabs.em.TDEMHorizontalLoopCylWidget import TDEMHorizontalLoopCylWidget
APP = TDEMHorizontalLoopCylWidget()
from matplotlib import rcParams
rcParams['font.size'] = 16
Here, we show the transient fields and fluxes that result from placing a vertical magnetic dipole (VMD) source over a layered Earth. The transient response in this case refers to the fields and fluxes that are produced once a long-standing primary magnetic field is removed.
There are two commonly used models for describing the VMD source that produces a transient response: 1) as an infinitessimally small bar magnet that experiences a long-standing vertical magnetization which is then instantaneously removed at $t=0$, and 2) as an infinitessimally small horizontal loop of wire carrying a constant current which is then instantaneously shut off at $t=0$ (step-off current waveform).
True dipole sources do not exist in nature however they can be approximated in practice. For geophysical applications, we use small current loops to approximate transient VMD sources. These EM sources may be placed on the Earth's surface (ground-based surveys) or flown through the air (airborne surveys). According to the Biot-Savart law, a primary magnetic field is produced whenever there is current in the loop. When the current is shut-off, the sudden change in magnetic flux induces anomalous currents in the Earth which propagate and diffuse over time. The distribution and propagation of the induced currents depends on the subsurface conductivity distribution and how much time has passed since the current in the VMD source was shut off. The induced currents ultimately produce secondary magnetic fields which can be measured by a receiver.
In this app, we explore the following:
The geological scenario being modeled is shown in the figure below. Here, we assume the Earth is comprised of 3 layers. Each layer can have a different electrical conductivity ($\sigma$). However, a constant magnetic susceptibility ($\chi$) is used for all layers; where $\mu_0$ is the magnetic permeability of free space and $\mu = \mu_0 (1 +\chi)$. The thicknesses of the top two layers are given by $h_1$ and $h_2$, respectively.
In this case, a transient VMD source (Tx) is used to excite the Earth, and the Earth's TEM response (secondary magnetic field) is measured by a receiver (Rx). In practice, the transmitter and receiver may be placed near the Earth's surface or in the air. The receiver may also measure secondary fields at a variety of times after the source is shut off.
To understand the fields and currents resulting from a transient VMD source we have two apps:
Follow the exercise in a linear fashion. Some questions may use parameters set in a previous question.
Q1: Set $\sigma_1$, $\sigma_2$ and $\sigma_3$ to arbitrary conductivity values. Based on the geometry of the problem, which components (x, y, z) of each field (E, B, dBdt or J) are zero? Run the Fields app and set AmpDir = None. Next, try different combinations of Field and Comp. Does the simulation match what you initially thought?
Q2: Re-run the Fields app to set parameters back to default. What happens to the Ey and Jy as you increase time index starting at 1? How does the diffusion and propagation of the EM signal change if the conductivity of all the layers is increased to 1 S/m?
Q3: Re-run the Fields app to set parameters back to default. Set $\sigma_1 = 0.01$ S/m, $\sigma_2 = 1$ S/m and $\sigma_3 = 0.01$ S/m. Now increase time index starting at 1. Is the signal able to effectively penetrate the conductive layer? Why/why not? What if the layer was resistive (i.e. $\sigma_2 = 0.0001$ S/m) instead?
Q4: Repeat Q3 but examine the current density. Where is the highest concentration of current density at late time channels? Does this support your answer to Q3?
Q5: Re-run the Fields app to set parameters back to default. Set Field = B, AmpDir = Direction. What happens to the magnetic flux density as the time index is increased starting at 1? At (x,z)=(0,0), what is the vector direction of the magnetic flux density? Repeat Q5 for dBdt.
Q6: Re-run the Fields app to set parameters back to default. Set $\sigma_1 = 0.01$ S/m, $\sigma_2 = 1$ S/m and $\sigma_3 = 0.01$ S/m. Examine how B and dBdt are impacted by the conductive layer.
We use this app to simulate the fields and currents everywhere due to a transient VMD source. The fields and induced currents depend on time and the subsurface conductivity distribution. You will use the app to change various parameters in the model and see how the fields and currents change.
APP.InteractivePlane_Layer()
Using this app, we show how the fields observed at the receiver location depend on the parameters set in the previous app. Note that if you want to see changes in the data due to changes in the model, you MUST re-run the previous app.
APP.InteractiveData_Layer()
EM fields will be depenent upon a number of parameters, using a simple half-space model ($\sigma_1=\sigma_2=\sigma_3$) explore how EM fields and data changes upon below four parameters.
E1: Effects of frequency?
E2: Effects of Tx height?
E3: Effects of Conductivity?
E4: Effects of Susceptibility?