A laboratory study of seismic velocity and attenuation anisotropy in near-surface sedimentary rocks

The laboratory ultrasonic pulse‐echo method was used to collect accurate P‐ and S‐wave velocity (±0.3%) and attenuation (±10%) data at differential pressures of 5–50 MPa on water‐saturated core samples of sandstone, limestone and siltstone that were cut parallel and perpendicular to the vertical borehole axis. The results, when expressed in terms of the P‐ and S‐wave velocity and attenuation anisotropy parameters for weakly transversely isotropic media (ɛ, γ, ɛ_Q, γ_Q) show complex variations with pressure and lithology. In general, attenuation anisotropy is stronger and more sensitive to pressure changes than velocity anisotropy, regardless of lithology. Anisotropy is greatest (over 20% for velocity, over 70% for attenuation) in rocks with visible clay/organic matter laminations in hand specimens. Pressure sensitivities are attributed to the opening of microcracks with decreasing pressure. Changes in magnitude of velocity and attenuation anisotropy with effective pressure show similar trends, although they can show different signs (positive or negative values of ɛ, ɛ_Q, γ, γ_Q). We conclude that attenuation anisotropy in particular could prove useful to seismic monitoring of reservoir pressure changes if frequency‐dependent effects can be quantified and modelled.     

Attenuation of P- and S-waves in the Chamoli Region, Himalaya, India

The attenuation properties of the crust in the Chamoli region of Himalaya have been examined by estimating the frequency-dependent relationships of quality factors for P waves (Q_α) and for S waves (Q_β) in the frequency range 1.5–24 Hz. The extended coda normalization method has been applied on the waveforms of 25 aftershocks of the 1999 Chamoli earthquake (M 6.4) recorded at five stations. The average value of Q_α is found to be varied from 68 at 1.5 Hz to 588 at 24 Hz while it varies from 126 at 1.5 Hz to 868 at 24 Hz for Q_β. The estimated frequency-dependent relations for quality factors are Q_α = (44 ± 1)f^(0.82±.04) and Q_β = (87 ± 3)f^(0.71±.03). The rate of increase of Q(f) for P and S waves in the Chamoli region is comparable with the other regions of the world. The ratio Q_β/Q_α is greater than one in the region which along with the frequency dependence of quality factors indicates that scattering is an important factor contributing to the attenuation of body waves in the region. A comparison of attenuation relation for S wave estimated here (Q_β = 87f^0.71) with that of coda waves (Q_c = 30f^1.21) obtained by Mandal et al. (2001) for the same region shows that Q_c > Q_β for higher frequencies (>8 Hz) in the region. This indicates a possible high frequency coda enrichment which suggests that the scattering attenuation significantly influences the attenuation of S waves at frequencies >8 Hz. This observation may be further investigated using multiple scattering models. The attenuation relations for quality factors obtained here may be used for the estimation of source parameters and near-source simulation of earthquake ground motion of the earthquakes, which in turn are required for the assessment of seismic hazard in the region.     

VERTICAL SEISMIC PROFILING IN COAL*

The intellection of seismic wave propagation in coal measures demands direct observation of the wavefield progression. Two vertical seismic profiles with high spatial and temporal sampling, were recently recorded in the Sydney Basin coalfields as part of an experimental coal seismic program. Static corrections and interval velocities were obtained by an automated system to determine first kicks and pulse rise times. Upgoing and downgoing waves were separated in the f—k-plane using a novel technique of contour slice filtering. The isolated upgoing waves clearly display reflections from the major coal seams within the stratigraphic sequence. The downgoing wave spectra were subjected to attenuation analysis. The deduced specific quality factor Q for Permian coal measure rocks lies in the range 20–70. Similar estimates were obtained in the time domain from measurements of pulse broadening. Synthetic VSP seismograms, computed using an exact recursive formulation, are an indispensable aid to interpretation. They illustrate the filtering effects of coal seams and sequences, and the effects of the contribution of internal and free-surface multiple reflections in the recorded wavetrains.     

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.     

Seismic attenuation in Faroe Islands basalts

We analysed vertical seismic profiling (VSP) data from two boreholes at Glyvursnes and Vestmanna on the island of Streymoy, Faroe Islands, to determine the magnitude and causes of seismic attenuation in sequences of basalt flows. The work is part of SeiFaBa, a major project integrating data from vertical and offset VSP, surface seismic surveys, core samples and wireline log data from the two boreholes. Values of effective seismic quality factor (Q) obtained at Glyvursnes and Vestmanna are sufficiently low to significantly degrade the quality of a surface reflection seismic image. This observation is consistent with results from other VSP experiments in the North Atlantic region. We demonstrate that the most likely cause of the low values of effective Q at Glyvursnes and Vestmanna is a combination of 1D scattering and intrinsic attenuation due to seismic wave‐induced fluid flow within pores and micro‐cracks. Tests involving 3D elastic wave numerical modelling with a hypothetical basalt model based on field observations, indicate that little scattering attenuation is caused by lateral variations in basalt structure.     

Attenuation tomography: An application to gas-hydrate and free-gas detection

We estimate the quality factor (Q) from seismic reflections by using a tomographic inversion algorithm based on the frequency‐shift method. The algorithm is verified with a synthetic case and is applied to offshore data, acquired at western Svalbard, to detect the presence of bottom‐simulating reflectors (BSR) and gas hydrates. An array of 20 ocean‐bottom seismographs has been used. The combined use of traveltime and attenuation tomography provides a 3D velocity–Q cube, which can be used to map the spatial distribution of the gas‐hydrate concentration and free‐gas saturation. In general, high P‐wave velocity and quality factor indicate the presence of solid hydrates and low P‐wave velocity and quality factor correspond to free‐gas bearing sediments. The Q‐values vary between 200 and 25, with higher values (150–200) above the BSR and lower values below the BSR (25–40). These results seem to confirm that hydrates cement the grains, and attenuation decreases with increasing hydrate concentration.     

An attenuation study using earthquakes from the 1997 Umbria-Marche sequence

We studied spatial and temporal characteristics of seismic attenuation inCentral Italy using S- and coda- waves recorded by the MarchesanSeismograph Network from earthquakes located in the epicentral area ofthe 1997 Umbria-Marche sequence. The amplitude decay of the S waveswith distance was defined calculating empirical attenuation functions at 15frequencies between 1 and 25 Hz. We analyzed separately foreshocks andaftershocks and we found the same attenuation functions, suggesting thatthe possible temporal variations could be confined in a small area. Thefrequency dependence of Q _S was approximated by the equation Q _S=18 · f ^2.0between 1 and 10 Hz. At higher frequencies (10–25 Hz), the frequencydependence of Q _s weakens, having an average value of Q _S=990. We also estimated Q from coda waves (Q _C) using the single-scattering models of Aki andChouet (1975) and Sato (1977). We found that Q _C=77 · f ^0.6, (between 2 and 20Hz) at the western side of the mountain chain, using either foreshocks oraftershocks. This relation is consistent with previous estimates of Q _Creported for the Central Apennines. For a volume sampling the Colfioritobasin, the Apennines and the Marche region we found that Q _C=55 · f ^0.8,indicating highattenuation below the mountain belt. To detect small temporal changes ofQ, we calculated spectral ratios of 5 temporal doublets located in theepicentral area and recorded at the closest station. We found temporalchanges of Q that vary from 27% to 56%, depending on the locationof the doublets. This variability suggests that the temporal change ofattenuation may depend on the spatial variation of Q and perhaps on thespatial distribution of tectonic stress in the epicentral area.     

Intrinsic Qp seismic attenuation from the rise time of microearthquakes: a local scale application at Rio‐Antirrio,Western Greece

We determine the attenuation structure of a three‐dimensional medium based on first pulse‐width measurements from microearthquake data. Ninety‐five microearthquakes from a seventy stations local network were considered in this study. Measurements of the first half cycle of the wave, the so‐called rise time τ were carried out on high quality velocity seismograms and inverted to estimate the P‐waves intrinsic quality factor Q_p. The results of this investigation indicate that first pulse width data from a local microearthquake network permit retrieval with sufficient accuracy of the heterogeneous Q_p structure. The inferred attenuation variability corresponds to the known geological formations in the region.     

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.     

Comparison of the attenuation properties for two different areas in eastern Anatolia,Turkey

In this study, the attenuation properties of the crust and the quality factor of S wave in eastern Anatolia (Turkey) were determined by local earthquakes for two different areas, Oltu and Erzurum. Seismic wave attenuation can be changed with high pressure or structural effects. Therefore, we argued that the estimation of attenuation coefficient in seismic active zones in Eastern Anatolia is a very useful tool to determine seismic activities. It uses regional waveform data set from two stations, OLT and ERZ, for 95 events that occurred in these regions between 2001 and 2005. The attenuation has been determined using the Chobra–Alexeev model based on the epicenter distance–amplitude relations. This model allows for investigation of the effects of variations in attenuation properties for different areas. We introduced a new magnitude formula for these areas using the amplitude normalization methods for reference values M_L=4, so as to correct effects of the magnitudes. We also determined velocity of seismic waves. The average attenuation coefficient (α), average quality factor (Q_s) and P and S waves velocities were obtained with normalized amplitude values for Erzurum (ERZ) and Oltu (OLT) as 0.0135 km⁻¹, 37, 6.20 km/s and 3.38 km/s and 0.0151, 34, 6.13 and 3.48.     

Source/receiver array directivity effects on marine seismic attenuation measurements

Attenuation of seismic waves, quantified by the seismic quality factor Q, holds important information for seismic interpretation, due to its sensitivity to rock and fluid properties. A recently published study of Q, based on surface seismic reflection data, used a modified spectral ratio approach (QVO), but both source and receiver responses were treated as isotropic, based on simple raypath arguments. Here, this assumption has been tested by computing apparent attenuation generated by frequency-dependent directivity of typical marine source and receiver arrays and acquisition geometries. Synthetic wavelet spectra were computed for reflected rays, summed over the first Fresnel zone, from the base of a single interval, 50–3000 m thick and velocity 2000 m/s, overlying a 2200 m/s half-space, and for offsets of 71–2071 m. The source and receiver geometry were those of an actual survey. The modelled spectra are clearly affected by directivity, most strongly because of surface ghosts. In general, the strong high-frequency component, produced by the array design, leads to apparently negative attenuation in individual reflection events, though this is dependent on offset and target depth. For shallow targets (less than 400–500 ms two-way traveltime (TWT) depth), apparent Q-values as extreme as ?50 to ?100 were obtained. For deeper target depths, the directivity effect is far smaller. The implications of the model study were tested on real data. QVO was applied to 20 true-spectrum-processed CMPs, in a shallow (405–730 ms TWT) and a deeper (1000–1300 ms TWT) interval, firstly using a measured far-field source signature (effectively isotropic), and secondly using computed directivity effects instead. Mean interval Q^?1-values for the deeper interval, 0.029 ± 0.011 and 0.027 ± 0.018 for conventional and directional processing, respectively, suggested no directivity influence on attenuation estimation. For the shallow interval (despite poor spectral signal-to-noise ratios and hence scattered attenuation estimates), directional processing removed directivity-generated irregularities from the spectral ratios, resulting in an improvement from Q^?1_int = ?0.036 ± 0.130 to a realistic Q^?1_int = 0.012 ± 0.030: different at 94% confidence level. Equivalent Q-values are: for the deeper interval, 35 and 37 for conventional and directional processing, respectively, and ?28 and 86 for the shallow interval. These results support the conclusions of the model studies, i.e. that source/receiver directivity has a negligible effect except for shallow targets (e.g. TWT depth ≤ 500 ms) imaged with conventional acquisition geometry. In such cases directivity corrections to spectra are strongly recommended.     

Attenuation of High Frequency P and S Waves in the Gujarat Region,India

The local earthquake waveforms recorded on broadband seismograph network of Institute of Seismological Research in Gujarat, India have been analyzed to understand the attenuation of high frequency (2–25 Hz) P and S waves in the region. The frequency dependent relationships for quality factors for P (Q _P) and S (Q _S) waves have been obtained using the spectral ratio method for three regions namely, Kachchh, Saurashtra and Mainland Gujarat. The earthquakes recorded at nine stations of Kachchh, five stations of Saurashtra and one station in mainland Gujarat have been used for this analysis. The estimated relations for average Q _P and Q _S are: Q _P = (105 ± 2) f ^{0.82 ± 0.01}, Q _S = (74 ± 2) f ^{1.06 ± 0.01} for Kachchh region; Q _P = (148 ± 2) f ^{0.92 ± 0.01}, Q _S = (149 ± 14) f ^{1.43 ± 0.05} for Saurashtra region and Q _P = (163 ± 7) f ^{0.77 ± 0.03}, Q _S = (118 ± 34) f ^{0.65 ± 0.14} for mainland Gujarat region. The low Q (<200) and high exponent of f (>0.5) as obtained from present analysis indicate the predominant seismic activities in the region. The lowest Q values obtained for the Kachchh region implies that the area is relatively more attenuative and heterogeneous than other two regions. A comparison between Q _S estimated in this study and coda Q (Qc) previously reported by others for Kachchh region shows that Q _C > Q _S for the frequency range of interest showing the enrichment of coda waves and the importance of scattering attenuation to the attenuation of S waves in the Kachchh region infested with faults and fractures. The Q _S/Q _P ratio is found to be less than 1 for Kachchh and Mainland Gujarat regions and close to unity for Saurashtra region. This reflects the difference in the geological composition of rocks in the regions. The frequency dependent relations developed in this study could be used for the estimation of earthquake source parameters as well as for simulating the strong earthquake ground motions in the region.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.     

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.     

Lateral variations of coda Q and attenuation of seismic waves in Southwest Anatolia

The seismic quality factor (Q _c) and the attenuation coefficient (δ) in the earth’s crust in southwest (SW) Anatolia are estimated by using the coda wave method based on the decrease of coda wave amplitude by time on the seismogram. The quality factor Q _o, the value of Q _c at 1 Hz, and its frequency dependency η are determined from this method depending on the attenuation properties of scattered coda waves. δ is determined from the observations of amplitude variations of seismic waves. In applying the coda wave method, firstly, a type curve representing the average pattern of the individual coda decay curves for 0.75, 1.5, 3.0, 6.0, 12.0, and 24.0 Hz values was estimated. Secondly, lateral variation of coda Q and the attenuation coefficients for three main tectonic patterns are estimated. The shape of the type curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The Q _o and η values vary from 30 to 180 and from 0.55 to 1.25, respectively for SW Anatolia. In SW Anatolia, coda Q–f relation is described by and δ = 0.008 km⁻¹. These results are expected to help in understanding the degree of tectonic complexity of the crust in SW Anatolia.     

Attenuation of Coda Waves at the Changbaishan Tianchi Volcanic Area in Northeast China

21 earthquakes recorded by a temporary seismic network in the Changbaishan Tianchi volcanic area in Northeast China operated during the summer of 2002 and 2003 were analyzed to estimate the S coda attenuation. The attenuation quality factor Q_c was estimated using the single scattering attenuation model of Sato (1977) in the frequency band from 4 to 24 Hz. All the events studied in this paper occurred at depths from 2 to 6 km with ML of 1.4–2.8. The epicentral distances are less than 25 km. For all events which occurred near the Tianchi Lake (caldera), the Q_c patterns obtained at the stations near the lake are similar, and the Q_c values are relatively small. At the stations located about 15 km east of the Tianchi Lake, however, the average Q_c is significantly higher. For an event which occurred 25km from the lake to the west, Q_c patterns derived at the stations near the lake are quite similar to the above mentioned Q_c for stations located in the east. Further study shows that Q_c value in the north and central areas of the volcano is relatively lower than that in the surrounding area. Compared to other volcanic areas in the world, the average Q_c of the Changbaishan Tianchi volcanic area is obviously lower. The deep seismic sounding and teleseismic receiver function studies indicated more than one lower velocity layer in the crust. The MT studies suggested the presence of high conductive bodies beneath the area. We interpret the strong attenuation of coda waves near the Changbaishan Tianchi volcano as being possibly related to high temperature medium caused by shallow magma chambers.     

Joint cross-well and single-well seismic studies of CO2 injection in an oil reservoir

A series of time‐lapse seismic cross‐well and single‐well experiments were conducted in a diatomite reservoir to monitor the injection of CO₂ into a hydrofracture zone, based on P‐ and S‐wave data. A high‐frequency piezo‐electric P‐wave source and an orbital‐vibrator S‐wave source were used to generate waves that were recorded by hydrophones as well as 3‐component geophones. During the first phase the set of seismic experiments was conducted after the injection of water into the hydrofractured zone. The set of seismic experiments was repeated after a time period of seven months during which CO₂ was injected into the hydrofractured zone. The questions to be answered ranged from the detectability of the geological structure in the diatomic reservoir to the detectability of CO₂ within the hydrofracture. Furthermore, it was intended to determine which experiment (cross‐well or single‐well) is best suited to resolve these features. During the pre‐injection experiment, the P‐wave velocities exhibited relatively low values between 1700 and 1900 m/s, which decreased to 1600–1800 m/s during the post‐injection phase (?5%). The analysis of the pre‐injection S‐wave data revealed slow S‐wave velocities between 600 and 800 m/s, while the post‐injection data revealed velocities between 500 and 700 m/s (?6%). These velocity estimates produced high Poisson's ratios between 0.36 and 0.46 for this highly porous (～50%) material. Differencing post‐ and pre‐injection data revealed an increase in Poisson's ratio of up to 5%. Both velocity and Poisson's ratio estimates indicate the dissolution of CO₂ in the liquid phase of the reservoir accompanied by an increase in pore pressure. The single‐well data supported the findings of the cross‐well experiments. P‐ and S‐wave velocities as well as Poisson's ratios were comparable to the estimates of the cross‐well data. The cross‐well experiment did not detect the presence of the hydrofracture but appeared to be sensitive to overall changes in the reservoir and possibly the presence of a fault. In contrast, the single‐well reflection data revealed an arrival that could indicate the presence of the hydrofracture between the source and receiver wells, while it did not detect the presence of the fault, possibly due to out‐of‐plane reflections.     

The relationships between the velocities,attenuations and petrophysical properties of reservoir sedimentary rocks1

We have measured the velocities and attenuations of compressional and shear waves in 29 water-saturated samples of sandstones and shales at a confining pressure of 60 MPa and at frequencies of about 0.85 MHz. The measurements were made using a pulse echo method in which the samples (diameter 5 cm, length 1.5 cm to 2.5 cm) were placed between perspex buffer rods inside a high-pressure cell. The velocity of each seismic wave was determined from the traveltime difference of equivalent phase points (corrected for diffraction effects) of the signals reflected from the top and from the base of each sample. Attenuation was determined in a similar way by comparison of the diffraction corrected amplitudes of the signals. The attenuation data are presented as ‘quality factors’: Q_p and Q_s for compressional and shear waves respectively. The results show that Q_s is strongly correlated with V_s, that Q_p is weakly correlated with V_p, and that Q_p is strongly correlated with Q_s. Q_p is strongly dependent on the volume percentage of the assemblage of intra-pore minerals, whether they are clays or carbonates. It is concluded that the attenuation mechanism is due to the local fluid flow arising from the differential dilation of the solid rock frame and the intra-pore mineral assemblage, which is a result of their very different elastic moduli.     

Q P and Q S in Beijing and Jianchuan, Yunnan areas

Events from the December 1982 Huairou County, Beijing, and the July 1982 Jianchuan, Yunnan earthquake series were recorded at one station in Beijing City and at four stations in Eryuan area, Yunnan, respectively. Dividing the spectra (for P and S waves) from the smaller events in the series by spectra (for corresponding P and S waves) from the larger events in cach region, we have determined the high frequency source spectral decay rate to be ω⁻¹in both cases. Through trial and error method, we are able to determine the appropriate constantQ _P andQ _S that corrects the individual spectra to the proper high frequency decay rate. It is found that aQ _P of 800 and aQ _S of 550 can adequately compensate for the attenuation of the waves in Beijing area. For Jianchuan, Yunnan area the corresponding values are 900 and 400. This subject is supported by Chinese Joint Seismological Science Foundation.