BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Andrew Curtis

Andrew Curtis contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.geo-phphysics.comp-phHuman-Computer Interaction

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2026arXiv

Direct-3D Variational Bayesian Surface Wave Inversion and Its Application to Ambient Noise Tomography beneath Great Britain

We present a new, variational, fully nonlinear, probabilistic ambient noise tomography method, which estimates subsurface structure and quantifies the corresponding uncertainties directly in three dimensions (3D) from inter-receiver seismic surface wave dispersion data. We use the method to invert for high resolution 3D seismic velocity models of the upper crust beneath Great Britain using seismic ambient noise data recorded around the region - a task that proved too high-dimensional and hence computationally demanding for Monte Carlo sampling to converge to a stable solution. We compare the inversion results from the new method to those obtained from two standard, indirect inversion methods, in which 2D (geographical) surface wave velocity maps and 1D (depth) shear velocity profiles are estimated in two separate, consecutive steps. The results show that the direct-3D scheme preserves better lateral continuity and produces better data fit than the two-step methods, and provides information about lateral correlations that is absent from the two-step solutions. The inversion results are consistent with large-scale geology of Great Britain, and for the first time provide seismologically-imaged evidence of the Great Glen Fault and other major tectonic faults. We therefore propose that direct-3D inversion schemes should be used where possible for surface wave inversion as they provide improved results at little additional computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Bayesian Inversion, Uncertainty Analysis and Interrogation using Boosting Variational Inference

Geoscientists use observed data to estimate properties of the Earth's interior. This often requires non-linear inverse problems to be solved and uncertainties to be estimated. Bayesian inference solves inverse problems under a probabilistic framework, in which uncertainty is represented by a so-called posterior probability distribution. Recently, variational inference has emerged as an efficient method to estimate Bayesian solutions. By seeking the closest approximation to the posterior distribution within any chosen family of distributions, variational inference yields a fully probabilistic solution. It is important to define expressive variational families so that the posterior distribution can be represented accurately. We introduce boosting variational inference (BVI) as a computationally efficient means to construct a flexible approximating family comprising all possible finite mixtures of simpler component distributions. We use Gaussian mixture components due to their fully parametric nature and the ease with which they can be optimised. We apply BVI to seismic travel time tomography and full waveform inversion, comparing its performance with other methods of solution. The results demonstrate that BVI achieves reasonable efficiency and accuracy while enabling the construction of a fully analytic expression for the posterior distribution. Samples that represent major components of uncertainty in the solution can be obtained analytically from each mixture component. We demonstrate that these samples can be used to solve an interrogation problem: to assess the size of a subsurface target structure. To the best of our knowledge, this is the first method in geophysics that provides both analytic and reasonably accurate probabilistic solutions to fully non-linear, high-dimensional Bayesian full waveform inversion problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An introduction to variational inference in Geophysical inverse problems

In a variety of scientific applications we wish to characterize a physical system using measurements or observations. This often requires us to solve an inverse problem, which usually has non-unique solutions so uncertainty must be quantified in order to define the family of all possible solutions. Bayesian inference provides a powerful theoretical framework which defines the set of solutions to inverse problems, and variational inference is a method to solve Bayesian inference problems using optimization while still producing fully probabilistic solutions. This chapter provides an introduction to variational inference, and reviews its applications to a range of geophysical problems, including petrophysical inversion, travel time tomography and full-waveform inversion. We demonstrate that variational inference is an efficient and scalable method which can be deployed in many practical scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modular interface for managing cognitive bias in experts

Expert knowledge is required to interpret data across a range of fields. Experts bridge gaps that often exists in our knowledge about relationships between data and the parameters of interest. This is especially true in geoscientific applications, where knowledge of the Earth is derived from interpretations of observable features and relies on predominantly unproven but widely accepted theories. Thus, experts facilitate solutions to otherwise unsolvable problems. However, experts are inherently subjective, and susceptible to cognitive biases and adverse external effects. This work examines this problem within geoscience. Three compelling examples are provided of the prevalence of cognitive biases from previous work. The problem is then formally defined, and a set of design principles which ensure that any solution is sufficiently flexible to be readily applied to the range of geoscientific problems. No solutions exist that reliably capture and reduce cognitive bias in experts. However, formal expert elicitation methods can be used to assess expert variation, and a variety of approaches exist that may help to illuminate uncertainties, avoid misunderstandings, and reduce herding behaviours or single-expert over-dominance. This work combines existing and future approaches to reduce expert suboptimality through a flexible modular design where each module provides a specific function. The design centres around action modules that force a stop-and-perform step into interpretation tasks. A starter-pack of modules is provided as an example of the conceptual design. This simple bias-reduction system may readily be applied in organisations and during everyday interpretations through to tasks for major commercial ventures.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Surface wave dispersion inversion using an energy likelihood function

Seismic surface wave dispersion inversion is used widely to study the subsurface structure of the Earth. The dispersion property is usually measured by using frequency-phase velocity (f-c) analysis and by picking phase velocities from the obtained f-c spectrum. However, because of potential contamination the f-c spectrum often has multimodalities at each frequency for each mode. These introduce uncertainty and errors in the picked phase velocities, and consequently the obtained shear velocity structure is biased. To overcome this issue, in this study we introduce a new method which directly uses the spectrum as data. We achieve this by solving the inverse problem in a Bayesian framework and define a new likelihood function, the energy likelihood function, which uses the spectrum energy to define data fit. We apply the new method to a land dataset recorded by a dense receiver array, and compare the results to those obtained using the traditional method. The results show that the new method produces more accurate results since they better match independent data from refraction tomography. This real-data application also shows that it can be applied efficiently since it removes the need to pick phase velocities, and with relatively little adjustment to current practice since it uses standard f-c panels to define the likelihood. We therefore recommend using the energy likelihood function rather than explicitly picking phase velocities in surface wave dispersion inversion.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Interrogating probabilistic inversion results for subsurface structural information

The goal of a scientific investigation is to find answers to specific questions. In geosciences this is typically achieved by solving an inference or inverse problem and interpreting the solution. However, the answer obtained is often biased because the solution to an inverse problem is nonunique and human interpretation is a biased process. Interrogation theory provides a systematic way to find optimal answers by considering their full uncertainty estimates, and by designing an objective function that defines desirable qualities in the answer. In this study we demonstrate interrogation theory by quantifying the size of a particular subsurface structure. The results show that interrogation theory provides an accurate estimate of the true answer, which cannot be obtained by direct, subjective interpretation of the solution mean and standard deviation. This demonstrates the value of interrogation theory. It also shows that fully nonlinear uncertainty assessments may be critical in order to address real-world scientific problems, which goes some way towards justifying their computational expense.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Probabilistic Neural-Network Based 2D Travel Time Tomography

Travel time tomography for the velocity structure of a medium is a highly non-linear and non-unique inverse problem. Monte Carlo methods are becoming increasingly common choices to provide probabilistic solutions to tomographic problems but those methods are computationally expensive. Neural networks can often be used to solve highly non-linear problems at a much lower computational cost when multiple inversions are needed from similar data types. We present the first method to perform fully non-linear, rapid and probabilistic Bayesian inversion of travel time data for 2D velocity maps using a mixture density network. We compare multiple methods to estimate probability density functions that represent the tomographic solution, using different sets of prior information and different training methodologies. We demonstrate the importance of prior information in such high dimensional inverse problems due to the curse of dimensionality: unrealistically informative prior probability distributions may result in better estimates of the mean velocity structure, however the uncertainties represented in the posterior probability density functions then contain less information than is obtained when using a less informative prior. This is illustrated by the emergence of uncertainty loops in posterior standard deviation maps when inverting travel time data using a less informative prior, which are not observed when using networks trained on prior information that includes (unrealistic) a priori smoothness constraints in the velocity models. We show that after an expensive program of training the networks, repeated high-dimensional, probabilistic tomography is possible on timescales of the order of a second on a standard desktop computer.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Probabilistic Neural Network Tomography across Grane field (North Sea) from Surface Wave Dispersion Data

Surface wave tomography uses measured dispersion properties of surface waves to infer the spatial distribution of subsurface properties such as shear-wave velocities. These properties can be estimated vertically below any geographical location at which surface wave dispersion data are available. As the inversion is significantly non-linear, Monte Carlo methods are often used to invert dispersion curves for shear-wave velocity profiles with depth to give a probabilistic solution. Such methods provide uncertainty information but are computationally expensive. Neural network based inversion provides a more efficient way to obtain probabilistic solutions when those solutions are required beneath many geographical locations. Unlike Monte Carlo methods, once a network has been trained it can be applied rapidly to perform any number of inversions. We train a class of neural networks called mixture density networks, to invert dispersion curves for shear-wave velocity models and their non-linearised uncertainty. Mixture density networks are able to produce fully probabilistic solutions in the form of weighted sums of multivariate analytic kernels such as Gaussians, and we show that including data uncertainties in the mixture density network gives more reliable mean velocity estimates when data contains significant noise. The networks were applied to data from the Grane field in the Norwegian North sea to produce shear-wave velocity maps at several depth levels. Post-training we obtained probabilistic velocity profiles with depth beneath 26,772 locations to produce a 3D velocity model in 21 seconds on a standard desktop computer. This method is therefore ideally suited for rapid, repeated 3D subsurface imaging and monitoring.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Seismic tomography using variational inference methods

Seismic tomography is a methodology to image the interior of solid or fluid media, and is often used to map properties in the subsurface of the Earth. In order to better interpret the resulting images it is important to assess imaging uncertainties. Since tomography is significantly nonlinear, Monte Carlo sampling methods are often used for this purpose, but they are generally computationally intractable for large datasets and high-dimensional parameter spaces. To extend uncertainty analysis to larger systems we use variational inference methods to conduct seismic tomography. In contrast to Monte Carlo sampling, variational methods solve the Bayesian inference problem as an optimization problem, yet still provide probabilistic results. In this study, we applied two variational methods, automatic differential variational inference (ADVI) and Stein variational gradient descent (SVGD), to 2D seismic tomography problems using both synthetic and real data and we compare the results to those from two different Monte Carlo sampling methods. The results show that variational inference methods can produce accurate approximations to the results of Monte Carlo sampling methods at significantly lower computational cost, provided that gradients of parameters with respect to data can be calculated efficiently. We expect that the methods can be applied fruitfully to many other types of geophysical inverse problems.

0 citations0 reviewsNo code linkedOpen access