We establish a patchy saturation model and derive the seismic wave equations for patchy saturated porous media on the basis of Biot's equations and Johnson's bulk modulus. We solve the equations, obtain the attenuat...We establish a patchy saturation model and derive the seismic wave equations for patchy saturated porous media on the basis of Biot's equations and Johnson's bulk modulus. We solve the equations, obtain the attenuation coefficients, and analyze the characteristics of wave attenuation in the seismic frequency range. The results suggest that seismic waves show attenuation and dispersion in partially saturated rocks in the low frequency range. With frequency increasing, attenuation increases. The attenuation of P-waves of the second kind is more pronounced in agreement with Biot's theory. We also study the effect of porosity, saturation, and inner sphere radius on the attenuation of the P-waves of the first kind and find that attenuation increases with increasing frequency and porosity, and decreases with increasing frequency and degree of saturation. As for the inner sphere radius, wave attenuation is initially increasing with increasing frequency and inner sphere radius less than half the outer radius. Subsequently, wave attenuation decreases with increasing frequency and inner sphere radius is higher than half the outer sphere radius.展开更多
The constant Q property in viscoelastic media assumes that the quality factor Q does not change with frequency(i.e.,the Q value is independent of the frequency).For seismic waves propagating in viscoelastic media,the ...The constant Q property in viscoelastic media assumes that the quality factor Q does not change with frequency(i.e.,the Q value is independent of the frequency).For seismic waves propagating in viscoelastic media,the wave equation is determined by the viscoelastic media model.Equivalence relations exist between various frequency domain mathematical models and physical rheological models for the constant Q property.Considering two elastic moduli and three attenuation variables,24 kinds of wave equations based on diff erent generalized rheological models are divided into six classes in this study,and the 12 kinds of specifi c representation for the wave equations in the time domain are derived.On the basis of the equivalence relations between the generalized rheological models,the diff erence and equivalence relation between diff erent wave equations are proven and clarifi ed.Results show that the high-order generalized rheological model can accurately characterize the attenuation characteristics of seismic waves and has advantages in characterizing the dispersion characteristics in viscoelastic media.Lastly,the seismic refl ection characteristics caused by the diff erence of Q value are verifi ed by the forward modeling of the constant Q wave equation in this study,thereby providing a theoretical basis for the analysis and inversion of the formation Q value from refl ection seismic data.展开更多
3-D S-wave Q structure in Jiashi earthquake region is inverted based on the attenuation of seismic waves recorded from earthquakes in this region in 1998 by the Research Center of Exploration Geophysics (RCEG), CSB, a...3-D S-wave Q structure in Jiashi earthquake region is inverted based on the attenuation of seismic waves recorded from earthquakes in this region in 1998 by the Research Center of Exploration Geophysics (RCEG), CSB, and a rough configuration of deep crustal faults in the earthquake region is presented. First, amplitude spectra of S-waves are extracted from 450 carefully-chosen earthquake records, called observed amplitude spectra. Then, after instrumental and site effect correction, theoretical amplitude spectra are made to fit observed amplitude spectra with nonlinear damped least-squares method to get the observed travel time over Q, provided that earthquake sources conform to Brune's disk dislocation model. Finally, by 3-D ray tracing method, theoretical travel time over Q is made to fit observed travel time over Q with nonlinear damped least-squares method. In the course of fitting, the velocity model, which is obtained by 3-D travel time tomography, remains unchanged, while only Q model is modified. When fitting came to the given accuracy, the ultimate Q model is obtained. The result shows that an NE-trending low Q zone exists at the depths of 10-18 km, and an NW-trending low Q zone exists at the depths of 12-18 km. These roughly coincide with the NE-trending and the NW-trending low velocity zones revealed by other scientists. The difference is that the low Q zones have a wider range than the low velocity zones.展开更多
基金supported by the National Natural Science Foundation of China(Nos.41204089 and 41174087)the National Science and Technology Major Project(Nos.2011ZX05035-001 and 2011ZX05005-005)the National 863 Program(No.2013AA064201)
文摘We establish a patchy saturation model and derive the seismic wave equations for patchy saturated porous media on the basis of Biot's equations and Johnson's bulk modulus. We solve the equations, obtain the attenuation coefficients, and analyze the characteristics of wave attenuation in the seismic frequency range. The results suggest that seismic waves show attenuation and dispersion in partially saturated rocks in the low frequency range. With frequency increasing, attenuation increases. The attenuation of P-waves of the second kind is more pronounced in agreement with Biot's theory. We also study the effect of porosity, saturation, and inner sphere radius on the attenuation of the P-waves of the first kind and find that attenuation increases with increasing frequency and porosity, and decreases with increasing frequency and degree of saturation. As for the inner sphere radius, wave attenuation is initially increasing with increasing frequency and inner sphere radius less than half the outer radius. Subsequently, wave attenuation decreases with increasing frequency and inner sphere radius is higher than half the outer sphere radius.
基金This work was supported by National Natural Science Foundation of China(No.41774137)111 project(No.B18055),and the Fundamental Research Funds for the Central Universities(No.19CX02002A).
文摘The constant Q property in viscoelastic media assumes that the quality factor Q does not change with frequency(i.e.,the Q value is independent of the frequency).For seismic waves propagating in viscoelastic media,the wave equation is determined by the viscoelastic media model.Equivalence relations exist between various frequency domain mathematical models and physical rheological models for the constant Q property.Considering two elastic moduli and three attenuation variables,24 kinds of wave equations based on diff erent generalized rheological models are divided into six classes in this study,and the 12 kinds of specifi c representation for the wave equations in the time domain are derived.On the basis of the equivalence relations between the generalized rheological models,the diff erence and equivalence relation between diff erent wave equations are proven and clarifi ed.Results show that the high-order generalized rheological model can accurately characterize the attenuation characteristics of seismic waves and has advantages in characterizing the dispersion characteristics in viscoelastic media.Lastly,the seismic refl ection characteristics caused by the diff erence of Q value are verifi ed by the forward modeling of the constant Q wave equation in this study,thereby providing a theoretical basis for the analysis and inversion of the formation Q value from refl ection seismic data.
基金Joint Seismological Science Foundation of China (957-07-414)State Key Basic Research Development and Pro-gramming Project (95-13-02-02).
文摘3-D S-wave Q structure in Jiashi earthquake region is inverted based on the attenuation of seismic waves recorded from earthquakes in this region in 1998 by the Research Center of Exploration Geophysics (RCEG), CSB, and a rough configuration of deep crustal faults in the earthquake region is presented. First, amplitude spectra of S-waves are extracted from 450 carefully-chosen earthquake records, called observed amplitude spectra. Then, after instrumental and site effect correction, theoretical amplitude spectra are made to fit observed amplitude spectra with nonlinear damped least-squares method to get the observed travel time over Q, provided that earthquake sources conform to Brune's disk dislocation model. Finally, by 3-D ray tracing method, theoretical travel time over Q is made to fit observed travel time over Q with nonlinear damped least-squares method. In the course of fitting, the velocity model, which is obtained by 3-D travel time tomography, remains unchanged, while only Q model is modified. When fitting came to the given accuracy, the ultimate Q model is obtained. The result shows that an NE-trending low Q zone exists at the depths of 10-18 km, and an NW-trending low Q zone exists at the depths of 12-18 km. These roughly coincide with the NE-trending and the NW-trending low velocity zones revealed by other scientists. The difference is that the low Q zones have a wider range than the low velocity zones.