We are trying to model and invert the data set and tried to solve but not sure about it. need some better
solution for TDEM 1D inversion.
In attachment, a data file on which we are working right now as a test to better manage the code for both forward modeling and standard inversion. It is a 1D sounding . This data was collected along an ERT profile. The first test we want to do is to be able to model the data that should be recorded if we consider the ERT inversion as the true model, simplified as a 4-layer model (0-7 m : 2000 Ohm.m, 7-20 m: 120 Ohm.m, 20-45 m: 300 Ohm.m, 45- : 15 Ohm.m). We used a 25 x 25 square loop as transmitter/receiver. The theoretical input signal from the device is a square wave (thus a step-like shape) with ratio pulse/pause = 3/1, the dead-time is 4 µs, the intensity of the current was 3.4 A. but the true signal during turning off has a ramp. We would like to see the influence of this on the data and thus model both. We have seen that SIMPEG can model this type of signal, but are not sure how to proceed exactly.
We guess that if we are able to model this, this should get us started.
In attachment, there are some tentative code based on SimPEG published examples. In the first one, we just generate synthetic data that we invert. In the second, we are trying to invert the collected real field data. https://docs.simpeg.xyz/content/examples/09-tdem/plot_inversion_1D.html#sphx-glr-content-examples-09-tdem-plot-inversion-1d-py and https://docs.simpeg.xyz/content/examples/00_published/plot_heagyetal2017_cyl_inversions.html#sphx-glr-content-examples-00-published-plot-heagyetal2017-cyl-inversions-py
We have a few questions related to those:
- We are not sure how we can easily generate a grid which would fit specific boundaries as expected in our case. In the current example, we define 4 homogeneous layers but whose limits do not coincide exactly with the grid.
- There are a series of available source in SIMPEG, but we are not sure how to simulate correctly our source, which is a square loop 25 m x 25 m (for both transmitter and receiver). Here we used a circular loop of similar area. Similarly, we are not sure how to specify the receiver loop size.
- Our original data (see 1107_4.tem file) are expressed in V/A as recorded by the device (recorded voltage divided by input current intensity) as most TEM device do. However, it seems that SIMPEG is generating the vertical component of the magnetic field as calculated data. Therefore we tried to convert those back to V and Magnetic Field. However, it feels like this can probably be done within SIMPEG ? Is it possible to specify the type of data we want as input/output of the forward model or input of the inversion? Or should we always convert those as we tried here ?
In summary, we could use some help to get started with the proper use of the code.
Thanks in advance