Efficient 3D Gravity and Magnetic Modeling

Size: px
Start display at page:

Download "Efficient 3D Gravity and Magnetic Modeling"

Transcription

1 Efficient 3D Gravity and Magnetic Modeling X. Li Fugro Gravity & Magnetic Services Inc., Houston, Texas, USA Summary There are many different spatial-domain algorithms for 3D gravity and magnetic forward and inverse modeling. However, an efficient modeling tool for petroleum exploration needs algorithms in the wavenumber domain. We have developed, tested and applied such a tool, for both forward and inverse modeling, over two decades. Our algorithms work for density and susceptibility variations of any complex form; compute gravity, gravity gradient and total magnetic intensity (TMI) responses; and invert one field or gradient alone or any combination of these field and gradient components simultaneously. It takes only minutes, not many hours, to complete a joint inversion for structure of a practically sized project on a personal computer. Introduction Spatial-domain closed-form or numerical computation formulae for forward gravity and magnetic modeling have been extensively studied. People often represent an isolated body by simple geometries: an ellipsoid, sphere, cylinder, thin sheet, etc. A complex body or structure is expressed by a combination of right rectangular prisms, polygonal prisms, or polyhedrons (Li and Chouteau, 1998). Researchers have also developed formulae that allow variable densities and susceptibilities within a prism or polyhedron. These formulae are accurate but inefficient. For example, a widely used formula for computation of the gravity due to a right rectangular prism with a constant density contains 24 terms: 16 logarithms and 8 arctangents (Li and Chouteau, 1998, p. 344). Gravity or magnetic data may be inverted for either physical property or structure. In a property inversion, we often divide the subsurface into cells and invert for the constant density or susceptibility values of the cells. This is a linear inversion as in seismic tomography. However, the structure inversion is a much-preferred choice in petroleum exploration uses of gravity and magnetic data. Explorationists use such an inversion to determine a structure, i.e., the depth variation of a boundary such as the basement or base of salt, when there is not enough seismic data or its quality is poor. The structure inversion is nonlinear as in seismic depth imaging. The popular approach is to linearize a nonlinear problem and then solve the linear system in a least-squares sense. This approach requires many iterations of forward computation and solution of a linear system. The size of this system often gets huge for a field project: its number of rows is the total number of data points and its number of columns is the total number of unknown parameters. Researchers have thus designed many advanced mathematical solutions, with a focus on two aspects: (a) transforming a dense matrix into a sparse one (with many zero elements in the matrix), e.g., by the wavelet compression technique, and (b) using an effective solver of a sparse matrix, e.g., the conjugate gradient or LSQR method. Unfortunately, products from such great efforts are far from efficient, and it is still common to take many hours to run a 3D inversion. Gravity and magnetics are low-cost methods. A computer cluster is now widely used for seismic processing and interpretation while a laptop computer is still a popular machine for gravity and

2 magnetic modeling. At the same time, users don t want to wait hours for a modeling result. For this reason, we have to deviate from the conventional approaches to forward and inverse modeling and seek solutions in the wavenumber domain. Methodology In petroleum exploration, we deal with layered structures as well as isolated bodies such as igneous intrusions or salt emplacements. The density or susceptibility within a body, particularly a layer, may vary both vertically and horizontally. An advanced modeling tool should allow different forms of density variations within a layer. First, the simplest form is a constant and this is widely used when variable density information is not available. Second, densities of sediments increase with depth due to compaction. Density changes have thus been represented by a linear, exponential, quadratic, parabolic, or hyperbolic function of depth. A depth-density pair may be used to describe these functions. Third, the density in a sedimentary basin may also change horizontally, e.g., the density may decrease basinward from a major fault. A density grid is a general form for describing such a horizontal variation. Fourth, a combination of a horizontal density grid and a depth-density pair allows both vertical and horizontal variations. However, this combination involves one horizontal grid and one single depth-density pair and is not applicable to the most complex and the most general cases. Fifth, such cases can be defined only by a voxet. In practice, a velocity voxet may be available when gravity and/or gravity gradient data is used to help interpret a horizon that seismic encounters difficulties (e.g., base salt). We can convert this velocity voxet into a density voxet using a simple velocity-density relationship or a geostatistical analysis. Evidently, the first four forms are special cases of this last one. A voxet may be divided vertically into a number of laminas (thin sheets). Each lamina is equivalent to a surface density distribution σ on a horizontal plane. The gravity effect on the top surface of a lamina can be written in the wavenumber domain as F[ g( K )] = 2πγF[ σ ( K) ] (1) where γ is the gravitational constant, K the wavenumber vector, and F stands for the Fourier transform. Applying an upward continuation, we obtain the gravity effect on a horizontal plane at altitude z: K z F[ g( K) ] = 2 πγe F[ σ ( K) ] (2) Summing up the effects of all laminas produces the total gravity effect. Applying operations in the wavenumber domain further computes the gravity gradient components TXX, TXY, TXZ, TYY, TYZ, T, T = ( T T 2 and the magnetic effect Δ T. ZZ UV YY XX ) For inversion, Bott (1960) proposed a very simple update formula for interpretation of gravity anomalies over a sedimentary basin. The gravity effect due to an infinite horizontal slab is a constant independent of the vertical position (i.e., burial depth) and the horizontal location of observation. g = 2πγρt (3) where ρ is the bulk density (contrast) of the slab, and t the thickness of the slab. Rearranging equation (3) produces g t = (4) 2πγρ This equation can be used to compute the depth change of a horizon, point by point. The depth inversion is accomplished with an iterative process.

3 More sophisticated algorithms for forward and inverse modeling are needed in practice. For example, equation (4) cannot be used to update depth when the magnetic anomaly or a gravity gradient component is inverted. We have developed advanced and efficient algorithms that can invert not only a field or gradient component individually but also a selection of any number of field and gradient components simultaneously. Below I present a joint inversion example. Figure 1. (a) The base of salt and (b) the density distribution of the starting model used to test a joint inversion of gravity gradients T and T. XY UV Figure 2. The model (a) before and (b) after a joint inversion of TXY. Only TXY grids are shown. One example Our test model consists of a seawater layer, sediments, a salt body and a basement (see Figure 3a). Each of the depth grids is quite large: 960 rows by 960 columns. The sediments have a complex density variation, defined by a voxet. Figure 1a shows the depth to the base of salt and Figure 1b depicts the density distribution of the initial 3D model. Figure 2a displays maps and east-west crosssections of this starting model. We simulate TXY, the curvature gravity gradient components that are directly measured by the Falcon TM airborne gravity gradiometer system. We display T XY

4 responses only although a joint inversion of TXY is performed to make changes to the base of salt so that the calculated gradients match observed gradients. Figure 2b demonstrates the inverted results. We have a seismic section in this test and its location is shown as the NW-SE diagonal line in map windows of Figure 2b. Figure 3 illustrates the cross-section along this diagonal profile. In particular, density variations overlie seismic section in Figure 3b. The algorithms converge fast: it takes 10 iterations (13 minutes) to complete this joint inversion of TXY on a desktop computer (Dell Precision T3400n Workstation with 64-bit Intel Core TM 2 Duo E8500 processor, 3.16 GHz and 8 GB RAM). Figure 3. A cross-section along a seismic profile. The NW-SE diagonal line in maps of Figure 2b indicates the profile location. (a) Structure and (b) densities (in color) overlaid on seismic. Conclusions We have developed wavenumber-domain algorithms for 3D gravity and magnetic modeling. It works effectively with the complexity of petroleum exploration problems and efficiently for the size of a field project. The inversion is always nonunique and depends strongly on the starting model and constraints. An efficient inversion allows users to test different reasonable assumptions and thus to derive a reliable final model. Acknowledgments Many people at Fugro Gravity & Magnetic Services Inc. have contributed to the development of this efficient 3D modeling tool. I also thank Fugro Gravity & Magnetic Services Inc. for permission to publish this work.

5 References Bott, M. H. P., 1960, The use of rapid digital computing methods for direct gravity interpretation of sedimentary basins: Geophysical Journal of the Royal Astronomical Society, 3, Li, X., and Chouteau, M., 1998, Three-dimensional gravity modeling in all space: Surveys in Geophysics, 19,

Geophysics 224 B2. Gravity anomalies of some simple shapes. B2.1 Buried sphere

Geophysics 224 B2. Gravity anomalies of some simple shapes. B2.1 Buried sphere Geophysics 4 B. Gravity anomalies of some simple shapes B.1 Buried sphere Gravity measurements are made on a surface profile across a buried sphere. The sphere has an excess mass M S and the centre is

More information

prismatic discretization of the digital elevation model, and their associated volume integration problems. Summary

prismatic discretization of the digital elevation model, and their associated volume integration problems. Summary A new method of terrain correcting airborne gravity gradiometry data using 3D Cauchy-type integrals Michael S. Zhdanov*, University of Utah and TechnoImaging, Glenn A. Wilson, TechnoImaging, and Xiaojun

More information

(x, y, z) m 2. (x, y, z) ...] T. m 2. m = [m 1. m 3. Φ = r T V 1 r + λ 1. m T Wm. m T L T Lm + λ 2. m T Hm + λ 3. t(x, y, z) = m 1

(x, y, z) m 2. (x, y, z) ...] T. m 2. m = [m 1. m 3. Φ = r T V 1 r + λ 1. m T Wm. m T L T Lm + λ 2. m T Hm + λ 3. t(x, y, z) = m 1 Class 1: Joint Geophysical Inversions Wed, December 1, 29 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical, statistical Deal with data in the same

More information

GEOPHYS 242: Near Surface Geophysical Imaging. Class 8: Joint Geophysical Inversions Wed, April 20, 2011

GEOPHYS 242: Near Surface Geophysical Imaging. Class 8: Joint Geophysical Inversions Wed, April 20, 2011 GEOPHYS 4: Near Surface Geophysical Imaging Class 8: Joint Geophysical Inversions Wed, April, 11 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical,

More information

Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah

Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation

More information

SUMMARY. method to synthetic datasets is discussed in the present paper.

SUMMARY. method to synthetic datasets is discussed in the present paper. Geophysical modeling through simultaneous Joint Inversion of Seismic, Gravity and Magnetotelluric data Michele De Stefano (1), Daniele Colombo (1) WesternGeco EM - Geosystem, via Clericetti 42/A, 20133

More information

B15 Enhancement of Linear Features from Gravity Anomalies by Using Curvature Gradient Tensor Matrix

B15 Enhancement of Linear Features from Gravity Anomalies by Using Curvature Gradient Tensor Matrix B5 Enhancement of Linear Features from Gravity Anomalies by Usin Curvature Gradient Tensor Matrix B. Oruç* (Kocaeli University) SUMMARY In this study, a new ede enhancement technique based on the eienvalues

More information

Robust 3D gravity gradient inversion by planting anomalous densities

Robust 3D gravity gradient inversion by planting anomalous densities Robust 3D gravity gradient inversion by planting anomalous densities Leonardo Uieda Valéria C. F. Barbosa Observatório Nacional September, 2011 Outline Outline Forward Problem Outline Forward Problem Inverse

More information

Application of wavelet theory to the analysis of gravity data. P. Hornby, F. Boschetti* and F. Horowitz, Division of Exploration and Mining, CSIRO,

Application of wavelet theory to the analysis of gravity data. P. Hornby, F. Boschetti* and F. Horowitz, Division of Exploration and Mining, CSIRO, Application of wavelet theory to the analysis of gravity data. P. Hornby, F. Boschetti* and F. Horowitz, Division of Exploration and Mining, CSIRO, Australia. Summary. The fundamental equations of potential

More information

Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers

Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers J.W.C. Sherwood* (PGS), K. Sherwood (PGS), H. Tieman (PGS), R. Mager (PGS) & C. Zhou

More information

Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction

Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to

More information

Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines

Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic

More information

Downloaded 18 Jul 2011 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 18 Jul 2011 to Redistribution subject to SEG license or copyright; see Terms of Use at Rapid gravity and gravity gradiometry terrain correction via adaptive quadtree mesh discretization Kristofer Davis, M. Andy Kass, and Yaoguo Li, Center for Gravity, Electrical and Magnetic Studies, Colorado

More information

Data Acquisition. Chapter 2

Data Acquisition. Chapter 2 Data Acquisition Chapter 2 1 st step: get data Data Acquisition Usually data gathered by some geophysical device Most surveys are comprised of linear traverses or transects Typically constant data spacing

More information

EMIGMA V9.x Premium Series April 8, 2015

EMIGMA V9.x Premium Series April 8, 2015 EMIGMA V9.x Premium Series April 8, 2015 EMIGMA for Gravity EMIGMA for Gravity license is a comprehensive package that offers a wide array of processing, visualization and interpretation tools. The package

More information

A comparison between time domain and depth domain inversion to acoustic impedance Laurence Letki*, Kevin Darke, and Yan Araujo Borges, Schlumberger

A comparison between time domain and depth domain inversion to acoustic impedance Laurence Letki*, Kevin Darke, and Yan Araujo Borges, Schlumberger Laurence Letki*, Kevin Darke, and Yan Araujo Borges, Schlumberger Summary Geophysical reservoir characterization in a complex geologic environment remains a challenge. Conventional amplitude inversion

More information

Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited

Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears

More information

ANALYSIS OF GEOPHYSICAL POTENTIAL FIELDS A Digital Signal Processing Approach

ANALYSIS OF GEOPHYSICAL POTENTIAL FIELDS A Digital Signal Processing Approach ADVANCES IN EXPLORATION GEOPHYSICS 5 ANALYSIS OF GEOPHYSICAL POTENTIAL FIELDS A Digital Signal Processing Approach PRABHAKAR S. NAIDU Indian Institute of Science, Bangalore 560012, India AND M.P. MATHEW

More information

3D inversion of marine CSEM data: A feasibility study from the Shtokman gas field in the Barents Sea

3D inversion of marine CSEM data: A feasibility study from the Shtokman gas field in the Barents Sea 3D inversion of marine CSEM data: A feasibility study from the Shtokman gas field in the Barents Sea M. S. Zhdanov 1,2, M. Čuma 1,2, A. Gribenko 1,2, G. Wilson 2 and N. Black 2 1 The University of Utah,

More information

EOSC 454 Lab #3. 3D Magnetics. Date: Feb 3, Due: Feb 10, 2009.

EOSC 454 Lab #3. 3D Magnetics. Date: Feb 3, Due: Feb 10, 2009. EOSC 454 Lab #3 3D Magnetics Date: Feb 3, 2009. Due: Feb 10, 2009. 1 Introduction In this exercise you will perform both forward models and inversions of magnetic data. You will compare the response of

More information

The Immersed Interface Method

The Immersed Interface Method The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial

More information

2. Give an example of a non-constant function f(x, y) such that the average value of f over is 0.

2. Give an example of a non-constant function f(x, y) such that the average value of f over is 0. Midterm 3 Review Short Answer 2. Give an example of a non-constant function f(x, y) such that the average value of f over is 0. 3. Compute the Riemann sum for the double integral where for the given grid

More information

High Resolution Imaging by Wave Equation Reflectivity Inversion

High Resolution Imaging by Wave Equation Reflectivity Inversion High Resolution Imaging by Wave Equation Reflectivity Inversion A. Valenciano* (Petroleum Geo-Services), S. Lu (Petroleum Geo-Services), N. Chemingui (Petroleum Geo-Services) & J. Yang (University of Houston)

More information

APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA.

APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA. APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA. BY: ISAAC KUMA YEBOAH. Department of Engineering, Regent University

More information

Introduction to Geophysical Inversion

Introduction to Geophysical Inversion Introduction to Geophysical Inversion Goals Understand the non-uniqueness in geophysical interpretations Understand the concepts of inversion. Basic workflow for solving inversion problems. Some important

More information

Dynamic data integration and stochastic inversion of a two-dimensional confined aquifer

Dynamic data integration and stochastic inversion of a two-dimensional confined aquifer Hydrology Days 2013 Dynamic data integration and stochastic inversion of a two-dimensional confined aquifer Dongdong Wang 1 Ye Zhang 1 Juraj Irsa 1 Abstract. Much work has been done in developing and applying

More information

26257 Nonlinear Inverse Modeling of Magnetic Anomalies due to Thin Sheets and Cylinders Using Occam s Method

26257 Nonlinear Inverse Modeling of Magnetic Anomalies due to Thin Sheets and Cylinders Using Occam s Method 26257 Nonlinear Inverse Modeling of Anomalies due to Thin Sheets and Cylinders Using Occam s Method R. Ghanati* (University of Tehran, Insitute of Geophysics), H.A. Ghari (University of Tehran, Insitute

More information

Qualitative Depth Estimation by Differencing Upward Continuations

Qualitative Depth Estimation by Differencing Upward Continuations Qualitative Depth Estimation by Differencing Upward Continuations Jacobsen (1987) made a strong case for using upward continuation filtering as a method for separating causative sources from various depths.

More information

Objectives for Terrain Week

Objectives for Terrain Week The Terrain Correction Chuck Connor, Laura Connor Potential Fields Geophysics: Terrain Week for Terrain Week Learn about the terrain correction The inner terrain correction Learn about The outer terrain

More information

Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak

Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society

More information

We G High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation

We G High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation We G103 05 High-resolution Tomography Using Offsetdependent Picking and Inversion Conditioned by Image-guided Interpolation G. Hilburn* (TGS), Y. He (TGS), Z. Yan (TGS) & F. Sherrill (TGS) SUMMARY An approach

More information

G042 Subsalt Imaging Challenges - A Deepwater Imaging Analysis

G042 Subsalt Imaging Challenges - A Deepwater Imaging Analysis G042 Subsalt Imaging Challenges - A Deepwater Imaging Analysis M. Cogan* (WesternGeco), J. Gardner (WesternGeco) & N. Moldoveanu (WesternGeco) SUMMARY Upon completion of the final reverse-time migration

More information

Computation of the gravity gradient tensor due to topographic masses using tesseroids

Computation of the gravity gradient tensor due to topographic masses using tesseroids Computation of the gravity gradient tensor due to topographic masses using tesseroids Leonardo Uieda 1 Naomi Ussami 2 Carla F Braitenberg 3 1. Observatorio Nacional, Rio de Janeiro, Brazil 2. Universidade

More information

Least-squares Wave-Equation Migration for Broadband Imaging

Least-squares Wave-Equation Migration for Broadband Imaging Least-squares Wave-Equation Migration for Broadband Imaging S. Lu (Petroleum Geo-Services), X. Li (Petroleum Geo-Services), A. Valenciano (Petroleum Geo-Services), N. Chemingui* (Petroleum Geo-Services),

More information

Integral equation method for anisotropic inversion of towed streamer EM data: theory and application for the TWOP survey

Integral equation method for anisotropic inversion of towed streamer EM data: theory and application for the TWOP survey Integral equation method for anisotropic inversion of towed streamer EM data: theory and application for the TWOP survey Michael S. Zhdanov 1,2, Masashi Endo 1, Daeung Yoon 1,2, Johan Mattsson 3, and Jonathan

More information

P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies

P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies P257 Transform-domain Sparsity Regularization in Reconstruction of Channelized Facies. azemi* (University of Alberta) & H.R. Siahkoohi (University of Tehran) SUMMARY Petrophysical reservoir properties,

More information

Iterative Methods for Solving Linear Problems

Iterative Methods for Solving Linear Problems Iterative Methods for Solving Linear Problems When problems become too large (too many data points, too many model parameters), SVD and related approaches become impractical. Iterative Methods for Solving

More information

U043 3D Prestack Time Domain Full Waveform Inversion

U043 3D Prestack Time Domain Full Waveform Inversion U043 3D Prestack Time Domain Full Waveform Inversion D.V. Vigh* (WesternGeco), W.E.S. Starr (WesternGeco) & K.D. Kenneth Dingwall (WesternGeco) SUMMARY Despite the relatively high computational demand,

More information

Large-scale workflows for wave-equation based inversion in Julia

Large-scale workflows for wave-equation based inversion in Julia Large-scale workflows for wave-equation based inversion in Julia Philipp A. Witte, Mathias Louboutin and Felix J. Herrmann SLIM University of British Columbia Motivation Use Geophysics to understand the

More information

Noise in FTG Data and its Comparison With Conventional Gravity Surveys

Noise in FTG Data and its Comparison With Conventional Gravity Surveys Noise in FTG Data and its Comparison With Conventional Gravity Surveys G. Barnes, J. umley, P. Houghton, R. Gleave ARKeX td. Cambridge, U.K. Summary Employing various assumptions and approximations, we

More information

Estimation of Depth and Shape Factor of Buried Structure From Residual Gravity Anomaly Data

Estimation of Depth and Shape Factor of Buried Structure From Residual Gravity Anomaly Data Australian Journal of Basic and Applied Sciences, 5(): 2-25, 2 ISSN 99-878 Estimation of Depth and Shape Factor of Buried Structure From Residual Gravity Anomaly Data MojtabaBabaee, Ahmad Alvandi, 2 Hosein

More information

Sparse wavelet expansions for seismic tomography: Methods and algorithms

Sparse wavelet expansions for seismic tomography: Methods and algorithms Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves 24 28 July 2011 (Joint

More information

Digital Image Processing. Prof. P. K. Biswas. Department of Electronic & Electrical Communication Engineering

Digital Image Processing. Prof. P. K. Biswas. Department of Electronic & Electrical Communication Engineering Digital Image Processing Prof. P. K. Biswas Department of Electronic & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Lecture - 21 Image Enhancement Frequency Domain Processing

More information

Inversion after depth imaging

Inversion after depth imaging Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.

More information

EE795: Computer Vision and Intelligent Systems

EE795: Computer Vision and Intelligent Systems EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 FDH 204 Lecture 14 130307 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Stereo Dense Motion Estimation Translational

More information

Engineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering

Engineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering Engineering Effects of Boundary Conditions (Fixtures and Temperatures) J.E. Akin, Rice University, Mechanical Engineering Here SolidWorks stress simulation tutorials will be re-visited to show how they

More information

F020 Methods for Computing Angle Gathers Using RTM

F020 Methods for Computing Angle Gathers Using RTM F020 Methods for Computing Angle Gathers Using RTM M. Vyas* (WesternGeco), X. Du (WesternGeco), E. Mobley (WesternGeco) & R. Fletcher (WesternGeco) SUMMARY Different techniques can be used to compute angle-domain

More information

Question: What are the origins of the forces of magnetism (how are they produced/ generated)?

Question: What are the origins of the forces of magnetism (how are they produced/ generated)? This is an additional material to the one in the internet and may help you to develop interest with the method. You should try to integrate some of the discussions here while you are trying to answer the

More information

The determination of the correct

The determination of the correct SPECIAL High-performance SECTION: H i gh-performance computing computing MARK NOBLE, Mines ParisTech PHILIPPE THIERRY, Intel CEDRIC TAILLANDIER, CGGVeritas (formerly Mines ParisTech) HENRI CALANDRA, Total

More information

PGS hyperbeam - rapid scenario-testing of velocity models to optimize depth imaging

PGS hyperbeam - rapid scenario-testing of velocity models to optimize depth imaging A Publication of Petroleum Geo-Services Vol. 10 No. 4 April 2010 PGS hyperbeam - rapid scenario-testing of velocity models to optimize depth imaging Introduction Depth imaging is now established as a key

More information

An introduction to interpolation and splines

An introduction to interpolation and splines An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve

More information

Dynamic Geometry Processing

Dynamic Geometry Processing Dynamic Geometry Processing EG 2012 Tutorial Will Chang, Hao Li, Niloy Mitra, Mark Pauly, Michael Wand Tutorial: Dynamic Geometry Processing 1 Articulated Global Registration Introduction and Overview

More information

Oasis montaj How-To Guide. VOXI Earth Modelling - Running an AGG Unconstrained Inversion

Oasis montaj How-To Guide. VOXI Earth Modelling - Running an AGG Unconstrained Inversion Oasis montaj How-To Guide VOXI Earth Modelling - Running an AGG Unconstrained Inversion The software described in this manual is furnished under license and may only be used or copied in accordance with

More information

Joint seismic traveltime and TEM inversion for near surface imaging Jide Nosakare Ogunbo*, Jie Zhang, GeoTomo LLC

Joint seismic traveltime and TEM inversion for near surface imaging Jide Nosakare Ogunbo*, Jie Zhang, GeoTomo LLC Jide Nosaare Ogunbo*, Jie Zhang, GeoTomo LLC Summary For a reliable interpretation of the subsurface structure, the joint geophysical inversion approach is becoming a viable tool. Seismic and EM methods

More information

Directions in 3-D imaging - Strike, dip, both?

Directions in 3-D imaging - Strike, dip, both? Stanford Exploration Project, Report 113, July 8, 2003, pages 363 369 Short Note Directions in 3-D imaging - Strike, dip, both? Marie L. Clapp 1 INTRODUCTION In an ideal world, a 3-D seismic survey would

More information

Progress Report on: Interferometric Interpolation of 3D SSP Data

Progress Report on: Interferometric Interpolation of 3D SSP Data Progress Report on: Interferometric Interpolation of 3D SSP Data Sherif M. Hanafy ABSTRACT We present the theory and numerical results for interferometrically interpolating and extrapolating 3D marine

More information

GeoProbe Geophysical Interpretation Software

GeoProbe Geophysical Interpretation Software DATA SHEET GeoProbe Geophysical Interpretation Software overview DecisionSpace Geosciences key features Integrated building, editing and interactive deformation of sealed multi-z bodies extracted from

More information

3D pyramid interpolation

3D pyramid interpolation 3D pyramid interpolation Xukai Shen ABSTRACT Seismic data geometries are not always as nice and regular as we want due to various acquisition constraints. In such cases, data interpolation becomes necessary.

More information

F031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data

F031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data F031 Application of Time Domain and Single Frequency Waveform Inversion to Real Data D Yingst (ION Geophysical), C. Wang* (ION Geophysical), J. Park (ION Geophysical), R. Bloor (ION Geophysical), J. Leveille

More information

Modern Methods of Data Analysis - WS 07/08

Modern Methods of Data Analysis - WS 07/08 Modern Methods of Data Analysis Lecture XV (04.02.08) Contents: Function Minimization (see E. Lohrmann & V. Blobel) Optimization Problem Set of n independent variables Sometimes in addition some constraints

More information

Magnetics. Introduction to Filtering using the Fourier Transform. Chuck Connor, Laura Connor. Potential Fields. Magnetics.

Magnetics. Introduction to Filtering using the Fourier Transform. Chuck Connor, Laura Connor. Potential Fields. Magnetics. Introduction to using the Chuck Connor, Laura Connor Potential Fields for this week Nm Sm Schematic Earth dipolar magnetic field. The field lines placed in the page plane are drawn as thick lines, those

More information

3-D vertical cable processing using EOM

3-D vertical cable processing using EOM Carlos Rodriguez-Suarez, John C. Bancroft, Yong Xu and Robert R. Stewart ABSTRACT Three-dimensional seismic data using vertical cables was modeled and processed using equivalent offset migration (EOM),

More information

Simultaneous joint inversion of refracted and surface waves Simone Re *, Claudio Strobbia, Michele De Stefano and Massimo Virgilio - WesternGeco

Simultaneous joint inversion of refracted and surface waves Simone Re *, Claudio Strobbia, Michele De Stefano and Massimo Virgilio - WesternGeco Simultaneous joint inversion of refracted and surface waves Simone Re *, Claudio Strobbia, Michele De Stefano and Massimo Virgilio - WesternGeco Summary In this paper, we review the near-surface challenges

More information

MAPPING POISSON S RATIO OF UNCONSOLIDATED MATERIALS FROM A JOINT ANALYSIS OF SURFACE-WAVE AND REFRACTION EVENTS INTRODUCTION

MAPPING POISSON S RATIO OF UNCONSOLIDATED MATERIALS FROM A JOINT ANALYSIS OF SURFACE-WAVE AND REFRACTION EVENTS INTRODUCTION MAPPING POISSON S RATIO OF UNCONSOLIDATED MATERIALS FROM A JOINT ANALYSIS OF SURFACE-WAVE AND REFRACTION EVENTS Julian Ivanov, Choon B. Park, Richard D. Miller, and Jianghai Xia Kansas Geological Survey

More information

Potent can model multi-component data, such as down-hole Fx, Fy, Fz magnetic data, or tensor gravity data.

Potent can model multi-component data, such as down-hole Fx, Fy, Fz magnetic data, or tensor gravity data. Getting Started What is Potent? Main concepts A simple tutorial What is Potent? Potent is a program for modelling the magnetic and gravitational effects of subsurface structures. It provides a highly interactive

More information

Provides New Tech Solution For Seismic Velocity Modeling

Provides New Tech Solution For Seismic Velocity Modeling JNURY 2009 The etter usiness Publication Serving the Exploration / Drilling / Production Industry Provides New Tech Solution For Seismic Velocity Modeling y Moritz M. Fliedner & Dimitri evc SNT CLR, C.

More information

Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples

Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples A. Ordoñez* (PGS), W.F. Sollner (PGS), T. Klüver (PGS) & L.G. Gelius (UiO) SUMMARY A joint

More information

Issues During the Inversion of Crosshole Radar Data: Can We Have Confidence in the Outcome?

Issues During the Inversion of Crosshole Radar Data: Can We Have Confidence in the Outcome? Boise State University ScholarWorks Geosciences Faculty Publications and Presentations Department of Geosciences 12-1-2006 Issues During the Inversion of Crosshole Radar Data: Can We Have Confidence in

More information

GeOlympus. GeoPACS GeoTR GeoStaR GeoWZ

GeOlympus. GeoPACS GeoTR GeoStaR GeoWZ GeOlympus GeoPACS GeoTR GeoStaR GeoWZ GeOlympus Workstation Work station Reflection time Refracted waves first arrivals Additional information Static and NMO corrections Layer by layer delay and statics

More information

H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography

H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography M. Cavalca* (Schlumberger - Westerngeco) & R.P. Fletcher (Schlumberger - Westerngeco) SUMMARY Estimation of the

More information

3D modeling of the Quest Projects Geophysical Datasets. Nigel Phillips

3D modeling of the Quest Projects Geophysical Datasets. Nigel Phillips 3D modeling of the Quest Projects Geophysical Datasets Nigel Phillips Advanced Geophysical Interpretation Centre Undercover Exploration workshop KEG-25 April 2012 Mineral Physical Properties: density sus.

More information

Downloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at Recovering the Reflectivity Matrix and Angle-dependent Plane-wave Reflection Coefficients from Imaging of Multiples Alba Ordoñez PGS/UiO*, Walter Söllner PGS, Tilman Klüver PGS and Leiv J. Gelius UiO Summary

More information

Multi-azimuth velocity estimation

Multi-azimuth velocity estimation Stanford Exploration Project, Report 84, May 9, 2001, pages 1 87 Multi-azimuth velocity estimation Robert G. Clapp and Biondo Biondi 1 ABSTRACT It is well known that the inverse problem of estimating interval

More information

GRAPRISM. Version Markku Pirttijärvi

GRAPRISM. Version Markku Pirttijärvi GRAPRISM Version 1.1 2003 Markku Pirttijärvi Introduction: Gravity, the attraction force between masses, is one of the basic forces of the physical world. The objective of geophysical gravity method is

More information

Introduction to Digital Image Processing

Introduction to Digital Image Processing Fall 2005 Image Enhancement in the Spatial Domain: Histograms, Arithmetic/Logic Operators, Basics of Spatial Filtering, Smoothing Spatial Filters Tuesday, February 7 2006, Overview (1): Before We Begin

More information

Gravity Gradients in PreStack Depth Migration

Gravity Gradients in PreStack Depth Migration Index Table of contents Gravity Gradients in PreStack Depth Migration Ed. K. Biegert Shell International Exploration and Production Summary Subsurface de-risking is presently almost exclusively done by

More information

WinGLink. Geophysical interpretation software VERSION

WinGLink. Geophysical interpretation software VERSION WinGLink Geophysical interpretation software VERSION 2.21.08 WinGLink geophysical interpretation software WinGLink Geophysical interpretation software Essential tools for data QC, processing, and modeling

More information

Archbold Area Schools Math Curriculum Map

Archbold Area Schools Math Curriculum Map Math 8 August - May Mathematical Processes Formulate a problem or mathematical model in response to a specific need or situation, determine information required to solve the problem, choose method for

More information

Source Estimation for Wavefield Reconstruction Inversion

Source Estimation for Wavefield Reconstruction Inversion Source Estimation for Wavefield Reconstruction Inversion Zhilong Fang * and Felix J. Herrmann * * Seismic Laboratory for Imaging and Modeling SLIM), University of British Columbia Abstract Wavefield reconstruction

More information

Stochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston

Stochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston Stochastic conjugate gradient method for least-square seismic inversion problems Wei Huang*, Hua-Wei Zhou, University of Houston Summary With the development of computational power, there has been an increased

More information

We N Depth Domain Inversion Case Study in Complex Subsalt Area

We N Depth Domain Inversion Case Study in Complex Subsalt Area We N104 12 Depth Domain Inversion Case Study in Complex Subsalt Area L.P. Letki* (Schlumberger), J. Tang (Schlumberger) & X. Du (Schlumberger) SUMMARY Geophysical reservoir characterisation in a complex

More information

Downloaded 09/16/13 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 09/16/13 to Redistribution subject to SEG license or copyright; see Terms of Use at Time-domain incomplete Gauss-Newton full-waveform inversion of Gulf of Mexico data Abdullah AlTheyab*, Xin Wang, Gerard T. Schuster, King Abdullah University of Science and Technology Downloaded 9// to

More information

Chapter 3 Path Optimization

Chapter 3 Path Optimization Chapter 3 Path Optimization Background information on optimization is discussed in this chapter, along with the inequality constraints that are used for the problem. Additionally, the MATLAB program for

More information

Graphical Analysis TEACHER NOTES SCIENCE NSPIRED. Science Objectives. Vocabulary. About the Lesson. TI-Nspire Navigator. Activity Materials

Graphical Analysis TEACHER NOTES SCIENCE NSPIRED. Science Objectives. Vocabulary. About the Lesson. TI-Nspire Navigator. Activity Materials Science Objectives Students will interpret a graph. Students will linearize data to find the mathematical relationship between two variables. Vocabulary data directly proportional inverse inversely proportional

More information

Unit Maps: Grade 8 Math

Unit Maps: Grade 8 Math Real Number Relationships 8.3 Number and operations. The student represents and use real numbers in a variety of forms. Representation of Real Numbers 8.3A extend previous knowledge of sets and subsets

More information

1.2 Numerical Solutions of Flow Problems

1.2 Numerical Solutions of Flow Problems 1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian

More information

EULER DECONVOLUTION OF POTENTIAL FIELD MAGNETOMETER AND ITS USE IN DELINEATION OF REGIONAL FAULTS WITH DATA FROM NORTHERN PERU

EULER DECONVOLUTION OF POTENTIAL FIELD MAGNETOMETER AND ITS USE IN DELINEATION OF REGIONAL FAULTS WITH DATA FROM NORTHERN PERU Euler Deconvolution 1 EULER DECONVOLUTION OF POTENTIAL FIELD MAGNETOMETER AND ITS USE IN DELINEATION OF REGIONAL FAULTS WITH DATA FROM NORTHERN PERU C.J. Solano Salinas a,b and E. Borda c a Centro de Tecnologias

More information

Equivalence of source-receiver migration and shot-profile migration

Equivalence of source-receiver migration and shot-profile migration Stanford Exploration Project, Report 112, November 11, 2002, pages 109 117 Short Note Equivalence of source-receiver migration and shot-profile migration Biondo Biondi 1 INTRODUCTION At first glance, shot

More information

Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi

Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi Image Transformation Techniques Dr. Rajeev Srivastava Dept. of Computer Engineering, ITBHU, Varanasi 1. Introduction The choice of a particular transform in a given application depends on the amount of

More information

G017 Beyond WAZ - A Modeling-based Evaluation of Extensions to Current Wide Azimuth Streamer Acquisition Geometries

G017 Beyond WAZ - A Modeling-based Evaluation of Extensions to Current Wide Azimuth Streamer Acquisition Geometries G017 Beyond WAZ - A Modeling-based Evaluation of Extensions to Current Wide Azimuth Streamer Acquisition Geometries M. Cvetkovic* (ION Geophysical), Z. Zhou (ION Geophysical / GXT Imaging Solutions) &

More information

CoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction

CoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction CoE4TN4 Image Processing Chapter 5 Image Restoration and Reconstruction Image Restoration Similar to image enhancement, the ultimate goal of restoration techniques is to improve an image Restoration: a

More information

Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method

Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,

More information

Variogram Inversion and Uncertainty Using Dynamic Data. Simultaneouos Inversion with Variogram Updating

Variogram Inversion and Uncertainty Using Dynamic Data. Simultaneouos Inversion with Variogram Updating Variogram Inversion and Uncertainty Using Dynamic Data Z. A. Reza (zreza@ualberta.ca) and C. V. Deutsch (cdeutsch@civil.ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta

More information

Geostatistics 2D GMS 7.0 TUTORIALS. 1 Introduction. 1.1 Contents

Geostatistics 2D GMS 7.0 TUTORIALS. 1 Introduction. 1.1 Contents GMS 7.0 TUTORIALS 1 Introduction Two-dimensional geostatistics (interpolation) can be performed in GMS using the 2D Scatter Point module. The module is used to interpolate from sets of 2D scatter points

More information

Determination of 2D shallow S wave velocity profile using waveform inversion of P-SV refraction data

Determination of 2D shallow S wave velocity profile using waveform inversion of P-SV refraction data Determination of 2D shallow S wave velocity profile using waveform inversion of P-SV refraction data Mohamed Amrouche and Hiroaki Yamanaka Dept. of Environmental Science and Technology, Tokyo Institute

More information

product description product capabilities

product description product capabilities product description Geosoft s Airborne Geophy s- ics application for the Oasis montaj software platform provides field geophysicists with the ability to process, filter, grid, and map data from airborne

More information

Anisotropy-preserving 5D interpolation by hybrid Fourier transform

Anisotropy-preserving 5D interpolation by hybrid Fourier transform Anisotropy-preserving 5D interpolation by hybrid Fourier transform Juefu Wang and Shaowu Wang, CGG Summary We present an anisotropy-preserving interpolation method based on a hybrid 5D Fourier transform,

More information

Anisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS

Anisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS Anisotropic model building with well control Chaoguang Zhou*, Zijian Liu, N. D. Whitmore, and Samuel Brown, PGS Summary Anisotropic depth model building using surface seismic data alone is non-unique and

More information

2D Inversions of 3D Marine CSEM Data Hung-Wen Tseng*, Lucy MacGregor, and Rolf V. Ackermann, Rock Solid Images, Inc.

2D Inversions of 3D Marine CSEM Data Hung-Wen Tseng*, Lucy MacGregor, and Rolf V. Ackermann, Rock Solid Images, Inc. 2D Inversions of 3D Marine CSEM Data Hung-Wen Tseng*, Lucy MacGregor, and Rolf V. Ackermann, Rock Solid Images, Inc. Summary A combination of 3D forward simulations and 2D and 3D inversions have been used

More information

GPR DATA ANALYSIS OF AIRFIELD RUNWAY USING MIGRATION VELOCITY SIMULATION ALGORITHM

GPR DATA ANALYSIS OF AIRFIELD RUNWAY USING MIGRATION VELOCITY SIMULATION ALGORITHM GPR DATA ANALYSIS OF AIRFIELD RUNWAY USING MIGRATION VELOCITY SIMULATION ALGORITHM Ramya M and Krishnan Balasubramaniam Centre for Nondestructive Evaluation, Department of Mechanical Engineering, Indian

More information