数字观测时代的全球三维结构与地震定位研究

概述了地震学观测研究取得的地球内部结构的一些基本认识，并介绍了在数字地震观测时代，地震层析成像研究取得的一些进展，如敏感性分析、结构和定位联合反演、非线性反演和非规则的参量化方法等。同时就地震学的传统定位问题，概述了利用地震走时进行地球深部结构反演所涉及的一些基本问题，即数据资料、模型参量化、反演计算和敏感性分析中的具体细节问题。中就国际地学界近年来利用三维结构开展地震定位研究的动向，从地震精定位在地震成核、精细地质结构刻画、核查和震源破裂过程研究等方面的应用研究角度，希望在中国能尽快推进相应的地震精度定位工作。     

A low‐dimensional hillslope‐based catchment model for layered groundwater flow

Despite the strong interaction between surface and subsurface waters, groundwater flow representation is often oversimplified in hydrological models. For instance, the interplay between local or shallow aquifers and deeper regional‐scale aquifers is typically neglected. In this work, a novel hillslope‐based catchment model for the simulation of combined shallow and deep groundwater flow is presented. The model consists of the hillslope‐storage Boussinesq (hsB) model representing shallow groundwater flow and an analytic element (AE) model representing deep regional groundwater flow. The component models are iteratively coupled via a leakage term based on Darcy's law, representing delayed recharge to the regional aquifer through a low conductivity layer. Simulations on synthetic single hillslopes and on a two‐hillslope open‐book catchment are presented, and the results are compared against a benchmark three‐dimensional Richards equation model. The impact of hydraulic conductivity, hillslope plan geometry (uniform, convergent, divergent), and hillslope inclination (0.2%, 5%, and 30%) under drainage and recharge conditions are examined. On the single hillslopes, good matches for heads, hydrographs, and exchange fluxes are generally obtained, with the most significant differences in outflows and heads observed for the 30% slope and for hillslopes with convergent geometry. On the open‐book catchment, cumulative outflows are overestimated by 1–4%. Heads in the confined and unconfined aquifers are adequately reproduced throughout the catchment, whereas exchange fluxes are found to be very sensitive to the hillslope drainable porosity. The new model is highly efficient computationally compared to the benchmark model. The coupled hsB/AE model represents an alternative to commonly used groundwater flow representations in hydrological models, of particular appeal when surface–subsurface exchanges, local aquifer–regional aquifer interactions, and low flows play a key role in a watershed's dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Flow,bed topography,grain size and sedimentary structure in open channel bends: A three-dimensional model

Bed topography and grain size are predicted for steady, uniform flow in circular bends by consideration of the balance of fluid, gravity and frictional forces acting on bed load particles. Uniform flow pattern is adequately described by conventional hydraulic equations, with bed shear defined as that effectively acting on bed load grains. This analysis is used as a basis to predict bed topography and grain size for steady, non-uniform flow in non-circular bends (represented by a ‘sine-generated’ curve). The non-uniform flow pattern is calculated using the method of Engelund (1974a). Equilibrium bed form, hence sedimentary structure, is found by comparison of existing flow conditions with one of the schemes describing the hydraulic stability limits of the various bed forms. The model was compared with bankfull flow observations from a channel bend on the River South Esk, Scotland. Theoretical bed topography and velocity distribution were very close to the observed data. However, bed shear stress showed only a broad agreement, probably because of the use a constant friction coefficient value. Mean grain size distribution showed good agreement, but theory did not account adequately for gravel sizes in the talweg region and on the upstream, inner part of the bar, possibly due to theoretical underestimation of effective bed shear. Bed form and sedimentary structure are predicted well using the familiar stream power-grain size scheme. The behaviour of the model under unsteady uniform flow conditions in circular bends was analyzed, and suggests that any variation of grain size and bed topography with stage is likely to be limited to deeper parts of the channel.     

A stochastic model for initial movement of sand grains by wind

A stochastic model for entrainment of sand grains by wind is presented through analysis of the forces exerted on a single spherical grain, coupled with fluctuations of wind velocity and the change in grain position on the surface. The structure of the stochastic model is consistent with experimental data in the literature. The probability of initial motion increases first, and then decreases, with grain size. It reaches a maximum at diameters of about 0·9 mm. Some sand grains are still in motion at less than the conventional threshold velocity, even at very low velocities. The probability of sand grain movement reaches unity at twice the conventional threshold velocity. Considerable discrepancies amongst conventional threshold formulae may result from the different probabilities of initial movement implied in these formulae. Copyright © 2008 John Wiley & Sons, Ltd.     

A two‐dimensional model of hydraulic performance of stormwater infiltration systems

Stormwater infiltration systems are a popular method for urban stormwater control. They are often designed using an assumption of one‐dimensional saturated outflow, although this is not very accurate for many typical designs where two‐dimensional (2D) flows into unsaturated soils occur. Available 2D variably saturated flow models are not commonly used for design because of their complexity and difficulties with the required boundary conditions. A purpose‐built stormwater infiltration system model was thus developed for the simulation of 2D flow from a porous storage. The model combines a soil moisture–based model for unsaturated soils with a ponded storage model and uses a wetting front‐tracking approach for saturated flows. The model represents the main physical processes while minimizing input data requirements. The model was calibrated and validated using data from laboratory 2D stormwater infiltration trench experiments. Calibrations were undertaken using five different combinations of calibration data to examine calibration data requirements. It was found that storage water levels could be satisfactorily predicted using parameters calibrated with either data from laboratory soils tests or observed water level data, whereas the prediction of soil moistures was improved through the addition of observed soil moisture data to the calibration data set. Copyright © 2012 John Wiley & Sons, Ltd.     

三维模式约束的电离层掩星反演方法

目前电离层掩星数据反演是基于电离层电子密度分布局部球对称近似的Abel反演方法,实际电离层的非球对称性会给电子密度反演结果带来误差.本文研究利用三维电离层模式来提供电子密度水平变化的先验信息约束电离层掩星反演的方法,即三维模式约束法;并将该方法应用于模拟掩星观测数据和实测掩星数据的反演.模拟观测数据的反演结果表明,与Abel反演方法相比,三维模式约束法能够减小反演误差.采用IRI2001模式作为约束,对COSMIC电离层掩星实测数据反演,将反演结果与全球的垂测仪数据进行比较,结果表明,三维模式约束法和Abel反演方法都能很好地反演电离层掩星.     

Generation of three dimensional flow fields for statistically anisotropic heterogeneous porous media

A methodology for generating three dimensional (3D) flow fields for statistically anisotropic heterogeneous porous media is presented and demonstrated. The simulated flow fields are shown to exhibit the input spatial correlation structure and observe mass continuity. Sample flow fields are presented in the form of cross sectional slices of the 3D formation. These cross sections demonstrate visually the characteristics of subsurface flow. The method was found to be faster than traditional techniques in terms of its computational requirements. Given this method, it is possible to generate the large number of realizations of a velocity field necessary to compute high order statistics in transport problems.     

A laboratory model of surface crust formation and disruption on lava flows through non-uniform channels

Crust formation on basaltic lava flows dictates conditions of both flow cooling and emplacement. For this reason, flow histories are dramatically different depending on whether lava is transported through enclosed lava tubes or through open channels. Recent analog experiments in straight uniform channels (Griffiths et al. J Fluid Mech 496:33–62, 2003) have demonstrated that tube flow, dictated by a stationary surface crust, can be distinguished from a mobile crust regime, where a central solid crust is separated from channel walls by crust-free shear zones, by a simple dimensionless parameter ϑ, such that ϑ<25 produces tube flow and ϑ>25 describes the mobile crust regime. ϑ combines a previously determined parameter ψ, which describes the balance between the formation rate of surface solid and the shear strain that disrupts the solid crust, with the effects of thermal convection (described by the Rayleigh number Ra).Here we explore ways in which ϑ can be used to describe the behavior of basaltic lava channels. To do this we have extended the experimental approach to examine the effects of channel irregularities (expansions, contractions, sinuosity, and bottom roughness) on crust formation and disruption. We find that such changes affect local flow behavior and can thus change channel values of ϑ. For example, gradual widening of a channel results in a decrease in flow velocity that causes a decrease in ϑ and may allow a down-flow transition from the mobile crust to the tube regime. In contrast, narrowing of the channel causes an increase in flow velocity (increasing ϑ), thus inhibiting tube formation.We also quantify the fraction of surface covered by crust in the mobile crust regime. In shallow channels, variations in crust width (d _c) with channel width (W) are predicted to follow d _c∼W ^5/3. Analysis of channelized lava flows in Hawaii shows crustal coverage consistent with this theoretical result along gradually widening or narrowing channel reaches. An additional control on crustal coverage in both laboratory and basaltic flows is disruption of surface crust because of flow acceleration through constrictions, around bends, and over breaks in slope. Crustal breakage increases local rates of cooling and may cause local blockage of the channel, if crusts rotate and jam in narrow channel reaches. Together these observations illustrate the importance of both flow conditions and channel geometry on surface crust development and thus, by extension, on rates and mechanisms of flow cooling. Moreover, we note that this type of analysis could be easily extended through combined use of FLIR and LiDAR imaging to measure crustal coverage and channel geometry directly.Editorial responsibility: A. Harris     

台湾地区地壳形变的弹性块体位错模型

在经典的非震形变位错模型中,地壳形变被认为是活动块体刚性运动和上部断层锁定影响的叠加,本文对此模型进行了改进: (1) 用活动块体整体运动和内部线性应变、旋转的贡献代替活动块体刚性运动的贡献;(2) 用分层介质地壳模型代替半无限介质模型计算断层锁定的影响. 利用改进后的非震形变位错模型,拟合了台湾地区1990～1995年间GPS观测资料. 结果显示,在东部海岸山脉区,约有30 mm·a-1的汇聚率被奇美断层消耗掉,运动速度从奇美断层向北迅速衰减. 在西部平原地区,南部断层是岛内锁定最为强烈的断层,该地区相应的也是史上灾害性地震多发的地区. 根据反演结果计算出的应变率与旋转率分布与前人结果在大部分地区一致,主应变率场显示台湾大部分地区存在近NW-SE方向的主压应变,主压应变方向呈扇形分布. 旋转率场显示台湾东部和南部地区存在着逆时针旋转率,而西部和北部地区则为顺时针旋转率.     

球坐标系下多震相走时三参数同时反演成像

球坐标系下多震相走时三参数(速度、震源位置和反射界面)同时反演需要解决两个关键问题:(1)球坐标系下3D速度模型中多次透射、反射(折射)及转换波精确、快速的射线追踪;(2)同时反演时三种不同参数间的强耦合问题.为此,我们将直角坐标系下分区多步不规则最短路径算法推广至球坐标系中,进行区域或者全球尺度的多震相射线追踪.然后将其与适合多参数同时反演的子空间算法相结合,形成一种球坐标系下联合多震相走时三参数同时反演的方法技术.与双参数(速度和反射界面或速度和震源位置)同时反演的数值模拟对比分析显示:三参数与双参数的同时反演结果大体接近,并且它们对到时数据中可容许的随机噪声不太敏感.结果说明本文中的同时反演成像为一种提高成像分辨率,同时反演速度、震源位置和反射界面的有效方法.     

Numerical calculation for bed variation in compound-meandering channel using depth integrated model without assumption of shallow water flow

In a compound meandering channel, patterns of flow structures and bed variations change with increasing water depth owing to complex momentum exchange between high-velocity flow in a main channel and low-velocity flows in flood plains. We have developed a new quasi-three-dimensional model without the shallow water assumption, i.e., hydrostatic pressure distribution; our method is known as the general bottom velocity computation (BVC) method. In this method, a set of depth-integrated equations, including depth-integrated momentum and vorticity equations, are prepared for evaluating bottom velocity and vertical velocity distributions. The objective of this study is to develop a bed variation calculation method for both single and compound meandering channels by using the BVC method coupled with a sediment transport model. This paper shows that the BVC method can reproduce the pattern change of bed variation in a compound meandering channel flow with increasing relative depth. The variation in sediment transport rate due to overbank flow is explained by experimental and computational results.     

唐山响堂三维场地影响观测台阵

１９９４年７月中国地震局工程力学研究所在唐山余震区响堂镇建成了我国第一个三维场地影响观测台阵。该台阵目前有四个测点,分别布设在基岩地表、土层地表、地下１７ｍ和地下３２ｍ处。它安装了分辨率为１６位的井下数字观测系统,主机和从机同步运行,系统时钟采用Ｏｍｅｇａ导航信号自动校对（精度１ｍｓ）,地震数据采用固态方式（ＣＭＯＳ）存储,每个测点均布设一组 量力平衡式加速度计。该台阵自投入运行至１９９７年１２月     

A depth-averaged two-phase model for debris flows over fixed beds

A depth-averaged two-phase model is proposed for debris flows over fixed beds, explicitly incorporating interphase and particle-particle interactions, fluid and solid fluctuations and multi grain sizes. A first-order model based on the kinetic theory of granular flows is employed to determine the stresses due to solid fluctuations, while the turbulent kinetic energy - dissipation rate model is used to determine the stresses from fluid fluctuations. A well-balanced numerical algorithm is applied to solve the governing equations. The present model is benchmarked against USGS experimental debris flows over fixed beds. Incorporating the stresses due to fluid and solid fluctuations and properly estimating the bed shear stresses are shown to be crucial for reproducing the debris flows. Longitudinal particle segregation is resolved, demonstrating coarser sediments around the fronts and finer grains trailing the head. Based on extended modeling exercises, debris flow efficiency is shown to increase with initial volume, which is underpinned by observed datasets.     

BED SHEAR IN EVOLVING SCOUR AT A CYLINDER:A THEORETICAL APPROACH

ＢＥＤＳＨＥＡＲＩＮＥＶＯＬＶＩＮＧＳＣＯＵＲＡＴＡＣＹＬＩＮＤＥＲ：ＡＴＨＥＯＲＥＴＩＣＡＬＡＰＰＲＯＡＣＨＳｕｂｈａｓｉｓｈＤＥＹ１ａｎｄＳｒｉｊｉｂＫ．ＫＡＲ２Ａｂｓｔｒａｃｔ：Ｔｈｅｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｏｆｓｈｅａｒｓｔｒｅｓｓ...     

A constitutive model for water flow in unsaturated fractured rocks

A conceptual model for describing effective saturation in fractured hard rock is presented. The fracture network and the rock matrix are considered as an equivalent continuum medium where each fracture is conceptualized as a porous medium of granular structure and the rock matrix is assumed to be impermeable. The proposed model is based on the representation of a rough‐walled fracture by an equivalent porous medium, which is described using classical constitutive models. A simple closed‐form equation for the effective saturation is obtained when the van Genuchten model is used to describe saturation inside fractures and fractal laws are assumed for both aperture and number of fractures. The relative hydraulic conductivity for the fractured rock is predicted from a simple relation derived by Liu and Bodvarsson. The proposed constitutive model contains three independent parameters, which may be obtained by fitting the proposed effective saturation curve to experimental data. Two of the model parameters have physical meaning and can be identified with the reciprocal of the air entry pressure values in the fractures of minimum and maximum apertures. Effective saturation and relative hydraulic conductivity curves match fairly well the simulated constitutive relations obtained by Liu and Bodvarsson. Copyright © 2008 John Wiley & Sons, Ltd.     

A two-dimensional dam-break flood plain model

A simple two-dimensional dam-break model is developed for flood plain study purposes. Both a finite difference grid and an irregular triangle element integrated finite difference formulation are presented. The governing flow equations are approximately solved as a diffusion model coupled to the equation of continuity. Application of the model to a hypothetical dam-break study indicates that the approach can be used to predict a two-dimensional dam-break flood plain over a broad, flat plain more accurately than a one-dimensional model, especially when the flow can break-out of the main channel and then return to the channel at other downstream reaches.     

Analytical model for fully three‐dimensional radial dispersion in a finite‐thickness aquifer

Analytical solutions for contaminant transport in a non‐uniform flow filed are very difficult and relatively rare in subsurface hydrology. The difficulty is because of the fact that velocity vector in the non‐uniform flow field is space‐dependent rather than constant. In this study, an analytical model is presented for describing the three‐dimensional contaminant transport from an area source in a radial flow field which is a simplest case of the non‐uniform flow. The development of the analytical model is achieved by coupling the power series technique, the Laplace transform and the two finite Fourier cosine transform. The developed analytical model is examined by comparing with the Laplace transform finite difference (LTFD) solution. Excellent agreements between the developed analytical model and the numerical model certificate the accuracy of the developed model. The developed model can evaluate solution for Peclet number up to 100. Moreover, the mathematical behaviours of the developed solution are also studied. More specifically, a hypothetical convergent flow tracer test is considered as an illustrative example to demonstrate the three‐dimensional concentration distribution in a radial flow field. The developed model can serve as benchmark to check the more comprehensive three‐dimensional numerical solutions describing non‐uniform flow contaminant transport. Copyright © 2009 John Wiley & Sons, Ltd.     

A low‐dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes

Hillslopes have complex three‐dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope‐storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three‐dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr–Coulomb failure law. The resulting hillslope‐storage Boussinesq stability model (HSB‐SM) is able to simulate rain‐induced shallow landsliding on hillslopes with non‐constant bedrock slope and non‐parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB‐SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd.     

A nonlinear polar coordinate shallow water model for tsunami computation along North Sumatra and Penang Island

A nonlinear shallow water model in cylindrical polar coordinate system is developed, using an explicit finite difference scheme with a very fine resolution, to compute different aspects of tsunami at North Sumatra and the adjacent island Simeulue in Indonesia, and the Penang Island in Peninsular Malaysia. The pole of the frame is placed on the mainland of Penang (100.5°E) and the model area extends up to the west of Sumatra (87.5°E). The model is applied to simulate the propagation of tsunami wave towards North Sumatra, Simeulue and Penang Islands associated with Indonesian tsunami of 26 December 2004. The model is also applied to compute water levels along the coastal belts of those islands. Computed and observed water level data are found to be in good agreement and North Sumatra is found to be vulnerable for very high surges. The computed and observed arrival times of high surges are also in reasonable agreement everywhere. Further studies are carried out to investigate the effect of convective terms and it is found that their effects are insignificant in tsunami propagation and weakly significant for wave amplitude very near to the coast.