分形理论在放射性场异常晕圈圈定中应用.doc

分形理论在放射性场异常晕圈圈定中应用 　　[摘 要]分形理论是处理复杂地质地球化学问题的手段之一，文章通过MAPGIS平台，采用分形理论含量―面积法研究某地区土壤氡浓度异常晕圈的圈定，为寻找该地区找矿远景区提供依据 [关键词]分形理论 含量―面积 氡浓度 MAPGIS 数据处理 成矿预测 中图分类号：F270+1 文献标识码：A 文章编号：1009-914X（2017）11-0359-03 1 引言 氡及其子体测量是一种直接寻找隐伏铀矿的放射性勘查方法，氡浓度背景值及异常值的确定是放射性数据处理中的重要工作之一，也是铀矿资源预测及评价的重要依据。传统氡背景值及异常值的确定方法是采用数理统计方法求算术平均值与标准偏差，这种方法是基于氡浓度数据（近似）服从正态分布或对数正态分布。但新的研究表明元素在地球化学场中的分布并不局限于正态分布或对数正态分布[1]，而是具有不规则和自相似性，具有典型的分形特征。地质问题是复杂的非线性问题，分形方法是解决非线性问题的一种有效的新方法，因此分形几何在石油、核能、煤炭勘探中得到了广泛的应用。本文应用分形理论中的元素含量―面积法来研究某地区氡浓度异常下限的确定 2 分形理论与元素含量―面积法 2.1 分形理论 在欧氏空间中，人们习惯把空间看成三维的，平面或球面看成二维，而把直线或曲线看成一维，通常人们习惯于整数的维数。分形理论把维数视为分数，这类维数是物理学家在研究混沌吸引子等理论时需要引入的重要概念。为了定量地描述客观事物的“非规则”程度，1919年，数学家豪斯道夫（Hausdoff）提出了连续空间的概念，也就是空间维数是可以连续变化的，它可以是整数也可以是分数，称为豪斯道夫维数。记作D，一般的表达式为： 其中L为某客体沿其每个独立方向皆扩大的倍数，K为得到的新客体是原客体的倍数。显然，D在一般情况下是一个分数。在此基础上，1975年，Mandelbrot在其《自然界中的分形几何》一书中引入了分形（fractal）这一概念，在数学上分形有以下凡个特点： ⑴具有无限精细的结构； ⑵比例自相似性； ⑶一般它的分数维大子它的拓扑维数； ⑷可以由非常简单的方法定义，并由递归、迭代产生 分形维数是复杂形体的不规则性的量度[1]，对于一个有确定维数的几何形体，用与它维数相同的尺度q去度量，根据豪斯道夫维数定义，其测度大小S（q）与单位量度q之间存在如下关系： 为求出其分形模型中的分维数D，将数据（S（q1），S（q2）…S（qn））和（q1，q2…qn）绘制在双对数坐标图上，如果数据点大致分布在一条直线上，用最小二乘法求出的直线的斜率的估算量可近似看作是分维数D值，如数据大致分布在两段或多段直线上时，采用分段拟合，在各区间上用最小二乘法回归时用最优法确定分界点，即使区间拟合的直线与原始数据之间的剩余平方和Ei（i=1，2，…）在各个区间的总和最小： 让i取遍所有数据点，使E取最小值的点为分界点，D1和D2分别为相应区间的斜率分维数 2.2 元素含量―面积 目前在地球物理学中利用分形技术确定异常下限的方法主要有含量―总量法、含量―面积法、含量―距离法和含量―周长法等。陈秋明等（1994）从分形的观点认识到地球化学背景值和异常值具有各自独立的幂指数关系，由此导致了一种多重分形分布[1]，元素含量―面积多重分形模型可表示为： S（q≤c1）∝，S（c1≤q≤c2）∝，… 其中∝代表正比于，c1，c2表示各分形分布界限点元素含量值（区域异常下限等）。表明元素的含量在空间上分布存在层次性特征，即存在背景分形分布和异常分形分布等多重分形分布 本文采用的元素含量―面积法以野外测量氡浓度值qi（单位Bq/l，贝克每升）为尺度，以MapGIS为平台绘制平面等值线图，然后将等值线所圈闭的区面积导出到Excel中进行统计，以浓度大于qi的等值线圈闭的区域面积之和S（qi）为测度。显然，在分形模型中随着圈定等值线的氡浓度值的升高，等值线所包围的面积将逐渐减小，即S（q）是一个随着q增加而减小的函数 3 工作区成矿地质条件 工作区位于某复式花岗岩体北部。地史上经历了前泥盆纪裂谷海槽演化、泥盆纪―三叠纪中期大陆形成、晚三叠世―新生代陆缘活动等三大发展阶段。出露的岩浆岩为燕山早期第二阶段花岗岩，如图1所示，主要岩性为中粗粒斑状黑云母花岗岩，岩石的含铀量较高，普遍高出铀克拉克值数倍至十几倍，且钍铀比值小，活化铀多，属典型的产铀岩性 区域上经历了多期次构造活动，形成一系列北北东向、北西向、近南北向、近东西向断裂构造。工作区近东西向和北西向两组断裂构造规模较大，是岩体内的主要含矿构造，北东向构造多为规模较小的次级带。本次氡