3‐D imaging using higher order fast marching traveltimes

Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 604-609 ◽  
Author(s):  
Alexander Mihai Popovici ◽  
James A. Sethian

Recently, fast marching methods (FMM) beyond first order have been developed for producing rapid solutions to the eikonal equation. In this paper, we present imaging results for 3‐D prestack Kirchhoff migration using traveltimes computed using the first‐order and second‐order FMM on several 3‐D prestack synthetic and real data sets. The second order traveltimes produce a much better image of the structure. Moreover, insufficiently sampled first order traveltimes can introduce consistent errors in the common reflection point gathers that affect velocity analysis. First‐order traveltimes tend to be smaller than analytic traveltimes, which in turn affects the migration velocity analysis, falsely indicating that the interval velocity was too low.

Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. U21-U28 ◽  
Author(s):  
Weihong Fei ◽  
George A. McMechan

A new migration velocity analysis is developed by combining the speed of parsimonious prestack depth migration with velocity adjustments estimated within and across common-reflection-point (CRP) gathers. The proposed approach is much more efficient than conventional tomographic velocity analysis because only the traces that contribute to a series of CRP gathers are depth migrated at each iteration. The local interval-velocity adjustments for each CRP are obtained by maximizing the stack amplitude over the predicted (nonhyperbolic) moveout in each CRP gather; this does not involve retracing rays. At every iteration, the velocity in each pixel is updated by averaging over all the predicted velocity updates. Finally, CRP positions and orientations are updated by parsimonious migration, and rays are retraced to define new CRP gathers for the next iteration; this ensures internal consistency between the updated velocity model and the CRP gather. Because the algorithm has a gridded-model parameterization, no explicit representation or fitting of reflectors is involved. Strong lateral-velocity variations, such as those found at salt flanks, can be handled. Application to synthetic and field data sets show that the proposed algorithm works effectively and efficiently.


Axioms ◽  
2021 ◽  
Vol 10 (2) ◽  
pp. 102
Author(s):  
Maya Briani ◽  
Emiliano Cristiani ◽  
Paolo Ranut

In this paper, we propose two models describing the dynamics of heavy and light vehicles on a road network, taking into account the interactions between the two classes. The models are tailored for two-lane highways where heavy vehicles cannot overtake. This means that heavy vehicles cannot saturate the whole road space, while light vehicles can. In these conditions, the creeping phenomenon can appear, i.e., one class of vehicles can proceed even if the other class has reached the maximal density. The first model we propose couples two first-order macroscopic LWR models, while the second model couples a second-order microscopic follow-the-leader model with a first-order macroscopic LWR model. Numerical results show that both models are able to catch some second-order (inertial) phenomena such as stop and go waves. Models are calibrated by means of real data measured by fixed sensors placed along the A4 Italian highway Trieste–Venice and its branches, provided by Autovie Venete S.p.A.


Geophysics ◽  
2021 ◽  
pp. 1-50
Author(s):  
German Garabito ◽  
José Silas dos Santos Silva ◽  
Williams Lima

In land seismic data processing, the prestack time migration (PSTM) image remains the standard imaging output, but a reliable migrated image of the subsurface depends on the accuracy of the migration velocity model. We have adopted two new algorithms for time-domain migration velocity analysis based on wavefield attributes of the common-reflection-surface (CRS) stack method. These attributes, extracted from multicoverage data, were successfully applied to build the velocity model in the depth domain through tomographic inversion of the normal-incidence-point (NIP) wave. However, there is no practical and reliable method for determining an accurate and geologically consistent time-migration velocity model from these CRS attributes. We introduce an interactive method to determine the migration velocity model in the time domain based on the application of NIP wave attributes and the CRS stacking operator for diffractions, to generate synthetic diffractions on the reflection events of the zero-offset (ZO) CRS stacked section. In the ZO data with diffractions, the poststack time migration (post-STM) is applied with a set of constant velocities, and the migration velocities are then selected through a focusing analysis of the simulated diffractions. We also introduce an algorithm to automatically calculate the migration velocity model from the CRS attributes picked for the main reflection events in the ZO data. We determine the precision of our diffraction focusing velocity analysis and the automatic velocity calculation algorithms using two synthetic models. We also applied them to real 2D land data with low quality and low fold to estimate the time-domain migration velocity model. The velocity models obtained through our methods were validated by applying them in the Kirchhoff PSTM of real data, in which the velocity model from the diffraction focusing analysis provided significant improvements in the quality of the migrated image compared to the legacy image and to the migrated image obtained using the automatically calculated velocity model.


Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. VE161-VE171 ◽  
Author(s):  
J. Schleicher ◽  
J. C. Costa ◽  
A. Novais

Image-wave propagation or velocity continuation describes the variation of the migrated position of a seismic event as a function of migration velocity. Image-wave propagation in the common-image gather (CIG) domain can be combined with residual-moveout analysis for iterative migration velocity analysis (MVA). Velocity continuation of CIGs leads to a detection of those velocities in which events flatten. Although image-wave continuation is based on the assumption of a constant migration velocity, the procedure can be applied in inhomogeneous media. For this purpose, CIGs obtained by migration with an inhomogeneous macrovelocity model are continued starting from a constant reference velocity. The interpretation of continued CIGs, as if they were obtained from residual migrations, leads to a correction formula that translates residual flattening velocities into absolute time-migration velocities. In this way, the migration velocity model can be improved iteratively until a satisfactory result is reached. With a numerical example, we found that MVA with iterative image continuation applied exclusively to selected CIGs can construct a reasonable migration velocity model from scratch, without the need to build an initial model from a previous conventional normal-moveout/dip-moveout velocity analysis.


Author(s):  
J. Pegna ◽  
F.-E. Wolter

Abstract Computer Aided Geometric Design of surfaces sometimes presents problems that were not envisioned by mathematicians in differential geometry. This paper presents mathematical results that pertain to the design of second order smooth blending surfaces. Second order smoothness normally requires that normal curvatures agree along all tangent directions at all points of the common boundary of two patches, called the linkage curve. The Linkage Curve Theorem proved here shows that, for the blend to be second order smooth when it is already first order smooth, it is sufficient that normal curvatures agree in one direction other than the tangent to a first order continuous linkage curve. This result is significant for it substantiates earlier works in computer aided geometric design. It also offers simple practical means of generating second order blends for it reduces the dimensionality of the problem to that of curve fairing, and is well adapted to a formulation of the blend surface using sweeps. From a theoretical viewpoint, it is remarkable that one can generate second order smooth blends with the assumption that the linkage curve is only first order smooth. This property may be helpful to the designer since linkage curves can be constructed from low order piecewise continuous curves.


Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. S437-S447 ◽  
Author(s):  
Jean-Philippe Montel ◽  
Gilles Lambaré

Common-image gathers are a useful output of the migration process. Their kinematic behavior (i.e., the way they curve up or down) is an indicator of the quality of the velocity model used for migration. Traditionally, when used for migration velocity analysis, we pick structural dips in the common attribute panels (offset, angle, etc.) and residual moveout (RMO) in the gathers. The measured RMO will then tell us how much we need to update the velocity model to improve the gather’s flatness. Understanding the kinematics of the picked events is the key to an accurate model update. This point has been widely underestimated in many cases. For example, when dealing with angle gathers, there is a general assumption that the associated tomographic rays are fully defined by the picked structural dips and the gather opening and azimuth angle, and that if the velocity model is correctly updated down to a given horizon, it is not necessary to shoot the tomographic rays upward through this horizon. We find through an original theoretical analysis that both of these assumptions have to be modified when the gathers exhibit RMO. Using a kinematic analysis, we determine that knowledge of the RMO slopes is necessary to compute the tomographic rays.


2020 ◽  
Vol 34 (07) ◽  
pp. 11580-11587
Author(s):  
Haojie Liu ◽  
Han Shen ◽  
Lichao Huang ◽  
Ming Lu ◽  
Tong Chen ◽  
...  

Traditional video compression technologies have been developed over decades in pursuit of higher coding efficiency. Efficient temporal information representation plays a key role in video coding. Thus, in this paper, we propose to exploit the temporal correlation using both first-order optical flow and second-order flow prediction. We suggest an one-stage learning approach to encapsulate flow as quantized features from consecutive frames which is then entropy coded with adaptive contexts conditioned on joint spatial-temporal priors to exploit second-order correlations. Joint priors are embedded in autoregressive spatial neighbors, co-located hyper elements and temporal neighbors using ConvLSTM recurrently. We evaluate our approach for the low-delay scenario with High-Efficiency Video Coding (H.265/HEVC), H.264/AVC and another learned video compression method, following the common test settings. Our work offers the state-of-the-art performance, with consistent gains across all popular test sequences.


Sign in / Sign up

Export Citation Format

Share Document