FLAC数值模拟软件及其在地学中的应用

FLAC是以拉格朗日差分法为基础的一种数值模拟应用软件,主要用于岩土力学工程。近年地质学家将其应用到地质构造和块体运动学、动力学研究中,以证实或证伪地球学家们提出的相关地球动力学模型。本文初步介绍了FLAC应用软件的基本原理、使用方法以及它在国内外地质学上的几个应用实例。     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

湘中锑(金)矿床成矿物质来源——Ⅰ.微量元素及其实验地球化学证据

为了查明湘中盆地中-低温梯（金）矿床的成矿物质来源，文章应用高分辨率ICP-MS、ICP-AES等方法对岩石与矿石进行了系统测定，并结合作者以往分析数据和实验地球化学结果，对该区盆地及周边基底地层与矿床中的微量元素分布特征及其水/岩反应效应进行了综合分析和系统研究。结果表明：Sb、Au等成矿元素在元古界基底碎屑岩中具有高的背景含量，在近似成矿条件下的水/岩反应实验中，成矿元素的淋出率达20%～90%，元古界基底碎屑岩中的金（锑、钨）矿床的周围蚀变围岩中出现明显的金负异常区。而盆地内泥盆系统地层中的高锑含量仅局限于锑矿床（化）及其蚀变围岩中；湘中锑（金）矿床中的矿石与基底碎屑岩具有相似的稀土元素球粒陨石标准化配分曲线特征，以上微量元素地球化学特征充分证明该区锑（金）矿床锑、金等成矿物质主要来源于元古界底碎屑岩。     

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     

川滇黔桂76种元素地球化学图编制中分析方法与分析质量研究(一) 不同实验室产生地球化学图的相似性——以Ag、Cs、Ga、Ge为例

编制中国76种元素地球化学图将使我们首次了解元素周期表上除氢、氧及惰性气体之外所有元素在中国大陆表层的分布,川滇黔桂76种元素地球化学图的编制是这一研究项目的试点。本文以不同实验室取得的Ag、Cs、Ga、Ge分析数据为例,探索、衡量不同实验室地球化学图相似性的方法,以优选实验室进行川滇黔桂4省区约4000件组合样76种元素的大规模分析工作。     

Experimental Study on Mine Spoils During Superficial Geochemistry Processes

Based on the mineral and chemical compositions of spoils taken from the Shilu Cu-Mo deposit and Hetai gold deposit, the leaching and batch experiments were made on spoils taken from these two deposits using the leaching column designed by the authors. The experimental results showed that it is not always true that the mine drainage is acidic. Its acidity depends on acid-buffering capacity of gangue, host rock and its alteration mineral assemblage. The composition of the drainage water is related with the interaction between superficial or underground water and solid materials in the spoils, including minerals, hydroxides and amorphous substances. The leaching extent of the element is related with its occurrence form in the deposit. The preferential flow results in leaching-out of heavy metals in large amounts. So it is important to prevent the generating of preferential flow in the system. The results will provide very important grounds for the comprehensive management of land and ecological rehabilitation of the mine site.     

Computational modeling of cyclic mobility and post-liquefaction site response

Under seismic excitation, liquefied clean medium to dense cohesionless soils may regain a high level of shear resistance at large shear strain excursions. This pattern of response, known as a form of cyclic mobility, has been documented by a large body of laboratory sample tests and centrifuge experiments. A plasticity-based constitutive model is developed with emphasis on simulating the cyclic mobility response mechanism and associated pattern of shear strain accumulation. This constitutive model is incorporated into a two-phase (solid–fluid), fully coupled finite element code. Calibration of the constitutive model is described, based on a unique set of laboratory triaxial tests (monotonic and cyclic) and dynamic centrifuge experiments. In this experimental series, Nevada sand at a relative density of about 40% is employed. The calibration effort focused on reproducing the salient characteristics of dynamic site response as dictated by the cyclic mobility mechanism. Finally, using the calibrated model, a numerical simulation is conducted to highlight the effect of excitation frequency content on post-liquefaction ground deformations.     

地震前兆时空非均匀性指标Cv值的实验检验

在双向加载条件下，对两种构造和介质不同的岩石标本进行变形实验，测量了标本表面不同部位的应变．在此基础上，利用王晓青和陈学忠等提出的描述地震前兆非均匀分布的参量————Cv值，分析了标本变形失稳过程中应变异常的分布，以期从实验角度检验Cv值方法并探讨其物理意义．研究表明，Cv值的变化与岩石变形特征的变化有关，是描述前兆分布非均匀性的一种有效指标；Cv值在失稳发生前先异常上升，然后下降，并在Cv值恢复过程中或恢复后伴随着失稳事件的发生．因此，Cv值是一种有效的预报指标．      

Hybrid vibration experiments with a bridge foundation system model

In order to improve seismic design technology of bridges, it is necessary to evaluate the vibration characteristics of a bridge–soil system that consists of soil, foundation structure, pier and superstructure. However, there have been few experimental studies on seismic behavior of bridge–soil system. In this paper, we conducted the hybrid vibration experiment on seismic behavior of bridge–soil system, and examined the applicability of hybrid vibration experiment to study seismic response of bridge–soil system. Based on the experiment results, seismic response of bridge was quantitatively studied.     

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E -polarization induction problems

To improve the accuracy of the numerical evaluation through the 3-D finite difference method, the surface boundary conditions are added to modify the old program. The author has tested the new program by making calculations for the model constructed by Wanamaker, et al (1984). The comparison between the numerical results obtained from this paper and those by Wannamaker, et al (1984) indicates that a pronounced improvement is realized in the evaluation of the horizontal magnetic components. Moreover, better calculations for the vertical magnetic components are also obtainable by using the new program.