Attenuation of ground vibrations using in-filled wave barriers

Ground vibrations generated by construction activities can adversely affect the structural health of adjacent buildings and foundations supporting them. Therefore propagation and rate of attenuation of construction induced ground vibrations is important during construction activities, particularly in urban areas where constructions are carried out in the vicinity of existing structures. In practice wave barriers are installed in the ground to mitigate the ground vibration propagation and hence to minimise the effect of ground vibrations on surrounding structures. Different types of fill materials such as bentonite, EPS geofoam and concrete are used in constructing wave barriers. In this study, a three-dimensional finite element model is developed to study the efficiency of different fill materials in attenuating ground vibrations. The model is first verified using data from full scale field experiments, where EPS geofoam has been used as a fill material in wave barriers. Then the same model has been used to evaluate the efficiency of open trenches, water filled wave barriers and EPS geofoam filled wave barriers on attenuation of ground vibrations. EPS geofoam is found to be the most efficient fill material, providing attenuation efficiency closer to open trenches. The efficiency of EPS geofoam and water filled wave barriers can be significantly increased by increasing the depth of the wave barrier.     

Analysis of ground vibrations due to underground trains by 2.5D finite/infinite element approach

The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, are modeled by the finite elements, and the wave propagation properties of the far field extending to infinity are modeled by the infinite elements. One particular feature of the 2.5D approach is that it enables the computation of the three-dimensional response of the half-space, taking into account the load-moving effect, using only a two-dimensional profile. Although the 2.5D finite/infinite element approach shows a great advantage in studying the wave propagation caused by moving trains, attention should be given to the calculation aspects, such as the rules for mesh establishment, in order to avoid producing inaccurate or erroneous results. In this paper, some essential points for consideration in analysis are highlighted, along with techniques to enhance the speed of the calculations. All these observations should prove useful in making the 2.5D finite/infinite element approach an effective one.     

黏弹性介质中瑞利波频散曲线和衰减系数曲线的反演

瑞利波具有能量大、信噪比高等特点,可以用来反演介质内部的力学信息,近年来在浅层地球物理勘探、深层地震学研究以及超声波无损检测等多个领域都有较广泛的应用。目前大多数瑞利波的应用中都假设介质是弹性的,然而实际中岩石、土壤和金属等介质都在一定程度上体现出了黏弹性。当介质的黏弹性较强时仍然采用弹性假设研究其中瑞利波的反演将增大误差,因此有必要考虑黏弹性介质中的瑞利波反演,但是目前这方面的研究仍不够深入。本文研究黏弹性介质中瑞利波频散曲线和衰减系数曲线的反演问题,给出其在半空间中联合反演的方法,并对该方法的误差进行分析。     

Attenuation and dispersion ofSH waves due to scattering by randomly distributed cracks

The effect of randomly distributed cracks on the attenuation and dispersion ofSH waves is theoretically studied. If earthquake ruptures are caused by sudden coalescence of preexisting cracks, it will be crucial for earthquake prediction to monitor the temporal variation of the crack distribution. Our aim is to investigate how the property of crack distribution is reflected in the attenuation and dispersion of elastic waves.We introduce the stochastic property, in the mathematical analysis, for the distributions of crack location, crack size and crack orientation. The crack size distribution is assumed to be described by a power law probability density (p(a) a ^– fora _minaa _max according to recent seismological and experimental knowledge, wherea is a half crack length and the range 13 is assumed. The distribution of crack location is assumed to be homogeneous for the sake of mathematical simplicity, and a low crack density is assumed. The stochastic property of each crack is assumed to be independent of that of the other cracks. We assume two models, that is, the aligned crack model and the randomly oriented crack model, for the distribution of crack orientation. All cracks are assumed to be aligned in the former model. The orientation of each crack is assumed to be random in the latter model, and the homogeneous distribution is assumed for the crack orientation. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.We observe the following features common to both the aligned crack model and the randomly oriented crack model. The attenuation coefficientQ ^–1 decays in proportion tok ^–1 in the high frequency range and its growth is proportional tok ² in the low frequency range, wherek is the intrinsic wave number. This asymptotic behavior is parameter-independent, too. The attenuation coefficientQ ^–1 has a broader peak as increases and/ora _min/a _max decreases. The nondimensional peak wave numberk _p a _max at whichQ ^–1 takes the peak value is almost independent ofa _min/a _max for =1 and 2 while it considerably depends ona _min/a _max for =3. The phase velocity is almost independent ofk in the rangeka _max<1 and increases monotonically ask increases in the rangeka _max>1. While the magnitude ofQ ^–1 and the phase velocity considerably depend on the orientation of the crack in the aligned crack model, the above feature does not depend on the crack orientation.The accumulation of seismological measurements suggests thatQ ^–1 ofS waves has a peak at around 0.5 Hz. If this observation is combined with our theoretical results onk _p a _max, the probable range ofa _max of the crack distribution in the earth can be estimated for =1 or 2. If we assume 4 km/sec as theS wave velocity of the matrix medium,a _max is estimated to range from 2 to 5 km. We cannot estimatea _max in a narrow range for =3.     

Laboratory measurements of spatial fluctuation and attenuation of elastic waves by scattering due to random heterogeneities

To study the effects of strong scattering on elastic waves, spatial fluctuation and scattering attenuation ofP waves were examined by laboratory experiments for 2-D models of random media approximately characterized by a triangular correlation function in the range of 2<ka<33, wherek is the wave number anda is the correlation distance of the heterogeneities, i.e., the heterogeneity size. The results obtained are as follows: (1) Forka>10, both the intensity and the correlation distance of the amplitude fluctuation are approximate for any phase of theP-wave train. The correlation distance nearly agrees with the heterogeneity size. These fluctuation properties are quite consistent with the theoretical prediction by the forward-scattering approximation. (2) For 3<ka<6, the fluctuation intensity becomes stronger in later phases of theP-wave train. This shows that scattering is approximately isotropic, and therefore, the scattered energy increases with time within theP-wave train. The correlation distance of the amplitude fluctuation disagrees with the heterogeneity size, and it shows a frequency-dependent property decreasing from 7a to 4a with the increase ofka from 3 to 6. These properties for 3<ka<6 have not yet been predicted theoretically. (3) Forka<3, though the fluctuation is considerably smaller compared with that ofka>10 and 3<ka<6, the fluctuation property is considered similar to that of 3<ka<6. (4) The observed scattering attenuation,Q ^–1, increases withka forka<3, has a peak aroundka=35, and then decreases withka. (5) When _min = 15° and = 0.075, the theoreticalQ ^–1 curve, predicted by the approximate theory of Wu, roughly matches the observedQ ^–1 values, where _min is the minimum scattering angle measured from the propagation direction of theP waves and is the rms of fractional velocity fluctuation. This suggests that the energy scattered in the range of >15° is lost from theP waves, while the energy scattered in the range of <15° is retained; and that the approximate theory overestimates by about three times the value of the model media used owing to the neglect of multiple scattering. (6) When the size of velocity heterogeneities responsible for forward scattering at 3<ka<6 is estimated from the _min value of 15° on the basis of Wu's theory, it nearly agrees with the correlation distance for the initial phase of theP-wave train.     

Experimental investigation of railway train-induced vibrations of surrounding ground and a nearby multi-story building

In this paper, a field experiment was carried out to study train-induced environmental vibrations. During the field experiment, velocity responses were measured at different locations of a six-story masonry structure near the Beijing-Guangzhou Railway and along a small road adjacent to the building. The results show that the velocity response levels of the environmental ground and the building floors increase with train speed, and attenuate with the distance to the railway track. Heavier freight trains indu...     

Three-dimensional analysis of rows of piles as passive barriers for ground vibration isolation

This paper presents a theoretical study on the ground vibration isolation efficiency by multi-rows of piles as passive barrier in a three-dimensional context. Integral equations governing Rayleigh wave scattering are derived according to the Green's solution of Lamb problem. The integral equations are solved accurately and efficiently with an iteration technique. They are used to predict the complicated Rayleigh wave field that is generated by a number of irregular scatterers embedded in an elastic half-space solid. The method is verified with a numerical solution available in the literature for a simplified Rayleigh wave scattering problem. Passive isolation effectiveness of ground vibration with two or three rows of small piles is further studied in detail. Effects of relevant parameters on the effectiveness of vibration isolation are analyzed and presented.     

Local amplification of deep mining induced vibrations part.2: Simulation of ground motion in a coal basin

This work investigates the impact of deep coal mining induced vibrations on surface constructions using numerical tools. An experimental study of the geological site amplification and of its influence on mining induced vibrations has already been published in the previous paper (Part 1: Experimental evidence for site effects in a coal basin). Measurements have shown the existence of an amplification area in the southern part of the basin where drilling data have shown the presence of particularly fractured and soft stratigraphic units. The present study, using the boundary element method (BEM) in the frequency domain, first investigates canonical geological structures in order to get general results for various sites. The amplification level at the surface is given as a function of the shape of the basin and of the velocity contrast with the bedrock. Next, the particular coal basin previously studied experimentally (Driad-Lebeau et al. [1]) is modeled numerically by BEM. The amplification phenomena characterized numerically for the induced vibrations are found to be compatible with the experimental findings such as: amplification level, frequency range and location. Finally, the whole work was necessary to fully assess the propagation and amplification of mine induced vibrations. The numerical results quantifying amplification can also be used to study other coal basins or various types of alluvial sites.     

地下衬砌洞室对Rayleigh波的放大作用

本文利用波函数展开法,给出了半空间中圆形衬砌洞室引起的Rayle igh波散射问题的一个级数解,并通过对求解技术的改进,获得了高频解答。论文通过入射波长、洞室直径、衬砌刚度和埋置深度等参数分析,研究了地下洞室对Rayle igh波的放大作用。数值结果表明,地下洞室情况下地表位移幅值可以达到无洞室情况的2.6倍以上;多数情况下,柔性衬砌洞室情况地表位移幅值最大,无衬砌洞室情况次之,刚性衬砌情况最小。     

圆弧状多层沉积谷地在Rayleigh波入射下动力响应的解析解

利用波函数的Fourier-Bessel级数展开法,推导出了含任意多个圆弧状沉积层的谷地在Rayleigh波入射下二维稳态动力响应的解析解,分析了该级数形式解析解随截断项数的收敛问题,讨论了在获取该解析解的过程中,用于计算Rayleigh波自由场的有限Fourier级数的项数和用于模拟水平地表的大圆弧的半径对解的影响.研究结果表明,该解析解能够在一个非常宽的频带内收敛.最后,利用该解析解,在一个较宽的频带范围内,讨论了谷地中沉积介质的成层性,包括沉积层的层数、软弱夹层的存在及其厚度等因素对地面运动的影响.     

Determination of depths to centroids of three-dimensional sources of potential-field anomalies with examples from environmental and geologic applications

A method is developed for determining the depth to the centroid (the geometric center) of ‘semi-compact' sources. The method, called the anomaly attenuation rate (AAR) method, involves computing radial averages of AARs with increasing distances from a range of assumed source centers. For well-isolated magnetic anomalies from ‘semi-compact' sources, the theoretical AARs range from 2 (close to the sources) to 3 (in the far-field region); the corresponding theoretical range of AARs for gravity anomalies is 1 to 2. When the estimated source centroid is incorrect, the AARs either exceed or fall short of the theoretical values. The levelling-off of the far-field AARs near their theoretical maximum values indicates the upper (deeper) bound of the centroid location. Similarly, near-field AARs lower than the theoretical minimum indicate the lower (shallower) bound of the centroid location. It is not always possible to determine usable upper and lower bounds of the centroids because the method depends on characteristics of sources/anomalies and the noise level of the data. For the environmental magnetic examples considered in this study, the determined deeper bounds were within 4% of the true centroid-to-observation distance. For the case of the gravity anomaly from the Bloomfield Pluton, Missouri, USA, determination of only the shallower bound of the centroid location (7 km) was possible. This estimate agrees closely with the centroid of a previously determined three-dimensional model of the Bloomfield Pluton. For satellite magnetic anomalies, the method is appropriate only for high-amplitude, near-circular anomalies due to the inherent low signal-to-noise ratio of satellite magnetic anomalies. Model studies indicate that the AAR method is able to place depths within ±20–30 km of actual center locations from a 400-km observation altitude. Thus, the method may be able to discriminate between upper crustal, lower crustal, and mantle magnetic sources. The results from the prominent Kentucky anomaly are relatively well-resolved (centroid depth 30 km below the Earth's surface). For the Kiruna Magsat anomaly, the deleterious effects from neighboring anomalies make a determination difficult (possible depth could be between 20 and 30 km). The centroid depths are deeper for the Kursk anomaly (40–50 km). These depths may indicate that magnetic anomalies from the near-surface Kursk iron formations (a known contributor) and deep crustal magnetic sources could combine to form the Kursk Magsat anomaly.     

Prediction of free field vibrations due to pile driving using a dynamic soil–structure interaction formulation

This paper presents a numerical model for the prediction of free field vibrations due to vibratory and impact pile driving. As the focus is on the response in the far field where deformations are relatively small, a linear elastic constitutive behaviour is assumed for the soil. The free field vibrations are calculated by means of a coupled FE–BE model based on a subdomain formulation. First, the case of vibratory pile driving is considered, where the contributions of different types of waves are investigated for several penetration depths. In the near field, the soil response is dominated by a vertically polarized shear wave, whereas in the far field, body waves are importantly attenuated and Rayleigh waves dominate the ground vibration. Second, the case of impact pile driving is considered. A linear wave equation model is used to estimate the impact force during the driving process. Apart from the response of a homogeneous halfspace, it is also investigated how the soil stratification influences the ground vibration for the case of a soft layer on a stiffer halfspace. When the penetration depth is smaller than the layer thickness, the layered medium has no significant influence on ground vibrations. However, when the penetration depth is larger than the layer thickness, the influence of the layered medium becomes more significant. The computed ground vibrations are finally compared with field measurements reported in the literature.     

Sensitivity of PSHA results to ground motion prediction relations and logic-tree weights

Epistemic uncertainty in ground motion prediction relations is recognized as an important factor to be considered in probabilistic seismic hazard analysis (PSHA), together with the aleatory variability that is incorporated directly into the hazard calculations through integration across the log-normal scatter in the ground motion relations. The epistemic uncertainty, which is revealed by the differences in median values of ground motion parameters obtained from relations derived for different regions, is accounted for by the inclusion of two or more ground motion prediction relations in a logic-tree formalism. The sensitivity of the hazard results to the relative weights assigned to the branches of the logic-tree, is explored through hazard analyses for two sites in Europe, in areas of high and moderate seismicity, respectively. The analyses reveal a strong influence of the ground motion models on the results of PSHA, particularly for low annual exceedance frequencies (long return periods) and higher confidence levels. The results also show, however, that as soon as four or more relations are included in the logic-tree, the relative weights, unless strongly biased towards one or two relations, do not significantly affect the hazard. The selection of appropriate prediction relations to include in the analysis, therefore, has a greater impact than the expert judgment applied in assigning relative weights to the branches of the logic-tree.     

饱和地基中空沟近场隔振的现场试验与二维半解析边界元分析

动力机器运行和车辆行驶等会产生振动污染,危及邻近建筑物安全和干扰精密仪器设备正常运行等。这些振动污染可通过在地基中设置空沟的方式来降低或消除。针对饱和地基上明置动力机器基础的环境振动影响及空沟近场隔振问题,进行了饱和地基上空沟近场隔振的现场试验,并对试验结果进行了无量纲化分析;基于饱和土半解析边界元法,分别推导了动力机器基础环境振动影响和空沟近场隔振的边界元方程;在此基础上,详细研究了空沟对动力机器基础振动影响的隔振效果,分析了空沟深度、宽度和距振源距离对其隔振效果的影响。结果表明:空沟能够有效的降低动力机器基础的环境振动影响;空沟宽度对其隔振效果影响相对较小,而空沟深度对其隔振效果影响较大,为获得较好的隔振效果,空沟深度建议取1倍Rayleigh波波长;空沟距振源距离对其隔振效果也有较大影响,距离越远则隔振效果也越好,当被保护建筑距振源较远时,建议空沟在被保护建筑附近设置。此外,在某些特殊情况下,空沟隔振系统会由于共振现象而出现隔振效果劣化的现象,在工程设计中应予以注意。     

多道瞬态瑞雷波技术在公路采空塌陷区探测中应用

地下采空塌陷区岩体结构具有明显的结构不连续、松散、断裂、空洞等异常结构特征,本文在分析公路采空塌陷区地球物理场特征的基础上,详细讨论了应用多道瞬态瑞雷波技术探测采空塌陷区的可行性以及相关技术,深入分析基于瑞雷波波场特征、频谱结构、能量衰减、频散特征、相速度以及多阶模等特征探测和评价地层结构的均匀性、地下空洞及异常结构体等地下障碍物的存在及其空间分布的解释标志.滚动式多道瑞雷波探测技术可以有效地探测地下采空区、软弱结构体等异常结构体的赋存状态及空间分布.