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

Tariq Alkhalifah

Tariq Alkhalifah contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.geo-phMachine Learningphysics.comp-phmath.OCmath.PRphysics.app-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
6topics
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

12 published item(s)

preprint2022arXiv

$ω$-FWI: Robust full-waveform inversion with Fourier-based metric

Full-waveform inversion is a cutting-edge methodology for recovering high-resolution subsurface models. However, one of the main conventional full-waveform optimization problems challenges is cycle-skipping, usually leading us to an inaccurate local minimum model. A highly investigated track to alleviate this challenge involves designing a more global measure of misfit between the observed and modelled data beyond the sample-to-sample comparison. However, most of these approaches admit relatively smooth inversion results. Here, we introduce a novel misfit function based on the Fourier-based metric. This metric has been successfully applied in molecular physics for solving the Boltzmann equation, and we adapt it to full-waveform inversion. This misfit function exploits the power spectrum information between the modelled and observed data to provide low-wavenumber velocity model updates early, and more high resolution updates as we approach the solution. Thus, it also can be reformulated as a weighted $\ell_{2}$-norm in a quadratic case, which can be seen as a simple extension for conventional full-waveform inversion. Thus, despite its robustness to cycle skipping, it is capable of delivering high-resolution models synonymous to conventional FWI. Considering its frequency domain utilization, we refer to this inversion method as $ω$-FWI. Through the synthetic Marmousi model example, this method successfully recovers an accurate velocity model, starting from a linearly increasing model even for the case of noisy observed data and the lack of low frequencies below 3 Hz and 5Hz, in which the conventional $\ell_{2}$-norm full-waveform inversion suffers from cycle skipping.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Coherent noise suppression via a self-supervised blind-trace deep learning scheme

Coherent noise regularly plagues seismic recordings, causing artefacts and uncertainties in products derived from down-the-line processing and imaging tasks. The outstanding capabilities of deep learning in denoising of natural and medical images have recently spur a number of applications of neural networks in the context of seismic data denoising. A limitation of the majority of such methods is that the deep learning procedure is supervised and requires clean (noise-free) data as a target for training the network. Blindspot networks were recently proposed to overcome this requirement, allowing training to be performed directly on the noisy field data as a powerful suppressor of random noise. A careful adaptation of the blind-spot methodology allows for an extension to coherent noise suppression. In this work, we expand the methodology of blind-spot networks to create a blind-trace network that successfully removes trace-wise coherent noise. Through an extensive synthetic analysis, we illustrate the denoising procedure's robustness to varying noise levels, as well as varying numbers of noisy traces within shot gathers. It is shown that the network can accurately learn to suppress the noise when up to 60% of the original traces are noisy. Furthermore, the proposed procedure is implemented on the Stratton 3D field dataset and is shown to restore the previously corrupted direct arrivals. Our adaptation of the blind-spot network for self-supervised, trace-wise noise suppression could lead to other use-cases such as the suppression of coherent noise arising from wellsite activity, passing vessels or nearby industrial activity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deep learning unflooding for robust subsalt waveform inversion

Full-waveform inversion (FWI), a popular technique that promises high-resolution models, has helped in improving the salt definition in inverted velocity models. The success of the inversion relies heavily on having prior knowledge of the salt, and using advanced acquisition technology with long offsets and low frequencies. Salt bodies are often constructed by recursively picking the top and bottom of the salt from seismic images corresponding to tomography models, combined with flooding techniques. The process is time-consuming and highly prone to error, especially in picking the bottom of the salt (BoS). Many studies suggest performing FWI with long offsets and low frequencies after constructing the salt bodies to correct the miss-interpreted boundaries. Here, we focus on detecting the BoS automatically by utilizing deep learning tools. We specifically generate many random 1D models, containing or free of salt bodies, and calculate the corresponding shot gathers. We then apply FWI starting with salt flooded versions of those models, and the results of the FWI become inputs to the neural network, whereas the corresponding true 1D models are the output. The network is trained in a regression manner to detect the BoS and estimate the subsalt velocity. We analyze three scenarios in creating the training datasets and test their performance on the 2D BP 2004 salt model. We show that when the network succeeds in estimating the subsalt velocity, the requirement of low frequencies and long offsets are somewhat mitigated. In general, this work allows us to merge the top-to-bottom approach with FWI, save the BoS picking time, and empower FWI to converge in the absence of low frequencies and long offsets in the data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Laplace HypoPINN: Physics-Informed Neural Network for hypocenter localization and its predictive uncertainty

Several techniques have been proposed over the years for automatic hypocenter localization. While those techniques have pros and cons that trade-off computational efficiency and the susceptibility of getting trapped in local minima, an alternate approach is needed that allows robust localization performance and holds the potential to make the elusive goal of real-time microseismic monitoring possible. Physics-informed neural networks (PINNs) have appeared on the scene as a flexible and versatile framework for solving partial differential equations (PDEs) along with the associated initial or boundary conditions. We develop HypoPINN -- a PINN-based inversion framework for hypocenter localization and introduce an approximate Bayesian framework for estimating its predictive uncertainties. This work focuses on predicting the hypocenter locations using HypoPINN and investigates the propagation of uncertainties from the random realizations of HypoPINN's weights and biases using the Laplace approximation. We train HypoPINN to obtain the optimized weights for predicting hypocenter location. Next, we approximate the covariance matrix at the optimized HypoPINN's weights for posterior sampling with the Laplace approximation. The posterior samples represent various realizations of HypoPINN's weights. Finally, we predict the locations of the hypocenter associated with those weights' realizations to investigate the uncertainty propagation that comes from those realisations. We demonstrate the features of this methodology through several numerical examples, including using the Otway velocity model based on the Otway project in Australia.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Posterior sampling with CNN-based, Plug-and-Play regularization with applications to Post-Stack Seismic Inversion

Uncertainty quantification is crucial to inverse problems, as it could provide decision-makers with valuable information about the inversion results. For example, seismic inversion is a notoriously ill-posed inverse problem due to the band-limited and noisy nature of seismic data. It is therefore of paramount importance to quantify the uncertainties associated to the inversion process to ease the subsequent interpretation and decision making processes. Within this framework of reference, sampling from a target posterior provides a fundamental approach to quantifying the uncertainty in seismic inversion. However, selecting appropriate prior information in a probabilistic inversion is crucial, yet non-trivial, as it influences the ability of a sampling-based inference in providing geological realism in the posterior samples. To overcome such limitations, we present a regularized variational inference framework that performs posterior inference by implicitly regularizing the Kullback-Leibler divergence loss with a CNN-based denoiser by means of the Plug-and-Play methods. We call this new algorithm Plug-and-Play Stein Variational Gradient Descent (PnP-SVGD) and demonstrate its ability in producing high-resolution, trustworthy samples representative of the subsurface structures, which we argue could be used for post-inference tasks such as reservoir modelling and history matching. To validate the proposed method, numerical tests are performed on both synthetic and field post-stack seismic data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Time-lapse data matching using a recurrent neural network approach

Time-lapse seismic data acquisition is an essential tool to monitor changes in a reservoir due to fluid injection, such as CO$_2$ injection. By acquiring multiple seismic surveys in the exact location, we can identify the reservoir changes by analyzing the difference in the data. However, such analysis can be skewed by the near-surface seasonal velocity variations, inaccuracy in the acquisition parameters, and other inevitable noise. The common practice (cross-equalization) to address this problem uses the part of the data where changes are not expected to design a matching filter and then apply it to the whole data, including the reservoir area. Like cross-equalization, we train a recurrent neural network on parts of the data excluding the reservoir area and then infer the reservoir-related data. The recurrent neural network can learn the time dependency of the data, unlike the matching filter that processes the data based on the local information obtained in the filter window. We demonstrate the method of matching the data in various examples and compare it with the conventional matching filter. Specifically, we start by demonstrating the ability of the approach in matching two traces and then test the method on a pre-stack 2D synthetic data. Then, we verify the enhancements of the 4D signal by providing RTM images. We measure the repeatability using normalized root mean square and predictability metrics and show that in some cases, it overcomes the matching filter approach.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

PINNup: Robust neural network wavefield solutions using frequency upscaling and neuron splitting

Solving for the frequency-domain scattered wavefield via physics-informed neural network (PINN) has great potential in seismic modeling and inversion. However, when dealing with high-frequency wavefields, its accuracy and training cost limits its applications. Thus, we propose a novel implementation of PINN using frequency upscaling and neuron splitting, which allows the neural network model to grow in size as we increase the frequency while leveraging the information from the pre-trained model for lower-frequency wavefields, resulting in fast convergence to high-accuracy solutions. Numerical results show that, compared to the commonly used PINN with random initialization, the proposed PINN exhibits notable superiority in terms of convergence and accuracy and can achieve neuron based high-frequency wavefield solutions with a two-hidden-layer model.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Solving the Eikonal equation for compressional and shear waves in anisotropic media using peridynamic differential operator

The traveltime of compressional (P) and shear (S) waves have proven essential in many applications of earthquake and exploration seismology. An accurate and efficient traveltime computation for P and S waves is crucial for the success of these applications. However, solutions to the Eikonal equation with a complex phase velocity field in anisotropic media is challenging. The Eikonal equation is a first-order, hyperbolic, nonlinear partial differential equation (PDE) that represents the high-frequency asymptotic approximation of the wave equation. The fast marching and sweeping methods are commonly used due to their efficiency in numercally solving Eikonal equation. However, these methods suffer from numerical inaccuracy in anisotropic media with sharp heterogeneity, irregular surface topography and complex phase velocity fields. This study presents a new method to solving the Eikonal equation by employing the peridynamic differential operator (PDDO). The PDDO provides the nonlocal form of the Eikonal equation by introducing an internal length parameter (horizon) and a weight function with directional nonlocality. The operator is immune to discontinuities in the form sharp changes in field or model variables and invokes the direction of traveltime in a consistent manner. The weight function controls the degree of association among points within the horizon. Solutions are constructed in a consistent manner without upwind assumptions through simple discretization. The capability of this approach is demonstrated by considering different types of Eikonal equations on complex velocity models in anisotropic media. The examples demonstrate its unconditional numerical stability and results compare well with the reference solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Wavefield reconstruction inversion via physics-informed neural networks

Wavefield reconstruction inversion (WRI) formulates a PDE-constrained optimization problem to reduce cycle skipping in full-waveform inversion (FWI). WRI often requires expensive matrix inversions to reconstruct frequency-domain wavefields. Physics-informed neural network (PINN) uses the underlying physical laws as loss functions to train the neural network (NN), and it has shown its effectiveness in solving the Helmholtz equation and generating Green's functions, specifically for the scattered wavefield. By including a data-constrained term in the loss function, the trained NN can reconstruct a wavefield that simultaneously fits the recorded data and satisfies the Helmholtz equation for a given initial velocity model. Using the predicted wavefields, we rely on a small-size NN to predict the velocity using the reconstructed wavefield. In this velocity prediction NN, spatial coordinates are used as input data to the network and the scattered Helmholtz equation is used to define the loss function. After we train this network, we are able to predict the velocity in the domain of interest. We develop this PINN-based WRI method and demonstrate its potential using a part of the Sigsbee2A model and a modified Marmousi model. The results show that the PINN-based WRI is able to invert for a reasonable velocity with very limited iterations and frequencies, which can be used in a subsequent FWI application.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A data-driven choice of misfit function for FWI using reinforcement learning

In the workflow of Full-Waveform Inversion (FWI), we often tune the parameters of the inversion to help us avoid cycle skipping and obtain high resolution models. For example, typically start by using objective functions that avoid cycle skipping, like tomographic and image based or using only low frequency, and then later, we utilize the least squares misfit to admit high resolution information. We also may perform an isotropic (acoustic) inversion to first update the velocity model and then switch to multi-parameter anisotropic (elastic) inversions to fully recover the complex physics. Such hierarchical approaches are common in FWI, and they often depend on our manual intervention based on many factors, and of course, results depend on experience. However, with the large data size often involved in the inversion and the complexity of the process, making optimal choices is difficult even for an experienced practitioner. Thus, as an example, and within the framework of reinforcement learning, we utilize a deep-Q network (DQN) to learn an optimal policy to determine the proper timing to switch between different misfit functions. Specifically, we train the state-action value function (Q) to predict when to use the conventional L2-norm misfit function or the more advanced optimal-transport matching-filter (OTMF) misfit to mitigate the cycle-skipping and obtain high resolution, as well as improve convergence. We use a simple while demonstrative shifted-signal inversion examples to demonstrate the basic principles of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

ML-misfit: Learning a robust misfit function for full-waveform inversion using machine learning

Most of the available advanced misfit functions for full waveform inversion (FWI) are hand-crafted, and the performance of those misfit functions is data-dependent. Thus, we propose to learn a misfit function for FWI, entitled ML-misfit, based on machine learning. Inspired by the optimal transport of the matching filter misfit, we design a neural network (NN) architecture for the misfit function in a form similar to comparing the mean and variance for two distributions. To guarantee the resulting learned misfit is a metric, we accommodate the symmetry of the misfit with respect to its input and a Hinge loss regularization term in a meta-loss function to satisfy the "triangle inequality" rule. In the framework of meta-learning, we train the network by running FWI to invert for randomly generated velocity models and update the parameters of the NN by minimizing the meta-loss, which is defined as accumulated difference between the true and inverted models. We first illustrate the basic principle of the ML-misfit for learning a convex misfit function for travel-time shifted signals. Further, we train the NN on 2D horizontally layered models, and we demonstrate the effectiveness and robustness of the learned ML-misfit by applying it to the well-known Marmousi model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solving the acoustic VTI wave equation using physics-informed neural networks

Frequency-domain wavefield solutions corresponding to the anisotropic acoustic wave equations can be used to describe the anisotropic nature of the earth. To solve a frequency-domain wave equation, we often need to invert the impedance matrix. This results in a dramatic increase in computational cost as the model size increases. It is even a bigger challenge for anisotropic media, where the impedance matrix is far more complex. To address this issue, we use the emerging paradigm of physics-informed neural networks (PINNs) to obtain wavefield solutions for an acoustic wave equation for transversely isotropic (TI) media with a vertical axis of symmetry (VTI). PINNs utilize the concept of automatic differentiation to calculate its partial derivatives. Thus, we use the wave equation as a loss function to train a neural network to provide functional solutions to form of the acoustic VTI wave equation. Instead of predicting the pressure wavefields directly, we solve for the scattered pressure wavefields to avoid dealing with the point source singularity. We use the spatial coordinates as input data to the network, which outputs the real and imaginary parts of the scattered wavefields and auxiliary function. After training a deep neural network (NN), we can evaluate the wavefield at any point in space instantly using this trained NN. We demonstrate these features on a simple anomaly model and a layered model. Additional tests on a modified 3D Overthrust model and a model with irregular topography also show the effectiveness of the proposed method.

0 citations0 reviewsNo code linkedOpen access