广东地区地震动衰减和场地响应的研究

根据广东数字地震台网14个地震台的249条波形资料，研究了广东地区的衰减模型和各台站的场地响应. 采用三段几何衰减模型拟合，得到了广东地区的几何衰减函数. 同时也采用线性几何衰减模型和两段几何衰减模型对观测数据进行了拟合. 残差比较表明，线性几何衰减模型的残差最大，两段几何衰减模型的残差其次，而三段几何衰减模型的残差最小. 由此可以推论广东地区采用三段几何衰减模型具合理性和适用性. 得到广东地区非弹性衰减Q值随频率f的关系为Q(f)=481.5·f0.31；14个台站的场地响应均没有显示出明显的放大效应，这与它们均处于岩石地基相符.     

Can We Improve Estimates of Seismological Q Using a New “Geometrical Spreading” Model?

Precise measurements of seismological Q are difficult because we lack detailed knowledge on how the Earth’s fine velocity structure affects the amplitude data. In a number of recent papers, Morozov (Geophys J Int 175:239–252, 2008; Seism Res Lett 80:5–7, 2009; Pure Appl Geophys, this volume, 2010) proposes a new procedure intended to improve Q determinations. The procedure relies on quantifying the structural effects using a new form of geometrical spreading (GS) model that has an exponentially decaying component with time, e ^?γt·γ is a free parameter and is measured together with Q. Morozov has refit many previously published sets of amplitude attenuation data. In general, the new Q estimates are much higher than previous estimates, and all of the previously estimated frequency-dependence values for Q disappear in the new estimates. In this paper I show that (1) the traditional modeling of seismic amplitudes is physically based, whereas the new model lacks a physical basis; (2) the method of measuring Q using the new model is effectively just a curve fitting procedure using a first-order Taylor series expansion; (3) previous high-frequency data that were fit by a power-law frequency dependence for Q are expected to be also fit by the first-order expansion in the limited frequency bands involved, because of the long tails of power-law functions; (4) recent laboratory measurements of intrinsic Q of mantle materials at seismic frequencies provide independent evidence that intrinsic Q is often frequency-dependent, which should lead to frequency-dependent total Q; (5) published long-period surface wave data that were used to derive several recent Q models inherently contradict the new GS model; and (6) previous modeling has already included a special case that is mathematically identical to the new GS model, but with physical assumptions and measured Q values that differ from those with the new GS model. Therefore, while individually the previous Q measurements have limited precision, they cannot be improved by using the new GS model. The large number of Q measurements by seismologists are sufficient to show that Q values in the Earth are highly laterally variable and are often frequency dependent.     

Strongly Frequency-Dependent Intrinsic and Scattering Q for the Western Margin of the Sichuan Basin, China, as Revealed by a Phenomenological Model

Through applying the energy-flux model to teleseismic P-wave data, seismological scattering Q and intrinsic Q for the western margin of the Sichuan basin in southwestern China are estimated separately. The results exemplify the frequency-dependent Q values in the seismic frequency band. The scattering Q is heavily frequency-dependent and reaches a minimum at ~1.5 Hz. The intrinsic Q is nearly invariant for frequencies above 2 Hz and rapidly increases at lower frequencies. A major advantage of the energy-flux model is that it is a phenomenological model based on the conservation of energy and requires no geometrical spreading correction. The results are good counterexamples against the claim of a frequency-independent seismological Q.     

宁夏及邻区地震动衰减和场地响应的分区研究

根据宁夏数字地震台网19个地震台的294条波形资料,研究宁夏南区及北区的衰减模型和各台站的场地响应。采用三段几何衰减模型拟合,得到宁夏南区和北区非弹性衰减Q值随频率f的关系分别为：Q（f）=503．92f^0.32和Q（f）=267．6f^0.580;19个台站的场地响应除泾源地震台外,其余台站没有显示出明显的放大效应,与实际情况基本相符。     

Source and high-frequency decay parameters for the Azores region for stochastic finite-fault ground motion simulations

Strong ground motion prediction based on finite-fault simulation requires the identification of the fault (strike, dip, length and width), source kinematics parameters (stress drop, rupture velocity and slip distribution), regional crustal properties (geometrical spreading, anelastic structure, and upper crustal amplification and attenuation parameters) and the determination of amplification effects due to the local site geology. The general purpose of this study is to understand source and attenuation properties in the Azores, by the determination of stress drop, quality factor and kappa, through records obtained by the Portuguese digital seismic and accelerometer network. Source Spectra were obtained, for each record, after correcting observed spectra from geometrical spreading and anelastic attenuation effect: quality factor was estimated based on coda decay in the time domain and the kappa parameter was estimated by fitting the high-frequency decay of the acceleration spectrum with a straight line in a log-linear scale. Mean stress drop value was obtained considering that ω-squared model for the source spectra prevails. Parameters kappa, k, and quality factor, Q, have been estimated to be k = (0.075 ± 0.02)s and \(Q(f) = (76 \pm 11)f^{0.69 \pm 0.09}\), respectively. A mean value of stress drop was estimated to be around 90–130 bars for the earthquakes in analysis. The knowledge of source and path parameters, in association with other ground motion parameters, allows improving ground-motion estimates for the Azores and, consequently, will lead to more accurate seismic hazard assessment for the Azores and better characterization of seismic scenarios.     

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.     

紫坪铺水库库区介质衰减、台站响应和震源参数研究

利用紫坪铺水库数字地震台网和成都区域数字台网中的油榨坪台2004年8月到2008年5月12日汶川大地震前的紫坪铺库区2级以上地震波形资料,采用Atkinson方法计算了紫坪铺水库库区地震波非弹性衰减系数、几何扩散系数和介质品质因子;利用Moya等提出的方法计算了台站场地响应;对震源谱的低频水平和拐角频率进行了联合反演,并计算了紫坪铺库区287次1.6级以上地震的震源波谱参数,同时讨论了它们之间的关系。结果表明:紫坪铺库区Q值与频率的关系为Q(f)=47.8×f0.91;场地响应均表现为与频率相关;紫坪铺库区小震地震矩、震级与震源谱拐角频率的依赖关系不明确,这可能与水库诱发地震有关     

Compensating elastic transmission losses for P-wave attenuation estimation from sonic logs

The decay of seismic amplitude is caused by a variety of physical phenomena that can be divided broadly into elastic transmission losses (including geometrical spreading, interface transmission losses and scattering attenuation) and intrinsic attenuation, where wave energy is converted into heat due to viscous friction. The so-called statistical averaging method is currently considered as the most advanced sonic wave attenuation estimation method, and there exist various implementations of this method. But the way elastic transmission losses – that mask the true intrinsic attenuation – are compensated for appears to be an issue and in some cases this correction has been overlooked. In this paper, we revisit the statistical averaging method for intrinsic attenuation estimation with particular focus on the role of elastic transmission losses. Through synthetic examples, we demonstrate the importance of compensating for elastic transmission losses even if the variation of velocity and density with depth is not notable. Our implementation of the method uses finite-difference simulations thereby providing a versatile and accurate way to generate synthetic seismograms. We use a combination of elastic and viscoelastic finite-difference simulations to demonstrate the significant error without accurate compensation of the elastic transmission losses. We apply our implementation of the method to sonic waveforms acquired in an exploration well from Browse basin, Australia. The resulting intrinsic attenuation estimates are indeed indicative of gas-saturated zones identified from petrophysical analysis in which viscous friction are thought to be of importance.     

SEISMIC REFLECTION AMPLITUDES1

Seismic amplitude variations with offset contain information about the elastic parameters. Prestack amplitude analysis seeks to extract this information by using the variations of the reflection coefficients as functions of angle of incidence. Normally, an approximate formula is used for the reflection coefficients, and variations with offset of the geometrical spreading and the anelastic attenuation are often ignored. Using angle of incidence as the dependent variable is also computationally inefficient since the data are recorded as a function of offset. Improved approximations have been derived for the elastic reflection and transmission coefficients, the geometrical spreading and the complex travel-time (including anelastic attenuation). For a 1 D medium, these approximations are combined to produce seismic reflection amplitudes (P-wave, S-wave or converted wave) as a Taylor series in the offset coordinate. The coefficients of the Taylor series are computed directly from the parameters of the medium, without using the ray parameter. For primary reflected P-waves, dynamic ray tracing has been used to compute the offset variations of the transmission coefficients, the reflection coefficient, the geometrical spreading and the anelastic attenuation. The offset variation of the transmission factor is small, while the variations in the geometrical spreading, absorption and reflection coefficient are all significant. The new approximations have been used for seismic modelling without ray tracing. The amplitude was approximated by a fourth-order polynomial in offset, the traveltime by the normal square-root approximation and the absorption factor by a similar expression. This approximate modelling was compared to dynamic ray tracing, and the results are the same for zero offset and very close for offsets less than the reflector depth.     

Intrinsic attenuation and scattering of shear waves in the lithosphere of Kamchatka

Records of small local Kamchatka earthquakes were processed using Multiple Lapse Time Window Analysis (MLTWA). The method makes use of normalized integrals of 3D seismic energy density in several windows applied to an earthquake record that has been put through a bandpass filter. The intrinsic attenuation and scattering properties of the earth were estimated by choosing parameters that provide the best agreement between experimental and theoretical integrals as functions of hypocentral distance. The theoretical normalized integrals as functions of distance were computed using an analytical solution to the equation of seismic energy transport based on a simple model of isotropic scattering for scalar seismic waves excited by an impulse point source in an earth that is uniform as to scattering and intrinsic attenuation. The resulting estimates of attenuation and scattering parameters are similar to those derived in other studies of Kamchatka, as well as of other continents: America (California, US), Europe (the Canary Is., Spain, and Italy), and Asia (Japan).     

小波尺度域含气储层地震波衰减特征

黏弹性衰减因子Q的可靠估计可通过Q反褶积来提高地震资料的分辨率并有助于振幅分析. 本文从小波理论出发,结合地震波在黏弹性介质中的传播方程,推导出小波尺度域地震波能量衰减公式. 能量衰减公式具有下列性质：（1）Q值越大,能量衰减得越慢;Q值越小,能量衰减越严重;（2）尺度越小,信号中保留的能量越少;（3）对于脉冲源来说在理想的无衰减介质(即Q趋近于∞)中传播时,信号在不同尺度内的能量相同. 利用尺度能量公式,可从反射地震资料中直接估计品质因子Q（即衰减因子）,也可以提取不同尺度的能量衰减剖面作为储层描述的属性参数,用来进行岩性识别和指示气藏,与经典的谱比法相比,避免了谱比法所面临的双时窗问题以及进行谱估计的窗选择问题. 理论模型试验表明了本文方法的正确性和有效性.     

DISCUSSION ON ATTENUATION CHARACTERISTICS,SITE RESPONSE AND MAGNITUDE DETERMINATION IN SICHUAN

The attenuation characteristics and site response are calculated respectively for each individual tectonic unit in Sichuan (Sichuan Basin,west Sichuan plateau and Panzhihua-Xichang area),using digital waveform data recorded by regional seismic networks and relevant seismic phase data collected from China Seismograph Network.The frequency dependent Q(f) is obtained by the iterative grid-search technique described by Atkinson and Mereu based on trilinear geometrical spreading model.The source spectra are determined by the model of Brune and the site responses of seismic stations are derived by Moya's method using genetic algorithms.Comparison to conventional M_L estimates shows that the network local magnitude bias is quite significant at low and intermediate magnitudes.The bias at the j^th station for the i^th event is defined as ΔM_ij=M_ij-M_i, where ΔM_ij is the station magnitude and M_i the network-average value.For comparison,we mapped the spatial distribution of biases by digital seismograms recorded from 10535 earthquakes of magnitude 2.5≤M_L≤4.9 that occurred in Sichuan from January 1,2009 to June 30,2015.Based on the above data,the attenuation characteristics,site response and their effects on magnitude determination in Sichuan are analyzed.Our results demonstrate that the associated model for regional quality factor for frequencies can be expressed as Q₁(f)=450.6f^{0.513 4} for Sichuan Basin,Q₂(f)=136.6f^{0.581 3} for west Sichuan Plateau and Q₃(f)=101.9f^{0.666 3} for Panzhihua-Xichang area.Site response results indicate that different stations show different amplifications.Maps of biases appear to be different,but with similar dominant spatial distribution.For stations in Sichuan Basin,their greater magnitudes are functions of low attenuation in structure and amplification effects of both seismic stations and basin effects.For stations in west Sichuan Plateau,the possible causes of these lower magnitudes are severe dependence upon source region due to extreme lateral variations in either structure or path effect attenuation.For stations in Panzhihua-Xichang area,broken medium caused by strong tectonic activity or large earthquakes and heat flow up-welling along active faults may be the main reasons of low magnitude values when earthquakes occur in western Sichuan and eastern Tibetan region.And the greater magnitudes for earthquakes along the Longmen Mountains appear to be well correlated with edge effect of sedimentary basin on strong ground motion.In our study,stations magnitude biases appear to be extremely correlated with tectonic structures and different regions when seismic rays passing through,magnitudes are affected significantly by lateral variations in attenuation characters rather than site responses.     

Estimation of shear-wave attenuation characteristics for the Klyuchevskoi volcanic edifice

This paper reports a study of intrinsic attenuation and scattering in the Klyuchevskoi volcanic edifice. The data set consisted of small volcano-tectonic earthquakes occurring as deep as 30 km beneath Klyuchevskoi Volcano and recorded by radio telemetry seismograph stations installed on the edifice and near it. The digital seismograms of small volcano-tectonic earthquakes were processed by the multiple lapse time window analysis (MLTWA) method currently in use in seismology. The method uses experimental normalized integrals of 3D seismic energy density determined from several time windows applied to earthquake records that have been put through bandpass filters. The parameters that characterize the intrinsic attenuation and scattering of seismic energy of small volcano-tectonic earthquakes in Klyuchevskoi Volcano were estimated by adjusting the coefficients to ensure the best fit between experimental and theoretical integrals. An analytical solution to the seismic energy transfer equation was used to calculate the theoretical integrals of 3D seismic energy. Reliable estimates of the parameters that characterize the intrinsic attenuation and scattering in the material of the Klyuchevskoi volcanic edifice have been obtained. Our estimates of the Q for the edifice are below those derived by other workers for the Kamchatka lithosphere.     

Lg波衰减模型中建模误差的定量评估

对Lg波衰减模型中建模误差的统计特征进行了系统研究,并建立了地壳二维Lg波衰减模型。由于Lg波振幅可能受到几何扩散函数的强烈影响,合理评估反演过程的误差对于能否使用最小二乘意义下的反演非常重要。通过在川滇及其邻近地区收集的建模误差样本,使用K-S数值检验方法、Q-Q图和正态分布图形检验方法对Lg波衰减层析成像反演的输入数据中建模误差的分布特征进行了统计分析。采用奇异值分解(SVD)和反投影方法,分别获得了川滇地区的Q_Lg模型,定量计算模型的协方差矩阵和分辨率矩阵,定量评估了Q_Lg模型中每个格点的分辨率和误差。结果表明：在一阶近似条件下建模误差服从正态分布;通过开发的数据筛选程序,可以产生一个接近完美正态分布的数据集;与反投影方法相比,利用SVD方法获得的地壳Q值的分辨率更高;在射线覆盖较好的区域,Q_Lg模型的分辨率达到100 km,相对误差小于3%。     

Attenuation models inferred from intraplate earthquake recordings

Strong-motion recordings at 87 sites from 56 different intraplate earthquakes from North America, Europe, China and Australia have been used through a two-step regression analysis to develop new attenuation models for peak ground acceleration, and for pseudo-relative velocity for frequencies of 0.25, 0.5, 1.0, 2.0, 5.0 and 10.0 Hz, all for 5 per cent of critical damping. The estimates are obtained along with an analysis of residuals and scatter. A similar regression analysis has been performed also for Fourier spectra of acceleration, in which case the coefficient for the anelastic term has been interpreted in terms of a frequency dependent quality factor Q. The resulting Q-model shows a strong frequency sensitivity with values around 600–700 at 1 Hz, around 2000 at 10 Hz and around 5200 at 25 Hz. These PGA, PSV and Q results depend, however, on the underlying assumption for geometrical spreading, in particular for low frequencies.     

Monitoring viscosity changes from time‐lapse seismic attenuation: case study from a heavy oil reservoir

Heating heavy oil reservoirs is a common method for reducing the high viscosity of heavy oil and thus increasing the recovery factor. Monitoring of these viscosity changes in the reservoir is essential for delineating the heated region and controlling production. In this study, we present an approach for estimating viscosity changes in a heavy oil reservoir. The approach consists of three steps: measuring seismic wave attenuation between reflections from above and below the reservoir, constructing time‐lapse Q and Q^?1 factor maps, and interpreting these maps using Kelvin–Voigt and Maxwell viscoelastic models. We use a 4D relative spectrum method to measure changes in attenuation. The method is tested with synthetic seismic data that are noise free and data with additive Gaussian noise to show the robustness and the accuracy of the estimates of the Q‐factor. The results of the application of the method to a field data set exhibit alignment of high attenuation zones along the steam‐injection wells, and indicate that temperature dependent viscosity changes in the heavy oil reservoir can be explained by the Kelvin–Voigt model.     

基于小波包分解的地层吸收补偿

在分析地震反射波勘探信号小波包分解特性的基础上,提出了基于小波包分解的地层吸收补偿方法．如果没有地层吸收,浅层反射波和深层反射波具有相同的振幅谱而相位谱相差一线性相位．如果把地震记录分成不同的频段,不同频段对于时间的能量分布关系具有相似性,即对所有频段的深层反射的能量与同一频段浅层反射能量之比应该相同;只是不同频段的绝对能量大小不一样．由于地层的吸收造成各频段能量对时间的分布不同,如果给这些频段乘以时变权使各频段能量对时间的分布相同,就起到了对地层吸收的补偿作用．该方法首先用小波包对地震勘探信号进行分频,然后根据地层吸收的特点,在每一个频段上分别提取每一个时间点的吸收系数并用其倒数加权每一个小波包系数,再用加权后的小波包系数重构地震反射波记录消除地层吸收．该方法适合Ｑ值未知的情况,对变Ｑ和常Ｑ都能适应．从实际应用效果看,该方法能很好的补偿地层的吸收。     

随机弹性介质中地震波散射衰减分析(英文)

地震波衰减一直是许多学科研究的热点,因为可以反映介质的特性。导致地震波衰减的因素很多,如:传播过程中由于能量扩散导致的几何衰减,固体岩石内部晶粒间相对滑移导致的摩擦衰减,岩石结构不均匀引起的地震波散射衰减。本文主要从统计的观点出发,通过多次数值模拟的方法研究纵波散射在随机弹性介质中所引发的衰减。首先用随机理论建立了二维空间随机弹性介质模型,然后用错格伪谱法的数值方法模拟了波在随机介质中的传播,再通过波场中虚拟检波器的记录,用谱比法估计了弹性波在随机介质中的散射衰减。不同非均匀程度随机弹性介质中的数值结果表明:介质不均匀程度越高,散射衰减越大;在散射体尺寸小于波长的前提下,不同散射体尺寸的计算结果说明:散射体尺寸越大,弹性波衰减越明显。最后提出了一种不均匀孔隙介质中流体流动衰减的方法。通过对随机孔隙介质中地震波的总衰减和散射衰减分别进行了计算,并定量得出了随机孔隙介质中流体流动衰减,结果表明:在实际地震频段下,当介质不均匀尺度101米量级时,散射衰减比流体流动衰减要大,散射衰减是地震波在实际不均匀岩石孔隙介质中衰减的主要原因。     

2.5D modelling, inversion and angle migration in anisotropic elastic media

2.5D modelling approximates 3D wave propagation in the dip‐direction of a 2D geological model. Attention is restricted to raypaths for waves propagating in a plane. In this way, fast inversion or migration can be performed. For velocity analysis, this reduction of the problem is particularly useful. We review 2.5D modelling for Born volume scattering and Born–Helmholtz surface scattering. The amplitudes are corrected for 3D wave propagation, taking into account both in‐plane and out‐of‐plane geometrical spreading. We also derive some new inversion/migration results. An AVA‐compensated migration routine is presented that is simplified compared with earlier results. This formula can be used to create common‐image gathers for use in velocity analysis by studying the residual moveout. We also give a migration formula for the energy‐flux‐normalized plane‐wave reflection coefficient that models large contrast in the medium parameters not treated by the Born and the Born–Helmholtz equation results. All results are derived using the generalized Radon transform (GRT) directly in the natural coordinate system characterized by scattering angle and migration dip. Consequently, no Jacobians are needed in their calculation. Inversion and migration in an orthorhombic medium or a transversely isotropic (TI) medium with tilted symmetry axis are the lowest symmetries for practical purposes (symmetry axis is in the plane). We give an analysis, using derived methods, of the parameters for these two types of media used in velocity analysis, inversion and migration. The kinematics of the two media involve the same parameters, hence there is no distinction when carrying out velocity analysis. The in‐plane scattering coefficient, used in the inversion and migration, also depends on the same parameters for both media. The out‐of‐plane geometrical spreading, necessary for amplitude‐preserving computations, for the TI medium is dependent on the same parameters that govern in‐plane kinematics. For orthorhombic media, information on additional parameters is required that is not needed for in‐plane kinematics and the scattering coefficients. Resolution analysis of the scattering coefficient suggests that direct inversion by GRT yields unreliable parameter estimates. A more practical approach to inversion is amplitude‐preserving migration followed by AVA analysis. SYMBOLS AND NOTATION A list of symbols and notation is given in Appendix D .     

岩石的粘弹性谐振Q模型

实验表明,在应变小于10^-6范围内,砂岩对地震波的吸收主要由孔隙流体的局部运动引起,而且Q值随频率的变化出现谐振现象。据此,作者认为,地壳上部地震波能量的耗损可以表示为以下两种成分的线性组合:1.由滑动摩擦、热弹性驰豫、位错运动等引起的能量消散,它们主要与岩石的固相成分有关,可近似地用恒Q模型描述;2.由孔隙流体运动引起的能量消散,可以用谐振公式描述。根据这一认识,利用描述因果关系的Kramers-Krnig关系式可推导出表示地壳上部岩石粘弹性（复弹性模量,相速度频散和衰减函数）的公式,它们综合地描述了由各种机制引起的波的频散和吸收,并在谐振Q值等于参考常数Q值时退化为目前常用的Futterman模型。作为这种谐振Q模型的应用,介绍了它用于Q值测量结果外推和频散一吸收研究以及粘弹性介质中反射地震道合成的结果。