Coda wave analysis for distinguishing attenuation due to isotropic scattering from attenuation due to absorption

When the quality factorQ is taken into account in attenuation studies, it is necessary to know the relative losses of wave energy due to scattering and to anelastic absorption. The coda is the most important phenomenon now known which is related to elastic scattering of seismic waves. Utilizing coda, this study presents relationships which give theQ factors of the medium around the recording station and discriminate between attenuations arising from elastic scattering (under the assumption of isotropic scattering) and those arising from anelastic absorption. This work proposes a technique for separately determining the attenuation due to isotropic scattering and that due to absorption from the observed envelope of coda waves.     

Spatial variation of crustal codaQ in California

Coda wave data from California microearthquakes were studied in order to delineate regional fluctuations of apparent crustal attenuation in the band 1.5 to 24 Hz. Apparent attenuation was estimated using a single back scattering model of coda waves. The coda wave data were restricted to 30 s following the origin time; this insures that crustal effects dominate the results as the backscattered shear waves thought to form the coda would not have had time to penetrate much deeper. Results indicate a strong variation in apparent crustal attenuation at high frequencies between the Franciscan and Salinian regions of central California and the Long Valley area of the Sierra Nevada. Although the codaQ measurements coincide at 1.5 Hz (Q _c =100), at 24 Hz there is a factor of four difference between the measurements made in Franciscan (Q _c =525) and Long Valley (Q _c =2100) with the Salinian midway between (Q _c =900). These are extremely large variations compared to measures of seismic velocities of comparable resolution, demonstrating the exceptional sensitivity of the high frequency codaQ measurement to regional geology. In addition, the frequency trend of the results is opposite to that seen in a compilation of codaQ measurements made worldwide by other authors which tend to converge at high and diverge at low frequencies, however, the worldwide results generally were obtained without limiting the coda lengths and probably reflect upper mantle rather than crustal properties. Our results match those expected due to scattering in random media represented by Von Karman autocorrelation functions of orders 1/2 to 1/3. The Von Karman medium of order 1/3 corresponding to the Franciscan codaQ measurement contains greater amounts of high wavenumber fluctuations. This indicates relatively large medium fluctuations with wavelengths on the order of 100 m in the highly deformed crust associated with the Franciscan, however, the influence of scattering on the codaQ measurement is currently a matter of controversy.     

Near‐field seismic effects in a homogeneous medium and their removal in vertical seismic profile attenuation estimates

To better understand (and correct for) the factors affecting the estimation of attenuation (Q), we simulate subsurface wave propagation with the Weyl/Sommerfeld integral. The complete spherical wavefield emanating from a P‐wave point source surrounded by a homogeneous, isotropic and attenuative medium is thus computed. In a resulting synthetic vertical seismic profile, we observe near‐field and far‐field responses and a 90° phase rotation between them. Depth dependence of the magnitude spectra in these two depth regions is distinctly different. The logarithm of the magnitude spectra shows a linear dependence on frequency in the far‐field but not in those depth regions where the near‐field becomes significant. Near‐field effects are one possible explanation for large positive and even negative Q‐factors in the shallow section that may be estimated from real vertical seismic profile data when applying the spectral ratio method. We outline a near‐field compensation technique that can reduce errors in the resultant Q estimates.     

Multiwell study of seismic attenuation at the CO2CRC Otway project geosequestration site: Comparison of amplitude decay,centroid frequency shift and 1D waveform inversion methods

At the CO2CRC Otway geosequestration site, the abundance of borehole seismic and logging data provides a unique opportunity to compare techniques of Q (measure of attenuation) estimation and validate their reliability. Specifically, we test conventional time-domain amplitude decay and spectral-domain centroid frequency shift methods versus the 1D waveform inversion constrained by well logs on a set of zero-offset vertical seismic profiles. The amplitude decay and centroid frequency shift methods of Q estimation assume that a seismic pulse propagates in a homogeneous medium and ignore the interference of the propagating wave with short-period multiples. The waveform inversion explicitly models multiple scattering and interference on a stack of thin layers using high-resolution data from sonic and density logs. This allows for stable Q estimation in small depth windows (in this study, 150 m), and separation of the frequency-dependent layer-induced scattering from intrinsic absorption. Besides, the inversion takes into account band-limited nature of seismic data, and thus, it is less dependent on the operating frequency bandwidth than on the other methods. However, all considered methods of Q estimation are unreliable in the intervals where subsurface significantly deviates from 1D geometry. At the Otway site, the attenuation estimates are distorted by sub-vertical faults close to the boreholes. Analysis of repeated vertical seismic profiles reveals that 15 kt injection of the CO₂-rich fluid into a thin saline aquifer at 1.5 km depth does not induce detectable absorption of P-waves at generated frequencies 5–150 Hz, most likely because the CO₂ plume in the monitoring well is thin, <15 m. At the Otway research site, strong attenuation Q ≈ 30–50 is observed only in shaly formations (Skull Creek Mudstone, Belfast Mudstone). Layer-induced scattering attenuation is negligible except for a few intervals, namely 500–650 m from the surface, and near the injection interval, at around 1400–1550 m, where Q_scat ≈ 50–65.     

VTI介质中非线性对地震波频率和频宽变化的约束

采用改进型FCT有限差分法分别模拟各向同性、具有垂直对称轴的横向各向同性(VTI)以及非线性VTI介质中地震波场,随后采用小波分析探讨VTI介质中非线性对地震波频率和频宽变化的约束.数值模拟及分析结果表明：介质各向异性会引发地震记录中各波型的能量随炮检距的变化,而介质各向异性与非线性的耦合作用则会进一步加剧这种变化.在遵守能量守恒原理的前提下,介质的非线性会引发地震波场中各波型间能量的再分配.同时还会引发地震波速度的波动,且随着入射角的增大,波动幅度变小.这些现象在频率域内的响应表现为主频和频宽的降低或增大、地震记录中与各波型相对应的瞬时主频和频宽峰值出现时刻的偏移以及它们随炮检距变化而表现出的幅值波动.在本文给定的弱非线性介质中（非线性炽较薄且远离震源）,在震源能量较低的情况下,介质的非线性效应会造成地震波瞬时主频和频宽相对于线性介质情况分别为1/10～1/4和1/10～1/3左右的变化量.此外,随着入射波能量的下降,介质非线性所造成的地震波瞬时特征（振幅、主频、频宽）的变化量减小,但随着地震波在穿越非线性介质时路径的加长,非线性作用的累积效应使得波形畸变和频率分布特征的变化（尤其是频宽）等现象得以维持.     

Attenuation models of the earth

Free oscillation and body wave data are used to construct average Q models for the earth. The data set includes fundamental and overtone observations of the radial, spheroidal and toroidal modes, ScS observations and amplitudes of body waves as a function of distance. The preferred model includes a low-Q zone at both the top and the bottom of the mantle. In these regions the seismic velocities are likely to be frequency dependent in the “seismic” band. Absorption in the mantle is predominantly due to losses in shear. Compressional absorption may be important in the inner core.A grain-boundary relaxation model is proposed that explains the dominance of shear over compressional dissipation, the roughly frequency independent average values for Q and the variation of Q with depth. In the high-Q regions, the lithosphere and the midmantle (200–2000 km), Q is predicted to be frequency dependent. However, the low-Q regions of the earth, where Q is roughly frequency independent, dominate the observations of attenuation.     

Seismic inversion with generalized Radon transform based on local second-order approximation of scattered field in acoustic media

Sound velocity inversion problem based on scattering theory is formulated in terms of a nonlinear integral equation associated with scattered field. Because of its nonlinearity, in practice, linearization algorisms (Born/single scattering approximation) are widely used to obtain an approximate inversion solution. However, the linearized strategy is not congruent with seismic wave propagation mechanics in strong perturbation (heterogeneous) medium. In order to partially dispense with the weak perturbation assumption of the Born approximation, we present a new approach from the following two steps: firstly, to handle the forward scattering by taking into account the second-order Born approximation, which is related to generalized Radon transform (GRT) about quadratic scattering potential; then to derive a nonlinear quadratic inversion formula by resorting to inverse GRT. In our formulation, there is a significant quadratic term regarding scattering potential, and it can provide an amplitude correction for inversion results beyond standard linear inversion. The numerical experiments demonstrate that the linear single scattering inversion is only good in amplitude for relative velocity perturbation ( \( \delta_{c}/c_{0} \) ) of background media up to 10 %, and its inversion errors are unacceptable for the perturbation beyond 10 %. In contrast, the quadratic inversion can give more accurate amplitude-preserved recovery for the perturbation up to 40 %. Our inversion scheme is able to manage double scattering effects by estimating a transmission factor from an integral over a small area, and therefore, only a small portion of computational time is added to the original linear migration/inversion process.     

Estimation of frequency dependent coda wave attenuation structure at the vicinity of Cairo Metropolitan Area

Estimation of seismic wave attenuation in the shallow crust in terms of coda wave Q structure previously investigated in the vicinity of Cairo Metropolitan Area was improved using seismograms of local earthquakes recorded by the Egyptian National Seismic Network. The seismic wave attenuation was measured from the time decay of coda wave amplitudes on narrow bandpass filtered seismograms based on the single scattering theory. The frequency bands of interest are from 1.5 to 18 Hz. In general, the values obtained for various events recorded at El-Fayoum and Wadi Hagul stations are very similar for all frequency bands. A regional attenuation law Q _c = 85.66 f ^0.79 was obtained.     

CodaQ as a combination of scattering and intrinsic attenuation: Numerical simulations with the boundary integral method

Numerical modelling ofSH wave seismograms in media whose material properties are prescribed by a random distribution of many perfectly elastic cavities and by intrinsic absorption of seismic energy (anelasticity) demonstrates that the main characteristics of the coda waves, namely amplitude decay and duration, are well described by singly scattered waves in anelastic media rather than by multiply scattered waves in either elastic or anelastic media. We use the Boundary Integral scheme developed byBenites et al. (1992) to compute the complete wave field and measure the values of the direct waveQ and coda wavesQ in a wide range of frequencies, determining the spatial decay of the direct wave log-amplitude relation and the temporal decay of the coda envelope, respectively. The effects of both intrinsic absorption and pure scattering on the overall attenuation can be quantified separately by computing theQ values for corresponding models with (anelastic) and without (elastic) absorption. For the models considered in this study, the values of codaQ ^–1 in anelastic media are in good agreement with the sum of the corresponding scatteringQ ^–1 and intrinsicQ ^–1 values, as established by the single-scattering model ofAki andChouet (1975). Also, for the same random model with intrinsic absorption it appears that the singly scattered waves propagate without significant loss of energy as compared with the multiply scattered waves, which are strongly affected by absorption, suggesting its dominant role in the attenuation of coda waves.     

Borehole study of compressional and shear attenuation of basalt flows penetrated by the Brugdan and William wells on the Faroes shelf

We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective Q_S/Q_P ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small.     

Dispersion attenuation structure of the crust around the Weifang-Jiashan segment of the Tanlu fault zone in North China

On the condition that the velocity ratio of compressive wave to shear wave is stable during the propagating process of seismic wave, this paper develops an estimation method for medium quality factor based on single station's travel time difference between direct S and P waves and the first period signal of direct P wave. 8 774 high SNR wave data altogether recorded by 75 stations are analyzed. The results show that: (1) under the normal regional stress field, the quality factors Q_mSP in the stable tectonic area are higher than that in the active tectonic area around the Weifang-Jiashan section and its surrounding area in the Tanlu fault zone; (2) in the Juxian-Tancheng section the seismic wave attenuation is relatively quick, and the media is relatively broken, suggesting no tectonic stress accumulation; (3) the Xinyi-Sihong section is currently locking and in accumulating elastic strain energy stage, which has the deep environment and conditions of strong earthquake generation similar to those of the 1668 M8.5 Tancheng strong earthquake.     

Estimating permeability from field measurements of seismic attenuation in fractured chalk

Broadband (100–4000 Hz) cross‐hole seismic data have been acquired at a borehole test site where extensive hydrological investigations have previously been performed, including in situ estimates of permeability. The rock type is homogeneous chalk and fractures and bedding planes have been identified from well logs. High values of seismic attenuation, Q= 22 ≤ 27 ≤ 33, were observed over a 10 m depth interval where fracture permeability values of 20–50 darcy had been recorded. An attempt has been made to separate the attenuation due to scattering and intrinsic mechanisms. The estimated values of intrinsic attenuation, Q= 31 ≤ 43 ≤ 71, have been reproduced using a number of current theories of seismic‐wave propagation and fluid‐flow‐induced seismic attenuation in cracked and fractured media. A model that considers wavelength‐scale pressure gradients is the preferred attenuation mechanism. Model parameters were obtained from the hydro‐geological and seismic data. However, we conclude that it is not possible to use seismic Q to measure rock permeability remotely, principally because of the inherent uncertainties arising from model parameterisations.     

Effect of anelastic and scattering structures of the lithosphere on the shape of local earthquake coda

A simple model of single acoustic scattering is used to study the dependence of the shape of local earthquake coda on the anelastic and scattering structures of the lithosphere. The model is applied to the coda of earthquakes located near Stone Canyon, central California, and provides an explanation for the features observed in the data, which include an interesting temporal variation in the coda shape. A surficial layer with aQ of 50 and thickness of 10 or 25 km underlain by a zone with aQ of 1000 extending to the bottom of the lithosphere, together with a scattering scale length,a, that varies with depthz according to the relationa=0.3 exp[-(z/45)²] are found to constitute the simplest structure of the medium compatible with the coda data and with body and surface wave attenuation data. The profile of heterogeneity sizes implies that the scattering strength increases strongly with depth, a constraint required by the necessity to boost the energy of the later coda without forcing the intrinsicQ to be excessively high in the uppermost mantle. This constraint is viewed as an artifact of the single scattering model which overstimates the scattering coefficient due to the neglect of multiple scattering. The observed temporal variation of the signal is difficult to explain by a simple change of the intrinsicQ at some depth. Rather, it is suggested that the scattering properties at depth changed with time through a variation of the fractional rms velocity fluctuation on the order of one percent.     

Incorporation of the Effects of Absorption into Kirchhoff Migration

—The effects of absorption are incorporated into Kirchhoff migration. The aim is to reconstruct the structures with true-amplitude seismic data and to increase the resolution of migrated data. A complex wave velocity is introduced into the solution of the Helmholtz equation, the starting point of Kirchhoff migration. This leads to an additional filter, the antidissipation operator, which is convolved with the wave field. The general structure of Kirchhoff migration remains unchanged. The effects of the antidissipation operator are illustrated on synthetic data. The new operator is valid for complex media with varying velocity and varying quality factor Q. Moreover there is no limitation to constant-Q, frequency-dependent Q can also be handled. The success of anelastic migration depends on how well the Q macro model is known.     

Q c of three component seismograms of volcanic microearthquakes at Campi Flegrei volcanic area — Southern Italy

Digital recordings of three component microearthquake codas from shallow seismic events in the volcanic region of Campi Flegrei — Southern Italy — were used with an automatic technique to calculate the attenuation factorQ _c (codaQ) in the hypothesis of singleS toS backscattering.Results show the same value ofQ for each of the three components. This result is interpreted as due to isotropicS wave radiation pattern.A check of the coda method was performed using a single station method based on simple assumptions on the direct SH wave spectrum. Single stationQ was averaged over the stations and over the earthquakes. Results show that the two methods lead to comparable results.A frequency dependence quite different from that evaluated in active tectonic regions was found for coda attenuation, comparable to other volcanic areas throughout the world. This is interpreted as due to the presence of magma that affects anelasticity and scattering.     

Unified approach to amplitude attenuation and coda excitation in the randomly inhomogeneous lithosphere

A unified model is proposed for explaining the frequency dependent amplitude attenuation and the coda wave excitation on the basis of the single scattering process in the randomly inhomogeneous lithosphere. Adopting Birch's law and a direct proportion between density and wave velocity, we statistically describe the inhomogeneous medium by one random function characterized by the von Karman autocorrelation function. We calculate the amplitude attenuation from the solid angle integral of scattered wave energy on the basis of the Born approxiimation after subtracting the travel-time fluctuation effect caused by slowly varying velocity inhomogeneities. This subtraction is equivalent to neglect energy loss by scattering within a cone around the forward direction. The random inhomogeneity of the von Karman autocorrelation function of order 0.35 with the mean square fractional fluctuation of 7.2×10^–3 1.3×10^–2 and the correlation distance of 2.15.1 km well explains observed backward scattering coefficientg and the ratioQ _P ^–1 /Q _S ^–1 , and observed and partially conjecturedQ _S ^–1 for frequencies between 0.5 Hz and 30 Hz.     

Changes in the early part of the seismic coda due to localized scatterers: The estimation ofQ in a stope environment

A single scattering model was used to analyse the temporary changes in the mean density of scattered waves in a discrete random medium. The model of the mean energy density, originally proposed bySato (1977) for spherical radiation and isotropic scattering, has been modified and applied to a medium in which the scatterers are confined to a specified volume. The time variation of the early part of the mean energy density function for the different source durations was investigated. The dominant effect on the theoretical mean energy density is caused by the specified volume containing scatterers. The duration of the source pulse influences the early part of the coda fort/t ₀<1.2, wheret is the lapse time measured from the source origin time, andt ₀is arrival time of the body wave.The analysis of the coda signal of micro-events occurring immediately in front of the face enables us to estimate the size of the fracture zone induced by the stope. The model of the mean energy density of coda for a medium containing scatterers close to the seismic source was used to analyse a large number of events recorded close to an advancing mine face in a deep level gold mine in South Africa. The coda decay rate has two trends: the first, with a steep decay of coda, is produced by a larger deviation of rock parameters and/or larger size of the scatterers; the second trend, which decays more slowly, has the corresponding mean-free path ranging from 20 m to 200 m. The analysis indicates that the rock mass about 15–20 m from the stope contains a large proportion of fractured and blocked rock, which is the source of scattering. The scattering of theS-wave was much stronger and more stable, with the mean-free path varying from 11 m to 45 m. This is due to the shorter wavelength of theS wave in comparison with theP wave. The quality factor for theP coda wave varies from 30 to 100 in the fracture zone of stope and outside this zone it has a value of 300. The quality factor of theS wave varies from 20 to 78 in the equivalent volume. For rock surrounding the stope the ratioQ _sp ^–1 /Q _ss ^–1 varied from 0.31 to 0.69. This suggests that the radii of scatterers are smaller than 3.5 m.     

大同地震窗P波频散衰减的异常变化特征

通过测量介质品质因子来开展大同窗附近测震台站记录P波衰减频散特征的研究,并分析大同窗近场区和外围区震例前后的变化。对山自皂台记录的2001年以来大同窗地震序列的研究结果表明,近场区3次和外围区6次地震前,介质品质因子显示超出1倍均方差的高值变化特征,震后逐渐恢复到背景起伏状态;与大同窗传统的频次指标相比,介质品质因子的映震效果更好;阳原台的结果表明,仅近场区1次和外围区2次地震前,介质品质因子显示超出1倍均方差的高值变化特征,且有虚报、漏报。所处地质构造位置不同、传播路径介质差异可能是引起上述两个台站附近介质品质因子的地震对应效果明显不同的原因。     

Characteristics of coda wave attenuation in Yunnan area

The characteristic of seismic coda wave attenuation in Yunnan area in 7 frequency-bands range from 1 Hz to 20 Hz was estimated by using the local earthquake's waveform data recorded from 22 Yunnan digital seismic stations.Coda attenuation Q-c1 of each station was firstly calculated by single scattering method. Then, mean free path Le and seismic albedo Bo of each station were calculated, and scattering attenuation Q-1s and intrinsic attenuation Q-1i were separated from total attenuation Q-1t by multiple lapse time window analysis based on the multiple scattering model in uniform random isotropic scattering medium. The attenuating characteristics in Yunnan show that most value of Le are in 10～30 km, with maximal within 2～6 Hz;Bo are about 0.5 at 1～2 Hz, but less than 0.5at other frequency-bands, which means Q-1i is comparable with Q-1s at 1～2 Hz, and after 1～2 Hz, Q-1i is greater than Q-1s and dominates the attenuation process. Q-1c is close to Q-1i at other frequency bands except 1～2 Hz.Results show that Q-1 especially Qs-1 varies spatially, Q-1 in eastern Yunnan zone is a bit higher than in northwestern Yunnan zone;northwestern Yunnan zone higher than southwestern Yunnan zone. Comparing with other results in global, Qs-1 in Yunnan is lower than the global average value among these results, Q-1i is higher than the global average value, and Q-1t lies the middle among these results.