Stabilization of the Euler deconvolution algorithm by means of a two steps regularization approach
|
|
- Peter Garrison
- 5 years ago
- Views:
Transcription
1 Stabilization of the Euler deconvolution algorithm by means of a two steps regularization approach R. Pašteka ( 1,2 ), D. Kušnirák ( 1 ), H.-J. Götze ( 2 ) ( 1 )Department of Applied Geophysics, Comenius University, Bratislava, Slovakia ( 2 )Institute of Geosciences, Dept. of Geophysics, Christian-Albrechts University, Kiel, Germany Summary Euler deconvolution is one of the best developed methods for semi-automated interpretation of potential fields in the last two decades. This paper describes a two steps procedure based on: the stabilization of the input derivatives by regularization (1.step) and involving the regularized solution of the linear equation system (included in a 2.step). This can significantly improve the Euler deconvolution (ED) method which leads to better interpretations of potential fields. In this study the two steps stabilization approach will be applied to a synthetic model and a real data set which stems from an archaeological microgravity survey. These two case studies were calculated by the classical ED and the new derivative ED methods (DED). The latter is a modification of the classical approach and bases on initial vertical derivative evaluation for input data and it results in a direct estimate of structural index values. Presented new two steps procedure provides better estimations of Euler depths and eliminates significantly false solutions. Introduction Since long, Euler deconvolution is one of the best developed methods and thoroughly tried and tested for any semi-automated interpretation of potential fields (e.g. Reid 1995, Reid 2007). The reason is that it bases on a clearly formulated theoretical background (the Euler s homogenity theorem), is relatively simple and ensures good accessibility to professional and/or freely available interpretation software. On the other hand this interpretation method is very sensitive to noise and errors in interpreted data sets. Anyway, it is very effective in the interpretation of well shaped and sampled anomalies, which are caused by isolated geological structures and near-surface objects (e.g.: Bilgili et al. 2009, Pašteka et al. 2009). The deconvolution of noisy data often results in a large amount of erroneous and defocused solutions. In the past many techniques have been published to improve and select satisfying solutions (or reject bad ones). The majority of them is focused on clustering algorithms and/or fulfilling some of selected statistical criteria (standard deviation and/or condition number of the solution). In some cases these methods are very helpful, however, we think that such improvements are included to late in the numerical process. Therefore, stabilization of solutions should start much earlier. Toward this end we are suggesting a two folded methodological procedure: a) stabilized evaluation of the input data derivatives (Pašteka et al. 2009) and b) stabilized solution of the linear equation system (by means of the Levenberg-Marquardt method). This second regularization step is based on an introduction of a smoothing term in the equation matrix system, controlled by a dumping (regularization) parameter. The search for the optimum value of this parameter remains still an open problem. At this stage we are experimenting with clustering principle - well behaved and focused depth-solution clusters are accepted. Here we will present the results and our experience with the two steps stabilization of the Euler method. In the first example we get use of synthetic gravity data along a profile with additional noise. The classical (standard) Euler Deconvolution will be applied with a prescribed structural index). In the second one real world grided microgravity data were interpreted. They stem from a survey, which was
2 conducted in a baroque church in Western Slovakia and was focused on the detection of medieval crypts. In this example we tested the derivative Euler Deconvolution which included the direct estimation of structural index values. Both results show that the new two steps procedure provides two kinds of results (a) more focused estimates of Euler depths and (b) significant reductions of false solutions. Methodology The "two steps stabilization" procedure is based on: a) a Fourier domain low-pass filtering of the evaluated higher derivatives, entering into the Euler deconvolution method by means of regularization (Pašteka et al. 2009) and b) use of the Levenberg Marquardt (LM-) method (e.g. Aster et al. 2005) which enables a stable solution of the linear equation system. This numerical method modifies the classical Gauss Newton method for the solution of equation systems by means of the term λi (where I is the identity matrix and λ the regularization parameter), entering into the main matrix of the equation system, which forms then a nonsingular matrix. The problem which is related to the Levenberg Marquardt method is determination of the optimum λ value. We do not have a realistic solution of this problem at the moment. Application of the L-curve (Hansen, 1994) did not lead satisfying results because the best λ values did not necessarily point to the important breaking-point of the L-curve. The optimum λ was estimated by means of a clustering principle which will be described below. In Fig 1. we present results which were obtained by the classical Euler Deconvolution method with a prescribed structural index (Thompson, 1982). As mentioned earlier we used synthetic gravity data. They are calculated along a profile over a 2D horizontal cylinder (SI = 1) in the depth of 20 m. White noise was added with the range of 5% of the main anomaly amplitude. The influence of the added noise is crucial and without any smoothing or regularizing of the orthogonal x- and z-derivatives, all obtained solutions are obviously wrong (Fig.1a, black crosses). By means of regularization of derivatives (Pašteka et al. 2009) the quality increased (Fig. 1a, blue crosses are building a cluster close to the location of the center of the cylinder). However, there still exists a larger amount of false solutions. Only the application of the second regularization step the LM-method excludes these false solutions (Fig. 1b), but we realize that the solution clustering strongly depends of the selection of the regularization parameter λ (best solutions: Fig. 1b, red crosses, λ=10-6 ). To select an optimum value for λ, we used a "clustering principle": Well developed clusters are selected which corresponds with the most important anomalous features of the field. Simultaneously a large amount of irregularly distributed false solution points have to been eliminated. For very low λ values (e.g ) solutions are almost identical with that for λ = 0 (not shown here). For higher values (e.g and larger) solutions are defocused and moving upwards, forming an artificial and false cluster at z = 0 m; it is well developed - shown in Fig. 1b, by violet crosses, λ = However, the influence of regularized derivatives is more important if integrated in the method of the derivative Euler Deconvolution (DED) (Hsu 2002, Fedi and Florio 2002). This modification is based on an initial evaluation of vertical derivative and subsequent application of the standard Euler deconvolution methodology. This initial evaluation of vertical derivative minimize and remove the role of the background term and the structural index value can be estimated directly from the equation system together with the coordinates of the source. In Fig.2 we present the results of the DED method, applied on a selected part of residual Bouguer anomaly field (correction density: 1.80 t m -3 ), obtained in the framework of a microgravity and GPR survey which aimed in the detection of medieval crypts in the St. Nicholas baroque church in Trnava, Western Slovakia (Pašteka et al. 2007). In this survey, seven unknown crypts have been detected by the above mentioned methods and finally verified by video-inspection with a mini-camera. The selected part of the residual Bouguer anomaly field contains two small negative anomalies (min.: approx. 45 μgal), caused by two crypts with almost identical depth position and volume. In this case
3 the solutions without regularized derivatives are defocused and too shallow (Fig. 2a and 2d; the approximate cross-section of the crypt based on GPR interpretation is plotted in blue lines). Fig.1 Results of classical Euler Deconvolution with a prescribed structural index value (SI = 1), obtained from the synthetic gravity of a 2D horizontal cylinder (depth 20 m) with added synthetic white noise; a) the results with and without regularized derivatives, b) the results with regularized derivatives and Levenberg-Marquardt method (for various values of damping parameter ). Fig.2 Results of the DED method, which includes a direct estimation of structural index values. They are obtained from data interpretation of a microgravity survey. In the upper row residual Bouguer anomaly maps are displayed together with DED solutions (colored dots); in the lower row
4 depth sections from a GPR survey are shown in combination with DED solutions; a) & d) present results without regularized derivatives, and = 0; b) & e) portray results which were calculated by regularized derivatives and = 0; finally c) & f) shows the results with regularized derivatives and = 10-5 (the optimum value, selected by the clustering principle) In the case of regularized derivatives no bigger differences can be detected in the map between the results which based on a processing with and without the LM-method (Fig. 2b and c). On contrary to Fig. 2b only a small artificial cluster with coordinates 85.8m, 86.3m disappeared from the map (Fig. 2c). However, in the xz-projections (Fig. 3e, f) it is shown that the amount of wrong solutions with SI = 1 is much lower. Actually this SI value is not real for such isometric anomalies, because it is more reliable for contact structures. In the case without LM-introduction the number of these solutions achieves a relative level of 65% (λ = 0) and with it 49% (λ = 10-5 ), if they are compared to the total number of all solutions. Clusters with SI=1 are similar and show the approx. depth of the crypt center. Conclusions Based on our numerical experiments and results we described the role of regularized derivatives (1. stabilization step) as an important item and the introduction of the Levenberg-Marquardt (2. stabilization step) as an item with minor importance. Without stabilization of the derivatives we could not obtain any reliable Euler solutions in a great variety of synthetic and real data studies. The LM-method can help to cancel defocused and erroneous solutions, but it strongly depends on the correct selection of the optimum value of the dumping (regularization) parameter λ. It is important to mention that in the case of noisy input data, final clustering or statistical selection of received solutions can not basically affect the received shallower depths (resulting from the noise influence). However, the presented "two step stabilization" procedure results in Euler solution points which are located closer to the real positions of the bodies. Again, the key role plays the search for the optimum regularization parameters. In the case of the regularized derivatives, the concept of C-norm functions analysis (Pašteka et al. 2009) led to acceptable results in many cases. The L-curve concept (Hansen, 1994) included in the LM-method did not result in satisfying solutions till now the most acceptable solutions (clustered around a real source position) was not found close to the important breaking-point of the L-curve. This task/problem remains still open, in the actual version of the presented procedure a rather time-consuming trial and error method is applied for searching the optimum λ-parameter based on the clustering principle. Development and test of the presented "two-steps stabilization procedure" for improvement of Euler Deconvolution have been done under the framework of bilateral cooperation between the geophysical departments of Comenius University in Bratislava and Christian-Albrechts University in Kiel. The Matlab programs REGDER and REGDED have been completed by RP in the framework of his Mercator Professor ship in Kiel which was financially supported by the Deutsche Forschungsgemeinschaft, DFG (INST 5659/1-1) project. References Aster,R., Borchers, B., Thurber, C., [2005] Parameter Estimation and Inverse Problems. Elsevier Press. Bilgili F., Götze H.-J., Pašteka R., Schmidt S., Hackney R. [2009] Intrusion versus Inversion - a 3D Density Model of the Southern Rim of the Northwest German Basin. International Journal of Earth Sciences, Vol. 98, Nr. 3, Fedi, M., Florio, G., [2002] Euler Deconvolution with no a priori definition of Structural Index. Geoph. Research Abstracts, vol. 4. Europ. Geophys. Soc. 27th General Assembly, Nice, april,
5 Hansen, P. C., [1994] Regularization tools. A MatLab Package for Analysis and Solution of Discrete Ill-posed Problems. Numerical Algorithms 6, Hsu S., [2002] Imaging magnetic sources using Euler s equation. Geophysical Prospecting, 50, Pašteka R., Terray M., Hajach M., Pašiaková, M., [2007] Microgravity measurements and GPR technique in the search for medieval crypts: a case study from the St. Nicholas church in Trnava, SW Slovakia. Proc. of the Archaeol. Prosp. 7th conference, Nitra, Štúdijné zvesti Vol. 41, Pašteka R., Richter F.P., Karcol R., Brazda K., Hajach M. [2009] Regularized derivatives of potential fields and their role in semi-automated interpretation methods. Geophys. Prospecting, Vol. 57, Nr. 4, Reid, A. B. [1995] Euler Deconvolution: Past, Present and Future. A Review. SEG, Exp. Abstr., Reid, A. [2007] Semi-Automated Methods of Potential Field Interpretation Innovations, and Recent and Future Developments. Expaned Abstracts, EGM2007 Intern. Workshop., 4 p.
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 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 informationESTIMATION OF SUBSURFACE QANATS DEPTH BY MULTI LAYER PERCEPTRON NEURAL NETWORK VIA MICROGRAVITY DATA
Advances in Geosciences Vol. 20: Solid Earth (2008) Ed. Kenji Satake c World Scientific Publishing Company ESTIMATION OF SUBSURFACE QANATS DEPTH BY MULTI LAYER PERCEPTRON NEURAL NETWORK VIA MICROGRAVITY
More informationEstimation 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 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 informationDiffuse Optical Tomography, Inverse Problems, and Optimization. Mary Katherine Huffman. Undergraduate Research Fall 2011 Spring 2012
Diffuse Optical Tomography, Inverse Problems, and Optimization Mary Katherine Huffman Undergraduate Research Fall 11 Spring 12 1. Introduction. This paper discusses research conducted in order to investigate
More informationbrahim KARA and Nihan HOSKAN
Acta Geophysica vol. 64, no. 6, Dec. 2016, pp. 2232-2243 DOI: 10.1515/acgeo-2016-0097 An Easy Method for Interpretation of Gravity Anomalies Due to Vertical Finite Lines brahim KARA and Nihan HOSKAN Department
More informationB15 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 informationQualitative 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 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 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 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 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 informationAdaptive Filters Algorithms (Part 2)
Adaptive Filters Algorithms (Part 2) Gerhard Schmidt Christian-Albrechts-Universität zu Kiel Faculty of Engineering Electrical Engineering and Information Technology Digital Signal Processing and System
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 informationTotal variation tomographic inversion via the Alternating Direction Method of Multipliers
Total variation tomographic inversion via the Alternating Direction Method of Multipliers Landon Safron and Mauricio D. Sacchi Department of Physics, University of Alberta Summary Geophysical inverse problems
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 informationMonte Carlo Statics: The Last Frontier
Monte Carlo Statics: The Last Frontier David LeMeur* CGG Canada, Calgary, Alberta, Canada dlemeur@cgg.com and Sophie Merrer CGG Canada, Calgary, Alberta, Canada Abstract Summary Most surface-consistent
More information1. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems:
ILL-POSED LEAST SQUARES PROBLEMS. Techniques for ill-posed least squares problems. We write throughout this lecture = 2. Consider the following least squares problems: where A = min b Ax, min b A ε x,
More informationThree-dimensional inversion of borehole, time-domain, electromagnetic data for highly conductive ore-bodies.
Three-dimensional inversion of borehole, time-domain, electromagnetic data for highly conductive ore-bodies. Nigel Phillips, Doug Oldenburg, Eldad Haber, Roman Shekhtman. UBC-Geophysical Inversion Facility
More informationDeconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Deconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY We use the recently introduced multiscale and multidirectional curvelet transform to exploit the
More informationPARAMETERIZATION AND SAMPLING DESIGN FOR WATER NETWORKS DEMAND CALIBRATION USING THE SINGULAR VALUE DECOMPOSITION: APPLICATION TO A REAL NETWORK
11 th International Conference on Hydroinformatics HIC 2014, New York City, USA PARAMETERIZATION AND SAMPLING DESIGN FOR WATER NETWORKS DEMAND CALIBRATION USING THE SINGULAR VALUE DECOMPOSITION: APPLICATION
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 informationCOMBINING MULTIPLE GEOPHYSICAL DATA SETS INTO A SINGLE 3D IMAGE. Abstract
COMBINING MULTIPLE GEOPHYSICAL DATA SETS INTO A SINGLE 3D IMAGE Jeffrey J. Daniels*, Mark Vendl+, Jennifer Holt*, and Erich Guy* *OSU, Columbus, OH +USEPA, Chicago, IL Abstract Traditional geophysical
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 informationThe art of well tying with new MATLAB tools
Art of well tying The art of well tying with new MATLAB tools Heather J.E. Lloyd and Gary F. Margrave ABSTRACT To tie seismic data to well logs, four steps are required. The first step involves editing
More informationTimelapse ERT inversion approaches and their applications
Timelapse ERT inversion approaches and their applications THOMAS GÜNTHER Leibniz-Institute for Applied Geophysics (LIAG), Hannover, Germany. thomas.guenther@liag-hannover.de Introduction Aim of timelapse
More informationAnalysis of different interpolation methods for uphole data using Surfer software
10 th Biennial International Conference & Exposition P 203 Analysis of different interpolation methods for uphole data using Surfer software P. Vohat *, V. Gupta,, T. K. Bordoloi, H. Naswa, G. Singh and
More informationModel parametrization strategies for Newton-based acoustic full waveform
Model parametrization strategies for Newton-based acoustic full waveform inversion Amsalu Y. Anagaw, University of Alberta, Edmonton, Canada, aanagaw@ualberta.ca Summary This paper studies the effects
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 informationINTERNATIONAL JOURNAL OF GEOMATICS AND GEOSCIENCES Volume 5, No 2, 2014
INTERNATIONAL JOURNAL OF GEOMATICS AND GEOSCIENCES Volume 5, No 2, 2014 Copyright by the authors - Licensee IPA- Under Creative Commons license 3.0 Research article ISSN 0976 4380 Edge detection process
More informationImproving the Discrimination Capability with an Adaptive Synthetic Discriminant Function Filter
Improving the Discrimination Capability with an Adaptive Synthetic Discriminant Function Filter 83 J. Ángel González-Fraga 1, Víctor H. Díaz-Ramírez 1, Vitaly Kober 1, and Josué Álvarez-Borrego 2 1 Department
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 informationUsing Statistical Techniques to Improve the QC Process of Swell Noise Filtering
Using Statistical Techniques to Improve the QC Process of Swell Noise Filtering A. Spanos* (Petroleum Geo-Services) & M. Bekara (PGS - Petroleum Geo- Services) SUMMARY The current approach for the quality
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 informationMo 21P1 08 Comparison of Different Acquisition Patterns for 2D Tomographic Resistivity Surveys
Mo 21P1 08 Comparison of Different Acquisition Patterns for 2D Tomographic Resistivity Surveys R. Martorana* (University of Palermo), P. Capizzi (University of Palermo), A. D'Alessandro (INGV - Roma) &
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 informationN. Kazemi* (University of Alberta), A. Gholami (University of Tehran) & M. D. Sacchi (University of Alberta)
3 May 2 June 216 Reed Messe Wien We STZ1 9 Modified Sparse Multichannel Blind Deconvolution N. Kazemi* (University of Alberta), A. Gholami (University of Tehran) & M. D. Sacchi (University of Alberta)
More informationChapter 5. 3D data examples REALISTICALLY COMPLEX SYNTHETIC INVERSION. Modeling generation and survey design
Chapter 5 3D data examples In this chapter I will demonstrate the e ectiveness of the methodologies developed in the previous chapters using 3D data examples. I will first show inversion results of a realistically
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 informationA MODELING STUDY OF LOW-FREQUENCY CSEM IN SHALLOW WATER
A MODELING STUDY OF LOW-FREQUENCY CSEM IN SHALLOW WATER Michael A. Frenkel and Sofia Davydycheva EMGS ASA ABSTRACT The applicability and resolution power of low-frequency Controlled-Source Electromagnetic
More informationHow reliable is statistical wavelet estimation? Jonathan Edgar*, BG Group, and Mirko van der Baan, University of Alberta
Jonathan Edgar*, BG Group, and Mirko van der Baan, University of Alberta Summary Well logs are often used for the estimation of seismic wavelets. The phase is obtained by forcing a well-derived synthetic
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 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 informationParameter Estimation in Differential Equations: A Numerical Study of Shooting Methods
Parameter Estimation in Differential Equations: A Numerical Study of Shooting Methods Franz Hamilton Faculty Advisor: Dr Timothy Sauer January 5, 2011 Abstract Differential equation modeling is central
More informationENHANCED RADAR IMAGING VIA SPARSITY REGULARIZED 2D LINEAR PREDICTION
ENHANCED RADAR IMAGING VIA SPARSITY REGULARIZED 2D LINEAR PREDICTION I.Erer 1, K. Sarikaya 1,2, H.Bozkurt 1 1 Department of Electronics and Telecommunications Engineering Electrics and Electronics Faculty,
More informationMostafa Naghizadeh and Mauricio D. Sacchi
Ground-roll elimination by scale and direction guided curvelet transform Mostafa Naghizadeh and Mauricio D. Sacchi Summary We propose a curvelet domain strategy to subtract ground-roll from seismic records.
More informationImproving Levenberg-Marquardt Algorithm Inversion Result Using Singular Value Decomposition
20 Earth Science Research; Vol. 5, No. 2; 2016 ISSN 1927-0542 E-ISSN 1927-0550 Published by Canadian Center of Science and Education Improving Levenberg-Marquardt Algorithm Inversion Result Using Singular
More informationEULER 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 information1D internal multiple prediction in a multidimensional world: errors and recommendations
1D internal multiple prediction 1D internal multiple prediction in a multidimensional world: errors and recommendations Pan Pan and Kris Innanen ABSTRACT Internal multiples are more difficult to estimate
More informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationData 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 informationHigh 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 informationJournal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp. 59 74 SOLID MECHANICS WAVE PROPAGATION DUE TO AN EMBEDDED SEISMIC SOURCE IN A GRADED HALF-PLANE WITH RELIEF PECULIARITIES.
More informationTu A4 09 3D CSEM Inversion Of Data Affected by Infrastructure
Tu A4 09 3D CSEM Inversion Of Data Affected by Infrastructure J.P. Morten (EMGS), L. Berre* (EMGS), S. de la Kethulle de Ryhove (EMGS), V. Markhus (EMGS) Summary We consider the effect of metal infrastructure
More informationUMASIS, an analysis and visualization tool for developing and optimizing ultrasonic inspection techniques
17th World Conference on Nondestructive Testing, 25-28 Oct 2008, Shanghai, China UMASIS, an analysis and visualization tool for developing and optimizing ultrasonic inspection techniques Abstract Joost
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 informationFINITE ELEMENT MODELING OF IP ANOMALOUS EFFECT FROM BODIES OF ANY GEOMETRICAL SHAPE LOCATED IN RUGED RELIEF AREA
FINITE ELEMENT MODELING OF IP ANOMALOUS EFFECT FROM BODIES OF ANY GEOMETRICAL SHAPE LOCATED IN RUGED RELIEF AREA Alfred FRASHËRI, Neki FRASHËRI Geophysical Technologist Training Program Quantec Geophysics
More informationSurface Wave Suppression with Joint S Transform and TT Transform
Available online at www.sciencedirect.com Procedia Earth and Planetary Science 3 ( 011 ) 46 5 011 Xian International Conference on Fine Geological Exploration and Groundwater & Gas Hazards Control in Coal
More informationA simple method for depth determination from self-potential anomalies due to two superimposed structures
CSIRO PUBLISHING Exploration Geophysics, 16, 47, 38 314 http://dx.doi.org/1.171/eg151 A simple method for depth determination from self-potential anomalies due to two superimposed structures El-Sayed M.
More informationIncreasing/Decreasing Behavior
Derivatives and the Shapes of Graphs In this section, we will specifically discuss the information that f (x) and f (x) give us about the graph of f(x); it turns out understanding the first and second
More informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011
GEOPHYS 242: Near Surface Geophysical Imaging Class 4: First-Arrival Traveltime Tomography Mon, April 11, 2011 Wavefront tracing methods - speed versus accuracy Inversion algorithms - any magic approach?
More informationIntegral 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 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 informationFermat s principle and ray tracing in anisotropic layered media
Ray tracing in anisotropic layers Fermat s principle and ray tracing in anisotropic layered media Joe Wong ABSTRACT I consider the path followed by a seismic signal travelling through velocity models consisting
More informationROBUST LINE-BASED CALIBRATION OF LENS DISTORTION FROM A SINGLE VIEW
ROBUST LINE-BASED CALIBRATION OF LENS DISTORTION FROM A SINGLE VIEW Thorsten Thormählen, Hellward Broszio, Ingolf Wassermann thormae@tnt.uni-hannover.de University of Hannover, Information Technology Laboratory,
More informationP063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA
1 P063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA I. TRINKS 1, S. SINGH 2, C. CHAPMAN 3, P. BARTON 1, M. BOSCH 4, A. CHERRETT 5 Addresses 1 Bullard Laboratories,
More informationTarget-oriented wave-equation inversion with regularization in the subsurface-offset domain
Stanford Exploration Project, Report 124, April 4, 2006, pages 1?? Target-oriented wave-equation inversion with regularization in the subsurface-offset domain Alejandro A. Valenciano ABSTRACT A complex
More informationPoints Lines Connected points X-Y Scatter. X-Y Matrix Star Plot Histogram Box Plot. Bar Group Bar Stacked H-Bar Grouped H-Bar Stacked
Plotting Menu: QCExpert Plotting Module graphs offers various tools for visualization of uni- and multivariate data. Settings and options in different types of graphs allow for modifications and customizations
More informationHigh definition tomography brings velocities to light Summary Introduction Figure 1:
Saverio Sioni, Patrice Guillaume*, Gilles Lambaré, Anthony Prescott, Xiaoming Zhang, Gregory Culianez, and Jean- Philippe Montel (CGGVeritas) Summary Velocity model building remains a crucial step in seismic
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 informationWave-equation migration from topography: Imaging Husky
Stanford Exploration Project, Report 123, October 31, 2005, pages 49 56 Short Note Wave-equation migration from topography: Imaging Husky Jeff Shragge 1 INTRODUCTION Imaging land seismic data is wrought
More informationIncreasing/Decreasing Behavior
Derivatives and the Shapes of Graphs In this section, we will specifically discuss the information that f (x) and f (x) give us about the graph of f(x); it turns out understanding the first and second
More informationTERRESTRIAL LASER SYSTEM TESTING USING REFERENCE BODIES
TERRESTRIAL LASER SYSTEM TESTING USING REFERENCE BODIES Miriam Zámečníková and Alojz Kopáčik Department of Surveying, Faculty of Civil Engineering Slovak University of Technology Email: miriam.zamecnikova@stuba.sk
More informationD020 Statics in Magnetotellurics - Shift or Model?
D020 Statics in Magnetotellurics - Shift or Model? W. Soyer* (WesternGeco-Geosystem), S. Hallinan (WesternGeco- Geosystem), R.L. Mackie (WesternGeco-Geosystem) & W. Cumming (Cumming Geoscience) SUMMARY
More informationSponge boundary condition for frequency-domain modeling
GEOPHYSIS, VOL. 60, NO. 6 (NOVEMBER-DEEMBER 1995); P. 1870-1874, 6 FIGS. Sponge boundary condition for frequency-domain modeling hangsoo Shin ABSTRAT Several techniques have been developed to get rid of
More informationA simulated Simultaneous Source Experiment in Shallow waters and the Impact of Randomization Schemes Rolf Baardman, Roald van Borselen, PGS
simulated Simultaneous Source Experiment in Shallow waters and the Impact of Randomization Schemes Rolf aardman, Roald van orselen, PGS Summary In simultaneous source acquisition, seismic data can be recorded
More informationAn imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults
Earth Planets Space, 52, 3 11, 2000 An imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults Takumi Murakoshi 1, Hiroshi Takenaka
More informationUMASIS, AN ANALYSIS AND VISUALIZATION TOOL FOR DEVELOPING AND OPTIMIZING ULTRASONIC INSPECTION TECHNIQUES
UMASIS, AN ANALYSIS AND VISUALIZATION TOOL FOR DEVELOPING AND OPTIMIZING ULTRASONIC INSPECTION TECHNIQUES A.W.F. Volker, J. G.P. Bloom TNO Science & Industry, Stieltjesweg 1, 2628CK Delft, The Netherlands
More informationON THE VELOCITY OF A WEIGHTED CYLINDER DOWN AN INCLINED PLANE
ON THE VELOCITY OF A WEIGHTED CYLINDER DOWN AN INCLINED PLANE Raghav Grover and Aneesh Agarwal RG (Grade 12 High School), AA (Grade 11 High School) Department of Physics, The Doon School, Dehradun. raghav.503.2019@doonschool.com,
More informationIterative resolution estimation in Kirchhoff imaging
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 369?? Iterative resolution estimation in Kirchhoff imaging Robert G. Clapp, Sergey Fomel, and Marie Prucha 1 ABSTRACT We apply iterative
More informationA Simulation Study and Experimental Verification of Hand-Eye-Calibration using Monocular X-Ray
A Simulation Study and Experimental Verification of Hand-Eye-Calibration using Monocular X-Ray Petra Dorn, Peter Fischer,, Holger Mönnich, Philip Mewes, Muhammad Asim Khalil, Abhinav Gulhar, Andreas Maier
More informationSeismic-to-well ties by smooth dynamic time warping
Seismic-to-well ties by smooth dynamic time warping Tianci Cui and Gary F. Margrave ABSTRACT Well tying by SDTW Without knowledge of Q or a check-shot/vsp survey, synthetic seismogram has to be manually
More informationA Multiple-Line Fitting Algorithm Without Initialization Yan Guo
A Multiple-Line Fitting Algorithm Without Initialization Yan Guo Abstract: The commonest way to fit multiple lines is to use methods incorporate the EM algorithm. However, the EM algorithm dose not guarantee
More informationINTERSECTION OF CURVES FACTS, COMPUTATIONS, APPLICATIONS IN BLOWUP*
South Bohemia Mathematical Letters Volume 24, (2016), No. 1, 10-16. INTERSECTION OF CURVES FACTS, COMPUTATIONS, APPLICATIONS IN BLOWUP* PAVEL CHALMOVIANSKÝ abstrakt. We deal with application of intersection
More informationParameter estimation. Christiano Gava Gabriele Bleser
Parameter estimation Christiano Gava Christiano.Gava@dfki.de Gabriele Bleser gabriele.bleser@dfki.de Introduction Previous lectures: P-matrix 2D projective transformations Estimation (direct linear transform)
More informationSelection of an optimised multiple attenuation scheme for a west coast of India data set
P-391 Selection of an optimised multiple attenuation scheme for a west coast of India data set Summary R Pathak*, PC Kalita, CPS Rana, Dr. S. Viswanathan, ONGC In recent years a number of new algorithms
More information1498. End-effector vibrations reduction in trajectory tracking for mobile manipulator
1498. End-effector vibrations reduction in trajectory tracking for mobile manipulator G. Pajak University of Zielona Gora, Faculty of Mechanical Engineering, Zielona Góra, Poland E-mail: g.pajak@iizp.uz.zgora.pl
More informationEffects 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 informationLAB 2: DATA FILTERING AND NOISE REDUCTION
NAME: LAB TIME: LAB 2: DATA FILTERING AND NOISE REDUCTION In this exercise, you will use Microsoft Excel to generate several synthetic data sets based on a simplified model of daily high temperatures in
More informationSEMI-BLIND IMAGE RESTORATION USING A LOCAL NEURAL APPROACH
SEMI-BLIND IMAGE RESTORATION USING A LOCAL NEURAL APPROACH Ignazio Gallo, Elisabetta Binaghi and Mario Raspanti Universitá degli Studi dell Insubria Varese, Italy email: ignazio.gallo@uninsubria.it ABSTRACT
More informationP. Bilsby (WesternGeco), D.F. Halliday* (Schlumberger Cambridge Research) & L.R. West (WesternGeco)
I040 Case Study - Residual Scattered Noise Attenuation for 3D Land Seismic Data P. Bilsby (WesternGeco), D.F. Halliday* (Schlumberger Cambridge Research) & L.R. West (WesternGeco) SUMMARY We show that
More informationConsiderations in 3D depth-specific P-S survey design
Considerations in 3D depth-specific P-S survey design Don C. Lawton and Peter W. Cary 3D P-S survey design ABSTRACT A new sparse-shot design for 3D P-S surveys is introduced. In the sparse shot design
More informationCalibration of a Multi-Camera Rig From Non-Overlapping Views
Calibration of a Multi-Camera Rig From Non-Overlapping Views Sandro Esquivel, Felix Woelk, and Reinhard Koch Christian-Albrechts-University, 48 Kiel, Germany Abstract. A simple, stable and generic approach
More informationForward and Inverse Modeling of Gravity Data: Locating buried glacial channels and evaluating the results of published analysis
GRAVITY COMPUTER LAB Forward and Inverse Modeling of Gravity Data: Locating buried glacial channels and evaluating the results of published analysis During this lab your task will be to evaluate the accuracy
More informationGEOPHYSICS. Time to depth conversion and seismic velocity estimation using time migration velocities
Time to depth conversion and seismic velocity estimation using time migration velocities Journal: Manuscript ID: Manuscript Type: Date Submitted by the Author: Geophysics draft Letters n/a Complete List
More informationComparison of 1-D inversion schemes for TDEM data
Comparison of -D inversion schemes for DEM data ILIAS PAPADOPOULOS, PANAGIOIS SOURLOS and FILLIPOS VALLIANAOS 2 Department of Geophysics, School of Geology, Aristotle University of hessaloniki 2 Department
More informationDetermination 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 informationA Geostatistical and Flow Simulation Study on a Real Training Image
A Geostatistical and Flow Simulation Study on a Real Training Image Weishan Ren (wren@ualberta.ca) Department of Civil & Environmental Engineering, University of Alberta Abstract A 12 cm by 18 cm slab
More information