印度与欧亚大陆碰撞构造变形数值分析

在全面考虑了大陆构造变形三重非线性特征的基础上，用动态有限元法研究了藏东对印度与欧亚大陆碰撞作用的大变形响应。数值模拟结果显示了明显的纵向收缩变形，变形主轴多具有顺时针转动的趋势，但横向位移在总体上并不显著。其最大值正位于红河和鲜水河断裂之间，这一结果论证了印度板块对欧亚大陆的碰撞可能是川滇菱形块体向南东运动的动力来源。红河、鲜水河、昆仑及阿尔金活断层在碰撞作用下发生了顺时针的旋转变形；而且，计算结果模拟出了红河断裂由左旋变为右旋，另３条断裂始终左旋的运动特征     

Estimating groundwater recharge for a freshwater lens in an arid region: Formative and stability assessment

The formation of subsurface freshwater lenses on top of brackish groundwater is a fascinating hydrologic phenomenon that creates groundwater supplies of great potential value in arid regions. Information on the recharge quantity and mechanism of these lenses is both scarce and uncertain. This study examines the formation and macroscale stability of the Rawdatain freshwater lens in Kuwait, for which significant pre-development data are available. The Rawdatain is a large (150 million m³) subsurface freshwater lens overlying brackish groundwater compared to the other freshwater lenses in the Arabian Peninsula. In this study, a three-dimensional (3D) density-dependent groundwater flow model is tested against the following data targets to estimate long-term diffuse and focused groundwater recharge: (i) groundwater head, (ii) total dissolved solids (TDS) groundwater concentration, (iii) volume and vertical thickness of stored groundwater of three different water quality TDS ranges (0–700, 700–1000 and 1000–2000 mg/L) and (iv) geometrical shape features of the lens along cross-sections. To better represent the spatial variation in TDS, six different recharge zones were assigned to allocate diffuse and focused recharge conditions. Twelve recharge rate scenarios, encompassing a wide range of feasible long-term average annual recharge values (200,000–5,000,000 m³/year), were tested against the multiple targets and compared with the groundwater age of the Rawdatain lens. Based on comparison with data targets, the long-term average annual recharge is estimated to be 500,000 m³/year. Scenarios of reduced recharge, which may occur due to changes in land-use or climate, demonstrate the extremely slow response of the lens, which is in agreement with the slow development and formation of the lens (>2,000 years). Within a 100-year time frame, a 50% reduction in annual recharge reduces the lens volumes by 21, 17 and 9% for the three water quality categories, respectively. This study demonstrates the stability of freshwater lenses in arid regions and also provides methodology for similar focused rainfall recharge freshwater lenses.     

Site effects at Euroseistest—I. Determination of the valley structure and confrontation of observations with 1D analysis

This paper describes the process of construction of the 2D model of Volvi's geological structure and results of empirical and theoretical approaches to the evaluation of site response at Euroseistest. The construction of the 2D model is based on a re-interpretation of the available geophysical and geotechnical data in an effort to improve the definition of the subsoil structure at Euroseistest in terms of the most important parameters needed to model site response. The results of this re-interpretation are compared with a previous published 2D model of the same alluvial valley. Different analysis of the measurements and different criteria in the synthesis of data have led to a different model, even if both studies had access to the same field measurements. This underscores the fact that a model results of an interpretation and is not uniquely determined by the data, no matter how detailed they are. The well known subsoil structure opened the possibility to correlate the geometry and the dynamic properties of the 2D model with the results of site response determined from a detailed analysis of two events in frequency and time domains and 1D numerical modeling. The study of site response shows the important effect of the lateral variations on the ground motion and suggests that the contribution of locally generated surface waves to the resonant peak may be important. In the case of Volvi's graben, the limitations of the 1D approximation to simulate ground motion under complex soil conditions in both frequency and time domains are also shown. This paper lays the ground for a companion article dealing with 2D site effects in this basin.     

Three-dimensional elastic wave scattering and diffraction due to a rigid cylinder embedded in an elastic medium by a point source

The formal solutions of displacement field to the problem of elastic wave scattering and diffraction due to an infinitely long rigid cylinder embedded in an infinite elastic medium by an impulsive point source have been obtained in the integral form. The integrals for the reflected and the diffracted waves both in the shadow zone and in the illuminated zone are evaluated asymptotically for the early time motion by the Reisdue-Cagniard method and the Saddle-point-Cagniard method.Numerical results of the diffractedP, S andPS waves at a fixed circum-distance from the surface of the rigid cylinder show noticeably that (1) the energy partition for the diffractedS wave is small in comparison with that for the diffractedP wave, (2) the wave form of the diffractedS wave is broader and more diffused than that of the diffractedP wave, (3) the direction of the radial motions of the diffractedP andS waves varies as a function of the observational point, and (4) the energy partition for the diffractedP wave is much smaller than that for the direct or the reflectedP waves.This paper has been presented at the 46th Annual International Meeting of Society of Exploration Geophysicists in Houston, Texas, Oct. 28, 1976.     

Multi-scale modeling of three-phase–three-component processes in heterogeneous porous media

Flow and transport processes in porous media occur on different spatial and temporal scales and may also be locally different. Additionally, the structure of the porous medium itself generally shows a high dependence on the spatial scale.     

ACTIVE system for monitoring volcanic activity: A case study of the Izu-Oshima Volcano,Central Japan

A system is proposed for the monitoring of changes in the underground structure of an active volcano over time by applying a transient electromagnetic method. The monitoring system is named ACTIVE, which stands for Array of Controlled Transient-electromagnetics for Imaging Volcano Edifice. The system consists of a transmitter dipole used to generate a controlled transient electromagnetic (EM) field and an array of receivers used to measure the vertical component of the transient magnetic field at various distances, with automatic operation of both units. In order to verify the performance of the proposed system, numerical and field experiments were carried out by application of the system to the Izu-Oshima volcano, where a remarkable change in the apparent DC resistivity over time had been detected in association with the eruption in 1986.     

The effect of interpolation methods in surface definition: an experimental study

Mappings of the earth surface and their representation in 3D (three‐dimensional) models are commonly used in most recent research. Modeling research, which starts with classical surveying methods, acquires new dimensions matching the modern technologies. 3D models of any object or earth surface can be used in much visual and scientific research. A digital model of the landscape is an important part within creation of geo‐information systems used in the public administration and in the commercial sphere. It is an important tool in applications such as geomorphology, hydrology, geology, cartography, ecology, mining etc. Values of volume in terrains that do not have regular geometric structure can be obtained more accurately by using 3D models of surfaces with respect to developing technology. Basic data of 3D models must indicate 3D coordinates of the surveyed object in the reference frame. Distribution and intensity of points are important factors in modeling earth surfaces. A minimum number of points is desired in defining an object in 3D. Interpolation methods employing different mathematical models are used to obtain 3D models of terrain surfaces. In this study, the effect of interpolation methods in defining a terrain surface is investigated. For this purpose, a uniform surface, hill‐shaped artificial object with a known volume is employed. The 3D surface and volume are calculated by using 12 different interpolation methods. Point distribution, point intensity and accuracy of point measurements are not considered. The same data set was used for all the interpolation methods. The interpolation methods are compared and evaluated based on the results. Copyright © 2007 John Wiley & Sons, Ltd.     

Seismic traveltime inversion of 3D velocity model with triangulated interfaces

Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks, which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block, and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization, which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models. The damped least squares method is employed in seismic traveltime inversion, which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.     

一阶弹性波方程数值模拟中的混合吸收边界条件

Liu和Sen（2010和2012）在地震波场数值模拟中提出一种混合吸收边界条件. 该方法具有计算量小、容易实现及吸收效果好等优点. 但现有的混合吸收边界条件是针对二阶位移-应力方程设计的,存在稳定性问题. 本文首先推导了两种速度-应力单程波方程：二阶Higdon单程波方程和一阶Higdon单程波方程. 进而提出基于一阶弹性波方程的混合吸收边界条件方法. 在内部区域和边界之间引入一个过渡区域,通过单程波与双程波方程平滑过渡来消除人工边界反射. 为了改善混合吸收边界条件的吸收效果和稳定性,我们采用了能同时吸收纵、横波反射的一阶单程波方程和与变量位置有关的加权系数. 为了验证混合吸收边界条件的有效性,将其与常规分裂完全匹配层（PML）方法进行了比较. 数值模拟结果表明,与PML边界条件相比,混合吸收边界条件在耗用更小计算时间和存储量的前提下,可以获得更好的吸收效果. 另外,本文提出的两种混合吸收边界条件中,混合一阶Higdon吸收边界条件具有更好的稳定性.