Architecture and evolution of the Kushiro submarine canyon in the Kurile Trench forearc slope,North‐western Pacific

An air‐gun survey, conducted over a total distance of 4356 km in the western end of the Kurile Arc offshore, has revealed the architecture and evolution of the Kushiro submarine canyon and Tokachi submarine channels of the Tokachi‐oki forearc basin. The Kushiro submarine canyon, which runs along the eastern margin of the forearc basin, is characterized by an entrenchment of up to several hundred metres in depth. The Tokachi submarine channels, by contrast, occupy the centre of the basin and consist of small, branching and levéed channels. The Kushiro submarine canyon is not connected to the Tokachi River, which has the largest drainage area in eastern Hokkaido, with a catchment area of approximately 9010 km² that includes high mountains and a volcanic region. Instead, the Kushiro submarine canyon exhibits an offset connection/quasi‐connection (probably having been connected during a prior sea‐level lowstand) with the Kushiro River (drainage area of 2500 km²) which contains the Kushiro Swamp at its mouth. To understand this unusual arrangement of rivers and submarine channels, acoustic facies analysis was undertaken to establish the seismic stratigraphy of the area. Subsurface strata can be divided into six seismic units of Miocene to Recent age. Analyses of seismic facies and isopach maps indicate that: (i) the palaeo‐Kushiro submarine canyon, which was ancestral to the Kushiro submarine canyon, was an aggradational levéed channel; and (ii) the palaeo‐Tokachi submarine channel was much larger than the present‐day channel and changed its course several times. Both the palaeo‐Kushiro submarine canyon and palaeo‐Tokachi submarine channel were fed predominantly by the ancestral Tokachi River, whereas the present‐day channels are no longer connected or quasi‐connected to the Tokachi River. Entrenchment of the Kushiro submarine canyon began in its distal reaches during the Early Pleistocene and propagated landward over time, which was possibly caused by base‐level fall (i.e. subsidence of the trench floor) or uplift of the forearc basin. Entrenchment of the upper part of the Kushiro submarine canyon began during the Middle Pleistocene, which may have been related to: (i) depositional progradation; (ii) uplift of the coastal area; or (iii) a change in source area from the ancestral Tokachi River to the Kushiro River.     

非连续子母块体理论模型研究(Ⅰ):基本理论

对节理岩体的非连续变形与破坏失稳的数值模拟分析是目前岩土力学与工程领域中的前沿课题之一,岩体由节理裂隙以及由其切割而成岩石块体组成,介于连续与非连续之间,由于其内部结构的不确定性造成了其宏观力学行为的复杂性。在讨论了岩体破坏失稳的研究尺度的基础上,从多尺度耦合的角度出发,同时考虑连续问题、非连续问题以及连续问题向非连续问题转化,提出了非连续子母块体理论模型。在非连续子母块体理论模型中,假定岩体与岩体内部被层理、节理等结构面所切割而成的岩石块体均属于宏观尺度,把宏观的岩石块体定义为母块体,用假想的人工节理将其继续分为若干更小的块体,更小的块体称为子块体,子块体为各向同性的均质体,并且在计算中不再发生破坏分解为更小的块体,作为计算分析的基本单元,子块体的研究尺度属于细观尺度。模型也重新定义了子块体间的接触,并将其划分为连续接触与非连续接触两种类型,连续接触在一定条件下发生破坏转化为非连续接触;采用增广拉格朗日乘子法对子块体间的接触进行处理,能够计算出准确的弹簧接触力。     

Numerical analysis of block stone displacements in ancient masonry structures: A new method to estimate historic ground motions

An innovative approach is presented, in which the discontinuous deformation analysis (DDA) method is used to estimate historic ground motions by back analysis of unique structural failures in archaeological sites. Two archaeological sites in Israel are investigated using this new approach and results are presented in terms of displacement evolution of selected structural elements in the studied masonry structure. The response of the structure is studied up to the point of incipient failure, in a mechanism similar to the one observed in the field. Structural response is found to be very sensitive to dynamic parameters of the loading function such as amplitude and frequency. Prior to back analysis of case studies, two validations are presented. Both compare the performance of DDA with analytical solutions and present strong agreement between the two. Using comprehensive sensitivity analyses, the most likely peak ground acceleration (PGA) and frequency that must have driven the observed block displacements are found for the two case studies—the Nabatean city of Mamshit and the medieval fortress of Nimrod in southern and northern Israel, respectively. It is found that horizontal peak ground accelerations (HPGA) of 0.5g and 1g were required to generate the observed deformations in Mamshit and Nimrod, respectively. Although these might seem too high, considering structural and topographic amplifications it is concluded that the analyses suggest ground motions of 0.2g at a frequency of 1.5 Hz for Mamshit and up to 0.4g at a frequency of 1 Hz for Nimrod. These values provide constraints on the seismic risk associated with these regions as appears in the local building code using a completely independent approach. Copyright © 2007 John Wiley & Sons, Ltd.     

气候序列不连续点的小波分析和混合聚类分析

不论是单纯的小波分析还是单纯的非参数分析法都有其局限性.当样本量很大时,Fisher算法的计算量也是很大的,而且它不具有谱系结构,因此,为了保证谱系结构以及较快的计算速度,作者以小波变换和聚类分析相结合的方法来研究各类不连续点,给出了相应的算法和实验分析,并且实验结果给出了各类不连续点位置的比较精确的估计.     

Couplex Benchmark Computations Obtained with the Software Toolbox UG

This paper describes the numerical results for the COUPLEX benchmark obtained with the simulation software UG using vertex centered finite volume and higher order discontinuous Galerkin schemes. Multigrid solvers on unstructured grids, local mesh refinement and parallel computation are employed to yield very accurate solutions. Since the full range of results required in the benchmarks is too large to be displayed in this paper we focus on the comparison of discretization schemes, assessment of numerical errors and the presentation of parallel computations.     

Modelling Two-Phase Incompressible Flow in Porous Media Using Mixed Hybrid and Discontinuous Finite Elements

In this paper, we present a numerical model for simulating two-phase (oil–water and air–water) incompressible and immiscible flow in porous media. The mathematical model which is based on a fractional flow formulation is formed of two nonlinear partial differential equations: a mean pressure equation and a water saturation equation. These two equations can be solved in a sequential manner. Two numerical methods are used to discretize the equations of the two-phase flow model: mixed hybrid finite elements are used to treat the pressure equation, h-based Richards' equation and the diffusion term in the saturation equation, the advection term in the saturation equation is treated with the discontinuous finite elements. We propose a better way to calculate the nonlinear coefficients contained in our equations on each element of the discretized domain. In heterogeneous porous media, the saturation becomes discontinuous at the interface between two porous media. We show in this paper how to use the capillary pressure–saturation relationship in order to handle the saturation jump in the mixed hybrid finite element method. The two-phase flow simulator is verified against analytical solutions for some flow problems treated by other authors.     

不连续岩体-拱坝系统的动力过程数值模拟方法

提出了一种分析三维不连续变形块体系统动力学过程和建坝的数值方法,该方法不但能够模拟位移不连续面在地震载荷作用下的张开和大尺度的摩擦滑移,而且还能给出不连续面上安全系数的分布和随时间的变化.通过有2个块体组成的系统验证了所提出的方法的正确性,同时用它分析了一个具有断层和伸缩缝的拱坝系统模型在地震载荷下的动力反应.     

A discontinuous Galerkin method for seismic soil–structure interaction analysis in the time domain

A technique for modeling transient wave propagation in unbounded media is extended and applied to seismic soil–structure interaction analysis in the time domain. The technique, based on the discontinuous Galerkin method, requires lower computational cost and less storage than the boundary element method, and the time‐stepping scheme resulting from Newmark's method in conjunction with the technique is unconditionally stable, allowing for efficient and robust time‐domain computations. To extend the technique to cases characterized by seismic excitation, the free‐field motion is used to compute effective forces, which are introduced on the boundary of the computational domain containing the structure and the soil in the vicinity of the structure. A numerical example on a dam–foundation system subjected to seismic excitation demonstrates the performance of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

川滇地区现今块体水平运动变形特征分析

依据川滇地区1991，1999和2001年GPS测量数据，利用DDA正分析基本原理，结合现有的地质资料，对该地区近年来的水平构造运动特征进行了模拟分析，认为川滇地区水平运动整体速率在1999～2001年期间明显减弱；川滇决体主压应变以顺时针旋转为主，方向为SE—SEE；其南部是运动变形相对剧烈、易于积累能量的区域；未来一两年内，该地区地震活动水平不会太高。     

岩土力学数值分析方法研究

随着岩土工程领域的不断扩展与延伸,岩土力学数值分析方法得到了迅速发展,出现了各种各样的数值分析方法。将岩土力学常用的数值分析方法分为连续变形分析方法与非连续变形分析方法两大类,对各种方法的研究现状和最新进展进行了评述,着重介绍了各种方法的优缺点,并提出了解决问题的思路、方法和建议。