SimPEG Analytical Solutions


Does anyone have an analytical/numerical comparison of TDEM and FDEM modules for electric dipole? (Ward & Hohmann 2.40 and 2.50). I tried them but the analytical result was like on time and simpeg results were off time? Also, field values were really big.

Thank you.

Hi @Deniz, there is a python package called empymod:

where you can compute analytic solution for 1D layered earth.

Thank you, @sgkang09.

Hi @Deniz

I just stumbled upon this today. Did you found what you were looking for?

I actually have the figures around the equations you mentioned all reproduced, see here:

empymod has both the numerical and analytical solutions for all possible source-receiver combinations for fullspaces and halfspaces. However, the implementation does not follow Ward and Hohmann (which are derived from potentials, I think), but they follow Hunziker et al, 2015 (fullspaces) and Slob et al, 2010 (halfspaces), see for the full references.

I am happy to help with any questions if you have some.

(Thanks @sgkang09 for pointing to empymod!)

Hi @prisae

Thank you so much for the information. I was trying to compare analytical and numerical solutions for SimPEG. FDEM results were perfectly matched with Ward and Hohmann. However, whatever I do, TDEM result doesn’t match perfectly. Analytical and numerical results for TDEM are really close to each other but I was wondering if it is possible to get a better result.

Can you share your scripts? Then I could have a look, otherwise it is hard to tell…