Application of the numerical manifold method for stress wave propagation across rock masses

In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.     

Sensitivity of advective travel time of contaminants to correlated formation porosities

An approach is presented to quantify sensitivity of advective contaminant travel time to porosities of hydrogeologic units (HGUs) along the flowpaths when the porosities of different HGUs are correlated. The approach is an extension of the previous sensitivity analysis technique for independent input porosity cases. It is applicable in situations where a calibrated groundwater flow model exists, but a full contaminant transport model is not available. Three sensitivity indices are introduced based on the decomposition of covariance between the advective contaminant travel time and individual input porosities of HGUs. When the input HGU porosities are correlated, the three sensitivity indices quantify the total, intrinsic and correlated contributions from each individual HGU porosity, which should be considered in order to determine the relative importance to the uncertainty in advective contaminant travel time of the input HGU porosities contributing either independently or correlatedly. Two simple one-dimensional flow examples are presented to illustrate the applications of the approach to scenarios when each HGU has distinct porosity and situations of spatially variable porosity field. The approach is applicable to more complex multi-dimensional cases where advective contaminant travel time can be calculated based on simulated flow results from groundwater flow models.     

Iron isotopic systematics of UHP eclogites respond to oxidizing fluid during exhumation

The iron stable isotope compositions (δ⁵⁶Fe) and iron valence states of ultrahigh‐pressure eclogites from Bixiling in the Dabie orogen belt, China, were measured to trace the changes of geochemical conditions during vertical transportation of earth materials, for example, oxygen fugacity. The bulk Fe³⁺/ΣFe ratios of retrograde eclogites, determined by Mössbauer spectroscopy, are consistently higher than those of fresh eclogites, suggesting oxidation during retrograde metamorphism and fluid infiltration. The studied eclogites (five samples) display limited mid‐ocean ridge basalts (MORB)‐like (~0.10‰) δ⁵⁶Fe values, which are indistinguishable from their protoliths, that is, gabbro cumulates formed through differentiation of mantle‐derived basaltic magma. This suggests that Fe isotope fractionation during continental subduction is limited. Garnet separates display limited δ⁵⁶Fe variation ranging from ?0.08 ± 0.07 ‰ to 0.02 ± 0.07‰, whereas coexisting omphacite displays a large variation of δ⁵⁶Fe values from 0.15 ± 0.07‰ to 0.47 ± 0.07‰. Omphacite also has highly variable Fe³⁺/ΣFe ratios from 0.367 ± 0.025 to 0.598 ± 0.024, indicating modification after peak metamorphism. Omphacite from retrograde eclogites has elevated Fe³⁺/ΣFe ratios (0.54–0.60) compared to that from fresh eclogites (~0.37), whereas garnet displays a narrow range of ferric iron content with Fe³⁺/ΣFe ratios from 0.039 ± 0.013 to 0.065 ± 0.022. The homogenous δ⁵⁶Fe values and Fe³⁺/ΣFe ratios of garnet suggest that it survived the retrograde metamorphism and preserved its Fe‐isotopic features and ferric contents of peak metamorphism. Because of similar diffusion rates of Fe and Mg in garnet and omphacite, and constant Δ²⁶Mg_{omphacite‐garnet} values (1.14 ± 0.04‰), equilibrium iron isotope fractionation between garnet and omphacite was probably achieved during peak metamorphism. Elevated Fe³⁺/ΣFe ratios of omphacite from retrograde eclogites and variant Δ⁵⁶Fe_{omphacite‐garnet} values of the studied eclogites (0.13 ± 0.10‰ to 0.48 ± 0.10‰) indicate that oxidized geofluid infiltration resulted in the elevation of δ⁵⁶Fe values of omphacite during retrograde metamorphism.     

依兰-伊通断裂新构造活动规律分析

作为郯庐断裂带北段主干的依兰-伊通断裂, 其新构造活动性与活动规律仍然存在不同的认识.本次工作通过详细的野外调查, 发现该断裂内活断层广泛存在, 由东、西两支北东走向的主干活断层构成, 沿着古近纪地堑边界断层发育.这些活断层主要呈破碎型结构, 多为逆右行平移活动.通过对这些活断层一系列实测擦痕反演应力场, 显示它们多是在东西向挤压中活动的, 而现今应力场转变为北东东-南西西向区域性挤压.依据本次野外观察与¹⁴ C定年, 并结合前人定年结果与近代地震分布, 表明依兰-伊通西支活断层的最新活动时代为全新世与晚更新世相间, 而东支活断层的最新活动时代主要为早-中更新世.依兰-伊通断裂内活断层显示了明显的差异性活动, 表现为西支的活动强度明显大于东支, 西支的最新活动时代皆晚于东支, 沿走向上活动性强、弱相间与最新活动时代不断变化, 以及近代地震活动不均一分布.它们沿走向上的分段性、差异性活动主要是因为被一系列北西向断层切断所致.     

三塘湖盆地石炭系卡拉岗组火山岩岩相研究

三塘湖盆地石炭系卡拉岗组火山岩岩性揭示该区火山岩岩相主要包括火山爆发相、溢流相,火山沉积相不发育。为了明确不同岩相分布规律,本文系统利用测井、地质、地震等多种资料,建立测井岩性识别模板,以单井相作为约束,开展地震反演,明确火山岩岩相的平面分布规律。研究结果表明:该区火山喷发方式以裂隙式为主,火山口呈串珠状沿主断裂排列;溢流相平面呈层状分布,溢流相安山岩最为发育。研究结果可有效指导该区火山岩油藏的勘探评价工作,同时为火山岩岩相研究提供可借鉴的研究思路。     

Effects of low- to high-angle normal faults on sedimentary architectures in the Eocene Wenchang Formation,Enping Sag,Pearl River Mouth Basin,South China Sea

The effects of low- to high-angle (>30°) normal faults on sedimentary architectural units in the Eocene Wenchang Formation, Enping Sag, Pearl River Mouth Basin (PRMB), South China Sea were investigated utilising a high-quality 3D seismic data set and restored paleogeomorphology. It has been shown that sequence stratigraphic units and sedimentary architecture are significantly controlled by the low- to high-angle normal faults. The Wenchang Formation, a second-order sequence, can be subdivided into two para-second-sequences (the Lower and Upper Wenchang sequences, E₂W^L and E₂W^U) and seven third-order sequences (from base to top: SQ1～SQ7). The low-angle fault confined sequence architecture of the Wenchang Formation is mainly characterised by lateral stacking with the ratio of the vertical subsidence (V) to horizontal slip (H) being reduced from 1/2 for E₂W^L to 1/6 for E₂W^U. In contrast, the high-angle fault confined sequence is characterised by vertical stacking with the ratio of V/H close to 1 for sequences SQ1 to SQ7. In the 3D seismic area, the features of sediment-dispersal pattern were interpreted based on an integrated analysis of paleogeomorphology, seismic reflection characteristics, stratal thickness distribution and multiple attribute clustering. The results show that the large-scale fan delta, belt-shape lacustrine deposit and bird-foot braided delta systems mainly developed in the low-angle fault confined sequences, whereas small-scale fan delta, rhombus-shaped lacustrine deposit and lobe-shaped braided delta systems inherited tectono-sedimentary architectures in the high-angle fault confined sequences.     

FDR土壤水分传感器原状土壤标定方法研究

为获得测站准确的土壤水分传感器标定参数,利用自动土壤水分站点的大型原状土壤为样本,在实验室同时进行人工称重观测和仪器自动对比观测,得到一系列的人工自动测量数据,以人工称重数据为准对自动土壤水分传感器进行订正。由于采用原状土壤,样本土壤的质地、密度、土壤颗粒和结合紧密度基本没有受到破坏,订正后的自动土壤体积含水量值能更加真实反映台站实际土壤墒情,能够修订人工对比观测期间由于样本空间不足导致的不合理的田间标定方程。