Aurora Geosciences (AGL) based out of Whitehorse would like some assistance on inverting a public domain data set of Extremely Low Frequency (ELF) and gravity surveys .
The ELF surveys were conducted by an AGL ground geophysical crew. The data collected is Z-Axis Tipper Electromagnetics. We are not referencing an E- Field. We are just reading the inphase and quadrature tilt angle of the magnetic field at frequencies from 11Hz to 1440Hz.
There is no controlled source involved, and there are no electric dipoles being measured. It is a passive system.There is no base station in use either.
We can provide all technical specs, the data is similar to ZTEM tipper data
I have attached a written confirmation that the data set from the Yukon Geological Survey has released the data set to AGL and placed into the public domain.
A Gravity survey , conducted for the Yukon Geological Survey (YGS), was used to assist in the imaging of the Denali fault for Geothermal exploration.The survey was conducted in a 10 km X 7 km area located within NTS map sheet 115 G/4 with a station separation of 250 m and a line separation of 1000 m.In total 295 data points were collected in the 2019 program between March 15, 2019 and April 3, 2019.
An ELF survey consisting of 329 data points was conducted from June 12th to June 29th, 2019.
The survey consisted of eleven 7000 m lines with stations every 250 m. The in-phase & quadrature maps for the N/E component profile lines, tipper vectors, & divergence plots for all measured frequencies (11, 22, 45, 90, 180, 360, 720, & 1440 Hz) were plotted and modeled.
Logistical and technical Geophysical reports for the above surveys have been written by AGL and submitted to YGS during 2019.
I have attached a poster "Geothermal Studies by the Yukon Geological Survey " as a pdf. This poster was presented at the 2019 Yukon Geoscience Forum. Geothermal YGF 2019_lores.pdf (3.8 MB)
I have included the raw data from Line 5000 at the bottom of this post, and I have attached a pdf the contains a description of the columns
YGS-20191209 Description of the colums in the ELF observations.pdf (181.0 KB) Raw Data contains a CSV file of the observations and a pdf containing a description of the columns of the CSV.
Since I currently only have permissions to attache image files and pdfs, I attempted to attach a pdf of the raw data. The PDF was not able to preserve some formatting when i copied the contents back into csv.
For this reason, I have attached the raw data here. I was able to copy and paste the text below into a MS spreadsheet by using a white space deliminator. So I hope that you will be able to do the same.
Hi Alex, you are in luck, Z-Axis tipper is good to go out of the box with SimPEG. I have successfully inverted a few synthetic models using tipper only. The biggest initial hurdle will be writing an IO for your data file so that a SimPEG Survey class can be generated. I would suggest taking a look at the 3D MT forward model example. It includes both tipper setup and regular MT setup. As long as you can get your data into a SimPEG survey we should be able to guide you through to inversion. Let me know how it goes and we can discuss further.
Hi @AlexJ, I have been working away at NSEM problems and I have made a general data loader for databases like the one you provided here. This script sets up the SimPEG survey object that can be used in the inversion. It would still require a little modification in the code but most of the heavy lifting is done on that part. If you were interested I could share the code for you to develop further or at least I can guide you through it? Let me know what you think!
Perfect! Give me a few days to put together the script with a bit more detailed comments to help you along the way. I’m a little busy right now but I will get back to you by the end of the week.
@AlexJ, apologies for the delay, work projects have not let up all that much lately. I have attached my module for IO of NSEM data. I didn’t get to commenting as much as I’d like to help guide you but feel free to contact me on here and we can work this out. I have broken things up into a “Patch” class that hold “stations”. Also included is a tool named readDat. This is the function you will be modifying. What this does is determine the stations of the survey and then for each station assign the frequency data and whichever else you want to add. If there are some parameters that you would like include, you can modify the “mtStation” class to initiate them. The file format I have this setup for currently is a little different than your’s but not far off. Start by having a look and I will help the best I can. AMTtools.py (13.3 KB)