from em_examples.CondUtils import ColeColePelton, vizColeCole
from em_examples.FreqtoTime import transFilt
import matplotlib
from ipywidgets import interact, FloatText, FloatSlider, ToggleButtons
from IPython.display import display
%pylab inline
matplotlib.rcParams['font.size'] = 16
Populating the interactive namespace from numpy and matplotlib
Using a simple Cole-Cole model, we parameterize complex resistivity with four parameters: resistivity at zero frequency ($\rho_0$), chargeability($\eta$), time constant ($\tau$), and frequency dependence ($c$). Based upon those parameters, we understand how resistivity and conductivity changes when medium is chargeable both in frequency domain and time domain.
Pelton's Cole-Cole model for resistivity and conductivity can be written as
$$ \rho(\omega) = \rho_0 \Big[1 - \eta \Big(1-\frac{1}{1+(\imath\omega\tau)^c}\Big) \Big] $$and
$$ \sigma(\omega) = \sigma_{\infty}\Big(1-\frac{\eta}{1+(1-\eta)(\imath\omega\tau)^c} \Big) $$respectively.
$\sigma_1$: Conductivity of the first layer (S/m)
$\sigma_2$: Conductivity of the first layer (S/m)
$f$ (Hz): Frequency (Hz)
Type:
interact(vizColeCole, eta=FloatSlider(min=0.1, max=0.5, step=0.05, value=0.4),
tau=FloatText(value=0.1),
c=FloatSlider(min=0.1, max=1., step=0.1, value=0.5),
sigres = ToggleButtons(options=['sigma','resis']),
t1=FloatText(value=800),
t2=FloatText(value=1400),
)
<function em_examples.CondUtils.vizColeCole>