吉林省和龙市南坪地区钼矿综合信息成矿预测研究
|
摘要
南坪地区在吉林省延边朝鲜族自治州和龙市境内，大地构造上位于中朝准地台北缘东段之和龙地块，属古亚洲和濒西太平洋构造成矿域的结合部位。是中国地质调查局确定的16个重点成矿区带之一，地处辽东～吉南成矿带东部，在夹皮沟～和龙金铜镍成矿远景区与抚顺～和龙太古代铜金绿岩建造分布区的交叉部位，因而区域成矿条件良好。
目前众多学者对成矿预测做了大量研究且各有特色，但我国目前成矿规律和成矿预测研究工作仍处于经验性的总结阶段，没有形成一套统一的理论体系，随着各种新理论、新方法的出现，成矿预测理论体系需要不断完善。针对这一现状本文选定和龙市南坪地区作为研究区域，以综合信息预测理论为指导思想，全面利用物探、化探所显示的地质找矿信息，通过对与区域成矿有关的基础地质问题的研究，实现了研究区钼资源定量预测评价，证明和龙东部南坪地区依然拥有巨大的找矿潜力可以发掘，为后续的找矿靶区优选及其资源评价提供了科学依据。
本研究以综合成矿预测理论为指导，对南坪地区进行综合成矿预测，对研究区的成矿规律认识和找矿指导有所创新,对指导区域找矿勘查具有一定的实际意义。
作者主要做了以下的研究工作：
1、总结了南坪地区的成矿规律，初步查明研究区地层、岩浆岩、构造等基础地质特征及与成矿有关的主要地质因素。对南坪地区的典型矿床～石马洞矿床含矿岩体的规模、形态、产状及蚀变矿化特征进行了系统的分析和提炼。
2、通过改进谢别德法对二维经验模态分解的包络方式，提高了其运算速度，并利用其对区域地球物理、地球化学数据进行了分析，提出新的基于经验模态分解的物、化探异常提取方法。能够提取比常规方法更丰富的与成矿有关的地球物理、地球化学信息，可从多个层面认识地球物理场、地球化学场的特征，对成矿远景区的圈定有重要作用。
3、利用解析延拓等方法对研究区1：5万高精度磁测数据等资料进行了处理。对比向上延拓与二维经验模态分解后的结果可知，二维经验模态分解获得的IMF分量能更精确、客观地揭示和龙市南坪地区的地球物理（地质）结构特征。二维经验模态分解可以达到场分解的效果，不会出现信息混淆的情况。对IMF分量进行水平方向导数、垂向导数处理，通过研究区不同地质单元磁场空间形态的描述，对磁异常进行推断解释，重建研究区的构造框架，并在此基础上编制了磁性体分布图和构造解释图。
4、利用南坪地区1：5万水系沉积物地球化学数据，合理确定Au、Ag、Mo、Cu、W等元素的异常下限，通过因子分析、相关分析、聚类分析等多元统计方法研究区域内各元素的分布富集特征及共生组合规律，确定地球化学异常找矿标志。并结合地球化学数据特点编制了单元素地球化学异常图、组合元素地球化学异常图、多元素累加地球化学异常图。
5、研究工作基于GIS平台，结合石马洞典型矿床、成矿规律、地层、岩浆岩、构造等地质特征，以及物化探异常信息的研究，对地质统计单元进行了划分。提取了与矿床形成关系密切的各种能够直接或间接指示地下矿床存在的找矿标志作为预测变量，进一步通过地质—数学模型转换成为预测变量，对预测变量进行优化处理，建立了基于石马洞典型矿床的地、物、化综合找矿模型。
6、以成矿预测理论为指导思想，运用模糊神经网络方法，结合找矿模型编制了综合信息成矿预测图，确定和优选了找矿靶区。根据不同的成矿条件，将找矿预测单元划分为三级：A级（找矿有利度>0.7）、B级（0.7>找矿有利度>0.6）、C级（0.6>找矿有利度>0.5）。其中A级预测单元对于找矿最为有利，最有希望发现钼矿床（体）；B级预测单元找矿标志较为明显，有可能发现钼矿床（体）；C级预测单元具有一定的找矿标志，有发现钼矿床的可能性。
7、利用基于蒙特卡洛法的定量预测方法对矿产资源加以定量评价。经过模拟预测估算，吉林省和龙市南坪地区钼矿资源量在概率60%处大约为120000吨；80%处大约为65000吨。该区域已经探明的储量为48956吨，因此确定仍然有着大量的未知潜力矿体可供开发。
The Nanping area is located in the Yanbian Korean Autonomous Prefecture andHelong city.It is on the Helong block which is on the northern margin of Sino-Koreanparaplatform and plots, is the ancient Asian metallogenic domain. The Nanping area isone of the16Metallogenic belts which are identified by the China Geologicalsurvey.This area is in the eastern part of the Liaodong-Southern Jilin ore belt, andeast-and Jin Tongnie, the mineralization and Fushun-and too ancient bronzegreenstone formation distribution area of cross section, regional metallogeniccondition is good.
At present, many scholars have done a lot of research and have their owncharacteristics. But the metallogenic regularities and metallogenic prediction researchin our country are still limited to empirical summary stage, did not formed a unifiedtheory system. With the emergence of a variety of new theories,new methods still need to constantly improve. This paper selected the Nanping areain Helong city and with integrated information prediction theory as the guide,strengthen and mineralization related to basic geological research, strengthen theregional metallogenic regularity, full use of geophysical, geochemical shows thegeological prospecting information; the forecasting technology, qualitative predictionin delineation basis, realization study on prediction of resource potential evaluation ofdistrict molybdenum, as follow-up to the optimization of prospecting targets in andprovide scientific basis for resource evaluation. This paper mainly focuses on thefollowing research work:
1. At first, I summarize the Nanping regional metallogenic regularity.Preliminary I identify the basic geological characteristics of the stratigraphic,magmatic, tectonic and metallogenic geological’s factors. Papers on the Nanpingregion characteristic of ore deposits-ShiMa hole ore potentiality of rock mass scale,shape, form in altered mineralization and characteristics of the system are analyzedand refining.
2. The2d empirical mode decomposition is a kind of decomposition method ofnonlinear data. This article uses the improved Xie Biede method, envelope methodand operation speed improvement. Through the two-dimensional empirical modedecomposition to obtain rich geochemical information for from multiple aspects of geophysical, geochemical field characteristics, extraction and mineralization related tocontent, metallogenic prospective area of objective prediction has important role.
3. I collected1:50000geophysical (magnetic) data in the Nanping area, and usedjust like analytic continuation methods deal with the study area1:50000high precisionmagnetic survey data. I compared the upward continuation results with the empiricalmode decomposition, application of empirical mode decomposition for IMFcomponent more accurately, objectively reveals and city Nanping area geophysical(Geological) structural characteristics. And the application of empirical modedecomposition can reach field decomposition effect, there is no information confusedsituation, this article on the basis of research on the IMF components of the horizontaldirection, vertical direction derivative derivative processing, based on the differentgeological unit magnetic field space form of description, inference and interpretationof magnetic anomaly, reconstruction of study area tectonic framework, and on thebasis of the preparation of the magnetic body distribution and structure interpretationmap.
4. Collected the Nanping1:50000stream sediment geochemical data, anddetermined the lower limited of the elements such as Au, Ag, Mo, Cu, W. This articleuse the factor analysis,correlation analysis, cluster analysis, multivariate statisticalmethods the distribution and enrichment of various elements in the studyarea characteristics and paragenesis, and identify the geochemical anomalyprospecting criteria. Combined with the geochemical data on the characteristics of thepreparation of a single-element geochemical anomaly maps, anomaly maps of thecomposite element geochemical anomaly map, multi-element geochemicalaccumulation.
5. The reasearch work is based on the GIS platform,combined withthe Shima hole typical deposit, and metallogenic regularity of the known in the studyarea and stratigraphic, magmatic,tectonicand other geological features,and geologicalfeatures, and geophysical and geochemical anomalies, geostatistical unit has beendivided, extract and deposit close gelogicalsigns,be able to the geophysicalsign,geochemical,and other direct or indirectinstructions underground existenceof the deposit prospecting indicator as a predictor variable, and further convertedinto geological–mathematical model of the predictor variables to optimizethe predictor variable established based on the Shima hole typical deposits in theground,things of comprehensive prospecting model.
6. Metallogenic prediction theory, the use of fuzzy neural network, combinedwith prospecting model based on the preparation of comprehensive informationmineralization forecast map to determine the mineralization prospect area, andpreferably prospecting targets. Depending on the ore-forming conditions,prospectingprediction unit is divided into three levels: level A(prospecting favorable>0.7),GradeB(0.7> prospecting beneficial degree>0.6), C-level (0.6> favorable prospecting> 0.5).Class A prediction unit is most favorable for prospecting, most have foundthat themolybdenum deposit (body); more obvious prediction unit B prospectingcriteria, may find that the molybdenum deposit (body); prediction unit C has acertain prospectingsigns, there is the possibility to find molybdenum deposit.
7. Use of quantitative forecasting methods based on the Monte Carlo method tobequantitative evaluation of mineral resources, resource prediction, and then on theNanping area molybdenum reserves. Monte Carlo simulation and predictionestimates, Jilin Province and the Dragon Nanping molybdenum ore resources about120,000tons in the probability of60％;80％at about65,000tons, the region hasproven reserves of48,956tons, soto determine the study area is still a large number ofunknown potential ore body can be developed.
In summary, the Nanping area in the east of the Helong city still has a hugeprospecting potential to explore, this study metallogenic prediction comprehensiveexploration of the study area has played a guidingrole to provide practical theoreticalbasis for further detailed investigation of evaluation.
目前众多学者对成矿预测做了大量研究且各有特色，但我国目前成矿规律和成矿预测研究工作仍处于经验性的总结阶段，没有形成一套统一的理论体系，随着各种新理论、新方法的出现，成矿预测理论体系需要不断完善。针对这一现状本文选定和龙市南坪地区作为研究区域，以综合信息预测理论为指导思想，全面利用物探、化探所显示的地质找矿信息，通过对与区域成矿有关的基础地质问题的研究，实现了研究区钼资源定量预测评价，证明和龙东部南坪地区依然拥有巨大的找矿潜力可以发掘，为后续的找矿靶区优选及其资源评价提供了科学依据。
本研究以综合成矿预测理论为指导，对南坪地区进行综合成矿预测，对研究区的成矿规律认识和找矿指导有所创新,对指导区域找矿勘查具有一定的实际意义。
作者主要做了以下的研究工作：
1、总结了南坪地区的成矿规律，初步查明研究区地层、岩浆岩、构造等基础地质特征及与成矿有关的主要地质因素。对南坪地区的典型矿床～石马洞矿床含矿岩体的规模、形态、产状及蚀变矿化特征进行了系统的分析和提炼。
2、通过改进谢别德法对二维经验模态分解的包络方式，提高了其运算速度，并利用其对区域地球物理、地球化学数据进行了分析，提出新的基于经验模态分解的物、化探异常提取方法。能够提取比常规方法更丰富的与成矿有关的地球物理、地球化学信息，可从多个层面认识地球物理场、地球化学场的特征，对成矿远景区的圈定有重要作用。
3、利用解析延拓等方法对研究区1：5万高精度磁测数据等资料进行了处理。对比向上延拓与二维经验模态分解后的结果可知，二维经验模态分解获得的IMF分量能更精确、客观地揭示和龙市南坪地区的地球物理（地质）结构特征。二维经验模态分解可以达到场分解的效果，不会出现信息混淆的情况。对IMF分量进行水平方向导数、垂向导数处理，通过研究区不同地质单元磁场空间形态的描述，对磁异常进行推断解释，重建研究区的构造框架，并在此基础上编制了磁性体分布图和构造解释图。
4、利用南坪地区1：5万水系沉积物地球化学数据，合理确定Au、Ag、Mo、Cu、W等元素的异常下限，通过因子分析、相关分析、聚类分析等多元统计方法研究区域内各元素的分布富集特征及共生组合规律，确定地球化学异常找矿标志。并结合地球化学数据特点编制了单元素地球化学异常图、组合元素地球化学异常图、多元素累加地球化学异常图。
5、研究工作基于GIS平台，结合石马洞典型矿床、成矿规律、地层、岩浆岩、构造等地质特征，以及物化探异常信息的研究，对地质统计单元进行了划分。提取了与矿床形成关系密切的各种能够直接或间接指示地下矿床存在的找矿标志作为预测变量，进一步通过地质—数学模型转换成为预测变量，对预测变量进行优化处理，建立了基于石马洞典型矿床的地、物、化综合找矿模型。
6、以成矿预测理论为指导思想，运用模糊神经网络方法，结合找矿模型编制了综合信息成矿预测图，确定和优选了找矿靶区。根据不同的成矿条件，将找矿预测单元划分为三级：A级（找矿有利度>0.7）、B级（0.7>找矿有利度>0.6）、C级（0.6>找矿有利度>0.5）。其中A级预测单元对于找矿最为有利，最有希望发现钼矿床（体）；B级预测单元找矿标志较为明显，有可能发现钼矿床（体）；C级预测单元具有一定的找矿标志，有发现钼矿床的可能性。
7、利用基于蒙特卡洛法的定量预测方法对矿产资源加以定量评价。经过模拟预测估算，吉林省和龙市南坪地区钼矿资源量在概率60%处大约为120000吨；80%处大约为65000吨。该区域已经探明的储量为48956吨，因此确定仍然有着大量的未知潜力矿体可供开发。
The Nanping area is located in the Yanbian Korean Autonomous Prefecture andHelong city.It is on the Helong block which is on the northern margin of Sino-Koreanparaplatform and plots, is the ancient Asian metallogenic domain. The Nanping area isone of the16Metallogenic belts which are identified by the China Geologicalsurvey.This area is in the eastern part of the Liaodong-Southern Jilin ore belt, andeast-and Jin Tongnie, the mineralization and Fushun-and too ancient bronzegreenstone formation distribution area of cross section, regional metallogeniccondition is good.
At present, many scholars have done a lot of research and have their owncharacteristics. But the metallogenic regularities and metallogenic prediction researchin our country are still limited to empirical summary stage, did not formed a unifiedtheory system. With the emergence of a variety of new theories,new methods still need to constantly improve. This paper selected the Nanping areain Helong city and with integrated information prediction theory as the guide,strengthen and mineralization related to basic geological research, strengthen theregional metallogenic regularity, full use of geophysical, geochemical shows thegeological prospecting information; the forecasting technology, qualitative predictionin delineation basis, realization study on prediction of resource potential evaluation ofdistrict molybdenum, as follow-up to the optimization of prospecting targets in andprovide scientific basis for resource evaluation. This paper mainly focuses on thefollowing research work:
1. At first, I summarize the Nanping regional metallogenic regularity.Preliminary I identify the basic geological characteristics of the stratigraphic,magmatic, tectonic and metallogenic geological’s factors. Papers on the Nanpingregion characteristic of ore deposits-ShiMa hole ore potentiality of rock mass scale,shape, form in altered mineralization and characteristics of the system are analyzedand refining.
2. The2d empirical mode decomposition is a kind of decomposition method ofnonlinear data. This article uses the improved Xie Biede method, envelope methodand operation speed improvement. Through the two-dimensional empirical modedecomposition to obtain rich geochemical information for from multiple aspects of geophysical, geochemical field characteristics, extraction and mineralization related tocontent, metallogenic prospective area of objective prediction has important role.
3. I collected1:50000geophysical (magnetic) data in the Nanping area, and usedjust like analytic continuation methods deal with the study area1:50000high precisionmagnetic survey data. I compared the upward continuation results with the empiricalmode decomposition, application of empirical mode decomposition for IMFcomponent more accurately, objectively reveals and city Nanping area geophysical(Geological) structural characteristics. And the application of empirical modedecomposition can reach field decomposition effect, there is no information confusedsituation, this article on the basis of research on the IMF components of the horizontaldirection, vertical direction derivative derivative processing, based on the differentgeological unit magnetic field space form of description, inference and interpretationof magnetic anomaly, reconstruction of study area tectonic framework, and on thebasis of the preparation of the magnetic body distribution and structure interpretationmap.
4. Collected the Nanping1:50000stream sediment geochemical data, anddetermined the lower limited of the elements such as Au, Ag, Mo, Cu, W. This articleuse the factor analysis,correlation analysis, cluster analysis, multivariate statisticalmethods the distribution and enrichment of various elements in the studyarea characteristics and paragenesis, and identify the geochemical anomalyprospecting criteria. Combined with the geochemical data on the characteristics of thepreparation of a single-element geochemical anomaly maps, anomaly maps of thecomposite element geochemical anomaly map, multi-element geochemicalaccumulation.
5. The reasearch work is based on the GIS platform,combined withthe Shima hole typical deposit, and metallogenic regularity of the known in the studyarea and stratigraphic, magmatic,tectonicand other geological features,and geologicalfeatures, and geophysical and geochemical anomalies, geostatistical unit has beendivided, extract and deposit close gelogicalsigns,be able to the geophysicalsign,geochemical,and other direct or indirectinstructions underground existenceof the deposit prospecting indicator as a predictor variable, and further convertedinto geological–mathematical model of the predictor variables to optimizethe predictor variable established based on the Shima hole typical deposits in theground,things of comprehensive prospecting model.
6. Metallogenic prediction theory, the use of fuzzy neural network, combinedwith prospecting model based on the preparation of comprehensive informationmineralization forecast map to determine the mineralization prospect area, andpreferably prospecting targets. Depending on the ore-forming conditions,prospectingprediction unit is divided into three levels: level A(prospecting favorable>0.7),GradeB(0.7> prospecting beneficial degree>0.6), C-level (0.6> favorable prospecting> 0.5).Class A prediction unit is most favorable for prospecting, most have foundthat themolybdenum deposit (body); more obvious prediction unit B prospectingcriteria, may find that the molybdenum deposit (body); prediction unit C has acertain prospectingsigns, there is the possibility to find molybdenum deposit.
7. Use of quantitative forecasting methods based on the Monte Carlo method tobequantitative evaluation of mineral resources, resource prediction, and then on theNanping area molybdenum reserves. Monte Carlo simulation and predictionestimates, Jilin Province and the Dragon Nanping molybdenum ore resources about120,000tons in the probability of60％;80％at about65,000tons, the region hasproven reserves of48,956tons, soto determine the study area is still a large number ofunknown potential ore body can be developed.
In summary, the Nanping area in the east of the Helong city still has a hugeprospecting potential to explore, this study metallogenic prediction comprehensiveexploration of the study area has played a guidingrole to provide practical theoreticalbasis for further detailed investigation of evaluation.
引文
[1]黄崇轲,白冶,朱裕生等.中国铜矿床(下册)[M].地质出版社,2001:373-391.
[2]刘士倜.国外成矿预测工作简介[J].世界核地质科学,1978,1:33-35.
[3]赵鹏大.成矿定量预测与深部找矿[J].地学前缘,2007,5:1-10.
[4] Garrett, R.G.. Abundance model resource appraisa for Some CanadianCommodities [J]. Math. Geol,1978,10-14.
[5]赵鹏大,池顺都.当今矿产勘察问题的思考[J].地球科学,1998,23（1）:70-74.
[6]宋国耀,张晓华,肖克炎,朱裕生.矿产资源潜力评价的理论和GIS技术[J].物探化探计算技术,1999,3:199-205.
[7]肖克炎,丁建华,刘锐.美国“三步式”固体矿产资源潜力评价方法评述[J].地质论评,2006,52(6)：793-798.
[8]成秋明.成矿过程奇异性与矿产预测定量化的新理论与新方法[J].地学前缘,2007,14（5）：42-53.
[9]成秋明.非线性成矿预测理论:多重分形奇异性-广义自相似性-分形谱系模型与方法[J].地球科学,2006,3：338-348.
[10]雷万杉.内蒙古赤峰南部地区金矿综合信息矿产预测[D].长春：吉林大学,2009.
[11]程裕淇,陈毓川,赵一鸣.初论矿床的成矿系列问题[J].中国地质科学院院报,1979:32-58.
[12]贾大成,胡瑞忠,裴春雷,冯本智利用成矿系列对延边金苍地区金矿的成矿预测[J].地球与环境,2003,01：42-46.
[13]涂光炽.关于超大型矿床的寻找和理论研究[J].地球科学进展,1989,6:14-20.
[14] Xie Xuejing. New concept and new technique in searching for giant ore deposit[J]. Scientific Chinese,1995(5):14-16.
[15] Xie Xuejing. The surfacial geochemical expressions of Giant ore-deposits [J].Whiting, B H, Hodgson C J and Mason R (Editors), Giant Ore Deposits II,Queen’s University, Kingston,1995:479-492.
[16]王世称,王於天．综合信息解译原理与矿产预测图编制方法[M]．长春:吉林大学出版社,1989：3-18.
[17]王世称,许亚光,侯惠群.综合信息成矿系列预测的基本思路与方法[J].中国地质,1992,10:14-16.
[18]赵鹏大,池顺都.初论地质异常[J].地球科学,1991,3:3-10.
[19]成秋明.应用复杂性一非线性理论开展成矿预测——奇异性理论—广义自相似性—分形谱系多重分形与应用.第八届全国矿床会议论文集[C].北京：中国地质学会,2006.
[20]吉林省第六地质调查所.吉林省和龙市白石洞铁矿详查地质报告[R].2006.
[21]沈阳地质矿产研究所.和龙市幅K52C0020021/25万区域地质调查报告[R].吉林：吉林省地质调查研究院,2004.
[22]吉林地质局延边地区综合地质大队.吉林省和龙县北东矿区铜镍矿普查报告[R].吉林：吉林省地质调查研究院,1974.
[23]刘洪文,孟广才,王振贵等.和龙煤矿幅K-52-54-B1/5万区域地质调查说明书[R].吉林：地矿部吉林地勘局第二地质调查所,1995.
[24]吉林区域地质调查大队.1:20万白头山幅(K-52-XV)区域地质报告[R].吉林：吉林区域地质调查大队,1974.
[25]吉林区域地质调查大队.1:20万桦树林子幅(K-52-VⅡ)区域地质报告[R].1972.
[26]吉林区域地质调查大队.1:20万明月幅(K-52-ⅠX)区域地质报告[R].1973.
[27]王秀英,潘月栋,杨涛,张学谦.激电中梯和激电测深在吉林和龙市石马洞钼矿的应用[J].吉林地质,2011,3：78-83.
[28]贾大成.青龙村群及其地质时代[J].吉林地质,1994,1:31-34.
[29]张春艳,张兴洲,邱殿明.延边地区青龙村群斜长角闪岩中锆石U-Pb同位素年龄及地质意义[R].吉林大学学报（地球科学版）,2007,4:42-47.
[30]党兴彦,彭齐鸣,赵俊伟等.全国主要成矿远景区矿产资源评价重点选区研究[M].北京：中国地质调查局,2003,3
[31]李东津.吉林省岩石地层[M].北京：中国地质大学出版社,1997.
[32]张春华.河北太行山阜平期变质作用分析[J].中国区域地质,1987,2:27-34.
[33]张艳斌,吴福元,翟明国,路孝平.和龙地块的构造属性与华北地台北缘东段边界[R].中国科学D辑,2004,9:14-25.
[34]黄光南.希尔伯特变换及其在地震资料分析处理中的应用[D].青岛：中国海洋大学,2009.
[35]胡红英.局域波分解方法、特征剖析及应用研究[D].大连：大连理工大学,2006.
[36]钟佑明,秦树人.希尔伯特-黄变换的统一理论依据研究[J].振动与冲击,2006,25(3):40-43.
[37] Huang N E,Zhaohua Wu. A study of the characteristics of white noise using theempirical mode decomposition method[J]. Proceedings of the Royal Society,2004,460:1597-1611.
[38] Wu Zhaohua,Huang N E. A Study of the Characteristics of White Noise Usingthe Empirical Mode Decomposition Method[J]. Proceedings of the RoyalSociety,2003,1:1-27.
[39] Stephen C P, Robert J G, Jonathan W E. Application of the Hilbert-HuangTransform to the Analysis of Molecular Dynamics Simulations[J]. PhysChem,2003,107:4869-4876.
[40] Patrick Flandrin,Paulo Goncalves,Gabriel Rilling.Detrending and denoisingwith empirical mode decompositions[J].Pro.EUSIPCO,2004:20-30
[41] Tanaka T, Mandic D P. Complex empirical mode decomposition [J]. IEEESignal Processing Letters,2007,14(2):101―104.
[42] Rilling G, Flandrin P, Gon alves P,Lilly J M. Bivariate empirical modedecomposition [J]. IEEE Signal Processing Letters,2007,14(12):936-939.
[43] A.Linderhed.2D empirical mode decomposition in the spirit of imagecompression[J].Wavelet and Independent Component Analysis Application IXOrlando,2002,4.
[44] Zhu Xun, Shen Zheng, Eckermann S D, et al. Gravity wave characteristics inthe middle atmosphere derived from the Empirical Mode Decompositionmethod[J]. Geophysical Research,1997,102(14):545-561.
[45]周文纳,曾昭发,杜晓娟,等.基于经验模态分解的重力异常分离[J].世界地质,2010,29(3):495-502.
[46]曾琴琴,刘天佑.一种基于经验模态分解(EMD)的位场分离方法[J].石油地球物理勘探,2010,45(6):914-917.
[47]宋立新,高凤娇,郗朝晖.二维EMD分解方法的比较与改进[J].电子与信息学报,2008,12:2890-2893.
[48]孙艳争. EMD时频分析理论与应用研究[D].成都：电子科技大学,2007.
[49] Linderhed A.2D empirical mode decomposition in the spirit of imagecompression [J]. Wavelet and Independent Component Analysis Application IX,2002,4738:1-8.
[50] Nunes J C, Guyot S, and Del′echelle E. Texture analysis based on localanalysis of the Bidimensional Empirical Mode Decomposition [J]. MachineVision and Applications,2005,16:177-188.
[51]陈建国,肖凡,常韬.基于二维经验模态分解的重磁异常分离[J].地球科学,2011.2:327-335.
[52] Franke, R., Nielson, G.. Smooth Interpolation of Large Sets of Scattered Data[J].International Journal for Numerical Methods in Engineering,1980,15:1691-1704.
[53]王慧,刘正士.基于希尔伯特-黄变换的信号总失真度评价[J].华中科技大学学报,2008,10:302-305.
[54] Huang N E et al. A confidence limit for the empirical mode decomposition andHilbert spectral analysis[J]. Proceedings of the Royal Society of LondonSE-RIESA,2003,459:2317-2345.
[55] Flandrin P, Rilling G et al. Empirical mode decompo-sition as filter bank[j].IEEE Signal Process Lett,2003,11:112-114.
[56]何丽娟,张毅,张焱.二维ＥＭＤ方法在物化探数据处理中的应用研究——以云南个旧地区为例[J].地质找矿论丛,2011,9:311-315.
[57]黄静宁.深层次致矿异常信息提取及其找矿应用研究[D].北京：中国地质大学,2010.
[58] Xie Xuejing, Yin Binchuan. Geochemical patterns from local to global[J].Journal of Geochemical Exploration,1993,47:109-129.
[59] Bolviken B, Demetriades A, Ebbing J, et al. Geochemical mapping of westernEurope towards the year2000[R]. Trondheim: Geological Survey ofNorway,1990.
[60] N.E. Huang. Review of empirical mode decomposition analysis[J]. Proc. OfSPIE,2001,4391:71-79.
[61]滕彦国,庹先国,倪师军,张成江,马玉孝.地球化学基线的确定方法研究—以攀枝花地区为例[J].成都理工大学学报（自然科学版）,2008,4:422-427.
[62]滕彦国,倪师军等.地球化学基线的理论与实践[M].北京：化学工业出版社,2007.
[63] Bauer I,Spernger M, Bor J. Die Berechnung Lithogener und geonenerSchwermetallgehalte von Lobboden am Beispielen von Cu, Zn und Pb [J].Mainzer Geowiss Mitt,1992,21:47-70.
[64] Bauer I, Bor J. Vertikale Bilanzierung von Schwermetallen in Boden-Kennzeichnung der Empfindlichkeit der boden gegenuber Schwermetallenunter Berucksichtigung von lithogenem Grundgehalt, pedogener An-undAbreicherung some antheopogener Zusatzbelastung, Teil2[M].Berlin:Texte56,Umweltbundesam,1993.
[65] BauerI,Bor J.Lithogene,geonene und anthropogene SchwermeTallgeh-alte vonLobboden an den Beispielen von Cu,Zn,Ni,Pb,Hg und Cd[J].Mainzer GeowissMitt,1995,24:47-70.
[66]曹殿华,王安建,王高尚等.勘查地球化学异常多尺度分析方法—以赣东北德兴矿集区为例[J].现代地质,2008,6(22):1028-1033.
[67]罗孝宽,郭绍雍.应用地球物理教程[M].北京:地质出版社,1991.
[68]王谦身,安玉林,张赤军等.重力学[M].北京:地震出版社,2003.
[69]管志宁.地磁场与磁力勘探[M].北京：地质出版社,2005.
[70]赵俊猛.天山造山带与准噶尔盆地岩石圈结构及其动力学过程[D].北京：中国地震局地质研究所,1998.
[71]赵俊猛,马宗晋,姚长利,黎益仕,刘占坡.准噶尔盆地基底构造分区的重磁异常分析[J].新疆石油地质,2008,1:7-11.
[72]长春地质学院,吉林省地矿局第六地质调查所.吉林省延边地区（槽区）花岗岩成因系列及其成矿作用科研报告[R].吉林：吉林省地质矿产局,1989.
[73]陈化奇.多元素衬值累加晕在化探异常评价中的应用[J].甘肃地质,2008,4:90-93.
[74]赵鹏大.矿产勘查理论与方法[M].北京：中国地质大学出版社,2010.
[75]李新中.矿床统计预测单元划分的方法与程序[J].矿床地质,1998,17(4):369-375.
[76]陈永清,夏庆霖.应用地质异常单元圈定矿产资源体潜在地段[J].地球科学一中国地质大学学报,1999,24(5):459-463.
[77]陈石羡.地理信息系统在金属矿产预测中的应用[J].地质找矿论丛,1998,13(l):74-83.
[78]王于天.成矿预测单元的基本概念及其划分方法[J].地质论评,1990,36(6):24-29.
[79]王雄军.云南老君山矿集区多因复成成矿模式及空间信息成矿预测模型研究[D].湖南：中南大学,2008.
[80]武汉地质学院劫探教研室数学地质组．矿床统计预测简介_地质变量的选择、取值和变换[J]．地质与勘探,1980：32-37.
[81]朱裕生,肖克炎等.成矿预测方法通则之二——成矿预测方法[M].北京：地质出版社,1997.
[82]曹恩生,徐元才.特征分析及其应用—赣东北地区铜矿预测[J].物探与化探,1985,6:33-38.
[83]肖克炎.矿产资源评价系统（MRAS）用户使用指南[M].北京：中国地质科学院,2003.
[84]李堃,胡光道,段其发.基于MORPAS平台特征分析法的成矿远景区预测——以个旧西区锡多金属矿为例[J].地质科技情报,2009,4:65-70.
[85] Liu,L.M.,Peng,S.L. Prediction of hidden ore bodies by synthesis ofgeological,geophysical and geochemical information based on dynamic modelin Fenghuangshan ore field,tongling district, China.[J]. Geochem.Explor,2003,81：81–98.
[86]戚莹.基于模糊神经网络的森林虫害预测预报的应用研究[D].黑龙江：东北林业大学,2011.
[87]万中,梁文冬,卢宗娟.模糊数的隶属度区间分布函数[J].重庆理工大学学报（自然科学版）,2011,1:107-112.
[88]吴小俊,曹奇英,王士同,陈保香,刘同明.模糊神经系统及其应用研究[J].系统工程,1999,4:61-64.
[89]龚利,王云宽,宋英华,曹志强.基于模糊神经网络的固体粉末自动配料控制系统[J].计量学报,2008,3:284-288.
[90]孙祥,杨子荣,徐大地.基于模糊神经网络的义县萤石成矿预测[R].地质找矿论丛,2011,02:149-152.
[91]姜德义,蒋再文,章永武,任松.模糊神经网络预测公路隧道揭煤突出危险性[J].2000,S1:200-202.
[92]付海涛,王恩德.GIS技术在成矿预测中的应用实例——以杨家杖子矿田为例[J].矿床地质,2005,24(6):684-691.
[93]陈建平,王功文,侯昌波等.基于GIS技术的西南三江北段矿产资源定量预测与评价[J].矿床地质,2005,24(1):15-24.
[94]杜灵通,吕新彪.GIS在矿产资源评价中的应用[J].地质找矿论丛,2003,18(4):275-278.
[95]中华人民共和国国家质量监督检验检疫总局.中华人民共和国国家标准固体矿产地质勘查规范总则[M].北京:中国标准出版社,2002.
[96]宋国耀,张晓华,肖克炎,等.矿产资源潜力评价的若干问题[J].中国地质,1999(8):17-20.
[97]朱海宾.蒙特卡洛模型在矿产资源量预测中的应用[J].地质找矿论丛,2010,1:50-54.
[2]刘士倜.国外成矿预测工作简介[J].世界核地质科学,1978,1:33-35.
[3]赵鹏大.成矿定量预测与深部找矿[J].地学前缘,2007,5:1-10.
[4] Garrett, R.G.. Abundance model resource appraisa for Some CanadianCommodities [J]. Math. Geol,1978,10-14.
[5]赵鹏大,池顺都.当今矿产勘察问题的思考[J].地球科学,1998,23（1）:70-74.
[6]宋国耀,张晓华,肖克炎,朱裕生.矿产资源潜力评价的理论和GIS技术[J].物探化探计算技术,1999,3:199-205.
[7]肖克炎,丁建华,刘锐.美国“三步式”固体矿产资源潜力评价方法评述[J].地质论评,2006,52(6)：793-798.
[8]成秋明.成矿过程奇异性与矿产预测定量化的新理论与新方法[J].地学前缘,2007,14（5）：42-53.
[9]成秋明.非线性成矿预测理论:多重分形奇异性-广义自相似性-分形谱系模型与方法[J].地球科学,2006,3：338-348.
[10]雷万杉.内蒙古赤峰南部地区金矿综合信息矿产预测[D].长春：吉林大学,2009.
[11]程裕淇,陈毓川,赵一鸣.初论矿床的成矿系列问题[J].中国地质科学院院报,1979:32-58.
[12]贾大成,胡瑞忠,裴春雷,冯本智利用成矿系列对延边金苍地区金矿的成矿预测[J].地球与环境,2003,01：42-46.
[13]涂光炽.关于超大型矿床的寻找和理论研究[J].地球科学进展,1989,6:14-20.
[14] Xie Xuejing. New concept and new technique in searching for giant ore deposit[J]. Scientific Chinese,1995(5):14-16.
[15] Xie Xuejing. The surfacial geochemical expressions of Giant ore-deposits [J].Whiting, B H, Hodgson C J and Mason R (Editors), Giant Ore Deposits II,Queen’s University, Kingston,1995:479-492.
[16]王世称,王於天．综合信息解译原理与矿产预测图编制方法[M]．长春:吉林大学出版社,1989：3-18.
[17]王世称,许亚光,侯惠群.综合信息成矿系列预测的基本思路与方法[J].中国地质,1992,10:14-16.
[18]赵鹏大,池顺都.初论地质异常[J].地球科学,1991,3:3-10.
[19]成秋明.应用复杂性一非线性理论开展成矿预测——奇异性理论—广义自相似性—分形谱系多重分形与应用.第八届全国矿床会议论文集[C].北京：中国地质学会,2006.
[20]吉林省第六地质调查所.吉林省和龙市白石洞铁矿详查地质报告[R].2006.
[21]沈阳地质矿产研究所.和龙市幅K52C0020021/25万区域地质调查报告[R].吉林：吉林省地质调查研究院,2004.
[22]吉林地质局延边地区综合地质大队.吉林省和龙县北东矿区铜镍矿普查报告[R].吉林：吉林省地质调查研究院,1974.
[23]刘洪文,孟广才,王振贵等.和龙煤矿幅K-52-54-B1/5万区域地质调查说明书[R].吉林：地矿部吉林地勘局第二地质调查所,1995.
[24]吉林区域地质调查大队.1:20万白头山幅(K-52-XV)区域地质报告[R].吉林：吉林区域地质调查大队,1974.
[25]吉林区域地质调查大队.1:20万桦树林子幅(K-52-VⅡ)区域地质报告[R].1972.
[26]吉林区域地质调查大队.1:20万明月幅(K-52-ⅠX)区域地质报告[R].1973.
[27]王秀英,潘月栋,杨涛,张学谦.激电中梯和激电测深在吉林和龙市石马洞钼矿的应用[J].吉林地质,2011,3：78-83.
[28]贾大成.青龙村群及其地质时代[J].吉林地质,1994,1:31-34.
[29]张春艳,张兴洲,邱殿明.延边地区青龙村群斜长角闪岩中锆石U-Pb同位素年龄及地质意义[R].吉林大学学报（地球科学版）,2007,4:42-47.
[30]党兴彦,彭齐鸣,赵俊伟等.全国主要成矿远景区矿产资源评价重点选区研究[M].北京：中国地质调查局,2003,3
[31]李东津.吉林省岩石地层[M].北京：中国地质大学出版社,1997.
[32]张春华.河北太行山阜平期变质作用分析[J].中国区域地质,1987,2:27-34.
[33]张艳斌,吴福元,翟明国,路孝平.和龙地块的构造属性与华北地台北缘东段边界[R].中国科学D辑,2004,9:14-25.
[34]黄光南.希尔伯特变换及其在地震资料分析处理中的应用[D].青岛：中国海洋大学,2009.
[35]胡红英.局域波分解方法、特征剖析及应用研究[D].大连：大连理工大学,2006.
[36]钟佑明,秦树人.希尔伯特-黄变换的统一理论依据研究[J].振动与冲击,2006,25(3):40-43.
[37] Huang N E,Zhaohua Wu. A study of the characteristics of white noise using theempirical mode decomposition method[J]. Proceedings of the Royal Society,2004,460:1597-1611.
[38] Wu Zhaohua,Huang N E. A Study of the Characteristics of White Noise Usingthe Empirical Mode Decomposition Method[J]. Proceedings of the RoyalSociety,2003,1:1-27.
[39] Stephen C P, Robert J G, Jonathan W E. Application of the Hilbert-HuangTransform to the Analysis of Molecular Dynamics Simulations[J]. PhysChem,2003,107:4869-4876.
[40] Patrick Flandrin,Paulo Goncalves,Gabriel Rilling.Detrending and denoisingwith empirical mode decompositions[J].Pro.EUSIPCO,2004:20-30
[41] Tanaka T, Mandic D P. Complex empirical mode decomposition [J]. IEEESignal Processing Letters,2007,14(2):101―104.
[42] Rilling G, Flandrin P, Gon alves P,Lilly J M. Bivariate empirical modedecomposition [J]. IEEE Signal Processing Letters,2007,14(12):936-939.
[43] A.Linderhed.2D empirical mode decomposition in the spirit of imagecompression[J].Wavelet and Independent Component Analysis Application IXOrlando,2002,4.
[44] Zhu Xun, Shen Zheng, Eckermann S D, et al. Gravity wave characteristics inthe middle atmosphere derived from the Empirical Mode Decompositionmethod[J]. Geophysical Research,1997,102(14):545-561.
[45]周文纳,曾昭发,杜晓娟,等.基于经验模态分解的重力异常分离[J].世界地质,2010,29(3):495-502.
[46]曾琴琴,刘天佑.一种基于经验模态分解(EMD)的位场分离方法[J].石油地球物理勘探,2010,45(6):914-917.
[47]宋立新,高凤娇,郗朝晖.二维EMD分解方法的比较与改进[J].电子与信息学报,2008,12:2890-2893.
[48]孙艳争. EMD时频分析理论与应用研究[D].成都：电子科技大学,2007.
[49] Linderhed A.2D empirical mode decomposition in the spirit of imagecompression [J]. Wavelet and Independent Component Analysis Application IX,2002,4738:1-8.
[50] Nunes J C, Guyot S, and Del′echelle E. Texture analysis based on localanalysis of the Bidimensional Empirical Mode Decomposition [J]. MachineVision and Applications,2005,16:177-188.
[51]陈建国,肖凡,常韬.基于二维经验模态分解的重磁异常分离[J].地球科学,2011.2:327-335.
[52] Franke, R., Nielson, G.. Smooth Interpolation of Large Sets of Scattered Data[J].International Journal for Numerical Methods in Engineering,1980,15:1691-1704.
[53]王慧,刘正士.基于希尔伯特-黄变换的信号总失真度评价[J].华中科技大学学报,2008,10:302-305.
[54] Huang N E et al. A confidence limit for the empirical mode decomposition andHilbert spectral analysis[J]. Proceedings of the Royal Society of LondonSE-RIESA,2003,459:2317-2345.
[55] Flandrin P, Rilling G et al. Empirical mode decompo-sition as filter bank[j].IEEE Signal Process Lett,2003,11:112-114.
[56]何丽娟,张毅,张焱.二维ＥＭＤ方法在物化探数据处理中的应用研究——以云南个旧地区为例[J].地质找矿论丛,2011,9:311-315.
[57]黄静宁.深层次致矿异常信息提取及其找矿应用研究[D].北京：中国地质大学,2010.
[58] Xie Xuejing, Yin Binchuan. Geochemical patterns from local to global[J].Journal of Geochemical Exploration,1993,47:109-129.
[59] Bolviken B, Demetriades A, Ebbing J, et al. Geochemical mapping of westernEurope towards the year2000[R]. Trondheim: Geological Survey ofNorway,1990.
[60] N.E. Huang. Review of empirical mode decomposition analysis[J]. Proc. OfSPIE,2001,4391:71-79.
[61]滕彦国,庹先国,倪师军,张成江,马玉孝.地球化学基线的确定方法研究—以攀枝花地区为例[J].成都理工大学学报（自然科学版）,2008,4:422-427.
[62]滕彦国,倪师军等.地球化学基线的理论与实践[M].北京：化学工业出版社,2007.
[63] Bauer I,Spernger M, Bor J. Die Berechnung Lithogener und geonenerSchwermetallgehalte von Lobboden am Beispielen von Cu, Zn und Pb [J].Mainzer Geowiss Mitt,1992,21:47-70.
[64] Bauer I, Bor J. Vertikale Bilanzierung von Schwermetallen in Boden-Kennzeichnung der Empfindlichkeit der boden gegenuber Schwermetallenunter Berucksichtigung von lithogenem Grundgehalt, pedogener An-undAbreicherung some antheopogener Zusatzbelastung, Teil2[M].Berlin:Texte56,Umweltbundesam,1993.
[65] BauerI,Bor J.Lithogene,geonene und anthropogene SchwermeTallgeh-alte vonLobboden an den Beispielen von Cu,Zn,Ni,Pb,Hg und Cd[J].Mainzer GeowissMitt,1995,24:47-70.
[66]曹殿华,王安建,王高尚等.勘查地球化学异常多尺度分析方法—以赣东北德兴矿集区为例[J].现代地质,2008,6(22):1028-1033.
[67]罗孝宽,郭绍雍.应用地球物理教程[M].北京:地质出版社,1991.
[68]王谦身,安玉林,张赤军等.重力学[M].北京:地震出版社,2003.
[69]管志宁.地磁场与磁力勘探[M].北京：地质出版社,2005.
[70]赵俊猛.天山造山带与准噶尔盆地岩石圈结构及其动力学过程[D].北京：中国地震局地质研究所,1998.
[71]赵俊猛,马宗晋,姚长利,黎益仕,刘占坡.准噶尔盆地基底构造分区的重磁异常分析[J].新疆石油地质,2008,1:7-11.
[72]长春地质学院,吉林省地矿局第六地质调查所.吉林省延边地区（槽区）花岗岩成因系列及其成矿作用科研报告[R].吉林：吉林省地质矿产局,1989.
[73]陈化奇.多元素衬值累加晕在化探异常评价中的应用[J].甘肃地质,2008,4:90-93.
[74]赵鹏大.矿产勘查理论与方法[M].北京：中国地质大学出版社,2010.
[75]李新中.矿床统计预测单元划分的方法与程序[J].矿床地质,1998,17(4):369-375.
[76]陈永清,夏庆霖.应用地质异常单元圈定矿产资源体潜在地段[J].地球科学一中国地质大学学报,1999,24(5):459-463.
[77]陈石羡.地理信息系统在金属矿产预测中的应用[J].地质找矿论丛,1998,13(l):74-83.
[78]王于天.成矿预测单元的基本概念及其划分方法[J].地质论评,1990,36(6):24-29.
[79]王雄军.云南老君山矿集区多因复成成矿模式及空间信息成矿预测模型研究[D].湖南：中南大学,2008.
[80]武汉地质学院劫探教研室数学地质组．矿床统计预测简介_地质变量的选择、取值和变换[J]．地质与勘探,1980：32-37.
[81]朱裕生,肖克炎等.成矿预测方法通则之二——成矿预测方法[M].北京：地质出版社,1997.
[82]曹恩生,徐元才.特征分析及其应用—赣东北地区铜矿预测[J].物探与化探,1985,6:33-38.
[83]肖克炎.矿产资源评价系统（MRAS）用户使用指南[M].北京：中国地质科学院,2003.
[84]李堃,胡光道,段其发.基于MORPAS平台特征分析法的成矿远景区预测——以个旧西区锡多金属矿为例[J].地质科技情报,2009,4:65-70.
[85] Liu,L.M.,Peng,S.L. Prediction of hidden ore bodies by synthesis ofgeological,geophysical and geochemical information based on dynamic modelin Fenghuangshan ore field,tongling district, China.[J]. Geochem.Explor,2003,81：81–98.
[86]戚莹.基于模糊神经网络的森林虫害预测预报的应用研究[D].黑龙江：东北林业大学,2011.
[87]万中,梁文冬,卢宗娟.模糊数的隶属度区间分布函数[J].重庆理工大学学报（自然科学版）,2011,1:107-112.
[88]吴小俊,曹奇英,王士同,陈保香,刘同明.模糊神经系统及其应用研究[J].系统工程,1999,4:61-64.
[89]龚利,王云宽,宋英华,曹志强.基于模糊神经网络的固体粉末自动配料控制系统[J].计量学报,2008,3:284-288.
[90]孙祥,杨子荣,徐大地.基于模糊神经网络的义县萤石成矿预测[R].地质找矿论丛,2011,02:149-152.
[91]姜德义,蒋再文,章永武,任松.模糊神经网络预测公路隧道揭煤突出危险性[J].2000,S1:200-202.
[92]付海涛,王恩德.GIS技术在成矿预测中的应用实例——以杨家杖子矿田为例[J].矿床地质,2005,24(6):684-691.
[93]陈建平,王功文,侯昌波等.基于GIS技术的西南三江北段矿产资源定量预测与评价[J].矿床地质,2005,24(1):15-24.
[94]杜灵通,吕新彪.GIS在矿产资源评价中的应用[J].地质找矿论丛,2003,18(4):275-278.
[95]中华人民共和国国家质量监督检验检疫总局.中华人民共和国国家标准固体矿产地质勘查规范总则[M].北京:中国标准出版社,2002.
[96]宋国耀,张晓华,肖克炎,等.矿产资源潜力评价的若干问题[J].中国地质,1999(8):17-20.
[97]朱海宾.蒙特卡洛模型在矿产资源量预测中的应用[J].地质找矿论丛,2010,1:50-54.