Scattering and Diffraction of elastic in-plane P-and SV-waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and ...Scattering and Diffraction of elastic in-plane P-and SV-waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and strong-motion seismologists for over forty years. The case of out-of-plane SH waves on the same elastic canyon that is semi-circular in shape on the half-space surface is the first such problem that was solved by analytic closed form solutions over forty years ago by Trifunac. The corresponding case of in-plane P-and SV-waves on the same circular canyon is a much more complicated problem because, the in-plane P-and SV-scattered waves have different wave speeds and together they must have zero normal and shear stresses at the half-space surface. It is not until recently in 2014 that analytic solution for such problem is found by the author in the work of Lee and Liu. This paper uses the technique of Lee and Liu of defining these stress-free scattered waves to solve the problem of the scattered and diffraction of these in-plane waves on an almost-circular surface canyon that is arbitrary in shape.展开更多
Based on the measurement of the arrival time of maxima magnitude from band-pass filtering signals which were determined using a new Morlet wavelet multiple-filter method, we develop a method for measuring intrinsic an...Based on the measurement of the arrival time of maxima magnitude from band-pass filtering signals which were determined using a new Morlet wavelet multiple-filter method, we develop a method for measuring intrinsic and attenuative dispersion of the first cycle direct P-wave. We determine relative group delays of spectral components of direct P-waves for 984 ray paths from SML and ALS stations of the Taiwan Central Weather Bureau Seismic Network （CWBSN）. Using continuous relaxation model, we deduce a new transfer function that relates intrinsic dispersion to attenuation. Based on the genetic algorithm （GA）, we put forward a new inversion procedure for determining which is defined the flat part of quality factor Q（ω） spectrum, τ1 and τ2 parameters. The results indicate that ① The distribution of Om values versus epicentral distance and depth show that Qm values linearly increase with increasing of epicentral distance and depth, and Qm values is clearly independent of earthquakes magnitude; ② In the different depth ranges, Qm residual show no correlation with variations in epicentral distance. Some significant changes of Qm residual with time is likely caused by pre-seismic stress accumulation, and associated with fluid-filled higher density fractures rock volume in the source area of 1999 Chi-Chi Taiwan earthquake. We confirm that Qm residual with time anomaly appears about 2.5 years before the Chi-Chi earthquake; ③ A comparison of Qm residual for different depth range between SML and ALS stations show that the level of stress has vertical and lateral difference; ④ The area near observation station with both anomalously increasing and decreasing averaged Qm residual is likely an unstable environment for future strong earthquake occurrence. This study demonstrates the capability of direct P-waves dispersion for monitoring attenuation characteristics and its state changes of anelastic medium of the Earth at short propagation distance using seismograms recorded from very small events.展开更多
The indirect boundary element method (IBEM) is used to study the surface motion of an alluvial valley in layered half-space for incident plane P-waves based on Wolf’s theory. Firstly, the free field response can be s...The indirect boundary element method (IBEM) is used to study the surface motion of an alluvial valley in layered half-space for incident plane P-waves based on Wolf’s theory. Firstly, the free field response can be solved by the direct stiffness method, and the scattering wave response is calculated by Green’s functions of distributed loads acting on inclined lines in a layered half-space. The method is verified by comparing its results with literature and numerical analyses are performed by taking the amplification of incident plane P-waves by an alluvial valley in one soil layer resting on bedrock as an example. The results show that there exist distinct differences between the wave amplification by an alluvial valley embedded in layered half-space and that in homogeneous half-space and there is interaction between the valley and the soil layer. The amplitudes are relatively large when incident frequencies are close to the soil layer’s resonant frequencies.展开更多
Finite-frequency travel time tomography is a newly developing method.The main procedure in this new method is to compute the traveltime sensitive kernel.The travel time of the same scatterer needs to be used for compu...Finite-frequency travel time tomography is a newly developing method.The main procedure in this new method is to compute the traveltime sensitive kernel.The travel time of the same scatterer needs to be used for computing the traveltime sensitive kernel many times.It is a time-consuming task.It is easy and fast to get the travel time from analytic equations in a simple model such as a homogenous or linear velocity media.However,most of the earth models are layered.It is cumbersome to get the travel time from analytic equations.In order to enhance the computation efficiency,we used the table look-up method to compute the finite-frequency travel time sensitive kernel for P-waves in a layered structure model.We chose the AK135 earth model for the velocity model.The table look-up method saved about 50% of the computation time.We enhanced the computation speed by using the table lookup method in the same velocity model,which was very useful for enhancing the computation efficiency for the finite-frequency travel time tomography.展开更多
A P-wave velocity model was built in the central southern of the Tanlu Fault based on double-difference tomography.The results suggest the presence of a low-velocity anomaly extending from the surface to a depth of 25...A P-wave velocity model was built in the central southern of the Tanlu Fault based on double-difference tomography.The results suggest the presence of a low-velocity anomaly extending from the surface to a depth of 25 km around the Tanlu and Feixi Faults,representing fault-related fluids caused by partial melting.The relocated earthquakes indicate a significant concentration of seismic activity above 20 km around the Tanlu and Feixi Faults,suggesting that prominent fault systems possibly serve as conduits for the upward migration of deep minerals.The proposed geodynamic model,supported by geological and geophysical data,suggests that the migration of deep mineralized materials extends along the Tanlu Fault.The obtained results serve as a crucial foundation for elucidating the intricate process of mineralization in the central southern segment of the Tanlu Fault,thereby enhancing comprehension regarding the interaction among ore body formation,fault fluids,localized melting,and seismic activity.展开更多
We experimentally studied the effect of crack aperture on P-wave velocity, amplitude, anisotropy and dispersion. Experimental models were constructed based on Hudson＇s theory. Six crack models were embedded with equa...We experimentally studied the effect of crack aperture on P-wave velocity, amplitude, anisotropy and dispersion. Experimental models were constructed based on Hudson＇s theory. Six crack models were embedded with equal-radius penny-shaped crack inclusions in each layer. The P-wave velocity and amplitude were measured parallel and perpendicular to the layers of cracks at frequencies of 0.1 MHz to 1 MHz. The experiments show that as the crack aperture increases from 0.l mm to 0.34 mm, the amplitude of the P-waves parallel to the crack layers decreases linearly with increasing frequency and the P-wave velocity dispersion varies from 1.5% to 2.1%, whereas the amplitude of the P-wave perpendicular to the crack layers decreases quadratically with increasing frequency and the velocity dispersion varies from 1.9% to 4.7%. The variation in the velocity dispersion parallel and perpendicular to the cracks intensifies the anisotropy dispersion of the P-waves in the crack models （6.7% to 83%）. The P-wave dispersion strongly depends on the scattering characteristics of the crack apertures.展开更多
The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including...The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850＋ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth＇s deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.展开更多
It is generally accepted that the uniaxial compressive strength(UCS)and P-wave velocity of rocks tend to decrease simultaneously with increasing temperature.However,based on a great number of statistical data and syst...It is generally accepted that the uniaxial compressive strength(UCS)and P-wave velocity of rocks tend to decrease simultaneously with increasing temperature.However,based on a great number of statistical data and systematic analysis of the microstructure variation of rocks with temperature rising and corresponding propagation mechanism of elastic wave,the results show that(1)There are three different trends for the changes of UCS and P-wave velocity of sandstone when heated from room temperature(20C or 25C)to 800C:(i)Both the UCS and P-wave velocity decrease simultaneously;(ii)The UCS increases initially and then decreases,while the P-wave velocity decreases continuously;and(iii)The UCS increases initially and then fluctuates,while the P-wave velocity continuously decreases.(2)The UCS changes at room temperaturee400C,400Ce600C,and 600Ce800C are mainly attributed to the discrepancy of microstructure characteristics and quartz content,the transformation plasticity of clay minerals,and the balance between the thermal cementation and thermal damage,respectively.(3)The inconsistency in the trends of UCS and P-wave velocity changes is caused by the change of quartz content,phase transition of water and certain minerals.展开更多
Engineering rock mass classification,based on empirical relations between rock mass parameters and engineering applications,is commonly used in rock engineering and forms the basis for designing rock structures.The ba...Engineering rock mass classification,based on empirical relations between rock mass parameters and engineering applications,is commonly used in rock engineering and forms the basis for designing rock structures.The basic data required may be obtained from visual observation and laboratory or field tests.However,owing to the discontinuous and variable nature of rock masses,it is difficult for rock engineers to directly obtain the specific design parameters needed.As an alternative,the use of geophysical methods in geomechanics such as seismography may largely address this problem.In this study,25 seismic profiles with the total length of 543 m have been scanned to determine the geomechanical properties of the rock mass in blocks Ⅰ,Ⅲ and Ⅳ-2 of the Choghart iron mine.Moreover,rock joint measurements and sampling for laboratory tests were conducted.The results show that the rock mass rating（RMR） and Q values have a close relation with P-wave velocity parameters,including P-wave velocity in field(V;).P-wave velocity in the laboratory(V;) and the ratio of V;V;(i.e.K;= V;/V;.However,Q value,totally,has greater correlation coefficient and less error than the RMR,In addition,rock mass parameters including rock quality designation（RQD）,uniaxial compressive strength（UCS）,joint roughness coefficient（JRC） and Schmidt number（RN） show close relationship with P-wave velocity.An equation based on these parameters was obtained to estimate the P-wave velocity in the rock mass with a correlation coefficient of 91%.The velocities in two orthogonal directions and the results of joint study show that the wave velocity anisotropy in rock mass may be used as an efficient tool to assess the strong and weak directions in rock mass.展开更多
Fast and accurate P-wave arrival picking significantly affects the performance of earthquake early warning(EEW)systems.Automated P-wave picking algorithms used in EEW have encountered problems of falsely picking up no...Fast and accurate P-wave arrival picking significantly affects the performance of earthquake early warning(EEW)systems.Automated P-wave picking algorithms used in EEW have encountered problems of falsely picking up noise,missing P-waves and inaccurate P-wave arrival estimation.To address these issues,an automatic algorithm based on the convolution neural network(DPick)was developed,and trained with a moderate number of data sets of 17,717 accelerograms.Compared to the widely used approach of the short-term average/long-term average of signal characteristic function(STA/LTA),DPick is 1.6 times less likely to detect noise as a P-wave,and 76 times less likely to miss P-waves.In terms of estimating P-wave arrival time,when the detection task is completed within 1 s,DPick′s detection occurrence is 7.4 times that of STA/LTA in the 0.05 s error band,and 1.6 times when the error band is 0.10 s.This verified that the proposed method has the potential for wide applications in EEW.展开更多
Crack closure is one of the reasons inducing changes of P-wave velocity of rocks under compression.In this context,a method is proposed to investigate the relationships among P-wave velocity,pre-existing cracks,and co...Crack closure is one of the reasons inducing changes of P-wave velocity of rocks under compression.In this context,a method is proposed to investigate the relationships among P-wave velocity,pre-existing cracks,and confining pressure based on the discrete element method(DEM).Pre-existing open cracks inside the rocks are generated by the initial gap of the flat-joint model.The validity of the method is evaluated by comparing the P-wave velocity tested on a sandstone specimen with numerical result.As the crack size is determined by the diameter of particles,the effects of three factors,i.e.number,aspect ratio,and orientation of cracks on the P-wave velocity are discussed.The results show that P-wave velocity is controlled by the(i.e.number) of open micro-cracks,while the closure pressure is determined by the aspect ratio of crack.The reason accounting for the anisotropy of P-wave velocity is the difference in crack number in measurement paths.Both of the number and aspect ratio of cracks can affect the responses of P-wave velocity to the applied confining pressure.Under confining pressure,the number of open cracks inside rocks will dominate the lowest P-wave velocity,and the P-wave velocity of the rock containing narrower cracks is more sensitive to the confining pressure.In this sense,crack density is difficult to be back-calculated merely by P-wave velocity.The proposed method offers a means to analyze the effect of pre-existing cracks on P-wave velocity.展开更多
Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven method...Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven methods primarily use the limited frequency bandwidth information of seismic data and can invert P-wave impedance with high accuracy,but not high resolution.Conventional data-driven methods mainly employ the information from well-log data and can provide high-accuracy and highresolution P-wave impedance owing to the superior nonlinear curve fitting capacity of neural networks.However,these methods require a significant number of training samples,which are frequently insufficient.To obtain P-wave impedance with both high accuracy and high resolution,we propose a model-data-driven inversion method using Res Nets and the normalized zero-lag cross-correlation objective function which is effective for avoiding local minima and suppressing random noise.By using initial models and training samples,the proposed model-data-driven method can invert P-wave impedance with satisfactory accuracy and resolution.Tests on synthetic and field data demonstrate the proposed method’s efficacy and practicability.展开更多
As P-wave velocity is sensitive to the variations in coal reservoir parameters,it is possible to monitor the injected CO_(2)through P-wave velocity during CO_(2)sequestration in coal.However,the effects of CO_(2)on th...As P-wave velocity is sensitive to the variations in coal reservoir parameters,it is possible to monitor the injected CO_(2)through P-wave velocity during CO_(2)sequestration in coal.However,the effects of CO_(2)on the coal P-wave velocity under triaxial stress are not clearly discerned.In the present study,different boundary conditions and gases were utilised to investigate the factors affecting the P-wave velocity after the interaction of coal with CO_(2).Experiments with helium indicated that the pore pressure primarily affected the P-wave velocity by altering the effective stress.Experiments with CH4 and CO_(2)indicated that matrix swelling induced-cleats porosity decline significantly promoted P-wave velocity.Moreover,CO_(2)caused a wider scale and severe weakening of coal matrix than CH4,thereby significantly decreasing the P-wave velocity,and the decline in P-wave velocity increases with vitrinite content.Furthermore,experiments under different boundary conditions showed that with the boundary condition having more constraints,the decrement of pore pressure on P-wave velocity is more weaken,whereas the improvement of matrix swelling on P-wave velocity is more evident.This study contributes to understanding the mechanism of effect of CO_(2)on P-wave velocity under triaxial stress condition and provides guidance for monitoring CO_(2)sequestration in coal.展开更多
If the thickness of coal seams and the lithology of both roofs and floors of coal seams have not changed at all or only a little, then it is thought that the elastic anisotropy of coal seams depends mainly on fracture...If the thickness of coal seams and the lithology of both roofs and floors of coal seams have not changed at all or only a little, then it is thought that the elastic anisotropy of coal seams depends mainly on fractures and obeys the horizontally symmetric model of an azimuth anisotropy. For a fixed offset, the amplitude A of the reflection P-wave and the cosine of 2φ has an approximately linear relation, (φ is the source-detector azimuth with respect to the fracture strike. Based on this relationship, many things can be done, such as the extraction of macro bins, the correction of residual normal moveout, the formation of azimuth gather, the transformation and normalization of azimuth gathers and the extraction of reflection wave amplitudes of coal seams. The least squares method was used to inverse theoretically the direction and density of fractures of coal seams. The result is in good agreement with the regional geological structure, indicating that the azimuth anisotropic analysis of the P-wave is feasible in evaluating the density and direction of fractures in coal seams.展开更多
Laboratory tests were performed on Toyoura sand specimens to investigate the relationship between degree of saturation Sr, B-value and P-wave velocity Vp. Different types of pore water （de-aired water or tap water） ...Laboratory tests were performed on Toyoura sand specimens to investigate the relationship between degree of saturation Sr, B-value and P-wave velocity Vp. Different types of pore water （de-aired water or tap water） and pore gas （air or CO2） as well as different magnitudes of back pressure were used to achieve different Sr （or B-value）. The measured relationship between B-value and Vp was not consistent with the theoretical prediction. The measurement shows that the Vp value in the specimen flushed with de-aired water is independent of B-value （or St） and is always around the one in fully saturated condition. However, the Vp value in the specimen flushed with tap water increases with B-value, but the shape of the relationship between Vp and B-value is quite different from the theoretical prediction. The possible explanation for the discrepancy between laboratory measurement and theoretical prediction lies in that the air exists in the water as air bubbles and therefore the pore fluid （air-water mixture） is heterogeneous instead of homogenous assumed in the theoretical prediction.展开更多
Sinkholes and cavities can represent serious hazards to human safety and urban infrastructures,cause roadbed subsidence,and so on.It is therefore essential to evaluate various sinkholes in different depths and sizes t...Sinkholes and cavities can represent serious hazards to human safety and urban infrastructures,cause roadbed subsidence,and so on.It is therefore essential to evaluate various sinkholes in different depths and sizes to assess the risk of collapse.This paper evaluates the effect of different cavities on Rayleigh-wave propagation and body(P)wave refraction.Rayleigh(R)wave propagation is analyzed according to the classical multi-channel analysis of surface waves(MASW)method also considering the R-wave backscattering.Synthetic seismic traces are computed by means of finite element modeling(FEM)for cavity and intrusion at different depths and sizes.Furthermore,field acquisition data is used to verify the detection effect of a cavity on R-and P-waves.The results show that the presence of backscattered R-waves and the changes in the R-wave velocity spectrum can help in cavity identification.Additional possible evidence is represented by significant changes in the refraction travel times over the cavity location.It can be concluded that the field data are in good agreement with the synthetic,and it could be effective to consider the results of both R-and P-wave analysis in order to efficiently identify the cavities.展开更多
It is essential to study the porosity,thermal conductivity,and P-wave velocity of calcarenites,as well as the anisotropy coefficients of the thermal conductivity and P-wave velocity,for civil engineering,and conservat...It is essential to study the porosity,thermal conductivity,and P-wave velocity of calcarenites,as well as the anisotropy coefficients of the thermal conductivity and P-wave velocity,for civil engineering,and conservation and restoration of historical monuments.This study focuses on measuring the thermal conductivity using the thermal conductivity scanning(TCS)technique and measuring the P-wave ve-locity using portable equipment.This was applied for some dry and saturated calcarenite samples in the horizontal and vertical directions(parallel and perpendicular to the bedding plane,respectively).The calcarenites were selected from some historical monuments in Morocco.These physical properties were measured in the laboratory to find a reliable relationship between all of these properties.As a result of the statistical analysis of the obtained data,excellent linear relationships were observed between the porosity and both the thermal conductivity and porosity.These relationships are characterized by relatively high coefficients of determination for the horizontal and vertical samples.Based on the thermal conductivity and P-wave velocity values in these two directions,the anisotropy coefficients of these two properties were calculated.The internal structure and the pore fabric of the calcarenite samples were delineated using scanning electron microscopy(SEM),while their chemical and mineral compositions were studied using the energy dispersive X-ray analysis(EDXA)and X-ray diffraction(XRD)techniques.展开更多
A magnitude 5.5 earthquakes occurred in Eryuan County,Dali Bai Autonomous Prefecture,Yunnan Province,China,on March 3.And a magnitude 5.0 earthquake occurred in the same place on April 17,2013,i.e.,45 days later.Then,...A magnitude 5.5 earthquakes occurred in Eryuan County,Dali Bai Autonomous Prefecture,Yunnan Province,China,on March 3.And a magnitude 5.0 earthquake occurred in the same place on April 17,2013,i.e.,45 days later.Then,on May 21,2021,multiple earthquakes,one with magnitude 6.4 and several at 5.0 or above,occurred in Yangbi County,Dali Bai Autonomous Prefecture,Yunnan Province,China.All of these occurred in the Weixi-QiaohouWeishan fault zone.In this study,1,874 seismic events in Yangbi and Eryuan counties were identified by automatic micro-seismic identification technology and the first arrivals were picked up manually.Following this,a total of 11,968 direct P-wave absolute arrivals and 73,987 high-quality Pwave relative arrivals were collected for joint inversion via the double difference tomography method.This was done to obtain the regional three-dimensional fine crustal P-wave velocity structure.The results show that the travel time residuals before and after inversion decreased from the initial–0.1–0.1 s to–0.06–0.06 s.The upper crust in the study area,which exhibited a low-velocity anomaly,corresponded to the basin region;this indicated that the low-velocity anomaly in the shallow part of the study area was affected by the basin.Results also showed some correlation between the distribution of the earthquakes and velocity structure,as there was a lowvelocity body Lv1 with a wide distribution at depths ranging from 15–20 km in the Yangbi and Eryuan earthquake regions.In addition,earthquakes occurred predominantly in the highlow velocity abnormal transition zone.The low-velocity body in the middle and lower crust may be prone to concentrating upper crustal stress,thus leading to the occurrence of earthquakes.展开更多
Wave velocities in haloanhydrites are difficult to determine and significantly depend on the mineralogy. We used petrophysical parameters to study the wave velocity in haloanhydrites in the Amur Darya Basin and constr...Wave velocities in haloanhydrites are difficult to determine and significantly depend on the mineralogy. We used petrophysical parameters to study the wave velocity in haloanhydrites in the Amur Darya Basin and constructed a template of the relation between haloanhydrite mineralogy （anhydrite, salt, mudstone, and pore water） and wave velocities. We used the relation between the P-wave rnoduli ratio and porosity as constraint and constructed a graphical model （petrophysical template） for the relation between wave velocity, mineral content and porosity. We tested the graphical model using rock core and well logging data.展开更多
Conventional land vertical seismic profiling （VSP） exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite ...Conventional land vertical seismic profiling （VSP） exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite borehole shots and vertical vibrations from controllable seismic sources at the surface will produce relatively strong pure P-waves and weaker pure S-waves. Interfaces with a large Poisson＇s ratio difference have a positive influence on the formation of strong transmitted converted S-waves. By a comparative analysis of pure S-waves from sources and converted downgoing S-waves, we believe that the main frequency of pure S-waves is usually lower than pure P-waves while the main frequency of downgoing converted S-waves is close to that of P-waves. We have studied zero-offset and offset VSP data from land P-wave sources. Results show that pure S-waves commonly exist in these data with differences in wave intensity. S-wave velocity can be obtained from the P-wave source zero-offset VSP data. Finally, we discuss the bright future of joint application of VSP P-and S-waves and the full use of S-waves in P-wave source VSP data.展开更多
文摘Scattering and Diffraction of elastic in-plane P-and SV-waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and strong-motion seismologists for over forty years. The case of out-of-plane SH waves on the same elastic canyon that is semi-circular in shape on the half-space surface is the first such problem that was solved by analytic closed form solutions over forty years ago by Trifunac. The corresponding case of in-plane P-and SV-waves on the same circular canyon is a much more complicated problem because, the in-plane P-and SV-scattered waves have different wave speeds and together they must have zero normal and shear stresses at the half-space surface. It is not until recently in 2014 that analytic solution for such problem is found by the author in the work of Lee and Liu. This paper uses the technique of Lee and Liu of defining these stress-free scattered waves to solve the problem of the scattered and diffraction of these in-plane waves on an almost-circular surface canyon that is arbitrary in shape.
基金supported by the National Science Council of Taiwan granted NSC-94-2816-M-194-004
文摘Based on the measurement of the arrival time of maxima magnitude from band-pass filtering signals which were determined using a new Morlet wavelet multiple-filter method, we develop a method for measuring intrinsic and attenuative dispersion of the first cycle direct P-wave. We determine relative group delays of spectral components of direct P-waves for 984 ray paths from SML and ALS stations of the Taiwan Central Weather Bureau Seismic Network （CWBSN）. Using continuous relaxation model, we deduce a new transfer function that relates intrinsic dispersion to attenuation. Based on the genetic algorithm （GA）, we put forward a new inversion procedure for determining which is defined the flat part of quality factor Q（ω） spectrum, τ1 and τ2 parameters. The results indicate that ① The distribution of Om values versus epicentral distance and depth show that Qm values linearly increase with increasing of epicentral distance and depth, and Qm values is clearly independent of earthquakes magnitude; ② In the different depth ranges, Qm residual show no correlation with variations in epicentral distance. Some significant changes of Qm residual with time is likely caused by pre-seismic stress accumulation, and associated with fluid-filled higher density fractures rock volume in the source area of 1999 Chi-Chi Taiwan earthquake. We confirm that Qm residual with time anomaly appears about 2.5 years before the Chi-Chi earthquake; ③ A comparison of Qm residual for different depth range between SML and ALS stations show that the level of stress has vertical and lateral difference; ④ The area near observation station with both anomalously increasing and decreasing averaged Qm residual is likely an unstable environment for future strong earthquake occurrence. This study demonstrates the capability of direct P-waves dispersion for monitoring attenuation characteristics and its state changes of anelastic medium of the Earth at short propagation distance using seismograms recorded from very small events.
基金Supported by National Natural Science Foundation of China (No. 50978156 and No. 50908183)
文摘The indirect boundary element method (IBEM) is used to study the surface motion of an alluvial valley in layered half-space for incident plane P-waves based on Wolf’s theory. Firstly, the free field response can be solved by the direct stiffness method, and the scattering wave response is calculated by Green’s functions of distributed loads acting on inclined lines in a layered half-space. The method is verified by comparing its results with literature and numerical analyses are performed by taking the amplification of incident plane P-waves by an alluvial valley in one soil layer resting on bedrock as an example. The results show that there exist distinct differences between the wave amplification by an alluvial valley embedded in layered half-space and that in homogeneous half-space and there is interaction between the valley and the soil layer. The amplitudes are relatively large when incident frequencies are close to the soil layer’s resonant frequencies.
基金supported by the National Natural Science Foundation of China (Grant No. 90814013)
文摘Finite-frequency travel time tomography is a newly developing method.The main procedure in this new method is to compute the traveltime sensitive kernel.The travel time of the same scatterer needs to be used for computing the traveltime sensitive kernel many times.It is a time-consuming task.It is easy and fast to get the travel time from analytic equations in a simple model such as a homogenous or linear velocity media.However,most of the earth models are layered.It is cumbersome to get the travel time from analytic equations.In order to enhance the computation efficiency,we used the table look-up method to compute the finite-frequency travel time sensitive kernel for P-waves in a layered structure model.We chose the AK135 earth model for the velocity model.The table look-up method saved about 50% of the computation time.We enhanced the computation speed by using the table lookup method in the same velocity model,which was very useful for enhancing the computation efficiency for the finite-frequency travel time tomography.
基金supported by the National Natural Science Foundation of China(Nos.42574119,42274083,41974049)partly supported by the Urban Geological Survey Project of Linyi,Shandong Province,China(No.SDGP371300202102000468).
文摘A P-wave velocity model was built in the central southern of the Tanlu Fault based on double-difference tomography.The results suggest the presence of a low-velocity anomaly extending from the surface to a depth of 25 km around the Tanlu and Feixi Faults,representing fault-related fluids caused by partial melting.The relocated earthquakes indicate a significant concentration of seismic activity above 20 km around the Tanlu and Feixi Faults,suggesting that prominent fault systems possibly serve as conduits for the upward migration of deep minerals.The proposed geodynamic model,supported by geological and geophysical data,suggests that the migration of deep mineralized materials extends along the Tanlu Fault.The obtained results serve as a crucial foundation for elucidating the intricate process of mineralization in the central southern segment of the Tanlu Fault,thereby enhancing comprehension regarding the interaction among ore body formation,fault fluids,localized melting,and seismic activity.
基金supported by the Major National Project Program (No.2011ZX05007-006)
文摘We experimentally studied the effect of crack aperture on P-wave velocity, amplitude, anisotropy and dispersion. Experimental models were constructed based on Hudson＇s theory. Six crack models were embedded with equal-radius penny-shaped crack inclusions in each layer. The P-wave velocity and amplitude were measured parallel and perpendicular to the layers of cracks at frequencies of 0.1 MHz to 1 MHz. The experiments show that as the crack aperture increases from 0.l mm to 0.34 mm, the amplitude of the P-waves parallel to the crack layers decreases linearly with increasing frequency and the P-wave velocity dispersion varies from 1.5% to 2.1%, whereas the amplitude of the P-wave perpendicular to the crack layers decreases quadratically with increasing frequency and the velocity dispersion varies from 1.9% to 4.7%. The variation in the velocity dispersion parallel and perpendicular to the cracks intensifies the anisotropy dispersion of the P-waves in the crack models （6.7% to 83%）. The P-wave dispersion strongly depends on the scattering characteristics of the crack apertures.
基金supported by NSF grant EAR-063566(F.N.)National Natural Science Foundation of China grant 40774042(J.L.)
文摘The recently built China Digital Seismic Network consists of the China National Digital Seismic Network （CNDSN）, 31 regional seismic networks and several small aperture arrays with more than 1 000 stations including 850＋ broadband stations. It forms a gigantic seismic array that provides an unprecedented opportunity to study the Earth＇s deep interior besides its routine task of seismic monitoring. Many modern seismic studies rely on rotation of vertical and horizontal components in order to separate different types of seismic waves. Knowledge of the orientations of the two horizontal components thus is important to perform a correction rotation. We analyzed particle motions of teleseismic P waves recorded by the network and used them to estimate the northcomponent azimuth of each station. An SNR-weighted-multi-event method was introduced to obtain component azimuths that best explain the P-wave particle motions of all the events recorded at a station. The method provides robust estimates including a measurement error calculated from background noise levels. We found that about one third of the stations have some sort of problems, including misorientation of the two horizontal components, mislabeling and polarity reversal in one or more components. These problems need to be taken into account for any rotation based seismic studies.
基金This work was supported by the National Natural Science Foundation of China(Grant No.41772333)the program of State Key Laboratory of Frozen Soil Engineering(Grant No.SKLFSE201713)the Shaanxi Province New-Star Talents Promotion Project of Science and Technology(Grant No.2019KJXX-049).
文摘It is generally accepted that the uniaxial compressive strength(UCS)and P-wave velocity of rocks tend to decrease simultaneously with increasing temperature.However,based on a great number of statistical data and systematic analysis of the microstructure variation of rocks with temperature rising and corresponding propagation mechanism of elastic wave,the results show that(1)There are three different trends for the changes of UCS and P-wave velocity of sandstone when heated from room temperature(20C or 25C)to 800C:(i)Both the UCS and P-wave velocity decrease simultaneously;(ii)The UCS increases initially and then decreases,while the P-wave velocity decreases continuously;and(iii)The UCS increases initially and then fluctuates,while the P-wave velocity continuously decreases.(2)The UCS changes at room temperaturee400C,400Ce600C,and 600Ce800C are mainly attributed to the discrepancy of microstructure characteristics and quartz content,the transformation plasticity of clay minerals,and the balance between the thermal cementation and thermal damage,respectively.(3)The inconsistency in the trends of UCS and P-wave velocity changes is caused by the change of quartz content,phase transition of water and certain minerals.
文摘Engineering rock mass classification,based on empirical relations between rock mass parameters and engineering applications,is commonly used in rock engineering and forms the basis for designing rock structures.The basic data required may be obtained from visual observation and laboratory or field tests.However,owing to the discontinuous and variable nature of rock masses,it is difficult for rock engineers to directly obtain the specific design parameters needed.As an alternative,the use of geophysical methods in geomechanics such as seismography may largely address this problem.In this study,25 seismic profiles with the total length of 543 m have been scanned to determine the geomechanical properties of the rock mass in blocks Ⅰ,Ⅲ and Ⅳ-2 of the Choghart iron mine.Moreover,rock joint measurements and sampling for laboratory tests were conducted.The results show that the rock mass rating（RMR） and Q values have a close relation with P-wave velocity parameters,including P-wave velocity in field(V;).P-wave velocity in the laboratory(V;) and the ratio of V;V;(i.e.K;= V;/V;.However,Q value,totally,has greater correlation coefficient and less error than the RMR,In addition,rock mass parameters including rock quality designation（RQD）,uniaxial compressive strength（UCS）,joint roughness coefficient（JRC） and Schmidt number（RN） show close relationship with P-wave velocity.An equation based on these parameters was obtained to estimate the P-wave velocity in the rock mass with a correlation coefficient of 91%.The velocities in two orthogonal directions and the results of joint study show that the wave velocity anisotropy in rock mass may be used as an efficient tool to assess the strong and weak directions in rock mass.
基金National Natural Science Foundation of China under Grant Nos.51968016 and 5197083806the Guangxi Innovation Driven Development Project(Science and Technology Major Project,Grant No.Guike AA18118008).
文摘Fast and accurate P-wave arrival picking significantly affects the performance of earthquake early warning(EEW)systems.Automated P-wave picking algorithms used in EEW have encountered problems of falsely picking up noise,missing P-waves and inaccurate P-wave arrival estimation.To address these issues,an automatic algorithm based on the convolution neural network(DPick)was developed,and trained with a moderate number of data sets of 17,717 accelerograms.Compared to the widely used approach of the short-term average/long-term average of signal characteristic function(STA/LTA),DPick is 1.6 times less likely to detect noise as a P-wave,and 76 times less likely to miss P-waves.In terms of estimating P-wave arrival time,when the detection task is completed within 1 s,DPick′s detection occurrence is 7.4 times that of STA/LTA in the 0.05 s error band,and 1.6 times when the error band is 0.10 s.This verified that the proposed method has the potential for wide applications in EEW.
基金partially supported by the National Natural Science Foundation of China(Grant Nos.41872210 and 51674238)the Open Research Fund of the State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.Z017006)。
文摘Crack closure is one of the reasons inducing changes of P-wave velocity of rocks under compression.In this context,a method is proposed to investigate the relationships among P-wave velocity,pre-existing cracks,and confining pressure based on the discrete element method(DEM).Pre-existing open cracks inside the rocks are generated by the initial gap of the flat-joint model.The validity of the method is evaluated by comparing the P-wave velocity tested on a sandstone specimen with numerical result.As the crack size is determined by the diameter of particles,the effects of three factors,i.e.number,aspect ratio,and orientation of cracks on the P-wave velocity are discussed.The results show that P-wave velocity is controlled by the(i.e.number) of open micro-cracks,while the closure pressure is determined by the aspect ratio of crack.The reason accounting for the anisotropy of P-wave velocity is the difference in crack number in measurement paths.Both of the number and aspect ratio of cracks can affect the responses of P-wave velocity to the applied confining pressure.Under confining pressure,the number of open cracks inside rocks will dominate the lowest P-wave velocity,and the P-wave velocity of the rock containing narrower cracks is more sensitive to the confining pressure.In this sense,crack density is difficult to be back-calculated merely by P-wave velocity.The proposed method offers a means to analyze the effect of pre-existing cracks on P-wave velocity.
基金financially supported by the Important National Science&Technology Specific Project of China(Grant No.2017ZX05018-005)
文摘Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven methods primarily use the limited frequency bandwidth information of seismic data and can invert P-wave impedance with high accuracy,but not high resolution.Conventional data-driven methods mainly employ the information from well-log data and can provide high-accuracy and highresolution P-wave impedance owing to the superior nonlinear curve fitting capacity of neural networks.However,these methods require a significant number of training samples,which are frequently insufficient.To obtain P-wave impedance with both high accuracy and high resolution,we propose a model-data-driven inversion method using Res Nets and the normalized zero-lag cross-correlation objective function which is effective for avoiding local minima and suppressing random noise.By using initial models and training samples,the proposed model-data-driven method can invert P-wave impedance with satisfactory accuracy and resolution.Tests on synthetic and field data demonstrate the proposed method’s efficacy and practicability.
基金supported by the National Natural Science Foundation of China(No.51974304)the Natural Science Foundation of Hebei Province(No.E2020402075)+2 种基金the 2nd Xplorer Prize sponsored by the Tencent Foundationthe Program for Changjiang Scholars and Innovative Research Team in University(No.IRT 17R103)the Qinglan Project of Jiangsu Province.
文摘As P-wave velocity is sensitive to the variations in coal reservoir parameters,it is possible to monitor the injected CO_(2)through P-wave velocity during CO_(2)sequestration in coal.However,the effects of CO_(2)on the coal P-wave velocity under triaxial stress are not clearly discerned.In the present study,different boundary conditions and gases were utilised to investigate the factors affecting the P-wave velocity after the interaction of coal with CO_(2).Experiments with helium indicated that the pore pressure primarily affected the P-wave velocity by altering the effective stress.Experiments with CH4 and CO_(2)indicated that matrix swelling induced-cleats porosity decline significantly promoted P-wave velocity.Moreover,CO_(2)caused a wider scale and severe weakening of coal matrix than CH4,thereby significantly decreasing the P-wave velocity,and the decline in P-wave velocity increases with vitrinite content.Furthermore,experiments under different boundary conditions showed that with the boundary condition having more constraints,the decrement of pore pressure on P-wave velocity is more weaken,whereas the improvement of matrix swelling on P-wave velocity is more evident.This study contributes to understanding the mechanism of effect of CO_(2)on P-wave velocity under triaxial stress condition and provides guidance for monitoring CO_(2)sequestration in coal.
基金Projects 40574058 supported by the National Natural Science Foundation of China, 2005cb221500 the National Key Basic Research and Development(973) Program of China and 03(2007) the Scientific and Technological Project about Geology and Mineral Resources of Henan Land Resources Department
文摘If the thickness of coal seams and the lithology of both roofs and floors of coal seams have not changed at all or only a little, then it is thought that the elastic anisotropy of coal seams depends mainly on fractures and obeys the horizontally symmetric model of an azimuth anisotropy. For a fixed offset, the amplitude A of the reflection P-wave and the cosine of 2φ has an approximately linear relation, (φ is the source-detector azimuth with respect to the fracture strike. Based on this relationship, many things can be done, such as the extraction of macro bins, the correction of residual normal moveout, the formation of azimuth gather, the transformation and normalization of azimuth gathers and the extraction of reflection wave amplitudes of coal seams. The least squares method was used to inverse theoretically the direction and density of fractures of coal seams. The result is in good agreement with the regional geological structure, indicating that the azimuth anisotropic analysis of the P-wave is feasible in evaluating the density and direction of fractures in coal seams.
基金Foundation item： Project（2012CB719803） supported by the National Basic Research Program of China Project（201011159098） supported by the Seed Funding for Basic Research Scheme from The University of Hong Kong, China
文摘Laboratory tests were performed on Toyoura sand specimens to investigate the relationship between degree of saturation Sr, B-value and P-wave velocity Vp. Different types of pore water （de-aired water or tap water） and pore gas （air or CO2） as well as different magnitudes of back pressure were used to achieve different Sr （or B-value）. The measured relationship between B-value and Vp was not consistent with the theoretical prediction. The measurement shows that the Vp value in the specimen flushed with de-aired water is independent of B-value （or St） and is always around the one in fully saturated condition. However, the Vp value in the specimen flushed with tap water increases with B-value, but the shape of the relationship between Vp and B-value is quite different from the theoretical prediction. The possible explanation for the discrepancy between laboratory measurement and theoretical prediction lies in that the air exists in the water as air bubbles and therefore the pore fluid （air-water mixture） is heterogeneous instead of homogenous assumed in the theoretical prediction.
文摘Sinkholes and cavities can represent serious hazards to human safety and urban infrastructures,cause roadbed subsidence,and so on.It is therefore essential to evaluate various sinkholes in different depths and sizes to assess the risk of collapse.This paper evaluates the effect of different cavities on Rayleigh-wave propagation and body(P)wave refraction.Rayleigh(R)wave propagation is analyzed according to the classical multi-channel analysis of surface waves(MASW)method also considering the R-wave backscattering.Synthetic seismic traces are computed by means of finite element modeling(FEM)for cavity and intrusion at different depths and sizes.Furthermore,field acquisition data is used to verify the detection effect of a cavity on R-and P-waves.The results show that the presence of backscattered R-waves and the changes in the R-wave velocity spectrum can help in cavity identification.Additional possible evidence is represented by significant changes in the refraction travel times over the cavity location.It can be concluded that the field data are in good agreement with the synthetic,and it could be effective to consider the results of both R-and P-wave analysis in order to efficiently identify the cavities.
文摘It is essential to study the porosity,thermal conductivity,and P-wave velocity of calcarenites,as well as the anisotropy coefficients of the thermal conductivity and P-wave velocity,for civil engineering,and conservation and restoration of historical monuments.This study focuses on measuring the thermal conductivity using the thermal conductivity scanning(TCS)technique and measuring the P-wave ve-locity using portable equipment.This was applied for some dry and saturated calcarenite samples in the horizontal and vertical directions(parallel and perpendicular to the bedding plane,respectively).The calcarenites were selected from some historical monuments in Morocco.These physical properties were measured in the laboratory to find a reliable relationship between all of these properties.As a result of the statistical analysis of the obtained data,excellent linear relationships were observed between the porosity and both the thermal conductivity and porosity.These relationships are characterized by relatively high coefficients of determination for the horizontal and vertical samples.Based on the thermal conductivity and P-wave velocity values in these two directions,the anisotropy coefficients of these two properties were calculated.The internal structure and the pore fabric of the calcarenite samples were delineated using scanning electron microscopy(SEM),while their chemical and mineral compositions were studied using the energy dispersive X-ray analysis(EDXA)and X-ray diffraction(XRD)techniques.
基金funded by the general project of National Natural Science Foundation of China(No.41774072).
文摘A magnitude 5.5 earthquakes occurred in Eryuan County,Dali Bai Autonomous Prefecture,Yunnan Province,China,on March 3.And a magnitude 5.0 earthquake occurred in the same place on April 17,2013,i.e.,45 days later.Then,on May 21,2021,multiple earthquakes,one with magnitude 6.4 and several at 5.0 or above,occurred in Yangbi County,Dali Bai Autonomous Prefecture,Yunnan Province,China.All of these occurred in the Weixi-QiaohouWeishan fault zone.In this study,1,874 seismic events in Yangbi and Eryuan counties were identified by automatic micro-seismic identification technology and the first arrivals were picked up manually.Following this,a total of 11,968 direct P-wave absolute arrivals and 73,987 high-quality Pwave relative arrivals were collected for joint inversion via the double difference tomography method.This was done to obtain the regional three-dimensional fine crustal P-wave velocity structure.The results show that the travel time residuals before and after inversion decreased from the initial–0.1–0.1 s to–0.06–0.06 s.The upper crust in the study area,which exhibited a low-velocity anomaly,corresponded to the basin region;this indicated that the low-velocity anomaly in the shallow part of the study area was affected by the basin.Results also showed some correlation between the distribution of the earthquakes and velocity structure,as there was a lowvelocity body Lv1 with a wide distribution at depths ranging from 15–20 km in the Yangbi and Eryuan earthquake regions.In addition,earthquakes occurred predominantly in the highlow velocity abnormal transition zone.The low-velocity body in the middle and lower crust may be prone to concentrating upper crustal stress,thus leading to the occurrence of earthquakes.
基金supported by the National Major Scientific and Technological Special Project(No.2011ZX05029-003)the project of the Research Institute of Petroleum Exploration&Development(No.2012Y-058)
文摘Wave velocities in haloanhydrites are difficult to determine and significantly depend on the mineralogy. We used petrophysical parameters to study the wave velocity in haloanhydrites in the Amur Darya Basin and constructed a template of the relation between haloanhydrite mineralogy （anhydrite, salt, mudstone, and pore water） and wave velocities. We used the relation between the P-wave rnoduli ratio and porosity as constraint and constructed a graphical model （petrophysical template） for the relation between wave velocity, mineral content and porosity. We tested the graphical model using rock core and well logging data.
文摘Conventional land vertical seismic profiling （VSP） exploration usually uses P-wave sources and three-component geophones for receivers, emphasizing P- and converted S-waves. Previous studies show that both dynamite borehole shots and vertical vibrations from controllable seismic sources at the surface will produce relatively strong pure P-waves and weaker pure S-waves. Interfaces with a large Poisson＇s ratio difference have a positive influence on the formation of strong transmitted converted S-waves. By a comparative analysis of pure S-waves from sources and converted downgoing S-waves, we believe that the main frequency of pure S-waves is usually lower than pure P-waves while the main frequency of downgoing converted S-waves is close to that of P-waves. We have studied zero-offset and offset VSP data from land P-wave sources. Results show that pure S-waves commonly exist in these data with differences in wave intensity. S-wave velocity can be obtained from the P-wave source zero-offset VSP data. Finally, we discuss the bright future of joint application of VSP P-and S-waves and the full use of S-waves in P-wave source VSP data.
