一、区域重磁资料解释中的相关分析方法（论文文献综述）
楼海[1]（2001）在《重力方法在地壳结构研究中的应用》文中提出在地壳结构研究中重力学方法是一种最古老的地球物理方法。重力方法的理论已经发展的很成熟了。但随着技术的进步，重力方法在地壳结构研究中仍然发挥着重要的作用。本文对于深部构造研究中的重力学方法进行了仔细的分析和研究。对于重力方法的几个重要内容：重力资料的数据处理和反演方法，均衡理论和应用，卫星重力资料在深部构造解释中的应用等几个方面展开了讨论。在本文中对众多的数据处理方法的假设条件和地球物理意义进行了分析。提出在重力数据处理中那些地球物理意义明确，假设条件合理，计算稳定的方法最具有实用性。在对重力反问题基本公式及其多解性的讨论的基础上，提出交互模拟方法和多种资料的综合反演方法是克服重力反问题多解性的重要途径。文中给出了可以用于三维模拟反演的线性变密度四面体重力异常计算公式。为了进行地震-重力综合解释和联合反演，文中首先给出与连续速度模型一致的密度参数化方法，推导出相应的重力异常计算公式，并应用于大别山造山带的地震-重力综合解释。结果表明，综合解释对于认识三维地震层析得到的上地壳高速体的地质意义是很有帮助的。在速度与密度一致的模型参数化的基础上，本文提出了用于密度反演的地震-重力联合反演算法。算法以地震层析方法为基础，用速度-密度关系公式把地震走时残差与密度扰动联系在一起。通过联立求解走时方程和位场方程，构成地震-重力联合反演算法。实验计算表明，在速度-密度关系已知的情况下，增加地震资料的约束，可以改善重力反演的效果。在新疆西北部地区的地壳结构研究中，结合航磁资料和地热资料，对重力资料进行了多种处理和反演。得到了关于研究区地壳结构的很多信息，包括大地构造分区，盆地基底性质，大型断裂带的分布与深部构造意义，岩浆活动性，Moho面和地壳磁性层底面深度等方面的重要信息。对新疆西部地学断面的重力与航磁资料的模拟解释，给出了地壳的密度模型和磁性模型，这些模型反映地壳深部的结构和性质。为地壳深部结构的地质地球物理综合解释提供了依据。在新疆西北部地区，均衡异常显示出天山和准噶尔南缘地区没有达到重力均衡状态，而且，重力均衡作用力的方向与地壳垂直运动的方向相反，说明挤压构造作用与均衡作用相比明显占优势。重力均衡异常显示出天山山脉的均衡状态沿山脉走向方向有很大差异，表明了天山深部物质的不均匀性。同时还显示出西准噶尔与准噶尔盆地在统一的挤压构造作用下，有不同的响应。根据对新疆及其周围地区的卫星重力资料的解释，提出天山及其周围地区处于地幔对流形成的挤压沉降构造环境。天山造山带快速隆起是南北向不对称挤压作用的结果，准噶尔盆地则在挤压作用下弯曲而发生沉降，形成南部坳陷区。卫星重力资料还显示出东，西天山深部构造上的差异。卫星重力资料的解释结果与已有的地震面波反演结果及地质资料是相容的。在腾冲火山区地壳结构的研究中，对人工地震测深资料和重力资料进行了综合解释和联合反演。人工地震测深资料的二维射线追踪解释给出了从中山至自治的南北向穿过火山区的二维速度模型。速度模型显示火山区上地壳内P波速度明显偏低，中、下地壳速度略低。在腾冲至固东之间下地壳内有P4反射面，与<WP=3>其它部位明显不同。速度模型还显示龙陵断裂南北两侧地壳速度结构有较大差异，表明分属不同构造单元。地震资料的三维层析成像给出了火山区上地壳内部的P波低速体的位置和立体图像。低速体埋深约为7km。低速体不在火山区和热海地热区的正下方。参考其他地热，大地电磁，天然地震资料的解释结果，推测低速体可能代表一个岩浆囊或局部熔融体。地震层析成像显示热海地热区位于地表低速区内。地表低速区可能有断裂通道与上地壳低速区相连。可以推测热海热田中的高温热泉和沸泉是深部热流体沿断裂通道上升，与地表热水混合后形成的。重力资料与地震资料的联合反演给出腾冲火山区中下地壳内有一南北向的低密度区。这个低密度区可能是中下地壳受热膨胀，密度降低形成的。火山区地表的张性构造与中下地壳的膨胀有关。综合分析地震，重力及地质资料，对腾冲火山区地壳构造和火山喷发机制提出如下认识：火山喷发是印度板块向东俯冲造成的。俯冲板块的重熔作用在腾冲以下的上地幔中形成岩浆积累，使中下地壳受热膨胀，密度降低，产生地表的重力低异常。中下地壳受热膨胀在地表形成张性构造，同时，也为地幔岩浆上升提供通道，形成火山喷发。
胡斌[2]（2019）在《基于重磁数据的冈底斯成矿带深部结构分析与成矿预测》文中研究表明冈底斯成矿带形成于印度大陆和欧亚大陆碰撞的背景下,具有复杂的构造-岩浆-成矿过程,一直以来都是地学界的热点研究地区。多年来,地质学家对冈底斯成矿带的深部结构、地球动力学背景及区域成矿进行了深入的研究,极大地促进了人们对冈底斯成矿带的地质认识。作为该区仅有的两个覆盖面积全的地球物理资料,重磁数据对解释整体构造框架、圈定侵入岩体水平位置、描绘深部三维物性结构具有重要作用,并可为区域成矿远景区划分提供依据。为了从重磁数据中提取更多可靠的有效信息,本文首先研究了结合功率谱分析的二维经验模态分解位场分离方法,并通过模型试验验证了该方法对叠加异常分离的有效性;接着针对影响位场边界识别效果的各个因素,探讨各边界识别方法的应用效果,获得最佳的边界识别方法组合;最后针对重磁反演垂向分辨率低的问题,本文提出了一种合理的基于先验信息约束的重磁三维反演流程。将上述研究的方法应用到研究区重磁数据的处理、反演和解释中,得到如下主要结果:（1）研究区布格重力异常和化极磁异常均可分为由浅到深的3个等效层场源产生的异常组合;且深部重磁场以南北分区、东西分块为特征,南北以约31°N为界,东西以约88°E和93°95°E为界;深部重力异常解释的低密度体与MT反演结果以及地震接收函数结果综合解释的低速高导体基本吻合;（2）将位场综合边界识别方法应用于重磁异常边界识别,并结合地质资料圈定了基性、超基性（混杂）岩体和中、酸性岩体的分布情况,推断了主要断裂构造;（3）重磁三维反演结果显示研究区局部剩余密度主要由沉积层凹陷、结晶基底的起伏和岩浆岩的侵入等作用引起,强磁性体在南冈底斯带广泛连续分布。通过分析上述结果认为:由于印度板块不均匀地向北推进,使喜马拉雅地体、冈底斯地体的深部物质属性在东西方向上存在显着的差异;青藏高原南部区域中下地壳的低密度体东西向是不连续的;班公湖-怒江洋壳在北冈底斯西部约84°E88°E的范围内先后存在向北和向南俯冲的可能。最后通过综合分析上述成果与矿点分布的关系,总结重磁成矿要素,划分了21个成矿远景区。
许顺芳[3]（2015）在《克拉玛依后山地区重磁场及其与构造格架关系研究》文中提出克拉玛依后山地区是中亚造山带的一个重要组成部分——西准噶尔地区的东南部,该地区被认为在古生代期间经历了复杂的洋陆转换和地壳增生过程,备受学术界关注。论文利用不同尺度的重磁异常资料,以壳幔结构特征为主要目标,结合岩石学与构造地质学等的新成果与新认识,研究该地区现今构造格架,推断其演化过程,具有重要的科学意义。论文利用新近发布的卫星重磁数据、区域重磁数据和克拉玛依后山地区4条高精度重磁剖面测量数据,从不同空间尺度来解析该地区的重磁异常特征,并且利用新方法与新技术对重、磁异常进行了定性分析和定量解释,以及通过反演完成了地壳三维密度与磁化率成像工作,取得了丰富的成果。卫星重磁数据具有较高的观测高度、全球较高的覆盖度、观测数据处理的较统一性以及数据精度与分辨率的较一致性等突出特点,可以弥补传统地面与航空重磁测量的数据空白区以及长波长段的不可靠性。因此,论文选取了西准噶尔及其周边地区最新的卫星重、磁数据进行处理,分析了卫星重磁异常的分布特征,并且对其进行三维反演,构建了地壳与上地幔顶部的大尺度三维密度与磁化率结构,进而解释了其分布特征及其蕴含的大地构造格架和演化信息。针对克拉玛依后山地区的地球物理工作现状,即研究区现有面积性地面重力和航空磁力资料具有不同年代、不同比例尺的特点,论文选用不同小波进行了试验对比,最终采用db4 （Daubechies-4）小波进行多尺度分解,将小波分解5阶逼近作为本区的重力异常的区域背景,6阶逼近异常作为本区岩石圈磁性基底引起的航磁异常；选择了小波模极大值法、归一化导数偏差法、总水平梯度倾斜角法与Canny边缘检测等多种技术处理平面重、磁异常数据,并且对本区地质体边界及断裂构造进行了识别和推断,进而根据地质体边界识别结果,结合局部重力异常和垂向一次导数异常特征,识别了出露地表和隐伏的大型构造；利用分离之后的局部重、磁异常,进行3D反演,结果刻画了克拉玛依后山地区中上地壳的中尺度3D密度与磁化率结构。论文利用实测的4条高精度、大比例尺的剖面重力和磁力剖面数据,分别进行2.5维密度与磁化率成像反演,在分析其分布特征的基础上,结合相同测线的MT反演结果,进行了地球物理剖面对比研究。论文取得的主要认识和成果如下：（1）西准噶尔及邻区卫星重磁异常特征及其区域大地构造格架认识根据重磁异常特征,研究区（40°N～50°N,80°E～90°E）的阿尔泰、哈萨克斯坦-准噶尔、天山和塔里木北缘构造区中的重磁异常特征明显。区域重力场较好地反映了本区主要构造单元的岩石圈深部重力场背景,区域背景之上的局部重力异常则较为详细地刻画了各构造单元之间交汇处的岩石圈构造特征：（i）区内北部的阿尔泰褶皱带东南的山前重力高值带反映了该区域岩石圈高密度特征,明显受到NW向构造的控制,尤其富蕴—福海一带广阔的重力高对应着准噶尔晚古生代活动大陆边缘板段,被认为是中国境内最大的西伯利亚板块二级构造单元；（ii）区内南部的天山褶皱带是本区重力异常变化最剧烈的区域,正负异常条带展示了该区岩石圈深大断裂分布及其相互关系,其中南天山造山带上的重力高值带显示了塔里木板块与其衔接带上岩石圈高密度物质带的存在；（iii）区内天山褶皱带与准噶尔盆地衔接带以北天山北地带上的重力低值带为特征；（iv）哈萨克斯坦地块与准噶尔盆地交汇处呈现的NW及NE向重力高异常带位于准噶尔及其以北的塔尔巴哈台山南缘一带,被认为是洋内俯冲增生带,具有古洋盆、岛弧、后碰撞等多种不同类型构造背景,这或许是引起重力高异常的成因之一。岩石圈磁异常的构造分区特征也很明显：（i）与区域重力场相类似,在研究区东北部呈现大面积低磁性,而西部与东南部呈高磁性分布特征；（ii）区域磁异常在北部的阿尔泰褶皱带NW-SE向分布的低值带,反映了该区域明显受到NW向构造的控制；（ⅲ）区内南部的NW-SE向磁力异常低的区域是天山褶皱带,显示了该区岩石圈深大断裂分布及其相互关系；（vi）塔里木盆地与准噶尔盆地具有高磁性岩石圈；（v）较明显地,克拉玛依后山地区是西北部哈萨克斯坦板块与东南部准噶尔盆地磁性分布的分界或转换地区：（vi）磁异常模量Ta表明准噶尔盆地呈现为外部高、中部略低的圆盘状,向西北方向与哈萨克斯坦板块挤压与碰撞,接触区域即为克拉玛依后山地区,而且似乎由于受挤压,被撕扯的哈萨克斯坦板块东部具有由向东和向西南包围准噶尔盆地的迹象；（vii）在区内准噶尔盆地西、北缘存在一个弧形的低磁化率条带状异常,该异常位于卡拉麦里构造带的西部延伸段上,而卡拉麦里构造带被认为是西伯利亚板块与哈萨克斯坦板块的缝合带。为了能更好地了解该区域岩石圈结构,对卫星重磁资料进行了三维反演,获得了不同深度的密度和磁化率异常分布特征。（i）在上地壳深度上,区域密度变化与表层构造十分对应,尤其在天山褶皱带区域,南、北天山超深断裂两侧均出现密度变异,展示了复杂的构造格局；在阿尔泰褶皱带上,NW-NWW向“片状”密度异常特征十分明显；在西准噶尔区域,塔尔巴哈台山南缘高密度异常十分显着,而在克拉玛依后山一带高密度异常表现较弱,准噶尔盆地内密度变化相对平缓。（ii）在中地壳深度上,在天山褶皱带密度变化异常剧烈,巨大低密度异常区带主要对应北天山断裂带和南天山断裂带的北缘地带；在哈萨克斯坦-准噶尔区域,分布于它们之间地带的低密度带以达尔布特断裂西北的密度异常带最为显着。（iii）在下地壳、地幔顶部,哈萨克斯坦—准噶尔盆地区域以高密度异常为主,但是结构趋于简单,地壳构造痕迹依然可见：天山褶皱带上的超岩石圈断裂带密度异常有可能延伸至更大深度。从反演的不同深度磁化率可以看出：（i）伊利盆地、吐哈盆地、准噶尔盆地与塔里木盆地呈现为高磁性分布；（ii）准噶尔盆地以高磁性“蘑菇”状向西北方向顶入哈萨克斯坦板块的特征非常明显,碰撞接触区域即克拉玛依后山地区呈现北东-南西向的低磁性特征,与西北部哈萨克斯坦板块和东南部准噶尔盆地的高磁性分布相区别；（iii）区内北部的阿尔泰褶皱带NW-SE向分布的高值带反映了该区域明显受到NW向构造的控制,该高磁性异常与密度异常对应的晚古生代活动大陆边缘板段一致,为现今的额尔齐斯河流域发生陆缘地壳拉张——扩张脊的代表性生成物蛇绿岩和早泥盆世中酸性火山岩密切相关；（iv）在准噶尔盆地西北方向的低磁性与高磁性条带（呈北东-南西向）很有可能即为该碰撞与挤压增生的结果；（v）在准噶尔盆地西、北缘存在一个弧形的低磁化率条带状异常,该异常位于卡拉麦里构造带的西部延伸段上,该大型断裂带对幔源岩浆及其分异的残余熔体能够上升到上地壳侵位起着重要控制作用,而构造体制从碰撞期挤压变为后碰撞期拉张,则成为幔源岩浆及其分异产物能够在上地壳侵位的有利条件。从三维密度、磁化率成像结果中截取的三条典型剖面,分析了西准噶尔及周边的岩石圈密度和磁化率结构特征及其构造含义。（i）准噶尔盆地北段以中、下地壳—上地幔密度高、磁化率高为主要特征,高密度、高磁化率体具有向北下插的趋势,推测古洋板片有可能向北俯冲至阿尔泰陆缘岩石圈之下。（ii）在AA’剖面,准噶尔盆北段低密度、低磁化率异常位于卡拉麦里构造带的西部延伸段上：在印度板块与欧亚板块碰撞过程中,准噶尔地块岩石圈相对坚硬,保持了相对完整的结构；在天山褶皱带—塔里木盆地边缘段,天山褶皱带密度结构复杂,显示了在南北两侧的推挤下,地壳与岩石圈发生了褶皱、断裂、重叠和斜插过程,而磁化率结构则相对简单,为低磁化率反映,可能与现今较强烈的构造活动有关。（iii）在BB’和CC’两条横穿西准噶尔构造带的剖面上,均显示准噶尔盆地存在壳幔高密度、高磁化率体,且由SE向NW下插,在达尔布特断裂带附近,存在壳幔高密度体聚集,在达尔布特断裂带以北段存在低密度、低磁化率异常段,与AA’剖面中卡拉麦里构造带的磁化率结构相似。（2）克拉玛依后山地区区域重力与航磁异常特征及其地质构造意义选择了小波模极大值法、归一化导数偏差法、总水平梯度倾斜角法与Canny边缘检测等多种技术处理平面重、磁异常数据,并且就本区地质体边界及断裂构造进行了识别和推断,进而根据地质体边界识别结果,结合局部重力异常和垂向一次导数异常特征,识别了地表出露的达尔布特断裂、哈图断裂,及隐伏在盆地内覆盖区的克乌断裂和本区南部一条近南北向的区域性断裂F4,其向北可延伸至小拐以西,长达80 km以上,断裂两侧局部重力异常特征具有明显差异。应用磁性地质体边界识别结果和磁异常解析信号振幅,推测了具有较强磁性的地质体分布,地表发现的较强磁性蛇绿岩混杂岩（含基性岩）带均被预测,且推测这类磁性岩石很有可能在盆地内存在。原始重力异常与小波分解5阶逼近异常之差作为反演的局部（或剩余）重力异常,局部异常在克乌断裂以北,达尔布特断裂以南为全区异常幅值最高、范围最大的重力条带；铁厂沟及其以东具有两个重力低异常带,位于乌尔喀什尔山南侧山前的新生代覆盖区,可能是山前断陷盆地的反映。以下延500 m的航磁异常作为浅源磁异常,主要与区内岩浆岩、蛇绿岩混杂岩带的分布有关。利用地面重力与航空磁力异常,进行三维反演,刻画了克拉玛依后山地区中上地壳的中尺度密度与磁化率结构。克拉玛依后山地区中尺度三维密度和磁化率特征分析了该区的构造、岩体分布。密度和磁化率异常条带的边界与地表的断裂构造带,NE-SW向的哈图和达尔布特断裂,极为吻合,同时对隐伏的NE-SW向的克-乌断裂及SN向的F4断裂,均有明显的反映。哈图岩体、阿克巴斯陶岩体、洪山岩体及克拉玛依岩体周边均能见到密度和磁化率异常体,且发育深度较深,这类异常体均与线性构造密切相关,可以认为受到断裂的控制。本区最高密度、高磁化率的异常条带位于克-乌断裂带西北侧,推测很可能为蛇绿岩带产生,且与克拉玛依蛇绿岩带相关。位于铁厂沟南侧的低密度异常条带,浅层低磁化率、深层高磁化率的近EW向展布的的异常条带,控制了克拉玛依后山地区北部的构造格局。该区域受NW-SE向的碰撞挤压和拉伸,形成了一系列NE-SW向的达尔布特断裂、哈图断裂、克-乌断裂,很可能是俯冲增生楔环境下的区域性断裂,与附近的大规模花岗岩形成具有密切的关系,断裂体系控制了该区准噶尔泥盆-早石炭世岛弧（准噶尔中新生代上叠盆地）、克拉玛依后山奥陶纪—晚泥盆世蛇绿混杂岩带、克拉玛依后山早石炭世残留海盆、达尔布特中志留世—中泥盆世蛇绿混杂岩带、铁厂沟-白杨镇泥盆纪岛弧的地层和构造格局。（3）克拉玛依后山地区高精度重磁剖面的构造认识Ⅰ和Ⅱ线的密度成像剖面,均反映了地表出露的哈图断裂和达尔布特断裂带及隐伏的克乌断裂,同时对Ⅰ线上高密度的克拉玛依和达尔布特蛇绿岩带、低密度的洪山岩体具有较好的反映；Ⅱ线上高密度的克拉玛依和达尔布特蛇绿岩带,低密度的克拉玛依岩体具有较好的反映。Ⅲ线的密度成像剖面,也反映了地表出露的哈图断裂和达尔布特断裂带。Ⅰ、Ⅱ、Ⅲ和Ⅳ线的磁化率成像剖面对构造断裂反映不明显,但对测线上分布的高磁化率克拉玛依和达尔布特蛇绿岩带均有明显的显示。这与该区域构造两侧岩性结构相近有关,这该区的物性分析中的高、中等磁化率和高密度岩石以及高磁化率岩石均与蛇绿岩带相关的结果相吻合。
付光明,严加永,张昆,胡浩,罗凡[4]（2017）在《岩性识别技术现状与进展》文中研究说明自从1815年世界上第一幅地质图诞生以来,岩性识别即如何刻画、认识地下岩性分布一直是地质学中的重要问题.一个世纪以来,地质填图技术从手工填图发展到了数字填图,但鉴于地表地质填图受深度的限制,只能反映浅表岩性的分布,难于识别深部岩性,在覆盖区更不能给出下伏岩性的变化,限制了资源勘查的效果.为了克服地表地质填图的缺点,地球物理、遥感等技术逐渐引入到岩性识别中,并取得了长足进展.本文从岩性识别原理、方法、效果等方面入手,结合实例,对比分析了当前主要几种岩性识别技术的特点.综合分析认为,三维岩性识别将是资源勘查中重要的一个环节,如何准确识别地下岩性的三维分布特征是实现深部矿产、油气资源勘查突破的关键;从技术层面分析,认为采用钻孔等先验信息约束,开展重磁三维反演,获取地下密度和磁化率模型,根据物性与岩性的逻辑关系,识别岩性的三维分布特征是现阶段三维岩性识别最有可能成功并适合推广的方法.
张建民[5]（2019）在《金川铜镍矿区地球物理特征及深部成矿预测》文中认为加强矿产资源储备是国家经济持续稳定增长的重要保障,而勘查深部资源是实现资源保障的重要方向。金川铜镍矿床是我国最大的铜镍生产基地,在此开展深部勘查对增加资源储备具有重要意义。关于深部勘查,地质研究是基础条件,地球物理能够依据地下介质的物性差异进行找矿,方法手段多元、探测深度大、分辨率高,是重要的科学技术支撑条件。本文以地质研究为铺垫,通过对分辨率互补的地球物理位场勘探和电法勘探数据的处理和解释,对金川铜镍矿区深部的成矿有利部位进行了推断。1.通过系统收集和分析前人资料及卫星重力数据处理,野外踏勘,岩矿鉴定及物性测试,获得以下认识:（1）金川铜镍矿床位于区域北西向展布的重、磁异常梯级带上,处于华北板块与祁连地块的交接部位,与已知的龙首山陆缘带相对应,显示了其形成与深大断裂的发育有关。界面反演结果表明,龙首山处于地壳厚度的陡变带,地壳厚度约55 km,北东向的华北板块地壳厚度变化相对较小,反映了其稳定地块的属性;西南部的祁连地块地壳厚度较大,且不均一,说明该区构造运动复杂,曾受到强烈的推覆挤压作用,而导致不均一的变形。（2）铜镍矿床的产出与地层、构造和岩浆岩的关系密切,含矿岩体主要侵位于长城系白家嘴子组地层中,走向北西,倾向南西,倾角较大;与北西向构造关系最为密切,北西向断裂常起控矿作用,褶皱往往起改造富集作用;矿区超镁铁质岩体是成矿母岩,亦是围岩,矿体往往产出于岩体的下部。矿床成因不一,主要以深部熔离-贯入成矿为主。矿石矿物主要由雌黄铁矿、镍黄铁矿、黄铜矿、橄榄石、辉石等组成,完整矿物晶形较少,矿物裂隙发育,且常具有明显的蚀变现象,表明矿体形成具有多期次的特征,且曾受到强烈的后期改造作用。（3）基于密度、磁性和电阻率特征,矿区岩、矿石可大致分为三类,第一类是铜镍矿石,表现出高密度、高磁性、低电阻率的特征;第二类是超基性岩体,呈密度高、磁性强、电阻率中等的特征;第三类是超基性岩体的围岩,具有密度低、磁性弱、电阻率高的特征。这些物性特征的差异性,为研究区地球物理的勘查工作提供了前提条件。2.为获取更丰富的与矿相关的信息,对位场数据处理的方法进行了研究。（1）传统的位场数据分离方法仅强调区域场和局部场或多尺度分离的方法存在模态混叠的现象,为了更多且更有效地提取不同尺度的位场数据信息,一维和二维双树复小波被首次引入到位场数据的处理中,该方法具有近似的平移不变性、更多的方向选择性,以及限制数据冗余的特点,模型试验表明,该方法相比小波分析具有优越性,可以更好地实现多频率混叠信号的不同尺度信息分离。另外,基于数学形态学的基本算子,通过调整不同尺寸的形态学结构元素,提出了对一维和二维信号实现多尺度分解的形态学算子,模型试验表明,形态学算子能够很好地对不同尺度的信号进行分离,而且具有较高的计算效率。（2）针对传统位场数据边界识别方法常存在输出边界模糊、抗噪能力差、不能均衡强弱异常或在正负异常同时存在时具有虚假边界存在的问题,基于数学形态学良好的边缘刻画能力及抗噪性能,本文先构建了一种新的边界识别算子MMA,该算子可以利用其极大值有效地对位场数据的边缘进行识别,为了使其具有均衡强弱异常的能力,提出了MMAZ和MMAT两种均衡的边界识别滤波器,模型试验结果表明,相比倾斜角法、倾斜角的总水平导数、Theta图法等几种传统的边界识别滤波器,这两种滤波器不仅可以均衡不同幅值的异常,而且输出的边界更加收敛,同时在正负异常同时存在的情况下,没有虚假额外边界的产生,以及具有良好的抗噪能力。（3）反演是位场数据处理与解释中定位异常地质体空间位置、物性参数估计及地下构造推断的重要方法,对深部矿产勘查具有重要意义。针对位场数据反演常存在多解性、反演结果精度低以及计算效率低的缺点,本文采用光滑L0范数紧支撑聚焦方法及GPU并行计算,来快速地获取准确度高的位场数据反演结果。并通过不同埋深的两个地质体模型对该方法的效果进行了验证。3.利用这些位场处理方法对矿区的航磁异常数据进行了处理。边界识别结果显示,矿区北西向构造最为发育,北东向次之;已知含矿岩体的及矿区南延区正航磁异常的边界被确定。位场分离及反演结果均表明,矿区深部具有良好的找矿前景。4.为了更好地了解不同地球物理方法的深部探测能力及建立找矿模型,通过多个模型的重、磁模拟以及矿区以往实际的CSAMT数据反演结果,对深部异常体的位场异常特征及CSAMT的找矿指示进行了研究和总结。对于深部找矿,应该重视研究区航磁数据等值线宽缓的低正异常区域,但低正舒缓磁异常未必反映的是深部的超基性岩体,亦有可能是与超基性岩体具有几倍磁性差异的混合岩带或部分变质岩系等其他磁性略大的地质体,这是因为对于同等规模的异常体,由有效磁化强度不同而引起的磁异常差异在浅部表现的比较明显,差异较大,在深部则相反,差异较小。另外,鉴于研究区岩、矿体较大的倾角（50°80°）,该区磁异常变化基本可以忽略岩、矿体倾角的影响。在研究区异常幅值相同的不同区域,正异常水平距离的范围越大,则可能代表目标体埋藏越深,规模越大。相似特征的磁异常可能是分布复杂的磁异常体所引起的,对于深部找矿,要尽可能多地了解和结合丰富的已知地质、钻探等综合信息,从而能更为合理地给出有利的深部找矿位置。面对深部找矿,重力与磁法基本表现出相似的异常特征。电阻率差异反映出的破碎带亦是重要的找矿指示。结合地质、地球化学、遥感及其他找矿标志,建立了以地质-地球物理为主的找矿勘查模型。5.依据找矿勘查模型,结合矿区内实测的重、磁、电剖面数据的处理结果,对矿区及周边深部有利的找矿地段进行了分析,并圈定有利靶区6个。关于主矿区东部大范围的高磁异常区（东湾勘查区）,在其深部未找到有利的找矿地段。
韩兆红[6]（2011）在《利用重磁资料进行构造边界识别与弱异常提取的方法研究及应用》文中研究表明本文的研究内容主要包括利用重磁数据进行边界识别和寻找区域弱异常等项研究,同时通过理论模型试验加以证实,最后对具有复杂区域地质构造格局的齐齐哈尔地区重磁资料进行了相关处理,给出了初步地质解释,主要工作如下：1.提出了利用双方向区块扫描Lyapunov指数极小值确定物质区块分界线的方法,进行了有益的尝试,具体过程如下：（1）重力和航磁的网格化数据都满足行或列在同一直线上,以每条线为一个窗口,即以行或列作为单变量时间序列直接计算其Lyapunov指数。（2）把m×n数据矩阵按照行来划分,得到n行数列,可以计算出每一行对应的一个Lyapunov指数值；相似地按照列来划分,可以得到每一列所对应的一个Lyapunov指数值。这样实质上就是得到两个Lyapunov指数序列。（3）在所得到的两个Lyapunov指数序列中,行和列两个方向交汇处的每一个点都存在唯一的一对Lyapunov指数值,当满足双方向的Lyapunov指数值同时达到极小时,这个点的坐标就是所计算曲面拐点或极值的平面位置。2.通常重磁异常曲面的拐点或极值的连线往往对应地下地质体物性的分界面或轴线在地表的投影。为了检验上述Lyapunov指数方法的有效性,进行了台阶、水平圆柱体两个理论模型试验,并对齐齐哈尔地区实测重力数据进行了处理。理论模型试验结果与设计位置吻合,实际处理所得极小Lyapunov指数值点的连线与大兴安岭梯度带边界吻合。3.为了更有效地提取弱异常,文中总结了自相关滤波法在重磁数据处理中的优缺点。通过建立单球体和多球体模型试验发现,秦葆瑚1991年提出的自相关滤波因子具有基本消除了随机干扰的影响、能够在一定程度上反映区域异常的变化趋势、能显示异常有无的优点；但也存在滤波结果仍受区域场影响、提取的弱异常幅值太小且伴随有假异常、其滤波公式部分物理意义不明确等缺点。针对自相关滤波方法在提取局部弱异常中的缺陷,对一维和二维自相关滤波因子进行了适当改造,具体方法如下：（1）一维自相关滤波因子改造：在每条测线内自动寻找剩余局部异常最大点处的观测值做为参考模型,并沿整条测线进行滑动自相关处理。数据处理时,只需通过分析异常规模给出滤波窗口K,并让K等于最大异常范围半径,模型试验结果位置准确,减小了区域场的影响,异常相对强度明显变大,有效地突出了弱异常。当对含有随机干扰的数据进行处理时,不能完全消除随机干扰,虚假异常分布不对称,区域场大的一侧虚假异常强度较大,处理结果有窗口滑动痕迹。（2）二维自相关滤波因子改造：改变自相关滤波函数ρ的系数,以及四周测点的平均值Z（l）的表达式。二维自相关滤波只适用于线距等于点距的情况,对于其它类型数据必须先网格化成正方形网格；自相关半径选取在3～6个网格节点距为佳,异常零值线能较好反映异常体范围；对随机干扰压制效果有所改善,不仅消除了背景影响,有效地提取了弱异常,还可以通过零值线大体划出异常边界；高值突变点对本方法有较大影响,所以要求对原始资料先作剔除畸变点的预处理。4.为了检验改造后自相关滤波因子的效果,对齐齐哈尔地区实测航磁数据进行了处理,结果能很好地确定出各局部异常位置,边界划分清晰,基本提出了强区域场内的所有弱异常。但由于自相关滤波因子对随机干扰压制能力有限,不能完全消除它们,结果中可能还存在一定虚假异常。选择用匹配滤波和自相关滤波进行对比分析,发现处理结果有区域场影响小、局部异常独立性强的特点,更好的反映了该区局部异常的分布状况。按局部弱异常密集程度将该区大致分为两个区：第一区属大兴安岭华力西褶皱带,局部弱异常受区域地质构造控制,构造线多数沿NE～NNE向延伸；第二区为松嫩平原地带,地磁场相对平稳,无明显磁异常。其中A处具有较大研究价值,这里磁异常反应强烈,相对集中,对此地可做为找矿远景区进一步勘探。最后绘制了构造特征线图。5.介绍了分形几何的理论知识、分形与混沌成矿的关系,总结了分形理论在重磁数据处理中的应用。由于分形维数和Lyapunov指数都是描述非线性系统特征的量,所以把双方向Lyapunov指数极小值识别物质区块分界线的思想应用到分形理论中对地质体进行边界识别,这是将来本人要做的后续工作之一
马国庆[7]（2013）在《位场（重&磁）及其梯度异常自动解释方法研究》文中研究指明随着地球物理仪器勘探效率、探测精度、数据参数和容量的不断提高，很多传统的位场（重、磁）数据自动解释方法由于存在精度低、计算分析复杂、人为干预多等缺陷而难以满足当今地球物理勘探的需求，快速的、高精度的位场自动解释方法成为地球物理工作者研究的热点问题。本文提出多种位场及其梯度异常自动解释的新方法，主要涉及小子域滤波法、边界识别滤波器、欧拉反褶积法、解析信号法、局部波数法、物性反演及相关成像法等，获得了更加令人满意的结果。实测重磁异常是地表到深部所有密度和磁性不均匀分布的综合反映，为了提取出目标体的异常需对综合异常进行低通滤波处理。频谱分析技术是进行重磁异常场分离操作的常用方法之一，但该类方法易模糊掉不同形态异常之间的界线特征。小子域滤波法是一种可有效保留异常之间界线特征的场分离方法，但常规小子域滤波法未充分考虑计算点本身已处于区域场以及弱异常去除等问题。本文提出优化小子域滤波法，该方法通过对常规方法的子域剖分方式和选取准则进行调整可更有效地完成场分离工作。理论模型试验表明优化小子域滤波法的输出结果更加合理、准确，有效地改善了常规小子域滤波算法存在的不足。将其应用于实际磁测数据的处理，能很好地完成数据的场分离工作，且相对常规小子域滤波法能更有效地保留异常的界线特征。边界识别是位场数据自动解释中必不可少的任务之一，其可清晰地反映出地层之间的界线以及场源体的分布范围。现有的边界识别滤波器大多仅能识别出较浅地质体的边界，而较深地质体的边界则比较模糊。为了改善这一问题，本文提出水平与垂直导数的相关系数法进行地质体边界的识别。通过理论模型试验表明该方法能同时识别出浅部与深部地质体的边界，且相对其它方法所识别出的边界更加清晰。将相关系数法应用于实际数据的解释，获得了断裂的水平位置及地层之间的界线，与实际地质资料对应较好。垂直导数受噪声干扰较大，会影响识别结果的稳定性，为此提出归一化总水平导数法进行地质体边界的识别。理论模型试验表明归一化总水平导数法能同时识别出浅部与深部地质体的边界，且相对其它方法更加稳定，受噪音影响较小。将其应用于实测重力异常的解释，获得了断裂的分布特征。但相关系数法和归一化总水平导数法所识别出的边界均存在一定程度的发散，不易准确地划分出地质体的分布范围。为了获得更加清晰地边界识别结果，提出增强型均衡滤波器，其为不同阶导数之间的非线性组合，并引入一种基于Laplace方程的垂直导数稳定算法来降低噪声的干扰。模型试验表明该方法能更清晰和准确地识别出地质体的边界，且对弱异常有更好地识别效果。将其用于实际数据的解释，所划分断裂的位置及走向均与实际地质资料相吻合，并发现了更多的细微构造。增强型均衡滤波器需要计算异常的三阶导数，当点距较大时高阶导数的计算是不稳定的，为此提出增强型局部相位滤波器，其仅利用一阶水平导数与二阶水平导数的非线性组合进行地质体边界的识别。理论模型试验证明该方法依旧能很清晰地识别出不同深度地质体的边界，且相对增强型均衡滤波器其结果更加稳定。将其应用于实际数据的解释，清晰地给出了地层的标记及异常体的分布范围。经研究发现现有边界识别滤波器在进行由正负密度体所产生的叠加异常的边界识别时会产生多余的干扰边界，本文提出优化边界识别滤波器，其通过在边界识别滤波器公式中加入一较小常数来消除多余边界的干扰，通过理论模型和实际数据试验证明了方法的有效性。将本文所提出的高分辨率边界识别滤波器插入到软件平台-Geosoft，补充和增强数据处理与解释功能，实现了方法研究向应用模块的迅速转化，与国外高端软件平台实现了良好地对接。随着位场勘探数据量的不断增大，地球物理人员更加倾向于利用自动解释方法进行数据的解释。但现有位场自动解释方法普遍存在计算条理不清晰、精度低、计算复杂等缺点。本文对自动解释方法中的欧拉反褶积法、解析信号法和局部波数法进行改进，使反演结果更加准确、稳定。扩展（梯度）欧拉反褶积法具有不受背景异常干扰的优点，相对常规欧拉反褶积法具有更广的应用范围。该方法有两种不同的求解策略，一是位置参数和构造指数同时求解出来；二是先计算异常体的位置参数再计算构造指数。现今均采用第一种策略进行异常的反演，本文推导出了利用第二种策略进行求解的公式，并通过理论模型对比了不同求解策略反演结果的差异，反演结果表明第二种求解策略的精度较高，且反演结果更加稳定。将第二种求解策略应用于铁矿区航磁异常的解释，获得了铁矿的实际分布特征，该项研究对于扩展欧拉反褶积法的实际应用具有很强的指导意义。解析信号法具有不受磁化方向干扰的特性，被广泛地应用于磁异常的解释。现有的解析信号法在进行深度和构造指数计算时需计算异常的三阶导数，这会明显地增大噪声的干扰，从而降低反演结果的可信度。本文提出三种优化的解析信号法，前两种方法需计算异常的二阶导数，而第三种方法仅需要计算异常的一阶导数，能有效地提高结果的准确性和稳定性。通过理论模型和实际数据试验证明改进的解析信号法能成功地完成磁异常的解释，且精度高于其它同类解析信号法，尤其是第三种解析信号法，其计算结果非常稳定和准确。局部波数法是近几年应用较为广泛的一种自动解释方法，其大多利用局部波数的导数进行深度和构造指数的计算。本文提出利用不同局部波数曲线之间的组合进行异常体参数的计算，可有效地降低噪声的干扰，提高反演结果的准确性。理论模型试验和实际数据应用证明改进局部波数法的反演结果精度较高。此外，还推导出一种可将位置参数和构造指数同时反演出来的局部波数法，对于大批量数据解释更有优势。通过理论模型试验和实际应用可以看出，该方法的精度不低于其它局部波数法，具有良好的实际应用效果。上述三种自动解释方法的组合可形成针对复杂应用对象的多种解释技术选项，有效地降低了解释结果的不确定性和风险。本文提出利用快速模拟退火法进行位场异常的物性反演，有效地避免了大规模方程的求解，提高了效率，并以异常的均方差和迭代前后属性参数的变化同时作为迭代停止条件，提高了计算结果的准确性。理论模型试验表明该方法可准确地完成密度和磁化率的反演，与实际地质模型相一致。将该方法应用于辽源市煤矿采空区重力异常的解释，获得了采空区的分布范围，其与电法测量结果一致。位场相关成像法是根据实测异常与假定场源所产生异常的相关系数来快速地确定地质体的空间分布。现有的相关成像法是利用球体模型来模拟地下地质体的形状，当场源体的实际形状与球体相差较大时，计算结果会产生较大误差。为了改善这一缺点，本文提出增强型位场相关成像法，其分别计算不同类型场源体引起异常与实测异常的相关系数，使相关系数取得最大值的场源体类型与实际地质体类型相一致，相关系数的最大值对应于地质体的中心，可同时完成场源体深度和类型的估计。磁异常的相关成像采用异常的解析信号来进行，有效地避免了磁化方向的干扰，且计算公式更加简单。理论模型试验表明该方法可成功地完成位场数据的反演，且稳定性较高，最后将其应用于上海实测磁异常的解释，获得了地下未爆炸物的分布范围。张量探测技术是一种新兴的地球物理观测手段，其能提供地质体在不同方向上的导数反映，能更准确地描述地质体的特征。但现今全张量数据的解释方法缺乏系统性的研究。本文提出张量局部波数法和方向解析信号法进行全张量数据的解释。张量局部波数法是利用张量局部波数组成的线性矩阵求解场源体的位置参数。理论模型试验证明张量局部波数法能成功地完成张量数据的解释，且精度高于常规局部波数法，将其应用于美国实测重力数据的解释，获得了地质体的深度，与已有的解释结果相一致。方向解析信号法是利用x,y,z方向的解析信号进行磁张量异常的解释，并采用水平方向解释信号进行异常体边界的识别。理论模型试验表明方向解析信号法能有效地完成磁张量数据的解释工作，且反演结果和识别出的边界受磁化方向影响较小。将其应用于某矿区实测磁梯度数据的解释，获得了铁磁性物质的分布形式。为了完善Geosoft平台高精度重磁数据处理与解释流程，将所提出的张量数据解释方法插入到该软件平台。
李丽丽[8]（2013）在《地球物理（重磁震）勘探数据一体化解释技术》文中提出重力、磁力及地震（简称重磁震）等勘探数据的一体化反演是目前地球物理数据综合解释中的一个重要研究方向，是解决复杂地质问题最有效的途径之一。重磁勘探可快速地完成区域地球物理调查，了解地质构造（地层、断裂）分布特征；地震勘探是目前进行地层精细结构划分最常用的勘探手段之一，具有较高的垂向分辨率，但对于某些特殊构造（盐丘、高角度断裂等）反映不敏感而形成盲区。重磁数据具有提高陡倾角地层分辨率、分辨地震盲区（泥岩断层、盐丘下部）、改进地震速度模拟等作用，因此重磁震数据的计算机综合模拟可建立更高精度的地质-地球物理解释模型，其中模型参数和反演方法的精度是制约最终反演结果准确性的关键因素，因此本文针对模型参数计算方法和联合反演方法进行改进，以获得更加准确的反演结果。本项研究中重磁震数据一体化反演按如下流程进行：（1）利用高精度同相轴识别技术解释地震数据，获得层位的分布信息；利用改进的边界识别方法解释重磁数据，分析研究区断层和地层界线的分布特征；（2）利用改进的位场数据解释方法根据重磁数据来估算盐丘、陡倾角地层及浅部异常体的属性（深度、密度、磁化率）参数，用于建立地质-地球物理模型；（3）利用地震随机反演计算地层的波阻抗，进而实现地层密度分布函数的计算，并通过对随机反演结果进行不确定性分析来获得更加可靠的结果；（4）以地震和井数据为约束利用重磁和重力梯度异常计算整个测区的层位分布信息。根据上述信息建立地质-地球物理模型，并利用快速模拟退火法实现重磁、重力梯度及地震数据的联合反演，获得最终的地质-地球物理模型。同相轴识别是地震数据解释中最基本的任务，可清晰地划分出层位的分布情况。为了提高地震同相轴的分辨率，本文提出增强型数学形态法，该方法利用数学形态学中膨胀和腐蚀运算与水平导数相结合进行地震同相轴增强。通过理论模型试验证明该方法能成功地完成地震同相轴的增强工作，具有较高的分辨率，且相对其它方法能更清晰地识别出地震同相轴，受噪音影响较小，尤其对弱振幅同相轴也有较强的识别能力。将该方法应用于实际地震数据的解释，获得了层位的具体分布信息。同时，经研究发现，增强型数学形态学法经过相应的修改后可用于重磁数据的边界识别，即直接利用原始重磁数据总水平导数腐蚀与膨胀结果的比值进行地质体边界的识别，结果的最大值与地质体的边界相对应，该方法不需要垂直导数参与计算，受噪音和外界环境干扰小。通过理论模型试验证明，该方法能同时识别出深部与浅部地质体的边界，且抗干扰能力较强。将其应用于实际数据的解释，清晰地给出断裂的位置及分布范围。但数学形态学法所识别的位场数据的边界存在一定程度的发散，不利于准确地划分地质体的边界。为了获得更加令人满意的结果，对数学形态学计算结果进行幂次和指数变换来增强边界的清晰度。通过理论模型试验和实际数据证明幂次和指数变换能有效地增强边界识别结果的对比度，使地质体的边界得到更加清晰地显示。地质-地球物理综合模型建立的一项关键技术是研究区内地质体属性参数的计算。地震勘探对于盐丘构造、陡倾角地层及范围较小地质体反映不灵敏，难以准确地描述地质体的埋深、密度、磁化率等信息，为此，利用重磁异常辅助地震数据来完成参数的计算。现有重磁数据解释方法精度低、计算速度慢，严重地影响了多元数据融合的精度和效率，本文提出几种高精度解释方法来获得更加准确的地质体属性参数。利用重力异常水平导数与垂直导数非线性组合定量计算断层深度、密度、倾角等参数。通过理论模型试验证明，该方法反演得到的断层参数与模型理论值相接近，误差较小。将该方法应用于已知断层重力异常的解释，获得的断层几何及密度参数与其它资料解释结果吻合较好，表明该方法可有效地辅助地震数据完成断层的解释工作。采用解析信号法和局部波数法进行磁异常的解释，因为这两种方法具有不受磁化方向干扰的特性，其计算公式更加简单，但现有方法需要计算异常的三阶甚至更高阶导数，会明显地增大噪声的干扰，降低反演结果的准确性。本文对解析信号法和局部波数法进行改进，改进后的方法仅需要计算异常的二阶导数，有效地提高了反演结果的准确性。改进的解析信号法是利用解析信号水平与垂直导数之间的关系来估算磁源体位置和构造指数信息。通过理论模型试验证明改进的解析信号法能成功地完成磁异常的反演，具有较强的适应性。将其应用于实测磁异常的解释，反演结果与测井解释结果相一致。改进的局部波数法是利用局部波数在不同水平位置或不同高度值的组合进行磁源体参数的反演。通过理论模型试验证明改进的局部波数方法受噪音影响小，能更加稳定和准确地完成磁异常的反演,反演结果与理论值之间的差距小于5%。将其应用于实测磁异常的反演，其反演结果与解析信号欧拉反褶积法的反演结果相一致，表明本文方法具有良好的实际应用效果。为了真实地描述地层的密度变化，采用地震随机反演来计算地层的密度信息。随机反演是地质统计学理论与基于模型反演相结合的产物，亦称地质统计学反演。随机反演充分地利用了地震反演储层预测及随机建模储层预测的优势，有效地综合了不同尺度数据信息（地质、地球物理），突破了地震频带宽度的限制，能建立更高精度的地层波阻抗模型，且可获得储层非阻抗信息（孔隙度、伽玛等），具有更直接的岩石地球物理意义。随机反演会产生多个等概率结果，采用不确定性分析技术对反演结果进行评估可获得更为可靠的结果，在反演中起着至关重要的作用。通过理论模型证明随机反演可有效地完成地震数据的波阻抗反演，将其应用于实际数据的解释，获得了地层的波阻抗分布，并利用不确定性分析技术对反演结果进行评价，通过计算结果的均值、方差、概率密度等参数来获得更为可靠的结果。此外，通过与确定性反演的对比可以看出，随机反演具有更高的分辨率，能获得更多的细节特征。获得地层的波阻抗信息后，利用波阻抗与波速的关系可获得地层的密度分布。对于不存在三维地震数据的地区，以剖面地震和井数据为约束利用重磁和重力梯度数据来计算整个测区的层位分布信息，如存在三维地震数据，则直接利用同相轴识别技术来获得层位信息。重力梯度对浅部地质体、地层的细节特征及起伏变化反映敏感，但现今还没有重力梯度数据进行界面反演的详细介绍。首先推导出密度界面重力梯度异常的正演公式，并给出快速模拟退火法利用梯度数据进行密度界面反演的基本步骤。理论模型试验证明利用快速模拟退火法根据梯度数据可成功地完成密度界面的反演，且对界面的局部变化具有较高的分辨率，最后将其应用于中国南海海域重力梯度异常的解释，获得了海底地形分布，其结果与LDEO船测结果差距很小。利用获得的层位和参数信息建立地质-地球物理模型进行重磁震一体化反演，并将梯度异常引入联合反演中来提高对地层局部变化的分辨率。快速模拟退火法计算速度快，对于大数据解释具有优势，利用其进行重磁、重力梯度和地震数据的联合反演。通过试验证明快速模拟退火法能有效地完成重磁、重力梯度与地震数据的联合反演，且对界面细节特征有较高的分辨率。将该方法应用于地区实测数据的解释，揭示了地下层位的分布特征，为进一步勘探提供了决策依据。
黄金明[9]（2013）在《重磁数据处理解释技术在华南地区岩体圈定与形态反演中的应用研究》文中研究说明随着国民经济的高速发展，对矿产资源的消耗量迅速增长，过去探明矿床的资源储量满足不了日益增长的需求，迫使地质工作者转变思路，另辟蹊径，开拓新的找矿空间。向地层深部进军，寻找更多、更好的隐伏盲矿体，成为地质工作不可忽视的新领域。区域重磁资料解释是揭示深部岩体最为快速、有效的途径。区域重磁资料在划分大地构造及区域地质构造，研究岩浆岩带分布与圈定隐伏岩体等方面独特的作用。深部和隐伏岩体对于寻找深部隐伏矿具有极为重要的价值，通常在岩体顶盖、边部以及岩体外带、内带等都是金属矿床的赋存空间。所以隐伏岩体的发现与空间范围的定量化研究，不仅可以作为直接或间接的找矿手段，可以为深部找矿、矿产资源调查评价的选区和工作部署提供重要的地球物理依据，对于我国开展下一个深度的金属矿产潜力调查具有重要意义。本文开展并完成了以下研究工作：第一、分别从理论模型处理和实际资料处理结果两个方面对常用岩体的定性解释方法技术进行了评价，提出了一套深部岩体的圈定、反演和综合解释的技术流程。第二、基于国内现有的区域重力、航磁资料进行二次开发，基于传统的重磁位场转换方法，结合斜导数法等场源边界识别新技术，参考地质图、地球化学和其它资料，完成了华南地区深部隐伏岩体圈定。第三、在重磁数据处理技术方面，提出了空间域和频率域数据扩边、网格化数据空区插值技术，实现了2.5维剖面反演模型数据和三维自动反演模型数据转换为三维人机交互程序模型文件，并将这些技术应用到华南地区深部隐伏岩体圈定与形态研究工作中，取得了令人满意的应用效果。第四、基于研究区内岩石的密度和磁性资料，充分参考利用地质、钻探、化探资料和其它地球物理资料，对华南地区的部分岩体异常进行了三维定量反演，推断出了这些岩体的埋深、大小和空间形态。
周帅[10]（2017）在《重磁数据地质目标位置及物性高精度解释方法研究》文中指出随着卫星、航空和地面重磁及其张量勘探的综合应用,以及数据精度和勘探精度要求的不断提高,传统单一尺度的解释方法精度低、人为干预多等缺点日益突显,迫切需要研发地质目标体位置及物性的快速、高精度解释方法。本文从重磁数据正演模拟分析出发,采用数据分析工具定量分析勘探测量组合参数对目标体探测能力的影响,以及不同测量场类型和多个观测面高度获取的多维数据对目标体解释结果的影响规律;在此基础上,提出重磁及梯度异常的目标体位置、埋深及物性反演的新方法,获得准确的地质体三维空间分布及物性特征,从而为地球物理建模提供准确的基础数据和可靠的依据。首先从应用角度,针对地球物理数据采集阶段,引入数据分析工具深度分辨率图和离散Picard图分析测量组合参数对解释结果的影响,指导反演方法的建立及勘探测量组合参数的设计。另外,本文利用数据分析技术对重磁异常与梯度异常、多维测量数据与单一观测面数据目标体探测能力进行对比分析。空间域和频率域数据分析技术显示梯度异常对浅部地质体具有较高的分辨能力,而重磁异常包含更多的低频信息。多维重磁数据分析结果显示多维数据可用于反演图像重构的奇异值幅值增加,能够更好地刻画地质体的三维空间结构。通过量化不同测量组合参数、数据类型和空间分布对地质目标解释结果的影响分析,可以看出建立综合高精度、多维观测张量数据的反演方法在地质体解释中具有更大的优势。本文充分挖掘数据中包含的目标体信息,强调多维数据的充分利用,综合更多梯度分量参与反演计算,进而建立快速、高精度、人为干预少的地质体参数综合解释方法体系。目标体参数包括水平位置、埋深及物性分布等。本文从地质体参数解释展开新方法的研究,建立张量数据高精度边界识别方法;发展多维重磁数据场源位置快速反演及成像方法;研究几何参数优化约束的物性反演方法。边界识别方法广泛应用于地质体水平位置的确定中,传统方法存在识别位置模糊及深部构造水平位置偏差较大的问题。本文利用梯度数据能够更好地反映短波长信息的特点,提出几种基于梯度张量的地质体边界增强方法。首先对现有曲率梯度张量特征值边界识别方法进行深浅模型体异常试验,发现该方法无法同时识别不同幅值异常体的边界。为了解决这一问题,本文提出基于方向梯度张量矩阵特征值的边界识别方法,采用垂向导数均衡深浅地质体异常并应用到合成及实测数据中,新方法可同时识别强弱异常的目标体边界,但识别的边界位置模糊。为了改善这一问题,本文将用于图像处理的构造张量引入到重磁数据处理中,定义基于构造张量的方向总水平导数边界识别方法,可更清晰地增强边界位置。通过合成和实测重磁数据试验验证了方法的有效性。另外,为了改善现有方法在深部地质构造水平位置识别的偏移问题,建立了基于三维构造张量边界识别方法,综合利用三维构造张量的分量使得边界位置圈定的更加准确,实测数据的应用表明新方法可更准确的获得构造断裂的水平位置,发现更多可靠的细节构造信息。重磁勘探相比其它地球物理方法纵向分辨率不足,但快速、可靠地获取目标体的深度信息仍是重磁数据解释的重要任务。针对现有深度估计方法存在对先验信息依赖性强、计算精度低、计算复杂的问题,本文建立基于多维数据的高精度深度反演方法。通过理论公式推导提出传统斜导数深度估计方法计算的改进策略,有效地提高了水平分辨能力。为了实现不同场源类型地质体的深度计算,基于多维重力异常建立广义斜导数深度估计方法,理论模型和实测数据表明该方法可快速、准确地反演出场源深度和形体参数。针对现有解析信号法大多需要构造指数先验信息或计算高阶导数的不足,本文提出多维数据解析信号法,建立埋深与构造指数关于三个观测面高度的计算关系式,合成及实测磁异常反演结果证明新方法可在压制噪声干扰的同时获得可靠的场源几何参数。归一化总梯度方法实质是不同下延深度位置异常解析信号的归一化。为了同时拾取出地质体的水平及深度位置特征,本文对方向总水平导数滤波器归一化得到归一化边界识别方法。合成及实测数据应用效果表明新方法三维成像结果指示地质体空间位置特征。边界识别、深度估计及成像方法可获得地质目标的空间位置,通常将其称为几何参数反演方法,而另一重要的反演以获得物性参数为目的。只有综合了几何参数和物性参数才能完成完整的重磁异常解释。本文选择聚焦反演方法求解重磁数据物性反演问题,简单模型试验显示聚焦反演有效地反演出陡峭边界地质体的物性分布,而不同埋深模型体同时反演时较大埋深异常体反演结果误差较大。为了提高反演结果的分辨率,将几何参数反演获得的目标体水平位置、深度信息通过深度加权函数及水平梯度加权函数约束到物性反演中,有效地提高了反演的横向与深度分辨率。本文提出的优化约束物性反演方法,即将几何参数反演获得的空间位置信息作为先验条件约束到物性反演过程中,在其它地球物理或地质条件缺乏的情况下获得更可靠的物性分布。本文建立的重磁数据地质目标位置及物性高精度解释方法旨在充分挖掘数据本身包含的异常体信息,更好地描述了地质体的空间位置及物性分布特征,合成及实测数据的成功应用说明方法具有很好的实用价值。
二、区域重磁资料解释中的相关分析方法（论文开题报告）
（1）论文研究背景及目的
此处内容要求：
首先简单简介论文所研究问题的基本概念和背景，再而简单明了地指出论文所要研究解决的具体问题，并提出你的论文准备的观点或解决方法。
写法范例：
本文主要提出一款精简64位RISC处理器存储管理单元结构并详细分析其设计过程。在该MMU结构中,TLB采用叁个分离的TLB,TLB采用基于内容查找的相联存储器并行查找,支持粗粒度为64KB和细粒度为4KB两种页面大小,采用多级分层页表结构映射地址空间,并详细论述了四级页表转换过程,TLB结构组织等。该MMU结构将作为该处理器存储系统实现的一个重要组成部分。
（2）本文研究方法
调查法：该方法是有目的、有系统的搜集有关研究对象的具体信息。
观察法：用自己的感官和辅助工具直接观察研究对象从而得到有关信息。
实验法：通过主支变革、控制研究对象来发现与确认事物间的因果关系。
文献研究法：通过调查文献来获得资料，从而全面的、正确的了解掌握研究方法。
实证研究法：依据现有的科学理论和实践的需要提出设计。
定性分析法：对研究对象进行“质”的方面的研究，这个方法需要计算的数据较少。
定量分析法：通过具体的数字，使人们对研究对象的认识进一步精确化。
跨学科研究法：运用多学科的理论、方法和成果从整体上对某一课题进行研究。
功能分析法：这是社会科学用来分析社会现象的一种方法，从某一功能出发研究多个方面的影响。
模拟法：通过创设一个与原型相似的模型来间接研究原型某种特性的一种形容方法。
三、区域重磁资料解释中的相关分析方法（论文提纲范文）
（1）重力方法在地壳结构研究中的应用（论文提纲范文）
| 摘要 |
| 第一章 引言 |
| 第二章 地壳结构研究中的重力学方法 |
| 2.1 重力异常的数据处理 |
| 2.2 位场反问题解法 |
| 2.3 重力均衡原理与应用 |
| 2.4 卫星重力场与地球深部密度不均匀性 |
| 附录A线性变密度四面体重力异常计算 |
| 第三章 地震-重力综合解释与联合反演 |
| 3.1 速度-密度关系 |
| 3.2 地震-重力综合解释 |
| 3.3 地震-重力的联合反演 |
| 3.4 小结 |
| 第四章 新疆西北部地壳结构与构造动力学特征的重力解释 |
| 4.1 新疆西北部地区的区域地质与大地构造格架 |
| 4.2 重力资料和航磁资料 |
| 4.3 重力与航磁资料的处理和解释 |
| 4.4 卫星重力异常与新疆地区的构造动力学特征 |
| 4.5 地热资料与解释 |
| 4.6 讨论 |
| 4.7 小结 |
| 第五章 腾冲火山区地壳结构的地震-重力综合反演与解释 |
| 5.1 腾冲火山区地质构造概况与已有的深部地球物理研究成果 |
| 5.2 地震资料及其解释 |
| 5.3 重力资料的数据处理与反演 |
| 5.4 腾冲火山区地壳结构综合解释 |
| 5.5 小结 |
| 参考文献 |
| 致谢 |
| 攻读博士学位期间发表的论文 |
（2）基于重磁数据的冈底斯成矿带深部结构分析与成矿预测（论文提纲范文）
| 摘要 |
| abstract |
| 1 绪论 |
| 1.1 研究背景及意义 |
| 1.2 研究现状及存在的问题 |
| 1.2.1 青藏高原东西分区的地球物理特征研究现状 |
| 1.2.2 冈底斯岩浆作用的地球动力学背景研究现状 |
| 1.2.3 重磁数据处理方法研究现状 |
| 1.2.4 重磁三维反演研究现状 |
| 1.3 研究思路和研究内容 |
| 1.3.1 研究思路 |
| 1.3.2 研究内容 |
| 1.4 创新点 |
| 1.5 章节安排 |
| 2 研究区地质概况与地球物理特征 |
| 2.1 冈底斯带构造-岩浆演化 |
| 2.2 区域地质概况 |
| 2.2.1 区域地层 |
| 2.2.2 区域构造-岩浆岩 |
| 2.3 地球物理特征 |
| 2.3.1 岩石物理特征与重磁异常成因分析 |
| 2.3.2 重磁场特征 |
| 2.4 深部结构特征 |
| 2.4.1 地震和MT揭示的深部结构 |
| 2.4.2 莫霍面深度 |
| 2.4.3 藏南壳内低速高导体 |
| 2.5 区域矿产及时空分布 |
| 2.6 本章小结 |
| 3 研究区重磁异常二维经验模态分解 |
| 3.1 重磁异常二维经验模态分解理论 |
| 3.1.1 二维经验模态分解基本原理 |
| 3.1.2 功率谱分析与场源深度估计 |
| 3.1.3 模型试验 |
| 3.2 冈底斯带及邻区重磁异常二维经验模态分解及深度估计 |
| 3.3 结果讨论 |
| 3.3.1 东西分区的重磁证据 |
| 3.3.2 中下地壳低密度区分布 |
| 3.4 本章小结 |
| 4 研究区地质体边界与断裂识别 |
| 4.1 边界识别改进方法研究 |
| 4.2 研究区边界识别方法选取与实际效果 |
| 4.3 研究区断裂与地质体边界推断 |
| 4.3.1 研究区重磁异常断裂的识别特征 |
| 4.3.2 研究区断裂与地质体边界 |
| 4.3.3 研究区局部地质体推断解释 |
| 4.4 本章小结 |
| 5 研究区重磁三维反演 |
| 5.1 基于先验信息约束的重磁反演研究 |
| 5.1.1 基于先验信息约束的重磁反演流程 |
| 5.1.2 重磁2.5维人机交互式拟合 |
| 5.1.3 重磁三维反演 |
| 5.2 冈底斯带及邻区的重磁反演 |
| 5.2.1 冈底斯带及邻区的重磁2.5维交互式反演 |
| 5.2.2 冈底斯带及邻区的重磁三维反演 |
| 5.3 结果与讨论 |
| 5.3.1 岩浆岩三维形态与分布特征 |
| 5.3.2 冈底斯带岩浆活动的地球动力学过程 |
| 5.4 本章小结 |
| 6 成矿远景区划分 |
| 6.1 重磁成果与矿产分布关系分析 |
| 6.1.1 重磁异常特征与矿床分布的关系 |
| 6.1.2 地质体边界和断裂与矿床分布的关系 |
| 6.1.3 深部密度、磁化率结果与矿床分布的关系 |
| 6.2 区域成矿预测 |
| 6.3 本章小结 |
| 7 结论与建议 |
| 7.1 结论 |
| 7.2 存在问题及建议 |
| 7.2.1 主要问题 |
| 7.2.2 建议 |
| 致谢 |
| 参考文献 |
| 附录 |
（3）克拉玛依后山地区重磁场及其与构造格架关系研究（论文提纲范文）
| 作者简介 |
| 摘要 |
| ABSTRACT |
| 第一章 绪论 |
| §1.1 研究背景、目的和意义 |
| §1.2 国内外研究现状和存在的问题 |
| 1.2.1 西准噶尔地区的地质研究现状 |
| 1.2.2 西准噶尔地区的地球物理工作现状 |
| 1.2.3 重磁资料解释方法的研究进展 |
| 1.2.4 卫星重磁技术的研究进展 |
| 1.2.5 值得探讨的几个问题 |
| §1.3 研究目标、内容、方法和技术路线 |
| 1.3.1 研究目标和研究内容 |
| 1.3.2 研究方法和技术路线 |
| §1.4 论文主要成果与创新之处 |
| §1.5 论文章节安排 |
| 第二章 大地构造及区域地质背景 |
| §2.1 西准噶尔及邻区大地构造及地质背景 |
| 2.1.1 西准噶尔及邻区大地构造背景 |
| 2.1.2 西准噶尔及邻区区域地质背景 |
| §2.2 克拉玛依后山地区构造与地质背景 |
| 2.2.1 克拉玛依后山地区主要构造单元及属性 |
| 2.2.2 克拉玛依后山地区蛇绿混杂岩与地层分布 |
| §2.3 克拉玛依后山地区主要地质单元的岩石物性特征 |
| §2.4 小结 |
| 第三章 重磁数据处理与解释方法 |
| §3.1 重磁异常转换处理方法 |
| 3.1.1 磁异常模量 |
| 3.1.2 重磁异常梯度 |
| §3.2 重磁异常分离方法 |
| 3.2.1 小波多尺度分解 |
| 3.2.2 经验模态多尺度分解 |
| 3.2.3 基于多尺度分解分离重磁异常的对比分析 |
| §3.3 重磁异常三维反演方法 |
| §3.4 小结 |
| 第四章 西准噶尔及邻区卫星重磁异常特征与大地构造 |
| §4.1 卫星重磁资料来源与解算 |
| 4.1.1 卫星重力资料来源与解算 |
| 4.1.2 卫星磁力资料来源与解算 |
| §4.2 卫星重磁异常分布特征及其大地构造意义 |
| 4.2.1 卫星重力异常分布特征与大地构造 |
| 4.2.2 卫星磁力异常分布特征与大地构造 |
| §4.3 岩石圈三维密度与磁化率结构特征及其大地构造意义 |
| 4.3.1 卫星重磁异常三维反演 |
| 4.3.2 岩石圈三维密度与磁化率结构特征 |
| 4.3.3 岩石圈三维密度与磁性结构的大地构造含义及演化启示 |
| §4.4 小结 |
| 第五章 克拉玛依后山地区区域重磁异常处理与解释 |
| §5.1 区域重力与航磁资料来源与汇编 |
| §5.2 区域重力与航磁异常处理 |
| 5.2.1 重磁异常常规方法处理结果 |
| 5.2.2 多尺度分解方法的应用与评价 |
| 5.2.3 重磁异常线性信号的提取 |
| 5.2.4 浅源与深源重磁异常 |
| §5.3 上地壳三维密度与磁化率结构特征及其地质构造意义 |
| 5.3.1 重磁异常三维反演 |
| 5.3.2 上地壳三维密度与磁化率结构特征 |
| 5.3.3 克拉玛依后山地区的地质构造格架 |
| §5.4 小结 |
| 第六章 克拉玛依后山地区高精度重磁与MT剖面综合研究 |
| §6.1 地面高精度重磁剖面数据 |
| §6.2 地面高精度重磁剖面数据反演 |
| §6.3 上地壳密度与磁化率剖面结构特征 |
| §6.4 综合地球物理剖面解释 |
| 6.4.1. Ⅰ和Ⅱ线综合大剖面 |
| 6.4.2. Ⅲ线综合大剖面 |
| 6.4.3. Ⅳ线综合大剖面 |
| §6.5 小结 |
| 第七章 结论与建议 |
| §7.1 结论 |
| 7.1.1 西准噶尔及邻区卫星重磁异常特征与大地构造意义 |
| 7.1.2 克拉玛依后山地区区域重力与航磁异常特征与地质构造分布 |
| 7.1.3 克拉玛依后山地区高精度剖面重磁异常特征与综合解释 |
| §7.2 建议 |
| 致谢 |
| 参考文献 |
| 附录 |
| 附录A：克拉玛依后山地区主要地层的岩石物性特征 |
| 附录B：克拉玛依后山地区主要岩体的岩石物性特征 |
| 附录C：克拉玛依后山地区主要蛇绿岩带的岩石物性特征 |
（4）岩性识别技术现状与进展（论文提纲范文）
| 0 引言 |
| 1 岩性识别的主要方法 |
| 1.1 重磁技术 |
| 应用实例 |
| 1.2 地震技术 |
| 应用实例 |
| 1.3 测井技术 |
| 应用实例 |
| 1.4 遥感技术 |
| 应用实例 |
| 1.5 综合地球物理技术 |
| 2 岩性识别技术对比 |
| 3 结论与展望 |
（5）金川铜镍矿区地球物理特征及深部成矿预测（论文提纲范文）
| 摘要 |
| Abstract |
| 第1章 绪论 |
| 1.1 研究背景及意义 |
| 1.2 研究动态及发展现状 |
| 1.2.1 深部金属矿床地球物理探测进展 |
| 1.2.2 位场数据处理方法概述 |
| 1.2.3 金川铜镍矿床深部及外围找矿现状 |
| 1.3 研究内容及创新点 |
| 1.3.1 论文主要研究内容 |
| 1.3.2 论文主要创新点 |
| 第2章 区域地质背景 |
| 2.1 地层 |
| 2.2 岩浆活动 |
| 2.3 区域构造 |
| 2.3.1 构造特征 |
| 2.3.2 区域构造演化 |
| 2.4 地球物理场特征 |
| 2.4.1 区域重力场特征 |
| 2.4.2 区域磁场特征 |
| 第3章 矿床地质特征 |
| 3.1 矿区地质 |
| 3.1.1 地层及其与成矿的关系 |
| 3.1.2 构造及其与成矿的关系 |
| 3.1.3 岩浆岩及其与成矿的关系 |
| 3.1.4 铜镍矿成因模式 |
| 3.2 矿体特征 |
| 3.2.1 矿体类型、赋存位置及规模 |
| 3.2.2 矿石特征 |
| 3.3 岩、矿石物性特征研究 |
| 3.3.1 密度参数 |
| 3.3.2 磁性参数 |
| 3.3.3 电阻率特征 |
| 3.3.4 岩、矿石综合物性特征 |
| 第4章 位场数据处理方法与矿区航磁异常信息提取研究 |
| 4.1 位场分离方法研究 |
| 4.1.1 小波分析与经验模态分解多尺度分离方法 |
| 4.1.2 基于双树复小波的异常多尺度分离 |
| 4.1.3 基于形态学滤波的异常多尺度分离 |
| 4.1.4 理论模型试验 |
| 4.2 位场数据边界识别方法研究 |
| 4.2.1 基于形态学的基本边界识别算子 |
| 4.2.2 边界识别均衡滤波器 |
| 4.2.3 理论模型试验 |
| 4.3 基于SL0范数约束及GPU并行计算的聚焦反演 |
| 4.3.1 基于SL0范数紧支撑聚焦反演的原理 |
| 4.3.2 等效格架与GPU并行计算 |
| 4.3.3 理论模型试验 |
| 4.4 金川铜镍矿区航磁异常特征与信息提取研究 |
| 4.4.1 金川铜镍矿区航磁异常基本特征 |
| 4.4.2 矿区磁异常边界与磁源信息提取研究 |
| 4.5 本章小结 |
| 第5章 金川铜镍矿地质-地球物理找矿模型研究 |
| 5.1 地质找矿标志 |
| 5.2 地球物理深部探测能力研究及其找矿标志 |
| 5.3 找矿模型建立 |
| 第6章 铜镍矿深部及周边成矿预测 |
| 6.1 矿区周边航磁异常找矿意义分类 |
| 6.2 研究区深部成矿预测 |
| 6.3 小结 |
| 第7章 结论与展望 |
| 7.1 主要认识和结论 |
| 7.2 问题与展望 |
| 参考文献 |
| 作者简介及科研成果 |
| 致谢 |
（6）利用重磁资料进行构造边界识别与弱异常提取的方法研究及应用（论文提纲范文）
| 摘要 |
| Abstract |
| 第1章 绪论 |
| 1.1 选题依据和研究思路 |
| 1.2 研究的目的和意义 |
| 1.3 研究现状与进展 |
| 1.4 主要内容和创新点 |
| 1.4.1 主要研究内容 |
| 1.4.2 创新点 |
| 第2章 位场分离和边界增强基本理论 |
| 2.1 边界增强 |
| 2.1.1 水平方向导数 |
| 2.1.2 欧拉反褶积 |
| 2.1.3 小子域滤波 |
| 2.2 人工神经网络方法 |
| 2.2.1 神经网络模型 |
| 2.2.2 人工神经网络在重磁异常解释中的应用 |
| 2.3 小波分析 |
| 2.3.1 小波变换的概念 |
| 2.3.2 小波分析的基本理论 |
| 2.3.3 小波分析的应用 |
| 第3章 非线性理论在重磁边界识别中的应用 |
| 3.1 混沌动力学 |
| 3.1.1 混沌理论 |
| 3.1.2 混沌与成矿的关系 |
| 3.1.3 混沌的应用 |
| 3.2 Lyapunov指数 |
| 3.2.1 Lyapunov指数的定义 |
| 3.2.2 Lyapunov指数计算方法 |
| 3.3 Lyapunov指数与混沌、成矿的关系 |
| 3.3.1 Lyapunov指数与混沌的关系 |
| 3.3.2 Lyapunov指数与成矿的关系 |
| 3.4 Lyapunov指数在边界识别中的应用 |
| 3.4.1 双方向Lyapunov指数极小值方法 |
| 3.4.2 理论模型 |
| 3.4.3 实例 |
| 3.5 分形理论及其应用 |
| 3.5.1 分形几何基本理论 |
| 3.5.2 分形的应用 |
| 第4章 基于自相关滤波法提取重磁弱异常的方法研究 |
| 4.1 自相关滤波法研究现状及理论基础 |
| 4.1.1 自相关滤波法研究现状 |
| 4.1.2 自相关滤波法理论 |
| 4.2 一维自相关滤波 |
| 4.2.1 一维自相关滤波法的原理 |
| 4.2.2 一维自相关滤波法理论模型试验 |
| 4.2.3 一维自相关滤波法的改善及试验 |
| 4.3 维自相关滤波 |
| 4.3.1 二维自相关滤波法原理 |
| 4.3.2 二维自相关滤波法理论模型试验 |
| 4.3.3 二维自相关滤波法的改善及试验 |
| 4.4 实际资料处理 |
| 第5章 结论与建议 |
| 5.1 结论 |
| 5.2 建议 |
| 参考文献 |
| 攻读博士学位期间发表的学术论文 |
| 致谢 |
（7）位场（重&磁）及其梯度异常自动解释方法研究（论文提纲范文）
| 摘要 |
| Abstract |
| 第1章 绪论 |
| 1.1 研究目的及意义 |
| 1.2 重磁场自动解释方法概述 |
| 1.2.1 重磁场分离技术 |
| 1.2.2 异常场源边界识别方法 |
| 1.2.3 重磁异常场源几何参数自动解释方法 |
| 1.2.4 场源物性反演及相关成像法 |
| 1.2.5 重磁场全张量数据解释 |
| 1.3 现有解释方法所存在的主要问题 |
| 1.4 论文研究内容及思路 |
| 1.4.1 研究内容 |
| 1.4.2 研究思路 |
| 1.5 论文创新点 |
| 第2章 优化小子域滤波法 |
| 2.2 优化小子域滤波法 |
| 2.2.1 子域剖分方式 |
| 2.2.2 子域判别准则 |
| 2.3 理论模型试验 |
| 2.4 四川盆地航磁异常应用 |
| 2.5 小结 |
| 第3章 重磁异常场源边界识别方法研究 |
| 3.1 相关系数法 |
| 3.1.1 相关系数法基本原理 |
| 3.1.2 模型试验 |
| 3.1.3 相关系数法实际应用效果 |
| 3.2 归一化总水平导数法 |
| 3.2.1 归一化总水平导数法基本原理 |
| 3.2.2 模型重力异常试验 |
| 3.2.3 实际应用效果及讨论 |
| 3.3 增强型均衡滤波器 |
| 3.3.1 增强型均衡滤波器原理 |
| 3.3.2 理论模型试验 |
| 3.3.3 实际应用效果 |
| 3.4 增强型局部相位滤波器 |
| 3.4.1 增强型局部相位滤波器基本原理 |
| 3.4.2 理论模型试验 |
| 3.4.3 实际数据应用效果 |
| 3.5 优化场源边界识别方法 |
| 3.6 边界识别新方法的 Geosoft 独立插件实现 |
| 3.7 小结 |
| 第4章 重磁异常场源几何参数自动解释方法研究 |
| 4.1 扩展型欧拉反褶积法 |
| 4.1.1 扩展型欧拉反褶积法原理 |
| 4.1.2 理论模型试验 |
| 4.1.3 实际数据处理 |
| 4.1.4 小结 |
| 4.2 优化解析信号法 |
| 4.2.1 导数解析信号法 |
| 4.2.2 和解析信号法 |
| 4.2.3 直接解析信号法 |
| 4.2.4 小结 |
| 4.3 优化局部波数法 |
| 4.3.1 多重位置局部波数法 |
| 4.3.2 快速局部波数法 |
| 4.3.3 小结 |
| 第5章 场源物性反演及相关成像法研究 |
| 5.1 基于快速模拟退火法的物性反演 |
| 5.1.1 方法基本原理 |
| 5.1.2 理论模型试验 |
| 5.1.3 吉林辽源地区煤矿采空区应用 |
| 5.2 增强型位场相关成像法 |
| 5.2.1 方法基本原理 |
| 5.2.2 理论模型试验 |
| 5.2.3 实测磁异常数据应用 |
| 5.3 小结 |
| 第6章 重磁场全张量数据解释方法研究 |
| 6.1 全张量测量技术 |
| 6.1.1 全张量测量简介 |
| 6.1.2 全张量数据的信息优势 |
| 6.1.3 实测张量数据与计算张量数据对比分析 |
| 6.2 张量局部波数法 |
| 6.2.1 张量局部波数法基本原理 |
| 6.2.2 理论模型试验 |
| 6.2.3 实际重力张量数据应用 |
| 6.3 方向解析信号法 |
| 6.3.1 方向解析信号法原理 |
| 6.3.2 模型试验 |
| 6.3.3 实际数据处理 |
| 6.4 重磁全张量数据解释方法的 Geosoft 独立插件实现 |
| 6.5 小结 |
| 第7章 结论 |
| 参考文献 |
| 附录Ⅰ：快速局部波数法 |
| 附录Ⅱ：基于 Hartley 变换的重磁异常导数计算公式 |
| 作者简介及读博期间所取得的科研成果 |
| 后记和致谢 |
（8）地球物理（重磁震）勘探数据一体化解释技术（论文提纲范文）
| 摘要 |
| Abstract |
| 第1章 绪论 |
| 1.1 研究目的及意义 |
| 1.2 相关研究方法及发展进程 |
| 1.2.1 数学形态学在地球物理数据解释中的应用 |
| 1.2.2 地震同相轴及位场边界识别技术研究现状 |
| 1.2.3 地层属性参数获取方法研究现状 |
| 1.2.4 随机反演地球物理解释中的应用现状 |
| 1.2.5 梯度异常界面反演及地球物理数据联合反演 |
| 1.3 现有方法所存在的主要问题 |
| 1.4 本论文研究内容和技术路线 |
| 1.4.1 研究内容 |
| 1.4.2 技术路线 |
| 1.5 论文创新点 |
| 第2章 数学形态学法在地震数据同相轴增强及重磁数据边界识别中的应用 |
| 2.1 数学形态学基本运算 |
| 2.2 数学形态学在地震同相轴增强中的应用 |
| 2.2.1 增强型数学形态学法基本原理 |
| 2.2.2 理论模型实验 |
| 2.2.3 实际地震剖面同相轴增强 |
| 2.3 形态学方法在重磁数据边界识别中的应用 |
| 2.3.1 形态学边界识别滤波器基本原理 |
| 2.3.2 形态学边界识别滤波器模型试验 |
| 2.3.3 实测重磁数据应用效果 |
| 2.3.4 改进的形态学滤波器 |
| 2.4 小结 |
| 第3章 地层属性(密度、磁化率)及几何参数反演 |
| 3.1 一种利用重力异常定量解释断层的新方法 |
| 3.1.1 方法基本原理 |
| 3.1.2 模型实验 |
| 3.1.3 实际资料处理 |
| 3.2 改进的解析信号法在磁异常解释中的应用 |
| 3.2.1 方法基本原理 |
| 3.2.2 理论模型试验 |
| 3.2.3 实际磁异常解释 |
| 3.3 改进的局部波数法及其在磁测数据解释的应用 |
| 3.3.1 方法基本原理 |
| 3.3.2 模型试验 |
| 3.3.3 实际数据处理 |
| 3.4 小结 |
| 第4章 地震随机反演及不确定性分析 |
| 4.1 地震随机反演介绍 |
| 4.1.1 随机反演定义 |
| 4.1.2 确定性反演与随机反演区别 |
| 4.1.3 地震随机反演基本原理 |
| 4.2 随机反演相关参数概念 |
| 4.2.1 变差函数 |
| 4.2.2 克里金技术 |
| 4.2.3 随机模拟方法 |
| 4.3 不确定性分析 |
| 4.4 地震随机反演及不确定性分析应用 |
| 4.4.1 理论数据反演 |
| 4.4.2 实际地震剖面数据随机反演及不确定性分析 |
| 4.5 小结 |
| 第5章 利用快速模拟退火算法进行重震联合反演 |
| 5.1 梯度界面反演理论基础 |
| 5.1.1 梯度测量优势 |
| 5.1.2 界面起伏引起的重力梯度异常 |
| 5.2 梯度数据在界面反演中的应用研究 |
| 5.2.1 快速模拟退火算法重力梯度异常反演流程 |
| 5.2.2 模型试验 |
| 5.2.3 中国南海海底地形计算 |
| 5.3 快速模拟退火法进行重震数据联合反演 |
| 5.3.1 联合反演方法 |
| 5.3.2 理论模型试验 |
| 5.3.3 实际数据应用 |
| 5.4 小结 |
| 第6章 结论 |
| 参考文献 |
| 作者简介及在读期间所取得的科研成果 |
| 后记和致谢 |
（9）重磁数据处理解释技术在华南地区岩体圈定与形态反演中的应用研究（论文提纲范文）
| 摘要 |
| Abstract |
| 第1章 绪论 |
| 1.1 选题依据及目的意义 |
| 1.2 国内外研究现状 |
| 1.2.1 重磁数据处理和解释技术研究现状 |
| 1.2.2 基于重磁数据的岩浆岩体研究现状 |
| 1.2.3 存在问题 |
| 1.3 主要研究内容和创新点 |
| 1.3.1 主要研究内容 |
| 1.3.2 创新点 |
| 第2章 华南区域地质及地球物理特征 |
| 2.1 地层 |
| 2.2 岩浆岩 |
| 2.2.1 侵入岩 |
| 2.2.2 火山岩 |
| 2.3 构造 |
| 2.4 矿产 |
| 2.5 区域地球物理特征 |
| 2.5.1 密度参数与重力场特征 |
| 2.5.2 磁性参数与磁场特征 |
| 2.5.3 华南花岗岩与重磁场的关系 |
| 第3章 华南地区资料收集与预处理 |
| 3.1 资料收集 |
| 3.2 数据整理 |
| 3.3 投影转换 |
| 3.4 网格化 |
| 3.5 其它处理 |
| 第4章 重磁数据处理方法 |
| 4.1 异常分离 |
| 4.2 边缘识别 |
| 4.2.1 垂向导数 |
| 4.2.2 总水平导数 |
| 4.2.3 解析信号 |
| 4.2.4 倾斜角 |
| 4.2.5 图 |
| 4.2.6 小子域滤波 |
| 4.2.7 实测资料处理 |
| 4.2.8 方法总结 |
| 4.3 数据扩边和空白区填补 |
| 4.3.1 数据扩边 |
| 4.3.2 空白区填补 |
| 4.3.3 关于扩边数选择的讨论 |
| 第5章 三维重磁反演技术 |
| 5.1 随机子域三维自动反演 |
| 5.1.1 重磁场正演一般公式 |
| 5.1.2 规则长方体重磁场公式 |
| 5.1.3 格架分离 |
| 5.1.4 格架等效计算 |
| 5.1.5 随机子域反演策略 |
| 5.2 三维人机交互正反演 |
| 5.2.1 三角形多面体模型正演 |
| 5.2.2 三角形多面体模型几何形状反演 |
| 第6章 三维模型可视化和转换技术 |
| 6.1 三维可视化图形工具 |
| 6.2 物性反演模型数据可视化 |
| 6.2.1 三维图形显示和拾取变换 |
| 6.2.2 等值面重构算法 |
| 6.3 模型数据转换 |
| 6.3.1 轮廓线的预处理 |
| 6.3.2 二维轮廓线的三维形体重构算法 |
| 6.3.3 轮廓线分支问题 |
| 6.4 程序功能 |
| 6.4.1 三维显示 |
| 6.4.2 模型转换 |
| 第7章 华南地区岩体圈定 |
| 7.1 局部异常分离 |
| 7.2 异常识别 |
| 7.3 隐伏岩体预测 |
| 7.4 岩体圈定 |
| 第8章 华南地区岩体形态三维反演解释 |
| 8.1 岩体异常选取 |
| 8.2 三维定量反演流程 |
| 8.2.1 岩体异常求取 |
| 8.2.2 岩体异常选择 |
| 8.2.3 围岩影响消除 |
| 8.2.4 三维岩体建模 |
| 8.3 华南地区岩体形态三维反演结果 |
| 附图表 |
| 第9章 结论与建议 |
| 9.1 结论 |
| 9.2 建议 |
| 参考文献 |
| 致谢 |
| 个人简介及在读期间发表论文情况 |
（10）重磁数据地质目标位置及物性高精度解释方法研究（论文提纲范文）
| 摘要 |
| Abstract |
| 第一章 绪论 |
| 1.1 研究背景及意义 |
| 1.2 国内外研究现状 |
| 1.2.1 异常场源边界识别方法 |
| 1.2.2 重磁数据地质体深度估计方法 |
| 1.2.3 重磁数据快速成像技术 |
| 1.2.4 重磁异常物性反演方法 |
| 1.3 存在的主要问题 |
| 1.4 研究内容及创新点 |
| 1.4.1 论文主要研究内容 |
| 1.4.2 论文主要创新点 |
| 第二章 重磁数据目标体探测能力分析技术 |
| 2.1 数据分析工具DPP和DRP |
| 2.1.1 理论基础 |
| 2.1.2 理论模型计算分析 |
| 2.2 DPP和DRP测量组合参数对解释结果影响分析 |
| 2.2.1 测网间距及测点分布影响分析 |
| 2.2.2 飞行高度影响分析 |
| 2.2.3 测量精度影响分析 |
| 2.3 重磁异常与其梯度异常目标体发现能力分析 |
| 2.3.1 DPP和DRP测量场类型影响分析 |
| 2.3.2 基于频谱分析技术测量场类型影响分析 |
| 2.4 多维测量数据目标体探测能力分析 |
| 2.4.1 空间域数据分析 |
| 2.4.2 基于DRP和DPP的数据分析 |
| 2.5 本章小结 |
| 第三章 张量数据高精度边界识别方法研究 |
| 3.1 基于曲率梯度张量特征值的边界识别方法 |
| 3.1.1 曲率梯度张量矩阵及特征值 |
| 3.1.2 理论模型异常处理 |
| 3.2 基于张量方向特征值的边界识别方法 |
| 3.2.1 方向梯度张量矩阵及特征值 |
| 3.2.2 模型异常试验 |
| 3.2.3 实测重磁数据处理 |
| 3.3 构造张量方向总水平导数边界识别方法 |
| 3.3.1 理论基础 |
| 3.3.2 合成模型异常处理 |
| 3.3.3 实测重磁数据处理效果 |
| 3.4 基于三维构造张量的边界识别方法 |
| 3.4.1 三维构造张量滤波器理论建立 |
| 3.4.2 合成异常数据处理 |
| 3.4.3 实测重磁数据处理效果 |
| 3.5 本章小结 |
| 第四章 多维重磁数据场源深度快速反演方法 |
| 4.1 传统斜导数深度估计方法及改进策略 |
| 4.1.1 斜导数深度估计方法理论基础 |
| 4.1.2 理论模型应用分析 |
| 4.2 基于多维数据的广义斜导数深度估计方法 |
| 4.2.1 简单模型的深度计算 |
| 4.2.2 多维数据广义斜导数深度反演理论 |
| 4.2.3 合成异常数据试验 |
| 4.2.4 实测数据应用效果 |
| 4.3 多维数据解析信号法场源深度反演 |
| 4.3.1 解析信号法理论基础 |
| 4.3.2 合成异常模型试验 |
| 4.3.3 实际数据应用 |
| 4.4 本章小结 |
| 第五章 地质体特征拾取的归一化边界识别方法 |
| 5.1 归一化总梯度法 |
| 5.1.1 理论基础 |
| 5.1.2 理论模型试验 |
| 5.2 稳定的向下延拓方法 |
| 5.2.1 积分迭代下延方法理论 |
| 5.2.2 理论模型试验 |
| 5.3 高分辨率的归一化边界识别方法 |
| 5.3.1 归一化方向总水平导数方法理论 |
| 5.3.2 合成异常处理试验 |
| 5.3.3 Humble盐丘重力数据应用 |
| 5.4 本章小结 |
| 第六章 位场数据物性反演的优化约束方法 |
| 6.1 位场正演理论基础 |
| 6.2 梯度全张量数据物性反演方法 |
| 6.2.1 基于最小梯度支撑方程稳定函数的聚焦反演 |
| 6.2.2 深度加权函数及物性约束 |
| 6.2.3 理论模型试验及分析 |
| 6.3 基于几何参数的优化约束物性反演 |
| 6.3.1 优化约束物性反演理论建立 |
| 6.3.2 合成异常应用效果 |
| 6.3.3 美国文顿岩丘实测航空重力梯度全张量数据反演 |
| 6.4 本章小结 |
| 第七章 结论 |
| 参考文献 |
| 作者简介及在学期间取得的科研成果 |
| 致谢 |
四、区域重磁资料解释中的相关分析方法（论文参考文献）
- [1]重力方法在地壳结构研究中的应用[D]. 楼海. 中国地震局地球物理研究所, 2001(01)
- [2]基于重磁数据的冈底斯成矿带深部结构分析与成矿预测[D]. 胡斌. 中国地质大学(北京), 2019(02)
- [3]克拉玛依后山地区重磁场及其与构造格架关系研究[D]. 许顺芳. 中国地质大学, 2015(03)
- [4]岩性识别技术现状与进展[J]. 付光明,严加永,张昆,胡浩,罗凡. 地球物理学进展, 2017(01)
- [5]金川铜镍矿区地球物理特征及深部成矿预测[D]. 张建民. 吉林大学, 2019(02)
- [6]利用重磁资料进行构造边界识别与弱异常提取的方法研究及应用[D]. 韩兆红. 吉林大学, 2011(09)
- [7]位场（重&磁）及其梯度异常自动解释方法研究[D]. 马国庆. 吉林大学, 2013(08)
- [8]地球物理（重磁震）勘探数据一体化解释技术[D]. 李丽丽. 吉林大学, 2013(08)
- [9]重磁数据处理解释技术在华南地区岩体圈定与形态反演中的应用研究[D]. 黄金明. 中国地质大学（北京）, 2013(10)
- [10]重磁数据地质目标位置及物性高精度解释方法研究[D]. 周帅. 吉林大学, 2017(09)