Recent developments in curvelet-based seismic processing
|
|
- Asher Osborne
- 5 years ago
- Views:
Transcription
1 0000 Recent developments in curvelet-based seismic processing Felix J. Herrmann Seismic Laboratory for Imaging and Modeling Department of Earth and Ocean Sciences The University of British Columbia SUMMARY In this paper, we present recent developments in nonlinear curvelet-based sparsity-promoting formulations of problems in the seismic data processing flow. We present our latest work on a parallel curvelet transform and recent work on a curvelet-regularized formulation for the focal transform, the prediction of multiples and the computation of the inverse data space. We show that the curvelet s wavefront detection capability and invariance under wave propagation lead to a formulation of these problems that is stable under noise and missing data.
2 OUTLINE During this presentation an overview will be given of the latest developments concerning the application of the curvelet transform to seismic data processing. The following topics will be discussed extensively: the definition of a windowed parallel curvelet transform, which allows for the application of the curvelet transform to large data volumes (Thomson et al., 2006); combination of the curvelet transform with the (de)focal transform (= a data adaptive transform defined by a cross-convolution/correlation with the primaries) for the purpose of improved interpolation and multiple prediction (Berkhout and Verschuur, 2006) a curvelet-based formulation for the computation of the data inverse (Berkhout, 2006). Parallel Windowed FDCT Seismic data volumes are often too large to fit into the memory of single processor. At the expense of a low-frequency cutoff, a parallel curvelet transform can be constructed by applying a domain decomposition (see Fig. 1), implementing a partitioning of unity (the squares of the windows sum to 1). Consider equally sized windows with uniform overlap ɛ. The overlapping regions are tapered such that the energy of the system remains unity when points in the data set are duplicated at the overlaps. The tapering ensures that edges of the data goes smoothly to zero, eliminating potential edge artifacts. Once the data is split into overlapping windows and tapered the fast discrete curvelet transform (FDCT, Candes et al., 2006; Hennenfent and Herrmann, 2006), is performed on each window separately. The shape of the taper function and the overlap of neighboring windows are both shown in Fig. 1. The data in this example is split into sixteen overlapping windows. The dashed lines in the image represent the edges of the overlapping windows, with the region between nearby parallel lines being shared amongst the two windows. The taper function is shown, where one can see the value going to zero at the window edges.
3 THEORY This entire process of windowing can be considered as a linear operator acting on a vector containing the reordered to-be-transformed data volumes. For the parallel transform, it is useful to distinguish between global vectors and transforms (=matrix-vector multiplication), which comprise the entire data volume, and local vectors and transforms that exist only on individual nodes. Our notation will reflect this distinction. For instance, the distributed transform and vector C and x are distinguished from the local transform and vector C i and x i that exist on a node indexed by i. Note that there will be cases when we want to apply a local transform C i to the data on every node in a global vector x. This involves a block diagonal global operator where each block consists of the particular C i corresponding to a given x i, which is simply a subsection of x. We now consider the entire process in the context of linear operators We denote the block diagonal acting data vectors. global In the paralleloperator context, it is useful tofor distinguish a local operator as [C i ]. nodes. Our notation will reflect this distinction. For instance, the distributed data operator and vector A and x are expanded easily distinguished from into the overlapping sections An arbitrary global local operator and vector A i and x i that exist on a node indexed by i. using the global windowing operator W. The tapering operator T is then Note that there will be cases when we want to apply a local operator A i to the data on every node in a global vector x. This involves a block diagonal global operator where each block consists of the particular applied across the system. A i corresponding It tois a given easy x i, which is tosimply see a subsection that of x. We the tapering operator is denote the block diagonal global operator for a local operator as [A i]. diagonal, and will have An arbitrary repeating global data vector ispatterns expanded into overlapping in sections blocks representing each using the global windowing operator W. The tapering operator T is then applied across the system. It is easy to see that the tapering opera- node. These tapering operators is diagonal, and will have T i repeating arepatterns local in blockson representing each node. Once W and each node. Alternatively, one can simply look at the tapering operator T have been applied, the as a local linear T operator C i can be applied on each node. i on each node. Once W and T have been applied, any We propose a method to parallelize an arbitrary linear operator that avoids problems related to edge artifacts and preserves overall energy. It requires relatively little communication between parallel nodes, making it highly scalable. Our particular interest here is focused on a scalable parallel generalization of the FDCT, but we stress the fact that any linear operator could take the place of the FDCT in this framework. Parallel Windowed FDCT The basic structures of this framework are overlapping and tapered windows. This scheme has previously been used in various parallel processing applications (see e.g. J. Fan and Rector, 1997). The details of these structures are quite flexible. In general, the sizes of the windows and the overlapping regions do not have to be uniform throughout the system. Here, though, we will only consider equally sized windows with uniform overlap throughout the system. One can express the overlap between windows by the value ε, which represents the depth to which one window receives adjacent data from its neighbours (Mallat, 1998), as illustrated in Figure 1. The overlapping regions of the windows are tapered such that the energy of the system remains constant when points in the data set are duplicated due to the overlaps. The tapering also ensures that edges of the data goes smoothly to zero at the edges, eliminating the potential of creating edge artifacts. The taper is applied across the outer region of the overlapping windows, affecting points that are within a distance of 2ε from an edge. The tapering function b must satisfy the relationship b 2 + b 2 = 1, (1) where b and b signify the tapers in adjacent windows that cover the same region. For consistency, we consider the application of an identical tapering function b to all windows. It follows from Eq. 1 that b 2 n + b 2 2ε n+1 = 1, (2) where n is an integer on the interval [1,2ε], and subject to the boundary conditions b 1 = 0 and b 2ε = 1. This condition ensures that the energy of the system is conserved when summing values from adjacent windows that represent the data at the same location. There are many examples of functions that satisfy this relationship (Mallat, 1998), the simplest being: ( ) (n 1)π b n = sin. (3) 2(2ε 1) To taper the end of the data set along each dimension, the function b is translated to the points n = {N,N 1,, N 2ε +1}, where N is the total size of the overlapping windows. The tapering function is then ( ) (N n)π b n = sin 2(2ε 1). (4) Once the data is split into overlapping windows and tapered, operations, like the FDCT in our particular case, can be performed on each window independently. The shape of the taper function and the overlap of neighbouring windows are both shown in Figure 1. The data in this case is split into sixteen overlapping windows, four horizontally and four vertically. The dashed lines in the image represent the edges of the overlapping windows, with the region between nearby parallel lines being shared between the two windows. The taper function is shown, where one can see the value going to zero at the window edges. Since the image and all windows are square in this case, the taper functions along the vertical axis are the same as the horizontal ones that are shown. Also, note that tapering is not done at the edges of the image, since this would remove energy from the system. This issue can be successfully dealt with in a number of ways, but we omit this discussion. between global vectors and operators, which comprise the entire parallel system, and local vectors and operators that exist only on individual arbitrary linear operator A i can be applied on each node. Figure 1: 2D synthetic seismic data with vertical traces. Figure 1: 2D synthetic seismic data with 30% missing vertical traces. Dotted lines Dotted lines represent borders of overlapping windows. Taper function is shown above for clarity. represent borders of overlapping windows. Taper function is shown above for clarity. Perhaps the most useful part of considering windowing and tapering as linear operators in matrix form is that looking at these matrices makes understanding the adjoints of these operations simple. Since T is diagonal, it is its own adjoint. For W, the adjoint operation involves summing together overlapping regions. Since the forward operation was a scatter, the adjoint becomes a gather, where data that correspond to the same point are summed (Claerbout, 1992). In a parallel computing realization, this means that a given node sends its outer band to the nodes to which they belong, then gathers data related to its inner band from those same neigbours, and sums them together. Importantly, the combination of these processes satisfies the relation Since T is diagonal, it is its own adjoint. For W, the adjoint operation involves summing together overlapping regions. Since the forward operation was a scatter, the adjoint becomes a gather, where data that correspond to the same point are summed (Claerbout, 1992). In a parallel computing realization, this means that a given node sends its outer W H T H TW = I, (5) band to the nodes to which they belong, then gathers data related to its
4 inner band from those same neigbours, and sums them together. The combination of these processes satisfies the relation W H T H TW = I, (1) which ensures perfect reconstruction of the data after applying the operators followed by their adjoints, i.e. our parallel curvelet transform is a tight frame. This implies that energy is preserved through the entire process. So, a new transform, that we will call the Parallel Windowed FDCT (PWFDCT), can thus be defined. The parallel curvelet transform corresponds to a block-diagonal global operator that is defined by applying the FDCT, C i, on each node independently after applying the global windowing and tapering operators T and W. In other words, C := [C i ]TW. (2) Since the local C is numerically tight, it follows that C H C = I. (3) Other properties of C are similarly shared by C. It is possible, then, to make use of the PWFDCT in existing algorithms that include the FDCT. It should be noted, though, that curvelets at very large scales (or, equivalently, low frequencies) are not represented in the same way they would be if a single FDCT were performed on the global data. Since the benefits of curvelets are mostly found at the finer scales (higher frequencies), the lack of large scale curvelet representation is not a problem. This makes a wide variety of algorithms capable of handling data sets much larger than otherwise possible. An example of a seismic data interpolation algorithm where the PWFDCT takes the place of the FDCT is included in Fig. 2. From the results presented in that figure one can conclude that algorithms, such as curvelet-based sparsity-promoting data recovery, can readily be scaled up by using our parallel windowed curvelet transform based on a partition of unity. EXTENSIONS OF THE CURVELET TRANSFORM So far, curvelets have mainly found their applications in the approximation of migration operators (Douma and de Hoop, 2006; Chauris, 2006), amplitude recovery (Herrmann et al., 2006) and seismic data interpolation (Herrmann and Hennenfent, 2007) and separation (Herrmann
5 evident near window borders, as would be expected. In Figure 1, we show a subsection of the full output of interpolation runs using overlapping and non-overlapping windows. The same subsection is used for both, and contains one full window and parts of its nearest neighbours. In the example using non-overlapping windows (Figure 2(a)), the edges of the windows are obvious from the artifacts that appear. When overlapping windows are used (Figure 2(b)), the artifacts are no longer evident, and it is not clear at all where the window borders are. In essence, the interpolation performs just as well in the proximity of the window edges as it does in the middle of the window, which is clearly untrue for non-overlapping windows. The importance of overlapping and tapering is thus, as expected, clear from this example. We could potentially hold in memory. In order to run many algorithms in parallel, it is not sufficient to process data in separate pieces. At the same time, these same algorithms are often not scalable in their normal form due to exponential growth in the amount of communication between nodes that they require. For these reasons, we have developed a scheme that involves overlapping, tapered data windows that can be processed in parallel that is highly scalable since the communication costs do not grow with the number of processing nodes. We have used this method to define the PWFDCT, but we stress that this method is general and that the FDCT can be replaced by an arbitrary operator acting on each overlapping window independently as desired. We applied the PWFDCT to a seismic (a) data interpolation algorithm that is shown, in another presentation to this conference, to be successful using the FDCT. We have demonstrated that good results can be achieved with the PWFDCT in this algorithm, and shown the importance of the overlapping and tapering. Furthermore, we have demonstrated that it is possible to correct for the effect of tapering on threshold values in the transform domain. Much future work is expected to arise from the ideas described herein. In particular, we will apply the PWFDCT to other seismic data processing algorithms, enabling them to work on much larger data sets than is currently feasible. The parallelization method described here will also be applied to other operators besides the FDCT. Finally, we hope that this method will open up new possibilities in parallel signal processing in a variety of fields. ACKNOWLEDGMENTS The authors would like to thank the authors of the Fast Discrete Curvele Transform (FDCT) for making this code available at This work was in part financially supported by NSERC Discovery Grant 22R81254 of Felix J. Herrmann and was carried out as part of the SINBAD project with support, secured through ITF (the Industry Technology Facilitator), from the following organizations: BG Group, BP, Chevron, ExxonMobil and Shell. Much future work is expected to arise from the ideas described herein. In particular, we will apply the PWFDCT to other seismic data processing algorithms, enabling them to work on much larger data sets REFERENCES than is currently feasible. The parallelization method described here will also be applied to other operators besides the FDCT. Finally, we Claerbout, J. F., 1992, Earth soundings analysis: Processing versus hope that this method will open up new possibilities in parallel signal inversion: Blackwell Scientific publishing. processing in a variety of fields. Daubechies, I., M. Defrise, and C. de Mol, 2004, An iterative thresholding algorithm for linear inverse problems with a sparsity constraint: Comm. Pure Appl. Math., ACKNOWLEDGMENTS Gupta, A. and V. Kumar, 1993, The scalability of FFT on parallel computers: IEEE Transaction on Parallel and Distributed Systems. (a) (b) The authors would like to thank the authors of the Fast Discrete Curvelet Hennenfent, G. and F. J. Herrmann, 2006, Application of stable seismic signal recovery to seismic interpolation. (submitted for pre- (a) Transform (FDCT) for making(b) this code available at Figure 2: Interpolation output for (a) non-overlapping windows and This work was in part financially supported by NSERC Discovery sentation at 76th SEG Conference & Exhibition). for (b) overlapping, tapered windows with threshold correction. Arrows indicate locations of artifacts in (a) and show the improvement in Grant 22R81254 of Felix J. Herrmann and was carried out as part of Herrmann, F. J. and G. Hennenfent, 2006, Non-parametric seismic the SINBAD project with support, secured through ITF (the Industry data recovery with curvelet frames. (to be submitted). Figure 2: Interpolation output for (a) non-overlapping (b). Technology Facilitator), windows from the following and organizations: for (b) BG Group, overlapping, J. Fan, K.T. Nihei, L. M. N. C. and J. Rector, 1997, Overlap domain BP, Chevron, ExxonMobil and Shell. decomposition method for wave propagation: Presented at the 67th also compared different methods for correcting the threshold value in SEG Conference & Exhibition. tapered windows with threshold correction. Arrows indicate locations of artifacts in (a) the curvelet domain to account for the taper. When the taper was not L. Ying, L. D. and E. Candès, 2005, 3D discrete curvelet transform. taken REFERENCES into account in thresholding, errors were evident in the overlapping region. When correcting threshold values by using a Monte Carlo submitted for publication. and show the improvement in (b). Mallat, S., 1998, A wavelet tour of signal processing: Academic Press. sampling Claerbout, or by J. F., evaluating 1992, Earth the taper soundings function analysis: at curvelet Processing centroids, versus artifacts inversion: related toblackwell erroneousscientific thresholding publishing. are eliminated. We omit the figures Daubechies, demonstrating I., M. Defrise, this dueand to space C. de limitations. Mol, 2004, An iterative thresholding algorithm for linear inverse problems with a sparsity constraint: Comm. Pure Appl. Math., et al., 2007). During these applications, Gupta, A. andthe V. Kumar, invariance 1993, The scalability of FFTand on parallel computers: IEEE Transaction on Parallel and Distributed Systems. compressibility of curvelets is used. The curvelet mic signal properties recovery to seismic interpolation. can (submitted alsoforbe pre- used to (b) Hennenfent, G. and F. J. Herrmann, 2006, Application stable seis- Figure 2: Interpolation output for (a) non-overlapping windows and sentation at 76th SEG Conference & Exhibition). for (b) overlapping, tapered windows with threshold correction. Arrows seismic indicate locations of data artifacts in (a) recovery and show the improvement with in the data recovery focal with curvelet transform frames. (to be submitted). (Berkhout and Herrmann, F. J. and G. Hennenfent, 2006, Non-parametric seismic improve (b). J. Fan, K.T. Nihei, L. M. N. C. and J. Rector, 1997, Overlap domain Verschuur, 2006), multiple prediction decomposition method for wave propagation: Presented the 67th also compared different methods for correcting the threshold value in SEG and Conference the & Exhibition. computation of data inverse (Berkhout, taken into account in thresholding, 2006). errors werein evidenthis in the overlap- case, use is made of the property that submitted for publication. the curvelet domain to account for the taper. When the taper was not L. Ying, L. D. and E. Candès, 2005, 3D discrete curvelet transform. ping region. When correcting threshold values by using a Monte Carlo Mallat, S., 1998, A wavelet tour of signal processing: Academic Press. sampling or by evaluating the taper function curvelet centroids, ar- remain related to erroneous relatively thresholding are eliminated. sparse We omit theon data that is focused, defocused or curveletstifacts figures demonstrating this due to space limitations. inverted. Focused data in this context corresponds to data stripped of one propagation path interacting with the surface, while defocused data has a single propagation path added, mapping primaries into first-order multiples. The latter corresponds to SRME multiple predicted (Verschuur and Berkhout, 1997), while the former corresponds to the recently introduced focal transform (Berkhout and Verschuur, 2006). In all above described manifestations of seismic data, we are dealing with data that contains wavefronts whether these are primaries, multiples or the data inverse. Curvelets compress data with wavefronts and therefore serve as a regularization, where the primary operator, adjoint primary operator and data operator can stably be inverted with a curvelet sparsity norm. These different data manifestations can be calculated by solving of the following norm-one nonlinear program: P ɛ : { x = arg minx x 1 s.t. Ax y 2 ɛ f = ST x (4) in which y is the (incomplete) data, A the synthesis matrix and S T is the inverse sparsity transform and ɛ, a noise-dependent tolerance level. This constrained optimization problem is solved to within ɛ.
6 Recovery from incomplete data by focusing Combination of the non-adaptive curvelet transform with the data-adaptive focal transform (Berkhout and Verschuur, 2006) gives rise to an interesting extension, where data is focused by inverting the primary operator (= a multidimensional convolution with the primaries). In the focal domain, data is sparser and hence the curvelet coefficients are sparser, which leads to an improved recovery. Recovery with focusing is possible by defining the synthesis matrix as A := R PC T and the inverse sparsity transform as S T := PC T. The operator P stands for a multidimensional convolution with the primaries. The data vector y contains the incomplete measurements, y = Rd with R the restriction matrix that selects the rows from the curvelet matrix that correspond to active traces. The solution yielded by P ɛ corresponds to a curvelet-regularized inverse of the primary operator. The estimated coefficients represent an estimate for the focused data. The synthesis operator defocuses the focused data back to data domain. The results as shown Fig 3 show as expected an improvement over interpolation without focusing. Multiple prediction with defocusing By defining the synthesis matrix by the adjoint of the primary operator, the solution of P ɛ, with the data y containing the data including multiples, yields an estimate for the predicted multiples, i.e. by setting S T = C T. This estimate corresponds to the inverse of the adjoint of the primary operator, i.e., A := P T C T. As can be observed from Fig. 4, the frequency content of the predicted multiples is increased by this prediction through curvelet-regularized inversion. Computation of the data inverse The previous two examples showed that a curvelet regularization for the inversion of the primary operator and its inverse improves the stability and yields better results for the multiple prediction. Similarly, the computation of the data inverse can benefit from the same regularization by solving P ɛ for the synthesis matrix, A := PC T with P the multidimensional convolution with the total data, and the data vector, y = δ, with δ the discrete delta-function. As the example in Fig. 5 illustrates, the
7 (a) (b) (c) Figure 3: Comparison between curvelet-based recovery by sparsity-promoting inversion with and without focusing. (a) Randomly subsampled data with 80 % of the traces missing. (b) Curvelet-based recovery. (c) Curvelet-based recovery with focusing. Notice the significant improvement from the focusing with the primary operator. SMRE predicted multiples CSRI predicted multiples (a) (b) Figure 4: Comparison between convolution-based multiple prediction (a) and sparsity-based multiple prediction (b).
8 (a) (b) Figure 5: Computation of the data inverse. (a) shot record in the data space. (b) shot record in the inverse data space computed according to P ɛ. solution of P ɛ yields a stable estimate for the data inverse. DISCUSSION AND CONCLUSIONS Applications of the curvelet transform to seismic data processing bank on two favorable properties of curvelets, namely their ability to detect wavefronts and their approximate invariance under wave propagation. In this paper, we showed some exciting recent developments of this transform towards a large-scale parallelization and a combination with the physics of wave propagation. The successful application of curvelets, juxtaposed by sparsity-promoting inversion, opens a range of new perspectives on seismic data processing, wavefield extrapolation and imaging. Because of their singular wavefront detection capability, curvelets represent in our vision an ideal domain for future developments in exploration seismology. ACKNOWLEDGMENTS: The author would like to thank the authors of CurveLab for making their codes available, Darren Thomson for helping to implement the parallel curvelet transform and Eric Verschuur for his input. The examples presented in this paper were prepared with Madagascar (rsf.sourceforge. net/). SAGA and ExxonMobil are thanked for making test datasets available. This work was in part financially supported by the Natural Sciences and Engineering Research Council of Canada Discovery
9 Grant (22R81254) and Collaborative Research and Development Grant DNOISE ( ) of Felix J. Herrmann and was carried out as part of the SINBAD project with support, secured through ITF (the Industry Technology Facilitator), from the following organizations: BG Group, BP, Chevron, ExxonMobil and Shell. REFERENCES Berkhout, A. J., 2006, Seismic processing in the inverse data space: Geophysics, 71. Berkhout, A. J. and D. J. Verschuur, 2006, Focal transformation, an imaging concept for signal restoration and noise removal: Geophysics, 71. Candes, E. J., L. Demanet, D. L. Donoho, and L. Ying, 2006, Fast discrete curvelet transforms: SIAM Multiscale Model. Simul., 5, Chauris, H., 2006, Seismic imaging in the curvelet domain and its implications for the curvelet design: Presented at the 76th Ann. Internat. Mtg., SEG, Soc. Expl. Geophys., Expanded abstracts. Claerbout, J. F., 1992, Earth soundings analysis: Processing versus inversion: Blackwell Scientific publishing. Douma, H. and M. de Hoop, 2006, Leading-order seismic imaging using curvelets: Presented at the 76th Ann. Internat. Mtg., SEG, Soc. Expl. Geophys., Expanded abstracts. Hennenfent, G. and F. J. Herrmann, 2006, Seismic denoising with nonuniformly sampled curvelets: IEEE Comp. in Sci. and Eng., 8, Herrmann, F. J., U. Boeniger, and D.-J. E. Verschuur, 2007, Nonlinear primary-multiple separation with directional curvelet frames: Geoph. J. Int. To appear. Herrmann, F. J. and G. Hennenfent, 2007, Non-parametric seismic data recovery with curvelet frames. Submitted for publication. Herrmann, F. J., P. P. Moghaddam, and C. Stolk, 2006, Sparsety- and continuity-promoting seismic imaging with curvelet frames. In revision. Thomson, D., G. Hennenfent, H. Modzelewski, and F. Herrmann, 2006, A parallel windowed fast discrete curvelet transform applied to seismic processing: Presented at the SEG International Exposition and 76th Annual Meeting. Verschuur, D. J. and A. J. Berkhout, 1997, Estimation of multiple scat-
10 tering by iterative inversion, part II: practical aspects and examples: Geophysics, 62,
Incoherent noise suppression with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Incoherent noise suppression with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY The separation of signal and noise is a key issue in seismic data processing. By
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 informationSUMMARY. forward/inverse discrete curvelet transform matrices (defined by the fast discrete curvelet transform, FDCT, with wrapping
Seismic data processing with curvelets: a multiscale and nonlinear approach Felix J. Herrmann, EOS-UBC, Deli Wang, Jilin University and Gilles Hennenfent and Peyman Moghaddam SUMMARY In this abstract,
More informationSeismic wavefield inversion with curvelet-domain sparsity promotion Felix J. Herrmann, EOS-UBC and Deli Wang, Jilin University
Seismic wavefield inversion with curvelet-domain sparsity promotion Felix J. Herrmann, EOS-UBC and Deli Wang, Jilin University SUMMARY Inverting seismic wavefields lies at the heart of seismic data processing
More informationCurvelet-domain matched filtering Felix J. Herrmann and Peyman P. Moghaddam EOS-UBC, and Deli Wang, Jilin University
Curvelet-domain matched filtering Felix J. Herrmann and Peyman P. Moghaddam EOS-UBC, and Deli Wang, Jilin University SUMMARY Matching seismic wavefields and images lies at the heart of many pre-/post-processing
More informationRobust Seismic Image Amplitude Recovery Using Curvelets
Robust Seismic Image Amplitude Recovery Using Curvelets Peyman P. Moghaddam Joint work with Felix J. Herrmann and Christiaan C. Stolk UNIVERSITY OF BRITISH COLUMBIA University of Twente References Curvelets
More informationRandomized sampling strategies
Randomized sampling strategies Felix J. Herrmann SLIM Seismic Laboratory for Imaging and Modeling the University of British Columbia SLIM Drivers & Acquisition costs impediments Full-waveform inversion
More informationSEG/New Orleans 2006 Annual Meeting
and its implications for the curvelet design Hervé Chauris, Ecole des Mines de Paris Summary This paper is a first attempt towards the migration of seismic data in the curvelet domain for heterogeneous
More informationMain Menu. Summary. sampled) f has a sparse representation transform domain S with. in certain. f S x, the relation becomes
Preliminary study on Dreamlet based compressive sensing data recovery Ru-Shan Wu*, Yu Geng 1 and Lingling Ye, Modeling and Imaging Lab, Earth & Planetary Sciences/IGPP, University of California, Santa
More informationSUMMARY. In combination with compressive sensing, a successful reconstruction
Higher dimensional blue-noise sampling schemes for curvelet-based seismic data recovery Gang Tang, Tsinghua University & UBC-Seismic Laboratory for Imaging and Modeling (UBC-SLIM), Reza Shahidi, UBC-SLIM,
More informationSeismic Data Interpolation With Symmetry
Seismic Data Interpolation With Symmetry James Johnson joint work with Gilles Hennenfent 8 Consortium Meeting Outline Introduction Reciprocity Real Data Example Data Interpolation Problem Statement CRSI
More informationG009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain
G009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain M. Naghizadeh* (University of Alberta) & M. Sacchi (University of Alberta) SUMMARY We propose a robust interpolation
More informationTime-jittered ocean bottom seismic acquisition Haneet Wason and Felix J. Herrmann
Time-jittered ocean bottom seismic acquisition Haneet Wason and Felix J. Herrmann SLIM University of British Columbia Challenges Need for full sampling - wave-equation based inversion (RTM & FWI) - SRME/EPSI
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 informationCurvelet-based non-linear adaptive subtraction with sparseness constraints. Felix J Herrmann, Peyman P Moghaddam (EOS-UBC)
Curvelet-based non-linear adaptive subtraction with sparseness constraints Felix J Herrmann, Peyman P Moghaddam (EOS-UBC) Context Spiky/minimal structure deconvolution (Claerbout, Ulrych, Oldenburg, Sacchi,
More informationSUMMARY. Pursuit De-Noise (BPDN) problem, Chen et al., 2001; van den Berg and Friedlander, 2008):
Controlling linearization errors in l 1 regularized inversion by rerandomization Ning Tu, Xiang Li and Felix J. Herrmann, Earth and Ocean Sciences, University of British Columbia SUMMARY Linearized inversion
More informationOne-norm regularized inversion: learning from the Pareto curve
One-norm regularized inversion: learning from the Pareto curve Gilles Hennenfent and Felix J. Herrmann, Earth & Ocean Sciences Dept., the University of British Columbia ABSTRACT Geophysical inverse problems
More informationAmplitude and kinematic corrections of migrated images for non-unitary imaging operators
Stanford Exploration Project, Report 113, July 8, 2003, pages 349 363 Amplitude and kinematic corrections of migrated images for non-unitary imaging operators Antoine Guitton 1 ABSTRACT Obtaining true-amplitude
More informationCurvelet-domain multiple elimination with sparseness constraints. Felix J Herrmann (EOS-UBC) and Eric J. Verschuur (Delphi, TUD)
Curvelet-domain multiple elimination with sparseness constraints Felix J Herrmann (EOS-UBC) and Eric J. Verschuur (Delphi, TUD) Context SRME (Verschuur, Guitton, Berkhout, Weglein, Innanen) Sparse Radon
More informationInterpolation with pseudo-primaries: revisited
Stanford Exploration Project, Report 129, May 6, 2007, pages 85 94 Interpolation with pseudo-primaries: revisited William Curry ABSTRACT Large gaps may exist in marine data at near offsets. I generate
More informationA new sparse representation of seismic data using adaptive easy-path wavelet transform
A new sparse representation of seismic data using adaptive easy-path wavelet transform Jianwei Ma 1,3, Gerlind Plonka 2, Hervé Chauris 3 1 Institute of Seismic Exploration, School of Aerospace, Tsinghua
More informationSEG Houston 2009 International Exposition and Annual Meeting
Yu Geng* 1, Ru-Shan Wu and Jinghuai Gao 2 Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz, CA 95064 Summary Local cosine/sine basis is a localized version of cosine/sine basis
More informationSub-Nyquist sampling and sparsity: getting more information from fewer samples
Sub-Nyquist sampling and sparsity: getting more information from fewer samples Felix J. Herrmann SLIM Seismic Laboratory for Imaging and Modeling the University of British Columbia SLIM Drivers Our incessant
More informationRandomized sampling without repetition in timelapse seismic surveys
Randomized sampling without repetition in timelapse seismic surveys Felix Oghenekohwo SLIM University of British Columbia Randomized sampling without repetition in timelapse seismic surveys Felix Oghenekohwo
More informationSeismic data interpolation beyond aliasing using regularized nonstationary autoregression a
Seismic data interpolation beyond aliasing using regularized nonstationary autoregression a a Published in Geophysics, 76, V69-V77, (2011) Yang Liu, Sergey Fomel ABSTRACT Seismic data are often inadequately
More informationSource Estimation for Wavefield Reconstruction Inversion
Source Estimation for Wavefield Reconstruction Inversion Zhilong Fang * and Felix J. Herrmann * * Seismic Laboratory for Imaging and Modeling SLIM), University of British Columbia Abstract Wavefield reconstruction
More informationSeismic data interpolation using nonlinear shaping regularization a
Seismic data interpolation using nonlinear shaping regularization a a Published in Journal of Seismic Exploration, 24, no. 5, 327-342, (2015) Yangkang Chen, Lele Zhang and Le-wei Mo ABSTRACT Seismic data
More informationWave-equation inversion prestack Hessian
Stanford Exploration Project, Report 125, January 16, 2007, pages 201 209 Wave-equation inversion prestack Hessian Alejandro A. Valenciano and Biondo Biondi ABSTRACT The angle-domain Hessian can be computed
More informationLeast squares Kirchhoff depth migration: important details
Least squares Kirchhoff depth migration: important details Daniel Trad CREWES-University of Calgary Summary Least squares migration has been an important research topic in the academia for about two decades,
More informationTarget-oriented wave-equation inversion
Stanford Exploration Project, Report 120, May 3, 2005, pages 23 40 Target-oriented wave-equation inversion Alejandro A. Valenciano, Biondo Biondi, and Antoine Guitton 1 ABSTRACT A target-oriented strategy
More informationPractical implementation of SRME for land multiple attenuation
Practical implementation of SRME for land multiple attenuation Juefu Wang* and Shaowu Wang, CGGVeritas, Calgary, Canada juefu.wang@cggveritas.com Summary We present a practical implementation of Surface
More informationAngle-dependent reflectivity by wavefront synthesis imaging
Stanford Exploration Project, Report 80, May 15, 2001, pages 1 477 Angle-dependent reflectivity by wavefront synthesis imaging Jun Ji 1 ABSTRACT Elsewhere in this report, Ji and Palacharla (1994) show
More informationData dependent parameterization and covariance calculation for inversion of focusing operators
Stanford Exploration Project, Report 11, September 18, 21, pages 1 91 Data dependent parameterization and covariance calculation for inversion of focusing operators Barbara E. Cox 1 ABSTRACT The Common
More informationHierarchical scale curvelet interpolation of aliased seismic data Mostafa Naghizadeh and Mauricio Sacchi
Hierarchical scale curvelet interpolation of aliased seismic data Mostafa Naghizadeh and Mauricio Sacchi SUMMARY We propose a robust interpolation scheme for aliased regularly sampled seismic data that
More informationSUMMARY. In this paper, we present a methodology to improve trace interpolation
Reconstruction of seismic wavefields via low-rank matrix factorization in the hierarchical-separable matrix representation Rajiv Kumar, Hassan Mansour, Aleksandr Y. Aravkin, and Felix J. Herrmann University
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 informationSparse wavelet expansions for seismic tomography: Methods and algorithms
Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves 24 28 July 2011 (Joint
More informationCompressed sensing based land simultaneous acquisition using encoded sweeps
Compressed sensing based land simultaneous acquisition using encoded sweeps Rajiv Kumar 1, Shashin Sharan 1, Nick Moldoveanu 2, and Felix J. Herrmann 3 1 University of British Columbia, Vancouver 2 WesternGeco,
More informationShort Note. Non-stationary PEFs and large gaps. William Curry 1 INTRODUCTION
Stanford Exploration Project, Report 120, May 3, 2005, pages 247 257 Short Note Non-stationary PEFs and large gaps William Curry 1 INTRODUCTION Prediction-error filters (PEFs) may be used to interpolate
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 informationWriting Kirchhoff migration/modelling in a matrix form
Writing Kirchhoff migration/modelling in a matrix form Abdolnaser Yousefzadeh and John C. Bancroft Kirchhoff migration matrix ABSTRACT Kirchhoff prestack migration and modelling are linear operators. Therefore,
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 new take on FWI: Wavefield Reconstruction Inversion
A new take on FWI: Wavefield Reconstruction Inversion T. van Leeuwen 1, F.J. Herrmann and B. Peters 1 Centrum Wiskunde & Informatica, Amsterdam, The Netherlands University of British Columbia, dept. of
More informationFull-waveform inversion with randomized L1 recovery for the model updates
Full-waveform inversion with randomized L1 recovery for the model updates Xiang Li, Aleksandr Aravkin, Tristan van Leeuwen, Felix Herrmann March 5, 11 Abstract Full-waveform inversion (FWI) is a data fitting
More informationOptimized Compressed Sensing for Curvelet-based. Seismic Data Reconstruction
Optimized Compressed Sensing for Curvelet-based Seismic Data Reconstruction Wen Tang 1, Jianwei Ma 1, Felix J. Herrmann 2 1 Institute of Seismic Exploration, School of Aerospace, Tsinghua University, Beijing
More information= 0) is the 2-D Fourier transform of the field (1),, z = 0). The phase ) is defined in the dispersion relation as where
GEOPHYSICS, VOL. 61, NO. 5 (SEPTEMBER-OCTOBER 1996); P. 1412 1416, 4 FIGS. Short Note Prestack migration by split-step DSR. Alexander Mihai Popovici* INTRODUCTION The double-square-root (DSR) prestack
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 informationSeismic Noise Attenuation Using Curvelet Transform and Dip Map Data Structure
Seismic Noise Attenuation Using Curvelet Transform and Dip Map Data Structure T. Nguyen* (PGS), Y.J. Liu (PGS) Summary The curvelet transform is a known tool used in the attenuation of coherent and incoherent
More informationMissing trace interpolation and its enhancement of seismic processes
Missing trace interpolation Missing trace interpolation and its enhancement of seismic processes Wai-kin Chan and Robert R. Stewart ABSTRACT Many multi-channel seismic algorithms assume that the input
More information# of PDE solves flops conventional 2N s N s N hx N hy N hz this paper 2N x N r N s
AVA analysis and geological dip estimation via two-way wave-equation based extended images Rajiv Kumar, Tristan van Leeuwen and Felix J. Herrmann University of British Columbia, Department of Earth, Ocean
More informationSEG/San Antonio 2007 Annual Meeting
Imaging steep salt flanks by super-wide angle one-way method Xiaofeng Jia* and Ru-Shan Wu, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary The super-wide angle one-way
More informationSEG/New Orleans 2006 Annual Meeting
3-D tomographic updating with automatic volume-based picking Dimitri Bevc*, Moritz Fliedner, Joel VanderKwaak, 3DGeo Development Inc. Summary Whether refining seismic images to evaluate opportunities in
More informationApex Shifted Radon Transform
Apex Shifted Radon Transform Daniel Trad* Veritas DGC Inc., Calgary, Alberta, Canada dtrad@veritasdgc.com ABSTRACT The Apex Shifted Radon transform (ASRT) is an extension of the standard hyperbolic RT,
More informationApproximate Constant Density Acoustic Inverse Scattering Using Dip-Dependent Scaling Rami Nammour and William Symes, Rice University
Approximate Constant Density Acoustic Inverse Scattering Using Dip-Dependent Scaling Rami Nammour and William Symes, Rice University SUMMARY This abstract presents a computationally efficient method to
More informationExtending the search space of time-domain adjoint-state FWI with randomized implicit time shifts
Extending the search space of time-domain adjoint-state FWI with randomized implicit time shifts Mathias Louboutin * and Felix J. Herrmann Seismic Laboratory for Imaging and Modeling (SLIM), The University
More informationSummary. Introduction. Source separation method
Separability of simultaneous source data based on distribution of firing times Jinkun Cheng, Department of Physics, University of Alberta, jinkun@ualberta.ca Summary Simultaneous source acquisition has
More informationSeismic reflector characterization by a multiscale detection-estimation method
Seismic reflector characterization by a multiscale detection-estimation method Mohammad Maysami ABSTRACT Seismic reflector interfaces in the subsurface are typically idealized as zero-order discontinuities.
More informationEquivalent offset migration: the implementation and application update
Equivalent offset migration: the implementation and application update Xinxiang Li, Yong Xu and John C. Bancroft INTRODUCTION We have some improvements about equivalent offset migration (EOM) method. In
More informationData interpolation in pyramid domain
Data interpolation in pyramid domain Xukai Shen ABSTRACT Pyramid domain is defined as a frequency-space domain with different spatial grids for different frequencies. Data interpolation in pyramid domain
More informationLarge-scale workflows for wave-equation based inversion in Julia
Large-scale workflows for wave-equation based inversion in Julia Philipp A. Witte, Mathias Louboutin and Felix J. Herrmann SLIM University of British Columbia Motivation Use Geophysics to understand the
More informationStanford Exploration Project, Report 111, June 9, 2002, pages INTRODUCTION THEORY
Stanford Exploration Project, Report 111, June 9, 2002, pages 231 238 Short Note Speeding up wave equation migration Robert G. Clapp 1 INTRODUCTION Wave equation migration is gaining prominence over Kirchhoff
More informationSource-receiver migration of multiple reflections
Stanford Exploration Project, Report 113, July 8, 2003, pages 75 85 Source-receiver migration of multiple reflections Guojian Shan 1 ABSTRACT Multiple reflections are usually considered to be noise and
More informationAVA analysis and geological dip estimation via twoway wave-equation based extended images
AVA analysis and geological dip estimation via twoway wave-equation based extended images Rajiv Kumar, Tristan van Leeuwen and Felix J. Herrmann SLIM University of British Columbia Motivation Computation
More informationModel-space vs. data-space normalization for recursive depth migration
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Model-space vs. data-space normalization for recursive depth migration James Rickett 1 ABSTRACT Illumination problems caused by finite-recording
More informationTh C 02 Model-Based Surface Wave Analysis and Attenuation
Th C 02 Model-Based Surface Wave Analysis and Attenuation J. Bai* (Paradigm Geophysical), O. Yilmaz (Paradigm Geophysical) Summary Surface waves can significantly degrade overall data quality in seismic
More informationInterpolation using asymptote and apex shifted hyperbolic Radon transform
Interpolation using asymptote and apex shifted hyperbolic Radon transform Amr Ibrahim, Paolo Terenghi and Mauricio D. Sacchi Department of Physics, University of Alberta PGS Abstract The asymptote and
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 informationWe N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions
We N106 02 Converted-phase Seismic Imaging - Amplitudebalancing -independent Imaging Conditions A.H. Shabelansky* (Massachusetts Institute of Technology), A.E. Malcolm (Memorial University of Newfoundland)
More informationAdvances in radial trace domain coherent noise attenuation
Advances in radial trace domain coherent noise attenuation ABSTRACT David C. Henley* CREWES, Department of Geology and Geophysics University of Calgary, Calgary, AB henley@crewes.org The radial trace transform,
More informationPrestack residual migration in the frequency domain
GEOPHYSICS, VOL. 68, NO. (MARCH APRIL 3); P. 634 64, 8 FIGS. 1.119/1.156733 Prestack residual migration in the frequency domain Paul C. Sava ABSTRACT Prestack Stolt residual migration can be applied to
More informationGeneralized Network Flow Programming
Appendix C Page Generalized Network Flow Programming This chapter adapts the bounded variable primal simplex method to the generalized minimum cost flow problem. Generalized networks are far more useful
More information3D SRME: An Acquisition Independent Approach
Summary 3D SRME: An Acquisition Independent Approach A.Pica, P. Lanfranchi*, Ph. Herrmann CGG, Paris, France Conventional data-driven SRME techniques do not require any a priori knowledge of the subsurface
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 informationTarget-oriented computation of the wave-equation imaging Hessian
Stanford Exploration Project, Report 117, October 23, 2004, pages 63 77 Target-oriented computation of the wave-equation imaging Hessian Alejandro A. Valenciano and Biondo Biondi 1 ABSTRACT A target-oriented
More information(t x) domain, pattern-based multiple separation
Stanford Exploration Project, Report 103, April 27, 2000, pages 63?? (t x) domain, pattern-based multiple separation Robert G. Clapp and Morgan Brown 1 ABSTRACT Pattern-based signal/noise separation is
More informationSimultaneous source separation via robust time variant Radon operators
Simultaneous source separation via robust time variant Radon operators Summary Amr Ibrahim, University of Alberta, Edmonton, Canada amr.ibrahim@ualberta.ca and Mauricio D. Sacchi, University of Alberta,
More informationSparsity-promoting migration accelerated by message passing Xiang LI and Felix J.Herrmann EOS-UBC
Sparsity-promoting migration accelerated by message passing Xiang LI and Feli J.Herrmann EOS-UBC SUMMARY Seismic imaging via linearized inversion requires multiple iterations to the least-squares misfit
More informationEARTH SCIENCES RESEARCH JOURNAL
EARTH SCIENCES RESEARCH JOURNAL Earth Sci. Res. J. Vol. 10, No. 2 (December 2006): 117-129 ATTENUATION OF DIFFRACTED MULTIPLES WITH AN APEX-SHIFTED TANGENT- SQUARED RADON TRANSFORM IN IMAGE SPACE Gabriel
More informationEstimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples
GEOPHYSICS, VOL. 62, NO. 5 (SEPTEMBER-OCTOBER 1997); P. 1596 1611, 18 FIGS. Estimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples D. J. Verschuur* and A. J.
More informationRemoving blended sources interference using Stolt operator Amr Ibrahim, Department of Physics, University of Alberta,
Removing blended sources interference using Stolt operator Amr Ibrahim, Department of Physics, University of Alberta, amr.ibrahim@ualberta.ca Summary We use Stolt migration/de-migration operator to remove
More informationSeismic data reconstruction with Generative Adversarial Networks
Seismic data reconstruction with Generative Adversarial Networks Ali Siahkoohi 1, Rajiv Kumar 1,2 and Felix J. Herrmann 2 1 Seismic Laboratory for Imaging and Modeling (SLIM), The University of British
More informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 5: Refraction Migration Methods Wed, April 13, 2011
GEOPHYS 242: Near Surface Geophysical Imaging Class 5: Refraction Migration Methods Wed, April 13, 2011 Migration versus tomography Refraction traveltime and wavefield migration The theory of interferometry
More informationDesign and specifications for SLIM s software framework
Hello, World! Design and specifications for SLIM s software framework Progress and road map SLIM group at UBC EOS February 2008 About Our working environment Moving towards Python From serial to parallel
More informationStanford Exploration Project, Report 120, May 3, 2005, pages
Stanford Exploration Project, Report 120, May 3, 2005, pages 167 179 166 Stanford Exploration Project, Report 120, May 3, 2005, pages 167 179 Non-linear estimation of vertical delays with a quasi-newton
More informationDownloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at
Recovering the Reflectivity Matrix and Angle-dependent Plane-wave Reflection Coefficients from Imaging of Multiples Alba Ordoñez PGS/UiO*, Walter Söllner PGS, Tilman Klüver PGS and Leiv J. Gelius UiO Summary
More informationAttenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space
Attenuation of diffracted multiples with an apex-shifted tangent-squared radon transform in image space Gabriel Alvarez, Biondo Biondi, and Antoine Guitton 1 ABSTRACT We propose to attenuate diffracted
More informationAnisotropy-preserving 5D interpolation by hybrid Fourier transform
Anisotropy-preserving 5D interpolation by hybrid Fourier transform Juefu Wang and Shaowu Wang, CGG Summary We present an anisotropy-preserving interpolation method based on a hybrid 5D Fourier transform,
More informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationSimply denoise: wavefield reconstruction via jittered undersampling
Simply denoise: wavefield reconstruction via jittered undersampling Gilles Hennenfent 1 and Felix J. Herrmann 1 ABSTRACT In this paper, we present a new discrete undersampling scheme designed to favor
More informationSUMMARY. min. m f (m) s.t. m C 1
Regularizing waveform inversion by projections onto convex sets application to the D Chevron synthetic blind-test dataset Bas Peters *, Zhilong Fang *, Brendan Smithyman #, Felix J. Herrmann * * Seismic
More informationInverse Continuous Wavelet Transform Deconvolution Summary ICWT deconvolution Introduction
Marcilio Castro de Matos*, Sismo Research&Consulting and AASPI/OU, and Kurt J. Marfurt, The University of Oklahoma (OU) Summary Most deconvolution algorithms try to transform the seismic wavelet into spikes
More informationEquivalence of source-receiver migration and shot-profile migration
Stanford Exploration Project, Report 112, November 11, 2002, pages 109 117 Short Note Equivalence of source-receiver migration and shot-profile migration Biondo Biondi 1 INTRODUCTION At first glance, shot
More informationFull-waveform inversion using seislet regularization a
1 Full-waveform inversion using seislet regularization a a Published in Geophysics, 82, no. 5, A43-A49, (2017) Zhiguang Xue 1, Hejun Zhu 2 and Sergey Fomel 1 ABSTRACT Because of inaccurate, incomplete
More informationRecovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples
Recovering the Reflectivity Matrix and Angledependent Plane-wave Reflection Coefficients from Imaging of Multiples A. Ordoñez* (PGS), W.F. Sollner (PGS), T. Klüver (PGS) & L.G. Gelius (UiO) SUMMARY A joint
More informationOn the Scattering Effect of Lateral Discontinuities on AVA Migration
On the Scattering Effect of Lateral Discontinuities on AVA Migration Juefu Wang* and Mauricio D. Sacchi Department of Physics, University of Alberta, 4 Avadh Bhatia Physics Laboratory, Edmonton, AB, T6G
More informationSeismic reflection data interpolation with differential offset and shot continuation a
Seismic reflection data interpolation with differential offset and shot continuation a a Published in Geophysics, 68, 733-744 (2003) Sergey Fomel ABSTRACT I propose a finite-difference offset continuation
More informationAngle Gathers for Gaussian Beam Depth Migration
Angle Gathers for Gaussian Beam Depth Migration Samuel Gray* Veritas DGC Inc, Calgary, Alberta, Canada Sam Gray@veritasdgc.com Abstract Summary Migrated common-image-gathers (CIG s) are of central importance
More informationCurvelet Transform with Adaptive Tiling
Curvelet Transform with Adaptive Tiling Hasan Al-Marzouqi and Ghassan AlRegib School of Electrical and Computer Engineering Georgia Institute of Technology, Atlanta, GA, 30332-0250 {almarzouqi, alregib}@gatech.edu
More information3D pyramid interpolation
3D pyramid interpolation Xukai Shen ABSTRACT Seismic data geometries are not always as nice and regular as we want due to various acquisition constraints. In such cases, data interpolation becomes necessary.
More informationCHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT
CHAPTER 9 INPAINTING USING SPARSE REPRESENTATION AND INVERSE DCT 9.1 Introduction In the previous chapters the inpainting was considered as an iterative algorithm. PDE based method uses iterations to converge
More information