Directions and possibilities of mathematical geology

This paper considers the present state of mathematical geology. Three directions are recognized: applied, theoretical, and mathematical. Applied mathematical geology includes formal use of mathematics to solve problems and computer processing of data. Success is achieved by a correspondence of mathematical methods used to the nature of geological data. This correspondence can be demonstrated by purely mathematical means. Theoretical mathematical geology uses mathematics as a language of geology; however, a number of methodological problems must be solved: formalization of initial geological concepts and creation of a strict conceptual basis, substantiation of initial principles of mathematical simulation, creation of theoretical geological models, problems of elementary and coincidence in geology, and methodological substantiations of possibilities of any mathematical model to approximate geological models. The essense and significance of these problems are considered. The main task of mathematical geology is to prove its correspondence to the nature of the geological objects studied, geological data obtained, and geological problems solvable. Finally, the main problems of mathematical geology are not so much mathematical as geological and methodological.     

岩体三维节点偶对分析及其应用

本文讨论了节点偶对分析的三维扩展。节点偶对分析是针对摩擦滑动节理单元的约束边界进行的序列矢量判定。在三维节理单元的分析中,几何约束和相应的力学判据极为复杂,采用节点偶对等效结合力方法才能迎刃而解。文中对三维转向节理模型的分析得到合理的结果,它表明本文所提出的方法可成功地应用于节理岩体的有限元分析。     

地壳演化与成矿作用——以川滇地洼系“四层楼”铜矿床序列为例

提出川滇地洼系“四层楼”铜矿床序列的形成与陆壳演化的成生联系,是与本区陆壳由前地槽—地槽—地台—地洼演化各阶段与之相匹配的成矿作用的产物.与此同时,并总结了本区“四层楼”铜矿床序列的成矿作用具有明显的继承性、新生性、旋回性及层控性四大特点和多因复成矿床的成矿模式.     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

电感耦合等离子体发射光谱法同时测定玻璃中铝钙铁钾镁钠钛硫

通过选择分析谱线、处理样品方法和消除干扰因素等实验,建立了电感耦合等离子体发射光谱法同时测定玻璃中的Al2O3、CaO、Fe2O3、K2O、MgO、Na2O、TiO2和SO3的方法,克服了利用常规化学法测定玻璃中各氧化物步骤繁琐、耗时长、工作量大的不足。方法的回收率为95.0%~103.0%,精密度(RSD,n=10)为0.20%~1.72%。方法具有快速、简便、线性范围宽等优点,分析误差满足常规化学分析法的要求。用于钠钙硅玻璃及其制品的分析,结果令人满意。     

碱消解-高效液相色谱-电感耦合等离子体质谱法测定生物样品中的甲基汞和乙基汞

建立了碱消解-高效液相色谱-电感耦合等离子体质谱联用系统测定生物样品中甲基汞(MeHg)与乙基汞(EtHg)的分析方法。为提高灵敏度,选用微流量的PFA雾化器,在优化的检测条件下,MeHg及EtHg检出限可达到0.036μg/L和0.03μg/L;线性范围达到4个数量级,两条工作曲线线性相关系数为1。对1.78μg/L MeHg、1.65μg/L EtHg的混合标准溶液重复测定7次,色谱峰面积的相对标准偏差(RSD)分别为1.79%和1.44%。对标准物质BCR 464(金枪鱼)的分析结果表明,测定值与标准值基本吻合,但略低于标准值;甲基汞和乙基汞的加标回收率分别为85.9%和84.5%。高效液相色谱与质谱联用技术的高灵敏度和低检出限能够满足生物样品中汞形态定量分析的要求。     

等离子体质谱法测定玄武岩中微量元素三种样品预处理方法的比较

在测定玄武岩的微量元素时,为了保证得到的微量元素结果反映源区特征,需经预处理去除玄武岩的孔洞中充填的碳酸盐等物质。选取雷州半岛新生代玄武岩三个样品,每一个样品用三种预处理方法进行预处理,等离子体质谱法测定其中的微量元素。分析结果表明,用稀硝酸浸泡样品,会导致玄武岩中大多数微量元素严重丢失,得到的数据失真;而用稀盐酸浸泡去除碳酸盐的方法基本不改变玄武岩中微量元素的含量。     

Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China

MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology. Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December 2003. During this period, the total recharge of the groundwater system was 49,374 × 10⁶ m³ and the total discharge was 56,530 × 10⁶ m³ the budget deficit was −7,156 × 10⁶ m³. In this integrated system, the original data including graphs and attribution data could be stored in the database. When the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online.