A hybrid indirect boundary element—discrete wave number method applied to simulate the seismic response of stratified alluvial valleys

A hybrid indirect boundary element – discrete wavenumber method is presented and applied to model the ground motion on stratified alluvial valleys under incident plane SH waves from an elastic half-space. The method is based on the single-layer integral representation for diffracted waves. Refracted waves in the horizontally stratified region can be expressed as a linear superposition of solutions for a set of discrete wavenumbers. These solutions are obtained in terms of the Thomson–Haskell propagators formalism. Boundary conditions of continuity of displacements and tractions along the common boundary between the half-space and the stratified region lead to a system of equations for the sources strengths and the coefficients of the plane wave expansion. Although the regions share the boundary, the discretization schemes are different for both sides: for the exterior region, it is based on the numerical and analytical integration of exact Green's functions for displacements and tractions whereas for the layered part, a collocation approach is used. In order to validate this approach results are compared for well-known cases studied in the literature. A homogeneous trapezoidal valley and a parabolic stratified valley were studied and excellent agreement with previous computations was found. An example is given for a stratified inclusion model of an alluvial deposit with an irregular interface with the half-space. Results are displayed in both frequency and time domains. These results show the significant influence of lateral heterogeneity and the emergence of locally generated surface waves in the seismic response of alluvial valleys.     

Simplified vibratory characterization of alluvial basins

For the analysis of seismic wave amplification, modal methods are interesting tools to study the modal properties of geological structures. Modal approaches mainly lead to information on such parameters as fundamental frequencies and eigenmodes of alluvial basins. For a specific alluvial deposit in Nice (France), a simplified modal approach involving the Rayleigh method is considered. This approach assumes a set of admissible shape functions for the eigenmodes and allows a fast estimation of the fundamental frequency of the basin. The agreement with other numerical results (Boundary Element Method) and experimental ones is fully satisfactory. The simplified modal method then appears as an efficient mean for the global vibratory characterization of geological structures towards resonance. To cite this article: J.-F. Semblat et al., C. R. Geoscience 335 (2003).     

新疆阿拉尔花岗岩MFC过程的微量元素综合效应

新疆阿拉尔花岗岩的微量元素比相关图显示了线性样点列与幂函数曲线形式的样点列相 复杂图型。该图型用简单混合或者单一结果分异模型都不能解释。本文建立了混合－结晶分异复合过程的微量元素综合效应模型。模型１：从初始混合线引起害虫函数曲线束形式的结晶分异线。特殊情况有结晶分异线简化为通过原点的直线束形式或者与初始混合线重合。     

Application of an optimum design technique for determining the coefficient of consolidation by using piezocone test data

For normally consolidated clay, several researchers have developed a number of theoretical time factors to determine the coefficient of consolidation from piezocone test results. However, depending on assumptions and analytical techniques, it could vary considerably, even for a specific degree of consolidation. In this paper a method is proposed to determine a consistent coefficient of consolidation by applying the concept of an optimum design technique over all ranges of the degree of consolidation. Initial excess pore pressure distribution is assumed to be capable of being obtained by the successive spherical cavity expansion theory. The dissipation of pore pressure is simulated by means of a two-dimensional linear-uncoupled axi-symmetric consolidation analysis. The minimization of differences between measured and predicted excess pore pressure was carried out by the BFGS unconstrained optimum design algorithm with a one-dimensional golden section search technique. By analyzing numerical examples and in-situ test results, it was found that the adopted optimum design technique gives consistent and convergent results.     

Site effects at a vertical accelerometer array near Paducah, Kentucky

The downhole vertical accelerometer array VSAP near Paducah, KY, consists of three-component accelerometers at the surface, the top of the McNairy Formation (−41 m), and the top of the Paleozoic bedrock (—99 m). The array is at the northern end of the Mississippi Embayment, and it was installed to verify the ground-motion modeling for the site, assuming a significant earthquake in the New Madrid Seismic Zone. Accelerograms from 4.2 and 2.0 m_b earthquakes were used to check aspects of the modeling pertaining to linear behavior of the soil column, and to review the soil column models derived by drilling and geotechnical methods and through the use of high-resolution P- and SH-wave seismic refraction in reflection techniques. Results of the study indicate that for the linear case the soil column models derived by the two techniques are equivalent, and that the most important boundary in the soil column, with respect to amplification of the ground motions, is the interface between the limestone bedrock and soil.     

Estimation of Site Effects in Beijing City

— For the realistic modeling of the seismic ground motion in lateral heterogeneous anelastic media, the database of 3-D geophysical structures for Beijing City has been built up to model the seismic ground motion in the City, caused by the 1976 Tangshan and the 1998 Zhangbei earthquakes. The hybrid method, which combines the modal summation and the finite-difference algorithms, is used in the simulation. The modeling of the seismic ground motion, for both the Tangshan and the Zhangbei earthquakes, shows that the thick Quaternary sedimentary cover amplifies the peak values and increases the duration of the seismic ground motion in the northwestern part of the City. Therefore the thickness of the Quaternary sediments in Beijing City is the key factor controling the local ground effects. Four zones are defined on the base of the different thickness of the Quaternary sediments. The response spectra for each zone are computed, indicating that peak spectral values as high as 0.1 g are compatible with past seismicity and can be well exceeded if an event similar to the 1697 Sanhe-Pinggu occurs.     

Use of microtremor measurement for assessing site effects in Northern Belgium – interpretation of the observed intensity during the MS = 5.0 June 11 1938 earthquake

To provide quantitative information on site effects in the northern part of Belgium, forty-seven H/V microtremor measurements were performed with 5 second seismometers over an area of about 15.000 km². Most of the results show a northward regular increase of the fundamental period in agreement with the augmentation of the Mesozoic and Cainozoic soft sedimentsthickness from a few meters 40 km south of Brussels to 900 m at the Netherlands-Belgium border. The measured resonance frequency values wereconsistent with theoretical computations performed at different sites onthe basis of existing information and shallow seismic experiments. At oneparticular site (Uccle) where borehole data were available, microtremor measurements using an array of four seismological stations with differentapertures allowed to obtain the low frequency part of the Rayleigh wave dispersion curve, extending the range covered by the analysis of surfacewaves artificially generated. The Vs profile derived from the surface waveinversion corroborates the 1 Hz natural frequency of the site. Comparison of these results with the macroseismic information concerning the M_S = 5.0 1938 earthquake which occurred 50 km west of Brussels,confirmed the hypothesis that the geological structure of the Brabant massifis likely to control damage distribution during such an earthquake. Comparisonbetween the intensity map of the 1938 earthquake and the resonance period ofsediments obtained by our microtremor study shows a clear relation betweenthe two parameters. During the 1938 earthquake, site effects played a prominent role due to the dimension of the source whose corner frequency wasabout 1 Hz.     

Approximate modal decomposition of inelastic dynamic responses of wall buildings

Two approximate methods for decomposing complicated inelastic dynamic responses of wall buildings into simple modal responses are presented. Both methods are based on the equivalent linear concept, where a non‐linear structure is represented by a set of equivalent linear models. One linear model is used for representing only one vibration mode of the non‐linear structure, and its equivalent linear parameters are identified from the inelastic response time histories by using a numerical optimizer. Several theoretical relations essential for the modal decomposition are derived under the framework of complex modal analysis. Various numerical examinations have been carried out to check the validity of the proposed modal decomposition methods, and the results are quite satisfactory in all cases. Fluctuating bending moment and shear at any location along the wall height contributed by each individual vibration mode can be obtained. Modal contributions to shear and flexural strength demands, as well as the corresponding modal properties, under various seismic loading conditions can also be identified and examined in detail. Furthermore, the effects of higher vibration modes on seismic demands of wall buildings are investigated by using the modal decomposition methods. Several new insights into the complicated inelastic dynamics of multi‐story wall buildings are presented. Copyright © 2004 John Wiley & Sons, Ltd.     

岩石Kaiser效应测地应力原理中的若干问题研究

该文通过对几种不同类型岩石的实验研究，提出了关于在围压条件下影响岩石Kaiser效应的若干问题，分析了它们的声发射特征，为深部地应力测量提供了一种新的解决方法。     

澜沧江中游某崩塌堆积体变形空间效应研究

某崩塌堆积体范围内布置有公路、缆机平台、电站进水口等重要建筑物。其变形稳定性直接影响到治理工程设计、施工安全以及电站的运行。因此，文章在对其环境地质条件、物质结构、形成机制分析的基础上，对其变形空间效应作系统的地质分析；并建立了相应的三维地质力学模型，应用FLAC-3D快速拉各朗日差分程序对其空间变形特征进行了详细的分析。结果表明：收口转向及变形收敛特征使得堆积体下部形成一个“支撑拱”，从而阻碍了上部堆积体向下部的变形与位移，利于堆积体的稳定；同时，三维数值模拟结果也表明：堆积体的变形具有明显的空间效应。其位移矢量表现出一定的收敛性，且变形主要分布在1360～1650m高程；变形方式表现出一定的分区特征，表现出后缘压缩固结，中部滑移，前缘受阻压缩的位移特征。纵向具有明显的收口效应；横向上比较来看，也存在变形受阻的收口特征；收口范围主要集中在1330～1440m高程范围。