Finite difference modeling of contaminant transport using analytic element flow solutions

The two-dimensional implementation of the analytic element method (AEM) is commonly used to simulate steady-state saturated groundwater flow phenomena at regional and local scales. However, unlike alternative groundwater flow simulation methods, AEM results are not ordinarily used as the basis for simulation of reactive solute transport. The use of AEM-simulated flow fields is impeded by the discrepancy between a continuous representation of flow and a typically discrete representation of transport, and requires translation of the flow solution to a discrete analog. This paper presents a variety of methods for analytically calculating conservative discrete water fluxes and integrated components of the dispersion tensor across cell interfaces. An Eulerian finite difference method based on these AEM-derived parameters is implemented for use in simulation of 2D (vertically averaged) solute transport. This implementation is first benchmarked against existing methods that use standard finite difference flow solutions, then used to investigate the effects of an inaccurate discrete water balance. It is shown that improper translation of AEM fluxes leads to significant water balance errors and inaccurate simulation of contaminant transport.     

System identification of a concrete arch dam and calibration of its finite element model

A magnitude 4.3 earthquake occurred near Pacoima Dam on 13 January 2001. An accelerometer array that had been upgraded after the Northridge earthquake recorded the motion with 17 channels on the dam and the dam–foundation interface. Using this data, properties of the first two modes are found from a system identification study. Modal properties are also determined from a forced vibration experiment performed in 2002 and indicate a significantly stiffer system than is estimated from the 2001 earthquake records. The 2001 earthquake, although small, must have induced temporary nonlinearity. This has implications for structural health monitoring. The source of the nonlinear behaviour is believed to be loss of stiffness in the foundation rock. A finite element model of Pacoima Dam is constructed and calibrated to match modal properties determined from the system identification study. A dynamic simulation of the 2001 earthquake response produces computed motions that agree fairly well with the recorded ones. Copyright © 2006 John Wiley & Sons, Ltd.     

Transient analysis of wave propagation in non‐homogeneous elastic unbounded domains by using the scaled boundary finite‐element method

The scaled boundary finite‐element method has been developed for the dynamic analysis of unbounded domains. In this method only the boundary is discretized resulting in a reduction of the spatial dimension by one. Like the finite‐element method no fundamental solution is required. This paper extends the scaled boundary finite‐element method to simulate the transient response of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. To reduce the number of degrees of freedom and the computational cost, the technique of reduced set of base functions is applied. The scaled boundary finite‐element equation for an unbounded domain is reformulated in generalized coordinates. The resulting acceleration unit‐impulse response matrix is obtained and assembled with the equation of motion of standard finite elements. Numerical examples of non‐homogeneous isotropic and transversely isotropic unbounded domains demonstrate the accuracy of the scaled boundary finite‐element method. Copyright © 2006 John Wiley & Sons, Ltd.     

河南省固始县柳林-长湾化探Pb、Zn异常评价方法及效果

对河南固始县1989年1/20万区域化探扫面资料进行二次开发,将异常区中的Pb、Zn(A g)地球化学异常重新筛选,评序后提出以排序第二的13-乙1号Pb、Zn(A g)异常为重点评价对象。评价结果显示,在南、北、中三个异常带中,均找到了新矿体,取得了一定的地质找矿效果。     

Equatorial and low latitude ionosphere during intense geomagnetic storms

An investigation is made in order to analyse the role of neutral gas composition in the equatorial and low latitude ionosphere during intense geomagnetic storms. To this end data taken by the Dynamic Explorer 2 satellite at 280–300 km (molecular nitrogen N₂ and atomic oxygen O concentrations, electron density and vertical plasma drifts) are used. The sudden commencements of the events considered occurred at 11:38 UT on March 1, 1982, 18:41 UT on November 20, 1982 and 16:14 UT on February 4, 1983. Vertical plasma drifts are the most important contributor to the initial storm time response of the equatorial F region. Neutral composition changes (expressed as an increase in the molecular species, mainly N₂) possibly play a predominant role in the equatorial and low latitude (10–20°) decreases of electron density at heights near F2-region maximum during the main and recovery phases of intense geomagnetic storms. Delayed increases of electron density observed at daytime during the recovery phase may be also attributed to increases in atomic oxygen. At low latitudes possibly a combined effect of O increase and upward plasma drift due to enhanced equatorward winds is the responsible mechanism for the maintenance of enhanced electron density values.     

Static and dynamic response of massive caisson foundations with soil and interface nonlinearities—validation and results

The static, cyclic, and dynamic response of a massive caisson foundation embedded in nonlinear layered or inhomogeneous soil and loaded at its top is investigated. The caisson is supported against horizontal displacement and rotation by four types of inelastic springs and dashpots, described with the BWGG model that was developed in the preceding companion paper [Gerolymos N, Gazetas G. Development of winkler model for static and dynamic response of caisson foundations with soil and interface nonlinearities. Soil Dyn Earthq Eng, submitted companion paper]. The prediction of the model is satisfactorily compared with results from 3D-finite element analysis. Some experimental corroboration of the method is provided with the help of a 1/3-scale lateral load test that had been conducted in the field by EPRI. An illustrative example of a caisson embedded in linearly-inhomogeneous clay and subjected to static and dynamic loading is analysed. Characteristic results are presented highlighting the role of soil inelasticity and its interplay with the two dominant interface nonlinearities: separation (gapping) of the caisson shaft from the surrounding soil, and uplifting of the base from the underlying soil.     

可地浸砂岩型铀矿床航放弱铀异常研究及其找矿意义

认真研究新疆伊犁盆地南缘511、512两个已知可地浸砂岩型铀矿床的航放数据,发现矿体上方均存在微弱的航放铀偏高显示,由此提出了可地浸砂岩型铀矿床“弱铀异常”的新概念。基于地气运移机制和氡的接力迁移理论,探讨了“弱铀异常”的形成机理,介绍了航放弱铀异常的提取方法,指出高精度航放方法完全能探测到“弱铀异常”。弱铀异常可以作为一种直接标志信息运用于前期成矿预测。     

巴姆地震变形场和应力场：Ⅱ．用FEPG有限元方法求解

本文利用 FEPG 有限元程序对2003年12月的巴姆地震变形场和应力场进行了模拟计算。2003年12月26日,巴姆地区发生6.6级地震。Nakamura、Suzuki 等人利用余震资料和地震破坏程度的分布图分析了余震震中的统计分布形式,得到了此次地震的发震断层是巴姆断层以西大约5km 的 Arg-e-Bam 隐伏断层的结论。凌勇等人从差分干涉雷达(D-InSAR)出发,得到了此次地震的同震干涉条纹和雷达视线方向的位移场,并从雷达干涉的相干图上确定了发震断层的痕迹。他们还利用 Okada 程序模拟计算了此次地震的位移场,模拟计算的结果与 D-InSAR 测量的位移场比较符合。Okada 程序从断层的位错为出发点进行计算,从运动学角度对此次地震变形场和应力场进行研究;本文从动力学的角度出发,研究远场区域边界的应力状态以及断层的物性参数的变化对巴姆地震形变场的影响,将各种计算结果进行对比的同时也和前人的结果相比较。本文计算中采用了凌勇等人的断层参数,即 Arg-e-Bam 隐伏断层为近南北走向,向东倾斜,近直立的倾角,右旋走滑,其(Strike,dip,slip)为(175°,80°,173°)。断层未延伸到地表,自地面以下1km 处往下延伸16km 的宽度,断层长度为20km。对断层的杨氏模量、泊松比、断层附近区域的边界受力大小、两个主应力的比例以及主应力的方向的改变对计算结果的影响做了讨论。最后通过上述几组计算结果建立了一个与前人结果较为符合的巴姆地震参数模型。通过计算,一方面揭示了应力状态和物性参数的影响,另一方面也说明 Arg-e-Bam 隐伏断层是此次地震的发震断层这一结论是可信的。     

东沟特大型斑岩钼矿床的物化探找矿效果

对河南东沟钼矿床的区域物化探成果进行整理、研究,运用物化探资料循序渐进地揭示出钼矿床的空间分布规律。东沟特大型钼矿床的发现,不仅为我国的钼矿产业今后几十年的稳步发展将产生深远的影响,而且为在东秦岭—大别山钼成矿带进一步运用物化探作为充分必要手段,为寻找更多的钼矿资源提供了难得的成功经验。     

土壤污染地球化学标准及等级划分问题讨论

全国多目标区域地球化学调查发现有害元素广泛分布引起极大关注。土壤污染将长期危及生态系统和人类生存安全。制定土壤污染等级标准十分重要。国家土壤环境保护标准体系中,土壤污染等级标准是针对调查发现的污染进行量度划分,土壤环境质量标准是依据污染对环境影响程度确定临界值,土壤污染评价标准是针对存在问题进行科学研究的指标。土壤污染等级采用统一的分级标准,划分为国家级、流域级、省级和地市级。标准制定是一个由局部到全局的认识,又由全局到局部的决策过程。在多目标区域地球化学调查基础上首先制定土壤污染标准及划分污染等级是笔者的主导思想。