from em_examples.DCIP_overburden_PseudoSection import *
from IPython.display import display
For a direct current resistivity (DCR) or Induced Polarization (IP) survey, currents are injected to the earth, and flow. Depending upon the conductivity and chargeability contrasts current flow in the earth will be distorted, and these changes can be measurable on the sufurface electrodes. Here, we focus on a bloc target embedded in a halfspace below a highly conductive surface layer, and investigate what are happening in the earth when static currents are injected. The conductive layer will also impact the illumination of the target (conductor or resistor). By investigating in data upon different physical properties contrasts, we explore the sensitivity of these surveys to the target.
test = DCIP2DfwdWidget();
display(test)
2D profiles are often plotted as pseudo-sections by extending $45^{\circ}$ lines downwards from the A-B and M-N midpoints and plotting the corresponding $\Delta V_{MN}$, $\rho_a$, or misfit value at the intersection of these lines as shown below. For pole-dipole or dipole-pole surveys the $45^{\circ}$ line is simply extended from the location of the pole. By using this method of plotting, the long offset electrodes plot deeper than those with short offsets. This provides a rough idea of the region sampled by each data point, but the vertical axis of a pseudo-section is not a true depth.
In the widget below the red dot marks the midpoint of the current dipole or the location of the A electrode location in a pole-dipole array while the green dots mark the midpoints of the potential dipoles or M electrode locations in a dipole-pole array. The blue dots then mark the location in the pseudo-section where the lines from Tx and Rx midpoints intersect and the data is plotted. By stepping through the Tx (current electrode pairs) using the slider you can see how the pseudo section is built up.
The figures shown below show how the points in a pseudo-section are plotted for pole-dipole, dipole-pole, and dipole-dipole arrays. The color coding of the dots match those shown in the widget.
surveyType = 'DipoleDipole'
survey, xzlocs = DC2Dsurvey(np.ones(mesh.nC),surveyType)
PseudoSectionWidget(survey,surveyType)
Widget Javascript not detected. It may not be installed or enabled properly.