Transient S wave velocity rupture (TSVR) means the velocity of fault rupture propagation is between S wave velocity α and P wave velocity β . Its existing in the rupture of in plane ( i.e . strike slip...Transient S wave velocity rupture (TSVR) means the velocity of fault rupture propagation is between S wave velocity α and P wave velocity β . Its existing in the rupture of in plane ( i.e . strike slip) fault has been proved, but in 2 dimensional classical model, there are two difficulties in transient S wave velocity rupture, i.e ., initialization difficulty and divergence difficulty in interpreting the realization of TSVR. The initialization difficulty means, when v ↑ v R (Rayleigh wave velocity), the dynamic stress strength factor K 2(t) →+0, and changes from positive into negative in the interval ( v R, β ). How v transit the forbidden of ( v R, β )? The divergence difficulty means K 2(t) →+∞ when v ↓ β . Here we introduce the concept of fractal and tunnel effect that exist everywhere in fault. The structure of all the faults is fractal with multiple cracks. The velocity of fault rupture is differentiate of the length of the fault respect to time, so the rupture velocity is also fractal. The tunnel effect means the dynamic rupture crosses over the interval of the cracks, and the coalescence of the intervals is slower than the propagation of disturbance. Suppose the area of earthquake nucleation is critical or sub critical propagation everywhere, the arriving of disturbance triggers or accelerates the propagation of cracks tip at once, and the observation system cannot distinguish the front of disturbance and the tip of fracture. Then the speed of disturbance may be identified as fracture velocity, and the phenomenon of TSVR appears, which is an apparent velocity. The real reason of apparent velocity is that the mathematics model of shear rupture is simplified of complex process originally. The dual character of rupture velocity means that the apparent velocity of fault and the real velocity of micro crack extending, which are different in physics, but are unified in rupture criterion. Introducing the above mentioned concept to the calculation of K 2 (t) , the difficulty of initialization can be overcome, and the integral equation of triggering the initialization of TSVR is given quantitatively. By solving this integral equation, the lower limit of TSVR is 1.105 3 β , not β , and the divergence difficulty is overcome. TSVR is unstable solution, and may degenerate to sub Rayleigh wave velocity rupture immediately where the non critical condition can be measured. The results of this paper show that the initialization and continuum depends on the condition of earthquake nucleation in seismogenic area.展开更多
Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the ...Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the crust is thick in the studied region, reaching 40 km in thickness under 4 stations. The crust all alternatives high velocity layer with low velocity one. There appear varied velocity structures for different stations, and the stations around the same tectonic region exhibit similar structure characteristics. Combined with dominant depth distribution of many small-moderate earthquakes, the correlation between seismogenic layers and crustal structures of high and low velocity layers has been discussed.展开更多
Tengchong volcanic area is located near the impinging and underthrust margin of India and Eurasia plates. The volcanic activity is closely related to the tectonic environment. The deep structure characteristics are in...Tengchong volcanic area is located near the impinging and underthrust margin of India and Eurasia plates. The volcanic activity is closely related to the tectonic environment. The deep structure characteristics are inferred from the receiver function inversion with the teleseismic records in the paper. The results show that the low velocity zone is influenced by the NE-trending Dayingjiang fault. The S-wave low velocity structure occurs obviously in the southern part of the fault, but unobviously in its northern part. There are low velocity zones in the shallow po-sition, which coincides with the seismicity. It also demonstrates that the low velocity zone is directly related to the thermal activity in the volcanic area. Therefore, we consider that the volcano may be alive again.展开更多
This article presents a case study concerning a seismic characterization project.Full-wave sonic logging was used to characterize the shallow compressional wave and shear wave velocity profiles in the site.Anomalous v...This article presents a case study concerning a seismic characterization project.Full-wave sonic logging was used to characterize the shallow compressional wave and shear wave velocity profiles in the site.Anomalous values of the Poisson’s ratio derived from the velocity profiles suggested that the boreholes might have traversed slow formations(i.e.with shear wave velocity smaller than the borehole fluid compressional wave velocity or“mud-wave speed”)and that conventional processing of the sonic logs might have misinterpreted the direct arrivals of fluid acoustic waves as arrivals caused by shear wave propagation in the rock.Consequently,the shear wave velocity profiles provided by the contractor were considered to be unreliable by the project team.To address these problems,a non-conventional determination of the shear wave velocity was implemented,based on the relationship between the Poisson’s ratio of the rock formation and the shape of the first train of sonic waves which arrived to the receivers in the sonic probe.The relationship was determined based on several hundreds of finite element simulations of the acoustic wave propagation in boreholes with the same diameter as used in the perforations.The present article describes how this non-conventional approach was developed and implemented to obtain the shear wave velocity profiles from the raw sonic logs.The approach allows an extension of the range of applicability of full-wave sonic logging to determination of shear wave velocity profiles in formations with low compressional wave velocities.The method could be used to obtain shear wave velocity profiles where compressional wave velocity is as low as slightly larger than the mud-wave speed.A sample sonic log in Log ASCII Standard(LAS)format is provided as supplementary material to this paper via Mendeley Data,together with the FORTRAN source code used to process the log following the approach described in this study.展开更多
To deal with over-shooting and gouging in high speed machining, a novel approach for velocity smooth link is proposed. Considering discrete tool path, cubic spline curve fitting is used to find dangerous points, and a...To deal with over-shooting and gouging in high speed machining, a novel approach for velocity smooth link is proposed. Considering discrete tool path, cubic spline curve fitting is used to find dangerous points, and according to spatial geometric properties of tool path and the kinematics theory, maximum optimal velocities at dangerous points are obtained. Based on method of velocity control characteristics stored in control system, a fast algorithm for velocity smooth link is analyzed and formulated. On-line implementation results show that the proposed approach makes velocity changing more smoothly compared with traditional velocity control methods and improves productivity greatly.展开更多
Estimation of S-wave velocity using logging data has mainly been performed for sandstone,mudstone and oil and gas strata,while its application to hydrate reservoirs has been largely overlooked.In this paper we present...Estimation of S-wave velocity using logging data has mainly been performed for sandstone,mudstone and oil and gas strata,while its application to hydrate reservoirs has been largely overlooked.In this paper we present petrophysical methods to estimate the S-wave velocity of hydrate reservoirs with the P-wave velocity and the density as constraints.The three models used in this paper are an equivalent model(MBGL),a three-phase model(TPBE),and a thermo-elasticity model(TEM).The MBGL model can effectively describe the internal relationship among the components of the rock,and the estimated P-wave velocities are in good agreement with the measured data(2.8% error).However,in the TPBE model,the solid,liquid and gas phases are considered to be independent of each other,and the estimation results are relatively low(46.6% error).The TEM model is based on the sensitivity of the gas hydrate to temperature and pressure,and the accuracy of the estimation results is also high(3.6% error).Before the estimation,the occurrence patterns of hydrates in the Shenhu area were examined,and occurrence state one(the hydrate is in solid form in the reservoir) was selected for analysis.By using the known P-wave velocity and density as constraints,a reasonable S-wave velocity value(ranging from 400 to 1100 m s^(-1) and for a hydrate layer of 1100 m s^(-1)) can be obtained through multiple iterations.These methods and results provide new data and technical support for further research on hydrates and other geological features in the Shenhu area.展开更多
Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in ...Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in the interior of Ordos block is relatively stable, where no apparent change between high and low velocity layers exists and the shear wave velocity increases steadily with the depth. There is a 12km thick layer at the depth of 25km under this station, with an S wave velocity (V S=3.90km/s) lower than that at the same depth in its eastern and southern areas (V S≥4.00km/s). The crust under the eastern margin of Ordos block is thicker than that of the Yulin station, and the velocity structures alternate between the high and low velocity layers, with more low velocity layers. It has the same characteristic as having a 10km-thick low velocity layer (V S=3.80km/s) in the lower crust but buried at a depth of about 35km. Moreover, we studied the V P/V S ratio under each station in combination with the result of P wave velocity inversion. The results show that, the average velocity ratio of the Yulin station at the interior of Ordos block is only 1.68, with a very low ratio (about 1.60) in the upper crust and a stable ratio of about 1.73 in the mid and lower crust, which indicates the media under this station is homogenous and stable, being in a state of rigidity. But at the stations in the eastern and southern margins of the Ordos block, several layers of high velocity ratio (about 1.80) have been found, in which the average velocity ratio under Kelan and Lishi stations at the eastern margin is systemically higher than that of the general elastical body waves (1.732). This reflects that the crust under the marginal areas is more active relatively, and other materials may exist in these layers. Finally, we discussed the relationship among earthquakes, velocity structures beneath stations and faults.展开更多
In this article, we analyze the characters of SV-component receiver function of teleseismic body waves and its advantages in mapping the S-wave velocity structure of crust in detail. Similar to radial receiver functio...In this article, we analyze the characters of SV-component receiver function of teleseismic body waves and its advantages in mapping the S-wave velocity structure of crust in detail. Similar to radial receiver function, SV-component receiver function can be obtained by directly deconvolving the P-component from the SV-component of teleseismic recordings. Our analyses indicate that the change of amplitude of SV-component receiver function against the change of epicentral distance is less than that of radial receiver function. Moreover, the waveform of SV-component receiver function is simpler than the radial receiver function and gives prominence to the PS converted phases that are the most sensitive to the shear wave velocity structure in the inversion. The synthetic tests show that the convergence of SV-component receiver function inversion is faster than that of the radial receiver function inversion. As an example, we investigate the S-wave velocity structure beneath HIA sta-tion by using the SV-component receiver function inversion method.展开更多
Microtremors array observation for estimating S-wave velocity structure from phase velocities of Rayleigh and Love wave on two practical sites in Tangshan area by a China-US joint group are researched.The phase veloci...Microtremors array observation for estimating S-wave velocity structure from phase velocities of Rayleigh and Love wave on two practical sites in Tangshan area by a China-US joint group are researched.The phase velocities of Rayleigh wave are estimated from vertical component records and those of Love wave are estimated from three-component records of microtremors array using modified spatial auto-correlation method.Haskell matrix method is used in calculating Rayleigh and Love wave phase velocities,and the shallow S-wave velocity structure of two practical sites are estimated by means of a hybrid approach of Genetic Algorithm and Simplex.The results are compared with the PS logging data of the two sites,showing it is feasible to estimate the shallow S-wave velocity structure of practical site from the observation of microtremor array.展开更多
文摘Transient S wave velocity rupture (TSVR) means the velocity of fault rupture propagation is between S wave velocity α and P wave velocity β . Its existing in the rupture of in plane ( i.e . strike slip) fault has been proved, but in 2 dimensional classical model, there are two difficulties in transient S wave velocity rupture, i.e ., initialization difficulty and divergence difficulty in interpreting the realization of TSVR. The initialization difficulty means, when v ↑ v R (Rayleigh wave velocity), the dynamic stress strength factor K 2(t) →+0, and changes from positive into negative in the interval ( v R, β ). How v transit the forbidden of ( v R, β )? The divergence difficulty means K 2(t) →+∞ when v ↓ β . Here we introduce the concept of fractal and tunnel effect that exist everywhere in fault. The structure of all the faults is fractal with multiple cracks. The velocity of fault rupture is differentiate of the length of the fault respect to time, so the rupture velocity is also fractal. The tunnel effect means the dynamic rupture crosses over the interval of the cracks, and the coalescence of the intervals is slower than the propagation of disturbance. Suppose the area of earthquake nucleation is critical or sub critical propagation everywhere, the arriving of disturbance triggers or accelerates the propagation of cracks tip at once, and the observation system cannot distinguish the front of disturbance and the tip of fracture. Then the speed of disturbance may be identified as fracture velocity, and the phenomenon of TSVR appears, which is an apparent velocity. The real reason of apparent velocity is that the mathematics model of shear rupture is simplified of complex process originally. The dual character of rupture velocity means that the apparent velocity of fault and the real velocity of micro crack extending, which are different in physics, but are unified in rupture criterion. Introducing the above mentioned concept to the calculation of K 2 (t) , the difficulty of initialization can be overcome, and the integral equation of triggering the initialization of TSVR is given quantitatively. By solving this integral equation, the lower limit of TSVR is 1.105 3 β , not β , and the divergence difficulty is overcome. TSVR is unstable solution, and may degenerate to sub Rayleigh wave velocity rupture immediately where the non critical condition can be measured. The results of this paper show that the initialization and continuum depends on the condition of earthquake nucleation in seismogenic area.
基金funded by the National Natural Science Foundation of China (grant no.42074149)the Natural Science Foundation of Jiangsu Province (grant no.BK20201318)。
基金State Key Basic Development and Programming Project Mechanism and Prediction of Continental Strong Earthquakes (G1998040705).
文摘Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the crust is thick in the studied region, reaching 40 km in thickness under 4 stations. The crust all alternatives high velocity layer with low velocity one. There appear varied velocity structures for different stations, and the stations around the same tectonic region exhibit similar structure characteristics. Combined with dominant depth distribution of many small-moderate earthquakes, the correlation between seismogenic layers and crustal structures of high and low velocity layers has been discussed.
文摘Tengchong volcanic area is located near the impinging and underthrust margin of India and Eurasia plates. The volcanic activity is closely related to the tectonic environment. The deep structure characteristics are inferred from the receiver function inversion with the teleseismic records in the paper. The results show that the low velocity zone is influenced by the NE-trending Dayingjiang fault. The S-wave low velocity structure occurs obviously in the southern part of the fault, but unobviously in its northern part. There are low velocity zones in the shallow po-sition, which coincides with the seismicity. It also demonstrates that the low velocity zone is directly related to the thermal activity in the volcanic area. Therefore, we consider that the volcano may be alive again.
文摘This article presents a case study concerning a seismic characterization project.Full-wave sonic logging was used to characterize the shallow compressional wave and shear wave velocity profiles in the site.Anomalous values of the Poisson’s ratio derived from the velocity profiles suggested that the boreholes might have traversed slow formations(i.e.with shear wave velocity smaller than the borehole fluid compressional wave velocity or“mud-wave speed”)and that conventional processing of the sonic logs might have misinterpreted the direct arrivals of fluid acoustic waves as arrivals caused by shear wave propagation in the rock.Consequently,the shear wave velocity profiles provided by the contractor were considered to be unreliable by the project team.To address these problems,a non-conventional determination of the shear wave velocity was implemented,based on the relationship between the Poisson’s ratio of the rock formation and the shape of the first train of sonic waves which arrived to the receivers in the sonic probe.The relationship was determined based on several hundreds of finite element simulations of the acoustic wave propagation in boreholes with the same diameter as used in the perforations.The present article describes how this non-conventional approach was developed and implemented to obtain the shear wave velocity profiles from the raw sonic logs.The approach allows an extension of the range of applicability of full-wave sonic logging to determination of shear wave velocity profiles in formations with low compressional wave velocities.The method could be used to obtain shear wave velocity profiles where compressional wave velocity is as low as slightly larger than the mud-wave speed.A sample sonic log in Log ASCII Standard(LAS)format is provided as supplementary material to this paper via Mendeley Data,together with the FORTRAN source code used to process the log following the approach described in this study.
基金This project is supported by National Hi-tech Research and Development Program of China (863 Program, No. 2002AA421150)Specialized Re-search Fund for Doctor Program of Higher Education of China (No. 20030335091).
文摘To deal with over-shooting and gouging in high speed machining, a novel approach for velocity smooth link is proposed. Considering discrete tool path, cubic spline curve fitting is used to find dangerous points, and according to spatial geometric properties of tool path and the kinematics theory, maximum optimal velocities at dangerous points are obtained. Based on method of velocity control characteristics stored in control system, a fast algorithm for velocity smooth link is analyzed and formulated. On-line implementation results show that the proposed approach makes velocity changing more smoothly compared with traditional velocity control methods and improves productivity greatly.
基金supported by the National Natural Science Foundation of China(Nos.41304096 and 41176077)the National Science and Technology Major Project of China(No.2016ZX05024-001-002)+2 种基金the National High-tech R&D Program of China(863 ProgramNo.2013AA0925 01)the Fundamental Research Funds for the Central Universities(No.201762019)
文摘Estimation of S-wave velocity using logging data has mainly been performed for sandstone,mudstone and oil and gas strata,while its application to hydrate reservoirs has been largely overlooked.In this paper we present petrophysical methods to estimate the S-wave velocity of hydrate reservoirs with the P-wave velocity and the density as constraints.The three models used in this paper are an equivalent model(MBGL),a three-phase model(TPBE),and a thermo-elasticity model(TEM).The MBGL model can effectively describe the internal relationship among the components of the rock,and the estimated P-wave velocities are in good agreement with the measured data(2.8% error).However,in the TPBE model,the solid,liquid and gas phases are considered to be independent of each other,and the estimation results are relatively low(46.6% error).The TEM model is based on the sensitivity of the gas hydrate to temperature and pressure,and the accuracy of the estimation results is also high(3.6% error).Before the estimation,the occurrence patterns of hydrates in the Shenhu area were examined,and occurrence state one(the hydrate is in solid form in the reservoir) was selected for analysis.By using the known P-wave velocity and density as constraints,a reasonable S-wave velocity value(ranging from 400 to 1100 m s^(-1) and for a hydrate layer of 1100 m s^(-1)) can be obtained through multiple iterations.These methods and results provide new data and technical support for further research on hydrates and other geological features in the Shenhu area.
文摘Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in the interior of Ordos block is relatively stable, where no apparent change between high and low velocity layers exists and the shear wave velocity increases steadily with the depth. There is a 12km thick layer at the depth of 25km under this station, with an S wave velocity (V S=3.90km/s) lower than that at the same depth in its eastern and southern areas (V S≥4.00km/s). The crust under the eastern margin of Ordos block is thicker than that of the Yulin station, and the velocity structures alternate between the high and low velocity layers, with more low velocity layers. It has the same characteristic as having a 10km-thick low velocity layer (V S=3.80km/s) in the lower crust but buried at a depth of about 35km. Moreover, we studied the V P/V S ratio under each station in combination with the result of P wave velocity inversion. The results show that, the average velocity ratio of the Yulin station at the interior of Ordos block is only 1.68, with a very low ratio (about 1.60) in the upper crust and a stable ratio of about 1.73 in the mid and lower crust, which indicates the media under this station is homogenous and stable, being in a state of rigidity. But at the stations in the eastern and southern margins of the Ordos block, several layers of high velocity ratio (about 1.80) have been found, in which the average velocity ratio under Kelan and Lishi stations at the eastern margin is systemically higher than that of the general elastical body waves (1.732). This reflects that the crust under the marginal areas is more active relatively, and other materials may exist in these layers. Finally, we discussed the relationship among earthquakes, velocity structures beneath stations and faults.
基金State Key Basic Research Development and Programming Project (G199804070201) State Natural Science Foundation (40074008).
文摘In this article, we analyze the characters of SV-component receiver function of teleseismic body waves and its advantages in mapping the S-wave velocity structure of crust in detail. Similar to radial receiver function, SV-component receiver function can be obtained by directly deconvolving the P-component from the SV-component of teleseismic recordings. Our analyses indicate that the change of amplitude of SV-component receiver function against the change of epicentral distance is less than that of radial receiver function. Moreover, the waveform of SV-component receiver function is simpler than the radial receiver function and gives prominence to the PS converted phases that are the most sensitive to the shear wave velocity structure in the inversion. The synthetic tests show that the convergence of SV-component receiver function inversion is faster than that of the radial receiver function inversion. As an example, we investigate the S-wave velocity structure beneath HIA sta-tion by using the SV-component receiver function inversion method.
基金Supported by National Natural Science Foundation of China(No.50378032and No.50538030)Associated Foundation of Earthquake Science(No.201009)Foundation of Heilongjiang Institute of Science and Technology(No.04-15).
文摘Microtremors array observation for estimating S-wave velocity structure from phase velocities of Rayleigh and Love wave on two practical sites in Tangshan area by a China-US joint group are researched.The phase velocities of Rayleigh wave are estimated from vertical component records and those of Love wave are estimated from three-component records of microtremors array using modified spatial auto-correlation method.Haskell matrix method is used in calculating Rayleigh and Love wave phase velocities,and the shallow S-wave velocity structure of two practical sites are estimated by means of a hybrid approach of Genetic Algorithm and Simplex.The results are compared with the PS logging data of the two sites,showing it is feasible to estimate the shallow S-wave velocity structure of practical site from the observation of microtremor array.