重磁异常的线性反演

本文详细地讨论了磁化强度填图的方法原理,并且通过理论模型和实测异常计算了引起磁异常激发体的磁化强度分布,还讨论了该方法的效果。一、方法原理根据观测到的重磁异常,确定引起这一异常的激发体的各种参数,构成了重磁反问题。对反问题而言,一般说来唯一解是不可能的。这是因为,无限数目的不同密度或磁性分布,都可以引起一个给定的重力或磁力异常。因此,为了使反问题得到唯一解,必须事先附加某种假设。按照给定的附加条件,实际问题可分为两种类型。第一种类型是假设源的几何形状和空间位置已知,根据观测到的异常求解密度或磁化强度的分布;第二种类型则是假设密度或     

地球化学综合异常的圈定及找矿效果

能否正确有效地圈定地球化学综合异常是地球化学数据处理的重要内容,对找矿效果起决定作用。传统方法多在单元素异常的基础上圈定综合异常,而忽略了异常可能由地层本身高背景富集引起。笔者引用马氏距离识别地球化学离群点的方法圈定综合异常。由于马氏距离是基于多元正态分布理论,是多变量统计方法的直接推广和应用,其全面考虑了元素均值、方差及元素间协方差3个参数,是直接针对样本的计算过程,同时也是一元方法的直接推广,所以其具有独特的异常识别功能。经各种不同比例尺和不同采样介质化探数据处理应用显示,利用马氏距离方法求得的地球化学综合异常具边界唯一、界线清楚,减少了人为圈定异常的影响;异常强度突出,指标单一,可作为评价异常的一个重要参数使用;发现矿体(矿化)的概率程度高等优势。该方法可在实际中推广应用。     

超声电视测井图像中异常点剔除算法

针对超声电视测井图像中经常会出现的异常点问题,提出了一种采用全变分的修复算法（TV算法）来剔除异常点。该方法先在选定的区域里根据阈值来确定要剔除的异常点,然后利用全变分的图像修复算法去修复这些异常点。实验结果证明,该方法能够很好地对异常点进行剔除,保证测井图像的完整性,明显地改善图像质量,对实际测井资料的处理具有很重要实际意义。     

河北省铁矿分布与重磁场特征的关系

在河北省铁矿资源潜力评价及成矿预测中,对重力场、磁场特征进行了总结和归纳,认为规模较大的航磁正异常主要是出露或隐伏侵入岩基的反映;大面积分布的低缓正磁异常区则主要反映了太古宙变质岩系或基底隆起的边缘,沉积变质型铁矿主要密集分布在此区域;负磁场区主要出现在中上元古界、古生界碳酸盐岩和中生界酸性火山—沉积岩系;强度高、梯度大、或有明显负值伴生的局部磁异常则可能是磁性铁矿的反映。从重力场特征看,磁性铁矿主要分布在布格重力异常的梯级带或异常扭曲处,以及正、负异常的交界处,这些地方多与地质构造强烈变形密切相关。只有对重磁场进行深入地质剖析,结合典型矿床上建立的识别标志进行对比分析,才能得出一个符合客观实际的地质解释。     

数据断面中突变点的快速剔除方法

针对在高密度电阻率法测量过程中,数据断面易出现突变点,而造成电阻率拟断面图虚假异常,难以准确解释的情况提出了一种快速、高效的剔除突变点的方法。该方法以曲线的曲率公式为理论基础,用相邻点间的半二阶差分对突变点进行自动识别,并将其剔除。充分利用断面数据资料,然后用一元三点插值法或趋势面拟合加残差叠加算法对缺失点进行内插,以补齐数据断面中损失的数据点,使其便于进行定性或半定量解释。最后用该方法对含有突变点的剖面模型数据和野外实测数据断面进行处理,处理结果表明,突变点被完全剔除,但不损失有用信息,易于对其进行推断解释。通过实际应用,该方法具有一定的应用价值,值得推广使用。     

多元地球化学异常识别的核马氏距离方法

地球化学数据满足多元正态分布时,马氏距离是一种有效识别多元地球化学异常的综合指标。然而,由于地质系统的复杂性、成矿作用的多期多阶段性以及控矿因素的多重性常常导致多元地球化异常临界面是非线性的和模糊的,用马氏距离定义的平滑超椭球面不能准确表示这种复杂曲面。核函数能够将地球化学样品集非线性变换至特征空间,背景样品的映像集合在特征空间中构成一种流型,异常样品的映像则零散分布于流型的边缘及外围。计算和比较样品映像到样品映像总体的核马氏距离,可以识别异常样品。把该方法应用于白山地区多元地球化学异常识别,用核马氏距离、马氏距离和主成分得分识别金-银、金-银-砷-铋-汞、金-银-铜-铅-锌-锑-钴、金-银-铜-铅-锌-砷-锑-铋-汞-钴4种组合模式的多元地球化学异常。研究结果表明：复合核函数马氏距离的多元地球化学异常识别效果优于其他方法。     

柳江盆地浅层地下水硝酸盐背景值研究

为探索地下水硝酸盐背景值获取方法,文章以柳江盆地为研究对象,在对比分析国内外研究方法的基础上,首先采用绝对含量和毫克当量百分位数双因子法从宏观上剔除硝酸盐异常数据,然后再利用层次聚类分析结合主成分分析法,分析地下水水化学特征及识别异常分类,进一步剔除异常数据。最后剩余数据进行分布类型检验,采用浓度累计频率法确定地下水硝酸盐背景值范围。研究结果表明,绝对含量和毫克当量百分位数双因子法虽然不能完全剔除异常值,但可为后续层次聚类分析异常识别减少异常信息和子集;层次聚类分析法注重对各亚类的水化学特征分析来识别分析异常数据,具有识别人为异常和天然异常的优势。对比分析常用的数理学方法计算表明,2种方法结合,更能有效识别异常,计算出的地下水硝酸盐背景值更合理。异常数据剔除分析表明,柳江盆地浅层地下水硝酸盐异常与农业化肥的超量施用和居民生活污水与垃圾粪便的下渗污染具有密切的关系。     

航磁异常总梯度模反演方法的实用化改进及软件研制

通过采用网格数据文件按列与分块随机存取技术,使得大数据量网格位场异常总梯度模的计算及数据极值点的自动搜索和连续追踪反演成为可能;提出自动搜索出总梯度模网格数据极大值点后沿异常走向方向连续追踪相邻极大值点串的做法,所获极值点串的连线能够更直观和有效地反映磁性或密度地质体与构造边界(或中心)的水平位置.经实用化改进研制的反演软件,其模型及实际数据计算均获得较好的预期效果.     

相山地区地面高精度磁测数据

本文针对相山火山岩盆地的弱磁异常受地表磁性不均的强烈干扰形成磁场突变点多、幅值高的特点，提出了一种数据预处理方法──剖面三点插值切割法，能有效地剔除单个突变点，使磁异常曲线圆滑，并能保持异常特征，在实际应用中具有一定效果。     

遥感图像像素级异常识别的一种方法

遥感图像异常识别是遥感应用领域一个颇受关注的研究课题,在军事目标识别和自然环境保护等许多领域都有潜在应用价值。不妨假设遥感图像背景像素分布于随空间位置缓慢变化的一系列高斯超椭球体内,异常像素则分布于超椭球体之外。在这种假设前提下,首先应用Weiszfeld方法估算遥感图像中一系列高斯超椭球体的重心和波段协方差矩阵;然后,计算各像素到对应的超椭球体重心的马氏距离,并用直方图法确定马氏距离的异常下限;最后,把马氏距离高于异常下限的像素作为异常像素识别出来。在GDAL遥感图像数据输入输出函数库基础上,用VC++语言开发了遥感图像像素级异常识别的算法程序;用美国亚特兰大TM图像进行了方法的应用实验研究。结果表明,该方法对遥感图像中的局部异常具有很好的识别效果。     

判断正态样本异常值的两种方法及程序

物探、化探和数学地质中一般接触的是未知总体方差的样本。根据理论和经验,这些样本是符合正态分布或近似正态分布的。但是这些样本中,或许会发现个别的可疑数据,它们与样本中其它值比起来偏大或偏小。误差理论认为,这种个别的可疑数据往往是“粗大误差”,或叫“过失误差”、“坏值”,是由实验过程中读错、记错、仪器突然跳动、突然震动等异常情况引起的异常值。显然,这样的异常值若不剔除或改正,势必歪曲实验结果。但另一方面,如果在某地质单元范围内,从采集的数据看,某指示元素的背景值本来符合正态分布,却有个别的值过高,这样的异常值又正是工作的目标和工作继续深入的指示。     

赣东北金矿矿源岩的区域岩石密度特征及找矿方向

对研究区内金矿岩双桥山群的4336件岩石密度数据的统计分布、空间制图分析。得到岩石密度参数的地质制图，发现双桥山群岩石密度为偏态分布，推测岩石密度的增高可能是变质作用或成矿作用增强的反映，区内主要岩金矿床均分布在密度高异常区内，因此认为高密度异常区或重力正异常区，是区内岩金矿床的找矿远景区。     

多元正态总体假设检验在矿化带识别中的应用

根据北瓦地区岩心样品的地球化学数据处理结果,介绍多元正态总体假设检验方法与对应分析相结合在该区铅锌矿化带识别中的应用。从全体数据的对应分析结果出发,将见矿钻孔的样品得分值经稳健处理后作为铅锌矿化带的已知样本,并以马氏距离为统计量对未见矿钻孔的样品进行识别,在一定信度下作出了该区已知铅锌矿化带在深部存在较有意义延伸的统计推断。     

用SHRIMP测定锆石U-Pb年龄的工作方法

用SHRIMP测定锆石U-Pb年龄涉及锆石单矿物分选、样品靶制备、U-Pb同位素组成测定、原始数据处理、数据解释等多方面。锆石分选中要避免不同样品的交叉污染,要尽量将岩石中所有的锆石分选出来;淘洗岩石粉末时要避免丢失细小的锆石,因细小的锆石更重要。制作样品靶时,粘贴到样品靶的锆石应反映岩石中锆石的全貌。数据处理中,要修正或剔除异常数据、评估拟参加加权平均计算的每个数据权重的合理性。通常用测得的众数年龄的加权平均值作为侵位年龄或喷发年龄会老于实际年龄,应使用样品中最年轻的锆石年龄来限定岩浆侵位或喷发年龄。     

在复杂磁异常地段找矿的几点体会

在无产阶级文化大革命中,我队广大职工与兄弟队一起,以毛主度光辉哲学思想为指导,破除条条框框,大胆实践,在晋南一个被下了"矿体小而分散"结论的地区找到了富铁矿,使被修正主义路线砍下马的矿山得到新生.该区磁异常极为复杂,异常值有时为正,有时为负;有时明显,有时不明显.有人说该区"正异常有矿,负异常也有矿;有伽马有矿,无伽马也有矿".咋听起来使人莫明其妙,然而它正是该区复杂磁异常见矿的实际情况.这并不是说,复杂磁异常没有规律,不能认识.而在于我们的认识和工作是否如实地反映了客观实际,抓住了事物的本质.现将我们在晋南塔儿山地区复杂磁异常中找矿的几点体会叙述如下.     

可疑数据的取舍,格勒布斯(Grubbs)方法

一、引言由试验室提供的土的物理力学性质指标或野外原位测试所提供的测定数据中,有时在同一土层中的一个样本数据中,往往可能有一个或数个过大或过小的数据。过去地质人员根据主观的判断加以取舍。结果,数据的取舍因人而异,缺乏统一的准则。一般讲取自于同一土层的一个样本数据应该都来自同一土层,但是由于试验或测试上的原因,混入了另一个土层的数据;另一种可能是土样取自“过渡层”或二层交界面的情况确实是存在的,对于这种数据的取舍问题过去也往往采用人为的舍弃,这样的数据往往在值上的反映是过小或过大,我们     

地球化学异常识别的两种机器学习算法之比较

在Sklearn的Python语言代码基础上,开发了基于孤独森林和一类支持向量机的多元地球化学异常识别方法程序。选择吉林省和龙地区为实验区,从1∶5万水系沉积物资料中提取地球化学异常。把实验区已知矿点的空间分布位置作为"地真"数据,绘制两种机器学习算法的ROC曲线并计算AUC值,用来对比两种方法的多元地球化学异常识别效果。研究结果表明:两种机器学习算法都能够有效识别多元地球化学异常,所提取的多元地球化学异常与已知矿点具有显著的空间关联性;孤独森林算法在数据处理耗时和多元地球化学异常识别效果方面略优于一类支持向量机。     

单次流量测验中Ⅰ、Ⅱ型误差分布的研究

目前许多数据处理方法和统计推断都是以正态分布为其理论依据的。然而对水文观测数据和误差是否符合上述特性,国内外水文学者都曾有疑义,近年来数理统计学者也注意到当观测值不符合某些特性假定或含有异常点时,误差分布将不是正态,而是所谓的“长     

利用样本排序方法比较化探异常识别模型的效果

地球化学异常的有效识别是化探找矿成败的关键环节.利用样本排序方法对各种化探异常识别模型的处理效果和优劣性进行了比较.以云南个旧及其周边地区铜元素水系沉积物为例, 应用元素含量、异常衬度、局部奇异性分析三大类方法对数据进行了处理, 对化探异常指示变量的排序值进行了3个方面的对比: (1) 在高背景区和低背景区样品的排序特征; (2) 有矿出现位置处样品的样本排序特征; (3) 累积面积(等效于上侧样本排序值) 不同分组所圈定的异常范围和矿床(点) 之间的空间相关性.结果表明, 局部奇异性分析方法较按含量高低的全局性方法对弱缓异常识别更为有效, 也相对优于滑动衬值.按奇异性指数基于证据权法圈定的异常远景区更具预测意义, 显著优于用元素含量值所圈定的异常范围.局部奇异性方法原理清晰、方法简便、可操作性强, 在地球化学异常识别中可以用其替代滑动衬值方法.      

P点铅及其应用意义

在铅同位素研究中，异常铅是经常遇到并难于解释的问题。在一个矿区或矿带范围内，某些矿床的铅同位素资料中往往显示存在两种类型的铅：具有单阶段演化历史的正常铅和具有两个(或多个)阶段演化历史的异常铅。在常规的铅同位素组成图解上，这两种铅的数据点有时构成异常铅等时线，正常铅则位于等时线上含放射性成因铅最低的位置处。如果这种正常铅能给出合理的矿化年龄，并且该年龄与容矿围岩的成岩年龄基本一致，但显著老于异常铅瞬间增长模式年龄，那么，这种正常铅可能是异常铅等时线的起点，笔者称其为P点铅。显然，异常铅是后期放射性成因铅加入到P点铅形成的，P点铅的模式年龄为异常铅来源区的年龄：利用P点铅及异常铅等时线的资料，可以计算获得异常铅的矿化年龄。因此，P点铅这个概念的提出，为探讨矿化年龄问题开辟了一条新途径。文章还以加拿大和北欧地区某些矿床的铅同位素资料为例，讨论了在实际成矿过程中是否有P点铅存在的问题，讨论了P点铅在矿床成因研究中的重要意义。并给出了如何判断P点铅的具体条件。