Model-based Approach to Automatic 3D Seismic Horizon Correlation across Faults
|
|
- Marcia Moody
- 5 years ago
- Views:
Transcription
1 Model-based Approach to Automatic 3D Seismic Horizon Correlation across Faults Fitsum Admasu and Klaus Toennies Computer Vision Group Department of Simulation and Graphics Otto-von-Guericke Universität, Magdeburg, Germany Abstract Seismic data provide detailed information about subsurface structures. Reflection events visible in the seismic data are known as horizons, and indicate boundaries between different rock layers. A fault is a surface along which one side of rock layers has moved relative to the other in a direction parallel to the surface. Faults are recognized in seismic data by discontinuities of horizons. Fault interpretation is an important but time-consuming task of seismic interpretations. It is least supported by automatic tools due to seismic data distortions near fault regions. In this paper, we present an automatic method for correlations of horizons across faults in 3d seismic data. As automating horizons correlations using only seismic data features is not feasible, we find optimal correlations by imposing geological constraints. We formulate the correlation task as a non-rigid continuous point matching between the two sides of the fault. A fault surface is identified and seismic features on both sides of the fault are gathered. One side of the fault serves as the floating image while the other side is the reference image. A fault displacement model is constructed by performing initial discrete match of some prominent regions. Then the simulated annealing optimization technique is used to perform continuous point matching between the two sides according to seismic feature similarity and the fault model. A possible multiresolution scheme for the simulated annealing optimization is discussed. The method was applied to real 3D seismic data, and has shown to produce geologically acceptable results. Key Words: automatic seismic interpretation, model-based correspondence analysis. 1 Introduction Faults are breaks or planar surfaces in rocks across which there are displacements of rock layers. They exert control over oil formation by acting as seals or conduits of hydrocarbons. As result accurate assessment of fault geometry and displacements are crucial for the petroleum industry. Since 1920s seismic data have been used extensively for hypothesizing subsurface structures. Three-dimensional seismic data consist of numerous closely-spaced seismic lines that provide measures of subsurface reflectivity. Different subsurface s rock layers have 1
2 different acoustic impedances. As result when seismic waves are sent to underground structures, the changes in the seismic wave velocities give strong reflections being visible in the seismic images. These strong reflection events are known as horizons. Faults rarely give reflection events rather they are recognized in seismic data by the discontinuities of horizons events [5]. Structural interpretation of seismic data attempts to create a 3D subsurface model. Some of the tasks involved in structural interpretation are tracking of uninterrupted horizon segments, fault surfaces identification and correlating horizon segments across faults [2]. These interpretation tasks are done manually for some parts of the seismic data shown on figure 1. The horizon segments which are offset by the fault line are matched by the arrows. Horizon segment Fault line Timeline Crossline Inline Figure 1: 3-D seismic data with some interpretations. Our work aims at developing a computer-based methodology for correlation of horizons across faults. Seismic interpreters find displaced seismic features and connect them to each other on the basis of reflection characters and geological reasoning. By automating this task, we can save much time and avoid the inconsistencies existing among human interpreters. However automating the correlation task is very challenging. Seismic data contain a little image information, which is further obscured by noise. Since the two sides of the fault may also undergo different geological processes, such as compression and erosion, there are scale differences and some regions on the one side may not have matches on the other side. The methodology that we envisage is fusing the seismic data with information from a geological model in an iterative fashion. The method first finds and matches point regions that contain sufficiently distinctive structures on the two sides of the fault regions. Then the parameters required for the geological model are estimated based on displacements computed at those regions. Finally simulated annealing global search technique is used to find the optimal correlation between the two sides of the fault, optimal in sense of maximizing seismic features similarities while maintaining geologically valid solutions. The 2
3 method has been applied to fault patches extracted from real 3D seismic data and produces acceptable correspondence results. This paper is organized as follows. We review some of the existing automated structural interpretation tools in chapter 2. We formulate the horizons correlations task as correspondence analysis between two sides of a fault in chapter 3. Seismic similarity measures and the fault model with indications of how we perform the parameter estimations are described in chapter 4. The optimization problem and the possible multi-resolution simulated annealing scheme are explained in chapter 5. Our experimental results appear in chapter 6, followed by concluding remarks in chapter 7. 2 State of the Art Due to the complexity and size of seismic data sets, several automated tools are developed for structural interpretation of seismic data. Auto-picking or auto-trackers (reviewed in [2], [9]) have been commonly used to assist horizon tracking. Auto-picking tools are aimed at extending manually selected seismic traces based on local similarity measures. They perform well if there are uninterrupted horizon features. But horizon interruptions are very common. Steen et. al. [15] describe a semiautomatic method for fault tracking. They extract directive attributes in user- or data-steered directions, then use trained neural networks to enhance fault features from the directive attributes. Further the outputs of the neural networks are passed to the Canny non-maximum suppression edge detection algorithm [6] to generate the potential fault pixels. Bahorich et.al. [4] introduce a technique for fault highlighting through local coherence analysis. The coherence analysis measures the cross correlation of neighboring traces. Then the fault features are extracted as poorly correlated adjacent traces. Another method, reported by Gibson et.al [10], uses a coherency measure to detect points of significant horizon discontinuities. Then a highest confidence first merging strategy with a flexible surface model is applied to estimate the 3D fault surfaces iteratively. Alberts et.al. [1] explain a method for tracking horizons across discontinuities. They trained artificial neural networks (ANN) to track similar seismic intensity. However, horizon tracking across faults using solely seismic patterns is infeasible due to large seismic data distortion near faults. To alleviate this matter, Aurnhammer et.al [3] propose a model-based scheme for correlation of horizons at normal faults in 2D seismic images. They extract well-defined horizons segments on both sides of the fault and match the segments based on local correlation of seismic intensity and geological knowledge. Since exhaustive search for optimal solution of correlation is unfeasible, they suggest genetic algorithm as optimization technique. However, a pure two-dimensional approach lacks efficiency and is suitable only if the information of the 2D seismic slice is sufficient for evaluation of the geological constraints. We strive for exploiting existing three-dimensional spatial relationships in the data (such as continuity) directly for robust data analysis. 3
4 3 Problem Statement We formulate the horizons correlation task as a non-rigid continuous point matching between the two sides of the fault. Here continuous means each pixel on the left side of the fault is assigned a corresponding position on the right side. Our work concentrates on a single fault surface. Moreover, we assume that the fault surface has been already extracted as a plane with its associated patch. The seismic sediment features which are actually the grey (amplitude) values are averaged in five to ten pixels length along the horizons orientations starting at five pixels distance from the fault plane. The orientations of the horizons are defined by canny edge detector [6]. The average values are projected on the fault plane. These processes produce left and right fault feature 2D images (see Figure 2). The correspondence analysis problem is defined as finding displacement map for left fault image to overlap it with the right fault image so that corresponding positions in the two images are superposed. Fault Plane Left side of the fault mapped to the left plane Right side of the fault mapped to the right plane Horizons Fault patch Mapping from fault patch to fault plane Figure 2: A fault patch is mapped onto left and right fault planes. 4 Match Function We define a match function, ξ : Z 2 R, which measures the degree of match between the left and right side of the fault plane as they are overlaid on top of each other. The two images (left and right side of the fault plane) are given as two dimensional arrays and denoted by I l and I r where each of I l (x, y) and I r (x, y) maps to its respective averaged seismic features. The right image, I r, serves as reference image while I l (x, y) is the floating image which is displaced with displacement field T (x, y). T (x, y) is a 2D spatial coordinate transformation, such that T (x, y) =(x,y ). The match function, ξ, is then given as 4
5 ξ(t )=E s (T )+λ E g (T ) (1) where E s is the energy which is computed based on the seismic similarity, whereas E g is the energy that measures similarity of a given transformation and a geological model. They are described in details in section 4.1 and 4.2. λ is a negative scalar real value and used to balance E s and E g. 4.1 Computation of Seismic Similarity, E s Seismic similarity is defined as a mathematical measure of intensity similarity. In order to estimate the similarity, we choose the normalized cross-correlation technique. The main reasons for this choice are its potential accuracy and robustness with respect to noise as it is computed in some neighborhood of a point. For inputs I r, I l, and T, the normalized local cross-correlation function at point (x 0,y 0 ) with T (x 0,y 0 )=(x 0,y 0) is calculated as C(x 0,y 0 )= n x,y= n (I lxy) I rxy ) n x,y= n (I n lxy) x,y= n (I rxy) (2) where I lxy = I l (x 0 x, y 0 y) I l (x 0,y 0 ) I rxy = I r (x 0 x, y 0 y) I r (x 0,y 0 ) I l (x 0,y 0 ) and I r (x 0,y 0 ) are respectively the mean of I l and I r for n-neighborhood around point (x 0,y 0 ) and (x 0,y 0). The seismic similarity energy for given transformation is computed as the sum of the values of the normalized local cross-correlation at each point. 4.2 Computation of Geometrical Energy, E g The geometrical energy, E g, is computed as similarity between a given displacement map and a geological model. The geological model is based on pattern of displacement on a fault surface Fault Displacement Model Depending on the relative direction of displacement between the rocks on either side of the fault, faults are classified as normal, reverse or strike-slip [14](figure 3). The fault block above the fault surface is called the hanging wall, while the fault block below the fault is the footwall. We deal with normal faults. A normal fault is the most common fault type in which the hanging wall moves down relative to the the footwall. Thus normal fault displacements have only one direction, which means for any (x, y), T (x, y) =(x, y ). 5
6 Normal fault Footwall block Hanging wall block Thrust fault Strike-slip Fault Figure 3: Fault types. Furthermore, it is known from structural geology that horizons do not cross each other, that is for any T (x, a) =(x, a ) and T (x, b) =(x, b ), wehave a<b= a b (3) According to heuristics of Walsh et.al. [16] [17], a normalized displacement, D, at a point on a fault surface is given by D = 2(((1 + r)/2) 2 r 2 ) 2 (1 r) (4) where r is normalized radial distance from the fault center. The normalized displacement is D = d d max where d is the fault displacement at a point and d max is the maximum displacement on a fault surface Estimation of Fault Parameters The fault model of equation (4) assumes that the center of the fault, the width of the fault and the maximum displacement on the fault are known. But faults are often not contained to their complete extents in seismic data set. We have estimated the extent of the fault by computing landmark displacements. The simplest method for extraction of landmark corresponding points is to manually specify them. However, it is very difficult to specify accurate corresponding points and time-consuming. Thus, landmark displacements are extracted automatically in the following fashion. If the part of the fault ends in the seismic data, then there are regions with zero displacement. We take partially overlapping segments in those regions and propagate them to the 6
7 next slices. When any offset is found, the displacements are calculated and serve as landmarks. But if neither end of the fault is included in the seismic data, we identify particular prominent linear structures from the two sides of the fault features images. These prominent structures are extracted as segments by threshold of high contrasts. The segments can be matched with high confidence by maximizing the total cross-correlation values of the intensities around small neighborhood, as it is described in Aurnhammer [2]. The resulting landmark displacements are plugged into equation (4). Then the Marquardt-Levenberg constrained nonlinear optimization method, implemented in Matlab, is used to estimate the parameters for the fault displacement model. Finally, the theoretical transformation value for each point of the floating image is computed. The geometrical energy, E g, is computed as the least square error between any given transformation, T, and the theoretical transformation map, T computed. i.e E g (T )= i (T computed,i T i ) 2 (5) 5 Optimization After we have defined the transformation function and similarity measures, the next step is to find a suitable optimization procedure to generate the optimal transformation map, T max, which maximizes the value of the match function, ξ. T max = argmax T {T :Z2 Z 2 }(ξ(t)) (6) The general difficulty of the optimization of the match function in equation (6) is the nonlinearity existing in the search. It usually has many local maxima. Simple search strategies such as gradient decent are not appropriate. Therefore we use a simulated annealing (SA) [11], a stochastic non-linear optimization technique. 5.1 Simulated Annealing Optimization We set up the metropolis algorithm [13] to perform the minimization form of equation (6), based on the assumption that though the match function, ξ, might be noisy, its optimal would be surrounded by points whose values are also good. The metropolis algorithm incorporates ξ with the regularizer term which imposes a priori smoothness constraint on the solution (after [12]). It has to also make sure that all the current randomly generated candidate solutions satisfy the geological constraint defined at equation (3). We have used the geometric optimization schedule, [7]. The temperature is held fixed during each loop. At the end of each loop, k, the temperature,, is dropped according to the rule: k+1 = α k (7) The vales of the initial temperature ( 0 ), α and the number of iterations in each loop are to be determined experimentally. 7
8 In Simulation und Visualisierung 2004, Magdeburg, March Multi-resolution Optimization Though simulated annealing (SA) algorithms could find the global optimal results, it has high computational complexity [7]. We propose to take advantage of a multi-resolution analysis to increase the convergence rate of SA. The multi-resolution analysis is obtained by wavelet decomposition of the left and right side fault images. The decomposition is done by calculating the coefficient of a one dimensional continuous wavelet transform. Each column of the two dimensional image is fed to the one-dimensional continuous wavelet transform which computes the continuous wavelet coefficients of the input column vector at real, positive scales using a given type of wavelet. Figure 4 shows the result of the continuous 1-D wavelet coefficients using Daubechies wavelet [8]. The correspondence analysis between (a) and (b) on this figure is faster if we perform the analysis starting from the coarser-level and go to the finer-level. (a) (b) (c) (d) (e) (f) Figure 4: Wavelet decomposition: (a) and (b) show the left and the right-side fault feature images. (c) and (d) are respectively decompositions of (a) and (b) at coarse level. (e) and (f ) are respectively decompositions of (a) and (b) at coarser level. In this multi-resolution scenario, the simulated annealing optimization starts from lower 8
9 resolution level, then at higher resolution it searches with the convergence solutions reached at the lower resolution. Since the convergence solutions at lower resolution are close to the higher resolution optimum solutions, lower search temperature (that is, less number of iterations) can be used for searching at higher resolutions. But the multi-resolution analysis of the fault displacement constraint requires further investigations. This issues will be addressed in a paper to follow this. The results described in following chapter are done by single resolution optimization scheme. 6 Results and Discussion Our experimental data consist of seven fault patches taken from 3D seismic data which were surveyed from two different geographical locations. We have extracted the fault plane by linearly interpolating between manually provided seed points. The optimal solutions are defined as the would be solutions if experts performed the correlation. For the match function, we observed that the appropriate values of λ which compensates between the local seismic features and the global geological constraint range from 0.5 to 0.3. The initial temperature for the simulated annealing (SA) is set to values between 1000 and 1200 but still needs further experiments. We found that it is favorable for the optimization structure to cool down quickly and to stay more at lower temperatures. Thus the temperature is decremented at each step by 93% as cooling proceeds. The results of SA for a sample fault patch are shown on figure 5. To verify the results, we have restored the feature images to the 3D fault patch by inverting the feature mapping process (see figure 2). Some optimally correlated horizons are shown on two seismic slices on figure 6. These slices are taken from the restored 3D fault patch at locations I and II of figure 5 (a) and (b). The correlations are done according to the displacement map of figure 5 (d). The correlations have been verified by comparing them with manual interpretation. Our method fails in two test cases. The reason for failure appears to be insufficient information in seismic data to make correct initial discrete segments match. Interactions from nearby faults distort the fault displacement model and lead to incorrect global constraint (figure 7). However an overall analysis of the test cases reveals that the resulting correspondence analysis of our method is satisfactory with very good performance. 7 Summary and Conclusions We have presented an automatic method, which aligns locations of the left-side of the fault onto those of the right-side of the fault. The 3D spatial information of the seismic data is exploited in designing the fault displacement model and the smoothness constraints. Optimal displacement vector fields for correlating the displaced horizons are searched for maximizing the similarity of seismic image information in both sides while maintaining the fault model constraints. Since the search space is large and has many local maxima, we have utilized the stochastic stimulated annealing global search procedure. To speed up the convergence rate of the simulated annealing, the possible multi-resolution framework 9
10 In Simulation und Visualisierung 2004, Magdeburg, March I II (a) (c) I II (b) (d) Figure 5: (a) and (b) show respectively left and right fault images. (c) is the deformed image of (b) after aligning it with (a). (d) is displacement map of the alignment. is described. While the results are somewhat preliminary, they clearly demonstrate the applicability of our approach to real seismic data. Besides reducing the seismic interpreters time, the presented method, based on a quantitative model, provides repeatable data analysis tool. Moreover, the method can be used not only as an automatic verification tool for manually interpreted faults but also as a visualization tool for studying the full extent of the fault when the fault is not fully included in the seismic data set. However, the employed fault model needs to be modified to include fault interactions and lateral displacements. 8 Acknowledgements We would like to acknowledge Shell for the seismic data and stimulating discussions. We also thank Stefan Back and Janos Urai for their expertise advices regarding geology. The work presented here was supported by DFG Grant TO-166/
11 (a) (b) Figure 6: (a) and (b) seismic slices extracted at positions I and II of figure 5. The black arrows for some horizons indicate optimal correlation results of SA. The white curves show the fault lines. (a) (b) Figure 7: Seismic slices. The black arrows indicate correlation for horizons found by SA algorithm, while the white arrows show experts manual correlations for incorrect correlation of the black arrows. The white curves show fault lines. References [1] P. Alberts, M. Warner and D. Lister, Artificial Neural Networks for Simultaneous Multi Horizon Tracking across Discontinuities, 70th Annual International Meeting, Expanded Abstracts, Society of Exploration Geophysicists, Calgary, Canada, [2] M. Aurnhammer, Model-based Image Analysis for Automated Horizon Correlation across Faults in Seismic Data, PhD Thesis, University of Magdeburg, [3] M. Aurnhammer and K. Tönnies, Model-based analysis of geological structures in seismic images, In Simulation und Visualisierung, pp , Magdeburg, March, [4] M. Bahorich and S. Farmer, 3-D seismic Discontinuity for faults and stratigraphic features, The leading edge, Vol. 14, No. 10, pp , October,
12 [5] A. Brown, Interpretation of Three-Dimensional Seismic Data, American Association of Petroleum Geologists, 5th edition, December, [6] J. Canny, A computational approach to edge detection, IEEE Trans. patt. anal. mach. intell., Vol. 8, No. 6, pp , [7] H. Cohn and M. Fielding, Simulated annealing: searching for an optimal temperature schedule, Society for Industrial and Applied Mathematics, Journal of Optimization, Vol. 9, No. 3, pp , [8] I. Daubechies, Ten lectures on wavelets, Society for Industrial and Applied Mathematics, Philadelphia, [9] G. Dorn, Modern 3-D Seismic Interpretation, The Leading Edge, Vol. 17, No. 9, pp , [10] D. Gibson, M. Spann and J. Turner, Automatic Fault Detection for 3D Seismic Data, Proc. 7th Digital Image Computing: Techniques and Applications, pp , Sydney, December, [11] S. Kirkpatrick, J. Gelatt and M. Vecchi, Optimization by Simulated Annealing, Science, Vol. 220, No. 4598, pp , [12] S. Li, Markov Random Field Modeling in Computer Vision, Springer-Verlag, [13] N. Metropolis, A. Rosenbluth, M. Rosenbluth, A. Teller and E. Teller, Equations of state calculations by fast computing machines, J. Chemical Physics, Vol. 21, No. 6, pp , [14] E. Moores and R. Twiss, Structural Geology, W. H. Freeman and Company, [15] K. Tingdahl, Ø. Steen, P. Meldahl and J. Ligtenberg, Semi-automatic detection of faults in 3-D seismic signals, Society of Exploration Geophysicists, 71st Annual Meeting, Expanded Abstracts, pp , San Antonio, [16] J. Walsh and J. Watterson, Distributions of cumulative displacement and seismic slip on a single normal fault surface, Journal of Structural Geology, Vol. 9, No. 8, pp , [17] J. Walsh and J. Watterson, Analysis of the relationship between displacements and dimensions of faults, Journal of Structural Geology, Vol. 10, No. 3, pp ,
Model-based analysis of geological structures in seismic images
Model-based analysis of geological structures in seismic images Melanie Aurnhammer and Klaus Tönnies Computer Vision Group, Department of Simulation and Graphics Otto-von-Guericke University Magdeburg,
More informationMultiple horizons mapping: A better approach for maximizing the value of seismic data
Multiple horizons mapping: A better approach for maximizing the value of seismic data Das Ujjal Kumar *, SG(S) ONGC Ltd., New Delhi, Deputed in Ministry of Petroleum and Natural Gas, Govt. of India Email:
More informationThe coherence cube. MIKE BAHORICH Amoco Corporation Denver, CO. Faults parallel to strike. Conventional amplitude time
3-D seismic discontinuity for faults and stratigraphic features: The coherence cube MIKE BAHORICH Amoco Corporation Denver, CO STEVE FARMER Amoco Corporation Tulsa, OK Seismic data are usually acquired
More informationtechnical article Satinder Chopra 1*, Kurt J. Marfurt 2 and Ha T. Mai 2
first break volume 27, October 2009 technical article Using automatically generated 3D rose diagrams for correlation of seismic fracture lineaments with similar lineaments from attributes and well log
More informationConstrained Fault Construction
Constrained Fault Construction Providing realistic interpretations of faults is critical in hydrocarbon and mineral exploration. Faults can act as conduits or barriers to subsurface fluid migration and
More informationCULTURAL EDITING OF HRAM DATA COMPARISON OF TECHNIQUES. Canadian Journal of Exploration Geophysics, no. 1&2, vol. 34, 1998, pp.
CULTURAL EDITING OF HRAM DATA COMPARISON OF TECHNIQUES H. H. Hassan 1, J. W. Peirce 1, W. C. Pearson 2 and M. J. Pearson 3 Canadian Journal of Exploration Geophysics, no. 1&2, vol. 34, 1998, pp. 16-22
More informationAn Integrated approach for faults and fractures delineation with dip and curvature attributes
10 th Biennial International Conference & Exposition P 265 An Integrated approach for faults and fractures delineation with dip and curvature attributes Santosh, D.*, Aditi, B., Poonam, K., Priyanka S.,
More informationMultimedia Databases. Previous Lecture. 4.1 Multiresolution Analysis. 4 Shape-based Features. 4.1 Multiresolution Analysis
Previous Lecture Multimedia Databases Texture-Based Image Retrieval Low Level Features Tamura Measure, Random Field Model High-Level Features Fourier-Transform, Wavelets Wolf-Tilo Balke Silviu Homoceanu
More informationMultimedia Databases. Wolf-Tilo Balke Philipp Wille Institut für Informationssysteme Technische Universität Braunschweig
Multimedia Databases Wolf-Tilo Balke Philipp Wille Institut für Informationssysteme Technische Universität Braunschweig http://www.ifis.cs.tu-bs.de 4 Previous Lecture Texture-Based Image Retrieval Low
More informationThe SPE Foundation through member donations and a contribution from Offshore Europe
Primary funding is provided by The SPE Foundation through member donations and a contribution from Offshore Europe The Society is grateful to those companies that allow their professionals to serve as
More informationIntroduction Faults blind attitude strike dip
Chapter 5 Faults by G.H. Girty, Department of Geological Sciences, San Diego State University Page 1 Introduction Faults are surfaces across which Earth material has lost cohesion and across which there
More informationMultimedia Databases. 4 Shape-based Features. 4.1 Multiresolution Analysis. 4.1 Multiresolution Analysis. 4.1 Multiresolution Analysis
4 Shape-based Features Multimedia Databases Wolf-Tilo Balke Silviu Homoceanu Institut für Informationssysteme Technische Universität Braunschweig http://www.ifis.cs.tu-bs.de 4 Multiresolution Analysis
More informationCourse Title: Discipline: Geology Level: Basic-Intermediate Duration: 5 Days Instructor: Prof. Charles Kluth. About the course: Audience: Agenda:
Course Title: Structural Geology Discipline: Geology Level: Basic-Intermediate Duration: 5 Days Instructor: Prof. Charles Kluth About the course: This course covers the basic ideas of structural geometry
More information5. Simulated Annealing 5.1 Basic Concepts. Fall 2010 Instructor: Dr. Masoud Yaghini
5. Simulated Annealing 5.1 Basic Concepts Fall 2010 Instructor: Dr. Masoud Yaghini Outline Introduction Real Annealing and Simulated Annealing Metropolis Algorithm Template of SA A Simple Example References
More informationWe A1 12 Method to Generate a Watertight Geological Model Directly From a Seismic Volume
We A1 12 Method to Generate a Watertight Geological Model Directly From a Seismic Volume F. Pauget (ELIIS), S. Lacaze* (ELIIS) Summary This paper presents a novel method to compute a watertight geological
More informationMerging chronostratigraphic modeling and global interpretation Emmanuel Labrunye, Stanislas Jayr, Paradigm
Emmanuel Labrunye, Stanislas Jayr, Paradigm SUMMRY This paper proposes to combine the automatic interpretation of horizons in a seismic cube and the power of a space/time framework to quickly build an
More information2011 SEG SEG San Antonio 2011 Annual Meeting 1134
Seismic Stratigraphic Interpretation from a Geological Model A North Sea Case Study Sébastien Lacaze* and Fabien Pauget, Eliis, Michel Lopez and Aurélien Gay, University of Montpellier II, Marion Mangue,
More information23855 Rock Physics Constraints on Seismic Inversion
23855 Rock Physics Constraints on Seismic Inversion M. Sams* (Ikon Science Ltd) & D. Saussus (Ikon Science) SUMMARY Seismic data are bandlimited, offset limited and noisy. Consequently interpretation of
More informationDownloaded 10/02/18 to Redistribution subject to SEG license or copyright; see Terms of Use at
Multi-scenario, multi-realization seismic inversion for probabilistic seismic reservoir characterization Kester Waters* and Michael Kemper, Ikon Science Ltd. Summary We propose a two tiered inversion strategy
More informationKinematic structural forward modeling for fault trajectory prediction in seismic interpretation
Fault prediction by forward modeling Kinematic structural forward modeling for fault trajectory prediction in seismic interpretation Mohammed Alarfaj and Don C. Lawton ABSTRACT The unique relationship
More informationUsing Curvature to Map Faults, Fractures
Using Curvature to Map Faults, Fractures by SATINDER CHOPRA and KURT J. MARFURT Editor s note: Chopra is with Arcis Corp., Calgary, Canada; Marfurt is with the University of Oklahoma. Both are AAPG members.
More informationSeismic attributes for fault/fracture characterization
Satinder Chopra* and Kurt J. Marfurt + *Arcis Corporation, Calgary; + University of Houston, Houston Summary Seismic attributes have proliferated in the last three decades at a rapid rate and have helped
More information3D curvature attributes: a new approach for seismic interpretation
first break volume 26, April 2008 special topic 3D curvature attributes: a new approach for seismic interpretation Pascal Klein, Loic Richard, and Huw James* (Paradigm) present a new method to compute
More informationSeismic tomography with co-located soft data
Seismic tomography with co-located soft data Mohammad Maysami and Robert G. Clapp ABSTRACT There is a wide range of uncertainties present in seismic data. Limited subsurface illumination is also common,
More information3-D seismic continuity attribute for mapping discontinuities in a target zone: optimum parameters and evaluation
JOURNAL OF THE BALKAN GEOPHYSICAL SOCIETY, Vol. 2, No 4, November 1999, p. 112-119, 7 figs. 3-D seismic continuity attribute for mapping discontinuities in a target zone: optimum parameters and evaluation
More informationFAULT SLICING IN THE INTERPRETATION OF FAULTS IN SEISMIC DATA PROC- ESSING IN ATALA PROSPECT OF RIVER STATE, NIGERIA
FAULT SLICING IN THE INTERPRETATION OF FAULTS IN SEISMIC DATA PROC- ESSING IN ATALA PROSPECT OF RIVER STATE, NIGERIA Egbai, J.C. Department of Physics, Delta State University, Abraka ABSTRACT e-mail: jamesegbai@yahoo.com
More informationWe Improved Salt Body Delineation Using a new Structure Extraction Workflow
We-08-08 Improved Salt Body Delineation Using a new Structure Extraction Workflow A. Laake* (WesternGeco) SUMMARY Current salt imaging workflows require thorough geological understanding in the selection
More informationTemporal analysis for implicit compensation of local variations of emission coefficient applied for laser induced crack checking
More Info at Open Access Database www.ndt.net/?id=17661 Abstract Temporal analysis for implicit compensation of local variations of emission coefficient applied for laser induced crack checking by G. Traxler*,
More informationImproving the quality of Velocity Models and Seismic Images. Alice Chanvin-Laaouissi
Improving the quality of Velocity Models and Seismic Images Alice Chanvin-Laaouissi 2015, 2015, PARADIGM. PARADIGM. ALL RIGHTS ALL RIGHTS RESERVED. RESERVED. Velocity Volumes Challenges 1. Define sealed
More informationEstimators for Orientation and Anisotropy in Digitized Images
Estimators for Orientation and Anisotropy in Digitized Images Lucas J. van Vliet and Piet W. Verbeek Pattern Recognition Group of the Faculty of Applied Physics Delft University of Technolo Lorentzweg,
More informationQuantitative Seismic Interpretation An Earth Modeling Perspective
Quantitative Seismic Interpretation An Earth Modeling Perspective Damien Thenin*, RPS, Calgary, AB, Canada TheninD@rpsgroup.com Ron Larson, RPS, Calgary, AB, Canada LarsonR@rpsgroup.com Summary Earth models
More informationDelineating Karst features using Advanced Interpretation
P-152 Asheesh Singh, Sibam Chakraborty*, Shafique Ahmad Summary We use Amplitude, Instantaneous Phase, Trace Envelope and Dip of Maximum Similarity Attributes as a tool to delineate Karst induced features
More information11301 Reservoir Analogues Characterization by Means of GPR
11301 Reservoir Analogues Characterization by Means of GPR E. Forte* (University of Trieste) & M. Pipan (University of Trieste) SUMMARY The study of hydrocarbon reservoir analogues is increasing important
More informationKurt Marfurt Arnaud Huck THE ADVANCED SEISMIC ATTRIBUTES ANALYSIS
The Society of Exploration Geophysicists and GeoNeurale announce Kurt Marfurt Arnaud Huck THE ADVANCED SEISMIC ATTRIBUTES ANALYSIS 3D Seismic Attributes for Prospect Identification and Reservoir Characterization
More informationKMS Technologies KJT Enterprises Inc. Publication
KMS Technologies KJT Enterprises Inc. Publication Yu, G., Ma, Y., Zhu, X., Guo, T., Rebec, T. & Azbel, K. 2006 Reservoir characterization with high frequency bandwidth seismic data and coherence processing
More informationFault seal analysis: a regional calibration Nile delta, Egypt
International Research Journal of Geology and Mining (IRJGM) (2276-6618) Vol. 3(5) pp. 190-194, June, 2013 Available online http://www.interesjournals.org/irjgm Copyright 2013 International Research Journals
More informationEnhancing the resolution of CSEM inversion using seismic constraints Peter Harris*, Rock Solid Images AS Lucy MacGregor, OHM Surveys Ltd
Enhancing the resolution of CSEM inversion using seismic constraints Peter Harris*, Rock Solid Images AS Lucy MacGregor, OHM Surveys Ltd Summary The vertical resolution of inverted CSEM data remains an
More informationSeismic attributes of time-vs. depth-migrated data using self-adaptive window
Seismic attributes of time-vs. depth-migrated data using self-adaptive window Tengfei Lin*, Bo Zhang, The University of Oklahoma; Zhifa Zhan, Zhonghong Wan, BGP., CNPC; Fangyu Li, Huailai Zhou and Kurt
More informationMapping Magnetic Lineaments in the Foothills of Northeastern British Columbia using 2-D Wavelet Transform
Mapping Magnetic Lineaments in the Foothills of Northeastern British Columbia using 2-D Wavelet Transform Hassan Hassan* GEDCO, Calgary, Alberta, Canada hassan@gedco.com Abstract Summary This work describes
More informationUnravel Faults on Seismic Migration Images Using Structure-Oriented, Fault-Preserving and Nonlinear Anisotropic Diffusion Filtering
PROCEEDINGS, 44th Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2019 SGP-TR-214 Unravel Faults on Seismic Migration Images Using Structure-Oriented,
More informationWhat do you think this is? Conceptual uncertainty in geoscience interpretation.
Bond et al., GSA Today v. 17, no. 11, doi: 10.1130/GSAT01711A.1 What do you think this is? Conceptual uncertainty in geoscience interpretation. 1km How would you interpret this section? Midland Valley
More informationPETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR
PETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR APPLIED GRADUATE STUDIES Geology Geophysics GEO1 Introduction to the petroleum geosciences GEO2 Seismic methods GEO3 Multi-scale geological analysis GEO4
More informationNeural Networks biological neuron artificial neuron 1
Neural Networks biological neuron artificial neuron 1 A two-layer neural network Output layer (activation represents classification) Weighted connections Hidden layer ( internal representation ) Input
More informationLecture 6: Edge Detection. CAP 5415: Computer Vision Fall 2008
Lecture 6: Edge Detection CAP 5415: Computer Vision Fall 2008 Announcements PS 2 is available Please read it by Thursday During Thursday lecture, I will be going over it in some detail Monday - Computer
More informationFinite difference elastic modeling of the topography and the weathering layer
Finite difference elastic modeling of the topography and the weathering layer Saul E. Guevara and Gary F. Margrave ABSTRACT Finite difference 2D elastic modeling is used to study characteristics of the
More informationMotivation, Basic Concepts, Basic Methods, Travelling Salesperson Problem (TSP), Algorithms
Motivation, Basic Concepts, Basic Methods, Travelling Salesperson Problem (TSP), Algorithms 1 What is Combinatorial Optimization? Combinatorial Optimization deals with problems where we have to search
More informationNEW GEOLOGIC GRIDS FOR ROBUST GEOSTATISTICAL MODELING OF HYDROCARBON RESERVOIRS
FOR ROBUST GEOSTATISTICAL MODELING OF HYDROCARBON RESERVOIRS EMMANUEL GRINGARTEN, BURC ARPAT, STANISLAS JAYR and JEAN- LAURENT MALLET Paradigm Houston, USA. ABSTRACT Geostatistical modeling of reservoir
More informationPitfalls of seismic interpretation in prestack time- vs. depthmigration
2104181 Pitfalls of seismic interpretation in prestack time- vs. depthmigration data Tengfei Lin 1, Hang Deng 1, Zhifa Zhan 2, Zhonghong Wan 2, Kurt Marfurt 1 1. School of Geology and Geophysics, University
More informationLecture 8: Interest Point Detection. Saad J Bedros
#1 Lecture 8: Interest Point Detection Saad J Bedros sbedros@umn.edu Last Lecture : Edge Detection Preprocessing of image is desired to eliminate or at least minimize noise effects There is always tradeoff
More informationThin-Bed Reflectivity An Aid to Seismic Interpretation
Thin-Bed Reflectivity An Aid to Seismic Interpretation Satinder Chopra* Arcis Corporation, Calgary, AB schopra@arcis.com John Castagna University of Houston, Houston, TX, United States and Yong Xu Arcis
More informationLecture 8: Interest Point Detection. Saad J Bedros
#1 Lecture 8: Interest Point Detection Saad J Bedros sbedros@umn.edu Review of Edge Detectors #2 Today s Lecture Interest Points Detection What do we mean with Interest Point Detection in an Image Goal:
More informationFractured Volcanic Reservoir Characterization: A Case Study in the Deep Songliao Basin*
Fractured Volcanic Reservoir Characterization: A Case Study in the Deep Songliao Basin* Desheng Sun 1, Ling Yun 1, Gao Jun 1, Xiaoyu Xi 1, and Jixiang Lin 1 Search and Discovery Article #10584 (2014) Posted
More informationImproving Resolution with Spectral Balancing- A Case study
P-299 Improving Resolution with Spectral Balancing- A Case study M Fatima, Lavendra Kumar, RK Bhattacharjee, PH Rao, DP Sinha Western Offshore Basin, ONGC, Panvel, Mumbai Summary: The resolution limit
More informationMachine Learning in Modern Well Testing
Machine Learning in Modern Well Testing Yang Liu and Yinfeng Qin December 11, 29 1 Introduction Well testing is a crucial stage in the decision of setting up new wells on oil field. Decision makers rely
More informationConfidence and curvature estimation of curvilinear structures in 3-D
Confidence and curvature estimation of curvilinear structures in 3-D P. Bakker, L.J. van Vliet, P.W. Verbeek Pattern Recognition Group, Department of Applied Physics, Delft University of Technology, Lorentzweg
More information2011 SEG SEG San Antonio 2011 Annual Meeting 771. Summary. Method
Geological Parameters Effecting Controlled-Source Electromagnetic Feasibility: A North Sea Sand Reservoir Example Michelle Ellis and Robert Keirstead, RSI Summary Seismic and electromagnetic data measure
More informationSynthetic seismic modelling and imaging of an impact structure
Modelling and imaging of an impact structure Synthetic seismic modelling and imaging of an impact structure Matteo Niccoli ABSTRACT A geologic depth and velocity model of an impact crater was created and
More informationA Petroleum Geologist's Guide to Seismic Reflection
A Petroleum Geologist's Guide to Seismic Reflection William Ashcroft WILEY-BLACKWELL A John Wiley & Sons, Ltd., Publication Contents Preface Acknowledgements xi xiii Part I Basic topics and 2D interpretation
More informationA knowledge-based approach of seismic interpretation : horizon and dip-fault detection by means of cognitive vision.
A knowledge-based approach of seismic interpretation : horizon and dip-fault detection by means of cognitive vision. Philippe Verney, Jean-Francois Rainaud, Michel Perrin, Monique Thonnat To cite this
More informationRegularizing seismic inverse problems by model reparameterization using plane-wave construction
GEOPHYSICS, VOL. 71, NO. 5 SEPTEMBER-OCTOBER 2006 ; P. A43 A47, 6 FIGS. 10.1190/1.2335609 Regularizing seismic inverse problems by model reparameterization using plane-wave construction Sergey Fomel 1
More informationLaplacian Filters. Sobel Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters
Sobel Filters Note that smoothing the image before applying a Sobel filter typically gives better results. Even thresholding the Sobel filtered image cannot usually create precise, i.e., -pixel wide, edges.
More informationConjugate Directions for Stochastic Gradient Descent
Conjugate Directions for Stochastic Gradient Descent Nicol N Schraudolph Thore Graepel Institute of Computational Science ETH Zürich, Switzerland {schraudo,graepel}@infethzch Abstract The method of conjugate
More informationSummary. We present the results of the near-surface characterization for a 3D survey in thrust belt area in Sharjah, United Arab Emirates.
Near-surface characterization, challenges, and solutions for high-density, high-productivity, Alexander Zarkhidze*, Claudio Strobbia, Abdallah Ibrahim, WesternGeco; Luis Viertel Herrera, Abdulla Al Qadi,
More informationIdentifying faults and gas chimneys using multiattributes and neural networks
CORNER INTERPRETER S Coordinated by Linda R. Sternbach Identifying faults and gas chimneys using multiattributes and neural networks PAUL MELDAHL and ROAR HEGGLAND, Statoil, Stavanger, Norway BERT BRIL
More informationRock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E.
Rock Physics and Quantitative Wavelet Estimation for Seismic Interpretation: Tertiary North Sea R.W.Simm 1, S.Xu 2 and R.E.White 2 1. Enterprise Oil plc, Grand Buildings, Trafalgar Square, London WC2N
More informationEnhanced Subsurface Interpolation by Geological Cross-Sections by SangGi Hwang, PaiChai University, Korea
Enhanced Subsurface Interpolation by Geological Cross-Sections by SangGi Hwang, PaiChai University, Korea Abstract Subsurface geological structures, such as bedding, fault planes and ore body, are disturbed
More informationStratimagic. Seismic Facies Classification
Stratimagic Seismic Facies Classification 1 Stratimagic A truly geological interpretation of seismic data Stratimagic seismic facies maps allow geoscientists to rapidly gain insight into the depositional
More informationRecent developments in object modelling opens new era for characterization of fluvial reservoirs
Recent developments in object modelling opens new era for characterization of fluvial reservoirs Markus L. Vevle 1*, Arne Skorstad 1 and Julie Vonnet 1 present and discuss different techniques applied
More informationSeismic Inversion on 3D Data of Bassein Field, India
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-2004, India PP 526-532 Seismic Inversion on 3D Data of Bassein Field, India K.Sridhar, A.A.K.Sundaram, V.B.G.Tilak & Shyam Mohan Institute
More informationFUNDAMENTALS OF SEISMIC EXPLORATION FOR HYDROCARBON
FUNDAMENTALS OF SEISMIC EXPLORATION FOR HYDROCARBON Instructor : Kumar Ramachandran 10 14 July 2017 Jakarta The course is aimed at teaching the physical concepts involved in the application of seismic
More informationColombia s Offshore*
PS A Seismic-Structural Interpretation, on the Identification of Possible Causes in the Formation of Gas Chimneys in Colombia s Offshore* Tatiana Mayorga 1, Andrés E. Calle 2, Freddy M. Niño 2, Jorge Rubiano
More informationMETA ATTRIBUTES: A NEW CONCEPT DETECTING GEOLOGIC FEATURES & PREDICTING RESERVOIR PROPERTIES
META ATTRIBUTES: A NEW CONCEPT DETECTING GEOLOGIC FEATURES & PREDICTING RESERVOIR PROPERTIES FredAminzadeh*,FrisoBrouwer*,DavidConnolly* and Sigfrido Nielsen** *dgb-usa,onesugarcreekcenterblvd.,suite935,sugarland,tx,77478usa,tel.281-2403939,info.usa@dgb-group.com
More informationChurning seismic attributes with principal component analysis
Satinder Chopra + * and Kurt J. Marfurt + Arcis Seismic Solutions, Calgary; The University of Oklahoma, Norman Summary Seismic attributes are an invaluable aid in the interpretation of seismic data. Different
More informationLab 7: STRUCTURAL GEOLOGY FOLDS AND FAULTS
Lab 7: STRUCTURAL GEOLOGY FOLDS AND FAULTS This set of labs will focus on the structures that result from deformation in earth s crust, namely folds and faults. By the end of these labs you should be able
More informationSemi-Automated Stratigraphic Interpretation The HorizonCube
Fig. 2: A horizon slice for (a) non-steered similarity and (b) dipsteered similarity. Note that the dip-steered similarity is essentially showing correct definitions of the fault patterns. of 1 Km (40
More informationSalt Dome Detection and Tracking Using Texture Analysis and Tensor-based Subspace Learning
Salt Dome Detection and Tracking Using Texture Analysis and Tensor-based Subspace Learning Zhen Wang*, Dr. Tamir Hegazy*, Dr. Zhiling Long, and Prof. Ghassan AlRegib 02/18/2015 1 /42 Outline Introduction
More informationPart I. PRELAB SECTION To be completed before labs starts:
Student Name: Physical Geology 101 Laboratory #13 Structural Geology II Drawing and Analyzing Folds and Faults Grade: Introduction & Purpose: Structural geology is the study of how geologic rock units
More informationFOLDS AND THRUST SYSTEMS IN MASS TRANSPORT DEPOSITS
FOLDS AND THRUST SYSTEMS IN MASS TRANSPORT DEPOSITS G.I Aslop, S. Marco, T. Levi, R. Weinberger Presentation by Aaron Leonard INTRODUCTION Examine fold and thrust geometries associated with downslope movement
More informationEstimating Fault Seal and Capillary Sealing Properties in the Visund Field, North Sea A study carried out for Norsk Hydro
Estimating Fault Seal and Capillary Sealing Properties in the Visund Field, North Sea A study carried out for Norsk Hydro Abstract This study investigates the difference in seal/leakage mechanisms across
More informationObservation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story
Observation of shear-wave splitting from microseismicity induced by hydraulic fracturing: A non-vti story Petr Kolinsky 1, Leo Eisner 1, Vladimir Grechka 2, Dana Jurick 3, Peter Duncan 1 Summary Shear
More informationMultiple realizations: Model variance and data uncertainty
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Multiple realizations: Model variance and data uncertainty Robert G. Clapp 1 ABSTRACT Geophysicists typically produce a single model,
More informationTHE HETEROGENEOUS STRUCTURE OF FAULT ZONES WITHIN CARBONATE ROCKS: EVIDENCE FROM OUTCROP STUDIES AND IMPLICATIONS FOR FLUID FLOW
THE HETEROGENEOUS STRUCTURE OF FAULT ZONES WITHIN CARBONATE ROCKS: EVIDENCE FROM OUTCROP STUDIES AND IMPLICATIONS FOR FLUID FLOW C.G. Bonson*, J.J. Walsh, C. Childs, M.P.J. Schöpfer & V. Carboni Fault
More informationEarth models for early exploration stages
ANNUAL MEETING MASTER OF PETROLEUM ENGINEERING Earth models for early exploration stages Ângela Pereira PhD student angela.pereira@tecnico.ulisboa.pt 3/May/2016 Instituto Superior Técnico 1 Outline Motivation
More informationSound Recognition in Mixtures
Sound Recognition in Mixtures Juhan Nam, Gautham J. Mysore 2, and Paris Smaragdis 2,3 Center for Computer Research in Music and Acoustics, Stanford University, 2 Advanced Technology Labs, Adobe Systems
More informationKurt Marfurt Arnaud Huck THE ADVANCED SEISMIC ATTRIBUTES ANALYSIS
The Society of Exploration Geophysicists and GeoNeurale announce Kurt Marfurt Arnaud Huck THE ADVANCED SEISMIC ATTRIBUTES ANALYSIS 3D Seismic Attributes for Prospect Identification and Reservoir Characterization
More informationWe Prediction of Geological Characteristic Using Gaussian Mixture Model
We-07-06 Prediction of Geological Characteristic Using Gaussian Mixture Model L. Li* (BGP,CNPC), Z.H. Wan (BGP,CNPC), S.F. Zhan (BGP,CNPC), C.F. Tao (BGP,CNPC) & X.H. Ran (BGP,CNPC) SUMMARY The multi-attribute
More informationAppendix B: Geophysical Data (Thesis Appendix, 2013)
Utah State University From the SelectedWorks of David J Richey 2013 Appendix B: Geophysical Data (Thesis Appendix, 2013) David J Richey, Utah State University Available at: https://works.bepress.com/david_richey/2/
More informationMultiple-Point Geostatistics: from Theory to Practice Sebastien Strebelle 1
Multiple-Point Geostatistics: from Theory to Practice Sebastien Strebelle 1 Abstract The limitations of variogram-based simulation programs to model complex, yet fairly common, geological elements, e.g.
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2017) 210, 184 195 GJI Marine geosciences and applied geophysics doi: 10.1093/gji/ggx150 Structure-, stratigraphy- and fault-guided regularization in
More informationUsing structural validation and balancing tools to aid interpretation
Using structural validation and balancing tools to aid interpretation Creating a balanced interpretation is the first step in reducing the uncertainty in your geological model. Balancing is based on the
More informationDATA ANALYSIS AND INTERPRETATION
III. DATA ANALYSIS AND INTERPRETATION 3.1. Rift Geometry Identification Based on recent analysis of modern and ancient rifts, many previous workers concluded that the basic structural unit of continental
More informationStatics preserving projection filtering Yann Traonmilin*and Necati Gulunay, CGGVeritas
Yann Traonmilin*and Necati Gulunay, CGGVeritas Summary Projection filtering has been used for many years in seismic processing as a tool to extract a meaningful signal out of noisy data. We show that its
More informationScale Space Smoothing, Image Feature Extraction and Bessel Filters
Scale Space Smoothing, Image Feature Extraction and Bessel Filters Sasan Mahmoodi and Steve Gunn School of Electronics and Computer Science, Building 1, Southampton University, Southampton, SO17 1BJ, UK
More informationAmplitude variation with offset AVO. and. Direct Hydrocarbon Indicators DHI. Reflection at vertical incidence. Reflection at oblique incidence
Amplitude variation with offset AVO and Direct Hydrocarbon Indicators DHI Reflection at vertical incidence Reflection coefficient R(p) c α 1 S wavespeed β 1 density ρ 1 α 2 S wavespeed β 2 density ρ 2
More informationMachine Learning Applied to 3-D Reservoir Simulation
Machine Learning Applied to 3-D Reservoir Simulation Marco A. Cardoso 1 Introduction The optimization of subsurface flow processes is important for many applications including oil field operations and
More informationDownloaded 10/25/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
Facies modeling in unconventional reservoirs using seismic derived facies probabilities Reinaldo J. Michelena*, Omar G. Angola, and Kevin S. Godbey, ireservoir.com, Inc. Summary We present in this paper
More informationDetermine the azimuths of conjugate fracture trends in the subsurface
Reconnaissance of geological prospectivity and reservoir characterization using multiple seismic attributes on 3-D surveys: an example from hydrothermal dolomite, Devonian Slave Point Formation, northeast
More informationSimulated Annealing for Constrained Global Optimization
Monte Carlo Methods for Computation and Optimization Final Presentation Simulated Annealing for Constrained Global Optimization H. Edwin Romeijn & Robert L.Smith (1994) Presented by Ariel Schwartz Objective
More informationAutomatic time picking and velocity determination on full waveform sonic well logs
Automatic time picking and velocity determination on full waveform sonic well logs Lejia Han, Joe Wong, John C. Bancroft, and Robert R. Stewart ABSTRACT Full waveform sonic logging is used to determine
More informationLearning Gaussian Process Models from Uncertain Data
Learning Gaussian Process Models from Uncertain Data Patrick Dallaire, Camille Besse, and Brahim Chaib-draa DAMAS Laboratory, Computer Science & Software Engineering Department, Laval University, Canada
More information