- Home
- Conferences
- Conference Proceedings
- Conferences
73rd EAGE Conference and Exhibition incorporating SPE EUROPEC 2011
- Conference date: 23 May 2011 - 27 May 2011
- Location: Vienna, Austria
- ISBN: 978-90-73834-12-5
- Published: 23 May 2011
1 - 50 of 799 results
-
-
Joint Inversion of Surface-wave Dispersion and P-wave Refraction Data for Laterally Varying Layered Models
Authors L. V. Socco, D. Boiero and C. CalzoniWe present here a joint inversion method to build P- and S-wave velocity models from surface-wave and P-wave refraction data, specifically designed to deal with laterally varying layered environments which can present strong velocity contrasts with depth. In this case, a smooth minimum-structure inversion produces smooth models even for geological models, which are overall layered. The proposed algorithm is also able to incorporate a-priori information available over the site and any physical law to link model parameters. This method presents advantages with respect to individual surface wave analysis and refraction tomography since it imposes internal consistency for all the model parameters, reducing the required a-priori assumptions and the ill-ness of the two methods. We describe the algorithm and we show its application to synthetic and field datasets.
-
-
-
Shallow Seismic on Top of the European Roof – A Vibroseismic Experiment on Firn and Ice at Colle Gnifetti
Authors U. Polom, C. Hofstede, A. Diez and O. EisenIn the summer season 2010, a small shallow reflection seismic experiment was carried out on the firn and ice cover of the Colle Gnifetti, Monte Rosa group, Swiss/Italian Alps. At this site, the physical properties of ice are similar to polar regions, wherefore it is widely used for testing. The challenging experiment in 4500 m asl was designed to explore the scope of shallow vibroseis for seismic targets within and below the firn and ice mass. The small ElVIS vibrator system was used to generate shear waves and compression waves for SH-wave and P-wave receiver setups of two profiles. The resulting sections clearly show the boundary from ice to rock and deeper structures. The deepest features are estimated to 150 m for the SH-wave and 220 m for the P-waves. Reflections could be detected within the ice overburden, which are preliminarily interpreted as change of the crystal structure in the ice column. Furthermore, elastic parameters could be derived from seismic velocities, due to the clear basement reflections. The results of this unique experiment enable new insights in the internal structure of ice formations, and open a promising new investigation method for sub-ice structures and properties, such as basal sediments.
-
-
-
Determination of Spatial Distribution of Topsoil Texture and Available Potassium Contents Using Gamma Ray Spectrometry
By A. M. FakeyeThere are often difficulties in adequately mapping spatial distribution of soil properties at small scale due to soil variations. Gamma-ray spectrometry is a new technique that can potentially address this while improving the mapping of soil texture and available potassium (avK) and minimising huge soil survey cost and maximizing economic returns. To achieve this, secondary spatial information is incorporated into the study. Using portable Exploranium GR 320 gamma ray spectrometer, concentrations of 40K, 238U and 232Th in sandstone washed materials were measured in three plots while representative soil samples were analysed. The result showed homogeneity in on-ground γ-rays concentrations derived from γK based on the geo-lithological nature of soils. The recorded γK values and total gamma ray intensities were related to the soil texture and avK properties. Strong linear relationships were identified between γK and avK (r2 = 0.78, P < 0.001); and soil texture (r2 = 0.48 to 0.67). The success of γ-ray spectrometry in predicting avK and soil texture classes relied on this strong relationship that transformed γK map into a 0.2m resolution soil property map by regression analysis, associated with 78% variance in avK and 51% clay.
-
-
-
Algorithm for the Restoration of Clipped GPR Amplitudes
Authors A. Gulati and R. J. FergusonIt is common in Ground Penetrating Radar (GPR) imagery to have missing or corrupted traces. This can be either due to obstacles, noise, technical problems or economic considerations. Antenna-ground coupling is another reason for clipped amplitudes in GPR data. Most commercially available software use the famous "rubber band interpolation", which uses the spline polynomial to undo the clippings. This method is a simple polynomial based interpolation which performs declipping without considering any prior knowledge about the signal. In this paper, a modified Projection on convex set (POCS) method hybird with Non uniform fast Fourier transform (NFFT) is adopted for reconstruction of clipped amplitudes. Restoration of bandlimited GPR data which has undergone amplitude clipping is studied. This algorithm is tested on synthetic and then real clipped GPR data. To study the effectiveness of the technique, results obtained are compared with industry standard spline interpolation.
-
-
-
Using GPR to Detect Bauxite Horizons in Laterite Deposits of Amazon Basin, Brazil
Authors M. A. Barsottelli-Botelho and J. G. LuizThe GPR was applied on areas with laterite deposits and bauxite occurrences to see its capability of discriminate the different weathered horizons and detect the presence of bauxite horizon and also measure its width. The GPR was tested on two bauxite occurrences in Pará State, Brazil: Paragominas Region and Trombetas River. The good quality of the radargrams was expected since the electric resistivity of the soil, evaluated by Vertical Electrical Sounding which reveals a resistivity varying from 450 to 3000 Ohm.m. The GPR data acquisition with constant-offset results in clear GPR sections where it was possible to separate four zones or horizons as well as its thicknesses. Each zone shows a different reflection pattern or radarfacies, which helps to differentiate them. The first horizon corresponds to the Belterra Clay, which is the most superficial layer, (AT); below, there is the laterite sequence composed by more three horizons, where the second one is a diffractions zone, which corresponds to a clay with concretionary laterites (ND); the third zone shows a homogeneous pattern and it corresponds to a massive ferruginous horizon (LT); an fourth zone shows discontinuous and small reflectors and/or chaotic pattern which corresponds to the main bauxite horizon (BX).
-
-
-
Seismic Radar Using Time Reversal Mirrors
By S. M. HanafyA new application of Time Reversal Mirrors (TRM) is introduced in this work. It is similar to radar surveillance, detecting and tracking moving object(s) in a remote area using seismic waves; this includes the exact location of the moving object and its speed. Recordings of high signal-to-noise ratio reference Green's functions are required, where they are used as calibration records or moving sensors. The proposed method is tested with a field experiment.
-
-
-
SEM Wave Propagation in Complex Media with Tetrahedral to Hexahedral Mesh
Authors M. Charara, A. Vershinin, D. Sabitov and G. Pekarated the ability to handle high-resolution simulations of seismic wave propagation in 3D complex domains. However, the exponential accuracy of the method and the reduction of the computational effort rely on the use of conformal hexahedral meshes. Generating an all-hexahedral mesh based upon the available meshing can be difficult. We propose to use a 3D unstructured tetrahedral mesh generator and then split the resulting mesh into hexahedral elements. This approach allows using the same modeling SEM program without any modifications, and it remains faster than dealing with tetrahedral spectral element due to the tensorisation property for spatial derivatives. With this approach, the SEM method doesn’t have restrictions to solving complex problems in sonic or seismic data with respect to methods using tetrahedral meshes. Two examples dealing, respectively, with sonic modeling and the other with surface seismic modeling illustrate the feasibility of the proposed method.
-
-
-
Solving the 3D Acoustic Wave Equation with Higher-order Mass-lumped Tetrahedral Finite Elements
Authors E. Zhebel, S. Minisini, A. Kononov and W. A. MulderPresent-day computers allow for realistic 3D simulations of seismic wave propagation, as well as migration and inversion of seismic data with numerical solutions of the full wave equation. The finite-difference method is popular because of its simplicity but suffers from accuracy degradation for complex models with sharp interfaces between large impedance contrasts and for models with rough topography. A tetrahedral mesh offers more flexibility and maintains its accuracy if element boundaries are aligned with sharp interfaces. Higher-order finite elements with mass lumping provide a fully explicit time-stepping scheme. We have implemented elements of degree one, two, and three for the 3D acoustic wave equation. Numerical tests confirm the accuracy of the mass-lumped elements. There are two different third-degree elements that have almost the same accuracy, but one has a more favourable stability limit than the other. Convergence analysis shows that the higher the order of the element, the better the computational performance is. A low-storage implementation with OpenMP shows good scaling on 4-and 8-node platforms.
-
-
-
3D Wave Propagation Modeling in Presence of Topography
Authors I. Tarrass, F. Aye, L. Giraud and P. ThoreWe studied and developed an efficient method to handle the wave simulation in presence of topography based on curvilinear finite difference. Foremost, we derived the modified equations on the continuous problem. We used afterwards optimized stencils and optimized selective filters adapted from aeroacoustics.
-
-
-
Time-varying Boundary Conditions in Simulation of Seismic Wave Propagation
Authors R. P. Fletcher and J. O. A. RobertssonWe propose two new boundary conditions to regulate coherent reflections from the model boundaries in numerical solutions of wave equations. Both boundary conditions have the common feature that the boundary condition is varied with respect to time. The first boundary condition expands or contracts the computational model during a modeling simulation. The effect is to cause a Doppler shift in the reflected wavefield that can be used to shift energy outside a frequency band of interest. Additionally, when the computational domain is expanding, the range of possible incidence angles on the boundary is restricted. This can be used to increase the effectiveness of many existing absorbing boundary conditions that are more effective for incidence angles close to normal. The second boundary condition is an extension of random boundaries. By carefully changing the realization of a random boundary over time, a more diffusive wavefield can be simulated. We show results with 2D numerical simulations of the scalar wave equation for both these boundary conditions. While the first boundary condition has application to modeling, both of these boundary conditions have potential application within algorithms that rely upon modeling kernels, such as reverse-time migration and full-waveform inversion.
-
-
-
Elastic Corrections to Acoustic Finite-difference Simulations – Plane-wave Analysis and Examples
Authors J. W. D. Hobro, C. H. Chapman and J. O. A. RobertssonAcoustic finite-difference modelling is playing an increasingly important role in seismic imaging (e.g. in reverse time migration) but the additional cost of elastic finite-difference modelling restricts its use in commercial imaging technology. The cost of full elastic finite-difference modelling can exceed the cost of acoustic modelling in the same velocity model by two orders of magnitude or more. A technique is discussed that corrects an acoustic finite-difference simulation for elastic effects. It is based on estimating the error incurred when using the acoustic wave equation as an approximate solution to the elastic wave equation. The errors are used to generate an effective source field for an additional acoustic simulation that calculates a correction to the wavefield produced in the original acoustic simulation. The cost of this approach is roughly twice that of an acoustic simulation and therefore much less than that of a full elastic simulation. Plane wave theory is used to examine the properties of the method and results from finite-difference simulations are presented to illustrate its accuracy.
-
-
-
Wavefield Extrapolation in the Pseudo-depth Domain
Authors X. Ma and T. AlkhalifahWe introduce the pseudo-depth as an extension to vertical time representation of acoustic wavefield. By formulating the mapping relations between conventional depth and pseudo-depth, we derive an acoustic wave equation for isotropic media in the pseudo-depth domain. Because of a uniform vertical wavelength in the pseudo-depth domain, wavefield extrapolation requires less number of vertical samples to avoid aliasing. Reduction of samples implies a reduction of run time for finite-difference time extrapolation. To highlight its application potentials, we performed an exploding reflector modeling and zero-offset migration using pseudo-depth wavefield extrapolation. The results agree well with conventional vertical-depth modeling and migration. The pseudo-depth wavefield extrapolation also showed its potential in eliminating multiple artifacts resulting from using a discontinuous velocity model in migration.
-
-
-
Seeing Coal Top Ahead of the Drill-bit Through Seismic-while-drilling – Results from Numerical Modeling
Authors B. Zhou, I. M. Mason and S. A. GreenhalghBlast damage to coal seams due to incorrect blast stand-off distances costs the equivalent of one unnecessary open cut mine for every 10 existing Australian mines. The current approach for mapping coal seam tops is through drilling and pierce-point logging. To provide appropriate depth control with ±0.2 m precision for blast hole drilling, a 50m x 50m grid is required in advance of overburden removal. Pierce-point mapping can be inaccurate and expensive as the seam is not flat. To date, there are no effective and economic techniques to map coal seam structures in the open cut environment where coal seams are shallow. In this paper, we propose two Seismic-While-Drilling (SWD)-based approaches to see the top of the coal seam ahead of the drill-bit: 1) reverse VSP recording, in which the geophones are planted on the surface near the drill rig; 2) in-seam seismic recording, in which the geophones are locked in the coal seam in a borehole. Numerical modeling demonstrates the feasibility of these two approaches for predicting the approach of coal seams during blasthole drilling. If these methods can be reliably realized in practice, significant coal resources can be saved and the mines can become more productive and profitable.
-
-
-
Benchmark of Amplitude Preserved Imaging Algorithms Using Large Scale Synthetic Data
Authors H. Allouche, E. Crouzy, P. Thore, M. Brucher, G. Bailly and C. GodefroyIn the work we present, we analyze the amplitude preservation issue of imaging algorithms. We use for that purpose synthetic data that are generated by Full-wave seismic modeling. A particular emphasize is placed on the model building step which is very challenging because there are few robust tools that are adapted to build large and complex models. A 3D full elastic Finite Difference scheme is used to generate the seismic data. These generated data are primarily used to test the amplitude preservation of imaging algorithms and for other geophysical applications as well. The main results will be presented.
-
-
-
Analysis of Patchy CO2 Saturation from Time-lapse Sonic Logs Using Rock Physics Modelling
Authors E. Caspari, T. M. Müller and B. GurevichWe compare time-lapse sonic and neutron porosity logs of the Nagaoka CO2 sequestration experiment against the uniform and patchy saturation models, which represent two end-members of the P-wave velocity and CO2 saturation relationship. Most of the data points fall between the two limits, suggesting that the relationship is somewhere between uniform and patchy saturation. The behaviour between these limits can be explained by the mechanism of wave-induced-fluid on mesoscopic fluid heterogeneities (porescale << mesoscale << P-wavelength), which causes wave attenuation and velocity dispersion. We model these fluid effects using the 1D and 3D continuous random media model (CRM). The log data approximately follow the predictions of the CRM models for fluid patch sizes of 1 to 5 mm. This heterogeneity scale is much larger than the porescale features of a reservoir thin section, indicating that the mechanism of wave-induced fluid flow on the mesoscale can occur in sonic log data and therefore controls the velocity-saturation relation.
-
-
-
Borehole Transient Electromagnetics for Monitoring CO2 Sequestration in Saline Aquifers
Authors K. Spitzer, J. Börner, M. Afanasjew, R. U. Börner, O. Ernst and M. EiermannWe present a technique for monitoring CO2 migration in an underground gas storage by using a surface-to-borehole transient electromagnetic method (BTEM). As the targeted reservoir consists of a saline aquifer, the main idea is to detect the changes in the electric conductivity that go along with the resistive CO2 displacing the conductive pore fluid thus decreasing the rock conductivity significantly. There are several major problems to tackle: 1. the three-dimensional simulation of the transient EM fields, 2. the inversion of BTEM data in a time-lapse mode, 3. the description of the petrophysical processes associated with the CO2 migration, and 4. the experimental design to increase the sensitivity at the target depth as well as its practical implementation. Here, we mainly focus on the experimental setup employing BTEM and its interpretation technique using efficient vector finite element simulation procedures combined with Krylov subspace methods. This work is supplemented by extensive laboratory studies. Due to the wide range of practical and theoretical problems, our group is strongly interdisciplinary and consists of numerical mathematicians as well as laboratory, field, and theoretical geophysicists.
-
-
-
Electromagnetic Characterization of CO2 Sequestration Sites – Feasibility Studies and First Field Results from Ketzin
Authors R. Streich, M. Becken and O. RitterThe electrical resistivity of sedimentary rocks depends strongly on pore fluid content and pore connectivity, and bulk resistivity increases where carbon dioxide displaces conductive saline pore fluid. Therefore, controlled-source electromagnetic (CSEM) methods are potentially a key geophysical technique for exploring and monitoring carbon sequestration sites. We investigate the feasibility of CSEM characterization of the CO2 sequestration pilot site at Ketzin, Germany. Simulation studies using newly developed, accurate 1D and 3D modelling tools indicate that surface-based CSEM surveys may well be able to image the background resistivity structure, but will probably not detect the small CO2 reservoir. Borehole-to-surface configurations exhibit much higher sensitivity to the reservoir and may thus be particularly well-suited for monitoring applications. A recently conducted large-scale CSEM field survey has demonstrated the feasibility of CSEM measurements at the site, and provided data that contain subsurface responses down to the target reservoir depth.
-
-
-
Data-driven Hydraulic Barriers Delineation and Migration Pathway Detection within the CO2 Geological Storage of Sleipner
More LessThe deep saline aquifers constitute the strongest storage capacity for CO2 storage among the portfolio of possible geological opportunities. The upward and lateral migration of CO2 after being injected into such formations depends (i) on permeability and porosity of the storage formation, (ii) on the potential hydraulic barriers. At Sleipner, the CO2 is injected inside the Utsira formation, where hydraulic barriers correspond to thin intra-shale layers identified by logs without information on their lateral continuity. Furthermore, conventional interpretations of 3D seismic data cannot provide better constraints because the intra-shale layer thickness is smaller than seismic resolution. The proposed method consists in refining intra-shale layer delineation using 4D seismic data in the area invaded by CO2, but also in the Utsira formation part covered by time-lapse seismics. The shale bottom delineation is consistent with 3D-detailed geometry of the CO2 plume and with independent observations of shale at a well. To our knowledge such a detailed 3D description of the Utsira hydraulic barriers has not been published before. Hence, the intra-shale layer discontinuities, allowing the upward migration of CO2, can be identified and could be integrated, in further steps, in the geological static model, unavoidable stage for an efficient dynamic flow simulation.
-
-
-
Traveltime and Attenuation Tomography of CO2 Plume at Sleipner
Authors G. Rossi, R. A. Chadwick and G. A. WilliamsWithin the European project CO2ReMoVe traveltime and attenuation tomography have been applied on one of the high resolution seismic line acquired in 2006 on the Sleipner CO2 geological storage site. The aim is to gain an accurate model in depth, of both velocities and attenuation, to constrain better the quantification studies. The results obtained from the tomographyc inversion of traveltimes and spectral-centroid frequency shift, respectively, confirm the cumulation of CO2 in the uppermost layer of the reservoir, whereas within the plume, lateral variations of both velocity and quality factor values are observed. The presence of fine shale lenses, and their impact on the CO2 distribution within the plume may be the reason for these lateral variations. Through the comparison of the velocity and Q tomographic values with the theoretical curves from petrophysics studies, it will be possible to interpretate them in terms of effective variations in CO2 saturation and heterogeneities within the plume.
-
-
-
Wettability Determination of HVBB Coal – Water System with Injection of Model Flue Gas and CO2
Authors N. Shojaikaveh, E. S. J. Rudolph, J. Bruining and S. N. AshrafizadehGeological sequestration of carbon dioxide (CO2) in coal is one possible method to sequester CO2. CO2 or flue-gas injection into coal layers provides the dual benefits of enhancing coal-bed methane production (ECBM) and sequestering CO2. The success of this combined process depends strongly on the wetting behavior of the coal as a function of coal rank and roughness, system pressure and temperature, and also the composition of the gas. The wetting behavior can be evaluated from the contact angle of a gas bubble on the coal surface. In this study, we determine contact angles of a model flue gas (80/20 N2/ CO2 mixture) and pure CO2on a wet Warndt Luisenthal coal (WL), representing a high volatile bituminous coal (hvbB coal). Experiments were performed in the pressure range between 0 and 150 bar using a modified pendant drop cell at a constant temperature of 318 K. In general, contact angles of flue gas on WL coal are smaller than of CO2 on WL coal. The contact angle of flue gas against water on WL coal changes from water-wet to intermediate-wet with increasing pressure. With CO2 injection the wettability of WL coal becomes gas-wet for pressures above 90 bar.
-
-
-
Data Reconstruction via Sparse Bi-focal Transformation
Authors H. Kutscha and D. J. VerschuurIn most reconstruction algorithms information about the subsurface can not be utilised, even if such is available. Focal transformation is a way to effectively incorporate prior knowledge of the subsurface in seismic data reconstruction. The basis functions of this transformation are the focal operators. They can be understood as one-way propagation operators from certain effective depth levels in a prior velocity model. A sparseness constraint is used to penalize aliasing noise. By using several depth levels simultaneously the data can be described with less parameters in the transform domain. This results in a better signal to noise separation and, therefore, improved reconstruction. In addition, we introduce a smart choice of focal operators, based on data-driven operator updating where we follow closely the reflection events in the data. This allows an even stronger compression of the data in the focal domain.
-
-
-
Parametric Matching Pursuit Methods to Reconstruct Seismic Data Acquired with Multichannel Sampling
Authors A. Ozbek, M. Vassallo, K. Ozdemir, D. J. van Manen, K. Eggenberger and J. O. A. RobertssonIn marine seismic data acquisition, multichannel sampling can occur in various scenarios, including over/under streamers, multicomponent streamer, and multicomponent ocean-bottom node. In these various multichannel sampling scenarios, each channel outputs samples of a filtered version of a seismic wavefield. These filters are different for each channel; the data are often coarsely and irregularly sampled. The objective is to reconstruct the unknown wavefield from the multichannel samples. In this paper, we discuss how parametric matching pursuit methods can be extended to solve such reconstruction problems. We give two examples based on data that a multicomponent streamer could record and show the form the reconstruction problem takes when two- and three-component data are available.
-
-
-
Matching Pursuit Methods Applied to Multicomponent Marine Seismic Acquisition – The Issue of Crossline Aliasing
Authors M. Vassallo, A. Ozbek, K. Eggenberger, K. Ozdemir, D. J. van Manen and J. O. A. RobertssonIn this work we analyze the theoretical aspects of spatial aliasing in the crossline direction in marine seismic acquisition. We also explain the benefits of the additional measurements acquired by a multicomponent towed streamer, able to measure the three components of the particle velocity vector in addition to the pressure wavefield. We propose matching pursuit based techniques to reconstruct a 3D full bandwidth seismic wavefield on a fine receiver grid. The techniques that we describe process multicomponent seismic data; they calculate the desired 3D wavefield with satisfactory quality despite the severe aliasing that affects each of the individual input measurements in the crossline direction.
-
-
-
Convergence Improvement and Noise Attenuation Considerations for POCS Reconstruction
Authors J. J. Gao, M. D. Sacchi and X. H. ChenIn recent years, a variety of algorithms for multidimensional seismic survey reconstruction have been proposed. Projection onto Convex Sets (POCS) seismic trace reconstruction appears to be an affective, uncomplicated and robust algorithm for the recovery of unrecorded seismic traces. Its slow convergence, however, could jeopardize its computational appeal. For this reason, we investigate the performance of POCS in terms of different threshold schedules. We show that a data driven threshold schedule leads to an efficient implementation of POCS where high-quality solutions can be obtained in a modest number of iterations (>30-40). In addition, we present a variant of POCS that permits concurrent reconstruction and denoising of seismic volumes. The latter is achieved by defining a weighted reinsertion trace strategy that alleviates the influence of the noisy traces in the final reconstruction of the seismic volume. We illustrate our findings with synthetic and real data examples.
-
-
-
5D Seismic Volume Reconstruction Using HOSVD
Authors N. Kreimer and M. D. SacchiWe propose a tensor completion method to reconstruct sparse prestack data. An algorithm based on the High-Order SVD (HOSVD) is used to estimate a low rank decomposition of the seismic volume. The latter allows us to use rank-reduction of tensors to estimate missing traces and increase the SNR of the seismic volume. A real data set from the Western Canadian Sedimentary Basin is used to test the proposed algorithm.
-
-
-
Minimizing Ellipsoidal Norm Seismic Data Interpolation with Radon Spectrum Constraints
Authors X. Wang, H. Z. Wang, D. H. Zhou and L. X. TianA data adaptive interpolation is designed under the Radon spectrum constraints for spatially aliased seismic data. The inverse process is formed through constructing an ellipsoidal norm, which can make full use of the information about the Fourier spectrum of the data. The ellipsoidal norm can enhance the relative distribution of the Fourier spectrum and is computed with the help of Radon spectrum. Radon spectrum can helpfully concentrate the energy distribution of the Fourier spectrum. Correspondingly, the aliasing components in the Fourier spectrum can be relatively attenuated. The solution of the interpolation inverse problem is sparser and the inverse process becomes more stable. The numerical examples show the effectiveness of the method. The method is applicable to 2D and 3D seismic survey data recorded sparsely in a horizontal plane. The most common pre-stack applications of the algorithm are common-mid-point and half-offset domain interpolation.
-
-
-
Fast Generalized Fourier Interpolation of Nonstationary Seismic Records
Authors M. Naghizadeh and K. A. InnanenWe propose a fast and efficient method for the interpolation of nonstationary seismic data. The method uses the fast generalized Fourier transform FGFT to identify the space-wavenumber evolution of nonstationary spatial signals at each temporal frequency. The nonredundant nature of FGFT renders a big computational advantage to this interpolation method. A least-squares fitting scheme is used next to retrieve the optimal FGFT coefficients representative of the ideal interpolated data. For randomly sampled data on a regular grid, we seek a sparse representation of FGFT coefficients to retrieve the missing samples. A synthetic seismic data example was used to examine the performance of the method.
-
-
-
How to Obtain a Seamless Dataset for a Pre-stack Multi-survey Merge Using 5D Data Reconstruction?
Authors J. L. Rivault and A. MotagallyPre-stack merging multi-survey datasets is often challenging in terms of processing effort as well as time and human resources. It is problematic not only due to issues arising from survey matching (due to differences in sources and receivers characteristics) but also from variations in sampling, which may relate to acquisition azimuth, bin size, shot/receiver line spacing. For such datasets, it is important to regularise the sampling of the data in order to harmonise fold and minimise migration artefacts. In this paper we demonstrate the use of 5D data reconstruction for the benefit of multi-survey merging by mapping all the different surveys to a single acquisition design. The result is a dataset with a single acquisition configuration which has a constant fold and a significantly reduced level of migration noise. It also allows the preservation of azimuthal information by migrating in the COV domain. Its shows that 5D data reconstruction fully preserves the original seismic aspect allowing to easily migrate pre-stack multi-survey vintages in a single pass, which reduces the work of the geophysicist.
-
-
-
Determination of Shear-wave Velocity in Seabed from Dispersion Data of Love Waves
Authors H. Dong, G. Ke and K. DuffautEstimation of the shear (S-) wave velocity or shear strength in the seabed is of importance for geotechnical applications, characterization of subsurface, seismic exploration studies and geohazard investigation. S-wave velocities in the sediment can be estimated by invert the dispersion data of seismic interface waves, which will be considered in this study. Numerical study shows that Love wave sensitivity and inversion stability in linearized inversion is higher than that of the Rayleigh wave and the inversion of Love wave dispersion curve provides a more accurate results. However, in conventional underwater seismic experiments only Scholte wave can be generated. In a testing experiment S-wave sources were used and both Scholte wave and Love wave were generated and recorded. In this study the dispersion of Love waves are analyzed using three time-frequency analysis methods. The S-wave velocities are estimated by inverting the dispersion data. The inversion results show the maximum penetration depth of the estimates is around 25 m in the seafloor, where the estimated S-wave velocity is around 220m/s. This approach provides an alternate and supplementary means to estimate the seismic velocity structure that is valuable in both seafloor geotechnical engineering and the static correction in seismic exploration using reflection methods.
-
-
-
PZ Summation of 3D WAZ OBS Receiver Gathers
Authors P. Hugonnet, J. L. Boelle, P. Herrmann, F. Prat and A. LaframPZ summation is a typical processing flow for 2-Component Ocean Bottom Seismic (OBS) data, which takes advantage of the simultaneous recording of the pressure and vertical velocity wavefields to separate the up-going and down-going wavefields and to attenuate the water bottom receiver-side peg-legs. However, practical implementations are often either monochannel/1D or designed for 2D receiver gathers, while modern acquisition geometries often lead to 3D, wide-azimuth (WAZ), densely sampled receiver gathers. Such 3D gathers can be processed in the tau-px-py domain with a monochannel algorithm. In this paper we extend the method to directly process 3D gathers, without the need for tau-px-py transforms. Some current implementations also attenuate the source-side multiples, but they rely on a first order approximation for the amplitudes and on the 1D geology (horizontally layered medium) assumption. By introducing a simultaneous 3D predictive deconvolution operator in the peg-leg attenuation formulation, we remove the first order approximation and we improve the attenuation in the case of deviations from the 1D assumption.
-
-
-
Enhanced Wavefield Separation of OBS Data
Authors E. Zabihi Naeini, L. Baboulaz and S. GrionIn ocean-bottom seismic (OBS) data processing, wavefield separation results are sometimes affected by high levels of noise on the vertical component Z, while the pressure component P is in general of good quality. Nonetheless, Z is needed to achieve complete pre-stack wavefield separation and also to drive processes such as mirror imaging and up-down deconvolution. To address the problem of noise on Z affecting wavefield separation results, we propose a new method which first estimates a multiple model from the downgoing wavefield in a least-squares fashion. Next, this multiple model is used as a measure of seismic signal coherency to calculate an enhanced upgoing wavefield with minimal noise degradation. We show the benefits of this new method on a real data example from the Caspian Sea.
-
-
-
Velocity-based Wavelet Corrections for Domain Transformation
By J. E. GaiserIt is well established that similar vertical wavelength ranges must be preserved in multicomponent data and that wavelengths of P- and S-waves must match in order to sample reflectivity in an equivalent manner. Conversion of a wavefield to another time or depth domain is described by transformation functions that depend on average VP/VS ratios and velocity. Although these functions align corresponding stratigraphic events of different wavefields, they distort the seismic wavelet because global average velocity properties are independent of local interval properties that define wavelength. In this study we develop a theory of velocity-based wavelet corrections for domain transformations, which are expressed as functions of interval and average VP/VS and velocity, to match wavelength of multicomponent wavefields. We examine the effects for both land and marine data examples and find that land surveys are affected more than marine, and may require spectral broadening of the wavelet. Data from the Marcellus shale in northeast Pennsylvania, USA, shows significant bandwidth improvements for C-waves when wavelet corrections based on velocity match their wavelengths with P-waves. Application of these wavelet corrections should benefit registration fidelity, joint AVO/A (offset and azimuth) inversions and attribute analyses.
-
-
-
Polarity Reversal Correction for Multicomponent Joint Elastic Reverse Time Migration
Authors Q. Z. Du, Y. T. Zhu, G. Fang and X. F. GongMulticomponent joint elastic reverse time migration is one of important parts in multicomponent seismic data processing. We analyse the implementation of elastic reverse time migration and the polarity reversal in PS and SP images. We present a method using Poynting vector to correct the polarity reversal in elastic reverse time migration. The polarity reversal correction using Poynting vector can be done with the imaging conditions at every time step of the elastic RTM. And the whole correction procedure can be carried out in common shot domain rather than angle domain, which is low cost and convenient. Examples show this method can accurately locate the normal incidence point and correct the polarity reversal for both simple and complex models.
-
-
-
Using 2D–3C Seismic Data to Discriminate the Sandstone Reservoir in Guang’an Gas Field
More LessThis paper describes the interpretation of 2D-3C seismic survey acquired in the Guang’an Gas Field of Sichuan Basin, which is located at the southwest part of China. Estimating reservoir parameters from geophysical data is vitally important in hydrocarbon exploration and production. We present a joint-inversion algorithm to estimate reservoir parameters directly, using both PP and PS seismic amplitude variation with offset (AVO) data. Reservoir parameters are linked to geophysical parameters through a rock-properties model. The joint inversion method for PP and PS seismic data solves approximated linear expressions of PP and PS reflection coefficients simultaneously using simulated annealing algorithm. By combining the two independent measurements (PP and PS seismic data) together, we stabilize the system of equations for PP and PS seismic data separately, leading to more robust parameter estimation.
-
-
-
High-effort 3C 3D Seismic Survey in Northern Piceance Basin, Colorado, USA
Authors K. T. Lewallen, T. Johns and J. HeftiThe Upper Cretaceous Mesaverde Group in the Piceance Basin, northwest Colorado, USA, contains a large gas accumulation in a tight sand reservoir. Previous work using well log, borehole and surface seismic data demonstrated potential to utilize converted-wave images and attributes recorded on multicomponent data for fracture identification and lithology discrimination. A high-effort multicomponent 3D surface seismic survey was acquired that pushed the limits of available recording technologies. Preliminary processing results show PP images that overcome significant field noise issues and PS attributes consistent with previous assessments. The measured reservoir HTI attributes after layer stripping through the overburden indicate a principal axis direction of 130º and PS1 (fast) and PS2 (slow) traveltime delays between 10 and 40 ms consistent with pre-survey predictions. Processing and research continues to improve PP/PS images and attributes to support business decisions as well as to investigate the effects of fold, offset, azimuth and spatial sampling on data quality for future surveys in tight gas and other unconventional resource plays.
-
-
-
3D-3C Geomechanical Study of In-situ Bitumen Recovery, Alberta, Canada
By K. Wikel3D converted wave data have been used in the past as an indicator of fractures and differential stress, although the emphasis has historically been on fractures. Recently, industry has been analyzing stress directions and stress changes based on the direction of fast converted mode (PS1) and the time lag between the fast and slow converted modes (PS1/PS2) in the near surface. In the compliant near surface of NE Alberta the data provided from converted wave seismic data rotated into fast and slow shear directions allows us to evaluate cap rock integrity and areas where the cap rock may deviate from the regional faulting regime. This information is extremely useful in discerning the cap rocks ability to maintain integrity during shallow high pressure bitumen and heavy oil recovery. In addition to this, the changes in PS1 direction and PS1/PS2 time lag with time can assist us in monitoring reservoir changes and recovery of bitumen. This paper is an overview of a case study from Petrobanks pilot THAI facility near Conklin, Alberta, and the results of a converted wave cap rock integrity study.
-
-
-
Preserved-traveltime Smoothing
Authors V. Vinje, A. Stovas and D. ReynaudMost ray-based migration- and tomography methods require some degree of smoothness in the depth velocity models. A drawback of the conventional smoothing is errors in the traveltimes at the discontinuities of the velocity models. These errors are offset-dependent and they cause errors in both depth and RMO. Here we propose a new method, Preserved-Traveltime Smoothing (PTS) that preserves traveltimes (and hence depths) at these discontinuities. This is done by smoothing of composite anisotropic velocity parameters using a specific convolution filter.
-
-
-
High Resolution Velocity Model for Imaging Complex Structures
More LessWith the improvement of automatic CIG picking tool, the accuracy of tomography is mainly limited by the model representation (e.g. model grid size) and the regularization of the inversion solver. In practice, a dense model grid with the grid size less than 20 meters and a 100 meters’ smoothing are applied to achieve high resolution velocity updates in some well-conditioned areas (e.g. shallow sediment areas). The small grid size and small model smoothing may cause inversion artefacts in poorly conditioned areas (deep part) and increase the computational cost dramatically. In this abstract, we introduce a non uniform grid, which has a dense sampling in shallow and a sparse sampling in deep to reduce the computing cost. Also we propose a scheme for smoothing along geological structures to stabilize the inversion. We demonstrate with a production example in the Gulf of Mexico that our high resolution tomography is able to improve the subsurface images.
-
-
-
RTM and Kirchhoff Angle Domain Common-image Gathers for Migration Velocity Analysis
Authors J. P. Montel and G. LambaréAngle domain common image gathers are recommended in Kirchhoff and reverse time migration for velocity model building in complex area. For these both approaches there is a general agreement that the tomographic ray pairs are fully defined by the reflection and azimuth angle information and the reflection dip and that if the velocity model is correctly updated down to a given horizon, it is not necessary to shoot the tomographic ray pairs upwards through this horizon. We show here through examples and a theoretical analysis that these both statements have to be mitigated when the common image gathers exhibit a significant residual move out on. We also show how to accurately compute the tomographic ray pairs allowing then for an accurate angle domain migration velocity analysis.
-
-
-
Automatic Migration Velocity Analysis Using Reverse Time Migration
Authors W. W. Weibull and B. ArntsenThe objective of this paper is to describe an automatic velocity analysis method based on Reverse Time Migration and Differential Semblance Optimization. The velocity analysis is based on the solution of a nonlinear least squares problem aiming at the focusing of offset domain common image point gathers constructed by Reverse Time Migration. Because the method is based on the solution of the two-way wave equation, it can deal with strong and sharp velocity contrasts both in a stable and accurate manner. It is therefore expected that this method will help improve seismic imaging over complex geological settings. We illustrate the method with a simple synthetic 2D seismic example.
-
-
-
Migration-velocity Analysis Using Image-space Generalized Wavefields
Authors C. Cardoso and B. BiondiIn areas of complex geology, migration-velocity analysis (MVA) should use methods that describe the complexity of wavefield propagation, such as focusing and defocusing, multipathing, and frequency-dependent velocity sensitivity. MVA by wavefield extrapolation has the ability to address these issues because it uses wavefields as carriers of information. However, its high cost and lack of flexibility with respect to model parametrization have prevented its routine industrial use. We overcome those limitations by using new wavefields: the image-space generalized wavefields. These wavefields are synthesized from a prestack image computed with wavefield-extrapolation methods, using the pre-stack exploding-reflector model. Cost of MVA by wavefield extrapolation is decreased because only a small number of image-space generalized wavefields are necessary to accurately describe the kinematics of velocity errors and because these wavefields can be easily used in a target-oriented way. Flexibility is naturally incorporated because modeling these wavefields have selected reflectors as the initial conditions, allowing the use of a horizon-based parametrization of the model space. In a 3D real data example, we show that using wavefields synthesized by the prestack exploding-reflector model greatly improves efficiency of migration-velocity analysis by wavefield-extrapolation, while yielding a final accurate migration-velocity model evidenced by good quality and reliable depth images.
-
-
-
Well-based Interval Velocity Analysis Using Multi-parameter CIG
More LessThis study presents a unique approach of performing interval velocity analysis using multi-parameter CIGs for the 2D case. During the analysis, summation-free CIGs are generated and are analyzed as constant depth slices at every image location. When formation tops information is available, an efficient procedure for updating the velocity can be utilized. A specific pattern which can be used to identify correct and wrong velocity within each of the depth slices is used as the core of the analysis. The ability to analyze the structural parameters of each image point improves the accuracy of the velocity and at the same time provides valuable structural information at the analysis location. The procedure and the potential accuracy of proposed method are demonstrated for the 2D case.
-
-
-
Anisotropic Tomography Using Rock Physics Contraints
Authors Y. Li, D. Nichols, K. Osypov and R. BachrachAnisotropic model building is a well-known challenging problem for its non-linear, ambiguous nature. To reduce the null-space and stabilize the inversion, we propose a new preconditioning scheme in the linearized tomography to include the rock physics prior information. The rock physics information in terms of covariance among P-wave vertical velocity (v0), epsilon and delta is generated by stochastic realizations of a compacting shale model. We design a VSP synthetic survey with the common industry geometry on two different examples to study the effect of applying prior rock physics prior information on proper and improper media, respectively. The results show that by utilizing the proper prior rock physics information, tomography can better resolve the anisotropy parameters, especially in the area where inversion is poorly constrained by the data. However, cautions should be taken when the lithology of the subsurface is unknown. Finally, a posterior uncertainty analysis is performed to evaluate the contribution of the rock physics prior information. The results show that the null-space is greatly reduced by introducing the prior information.
-
-
-
Building Anisotropic Velocity Model Using Rock Physics and Local Tomography – A North Sea Case Study
Authors R. Bachrach and J. MathewsonAnisotropic velocity models are poorly constrained by surface seismic data alone. Rock physics anisotropic modelling of shale compaction and digenesis provides information which can be used to constrain anisotropic velocity model building. We demonstrate this new anisotropic model building methodology which incorporates shale compaction model, check shot data and tomography to build a 3D spatially varying anisotropic velocity model in the North Sea. The initial rock model is built by analysing well log data, temperature gradient and additional geologic information. Tomography with checkshots (local tomography) which takes into account both residual moveout of surface seismic data and vertical traveltime from checkshot data was used to calibrate rock model to predict a similar anisotropy profile. After calibration a 3D anisotropic model is built to produce spatially varying anisotropic parameters consistent with both rock model and checkshot data.
-
-
-
Nonlinear, Full Waveform Modeling in Terms of Forward Multiple Scattering
Authors A. J. Berkhout and D. J. VerschuurThe WRW-model shows how to decompose seismic measurements in terms of responses of individual subsurface grid points that include multiple scattering and, optionally, wave conversion. These Grid Point Responses (GPRs) can be used in full waveform seismic modeling algorithms, allowing a better understanding of the complex (blended) codas of seismic reflections and, above all, utilizing these codas for imaging purposes.
-
-
-
Seismic Diffraction Modeling with the Tip-wave Superposition Method
Authors A. M. Aizenberg, M. A. Ayzenberg and K. D. Klem-MusatovThe transmission-propagation operator theory provides mathematical tools for describing wave propagation in laterally inhomogeneous layered media. The theory is a hybrid of the integral representation (integral-equation) method and the spatial-frequency decomposition (slowness) method. The theory introduces a new wave-motion statement of the forward problem using two operators; the convolutional transmission operator and the feasible propagation operator. In this theory, the exact solution of the initial-boundary value problem is written as the sum of sequentially reflected and transmitted wave events. We further develop the tip-wave superposition method based on the approximations of the two operators in the seismic-frequency range. The feasible propagation operators inside the layers are approximated in the form of layer matrices. The transmission operators at the interfaces are approximated with effective reflection/transmission coefficients, which generalize plane-wave reflection/transmission coefficients to curved interfaces, non-planar wavefronts and finite frequencies. Each element of the layer matrix contains a feasible beam of the tip waves diverging from a small interface element with diffraction at concave parts of the layer boundary and reflecting/transmitting at the interface. We illustrate the potential of the tip-wave superposition method through diffraction modeling of the Green’s function in a layered medium.
-
-
-
3D CRS-attributes Based Diffraction Imaging
Authors S. Dell and D. GajewskiImaging of diffractions is a challenge in seismic processing. Standard seismic processing is tuned to enhance reflections. Separation of diffracted from reflected events is frequently used to achieve an optimized image of diffractions. We present a method to separate diffracted events in the time domain based on the Common-Reflection-Surface approach. The Common-Reflection-Surface method provides a simulated zero offset section as well as kinematic wavefield attributes. Using synthetic data we show, how these attributes can be used to effectively separate reflected and diffracted energy in the 3D case. The objectives of this work are seismic event separation in 3D case and time migration velocity analysis for diffractions implemented in the poststack domain.
-
-
-
Edge and Tip Diffraction Imaging in Three Dimensions
By T. J. MoserThe dimensionality of the diffracted wave is related to the dimensionality of the diffractor: a surface reflected ray is uniquely determined, edge diffracted rays strike out a one-parameter pencil, tip diffracted rays a two-parameter pencil of rays. These topological notions are important to analyse the effect of diffraction imaging by specularity suppression. In diffraction imaging reflectivity is fully suppressed, edge diffractivity partly suppressed, while tip diffractivity stands out unsuppressed. These observations have a significant impact on the structural interpretation of three-dimensional diffraction images. Edges constitute the skeleton of a structural model and tips indicate the joints of the skeleton.
-
-
-
Diffraction Imaging Based on the Diffraction Operator
Authors J. J. S. de Figueiredo, F. Oliveira, E. Esmi, L. Freitas, A. Novais and J. SchleicherSeismic diffractions carry detailed structure information. Some of these structures are well-known as common hydrocarbon traps, such as faults, pinch-outs, unconformities and structures the size of which is in the order of the wavelength. Therefore, the development of computational resource capable of detecting diffractor points with a good resolution is desirable, but has been a challenge in the area of seismic processing. In this work, we present an approach to seismic diffraction imaging based on the diffraction operator, which can be used in both the time and depth domains, in accordance with the complexity of the area. This method that does not require any knowledge apart from the migration velocity field, i.e., rms velocities or interval velocities consists in application of pattern recognition to the amplitudes along the diffraction operator.
-
-
-
Seismic Diffraction Separation in 2D and 3D Space
More LessSeismic diffractions contain valuable information of small objects underground. The main porosities in carbonate reservoirs with secondary storage space are dissolution caves and fractures, which generate diffraction waves. Therefore, it is valuable to research on diffraction separately. There exists a significant distinction between seismic diffractions and reflections in the post-migration dip-angle domain CRP. The difference response in dip-angle gather between the two types of waves is analyzed, and a flexible approach, which can separate the diffraction from refection efficiently, is described and illustrated in this abstract. The corresponding responses of diffraction and reflection in dip-angle gather for 2D and 3D cases are derived. This approach is proven to be effective and stable by applying it to two physical models data.
-