双时间步法在二维浅水方程求解中的应用

针对传统算法效率低的问题,将隐式双时间步法应用于求解二维浅水方程,建立了非结构网格下高效的有限体积模型。在应用双时间步法时,虚拟时间层中的定常问题采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法进行迭代求解。通过模拟计算4个典型算例以及与传统显式算法进行比较,对模型精度、效率及处理实际问题能力进行检验,分析了时间步长、内迭代次数对模型性能的影响。结果表明,双时间步法放宽了稳定性对时间步长的限制,时间步长可取到显式格式10倍以上,计算耗时减少了50%以上,模型具有良好精度与适应性,具有较好的推广应用价值。     

辅助微分方程完全匹配层在声波方程数值模拟中的应用

在地震波数值模拟中,为提高算法精度,需要使用高阶时间更新格式,而普通的非分裂完全匹配层（PML）吸收边界局限于低阶时间格式。辅助微分方程完全匹配层（ADE PML）是一种可以适应任意阶时间格式的非分裂完全匹配层技术,且可以直接应用复频移拉伸算子以提高PML在高角度入射时的效果。作者将ADE PML应用于声波方程四阶Runge Kutta时间格式的数值模拟中,对其吸收效能进行了检验。数值模拟表明,复频移ADE PML在高角度入射时表现优于非复频移ADE PML。另外,不同辅助变量更新格式的吸收效果存在微小差异,显格式下计算结果与解析解吻合较好。长时间能量衰减计算表明ADE PML可以稳定至2 × 10⁵时间步。     

Palaeotopography of the northern portion of the Telychian (Silurian) turbidite basin in central Wales

New palaeocurrent data from the country northeast and southwest of Corris (central Wales) indicate that Telychian sandstones in the Devil's Bridge Formation hereabouts were transported from both the northeast and northwest, consistent with topographical control by down‐to‐southeast movement on the Bala Lineament. This conclusion is supported by considerations of thickness variation and by evidence for listric detachment faulting in the formation at Bwlch y Groes. Towards Plynlimon, flow is broadly southerly and suggests confinement by the northward extension of the Bronnant Fault. A facies and palaeocurrent map is presented for the utilis sub‐Biozone of the Telychian which links with work by the British Geological Survey between Aberyswyth and Rhayader where, by contrast, palaeoflow is to the NNW. Copyright © 2005 John Wiley & Sons, Ltd.     

Effective block diagonal preconditioners for Biot's consolidation equations in piled‐raft foundations

The finite element (FE) simulation of large‐scale soil–structure interaction problems (e.g. piled‐raft, tunnelling, and excavation) typically involves structural and geomaterials with significant differences in stiffness and permeability. The symmetric quasi‐minimal residual solver coupled with recently developed generalized Jacobi, modified symmetric successive over‐relaxation (MSSOR), or standard incomplete LU factorization (ILU) preconditioners can be ineffective for this class of problems. Inexact block diagonal preconditioners that are inexpensive approximations of the theoretical form are systematically evaluated for mitigating the coupled adverse effects because of such heterogeneous material properties (stiffness and permeability) and because of the percentage of the structural component in the system in piled‐raft foundations. Such mitigation led the proposed preconditioners to offer a significant saving in runtime (up to more than 10 times faster) in comparison with generalized Jacobi, modified symmetric successive over‐relaxation, and ILU preconditioners in simulating piled‐raft foundations. Copyright © 2012 John Wiley & Sons, Ltd.     

基于反褶积算法的地震干涉技术被动源成像

地震干涉技术可以将任意2个检波器接收到的数据合成为在若干检波器之间传播的波,就好像其中的一个检波器作为一个虚拟震源来发挥作用。它可以从混沌无序的地震信号中发现有用信息,从地震噪声中提取有用信号以此推断地震波穿过介质的地质构造。基于反褶积算法,对其理论公式进行了较详细的推导,实现被动源地震干涉成像,证明了反褶积算法的可行性;并将其结果与互相关算法的结果进行对比,分析了2种方法在信噪比和分辨率方面的差异。数值计算表明,反褶积算法的纵向分辨率比互相关算法的高。对其进行的加噪试算表明,震源叠加后的反褶积算法呈现出高信噪比的特点。     

地表的分形测量及其大地构造学意义

以湖北红安地区为例,采用投影覆盖法(projectivecoveringmethod)对地表进行了二维分形测量,结果表明,地表面积具有双分形(bifractal)关系,即具有小尺度的结构分形(texturalfractal)和大尺度的构造分形(structuralfractal),分叉点(breakpoint)的尺度为3610m,分维值都在2～3之间且结构分维值大于其构造分维值。可见,地表形态具有分形性质,分维值可以指示地表形态的复杂程度。构造分维值可作为构造活动强度的一个指标,可为大地构造单元的划分提供定量依据。复杂地表形态主要是由构造活动(内营力作用)和各种复杂表生地质作用(外营力作用)引起的,前者主要控制大尺度的地形起伏,后者则塑造小尺度的地表形态。地表分维值可以指示地表的发育成熟度,该地区小尺度的结构分维值大于大尺度的构造分维值表明其处于地表形态的发育晚期。此外,地表的分形尺度可以来用确定构造活动尺度,从而指导构造地质与找矿勘探研究。     

Genesis of cristobalite and tridymite at low temperatures

It is suggested that an explanation of the observed metastable formation of cristobalite and tridymite in the stability field of quartz lies in the topology of these silica polymorphs and the resulting energies at the time of formation.     

Profondeur d'extraction racinaire et signature isotopique de l'eau prélevée par les racines des couverts végétaux

We seek to identify the depth to which water is extracted by the roots in the soil. Indeed, in an isotopic steady-state condition of leaf water, transpiration introduces into the atmosphere a vapour whose isotopic signature is identical to that of root water. In the isotopic models of atmospheric general circulation, it is classically allowed that the signature of transpiration belongs to the meteoric water line. This supposes that the water taken by the roots has escaped with the evaporation of the soil and comes thus from the deep layers of the soil. At the time of experimentation carried out on maize plants (Nemours, Seine-et-Marne, France), this extraction depth was inferred from the comparison between the signature of the water measured on the level of the first internode of the stems of the plants and the isotopic profile of water in the soil. When the flow of transpiration reaches a maximum value, the plant uptakes water resulting from precipitations and which preserves its non-evaporating character after having quickly infiltrated in the deep layers of the soil. This relates to only 55% of the flux transpired by the canopy, the remainder presenting an evaporating character more or less marked according to ambient conditions. This experiment invalidates the classical hypothesis used in isotopic models of general atmospheric circulation in temperate regions. In fact, only half the amount of water vapour transpired by the canopy during the day presents a signature similar to that of the rainwater sampled in deep soil layers. To cite this article: Z. Boujamlaoui et al., C. R. Geoscience 337 (2005).     

考虑Hansbo渗流的二维Biot固结有限元分析

为进一步深入研究弹性饱和黏性土地基的二维固结机制,引入Hansbo渗流方程描述固结过程中的非达西渗流,修正Biot二维固结方程。基于加权残数法,给出相应的有限元数值求解格式。通过和饱和黏性土一维非达西渗流固结理论有限体积法数值结果的对比,证明数值计算方法的有效性。在此基础上,探讨Hansbo渗流参数对二维地基固结进程的影响。计算结果表明,在固结初期,Hansbo渗流将增强Mandel-Cryer效应,增大孔压的峰值,并延长孔隙水压力达到峰值的时间;在固结中后期,整个土层存在孔隙水压力滞后现象。同时,Hansbo渗流将阻碍地基沉降的发展。而且,上述影响会随着Hansbo渗流参数的增大而更加明显。     

贾汪煤田为何未能向北东继续延伸

贾汪煤田的基本构造形态为贾汪向斜。通过对贾汪向斜北东延长段原地系统和外来系统的分析研究，认为在贾汪石炭——二叠纪煤田北界伏山前断裂的北东盘上，受构造应力的影响，奥陶系沿着其与寒武系界面间的软弱层被推覆，在原地代之以大面积分布的是青白口——震旦系及部分寒武系而非石炭——二叠系，从而造成了贾汪煤田没有向北东继续延伸。