23 Aug 2023
Weekly SimPEG meeting from August 23rd.
A discussion on solvers and bringing back MUMPS.
A discussion on solvers and bringing back MUMPS.
- 9 participants
- 1:00 hours
16 Aug 2023
Weekly SimPEG meeting from August 16th.
A discussion on TDEM implementations.
A discussion on TDEM implementations.
- 6 participants
- 41 minutes
11 Jan 2023
Recording of the weekly meeting from January 11th.
Some new faces, and discussions on upcoming TEM developments for the codebase.
Some new faces, and discussions on upcoming TEM developments for the codebase.
- 8 participants
- 30 minutes
23 Jun 2022
The salinization map of the region of Flanders, Belgium shows the depth of the interface between fresh and salt groundwater in the coastal and polder area. It serves as an exploratory tool to examine the potential of groundwater projects that improve freshwater availability in the shallow subsurface. Flanders environment agency published an updated map in 2019, based on airborne time-domain electromagnetic induction data. The result of the inversion is, however, overly smooth, which potentially conceals interesting features.
Via an inverse problem, the electromagnetic induction data can be mapped onto a conductivity profile, which serves as a proxy for salinity via petrophysical laws. The inverse problem is ill-posed and regularization improves the stability of the inversion. Based on Occam’s razor principle, a smoothing constraint is typically used with a very large number of thin layers. However, the salinity profiles in the Belgian coastal plains are sometimes sharp, impeding the correct estimation of the fresh-saltwater interface. In practice, the real underground might be either blocky or smooth, or somewhere in between. Standard constraints are thus not always appropriate.
With our novel wavelet-based inversion scheme, the original time-domain AEM data can be re-interpreted in a flexible fashion, meaning that we can easily generate an ensemble of inversion models with different types of Occam's razor minimum-structure. The flexibility is due to the wavelet basis. In simple terms, a wavelet function can be seen as a building block and a simple model is one that can be built with few building blocks of various sizes. Our proposed inversion scheme adds a regularization term that limits the number of building blocks in the wavelet-domain to make sure only the necessary complexity is retrieved. The scale-dependency makes the use of small blocks more expensive, as it corresponds to adding detail. The scheme is tuned by only one additional parameter (which determines the wavelet basis function) and can recover blocky, intermediate and smooth structures. It is also capable to recover high amplitude anomalies in combination with globally smooth profiles, a common problem for smooth inversion, and an essential feature to accurately predict the salinity.
We first demonstrated this alternative inversion scheme on 1D FDEM data and now extend it to 2D AEM data a in a saltwater intrusion context . The flexibility of the method allows for choosing the appropriate sharpness for each orientation. In Figure 1, an inversion model is shown with a relatively sharp wavelet basis (db3) in the vertical orientation, while a smoother wavelet basis (db8) is used along with the lateral orientation. Imposing the appropriate sharpness on the inversion model is crucial to obtaining a reliable estimation of the fresh-saltwater interface.
The code behind the 1D regularization has been adjusted to the SimPEG framework and is open to everyone . This allows to use this flexible regularization term for many types of inversion problems and promotes geoscience reproducibility!
Bio
Wouter Deleersnyder is a passionate Ph.D. candidate at the physics department of KU Leuven and Geology department of Ghent University. With his technical skills from his physics education (Master in Physics, KU Leuven, Belgium in 2019), he wants to tackle the problems of today. Stimulated by the importance of water for man and nature, he came into contact with groundwater and geophysics. He is working on new methods to get a better picture of what is under the ground. He is specializing in electromagnetic induction surveys, forward modelling, inversion, and uncertainty quantification.
Via an inverse problem, the electromagnetic induction data can be mapped onto a conductivity profile, which serves as a proxy for salinity via petrophysical laws. The inverse problem is ill-posed and regularization improves the stability of the inversion. Based on Occam’s razor principle, a smoothing constraint is typically used with a very large number of thin layers. However, the salinity profiles in the Belgian coastal plains are sometimes sharp, impeding the correct estimation of the fresh-saltwater interface. In practice, the real underground might be either blocky or smooth, or somewhere in between. Standard constraints are thus not always appropriate.
With our novel wavelet-based inversion scheme, the original time-domain AEM data can be re-interpreted in a flexible fashion, meaning that we can easily generate an ensemble of inversion models with different types of Occam's razor minimum-structure. The flexibility is due to the wavelet basis. In simple terms, a wavelet function can be seen as a building block and a simple model is one that can be built with few building blocks of various sizes. Our proposed inversion scheme adds a regularization term that limits the number of building blocks in the wavelet-domain to make sure only the necessary complexity is retrieved. The scale-dependency makes the use of small blocks more expensive, as it corresponds to adding detail. The scheme is tuned by only one additional parameter (which determines the wavelet basis function) and can recover blocky, intermediate and smooth structures. It is also capable to recover high amplitude anomalies in combination with globally smooth profiles, a common problem for smooth inversion, and an essential feature to accurately predict the salinity.
We first demonstrated this alternative inversion scheme on 1D FDEM data and now extend it to 2D AEM data a in a saltwater intrusion context . The flexibility of the method allows for choosing the appropriate sharpness for each orientation. In Figure 1, an inversion model is shown with a relatively sharp wavelet basis (db3) in the vertical orientation, while a smoother wavelet basis (db8) is used along with the lateral orientation. Imposing the appropriate sharpness on the inversion model is crucial to obtaining a reliable estimation of the fresh-saltwater interface.
The code behind the 1D regularization has been adjusted to the SimPEG framework and is open to everyone . This allows to use this flexible regularization term for many types of inversion problems and promotes geoscience reproducibility!
Bio
Wouter Deleersnyder is a passionate Ph.D. candidate at the physics department of KU Leuven and Geology department of Ghent University. With his technical skills from his physics education (Master in Physics, KU Leuven, Belgium in 2019), he wants to tackle the problems of today. Stimulated by the importance of water for man and nature, he came into contact with groundwater and geophysics. He is working on new methods to get a better picture of what is under the ground. He is specializing in electromagnetic induction surveys, forward modelling, inversion, and uncertainty quantification.
- 5 participants
- 45 minutes
27 May 2022
With modern instrumentation time-series data is readily available heightening the level of quality control. DCIP acquisition for example, resistivities, chargeabilities and decay curves are no longer the base of validating data points. The building blocks of the calculated properties can now be analyzed to uphold a better integrity of the data itself. The challenge now is analyzing numerous time-series manually. Each point has transmit and receiver data. Particularly with distributed systems, an overwhelming amount data is available. Classically used statistical methods can be employed but can cost a lot in time. Next best would be to simulate the geophysicist looking at the data. Borrowing techniques from computer vision, convolution neural networks can be trained to “visually” inspect time-series. With most of the heavy lifting of time-series analysis alleviated, the geophysicist can have more time and information to make better informed decisions. Though DCIP time-series is the primary focus here, the techniques presented can certainly extend into time-domain controlled source electromagnetics.
Bio
John Kuttai is a senior geophysicst for DIAS Geophysical’s research and development team (also, soon to be UBC Graduate student!!). While primarily servicing the mineral exploration sector, his contributions maintain the processing and inversion needs for DC resistivity, induced polarization, magnetic gradiometry, natural and controlled source frequency domain methods. Through signal processing to inversion, a complete comprehensive work flow to manage big data and inversion for drill target identification is the main focus of his work.
Bio
John Kuttai is a senior geophysicst for DIAS Geophysical’s research and development team (also, soon to be UBC Graduate student!!). While primarily servicing the mineral exploration sector, his contributions maintain the processing and inversion needs for DC resistivity, induced polarization, magnetic gradiometry, natural and controlled source frequency domain methods. Through signal processing to inversion, a complete comprehensive work flow to manage big data and inversion for drill target identification is the main focus of his work.
- 8 participants
- 54 minutes
21 Apr 2022
Abstract
Scientific communication today is designed around print documents and paywalled access to content. Over the last decade, the open-science movement has accelerated the use of pre-print services and data archives that are vastly improving the accessibility of scientific content. However, these systems are not designed for communicating modern scientific outputs, which encompasses much more than a “paper-centric model of the scholarly literature”.
In this presentation, we will give a background on challenges with today’s tools for research communication & collaboration, and present a vision for the future that follows the FORCE11 recommendations (Bourne et al., 2012). Specifically: (1) rethink the unit and form of scholarly publication; (2) develop tools and technologies to better support the scholarly lifecycle; and (3) add data, software, and workflows as first-class research objects.
We will discuss these recommendations in the context of (a) the ExecutableBooks community and new markup languages for scientific communication; (b) a collaborative writing tool called Curvenote that integrates with Jupyter Notebooks, and (c) new publishing tools that support networked scientific communication throughout the scholarly lifecycle.
Excerpts from FORCE11:
“We see a future in which scientific information and scholarly communication more generally become part of a global, universal and explicit network of knowledge; where every claim, hypothesis, argument—every significant element of the discourse—can be explicitly represented, along with supporting data, software, workflows, multimedia, external commentary, and information about provenance.”“This vision moves away from the paper-centric model of the scholarly literature, towards a more distributed network-centric model” that “vastly improves knowledge transfer and [has a] far wider impact.” (Bourne et al., 2012)
Bio
Rowan completed his PhD in Geophysics at the University of British Columbia, Canada, in 2017. His PhD laid out a framework for geophysical inversions (SimPEG!) and applied it to large-scale, vadose-zone groundwater flow. Currently, Rowan works at Curvenote, a company that sits at the intersection of scientific collaboration, publishing, and technology and provides tools for writing and collaborating on technical documents. Rowan founded Curvenote in 2020, with the mission to help get scientific communication out of PDFs and onto the web — the way we communicate scientific knowledge should evolve past the status quo of print-based publishing and all the limitations of paper. He is a core team member of the ExecutableBooks community, which supports MyST Markdown and JupyterBook.
Scientific communication today is designed around print documents and paywalled access to content. Over the last decade, the open-science movement has accelerated the use of pre-print services and data archives that are vastly improving the accessibility of scientific content. However, these systems are not designed for communicating modern scientific outputs, which encompasses much more than a “paper-centric model of the scholarly literature”.
In this presentation, we will give a background on challenges with today’s tools for research communication & collaboration, and present a vision for the future that follows the FORCE11 recommendations (Bourne et al., 2012). Specifically: (1) rethink the unit and form of scholarly publication; (2) develop tools and technologies to better support the scholarly lifecycle; and (3) add data, software, and workflows as first-class research objects.
We will discuss these recommendations in the context of (a) the ExecutableBooks community and new markup languages for scientific communication; (b) a collaborative writing tool called Curvenote that integrates with Jupyter Notebooks, and (c) new publishing tools that support networked scientific communication throughout the scholarly lifecycle.
Excerpts from FORCE11:
“We see a future in which scientific information and scholarly communication more generally become part of a global, universal and explicit network of knowledge; where every claim, hypothesis, argument—every significant element of the discourse—can be explicitly represented, along with supporting data, software, workflows, multimedia, external commentary, and information about provenance.”“This vision moves away from the paper-centric model of the scholarly literature, towards a more distributed network-centric model” that “vastly improves knowledge transfer and [has a] far wider impact.” (Bourne et al., 2012)
Bio
Rowan completed his PhD in Geophysics at the University of British Columbia, Canada, in 2017. His PhD laid out a framework for geophysical inversions (SimPEG!) and applied it to large-scale, vadose-zone groundwater flow. Currently, Rowan works at Curvenote, a company that sits at the intersection of scientific collaboration, publishing, and technology and provides tools for writing and collaborating on technical documents. Rowan founded Curvenote in 2020, with the mission to help get scientific communication out of PDFs and onto the web — the way we communicate scientific knowledge should evolve past the status quo of print-based publishing and all the limitations of paper. He is a core team member of the ExecutableBooks community, which supports MyST Markdown and JupyterBook.
- 6 participants
- 54 minutes
24 Mar 2022
Abstract
The use of airborne electromagnetic (AEM) data for geoscience applications is rapidly increasing. For instance, in California USA, there is an ongoing AEM project led by the California Department of Water Resources (CDWR), which plans to map out most of the Central Valley of California and some water basins in California. All acquired AEM data and resulting interpretation of the data, which are resistivity models of the subsurface will be publicly available. There are many big initiatives like this project throughout the world. Therefore, it is critical to understand how the resulting resistivity models are obtained from the acquired AEM data, and further equipped with an ability to download and explore the available resistivity data. In this talk, I will first introduce how a resistivity model is obtained from the AEM data then introduce open-source tools that can be used to explore this resistivity model.
Bio
Seogi completed his PhD in Geophysics at the University of British Columbia, Canada, in 2018. His thesis work focused on computational electromagnetics and its application to mining problems. Currently, he is a Postdoctoral Researcher in the Geophysics Department at Stanford. His research focus is on advancing the use of airborne electromagnetic and other remote sensing methods for groundwater management and groundwater science. He continues to contribute to the development of open source software, SimPEG, and educational resources, GeoSci.xyz, for geophysics.
The use of airborne electromagnetic (AEM) data for geoscience applications is rapidly increasing. For instance, in California USA, there is an ongoing AEM project led by the California Department of Water Resources (CDWR), which plans to map out most of the Central Valley of California and some water basins in California. All acquired AEM data and resulting interpretation of the data, which are resistivity models of the subsurface will be publicly available. There are many big initiatives like this project throughout the world. Therefore, it is critical to understand how the resulting resistivity models are obtained from the acquired AEM data, and further equipped with an ability to download and explore the available resistivity data. In this talk, I will first introduce how a resistivity model is obtained from the AEM data then introduce open-source tools that can be used to explore this resistivity model.
Bio
Seogi completed his PhD in Geophysics at the University of British Columbia, Canada, in 2018. His thesis work focused on computational electromagnetics and its application to mining problems. Currently, he is a Postdoctoral Researcher in the Geophysics Department at Stanford. His research focus is on advancing the use of airborne electromagnetic and other remote sensing methods for groundwater management and groundwater science. He continues to contribute to the development of open source software, SimPEG, and educational resources, GeoSci.xyz, for geophysics.
- 7 participants
- 55 minutes
18 Feb 2022
Randomize-then-optimize: nonlinear uncertainty quantification for regularized inversion
Much of our understanding of the Earth's subsurface comes from physical models produced by inverting geophysical observations made at the surface. Inversion is non-unique and nonlinear, however, meaning there is significant uncertainty in the inverted model parameters. The standard inversion method in geophysics remains regularized inversion despite the fact that it produces single model estimates without a meaningful way to quantify the uncertainty. In this talk I will present 'randomize-then-optimize' (RTO), an uncertainty quantification (UQ) strategy for regularized inversion. This method reinterprets regularized inversion in a Bayesian context, turning these familiar algorithms into Bayesian samplers capable of producing model uncertainty efficiently, even for large geophysical problems. I will discuss the basic theory behind RTO, describe our extension of it to hierarchically sample the regularization strength (which we call RTO-TKO), and show results on field data examples from electromagnetic geophysics.
Bio
Daniel is the John W Miles postdoctoral scholar in computational and theoretical geophysics at the Scripps Institution of Oceanography, UC San Diego. His research interests are primarily focused in two fields: (1) the lithosphere-asthenosphere system, and the role of fluids in particular, and their relation to plate tectonics; and (2) the development of algorithms capable of quantifying uncertainty in inverted subsurface models inferred from geophysical data. He studied computational mathematics at Stanford University, receiving a Masters in 2015. He did his doctoral studies in geophysics with professor Kerry Key at Columbia University, receiving a PhD in 2020.
Much of our understanding of the Earth's subsurface comes from physical models produced by inverting geophysical observations made at the surface. Inversion is non-unique and nonlinear, however, meaning there is significant uncertainty in the inverted model parameters. The standard inversion method in geophysics remains regularized inversion despite the fact that it produces single model estimates without a meaningful way to quantify the uncertainty. In this talk I will present 'randomize-then-optimize' (RTO), an uncertainty quantification (UQ) strategy for regularized inversion. This method reinterprets regularized inversion in a Bayesian context, turning these familiar algorithms into Bayesian samplers capable of producing model uncertainty efficiently, even for large geophysical problems. I will discuss the basic theory behind RTO, describe our extension of it to hierarchically sample the regularization strength (which we call RTO-TKO), and show results on field data examples from electromagnetic geophysics.
Bio
Daniel is the John W Miles postdoctoral scholar in computational and theoretical geophysics at the Scripps Institution of Oceanography, UC San Diego. His research interests are primarily focused in two fields: (1) the lithosphere-asthenosphere system, and the role of fluids in particular, and their relation to plate tectonics; and (2) the development of algorithms capable of quantifying uncertainty in inverted subsurface models inferred from geophysical data. He studied computational mathematics at Stanford University, receiving a Masters in 2015. He did his doctoral studies in geophysics with professor Kerry Key at Columbia University, receiving a PhD in 2020.
- 9 participants
- 1:20 hours
24 Jan 2022
With distributed array 3D DCIP aquisition and DIAS’s Common Voltage Reference method, the number of possible dipole collection grows signficantly compared to traditional 2D. Data sets can easily grow into a million or more data points. Being distributed, aquisition over extreme topography is done frequently. This often leads to large mesh cell counts for the fine discretization of the topography. The combination of numerous transmits and large cell counts heavily consume computational resources and take significant time for inversion. To improve performance, parallelizing the inversion is required. Configuring a cluster with the right combination of hardware and software can provide significant gains. SimPEG code is then modified to divide the inversion problem into smaller simulations which are assigned to each cluster node. Here I explore inverting datasets from small to large data and cell counts. The largest requiring more than 1TB of RAM, which is larger than most workstation’s capacity. The results and the work are based on the SimPEG framework parallelized for effective cluster usage.
Bio
John Kuttai is a senior geophysicst for DIAS Geophysical’s research and development team. While primarily servicing the mineral exploration sector, his contributions maintain the processing and inversion needs for DC resistivity, induced polarization, magnetic gradiometry, natural and controlled source frequency domain methods. Through signal processing to inversion, a complete comprehensive work flow to manage big data and inversion for drill target identification is the main focus of his work.
Bio
John Kuttai is a senior geophysicst for DIAS Geophysical’s research and development team. While primarily servicing the mineral exploration sector, his contributions maintain the processing and inversion needs for DC resistivity, induced polarization, magnetic gradiometry, natural and controlled source frequency domain methods. Through signal processing to inversion, a complete comprehensive work flow to manage big data and inversion for drill target identification is the main focus of his work.
- 7 participants
- 49 minutes
12 Aug 2021
August 2021 SimPEG Seminar. Implementing geological rules within geophysical inversion: A PGI perspective
Inferring geologically meaningful information from a geophysical inversion is a challenging task. Moreover, prior knowledge about the petrophysical contrasts or the relationships between various geological units can also prove difficult to translate into quantitative input for the inversion. In previous works, we developed a Petrophysically and Geologically guided Inversion (PGI) framework that enables desired petrophysical characteristics to be reproduced. This information is encoded into the objective function's smallness through a Gaussian Mixture Model (GMM). The resulting discrete geological representation of the subsurface thus fits both the geophysical and petrophysical information. The way we included geological information was limited to favouring the occurrence of chosen rock units in user-defined areas on a cell-by-cell approach. Transferring geological information from one area to another, such as an expected stratigraphy, was not easily done. Moreover, structural information (dip orientation, etc.), which by definition depends on multiple cells at once, was left to the objective function's smoothness, which acts on the physical property models rather than on the geological representation itself. We improve upon the existing PGI framework to make the inversion result geologically realistic by including geological rules as part of the process that builds the geological representation throughout the inversion's iterations. For this purpose, we incorporate image segmentation tools using Markov Random Field (MRF) as part of the PGI framework. The final recovered model fits geophysical and petrophysical information while reproducing geological characteristics, thus providing a more faithful and informed representation of the underground.
Thibaut Astic received a Ph.D. (2020) from the Geophysical Inversion Facility (GIF) at the University of British Columbia (Vancouver, BC, Canada), where he is currently working as a postdoctoral researcher. His research focuses on joint inversion coupled by petrophysical and geological information and the development of open-source tools for the geosciences community, mostly through the Python package SimPEG. Before his Ph.D., Thibaut worked in geological mapping and geophysical data acquisition and processing at the Quebec Department of Natural Resources and various geophysics companies.
Inferring geologically meaningful information from a geophysical inversion is a challenging task. Moreover, prior knowledge about the petrophysical contrasts or the relationships between various geological units can also prove difficult to translate into quantitative input for the inversion. In previous works, we developed a Petrophysically and Geologically guided Inversion (PGI) framework that enables desired petrophysical characteristics to be reproduced. This information is encoded into the objective function's smallness through a Gaussian Mixture Model (GMM). The resulting discrete geological representation of the subsurface thus fits both the geophysical and petrophysical information. The way we included geological information was limited to favouring the occurrence of chosen rock units in user-defined areas on a cell-by-cell approach. Transferring geological information from one area to another, such as an expected stratigraphy, was not easily done. Moreover, structural information (dip orientation, etc.), which by definition depends on multiple cells at once, was left to the objective function's smoothness, which acts on the physical property models rather than on the geological representation itself. We improve upon the existing PGI framework to make the inversion result geologically realistic by including geological rules as part of the process that builds the geological representation throughout the inversion's iterations. For this purpose, we incorporate image segmentation tools using Markov Random Field (MRF) as part of the PGI framework. The final recovered model fits geophysical and petrophysical information while reproducing geological characteristics, thus providing a more faithful and informed representation of the underground.
Thibaut Astic received a Ph.D. (2020) from the Geophysical Inversion Facility (GIF) at the University of British Columbia (Vancouver, BC, Canada), where he is currently working as a postdoctoral researcher. His research focuses on joint inversion coupled by petrophysical and geological information and the development of open-source tools for the geosciences community, mostly through the Python package SimPEG. Before his Ph.D., Thibaut worked in geological mapping and geophysical data acquisition and processing at the Quebec Department of Natural Resources and various geophysics companies.
- 7 participants
- 1:14 hours
15 Jul 2021
July 2021 SimPEG Seminar. On recovering changes of water head from satellite ground deformation data
Population growth and climate change in 21 century increase the demand for groundwater threatening the sustainability of groundwater resources. Therefore, there is an urgent need to understand the groundwater system. Major drivers of the system are the spatial and temporal changes of the water head, and therefore obtaining these changes is an essential task. Although head data can be measured from water wells (e.g. monitoring, irrigation), the spatial coverage and density of the wells are often poor resulting in a large data gap between the wells. The ground deformation signals, which can be measured from satellites (e.g. Sentinel), contain information about the changes of the head providing great potential to fill in this data gap. Working on the ground deformation data and head data acquired at the Central Valley of California, which is one of the most productive farmland in the world resulting in a large amount of water demand, I present a developed methodology that can 1) simulate the deformation signals from given head data and other aquifer parameters and 2) invert the observed deformation signals to estimate head data and other aquifer parameters.
Seogi completed his PhD in Geophysics at the University of British Columbia, Canada, in 2018. His thesis work focused on computational electromagnetics and its application to mining problems. Currently, he is a Postdoctoral Researcher in the Geophysics Department at Stanford. His research focus is on advancing the use of airborne electromagnetic and other remote sensing methods for groundwater management and groundwater science. He continues to contribute to the development of open source software, SimPEG, and educational resources, GeoSci.xyz, for geophysics.
Population growth and climate change in 21 century increase the demand for groundwater threatening the sustainability of groundwater resources. Therefore, there is an urgent need to understand the groundwater system. Major drivers of the system are the spatial and temporal changes of the water head, and therefore obtaining these changes is an essential task. Although head data can be measured from water wells (e.g. monitoring, irrigation), the spatial coverage and density of the wells are often poor resulting in a large data gap between the wells. The ground deformation signals, which can be measured from satellites (e.g. Sentinel), contain information about the changes of the head providing great potential to fill in this data gap. Working on the ground deformation data and head data acquired at the Central Valley of California, which is one of the most productive farmland in the world resulting in a large amount of water demand, I present a developed methodology that can 1) simulate the deformation signals from given head data and other aquifer parameters and 2) invert the observed deformation signals to estimate head data and other aquifer parameters.
Seogi completed his PhD in Geophysics at the University of British Columbia, Canada, in 2018. His thesis work focused on computational electromagnetics and its application to mining problems. Currently, he is a Postdoctoral Researcher in the Geophysics Department at Stanford. His research focus is on advancing the use of airborne electromagnetic and other remote sensing methods for groundwater management and groundwater science. He continues to contribute to the development of open source software, SimPEG, and educational resources, GeoSci.xyz, for geophysics.
- 9 participants
- 1:05 hours
24 Feb 2021
Weekly SimPEG Meeting from February 24th, 2021.
Includes, discussions on next steps for merging the EM1D subpackge into SimPEG, and some interpolation work by Seogi Kang.
Includes, discussions on next steps for merging the EM1D subpackge into SimPEG, and some interpolation work by Seogi Kang.
- 10 participants
- 50 minutes
17 Feb 2021
Weekly SimPEG Meeting from February 17th, 2021
Show preview of new boundary condition implementation coming to discretize!
Show preview of new boundary condition implementation coming to discretize!
- 10 participants
- 1:13 hours
26 Aug 2020
Thibaut Astic presents the preliminary version of his Ph.D. defence: "A framework for joint petrophysically and geologically guided geophysical inversion".
Presentation's slides: http://bit.ly/astic_phd_defence; Dissertation: http://bit.ly/astic_phd_thesis
Presentation's slides: http://bit.ly/astic_phd_defence; Dissertation: http://bit.ly/astic_phd_thesis
- 6 participants
- 1:03 hours
13 Aug 2019
SimPEG meeting on August 13 - updates on including dask and Zarr in the DC code
- 5 participants
- 32 minutes
6 Aug 2019
- Discussion on a trajectory for generating examples and answering user-questions on SimPEG
- John presents on using SimPEG to invert DCIP data / Mag Gradiometry data and progress on implementing new receiver types in the natural source electromagnetic code
- John presents on using SimPEG to invert DCIP data / Mag Gradiometry data and progress on implementing new receiver types in the natural source electromagnetic code
- 8 participants
- 1:11 hours
30 Jul 2019
This is the recording of the first hour of the SimPEG meeting in Montreal where we discussed the roadmap for SimPEG 1.0
- 2 participants
- 1:10 hours
30 Jul 2019
Progress on incorporating the simulation class into the potential fields codes, discussion on receiver classes, parallelizing the natural source EM codes and generating an example that uses GemPy to construct a geologic model.
- 6 participants
- 44 minutes
24 Jul 2019
Presentation by Jae on his work doing joint inversions with SimPEG
- 7 participants
- 39 minutes
9 Jul 2019
Discussion of high priority development items when the team is together in Montreal next week.
- 5 participants
- 50 minutes
2 Jul 2019
Seogi Kang provides an overview of inversions that consist of a parametric and a voxel model
Lindsey Heagy gives a demo of the simulation-refactor
Lindsey Heagy gives a demo of the simulation-refactor
- 9 participants
- 1:02 hours
22 Mar 2019
Friday Afternoon Talk Series - UCalgary - March 22, 2019
https://docs.google.com/presentation/d/1cWx__BqagqWZFEysZAlnwaSu_5jZe5wek8rPe3Vbhb8/edit?usp=sharing
https://docs.google.com/presentation/d/1cWx__BqagqWZFEysZAlnwaSu_5jZe5wek8rPe3Vbhb8/edit?usp=sharing
- 4 participants
- 36 minutes
8 Jan 2019
First SimPEG meeting of 2019! Dom provides an overview of developments in the potential fields codes and Lindsey introduces some ideas for how to maintain notebook-apps for multiple courses (and instructors).
- 7 participants
- 47 minutes
20 Nov 2018
Seogi Kang leads the meeting this week. To characterize large-scale hydrological structures at Edmonton to Calgary (ECC) corridor, time-domain airborne EM data are inverted, and a conductivity model is recovered (30km x 30 km). From the obtained conductivity model, we extract important hydrogeologic information about this region such as map of potential aquifer and 3D lithology model.
- 7 participants
- 48 minutes
13 Nov 2018
Bane Sullivan (Colorado School of Mines) shares his work on PVGeo, an open-source framework for visualizing geoscientific data and models in VTK powered platforms like ParaView.
The PVGeo Python package contains VTK powered tools for data visualization in geophysics which are wrapped for direct use within the application ParaView by Kitware. These tools are tailored to data visualization in the geosciences with a heavy focus on structured data sets like 2D or 3D time-varying grids. As an effort to bring interoperability between SimPEG projects and PVGeo, Bane will be demoing his work in SimPEG/discretize MR#114 to create an interface between discretize meshes and VTK data objects.
The PVGeo Python package contains VTK powered tools for data visualization in geophysics which are wrapped for direct use within the application ParaView by Kitware. These tools are tailored to data visualization in the geosciences with a heavy focus on structured data sets like 2D or 3D time-varying grids. As an effort to bring interoperability between SimPEG projects and PVGeo, Bane will be demoing his work in SimPEG/discretize MR#114 to create an interface between discretize meshes and VTK data objects.
- 5 participants
- 44 minutes
6 Nov 2018
Matthias Bücker gives us an overview of his work modelling IP and we explore trajectories for collaborating on IP modelling in SimPEG
- 11 participants
- 1:14 hours
23 Oct 2018
Catching up on SEG conversations and usage of SimPEG beyond UBC
- 11 participants
- 55 minutes
10 Jul 2018
As the SimPEG community grows, it will become helpful to document our processes (and have some processes in place!) for how we operate as a community.
There is a community repo that has been started, including two pull requests along with an issue that we can use as a starting point for the discussion.
There is a community repo that has been started, including two pull requests along with an issue that we can use as a starting point for the discussion.
- 10 participants
- 1:02 hours
27 Mar 2018
Work refactoring the DC code into the new Simulation framework
- 8 participants
- 49 minutes
20 Mar 2018
Johnathan Kuttai presented on “Time-Series Processing — From Raw Signal to IP Decay”
- 12 participants
- 1:05 hours
20 Feb 2018
Walk through of the FDEM upgrades as we move from the problem + survey to a simulation class
- 19 participants
- 1:16 hours
6 Feb 2018
Dikun Yang presents how he would like to use SimPEG for Marine CSEM inversions
- 8 participants
- 46 minutes
30 Jan 2018
- Plugging in empymod to SimPEG
- creating a Simulation class to replace the Problem + Survey
- parallelizing the code
- contributions to processing
- creating a Simulation class to replace the Problem + Survey
- parallelizing the code
- contributions to processing
- 11 participants
- 1:38 hours
6 Dec 2017
A framework for geophysical inversions with application to vadose zone parameter estimation.
- 1 participant
- 25 minutes
5 Sep 2017
Updates from the SimPEG team on current projects. Discussion on injecting properties into discretize and serializing meshes
- 11 participants
- 56 minutes
15 Aug 2017
updates from SciPy 2017 and the DISC presentations in Australia
- 7 participants
- 49 minutes
27 Jun 2017
Meeting on June 27.
Sustaining SimPEG. Ideas on how to keep momentum and grow SimPEG in the future.
Sustaining SimPEG. Ideas on how to keep momentum and grow SimPEG in the future.
- 9 participants
- 58 minutes
20 Jun 2017
Meeting on June 20
- updates on the EM paper
- fixes in the testing for deploying the docs
- discussion on a potential Mag project with Geoscience BC
- updates on the EM paper
- fixes in the testing for deploying the docs
- discussion on a potential Mag project with Geoscience BC
- 4 participants
- 38 minutes
13 Jun 2017
We chatted through a few issues with the usability of the regularization implementation and chatted through some updates coming in through em/dev.
- 7 participants
- 48 minutes
30 May 2017
May 30, 2017
- Releasing SimPEG 0.6.0
- A new example being added to the EM Paper (Seogi)
- Updates to em_examples and how we are using ipywidgets in EM Geosci (Thibaut)
- Releasing SimPEG 0.6.0
- A new example being added to the EM Paper (Seogi)
- Updates to em_examples and how we are using ipywidgets in EM Geosci (Thibaut)
- 10 participants
- 43 minutes
23 May 2017
We discuss the use of the Properties library for making SimPEG objects serializable and having validation on properties and classes.
- 11 participants
- 1:02 hours
16 May 2017
We do a dry run of the presentation that we will be giving at the BCGS.
- 10 participants
- 1:25 hours
2 May 2017
Dom and Seogi update us on their progress constructing the synthetic Kevista model.
We discuss the structure of GeoSci: computation and templates for the tutorials and case studies it contains.
We discuss the structure of GeoSci: computation and templates for the tutorials and case studies it contains.
- 10 participants
- 1:01 hours
18 Apr 2017
Seogi discusses the theory and implementation of the line-source formulation, and some of the challenges and opportunities.
- 4 participants
- 52 minutes
1 Mar 2017
Showing SimPEG (http://simpeg.xyz) and motivation behind quantitative geoscience communication.
- 1 participant
- 13 minutes
28 Feb 2017
SimPEG meeting on Feb 28. Thibaut Astic presents on Designing objective functions: a probability approach
- 5 participants
- 52 minutes
10 Jan 2017
SimPEG meeting on January 10, 2017. Rowan introduces the meshing package discretize and updates us on pymatsolver changes
- 5 participants
- 34 minutes
15 Jul 2016
Geophysical inversions are tools for constructing models of the subsurface (images) given a finite amount of data. SimPEG (http://simpeg.xyz) is an effort to synthesize geophysical forward and inverse methodologies into a consistent framework. We will show seven geophysical methods based around a diamond exploration case study, combining the results to drive a more informed decision. Slides may be found here: https://docs.google.com/presentation/d/1O6A85QwnnAibm7CsV2_VZ95HcRsnLyvPQIRYIz_K5xg/edit#slide=id.g15d5208fb1_2_246
- 2 participants
- 25 minutes