Source signature determination from ministreamer data

Geophysics ◽  
1994 ◽  
Vol 59 (8) ◽  
pp. 1261-1269 ◽  
Author(s):  
Martin Landrø ◽  
Jan Langhammer ◽  
Roger Sollie ◽  
Losse Amundsen ◽  
Eivind Berg
Keyword(s):  
Point Sources ◽  
Inversion Method ◽  
Far Field ◽  
Nonlinear Inversion ◽  
Inversion Algorithm ◽  
Linear Inversion ◽  
Air Bubble ◽  
Air Gun ◽  
Error Energy ◽  
Effective Source

Two methods for estimating the pressure wavefield generated by a marine airp‐gun array are tested. Data have been acquired at a ministreamer located below the source array. Effective source signatures for each air gun are estimated. In the first method a nonlinear inversion algorithm is used, where the forward modeling scheme is based upon a physical modeling of the air bubble generated by each air gun. In the second method a linear inversion method is used, with the assumption that the physics in the problem can be described by the acoustic wave equation with explosive point sources as the driving term. From the estimated effective source signatures, far‐field signatures have been calculated for both methods and compared with measured far‐field signatures. The error energy between the measured and estimated far‐field signatures was approximately 8 percent for both methods.

Source signature determination by inversion

Geophysics ◽  
1992 ◽  
Vol 57 (12) ◽  
pp. 1633-1640 ◽  
Author(s):  
M. Landrø ◽  
R. Sollie
Keyword(s):  
Reflection Coefficient ◽  
Forward Modeling ◽  
Sea Surface ◽  
Inversion Algorithm ◽  
Air Bubble ◽  
Damping Coefficients ◽  
Air Gun ◽  
Marine Air ◽  
Source Array ◽  
Effective Source

A new method for estimating the pressure wavefield generated by a marine air‐gun array is presented. It is assumed that data is acquired at a ministreamer located below the source array. Effective source signatures for each air gun are estimated by an inversion algorithm. The forward modeling scheme used in the inversion algorithm is based upon a physical modeling of the air bubble generated by each air gun. This means that typical inversion parameters are: gun depths, empirical damping coefficients, and reflection coefficient of the sea surface. Variations in streamer depth are also taken into account by the inversion scheme. The algorithm has been successfully tested on examples with unknown streamer positions, gun parameters, reflection coefficient of sea surface, and ministreamer data contaminated with white noise.

2-D and 3-D resistivity image reconstruction using crosshole data

Geophysics ◽  
1992 ◽  
Vol 57 (10) ◽  
pp. 1270-1281 ◽  
Author(s):  
Hiromasa Shima
Keyword(s):  
Field Test ◽  
Numerical Models ◽  
Electrical Potential ◽  
Nonlinear Least Squares ◽  
Three Dimensions ◽  
Inversion Method ◽  
Good Resolution ◽  
Inversion Algorithm ◽  
Linear Inversion ◽  
Resistivity Model

Theoretical changes in the distribution of electrical potential near subsurface resistivity anomalies have been studied using two resistivity models. The results suggest that the greatest response from such anomalies can be observed with buried electrodes, and that the resistivity model of a volume between boreholes can be accurately reconstructed by using crosshole data. The distributive properties of crosshole electrical potential data obtained by the pole‐pole array method have also been examined using the calculated partial derivative of the observed apparent resistivity with respect to a small cell within a given volume. The results show that for optimum two‐dimensional (2-D) and three‐dimensional (3-D) target imaging, in‐line data and crossline data should be combined, and an area outside the zone of exploration should be included in the analysis. In this paper, the 2-D and 3-D resistivity images presented are reconstructed from crosshole data by the combination of two inversion algorithms. The first algorithm uses the alpha center method for forward modeling and reconstructs a resistivity model by a nonlinear least‐squares inversion. Alpha centers express a continuously varying resistivity model, and the distribution of the electrical potential from the model can be calculated quickly. An initial general model is determined by the resistivity backprojection technique (RBPT) prior to the first inversion step. The second process uses finite elements and a linear inversion algorithm to improve the resolution of the resistivity model created by the first step. Simple 2-D and 3-D numerical models are discussed to illustrate the inversion method used in processing. Data from several field studies are also presented to demonstrate the capabilities of using crosshole resistivity exploration techniques. The numerical experiments show that by using the combined reconstruction algorithm, thin conductive layers can be imaged with good resolution for 2-D and 3-D cases. The integration of finite‐element computations is shown to improve the image obtained by the alpha center inversion process for 3-D applications. The first field test uses horizontal galleries to evaluate complex 2-D features of a zinc mine. The second field test illustrates the use of three boreholes at a dam site to investigate base rock features and define the distribution of an altered zone in three dimensions.

Turning‐ray tomography

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1917-1929 ◽  
Author(s):  
Joseph P. Stefani
Keyword(s):  
Velocity Structure ◽  
Initial Point ◽  
Real Data ◽  
Point Sources ◽  
Surface Velocity ◽  
Inversion Algorithm ◽  
Low Velocity ◽  
Linear Inversion ◽  
Near Surface ◽  
Velocity Inversion

Turning‐ray tomography is useful for estimating near‐surface velocity structure in areas where conventional refraction statics techniques fail because of poor data or lack of smooth refractor/velocity structure. This paper explores the accuracy and inherent smoothing of turning‐ray tomography in its capacity to estimate absolute near‐surface velocity and the statics times derived from these velocities, and the fidelity with which wavefields collapse to point diffractors when migrated through these estimated velocities. The method comprises nonlinear iterations of forward ray tracing through triangular cells linear in slowness squared, coupled with the LSQR linear inversion algorithm. It is applied to two synthetic finite‐ difference data sets of types that usually foil conventional refraction statics techniques. These models represent a complex hard‐rock overthrust structure with a low‐velocity zone and pinchouts, and a contemporaneous near‐shore marine trench filled with low‐ velocity unconsolidated deposits exhibiting no seismically apparent internal structure. In both cases velocities are estimated accurately to a depth of one‐ fifth the maximum offset, as are the associated statics times. Of equal importance, the velocities are sufficiently accurate to correctly focus synthetic wavefields back to their initial point sources, so migration/datuming applications can also use these velocities. The method is applied to a real data example from the Timbalier Trench in the Gulf of Mexico, which exhibits the same essential features as the marine trench synthetic model. The Timbalier velocity inversion is geologically reasonable and yields long and short wavelength statics that improve the CMP gathers and stack and that correctly align reflections to known well markers. Turning‐ray tomography estimates near‐surface velocities accurately enough for the three purposes of lithology interpretation, statics calculations, and wavefield focusing for shallow migration and datuming.

Refinements to the linear velocity inversion theory

Geophysics ◽  
1984 ◽  
Vol 49 (2) ◽  
pp. 112-118 ◽  
Author(s):  
Frank G. Hagin ◽  
Jack K. Cohen
Keyword(s):  
Synthetic Data ◽  
Scattered Wave ◽  
Scattering Data ◽  
Inversion Method ◽  
Impedance Change ◽  
Inversion Algorithm ◽  
Linear Inversion ◽  
Velocity Inversion ◽  
Linear Algorithm ◽  
Inversion Theory

The linear inversion method presented by Cohen and Bleistein in 1979 gives seriously degraded results when large reflectors are encountered. Obviously there is an irrecoverable loss of information when such a linear algorithm is applied to a nonlinear world. However, in many cases, excellent results can be achieved by suitable postprocessing of the output of the basic linear inversion algorithm. Although a certain degree of helpful postprocessing can be and has been performed by straightforward consideration of the linearization process, we present here a substantially improved postprocessing algorithm. The basis for these improvements is a more accurate scattering model due to Lahlou et al where, among other things, a WKB analysis of the wave equation led to a much more accurate accounting of the geometric spreading of the scattered wave. These notions plus an effective use of traveltime are used in the new algorithm to improve both the estimate of the reflector locations and the estimate of amplitude (velocity or acoustic impedance) change across the reflectors. The basic idea is to insert this idealized scattering data into the original linear algorithm, and then use the result of this computation as a guide in the interpretation of the numerical output of the algorithm. We demonstrate the result of computer implementation of this algorithm on synthetic data, with and without noise, and verify that the postprocessing algorithm produces dramatically improved reflector locations and speed estimates. Moreover, the new algorithm adds only very modest cost to the basic processing, which is, in turn, very competitive in cost to other multidimensional algorithms.

Air‐gun bubble damping by a screen

Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1765-1772 ◽  
Author(s):  
Jan Langhammer ◽  
Martin Landrø ◽  
James Martin ◽  
Eivind Berg
Keyword(s):  
Near Field ◽  
Field Measurements ◽  
Far Field ◽  
Pressure Amplitude ◽  
Air Bubble ◽  
Air Gun ◽  
Pressure Peak ◽  
Bubble Oscillations ◽  
Measured Pressure ◽  
Optimal Radius

A method for damping unwanted bubble oscillations from a seismic air gun is presented. The method exploits the fact that the primary pressure peak generated by an air gun is produced during the first 5–10 ms after firing. The air bubble is destroyed by mounting a perforated screen with an optimal radius about the gun. Once the primary pressure peak has been generated by the bubble, the bubble is destroyed by the screen, leading to a corresponding decrease in the measured pressure amplitude of the secondary bubble oscillations. Controlled near‐field measurements of 40‐cubic inch and 120‐cubic inch air guns with and without damping screens are used. The primary to bubble ratio improves from 1.4 without a screen to 4.4 with a screen in the near‐field. The corresponding values for estimated far‐field signatures are 1.8 to 9.0 when the signatures are filtered with an out‐128 Hz (72 dB/Oct) DFS V filter.

Nonlinear amplitude-variation-with-offset inversion for Lamé parameters using a direct inversion method

2017 ◽  
Vol 5 (3) ◽  
pp. SL57-SL67 ◽  
Author(s):  
Guangsen Cheng ◽  
Xingyao Yin ◽  
Zhaoyun Zong
Keyword(s):  
Nonlinear Equation ◽  
Real Data ◽  
Seismic Inversion ◽  
Inversion Method ◽  
Nonlinear Inversion ◽  
Reservoir Prediction ◽  
Amplitude Variation ◽  
Linear Inversion ◽  
Amplitude Variation With Offset ◽  
Avo Inversion

Prestack seismic inversion is widely used in fluid indication and reservoir prediction. Compared with linear inversion, nonlinear inversion is more precise and can be applied to high-contrast situations. The inversion results can be affected by the parameters’ sensitivity, so the parameterization of nonlinear equations is very significant. Considering the poor nonlinear amplitude-variation-with-offset (AVO) inversion results of impedance and velocity parameters, we adjust the parameters of the nonlinear equation, avoid the inaccuracy caused by parameters sensitivity and get the ideal nonlinear AVO inversion results of the Lamé parameters. The feasibility and stability of the nonlinear equation based on the Lamé parameters and method are verified by the model and the real data examples. The resolution and the lateral continuity of nonlinear inversion results are better compared with the linear inversion results.

A simultaneous inversion for background velocity and impedance maps

Geophysics ◽  
1990 ◽  
Vol 55 (4) ◽  
pp. 458-469 ◽  
Author(s):  
D. Cao ◽  
W. B. Beydoun ◽  
S. C. Singh ◽  
A. Tarantola
Keyword(s):  
Traditional Approach ◽  
Waveform Inversion ◽  
Synthetic Data ◽  
Inversion Method ◽  
Nonlinear Inversion ◽  
Seismic Velocities ◽  
Inversion Algorithm ◽  
Seismic Reflection Data ◽  
Simultaneous Inversion ◽  
Highly Nonlinear

Full‐waveform inversion of seismic reflection data is highly nonlinear because of the irregular form of the function measuring the misfit between the observed and the synthetic data. Since the nonlinearity results mainly from the parameters describing seismic velocities, an alternative to the full nonlinear inversion is to have an inversion method which remains nonlinear with respect to velocities but linear with respect to impedance contrasts. The traditional approach is to decouple the nonlinear and linear parts by first estimating the background velocity from traveltimes, using either traveltime inversion or velocity analysis, and then estimating impedance contrasts from waveforms, using either waveform inversion or conventional migration. A more sophisticated strategy is to obtain both the subsurface background velocities and impedance contrasts simultaneously by using a single least‐squares norm waveform‐fit criterion. The background velocity that adequately represents the gross features of the medium is parameterized using a sparse grid, whereas the impedance contrasts use a dense grid. For each updated velocity model, the impedance contrasts are computed using a linearized inversion algorithm. For a 1-D velocity background, it is very efficient to perform inversion in the f-k domain by using the WKBJ and Born approximations. The method performs well both with synthetic and field data.

Damping of secondary bubble oscillations for towed air guns with a screen

Geophysics ◽  
1997 ◽  
Vol 62 (2) ◽  
pp. 533-539
Author(s):  
Martin Landrø ◽  
Jan Langhammer ◽  
James Martin
Keyword(s):  
Near Field ◽  
Far Field ◽  
Pressure Amplitude ◽  
Air Bubble ◽  
Air Gun ◽  
Pressure Peak ◽  
Bubble Oscillations ◽  
Measured Pressure ◽  
Optimal Radius ◽  
Boat Speed

A method for damping unwanted bubble oscillations from a horizontally towed seismic air gun is presented. The air bubble is destroyed by a perforated screen mounted at an optimal radius about the gun. Once the primary pressure peak has been generated by the emerging bubble, the bubble continues to expand and is destroyed by the screen, leading to a corresponding decrease in the measured pressure amplitude of the secondary bubble oscillations. For a stationary gun fired first without, and then with, the screen fitted, the primary‐to‐ bubble ratio improves in the near field from 1.7 to 5.2, respectively, at a firing depth of 3 m and from 1.5 to 5.5, respectively, at 5 m depth. The primary‐to‐bubble ratio for a towed air gun in the quasi‐far‐field improves from 2.0 to 11.0 at 4 m depth and from 1.5 to 8.7 at 7 m depth when the screen is fitted. The boat speed was 1.6 knots and the signatures were filtered with an out‐128 Hz (72 dB/Oct) DFS V filter.

scholarly journals Application of an Improved Ultrasound Full-Waveform Inversion in Bone Quantitative Measurement

Symmetry ◽  
2021 ◽  
Vol 13 (2) ◽  
pp. 260
Author(s):  
Meng Suo ◽  
Dong Zhang ◽  
Yan Yang
Keyword(s):  
Nonlinear Equations ◽  
Waveform Inversion ◽  
Optimal Solution ◽  
Inversion Method ◽  
Full Waveform Inversion ◽  
Inversion Algorithm ◽  
Elastic Wave Equation ◽  
Full Waveform ◽  
Ultrasonic Propagation ◽  
Joint Velocity

Inspired by the large number of applications for symmetric nonlinear equations, an improved full waveform inversion algorithm is proposed in this paper in order to quantitatively measure the bone density and realize the early diagnosis of osteoporosis. The isotropic elastic wave equation is used to simulate ultrasonic propagation between bone and soft tissue, and the Gauss–Newton algorithm based on symmetric nonlinear equations is applied to solve the optimal solution in the inversion. In addition, the authors use several strategies including the frequency-grid multiscale method, the envelope inversion and the new joint velocity–density inversion to improve the result of conventional full-waveform inversion method. The effects of various inversion settings are also tested to find a balanced way of keeping good accuracy and high computational efficiency. Numerical inversion experiments showed that the improved full waveform inversion (FWI) method proposed in this paper shows superior inversion results as it can detect small velocity–density changes in bones, and the relative error of the numerical model is within 10%. This method can also avoid interference from small amounts of noise and satisfy the high precision requirements for quantitative ultrasound measurements of bone.

The shape of things to come — Development and testing of a new marine vibrator source

2019 ◽  
Vol 38 (9) ◽  
pp. 680-690 ◽  
Author(s):  
Benoît Teyssandier ◽  
John J. Sallas
Keyword(s):  
Seismic Source ◽  
Test Facility ◽  
Data Sets ◽  
Far Field ◽  
Ocean Bottom ◽  
Air Gun ◽  
Seismic Image ◽  
Field Radiation ◽  
Technical Issues ◽  
To Come

Ten years ago, CGG launched a project to develop a new concept of marine vibrator (MV) technology. We present our work, concluding with the successful acquisition of a seismic image using an ocean-bottom-node 2D survey. The expectation for MV technology is that it could reduce ocean exposure to seismic source sound, enable new acquisition solutions, and improve seismic data quality. After consideration of our objectives in terms of imaging, productivity, acoustic efficiency, and operational risk, we developed two spectrally complementary prototypes to cover the seismic bandwidth. In practice, an array composed of several MV units is needed for images of comparable quality to those produced from air-gun data sets. Because coupling to the water is invariant, MV signals tend to be repeatable. Since far-field pressure is directly proportional to piston volumetric acceleration, the far-field radiation can be well controlled through accurate piston motion control. These features allow us to shape signals to match precisely a desired spectrum while observing equipment constraints. Over the last few years, an intensive validation process was conducted at our dedicated test facility. The MV units were exposed to 2000 hours of in-sea testing with only minor technical issues.

Sign in / Sign up

Export Citation Format

Share Document