基于局部采样MCMC方法的时移探地雷达反演

马尔科夫链蒙特卡洛方法(MCMC)是一种启发式的全局寻优算法,可以用来解决概率反演的问题.基于MCMC方法的反演不依赖于准确的初始模型,可以引入任意复杂的先验信息,通过对先验概率密度函数的采样来获得大量的后验概率分布样本,在寻找最优解的过程中可以跳出局部最优得到全局最优解.MCMC方法由于计算量巨大,应用难度较高,在地...     

Weathering,erosion and sediment composition in a high‐gradient river,Calabria, Italy

Source rock lithology and immediate modifying processes, such as chemical weathering and mechanical erosion, are primary controls on fluvial sediment supply. Sand composition and Chemical Index of Alteration (CIA) of parent rocks, soil and fluvial sand of the Savuto River watershed, Calabria (Italy), were used to evaluate the modifications of source rocks through different sections of the basin, characterized by different geomorphic processes, in a sub‐humid Mediterranean climate. The headwaters, with gentle topography, produce a coarse‐grained sediment load derived from deeply weathered gneiss, having sand of quartzofeldspathic composition, compositionally very different from in situ degraded bedrock. Maximum estimated CIA values suggest that source rock has been affected significantly by weathering, and it testifies to a climatic threshold on the destruction of the bedrock. The mid‐course has steeper slopes and a deeply incised valley; bedrock consists of mica‐schist and phyllite with a very thin regolith, which provides large cobble to very coarse sand sediments to the main channel. Slope instability, with an areal incidence of over 40 per cent, largely supplies detritus to the main channel. Sand‐sized detritus of soil and fluvial sand is lithic. Estimated CIA value testifies to a significant weathering of the bedrock too, even if in this part of the drainage basin steeper slopes allow erosion to exceed chemical weathering. The lower course has a braided pattern and sediment load is coarse to medium–fine grained. The river cuts across Palaeozoic crystalline rocks and Miocene siliciclastic deposits. Sand‐sized detritus, contributed from these rocks and homogenized by transport processes, has been found in the quartzolithic distal samples. Field and laboratory evidence indicates that landscape development was the result of extensive weathering during the last postglacial temperature maximum in the headwaters, and of mass‐failure and fluvial erosional processes in the mid‐ and low course. Copyright © 2000 John Wiley & Sons, Ltd.     

Effect of viscous,viscoplastic and friction damping on the response of seismic isolated structures

In this paper the efficiency of various dissipative mechanisms to protect structures from pulse‐type and near‐source ground motions is examined. Physically realizable cycloidal pulses are introduced, and their resemblance to recorded near‐source ground motions is illustrated. The study uncovers the coherent component of some near‐source acceleration records, and the shaking potential of these records is examined. It is found that the response of structures with relatively low isolation periods is substantially affected by the high‐frequency fluctuations that override the long duration pulse. Therefore, the concept of seismic isolation is beneficial even for motions that contain a long duration pulse which generates most of the unusually large recorded displacements and velocities. Dissipation forces of the plastic (friction) type are very efficient in reducing displacement demands although occasionally they are responsible for substantial permanent displacements. It is found that the benefits by hysteretic dissipation are nearly indifferent to the level of the yield displacement of the hysteretic mechanism and that they depend primarily on the level of the plastic (friction) force. The study concludes that a combination of relatively low friction and viscous forces is attractive since base displacements are substantially reduced without appreciably increasing base shears and superstructure accelerations. Copyright © 2000 John Wiley & Sons, Ltd.     

Dynamic analysis of sliding structures with unsynchronized support motions

The dynamic analysis of sliding structures is complicated due to the presence of friction. Synchronization of the kinematics of all the isolation bearings is often granted to simplify the task. This, however, may lead to inaccurate prediction of the structural responses under certain circumstances. Stepped structures or continuous bridges with seismic isolation are notable examples where unsynchronized bearing motions are expected. In this paper, a logically simple and numerically efficient procedure is proposed to solve the dynamic problem of sliding systems with unsynchronized support motions. The motion equations for the sliding and non‐sliding modes of the isolated structure are unified into a single equation that is represented as a difference equation in a discrete‐time state‐space form and the base shear forces between the sliding interfaces can be determined through simple matrix algebraic analysis. The responses of the sliding structure can be obtained recursively from the discrete‐time version of the motion equation with constant integration time step even during the transitions between the non‐sliding and sliding phases. Therefore, both accuracy and efficiency in the dynamic analysis of the highly non‐linear system can be enhanced to a large extent. Rigorous assessment of seismic structures with unsynchronized support motions has been carried out for both a stepped structure and a continuous bridge. Effectiveness of friction pendulum bearings for earthquake protection of such structures has been verified. Moreover, evident unsynchronized sliding motions of the friction bearings have been observed, confirming the necessity to deal with each of the bearings independently in the analytical model. Copyright © 2000 John Wiley & Sons, Ltd.     

Ground motion characteristics of the Chi‐Chi earthquake of 21 September 1999

The purpose of this paper is to investigate the ground motion characteristics of the Chi‐Chi earthquake (21 September 1999) as well as the interpretation of structural damage due to this earthquake. Over 300 strong motion records were collected from the strong motion network of Taiwan for this earthquake. A lot of near‐field ground motion data were collected. They provide valuable information on the study of ground motion characteristics of pulse‐like near‐field ground motions as well as fault displacement. This study includes: attenuation of ground motion both in PGA and spectral amplitude, principal direction, elastic and inelastic response analysis of a SDOF system subjected to near‐field ground motion collected from this event. The distribution of spectral acceleration and spectral velocity along the Chelungpu fault is discussed. Based on the mode decomposition method the intrinsic mode function of ground acceleration of this earthquake is examined. A long‐period wave with large amplitude was observed in most of the near‐source ground acceleration. The seismic demand from the recorded near‐field ground motion is also investigated with an evaluation of seismic design criteria of Taiwan Building Code. Copyright © 2000 John Wiley & Sons, Ltd.     

Modified state–space procedures for pseudodynamic testing

The existing on‐line numerical integration algorithms are derived from the Newmark method, which is based on an approximation of derivatives in the differential equation. The state–space procedure (SSP), based on an interpolation of the discrete excitation signals for piecewise convolution integral, has been confirmed as more reliable than the Newmark method in terms of numerical accuracy and stability. In an attempt to enhance the pseudodynamic test, this study presents an on‐line integration algorithm (referred to as the OS–SSP method) via an integration of the state–space procedure with Nakashima's operator‐splitting concept. Numerical stability and accuracy assessment of the proposed algorithm in addition to the explicit Newmark method and the OS method were investigated via an eigenvalue, frequency‐domain and time‐domain analysis. Of the on‐line integration algorithms investigated, the OS–SSP method is demonstrated as the most accurate method with an acceptable stability (although not unconditionally stable) characteristic. Therefore, the OS–SSP method is the most desirable method for pseudodynamic testing if the numerical stability criterion (Δt/T⩽0.5) is ensured for every vibration mode involved. Copyright © 2001 John Wiley & Sons, Ltd.     

龙门山地区历史强震数值模拟

考虑区域地质构造差异、主要活动断裂分布特征和地表附加重力影响,建立反映龙门山地区地表起伏和岩石圈分层的三维粘弹性有限元模型。以GPS为约束重建研究区现今构造应力场,依次模拟龙门山地区1900年以来发生的5次M_S 7.0以上地震,从库仑应力和等效应力角度,分析应力场演化对强震的影响以及强震间的相互作用关系。研究结果表明:从库仑应力角度,有3次地震对后续地震有促进作用,其中汶川地震对芦山地震有触发作用;从等效应力角度,有4次地震对后续地震的发生具有加速作用。     

2016年3月21日金塔县M_S 4.7地震强震动记录分析

2016年3月21日11时29分甘肃省酒泉市金塔县发生M_S 4.7地震,甘肃省强震动台网中心收集到5组三分向加速度记录,常规处理后得到15条零基线校正后的三分向峰值加速度(PGA),选取土层台基类型的62XID台站加速度记录(PGA数值最大)进行分析,绘制校正后的加速度、速度和位移时程曲线及加速度反应谱,为中国西部工程震害调查提供参考资料。     

The influence of tectonic subsidence on volcaniclastic sedimentation: The Cretaceous upper Daeri Member,Wido Island,Korea

The Daeri Member, a Cretaceous volcanic–sedimentary succession, can be divided into lower, middle, and upper parts based on vertical changes in its lithologic characters. The lower Daeri Member is composed of siliciclastic deposits formed in a semi‐arid floodplain environment, which is overlain by the middle Daeri Member consisting mainly of andesite lava flow. After the emplacement of the andesite, activities of intrabasinal normal faults created accommodation on hanging wall blocks together with the development of intrabasinal topographic relief. The upper Daeri Member occurs only in hanging wall blocks and is composed of rhyolitic volcaniclastic sediments formed during an explosive volcanic eruption. Following the eruption, owing to semi‐arid climatic conditions and the destruction of vegetation, the eruptive materials were easily remobilized and deposited by episodic sediment gravity flows, resulting in deposition of the resedimented volcaniclastic deposits with sheet‐like geometry. Away from the intrabasinal normal faults, the resedimented volcaniclastic deposits show a decrease in grain size together with changes in inferred depositional processes from debris flows to hyperconcentrated flows and supercritical sheetfloods. This suggests that the resedimented volcaniclastic deposits were stacked on alluvial fan environments induced by intrabasinal topographic relief associated with normal fault activities. In addition, episodic movement of the faults gave rise to periodic fluctuation of the accommodation and an increase in gradient of the alluvial fan surface, resulting in the development of coarsening‐upward trends in the resedimented volcaniclastic deposits. The development of the alluvial fan and the coarsening‐upward trends indicate that dynamic tectonic subsidence and concomitant changes in the intrabasinal physiographic relief influenced the depositional processes and sizes of the transported volcaniclastic sediments of the upper Daeri Member. Thus, it is necessary to carefully observe tectonic signatures in volcaniclastic successions, particularly the syneruptive lithofacies, in order to reconstruct the tectonic and volcanic histories of receiving basins.     

Numerical analysis of multicomponent responses of surface-hole transient electromagnetic method

We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.