遥感地质论文

近40年来，航天遥感对地观测的宏伟发现和贡献就是直观形象地揭示并证实了地球表面到一定深度处存在不同类型、尺度规模、新老世代重叠、稠密交织、成群布列的圆环形地质结构。滤却了由于风水外动力作用、浮土、植被的干扰之后，使人们吃惊地发现，地球也如同月球、水星、火星一样满目尽是环形山的地质景观。其实表面满布环形结构的还有很多固体行星和卫星，不少表面为液态、气态的星体，同样也还存在流体旋涡，看来圆涌结构是太空星体的共有的动力物质形态。成因众说纷纭，如星体内热核爆炸圆筒涌射、火山爆发、小行星撞击等。

就地球地质科学而论，基于大量的遥感数字、影像信息，已有的地质矿产资料，现代地质学知识，参照超高能喷射流变力学原理、地下核试验的地质动力特征，对地球的圆环形喷涌结构进行了长期系统的分析对比研究，包括对圆环(筒)状构造的地质学研究、专题编图，青藏高原、长江中下游地带、哈密地区、白银厂、中条山等矿区的野外验证调查，于1977年提出“地球的圆涌构造与深源高能强爆热动力学”的新理论假设。认为这些源于地幔、下地壳岩石圈的竹笋形垂直筒状、饼状地质体，可划分为根基热核动力部，修长的管颈部，顶盖头部，以强大的热力、压力和气液固混合态，按流体力学规律，喷涌到地壳-岩石圈的不同深度，其中有部分头盖出露到地表，保有现实热动力源的，在海水中、大气圈中仍保持有明显的热动力晕梢，如同休眠火山。地球早期形成圆饼状原始地壳和结晶基底，后期根源逐步加深，直径越来越小，射喷动力越来越强，包括根基较浅，规模较小的金伯利型岩筒、强爆角砾岩筒、基性-超基性岩株、碱性和碳酸岩小岩颈，圆涌的定义范围，要比地幔柱(Mantle plume)更深广些。这一群地质体密集林立，新老包容，构成比地壳总质量大50多倍的岩石圈和上、下地幔，巨大的物质能量构成了地球的地质主力层。其注入形式为高能态的固液气混融的喷射流体状态，而不是中能态的固塑状态，更不属低能区的表浅固体状态。多数产生于急剧热核强爆灾变，是否有缓慢形成的，则有待考究。形成期内它们在岩石圈-地幔中吞蚀一切，自形程度和热力学独立性很强，根基部能量耗尽后停止喷涌，失去活力的躯体固化、消亡。不少个体长期处于休眠状态，后期还持续热动力效能。这类地质体对周边特别是顶盖部地质形态有极强的统领控制作用，即主力地质体，也是解析受控地质环境的强力核心。

这些遥感发现的圆形结构，多数在已有地质图件、资料中都能找到一定反映。地球物理的磁力场、重力场、电磁场、热力场、放射性场、地球化学场等资料中也有对应的显示。

深源强热动力的矿产聚集过程也很特殊，超高的热力可将地内根基部几立方千米，几十乃至几百立方千米的岩石充分熔融、析离、聚合，在强大的喷射动力推动下，运移到表浅部的特定地段，一次或多次地形成元素高度浓集，有用矿物高纯度的富矿，形成大矿和叠加矿床，如同一个高能的热分馏烟囱出矿口，众矿喷涌的“宝瓶口”，矿床成串的“烟囱”、百宝聚集的矿产树、矿产“百果圣诞树”。矿床可能是固态、液态和气态并存的多态系列结构，其中包括石油天然气的深源无机聚合，地下热能、深层地下水的上涌等。

无数实例证明，越是大型、超大型的多组分的富矿越与这种圆涌结构的特定部位高相关，而且成矿模式几何对称、对应性极强，预测的可靠性很高。

世界上很多知名的大矿、富矿都和特定的圆环形结构密切相关，与地内极高热强爆有关，与隐爆地域相关，例如加拿大的肖德贝里镍铜矿，澳大利亚的蒙特艾萨铅锌矿，奥林匹克坝铀矿，朝鲜的检德铜矿，中国的白云鄂博稀土-铁矿，中条山、白银厂的铜矿，长江中下游的铜铁矿床，金川白家咀子的镍铜铂矿，大庆油气田等。

环境灾害也发生在现活动的强热动力圆涌结构顶部，如百慕大海空罹难区，1908年的通古斯大爆炸，1662年的北京王恭厂大爆炸，印度的死丘事件，1556年的华县地震，三河-平谷地震，山丹地震，唐山大地震，印度尼西亚地震海啸，2008年5月的四川汶川地震，多次火山喷发，大型、特大型的崩塌、滑坡、泥石流，每年多次的台风，异常气象灾难等。

在后来的研究中，我们还陆续推衍出:岩石圈-水圈-大气圈-生态环境圈强热动力共驱链;矿产-环境-灾害垂直强热动力贯通协动链;深源强热动力成矿系列;石油-天然气深源无机聚合;地下热能、地内深水资源的上涌;地气控制天气，生物-人类兴衰背景强热动力谱，地球与天文背景的热共振谱等次一级理论，在这条理论链上有许多待探索的无穷奥秘。

30年来，我们结合遥感地质的推广，以培训、宣讲、座谈、会议论文交流、专题讨论咨询、对口探索等方式向社会和地质界宣传、介绍了新的深源高能强爆热动力学理论的科学原理和应用技巧，不少单位结合自己项目进行了试探性应用，取得了很好的效果，例如河北省、云南省的地质工作者已应用本理论在进行矿床研究预测。

但是理论推广和提升工作却不尽人意，可谓鸡鸣早而起舞缓慢。理论假说的提出也已30多年了，当时可谓很先锋了，至少尚未发现类似完备的说法，但是由于种种原因，一直未设立理论研究专题。

现在国内外地质界都在探索类似的地幔柱问题，但常规地质知识、物化探技术和深部钻探尚无法攻克这一难关。从研究策略上讲，以遥感圆环形结构为切入点，组织遥感、地质、物化探专业的多学科攻关项目，尚不失为一条快捷途径。围绕矿产开发、环境变迁、灾害研究等实用问题，将圆环形地质结构与已知地质矿产资料系统对比分析，提出一批预测验证方案，动员有关勘探、环保、防灾部门经过10年左右的探索，争取矿产资源储量翻番，重点解决环境灾害防治任务，还是有可能的。同时在理论研究上将得到进一步的充实提高。笔者深信本理论一定能快速发展起来，获得伟绩奇效，成为近期深部地质构造的统领性理论，在地质资源开发利用、生态环境优化、防灾减害，推动地质理论、勘查科技发展中为人类文明作出新贡献，不要再使原本创新的地质科学理论生长点继续停滞萎缩，错过时机，丧失机遇，给地质科学新发展带来损失。

(内部资料，2007-03-20)

［美］阿拉斯加卡特迈火山

杨武年朱章森

（成都理工大学遥感与GIS研究所，成都610059）

摘要遥感图像显现的地质构造形迹是地史以来地壳运动的综合结果，包含有历次构造变动和有关矿化的特征标识信息，因此，遥感信息场分层解析是无模型矿床预测中有效地进行成矿信息分离的重要手段。无模型矿床预测法的理论依据是，矿床的形成是多种地质因素发生变异的综合结果，特殊地质环境中多期次构造运动复合叠加区具备形成大型、超大型矿床的条件，巨大矿源体的存在势必导致其与周围环境（背景）在物质成分、结构构造等方面的显著差异，包括诸如成矿岩石组合、构造式样、矿化蚀变以及地球物理、地球化学场甚至由深部矿源异常体引起其上方地球生物圈、大气圈等发生变异的信息特征。因而该法不需要事先有已知有矿模型或已知无矿模型作为预测的基础，而以分解和提取上述致矿异常信息建立控矿地质体的自然源模型作为科学找矿的重要途径。由于此法不需要事先有已知有矿模型或已知无矿模型作为预测的基础，因此它不仅对研究程度较低的新区矿床预测，而且对研究程度较高的老区寻找新类型和大型、超大型矿床都具有重大理论意义和实用价值。

关键词遥感地质数学地质遥感信息场分层解析遥感图像分层解析法无模型矿床预测法

1引言

当前国内外成矿预测研究中，传统的基于相似类比理论的矿床模型预测法，通过数十年的实践日趋完善，在地质找矿中取得显著成效。然而这类方法的应用前提是必须有一定数量的已知有矿模型和无矿模型作为类比的基础，因此，它受资料水平、已知矿床模型及研究人员水平的限制，主要适用于地质研究程度较高、矿床模型易于建立的地区；而且因受相似类比理论限制，它只能找到与已知矿床类同的具相当规模的矿床，而难于发现新类型矿床或点状大型、超大型矿床。因此，在找矿难度日益增大的情况下，如何使多学科结合，应用新技术新方法和新思维开拓新区，寻找隐伏矿床和新类型矿床，特别是大型、超大型矿床，已成为现在和今后地质找矿工作中急需解决的重大理论和方法问题。为此，我们通过长期的找矿实践，提出了在求异理论和求同理论指导下的无模型矿床预测法的新思路[1,2]，后进一步发展，提出了遥感信息场分层解析与无模型矿床预测的方法理论[3]，经在云南、四川、新疆等地找矿实践，取得明显效果[1～5.11～14]。

2理论依据

遥感图像显现的地质构造形迹是地史以来地壳运动的综合结果，包含有历次构造变动和有关矿化的特征标识信息，因此，遥感信息场分层解析是有效地进行成矿信息分离的重要手段[2～6]。无模型矿床预测法的理论依据是，矿床的形成是多种地质因素发生变异的综合结果，特殊地质环境中多期次构造运动复合叠加区具备形成大型、超大型矿床的条件，巨大矿源体的存在势必导致其与周围环境（背景）在物质成分、结构构造等方面的显著差异[15]，包括诸如成矿岩石组合、构造式样、矿化蚀变以及地球物理、地球化学场甚至由深部矿源异常体引起其上方地球生物圈、大气圈等发生变异的信息特征。因此，应用无模型矿床预测法不需要事先有已知有矿模型或已知无矿模型作为预测的基础，而以分解和提取上述致矿异常信息、建立控矿地体的自然源模型作为科学找矿的重要途径。由于此法不需要事先有已知有矿模型或已知无矿模型作为预测的基础，所以它不仅对研究程度较低的新区矿床预测，而且对研究程度较高的老区寻找新类型和大型、超大型矿床都具重大理论意义和实用价值。

3方法技术

3.1遥感图像处理和信息提取

遥感图像是地面景观物体按一定比例缩小了的立体模型，地质构造形迹及地质现象的总体和个体的地表几何形态（纹形图案）以及物理特征（电磁波辐射特征）被真实、客观、连续、全面地记录，具有高度的概括性；同时含有大量地下一定深度隐伏地质构造的信息特征。遥感图像上的这些信息特征（解译标志）反映的是地物在内外动力作用过程中，在一定地质、地理条件下物质成分、结构构造、物理性质等方面的差异。隐伏构造信息则通过地壳的机械变形以及地表的地球物理场和地球化学场改造乃至生物圈、大气圈发生异常等方式而显示出来。而且，多波段遥感图像信息特征中除可见光波谱部分外，许多波段的信息特征位于人的肉眼光敏区以外，如红外、微波图像信息等。这类图像揭示了大量肉眼看不见的地物信息特征，大大增加了鉴别地物属性特征的信息量。因此，采用遥感图像处理能有效地提取成矿信息，尤其是提取那些微弱的矿化信息，通过地质解译、分析构造形迹特征及空间分布规律和应力状态，不但真实、客观，而且克服了常规地质方法有时由于点线观测的局限性，大大开拓视野，获得连续、系统、大量的信息特征，有助于将破裂系统与区域构造变形乃至地质建造等有机地联系起来进行深入的分析研究，得出与客观实际相吻合的结论。

3.2遥感构造信息场分“层”解析与致矿异常信息提取

遥感图像显现的地质构造形迹是地史以来地壳运动的综合结果，包含有历次构造变动和有关矿化的标识信息特征。因此，如何通过构造“分层”解析和应用数学方法去描述，实现线性体场的演译和分解，研究线性体场与其它地质异常的关系，找出其控岩控矿规律和提取矿化信息也就成了线性体场研究和成矿预测的主要内容。需要指出的是，以前大多数研究者注意到了线性构造的规模大小及空间尺度在构造解析和成矿预测中的重要性，但并未形成有效解决问题的方法。特别是在对线性构造量化处理时往往主次不分，把不同规模尺度和不同构造层次和期次的线性、环形构造混为一体进行处理，以致不仅使问题复杂化，而且很难确切找到线性、环形构造与相关地质体及有关矿产的内在联系，所得结论令人质疑。由控矿理论可知，特殊地质环境多期次构造复合叠加区具备形成大型、超大型矿床的条件，而这些地区线性体场的研究与分解，各期次构造变形场的建立及其构造应力场反演，在一定程度上是成矿预测的关键。因为构造应力场不仅控制着地壳岩石的机械变形、构造的成生演化及其组合形式，而且在一定程度上控制着沉积建造、岩浆活动及有关矿床的形成与分布。笔者经多年的科研攻关，根据遥感图像信息量大，对地质构造形迹显示的宏观性、信息连续性以及遥感大节理与宏观构造具有成生联系的性质；在理论上提出遥感大节理总是在野外构造小节理优势方位上发育发展的新论点，并提出了将宏观构造和相对微观构造（特别是遥感横张大节理）有机结合进行地质构造的“分层”解析，根据构造变形场特征反演应力场的新方法——“遥感图像分层解析法”及其相应理论[6,7]。该方法经在不同构造复合区试用，取得显著效果[5～7,10～13]。

例如，笔者在四川、贵州等地的构造复合变形区研究中，应用“遥感图像分层解析法”，通过遥感图像处理和褶皱断裂特别是横张大节理系统分层解析，结合野外调研，建立了研究区多期次构造变形场并反演了各时期构造应力场[10～12]，获得了与众不同的新认识。

笔者在黔西郎岱地区构造解析和成矿预测中[12]，应用航、卫片选择性地解译了该区主要褶皱及与其有关的横张大节理系统，结合野外调研，对该区构造变形特征进行了深入分析。根据该区宏观构造与横张大节理空间分布规律等特征建立了郎岱三角形构造格局的变形场，并据此反演得出该区燕山期在三角形边界条件（大断裂）控制下三边同时受力的联合应力场。其特点是应力和变形强度自三角形边部和顶点向三角形中心呈弧形递减，岩石变形严格按所得应力网络发生，在不同部位形成不同量级的褶皱并伴生相应的节理；三角形边部和顶点褶皱变形强烈，扩展到三角形中部褶皱趋于平缓。九条不同方位的褶皱相互嵌接，总体联系，局部复合，协调统一，控制着本区不同类型的矿床形成与分布。例如，在该三角形构造的北西和南西两个角点附近，集中分布了一系列铁矿床和铅锌矿床及矿化点（成因上与热液活动有关）；而煤矿则相对集中分布在三角形构造的内部。这是因为三角形顶点处应力强度最大，断裂发育，为热液活动和矿化富集提供了有利的空间条件，从而形成了一系列金属矿床，但对煤和天然气富集起了破坏作用，未形成具工业价值的矿床。而在三角形构造内部，由于应力减弱，褶皱平缓，断裂不发育，缺乏导矿储矿构造，对金属矿化不利，但褶皱变形有利于煤矿富集。从区域含矿层分布分析，铁矿和铅锌矿的矿源层（泥盆系和石炭系）在区内广泛分布。根据本区三角形构造变形场和应力场特征及北西、南西两个角点处成矿富集规律，可以预测，在三角形构造东部顶点处有望找到相同类型的工业矿床。其次，区域上黔西古生界具备生油气的地质条件，测区南西侧盘县三角形构造内已有煤成气发现，本区地质和构造条件与之类似，因而有望在该三角形构造中部构造平缓处或穹窿构造中找到具工业价值的气田。

3.3遥感线性体信息场的分“层”和相关信息组合量化解析与成矿预测

线性体场的分“层”量化即在前述线性构造分层解析结果基础上，采用网格法等方法分别对各期次各层次线性构造（成矿分析时结合环形构造等组合）进行量化处理和统计分析，以找出各层次各期次线性、环形构造的统计规律性及其与成矿作用的内在联系，定量、半定量地提取有用信息。不同的研究目的采用不同的量化处理方法。笔者在西昌地区有利含油气构造遥感解析和油气远景预测中[13]，在全区线性构造和环形构造系统解译基础上，首先对区域性断裂系统进行量化处理，建立了该区区域线性体场（断裂系统）的二维和三维彩色量化解析模式，结合其它资料，在构造活动区中确定了有利于油气生成环境的相对稳定的构造区块（小相岭—米市断块中的二级构造单元——米市断陷）；然后再对所选定区进行专题图像处理，提取有关信息，在此基础上，采用环形构造及与其密切相关的环状、放射状大节理有机结合的组合量化处理方法（环形构造采用自环心向环外缘反距离测度加权），结合物探地震、航磁重力及野外调研资料选定有关参数，通过图像图形处理并结合GIS建立了重要构造区块局部构造的二维和三维彩色定量解析自然模型。通过信息场的特征分析和求异，确定了局部构造隆起（异常）和背景的界线，定量提取了地腹隐伏构造的有关信息，作为评价局部构造和油气远景预测的依据。最后，通过综合分析，进行了油气远景预测。

3.4遥感线性体场的背景与异常数学解析

遥感线性体场中包含有许多非地质信息及其它无关的信息。数学分析的目的在于采用数学地质的方法确定线性体场的背景，提取和增强微弱的有用信息，压抑噪声干扰，找出具有地质找矿指示意义的（异常）标志。

异常是针对背景而言的。所谓背景是指研究区空间变量在其分布区域内的空间趋势，它代表了该变量在研究区内的总体分布特征。对线性体场来说，线性构造的区域变化趋势就是背景，反映局部构造存在的环形、线性体是异常，如环状、放射状大节理密集区。确定背景的方法很多，如常用的滑动平均法、趋势分析法、克里格法等。趋势分析将观测值分解为趋势和剩余两部分，前者即为背景，后者代表局部变化的异常特征。也可用矿化标志变量与多个其它变量的自相关分析研究背景场，然后分解不同背景场中的异常。

背景条件按地质总体的多少可以分为单一背景和多背景两种，前者指那些只经历了一次地质作用的地区，各种因素的取值来自于一个总体；后者指经历了多次地质作用的地区，各种因素取值来自于多个总体。多背景的确定方法复杂，但对于多期次构造运动造成的线性体场而言，若采用前述分层解析的方法，先将多期次线性体按“层”分解，那么确定每个层次的线性体场的背景及其求异就可按单背景处理，然后，对各层次（期次）的背景进行综合分析，这样可使问题大大简化，且能有效提高分析精度。

遥感信息场的空域和频域滤波处理是分解和提取线性构造的有效方法。一般而言，低通滤波得出的低频信息能反映遥感信息场中的深部（隐伏地质构造）的信息特征，高通滤波得到的高频信息则能反映地表及浅部有关地质构造的特征信息。不同方位的线性体场可以分别采用不同的定向滤波模板进行处理分解[8]。经滤波处理的图像再经锐化和二值处理，经统计分析或采用特殊的表达方式，如线性体的空间频数图和密度图、中心对称度图、优益度图、变异系数图、空间距离测度图以及各种参数类型的直方图等[9]，能减少人为干扰因素，有效地反映、提取地质构造和成矿异常具指示意义的特征参数。

信息熵是复杂程度和非均一化的数学度量，用它来查明或定量表征遥感信息场中的某种控矿因素的地质变异特征和统一性。以相对熵研究线性体场特征，若地质变量取线性体的长度或频数（或长度/频数），则所得空间密度熵能表征线性体密度场在空间内分布上的差异，若变量为线性体的方位，则所得方位熵可以揭示线性体场在方位分布上的非均匀变化程度。线性体的熵异常可以反映构造的多期次活动性和空间叠加等特征。

其次，分形理论等方法也是研究遥感信息场的有用工具。用分形理论研究线性构造和环形构造的空间特征，有助于查明构造活动性及与成矿作用关系，定量描述有关的自然地质指示标志空间变化的规律性。

3.5三“S”支持下的多源信息综合与复合图像处理和找矿靶区快速评价优选

在大区域范围的成矿预测研究中，对于大量地、物、化、遥资料，采用多源信息综合与复合图像处理技术，能快速提取有用信息并得到形象直观的成果图像表达，已在地质找矿中取得较好效果。然而，这一方法存在的问题是，目前一般应用的遥感图像系粗纠正产品，地物点在空间上存在一定畸变（像点位移）；而且，物化探信息本身存在一定的空间漂移。如果将这种变了形的图像和移了位的信息异场叠合在一起，将产生虚假的综合异常，有时会致人误入歧途。解决的办法是，采用三“S”支持下的多源信息综合与复合图像处理技术，先在GIS和GPS控制点控制下对遥感图像进行精纠正（尤其是投影差改正），制做正射投影遥感影像图，尔后通过区域化变量的结构分析，找出各种地学信息在不同尺度与水平的结构关联性，建立统一的随机场概念模型，并通过GIS实现各类信息在空间上精确匹配，充分利用计算机的可视化技术建立矿源体的自然信息源模型，定量提取综合找矿信息，达到找矿靶区快速评价优选的目的。笔者利用这一方法在南江区调成矿（金）预测中取得显著效果[14]。

南江地区地处特殊的地质构造环境，成矿地质条件较好，矿种较多，矿点星罗棋布。但相对而言，对金矿床的研究，前人研究程度很低，只在西部万金山和坪河发现沙金矿化点两处。本次研究中利用所研制的遥感正射影像，通过信息提取及金、银、铜等14个化探数据的数据处理、地、物、化、遥多源信息综合、复合图像处理有所新发现，并圈定了找矿远景区，经野外调研发现沿图像处理中识别出的构造带河谷中沙金矿化现象普遍，在金化探异常与遥感图像复合信息图上的较高异常地段采样化验分析，金矿化现象明显，已有数个样品达到金的工业品位，将来通过进一步工作，有望取得较大的进展。

参考文献

[1]朱章森，温世明，杨龙.来利山锡矿盲矿统计预测.地质科学，1987,（2）:131～147.

[2]朱章森.无模型预测法刍议.物化探计算技术，1987,14（1）:60～62.

[3]杨武年，朱章森.遥感构造信息场“分层”解析与无模型矿产预测的理论、方法及其意义.中国数学地质（6），北京：地质出版社，1995,63～67.

[4]杨武年，朱章森.遥感信息场分层解析与构造应力场定量研究.地质学报，1997,71（1）:87～96.

[5]Yang Wunian,Zhu Zhangsen.The theory and method of Phase-separation analysis of remote sensing information field of metallogenetic environment and nonmodel ore-deposit prediction.Proceedings of the 30th International Geological Conngress.VSP,International Science Publishers,The Netherlands.1997.

[6]Yang Wunian.Phase-separation analysis of remote sensing images,a new method for determining regional tectonic stress fields.Proceedings of the 30th International Geological Congress.VSP,International Science Publishers,The Netherlands.1997.

[7]杨武年，乐光禹等.恢复区域构造应力场的新方法.科学通报，1991,36（12）:391～934.

[8]杨武年.遥感图像线性特征检测——匹配滤波模板最佳参数选择的定量公式.成都地质学院学报，1985,1（总34）:101～108.

[9]王润生，杨为久.遥感线性体的数量化分析.国土资源遥感，1992,3（总13）：49～54.

[10]杨武年，乐光禹等.金佛山菱形构造格局区域变形场和应力场遥感图像解析.成都理工学院学报，1994,21（1）:102～109.

[11]杨武年，乐光禹.贵州西部六枝地区复合褶皱及其构造应力场遥感图像解析.见：第八届全国遥感技术学术会议论文集，1993，南宁，441～443.

[12]杨武年，黔西六枝—郎岱地区的构造格局及其构造应力场遥感图像解析.国土资源遥感，1996,2（总28）：21～28.

[13]杨武年，李永颐，易显志等.遥感信息量化技术在西昌地区含油气构造地质研究及远景预测中的应用.国土资源遥感，1994,3（总21）：63～70.

[14]杨武年，丁纯勤，王大可等.TM遥感正射影像地图在四川南江地区区调研究及成矿预测中的应用，地球科学（数学地质专刊），1997.

[15]赵鹏大，池顺都.初论地质异常.地球科学，1991,16（3）:241～248.