虚拟样机技术在地质工程机械领域的应用研究

基于对国内地质工程机械的研发现状的考虑，根据虚拟样机技术的特点，讨论了在地质工程机械领域实施虚拟样机技术的必要件，并提出了建立地质工程机械虚拟样机的策略。联合PRO／E、ANSYS及ADAMS建立了气动潜孔锤的虚拟样机模型，对气动潜孔锤虚拟样机进行了仿真，阐明了虚拟样机技术在地质工程机械领域的应用前景。     

A multi-method approach to study the stability of natural slopes and landslide susceptibility mapping

In this paper, a multi-method approach for the assessment of the stability of natural slopes and landslide hazard mapping applied to the Dakar coastal region is presented. This approach is based on the effective combination of geotechnical field and laboratory works, of GIS, and of mechanical (deterministic and numerical) stability analysis. By using this approach, valuable results were gained regarding instability factors, landslide kinematics, simulation of slope failure and coastal erosion. This led to a thorough assessment and strong reduction in the subjectivity of the slope stability and hazard assessment and to the development of an objective landslide danger map of the SW coast of Dakar. Analysis of the results shows that the slides were influenced by the geotechnical properties of the soil, the weathering, the hydrogeological situation, and the erosion by waves. The landslide susceptibility assessment based on this methodological approach has allowed for an appropriate and adequate consideration of the multiple factors affecting the stability and the optimization of planning and investment for land development in the city.     

Atomistic calculations of structural and elastic properties of serpentine minerals: the case of lizardite

The physical properties of the hydrous phyllosilicate lizardite have been investigated by atomistic simulation using the GULP code based on transferable semi-empirical interatomic potentials. Lizardite behavior was first investigated during structure relaxation at room temperature. The Helmholtz free energy is minimum for an equilibrium structure that is in agreement with experiment. The bulk, shear, and Young modulii for lizardite were calculated along with the Poisson ratio. From the shear and bulk modulii, we also calculated translational and longitudinal acoustic wave velocities that are important quantities for tectonophysics models. As expected, lizardite is stiffer in the a direction parallel to the layers than in the c perpendicular direction; the variation of the unit cell parameters with pressure is in good agreement with experiment. The cohesive energy between two successive layers along c direction was calculated at 0.33 eV (i.e., 0.11 eV per OH bond) in good agreement with recent ab initio calculations. Upon pressure and temperature variations, we evidenced that structural changes are mainly pressure induced; pressure being accommodated by a decrease of the c parameter up to 10 GPa. We also found that the change of slope in the derivative of the c cell parameter with respect to pressure occurring around 2 GPa originates from the bending of the interlayer hydroxyl groups with respect to the layer normal direction.     

Occurrence of salt water above fresh water in dynamic equilibrium in a coastal groundwater flow system near De Panne, Belgium

A salt water lens is found above fresh water under the shore between Dunkerque (France) and Nieuwpoort (Belgium). This inverse density distribution is in a dynamic equilibrium. It develops due to the infiltration of salt water on the back shore during high tide. Under this salt water lens, water infiltrated in the adjacent dune area flows towards the sea and discharges at the seabed. This water quality distribution differs from the classic salt water wedge under fresh water described in the literature. Here, the evolution to this water quality distribution is simulated with a density dependent numerical model. A large tidal range, shore morphology and a permeable groundwater reservoir are the main conditions for the observed water quality distribution. By altering these conditions, intermediate water quality distributions between the classic salt water wedge and the one discussed here develop. Based on these simulations, it is expected that similar kinds of inverse density distribution could be present in a number of coastal areas, which have tides, a gently sloping shore and a permeable substratum.     

3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.     

库车坳陷侏罗系煤成气动力学模拟研究

应用黄金管—高压釜封闭体系热模拟实验与GC、GC-IRMS分析技术,结合KINETICS专用软件,对库车坳陷侏罗系煤成气进行了动力学模拟研究。库车坳陷侏罗系煤具有高的产气性,在高演化阶段主要产甲烷气;侏罗系煤热解气甲烷碳同位素为-36‰～-25‰,乙烷碳同位素为-28‰～-16‰;甲烷、C2-C5气态烃的生成活化能分别为(47～64k)calm/ol、(55～72k)calm/ol,频率因子各为5.265×1013s-1、5.388×1018s-1。在此基础上,进一步探讨了克拉2气田天然气的成因。研究认为,克拉2气田天然气属阶段捕获的煤成气,主要聚集了5～1Ma时期的天然气,其成熟度Ro分布范围为1.3%～2.5%。     

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

含预制软弱带的岩石破裂过程的数值模拟

考虑到岩石脆性破坏过程中介质的不均匀性之特点，对岩石样品中含预制弱介质条带的岩石样品破坏过程进行二维有限元数值模拟，并对弱介质带的破坏贯通过程、新的断层的产生和有关的地震活动进行了研究．数值模型展示了岩石从变形、微观破坏到整体破坏的全过程以及微震活动的时空分布特征．应力、应变和微震活动的时空分布形象地描绘了岩石变形的局部化和时空迁移等现象，这与实际地壳中所观测到的现象是一致的．此外，模拟结果与实验观测结果也是一致或相似的．      

汾河口傍河水源地激发补给量研究

该文阐述傍河水源地地下水资源构成和补给特性。利用野外试验与室内数值模拟相结合方法，分析汾河口傍河水源地可采水资源量，数值模型刻画黄河实际摆动边界，考虑潜水延迟给水和极端干旱气候条件等多种方案；利用环境同位素技术追踪傍河水源地开采条件下的河水激发补给过程和补给强度。研究表明，在渗透性较好的常年性河流旁建立的水源地，具有抵御极端气象、水文环境以丰补歉的能力，枯水年动用的含水层储存量。丰水年恢复性强。