An improved method to estimate Q based on the logarithmic spectrum of moving peak points

2019 ◽  
Vol 7 (2) ◽  
pp. T255-T263 ◽  
Author(s):  
Yanli Liu ◽  
Zhenchun Li ◽  
Guoquan Yang ◽  
Qiang Liu

The quality factor ([Formula: see text]) is an important parameter for measuring the attenuation of seismic waves. Reliable [Formula: see text] estimation and stable inverse [Formula: see text] filtering are expected to improve the resolution of seismic data and deep-layer energy. Many methods of estimating [Formula: see text] are based on an individual wavelet. However, it is difficult to extract the individual wavelet precisely from seismic reflection data. To avoid this problem, we have developed a method of directly estimating [Formula: see text] from reflection data. The core of the methodology is selecting the peak-frequency points to linear fit their logarithmic spectrum and time-frequency product. Then, we calculated [Formula: see text] according to the relationship between [Formula: see text] and the optimized slope. First, to get the peak frequency points at different times, we use the generalized S transform to produce the 2D high-precision time-frequency spectrum. According to the seismic wave attenuation mechanism, the logarithmic spectrum attenuates linearly with the product of frequency and time. Thus, the second step of the method is transforming a 2D spectrum into 1D by variable substitution. In the process of transformation, we only selected the peak frequency points to participate in the fitting process, which can reduce the impact of the interference on the spectrum. Third, we obtain the optimized slope by least-squares fitting. To demonstrate the reliability of our method, we applied it to a constant [Formula: see text] model and the real data of a work area. For the real data, we calculated the [Formula: see text] curve of the seismic trace near a well and we get the high-resolution section by using stable inverse [Formula: see text] filtering. The model and real data indicate that our method is effective and reliable for estimating the [Formula: see text] value.

Geophysics ◽  
2000 ◽  
Vol 65 (2) ◽  
pp. 377-389 ◽  
Author(s):  
Paul J. Hatchell

Transmission distortions are observed on prestack seismic data at two locations in the Gulf of Mexico. These distortions produce anomalous amplitude versus offset (AVO) signatures. The locations of the distortion zones are determined using acquisition geometry and ray tracing. No obvious reflection events, such as shallow gas zones, are observed at the predicted locations of the distortion zones. Instead, the distortion zones correlate with buried faults and unconformities. It is postulated that the distortions are produced by velocity changes across buried faults and unconformities. The distortions result from an interference pattern resulting from seismic waves arriving from different sides of the faults. A simple model is developed to explain many of the characteristics of the distortion pattern.


2018 ◽  
Vol 477 (1) ◽  
pp. 413-426 ◽  
Author(s):  
Joshu J. Mountjoy ◽  
Xiaoming Wang ◽  
Susi Woelz ◽  
Sean Fitzsimons ◽  
Jamie D. Howarth ◽  
...  

AbstractLacustrine-tsunami risk from landslides can be significant yet for most locations globally the hazard remains unquantified. Lake Tekapo, in the tectonically active mountain belt of New Zealand's South Island, has been chosen to develop surveying and modelling techniques to assess the hazard from landslide tsunamis. Lake Tekapo is ideal for this study due to the high sedimentation rates, steep surrounds and the proximity to active faulting that indicate a high landslide potential. The shoreline tourist settlement and hydropower infrastructure mean the impact of any tsunami could be significant. In 2016 a survey was carried out to collect high-resolution (1 m grid) EM2040 multibeam bathymetry, high-resolution seismic reflection data (Boomer and chirp) and 6 m long sediment cores. These data reveal a diverse range of sedimentary processes in response to high sediment input and numerous landslides with varied styles of emplacement. For example, a one-off landslide initiated 40 m above the shoreline with debris deposits that have runout onto the lake floor to 100 m water depth contrasts with the Cass River delta on the western shore that has failed multiple times during the lake-basin infilling history. Landslide-generated tsunami scenarios are used to determine the relative hazard at different regions of the lake to guide development of a probabilistic tsunami model.


Solid Earth ◽  
2019 ◽  
Vol 10 (3) ◽  
pp. 851-869 ◽  
Author(s):  
Salomé Mignard ◽  
Thierry Mulder ◽  
Philippe Martinez ◽  
Thierry Garlan

Abstract. The effects of changes in slope gradient on deposition processes and architecture have been investigated in different deep-sea systems both in modern and ancient environments. However, the impact of subtle gradient changes (< 0.3∘) on sedimentary processes along deep-sea fans still needs to be clarified. The Ogooue Fan, located in the northeastern part of the Gulf of Guinea, extends over more than 550 km westwards of the Gabonese shelf and passes through the Cameroon volcanic line. Here, we present the first study of acoustic data (multibeam echosounder and 3.5 kHz, very high-resolution seismic data) and piston cores covering the deep-sea part of this West African system. This study documents the architecture and sedimentary facies distribution along the fan. Detailed mapping of near-seafloor seismic-reflection data reveals the influence of subtle slope gradient changes (< 0.2∘) along the fan morphology. The overall system corresponds to a well-developed deep-sea fan, fed by the Ogooue River sedimentary load, with tributary canyons, distributary channel–levee complexes and lobe elements. However, variations in the slope gradient due to inherited salt-related structures and the presence of several seamounts, including volcanic islands, result in a topographically complex slope profile including several ramps and steps. In particular, turbidity currents derived from the Gabonese shelf deposit cross several interconnected intra-slope basins located on the low gradient segments of the margin (< 0.3∘). On a higher gradient segment of the slope (0.6∘), a large mid-system valley developed connecting an intermediate sedimentary basin to the more distal lobe area. Distribution and thickness of turbidite sands is highly variable along the system. However, turbidite sands are preferentially deposited on the floor of the channel and the most proximal depositional areas. Core description indicates that the upper parts of the turbidity flows, mainly composed of fine-grained sediments, are found in the most distal depocenters.


Geophysics ◽  
2017 ◽  
Vol 82 (1) ◽  
pp. N1-N12 ◽  
Author(s):  
Francisco de S. Oliveira ◽  
Jose J. S. de Figueiredo ◽  
Andrei G. Oliveira ◽  
Jörg Schleicher ◽  
Iury C. S. Araújo

Quality factor estimation and correction are necessary to compensate the seismic energy dissipated during acoustic-/elastic-wave propagation in the earth. In this process, known as [Formula: see text]-filtering in the realm of seismic processing, the main goal is to improve the resolution of the seismic signal, as well as to recover part of the energy dissipated by the anelastic attenuation. We have found a way to improve [Formula: see text]-factor estimation from seismic reflection data. Our methodology is based on the combination of the peak-frequency-shift (PFS) method and the redatuming operator. Our innovation is in the way we correct traveltimes when the medium consists of many layers. In other words, the correction of the traveltime table used in the PFS method is performed using the redatuming operator. This operation, performed iteratively, allows a more accurate estimation of the [Formula: see text] factor layer by layer. Applications to synthetic and real data (Viking Graben) reveal the feasibility of our analysis.


Geophysics ◽  
1985 ◽  
Vol 50 (6) ◽  
pp. 903-923 ◽  
Author(s):  
T. N. Bishop ◽  
K. P. Bube ◽  
R. T. Cutler ◽  
R. T. Langan ◽  
P. L. Love ◽  
...  

Estimation of reflector depth and seismic velocity from seismic reflection data can be formulated as a general inverse problem. The method used to solve this problem is similar to tomographic techniques in medical diagnosis and we refer to it as seismic reflection tomography. Seismic tomography is formulated as an iterative Gauss‐Newton algorithm that produces a velocity‐depth model which minimizes the difference between traveltimes generated by tracing rays through the model and traveltimes measured from the data. The input to the process consists of traveltimes measured from selected events on unstacked seismic data and a first‐guess velocity‐depth model. Usually this first‐guess model has velocities which are laterally constant and is usually based on nearby well information and/or an analysis of the stacked section. The final model generated by the tomographic method yields traveltimes from ray tracing which differ from the measured values in recorded data by approximately 5 ms root‐mean‐square. The indeterminancy of the inversion and the associated nonuniqueness of the output model are both analyzed theoretically and tested numerically. It is found that certain aspects of the velocity field are poorly determined or undetermined. This technique is applied to an example using real data where the presence of permafrost causes a near‐surface lateral change in velocity. The permafrost is successfully imaged in the model output from tomography. In addition, depth estimates at the intersection of two lines differ by a significantly smaller amount than the corresponding estimates derived from conventional processing.


Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. V235-V247 ◽  
Author(s):  
Duan Li ◽  
John Castagna ◽  
Gennady Goloshubin

The frequency-dependent width of the Gaussian window function used in the S-transform may not be ideal for all applications. In particular, in seismic reflection prospecting, the temporal resolution of the resulting S-transform time-frequency spectrum at low frequencies may not be sufficient for certain seismic interpretation purposes. A simple parameterization of the generalized S-transform overcomes the drawback of poor temporal resolution at low frequencies inherent in the S-transform, at the necessary expense of reduced frequency resolution. This is accomplished by replacing the frequency variable in the Gaussian window with a linear function containing two coefficients that control resolution variation with frequency. The linear coefficients can be directly calculated by selecting desired temporal resolution at two frequencies. The resulting transform conserves energy and is readily invertible by an inverse Fourier transform. This modification of the S-transform, when applied to synthetic and real seismic data, exhibits improved temporal resolution relative to the S-transform and improved resolution control as compared with other generalized S-transform window functions.


2013 ◽  
Vol 838-841 ◽  
pp. 1256-1262
Author(s):  
Yun Wen Zheng ◽  
Quan Mei Gong ◽  
Yun Zhuang Zheng

The excavation of foundation pits near existing railway will cause innegligible deformation for the railway line. This paper first conducts numerical simulation for the impact of the excavation of foundation pit for a building close to existing railway lines on the deformation of the railway, and arrives at such conclusions: the excavation of foundation pits has a great impact on the railway lines. The amount of vertical deformation is beyond the requirements of relevant regulations. Accordingly, the foundation needs stabilizing treatment. After the reinforcement of high pressure jet-grouting pile, the maximum vertical deformation decreases obviously to meets the requirements. On this basis, the paper analyzes the reduction amount on different phases of the excavation to work out their respective reinforcement effect. After that, it compares the results with the real data collected during the construction to test the correctness of the model, finds out the major periods during which the vertical deformation occurs as well as the amount of deformation in each period, and points out the worst period for railway deformation during the excavation of foundation pit.


Geophysics ◽  
1986 ◽  
Vol 51 (6) ◽  
pp. 1266-1276 ◽  
Author(s):  
Luc T. Ikelle ◽  
Jean Paul Diet ◽  
Albert Tarantola

In the acoustic approximation, the Earth is described using only density and bulk modulus. Assuming smooth density variations, reflections can be described using a single function—the velocity of compressional waves. If a reference model which is close enough to the actual Earth is known, the problem of estimating the medium velocity from the observed data can be linearized. Using a least‐squares formulation and working in the ω-k domain, the linearized inverse problem for a homogeneous reference medium can be solved by a noniterative algorithm which is economically competitive with prestack migration. Numerical tests with synthetic and real data demonstrate the feasibility and the numerical stability of the method. The numerical results compare well with those obtained by migration of unstacked data, although superior results will only be obtained when the physics of the problem (including elastic versus acoustic effects, three‐dimensional propagation, and accurate source estimation) will realistically be taken into account.


Geophysics ◽  
2003 ◽  
Vol 68 (5) ◽  
pp. 1714-1730 ◽  
Author(s):  
Bertrand Iooss ◽  
David Geraets ◽  
Tapan Mukerji ◽  
Yann Samuelides ◽  
Mustafa Touati ◽  
...  

Understanding the internal heterogeneities of reservoirs is one of the key issues in better recovery and efficient reservoir management. Seismic data are widely used to map subsurface heterogeneities. These heterogeneities can include variations in wave velocity and rock density, which can be used to interpret variations in reservoir properties such as porosity, lithofacies, and fluids. This paper describes a statistical tomography method to infer the spatial statistics of subsurface velocity heterogeneities from seismic data. We consider an acoustic wave propagating in a medium represented as a single macromodel superimposed on statistically stationary random velocity perturbations. While the macromodel is retrieved by classical seismic methods, the picked traveltimes and their fluctuations are used to estimate the covariance function of the spatially varying velocity perturbations. We present a formulation based on ray‐theoretical results and describe two algorithms: one using the prestack traveltimes and the other using the stacking velocities. The methods are tested with synthetic seismic reflection data in an idealized medium (with a Gaussian spatial covariance) and with synthetic transmission data in a more geologically realistic medium. Then, the two algorithms are applied on real data. The estimates of the spatial statistics obtained from inverting the traveltime statistics match reasonably well with the true parameters of the heterogeneous media.


Sign in / Sign up

Export Citation Format

Share Document