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.     

Chaotic Diffusion And Effective Stability of Jupiter Trojans

It has recently been shown that Jupiter Trojans may exhibit chaotic behavior, a fact that has put in question their presumed long term stability. Previous numerical results suggest a slow dispersion of the Trojan swarms, but the extent of the ‘effective’ stability region in orbital elements space is still an open problem. In this paper, we tackle this problem by means of extensive numerical integrations. First, a set of 3,200 fictitious objects and 667 numbered Trojans is integrated for 4 Myrs and their Lyapunov time, T_L, is estimated. The ones following chaotic orbits are then integrated for 1 Gyr, or until they escape from the Trojan region. The results of these experiments are presented in the form of maps of T_Land the escape time, T_E, in the space of proper elements. An effective stability region for 1 Gyr is defined on these maps, in which chaotic orbits also exist. The distribution of the numbered Trojans follows closely the T_E=1 Gyr level curve, with 86% of the bodies lying inside and 14% outside the stability region. This result is confirmed by a 4.5 Gyr integration of the 246 chaotic numbered Trojans, which showed that 17% of the numbered Trojans are unstable over the age of the solar system. We show that the size distributions of the stable and unstable populations are nearly identical. Thus, the existence of unstable bodies should not be the result of a size-dependent transport mechanism but, rather, the result of chaotic diffusion. Finally, in the large chaotic region that surrounds the stability zone, a statistical correlation between T_LandT_E is found.     

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.     

地球化学景观吸引子的发现及其意义

勘查地球化学对异常评价的最大不确定性源于对地球化学景观特性的认识不够。我们应用构造相空间的方法分析了浙江儒岙幅（1:5万）水系沉积物测量的Ag含量数据，揭示了地球化学景观是具有低维吸引子（DZ≈2．9）的混饨系统。这一发现将对勘查地球化学中采样密度的选择、地球化学场模拟、异常评价等根本性问题产生重大影响。     

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.     

Room-temperature equation of state of Li<Subscript>2</Subscript>VOSiO<Subscript>4</Subscript> up to 8.5 GPa

High-pressure single-crystal X-ray diffraction measurements of lattice parameters of the compound Li₂VOSiO₄, which crystallises with a natisite-type structure, has been carried out to a pressure of 8.54(5) GPa at room temperature. Unit-cell volume data were fitted with a second-order Birch-Murnaghan EoS (BM-EoS), simultaneously refining V ₀ and K ₀ using the data weighted by the uncertainties in V. The bulk modulus is K ₀ = 99(1) GPa, with K′ fixed to 4. Refinements of third order equations-of-state yielded values of K′ that did not differ significantly from 4. The compressibility of the unit-cell is strongly anisotropic with the c axis (K ₀(c) = 49.7 ± 0.5 GPa) approximately four times more compressible than the a axis (K ₀(a) = 195 ± 3 GPa).     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.