Walk‐away VSP using drill noise as a source

Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 978-997 ◽  
Author(s):  
Jacob B. U. Haldorsen ◽  
Douglas E. Miller ◽  
John J. Walsh

We describe a method for extracting and deconvolving a signal generated by a drill bit and collected by an array of surface geophones. The drill‐noise signature is reduced to an effective impulse by means of a multichannel Wiener deconvolution technique, producing a walk‐away reverse vertical seismic profile (VSP) sampled almost continuously in depth. We show how the multichannel technique accounts for noise and for internal drill‐string reflections, automatically limiting the deconvolved data to frequencies containing significant energy. We have acquired and processed a data set from a well in Germany while drilling at a depth of almost 4000 m. The subsurface image derived from these data compares well with corresponding images from a 3-D surface seismic survey, a zero‐offset VSP survey, and a walk‐away VSP survey acquired using conventional wireline techniques. The effective bandwidth of the deconvolved drill‐noise data is comparable to the bandwidth of surface seismic data but significantly smaller than what can be achieved with wireline VSP techniques. Although the processing algorithm does not require the use of sensors mounted on the drill string, these sensors provide a very economic way to compress the data. The sensors on the drill string were also used for accurate timing of the deconvolved drill‐noise data.

Geophysics ◽  
1994 ◽  
Vol 59 (10) ◽  
pp. 1500-1511 ◽  
Author(s):  
Jakob B. U. Haldorsen ◽  
Douglas E. Miller ◽  
John J. Walsh

We describe a technique for performing optimal, least‐squares deconvolution of vertical seismic profile (VSP) data. The method is a two‐step process that involves (1) estimating the source signature and (2) applying a least‐squares optimum deconvolution operator that minimizes the noise not coherent with the source signature estimate. The optimum inverse problem, formulated in the frequency domain, gives as a solution an operator that can be interpreted as a simple inverse to the estimated aligned signature multiplied by semblance across the array. An application to a zero‐offset VSP acquired with a dynamite source shows the effectiveness of the operator in attaining the two conflicting goals of adaptively spiking the effective source signature and minimizing the noise. Signature design for seismic surveys could benefit from observing that the optimum deconvolution operator gives a flat signal spectrum if and only if the seismic source has the same amplitude spectrum as the noise.


Geophysics ◽  
2021 ◽  
pp. 1-56
Author(s):  
Flavio Poletto ◽  
Alex Goertz ◽  
Cinzia Bellezza ◽  
Endre Vange Bergfjord ◽  
Piero Corubolo ◽  
...  

Seismic while drilling (SWD) by drill-bit source has been successfully used in the past decades and is proven using variable configurations in onshore applications. The method creates a reverse vertical seismic profile (RVSP) dataset from surface sensors deployed as arrays in the proximity of the monitored wells. The typical application makes use of rig-pilot reference (pilot) sensors at the top of the drill-string and also downhole. This approach provides while-drilling checkshots as well as multioffset RVSP for 2-D and 3-D imaging around the well and prediction ahead of the bit. For logistical (sensor deployment) and cost (rig time related to technical installation) reasons the conventional drill-bit SWD application is typically much easier onshore than offshore. We present a novel approach that uses a network of passive-monitoring sea bottom nodes pre-deployed for microseismic monitoring to simultaneously and effectively record offshore SWD data. We study the results of a pilot test where we passively monitored the drilling of an appraisal well at the Wisting discovery in the Barents Sea with an ocean-bottom cable deployed temporarily around the drilling rig. The continuous passive recording of vibration signals emitted during the drilling of the well provides the SWD data set, which is treated as a reverse vertical seismic profile. The study is performed without rig-pilot signal. The results are compared with legacy data and demonstrate the effectiveness of the approach and point to future applications for real-time monitoring of the drilling progress, both in terms of geosteering the drill bit and predicting formation properties ahead of the bit by reflection imaging.


2016 ◽  
Vol 4 (3) ◽  
pp. SJ55-SJ65 ◽  
Author(s):  
Pascal Edme ◽  
David F. Halliday

We have introduced a workflow that allows subsurface imaging using upcoming body-wave arrivals extracted from ambient-noise land seismic data. Rather than using the conventional seismic interferometry approach based on correlation, we have developed a deconvolution technique to extract the earth response from the observed periodicity in the seismic traces. The technique consists of iteratively applying a gapped spiking deconvolution, providing multiple-free images with higher resolution than conventional correlation. We have validated the workflow for zero-offset traces with simple synthetic data and real data recorded during a small point-receiver land seismic survey.


Geophysics ◽  
2003 ◽  
Vol 68 (4) ◽  
pp. 1132-1139 ◽  
Author(s):  
Jérôme H. Le Rousseau ◽  
Henry Calandra ◽  
Maarten V. de Hoop

We illustrate the performance of the generalized screen propagator on real seismic data for 3D zero‐offset and prestack depth imaging. We use TotalFinaElf's L7D data set from the North Sea, a 3D marine seismic survey that contained limited azimuthal coverage. The subsurface shows significant tectonic deformation, including an intrusive salt body in sedimentary sequences. A transformation to common azimuth is applied prior to the 3D prestack depth imaging procedure. We compare the performance of the generalized screen propagator with that of a hybrid phase shift plus interpolation (PSPI)/split‐step Fourier method. Three‐dimensional prestack results confirm the generalized screen method handles multipathing more accurately. Comparisons are also made with Kirchhoff migration results. The results differ mainly in the fine‐scale irregularities of the image and not in the wavefront set of the image. Using synthetic models of similar structure (the SEG/EAGE salt model), we further illustrate the importance of multipathing and multiple scattering. Overall, our results show that our wave‐equation approach produces better images than the Kirchhoff approach to prestack depth migration; we attribute this mainly to a more complete handling of wave diffraction in the generalized screen expansion, which becomes important in strongly heterogeneous and irregular velocity models such as the ones containing salt bodies.


Geophysics ◽  
1993 ◽  
Vol 58 (11) ◽  
pp. 1634-1645 ◽  
Author(s):  
Børge Arntsen ◽  
Bjørn Ursin

The classical one‐dimensional (1-D) inverse problem consists of estimating reflection coefficients from surface seismic data using the 1-D wave equation. Several authors have found stable solutions to this problem using least‐squares model‐fitting methods. We show that the application of these plane‐wave solutions to seismic data generated with a point source can lead to errors in estimating reflection coefficients. This difficulty is avoided by using a least‐squares model fitting scheme describing vertically traveling waves originating from a point source. It is shown that this method is roughly equivalent to deterministic deconvolution with built‐in multiple removal and compensation for spherical spreading. A true zero‐offset field data set from a specially designed seismic experiment is then used as input to estimate reflection coefficients. Stacking velocities from a conventional seismic survey were used to estimate spherical spreading. The resulting reflection coefficients are shown to correlate well with an available well log.


2020 ◽  
Vol 223 (3) ◽  
pp. 1672-1687
Author(s):  
Philipp Prasse ◽  
James Wookey ◽  
J-Michael Kendall ◽  
Daniel Roberts ◽  
Martin Dutko

SUMMARY We present unambiguous evidence that the Mahogany salt body, located in the Northern part of the Gulf of Mexico, is seismically anisotropic. Evidence of anisotropy comes from shear wave splitting data obtained from a vertical seismic profile VSP. The data set consists of 48 vertically aligned receivers in a borehole drilled through the salt body. Splitting analysis is performed on shear wave phases that are converted from compressional waves at the top and bottom of the salt body. The phase converted at the top of the salt layer shows a clear signature of seismic anisotropy, while the phase at the base of the salt layer shows negligible splitting. We investigate the possibility of rock salt halite LPO as a cause of the observed anisotropy. A finite element geomechanical salt deformation model of the Mahogany salt body is developed, where deformation history is used as an input to the texture plasticity simulation program VPSC. Assuming a halite salt body, a full elasticity model is then calculated and used to create a synthetic VSP splitting data set. The comparison between the synthetic and real VSP data set shows that LPO of rock salt can explain the observed anisotropy remarkably well. This is the strongest evidence to date of seismic anisotropy in a deforming salt structure. Furthermore, for the first time, we are able to demonstrate clear evidence that deforming halite is the most likely cause of this anisotropy, combining data set analysis and synthetic full wave form modelling based on calculated rock salt elasticities. Neglecting anisotropy in seismic processing in salt settings could lead to potential imaging errors, for example the deformation models show an averaged delta parameter of δ = –0.06, which would lead in a zero offset reflection setting to a depth mismatch of 6.2 per cent. Our work also show how observations of salt anisotropy can be used to probe characteristics of salt deformation.


Geophysics ◽  
2002 ◽  
Vol 67 (4) ◽  
pp. 1028-1037 ◽  
Author(s):  
R. James Brown ◽  
Robert R. Stewart ◽  
Don C. Lawton

This paper proposes a multicomponent acquisition and preprocessing polarity standard that will apply generally to the three Cartesian geophone components and the hydrophone or microphone components of a 2‐D or 3‐D multicomponent survey on land, at the sea bottom, acquired as a vertical seismic profile, vertical‐cable, or marine streamer survey. We use a four‐component ocean‐bottom data set for purposes of illustration and example. A primary objective is a consistent system of polarity specifications to facilitate consistent horizon correlation among multicomponent data sets and enable determination of correct reflectivity polarity. The basis of this standard is the current SEG polarity standard, first enunciated as a field‐recording standard for vertical geophone data and hydrophone streamer data. It is founded on a right‐handed coordinate system: z positive downward; x positive in the forward line direction in a 2‐D survey, or a specified direction in a 3‐D survey, usually that of the receiver‐cable lines; and y positive in the direction 90° clockwise from x. The polarities of these axes determine the polarity of ground motion in any component direction (e.g., downward ground motion recording as positive values on the vertical‐geophone trace). According also to this SEG standard, a pressure decrease is to be recorded as positive output on the hydrophone trace. We also recommend a cyclic indexing convention, [W, X, Y, Z] or [0, 1, 2, 3], to denote hydrophone or microphone (pressure), inline (radial) geophone, crossline (transverse) geophone, and vertical geophone, respectively. We distinguish among three kinds of polarity standard: acquisition, preprocessing, and final‐display standards. The acquisition standard (summarized in the preceding paragraph) relates instrument output solely to sense of ground motion (geophones) and of pressure change (hydrophones). Polarity considerations beyond this [involving, e.g., source type, wave type (P or S), direction of arrival, anisotropy, tap‐test adjustments, etc.] fall under preprocessing polarity standards. We largely defer any consideration of a display standard.


2021 ◽  
Author(s):  
Olaf Hellwig ◽  
Stefan Buske

<p>The polymetallic, hydrothermal deposit of the Freiberg mining district in the southeastern part of Germany is characterised by ore veins that are framed by Proterozoic orthogneiss. The ore veins consist mainly of quarz, sulfides, carbonates, barite and flourite, which are associated with silver, lead and tin. Today the Freiberg University of Mining and Technology is operating the shafts Reiche Zeche and Alte Elisabeth for research and teaching purposes with altogether 14 km of accessible underground galleries. The mine together with the most prominent geological structures of the central mining district are included in a 3D digital model, which is used in this study to study seismic acquisition geometries that can help to image the shallow as well as the deeper parts of the ore-bearing veins. These veins with dip angles between 40° and 85° are represented by triangulated surfaces in the digital geological model. In order to import these surfaces into our seismic finite-difference simulation code, they have to be converted into bodies with a certain thickness and specific elastic properties in a first step. In a second step, these bodies with their properties have to be discretized on a hexahedral finite-difference grid with dimensions of 1000 m by 1000 m in the horizontal direction and 500 m in the vertical direction. Sources and receiver lines are placed on the surface along roads near the mine. A Ricker wavelet with a central frequency of 50 Hz is used as the source signature at all excitation points. Beside the surface receivers, additional receivers are situated in accessible galleries of the mine at three different depth levels of 100 m, 150 m and 220 m below the surface. Since previous mining activities followed primarily the ore veins, there are only few pilot-headings that cut through longer gneiss sections. Only these positions surrounded by gneiss are suitable for imaging the ore veins. Based on this geometry, a synthetic seismic data set is generated with our explicit finite-difference time-stepping scheme, which solves the acoustic wave equation with second order accurate finite-difference operators in space and time. The scheme is parallelised using a decomposition of the spatial finite-difference grid into subdomains and Message Passing Interface for the exchange of the wavefields between neighbouring subdomains. The resulting synthetic seismic shot gathers are used as input for Kirchhoff prestack depth migration as well as Fresnel volume migration in order to image the ore veins. Only a top mute to remove the direct waves and a time-dependent gain to correct the amplitude decay due to the geometrical spreading are applied to the data before the migration. The combination of surface and in-mine acquisition helps to improve the image of the deeper parts of the dipping ore veins. Considering the limitations for placing receivers in the mine, Fresnel volume migration as a focusing version of Kirchhoff prestack depth migration helps to avoid migration artefacts caused by this sparse and limited acquisition geometry.</p>


Author(s):  
A. Ogbamikhumi ◽  
T. Tralagba ◽  
E. E. Osagiede

Field ‘K’ is a mature field in the coastal swamp onshore Niger delta, which has been producing since 1960. As a huge producing field with some potential for further sustainable production, field monitoring is therefore important in the identification of areas of unproduced hydrocarbon. This can be achieved by comparing production data with the corresponding changes in acoustic impedance observed in the maps generated from base survey (initial 3D seismic) and monitor seismic survey (4D seismic) across the field. This will enable the 4D seismic data set to be used for mapping reservoir details such as advancing water front and un-swept zones. The availability of good quality onshore time-lapse seismic data for Field ‘K’ acquired in 1987 and 2002 provided the opportunity to evaluate the effect of changes in reservoir fluid saturations on time-lapse amplitudes. Rock physics modelling and fluid substitution studies on well logs were carried out, and acoustic impedance change in the reservoir was estimated to be in the range of 0.25% to about 8%. Changes in reservoir fluid saturations were confirmed with time-lapse amplitudes within the crest area of the reservoir structure where reservoir porosity is 0.25%. In this paper, we demonstrated the use of repeat Seismic to delineate swept zones and areas hit with water override in a producing onshore reservoir.


Sign in / Sign up

Export Citation Format

Share Document