中国与美国西部平均构造剪应力值的比较

提出了一种估计一个地区平均构造剪应力值的方法．用这种方法，选取哈佛大学公布的1977～1999年共15 993次地震的矩张量数据，估计了美国西部19个地区、中国和邻区43个地区(每个地区为1010范围)的平均环境剪应力值．结果表明，美国西部南加州西部海域和南加州应力最高，达13.7和12.0 MPa，然后向北、向南和向东逐渐递减，但最小也达8.7 MPa，是最高值的63%．中国新疆西北部地区和西藏察隅地区应力水平最高，达17.2和12.9 MPa，比美国高．中国的华北、云南、四川、台湾和美国南加州的应力水平差不多．中国南北地震带的应力水平为13 MPa左右，比南加州略高．两个重要地区的平均剪应力值分布图，提供了地学的基本数据．这些结果可为研究地震活动的大背景提供依据，对研究强地面运动参数(如峰值加速度及反应谱等)的衰减关系也是有用的．      

用应变积累释放模型研究中国大陆地块分区地震活动

依据活动地块假说及活动边界的划分，对中国大陆地震活动进行分区．在此基础上，应用应变积累释放模型，通过地震应变积累释放的计算，研究了中国大陆各活动地块分区的地震活动性，初步探讨了各地块的地震活动阶段划分．此外，还对应变积累释放模型的应用前提条件进行了必要的讨论，并对模型结果中可能出现的问题给予一定分析．      

流动成像测井研究进展

流动成像测井属于对非均匀介质动态的非线性测量，通过对油气井内多相流体信息的投影扫描和反演处理，实时成像显示流体分布及流动状况，国外研制的电导法，电容法仪器采用电导探针或电容元件构成阵列测量探头，分别利用油气与水的导电特性和介电特性差异辨识井内流体，对流动截面的测量局限于个别点上，已在水平井中见到应用效果，但对物场信息投影测量的数据量和分辨率未能满足成像要求。我国提出的电磁法采用环状阵列电极构成测量探头，综合利用油气与水的导电和介电特性差异辨识井内流体，可以测量获得多相流体流动截面的清晰图像，为研究开发全新的流动成像测井系统奠定了良好基础。     

Seismic analysis of base-isolated liquid storage tanks using the BE–FE–BE coupling technique

A three-dimensional soil–structure–liquid interaction problem is numerically simulated in order to analyze the dynamic behavior of a base-isolated liquid storage tank subjected to seismic ground motion. A dynamic analysis of a liquid storage tank is carried out using a hybrid formulation, which combines the finite shell elements for structures and the boundary elements for liquid and soil. The system is composed of three parts: the liquid–structure interaction part, the soil–foundation interaction part, and the base-isolation part. In the liquid–structure interaction part, the tank structure is modeled using the finite elements and the liquid is modeled using the internal boundary elements, which satisfy the free surface boundary condition. In the soil–foundation interaction part, the foundation is modeled using the finite elements and the half-space soil media are modeled using the external boundary elements, which satisfy the radiation condition in the infinite domain. Finally, above two parts are connected with the base-isolation system to solve the system's behavior. Numerical examples are presented to demonstrate the accuracy of the developed method, and an earthquake response analysis is carried out to demonstrate the applicability of the developed technique. The properties of a real LNG tank located in the west coast of Korea are used. The effects of the ground and the base-isolation system on the behavior of the tank are analyzed.     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

Screening model for volatile pollutants in dual porosity soils

This paper develops mass fraction models for transport and fate of agricultural pollutants in structured two-region soils. Mass fraction index models, based on a semi-infinite domain solution, are derived that describe leaching at depth, vapor losses through soil surface, absorption, and degradation in the dynamic- and stagnant-water soil regions. The models predict that leaching is the result of the combined effect of the upward vapor-phase transport relative to downward advection, residence time relative to half-life, dispersion, and lateral diffusive mass transfer. Simulations show that leached fraction of volatile compounds does not always decrease monotonically with increased residence time relative to the pollutant half-life, as a result of complex interactions among the different physical and biochemical processes. The results show that leaching, volatilization, and degradation losses can be affected significantly by lateral diffusive mass transfer into immobile-water regions and advection relative to dispersion (i.e. Peclet number) in the mobile-water regions. It is shown that solute diffusion into the immobile phase and subsequent biochemical decay reduces leaching and vapor losses through soil surface. Potential use of the modified leaching index for the screening of selected pesticides is illustrated for different soil textures and infiltration rates. The analysis may be useful to the management of pesticides and the design of landfills.     

Predicting average annual groundwater levels from climatic variables: an empirical model

On the basis of one-dimensional theoretical water flow model, we demonstrate that the groundwater level variation follows a pattern similar to recharge fluctuation, with a time delay that depends on the characteristics of aquifer, recharge pattern as well as the distance between the recharge and observation locations. On the basis of a water budget model and the groundwater flow model, we propose an empirical model that links climatic variables to groundwater level. The empirical model is tested using a partial data set from historical records of water levels from more than 80 wells in a monitoring network for the carbonate rock aquifer, southern Manitoba, Canada. The testing results show that the predicted groundwater levels are very close to the observed ones in most cases. The overall average correlation coefficient between the predicted and observed water levels is 0.92. This proposed empirical statistical model could be used to predict variations in groundwater level in response to different climate scenarios in a climate change impact assessment.     

Interpretation of charging on fracture or frictional slip surface of rocks

Fracturing and frictional sliding of quartz and granite under dry condition generates fractoluminescence, charged particle emission and electromagnetic radiation. Various kinds of experiments indicate that surface charge density on fracture or frictional slip surface of quartz and granite is 10⁻⁴ to 10⁻² C/m² which is larger than bound charges induced by the disappearance of piezoelectricity due to the release of stress. Hole and electron trapping centers, which is found in semiconductor devices with the Si–SiO₂ system, are causes of surface charging on fracture or frictional slip surface of quartz crystal. The quantity of the surface charge is enough to cause corona discharge that can generate earthquake lights. The mechanism considering the hole and electron trapping centers has a probability to explain why non-piezoelectric minerals or rocks generate electromagnetic phenomena. It can be one of origins of seismo-electromagnetic phenomena (SEP).     

结构动力反应分析的三阶显式方法

本通过对传统动力反应分析方法的总结，阐明了建立隐式和显式方法的一般思路及数学本质，提出了使用系统位移反向向量三阶导数的隐工和实用显式积分方法－3阶显式方法，分析了该显式方法的精度和稳定性，并对建立更高阶隐式和显式方法以及方法的精度和稳定性作了初步讨论。最后，通过算例对本方法、献[1]方法和经典的常平均加速度法（隐式方法、视为精确解）的精度和稳定性进行了比较分析。结果表明，本方法具有明显的优点。     

液化场地桩-土-桥梁结构动力相互作用振动台试验研究进展

本文在全面归纳与总结液化场地桩－土－桥梁结构动力相互作用振动台试验及与之相关领域的国内外研究进展基础上，直接针对我国桥梁工程中的主要震害问题，提出在我国开展液化场地桩－土－桥梁结构动力相互作用振动台试验研究的必要性，并阐述作者对液化场地桩－土－张桥梁结构动力相互作用振动台试验中若干问题的认识。