Downloaded 18 Jul 2011 to Redistribution subject to SEG license or copyright; see Terms of Use at
|
|
- Sherman Leon Warner
- 6 years ago
- Views:
Transcription
1 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 School of Mines SUMMARY We present a method for modeling the terrain response in gravity and gravity gradiometry surveys utilizing an adaptive quadtree mesh discretization. The data-dependent method is tailored to provide rapid and accurate spatial terrain corrections for ground or draped airborne surveys. The surface used in the modeling of terrain effect at each datum is discretized automatically to the largest cell size that will yield the desired accuracy, resulting in much faster modeling than fullresolution calculations. The largest cell sizes within the model occur in areas of minimal correction and at large distances away from the data location. We show synthetic and field examples for proof of concept. The adaptive quadtree method reduces the computational cost of the field example by performing 99.7% fewer calculations than the full model would require while retaining an accuracy of one Eötvös for the gradient data, and runs in.9% of the time for the gravity field terrain corrections compared with a commercial software package. INTRODUCTION The effect of terrain in gravity gradiometry surveys for petroleum and mining applications is often the predominant signal measured, as the field decays as distance-cubed (Chen and Macnae, 1997). Therefore accurate terrain corrections are critical to proper interpretation of the dataset. However, the accuracy of the terrain correction is most often a function of the resolution of the available digital elevation model (DEM). As a result, higher and higher resolution DEMs are acquired during gravity gradient surveys. While this results in more accurate terrain corrections, these high-resolution DEMs have the unfortunate consequence of creating extremely computationally intensive terrain corrections. Current methods of terrain correction often forward model the terrain effect at full resolution across the dataset, taking days and sometimes weeks to complete. Much faster algorithms are needed if field interpretation of data is required or if results are simply needed quickly. Gravity data also require high resolution DEMs, especially for ground-based surveys. In addition, a much larger spatial extent of the DEM is required as the field decays as distance-squared. While the minimum DEM extent and cell size for a groundbased gravity survey is well defined (Hammer, 1939), airborne gravity terrain corrections are more troublesome and less defined (Hammer, 1974). While extremely high resolution digital elevation models are required in areas with high relief and/or near observation locations, the resolution requirements can be relaxed as the distance from the observation point increases. Also, if the characteristics of the terrain change over the survey, the resolution requirements of the DEM may change as well. Moreover, the required resolution is also a function of the flight height (Reid, 198; Kass and Li, 28). Thus a data-adaptive discretization of the DEM can result in faster computations and smaller storage requirements while ensuring desired accuracy. In this abstract, we show that the adaptive quadtree method is appropriate for data-dependent discretization of DEMs for terrain corrections. The method allows for user-defined error bounds and parallelization, and reduces the number of required calculations by three to five orders of magnitude. We start with the forward modeling calculation of the prismatic model and then discuss the discretization method. We show a synthetic case for an illustration as well as a field example. FORWARD MODEL The terrain effect in gravity gradiometry is calculated as the response due to the prism contribution bounded above by the terrain and below by a plane that passes through the lowest terrain point in the survey area. The forward model of the gravity gradient follows a similar construction as described in Li (21). As gravitational potential satisfies Laplace s equation, the potentials due to each source superimpose at each observation point and thus the contribution of all prisms to the gravity gradient measured at a point is the sum of their individual responses. For a volume source at location r with density ρ, the gravity gradient at an observation point (denoted by vector r ) is a tensor T with individual components represented by T kl = γρ 1 dv V x k x l r r o (1) where V is the volume of the source. The integration follows the magnetic derivation by Sharma (1966) or Zhang et al. (2), with treatment to avoid singularities at observation locations coplanar with the faces of the prisms. As the equations are lengthy and not required for understanding of the method, we direct the reader to Sharma s paper for a detailed discussion of the discretization of the integral equations. For source-free regions, this tensor is symmetric and traceless, and commonly is written in the form: [ Txx T xy T xz ] T = T yy T yz T zz Here we have adopted a Cartesian coordinate system having its origin at the Earth s surface, x-axis pointing towards grid north, y-axis pointing towards grid east, and z-axis pointing vertically downward. Thus we can write the expression for the gravity gradient at one data point, d i, as d i T kl (r oi ) { M } 1 = γρ j V j x k x l r r oi dv (2) 1789
2 M ρ j G i j. In this construction, we are summing the response due to many prisms V j with constant density ρ j within the prism. In a compact notation for a linear system, we can show Equation 2 as Gρ = d, (3) where d is the data response from the prisms based on their density, ρ. It is important to note that the matrix, G, can be quite large. However, the forward modeling algorithm does not require the formation of the complete matrix at one time, but rather only the computation for one datum at a time. The forward model of the gravity field follows a similar approach. The field due to a volume source is represented as z z o g z (r o ) = γρ dv, (4) V r r o 3 leading to the gravity field at one data point d i due to the sum of responses: d i = { M } z z oi γρ j V j r r oi 3 dv (5) M ρ j G i j. DISCRETIZATION APPROACH In order to reduce the number of computations for the terrain correction, we employ an adaptive quadtree mesh discretization of the supplied digital elevation model. This results in fewer prisms for which to calculate the gravity gradient response. The derivation of the discretization process for the gravity field is identical to that of the gravity gradient. The adaptive quadtree mesh design (Frey and Marechal, 1998) is one that places larger cells where no or little signal is present (in this case response from topographic relief) and smaller cells to increase resolution where sources are located. The quadtree method is useful in large scale problems by reducing the total number of model cells (Ascher and Haber, 21). Quadtree mesh discretization has been used primarily in geophysics in remote sensing applications (e.g. Gerstner, 1999), but recently has found application in use with mesh padding for DC resistivity (Eso and Oldenburg, 27) and magnetic processing (Davis and Li, 27). The mesh design is adaptive due to the data dependence from dataset to dataset as well as region to region within the same data set, as opposed to a quadtree mesh pre-defined by a known function or distribution. In general, the mesh starts with cells of a large, user-defined width; if a property has a value higher than a given threshold, it is split into quarters. This process continues until all cell properties are within the threshold. The quadtree mesh threshold is chosen as a function of a userdefined error level (usually at or below the aggregate survey noise). The error is defined by the difference between the examined large prism and the smaller prisms (e.g. higher resolution) within the same area. A quadtree algorithm is then used to discretize the terrain for each data point. At each stage, the binned prisms have a height of average terrain within the respective cell boundaries. The gravity gradient response of each cell is calculated and then the cell is split into four. The terrain effect of each new cell with a new average height is calculated again (five total forward calculations). The four cells are kept if the difference between the gravity gradient components of the four cells is more than the error threshold. Those are then examined and split again, if necessary. As the process approaches the optimal resolution, splitting the cells will provide less and less benefit until further splitting no longer has an effect greater than the accuracy threshold or the split cell size is equal to the DEM resolution. The error level for the quadtree algorithm is normalized by the number of cells in the model at each iteration. Thus if the maximum error is allowed in each cell it will superimpose and the upper bound of the error between the full and quadtree model is the defined error level. Then, if a cell is split (into four smaller cells), the error threshold for that cell is divided by four in order to preserve the upper bound. For example, if the initial model is split into 256 cells, then the first error level is the user-defined error divided by 256. Though the error is an upper bound, in practice the true error is much smaller as shown in the examples. Selection of the initial discretization has no effect on final accuracy, but is critical for the overall computation speed. In order to quantitatively determine a data-dependent initial discretization we examine the largest radial wavelength of the digital elevation model. The radial wavelength is used because of the nature of the two-dimensional DEM. Moreover, this simple step is fast to calculate because it depends on the number of data points in each direction and their respective spacing and the Fourier transform is unnecessary. Once the largest wavelength is calculated, it becomes the width of the initial discretized prisms for modeling. The splitting of these cells into smaller ones then occurs to retain resolution of the DEM. Further discretization though the adaptive mesh is analogous to adding higher frequencies in the Fourier domain. SYNTHETIC EXAMPLE As an example for illustrating the accuracy of the method, we use the Gaussian hill topography model presented in Figure 1. The hill extends from meters at ground level to 3 meters above the ground at its peak. The total number of observations locations was 625. We calculated the gravity gradient response of this terrain at 35 meters above ground level for each component and allowed a 1-Eötvös error. The starting cell size is 16 m by 16 m based on the longest radial wavelength. We show the adaptive discretization of a Gaussian hill for a single observation location in Figure 1. The asymmetry seen of the discretization is a function of the cen- 179
3 meters Northing (m) Easting (m) Figure 1: Gaussian hill with height equal to 3 meters. Automated adaptive quadtree discretization of terrain to achieve a maximum of 1 Eo tvo s accuracy 5 meters immediately above the peak for all tensor components. ter of the initial cells not being coincident with the peak of the model. We compared the response of the Gaussian hill calculated with an adaptive quadtree discretized DEM (Figure 2) to the regularlygridded, surface-based forward modeling. Though we allowed Figure 2: Gravity gradient response for each component at the algorithm as much as 1-Eo tvo s error, the maximum differa plane 5 meters above the maximum height of a 3 meterence between the two datasets was.17 Eo tvo s (Figure 2). high Gaussian hill. Units in Eo tvo s. Difference between We next examine the number of calculations performed in orthe terrain effect calculated with the adaptive quadtree method der to evaluate the efficiency of the adaptive algorithm. The versus a full-resolution surface-based method. The maximum full model contains 729,88 cells which means that there are error threshold given to the adaptive mesh was 1 Eo tvo s. Units 455,68, calculations. Using the 1 Eo tvo s error level exare in 1/1 of an Eo tvo s (-1 to 1 Eo tvo s color bar). ample, our method made 1,468,486 calculations, 31 times less (or.3%) of the full model. This is an extreme case of efficiency because of the nature of the of the terrain model itself. FIELD EXAMPLE To demonstrate the efficiency of the algorithm, we calculated the terrain corrections for a draped airborne gravity gradient and gravity survey from an AGG test site in Australia. The gravity gradient terrain corrections used a density contrast of 2.67 g/cc and allowed a 1 Eo tvo s error while the gravity corrections used the same density contrast and 1 mgal error tolerance. Figure?? shows the discretization of the DEM for one observation location. The terrain corrections calculated in this manner had a compression ratio of approximately 351, Figure 3: Vertically-exaggerated terrain with survey locations shown in white. 1791
4 meaning there were approximately.3% of the calculations necessary to produce the corrections as compared to using fullresolution. Figure 4 shows the calculated Tzz component of the gravity gradient. The airborne terrain correction for the vertical component of the gravity field was calculated in the same manner. Figure?? shows the corrections which have the expected features; these appear smoothed due to the observation elevation (Davis et al., 28). In comparison with an industry-standard commercial software package, the terrain correction took less than.9% of the time to run. CONCLUSIONS Terrain calculations with adaptive DEM gridding are faster than full-resolution calculations and allow user-specified accuracy. Depending on the complexity and relative relief of the terrain model, the adaptive quadtree correction method can reduce the number of model cell calculations by at least a magnitude which, depending on the problem, results in significant computational savings. The algorithm begins with a starting model with cell discretized the largest wavelength of observed terrain. It then proceeds to split cells to increase resolution and accuracy until a pre-defined error threshold is achieved. ACKNOWLEDGMENTS The authors would like to thank the Korea Resources Corporation and the sponsors of the Gravity and Magnetics Research Consortium: Anadarko, BGP, BP, ConocoPhillips, Marathon, and Vale for funding support. We thank Theo Aravanis and Rio Tinto for providing the Lidar data, which were acquired by Fugro Spatial. We also thank Trevor Irons for helpful discussions. (c) Figure 4: Discretization for one gravity gradient observation location for the calculation of terrain effect. Calculated Tzz component of the calculated terrain correction. (c) Calculated terrain correction for the vertical component of the gravity field. Units in mgal. 1792
5 EDITED REFERENCES Note: This reference list is a copy-edited version of the reference list submitted by the author. Reference lists for the 21 SEG Technical Program Expanded Abstracts have been copy edited so that references provided with the online metadata for each paper will achieve a high degree of linking to cited sources that appear on the Web. REFERENCES Ascher, U. M., and E. Haber, 21, Grid refinement and scaling for distributed parameter estimation problems: Inverse Problems, 17, no. 3, , doi:1.188/ /17/3/314. Chen, J., and J. Macnae, 1997, Terrain corrections are critical for airborne gravity gradiometer data: Exploration Geophysics, 28, no. 2, 21 25, doi:1.171/eg Davis, K., M. Kass, R. Krahenbuhl, and Y. Li, 28, Survey design and model appraisal based on resolution analysis for 4D gravity monitoring: Presented at the 78th Annual International Meeting, SEG. Davis, K., and Y. Li, 27, A fast approach to magnetic equivalent source processing using an adaptive quadtree mesh discretization: Presented at the 19th Annual Conference and Exhibition, Australian Society of Exploration Geophysicists. Eso, R., and D. Oldenburg, 27, Efficient 2.5D resistivity modeling using a quadtree discretization: Conference Proceedings, Symposium of Applied Geophysics on Engineering and Environmental Problems, Frey, P. J., and L. Marechal, 1998, Fast adaptive quadtree mesh generation: 7th International Meshing Roundtable, Sandia National Laboratories, Gerstner, T., 1999, Adaptive hierarchical methods for landscape representation and analysis : Lecture Notes in Earth Sciences, 78, 75 92, doi:1.17/bfb972. Hammer, S., 1939, Terrain corrections for gravimeter stations: Geophysics, 4, 184., 1974, Topographic and terrain correction for airborne gravity: Geophysics, 39, 537. Kass, M., and Y. Li, 28, Practical aspects of terrain correction in airborne gravity gradiometry surveys: Exploration Geophysics, 39, no. 4, , doi:1.171/eg821. Li, Y., 21, 3D inversion of gravity gradiometry data: Presented at the 71st Annual International Meeting, SEG. Reid, A. B., 198, Aeromagnetic survey design: Geophysics, 45, , doi:1.119/ Sharma, P. V., 1966, Rapid computation of magnetic anomalies and demagnetization effects caused by bodies of arbitrary shape: Pure and Applied Geophysics, 64, no. 1, 89 19, doi:1.17/bf Zhang, C., M. F. Mushayandebvu, A. B. Reid, J. D. Fairhead, and M. E. Odegard, 2, Euler deconvolution of gravity tensor gradient data : Geophysics, 65, , doi:1.119/
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 informationEfficient 3D Gravity and Magnetic Modeling
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
More informationDownloaded 23 May 2011 to Redistribution subject to SEG license or copyright; see Terms of Use at
Efficient 3D inversion of magnetic data via octree mesh discretization, space-filling curves, and wavelets Kristofer Davis and Yaoguo Li, Center for Gravity, Electrical, and Magnetics, Colorado School
More informationRapid gravity and gravity gradiometry terrain corrections via an adaptive quadtree mesh discretization
CSIRO PUBLISHING www.publish.csiro.au/journals/eg Copublished paper, please cite all three: Exploration Geophysics, 211, 42, 88 97; Butsuri-Tansa, 211, 64, 88 97; Jigu-Mulli-wa-Mulli-Tasa, 211, 14, 88
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
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 informationGEOPHYS 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 informationComputation 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 information3D 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 informationNoise 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 informationObjectives 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 informationApplication 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 informationEMIGMA 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 informationIntroduction 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 informationOasis 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 informationUnstructured grid modelling to create 3-D Earth models that unify geological and geophysical information
Unstructured grid modelling to create 3-D Earth models that unify geological and geophysical information Peter Lelièvre, Angela Carter-McAuslan, Cassandra Tycholiz, Colin Farquharson and Charles Hurich
More information3D Inversion of Time-Domain Electromagnetic Data for Ground Water Aquifers
3D Inversion of Time-Domain Electromagnetic Data for Ground Water Aquifers Elliot M. Holtham 1, Mike McMillan 1 and Eldad Haber 2 (1) Computational Geosciences Inc. (2) University of British Columbia Summary
More informationElectromagnetic 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 informationNoise and repeatability of airborne gravity gradiometry
Noise and repeatability of airborne gravity gradiometry Asbjorn Norlund Christensen 1*, Mark H. Dransfield 2 and Christopher Van Galder 2 present various noise estimates derived from surveys over the R.J.
More informationSummary. Introduction
. Tony Martin*, Cristiano Saturni and Peter Ashby, ION Geophysical Summary Modern marine seismic surveys may contain many Terabytes of data. In seismic processing terms, understanding the impact and effectiveness
More informationHow Deep Can My Magnetometer See?
A common question when using magnetometers or gradiometers is, How deep is my instrumentation seeing? This Magnetic Moment provides some answers to this question while emphasizing quick analysis methods
More informationDIGITAL TERRAIN MODELS
DIGITAL TERRAIN MODELS 1 Digital Terrain Models Dr. Mohsen Mostafa Hassan Badawy Remote Sensing Center GENERAL: A Digital Terrain Models (DTM) is defined as the digital representation of the spatial distribution
More informationNumerical simulation of 3-D seepage field in tailing pond and its practical application
Available online at www.sciencedirect.com Procedia Engineering (0) 70 76 0 SREE Conference on Engineering Modelling and Simulation Numerical simulation of -D seepage field in tailing pond and its practical
More informationMulticomponent land data pre-processing for FWI: a benchmark dataset
Multicomponent land data pre-processing for FWI: a benchmark dataset Raul Cova, Bernie K. Law and Kris Innanen CRWES/University of Calgary Summary Successful full-waveform inversion (FWI) studies using
More informationGeophysics 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 informationTensor Based Approaches for LVA Field Inference
Tensor Based Approaches for LVA Field Inference Maksuda Lillah and Jeff Boisvert The importance of locally varying anisotropy (LVA) in model construction can be significant; however, it is often ignored
More informationSUMMARY. solve the matrix system using iterative solvers. We use the MUMPS codes and distribute the computation over many different
Forward Modelling and Inversion of Multi-Source TEM Data D. W. Oldenburg 1, E. Haber 2, and R. Shekhtman 1 1 University of British Columbia, Department of Earth & Ocean Sciences 2 Emory University, Atlanta,
More informationGenerate Digital Elevation Models Using Laser Altimetry (LIDAR) Data. Christopher Weed
Generate Digital Elevation Models Using Laser Altimetry (LIDAR) Data Christopher Weed Final Report EE 381K Multidimensional Digital Signal Processing December 11, 2000 Abstract A Laser Altimetry (LIDAR)
More informationMagnetics. 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 informationNEW MONITORING TECHNIQUES ON THE DETERMINATION OF STRUCTURE DEFORMATIONS
Proceedings, 11 th FIG Symposium on Deformation Measurements, Santorini, Greece, 003. NEW MONITORING TECHNIQUES ON THE DETERMINATION OF STRUCTURE DEFORMATIONS D.Stathas, O.Arabatzi, S.Dogouris, G.Piniotis,
More information3. Data Structures for Image Analysis L AK S H M O U. E D U
3. Data Structures for Image Analysis L AK S H M AN @ O U. E D U Different formulations Can be advantageous to treat a spatial grid as a: Levelset Matrix Markov chain Topographic map Relational structure
More informationBiometrics Technology: Image Processing & Pattern Recognition (by Dr. Dickson Tong)
Biometrics Technology: Image Processing & Pattern Recognition (by Dr. Dickson Tong) References: [1] http://homepages.inf.ed.ac.uk/rbf/hipr2/index.htm [2] http://www.cs.wisc.edu/~dyer/cs540/notes/vision.html
More informationGRAV3D Version 3.0. A Program Library for Forward Modelling and Inversion of Gravity Data over 3D Structures.
GRAV3D Version 3.0 A Program Library for Forward Modelling and Inversion of Gravity Data over 3D Structures. UBC-Geophysical Inversion Facility Department of Earth and Ocean Sciences University of British
More informationGravity Methods (VII) wrap up
Environmental and Exploration Geophysics II Gravity Methods (VII) wrap up tom.h.wilson tom.wilson@mail.wvu.edu Department of Geology and Geography West Virginia University Morgantown, WV Items on the list
More informationQUEST Project: 3D inversion modelling, integration, and visualization of airborne gravity, magnetic, and electromagnetic data, BC, Canada.
Mira Geoscience Limited 409 Granville Street, Suite 512 B Vancouver, BC Canada V6C 1T2 Tel: (778) 329-0430 Fax: (778) 329-0668 info@mirageoscience.com www.mirageoscience.com QUEST Project: 3D inversion
More informationFOOTPRINTS EXTRACTION
Building Footprints Extraction of Dense Residential Areas from LiDAR data KyoHyouk Kim and Jie Shan Purdue University School of Civil Engineering 550 Stadium Mall Drive West Lafayette, IN 47907, USA {kim458,
More information2D 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 information26257 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 informationRobust 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 informationproduct 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 informationEXAM SOLUTIONS. Image Processing and Computer Vision Course 2D1421 Monday, 13 th of March 2006,
School of Computer Science and Communication, KTH Danica Kragic EXAM SOLUTIONS Image Processing and Computer Vision Course 2D1421 Monday, 13 th of March 2006, 14.00 19.00 Grade table 0-25 U 26-35 3 36-45
More informationEngineering and Environmental Geophysics with terratem
Engineering and Environmental Geophysics with terratem (1) Geophysical Consultant Alpha Geosciences Pty. Ltd. Sydney NSW 2210 Australia Email: rhenderson@terratem.com Roger Henderson (1) ABSTRACT terratem
More information[Youn *, 5(11): November 2018] ISSN DOI /zenodo Impact Factor
GLOBAL JOURNAL OF ENGINEERING SCIENCE AND RESEARCHES AUTOMATIC EXTRACTING DEM FROM DSM WITH CONSECUTIVE MORPHOLOGICAL FILTERING Junhee Youn *1 & Tae-Hoon Kim 2 *1,2 Korea Institute of Civil Engineering
More informationEOSC 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 information3D 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 informationEMIGMA Premium and Professional editions
EMIGMA Premium and Professional editions EMIGMA Overview EMIGMA Database Version Database Design Objectives Reduce user s time required for modelling and data analyses Enhance data analyses capabilities
More informationChapters 1-4: Summary
Chapters 1-4: Summary So far, we have been investigating the image acquisition process. Chapter 1: General introduction Chapter 2: Radiation source and properties Chapter 3: Radiation interaction with
More informationGenerate Digital Elevation Models Using Laser Altimetry (LIDAR) Data
Generate Digital Elevation Models Using Laser Altimetry (LIDAR) Data Literature Survey Christopher Weed October 2000 Abstract Laser altimetry (LIDAR) data must be processed to generate a digital elevation
More informationReality Check: Processing LiDAR Data. A story of data, more data and some more data
Reality Check: Processing LiDAR Data A story of data, more data and some more data Red River of the North Red River of the North Red River of the North Red River of the North Introduction and Background
More informationPresentation and analysis of multidimensional data sets
Presentation and analysis of multidimensional data sets Overview 1. 3D data visualisation Multidimensional images Data pre-processing Visualisation methods Multidimensional images wavelength time 3D image
More informationGEON Points2Grid Utility Instructions By: Christopher Crosby OpenTopography Facility, San Diego Supercomputer Center
GEON Points2Grid Utility Instructions By: Christopher Crosby (ccrosby@sdsc.edu) OpenTopography Facility, San Diego Supercomputer Center (Formerly: GEON / Active Tectonics Research Group School of Earth
More informationMasters Thesis. Effects of low-pass filtering on inversion of airborne gravity gradient data. Center for Gravity, Electrical & Magnetic Studies
CGEM Masters Thesis Colorado School of Mines Center for Gravity, Electrical & Magnetic Studies Effects of low-pass filtering on inversion of airborne gravity gradient data Jeongmin Lee Department of Geophysics
More informationDirectional Derivatives. Directional Derivatives. Directional Derivatives. Directional Derivatives. Directional Derivatives. Directional Derivatives
Recall that if z = f(x, y), then the partial derivatives f x and f y are defined as and represent the rates of change of z in the x- and y-directions, that is, in the directions of the unit vectors i and
More informationSchemes for improving efficiency of pixel-based inversion algorithms for electromagnetic loggingwhile-drilling
Schemes for improving efficiency of pixel-based inversion algorithms for electromagnetic loggingwhile-drilling measurements Y. Lin, A. Abubakar, T. M. Habashy, G. Pan, M. Li and V. Druskin, Schlumberger-Doll
More informationDECONFLICTION AND SURFACE GENERATION FROM BATHYMETRY DATA USING LR B- SPLINES
DECONFLICTION AND SURFACE GENERATION FROM BATHYMETRY DATA USING LR B- SPLINES IQMULUS WORKSHOP BERGEN, SEPTEMBER 21, 2016 Vibeke Skytt, SINTEF Jennifer Herbert, HR Wallingford The research leading to these
More informationTopographic Lidar Data Employed to Map, Preserve U.S. History
OCTOBER 11, 2016 Topographic Lidar Data Employed to Map, Preserve U.S. History In August 2015, the National Park Service (NPS) contracted Woolpert for the Little Bighorn National Monument Mapping Project
More informationPreparation for flight testing the VK1 gravity gradiometer
Preparation for flight testing the VK1 gravity gradiometer J. Anstie 1, T. Aravanis 2, P. Johnston 3, A. Mann 4, M. Longman 5, A. Sergeant 6, R. Smith 7, F. Van Kann 8, G. Walker 9, G. Wells 10, and J.
More informationRobotics Programming Laboratory
Chair of Software Engineering Robotics Programming Laboratory Bertrand Meyer Jiwon Shin Lecture 8: Robot Perception Perception http://pascallin.ecs.soton.ac.uk/challenges/voc/databases.html#caltech car
More informationSite map 1 of 1. Complete 3D mag/grav workflow outline. 1. Setting up for inversion. 1. Introduction
Site map 1 of 1 Complete 3D mag/grav workflow outline 0. Introduction Goals Using the workflow Best practice 4. Relevant assumptions 5. Codes needed 6. Examples Setting up for inversion Clarify the problem
More informationOasis montaj How-To Guide. VOXI Earth Modelling - Running an Inversion
Oasis montaj How-To Guide VOXI Earth Modelling - Running an Inversion The software described in this manual is furnished under license and may only be used or copied in accordance with the terms of the
More informationAdvances in Airborne Gravity and Magnetics
Advances in Airborne Gravity and Magnetics J Derek Fairhead* JD GeoConsultancy Ltd., & School of Earth and Environment, University of Leeds, UK Gordon R J Cooper, School of Geosciences, University of the
More informationAdaptive 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 informationSpatial Enhancement Definition
Spatial Enhancement Nickolas Faust The Electro- Optics, Environment, and Materials Laboratory Georgia Tech Research Institute Georgia Institute of Technology Definition Spectral enhancement relies on changing
More informationWinGLink. 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 informationSummary. Introduction
Leif H. Cox*, Montana Tech and TechnoImaging, Glenn A. Wilson, TechnoImaging, and Michael S. Zhdanov, The University of Utah and TechnoImaging Summary We show that it is practical to invert entire airborne
More informationApplied Cartography and Introduction to GIS GEOG 2017 EL. Lecture-7 Chapters 13 and 14
Applied Cartography and Introduction to GIS GEOG 2017 EL Lecture-7 Chapters 13 and 14 Data for Terrain Mapping and Analysis DEM (digital elevation model) and TIN (triangulated irregular network) are two
More informationEffect of suction pipe leaning angle and water level on the internal flow of pump sump
IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Effect of suction pipe leaning angle and water level on the internal flow of pump sump To cite this article: Z-M Chen et al 216
More informationA METHOD TO PREDICT ACCURACY OF LEAST SQUARES SURFACE MATCHING FOR AIRBORNE LASER SCANNING DATA SETS
A METHOD TO PREDICT ACCURACY OF LEAST SQUARES SURFACE MATCHING FOR AIRBORNE LASER SCANNING DATA SETS Robert Pâquet School of Engineering, University of Newcastle Callaghan, NSW 238, Australia (rpaquet@mail.newcastle.edu.au)
More informationAPPLICABILITY ANALYSIS OF CLOTH SIMULATION FILTERING ALGORITHM FOR MOBILE LIDAR POINT CLOUD
APPLICABILITY ANALYSIS OF CLOTH SIMULATION FILTERING ALGORITHM FOR MOBILE LIDAR POINT CLOUD Shangshu Cai 1,, Wuming Zhang 1,, Jianbo Qi 1,, Peng Wan 1,, Jie Shao 1,, Aojie Shen 1, 1 State Key Laboratory
More informationLab 12: Sampling and Interpolation
Lab 12: Sampling and Interpolation What You ll Learn: -Systematic and random sampling -Majority filtering -Stratified sampling -A few basic interpolation methods Data for the exercise are in the L12 subdirectory.
More informationRipplet: a New Transform for Feature Extraction and Image Representation
Ripplet: a New Transform for Feature Extraction and Image Representation Dr. Dapeng Oliver Wu Joint work with Jun Xu Department of Electrical and Computer Engineering University of Florida Outline Motivation
More informationLidar and GIS: Applications and Examples. Dan Hedges Clayton Crawford
Lidar and GIS: Applications and Examples Dan Hedges Clayton Crawford Outline Data structures, tools, and workflows Assessing lidar point coverage and sample density Creating raster DEMs and DSMs Data area
More informationANALYSIS 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 informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More informationSUMMARY. 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 information6. Object Identification L AK S H M O U. E D U
6. Object Identification L AK S H M AN @ O U. E D U Objects Information extracted from spatial grids often need to be associated with objects not just an individual pixel Group of pixels that form a real-world
More informationDynamic 3D representation of information using low cost Cloud ready Technologies
National Technical University Of Athens School of Rural and Surveying Engineering Laboratory of Photogrammetry Dynamic 3D representation of information using low cost Cloud ready Technologies George MOURAFETIS,
More informationiric Software Changing River Science River2D Tutorials
iric Software Changing River Science River2D Tutorials iric Software Changing River Science Confluence of the Colorado River, Blue River and Indian Creek, Colorado, USA 1 TUTORIAL 1: RIVER2D STEADY SOLUTION
More informationOasis montaj How-To Guide. VOXI Earth Modelling - Running an Inversion Using Gradient Weighting
Oasis montaj How-To Guide VOXI Earth Modelling - Running an Inversion Using Gradient Weighting The software described in this manual is furnished under license and may only be used or copied in accordance
More informationMETRIC PLANE RECTIFICATION USING SYMMETRIC VANISHING POINTS
METRIC PLANE RECTIFICATION USING SYMMETRIC VANISHING POINTS M. Lefler, H. Hel-Or Dept. of CS, University of Haifa, Israel Y. Hel-Or School of CS, IDC, Herzliya, Israel ABSTRACT Video analysis often requires
More informationIntroduction to Geophysical Modelling and Inversion
Introduction to Geophysical Modelling and Inversion James Reid GEOPHYSICAL INVERSION FOR MINERAL EXPLORERS ASEG-WA, SEPTEMBER 2014 @ 2014 Mira Geoscience Ltd. Forward modelling vs. inversion Forward Modelling:
More informationLecture 06. Raster and Vector Data Models. Part (1) Common Data Models. Raster. Vector. Points. Points. ( x,y ) Area. Area Line.
Lecture 06 Raster and Vector Data Models Part (1) 1 Common Data Models Vector Raster Y Points Points ( x,y ) Line Area Line Area 2 X 1 3 Raster uses a grid cell structure Vector is more like a drawn map
More informationSPECIAL TECHNIQUES-II
SPECIAL TECHNIQUES-II Lecture 19: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Method of Images for a spherical conductor Example :A dipole near aconducting sphere The
More informationWe N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions
We N106 02 Converted-phase Seismic Imaging - Amplitudebalancing -independent Imaging Conditions A.H. Shabelansky* (Massachusetts Institute of Technology), A.E. Malcolm (Memorial University of Newfoundland)
More informationMulticomponent f-x seismic random noise attenuation via vector autoregressive operators
Multicomponent f-x seismic random noise attenuation via vector autoregressive operators Mostafa Naghizadeh and Mauricio Sacchi ABSTRACT We propose an extension of the traditional frequency-space (f-x)
More informationMore About Factoring Trinomials
Section 6.3 More About Factoring Trinomials 239 83. x 2 17x 70 x 7 x 10 Width of rectangle: Length of rectangle: x 7 x 10 Width of shaded region: 7 Length of shaded region: x 10 x 10 Area of shaded region:
More informationMulti-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 informationAPPENDIX E2. Vernal Pool Watershed Mapping
APPENDIX E2 Vernal Pool Watershed Mapping MEMORANDUM To: U.S. Fish and Wildlife Service From: Tyler Friesen, Dudek Subject: SSHCP Vernal Pool Watershed Analysis Using LIDAR Data Date: February 6, 2014
More informationTHE VIEWING TRANSFORMATION
ECS 178 Course Notes THE VIEWING TRANSFORMATION Kenneth I. Joy Institute for Data Analysis and Visualization Department of Computer Science University of California, Davis Overview One of the most important
More informationDigital Image Processing
Digital Image Processing Part 9: Representation and Description AASS Learning Systems Lab, Dep. Teknik Room T1209 (Fr, 11-12 o'clock) achim.lilienthal@oru.se Course Book Chapter 11 2011-05-17 Contents
More informationIntroduction to Computer Graphics. Modeling (3) April 27, 2017 Kenshi Takayama
Introduction to Computer Graphics Modeling (3) April 27, 2017 Kenshi Takayama Solid modeling 2 Solid models Thin shapes represented by single polygons Unorientable Clear definition of inside & outside
More informationLecture 15: Segmentation (Edge Based, Hough Transform)
Lecture 15: Segmentation (Edge Based, Hough Transform) c Bryan S. Morse, Brigham Young University, 1998 000 Last modified on February 3, 000 at :00 PM Contents 15.1 Introduction..............................................
More informationRECOMMENDATION ITU-R P DIGITAL TOPOGRAPHIC DATABASES FOR PROPAGATION STUDIES. (Question ITU-R 202/3)
Rec. ITU-R P.1058-1 1 RECOMMENDATION ITU-R P.1058-1 DIGITAL TOPOGRAPHIC DATABASES FOR PROPAGATION STUDIES (Question ITU-R 202/3) Rec. ITU-R P.1058-1 (1994-1997) The ITU Radiocommunication Assembly, considering
More informationMetroPro Surface Texture Parameters
MetroPro Surface Texture Parameters Contents ROUGHNESS PARAMETERS...1 R a, R q, R y, R t, R p, R v, R tm, R z, H, R ku, R 3z, SR z, SR z X, SR z Y, ISO Flatness WAVINESS PARAMETERS...4 W a, W q, W y HYBRID
More informationSegmentation and Tracking of Partial Planar Templates
Segmentation and Tracking of Partial Planar Templates Abdelsalam Masoud William Hoff Colorado School of Mines Colorado School of Mines Golden, CO 800 Golden, CO 800 amasoud@mines.edu whoff@mines.edu Abstract
More informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationClustering Part 4 DBSCAN
Clustering Part 4 Dr. Sanjay Ranka Professor Computer and Information Science and Engineering University of Florida, Gainesville DBSCAN DBSCAN is a density based clustering algorithm Density = number of
More informationGravity 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 informationComputation of Slope
Computation of Slope Prepared by David R. Maidment and David Tarboton GIS in Water Resources Class University of Texas at Austin September 2011, Revised December 2011 There are various ways in which slope
More informationThe AIRGrav airborne gravity system
The AIRGrav airborne gravity system Stephan Sander, Malcolm Argyle, Stefan Elieff, Stephen Ferguson, Veronique Lavoie and Luise Sander Sander Geophysics Ltd. stephan@sgl.com Abstract An overview of the
More informationHigh-Precision Positioning Unit 2.2 Student Exercise: Calculating Topographic Change
High-Precision Positioning Unit 2.2 Student Exercise: Calculating Topographic Change Ian Lauer and Ben Crosby (Idaho State University) Change is an inevitable part of our natural world and varies as a
More information