CN106126882B - 基于肯德尔和谐系数排序的元素地球化学数据优选方法 - Google Patents

基于肯德尔和谐系数排序的元素地球化学数据优选方法 Download PDF

Info

Publication number
CN106126882B
CN106126882B CN201610425349.9A CN201610425349A CN106126882B CN 106126882 B CN106126882 B CN 106126882B CN 201610425349 A CN201610425349 A CN 201610425349A CN 106126882 B CN106126882 B CN 106126882B
Authority
CN
China
Prior art keywords
mine
concordance
geochemistry
coefficient
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610425349.9A
Other languages
English (en)
Other versions
CN106126882A (zh
Inventor
王新华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN201610425349.9A priority Critical patent/CN106126882B/zh
Publication of CN106126882A publication Critical patent/CN106126882A/zh
Application granted granted Critical
Publication of CN106126882B publication Critical patent/CN106126882B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了属于找矿勘探技术领域的一种基于肯德尔和谐系数排序的勘查区元素地球化学数据优选方法;该方法是收集或采集元素地球化学数据;将地球化学数据插值处理后投影到规则网格矩阵内并网格化;计算各元素与成矿之间的肯德尔和谐系数;按肯德尔和谐系数高低进行排序,编制优选后元素地球化学异常图;然后对勘查区元素地球化学数据进行评价并排序,提取出与成矿相关性强的元素数据,作为找矿有利元素,应用于勘查区的找矿预测实践,通过对地球化学数据的优选,提高找矿预测效率与质量。

Description

基于肯德尔和谐系数排序的元素地球化学数据优选方法
技术领域
本发明属于找矿勘探技术领域,特别涉及一种基于肯德尔和谐系数排序的元素地球化学数据优选方法
背景技术
元素地球化学方法是金矿地球化学找矿的一种重要的、有效的方法,该方法通过区域成矿元素和伴生元素的含量异常对金矿成矿潜力进行分析。在找矿勘探工作的元素地球化学数据选择工作中,以往的方法是根据研究者/工作者的经验对区域内找矿有利元素进行选择,主观性较强,同时因为不同工作者的经验与知识背景不同往往造成找矿实践中认识分歧。以往并没有一种明确有效的评价勘查区成矿有利元素的指标,也没有明确的依据指标成矿有利元素优选方法流程,直接影响矿产勘查评价效果,因此亟需一种新的勘查区找矿有利地球化学元素优选方法。
发明内容
本发明的目的是提供一种基于肯德尔和谐系数排序的勘查区元素地球化学数据优选方法,其特征在于,使用肯德尔和谐系数对勘查区元素地球化学数据进行评价并排序,提取出与成矿相关性强的元素数据,作为找矿有利元素,应用于勘查区的找矿预测实践,通过对地球化学数据的优选,提高找矿预测效率与质量;包括以下步骤:
1)收集或采集元素地球化学数据;
2)将地球化学数据插值处理后投影到规则网格矩阵内;
3)矿点投影变换并网格化;
4)计算各元素与成矿之间的肯德尔和谐系数;
5)将各元素按照与成矿的肯德尔和谐系数高低进行排序,编制优选后元素地球化学异常图并用于找矿预测。
所述步骤1)采集元素地球化学数据是土壤、水系、水文、生物或气体地球化学采样数据。
所述步骤2)是将地球化学数据进行插值处理并投影到100米x100米的网格内,网格矩阵需覆盖整个调查区。
所述步骤3)是将矿点坐标进行投影变换并投影到上述网格矩阵,含矿点标记为1,不含矿点标记为0。
所述步骤4)计算各元素与成矿之间的肯德尔和谐系数是针对不同矿产种类,根据肯德尔和谐系数对勘查区内地球化学元素进行排序,并优选出排序前百分之二十到百分之三十的元素类别作为找矿有利元素和指示元素。
所述步骤5)为选择前六位的W、Hg、Au、Sb、Sn、Ag,编制优选后的元素地球化学异常图;并用于找矿预测。
本发明的有益效果是使用肯德尔和谐系数作为元素地球化学数据评价标准的优势是不要求数据绝对服从高斯分布,在边界处性能更好,鲁棒性更强。为我国的找矿勘探工作提供了涵盖面广、有效性高、实用性强、准确性好的方法。
附图说明
图1为找矿勘探操作流程示意图。
具体实施方式
本发明提供一种基于肯德尔和谐系数排序的勘查区元素地球化学数据优选方法,下面结合附图和实施例对本发明予以进一步说明。
如图1所示的找矿勘探操作流程示意图。
包括以下步骤:
1)收集或采集元素地球化学数据;收集调查区1:200000-1:10000比例尺元素地球化学数据,或以1:200000-1:50000比例尺在调查区布置测网,采集地球化学样品,得到分析数据;数据可以是土壤、水系、水文、生物或气体地球化学采样数据。
2)将地球化学数据插值处理后投影到规则网格矩阵内;将地球化学数据进行插值处理并投影到100米x100米(或根据比例尺选择其他网格度)的网格内,网格矩阵需覆盖整个调查区。
3)矿点投影变换并网格化;将矿点坐标进行投影变换并投影到上述网格矩阵,含矿点标记为1,不含矿点标记为0.
4)计算各元素与成矿之间的肯德尔和谐系数;根据肯德尔和谐系数对勘查区内地球化学元素进行排序,并优选出排序前百分之二十到百分之三十的元素类别,作为找矿有利元素(此专利可针对不同矿产种类,选取相应的肯德尔和谐系数较高的元素作为该矿种找矿有利元素/指示元素)。
5)按照肯德尔和谐系数高低进行排序,此例中则选择前六位的W、Hg、Au、Sb、Sn、Ag,编制优选后的元素地球化学异常图;并用于找矿预测,如表1、表2所示。
表1 某勘查区18中地球化学元素与成矿肯德尔和谐系数表
表2 根据肯德尔和谐系数对各单元素进行排序
实施例
计算各元素与成矿之间的肯德尔和谐系数W的计算公式:
式中:为计算过程中的一个系数;
m为第i个评价者的评定结果中有重复的个数;
pij为第i个评价者的评定结果中第j个重复等级的相同等级数;
N为被评分的对象数;k为评分者人数或评分所依据的标准数;
S为每个被评对象所评等级之和与所有这些和的平均数的离差平方和,当评分者意见完全一致时,s取得最大值,可见和谐系数是实际求得的s与其最大可能取值的比值,故0≤w≤1。
这里给出的是肯德尔和谐系数和相关系数的一种基本的计算方法,根据不同应用情况也可能使用特定的经过更改的公式进行计算,根据数据特征此处肯德尔系数亦可以更换成皮尔森相关系数或斯皮尔斯曼相关系数或其他表示相关性的量度标准。

Claims (1)

1.一种基于肯德尔和谐系数排序的勘查区元素地球化学数据优选方法,其特征在于,使用肯德尔和谐系数对勘查区元素地球化学数据进行评价并排序,提取出与成矿相关性强的元素数据,作为找矿有利元素,应用于勘查区的找矿预测实践,通过对地球化学数据的优选,提高找矿预测效率与质量;
包括以下步骤:
1)收集或采集元素地球化学数据;
2)将地球化学数据插值处理后投影到规则网格矩阵内;
3)矿点投影变换并网格化是将矿点坐标进行投影变换并投影到上述网格矩阵,含矿点标记为1,不含矿点标记为0;
4)计算各元素与成矿之间的肯德尔和谐系数是针对不同矿产种类,根据肯德尔和谐系数对勘查区内地球化学元素进行排序,并优选出排序前百分之二十到百分之三十的元素类别作为找矿有利元素和指示元素;其中计算各元素与成矿之间的肯德尔和谐系数W的计算公式:
式中:为计算过程中的一个系数;
m为第i个评价者的评定结果中有重复的个数;
pij为第i个评价者的评定结果中第j个重复等级的相同等级数;
n 为被评分的对象数;k为评分者人数或评分所依据的标准数;
s 为每个被评对象所评等级之和与所有这些和的平均数的离差平方和,当评分者意见完全一致时,s取得最大值,可见和谐系数是实际求得的s与其最大可能取值的比值,故0≤w≤1;
5)将各元素按照与成矿的肯德尔和谐系数高低进行排序,编制优选后元素地球化学异常图并用于找矿预测;具体选择前六位的W、Hg、Au、Sb、Sn、Ag,编制优选后的元素地球化学异常图;并用于找矿预测。
CN201610425349.9A 2016-06-15 2016-06-15 基于肯德尔和谐系数排序的元素地球化学数据优选方法 Expired - Fee Related CN106126882B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610425349.9A CN106126882B (zh) 2016-06-15 2016-06-15 基于肯德尔和谐系数排序的元素地球化学数据优选方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610425349.9A CN106126882B (zh) 2016-06-15 2016-06-15 基于肯德尔和谐系数排序的元素地球化学数据优选方法

Publications (2)

Publication Number Publication Date
CN106126882A CN106126882A (zh) 2016-11-16
CN106126882B true CN106126882B (zh) 2019-01-29

Family

ID=57469650

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610425349.9A Expired - Fee Related CN106126882B (zh) 2016-06-15 2016-06-15 基于肯德尔和谐系数排序的元素地球化学数据优选方法

Country Status (1)

Country Link
CN (1) CN106126882B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107291798A (zh) * 2017-05-09 2017-10-24 吉林大学 一种地球化学数据异常信息的提取方法及系统
CN107451226B (zh) * 2017-07-19 2019-12-17 中国地质大学(武汉) 一种基于gis的网格系统及用其制作综合异常剖析图的方法
CN110263791B (zh) * 2019-05-31 2021-11-09 北京京东智能城市大数据研究院 一种识别功能区的方法和装置
CN112990322A (zh) * 2021-03-23 2021-06-18 中国地质调查局发展研究中心 一种确定地球化学异常识别临界值的方法、装置及设备
CN113408019B (zh) * 2021-06-23 2022-04-22 河北地质大学 一种基于bme-gwr的水系沉积物地球化学异常制图方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101694671A (zh) * 2009-10-27 2010-04-14 中国地质大学(武汉) 一种基于地学栅格图像的空间加权主成分分析的方法
CN101739396A (zh) * 2008-11-14 2010-06-16 电子科技大学 基于不确定性空间数据挖掘的区域成矿预测方法
CN103529486A (zh) * 2013-07-02 2014-01-22 核工业北京地质研究院 一种地球化学异常圈定方法
WO2014179832A1 (en) * 2013-05-08 2014-11-13 Technological Resources Pty Ltd A method of, and a system for, controlling a drilling operation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101739396A (zh) * 2008-11-14 2010-06-16 电子科技大学 基于不确定性空间数据挖掘的区域成矿预测方法
CN101694671A (zh) * 2009-10-27 2010-04-14 中国地质大学(武汉) 一种基于地学栅格图像的空间加权主成分分析的方法
WO2014179832A1 (en) * 2013-05-08 2014-11-13 Technological Resources Pty Ltd A method of, and a system for, controlling a drilling operation
CN103529486A (zh) * 2013-07-02 2014-01-22 核工业北京地质研究院 一种地球化学异常圈定方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
地球化学异常评价及成矿远景预测;吕军阳;《中国优秀硕士学位论文全文数据库 基础科学辑》;20150815;摘要,第一章,第三章,图3-8,表3.1
水热因子对退化草原羊草恢复演替群落土壤呼吸的影响;王风玉等;《植物生态学报》;20031231;第27卷(第5期);第644-649
西秦岭寨上金矿原生晕地球化学特征及成矿预测;刘新会等;《黄金科学技术》;20120831;第20卷(第4期);第7-15页
西部高寒山区遥感与化探信息综合找矿定位预测研究;杨自安;《中国优秀博硕士学位论文全文数据库(博士) 基础科学辑》;20051015;第四章,表4-4-6

Also Published As

Publication number Publication date
CN106126882A (zh) 2016-11-16

Similar Documents

Publication Publication Date Title
CN106126882B (zh) 基于肯德尔和谐系数排序的元素地球化学数据优选方法
Zhou et al. Application of cluster analysis to geochemical compositional data for identifying ore-related geochemical anomalies
Zheng et al. Analysis of stream sediment data for exploring the Zhunuo porphyry Cu deposit, southern Tibet
Baldeck et al. Environmental drivers of tree community turnover in western Amazonian forests
Chork et al. Integrating a high-breakdown option into discriminant analysis in exploration geochemistry
CN111639067B (zh) 多特征融合卷积自编码的多元地球化学异常识别方法
Liu et al. Maximum entropy modeling for orogenic gold prospectivity mapping in the Tangbale-Hatu belt, western Junggar, China
Kord et al. Spatial analysis of Ardabil plain aquifer potable groundwater using fuzzy logic
CN110544286B (zh) 一种大尺度范围的场地快速分类方法及构建场地分类图的装置
Huang et al. Vertical structure-based classification of oceanic eddy using 3-D convolutional neural network
Haddadchi et al. Grid‐based sediment tracing approach to determine sediment sources
CN108595749A (zh) 一种使用变异函数单一方向结构分析的资源储量评估方法
Pelánková et al. The relationships of modern pollen spectra to vegetation and climate along a steppe–forest–tundra transition in southern Siberia, explored by decision trees
Yates et al. 1.5 Geostatistics
CN103903057B (zh) 基于元素地球化学异常的碳硅泥岩型铀矿靶区优选方法
Nazzal et al. A combined hydrochemical-statistical analysis of Saq aquifer, northwestern part of the Kingdom of Saudi Arabia
Sadeghi Evaluation of geochemical anomaly classification models based on the relevant uncertainties and error propagation per class to select the most robust model (s) for the follow-up exploration
Hani et al. Spatial variability of heavy metals in the soils of Ahwaz using geostatistical methods
DE VARGAS et al. HYDROGEOLOGY OF FRACTURED AQUIFERS: APPLICATION OF CONSISTENCY INDEXES FOR THE VALIDATION OF GEOSPATIAL MATHEMATICAL MODELS: Hidrogeologia de aquíferos fraturados: aplicação de índices de consistência para validação de modelos matemáticos geoespaciais
Mokhtari et al. Geochemical prospecting for Cu mineralization in an arid terrain-central Iran
CN114578448A (zh) 基于多地源地质参数的考察点定位方法
Keykha Hoseinpoor et al. The use of robust factor analysis of compositional geochemical data for the recognition of the target area in Khusf 1: 100000 sheet, South Khorasan, Iran
CN113779888A (zh) 地面沉降危险性评估方法、装置、设备及存储介质
Vergiev et al. Key parameters for landscape evolution and anthropogenisation estimation in the Kamchia River downstream region (Eastern Bulgaria)
CN114428290A (zh) 利用地震资料估算压扭构造运动学涡度的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190129

Termination date: 20210615

CF01 Termination of patent right due to non-payment of annual fee