The forward problem of electromagnetic induction: accurate finite-difference approximations for two-dimensional discrete boundaries with arbitrary geometry

We describe finite-difference approximations to the equations of 2-D electromagnetic induction that permit discrete boundaries to have arbitrary geometrical relationships to the nodes. This allows finite-difference modelling with the flexibility normally ascribed to finite-element modelling. Accuracy is demonstrated by comparison with finite-element computations. We also show that related approximations lead to substantially improved accuracy in regions of steep, but not discontinuous, conductivity gradient.     

Planform of convection with strongly temperature-dependent viscosity

Abstract

This paper experimentally investigates the convective planform near critical in a fluid layer whose temperature-dependent viscosity varies from top to bottom by up to a factor of 1500. Convection occurs in three different planforms: rolls, hexagons and squares. The square planform, which appears only for fluids with viscosity variation greater than about 50, replaces the hexagonal convection pattern as the Rayleigh number increases much above critical. The large amplitude of hexagonal convection with strong viscosity variation precludes studying the hexagon-square transition with perturbation methods of the type used to study the hexagon-roll transitions at smaller viscosity variations.     

Finite element analysis of excavation

The finite element method has often been used to simulate excavation. When the soil is linearly elastic, the results of excavation should be independent of the number of stages in the excavation process, and lack of such independence indicates an incorrect procedure. The simple direct method described in this paper provides the required independence in the case of linearly elastic materials, and hence can be used for multi-stage excavation in non-linear problems without excessive errors. However methods whose errors increase with the number of stages of excavation are quite unsuitable for non-linear problems. Alternative methods of analysis, errors arising from the inability of the elements to model adequately the stress gradients near the toe of the excavation and excavation adjacent to a diaphragm wall are discussed.     

西藏日喀则地区德村-昂仁蛇绿岩内基性岩的元素与Sr-Nd-Pb同位素地球化学及其揭示的特提斯地幔域特征

系统研究了西藏雅鲁藏布江蛇绿岩带中部日喀则地区德村、吉丁和昂仁蛇绿岩中基性岩石的元素与 Sr-Nd-Pb 同位素地球化学特征。这些基性岩石,包括玄武岩、辉长岩和辉绿岩,属于低钾拉斑玄武岩系列,球粒陨石标准化稀土元素分配模式为轻稀土元素亏损的 N-MORB 型,(La/Yh)_N=0.31～0.65(除样品 DC993为1.17)。在原始地幔标准化微量元素图上,亏损高度不相容元素,与 N-MORB 配分模式一致。相对于 Th,无 Nb、Ta的亏损,显示样品不是产于 SSZ 环境。经构造环境图解判别,样品落入了 N-MORB 区域内;这些元素成分特征表明样品具有洋中脊环境或成熟的弧后盆地环境属性。Sr、Nd 和 Pb同位素组成特征表明特提斯地幔源区以 DM(亏损地幔)为主,同时存在少量 EMⅡ(富集地幔类型Ⅱ)、Sr,Nd 和 Pb 同位素组成特征还表明特提斯地幔域具有印度洋 MORB 型的 Sr-Nd-Pb 同位素组成特征。本文的结果进一步支持 Zhang et al.(2005)的研究结果,现今印度洋不仅在地理位置上占据了曾经是特提斯洋的大部分,而且它的地幔域还继承了曾经特提斯的地幔域的地球化学特征。     

内蒙古温都尔庙和巴彦敖包-交其尔蛇绿岩的元素与同位素地球化学：对古亚洲洋东部地幔域特征的限制

内蒙古中部发育的三条蛇绿岩带是华北板块和西伯利亚板块之间的缝合带。本文系统研究了其中的温都尔庙和巴彦敖包-交其尔两个蛇绿岩带中变质玄武岩的元素和 Sr、Nd、Pb 同位素地球化学。苏右旗温都尔庙碱性玄武岩为轻稀土富集型;岩石具有板内和大陆裂谷区玄武岩的特征,可能代表了600Ma 左右,温都尔庙地区开始发育的新洋盆。采自苏左旗的巴彦敖包-交其尔玄武岩分为两类,一类呈现轻稀土富集型,呈洋岛玄武岩特征;另一类具有明显的 Nb、Ta 负异常,显示大洋岛弧玄武岩特征,洋岛玄武岩的存在表明古亚洲洋曾经发育洋盆,大洋岛弧玄武岩的存在表明古亚洲洋内部有大洋岩石圈之间的俯冲。将本文的古亚洲洋洋岛玄武岩与中国西南地区的特提斯洋岛玄武岩进行系统的元素和同位素地球化学特征对比表明,古亚洲洋的洋岛玄武岩显示高 U/Pb(HU)和北大西洋和太平洋省的特征,而特提斯洋岛玄武岩属于印度洋省。这些说明古亚洲洋地幔域与特提斯地幔域是两个独立的构造域,它们代表了长期演化的两个不同的地幔地球化学域。     

藏北高原地壳及上地幔导电性结构——超宽频带大地电磁测深研究结果

为了研究西藏中、北部壳、幔导电性结构，讨论高原中、北部岩石圈热状态，1998年和1999年（INDEPTH（Ⅲ） MT）在西藏中、北部完成了德庆—龙尾错（500线）和那曲—格尔木（600线）超宽频带大地电磁深探测剖面的研究.研究结果表明，西藏中、北部以昆仑山断裂为界，其南北壳、幔电性结构有很大差异.昆仑山断裂以北地壳和上地幔为高阻区.而昆仑山以南，地壳和上地幔的导电性有明显的分层结构：地壳上部以不连续的高阻体为主，夹有局部低阻异常体，沿南北方向上地壳的电性结构复杂，具有不连续、分块的特点；但中、下地壳为大范围的高导异常区，区内发育有大规模、不相连续、产状各异的高导体，其电阻率均小于4Ωm；在班公—怒江和金沙江缝合带之下，壳内高导体都具有向上地幔延伸的趋势，存在连通壳、幔的低阻通道.根据西藏高原中、北部壳、幔电性结构的研究推断：如同藏南一样，这里也普遍存在部分熔融体和热流体，它们的成因主要与班公—怒江和金沙江缝合带的壳-幔热交换、热活动有关，这是两期形成的壳-幔热交换通道.其中，班公—怒江缝合带的壳-幔热交换通道形成时间比金沙江缝合带早.因此，研究区壳、幔的热活动是从南边和西边开始，向北、向东扩展，导致现今西藏中、北部地壳和上地幔的热流分布由西向东、由南向北增大.     

Contaminant migration in a regular two- or three-dimensional fractured network: Reactive contaminants

A semi-analytic solution for the migration of dilute reactive contaminant species in a fractured medium having a one-, two- or three-dimensional fracture network is presented. The application of the technique to practical problems is discussed and illustrated by examples. The effects of sorption in the matrix and on the surface of the fractures are considered for a range of parameters. It is demonstrated that the synergistic interaction of the various parameters is such that the effect of changing a single variable may not be immediately obvious and a designer should perform computations for the range of expected parameters to determine the maximum impact. For example, it is shown that it is not necessarily conservative to consider only one-dimensional or only three-dimensional diffusion into the matrix; which of these cases is most critical will depend on the values of other parameters. It is also shown that an effective leachate collection system may significantly reduce the impact of those contaminants which do escape into the fractured media.