Residual curvature migration velocity analysis for angle domain common imaging gathers

Pre-stack depth migration velocity analysis is one of the keys to influencing the imaging quality of pre-stack migration.In this paper we cover a residual curvature velocity analysis method on angle-domain common image gathers（ADCIGs） which can depict the relationship between incident angle and migration depth at imaging points and update the migration velocity.Differing from offset-domain common image gathers（ODCIGs）,ADCIGs are not disturbed by the multi-path problem which contributes to imaging artifacts,thus influencing the velocity analysis.On the basis of horizontal layers,we derive the residual depth equation and also propose a velocity analysis workflow for velocity scanning.The tests to synthetic and field data prove the velocity analysis methods adopted in this paper are robust and valid.

Multiwave velocity analysis based on Gaussian beam prestack depth migration

Prestack depth migration of multicomponent seismic data improves the imaging accuracy of subsurface complex geological structures. An accurate velocity field is critical to accurate imaging. Gaussian beam migration was used to perform multicomponent migration velocity analysis of PP- and PS-waves. First, PP- and PS-wave Gaussian beam prestack depth migration algorithms that operate on common-offset gathers are presented to extract offsetdomain common-image gathers of PP- and PS-waves. Second, based on the residual moveout equation, the migration velocity fields of P- and S-waves are updated. Depth matching is used to ensure that the depth of the target layers in the PP- and PS-wave migration profiles are consistent, and high-precision P- and S-wave velocities are obtained. Finally, synthetic and field seismic data suggest that the method can be used effectively in multiwave migration velocity analysis.

Numerical Simulation for Accuracy of Velocity Analysis in Small-Scale High-Resolution Marine Multichannel Seismic Technology

When used with large energy sparkers, marine multichannel small-scale high-resolution seismic detection technology has a high resolution, high-detection precision, a wide applicable range, and is very flexible. Positive results have been achieved in submarine geological research, particularly in the investigation of marine gas hydrates. However, the amount of traveltime difference information is reduced for the velocity analysis under conditions of a shorter spread length, thus leading to poorer focusing of the velocity spectrum energy group and a lower accuracy of the velocity analysis. It is thus currently debatable whether the velocity analysis accuracy of short-arrangement multichannel seismic detection technology is able to meet the requirements of practical application in natural gas hydrate exploration. Therefore, in this study the bottom boundary of gas hydrates(Bottom Simulating Reflector, BSR) is used to conduct numerical simulation to discuss the accuracy of the velocity analysis related to such technology. Results show that a higher dominant frequency and smaller sampling interval are not only able to improve the seismic resolution, but they also compensate for the defects of the short-arrangement, thereby improving the accuracy of the velocity analysis. In conclusion, the accuracy of the velocity analysis in this small-scale, high-resolution, multi-channel seismic detection technology meets the requirements of natural gas hydrate exploration.

Four-parameter velocity analysis and its applications to the Ken-71 area converted-wave data

3-D converted-wave data were acquired using digital MEMS （micro-electromechanical system） three component （3C） sensors in the alternating sand and shale sequence in the overburden of the Shengli Ken-71 area. This gives rise to serious non-hyperbolic moveout effects in the converted-wave data due to both the asymmetrical ray path and anisotropic effects. Conventional velocity analysis and moveout correction based on isotropic methods do not flatten reflections events. Here, we use a four-parameter theory to evaluate these effects and process the data. These four parameters include the PS converted wave stacking velocity （Vc2）, the vertical velocity ratio （Y0）, the effective velocity ratio （Yeff）, and the anisotropy parameter （xoff）, The method utilizes the moveout information at different offsets to estimate the different parameters and ensures that the events are properly aligned for stacking, As a result, this four-parameter theory leads to an improvement in imaging quality and correlation between the P-waves and converted-waves.

Differential interformational velocity analysis as an effective direct hydrocarbon indicator under velocity reversal conditions,an example from the anomalously high temperature and over-pressured DF1-1 Gas Field in the Yinggehai Basin,South China Sea

In the DF1-1 Gas Field in the Yinggehai Basin, South China Sea, the velocity-depth plot and velocity spectra show significant variations from a linear trend, exhibiting a distinct reversal phenomenon. Velocity parameters derived from velocity spectral analysis of the seismic data and sonic logs indicate that the interval velocity reverses below 2,100 m (2.2 s two-way time (TWT)) in the DF1-1 Gas Field. Some direct hydrocarbon indicators (DHIs) models developed for the shallow strata cannot be simply applied to the moderately to deeply buried strata for direct exploration target recognition because the velocity reversal has caused the middle-deep gas reservoirs to exhibit a moderate or weak seismic amplitude. The hydrocarbon indicator method of “Differential Interformational Velocity Analysis (DIVA)” with the aid of other hydrocarbon indicating techniques was employed to effectively identify DHIs in the middle-deep strata under velocity inversion. The result has shown that the DIVA technique can be effectively used as a DHI in both the shallow and the middle-deep strata in the study area with the shallow strata characterized by Type I DIVA anomaly and the middle-deep strata characterized by the Type II DIVA anomaly.

VELOCITY ANALYSIS OF A PLANAR MECHANISM USING A METHOD FROM STATICS

A new method for solving the velocity analysis of a mechanism is presented. Central to the method is how to convert the probem of velocity analysis of a mechanism to one of static analysis. Application of the method to certain practical problems has advantages compared to conventional methods for both graphical and analytical solutions. For brevity an example of a planar mechanism only is presented.

VELOCITY ANALYSIS OF A SPATIALM ECHANISM USING THE PRINCIPLEOF VIRTUAL FORCES

A new method of both analysis and graphics is presented for solving the problem of velocity analysis of a spatial four-bar mechanism. Central to the method is how to use the principle of virtual forces equilibrium system connected with the principle of virtual velocity to solve the velocity analysis of a mechanism. The method is accurate in principle and much simpler than the conventional method. It can be applied to both planar and spatial mechanisms. For brevity an example of a spatial mechanism only is presented.

Velocity anomaly analysis on low-amplitude structure and its application in Bohai Sea:a case study on QHD33-1S structure

Varying degree of velocity anomalies has been appeared in Shijiutuo uplift in Bohai Sea,which is mainly demonstrated in the inconsistent between seismic interpretation and the actual drilling depth.In this paper, QHD33-1S area is taken as the example.First of all,the main reason that causes the velocity anomaly phenomena is the effect of sand-formation ratio by forward model analysis,and then technical approaches to improve the accuracy of the velocity field and the degree of understanding of anomalies are further explored,resulting in more precise determination of horizontal variation trend of the velocity in QHD33-1S area.Consequences of actual applications indicate that through the accurate analysis of the velocity anomaly,we can not only realize the fine description of low-amplitude structure,but also the effective prediction of the hydrocarbon-bearing properties of the reservoir.Meanwhile,the results also have a certain significance for the other low-amplitude structures in Bohai Sea.

THE APPLICATION OF COMPATIBLE STRESS ITERATIVE METHOD IN DYNAMIC FINITE ELEMENT ANALYSIS OF HIGH VELOCITY IMPACT

There is a common difficulty in elastic-plastic impact codes such as EPIC[2,3] NONSAP[4], etc.. Most of these codes use the simple linear functions usually taken from static problem to represent the displacement components. In such finite element formulation, the stress components are constant in each element and they are discontinuous in any two neighboring elements. Therefore, the bases of using the virtual work principle in such elements are unreliable. In this paper, we introduce a new method, namely, the compatible stress iterative method, to eliminate the above-said difficulty. The calculated examples show that the calculation using the new method in dynamic finite element analysis of high velocity impact is valid and stable, and the element stiffness can be somewhat reduced.

Combined migration velocity model-building and its application in tunnel seismic prediction

We propose a combined migration velocity analysis and imaging method based on Kirchhoff integral migration and reverse time migration,using the residual curvature analysis and layer stripping strategy to build the velocity model.This method improves the image resolution of Kirchhoff integral migration and reduces the computations of the reverse time migration.It combines the advantages of efficiency and accuracy of the two migration methods.Its application in tunnel seismic prediction shows good results.Numerical experiments show that the imaging results of reverse time migration are better than the imaging results of Kirchhoff integral migration in many aspects of tunnel prediction.Field data show that this method has efficient computations and can establish a reasonable velocity model and a high quality imaging section.Combination with geological information can make an accurate prediction of the front of the tunnel geological structure.

A Quadratic precision generalized nonlinear global optimization migration velocity inversion method

An important research topic for prospecting seismology is to provide a fast accurate velocity model from pre-stack depth migration. Aiming at such a problem, we propose a quadratic precision generalized nonlinear global optimization migration velocity inversion. First we discard the assumption that there is a linear relationship between residual depth and residual velocity and propose a velocity model correction equation with quadratic precision which enables the velocity model from each iteration to approach the real model as quickly as possible. Second, we use a generalized nonlinear inversion to get the global optimal velocity perturbation model to all traces. This method can expedite the convergence speed and also can decrease the probability of falling into a local minimum during inversion. The synthetic data and Mamlousi data examples show that our method has a higher precision and needs only a few iterations and consequently enhances the practicability and accuracy of migration velocity analysis （MVA） in complex areas.

Movement analysis on steel wire rope of continuous conveyor with disc-tube assembly

The steel wire rope of continuous conveyor with disc-tube assembly is droved by the driving wheel. When the driving wheel rotates, the gear is combined to the con- nection disc in turn, promoting the connection disc to move in succession. Turning the whirling torque of driving wheel into the straight-line traction force. When the steel wire rope is winded by the driving wheel some winded along the circumference, others winded along the straight line. Used motion subject law, this article analyses the change of the velocity and the acceleration of the steel wire rope in the straight movement, and observe the mathematics’ model of velocity and acceleration.

Parametrical analysis of the railways dynamic response at high speed moving loads

The paper introduces some findings about a sensitivity analysis conducted on every geometrical and mechanical parameters which characterize the use of a railway superstructure at the high velocity. This analysis was carried out by implementing a forecast model that is derived from the simplified Gazetas and Dobry one. This model turns out to be particularly appropriate in the explication of problems connected to high velocity, since it evaluates both inertial and viscous effects activated by the moving load speed. The model implementation requires the transfer function determination that represents the action occurred by the bed surfaces on the railway and it therefore contains information concerning the geometrical and the mechanical characteristics of the embankment, of the ballast and of the sub-ballast. The transfer function H has been evaluated with the finite elements method and particularly, by resorting the ANSYS code with a harmonic structural analysis in the frequencies field. The authors, from the critic examination of the system＇s dynamics response in its entirety, glean a series of observations both of a general and a specific character, finally attaining a propose of a design modification of the standard railway superstructure at the high velocity of train operation adopted today especially in Italy.

Multi-decadal variations in glacier flow velocity and the influencing factors of Urumqi Glacier No.1 in Tianshan Mountains, Northwest China

Urumqi Glacier No. 1 is a representative glacier in the inland areas of Central Asia and is the only Chinese reference glacier in the World Glacier Monitoring Service. In this study, we explored multi-decadal variations in the flow velocity of the glacier and the influencing factors based on continuous field observations and path coefficient analysis. Results show that the glacier flow velocity decreased from 5.5 m/a in 1980/1981 to 3.3 m/a in 2010/2011. The annual variation in the direction of glacier flow velocity in the western branch and eastern branch was less than 1°–3°, and the change of glacier flow velocity in the western branch was more dramatic than that in the eastern branch. Glacier flow velocity was influenced by glacier morphology（including glacier area, glacier length, and ice thickness）, glacier mass balance and local climate conditions（air temperature and precipitation）, the glacier morphology being the leading factor. The long-term flow velocity data set of Urumqi Glacier No. 1 contributes to a better understanding of glacier dynamics within the context of climatic warming.

Kirchhoff PSDM angle-gather generation based on the traveltime gradient

Angle-domain common-image gathers（ADCIGs） are the basic data in migration velocity analysis（MVA） and amplitude variation with angle（AVA） analysis. We propose a common-angle gather-generating scheme using Kirchhoff PSDM based on the traveltime gradient field. The scheme includes three major operations：（1） to calculate the traveltime field of the source and the receiver based on the dynamic programming approach;（2） to obtain the refl ection angle according to the traveltime gradient field in the image space; and（3） to generate the ADCIGs during the migration process. Because of the computation approach, the method for generating ADCIGs is superior to conventional ray-based methods. We use the proposed ADCIGs generation method in 3D large-scale seismic data. The key points of the method are the following.（1） We use common-shot datasets for migration,（2） we load traveltimes based on the shot aperture, and（3） we use the MPI and Open Mp memory sharing to decrease the amount of input and output（I/O）. Numerical examples using synthetic data suggest that the ADCIGs improve the quality of the velocity and the effectiveness of the 3D angle-gather generation scheme.

Tunnel seismic tomography method for geological prediction and its application

Typical existing methods of tunnel geological prediction include negative apparent velocity, horizontal seismic profile, and the Tunnel Seismic Prediction （TSP） method as this technology is under development at home and abroad. Considering simpler observational methods and data processing, it is hard to accurately determine the seismic velocity of the wall rock in the front of the tunnel face. Therefore, applying these defective methods may result in inaccurate geological inferences which will not provide sufficient evidence for classifying the wall rock characteristics. This paper proposes the Tunnel Seismic Tomography （TST） method using a spatial observation arrangement and migration and travel time inversion image processing to solve the problem of analyzing the velocity structure of wall rock in the front of the tunnel face and realize accurate imaging of the geological framework of the tunnel wall rock. This method is very appropriate for geological prediction under complex geological conditions.

Datum Level and NMO Corrections to Shallow Seismic Reflection Data on Rugged Topography

The application of the seismic reflection method is often limited in complex terrain areas. The problem is the incorrect correction of time-shift caused by topography. To apply normal moveout (NMO) correction to reflection data correctly, static corrections are necessary to be applied in advance for the compensation of the time distortions of topography and the time-delays from near-surface weathered layers. For environment and engineering investigation, weathered layers are our targets so that the static correction mainly serves the adjustment of time-shift due to an undulating surface. In practice, seismic reflected raypaths are assumed to be almost vertical through the near-surface layers because they have much lower velocities than layers below. This assumption is typically acceptable in most cases since it results in little residual error for small elevation changes and small offsets in reflection events. Although static algorithms based on choosing a floating datum related to common midpoint gathers or residual surface-consistent functions are available and effective, errors caused by the assumption of vertical raypaths often generate pseudo indications of structures. This paper presents the comparison of applying corrections based on the vertical raypaths and bias (non-vertical) raypaths. It also provides an approach of combining elevation and NMO corrections. The advantages of the approach are demonstrated by a synthetic example of multi-coverage seismic reflection surveys on rough topography.

Automatic estimation of traveltime parameters in VTI mediausing similarity-weighted clustering

Compared with hyperbolic velocity estimation methods,nonhyperbolic methods(such as shifted hyperbola)are better choices for large offsets or vertical transverse isotropy(VTI)media.Since local seismic event slope contains subsurface information,they can be used to estimate zero-offset two-way traveltime and normal moveout velocity.The traditional velocity estimation methods require a great deal of manual work and are also prone to human error.In order to estimate the traveltime parameters for VTI media automatically,in this paper,we propose to use predictive painting and similarity-weighted clustering to obtain traveltime parameters.The predictive painting is used to estimate zero-offset two-way traveltime,and the shifted-hyperbola traveltime equation is used to obtain velocity and anisotropy attributes.We first map local slopes to zero-offset two-way traveltime and moveout-parameters domain and then use similarity-weightedk-means clustering to find the maximum likelihood anisotropy parameters of the main subsurface structures.In order to demonstrate that,we apply the similarity-weighted clustering method to synthetic and field data examples and the results are of higher accuracy when compared to the ones obtained using multiparameter semblance-based method.From estimation error section,it can be seen that the estimation error of multiparameter semblance-based method is about 3-5 times that of the proposed method.

Structural Effect on Hydrocarbon Trapping, North Western Nile Delta, Egypt

The present study aims to interpret the available data, that derived from well logs and 2-D seismic reflection data, in the area located northeast onshore Nile Delta, in order to delineate its structural styles and its relation with the basement inferences and regional tectonic, and also to detect the possible existence of hydrocarbon accumulation at the Tertiary-Quaternary sedimentary section. The study area lies between latitudes 30°31' and 31°17'N, and longitudes 30°38' and 31°21'E. Using the available well logs, the average velocity gradients are calculated at El-Wastani, Kafr El Sheikh, Abu Madi, Qawasim and Sidi Salem Formations and velocity contour maps are constructed. The available 2-D seismic sections were interpreted and utilized to construct the isochronous reflection maps to the tops of the fore-mentioned five formations. These maps were converted into depth maps. The depth maps show the predominance of three tectonic trends. The structure contour maps were used to deduce the structural elements affecting the fore-mentioned formations in the area under study. These maps reflect the abundance of the E-W trending faults (sometimes ENE-WSW and in other times WNW-ESE), which characterize the general Mediterranean Sea trend. The second system of structural deformations is trending NE-SW (Syrian arc trend). The third type of faults is the NW-SE trending features (Red Sea faults trend). By integrating the structure contour maps together with the different kinds of maps have been done using different sources of data such as velocity data, the possible locations that can be oil traps for hydrocarbon accumulation have been detected.

Noise reduction in LOS wind velocity of Doppler lidar using discrete wavelet analysis

The line of sight (LOS) wind velocity can be determined from the incoherent Doppler lidar backscattering signals. Noise and interference in the measurement greatly degrade the inversion accuracy. In this paper, we apply the discrete wavelet denoising method by using biorthogonal wavelets and adopt a distance-dependent thresholds algorithm to improve the accuracy of wind velocity measurement by incoherent Doppler lidar. The noisy simulation data are processed and compared with the true LOS wind velocity. The results are compared by the evaluation of both the standard deviation and correlation coefficient. The results suggest that wavelet denoising with distance-dependent thresholds can considerably reduce the noise and interfering turbulence for wind lidar measurement.