Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

青藏块体东北缘地壳水平运动状态

应用青藏块体东北缘1999～2003年多期GPS观测资料，计算了不同时段GPS点水平运动速率。通过分析发现：甘青块体可分为东部块体和西部块体，东、西部块体的运动状态存在明显的差异；受2001年11月14日昆仑Ms8．1地震的影响，震后地壳运动状态发生了明显的改变。     

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Reliability assessment of serviceability performance of braced retaining walls using a neural network approach

In urban environments, one major concern with deep excavations in soft clay is the potentially large ground deformations in and around the excavation. Excessive movements can damage adjacent buildings and utilities. There are many uncertainties associated with the calculation of the ultimate or serviceability performance of a braced excavation system. These include the variabilities of the loadings, geotechnical soil properties, and engineering and geometrical properties of the wall. A risk‐based approach to serviceability performance failure is necessary to incorporate systematically the uncertainties associated with the various design parameters. This paper demonstrates the use of an integrated neural network–reliability method to assess the risk of serviceability failure through the calculation of the reliability index. By first performing a series of parametric studies using the finite element method and then approximating the non‐linear limit state surface (the boundary separating the safe and ‘failure’ domains) through a neural network model, the reliability index can be determined with the aid of a spreadsheet. Two illustrative examples are presented to show how the serviceability performance for braced excavation problems can be assessed using the reliability index. Copyright © 2005 John Wiley & Sons, Ltd.     

Relationship between footwall composition,crustal contamination,and fluid–rock interaction in the Platreef,Bushveld Complex,South Africa

In the northern limb of the 2.06-Ga Bushveld Complex, the Platreef is a platinum group elements (PGE)-, Cu-, and Ni-mineralized zone of pyroxenite that developed at the intrusion margin. From north to south, the footwall rocks of the Platreef change from Archaean granite to dolomite, hornfels, and quartzite. Where the footwall is granite, the Sr-isotope system is more strongly perturbed than where the footwall is Sr-poor dolomite, in which samples show an approximate isochron relationship. The Nd-isotope system for samples of pyroxenite and hanging wall norite shows an approximate isochron relationship with an implied age of 2.17 ± 0.2 Ga and initial Nd-isotope ratio of 0.5095. Assuming an age of 2.06 Ga, the ɛNd values range from −6.2 to −9.6 (ave. −7.8, n = 17) and on average are slightly more negative than the Main Zone of the Bushveld. These data are consistent with local contamination of an already contaminated magma of Main Zone composition. The similarity in isotope composition between the Platreef pyroxenites and the hanging wall norites suggests a common origin. Where the country rock is dolomite, the Platreef has generally higher plagioclase and pyroxene δ ¹⁸O values, and this indicates assimilation of the immediate footwall. Throughout the Platreef, there is considerable petrographic evidence for sub-solidus interaction with fluids, and the Δ _{plagioclase–pyroxene} values range from −2 to +6, which indicates interaction at both high and low temperatures. Whole-rock and mineral δD values suggest that the Platreef interacted with both magmatic and meteoric water, and the lack of disturbance to the Sr-isotope system suggests that fluid–rock interaction took place soon after emplacement. Where the footwall is granite, less negative δD values suggest a greater involvement of meteoric water. Consistently higher values of Δ _{plagioclase–pyroxene} in the Platreef pyroxenites and hanging wall norites in contact with dolomite suggest prolonged interaction with CO₂-rich fluid derived from decarbonation of the footwall rocks. The overprint of post crystallization fluid–rock interaction is the probable cause of the previously documented lack of correlation between PGE and sulfide content on the small scale. The Platreef in contact with dolomite is the focus of the highest PGE grades, and this suggests that dolomite contamination played a role in PGE concentration and deposition, but the exact link remains obscure. It is a possibility that the CO₂ produced by decarbonation of assimilated dolomite enhanced the process of PGE scavenging by sulfide precipitation.     

Possible chemical structures and biological precursors of different vitrinites in coal measure in Northwest China

Chemical structure of Jurassic vitrinites isolated from the coals in basins in NW China have been checked using solid state ¹³C NMR and flash pyrolysis-GC/MS. Study shows some Jurassic collodetrinites are rich in aliphatic products in pyrolysates, consisting with the high amount of methylene carbon in ¹³C NMR spectra. In contrast, pyrolysates of Jurassic collotelinites are rich in phenols and alkylbenzenes. Also one Pennsylvanian and one Permian vitrinite selected from the Ordos basin, NW China have been checked for comparison. The proportion of aliphatics is low in pyrolysates, and aliphatic carbon peak in ¹³C NMR spectrum of Permian vitrinite is mostly composed of gas-prone carbons compared with collodetrinites in those Jurassic basins. But both pyrolysis and ¹³C NMR data shows the Pennsylvanian vitrinite is not only gas-prone but also oil-prone. Relatively high proportion of long chain aliphatic structure of some Jurassic vitrinite in Junggar, Turpan-Hami basins may be due to the contribution of liptodetrinites, which may be included during the formation of vitrinites. And it seems that suberinite is the most possible precursor of long chain aliphatics in the structure of Jurassic collodetrinite.     

A preliminary stable isotope study on Mogok Ruby, Myanmar

The primary occurrence of ruby in the Mogok area, northern Myanmar is exclusively found in marble along with spinel–forsterite-bearing marble and phlogopite–graphite marble. These marble units are enclosed within banded biotite–garnet–sillimanite–oligoclase gneisses. Samples of these marbles collected for C–O stable isotope analysis show two trends of δ¹³C–δ¹⁸O variation resulting most likely from fluid–rock interactions. Ruby-bearing marble and phlogopite–graphite marble follow a trend with coupled C–O depletion, whereas spinel–forsterite-bearing marble follows a δ¹⁸O depletion trend with relatively constant δ¹³C values. Ruby formation might have resulted from CO₂-rich fluid–rock interaction, while spinel–forsterite-bearing marble was genetically related to CO₂-poor fluid–rock interaction. Both fluids may have arisen from external sources. Based on graphite Raman spectral thermometry, the estimated temperature for phlogopite–graphite marble, and probably ruby-bearing marble, was lower than 607 °C, and for spinel–forsterite-bearing marble, lower than 710 °C. Contrasting C/O diffusion between graphite/ruby/spinel/forsterite and calcite, local variations of isotopic compositions of newly formed minerals as a result of non-pervasive fluid infiltration, and open-system isotopic disturbance during cooling may have affected C-/O-isotopic fractionations between minerals. The estimated high formation temperatures for ruby and spinel/forsterite imply that the parental fluids may have been related to nearby igneous intrusions and/or metamorphic processes. Whether these two types of fluid were genetically related is unclear based on the present data.     

Biostratigraphic and isotopic record of the Cenomanian–Turonian deposits in the Ohaba-Ponor section (SW Haţeg, Romania)

The middle Cenomanian–lower Turonian deposits of Ohaba-Ponor section (Southern Carpathians) were studied from biostratigraphic and isotopic points of view. Both the qualitative and semiquantitative nannofloral analyses, as well as the stable isotope (δ¹³C and δ¹⁸O) data support significant palaeoenvironmental changes in the investigated interval. Two δ¹³C positive excursions were recognized: (1) an excursion up to 1.8‰ (PDB) within the middle/late Cenomanian boundary; (2) an excursion up to 2.2‰ (PDB) in the Cenomanian/Turonian boundary interval. The oldest δ¹³C positive excursion recorded (placed within the Acanthoceras jukes-brownei/Eucalycoceras pentagonum Ammonite Zone boundary interval, and in the NC11 Calcareous Nannofossil Zone respectively) could be assigned to the middle Cenomanian Event II (MCEII). During the above-mentioned event, significant increase in abundance of Watznaueria barnesae, followed by successive blooms of Biscutum constans and Eprolithus floralis, were observed. The youngest δ¹³C positive excursion was identified in the Cenomanian/Turonian boundary interval (in the NC12 and lower part of the NC13 Calcareous Nannofossil Zones). Even the amplitude of this δ¹³C positive excursion is lower in the Ohaba-Ponor section, as generally reported, this may represent the regional record of the OAE2. The successive peaks of the nannofossils Biscutum constans, Zeugrhabdotus erectus and Eprolithus floralis indicate episodes of cooler surface water and high fertility, which preceded and lasted the Cenomanian/Turonian boundary event. Additionally, fluctuations of δ¹⁸O values between −2 and −6‰ suggest also cooler conditions within the Cenomanian/Turonian boundary interval.     

Diagenesis, palaeoclimate and tectono-sedimentary influences on clay mineralogy and stable isotopes from Upper Cretaceous marine successions of the Basque-Cantabrian Basin (N Spain)

Clay mineralogy and whole-rock stable isotopes (δ¹⁸O and δ¹³C) of Upper Cretaceous marly sediments on the Basque-Cantabrian Basin have been integrated to determine the main effects of diagenesis, palaeoclimate and tectono-sedimentary factors in sections belonging to deep- (Barrika) and platform-marine (Isla de Castro, Villamartín and Olazagutía) settings.The mean values for the clay assemblages and δ¹⁸O exhibit notable differences among the sections, partially explainable by the influence of diagenesis. The Barrika sediments, with more diagenetically advanced illite-smectite (I-S) mixed-layer (R1, 70% illite), authigenic chlorite, and low δ¹⁸O (−4.05‰ PDB), experienced higher diagenetic grade than Isla de Castro and Olazagutía, which have R0 I-S (20% illite) and heavier δ¹⁸O. Villamartín was also affected by higher diagenesis than Isla de Castro and Olazagutía, given the occurrence of R1 I-S (60% illite) and low δ¹⁸O (−4.11‰ PDB). However, the absence of other clays in Villamartín (e.g. authigenic chlorite) is indicative of less diagenetic grade than Barrika. These results show the useful integration of clay mineralogy and stable isotopes to detect different diagenetic grades in distinct marine successions of the same basin.Despite being influenced by diagenesis, the clay mineralogy partially preserves its inherited signature. This allows detection of major contents of I-S and mica, and minor kaolinite, interpreted as indicative of warm palaeoclimatic conditions. High kaolinite content in Villamartín and absence of kaolinite in Isla de Castro, though, are considered to be a product of neither diagenesis nor palaeoclimatic influences. Instead, tectono-sedimentary causes, related to unsuitable conditions for clay formation and transport from the local source areas, contributed to original clay differences. The inferred effects of diagenesis, palaeoclimate and tectono-sedimentary factors make this work important to show the potentially great variety of controls on the clay mineralogy of marine sections, which are often uncritically treated in studies concerning the Late Cretaceous.