Application of Multi-Fractal Filtering to Extraction of Geochemical Anomalies from Multi-Geochemical Backgrounds: A Case Study of the Southern Section of "Sanjiang Ore-Forming Zone", Southwestern China
-
摘要: 西南“三江”南段位于印度板块与扬子板块的结合带, 是特提斯-喜马拉雅成矿域的重要组成部分.在多旋回的构造岩浆活动过程中, 形成了复杂的成矿地质地球化学背景和丰富的有色金属和贵金属矿产.正是这种复杂的地质地球化学背景, 使得应用地质统计学方法提取致矿地球化学异常遇到了挑战.为此, 多重分形滤波技术被引进, 试图从复杂的多重地球化学背景中有效地分离与成矿有关的Cu、Zn异常.研究表明: (1) 由于研究区存在多重地球化学背景, 克里格方法揭示的Cu、Zn剩余异常通常是受区域线形构造控制的区域异常, 多重分形滤波技术则揭示了与矿化有关的Cu、Zn局部致矿异常; (2) Cu致矿异常的提取表明, 多重分形滤波技术能够克服某些岩性(譬如玄武岩) 引起的Cu高背景的影响, 并从多重地球化学背景中有效提取致矿异常, 包括低背景中的弱异常和隐蔽异常; (3) 运用多重分形滤波技术清晰揭示了一系列Zn原始数据图上没有显示的、规模不等的NNE向Zn矿化异常带, 其中那条穿越整个研究区的Zn矿化异常带含有金顶巨型Pb-Zn矿床和一系列矿床、矿点、矿化点.整个研究工作是基于MORPAS3.0软件系统实现的.Abstract: The southern section of "Sanjiang ore-forming zone" (the western Yunnan terrain), southwestern China, located at the joined belt between the Indian plate and Yangtze plate, is an important area of eastern Tethyan tectonic ore-forming domain. The polycyclic tectono-magmatic activities resulted both in the diversity of ore-forming geochemical background, and in the numerous base metal deposits and precious metal deposits, including a few such world-class deposits as the Jinding Pb-Zn deposit and the Laowangzhai Au deposit within the section. Such a complex geochemical background challenges us with the geostatistical extraction of geochemical anomalies for the mineralization purpose. For this, the multi- fractal filtering (S-A method) technique is introduced to isolate effectively the Cu and Zn anomalies related to mineralization in the complex multi-geochemical background. The following three research results are obtained in this study: (1) In the complicated geochemical background, the residual anomalies of Cu and Zn uncovered by Kriging are typical of the regional anomalies controlled by the regional linear fractural pattern; but the S-A method reveals the local anomalies of Cu and Zn associated with mineralization; (2) The separation of the Cu anomaly from its background demonstrates that the S-A method could effectively identify the high geochemical background of Cu from basalts and extract the Cu anomalies from multi-geochemical background, including weak and hidden Cu anomalies from their low background; (3) A few NNE-trending Zn mineralized and anomalous belts containing the super-large Jinding Pb-Zn deposit and a series of some other Pb-Zn deposits, that are absent at the Zn original data map, are clearly shown at the S-A anomalous map. This discovery presents us with not only new areas for prospecting Cu-Pb-Zn deposits but also new perspectives for the research into new ore-controlling factors, for the regional distribution of Pb-Zn deposits, which may include NNE-trending hidden fractures. All the research work is completed with the help of MORPAS 3.0.
-
西南“三江”南段位于印度板块与扬子板块的结合带, 是特提斯-喜马拉雅成矿域的重要组成部分.在多旋回的构造岩浆活动过程中, 形成了复杂的成矿地质地球化学背景和丰富的有色金属和贵金属矿产(图 1).成矿地质地球化学背景的复杂性和矿床类型的多样性, 尤其是二叠纪峨嵋山玄武岩引起的Cu等元素的高背景(陈永清等, 2003), 使得用传统的数据处理方法, 譬如地质统计学方法提取致矿异常遇到了挑战.这里引进幂率谱-面积分形滤波技术(the power spectrum-area fractal method, 简称S-A) 从其多重地球化学背景中提取Cu、Zn异常, 取得了预期的效果.这种新开发的滤波技术已被成功地应用于加拿大多个地区地球化学、地球物理和遥感致矿异常信息的提取(Cheng et al., 2000; Cheng, 2001, 2004; Xu and Cheng, 2001).该方法的特点是在从复杂的地质矿化总体中提取致矿异常时既考虑矿化地质体幂率谱各向异性尺度不变性特征, 又考虑其奇异性特征.
图 1 西南“三江”南段地质矿产简图(据云南省地质调查院 1∶50 万数字地质图编绘, 2003)Kz.新生代碎屑岩、局部夹火山岩; Mz.中生代碎屑岩、碳酸盐岩、局部夹火山岩; Pz.古生代碳酸盐岩、碎屑岩、夹火山岩; Pt.元古宇变质岩; 1.新近纪-古近纪安山岩; 2.二叠纪玄武岩; 3.中酸性岩体; 4.基性岩体; 5.断层; 6.地质界线; 7.国界; 8.省界; 9.湖泊Fig. 1. Geological map and the distribution of the known mineral deposits in southern region of "Sanjiang ore-forming zone", southwestern China1. 多重分形滤波基本原理
地壳在纵向和横向的非均质性, 其物质结构的各向异性通常具有分形结构和自相似特征(Herzfeld and Overbeck, 1999; 成秋明, 2006), 这种特征在频率域中能够表达为下列幂律关系:
(1) 这里S是能谱密度(energy spectrum density), A是大于能谱密度某一临界值(S0) 的面积.不同的β值在logA (> S) -log (S) 图上能够获取, 这取决于分形滤波器的构置.通常在log-log图上, 所有直线段服从关系式(1).不同的直线段代表了不同的分形关系, 两条直线的交点所对应的横坐标值(能谱密度值) 被视为确定分形滤波器的阈值.借助这些阈值, 可以构造异常滤波器、背景滤波器.在频率域上运用不同的滤波器, 并将其变换到空间域上以实现将异常和背景进行分离.这种分离异常和背景的S-A方法能通过GeoDas (Cheng, 2002) 和MORPAS3.0 (陈永清等, 2006) 软件系统加以实现, 本次研究工作是基于后者完成的.
2. Cu地球化学异常提取
为了对比应用不同的方法提取Cu异常的效果, 首先编制了Cu原始含量图(图 2a).该图表明Cu高含量主要分布在二叠纪峨嵋山玄武岩区域(对比图 2a和图 1).玄武岩中Cu的平均含量为196×10-6, 是世界玄武岩Cu平均含量的2.25倍; 而在保山地块的上寒武统Cu的平均含量仅为25×10-6 (陈永清等, 2003,2005 ).这表明位于玄武岩区的某些大规模的所谓Cu异常是由玄武岩引起的.然后, 在5×5 km2和25×25 km2的滑动窗口内分别使用反距离加权插值法(IDW) 计算插值, 再将插值结果相减得到Cu剩余异常图(图 2b).该图反映了受区域线形构造控制的Cu区域地球化学异常, 因为这些异常的分布趋势与研究区断裂带在空间上具有高度的一致性, 皆呈SN向(研究区西部怒江、澜沧江流域) 和NW-SE向(哀牢山断裂带和红河流域) 分布(对比图 2b和图 1).最后, 应用S-A法获取了lnA (> S) -lnS关系图(图 3).该方法确保了如图 3中所表明的S和A之间的幂律关系, 不同斜率的直线段代表了不同的自相似性特征, 它们通常对应空间域中不同的地球化学背景和异常模式.例如, 图 3中, 用最小二乘法模拟的右边两条线段的交点确立的阈值lnS0=9.05.S < S0通常代表异常, S > S0代表背景.通常, 出于圈定和评价靶区的目的, 我们只对异常感兴趣.通过逆Fourier变换, 并取lnS0=9.05作为阈值绘制Cu异常图(图 4).该图表明, 大多数已知Cu矿床分布在Cu异常区, 同时提供了一些新的Cu异常区(找矿远景区).与图 2a含量相比, 在图 4中, 与Cu矿化有关的弱Cu异常被增强, 图 2a中没有显示的许多与Cu矿化有关的所谓隐蔽Cu异常被揭示; 而在图 2a中由玄武岩引起的所谓Cu的区域性异常, 在图 5a中呈现Cu的高背景.这表明, 多重分形滤波技术能够从多重地球化学背景中有效地提取致矿地球化学异常.
3. Zn地球化学异常提取
类似于研究Cu的程序, 绘制Zn原始含量图(图 5a)、Zn剩余异常图(图 5b)、lnA (> S) -ln (S) 幂律关系图(图 6) 和应用S-A方法获取的Zn异常图(图 7).和Cu剩余异常图一样, 图 5b表明Zn区域异常与Cu区域异常具有类似的空间分布, 皆受SN向和NW-SE向区域大断裂控制.在图 6上, 取lnS0=12.52作为阈值, 应用S-A方法绘制Zn异常分布图(图 7).图 7表明: (1) 已知Zn矿床空间分布与Zn异常具有较高的一致性; (2) 更有意义的是, 新发现一系列NNE向的规模不等的Pb-Zn矿化异常带; 其中最重要的是那条穿越全区的NNE向的Zn矿化异常带, 在该带上, 分布有包括金顶巨型矿床在内的一系列Pb-Zn多金属矿床、矿点和矿化点.Zn矿床和Zn异常具有等距分布、分段集中的特点.这一发现不仅为研究区Pb-Zn多金属矿床的勘查开辟了新的区域, 而且为研究该类矿床区域分布的控制因素(NNE向构造控矿) 提出了新的课题.
4. 结论
成矿作用是地壳局域范围内发生的一种物质和结构上的非均一化过程, 其产物(矿床) 是典型的具有各向异性的异常地质体.大型、超大型矿床通常产在地壳结构和物质组成复杂的区域(Gorelov, 1982), 其地球物理场和地球化学场在时空尺度上都是高度变化的.这就使我们很难用经典的地质统计学去准确刻画这种背景和异常的复杂关系.这是因为传统的地质统计学方法只能描述或模拟其变化性在一个数量级、至多不超过二个数量级的变化范围(Lovejoy et al., 2005), 因此, 我们必须求助谱分析和其他分形以及多重分形模型等强有力工具去客观揭示这种复杂规律.
本次研究表明: (1) 在复杂的多重地球化学背景下, 地质统计学方法揭示的Cu、Zn剩余异常是受区域线形构造控制的区域异常; 而多重分形滤波技术则揭示了与矿化相关的Cu、Zn局部异常, 其异常与已知矿床在空间分布上具有较高的一致性. (2) Cu异常的有效提取表明, 多重分形滤波技术能够克服玄武岩引起的高背景, 从而有效提取不同背景上的异常, 包括低背景区的弱异常和隐蔽异常; 从而为该区Cu矿勘查提供新的靶区. (3) 多重分形滤波技术揭示了一系列规模不等的NNE向的Zn矿化异常带; 这一发现不仅为研究区Pb-Zn多金属矿床的勘查开辟了新的区域, 而且为研究该类矿床区域分布的控制因素(NNE向构造控矿) 提供了新的启示.
-
图 1 西南“三江”南段地质矿产简图(据云南省地质调查院 1∶50 万数字地质图编绘, 2003)
Kz.新生代碎屑岩、局部夹火山岩; Mz.中生代碎屑岩、碳酸盐岩、局部夹火山岩; Pz.古生代碳酸盐岩、碎屑岩、夹火山岩; Pt.元古宇变质岩; 1.新近纪-古近纪安山岩; 2.二叠纪玄武岩; 3.中酸性岩体; 4.基性岩体; 5.断层; 6.地质界线; 7.国界; 8.省界; 9.湖泊
Fig. 1. Geological map and the distribution of the known mineral deposits in southern region of "Sanjiang ore-forming zone", southwestern China
-
[1] Chen, Y. Q., Chen, J. G., Wang, X. Q., 2006. GIS based integrated quantitative assessment of mineral resources. Geological Bulletin of China (in Press, in Chinese with English abstract). [2] Chen, Y. Q., Lu, Y. X., Xia, Q. L., 2005. Geochemical characteristics of the Hetaoping Pb-Zn deposit, Baoshan area, Yunnan Province, and its genetic model and ore prospecting pattern. Geology in China, 32 (1): 90-99 (in Chinese with English abstract). [3] Chen, Y. Q., Xia, Q. L., Liu, H. G., 2003. Geochemical characteristics of Pt-Pd-Cu ore-bearing formations in eastern Yunnan and analysis of their ore potentiality. Geology in China, 30 (3): 225-234 (in Chinese with English abstract). [4] Cheng, Q. M., 2001. Selection of multifractal scaling breaks and separation of geochemical and geophysical anomalies. Journal of China University of Geosciences, 12 (1): 54-59. [5] Cheng, Q. M., 2002. Geodata analysis system (Geo Das) for mineral exploration: User 'guide and exercise manual. Material for thetraining workshop on Geo Das held at York University, November1-3.289 (http://www.gisworld.org/geodas). [6] Cheng, Q. M., 2004. A new model for quantifying anisotropic scale invariance and for decomposition of mixing patterns. Mathematical Geology, 36 (3): 345-360. doi: 10.1023/B:MATG.0000028441.62108.8a [7] Cheng, Q. M., 2006. Singularity-generalized self-similarity-fractal spectrum (3S). Earth Science—Journal of China University of Geosciences, 31 (3): 337-348 (in Chinese with English abstract). [8] Cheng, Q. M., Xu, Y. G., Grunsky, E., 2000. Integrated spatial and spectrum analysis for geochemical anomaly separation. Nature Resources Research, 9 (1): 43-56. doi: 10.1023/A:1010109829861 [9] Gorelov, D. A., 1982. Quantitative characteristics of geologic anomalies in assessing ore capacity. International Geology Review, 24 (4): 457-465. doi: 10.1080/00206818209452430 [10] Herzfeld, U. C., Overbeck, C., 1999. Analysis and simulation of scale-dependent fractal surfaces with application to seafloor morphology. Computer & Geosciences, 25 (9): 963-978. [11] Lovejoy, S., Schertzer, D., Gagnon, J. S., et al., 2005. Multifractal simulation of the earth's surface and interior: Anisotropic singularities and morphology. In: Cheng, Q. M., Bonham-Carter, G., eds., Proceeding of IAMG: GIS and spatial analysis, 1: 37-54. [12] Xu, Y. G., Cheng, Q. M., 2001. A multifractal filter technique for geochemical data analysis from Nova Scotia, Canada. Geochemistry: Exploration, Analysis and Environment, 1 (2): 1-12. [13] 陈永清, 陈建国, 汪新庆, 2006. 基于GIS的矿产资源综合定量评价. 地质通报(待刊). [14] 陈永清, 卢映祥, 夏庆霖, 2005. 云南保山核桃坪铅锌矿床地球化学特征及其成矿模式与找矿模型. 中国地质, 32 (1): 90-99. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI200501012.htm [15] 陈永清, 夏庆霖, 刘红光, 2003. 滇东Pt-Pd-Cu含矿建造地球化学特征及其含矿性分析. 中国地质, 30 (3): 225-234. doi: 10.3969/j.issn.1000-3657.2003.03.001 [16] 成秋明, 2006. 非线性成矿预测理论: 多重分形奇异性-广义自相似性-分形谱系模型与方法. 地球科学——中国地质大学学报, 31 (3): 337-348. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200603008.htm -