scholarly journals Time-lapse full-waveform inversion with ocean-bottom-cable data: Application on Valhall field

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. R225-R235 ◽  
Author(s):  
Di Yang ◽  
Faqi Liu ◽  
Scott Morton ◽  
Alison Malcolm ◽  
Michael Fehler

Knowledge of changes in reservoir properties resulting from extracting hydrocarbons or injecting fluid is critical to future production planning. Full-waveform inversion (FWI) of time-lapse seismic data provides a quantitative approach to characterize the changes by taking the difference of the inverted baseline and monitor models. The baseline and monitor data sets can be inverted either independently or jointly. Time-lapse seismic data collected by ocean-bottom cables (OBCs) in the Valhall field in the North Sea are suitable for such time-lapse FWI practice because the acquisitions are of a long offset, and the surveys are well-repeated. We have applied independent and joint FWI schemes to two time-lapse Valhall OBC data sets, which were acquired 28 months apart. The joint FWI scheme is double-difference waveform inversion (DDWI), which inverts differenced data (the monitor survey subtracted by the baseline survey) for model changes. We have found that DDWI gave a cleaner and more easily interpreted image of the reservoir changes compared with that obtained with the independent FWI schemes. A synthetic example is used to demonstrate the advantage of DDWI in mitigating spurious estimates of property changes and to provide cross validations for the Valhall data results.

Geophysics ◽  
2021 ◽  
pp. 1-60
Author(s):  
Wei Zhou ◽  
David Lumley

Full waveform inversion (FWI) can be applied to time-lapse (4D) seismic data for subsurface reservoir monitoring. However, non-repeatability (NR) issues can contaminate the data and cause artifacts in the estimation of 4D rock and fluid property changes. Therefore, evaluating and studying the NR effects on the 4D data and FWI results can help, for instance, discriminate inversion artifacts from true changes, guide seismic survey design and processing workflow. Using realistic reservoir models, data and field measurements of NR, we show the effects of NR source-receiver position and seawater velocity changes on the data and the 4D FWI results. We find that ignoring these NR effects in the inversion can cause strong artifacts in the estimated velocity change models, and thus should be addressed before or during inversion. The NR source-receiver positioning issue can be addressed by 4D FWI successfully, whereas the NR water velocity issue requires measurements or estimations of water velocities. Furthermore, we compare the accuracy and robustness of the parallel, double-difference and central-difference 4D FWI methods to realistic NR ocean-bottom node data in a quantitative way. Parallel 4D FWI fails to capture geomechanical changes and also overestimates the aquifer layer changes with NR data. Double-difference 4D FWI is capable of recovering the geomechanical changes, but is also sensitive to NR noises, generating more artifacts in the overburden. By averaging the forward and reverse bootstrap 4D estimates, central-difference 4D FWI is more robust to NR noises, and also produces the most accurate 4D estimates.


Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. B311-B324 ◽  
Author(s):  
Laura Gassner ◽  
Tobias Gerach ◽  
Thomas Hertweck ◽  
Thomas Bohlen

Evidence for gas-hydrate occurrence in the Western Black Sea is found from seismic measurements revealing bottom-simulating reflectors (BSRs) of varying distinctness. From an ocean-bottom seismic data set, low-resolution traveltime-tomography models of P-wave velocity [Formula: see text] are constructed. They serve as input for acoustic full-waveform inversion (FWI), which we apply to derive high-resolution parameter models aiding the interpretation of the seismic data for potential hydrate and gas deposits. Synthetic tests indicate the applicability of the FWI approach to robustly reconstruct [Formula: see text] models with a typical hydrate and gas signature. Models of S-wave velocity [Formula: see text] containing a hydrate signature can only be reconstructed when the parameter distribution of [Formula: see text] is already well-known. When we add noise to the modeled data to simulate field-data conditions, it prevents the reconstruction of [Formula: see text] completely, justifying the application of an acoustic approach. We invert for [Formula: see text] models from field data of two parallel profiles of 14 km length with a distance of 1 km. Results indicate a characteristic velocity trend for hydrate and gas occurrence at BSR depth in the first of the analyzed profiles. We find no indications for gas accumulations below the BSR on the second profile and only weak indications for hydrate. These differences in the [Formula: see text] signature are consistent with the reflectivity behavior of the migrated seismic streamer data of both profiles in which a zone of high-reflectivity amplitudes is coincident with the potential gas zone derived from the FWI result. Calculating saturation estimates for the potential hydrate and gas zones yields values of up to 30% and 1.2%, respectively.


Solid Earth ◽  
2019 ◽  
Vol 10 (3) ◽  
pp. 765-784 ◽  
Author(s):  
Andrzej Górszczyk ◽  
Stéphane Operto ◽  
Laure Schenini ◽  
Yasuhiro Yamada

Abstract. Imaging via pre-stack depth migration (PSDM) of reflection towed-streamer multichannel seismic (MCS) data at the scale of the whole crust is inherently difficult. This is because the depth penetration of the seismic wavefield is controlled, firstly, by the acquisition design, such as streamer length and air-gun source configuration, and secondly by the complexity of the crustal structure. Indeed, the limited length of the streamer makes the estimation of velocities from deep targets challenging due to the velocity–depth ambiguity. This problem is even more pronounced when processing 2-D seismic data due to the lack of multi-azimuthal coverage. Therefore, in order to broaden our knowledge about the deep crust using seismic methods, we present the development of specific imaging workflows that integrate different seismic data. Here we propose the combination of velocity model building using (i) first-arrival tomography (FAT) and full-waveform inversion (FWI) of wide-angle, long-offset data collected by stationary ocean-bottom seismometers (OBSs) and (ii) PSDM of short-spread towed-streamer MCS data for reflectivity imaging, with the former velocity model as a background model. We present an application of such a workflow to seismic data collected by the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) and the Institut Français de Recherche pour l'Exploitation de la Mer (IFREMER) in the eastern Nankai Trough (Tokai area) during the 2000–2001 Seize France Japan (SFJ) experiment. We show that the FWI model, although derived from OBS data, provides an acceptable background velocity field for the PSDM of the MCS data. From the initial PSDM, we refine the FWI background velocity model by minimizing the residual move-outs (RMOs) picked in the pre-stack-migrated volume through slope tomography (ST), from which we generate a better-focused migrated image. Such integration of different seismic datasets and leading-edge imaging techniques led to greatly improved imaging at different scales. That is, large to intermediate crustal units identified in the high-resolution FWI velocity model extensively complement the short-wavelength reflectivity inferred from the MCS data to better constrain the structural factors controlling the geodynamics of the Nankai Trough.


2020 ◽  
Vol 224 (1) ◽  
pp. 306-311
Author(s):  
Dirk Philip van Herwaarden ◽  
Michael Afanasiev ◽  
Solvi Thrastarson ◽  
Andreas Fichtner

SUMMARY We present a new approach to full-waveform inversion (FWI) that enables the assimilation of data sets that expand over time without the need to reinvert all data. This evolutionary inversion rests on a reinterpretation of stochastic Limited-memory Broyden–Fletcher–Goldfarb–Shanno (L-BFGS), which randomly exploits redundancies to achieve convergence without ever considering the data set as a whole. Specifically for seismological applications, we consider a dynamic mini-batch stochastic L-BFGS, where the size of mini-batches adapts to the number of sources needed to approximate the complete gradient. As an illustration we present an evolutionary FWI for upper-mantle structure beneath Africa. Starting from a 1-D model and data recorded until 1995, we sequentially add contemporary data into an ongoing inversion, showing how (i) new events can be added without compromising convergence, (ii) a consistent measure of misfit can be maintained and (iii) the model evolves over times as a function of data coverage. Though applied retrospectively in this example, our method constitutes a possible approach to the continuous assimilation of seismic data volumes that often tend to grow exponentially.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. R109-R119 ◽  
Author(s):  
Timothy J. Sears ◽  
Penny J. Barton ◽  
Satish C. Singh

Elastic full waveform inversion of multichannel seismic data represents a data-driven form of analysis leading to direct quantification of the subsurface elastic parameters in the depth domain. Previous studies have focused on marine streamer data using acoustic or elastic inversion schemes for the inversion of P-wave data. In this paper, P- and S-wave velocities are inverted for using wide-angle multicomponent ocean-bottom cable (OBC) seismic data. Inversion is undertaken using a two-dimensional elastic algorithm operating in the time domain, which allows accurate modeling and inversion of the full elastic wavefield, including P- and mode-converted PS-waves and their respective amplitude variation with offset (AVO) responses. Results are presented from the application of this technique to an OBC seismic data set from the Alba Field, North Sea. After building an initial velocity model and extracting a seismic wavelet, the data are inverted instages. In the first stage, the intermediate wavelength P-wave velocity structure is recovered from the wide-angle data and then the short-scale detail from near-offset data using P-wave data on the [Formula: see text] (vertical geophone) component. In the second stage, intermediate wavelengths of S-wave velocity are inverted for, which exploits the information captured in the P-wave’s elastic AVO response. In the third stage, the earlier models are built on to invert mode-converted PS-wave events on the [Formula: see text] (horizontal geophone) component for S-wave velocity, targeting first shallow and then deeper structure. Inversion of [Formula: see text] alone has been able to delineate the Alba Field in P- and S-wave velocity, with the main field and outlier sands visible on the 2D results. Inversion of PS-wave data has demonstrated the potential of using converted waves to resolve shorter wavelength detail. Even at the low frequencies [Formula: see text] inverted here, improved spatial resolution was obtained by inverting S-wave data compared with P-wave data inversion results.


Sign in / Sign up

Export Citation Format

Share Document