首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The purpose of this paper is the comparison of P‐wave velocity and velocity anisotropy, measured at different scales under laboratory and field conditions. A shallow seismic refraction survey with shot/receiver spacing of up to 10 m was carried out on a flat outcrop of lhertzolite in the southern part of the Balmuccia massif. Oriented rock samples were also obtained from the locality. The particular advantage of the laboratory method used is the possibility of measuring velocity in any direction under controlled conditions. Laboratory tests were made on spherical peridotite samples, 50 mm in diameter, by ultrasonic velocity measurements in 132 directions (meridian and parallel networks) under confining stress ranging from atmospheric to 400 MPa. The mean P‐wave velocity of the field and laboratory data differed by between 20–30%. In addition, P‐wave velocity anisotropy of 25% was detected in the field data. Whereas the anisotropy in the laboratory samples in the same orientation as the field surveys was less than 2%. This observed scaling factor is related to the different sampling sizes and the difference in frequencies of applied elastic waves. With an ultrasonic wavelength of 10 mm, laboratory samples represent a continuum. The field velocities and velocity anisotropy reflect the presence of cracks, which the laboratory rock samples do not contain. Three sub‐vertical fracture sets with differing strikes were observed in the field outcrop. Estimates of fracture stiffness from the velocity anisotropy data are consistent with other published values. These results highlight the difficulty of using laboratory velocity estimates to interpret field data.  相似文献   

2.
A criterion is developed for the simulation of realistic artificial ground motion histories at soft‐soil sites, corresponding to a detailed ground motion record at a reference firm‐ground site. A complex transfer function is defined as the Fourier transform of the ground acceleration time history at the soft‐soil site divided by the Fourier transform of the acceleration record at the firm‐ground site. Working with both the real and the imaginary components of the transfer function, and not only with its modulus, serves to keep the statistical information about the wave phases (and, therefore, about the time variation of amplitudes and frequencies) in the algorithm used to generate the artificial records. Samples of these transfer functions, associated with a given pair of soft‐soil and firm‐ground sites, are empirically determined from the corresponding pairs of simultaneous records. Each function included in a sample is represented as the superposition of the transfer functions of the responses of a number of oscillators. This formulation is intended to account for the contributions of trains of waves following different patterns in the vicinity of both sites. The properties of the oscillators play the role of parameters of the transfer functions. They vary from one seismic event to another. Part of the variation is systematic, and can be explained in terms of the influence of ground motion intensity on the effective values of stiffness and damping of the artificial oscillators. Another part has random nature; it reflects the random characteristics of the wave propagation patterns associated with the different events. The semi‐empirical model proposed recognizes both types of variation. The influence of intensity is estimated by means of a conventional one‐dimensional shear wave propagation model. This model is used to derive an intensity‐dependent modification of the values of the empirically determined model parameters in those cases when the firm‐ground earthquake intensity used to determine these parameters differs from that corresponding to the seismic event for which the simulated records are to be obtained. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

3.
The depth to the surface of a refractor and the seismic velocity within the refractor are very often intimately related. In the shallow environment, increased thicknesses of weathering occur in areas of jointing, shearing or lithological variations, and these zones of deeper weathering can have lower subweathering refractor velocities. This association is important in geotechnical investigations and in the measurement of weathering thicknesses and sub-weathering velocities for statics corrections for reflection seismic surveys. Algorithms, which employ forward and reverse traveltime data and which explicitly accommodate the offset distance through the process known as refraction migration, are necessary if detailed structure on a refractor and rapid lateral variations of the seismic velocity within it are to be resolved. These requirements are satisfied with wavefront construction techniques, Hales’ method and the generalized reciprocal method (GRM). However, these methods employ refraction migration in fundamentally different manners. Most methods compute an offset distance with an often imprecise knowledge of the seismic velocities of the overlying layers. In contrast, the GRM uses a range of offset distances from less than to greater than the optimum value, with the optimum value being selected with a minimum-variance criterion. The approach of the GRM is essential where there are undetected layers and where there are rapid variations in the depth to a refractor and the seismic velocity within it. In the latter situations the offset distance necessary to define the seismic velocities can differ considerably from the value required to define depths. The efficacy of the GRM in resolving structure and seismic velocity is demonstrated with three model studies and two field examples.  相似文献   

4.
The seismic response spectrum defines the amplitude of the load, but it does not specify the number of cycles for which the load must be resisted by the structure. The amplitude by itself is not sufficient to evaluate the seismic resistance of a structure, because the structure's strength, stiffness and energy‐dissipation capacity reduce with an increase in the number of load cycles. This paper presents a cyclic‐demand spectrum, which, in conjunction with the amplitude spectrum, provides a more complete definition of the seismic load, hence a way to consider the degradation in strength, stiffness and energy‐dissipation capacity in a rational manner. Similarly to three amplitude parameters (peak ground acceleration, peak ground velocity, and peak ground displacement), three cyclic‐demand parameters are introduced for stiff, moderately stiff, and flexible systems. A design example is presented to illustrate the use of the cyclic‐demand spectrum. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

5.
The Binchuan region of Yunnan is a structurally complex region with mountains, basins, and active faults. In this situation,seismic wave propagation exhibits complex characteristics due to strong heterogeneity of underground media instead of following the great-circle path. In order to obtain a high-resolution shallow crustal structure, a dense seismic array was deployed during March 21 to May 30, 2017 in this area. To better understand the complexities of seismic wave propagation in this region, we perform array-based frequency-domain beamforming analysis and single-station based polarization analysis to investigate the characteristics of seismic wave propagation, using airgun-generated P-wave signals recorded by dense array stations in this experiment. The results from these two methods both reveal similar but complex characteristics of seismic wave propagation in the Binchuan basin. The azimuth anomalies off the great-circle path are quite large with values up to 30°, which is caused by strong structural heterogeneity in the very shallow crust. Our research provide a better understanding of the complex geologic structures in this area and provide guidance for detecting concealed faults and distribution of velocity anomalies.  相似文献   

6.
Semi‐active variable stiffness resettable devices can reduce seismic demands and damages in structures. Despite their advantages, variable stiffness resettable devices are under‐utilized mainly because of the shortage of fundamental research in quantifying the sensitivity of key seismic response parameters, and losses, in structures that use such systems for seismic hazard mitigation. Within this setting, the research summarized herein measures the effectiveness of semi‐active resettable energy dissipating devices in the Single‐Degree‐of‐Freedom domain aiming at quantifying the sensitivity of their seismic response to variation in control parameters and generating the required knowledge to utilize such semi‐active devices in the Multi‐Degree‐of‐Freedom domain. The performance (i.e. maximum relative displacement and peak absolute acceleration demands) of Single‐Degree‐of‐Freedom systems with an array of semi‐active control logics under various dynamic excitation regimes is studied. Two sets of 40 ground motions representing various seismic loading conditions (i.e. pulse‐like and rock‐site ground motions) are used, and an efficient control logic for mitigating these seismic demands is proposed. Numerical results show that proposed control logic enables a decrease of 40–60% for both maximum relative displacement and seismic base shear and 15–25% decrease for peak absolute acceleration. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

7.
We develop a semi‐empirical model which combines the theoretical model of Xu and White and the empirical formula of Han, Nur and Morgan in sand–clay environments. This new model may be used for petrophysical interpretation of P‐ and S‐wave velocities. In particular, we are able to obtain an independent estimation of aspect ratios based on log data and seismic velocity, and also the relationship between velocities and other reservoir parameters (e.g. porosity and clay content), thus providing a prediction of shear‐wave velocity. To achieve this, we first use Kuster and Toksöz's theory to derive bulk and shear moduli in a sand–clay mixture. Secondly, Xu and White's model is combined with an artificial neural network to invert the depth‐dependent variation of pore aspect ratios. Finally these aspect ratio results are linked to the empirical formula of Han, Nur and Morgan, using a multiple regression algorithm for petrophysical interpretation. Tests on field data from a North Sea reservoir show that this semi‐empirical model provides simple but satisfactory results for the prediction of shear‐wave velocities and the estimation of reservoir parameters.  相似文献   

8.
采用Shake 2000程序,以Turkey Flat试验场地为模型,通过输入不同类型下多种强度的地震波,计算研究多工况下剪切波速测试标准差对地表加速度反应谱和峰值加速度的影响。结论为:(1)浅硬场地上剪切波速测试标准差对地震动的影响很大,影响程度与输入地震波的强度和频率以及场地剪切波速计算值有关;(2)如果将反应谱残差大于20%或加速度峰值差别大于20%定义为统计意义上的不可忽略,那么剪切波速测试标准差对计算结果的影响在大多数情况下均不可忽略;(3)当输入波的卓越周期与场地特征周期接近时,浅硬场地上剪切波速测试标准差引起的反应谱变化非常显著;(4)只有当输入波的卓越周期与场地特征周期相差较大且输入波强度偏小时,剪切波速测试标准差引起的反应谱变化才可略去;(5)当浅硬场地上剪切波速实测结果低于统计均值时,地震动计算结果的偏差一般明显大于剪切波速实测结果,高于统计均值时引起的偏差,且地震输入越强表现越明显。  相似文献   

9.
Time‐lapse refraction can provide complementary seismic solutions for monitoring subtle subsurface changes that are challenging for conventional P‐wave reflection methods. The utilization of refraction time lapse has lagged behind in the past partly due to the lack of robust techniques that allow extracting easy‐to‐interpret reservoir information. However, with the recent emergence of the full‐waveform inversion technique as a more standard tool, we find it to be a promising platform for incorporating head waves and diving waves into the time‐lapse framework. Here we investigate the sensitivity of 2D acoustic, time‐domain, full‐waveform inversion for monitoring a shallow, weak velocity change (?30 m/s, or ?1.6%). The sensitivity tests are designed to address questions related to the feasibility and accuracy of full‐waveform inversion results for monitoring the field case of an underground gas blowout that occurred in the North Sea. The blowout caused the gas to migrate both vertically and horizontally into several shallow sand layers. Some of the shallow gas anomalies were not clearly detected by conventional 4D reflection methods (i.e., time shifts and amplitude difference) due to low 4D signal‐to‐noise ratio and weak velocity change. On the other hand, full‐waveform inversion sensitivity analysis showed that it is possible to detect the weak velocity change with the non‐optimal seismic input. Detectability was qualitative with variable degrees of accuracy depending on different inversion parameters. We inverted, the real 2D seismic data from the North Sea with a greater emphasis on refracted and diving waves’ energy (i.e., most of the reflected energy was removed for the shallow zone of interest after removing traces with offset less than 300 m). The full‐waveform inversion results provided more superior detectability compared with the conventional 4D stacked reflection difference method for a weak shallow gas anomaly (320 m deep).  相似文献   

10.
Seismic wave propagation in viscoelastic and saturated rock   总被引:1,自引:0,他引:1  
Considering engineering practice, the viscoelastic two-phase model is adopted, seismic wave propagation in saturated rock is studied. Not only the effect of the viscosity of rock skeleton but also the effect of ground water on the propagation of the seismic wave can be considered by this model, the propagation characteristics of seismic wave in saturated rock can be understood comprehensively and the model is more reasonable than other model by which seismic wave propagation is studied. The effect of frequency, water content and viscosity constant on the wave velocity and attenuation are studied by numerical examples and some valuable conclusions are drawn.  相似文献   

11.
Experimental measurements of fracture-induced seismic waves velocity variations at frequencies ~ 1 kHz, ~ 40 kHz and ~ 1 MHz were performed directly in the field at the rocky outcrop and in the laboratory on specific rock samples collected from the outcrops. The peridotite–lherzolite outcrop appeared macroscopically uniform and contained three systems of visible parallel sub-vertical fractures. This rock has substantial bulk density and higher than average value of seismic wave velocity. The presence of fracture systems gives rise to its velocity anisotropy. The seismic waves passing through the rock fractures are subject to velocity dispersion and frequency dependent attenuation. Our data, obtained from field and laboratory measurements, were compared with theoretical model predictions. In this model we successfully used displacement discontinuity approach. For the velocity dispersion evaluation we used multi-frequency measurements. The a priori observation of orientations and densities of fracture sets allowed evaluation of their stiffness. Our approach revealed that the first arrivals of seismic waves can be used for evaluation of P-wave group velocities, the specific case, in which we expect anomalous velocity dispersion. Our observations contribute to the issue of up-scaling of well-log derived velocities in fractured rock to the scale of standard seismic exploration frequencies.  相似文献   

12.
This paper focuses on the effects of long‐period pulse of near‐fault ground motions on the structural damage potential. Two sets of near‐fault ground motion records from Chi‐Chi, Taiwan earthquake and Northridge earthquake with and without distinct pulse are selected as the input, and the correlation analysis between 30 non‐structure‐specific intensity measure parameters and maximum inelastic displacements and energy responses (input energy and hysteretic energy) of bilinear single degree of freedom systems are conducted. Based on the frequency characteristic of near‐fault ground motions with remarkable long‐period components, two intensity indices are proposed, namely, the improved effective peak acceleration (IEPA) and improved effective peak velocity (IEPV). In addition a new characteristic period of these ground motions is defined based on IEPA and IEPV. Numerical results illustrate that the intensity measure parameters related to ground acceleration present the best correlation with the seismic responses for rigid systems; the velocity‐related and displacement‐related parameters are better for medium‐frequency systems and flexible systems, respectively. The correlation curves of near‐fault ground motions with velocity pulse differ from those of ground motions without pulse. Moreover, the improved parameters IEPA and IEPV of near‐fault impulsive ground motions enhance the performance of intensity measure of corresponding conventional parameters, i.e. EPA and EPV. The new characteristic period based on IEPA and IEPV can better reflect the frequency content of near‐fault ground motions. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

13.
本文选取山东地区2个Ⅲ类场地的工程地质勘探及土层剪切波速等资料,将土层厚度按5个深度段,每个分段给出了4个土层剪切波速的改变量,通过改变不同深度段土层剪切波速,建立了19种土层地震反应分析模型,分析了不同深度段,不同概率水平下土层剪切波速的变化对场地地震动参数的影响。研究表明,不同深度段土层剪切波速的变化对场地地震动参数的影响有差异。具体表现为,土层剪切波速的改变在1—10m、11—40m和地震输入界面处三个深度段对地震动加速度峰值影响较大;其中,41—70m和71—100m两个深度段剪切波速的改变对地震动加速度峰值影响小;在土层深度1—10m时,剪切波速降低,峰值变大,剪切波速的改变与峰值的改变呈负相关;在其它深度段,剪切波速降低,峰值变小,剪切波速的改变与峰值的改变呈正相关。剪切波速的改变在1—10m和11—40m两个深度段对地震加速度反应谱影响较大;在41—70m、71—100m和地震输入界面三个深度段对地震加速度反应谱影响很小。  相似文献   

14.
Three‐dimensional seismic survey design should provide an acquisition geometry that enables imaging and amplitude‐versus‐offset applications of target reflectors with sufficient data quality under given economical and operational constraints. However, in land or shallow‐water environments, surface waves are often dominant in the seismic data. The effectiveness of surface‐wave separation or attenuation significantly affects the quality of the final result. Therefore, the need for surface‐wave attenuation imposes additional constraints on the acquisition geometry. Recently, we have proposed a method for surface‐wave attenuation that can better deal with aliased seismic data than classic methods such as slowness/velocity‐based filtering. Here, we investigate how surface‐wave attenuation affects the selection of survey parameters and the resulting data quality. To quantify the latter, we introduce a measure that represents the estimated signal‐to‐noise ratio between the desired subsurface signal and the surface waves that are deemed to be noise. In a case study, we applied surface‐wave attenuation and signal‐to‐noise ratio estimation to several data sets with different survey parameters. The spatial sampling intervals of the basic subset are the survey parameters that affect the performance of surface‐wave attenuation methods the most. Finer spatial sampling will reduce aliasing and make surface‐wave attenuation easier, resulting in better data quality until no further improvement is obtained. We observed this behaviour as a main trend that levels off at increasingly denser sampling. With our method, this trend curve lies at a considerably higher signal‐to‐noise ratio than with a classic filtering method. This means that we can obtain a much better data quality for given survey effort or the same data quality as with a conventional method at a lower cost.  相似文献   

15.
选取某核电场地控制性钻孔的厚度、剪切波速、密度等实际勘探数据,通过改变回填土剪切波速,分析了回填土不确定性对场地地震动参数的影响。研究结果表明:在回填土层厚度不变和模型总厚度不变的情况下,地表的水平向峰值加速度随着回填土剪切波速的增大而减小,但水平向峰值加速度增幅逐渐减小;回填土剪切波速到达一定的波速就不再影响地表水平峰值加速度;随着回填土剪切波速的增加,整个反应谱的谱值都普遍减小。  相似文献   

16.
Surface waves in seismic data are often dominant in a land or shallow‐water environment. Separating them from primaries is of great importance either for removing them as noise for reservoir imaging and characterization or for extracting them as signal for near‐surface characterization. However, their complex properties make the surface‐wave separation significantly challenging in seismic processing. To address the challenges, we propose a method of three‐dimensional surface‐wave estimation and separation using an iterative closed‐loop approach. The closed loop contains a relatively simple forward model of surface waves and adaptive subtraction of the forward‐modelled surface waves from the observed surface waves, making it possible to evaluate the residual between them. In this approach, the surface‐wave model is parameterized by the frequency‐dependent slowness and source properties for each surface‐wave mode. The optimal parameters are estimated in such a way that the residual is minimized and, consequently, this approach solves the inverse problem. Through real data examples, we demonstrate that the proposed method successfully estimates the surface waves and separates them out from the seismic data. In addition, it is demonstrated that our method can also be applied to undersampled, irregularly sampled, and blended seismic data.  相似文献   

17.
Ultrasonic (500 kHz) P‐ and S‐wave velocity and attenuation anisotropy were measured in the laboratory on synthetic, octagonal‐shaped, silica‐cemented sandstone samples with aligned penny‐shaped voids as a function of pore fluid viscosity. One control (blank) sample was manufactured without fractures, another sample with a known fracture density (measured from X‐ray CT images). Velocity and attenuation were measured in four directions relative to the bedding fabric (introduced during packing of successive layers of sand grains during sample construction) and the coincident penny‐shaped voids (fractures). Both samples were measured when saturated with air, water (viscosity 1 cP) and glycerin (100 cP) to reveal poro‐visco‐elastic effects on velocity and attenuation, and their anisotropy. The blank sample was used to estimate the background anisotropy of the host rock in the fractured sample; the bedding fabric was found to show transverse isotropy with shear wave splitting (SWS) of 1.45 ± 1.18% (i.e. for S‐wave propagation along the bedding planes). In the fractured rock, maximum velocity and minimum attenuation of P‐waves was seen at 90° to the fracture normal. After correction for the background anisotropy, the fractured sample velocity anisotropy was expressed in terms of Thomsen's weak anisotropy parameters ε, γ & δ. A theory of frequency‐dependent seismic anisotropy in porous, fractured, media was able to predict the observed effect of viscosity and bulk modulus on ε and δ in water‐ and glycerin‐saturated samples, and the higher ε and δ values in air‐saturated samples. Theoretical predictions of fluid independent γ are also in agreement with the laboratory observations. We also observed the predicted polarisation cross‐over in shear‐wave splitting for wave propagation at 45° to the fracture normal as fluid viscosity and bulk modulus increases.  相似文献   

18.
The predominant period and corresponding equivalent viscous damping ratio, also known in various loading codes as effective period and effective damping coefficient, are two important parameters employed in the seismic design of base‐isolated and conventional building structures. Accurate determination of these two parameters can reduce the uncertainty in the computation of lateral displacement demands and interstory drifts for a given seismic design spectrum. This paper estimates these two parameters from data sets recorded from a full‐scale five‐story reinforced concrete building subjected to seismic base excitations of various intensities in base‐isolated and fixed‐base configurations on the outdoor shake table at the University of California, San Diego. The scope of this paper includes all test motions in which the yielding of the reinforcement has not occurred and the response can still be considered ‘elastic’. The data sets are used with three system identification methods to determine the predominant period of response for each of the test configurations. One of the methods also determines the equivalent viscous damping ratio corresponding to the predominant period. It was found that the predominant period of the fixed‐base building lengthened from 0.52 to 1.30 s. This corresponded to a significant reduction in effective system stiffness to about 16% of the original stiffness. The paper then establishes a correlation between predominant period and peak ground velocity. Finally, the predominant periods and equivalent viscous damping ratios recommended by the ASCE 7‐10 loading standard are compared with those determined from the test building. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

19.
A method for generating a suite of synthetic ground motion time‐histories for specified earthquake and site characteristics defining a design scenario is presented. The method employs a parameterized stochastic model that is based on a modulated, filtered white‐noise process. The model parameters characterize the evolving intensity, predominant frequency, and bandwidth of the acceleration time‐history, and can be identified by matching the statistics of the model to the statistics of a target‐recorded accelerogram. Sample ‘observations’ of the parameters are obtained by fitting the model to a subset of the NGA database for far‐field strong ground motion records on firm ground. Using this sample, predictive equations are developed for the model parameters in terms of the faulting mechanism, earthquake magnitude, source‐to‐site distance, and the site shear‐wave velocity. For any specified set of these earthquake and site characteristics, sets of the model parameters are generated, which are in turn used in the stochastic model to generate the ensemble of synthetic ground motions. The resulting synthetic acceleration as well as corresponding velocity and displacement time‐histories capture the main features of real earthquake ground motions, including the intensity, duration, spectral content, and peak values. Furthermore, the statistics of their resulting elastic response spectra closely agree with both the median and the variability of response spectra of recorded ground motions, as reflected in the existing prediction equations based on the NGA database. The proposed method can be used in seismic design and analysis in conjunction with or instead of recorded ground motions. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

20.
In this case study we consider the seismic processing of a challenging land data set from the Arabian Peninsula. It suffers from rough top‐surface topography, a strongly varying weathering layer, and complex near‐surface geology. We aim at establishing a new seismic imaging workflow, well‐suited to these specific problems of land data processing. This workflow is based on the common‐reflection‐surface stack for topography, a generalized high‐density velocity analysis and stacking process. It is applied in a non‐interactive manner and provides an entire set of physically interpretable stacking parameters that include and complement the conventional stacking velocity. The implementation introduced combines two different approaches to topography handling to minimize the computational effort: after initial values of the stacking parameters are determined for a smoothly curved floating datum using conventional elevation statics, the final stack and also the related residual static correction are applied to the original prestack data, considering the true source and receiver elevations without the assumption of nearly vertical rays. Finally, we extrapolate all results to a chosen planar reference level using the stacking parameters. This redatuming procedure removes the influence of the rough measurement surface and provides standardized input for interpretation, tomographic velocity model determination, and post‐stack depth migration. The methodology of the residual static correction employed and the details of its application to this data example are discussed in a separate paper in this issue. In view of the complex near‐surface conditions, the imaging workflow that is conducted, i.e. stack – residual static correction – redatuming – tomographic inversion – prestack and post‐stack depth migration, leads to a significant improvement in resolution, signal‐to‐noise ratio and reflector continuity.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号