Directional depth migration

Geophysics ◽  
1985 ◽  
Vol 50 (11) ◽  
pp. 1784-1789 ◽  
Author(s):  
J. H. Higginbotham ◽  
Y. Shin ◽  
D. V. Sukup

Complicated geologic structures such as folds, overthrusts, and salt domes can produce reflectors with dips as great as 90 degrees. Because oil and gas accumulations are often associated with these steeply dipping interfaces, accurate processing of reflection seismic information from such areas becomes an important though challenging task. The proper imaging of steeply dipping reflectors requires accurate knowledge of the velocity field through which the wavefronts propagate. Thus, velocity analysis becomes extremely important. In addition to this problem, most migration algorithms have serious difficulties when dip is greater than about 50 degrees. In this discussion, we assume the velocity field is known, the data may be stacked correctly before migration, and the chief concern is migration accuracy. We describe a method for depth migration of very steeply inclined wavefronts through inhomogeneous velocity fields. The extension of the proposed technique to migration before stack is obvious.

Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. WCA19-WCA34 ◽  
Author(s):  
Christiaan C. Stolk ◽  
Maarten V. de Hoop ◽  
William W. Symes

Recent analysis and synthetic examples have shown that many prestack depth migration methods produce nonflat image gathers containing spurious events, even when provided with a kinematically correct migration velocity field, if this velocity field is highly refractive. This pathology occurs in all migration methods that produce partial images as independent migrations of data bins. Shot-geophone prestack depth migration is an exception to this pattern: each point in the prestack image volume depends explicitly on all traces within the migration aperture. Using a ray-theoretical analysis, we have found that shot-geophone migration produces focused (subsurface-offset domain) or flat (scattering-angle domain) image gathers, provided there is a curvilinear coordinate system defining pseudodepth with respect to which the rays carrying significant energy do not turn, and that the acquisition coverage is sufficient to determine all such rays. Although the analysis is theoretical and idealized, a synthetic example suggests that its implications remain valid for practical implementations, and that shot-geophone prestack depth migration could be a particularly appropriate tool for velocity analysis in a complex structure.


Geophysics ◽  
2001 ◽  
Vol 66 (4) ◽  
pp. 1208-1216 ◽  
Author(s):  
H. Grubb ◽  
A. Tura ◽  
C. Hanitzsch

Estimating a suitable velocity field for use in prestack depth migration is inherently uncertain because of limitations on the available data and estimation techniques. This uncertainty affects both the migrated depth of structures and their amplitudes in the inverted images. These effects can be estimated by performing multiple migrations with a set of velocity fields and colocating features in the migrated images. This lets us examine the imaging procedure’s sensitivity to changes in the velocity field so we can assess both structural and amplitude uncertainties in migrated images. These two types of uncertainties affect interpretation in different ways. For instance, with structural uncertainty interpretation we consider the change in migrated location of structures when deciding on drilling locations, optimizing well trajectories, or computing uncertainty in volumetric calculations. With amplitude uncertainty or amplitude versus offset (AVO) uncertainty interpretation, we consider (1) uncertainty in crossplots of pairs of AVO attributes at a point of interest or (2) uncertainty of the attribute values along identified structures. For any interpretation informing a decision, the uncertainty can help estimate risk. Our data processing approach is based on amplitude‐preserving prestack depth migration followed by AVO inversion, or AVO migration/inversion. It is valid for estimating AVO attributes in simple to moderately complex structural settings. Our methods of assessing the effect of velocity uncertainty can also be applied when obtaining structural uncertainties for a complex overburden geology or amplitude uncertainties in conventional NMO‐based AVO analysis. They may also be applied straightforwardly to any poststack attribute analysis. Key to the approach is the availability of multiple velocity fields to generate multiple migrated images. In our application, an automatic algorithm samples possible fields, but the set of fields to consider could be generated from another source, such as interpretation.


Polymers ◽  
2021 ◽  
Vol 13 (8) ◽  
pp. 1205
Author(s):  
Ruiqi Wang ◽  
Riqiang Duan ◽  
Haijun Jia

This publication focuses on the experimental validation of film models by comparing constructed and experimental velocity fields based on model and elementary experimental data. The film experiment covers Kapitza numbers Ka = 278.8 and Ka = 4538.6, a Reynolds number range of 1.6–52, and disturbance frequencies of 0, 2, 5, and 7 Hz. Compared to previous publications, the applied methodology has boundary identification procedures that are more refined and provide additional adaptive particle image velocimetry (PIV) method access to synthetic particle images. The experimental method was validated with a comparison with experimental particle image velocimetry and planar laser induced fluorescence (PIV/PLIF) results, Nusselt’s theoretical prediction, and experimental particle tracking velocimetry (PTV) results of flat steady cases, and a good continuity equation reproduction of transient cases proves the method’s fidelity. The velocity fields are reconstructed based on different film flow model velocity profile assumptions such as experimental film thickness, flow rates, and their derivatives, providing a validation method of film model by comparison between reconstructed velocity experimental data and experimental velocity data. The comparison results show that the first-order weighted residual model (WRM) and regularized model (RM) are very similar, although they may fail to predict the velocity field in rapidly changing zones such as the front of the main hump and the first capillary wave troughs.


1991 ◽  
Vol 58 (3) ◽  
pp. 820-824 ◽  
Author(s):  
A. Bogobowicz ◽  
L. Rothenburg ◽  
M. B. Dusseault

A semi-analytical solution for plane velocity fields describing steady-state incompressible flow of nonlinearly viscous fluid into an elliptical opening is presented. The flow is driven by hydrostatic pressure applied at infinity. The solution is obtained by minimizing the rate of energy dissipation on a sufficiently flexible incompressible velocity field in elliptical coordinates. The medium is described by a power creep law and solutions are obtained for a range of exponents and ellipse eccentricites. The obtained solutions compare favorably with results of finite element analysis.


PLoS ONE ◽  
2021 ◽  
Vol 16 (3) ◽  
pp. e0247826
Author(s):  
Bård A. Bendiksen ◽  
Gary McGinley ◽  
Ivar Sjaastad ◽  
Lili Zhang ◽  
Emil K. S. Espe

Myocardial velocities carry important diagnostic information in a range of cardiac diseases, and play an important role in diagnosing and grading left ventricular diastolic dysfunction. Tissue Phase Mapping (TPM) Magnetic Resonance Imaging (MRI) enables discrete sampling of the myocardium’s underlying smooth and continuous velocity field. This paper presents a post-processing framework for constructing a spatially and temporally smooth and continuous representation of the myocardium’s velocity field from TPM data. In the proposed scheme, the velocity field is represented through either linear or cubic B-spline basis functions. The framework facilitates both interpolation and noise reducing approximation. As a proof-of-concept, the framework was evaluated using artificially noisy (i.e., synthetic) velocity fields created by adding different levels of noise to an original TPM data. The framework’s ability to restore the original velocity field was investigated using Bland-Altman statistics. Moreover, we calculated myocardial material point trajectories through temporal integration of the original and synthetic fields. The effect of noise reduction on the calculated trajectories was investigated by assessing the distance between the start and end position of material points after one complete cardiac cycle (end point error). We found that the Bland-Altman limits of agreement between the original and the synthetic velocity fields were reduced after application of the framework. Furthermore, the integrated trajectories exhibited consistently lower end point error. These results suggest that the proposed method generates a realistic continuous representation of myocardial velocity fields from noisy and discrete TPM data. Linear B-splines resulted in narrower limits of agreement between the original and synthetic fields, compared to Cubic B-splines. The end point errors were also consistently lower for Linear B-splines than for cubic. Linear B-splines therefore appear to be more suitable for TPM data.


2019 ◽  
Vol 10 (2) ◽  
pp. 459-470
Author(s):  
V. A. Kontorovich ◽  
В. V. Lunev ◽  
V. V. Lapkovsky

The article discusses the geological structure, oil‐and‐gas‐bearing capacities and salt tectogenesis of the Anabar‐Khatanga saddle located on the Laptev Sea shore. In the study area, the platform sediments are represented by the 14‐45 km thick Neoproterozoic‐Mesozoic sedimentary complexes. The regional cross‐sections show the early and middle Devonian salt‐bearing strata and associated salt domes in the sedimentary cover, which may be indicative of potential hydrocarbon‐containing structures. Diapirs reaching the ground surface can be associated with structures capable of trapping hydrocarbons, and typical anticline structures can occur above the domes buried beneath the sediments. In our study, we used the algorithms and software packages developed by A.A. Trofimuk Institute of Petroleum Geology and Geophysics (IPGG SB RAS). Taking into account the structural geological features of the study area, we conducted numerical simulation of the formation of salt dome structures. According to the numerical models, contrasting domes that reached the ground surface began to form in the early Permian and developed most intensely in the Mesozoic, and the buried diapirs developed mainly in the late Cretaceous and Cenozoic.


2011 ◽  
Vol 4 (3) ◽  
pp. 1569-1610
Author(s):  
J. J. Fürst ◽  
O. Rybak ◽  
H. Goelzer ◽  
B. De Smedt ◽  
P. de Groen ◽  
...  

Abstract. We present a novel finite difference implementation of a three-dimensional higher-order ice sheet model that performs well both in terms of convergence rate and numerical stability. In order to achieve these benefits the discretisation of the governing force balance equation makes extensive use of information on staggered grid points. Using the same iterative solver, an existing discretisation that operates exclusively on the regular grid serves as a reference. Participation in the ISMIP-HOM benchmark indicates that both discretisations are capable of reproducing the higher-order model inter-comparison results. This allows a direct comparison not only of the resultant velocity fields but also of the solver's convergence behaviour which holds main differences. First and foremost, the new finite difference scheme facilitates convergence by a factor of up to 7 and 2.6 in average. In addition to this decrease in computational costs, the precision for the resultant velocity field can be chosen higher in the novel finite difference implementation. For high precisions, the old discretisation experiences difficulties to converge due to large variation in the velocity fields of consecutive Picard iterations. Finally, changing discretisation prevents build-up of local field irregularites that occasionally cause divergence of the solution for the reference discretisation. The improved behaviour makes the new discretisation more reliable for extensive application to real ice geometries. Higher precision and robust numerics are crucial in time dependent applications since numerical oscillations in the velocity field of subsequent time steps are attenuated and divergence of the solution is prevented. Transient applications also benefit from the increased computational efficiency.


2013 ◽  
Vol 2013 ◽  
pp. 1-19 ◽  
Author(s):  
Daniel Carbunaru ◽  
Sabina Stefan ◽  
Monica Sasu ◽  
Victor Stefanescu

The mesoscale configurations are analysed associated withthesplitting process of convective cells responsible for severe weather phenomena in the south-eastern part of Romania. The analysis was performed using products from the S-band Doppler weather radar located in Medgidia. The cases studied were chosen to cover various synoptic configurations when the cell splitting process occurs. To detect the presence and intensity of the tropospheric jet, the Doppler velocity field and vertical wind profiles derived from radar algorithms were used. The relative Doppler velocity field was used to study relative flow associated with convective cells. Trajectories and rotational characteristics associated with convective cells were obtained from reflectivity and relative Doppler velocity fields at various elevations. This analysis highlights the main dynamic features associated with the splitting process of convective cells: the tropospheric jet and vertical moisture flow associated with the configuration of the flow relative to the convective cells for the lower and upper tropospheric layers. These dynamic characteristics seen in the Doppler based velocity field and in the relative Doppler velocity field to the storm can indicate further evolution of convective developments, with direct implications to very short range forecast (nowcasting).


Author(s):  
X. Zhu ◽  
J. Langhammer ◽  
D. King ◽  
E. Madtson ◽  
H. K. Helgesen ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document