内蒙古东部近54年气候变化对生态环境演变的影响

利用4个代表站1951—2004年温度和降水资料，通过气候统计诊断分析方法，从气候要素变化趋势、基本气候状态特征以及生态环境演化趋势等方面，研究气候变化对生态环境演变的影响。结果表明：（1）近54年来东部区温度呈显著的升高态势。降水量波动性较大，总体上呈缓慢的增加趋势，但趋势不明显，属于气候自然波动的范围，而1998年至本世纪降水量呈下降趋势，并存在11年和22年的周期性变化。（2）无论是温度还是降水，目前都处于高的气候基本态和高气候变率时期，致使极端气候事件呈增加趋势。气候的变干、变暖以及变率的增大，加剧了内蒙古东部地区干旱化的程度。由此引起一系列生态环境问题，如湿地萎缩、草场退化等。（3）从示范区卫星遥感资料获悉，固定沙地面积的增加，说明科尔沁沙地沙漠化整体出现生态良性好转，但各类沙地图斑数的增加，表明在目前的气候条件下，沙地破碎化程度加剧，沙漠化局部有所恶化。     

Metallogenic age and metallogenic dynamic setting of the Boziguoer super-large rare-rare earth metal deposit in Xinjiang: Insights from monazite U-Th-Pb dating and graphite Laser Raman spectroscopySCIEI北大核心CSCD

波孜果尔碱性岩型Nb-Ta-Rb-Zr-REE矿床位于塔里木克拉通北缘与南天山造山带的接合地段,是近年来在塔里木北缘碱性岩带中发现的规模最大的稀有-稀土综合性矿床。本文通过对含矿岩体中的独居石和晶质石墨包体开展成矿年代学和岩浆形成条件研究,结合含矿碱性岩全岩成分,揭示波孜果尔碱性岩型稀有-稀土金属矿床的岩石地球化学特征、晶质石墨对成岩/成矿的指示、稀土成矿时代和成矿动力学背景,并进一步总结区域成矿指示意义。研究发现,波孜果尔碱性岩以富含碱性暗色矿物(钠铁闪石和霓石)为特征,岩性以正长岩为主,有较高的稀土总量,富集轻稀土元素和高场强元素(Nb、Ta、Zr、Hf、Th、U等),亏损重稀土元素和Ba、Sr等大离子亲石元素,显示强烈的负Eu异常。通过对晶质石墨的研究发现,不同岩相组中石墨的拉曼谱峰相似,显示具有较高的结晶程度,岩浆演化从细晶到伟晶阶段,石墨的无序度增大,可能与烧绿石和锆石等含放射性元素矿物增多使得辐射加强的影响有关。估算石墨的结晶温度介于700~800℃之间,其形成早于碱性岩浆,指示碱性岩浆早期为还原环境,晚期趋于氧化。独居石LA-ICP-MS U-Th-Pb结果显示波孜果尔稀土矿化时代约为280Ma,结合前人已获得的锆石U-Pb年龄,表明波孜果尔碱性岩浆活动时间从~290Ma到280Ma,持续了至少10Myr。进一步研究发现,从时间和空间上,波孜果尔含矿碱性岩所在的碱性岩带形成与塔里木地幔柱有关。地幔柱就位诱发被先前俯冲物质交代富集稀有-稀土元素和挥发性成分的岩石圈地幔低程度部分熔融形成富成矿元素的碱性岩母岩浆,随后在岩浆演化晚期富集成矿。     

Spatial distribution pattern of ore forming elements in the Bayan Obo deposit and exploration implicationsSCIEI北大核心CSCD

白云鄂博矿床是全球最大的稀土-铌-铁矿床。本文对白云鄂博矿床主矿、东矿赋矿白云岩、不同类型稀土-铌-铁矿石和富钾板岩进行详细的野外地质调查和地球化学成分分析,揭示了成矿元素稀土、铌和主要蚀变元素钠、钾、氟、磷的空间分布规律。结果表明,白云鄂博主采坑比东采坑具有更高的稀土含量,矿体由外围向中心、由浅部向深部有明显的稀土元素(尤其是中-重稀土元素)富集的现象,矿体深部具有更大的中-重稀土元素找矿潜力;成矿元素铌的异常范围与重稀土的异常范围相吻合,表明二者具有密切的共生关系;蚀变元素氟的异常范围与稀土-铌-铁矿体的分布范围相一致,是重要的沉淀成矿元素;蚀变元素磷的异常要明显超过稀土-铌-铁矿体的分布范围,是重要的迁移载体元素;钠化蚀变主要出现在矿体下盘含稀土白云岩中,代表早期高温蚀变过程;钾化蚀变主要发育在矿体上盘的富钾板岩中,代表晚期低温蚀变过程;钾化、萤石化及磷的异常是碳酸岩型稀土矿床最主要的浅部和外围示矿指标。     

Sharp and laterally constrained multitrace impedance inversion based on blocky coordinate descent

Seismic impedance inversion is a well-known method used to obtain the image of subsurface geological structures. Utilizing the spatial coherence among seismic traces, the laterally constrained multitrace impedance inversion (LCI) is superior to trace-by-trace inversion and can produce a more realistic image of the subsurface structures. However, when the traces are numerous, it will take great computational cost and a lot of memory to solve the large-scale matrix in the multitrace inversion, which restricts the efficiency and applicability of the existing multitrace inversion algorithm. In addition, the multitrace inversion methods are not only needed to consider the lateral correlation but also should take the constraints in temporal dimension into account. As usual, these vertical constraints represent the stratigraphic characteristics of the reservoir. For instance, total-variation regularization is adopted to obtain the blocky structure. However, it still limits the magnitude of model parameter variation and therefore somewhat distorts the real image. In this paper, we propose two schemes to solve these issues. Firstly, we introduce a fast algorithm called blocky coordinate descent (BCD) to derive a new framework of laterally constrained multitrace impedance inversion. This new BCD-based inversion approach is fast and spends fewer memories. Next, we introduce a minimum gradient support regularization into the BCD-based laterally constrained inversion. This new approach can adapt to sharp layer boundaries and keep the spatial coherence. The feasibility of the proposed method is illustrated by numerical tests for both synthetic data and field seismic data.     

Connectivity constrained segmentation of geologic features

Segmentation of geologic features plays a significant role in seismic interpretation. Based on the segmentation results, interpreters can readily recognize the shape and distribution of geologic features in three-dimensional space and conduct further quantitative analysis. Usually, there are mainly two steps for the segmentation of geologic features: the first step is to extract seismic attributes that can highlight the occurrence of geologic features, and the second step is to apply the segmentation algorithm on the seismic attribute volumes. However, the occurrence of geologic features is not always corresponding to the anomaly value on the seismic attribute volumes and vice versa because of several factors, such as noise in the seismic data, the limited resolution of seismic images and the limited effectiveness of the seismic attribute. Therefore, the segmentation results, which are generated solely based on seismic attributes, are not sufficient to give an accurate depiction of geologic features. Aiming at this problem, we introduce the connectivity constraint into the process of segmentation based the assumption that for one single geologic feature all of its components should be connected to each other. Benefiting from this global constraint, the segmentation results can precisely exclude the interference by false negatives on seismic attribute volumes. However, directly introducing the connectivity constraint into segmentation would face the risk that the segmentation results would deteriorate significantly because of false positives with relatively large area when the connectivity constraints are enforced. Therefore, based on the seismic attribute that highlights the boundary of geologic feature, we further propose a post-processing technique, called pruning, to refine the segmentation results. By taking the segmentation of the channel as an example, we demonstrate that the proposed method is able to preserve the connectivity in the process of segmentation and generate better segmentation results on the field data.     

基于微分响应的流域产流分单元修正方法

为考虑洪水预报误差的空间变化,提出一种基于微分响应的流域产流分单元修正方法.该方法建立了各单元流域产流与流域出口流量之间的微分响应关系,采用正则化最小二乘法结合逐步迫近进行反演求解,将产流误差估计量分配给相应单元流域实现流域产流分单元修正.将构建的方法应用于大坡岭流域和七里街流域进行新安江模型产流修正,比较分析了流域产流分单元修正、流域面平均产流修正和自回归修正的效果.结果表明：流域产流分单元修正效果优于流域面平均产流修正;随着预见期的增大,产流微分响应修正效果优于自回归修正.该方法通过汇流系统将流域出口断面流量信息进行分解用于修正各单元流域产流,有利于提高实时洪水预报精度.     

内蒙古柴胡栏子金矿田燕山期岩浆岩与金矿床

岩相学、岩石化学及稀土元素和微量元素研究证明, 柴胡栏子金矿田的燕山期火成岩组合具有活动大陆边缘弧火成岩组合的特点, 介于典型的靠近海沟的钙碱性弧火成岩与陆内造山火成岩之间的过渡类型。流体包裹体测温 （１８０℃～２９０℃）, 硫同位素 （δ３４ Ｓ＝ ３１‰～７８‰）, 氢氧同位素 （成矿流体以岩浆水为主, 有大气水参与）, 岩石、矿石稀土元素配分特点的相似性等说明金的成矿主要与岩浆作用有关。燕山期岩浆构造作用引起原始矿源层太古界建平群中的金发生活化、迁移, 并在岩浆热液的进一步作用和大气水的参与下, 发生富集沉淀, 形成金矿床。柴胡栏子金矿田的北矿区金矿成矿主要与闪长岩有关, 以蚀变岩型为主。王府矿区金矿化与闪长玢岩脉关系最密切, 以石英脉型为主。岩浆作用为金矿床的形成提供了热源和大部分热液来源, 并且是成矿物质活化和转移的载体。本区金矿床应属与岩浆作用有关的中温热液矿床。     

一种新的断裂计算方法

利用单元中嵌合无厚度不连续面的方法 ,模拟断裂面的形成过程。用虚功原理建立包含断裂面影响的有限元平衡方程 ;利用余虚功原理推导了断裂面的矩阵的具体形式 ; 以第一类断裂问题 -开裂问题为例 ,建立了开裂面粘塑性开裂本构关系,进一步得到了以弹粘塑性形式表示的有限元计算格式。 通过简单算例证明了该方法的正确性。     

物化探方法对古墓完整性的研究

根据湖南春秋战国—秦汉古墓葬地质条件、埋葬特征与周边环境,遵循准确、简便、快捷和经济的原则,选择并运用放射性氡测量、汞测量和探地雷达等方法,对湘西北某县AB战国楚墓群4个大中型古墓和长沙杨家山古墓开展物化探综合勘查试验,取得了古墓完整性勘查研究的初步成功。     

一个大气边界层非线性三维数值模式及其在太湖流域的应用

本模式用以K半经验理论为基础的闭合方程组求解,引进了地形作用、动量、热量和水汽的侧向扩散和垂直湍流交换以及凝结加热等物理过程,其中湍流交换系数K是地表粗糙度、几何高度、大气稳定度及风速切变等因子的函数。本模式还考虑了山地、平原和水面在地形高度、地表粗糙度、辐射、蒸发及热交换等方面的差异,在下垫面上建立热平衡方程与边界层内控制方程耦合。将模式应用于太湖流域,计算得到的该流域边界层内温、压、湿、风的分布特征与实际情况相似,模式具有一定的实用性。