化探单元素异常统计内容参数公式
地质学中一些公式

地学中常用公式一、平均品位的计算公式: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—保存—确定。
- 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 最贴近。