矿床保存变化研究的热年代学技术方法

成矿后矿床的保存与变化是矿床地质研究的重要组成部分,但研究程度不高,至今仍属于薄弱环节。矿床的保存与其隆升剥露密切相关,深入细致研究恢复矿区矿床的隆升和剥露历史,是揭示矿床保存变化的一种重要途径。本文重点论述研究矿床保存变化的技术方法,通过综合应用裂变径迹、(U-Th)/He、~(40)Ar-~(39)Ar年代学以及La-ICP-MS、锆石U-Pb定年等多种技术手段,研究成矿期次、构造活动期次以及二者间的联系,定量计算不同矿区、不同矿体、不同部位、不同时间的冷却隆升速率、隆升幅度、剥蚀速率和剥蚀量,探讨矿床保存深度与剥蚀量间的关系,总结不同时代、不同矿区矿床保存-变化过程,建立新的矿床地质-保存环境模型,完善矿床预测的综合示踪标志,最终预测不同矿区、不同地段隐伏矿床可能的产出深度,给出矿床可能已经被剥蚀殆尽的区段或地段,为深部找矿和区域成矿潜力评价提供依据。本文可为地质工作者的相关应用提供借鉴和手段。     

MORPHOMETRIC ANALYSIS OF KARST DEPRESSIONS ON A MEDITERRANEAN KARST MASSIF

Dolines are the most typical landform in karst landscapes. They play an important role in karst hydrogeology; they are effective sediment traps that provide information about past climates and karst evolution, and can also be indicators of tectonic activity. The morphometic analysis of those depressions can help in their study and interpretation. However, detailed identification and delineation of dolines becomes increasingly difficult when the study area is very large, of limited access or covered by vegetation. Here, we have used a method developed by the authors for automatic identification and delineation of karst depressions based on the use of a pit removal algorithm. Next a complete morphometric analysis of the karst depressions has been carried out. The main advantage of the methodology is that by delineating and identifying the depth of each doline from its rim, many morphometric statistics can be obtained, including the relation of doline depth and area with respect to altitude, the main doline orientations, doline density, the fractal dimension of doline area and frequency, and point field analysis, among others. The results are interpreted from a geomorphological point of view and the mapped dolines can be easily integrated into geomorphological mapping. A case study is presented from the large Sierra Gorda karst massif in southern Spain, where 3100 depressions have been identified. A strong structural control of doline shape and orientation is demonstrated; the main areas of doline density are highlighted and doline size is shown to follow a fractal law.     

Hybrid N-order Lagrangian Interpolation Eulerian?Lagrangian Method for Salinity Calculation

The Eulerian?Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3?8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

基于非均质地层参数表征的颗粒流模型构建

非均质地层参数的获取与表征方法是浅层地质体改造的重要依据。基于非均质地层的地质建模,表征具有随机不确定性和模糊不确定性的土体表观参数和物理力学参数。利用 Weibull 分布统计描述局部地层信息,并应用Diamond-Square分形插值方法,结合差分盒子维计算非均质地层分形维数,进一步合理演算整体地层参数。采用PFC3D数值模拟软件,结合马尔科夫链蒙特卡洛方法采样数据,对非均质地层块石颗粒及土颗粒进行建模。通过非均质地层模拟结果,发现其破坏过程与应力应变曲线要比均质土体复杂得多,由于块石大颗粒的含量、分布特性对土样力学性质影响显著。本文建模方法有助于非均质地层参数表征研究。     

2011~2020年勘查地球化学数据处理研究进展

我国建立了包含海量数据的高质量的勘查地球化学数据库,为矿产勘查、环境评价和地质调查等提供了重要的数据支撑。如何高效处理勘查地球化学数据,并从中发掘和识别深层次信息一直是勘查地球化学学科研究的热点和前沿领域。本文在系统调研国内外学者过去十年发表的论著基础上,对勘查地球化学数据处理方法进行分析与对比,从勘查地球化学数据库建设、地球化学异常识别及其不确定性评价等方面概述了我国近十年来在该领域取得的主要研究进展,包括:(1)分形与多重分形模型由于考虑了地球化学空间模式的复杂性和尺度不变性,在全球范围内得到极大的发展和推广,我国学者引领了基于分形与多重分形的勘查地球化学数据处理;(2)机器学习和大数据思维开始在该领域启蒙,并迅速得到关注,正在成为研究热点和前沿领域,我国学者率先开展基于机器学习算法的勘查地球化学大数据挖掘研究;(3)我国学者需要进一步加强勘查地球化学数据缺失值处理以及成分数据闭合效应研究。今后该领域应进一步加强对弱缓地球化学异常识别、异常不确定性评价以及异常识别与其形成机理相结合等方面的研究。     

利用毛管压力曲线分形分维方法研究流动单元

利用取心井铸体薄片获得的图像资料和毛管压力曲线,通过图像分形几何学方法以分维数的形式定量地表征出了复杂的微观孔隙喉道结构特征,发现能够很好地划分和评价孔隙岩石中油、气、水的渗流差异,可以用于储层微观流动单元表征。文中阐述了岩石微观孔隙喉道结构分形的理论基础、计算方法和应用于表征流动单元的依据。建立了中国西部砾岩低渗透油藏微观孔隙喉道分维数与孔隙度、渗透率之间计算图版,据此在油藏中利用常规测井资料获得的孔隙度、渗透率参数计算微观孔隙喉道分维数,开展全油藏流动单元划分与评价,取得了良好的效果。研究结果表明,利用毛管压力曲线分形分维方法研究储层微观流动单元是一种很有效的途径。     

陕西关中地区城镇等级规模结构分形研究

文章从分形理论在城镇体系研究中的应用入手,根据豪斯道夫（Hausdorff）维数,运用Zipf公式,对陕西省关中地区的城镇等级规模结构进行了分维测算。结果表明,该区域城镇体系具有分形特征,说明分形理论适合关中地区城镇等级规模结构的测算分析。并在分形研究的基础上分析了该区城镇规模结构的分布特点及其存在的问题。最后基于分形理论的测算与分析结果,提出了关中地区城镇体系结构和功能优化的建议。     

长江三峡巴东复杂斜坡系统成因研究

三峡库区巴东县新城区所在的扇形大斜坡是一个复杂斜坡系统。斜坡后缘边界的挤压变形带是斜坡整体沿其底界面（T1j^3-T2b^1）之间岩层界面向长江发生重力蠕滑作用的结果,不存在“巴东断裂”问题。斜坡系统内部未发现成规模的褶皱,也不存在数量众多的内动力成因的断裂带。巴东斜坡区浅表生变形破坏类型可概括为侧向拉裂、“雪崩式”垮塌、表层风化、侧向滑移张剪破裂、岩层折断变位、软岩膝折剪破裂和压扭破裂又重胶结7种重力卸荷破坏力学模式。斜坡系统内部发育分期分区（块）滑动（滑坡）现象,但不存在构造地质意义上的“断层”形迹。在地貌形态上,巴东大斜坡被4条冲沟分割成5个斜坡单元。在斜坡地质结构、成因类型与空间分布上,巴东大斜坡可划分为3个层次,即表层崩塌滑坡成因为主的堆积层（第一层次）,冲沟分割且浅表生地质形迹发育的层状顺倾的中间基岩层（第二层次）和整体连续顺倾的深层基岩层（第三层次）。分形几何计算证明,巴东大斜坡的地貌形态尚处于侵蚀发育的青（幼）年期。FLAC^3D数值模拟发现,在长江侵蚀下切的不同阶段,巴东斜坡体前缘和后缘接近底界面位置塑性变形区分布集中,但不具备沿深层界面发生整体滑动剪出的可能性。基本认识是,“巴东复杂斜坡系统”是在官渡口-东瀼口向斜南翼（单斜山）的地质背景下,持续经受长江快速侵蚀下切外动力作用,河谷岸坡快速临空导致其自身重力产生强烈的侧向卸荷与滑移等浅表生地质改造作用过程而形成的,可简单地概括为“重力成因论”。     

基于GM(1,1)模型的塔里木盆地南缘孤立绿洲降水灾变预测

我们选取了具有一定代表性的塔里木盆地南缘的孤立绿洲民丰县作为研究区,在对数据有效性的验证之后,利用研究区1980~1999年的降水资料,运用灰色系统理论的季节灾变预测方法,对该地区未来年份的降水量进行预测,并建立灰色季节灾变GM(1,1)预测模型,预测未来旱季到来的年份,为预测未来降水趋势提供了信息。结果表明:该研究区下一个旱灾灾变年大致指向2006年至2007年左右。这对合理利用本地气候资源具有实际意义,并且能够为研究干旱区气候变化提供一定的帮助。     

Hydrothermal alteration processes in the giant Dahutang tungsten deposit,South China: Implications from litho-geochemistry and mass balance calculation

The giant Dahutang tungsten (W) deposit has a total reserve of more than 1.31 Mt WO₃. Veinlet-disseminated scheelite and vein type wolframite mineralization are developed in this deposit, which are related to Late Mesozoic biotite granite. Four major types of alterations, which include albitization, potassic-alteration, and greisenization, and overprinted silicification developed in contact zone. The mass balance calculate of the four alteration types were used to further understanding of the mineralization process. The fresh porphyritic biotite granite has high Nb, Ta, and W, but low Ca and Sr while the Jiuling granodiorite has high Ca and Sr, but low Nb, Ta, and W concentrations. The altered porphyritic biotite granite indicated that the Nb, Ta, and W were leached out from the fresh porphyritic biotite granite, especially by sodic alteration. The low Ca and Sr contents of the altered Neoproterozoic Jiuling granodiorite indicate that Ca and Sr had been leached out from the fresh granodiorite by the fluid from Mesozoic porphyritic biotite granites. The metal W of the Dahutang deposit was mainly derived from the fluid exsolution from the melt and alteration of W-bearing granites. This study of alteration presents a new hydrothermal circulation model to understand tungsten mineralization in the Dahutang deposit.