地球化学数据处理软件系统设计
|
摘要
地球化学元素分布规律的研究是揭示元素矿化富集及空间变化规律的重要途径之一。勘查地球化学数据的统计特征常被用来描述和刻画地球化学元素的分布规律，然而人们早已注意到普通的统计方法并不考虑样品的空间分布和统计特征随空间度量尺度的变化性，而分形理论能有效地克服统计方法的不足。随着分形方法在地球化学数据处理中应用日益广泛，将各种方法集成到软件系统中对生产和科学研究都有很大的意义。
笔者总结了国内外将分形理论应用到地球化学数据处理的几个大的方向，在Windows平台下用C++编程语言开发出了一套集成多种分形算法的软件系统。软件系统有五大核心模块：数据预处理；分维数计算；地球化学异常下限的求取；数据网格化；多重分形谱。
在开发出的软件平台上，笔者制作了14种元素富集模型的多重分形谱图，对谱图的各种指标做了详细分析，并针对多重分形谱图的特点制订了两种新的分析指标：富集效率指数(Enrichment Efficiency Index，EEI)和富集强度指数(Enrichment Intensity Index，EII)。还用南岭地区1∶5万、1∶20万勘查地球化学数据、江西大吉山地区1∶5万勘查地球化学数据共三十多种元素的多重分形谱图和理论模型对比，根据多重分形谱的各指标总结了三个地区的元素富集规律，效果良好。
从软件使用及结果分析可以看出：多重分形谱能较好地描述地球化学数据的分布特征；在开发出的软件平台上能方便、快捷地制作出多重分形谱图、频率直方图和各种分维散点图，以及快速网格化数据和确定异常下限等，有较好的实用价值。
The research of the law of geochemical elements distribution is important to realize the law of enrichment and spatial variation of elements. The statistic character of geochemical data is often used to describe the law of geochemical elements distribution. However, the general statistic methods can not distinguish the spatial distribution of the elements and the variation of statistic character with different scale, but the fractal methods can overcome those defects. With the increasing number of fractal methods being used in the geochemical data processing, a software system which is integrated many useful fractal methods is definitely significant to scientific research and practice.
In this paper, the main fractal methods in geochemical data processing are generalized, and the specific design of a software system which is designed by C++ programming language in Windows operation system is introduced. There are five cores in this software: data preprocessing, fractal dimension calculation, geochemical anomalous threshold, grid data and multifractal spectrum.
Fourteen multifractal spectrum graphs about element enrichment model are produced by the software introduced above, and every parameter of the multifractal spectrum is analyzed specifically, and two new indexes EEI(Enrichment Efficiency Index) and EII(Enrichment Intensity Index) are invented from the characters of multifractal spectrum graph. The multifractal spectrums of geochemical exploration data in Nanling region with the scale of 1:50000 and 1:200000 and in Dajishan mining district in Jiangxi province with the scale of 1:50000 are calculated and compared with the spectrums of the numerical models. In the end, the law of element enrichment in three regions is summarized from the indexes of the multifractal spectrum, and the outcome is satisfying.
To sum up, the multifractal spectrum can describe the distribution character of geochemical data finely. By using the software introduced in this paper, we can produce multifractal spectrum graph, the frequency distribution graph and disperse dot graph of fractal methods and obtain grid data and geochemical anomalous
笔者总结了国内外将分形理论应用到地球化学数据处理的几个大的方向，在Windows平台下用C++编程语言开发出了一套集成多种分形算法的软件系统。软件系统有五大核心模块：数据预处理；分维数计算；地球化学异常下限的求取；数据网格化；多重分形谱。
在开发出的软件平台上，笔者制作了14种元素富集模型的多重分形谱图，对谱图的各种指标做了详细分析，并针对多重分形谱图的特点制订了两种新的分析指标：富集效率指数(Enrichment Efficiency Index，EEI)和富集强度指数(Enrichment Intensity Index，EII)。还用南岭地区1∶5万、1∶20万勘查地球化学数据、江西大吉山地区1∶5万勘查地球化学数据共三十多种元素的多重分形谱图和理论模型对比，根据多重分形谱的各指标总结了三个地区的元素富集规律，效果良好。
从软件使用及结果分析可以看出：多重分形谱能较好地描述地球化学数据的分布特征；在开发出的软件平台上能方便、快捷地制作出多重分形谱图、频率直方图和各种分维散点图，以及快速网格化数据和确定异常下限等，有较好的实用价值。
The research of the law of geochemical elements distribution is important to realize the law of enrichment and spatial variation of elements. The statistic character of geochemical data is often used to describe the law of geochemical elements distribution. However, the general statistic methods can not distinguish the spatial distribution of the elements and the variation of statistic character with different scale, but the fractal methods can overcome those defects. With the increasing number of fractal methods being used in the geochemical data processing, a software system which is integrated many useful fractal methods is definitely significant to scientific research and practice.
In this paper, the main fractal methods in geochemical data processing are generalized, and the specific design of a software system which is designed by C++ programming language in Windows operation system is introduced. There are five cores in this software: data preprocessing, fractal dimension calculation, geochemical anomalous threshold, grid data and multifractal spectrum.
Fourteen multifractal spectrum graphs about element enrichment model are produced by the software introduced above, and every parameter of the multifractal spectrum is analyzed specifically, and two new indexes EEI(Enrichment Efficiency Index) and EII(Enrichment Intensity Index) are invented from the characters of multifractal spectrum graph. The multifractal spectrums of geochemical exploration data in Nanling region with the scale of 1:50000 and 1:200000 and in Dajishan mining district in Jiangxi province with the scale of 1:50000 are calculated and compared with the spectrums of the numerical models. In the end, the law of element enrichment in three regions is summarized from the indexes of the multifractal spectrum, and the outcome is satisfying.
To sum up, the multifractal spectrum can describe the distribution character of geochemical data finely. By using the software introduced in this paper, we can produce multifractal spectrum graph, the frequency distribution graph and disperse dot graph of fractal methods and obtain grid data and geochemical anomalous
引文
[1] 於崇文．固体地球系统的复杂性与自组织临界性．地学前缘，1998，5(3)：159～174
[2] 於崇文，岑况，鲍征宇等．成矿作用动力学．北京：地质出版社，1998
[3] 於崇文．大型矿床和成矿区(带)在混沌边缘．地学前缘，1999，6(1)：85～102
[4] 於崇文．地质系统的复杂性——地质科学的基本问题(Ⅰ)．地球科学——中国地质大学学报，2002，27(5)：509～519
[5] 於崇文．地质系统的复杂性——地质科学的基本问题(Ⅱ)．地球科学——中国地质大学学报，2003，28(1)：31～40
[6] 孟宪伟，余先川．地球化学场分形景观与地球化学分级预测．物探与化探，1994，18(5)：393～397
[7] 孟宪伟，张晓华．多标度分形与地球化学场分解．地质与勘探，1996，32(4)：47～49
[8] 李长江，麻土华．矿产勘查中的分形、混沌与ANN．北京：地质出版社，1999
[9] 成秋明．多维分形理论和地球化学元素分布规律．地球科学——中国地质大学学报，2000，25(3)：311～318
[10] 施俊法，向运川，王春宁．区域地球化学异常空间分形结构及其意义——以浙江省诸暨地区区域地球化学数据为例．矿物学报，2000，20(1)：68～72
[11] 施俊法．浙江省诸暨地区元素地球化学分布与标度律．地球科学——中国地质大学学报，2001，26(2)：167～171
[12] 成秋明．多重分形与地质统计学方法用于勘查地球化学异常空间结构和奇异性分析．地球科学——中国地质大学学报，2000，26(2)：161～166
[13] 龚庆杰，张德会，韩东昱．一种确定地球化学异常下限的简便方法．地质地球化学，2001，29(3)：215～220
[14] 谢淑云，鲍征宇．地球化学场的连续多重分形模式．地球化学，2002，31(2)：191～200
[15] 李庆谋，成秋明．分形奇异(特征)值分解方法与地球物理和地球化学异常重建．地球科学——中国地质大学学报，2004，29(1)：109～118
[16] 韩东昱，龚庆杰等．区域化探数据处理的几种分形方法．地质通报 2004，23(7)：714～719
[17] 邹林，彭省临等．青海阿尔茨托山地区地球化学异常场的多重分形研究．中国地质，2004，31(4)：436～441
[18] 申维．n维自仿射分形及其在地球化学中的应用．地质论评，2005，51(2)：208～211
[19] 谢淑云，鲍征宇等．新疆塔里木盆地艾协克-桑塔木南地区油气地球化学指标的多重分形性分析．地质论评 2005，51(1)：107～112
[20] 杨茂森，黎清华，杨海巍．分形方法在地球化学异常分析中的运用研究——以胶东矿集区为例．地球科学进展，2005，20(7)：809～814
[21] 倪秀静，李健，韩泽华．用MATLAB实现成矿区化学元素的第L次分形插值逼近曲面．物探化探计算技术，2005，27(2)：150～153
[22] 李随民，姚书振．基于MapGIS的分形方法确定化探异常．地球学报，2005，26(2)：187～190
[23] 陈志海，戴勇，郎兆新．缝洞性油藏油井产量变化分形特征及其意义．石油勘探与开发，2006，33(1)：95～98
[24] 李庆谋，刘少华，梁光和．ArcGeophysics1.0∶GIS环境下综合处理、解释地球物理与地质、地球化学、遥感等资料的组件(COM)包与应用实例．地球物理学进展，2006，21(1)：194～201
[25] 成秋明．非线性成矿预测理论：多重分形奇异性-广义自相似性-分形谱系模型与方法．地球科学—中国地质大学学报，2006，31(3)：337～348
[26] Qiu ming Cheng, Agterberg F P, Ballantyne S B. The separation of geochemical anomalies from background by fractal methods. Journal of Geochemical Exploration, 1994, 51:109～130
[27] Qiu ming Cheng, Agterberg F P ,etc.A spatial analysis method for geochemical anomaly separation. Journal of Geochemical Exploration, 1996, 56:183～195
[28] Qiu ming Cheng.Spatial and scaling modelling for geochemical anomaly separation. Journal of Geochemical exploration,1999,65:175～194
[29] A.Lima, B.De Vivo,etc. Multifractal IDW interpolation and fractal filtering method in environmental studies: an application on regional stream sediments of (Italy), Campania region. Applied GeoChemistry, 2003,18:1853～1865
[30] Changjiang Li,Tuhua Ma, Junfa Shi.Application of a fractal method relating concentrations and distances for separation of geochemical anomaliesfrom background. Journal of Geochemical Exploration, 2003, 77:167～175
[31] 孙霞，傅竹西，吴自勤．薄膜生长的多重分形谱的计算．计算物理，2001，18(3)：247～252
[32] Frederik P．Agterberg．地球化学图纹理的多重分形模拟．地球科学—中国地质大学学报，2001，26(2)：142～151
[33] 姜志强．R~2空间上分布数据的多重分形维数谱计算和分形特征提取．小型微型计算机系统，2004，25(10)：1844～1846
[34] 李平，胡可乐，汪秉宏．多重分形谱在材料分析中的应用研究．南京航空航天大学学报，2004，36(1)：77～81
[35] 刘晓冬，徐景祯，王新峰．断裂体系多重分形谱特征及其应用．石油地球物理勘探，1999，34(3)：285～302
[36] 王金安，谢和平．岩石断裂面的各向异性分形和多重分形研究．岩土工程学报，1998，20(6)：16～21
[37] Qiuming Cheng.The gliding box method ofr multifractal modeling.Computer&Geoscienees, 1999(25):1073～1079
[38] Qiuming Cheng.Multifraetality and spatial statistics.Computer&Geosciences, 1999(25):949～961
[39] 周宇欢．多分形理论及其在复杂系统中的应用：[硕士学位论文]．南京：中国人民解放军理工大学，2006
[40] 申维．分形混沌与矿产预测．北京：地质出版社，1999
[41] 李方林，鲍征宇，裴韬，等．地球化学空间数据处理原理及软件系统．地球科学——中国地质大学学报，1999，24(3)：316～319
[42] 沈步明，沈远超．新疆某金矿的分数维特征及其地质意义．中国科学(B辑)，1993，23(3)：297～302
[43] 程小久，卢建杭，宋亮明．铅锌品位分维D值的意义和计算程序．地质与勘探，1994，30(5)：32～35
[44] 王新社，郑亚东，杨崇辉，等．用动态重结晶石英颗粒的分形确定变形温度及应变速率．岩石矿物学杂志，2001，20(1)：36～40
[45] Keith C.Clarke. Computation of the fractal dimension of topographic surfaces using the triangular prism surface area, Computer & Geosciences.1986, 12(5):713～722
[46] 王喜生，杨永华．地球化学分散模式分维计算方法的探讨．地质地球化学，2000，28(3)：75～81
[47] 王崇云．地球化学找矿基础．北京：地质出版社，1999
[48] 谢淑云，鲍征宇．多重分形方法在金属成矿潜力评价中的应用．成都理工大学学报(自然科学版)，2004，30(1)：28～33
[49] 韩吟文，马振东．地球化学．北京：地质出版社，2003
[2] 於崇文，岑况，鲍征宇等．成矿作用动力学．北京：地质出版社，1998
[3] 於崇文．大型矿床和成矿区(带)在混沌边缘．地学前缘，1999，6(1)：85～102
[4] 於崇文．地质系统的复杂性——地质科学的基本问题(Ⅰ)．地球科学——中国地质大学学报，2002，27(5)：509～519
[5] 於崇文．地质系统的复杂性——地质科学的基本问题(Ⅱ)．地球科学——中国地质大学学报，2003，28(1)：31～40
[6] 孟宪伟，余先川．地球化学场分形景观与地球化学分级预测．物探与化探，1994，18(5)：393～397
[7] 孟宪伟，张晓华．多标度分形与地球化学场分解．地质与勘探，1996，32(4)：47～49
[8] 李长江，麻土华．矿产勘查中的分形、混沌与ANN．北京：地质出版社，1999
[9] 成秋明．多维分形理论和地球化学元素分布规律．地球科学——中国地质大学学报，2000，25(3)：311～318
[10] 施俊法，向运川，王春宁．区域地球化学异常空间分形结构及其意义——以浙江省诸暨地区区域地球化学数据为例．矿物学报，2000，20(1)：68～72
[11] 施俊法．浙江省诸暨地区元素地球化学分布与标度律．地球科学——中国地质大学学报，2001，26(2)：167～171
[12] 成秋明．多重分形与地质统计学方法用于勘查地球化学异常空间结构和奇异性分析．地球科学——中国地质大学学报，2000，26(2)：161～166
[13] 龚庆杰，张德会，韩东昱．一种确定地球化学异常下限的简便方法．地质地球化学，2001，29(3)：215～220
[14] 谢淑云，鲍征宇．地球化学场的连续多重分形模式．地球化学，2002，31(2)：191～200
[15] 李庆谋，成秋明．分形奇异(特征)值分解方法与地球物理和地球化学异常重建．地球科学——中国地质大学学报，2004，29(1)：109～118
[16] 韩东昱，龚庆杰等．区域化探数据处理的几种分形方法．地质通报 2004，23(7)：714～719
[17] 邹林，彭省临等．青海阿尔茨托山地区地球化学异常场的多重分形研究．中国地质，2004，31(4)：436～441
[18] 申维．n维自仿射分形及其在地球化学中的应用．地质论评，2005，51(2)：208～211
[19] 谢淑云，鲍征宇等．新疆塔里木盆地艾协克-桑塔木南地区油气地球化学指标的多重分形性分析．地质论评 2005，51(1)：107～112
[20] 杨茂森，黎清华，杨海巍．分形方法在地球化学异常分析中的运用研究——以胶东矿集区为例．地球科学进展，2005，20(7)：809～814
[21] 倪秀静，李健，韩泽华．用MATLAB实现成矿区化学元素的第L次分形插值逼近曲面．物探化探计算技术，2005，27(2)：150～153
[22] 李随民，姚书振．基于MapGIS的分形方法确定化探异常．地球学报，2005，26(2)：187～190
[23] 陈志海，戴勇，郎兆新．缝洞性油藏油井产量变化分形特征及其意义．石油勘探与开发，2006，33(1)：95～98
[24] 李庆谋，刘少华，梁光和．ArcGeophysics1.0∶GIS环境下综合处理、解释地球物理与地质、地球化学、遥感等资料的组件(COM)包与应用实例．地球物理学进展，2006，21(1)：194～201
[25] 成秋明．非线性成矿预测理论：多重分形奇异性-广义自相似性-分形谱系模型与方法．地球科学—中国地质大学学报，2006，31(3)：337～348
[26] Qiu ming Cheng, Agterberg F P, Ballantyne S B. The separation of geochemical anomalies from background by fractal methods. Journal of Geochemical Exploration, 1994, 51:109～130
[27] Qiu ming Cheng, Agterberg F P ,etc.A spatial analysis method for geochemical anomaly separation. Journal of Geochemical Exploration, 1996, 56:183～195
[28] Qiu ming Cheng.Spatial and scaling modelling for geochemical anomaly separation. Journal of Geochemical exploration,1999,65:175～194
[29] A.Lima, B.De Vivo,etc. Multifractal IDW interpolation and fractal filtering method in environmental studies: an application on regional stream sediments of (Italy), Campania region. Applied GeoChemistry, 2003,18:1853～1865
[30] Changjiang Li,Tuhua Ma, Junfa Shi.Application of a fractal method relating concentrations and distances for separation of geochemical anomaliesfrom background. Journal of Geochemical Exploration, 2003, 77:167～175
[31] 孙霞，傅竹西，吴自勤．薄膜生长的多重分形谱的计算．计算物理，2001，18(3)：247～252
[32] Frederik P．Agterberg．地球化学图纹理的多重分形模拟．地球科学—中国地质大学学报，2001，26(2)：142～151
[33] 姜志强．R~2空间上分布数据的多重分形维数谱计算和分形特征提取．小型微型计算机系统，2004，25(10)：1844～1846
[34] 李平，胡可乐，汪秉宏．多重分形谱在材料分析中的应用研究．南京航空航天大学学报，2004，36(1)：77～81
[35] 刘晓冬，徐景祯，王新峰．断裂体系多重分形谱特征及其应用．石油地球物理勘探，1999，34(3)：285～302
[36] 王金安，谢和平．岩石断裂面的各向异性分形和多重分形研究．岩土工程学报，1998，20(6)：16～21
[37] Qiuming Cheng.The gliding box method ofr multifractal modeling.Computer&Geoscienees, 1999(25):1073～1079
[38] Qiuming Cheng.Multifraetality and spatial statistics.Computer&Geosciences, 1999(25):949～961
[39] 周宇欢．多分形理论及其在复杂系统中的应用：[硕士学位论文]．南京：中国人民解放军理工大学，2006
[40] 申维．分形混沌与矿产预测．北京：地质出版社，1999
[41] 李方林，鲍征宇，裴韬，等．地球化学空间数据处理原理及软件系统．地球科学——中国地质大学学报，1999，24(3)：316～319
[42] 沈步明，沈远超．新疆某金矿的分数维特征及其地质意义．中国科学(B辑)，1993，23(3)：297～302
[43] 程小久，卢建杭，宋亮明．铅锌品位分维D值的意义和计算程序．地质与勘探，1994，30(5)：32～35
[44] 王新社，郑亚东，杨崇辉，等．用动态重结晶石英颗粒的分形确定变形温度及应变速率．岩石矿物学杂志，2001，20(1)：36～40
[45] Keith C.Clarke. Computation of the fractal dimension of topographic surfaces using the triangular prism surface area, Computer & Geosciences.1986, 12(5):713～722
[46] 王喜生，杨永华．地球化学分散模式分维计算方法的探讨．地质地球化学，2000，28(3)：75～81
[47] 王崇云．地球化学找矿基础．北京：地质出版社，1999
[48] 谢淑云，鲍征宇．多重分形方法在金属成矿潜力评价中的应用．成都理工大学学报(自然科学版)，2004，30(1)：28～33
[49] 韩吟文，马振东．地球化学．北京：地质出版社，2003