Gondwanan magmatism with adakite-like signature linked to Cu (Mo)-porphyry deposits from the San Rafael Massif,Mendoza Province,Argentina

The gondwanan magmatism in the San Rafael Massif, known as Choiyoi Magmatic Cycle, was emplaced during the inception of a magmatic arc setting during the early Permian. Two different sections can be differentiated in this volcanic sequence. The lower section (∼281 up to ∼265 Ma) consisting of andesites and dacitic to low-silica rhyolitic ignimbrites has geochemical characteristics that indicate a subduction zone setting. The upper section (∼265 up to ∼252 Ma) composed of rhyolitic ignimbrites and lava flows, dacitic to rhyolitic subvolcanics and alkalic basaltic andesites has geochemical characteristics transitional between subduction and continental intraplate settings. Several Cu–(Mo) porphyry deposits are genetically linked to the lower section (Infiernillo, San Pedro and La Chilca-Zanjón del Buitre). In this paper, we discussed the petrogenesis of the magmatism linked to the porphyry deposits from the San Rafael Massif. The petrogenetic analysis suggests that the lower section was produced in a thickened crust resulting in an adakite-like signature magmatism. The U/Pb LA-ICP-MS age of magmatic zircons from an intrusive associated to the San Pedro porphyry (263.1 ± 4.2 Ma) allowed confirming that the emplacement of Cu–Mo porphyry deposits in the San Rafael Massif occurred during the change in the geodynamical conditions from a transpressive to a transtensive tectonic regime.     

The San Jorge porphyry copper deposit, Mendoza, Argentina: a combination of orthomagmatic and hydrothermal mineralization

The San Jorge porphyry copper deposit (SJPCD) is hosted by Carboniferous clastic sedimentary rocks and Permian intrusions located within the Permo-Triassic belt of Chile and Argentina. Its hypogene mineralization and alteration are products of superposed orthomagmatic and hydrothermal events that were strongly fault controlled. Copper related to orthomagmatic processes includes disseminated chalcopyrite in the matrix of porphyritic granodiorite and andesite, and chalcopyrite with tourmaline and quartz in breccias, both of which have accompanying potassic alteration. Soon thereafter, disseminated chalcopyrite is associated with a structurally controlled silicification of the sedimentary sequence. Finally, multiple episodes of hydrofracturing, probably driven by a deep-seated intrusion, deposited sulfide minerals in veinlets throughout the sedimentary sequence; the centers of these systems are characterized by potassic alteration. Total sulfides, which include chalcopyrite, pyrite, arsenopyrite, and pyrrhotite, and pyrite:chalcopyrite form a linear NNE trend, parallel to the main faults. Quartz–sericite is the dominant alteration and is ubiquitous. Zones of potassic alteration can be delineated even though phyllic alteration can be superposed. Much of the system evolved under reducing conditions. Despite uplift along a reverse fault during the Tertiary, and subsequent erosion, the system is preserved at high levels. Supergene processes redistributed copper in secondary oxides and sulfides. These processes were more effective where the deposit is covered by unconsolidated alluvial sediments. The unique history of the San Jorge deposit renders it an important variation of porphyry copper-style mineralization.     

A statistical assessment of the uncertainty in a 3-D geological framework model

Three-dimensional framework models are the state of the art to present geologists’ understanding of a region in a form that can be used to support planning and decision making. However, there is little information on the uncertainty of such framework models. This paper reports an experiment in which five geologists each produced a framework model of a single region in the east of England. Each modeller was provided with a unique set of borehole observations from which to make their model. Each set was made by withholding five unique validation boreholes from the set of all available boreholes. The models could then be compared with the validation observations. There was no significant between-modeller source of variation in framework model error. There was no evidence of systematic bias in the modelled depth for any unit, and a statistically significant but small tendency for the mean error to increase with depth below the surface. The confidence interval for the predicted height of a surface at a point ranged from ±5.6 m to ±6.4 m. There was some evidence that the variance of the model error increased with depth, but no evidence that it differed between modellers or varied with the number of close-neighbouring boreholes or distance to the outcrop. These results are specific to the area that has been modelled, with relatively simple geology, and reflect the relatively dense set of boreholes available for modelling. The method should be applied under a range of conditions to derive more general conclusions.     

基于线框架模型的三维地质断层结构模型及其构建技术

基于地质断层构造的特点,文中提出了一种适合三维断层模拟的数据模型——线框架模型,并基于该数据模型详细描述三维断层结构模型的构建方法,实现了从层面结构的地质模型向块体结构的地质模型进行自动转换的关键算法。线框架模型以方向线框架来表达地质模型的骨架,以方向TIN网来表达地质模型的血肉,具备地质断层结构模型所要表达的各种几何元素与拓扑关系。基于线框架模型的三维地质断层构模技术将地质界面之间的交切在地质界面三角网形成之前就进行处理,将面与面之间的切割操作简化为线与线之间的操作,这种方法不仅能有效地减轻断层地质体构模的复杂程度,而且可以将线与线之间的交线作为地质界面三角网构建时的约束边,从而有效地保证了地质界面之间交切关系的正确性。     

An index to represent lateral variation of the confidence of experts in a 3-D geological model

A Confidence Index is proposed that expresses the confidence of experts in the quality of a 3-D model as a representation of the subsurface at particular locations. The Confidence Index is based on the notion that the variation of the height of a particular geological surface represents general geological variability and local variability. The general variability comprises simple trends which allow the modeller to project surface structure at locations remote from direct observations. The local variability limits the extent to which borehole observations constrain inferences which the modeller can make concerning local fluctuations around the broad trends. The general and local geological variability of particular contacts are modelled in terms of simple trend surfaces and variogram models. These are then used to extend measures of confidence that reflect expert opinion so as to assign a confidence value to any location where a particular contact is represented in a model. The index is illustrated with an example from the East Midlands region of the United Kingdom.     

三维质点追踪模型及其在胶南海域的应用

在已建立的变边界河口、陆架、海洋模型(ECOM)的基础上添加一个变边界质点追踪模块,研究胶东半岛南海域的物质输运特征,讨论了三维空间中,水质点在潮流作用下的迁移特点.该模型在近海海洋环境管理上具有实用价值.     

长江口斜压诊断模式三维流场数值模拟

建立了σ坐标系下长江口斜压诊断模式三维流场数学模型,采用k-kl二方程紊流闭合模型求解垂向涡粘系数,计算域内恒定、非均匀盐度场反映了计算域密度斜压效应。验证结果表明,斜压诊断模式的模拟成果较切合实际。模型成功复演了长江口水域由往复流向旋转流过渡的流场特性;潮波从四条汊道传入后在分汊口附近相遇形成多个会潮点,其位置随着径流的大小、潮汐的强弱、汊道的形态等因素而变化。     

基于三维地层模型的岩体质量可视化分级

岩体质量分级是岩石力学学科研究的热点和难点,在水电工程建设中具有重要的实际意义。在对岩级分类方法进行大量研究之后,工程人员对岩级结果的表达有了更高的要求,传统的方式已不能满足人们空间分析的需要,研究岩体质量可视化分级十分必要。首先,利用NURBS工具以一种新的数据结构和方法重建三维地层模型,然后,基于岩体质量分级结果,按照空间图形集合运算算法生成相应的三维岩级模型。将上述理论方法应用到锦屏一级水电工程岩体质量分级评价中,并进行了一系列与工程相关的岩级可视化分析,获得了很好的应用效果。     

岩体三维非稳定渗流模型及数值模拟

地下水是导致地质灾害发生的重要因素之一。岩体本身的复杂性及地下水在岩体内随时间而变化造成地下水在岩体内运动机理的复杂性。文章采用数值模拟方法探索地下水在岩体内的运动规律,了解其致灾机理。以岩体结构力学为基础,运用现代分形理论,建立主干裂隙分形网络;根据裂隙发育规模与工程尺度关系将岩体看作拟连续介质与块裂介质混合介质,并根据两类介质接触处水头相等及节点流量相等建立合理的三维非稳定渗流模型;给出渗流模型的有限元解法,开发出相应的有限元软件;给出算例,计算结果体现了主干裂隙在渗流中的强导水作用及网络状裂隙的贮水功能与渗流滞后效应;强调根据裂隙发育规模与工程尺度关系确定合理渗流模型的重要性。     

斜坡上异重流的三维数值模拟

针对异重流的流动特征,建立了适用于具有各向异性浮力紊动特征的三维异重流运动的数学模型,并模拟了异重流在15°斜坡底面上的潜行过程。计算结果准确地模拟了异重流的运动特征和形态,其前锋的潜行速度与实验结果相当吻合。该模型采用非结构同位网格上的SIMPLEC算法能适应复杂边界和地形,可应用于自然界实际环境中异重流的演进计算。     

3-D numerical models of viscous flow applied to fold nappes and the Rawil depression in the Helvetic nappe system (western Switzerland)

The Helvetic nappe system exhibits three-dimensional (3-D) features such as the lateral variation in geometry between the Morcles and Doldenhorn fold nappes or the Rawil depression. We perform 3-D finite element simulations of linear and power-law viscous flow to investigate fold nappe formation during shortening of a half graben with laterally varying thickness. 3-D ellipsoids and corresponding 2-D intersection ellipses are used to quantify finite strain. Fold nappes which formed above a thicker graben have (i) larger amplitudes, (ii) a less sheared and thinned overturned limb, and (iii) a larger thickness than fold nappes formed above a thinner graben. These results agree with observations for the Morcles and Doldenhorn nappes. We also perform 3-D simulations for a tectonic scenario suggested for the evolution of the Rawil depression. The basement is shortened and extended laterally and includes a graben which is oblique to the shortening direction and acts as mechanical weak zone. The graben causes laterally varying basement uplift generating a depression whose amplitude depends on the graben orientation and the stress exponent of basement and sediments. The axial plunge of the depression is smaller (approximately 10°) than the observed plunge (approximately 30°) indicating that additional processes are required to explain the geometry of the Rawil depression.     

滑坡体三维地质建模与可视化分析

针对滑坡地质研究的自身特点,提出了面向滑坡地质体三维建模的NURBS-BRep混合数据结构和地质结构单元实体构造技术.通过对滑坡基础地质数据的预处理、滑坡数字地形和滑动面的三维建模、地质结构面的系统构造、地质结构体的生成和显示,形成了一套完整的滑坡三维地质模型的构建方法.将该方法应用于某水库滑坡,建立了相应的三维滑坡地质模型,并基于此模型进行了三维剖切分析、数字钻孔、等值线生成、滑块自动剖分、滑坡失稳可视化动态模拟和滑坡体方量精确计算等一系列实用的可视化分析,为滑坡稳定性的准确计算和客观评价提供了有力的支持.     

Geodynamic evolution of the northern Molucca Sea area (Eastern Indonesia) constrained by 3-D gravity field inversion

In order to better understand the tectonic framework of the Northern Molucca Sea area, we inverted satellite and sea-surface gravity data into an iterative scheme including a priori seismological and geological data. The resulting 3-D density model images the various tectonic units from the surface down to 40 km. We proceed to various tests to assess the stability and robustness of our inversion. In particular, we performed an offset and average smoothing method to properly refine our results. The resulting model shows a striking vertical regularity of the structures through the different layers, whereas the density contrasts appear strongly uneven in the horizontal direction.The density model emphasizes the complexity of the upper lithospheric structure in the northern Molucca Sea, which is clearly dominated by the interaction between ophiolitic ridges, sedimentary wedges and rigid blocks of the Philippine Sea Plate. It also provides new, hard information that can be used in discussion of the evolution of the region.Large density variations are concentrated in the central part of northern Molucca Sea and dominate the upper lithospheric. North–south trending density structures along the Central Ridge and west dipping thrust faults on the western side of the region are clearly imaged. In the eastern part of the region, we distinguish several blocks, especially the Snellius Plateau which seems to be split into two parts. We interpret this as an oceanic plateau associated with thicker crust that previously belonged to the Philippine Sea Plate. This crust is now trapped between the Molucca Sea complex collision zone and the Philippine Trench, due to the development of a new subduction zone in its eastern side.     

Building a 3-D model of a flood basalt: an example from the Etendeka, NW Namibia

Flood basalts represent large outpourings of lavas which often cover, and interact with, sedimentary basins. For this reason areas with significant flood basalt cover are often targets for hydrocarbon exploration. Problems exist, however, when trying to image sediments and structures in offshore regions covered by basalts. Here we present preliminary 3-D models of the Etendeka flood basalt province from NW Namibia, which can act as an aid in understanding the internal and external architecture of the flood basalt cover. Satellite images, digital elevation models, measured geological logs, sections and maps, are used to create the geological model. Models are presented in 2 parts; 1) models created using topography with images such as Land Sat and geological maps draped over them, and 2) a 3-D model of key lava and sediment surfaces in the basin as defined by measured geological sections. Initial results show a palaeo-volcanic feature early in the flood basalt history which is onlapped by later Iavas. The modelling also allows a simple correction for post emplacement subsidence by assuming an original sub horizontal position for the first basin wide silicic flow unit.     

孤立源的密度、磁化强度及其埋深估算的3维位场方法研究

一个3维均匀源的密度、磁化强度及其埋深,可通过一种非线性、非迭代的3维位场方法反演获得。通过对该方法的研究,实现了对其编程,并将该方法应用于模型反演,取得了很好的反演效果。     

The San Blas Pluton: An example of Carboniferous plutonism in the Sierras Pampeanas, Argentina

New geochronological and geochemical data are reported for the San Blas Pluton (SBP), in the northwestern Sierra de Velasco, Sierras Pampeanas, which intrudes Ordovician granitoids developed during the Famatinian orogeny. A precise Carboniferous age of 340±3 Ma is established by U–Pb dating of zircon using a sensitive high-resolution ion microprobe (SHRIMP). The SBP illustrates several petrological and geochemical characteristics of previously reported Carboniferous granitoids in the Sierras Pampeanas. Their generation is consistent with a regional reheating of the crust at approximately 342 Ma, which resulted in the formation of relatively large amounts of granitic melts that were emplaced in higher crustal levels along master fractures (older master shear zones of Lower Paleozoic age). The SBP can be chemically defined as a typical A-type granitoid related to postcollisonal or postorogenic magmatism. Its high REE content and extraordinarily high U and Th concentrations may have economic significance. Many previously published Devonian and Carboniferous K–Ar dates are reset Ordovician ages, but the existence of other Carboniferous bodies in the Sierra de Velasco cannot be discounted until detailed mapping of the whole Sierra is completed.     

典型窝崩三维数值模拟

窝崩是冲积河流常见的一种崩岸类型,具有发展速度快、破坏力强等特点,且不易成功模拟。以长江下游扬中河段指南村窝崩为例,基于窝崩机理及其力学模式,采用考虑回流区紊动影响的挟沙力公式模拟窝塘内泥沙输移,初步建立了窝崩三维数值模拟方法。研究结果表明:窝崩发生初期窝塘内局部水流具有底层流速大于上层流速的特征,这为窝塘底部及坡脚淘刷提供了动力条件;窝塘内回流区水流挟沙力的计算对窝崩形态模拟结果存在显著影响,考虑回流区紊动影响的泥沙输移特性后,模拟结果与实测结果基本吻合,可较好地模拟窝崩的快速发展过程。研究成果可为窝崩机理的深入认识和窝崩治理提供科技支撑。     

拉萨地块西部亚热地区林子宗群火山岩地球化学特征及其构造意义

广泛分布于拉萨地块南部的林子宗群火山岩被普遍认为是新特提斯洋北向俯冲消减结束并由同碰撞向后碰撞构造背景转化过程中的岩浆作用产物,记录了从新特提斯洋俯冲结束到印度-欧亚大陆碰撞-后碰撞的丰富信息。本文对采集于拉萨地块西部亚热地区的林子宗群火山岩样品进行了年代学和岩石地球化学研究,结果显示:林子宗群火山岩为一套钾玄质系列的流纹岩,该岩石具有高硅(w(SiO₂)=73.90%～74.03%)、高钾(w(K₂O)=6.15%～6.28%)、高钾钠比(w(K₂O)/w(Na₂O)=2.14～2.65),富集K、Rb、Th、U、Hf等元素,亏损Nb、Ta、P、Ti、Sr等元素,具有明显的Eu负异常;林子宗群火山岩样品LA-ICP-MS锆石U-Pb年龄为50.1 Ma±0.3 Ma—51.3Ma±0.4Ma之间;结合前人研究,笔者认为亚热地区林子宗群火山岩为形成于印度-欧亚大陆同碰撞背景下的大陆地壳部分熔融的产物。     

典型窝崩三维数值模拟

窝崩是冲积河流常见的一种崩岸类型，具有发展速度快、破坏力强等特点，且不易成功模拟。以长江下游扬中河段指南村窝崩为例，基于窝崩机理及其力学模式，采用考虑回流区紊动影响的挟沙力公式模拟窝塘内泥沙输移，初步建立了窝崩三维数值模拟方法。研究结果表明：窝崩发生初期窝塘内局部水流具有底层流速大于上层流速的特征，这为窝塘底部及坡脚淘刷提供了动力条件；窝塘内回流区水流挟沙力的计算对窝崩形态模拟结果存在显著影响，考虑回流区紊动影响的泥沙输移特性后，模拟结果与实测结果基本吻合，可较好地模拟窝崩的快速发展过程。研究成果可为窝崩机理的深入认识和窝崩治理提供科技支撑。     

The geological and geodynamic evolution of the eastern Black Sea basin: insights from 2-D and 3-D tectonic modelling

Subsidence mechanisms that may have controlled the evolution of the eastern Black Sea have been studied and simulated using a numerical model that integrates structural, thermal, isostatic and surface processes in both two- (2-D) and three-dimensions (3-D). The model enables the forward modelling of extensional basin evolution followed by deformation due to subsequent extensional and compressional events. Seismic data show that the eastern Black Sea has evolved via a sequence of interrelated tectonic events that began with early Tertiary rifting followed by several phases of compression, mainly confined to the edges of the basin. A large magnitude (approximately 12 km) of regional subsidence also occurred in the central basin throughout the Tertiary. Models that simulate the magnitude of observed fault controlled extension (β=1.13) do not reproduce the total depth of the basin. Similarly, the modelling of compressional deformation around the edges of the basin does little to enhance subsidence in the central basin. A modelling approach that quantifies lithosphere extension according to the amount of observed crustal thinning and thickening across the basin provides the closest match to overall subsidence. The modelling also shows that deep crustal and mantle–lithosphere processes can significantly influence the rate and magnitude of syn- to post-rift subsidence and shows that such mechanisms may have played an important role in forming the anomalously thin syn-rift and thick Miocene–Quaternary sequences observed in the basin. It is also suggested that extension of a 40–45 km thick pre-rift crust is required to generate the observed magnitude of total subsidence when considering a realistic bathymetry.