Aseismic negative dislocation model and deformation analysis of crustal horizontal movement during 1999——2001 in the northeast margin of Qinghai-Xizang block

Through numerical simulation for GPS data, aseism/c negative dislocation model for crustal horizontal movement during 1999-2001 in the northeast margin of Qinghai-Xizang block is presented, combined with the spatial distri-bution of apparent strain field in this area, the characteristics of motion and deformation of active blocks and their boundary faults, together with the place and intensity of strain accumulation are analyzed. It is shown that: a) 9 active blocks appeared totally clockwise motion from eastward by north to eastward by south. Obvious sinistral strike-slip and NE-NEE relative compressive motion between the blocks separated by Qilianshan-Haiyuan fault zone was discovered; b) 20 fault segments (most of them showed compression) locked the relative motion between blocks to varying degrees, among the total, the mid-east segment of Qilianshan fault (containing the place where it meets Riyueshan-Lajishan fault) and the place where it meets Haiyuan fault and Zhuanglanghe fault, more favored accumulation of strain. Moreover, the region where Riyueshan-Lajishan fault meets north boundary of Qaidam block may have strain accumulation to some degree, c) Obtained magnitude of block velocities and locking of their boundaries were less than relevant results for observation in the period of 1993-1999.     

NNR constraint in ITRF2000 and the new global plate motion model NNR-ITRF2000VEL

There are 54 sites employed by ITRF2000 for ITRF2000 orientation. The deficiencies are obvious. First, these sites cannot well represent the rotation rate of the earth crust because there is no selected site in five out of fourteen tectonic plates and three of fourteen plates only have one site each. Second, the total angular momentum of the crust is non-vanishing in ITRF2000, even though it is declared that No Net Rotation (NNR) with respect to NNR-NUVEL1A is imposed on ITRF2000 construction according to the documentations of ITRF2000. So the NNR condition in conventional terrestrial reference system (CTRS) realization cannot be satisfied in ITRF2000. In this paper, the criteria of site selection for estimating the Euler vectors are suggested; the Tisserand system constraint equation in ITRF construction is derived; and as a product, the global plate motions can be obtained from the ITRF2000 construction.     

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.     

地下水及溶质运移数值模拟系统

介绍了地下水及溶质运移数值模拟的研究进展、数值模型的分类和数值模拟的步骤。根据数值模拟中涉及到的地质信息资源的特点，提出了基于地理信息系统(GIS)技术的地下水及溶质运移数值模拟系统的设计思路。该系统能够对模拟区域的基础地质资料、水文地质勘探资料以及数值模型的输入、输出资料进行动态管理，而且随着勘探、科研和生产的不断进行，通过对数据库的及时更新以及对数值模型的不断校正，可为不同的科研和生产方案的预测与评价提供信息资源。     

高填方路堤互锚式薄壁挡土墙土压力试验研究

在山区修建高等级公路工程中，挡土墙结构广泛应用于高填方路堤工程。在提出一种新型互锚式薄壁挡土墙结构的基础上，通过埋设土压力盒的方法对该种挡土墙的结构型式进行了力学模型试验，并进行了模型土压力强度特性分析。试验及分析表明，该挡土墙结构成本低廉，施工方便，结构合理，可以应用于高填方路堤实际工程。     

易门式大型铜矿床构造成矿动力学模型

针对金属矿床隐伏矿定位预测的关键问题,在提出构造成矿动力学的分支学术方向的基础上,以易门式大型铜矿床为例,从成矿地质背景、矿田构造、构造地球化学和成矿构造应力场等方面进行构造成矿动力学研究,阐述在构造应力场控制下成矿流体运移和聚集的规律。构造地球化学研究认为,构造地球化学异常反映矿体的原生晕;构造地球化学异常受构造控制,其分带特征可指示成矿流体的流向,并提供矿床成因和隐伏矿(化)体相对埋深的信息。构造应力场的控矿特征主要表现在:构造应力场导致控矿构造的形成,并驱动成矿流体的运移;构造应力场控制了成矿能量场;应力和能量的高值集中区分布特点可反映某些构造型式。在此基础上建立构造成矿动力学模型,概括隐伏矿成矿预测准则,据此提出若干重点找矿靶区和靶位,其中部分靶区得到工程验证,表明构造成矿动力学方法对隐伏矿定位预测和评价具有重要的指导意义。     

论层序地层学与含油气系统在油气勘探中的联系——以鄂尔多斯中生代盆地为例

虽然层序地层学与含油气系统理论的各自研究对象和研究方法不同，但由于层序地层学研究对象(沉积岩)是含油气系统研究对象(油气)的载体，二者可以由含油气盆地分析有机的结合起来。层序地层分析包含了对含油气系统的地质要素及成藏作用载体的分析，因而沉积层序的组成单元与含油气系统的地质要素有必然的联系：在一个沉积盆地的数个沉积层序中，成熟烃源岩往往是地史中具一定埋深、分布广、厚度大、有机质含量高的凝缩层，这个凝缩层往往是一个构造超层序的最大海(湖)泛面；储集岩往往是成熟烃源岩之上层序的低水位体系域或紧邻成熟烃源岩的高水位体系域；有效盖层为储集岩之上层序的水进体系域及凝缩层；成熟烃源岩之上的沉积层序为上覆岩层；低水位体系域储集岩有可能沿上倾方向尖灭，被层序界面及其上层序的水进体系域岩性圈闭；除构造裂缝外，低水位体系域下切谷可作为油气向上运移通道；一个含油气系统往往跨越不同的沉积层序甚至构造超层序。     

A neural network approach for structural identification and diagnosis of a building from seismic response data

This work presents a novel procedure for identifying the dynamic characteristics of a building and diagnosing whether the building has been damaged by earthquakes, using a back‐propagation neural network approach. The dynamic characteristics are directly evaluated from the weighting matrices of the neural network trained by observed acceleration responses and input base excitations. Whether the building is damaged under a large earthquake is assessed by comparing the modal parameters and responses for this large earthquake with those for a small earthquake that has not caused this building any damage. The feasibility of the approach is demonstrated through processing the dynamic responses of a five‐storey steel frame, subjected to different strengths of the Kobe earthquake, in shaking table tests. Copyright © 2002 John Wiley & Sons, Ltd.     

A procedure for evaluating seismic energy demand of framed structures

Energy serves as an alternative index to response quantities like force or displacement to include the duration‐related seismic damage effect. A procedure to evaluate the absorbed energy in a multistorey frame from energy spectra was developed. For low‐ to medium‐rise frames, it required a static pushover analysis of the structure to determine the modal yield force and ductility factor of an equivalent single‐degree‐of‐freedom system for the first two modes. The energy spectra were then used to determine the energy contribution of each mode. A procedure was also developed to distribute the energy along the frame height based on energy shapes. This study showed that the second‐mode response in some cases needs to be considered to reflect the energy (or damage) concentration in the upper floors. Copyright © 2002 John Wiley & Sons, Ltd.     

System identification of the Vincent Thomas suspension bridge using earthquake records

The Vincent Thomas Bridge in the Los Angeles metropolitan area, is a critical artery for commercial traffic flow in and out of the Los Angeles Harbor, and is at risk in the seismically active Southern California region, particularly because it straddles the Palos Verdes fault zone. A combination of linear and non‐linear system identification techniques is employed to obtain a complete reduced‐order, multi‐input–multi‐output (MIMO) dynamic model of the Vincent Thomas Bridge based on the dynamic response of the structure to the 1987 Whittier and 1994 Northridge earthquakes. Starting with the available acceleration measurements (which consists of 15 accelerometers on the bridge structure and 10 accelerometers at various locations on its base), an efficient least‐squares‐based time‐domain identification procedure is applied to the data set to develop a reduced‐order, equivalent linear, multi‐degree‐of‐freedom model. Although not the main focus of this study, the linear system identification method is also combined with a non‐parametric identification technique, to generate a reduced‐order non‐linear mathematical model suitable for use in subsequent studies to predict, with good fidelity, the total response of the bridge under arbitrary dynamic environments. Results of this study yield measurements of the equivalent linear modal properties (frequencies, mode shapes and non‐proportional damping) as well as quantitative measures of the extent and nature of non‐linear interaction forces arising from strong ground shaking. It is shown that, for the particular subset of observations used in the identification procedure, the apparent non‐linearities in the system restoring forces are quite significant, and they contribute substantially to the improved fidelity of the model. Also shown is the potential of the identification technique under discussion to detect slight changes in the structure's influence coefficients, which may be indicators of damage and degradation in the structure being monitored. Difficulties associated with accurately estimating damping for lightly damped long‐span structures from their earthquake response are discussed. The technical issues raised in this paper indicate the need for added spatial resolution in sensor instrumentation to obtain identified mathematical models of structural systems with the broadest range of validity. Copyright © 2003 John Wiley & Sons, Ltd.