鲁东造山带榴辉岩变质作用特征及其动力学机制

岩相学研究表明,广泛发育鲁东造山逞中的榴辉岩可以划分为四种类型：第一类为含柯石英及其假象榴辉岩;么二类为含蓝晶石 、黝帘石及多硅白云母等榴辉岩;第三类为石榴石－绿辉石－石英组合榴辉岩;第四类为（角闪石）石榴石－辉及有关岩石（非典型“榴辉岩”）。这四类榴辉岩峰期前进变质矿物共生组合及Ｐ－Ｔ条件估算结果均显示它们应属于高温与超高压变质作用产物,其变质作用ＰＴｔ轨迹表现为顺时针演化趋势,反映出板块俯冲碰     

A spatially and temporally adaptive solution of Richards’ equation

Efficient, robust simulation of groundwater flow in the unsaturated zone remains computationally expensive, especially for problems characterized by sharp fronts in both space and time. Standard approaches that employ uniform spatial and temporal discretizations for the numerical solution of these problems lead to inefficient and expensive simulations. In this work, we solve Richards’ equation using adaptive methods in both space and time. Spatial adaption is based upon a coarse grid solve and a gradient error indicator using a fixed-order approximation. Temporal adaption is accomplished using variable order, variable step size approximations based upon the backward difference formulas up to fifth order. Since the advantages of similar adaptive methods in time are now established, we evaluate our method by comparison with a uniform spatial discretization that is adaptive in time for four different one-dimensional test problems. The numerical results demonstrate that the proposed method provides a robust and efficient alternative to standard approaches for simulating variably saturated flow in one spatial dimension.     

Diffusive wave solutions for open channel flows with uniform and concentrated lateral inflow

The diffusive wave equation with inhomogeneous terms representing hydraulics with uniform or concentrated lateral inflow into a river is theoretically investigated in the current paper. All the solutions have been systematically expressed in a unified form in terms of response function or so called K-function. The integration of K-function obtained by using Laplace transform becomes S-function, which is examined in detail to improve the understanding of flood routing characters. The backwater effects usually resulting in the discharge reductions and water surface elevations upstream due to both the downstream boundary and lateral inflow are analyzed. With a pulse discharge in upstream boundary inflow, downstream boundary outflow and lateral inflow respectively, hydrographs of a channel are routed by using the S-functions. Moreover, the comparisons of hydrographs in infinite, semi-infinite and finite channels are pursued to exhibit the different backwater effects due to a concentrated lateral inflow for various channel types.     

Yielding oscillator subjected to simple pulse waveforms: numerical analysis & closed‐form solutions

Numerical and analytical solutions are presented for the elastic and inelastic response of single‐degree‐of‐freedom yielding oscillators to idealized ground acceleration pulses. These motions are typical of near‐fault earthquake recordings generated by forward rupture directivity and may inflict damage in the absence of substantial structural strength and ductility capacity. Four basic pulse waveforms are examined: (1) triangular; (2) sinusoidal; (3) exponential; and (4) rectangular. In the first part of the article, a numerical study is presented of the effect of oscillator period, strength, damping, post‐yielding stiffness and number of excitation cycles, on inelastic response. Results are presented in the form of dimensionless graphs and regression formulas that elucidate the salient features of the problem. It is shown that conventional R–µ relations may significantly underestimate ductility demand imposed by near‐fault motions. The second part of the article concentrates on elastic‐perfectly plastic oscillators. Closed‐form solutions are derived for post‐yielding response and associated ductility demand. It is shown that all three ground motion histories (i.e. acceleration, velocity, and displacement) control oscillator response—contrary to the widespread view that ground velocity alone is of leading importance. The derived solutions provide insight on the physics of inelastic response, which is often obscured by the complexity of numerical algorithms and actual earthquake motions. The model is evaluated against numerical results from near‐field recordings. A case study is presented. Copyright © 2006 John Wiley & Sons, Ltd.     

进化策略算法在地面沉降量计算参数反演中的应用

多数情况下,评估的研究区缺乏多年积累的地下水位及开采量资料。因此,通常采用土力学模型,根据太沙基固结理论,进行地面沉降量预测。在地面沉降量计算中,土层变形参数值的选取直接影响计算结果的精确性。深层含水层埋深通常大于100m。通过现场勘探,采取土样进行岩土物理力学性质试验获取参数值,成本过于昂贵。通过参考区域土层物理力学性质资料,结合工程经验获取参数值,预测结果的精度得不到保证。针对上述情况,作者提出用非线性优化算法-进化策略算法,根据实测地面沉降量反演土层变形参数值。通过这种方法确定参数值,既节约成本,又保证了计算结果的可靠性。进化策略算法通过模拟生物遗传及进化过程,利用转移概率来帮助指导搜索。搜索结果不依赖于初始点的选择,对于求解全局最优解有很强的鲁棒性。作者将进化策略算法用于某一工程实例土层变形参数的反演,结果表明了该算法的可行性及稳健性,值得在工程实践中推广应用。     

构造控制型天然气水合物矿藏及其特征

构造环境是天然气水合物富集成藏的重要控制因素,增生楔、断裂体系、褶皱、(泥)底辟、滑塌等特殊构造体是影响天然气水合物成藏的主要地质载体。通过对这些特殊构造体与天然气水合物成藏关系的研究,结合流体活动对水合物形成的影响,总结出陆缘地区有增生楔型、盆缘斜坡型、埋藏背斜型、断褶型、滑塌型及底辟型等六类构造控制型水合物矿藏。南海位于欧亚板块、太平洋板块及印澳板块的交汇处,早期为活动陆缘,晚期演化为被动陆缘,其构造活动具有早期张裂、后期挤压的特点,这既不同于被动陆缘,也有别于活动陆缘,可视为“复合型”大陆边缘,兼具了“被动陆缘沉积速率高、活动陆缘构造活跃”的优点,从而形成了“增生楔型、断褶型、底辟型、滑塌型、盆缘斜坡型”等多种构造控制型水合物矿藏,是“复合型”大陆边缘水合物成藏地质模式的典型代表。     

Study on Formation Mechanisms of Heavy Rainfall Within the Meiyu Along the Mid-Lower Yangtze River and Theories and Methods of Their Detection and Prediction

As the project of National Key Basic Research Development Program: Research on Formation Mechanisms and Predictive Theories of Major Weather Disasters in China has been fulfilled by 5-yr efforts of Chinese scientists, achieving results of great significance are as follows: 1) development of multi-scale physical models for Meiyu frontal heavy rainfall based on a range of real-time observations; 2) construction of synoptic models for such heavy rainfall; 3) the Meiyu front found to consist of multi-scale systems that represent a subtropical front, which shears structural features of an extratropical front and ITCZ, displaying sometimes a bi-front feature in the mid-lower Yangtze Basin (MLYB). The positive feedback between pre-frontal wet physical processes and over-front strong convective activities as well as interactions among multi-scale systems of the Meiyu front act as the important mechanism for the maintenance and development of the Meiyu front; 4) proposal of theories and methods for quantitative retrieval of multiple mesoscale torrential rains from satellite remote sensings, leading to a line of products; 5) investigation of applicable theories and techniques for retrieving the heavy rainfall system's 3D structure from dual-Doppler synchronous detectings; and 6) development of a system for meso heavy rainfall numerical prediction models with a 3D variational data assimilation scheme included, a tool that played an active role in flood combating and relief activities over the Huaihe River Basin (HRB) in 2003.     

2003年11月13日岷县ML5.5级地震震源机制的矩张量反演

利用2003年11月13日岷县5.5级地震震中周围7个台站的数字记录反演了本次地震的矩张量,对所得结果与用P波初动符号的结果进行了比较,二者具有较好的一致性。结合该地区的地震分布和本次地震的余震分布,对可能的发震断层进行了讨论。     

Scattering of Rayleigh waves by a shallow circular alluvial valley： High-frequency solution

Introduction The effect of local site conditions on wave propagation is one of the most attractive topics in engineering seismology. It may be resolved by either a numerical method or an analytical method. Numerical methods include the finite difference, finite element, boundary element method, etc. The analytical method is the wave function expansion method. The advantage of these numerical methods is that they can be applied to local inhomogeneity of arbitrary shapes, but analytical method …     

A tectonic stress map of China and adjacent regions from earthquakes occurred during 2000～2004

Earthquake focal mechanism solutions provide the basic information about the present-day regional tectonics stress distribution, which controls the activities of crustal faults. Therefore, continued efforts for updating the da-tabase of earthquake focal mechanism solutions are quite valuable and important. The international lithosphere program initiated the ″world stress map project″ (Zoback, 1992) since 1980′s. During this project, worldwide studies and measurements of present-day regional…