SUMMARY
Salt structures are widely distributed in many basins worldwide and play an important role in understanding tectonic movements, offering underground storage and sealing hydrocarbon traps. In ...addition to the acknowledged difficulties in mapping complex salt structures through seismic methods, when an evaporitic layer, such as anhydrite, forms over salt, it can introduce strong multimode conversions that can couple with the primary compressional wavefields and generate artefacts in resulting acoustic images. From two well-log suites from the Gulf of Mexico, we identify thin evaporitic caprocks on top of their salt bodies and analyse their elastic properties. Through controlled experiments, including physical and numerical modelling with a vertical seismic profiling survey geometry, we observe significant shear-mode conversions at the top of the ultra-high-velocity caprock, which further result in a family of prominent S and P (converted from S wave) modes across the top of the salt region. Similarly, in a field survey, we identify evident converted S waves and a multimode P wave (converted S wave in the anhydrite layer, and P wave elsewhere) following the primary P transmission inside the salt body. While separating the converted Smodes at the receiver end is more unambiguous, the multimode P waves could behave very similarly to the primary Pmodes and are more difficult to suppress. Under the common acoustic assumption of seismic velocity model building and imaging, complex mode conversions on top of the salt are generally ignored. Through controlled experiments and a field survey, we analyse the pitfalls associated with this omission. We emphasize the importance of understanding the physics of wave partitioning in the presence of a thin ultra-high-velocity layer on the top of the salt.
Variations in the topography and thickness of the near-surface low-velocity weathering zone negatively impact the quality of onshore seismic data and images. Such impact can be mitigated using ...accurate near-surface velocity models. We propose to estimate near-surface velocities using a layer-cell tomography method which sequentially combines the model parameterizations of multiscale deformable-layer tomography (DLT) and cell tomography. We first estimate the long-wavelength velocity variations using multiscale DLT, which maps the undulating geometry of the weathering base where the angular coverage of first-arrival raypaths is limited. Then we convert the solution model of multiscale DLT into a cell model, and further determine the short-wavelength velocity variations using cell tomography. We test the proposed layer-cell tomography method using synthetic datasets and a field dataset. The results indicate that the new method improves the resolution of the velocity model from the multiscale DLT, and reduces the artifacts in the velocity solution compared to those observed in the cell tomography method. Reverse time migration of the synthetic data shows that using the velocity model from the layer-cell tomography method gives more accurate images of reflectors. For the field dataset test, the layer-cell tomography method yields a high-resolution near-surface velocity model that could be validated by the improvements in the tomostatic correction results.
Full text
Available for:
DOBA, EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, IZUM, KILJ, KISLJ, MFDPS, NLZOH, NUK, OBVAL, OILJ, PILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UILJ, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
In complex near-surface geologic settings, first-arrival traveltime tomography suffers from strong data noise and uneven raypath coverage. To reduce the inversion non-uniqueness, model smoothing is ...necessary but often without incorporating structural constraints based on seismic images. The reasons include lack of shallow reflections, strong noise, and gradient velocity field with strong lateral variation. We propose to estimate the near-surface structural trend from the velocity solution of deformable-layer tomography (DLT), that inverts for velocity contours directly. Then the structural dips derived from the DLT solution are taken as the structural regularization in the subsequent traveltime tomographic inversion. Tests for the 2011 SAGEEP workshop model illustrate that this DLT-guided structural regularization steers the tomographic inversion toward more accurate near-surface models. The application of the DLT-guided structural regularization to a field data reveals the near-surface geology of a beach area in Qingdao, China. The DLT solution could complement seismic images or other structural information in regularizing tomographic inversion.
Reverse time migration (RTM) is a seismic imaging method to map the subsurface reflectivity using recorded seismic waveforms. The practice in exploration seismology has long established a two-fold ...approach of seismic imaging: Using velocity modeling building to establish the long-wavelength reference velocity models, and using seismic migration to map the short-wavelength reflectivity structures. Among various seismic migration methods for different situations, RTM is the only method that is capable to use all seismic wave types that can be computed numerically. Being initiated in early 1980's, RTM seeks an image of the subsurface reflectivity as the best match in an image space between the extrapolation of time-reversed waveform data and the prediction based on estimated velocity model and source parameters. Judging the image quality in the same space of forming the images is more advantageous than the approaches of modeling and inversion which seek the solution in the model space but judge its fitness in data space. Considering that most seismic migration applications today still use primary reflection as the only signal, the capability of RTM to use all computable wave types is unique and helpful reducing the imaging artifacts due to mistaking non-primary waves as primary reflections. Hence, we refer to those RTM algorithms using only primary reflections as the first-generation RTM methods, and the RTM algorithms making a full use of primary reflections, multiple reflections and other non-primary waveform data as the second-generation RTM methods. This paper reviews the development history of the RTM along with its major challenges, current solutions, and future perspectives.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK, ZRSKP
Seismic migrations are generally formulated as the adjoint operators of linear forward modeling and often lead to images with degraded resolution, unbalanced illumination and migration artifacts, ...especially in surveys with geologic complexity and irregular acquisition geometry. Least-squares migration (LSM) is able to mitigate these problems and produce better resolved images that are suitable for subsequent AVO/AVA inversion. However, no matter what domain LSM is implemented in, the computational cost is still several times or even one order of magnitude more than that of traditional migration. In this paper, we present an efficient image-domain least-square Kirchhoff depth migration (LSKDM), in which the Hessian matrix is approximated by a grid of point-spread-functions (PSFs). Traditional PSF computing algorithm requires a non-negligible cost caused by a successive operation of modeling and migration, and has to satisfy a sampling restriction to avoid interference between nearby PSFs. We present in this paper that, by using the ray-based Green's functions and the linear traveltime approximation, the PSFs can be constructed explicitly at a significantly reduced computational cost and is able to adapt flexible spatial sampling that is fine enough to detect small-scale illumination variation. With the constructed PSFs, we formulate an image-domain LSKDM to iteratively solve for the optimal reflectivities. Numerical tests on synthetic and field data examples demonstrate that the proposed LSKDM is highly efficient and is capable of producing images with enhanced spatial resolution and amplitude fidelity when compared with the Kirchhoff depth migration (KDM) image.
The pre-/subsalt fractured networks provide key clues to understanding the tectonic history and the subsurface reservoir evolution. Yet, they are among the most complicated targets for geophysical ...investigations. We develop a structure-oriented mapping technique to delineate the high-resolution dipping layers across the borehole region using the zero-offset vertical-seismic-profiling (ZVSP) survey, which is conventionally used to provide 1-D wave propagation information. The key information for the single-shot vertical-seismic-profiling (VSP) mapping, which is the structural dip across the borehole, is unreliable from the poor subsalt surface seismic image. Alternatively, we obtain the structural dips via reflection layer tomography. Taking advantage of the accurate interval velocities and the high-fidelity identifications from the reflection events, we manage to invert the geometry of the main reflectors identified across the wellbore. Based on such information, we propose an effective processing strategy and achieve a high-resolution image of the dipping structures around the borehole region from the ZVSP. The current result compares reasonably well with the corresponding surface seismic profile but supplies higher-resolution details of the dipping structures and fault networks below the thick evaporite caprock where the surface seismic image degrades sharply. The enhanced subsurface image encourages better structural evaluation, geologic interpretation, and future 2-D/3-D VSP survey design.
The inadequate resolution of cavity karst reservoir characterization is a key factor affecting the hydrocarbon production efficiency of ultra-deep marine carbonates. Full waveform inversion (FWI) ...using high-frequency seismic data can provide a higher resolution than conventional methods. However, computational efficiency limits its application. This paper proposed a target-oriented local FWI method for cavity karst reservoir characterization. Based on the wavefield injection and Marchenko redatuming methods, a local wavefield forward modeling operator is derived. It can obtain the local wavefield corresponding to the physical source at the surface. Based on this local wavefield reconstruction, an objective function of the local FWI is established. To enhance the reservoir boundaries in the inversion results, a stabilizing strategy that combines total variation and minimum support regularization is proposed. A synthetic data test demonstrates that the obtained model can well describe the cavity karst reservoir structures.
Accurate delineation of the salt flank is important for oil and gas exploration in areas with salt intrusion. We describe the application of deformable-layer tomography (DLT) to invert for the ...geometry of salt flank using travel times of transmitted P- and S-waves from surface sources to downhole receivers. The DLT allows us to take advantage of a common situation that the salt velocity is known and generally invariant, but we need to determine the variable geometry of salt flank. We demonstrate our new method using a physical model experiment mimicking the setup of a walkaway vertical-seismic-profiling (VSP) survey. We first use picked P-wave arrivals to invert for the salt flank geometry; the sinusoid shape of the salt flank is estimated fairly by the DLT due to the uneven P-wave raypath coverage. As an improvement, we further incorporate the S-wave arrivals in the DLT. The picking of S-wave arrivals is assisted with modeled S-wave arrivals based on the P-wave DLT model. The DLT solution using both P- and S-wave arrivals delineates the salt flank geometry more accurately than that using P-wave arrivals alone. The salt flank delineation using DLT can complement with the conventional salt proximity and migration methods. It could also serve as constraints or the initial model for the full-waveform inversion of salt geometry.
Using directional waves is advantageous in subsurface seismic imaging because such waves are localized in space and time. Hence, the recorded data contain information mostly about the illuminated ...targets. However, most seismic field data are physically generated by point sources that excite wavefield propagating into all directions and interact with all parts of the medium. How to convert point-source wavefields into a wavefield due to a directional wave packet without physically exciting it, is the subject of investigation here. In particular, we investigate the Gaussian wave packet (GWP) which is a directional wave packet localized in space and time. GWP is an exact solution of the wave equation and it differs from the widely used asymptotic solutions such as Gaussian Beams or Gaussian Packets. The spatial localization and propagation direction of GWP are controlled by parameters chosen by users. We propose a method to synthesize GWP field data in complex media using recorded shot records of point sources based on the reverse-time concept. To assess the quality of the synthesized GWP data, we study the influences from point-source wavelets and the spatial interval between point sources. Finally, we present the constructed GWP fields in two numerical examples. We also show one application of the GWP data in seismic imaging using multiples.
Full text
Available for:
DOBA, EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, IZUM, KILJ, KISLJ, MFDPS, NLZOH, NUK, OBVAL, OILJ, PILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UILJ, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ