dipole_k¶
-
empymod.model.
dipole_k
(src, rec, depth, res, freq, wavenumber, ab=11, aniso=None, epermH=None, epermV=None, mpermH=None, mpermV=None, **kwargs)[source]¶ Return electromagnetic wavenumber-domain field.
Calculate the electromagnetic wavenumber-domain field due to infinitesimal small electric or magnetic dipole source(s), measured by infinitesimal small electric or magnetic dipole receiver(s); sources and receivers are directed along the principal directions x, y, or z, and all sources are at the same depth, as well as all receivers are at the same depth.
- Parameters
- src, reclist of floats or arrays
Source and receiver coordinates [x, y, z] (m):
The x- and y-coordinates can be arrays, z is a single value.
The x- and y-coordinates must have the same dimension.
The x- and y-coordinates only matter for the angle-dependent factor.
Sources or receivers placed on a layer interface are considered in the upper layer.
- depthlist
Absolute layer interfaces z (m); #depth = #res - 1 (excluding +/- infinity).
- resarray_like
Horizontal resistivities rho_h (Ohm.m); #res = #depth + 1.
- freqarray_like
Frequencies f (Hz), used to calculate etaH/V and zetaH/V.
- wavenumberarray
Wavenumbers lambda (1/m)
- abint, optional
Source-receiver configuration, defaults to 11.
electric source
magnetic source
x
y
z
x
y
z
electric
receiver
x
11
12
13
14
15
16
y
21
22
23
24
25
26
z
31
32
33
34
35
36
magnetic
receiver
x
41
42
43
44
45
46
y
51
52
53
54
55
56
z
61
62
63
64
65
66
- anisoarray_like, optional
Anisotropies lambda = sqrt(rho_v/rho_h) (-); #aniso = #res. Defaults to ones.
- epermH, epermVarray_like, optional
Relative horizontal/vertical electric permittivities epsilon_h/epsilon_v (-); #epermH = #epermV = #res. If epermH is provided but not epermV, isotropic behaviour is assumed. Default is ones.
- mpermH, mpermVarray_like, optional
Relative horizontal/vertical magnetic permeabilities mu_h/mu_v (-); #mpermH = #mpermV = #res. If mpermH is provided but not mpermV, isotropic behaviour is assumed. Default is ones.
- verb{0, 1, 2, 3, 4}, optional
Level of verbosity, default is 2:
0: Print nothing.
1: Print warnings.
2: Print additional runtime and kernel calls
3: Print additional start/stop, condensed parameter information.
4: Print additional full parameter information
- Returns
- PJ0, PJ1array
Wavenumber-domain EM responses:
PJ0: Wavenumber-domain solution for the kernel with a Bessel function of the first kind of order zero.
PJ1: Wavenumber-domain solution for the kernel with a Bessel function of the first kind of order one.
See also
Examples
>>> import empymod >>> import numpy as np >>> src = [0, 0, 100] >>> rec = [5000, 0, 200] >>> depth = [0, 300, 1000, 1050] >>> res = [1e20, .3, 1, 50, 1] >>> freq = 1 >>> wavenr = np.logspace(-3.7, -3.6, 10) >>> PJ0, PJ1 = empymod.dipole_k(src, rec, depth, res, freq, wavenr, verb=0) >>> print(PJ0) [ -1.02638329e-08 +4.91531529e-09j -1.05289724e-08 +5.04222413e-09j -1.08009148e-08 +5.17238608e-09j -1.10798310e-08 +5.30588284e-09j -1.13658957e-08 +5.44279805e-09j -1.16592877e-08 +5.58321732e-09j -1.19601897e-08 +5.72722830e-09j -1.22687889e-08 +5.87492067e-09j -1.25852765e-08 +6.02638626e-09j -1.29098481e-08 +6.18171904e-09j] >>> print(PJ1) [ 1.79483705e-10 -6.59235332e-10j 1.88672497e-10 -6.93749344e-10j 1.98325814e-10 -7.30068377e-10j 2.08466693e-10 -7.68286748e-10j 2.19119282e-10 -8.08503709e-10j 2.30308887e-10 -8.50823701e-10j 2.42062030e-10 -8.95356636e-10j 2.54406501e-10 -9.42218177e-10j 2.67371420e-10 -9.91530051e-10j 2.80987292e-10 -1.04342036e-09j]