# MuLTI-TEM
Multimodal Layered Transdimensional Inversion of Transient Electromagnetic data with Depth Constraints
(Can be used for all ground-based TEM survey designs, see parameter list below)
The MuLTI-TEM algorithm was writted by Siobhan Killingbeck and Phil Livermore in 2019, School of Earth and Environemnt, The University of Leeds. It is based on Matlab code written by Thomas Bodin. MuLTI-TEM was developed in Matlab version 2017a, therefore all Matlab codes supplied will work with this version or subsequent Matlab versions. This code can be run on Linux or Windows based platforms. This github repository includes:
1. 1D and 2D MuLTI-TEM Matlab scripts
2. Two matlab functions called in the main scripts, “thicknesses_and_priors” and “whichnuclei”,
3. 1D and 2D example dataset used to test MuLTI-TEM 1D and 2D codes,
4. Adapted Leroi forward modelling code (called LEROI_TEM.F90), written in Fortran 95, and it’s corresponding LEROI_TEM mex files used to be called in Matlab within the MuLTI-TEM script, for both Linux and Windows based platforms. The original Leroi code can be downloaded from the CSIRO/AMIRA P223F project at: http://p223suite.sourceforge.net/.
The mex interface has been tested on CentOS 7, and Windows 7, with Matlab 2017a. Further details and instructions on how to compile and use the LEROI_TEM mex file can be found in the corresponding folders and at https://github.com/cemac/MEX-LEROI_TEM. If working on a Linux based platform the LEROI_TEM mex file created for Linux specifically must be compiled in the Matlab working directory and then the LEROI_TEM function can be called in Matlab. If working on a Windows platform the LEROI_TEM mex file created for windows specifically does not need compiled and can just be called from the Matlab working directory.
To run the MuLTI-TEM Matlab codes, .m files, the corresponding platform based LEROI_TEM mex file must be in the active Matlab working directory along with the Matlab functions ‘thicknesses_and_priors’ and ‘whichnuclei’ and the input TEM data files.
For the 1D code: The input TEM data files are .mat files with the received voltages saved as a column vector variable called “data” with column 1 and 2 being the start and end times of the time gates in ms, column 3 being the time gate centre in ms and column 4 being the received voltages in nanoVolts/(Amp.meter^2). The weighting is determined from the variance of the measured data points, and can be a percentage of the data or inputted as a column vector.
For the 2D code: The input TEM data files are .mat files with the received voltages saved as a cell array variable called “data_2D”. Each cell is populated with a 1D sounding, ordered by acquisition along the line e.g., the first sounding acquired at the start of the line will be in the first cell and the last sounding at the end of the line will be in the last cell. Within each cell, the 1D soundings are column vectors with column 1 being the time gate centre in ms and column 2 being the received voltages in nanoVolts/(Amp.meter^2).
IMPORTANT NOTE ABOUT THE FORWARD MODELLING CODE:
In MuLTI-TEM, the Leroi code has been adapted, from its original code, for ground-based TEM survey designs only. The TEM survey parameters are defined in MuLTI-TEM as follows (we try to keep each parameter name the same as the original in Leroi where possible):
NCHNL = Number of time gates (normally either 20 or 30, depending on repetition frequency)
REFTYM = Time (in ms) from which TMS or TOPN & TCLS are measured. For example, this could be signal off-time or start of downward ramp.
OFFTYM = time (milliseconds) between end of one pulse and the start of the next pulse (of opposite sign) since a bipolar waveform is assumed (lambda/4). For systems which have a signal which is always on, OFFTIME = 0.
TXON = digitised time (in milliseconds) of each point in the waveform (set at 4 points). In most cases, TXON(1) = 0, TXON(2) = pulse on-time, TXON(3) = pulse off-time, TXON(4) = REFTYM where TXON(4) - TXON(3) = turn off time e.g., [0.0, 0.001, 1.0492, 1.05];
TXAMP = transmitter current in amps at time TXON(J), normally signal is normalised e.g., [0.0, 1.0, 1.0, 0.0];
TOPN = Start times (in ms) of receiver windows, (1x20);
TCLS = End times (in ms) of receiver windows, (1x20);
SXE = east coordinate of vertex I for loop position J (m), fixed at 4 vertices e.g., [5, -5, -5, 5]; Note the transmitter is fixed on the ground (Z=0) in this adapted Leroi code, for airborne data more parameters will need to be passed through the mex file to model.
SXN = north coordinate of vertex I for loop position J (m), fixed at 4 vertices e.g., [5, 5, -5, -5]; Note the transmitter is fixed on the ground (Z=0) in this adapted Leroi code, for airborne data more parameters will need to be passed through the mex file to model.
RXE = receiver easting (m);
RXN = receiver northing (m);
RXZ = receiver z (always be 0 for ground based TEM);
However, the original Leroi code can model almost any EM system and survey design, therefore this main code should be downloaded and adapted to the user’s specific EM survey method. A corresponding mex file will need to be created for the users specific adapted Leroi code, our mex file can be used as a guideline for this process.
DOI 10.5281/zenodo.3471638
# MuLTI-TEMp
The MuLTI-TEMp algorithm was adapted from the original MuLTI TEM code. The petrophysical (p) version was created in 2021 by Siobhan Killingbeck, Department of Geography & Environmental Management, The University of Waterloo funded by the W. Garfield Weston Foundation as part of the Canadian SEARCHArtic project.
MuLTI-TEMp additionally derives salinity of subsurface fluids. It uses Archie’s mixing law and Fofonoff and Millard’s 1983 internationally recognized conductivity to salinity conversion method, outputting salinity in terms of its practical salinity from the Practical Salinity Scale 1978, (Perkin and Lewis, 1980). Archie's law and Fofonoff and Millard’s 1983 conversion are applied to each resistivity-depth model within the ensemble of the posterior distribution to output an additional PDF of practical salinity, provided that a prior distribution of porosity can be reliably estimated, for example, from nearby borehole measurements or detailed seismic velocity analysis. A prior distribution of pressure and temperature is also required for the conductivity to salinity conversion; this can be obtained from borehole measurements or modelling, for example, using a 1D steady-state advection-diffusion model. The code is set up to use a Gaussian distribution to randomly sample each petrophysical parameter (temperature, pressure, cementation factor and porosity) from a prior mean and standard deviation, however this can be changed to a uniform or skewed distribution with depth, depending on the observed data. Note the mixing law can also be easily changed if required.
ADDITIONAL FUNCTIONS TO BE USED WITH MuLTI-TEMp
* "salinity.m" matlab script written by Edward T Peltzer, MBARI in 2007
ADDITIONAL INPUTS TO BE USED WITH MuLTI-TEMp
* porosity = estimate of porosity model from known data, if you want to calcuate pore fluid salinity.
* sigma_porosity = estimate of porosity uncertainty
* Rfluid = estimate of pore fluid resistivity model from known data e.g., borehole or lake measurements, if you want to calcuate porosity.
* sigma_Rfluid = estimate of pore fluid uncertainty
* cementation = m (cementation factor) Archie's law
* sigma_cementation = estimate of uncertainty in m
* temp = temperature degrees C
* sigma_temp = temperature uncertainty degrees C
* pressure = pressure in dbars
* sigma_pressure = uncertainty in dbars
MULTIPLE BASE FREQUENCIES INVERSION
* MuLTI-TEMp has also been developed to invert for multiple base frequencies. The version: MuLTI-TEMp_III is currently written to invert for THREE base frequencies. In this version, the likelihood of a proposed mode
没有合适的资源?快使用搜索试试~ 我知道了~
具有深度约束 的瞬态电磁数据 的多模态分层跨维反演_Fortran_代码_下载
共41个文件
mat:15个
m:14个
md:3个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 19 浏览量
2022-07-04
14:27:25
上传
评论
收藏 3.79MB ZIP 举报
温馨提示
具有深度约束的瞬态电磁数据的多模态分层跨维反演 (可用于所有地面 TEM 测量设计,请参见下面的参数列表) MuLTI-TEM 算法由利兹大学地球与环境学院的 Siobhan Killingbeck 和 Phil Livermore 于 2019 年编写。它基于 Thomas Bodin 编写的 Matlab 代码。MuLTI-TEM 是在 Matlab 版本 2017a 中开发的,因此提供的所有 Matlab 代码都将适用于该版本或后续的 Matlab 版本。此代码可以在基于 Linux 或 Windows 的平台上运行。这个 github 存储库包括: 1D 和 2D MuLTI-TEM Matlab 脚本 在主脚本中调用了两个 matlab 函数,“thicknesses_and_priors”和“whichnuclei”, 用于测试 MuLTI-TEM 1D 和 2D 代码的 1D 和 2D 示例数据集, 改编的 Leroi 正向建模代码(称为 LEROI_TEM.F90) 更多详情、使用方法,请下载后阅读README.md文件
资源推荐
资源详情
资源评论
收起资源包目录
MuLTI-TEM-master.zip (41个子文件)
MuLTI-TEM-master
functions
salinity.m 2KB
thicknesses_and_priors.m 4KB
whichnuclei.m 1KB
LEROI_TEM_mex_Linux
LEROI_TEM.mexa64 448KB
test.m 2KB
Leroi.cfl 1KB
Leroi_FM_TEM_subs.f90 367KB
test_1lyr.m 1KB
Parameter_list.txt 2KB
LEROI_TEM.F90 394KB
README.md 4KB
tools
colourbars
cmap_blue.mat 682B
cmap_jet.mat 313B
salt_cmap.mat 393B
error_cmap.mat 2KB
Plotting_1D_responses_along_2D_line.m 3KB
MuLTI_TEMp
MuLTI_TEMp_1D.m 38KB
MuLTI_TEMp_2D.m 45KB
example_input_data
cmap_jet.mat 313B
1D_input_III.mat 2KB
2D_input_data_III.mat 5KB
Rfluid.mat 231B
1D_input.mat 690B
synthetic_data_2D_III.m 16KB
Rfluid2D.mat 759B
porosity.mat 225B
synthetic_data_2D.m 5KB
2D_input_data.mat 2KB
porosity2D.mat 713B
MuLTI_TEMp_III_1D.m 47KB
MuLTI_TEMp_III_2D.m 55KB
LICENSE 34KB
MuLTI_TEM
example_input_data
MuLTI_TEM_2D_input_data.mat 4KB
MuLTI_TEM_1D_input_data.mat 649B
MuLTI_TEM_1D.m 28KB
MuLTI_TEM_2D.m 32KB
LEROI_TEM mex file Windows
LEROI_TEM.mexw64 730KB
WindowsCompilationNotes.md 659B
README.md 9KB
Publications
Fofonoff_and_Millard_1983.pdf 2.02MB
Killingbeck_et_al_2020_solid_earth.pdf 1.23MB
共 41 条
- 1
资源评论
快撑死的鱼
- 粉丝: 1w+
- 资源: 9154
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功