The result of the data plot of the forward model of marmousi is as follows. For the parallel computing, we deal with the data of different shots simultaneously. Accuracy of first arrivals for a constant velocity gradient for a constant velocity gradient model, the analytical solution for the travel time is known sheriff and geldart, 1995. Module 2 velocity, speed and models in mechanics per wiki. Accelerating fullwaveform inversion with attenuation. The marmousi model was created in 1988 by the institut francais du petrole ifp. The first step is to load the velocity model and the observed data into julia, as well as setting up bound constraints for the inversion, which prevent too high or low velocities in the final result. Figure 1 shows a smoothed version of the marmousi model displayed in the xzdomain top and the domain bottom. The geometry of this model is based on a profile through the north quenguela trough in the cuanza basin. Gaussian filtered velocity models of the true model and a constant halfspace density model cdh are used as initial models. Building a good initial model for fullwaveform inversion. Application to waveequation imaging and tomography. The regional structure is constructed by extending the physical properties of most of the welldefined units to less wellconstrained areas in the mississippi embayment and the cus.
The agl elastic marmousi model, also commonly referred to as the marmousi2 model, is an elastic extension of the original marmousi model created by gary s. Cygwin executes the instruction, and you can see that the feedback information is color, so it is very convenient to analyze the information. The geometry and velocity model were created to produce complex seismic data which require advanced processing techniques to obtain a correct earth image. Any of these files may be downloaded to local machines using ftp protocols. Analysis of the problem of the velocity model determination for seismic imaging. The raw velocity model needs to be formatted in a similar fashion to the model section of this paper. Full waveform inversion lab king abdullah university of. A velocity model maps out the layers of rock in a monitoring. Distribution maps and fitting lines between velocity and density for marmousi2.
Any of these files may be downloaded to local machines. To run it yourself you can copytype each individual cell or directly download the full notebook, including all required files. Theonly way to make an excessively complex velocity model suitable for application of raybased methods, such as the gaussian beam or gaussian packet methods, is to smooth it. Gzipped tar file of marmousi data, velocity and density models, big endian ieee disk images of segy, approx. For the simulations, the minimum swave velocity was set at 600 ms by. Use velocity modeling to integrate measured seismic velocities, depth imaged seismic data, well velocity information and interpreted structural information into a calibrated velocity model. We used a model building technique presented by obrien et al. The repository contains several velocity models of varying smoothness as well. This model is an update to s362ani and s362wmani representing an anisotropic shear velocity model of the earths mantle using normal modes, body waves, surface waves and. We initiate inversion from a linear velocity model and search for updates of model parameters by lbfgs. Full waveform inversion in time and frequency domain of. Other seismological velocity models are avaiable for other regions. Layers of rock are created from different materials or under varying conditions and naturally each layer of rock will have different properties. A combined elastic waveform and gravity inversion for.
The grided values of velocity vary from 1500 ms 1 to 5500 ms 1. Image size 2000 px notice that very large and arbitrary shaped domains can be mapped. The 2004 bp velocityanalysis benchmark dataset from the velocityanalysis workshop at the 2004 eage meeting, paris, france. Visual explanation of python seismic data develop paper. The marmousi model, based on a real geological structure, is very complex, see figure 1. To request modeled data please contact michael welch at the agl. Marmousi inversion seismic 2 model the industry standard marmousi model was used as a test case, to validate our models ability to generalise for geometries that it has not seen before. Influence of error in estimating anisotropy parameters on. Citeseerx document details isaac councill, lee giles, pradeep teregowda. Download 20,390 velocity stock photos for free or amazingly low rates. Clearly, the velocity for this portion of the trip will be zero.
Marmousi synthetic dataset stanford exploration project. Circles, number 1, show misshaping of pinchout and a wrong dip for the reflector in the isotropic image. Open data on the seg wiki is a catalog of available open geophysical data online. The initial conditions for this leg are the same as the final conditions for the previous portion of the trip. Thus, we model it separately of the trip to the physics building. After the velocity modeling phase, domain conversion of various objects and how to model structural uncertainty is covered. Marmousi traveltime computation and imaging comparisons. Since velocity gradients account for compaction due to. Reverse time migration by interpolation and pseudonalytical methods within. The purpose of this paper is to analyze the generalized optimal velocity model proposed by nagatani7, 8. Velocity stock photos download 20,390 royalty free photos. We generated bornmodeled data with the smoothed velocity model and the true reflectivity model.
Information about the openaccess article full waveform inversion in time and frequency domain of velocity modeling in seismic imaging. In model parameterization the triangular cell size is 12 m and the metric distance for secondary node spacing is. An important aspect of microseismic monitoring is the creation of accurate velocity models for the monitoring region. Next, we test our algorithm under the marmousi velocity model, the sigsbee2a velocity model and the 2d canadian foothills velocity model. Prestack inversion outline 1 prestack inversion 2 full waveform inversion the basic idea 3 the math 4 marmousi example estimating the initial model fwi bee bednar panorama technologies seismic modeling, migration and velocity inversion may 30, 2014 3 53. The model contains many reflectors, steep dips, and strong velocity variations in both the lateral and the vertical direction with a minimum velocity of 1500 ms and a maximum of 5500 ms. The geometry and velocity model were created to produce a complex seismic data which requiere advanced processing techniques to obtain a. A single point particle, the center of mass of a rigid body or rotating body that is subject to zero net applied force or torque. The agl elastic marmousi model was created at the university of houston. For the purposes of this example a shot will be fired at 5 km along the horizontal coordinate and at a depth of 10 meters. Brochure attached to the cdrom formerly to the tape. In this section, we test the scalarpwave fwi with the wellknown 2d marmousi model.
The agl elastic marmousi model was created at the university of. Comparison of figures 5a and 5b demonstrates several differences. All data published is freely available to download where practical from. This notebook shows how to interpolate a given model defined on a regular grid to a salvus mesh using the marmousi model as an example. Due to multipathing, migration velocity analysis mva principle has to be revisited, requiring in particular the use of common angle gathers cags instead of the classical common offset panels. We have smoothed the marmousi model by choosing a coarser grid and by. Paleozoic, crystalline upper crust, lower crust, modified lower crust and upper mantle. For both methods we consider shot gathers independently and for each of them we compute two wavefields. Multiple arrival tracking within irregular triangular or. Model building the velocity model is 67km long and 12km deep, and was built on a 6. Community velocity models, within scec, are computer programs and data files that provide information about earths material properties from the surface to depths below 100km. Constant velocity and zero net force model applicable systems. Sedimentary rock forms through deposition over time.
Transform seismic time into subsurface depth as a critical step in all exploration, development and production project planning. Accurate initial model or available lowfrequency data is an important factor in the success of full waveform inversion fwi. Download scientific diagram the marmousi velocity model. The resulting forward and adjoint problems can be efficiently solved by using the fast sweeping method. We compare several migration traveltime computation methods in the complex marmousi velocity model.
In this study a twostep hierarchic joint inversion of seismic waveforms and gravity data is tested using the marmousiii model. Abstract we demonstrate a method for 2d velocity estimation in complex media characterized by triplicated ray fields. Hybrid kinematicdynamic approach to seismic waveequation. We demonstrate an application of fwi to marmousi ii benchmark velocity model. Tools for velocity analysis velocity spectrum or semblance best modern method industry standard define time window in which wavelet is contained calculate hyperbolic trajectory measure coherency of trace signal, e. Image domain leastsquares migration with a hessian matrix.
The 2004 bp velocity analysis benchmark dataset from the velocity analysis workshop at the 2004 eage meeting, paris, france. The data set includes subsurface models for density, pwave and swave velocity, as well as a series of processed data sets nmo stacks, pre and poststack time and depth migrated sections. Download scientific diagram marmousi model velocity field. The results of velocity model drawing are as follows. Scecs community velocity models cvms provide detailed 3d properties for southern california. Scons script that generates rsf formatted amoco velocity model.
The velocity model is shown in figure 7a, and the true normalincidence reflectivity model derived from the velocity model is shown in figure 7b. Furthermore, we define an 8 hertz ricker wavelet as the source function. Wavefield reconstruction by interpolating significantly decimated boundaries. We show the benefits of the proposed hybrid calculationandstoring windowed technique for speeding up the reversetime migration, 14 and the waveequation tomography. Receivers are spread at a depth of 25 meters every. Now we use a more complex model the marmousi model to further investigate the viability and validity of the multistage triangular shortestpath method algorithm for ray tracing. The above test is based on the simple velocity model. This is a good line to start experimenting with open data and su because there is a full processing sequence including data download, reformat, header loading, gain, prestack fk filter, brute velocity estimation, brute stack, residual statics, final velocity analysis, final stack, and two types of post stack migration phase shift and. In step 1 fwi is performed for all elastic parameters. We then demonstrate the capabilities of our lagrangian solver using the marmousi model from the 1996 inria workshop on multiarrival traveltimes see. Numerical examples using the marmousi model and a wide azimuth field data show that the proposed method can achieve reasonable accuracy and significantly enhances the. The lowfrequency helps determine the kinematical relevant components, lowwavenumber of the velocity model, which are in turn needed. Figure 5 illustrates the images of vti and isotropic psdm for marmousi model as well as the overlay of the images with the reflectivity model which was derived from the velocity model and its calculated acoustic impedance.
The dimensions of the model are 9200 metres length by 3000 metres depth. The purpose of this course is to introduce the student to velocity modeling, velocity handling and domain conversion functionality in petrel. Full waveform inversion for velocity and density with rock. Bandlimited greens functions, paraxial ray tracing, norsar wavefront construction, zhangs local ray tracing, van trier symes upwind finitedifference scheme, and podvins eikonal solver.
1622 528 658 717 1098 555 1095 315 1422 535 929 1508 643 30 977 961 406 302 917 1315 504 1318 4 843 589 1392 187 700 1356 873 773 980 1002 945 397 328 916 1499 277 1410 1374 560