This paper evaluates the seismic vulnerability of different classes of typical bridges in California when subjected to seismic shaking or liquefaction-induced lateral spreading. The detailed structural configurations ...This paper evaluates the seismic vulnerability of different classes of typical bridges in California when subjected to seismic shaking or liquefaction-induced lateral spreading. The detailed structural configurations in terms of superstructure type, connection, continuity at support and foundation type, etc. render different damage resistant capability. Six classes of bridges are established based on their anticipated failure mechanisms under earthquake shaking. The numerical models that are capable of simulating the complex soil-structure interaction effects, nonlinear behavior of columns and connections are developed for each bridge class. The dynamic responses are obtained using nonlinear time history analyses for a suite of 250 earthquake motions with increasing intensity. An equivalent static analysis procedure is also implemented to evaluate the vulnerability of the bridges when subjected to liquefaction-induced lateral spreading. Fragility functions for each bridge class are derived and compared for both seismic shaking （based on nonlinear dynamic analyses） and lateral spreading （based on equivalent static analyses） for different performance states. The study finds that the fragility functions due to either ground shaking or lateral spreading show significant correlation with the structural characterizations, but differences emerge for ground shaking and lateral spreading conditions. Structural properties that will mostly affect the bridges＇ damage resistant capacity are also identified.展开更多
Phase spectrum estimation of the seismic wavelet is an important issue in high-resolution seismic data processing and interpretation. On the basis of two patterns of constant-phase rotation and root transform for wave...Phase spectrum estimation of the seismic wavelet is an important issue in high-resolution seismic data processing and interpretation. On the basis of two patterns of constant-phase rotation and root transform for wavelet phase spectrum variation, we introduce six sparse criteria, including Lu’s improved kurtosis criterion, the parsimony criterion, exponential transform criterion, Sech criterion, Cauchy criterion, and the modified Cauchy criterion, to phase spectrum estimation of the seismic wavelet, obtaining an equivalent effect to the kurtosis criterion. Through numerical experiments, we find that when the reflectivity is not a sparse sequence, the estimated phase spectrum of the seismic wavelet based on the criterion function will deviate from the true value. In order to eliminate the influence of non-sparse reflectivity series in a single trace, we apply the method to the multi-trace seismogram, improving the accuracy of seismic wavelet phase spectrum estimation.展开更多
A mature mathematical technique called copula joint function is introduced in this paper, which is commonly used in the financial risk analysis to estimate uncertainty. The joint function is generalized to the n-dimen...A mature mathematical technique called copula joint function is introduced in this paper, which is commonly used in the financial risk analysis to estimate uncertainty. The joint function is generalized to the n-dimensional Frank’s copula. In addition, we adopt two attenuation models proposed by YU and Boore et al, respectively, and construct a two-dimensional copula joint probabilistic function as an example to illustrate the uncertainty treatment at low probability. The results show that copula joint function gives us a better prediction of peak ground motion than that resultant from the simple linear weight technique which is commonly used in the traditional logic-tree treatment of model uncertainties. In light of widespread application in the risk analysis from financial investment to insurance assessment, we believe that the copula-based technique will have a potential application in the seismic hazard analysis.展开更多
This study investigates seismic interferometry in which the Green＇s function is estimated between two receiv- ers by cross-correlation and integration over sources. For smoothly varying source strengths, the dominant...This study investigates seismic interferometry in which the Green＇s function is estimated between two receiv- ers by cross-correlation and integration over sources. For smoothly varying source strengths, the dominant contributions of the correlation integral come from the stationary phase directions in the forward and backward directions from the alignment of the two receivers. Gaussian beams can be used to evaluate the correlation integral and concentrate the amplitudes in a vicinity of the stationary phase regions instead of completely relying on phase interference. Several numerical examples are shown to illustrate how this process works. The use of Gaussian beams for the evaluation of the correlation integral results in stable estimates, and also provides physical insight into the estimation of the Green＇s function based on seismic interferometry.展开更多
Recovery of under-sampled seismic data is a critical problem,in oil and gas exploration,therefore recovery algorithms with iterative shrinkage based on compressed sensing have been recently proposed. However most of t...Recovery of under-sampled seismic data is a critical problem,in oil and gas exploration,therefore recovery algorithms with iterative shrinkage based on compressed sensing have been recently proposed. However most of these algorithms usually adopt a soft shrinkage function,which assumes that all of the sparse coefficients are independent of each other in curvelet or other domains,little attention has so far been devoted to the inter-dependencies of coefficients. In this paper,the dependencies of parent-child curvelet coefficients of seismic data are exploited by Bayesian estimation,moreover the new seismic data recovery algorithm via curvelet-based bivariate shrinkage function is proposed. First the respective parent-child curvelet coefficients joint distribution models of fully-sampled seismic data and noise signal caused by missing traces are established,then the bivariate shrinkage function according to the Bayesian maximum posterior probability estimation is obtained,finally the Landweber iterative shrinkage algorithm is used in the recovery process.When compared with existing recovery algorithms,it is proved that the proposed algorithm can obtain higher PSNR performance,and maintains the texture details better in events of seismic data展开更多
The S wave velocity structure in Changbaishan volcanic region was obtained from teleseismic receiver function modeling. The results show that there exist distinct low velocity layers in crust in volcano area. Beneath ...The S wave velocity structure in Changbaishan volcanic region was obtained from teleseismic receiver function modeling. The results show that there exist distinct low velocity layers in crust in volcano area. Beneath WQD station near to the Tianchi caldera the low velocity layer at 8 km depth is 20 km thick with the lowest S-wave velocity about 2.2 km/s At EDO station located 50 km north of Tianchi caldera, no obvious crustal low velocity layer is detected. In the volcanic region, the thickness of crustal low velocity layer is greater and the lowest velocity is more obvious with the distance shorter to the caldera. It indicates the existence of the high temperature material or magma reservoir in crust near the Tianchi caldera. The receiver functions and inversion result from different back azimuths at CBS permanent seismic station show that the thickness of near surface low velocity layer and Moho depth change with directions. The near surface low velocity layer is obviously thicker in south direction. The Moho depth shows slight uplifting in the direction of the caldera located. We con- sider that the special near surface velocity structure is the main cause of relatively lower prominent frequency of volcanic earthquake waveforms recorded by CBS station. The slight uplifting of Moho beneath Tianchi caldera indicates there is a material exchanging channel between upper mantle and magma reservoir in crust.展开更多
This paper presents a copula technique to develop time-variant seismic fragility curves for corroded bridges at the system level and considers the realistic time-varying dependence among component seismic demands. Bas...This paper presents a copula technique to develop time-variant seismic fragility curves for corroded bridges at the system level and considers the realistic time-varying dependence among component seismic demands. Based on material deterioration mechanisms and incremental dynamic analysis, the time-evolving seismic demands of components were obtained in the form of marginal probability distributions. The time-varying dependences among bridge components were then captured with the best fitting copula function, which was selected from the commonly used copula classes by the empirical distribution based analysis method. The system time-variant fragility curves at different damage states were developed and the effects of time-varying dependences among components on the bridge system fragility were investigated. The results indicate the time-varying dependence among components significantly affects the time-variant fragility of the bridge system. The copula technique captures the nonlinear dependence among component seismic demands accurately and easily by separating the marginal distributions and the dependence among them.展开更多
The concept of Initial Casualty Matrix is introduced. Using some probability distribution functions, the initial casualty matrix of masonry is determined. The dynamic method of seismic casualty assessment is establish...The concept of Initial Casualty Matrix is introduced. Using some probability distribution functions, the initial casualty matrix of masonry is determined. The dynamic method of seismic casualty assessment is established and then applied to the Tangshan earthquake data, with some conclusions drawn.展开更多
Because ambient seismic noise provides estimated Green’s function （EGF） between two sites with high accuracy, Rayleigh wave propagation along the path connecting the two sites is well resolved. Therefore, earthquak...Because ambient seismic noise provides estimated Green’s function （EGF） between two sites with high accuracy, Rayleigh wave propagation along the path connecting the two sites is well resolved. Therefore, earthquakes which are close to one seismic station can be well located with calibration extracting from EGF. We test two algorithms in locating the 1998 Zhangbei earthquake, one algorithm is waveform-based, and the other is traveltime-based. We first compute EGF between station ZHB （a station about 40 km away from the epicenter） and five IC/IRIS stations. With the waveform-based approach, we calculate 1D synthetic single-force Green’s functions between ZHB and other four stations, and obtain traveltime corrections by correlating synthetic Green’s functions with EGFs in period band of 10–30 s. Then we locate the earthquake by minimizing the differential travel times between observed earthquake waveform and the 1D synthetic earthquake waveforms computed with focal mechanism provided by Global CMT after traveltime correction from EGFs. This waveform-based approach yields a location which error is about 13 km away from the location observed with InSAR. With the traveltime-based approach, we begin with measuring group velocity from EGFs as well as group arrival time on observed earthquake waveforms, and then locate the earthquake by minimizing the difference between observed group arrival time and arrival time measured on EGFs. This traveltime-based approach yields accuracy of 3 km, Therefore it is feasible to achieve GT5 （ground truth location with accuracy 5 km） with ambient seismic noises. The less accuracy of the waveform-based approach was mainly caused by uncertainty of focal mechanism.展开更多
The noise data in vertical component records of 85 seismic stations in Fujian Province during 2012 is used as the research object in this paper. The noise data is divided into fiveminute segments to calculate the powe...The noise data in vertical component records of 85 seismic stations in Fujian Province during 2012 is used as the research object in this paper. The noise data is divided into fiveminute segments to calculate the power spectra. The high reference line and low reference line of station are then identified by drawing a probability density function graph( PDF)using the power spectral probability density function. Moreover, according to the anomalies of PDF graphs in 85 seismic stations,the abnormal noise is divided into four categories: dropped packet, low noise, high noise, and median noise anomalies.Afterwards,four selection methods are found by the high or low noise reference line of the stations,and the system of real-time monitoring of seismic noise is formed by combining the four selection methods. Noise records of 85 seismic stations in Fujian Province in July2013 are selected for verification,and the results show that the anomalous noise-recognition system could reach a 90% success rate at most stations and the effect of selection are very good. Therefore,it could be applied to the seismic noise real-time monitoring in stations.展开更多
The extraction of the Earth＇s Green＇s function from field fluctuations is a rapidly growing area of research. The principle of Green＇s function extraction is often related to the requirement of equipartitioning, wh...The extraction of the Earth＇s Green＇s function from field fluctuations is a rapidly growing area of research. The principle of Green＇s function extraction is often related to the requirement of equipartitioning, which stipulates that the energy of field fluctuations is distributed evenly in some sense. We show the meaning of equipartitioning for a variety of different formulations for Green＇s function retrieval. We show that equipartitioning is not a sufficient condition, and provide several examples that illustrate this point. We discuss the implications of lack of equipartitioning for various schemes for the reconstruction of the Green＇s function in seismology. The theory for Green＇s function extraction is usually based on a statistical theory that relies on ensemble averages. Since there is only one Earth, one usually replaces the ensemble average with a time average. We show that such a replacement only makes sense when attenuation is taken into account, and show how the theory for Green＇s function extraction for oscillating systems can be extended to incorporate attenuation.展开更多
基金Supported by:Pacific Earthquake Engineering Research Center Lifelines Program Under Project Task No.9C
文摘This paper evaluates the seismic vulnerability of different classes of typical bridges in California when subjected to seismic shaking or liquefaction-induced lateral spreading. The detailed structural configurations in terms of superstructure type, connection, continuity at support and foundation type, etc. render different damage resistant capability. Six classes of bridges are established based on their anticipated failure mechanisms under earthquake shaking. The numerical models that are capable of simulating the complex soil-structure interaction effects, nonlinear behavior of columns and connections are developed for each bridge class. The dynamic responses are obtained using nonlinear time history analyses for a suite of 250 earthquake motions with increasing intensity. An equivalent static analysis procedure is also implemented to evaluate the vulnerability of the bridges when subjected to liquefaction-induced lateral spreading. Fragility functions for each bridge class are derived and compared for both seismic shaking （based on nonlinear dynamic analyses） and lateral spreading （based on equivalent static analyses） for different performance states. The study finds that the fragility functions due to either ground shaking or lateral spreading show significant correlation with the structural characterizations, but differences emerge for ground shaking and lateral spreading conditions. Structural properties that will mostly affect the bridges＇ damage resistant capacity are also identified.
基金supported by the Major Basic Research Development Program of China (973 Project No. 2007CB209608)
文摘Phase spectrum estimation of the seismic wavelet is an important issue in high-resolution seismic data processing and interpretation. On the basis of two patterns of constant-phase rotation and root transform for wavelet phase spectrum variation, we introduce six sparse criteria, including Lu’s improved kurtosis criterion, the parsimony criterion, exponential transform criterion, Sech criterion, Cauchy criterion, and the modified Cauchy criterion, to phase spectrum estimation of the seismic wavelet, obtaining an equivalent effect to the kurtosis criterion. Through numerical experiments, we find that when the reflectivity is not a sparse sequence, the estimated phase spectrum of the seismic wavelet based on the criterion function will deviate from the true value. In order to eliminate the influence of non-sparse reflectivity series in a single trace, we apply the method to the multi-trace seismogram, improving the accuracy of seismic wavelet phase spectrum estimation.
基金Project of Institute of Crustal Dynamics, China Earthquake Administration (ZDJ2007-1)One Hundred Individual Program of Chinese Academy of Sciences (99M2009M02) National Natural Science Foundation of China (40574022)
文摘A mature mathematical technique called copula joint function is introduced in this paper, which is commonly used in the financial risk analysis to estimate uncertainty. The joint function is generalized to the n-dimensional Frank’s copula. In addition, we adopt two attenuation models proposed by YU and Boore et al, respectively, and construct a two-dimensional copula joint probabilistic function as an example to illustrate the uncertainty treatment at low probability. The results show that copula joint function gives us a better prediction of peak ground motion than that resultant from the simple linear weight technique which is commonly used in the traditional logic-tree treatment of model uncertainties. In light of widespread application in the risk analysis from financial investment to insurance assessment, we believe that the copula-based technique will have a potential application in the seismic hazard analysis.
基金supported by U.S. National Science Foundation EAR06-35611U.S. Air Force contract FA8718-08-C-002the members of the Geo-Mathematical Imaging Group (GMIG) at Purdue University
文摘This study investigates seismic interferometry in which the Green＇s function is estimated between two receiv- ers by cross-correlation and integration over sources. For smoothly varying source strengths, the dominant contributions of the correlation integral come from the stationary phase directions in the forward and backward directions from the alignment of the two receivers. Gaussian beams can be used to evaluate the correlation integral and concentrate the amplitudes in a vicinity of the stationary phase regions instead of completely relying on phase interference. Several numerical examples are shown to illustrate how this process works. The use of Gaussian beams for the evaluation of the correlation integral results in stable estimates, and also provides physical insight into the estimation of the Green＇s function based on seismic interferometry.
基金Sponsored by the National Natural Science Foundation of China(Grant o.61374127)
文摘Recovery of under-sampled seismic data is a critical problem,in oil and gas exploration,therefore recovery algorithms with iterative shrinkage based on compressed sensing have been recently proposed. However most of these algorithms usually adopt a soft shrinkage function,which assumes that all of the sparse coefficients are independent of each other in curvelet or other domains,little attention has so far been devoted to the inter-dependencies of coefficients. In this paper,the dependencies of parent-child curvelet coefficients of seismic data are exploited by Bayesian estimation,moreover the new seismic data recovery algorithm via curvelet-based bivariate shrinkage function is proposed. First the respective parent-child curvelet coefficients joint distribution models of fully-sampled seismic data and noise signal caused by missing traces are established,then the bivariate shrinkage function according to the Bayesian maximum posterior probability estimation is obtained,finally the Landweber iterative shrinkage algorithm is used in the recovery process.When compared with existing recovery algorithms,it is proved that the proposed algorithm can obtain higher PSNR performance,and maintains the texture details better in events of seismic data
基金supported by National Key Project of Scientific and Technical Supporting Programs Funded by Ministry of Science & Technology of China (2006BAC01B04)Joint Seismological Science Foundation of China (106023)Contribution No. is 09FE3006 of Institute of Geophysics,China Earthquake Administration
文摘The S wave velocity structure in Changbaishan volcanic region was obtained from teleseismic receiver function modeling. The results show that there exist distinct low velocity layers in crust in volcano area. Beneath WQD station near to the Tianchi caldera the low velocity layer at 8 km depth is 20 km thick with the lowest S-wave velocity about 2.2 km/s At EDO station located 50 km north of Tianchi caldera, no obvious crustal low velocity layer is detected. In the volcanic region, the thickness of crustal low velocity layer is greater and the lowest velocity is more obvious with the distance shorter to the caldera. It indicates the existence of the high temperature material or magma reservoir in crust near the Tianchi caldera. The receiver functions and inversion result from different back azimuths at CBS permanent seismic station show that the thickness of near surface low velocity layer and Moho depth change with directions. The near surface low velocity layer is obviously thicker in south direction. The Moho depth shows slight uplifting in the direction of the caldera located. We con- sider that the special near surface velocity structure is the main cause of relatively lower prominent frequency of volcanic earthquake waveforms recorded by CBS station. The slight uplifting of Moho beneath Tianchi caldera indicates there is a material exchanging channel between upper mantle and magma reservoir in crust.
基金Natural Science Foundation of China under Grant No.51808376
文摘This paper presents a copula technique to develop time-variant seismic fragility curves for corroded bridges at the system level and considers the realistic time-varying dependence among component seismic demands. Based on material deterioration mechanisms and incremental dynamic analysis, the time-evolving seismic demands of components were obtained in the form of marginal probability distributions. The time-varying dependences among bridge components were then captured with the best fitting copula function, which was selected from the commonly used copula classes by the empirical distribution based analysis method. The system time-variant fragility curves at different damage states were developed and the effects of time-varying dependences among components on the bridge system fragility were investigated. The results indicate the time-varying dependence among components significantly affects the time-variant fragility of the bridge system. The copula technique captures the nonlinear dependence among component seismic demands accurately and easily by separating the marginal distributions and the dependence among them.
文摘The concept of Initial Casualty Matrix is introduced. Using some probability distribution functions, the initial casualty matrix of masonry is determined. The dynamic method of seismic casualty assessment is established and then applied to the Tangshan earthquake data, with some conclusions drawn.
基金supported by Chinese Acadmy of Sciences Fund(No.KCZX-YW-116-1)Joint Seismological Science Fundation of China (Nos.20080878 and 200708035)
文摘Because ambient seismic noise provides estimated Green’s function （EGF） between two sites with high accuracy, Rayleigh wave propagation along the path connecting the two sites is well resolved. Therefore, earthquakes which are close to one seismic station can be well located with calibration extracting from EGF. We test two algorithms in locating the 1998 Zhangbei earthquake, one algorithm is waveform-based, and the other is traveltime-based. We first compute EGF between station ZHB （a station about 40 km away from the epicenter） and five IC/IRIS stations. With the waveform-based approach, we calculate 1D synthetic single-force Green’s functions between ZHB and other four stations, and obtain traveltime corrections by correlating synthetic Green’s functions with EGFs in period band of 10–30 s. Then we locate the earthquake by minimizing the differential travel times between observed earthquake waveform and the 1D synthetic earthquake waveforms computed with focal mechanism provided by Global CMT after traveltime correction from EGFs. This waveform-based approach yields a location which error is about 13 km away from the location observed with InSAR. With the traveltime-based approach, we begin with measuring group velocity from EGFs as well as group arrival time on observed earthquake waveforms, and then locate the earthquake by minimizing the difference between observed group arrival time and arrival time measured on EGFs. This traveltime-based approach yields accuracy of 3 km, Therefore it is feasible to achieve GT5 （ground truth location with accuracy 5 km） with ambient seismic noises. The less accuracy of the waveform-based approach was mainly caused by uncertainty of focal mechanism.
基金sponsored by the National Key Technology R&D Program of China(2009BAK55B00)the Earthquake Industry Research Project(201508012)
文摘The noise data in vertical component records of 85 seismic stations in Fujian Province during 2012 is used as the research object in this paper. The noise data is divided into fiveminute segments to calculate the power spectra. The high reference line and low reference line of station are then identified by drawing a probability density function graph( PDF)using the power spectral probability density function. Moreover, according to the anomalies of PDF graphs in 85 seismic stations,the abnormal noise is divided into four categories: dropped packet, low noise, high noise, and median noise anomalies.Afterwards,four selection methods are found by the high or low noise reference line of the stations,and the system of real-time monitoring of seismic noise is formed by combining the four selection methods. Noise records of 85 seismic stations in Fujian Province in July2013 are selected for verification,and the results show that the anomalous noise-recognition system could reach a 90% success rate at most stations and the effect of selection are very good. Therefore,it could be applied to the seismic noise real-time monitoring in stations.
文摘The extraction of the Earth＇s Green＇s function from field fluctuations is a rapidly growing area of research. The principle of Green＇s function extraction is often related to the requirement of equipartitioning, which stipulates that the energy of field fluctuations is distributed evenly in some sense. We show the meaning of equipartitioning for a variety of different formulations for Green＇s function retrieval. We show that equipartitioning is not a sufficient condition, and provide several examples that illustrate this point. We discuss the implications of lack of equipartitioning for various schemes for the reconstruction of the Green＇s function in seismology. The theory for Green＇s function extraction is usually based on a statistical theory that relies on ensemble averages. Since there is only one Earth, one usually replaces the ensemble average with a time average. We show that such a replacement only makes sense when attenuation is taken into account, and show how the theory for Green＇s function extraction for oscillating systems can be extended to incorporate attenuation.
