化探单元素异常统计内容参数公式
地质学中一些公式
地学中常用公式一、平均品位的计算公式:1、算术平均:(X1+X2-……+Xn)/n X1、X2、X n为样品品位2、加权平均:(X l×L l+X2×L2+……+ X n×Ln)/(L l+L2+……+L n) X1、X2……X n。
为样品品位,L l+L2+……+Ln为样品长度3、几何平均为Xn2⨯1 X1、X2、Xn为样品品位X⨯n⨯X注:品位为正态分布时,处理特高品位时,可用此公式。
二、矿体厚度(Vm)、品位(Vc)变化系数:—X=(X1+X2+……+Xn)/n 计算矿体厚度、品位的平均值∑-σ计算均方差X(2nXi/(-=)1)厚度、品位变化系数:Vm或Vc=⨯σ100%÷X三、地质剖面岩石厚度计算公式:y=sinα·cosβ·cosγ±cosα·sinβα--导线坡度角β--地层倾角γ --导线方向与地层倾角的夹角地层倾向与坡向相反取正号,地层倾向与坡向相同取负号;真厚度=L×y四、钻孔矿体厚度的确定矿体的厚度是根据矿体露头上、坑道中和从钻孔中所获得的资料进行的。
(一)坑道中矿体厚度的测定当坑道所揭露的矿体与围岩的接触界线清楚时,取样和编录时可在矿体上用钢尺直接捌量出来。
厚度测量的次数决定于坑道的布置情况,如矿体是用穿脉坑道圈定的,则测量次数与穿脉坑道的数量相符。
如果矿体是用沿脉坑道圈定的,则厚度的测定按一定间隔在取样的位置进行测量。
如果矿体与围岩的界线不清时,矿体厚度的测定必须根据取样结果来确定。
(二)钻孔中矿体厚度的测定因为钻孔中所截穿的矿体均在地下深处、只能间接地去测定矿体的厚度。
当钻孔是垂直矿层钻进时,且岩心采取率为100%,可直接丈量岩心,取得厚度的数据。
若岩心采取率不高,除用钢尺丈量岩心长度外,还要按下式进行换算:m nL(11-9)式中: m ——矿体的厚度(米); L ——实测矿心长度(米)I n ——矿心采取率(%)。
化探异常评序
二、异常评序方法地球化学综合异常评序采用相同主元素类内排序方法,即将同一主元素的异常分别置于一起,依其评序指数(JOI)由大到小进行评序;遇有某一异常内有2个以上主元素者,则分别计算各自的评序指数,参与同类主元素异常的评序。
评序指数(JOI)计算方法:JOI=D×K×C式中D—主元素规模;K—异常内(除主元素外)特征组合元素平均衬度,若是两个主元素在类内评序时,其余一个主元素与其它元素一起参加特征组合衬度计算;C—修正系数。
修正系数C由以下几个方面内容组成:1、主元素异常浓度分带的外、中、内带各给1、2、3分。
2、异常组合特征:组合复杂与主元素密切相关者给3分,组合单一者给1分,介于二者之间给2分。
3、异常套合程度:与主元素异常套合程度好或有一定规律分布者给2分,较好者给1分,不好者不给分。
4、异常与已知矿产吻合程度:反映了已知矿点、矿化点者给2分;不一致者给1分。
5、异常分类类别甲1、乙1、乙2、乙3、丙类分别给4、3、2、1、0.5分。
6、异常区的地层、岩浆岩、大地构造位置等条件,对成矿有利者,按其所处条件好、较好、一般者分别给3、2、1分。
最后将各项指标得分累加,以总分的高低给定修正系数值,按>15、14~11、10~7、<7分别给C赋值1.2、1、0.8、0.6参与运算,以评序指数的高低进行异常排序。
三、异常评序结果区内所有被选综合异常均参与了评序,评序指数(JOI)顺序排列前十位的元素如下:第一为Ni,第二为W,第三为Ag,第四为Sb,第五、第六为Au,第七为Mo,第八、第九为Pb,第十位为Mo,目标矿种显示良好,成矿潜力较大。
(一)Ni异常区内Ni异常主要出现于测区北端奥陶纪地层和东南角的花岗闪长岩脉中,元素组合以NiCr为主,其中北端异常呈NW-SE向分布,和区内主体断裂走向一致,伴生AuAgCu元素组合,具有一定的找矿意义;而东南角和花岗闪长岩体侵入有关。
进入地球化学异常特征组合者3处,成为主元素者为2处,评序结果见表4-4。
化探数据处理方法
内蒙古扎赉特旗东芒合矿和哈拉街吐矿化探数据处理及图件编制方法1 化探数据质量评价的数据处理(分矿区)⑴统计重采样和重分析抽查样所占样品总数的比例比例 = (重采样和重分析抽查样数/工作样总数)100%⑵作出SSPS数据文件将重采样和重分析样分别作成SSPS数据文件。
文件中列出项目为:①重采抽查样重采样号元素含量相应的工作样号元素含量②重分析抽查样重分析样号元素含量相应的工作样号元素含量⑶计算各元素相对误差重采样和重分析抽查样相对误差均按RE(%) = |C1-C2|/0.5×(C1+C2)×100%计算。
C1为重采样或重分析抽查样的分析含量C2为重采样或重分析抽查样的相应的工作样的分析含量| |为绝对值RE(%)≤30%为合格,>30为超差(不合格);(Au:RE(%)≤50%为合格,>50为超差)⑷计算各元素的合格率η= (抽查样品中合格的样品数/抽查样品的总数)100%合格率(η)应>80%,即这批样品的分析结果是可信的。
⑸列表表示检查或分析质量结果表××化探重采样抽查各元素的合格率(%)Cu Pb Zn Cr Ni Co Sn V Ag Ti2 矿区地球化学特征研究的数据处理(以哈拉街吐为例)⑴作出SSPS数据文件作出下列SSPS数据文件:①文件1:整个矿区数据文件;②文件2:矿区地层数据文件;③文件3:矿区岩浆岩数据文件;④文件4 :下二叠统大石寨组(P1d)数据文件;⑤文件5 :下白垩统大磨拐河含煤组(K1d)数据文件;⑥文件6 :华力西晚期侵入岩数据文件;⑦文件7 :燕山期早期侵入岩数据文件;⑧文件8 :燕山期晚期侵入岩数据文件;⑨文件9:已知矿附近一定范围数据文件每一数据文件的内容项目包括:序号野外号 X坐标 Y坐标各元素的含量⑵整个矿区和各地质单元(各地层、各岩浆岩)样品各元素含量特征统计统计的参数包括:①元素含量平均值;②最大值;③最小值;④标准离差;⑤变化系数(标准离差/含量平均值);⑥浓度克拉克值(元素含量平均值/该元素的克拉克值)整个矿区和各地质单元统计结果含量平均值、最小值、最大值用表表示。
地球化学异常下限确定方法
地球化学异常下限确定方法一、地球化学数据处理基础数据处理的意义是获得较为准确的平均值(背景)和异常下限。
1、地球化学数据处理归根结底仍属于统计学的范畴,所以要求数据应是正态分布的,不是拿来数据就能应用的,特别是用公式计算时更要注意这一点。
正态(μ =0, δ =1)----(偏态)。
大数定理:又称大数法则、大数率。
在一个随机事件中,随着试验次数的增加,事件发生的频率趋于一个稳定值;同时,在对物理量的测量实践中,测定值的算术平均也具有稳定性。
所以如果在计算时,数据中包含较多的野值时,实际获得的是一个不具稳定性的算术平均,它实际不能替代背景值。
2、异常是一个相对概念,有不同尺度上的要求,所以不要将其看作一个定值。
在悉尼国际化探会议上(1976),对异常下限定义:异常下限是地球化学工作者根据某种分析测试结果对样品所取定的一个数值,据此可以圈定能够识别出与矿化有关的异常。
并对异常下限提出了一个笼统的定义:凡能够划分出异常和非异常数据的数值即为异常下限。
据此,异常下限不能简单的理解为背景上限。
二、异常下限确定方法具体异常下限确定方法较多:地化剖面法、概率格纸法、直方图法、马氏距离法、单元素计算法、数据排序法、累积频率法……下面逐一介绍:1、地化剖面法:(可以不考虑野值)在已知区做地化剖面:要求剖面较长,穿过矿化区(含蚀变区)和正常地层(背景),能区分含矿区和非矿区就可确定为下限。
2、概率格纸法:(可以不考虑野值)以含量和频率作图15%--负异常50%--背景值85%--X+δ(高背景)98%-- ( X+2δ)异常下限3、直方图法:(可以不考虑野值)能分解出后期叠加的值就为异常下限4、马氏距离法:(在计算时已考虑野值)针对样本,实际为建立在多元素正态分布基础之上—多重样本的正态分布,超出椭球体时—异常样(如P3点)。
相似于因子得分的计算,最后为一个剔除异常样本时的计算值,实际计算出综合异常边界线。
当令m=1时,上式化解为Xa=Xo?KS,这是我们较为熟悉的单元素(一维)计算异常下限常用公式。
土壤化探中异常下限的确定(可编辑)
土壤化探中异常下限的确定土壤化探中异常下限的确定摘要土壤地球化学异常下限的确定是勘查地球化学的一个基本问题,也是勘查地球化学应用于矿产勘查时决定成败的一个关键性环节。
但由于地质背景和成矿模式的复杂多样,迄今为止仍然没有一种普遍适用的异常下限计算方法诞生,各种计算方法各有优势,同时又有假设条件的制约和使用的局限性。
为此,采取多种方法计算异常下限并根据地质背景进行综合比较以确定异常下限是当前圈定异常的一种有效途径。
地球化学异常下限值是区分背景区与异常区的基本指标,而计算异常下限值的准确性也直接关系到下一步探矿工作开展的关键。
本文分为三个部分论述土壤化探异常下限的确定。
首先介绍一些土壤化探异常下限的确定的相关概念;其次介绍各种方法,如:剖面图法、直方图解法、面积校正累积频率法、马氏距离法、单元素计算法、累积频率法、迭代法、传统统计方法、多重分形法分形、均值标准差法、含量-面积(C-A)分形方法、概率格纸图解法等);最后用一些矿床应用实例来验证及评价一些方法。
本文选取新疆西天山成矿带托逊地区1:50000土壤X荧光化探样品中Mn、Fe、Zn、As四种元素为例,使用传统统计方法、多重分形方法、85%累计频率法分别对化探数据进行处理后得出结论:传统统计方法计算出的异常范围小,且较为分散;多重分形方法对弱小异常的固定效果明显,但范围过大;85%累计频率法与传统方法所得异常下限值比较接近,但对弱小异常的识别效果相对于传统方法显著;对化探找金中背景值、异常下限的传统计算方法进行了讨论;土壤元素异常下限值的确定对环境地球化学评价具有重要意义。
传统异常下限值计算方法仅适用于元素含量数据呈正态分布的情况, 而事实上土壤元素含量的空间分布极其复杂, 很可能具有多重分形分布特征。
本文利用校正累积频率分形方法确定铜陵矿区土壤中的异常下限值为1.687 mg / kg , 并据此圈定了异常范围。
与传统方法所确定的异常下限值及相应异常区域对比, 分形方法圈定的异常区域范围更广, 更为合理、有效。
化探异常评序(精)
二、异常评序方法地球化学综合异常评序采用相同主元素类内排序方法,即将同一主元素的异常分别置于一起,依其评序指数(JOI)由大到小进行评序;遇有某一异常内有2个以上主元素者,则分别计算各自的评序指数,参与同类主元素异常的评序。
评序指数(JOI)计算方法:JOI=D×K×C式中D—主元素规模;K—异常内(除主元素外)特征组合元素平均衬度,若是两个主元素在类内评序时,其余一个主元素与其它元素一起参加特征组合衬度计算;C—修正系数。
修正系数C由以下几个方面内容组成:1、主元素异常浓度分带的外、中、内带各给1、2、3分。
2、异常组合特征:组合复杂与主元素密切相关者给3分,组合单一者给1分,介于二者之间给2分。
3、异常套合程度:与主元素异常套合程度好或有一定规律分布者给2分,较好者给1分,不好者不给分。
4、异常与已知矿产吻合程度:反映了已知矿点、矿化点者给2分;不一致者给1分。
5、异常分类类别甲1、乙1、乙2、乙3、丙类分别给4、3、2、1、0.5分。
6、异常区的地层、岩浆岩、大地构造位置等条件,对成矿有利者,按其所处条件好、较好、一般者分别给3、2、1分。
最后将各项指标得分累加,以总分的高低给定修正系数值,按>15、14~11、10~7、<7分别给C赋值1.2、1、0.8、0.6参与运算,以评序指数的高低进行异常排序。
三、异常评序结果区内所有被选综合异常均参与了评序,评序指数(JOI)顺序排列前十位的元素如下:第一为Ni,第二为W,第三为Ag,第四为Sb,第五、第六为Au,第七为Mo,第八、第九为Pb,第十位为Mo,目标矿种显示良好,成矿潜力较大。
(一)Ni异常区内Ni异常主要出现于测区北端奥陶纪地层和东南角的花岗闪长岩脉中,元素组合以NiCr为主,其中北端异常呈NW-SE向分布,和区内主体断裂走向一致,伴生AuAgCu 元素组合,具有一定的找矿意义;而东南角和花岗闪长岩体侵入有关。
进入地球化学异常特征组合者3处,成为主元素者为2处,评序结果见表4-4。
化探数据处理成图过程
化探数据处理成图的过程毕武1、2段新力1、2黄显义1、2袁小龙1、2彭仲秋1、2李永华1、21.乌鲁木齐金维图文信息科技有限公司,新疆,乌鲁木齐,8300912.新疆地矿局物化探大队计算中心,新疆,昌吉,8311000 前言GeoIPAS软件用户群不断扩大,由于各用户对系统的熟悉程度不同,对软件功能了解不够,有必要分专题将GeoIPAS处理数据及成图过程做一系统总结,下面就化探数据处理成图的过程做一总结。
1 处理步骤化探处理的成果包括:(1)参数统计表;(2)R型聚类分析-谱系图;(3)重复样三层套合方差分析或者重复样合格率计算结果;(4)点位数据图;(5)地球化学图;(6)直方图;(7)组合异常图;(8)综合异常图;(9)远景区划图;(10)单元素异常参数统计(附表册);(11)异常剖析(附图册);(12)综合异常登记卡(附表册)。
在GeoIPAS系统中,化探数据处理分为以下几个主要步骤:1.1 数据检查数值检查,坐标检查,重复样坐标检查。
1.2 分析处理重复样三层套合方差分析、重复样合格率计算、化探特征参数统计、化探背景值分析、R型聚类分析、因子分析。
1.3 数据分析数据变换;衬值、累加衬值;数据累加、累乘、比值;异常归一化。
1.4 网格化离散数据网格化、XYZ数据转网格数据1.5 成图点位数据图、彩色等量线图、直方图、组合异常图、单元素异常图、综合异常图、剖析图。
1.6 单元素异常参数统计1.7 综合异常登记卡图2 化探数据处理成图流程2 具体处理过程2.1 数据检查我们的数据处理工作从化验室提供的样品分析报告开始,项目要提供坐标和样品对应的分析数据,坐标我们一般取实际工作中的米单位,系统中默认东西向横坐标为X坐标,不加带号,南北向纵坐标为Y坐标,需要时还要提供样品对应的地质编码,我们拿到这个数据后首先进行数据检查,以确保数据中不出现写错、漏填、负数、0、>等字符,如果有这样的情况要找实验室给予纠正。
用MAPGIS作化探元素或物探参数异常图
用MAPGIS作化探元素或物探参数异常图一、改EXCEL数据格式为DET数据格式建立单元素EXCEL数据〔X、Y、单元素或物探参数符号〔不能带百分号和括号〕,此处的X、Y是数学意义上的横坐标、纵坐标,顺序不能颠倒〕—翻开单元素或物探参数EXCEL数据〔X、Y、单元素或物探参数含量值〕—另存为CSV格式〔逗号分隔〕—保存—是—是—是—翻开记事本—文件—翻开—查找围—××盘××夹—选文件(*.csv)—用记事本翻开CSV格式数据文件—加文件头NOTGRID〔非网格化数据〕—另存为—保存类型—所有文件—改文件名后缀CSV为DET格式。
另外,关于最后作出来的图形比例尺问题,如X、Y用实际坐标表示,即X为6位整数、Y为7位整数,作出来的图比例尺为1:1000,如X、Y用虚假坐标表示,即X为5位整数带1位小数、Y为6位整数带1位小数,作出来的图比例尺为1:10000二、高程点标注制图〔建点位图〕MAPGIS6.7主菜单—空间分析—DTM分析—关掉三角剖分显示窗口1—翻开DET格式文件,不能翻开GRD网格化数据〔并改目录为点位图〕—模型应用—高程点标注制图〔X、Y、标注对象〕—缺省符号—改符号〔点、颜色等〕—标注位置—对齐—右边—确认—左单击“文件〞—存数据于点数据文件—××.WT—保存类型—点文件〔*.WT〕—保存—存数据于线数据文件—××.WL —保存类型—线文件〔*.WL〕—保存—输入编辑—统改参数—保存。
关闭窗口—工作区数据已修改,建议先保存数据,先保存吗?—否。
实质就是建MAPGIS点位图。
三、等值线图MAPGIS6.7主菜单—空间分析—DTM分析—关闭三角剖分显示窗口1—翻开DET格式文件—GRD模型—离散数据网格化—网格参数设置—改起点、终点坐标—变小到适宜〔取整数〕—网格间距〔取适宜数据如0.05〕〔读者注:网格参数设置不知啥意思,建议不要动〕—网格方法〔克立格〕—文件换名—××.Grd—保存—确定。
化探工作方法简介2003版
化探工作方法简介主讲赵玉明化探工作在地、物、化、遥四大专业找矿工作中占有很重要位置。
它与其它三项找矿工作比较具有明显的直观性、有效性、经济和快捷性的特点,是在覆盖区域找矿不可缺少的重要找矿方法。
一、化探工作分类及其野外工作方法(一)按工作性质和测区范围的分类1、区域化探(战略踏勘性化探)。
目的是为了发现找矿远景区(带)、大型矿田、大中型矿床的大致区域和了解区域控矿因素(地层、构造和火山岩)。
面积数千平方千米或更大。
工作比例尺:1:10万、1:20万、1:50万。
采样密度分别是:2点/Km2、0.25-1点/Km2、0.4-0.8点/Km2。
2、地球化学普查(普查化探)。
目的是在区域化深成果基础上为了查明成矿有利地区和取得与成矿有关的地球化学特征资料等,在区域异常区内优选出一定面积的找矿靶区,进一步做化探工作。
工作面积数千平方千米或数百平方千米。
工作比例尺:1:2.5万、1:5万。
采样是采集水系沉积物质,密度为2-8点/Km23、地球化学详查或异常检查(详细化探)。
目的是在前两类化探异常区内,进一步查明异常与矿体之间的关系,或证实是不是矿化异常,为山地工程布设查证定位提供依据。
面积:0.几平方千米—几十平方千米。
比例尺:1:5千、1:1万。
采样主要是采取土壤。
密度100点----200点/Km2或>200点/Km2。
点线距为50X20~100X20m。
(二)按采样介质和野外工作方法不同进行的分类1、气体测量:寻找能产生气体的矿产。
如氡气、汞蒸气、甲烷气等2、植物测量:目前不采用。
用于沼泽、湿地、草原和森林区。
取植物样,灰化,测灰中金属含量。
3、水化学测量:目前不采用。
用于水网密集发育区,取水样,测定水中金属含量。
4、水系沉积物测量:用于区域化探和普查化探中。
比例尺:1:5万~20万之间。
测区应是水系发育(网状)或沟谷河道发育区。
采取表层以下(20cm)淤泥和粉砂为主。
采样一般粒径为-0.216mm(60目)或-0.172mm(80目)细粒物质。
化探中五个常用参数的应用
化探中五个常⽤参数的应⽤2019-08-17【摘要】地质找矿⼯作中,尤其是化探经常应⽤浓集克拉克值、变化系数、富集系数、异常下限、相对偏差五参数。
本⽂介绍了五参数的概念及含义、计算公式和⽅法、应⽤,以期地质化探⼯作者对五参数的认识⼀致,在有关地质化探的⽣产科研报告、书刊杂志中使⽤统⼀,有益于化探规范标准化。
【关键词】浓集克拉克值;变化系数;富集系数;异常下限;相对偏差我国的勘查地球化学,正处于全⾯发展时期,以全国多⽬标区域地球化学调查为主要标志;⽴⾜于地球化学填图及地球化学填图与矿产勘查⼀体化。
地质、物探、化探、遥感等综合⽅法找矿中,化探先⾏。
化探成果具有重要的指导找矿作⽤。
理解掌握浓集克拉克值、变化系数、富集系数、异常下限、相对偏差五参数,并应⽤于地质化探⼯作实践中,是取得化探成果的重要理论基础。
1 浓集克拉克值1.1 概念及含义浓集克拉克值,指在某⼀地区或某种地质体内元素含量平均值与该元素克拉克值之⽐值[1]。
⼜称浓度克拉克值[2]。
某⼀地区,指⾯积上千、⼏百、或⼏⼗平⽅千⽶的某区域、某测区、某⼯区;某种地质体,指某⼀地区内的某种地质单元、地层、岩体、岩性或某矿床、矿体、矿物等。
元素含量平均值,即元素丰度,未剔除异常含量(特⾼值或特低值)的所有样本参加计算的元素含量平均值。
元素克拉克值,即地壳的平均化学成分或元素在地壳中的丰度[3]。
元素克拉克值表[3]中的单位,先是采⽤重量百分含量(%),现地壳元素丰度表[4]中多采⽤相对含量(×10-6)。
两表中6~9⼈不同年份发表的地壳中元素丰度值,以维诺格拉多夫(1962)和泰勒(1964)、黎彤(1967)3⼈的数据较常⽤。
1.2 计算公式和⽅法1.3 应⽤1)⽤于描述某⼀地区或某种地质体内元素含量的分布特征当⼀元素的浓集克拉克值⼤于1时,则称该元素集中,其浓集克拉克值⼩于1时,则称该元素分散。
元素的集中和分散是元素迁移的结果[6]。
2)指⽰元素的地质找矿意义。
化探中五个常用参数的应用
化探中五个常用参数的应用作者:藏金生李诗言蔡新明来源:《科技视界》2013年第28期【摘要】地质找矿工作中,尤其是化探经常应用浓集克拉克值、变化系数、富集系数、异常下限、相对偏差五参数。
本文介绍了五参数的概念及含义、计算公式和方法、应用,以期地质化探工作者对五参数的认识一致,在有关地质化探的生产科研报告、书刊杂志中使用统一,有益于化探规范标准化。
【关键词】浓集克拉克值;变化系数;富集系数;异常下限;相对偏差我国的勘查地球化学,正处于全面发展时期,以全国多目标区域地球化学调查为主要标志;立足于地球化学填图及地球化学填图与矿产勘查一体化。
地质、物探、化探、遥感等综合方法找矿中,化探先行。
化探成果具有重要的指导找矿作用。
理解掌握浓集克拉克值、变化系数、富集系数、异常下限、相对偏差五参数,并应用于地质化探工作实践中,是取得化探成果的重要理论基础。
1 浓集克拉克值1.1 概念及含义浓集克拉克值,指在某一地区或某种地质体内元素含量平均值与该元素克拉克值之比值[1]。
又称浓度克拉克值[2]。
某一地区,指面积上千、几百、或几十平方千米的某区域、某测区、某工区;某种地质体,指某一地区内的某种地质单元、地层、岩体、岩性或某矿床、矿体、矿物等。
元素含量平均值,即元素丰度,未剔除异常含量(特高值或特低值)的所有样本参加计算的元素含量平均值。
元素克拉克值,即地壳的平均化学成分或元素在地壳中的丰度[3]。
元素克拉克值表[3]中的单位,先是采用重量百分含量(%),现地壳元素丰度表[4]中多采用相对含量(×10-6)。
两表中6~9人不同年份发表的地壳中元素丰度值,以维诺格拉多夫(1962)和泰勒(1964)、黎彤(1967)3人的数据较常用。
1.2 计算公式和方法1.3 应用1)用于描述某一地区或某种地质体内元素含量的分布特征当一元素的浓集克拉克值大于1时,则称该元素集中,其浓集克拉克值小于1时,则称该元素分散。
元素的集中和分散是元素迁移的结果[6]。
刘红杰关于异常下限的几种计算方法
地球化学元素含量的异常确定是勘查地球化学中最重要的工作之一,但迄今为止还没有找到一个完全令人满意的具有科学依据的方法。
长期以来,人们主要是使用经典的统计学方法,以样品数据呈正态分布为假设前提,通过计算数据的统计学参数(如均值、标准离差等)对异常进行筛选和评价。
一般是以平均值(X)与2倍(也有为1.5倍或3倍)的标准离差(δ)之和作为地球化学的异常下限值。
该方法仅适用于地球化学数据呈正态分布的情况,但实际上对于元素的地球化学分布而言正态分布并不是唯一的一种分布,人们已经发现许多元素,特别是微量元素并不遵循正态分布,而是呈明显的正向偏斜或表现为一种幂型的拖尾分布。
其他几种用来筛选和评价地球化学异常的方法,如移动平均法、趋势面法、克里格法以及概率格纸法等,除了概率格纸法仍是基于正态分布这一观点外,其他的几种方法虽然注意到了元素含量分布的空间信息,但都是以地球化学含量数据在空间上呈连续变化,且是一个光滑的连续曲面这一假设为基础建立的。
事实上,地球化学元素含量的空间分布是极其复杂、十分粗糙而并非处处可微的。
正如李长江等(1995)研究揭示的地球化学景观可能是一个具有低维(D=2.9)吸引子的混沌系统,是分形。
考虑到方法的实用性、有效性、易操作,通过几种方法在工作区的试验对比,叠代法确定的背景值及异常下限较低,更有利于突出弱异常。
因此,工作区背景值和异常下限的确定选用叠代法。
叠代法处理的步骤:①计算全区各元素原始数据的均值(X1)和标准偏差(S1);②按X1+3S1的条件剔除一批高值后获得一个新数据集,再计算此数据集的均值(X2)和标准偏差(S2);③重复第二步,直至无特高值点存在,求出最终数据集的均值(X)和标准偏差(S),则X做为背景值C0,X+nS(n根据情况选1.5或2,3)做为异常下限Ca。
化探异常圈定判别查证方法
异常圈定与判别
? “高、大、全”异常评价易,弱小异常评价难!浅部矿化异常评价易, 深部矿(包括隐伏矿和盲矿)异常评价难!
多元素规模
Σ NAP
序数
1.49
10
3.28
8
0.4
18
4.13
7
4.28
6
0.73
14
1.6
9
0.89
13
19.98
2
8.22
5
86.33
1
1.04
12
0.44
16
13.03
3
1.12
11
11.91
4
0.26பைடு நூலகம்
21
0.24
22
0.4
17
0.61
15
0.31
20
0.35
19
某地局部异常特征排序表
浓度分带
异常圈定与判别
? 元素异常下限确定 根据测区地质和景观特征,选择使用整个测区或划分子区来确定异 常下限。
? 对勘查区面积不大,地质情况比较简单,元素呈明显单峰或有一个 单一的背景全域和一个异常全域分布,就可以在区内剔除高、低值 含量(X±nS)根据工作区情况,n值一般为1.5-3,计算出一个统一 的背景平均值(X0)和异常下限(T);累积频率85%或95%对应的含 量,目估法经验判断。无论何种方法确定仅作为参考值,以能否客 观反映本测区(或子区)矿产和矿化分布特征为依据,作适当修正。
化探中五个常用参数的应用
【 K e y w o r d s ] C i a r k e v a l u e ; V a r i a t i o n c o e f f i c i e n t ; E n i r c h m e n t c o e f i f c i e n t ; T h r e s h o l d ; R e l a t i v e d e v i a t i o n
科技・ 探索・ 争鸣
S c 科 i e n c e & 技 T e c h 视 n o l o g y 界 V i s i o n
化探中五个常用参数的应用
藏 金生 李诗 言 蔡 新 明 ( 河 南省地质 矿产 勘查 开发局 第 三地质 矿产调 查院 , 河南 信阳 4 6 4 0 0 0 )
a p p l i e d d u in r g t h e p r o c e s s o f o r e - p r o s p e c t i n g . e s p e c i a l l y d u r i n g g e o c h e mi c a l e x p l o r a t i o n . Th i s t h e s i s i n t r o d u c e s t h e c o n c e p t i o n a n d d e in f i t i o n. d e s i g n
f o r mu l a s a n d me t h o d s a n d a p p l i c a t i o n o f t h e s e i f v e p ra a me t e s,S r O a s t o u n i f y g e o c h e mi c l a e x p l o r e s’c r o g n i t i o n o f t h e s e iv f e p a r a me t e r s ,a n d t h e i r u n i f o m r u s e i n t h e s c i e n t i i f c r e s e a r c h r e p o r t a n d ma g a z i n e s r e l a t e d t o g e o c h e mi c a l e x p l o r a t i o n, wh i c h b e n e it f s t o n o r ma l i z e t h e s t a n d a r d o f g e o c h e mi c a l
矿产勘查中化探异常下限的计算方法
第29卷 第1期物探化探计算技术 2007年1月基金项目:中国地质调查局科研项目(DK9902122)收稿日期:2005-12-28文章编号:1001—1749(2007)01—0054—04矿产勘查中化探异常下限的多重分形计算方法孙忠军(中国地质科学院物化探研究所,河北廊坊065000)摘 要:化探异常下限计算方法在矿产勘查和资源预测中非常重要。
传统的计算方法存在依据不充分的问题。
这里在研究地球化学场数据分布的基础上,提出了多重含量~频数分形计算方法。
研究表明,多重分形异常下限计算方法较之传统的方法具有理论基础扎实、适用性较广的特点。
关键词:异常下限;分形分布;多重含量~频数分形中图分类号:P 632 文献标识码:A0 前言化探异常下限的计算在矿产勘查和预测中占有十分重要的地位。
随着人们对地球化学数据分布特征的认识,已发展了多种计算方法。
传统的异常下限计算方法的依据是,假定地球化学场数据分布是一种光滑(处处可微)的连续曲面。
先后开发出移动平均法、趋势面法、克力格法、概率格纸法和均值加标准离差法等[1~5]。
国内、外近十几年的研究成果表明,区域地球化学场十分粗噪而并非处处可微,数据服从分形分布[6~12]。
研究区域地球化学场异常下限的分形方法已成为一种新的发展趋势。
本文作者研究了藏麻西孔和扎西尕日矿区地球化学场的数据分布特征,提出了多重含量~频数异常下限计算方法。
1 试验区地质特征试验区位于青海省风火山成矿带的藏麻西孔和扎西尕日矿区。
藏麻西孔是中型铅锌矿。
矿区地质构造背景是近东西向的背斜构造,主体断裂为近东西向的逆断层,近南北向的平移断层切割了背斜。
背斜核部出露的是白垩系上统风火山群下岩组中、下岩段的地层,背斜翼部由风火山群上岩段组成。
风火山群上岩段(K 2fn a -3)岩性为紫红色长石石英砂岩夹岩屑石英砂岩,中岩段(K 2fn a -2)为浅色岩屑石英砂岩夹长石石英砂岩,下岩段(K 2fn a -1)是紫红色长石石英粉砂夹岩屑石英砂岩。
化探常用表格全解
地球化学普查水系沉积物测量记录项目名称: 图幅号: 采样日期:年月日记录: 采样审核:A 样品号:N7。
图幅名拼音代码+采样大格编号+小格代码+小格样号,如:MP234B1。
该样品号中:MP-茅坪幅代码;234-大格号;B-小格号;1,B小格第一个样号)。
B 原始样号:被重复采样的样品号。
C样袋号:野外实地编号。
D 图幅代号:N10。
1:50000地形图图幅号,如H49E007008E 横坐标:N8。
统一确定为高斯6度带,记录带号+横坐标精确到m。
如20428303F 纵坐标:N7。
高斯6度带精确到m。
女口3395158G 海拔高程:N4。
采样点高程坐标,以米为单位。
从地形图等高线或通过GPS直接读取。
H 水系级别:C1。
记录:1、一级水系;2、二级水系;3、三级水系。
I 采样部位:C1。
采样点位于水系的位置,用代码表示:1:河底;2:水线附近;3:河漫滩上;4:水塘入口处J 样品组分:C3。
记录3位数:分别代表样品中粗砂(第1位)、细砂(第2位)和淤泥及有机物(第3位)含量。
此三项为样品的沉积物组分,以编码方式分级填写,分为:0:无; 1:少量(<30%); 2 :中量(30〜70%); 3:大量(>70 %),三者之和不能超过100 %。
K 样品颜色:C2。
1、灰黑色;2、灰色;3、褐色;4、灰黄色;5、红色;6、砖红色;7、灰绿色。
L 地质时代:C4。
记录所控汇水域内地质时代。
记录地质时代符号。
沉积地层按出露情况适当并层;侵入岩记录主要侵入期。
M 岩石类型:C4。
填写该点所控制汇水面积内占优势的基岩类型,参见“区域地球化学勘查规范”附录B表B2。
N 矿化蚀变:C1。
记录矿化蚀变程度。
0、无;1、弱;2、中等;3、强烈。
O 地貌类型:C1。
1、平原■准平原;2、低山-丘陵;3、山地-峡谷;4、高山-深谷;5、高原;6、高寒山地;7、盆地;&沼泽洼地;9、岩溶石山。
P 植被:C1。
0,无;1,稀疏,浅薄,覆盖度<1/3 ;2,中等,覆盖度在1/3〜2/3间;3,茂密,浓厚,覆盖度>2/3。
利用MAPGIS处理化探数据文档
地球化学探矿工作方法及注意事项异常下限计算公式:都采用Excle表格计算(原始数据参照表三)平均值:=A VERAGE(C2:C36)方差:STDEV(C2:C36)异常下限:平均值+K值*方差计算过程中自动剔除大于平均值加3倍离差的高值和小于平均值减3倍离差的低值。
K值自选(1-3)剔除次数的多少取决于没有大于和小于这些数据。
表一表二表三等值线图生成方法步骤:利用Excel先做一个三角部分文件,格式如图:A1:要打上notgrid.A2为生成点文件的X坐标 B2为Y坐标 C2为高程代表化学分析结果。
文件需要另存格式如图:点击保存,关闭表格。
打开MAPGIS6.7空间分析-DTM分析-点击文件打开三角部分文件。
如图:文件打开后在菜单栏中找到Grd模型打开里面有一栏离散数据网格化打开。
如图:网格化方法用:Kring泛克立格法网格化,文件换名保存文件。
“选择打开”,选三阶漂移如图:其它内容不用改动,改完后按确定。
“搜索打开”搜索类型改为八方向,其它内容不用改动,改完后按确定。
如图:全部改完后按确定。
再一次在Grd模型里找到平面等值线图绘制点击后打开离散数据网格化后保存的文件,如图:把文件打开如图:设置:等值线套区、绘制色阶、保留边界线全部选上。
光滑:等值线光滑处理选上,光滑度选择高程度。
等值线定层一般定四层,第一层为背景值,第二层为异常下限值,第三层为中带值(异常下限值的2倍),第四层为内带值(异常下限值的4倍)。
注记参数都改为YES。
区、线参数参照图例规范设定。
制图幅面:改为数据投影变换如图:原数据投影参数与结果数据投影参数都要改变,设置结果与投影点位图一致。
全部设置完后点击确定。
生成图像如图:生成图像后保存点、线、区。
保存好以后再一次打开MAPGIS6.7图形处理-输入编辑打开刚才已保存好的点、线、区文件进行图式修饰。
点位图制作方法步骤:利用Excel表格把所有化探数据先处理一下如图:注:ID代表总序号、YH代表样号、Y代表X轴坐标、X代表Y轴坐标、Au、Ag 代表元素化学分析结果。
化探异常圈定判别查证方法
概况
? 野外调查数据(记录卡信息) 野外工作中观察到和记录的数据,通过系统的整理和归类,使野外 现象能完整的在数据库体现出来。
? 地质数据 采样点岩石类型、矿化蚀变信息,成矿构造信息,成矿类型信息等。
? 综合数据
各种数据组合成工作数据文件,供数据处理和成图。
概况
? 基本化探参数及图解
? 单变量参数 地球化学区和子区、重要地质单元样品数(n)、算术平均值(X), 几何平均值(G),中位数(Me),众数(Mo)、均方差(S0),极差等;
? 分析质量评估
分析质量不合格,所得出的结论也是错误的,所以未做评估的数据无 可靠性可言。通过分析数据质量的评估,说明哪些元素数据满足化探 要求,哪些数据可作为参考或在哪些方面可以利用。比如报出率小于 50%显然不能满足编制地球化学图的要求,但仍可用于异常圈定。
概况
? 建立数据文件
化探数据量大、指标多,令人眼花缭乱,必须对所有取得的数据进 行认真严格整理。 ? 采样点坐标数据 主要对实际采样点位与分析测试样品号一一对应起来! ? 分析数据 对实验室提供的数据进行分析。按照样品性质进行分类整理,在 Excel电子表格中很方便地筛选出来所需的样品测试数据。
异常圈定与判别
某地采用剔除法确定有异常下线
异常圈定与判别
? 当勘查区内元素频率分布出现一 个单一的背景全域和一个异常全 域交迭出现双峰,或未出现双峰, 但频率分布曲线呈正向倾斜时, 一般可以用众值(M0)代替平均 值(X);
? 对于元素分配不服从正态或对数 正态分布,而是呈双峰或多峰态, 采用多重分形方法确定异常下限。
异常圈定与判别
多重分形确定异常下限
异常圈定与判别
? 特殊情况地球化学异常下限确定 ? 地表弱矿化地区,需要降低异常下限; ? 针对相邻高含量点相伴出现(圈定低缓异常)时,可根据高含量点数降
化探单元素异常统计内容参数公式
化探单元素异常统计内容1、异常ID ID2、样品个数 N3、异常面积 S4、样品最大值 Max5、样品最小值 Min6、异常下限 T7、算术平均值 nXiX n∑=18、几何平均值 ∑=ng Xi n X 1log 19、标准离差 1)X X(n1i 2__0--∑=n S i=10、异常衬度 TX A c =11、异常规模 ()T X S A d -⨯= 12、异常NAP 值 S A NAP c ⨯=浓集克拉克值(C)计算公式: 某元素的克拉克值XC =变化系数(Cv)计算公式: XS Cv =致矿系数(Z)计算公式:Z =Cv(全区)+10×Cv(剔高值后)+100×高值比例+C ,高值是大于3倍的标准离差。
化探背景分析 中位数:505050)50(f F H X Me -+=偏度:24)(1)(13231nX X f n X X f nR i i ii •⎪⎪⎭⎫ ⎝⎛--=∑∑ 峰度:963)(1)(14242n X X f n X X f nR i i i i •⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-⎪⎪⎭⎫ ⎝⎛--=∑∑正态检验: ∑•-=ifx F x F |)()(|max 1λXi :组中值或含量值;f i :Xi 所对应的频数;H :组距;X50:包括累计频率50%在内的所在组的组下限;F 50:累计频率50%所在组之前的累计频率;f 50:包括累计频率50%所在组的组频率;F(X):为经验累计频率;F 1(X):为理论累计频率。
R 型聚类分析iiij ij S X X X -='其中:∑==nj ij X n X 11;1)(12--=∑=n X (XS nj i iji∑∑∑===-⋅---=⋅=ni ni k kij jik kini j jikkjj jk jk X (X X (X X (X X (XS S S r 11221)())()(式中:r kj 为第j 个变量和第k 个变量的相关系数;X ji 为第j 个变量第i 个样品的观测值; X j 与X k 为第j 个和第k 个变量的平均值。
单元素异常特征统计方法
单元素异常特征统计方法四川省地质矿产勘查开发局202地质队 张应先单元素异常图绘制出来以后,需要统计单元素的异常特征参数,统计内容包括包括异常点数,异常面积、异常均值、几何平均值、异常极大值、异常均方差、异常衬度、异常规模等各项特征参数,传统的用mapgis进行统计过程比较繁琐,工作量较大,本文介绍由元素异常线来做异常特征参数的统计方法。
1、单元素异常图的绘制统计并计算出元素异常下限,并将原始数据进行网格化处理(注意选择正确的网格间距和搜索半径范围以及网格化方法),利用mapgis或者GeoiIPAS,surfer等软件绘制单元素异常图。
(图一) 推荐使用GeoIPAS,GeoMDSS,GeoExpl专业软件。
图一 金元素异常图2、元素异常统计方法在mapgi中调入异常线并对其做相关处理,具体包括:删除小的异常线,闭合不封闭的异常线等,然后对异常线进行异常编号。
修改异常线的时候最好放置在原始数据上参照修改。
(图二)修改好以后,就可以利用金维(GeoIPAS)进行单元素异常特征统计了。
(需要注意的是网格化方法和网格间距以及搜索半径是决定异常范围和形态的重要因素,一定要选中好正确的网格间距、搜索半径和网格化方法,并不是只要能把等值线图圈出来就行那么简单,圈出来的等值线需要和原始数据进行对比,多次实验。
)打开软件,并进入单元素异常特征参数统计菜单,然后装入已经修改好的异常线文件(1.64版本需要将线文件转换为明码文件,然后再装入)。
数据类型可以选中原始数据和网格化数据,文中以原始数据为例。
装入异常线文件和原始数据以后,正确设置原始数据和异常线文件的投影参数。
设置统计结果保存路径。
对应选择相关属性项目:X、Y,统计元素,填入统计元素异常下限值。
及单位。
2.8版本可以选择需要统计的项目。
设置原始数据投影参数:设置异常线文件投影参数:各项参数设置完成以后点击确定开始统计。
软件即自动统计出各项参数数据特征。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
化探单元素异常统计内容1、异常ID ID2、样品个数 N3、异常面积 S4、样品最大值 Max5、样品最小值 Min6、异常下限 T7、算术平均值 nXiX n∑=18、几何平均值 ∑=ng Xi n X 1log 19、标准离差 1)X X(n1i 2__0--∑=n S i=10、异常衬度 TX A c =11、异常规模 ()T X S A d -⨯= 12、异常NAP 值 S A NAP c ⨯=浓集克拉克值(C)计算公式: 某元素的克拉克值XC =变化系数(Cv)计算公式: XS Cv =致矿系数(Z)计算公式:Z =Cv(全区)+10×Cv(剔高值后)+100×高值比例+C ,高值是大于3倍的标准离差。
化探背景分析 中位数:505050)50(f F H X Me -+=偏度:24)(1)(13231nX X f n X X f nR i i ii ∙⎪⎪⎭⎫ ⎝⎛--=∑∑ 峰度:963)(1)(14242n X X f n X X f nR i i i i ∙⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-⎪⎪⎭⎫ ⎝⎛--=∑∑正态检验: ∑∙-=ifx F x F |)()(|max 1λXi :组中值或含量值;f i :Xi 所对应的频数;H :组距;X50:包括累计频率50%在内的所在组的组下限;F 50:累计频率50%所在组之前的累计频率;f 50:包括累计频率50%所在组的组频率;F(X):为经验累计频率;F 1(X):为理论累计频率。
R 型聚类分析iiij ij S X X X -='其中:∑==nj ij X n X 11;1)(12--=∑=n X (XS nj i iji∑∑∑===-⋅---=⋅=ni ni k kij jik kini j jikkjj jk jk X (X X (X X (X X (XS S S r 11221)())()(式中:r kj 为第j 个变量和第k 个变量的相关系数;X ji 为第j 个变量第i 个样品的观测值; X j 与X k 为第j 个和第k 个变量的平均值。
Q 型聚类分析(min)(max)(min)'i i i ij ij X X X X X --=(i=1,2,…,P ;j=1,2,…,n)其中:P ,n 分别为变量数和样品数;Xi(max)及Xi(min)分别为数据中第i 个指标的极大值与极小值∑=-=Pi ik ijjk P X (XD 12/)( (j ,k=1,2,…,n ;j ≠k)式中:Djk 为第j 个样品与第k 个样品的距离系数;Xij 为第i 个变量第j 个样品的观测值。
因子分析数学原理:设有一批含p 个变量,n 个样品的观测数据,如果其变量为X1、X2、……、Xp ,它们的综合变量记为F1、F2、……、Fm(m ≤p),其数学表达式为⎝⎛+++=+++=+++=p mp m m mpp p p Xa X a X a F X a X a X a F X a X a X a F ........................22112222121212121111 要求1 (2)1212211=+++p a a a (k=1,2,……,m )系数{}ija 由下列原则决定:1、 Fi 与Fj (ij ,i 、j=1,2,……,m )互相无关2、 F1是X1、X2、……,Xp 的一切线性组合中方差最大;F2是与F1不相关的X1、X2、……、Xp 所有线性组合中方差最大的; ……;Fm 是与F1、F2、……、Fm-1都不相关的X1、X2、……、Xp 所有线性组合中方差最大。
回归系数:()222XX Y X XY XX X X Y XY b -⋅-=⋅--=∑∑∑∑标准差:()22XX X -=σ()22XY Y -=σ相关系数: YX YX XY r σσ⋅⋅-=Y X b r σσ⨯= XY r b σσ⨯=回归分析1、多元线性回归分析研究某一变量与多个变量之间的线性关系(某两变量之间的非线性关系有时也可以转化为线性关系)。
这是以大量收集到的观测数据为基础,找出相关变量之间的内部规律性,以定量形式建立一个因变量与另一个自变量(或另几个自变量)之间关系的数学表达式,从而可以根据一个或几个变量的观测值来预测(预报)另一个变量的估计值,并能从多个指标变量中找出对所研究的问题起重要作用的某些指标。
2、多元逐步回归分析是以大量收集到的观测数据为基础,建立某一个变量与另一个变量(或几个变量)之间关系的数学表达式。
逐步回归是在多元回归基础是派生出来的一种算法,它能从众多的变量(或预先尽可能多地考虑一些变量)中自动挑选重要变量(指标或因子),并确定其数学表达式的一种统计方法。
3、正交化回归分析与线性回归不同,它不一定建立随机变量与全体变量之间的关系,而是一种具有挑选因子的方案。
线性回归和逐步回归仅考虑因变量(一个或多个)与自变量(多外)之间的关系,而不考虑自变量本身之间的关系,因而信息有重复的可能,而降低某一独立因素的相对权系数。
这种重复也因具体问题的不同,对回归效果就有一定的影响。
正交化回归不仅考虑自变量与因变量之间的关系,而且还要考虑自变量本身之间的相互影响。
即是将自变量因子进行正交化,排除自变量之间的相互影响,得到一组新的正交化了的因子,再建立新因子与因变量之间的回归方程,并还原到非正交化、非正规化的回归方程,作为建立预报模型之用。
对应分析对应分析是在R-型与Q-型因子分析基础上发展起来和一种多元统计方法,它把R-型和Q-型因子分析结合起来,综合考虑变量之间、样品之间及变量与样品之间的关系。
它揭示了R-型与Q-型分析之间的两重性,以较少重要的几个公共因子的综合指标去研究对象在成因上或空间上的联系,应用载荷平面投影图进行地质解释与推断,在地质学的应用可以包括:对矿床成因的解释;成岩(成矿)的物质来源、作用方式、作用因素;含矿岩体的预测、评价;地球化学的研究等。
根据R-型与Q-型分析具有的对偶性,由R-型分析结果很容易地得Q-型分析结果,它可以克服由于样品容量大而对Q-型因子分析所带来计算上的困难,把变量和样品同时反映到同一坐标轴(因子轴)的一张图形上,更便于地质解释与推断。
相关分析1、典型相关分析是研究两组地质指标(或变量)之间相关关系的一种统计方法。
它揭示了两个因素“集团”之间的内部联系,而两个因素“集团”的内容和变量数目又可以不同,这就决定了它在解决地质问题上的许多特点,在地质上用来研究二组地质特征之间的关系。
用以研究地质成因和地质体对比。
如研究某种矿物成分与其围岩成分的关系;研究古生物群同古地理环境之间依存关系;两条区域剖面的化学成分、岩石成分的对比等等。
典型相关分析的目的,是在两组众多地质变量中分别寻找若干对(每对为若干个变数的线性组合)有代表性、且具有相关关系的综合指标,通过研究两组综合指标间的关系来反映两组地质变量间的主要作用因素和作用形式。
2、秩相关分析是研究一组观测值与另外一组观测序列之间的相关关系,在矿床统计预测中,研究各种控矿地质因素和找矿标志与矿化的相关关系,查明这样的因素与标志最有利的数值范围,即查明找矿有利的统计标志。
两个观测序列的秩相关系数:)1(6122--=∑n n d P其中d 为对比序列的序差;n 为对比序列的等级数。
3、相关频数比值法:筛选因子(变量)的目的在于从大量的预报因子中选出与预报量相关较好,而因子之间相关较差,即因子独立性强的若干较优因子,由它们组成数学模型使预报效果达到最佳。
计算公式化:第i 个因子的相关频数比iii n n m '=n i :变量Xi 报对的频数,指预报量Y=1时Xi=1,Y=0时Xi=0n ’i :因子间相关频数,指Xi 因子在各样本个体中的报错时,相应样本个体其它因子重复报错的总频数。
特征分析特征分析是一种对定性描述资料的多元统计方法,其理论简单,通过对已知模型区内各变量间相互关系的定量考察,来确定每个变量的权系数,用于预测的原则是类比法。
可综合处理各种数据,尤其适用于多类型定性描述资料的综合处理。
在地质找矿过程中常用于预测与评价找矿远景区,亦可用来挑选控矿的重要变量。
方法原理:从乘积矩阵Z ’Z (Z 为原始数据矩阵)和匹配矩阵T 出发,采用主分量分析的雅可比法(Jacobi )寻找最大特征值(λ1)及其相应的一组特征向量(b ij ),来确定各标志的权系数,建立模型与变量线性关系,同时得到预测区的综合指标——关联值,提供定量评价、预测。
点聚图0-1法%100⨯-=NN 判错的样品数区分率隶属度和贴近度1、根据在集合上的隶属函数,按隶属原则识别对象,判定其属于哪一个类型;2、根据各弗晰两两间的贴近度,按择近原则,确定哪两个弗晰集最贴近。
1、隶属度设有n 种类型,m 种指标,则第i 种类型在第j 种指标上的隶属函数为:Xb a b a X a a X a a X b a b a X b a X b a X A ij ij ijijijijijijij ijijijijij ijij ij ≤++<<≤≤<<--≤⎝⎛⎪⎪⎭⎫ ⎝⎛--⎪⎪⎭⎫ ⎝⎛--=)2()2()2()2()1()1()1()1(2)2(2)1(01110其中a ij (1)和a ij (2)分别是第i 类元素第j 种指标的最小值和最大值;b ij 2=2σij ,而σij 2是第i 类元素第j 种指标上的方差,给定一具体对象X *,设它的m 个指标为X 1*、X 2*、…、X m *,令m j X A S j ij i ≤≤=1*)(min 又若n i S S i i ≤≤=1)max(0则认为此对象属第i 0类A i0 。
2、贴近度已知有n 种类型(A 1,A 2,…,A n ),它们都有m 种指标,均为正态型弗晰变量,相应的参数分别为(a ij (1)、a ij (2)、b ij )(i=1,2,…,n ;j=1,2,…,m ),其中a ij (1)=min (X ij ),a ij (2)=max (X ij ),b ij 2=2σij 2,σij 2是X ij 的方差,待判断对象B 的m 个指标分别是具有参数(a j ,b j )(j=1,2,…,m )的正态型弗晰变量,则B 与各类型的贴近度为:jij j ijij j ij j ijijj ij ijj ij j ijij j ij j ij j ij j ij j ij jij a b b a b b a a a a a a a a b b a b b a a b b a a b b a a B A ≤++++<<≤≤<<+-+-≤⎝⎛⎪⎪⎭⎫ ⎝⎛+--⎪⎪⎭⎫ ⎝⎛+--=)()()()(021112110),()2()2()2()2()1()1()1()1(2)2(2)1(又记 ),min(B A S ij i = 若有 n i S S i i ≤≤=1)max(0则按贴近原则,可以认为此B 与A i0 最贴近。