磁性界面反演方法
- 格式:docx
- 大小:117.66 KB
- 文档页数:2
地球物理反演方法的综述地球物理反演是一种利用地球物理方法来推断地下构造和物质分布的技术。
通过观测和测量地球物理场,如重力、地磁、电磁、地震等,结合数理统计和计算机模拟方法,可以对地下的地质构造、岩石性质和地下水资源等进行精确的推断。
本文将综述地球物理反演方法的原理、分类及应用。
一、地球物理反演方法的原理地球物理反演方法的原理在于根据地球物理场的观测数据,通过数学模型和计算方法,将地球物理场与地下介质属性之间的关系联系起来。
根据电磁波传播、物质密度、电阻率、磁化率等反演参数的变化规律,推断地下介质的结构和成分。
其中常用的地球物理反演方法包括重力法、磁法、电磁法、地电法和地震法等。
不同的反演方法适用于不同的地质介质和研究目标,各有其优势和限制。
二、地球物理反演方法的分类1. 重力反演法:利用重力场观测数据,通过计算物质的密度分布,来推断地下构造的方法。
重力反演法在石油勘探、地质灾害分析、水资源评价等领域具有广泛应用。
2. 磁法反演法:通过磁场观测数据,推断地下磁化率和磁性物质的空间分布。
磁法反演在矿产勘探、地震预测等方面发挥重要作用。
3. 电磁法反演法:通过电磁场观测数据,推断地下电阻率分布,来研究地下水资源、矿产和工程勘探。
电磁法反演在地下水资源评价、油气勘探、环境地球物理和岩土工程等方面有广泛应用。
4. 地电法反演法:通过电场和电位观测数据,推断地下电阻率分布,用于研究地下水位、地下水性质、污染监测和地下工程等。
地电法反演在工程地球物理勘探和水文地球物理领域具有广泛应用。
5. 地震法反演法:通过地震波在地下的传播与变化,推断地下介质的速度和密度分布,用于研究地质构造、地震预测和石油勘探等。
地震法反演是地球物理反演方法中应用最广泛的方法之一。
三、地球物理反演方法的应用地球物理反演方法广泛应用于地质探测、资源勘探、环境监测和工程勘察等领域。
以下是几个常见的应用领域:1. 石油勘探:地震反演方法可用于确定油气藏的位置、大小和分布,辅助油田开发和管理。
AbstractIn the early period of evaluation about basic structures, magnetic prospecting provided more convenient way in geophysical prospecting work, with taking advantage of its own characteristics of Less input, easy access to the field, wide range covering for unusual and effective information, data processing and shorter interpretation time. The processing and interpretation of magnetic data, the determination of the fluctuations of the magnetic geological interface and the depth of the magnetic basis, which is one of the important tasks for magnetic prospecting in in solving basin deep tectonic information. Studying the interface of magnetic strata has extremely important theoretical and practical value both in study the formation and evolution of the basin and the distribution of mineral, oil and gas resources, such as oil and gas formation, migration and occurrenceThe method of inversion magnetic interface depth can be divided into two categories: frequency domain and space domain,The spatial domain magnetic interface iterative inversion method that we studied in this paper is greatly influenced by the interface average depth of the proper selection, the deeper average interface you chose, the deeper entire interface the study will get, otherwise it will be more shallow. Therefore, it is important to determine the average depth of the interface before the inversion of the magnetic interface depth. The method of using power spectrum to determine the average depth of the interface are given in this paper, which has a big advantage of getting the result without magnetic interface material before the inversion. But before make confirmation of the initial interface value, We need to extend the observation surface field downward, Based on the iterative method proposed by the academician Xu Shizhe, this paper avoids the problem of divergence of field value effectively.Based on the theory of model test, which indicates that the magnetic interface in the spatial domain iterative inversion method is reliable, and It has certain pratical application value. At the same time, this method is used to deal with the magnetic survey data of northern south China sea, and obtains the characteristics of the magnetic basement and the fluctuation of the curie-surface in south China sea.Key words:Magnetic interface, the iterative inversion method, Power spectrum inversion method, the north of the South China Sea目录第一章绪论 (1)1.1选题背景及研究意义 (1)1.2国内外研究现状 (2)1.3主要研究内容及成果 (4)1.4技术路线 (5)第二章磁性界面正反演方法原理 (6)2.1磁性界面有限单元法正演 (6)2.2空间域磁性界面迭代反演法 (9)2.3磁性界面迭代反演的关键技术 (14)2.3.1界面平均深度的确定 (14)2.3.2迭代法位场向下延拓 (17)第三章磁性界面迭代反演方法程序设计 (19)3.1输入输出数据格式设计 (19)3.1.1平面位场数据输入格式 (19)3.1.2平面位场数据输出格式 (19)3.2磁性界面迭代反演方法程序设计 (20)3.3求取平均深度方法程序设计 (21)3.4位场向下延拓迭代法程序设计 (22)第四章理论模型试算 (24)4.1单一磁性界面模型反演 (24)4.2多磁性界面模型反演 (34)第五章实际资料处理及解释 (42)5.1研究区概况 (42)5.2岩石物性特征 (43)5.3磁性界面起伏反演 (46)第六章结论与建议 (53)6.1结论 (53)6.2建议 (53)参考文献 (55)攻读学位期间取得的研究成果 (58)致谢 (59)第一章绪论1.1 选题背景及研究意义自20世纪30年代起,磁法勘探开始应用于我国找矿试验工作中,此后随着地质工作的不断深入开展,以及数学物理理论和计算机技术的飞速发展,使得磁法勘探在方法技术、理论创新以及实际应用等各方面得到了快速发展,已成为现代地球物理方法中重要的技术手段之一。
第九节 磁性界面反演方法具有一定磁性差异的地质界面,如结晶基底面、大岩体的上顶面等,是找矿勘探与基础地质研究中常见的地质现象。
磁性界面反演方法是确定这一类地质模型界面深度的方法。
磁性界面反演方法有空间域和频率域两大类。
在磁性界面反演中,常常把磁性界面划分成大量的离散二度水平棱柱体或三度直立棱柱体组合模型,由于未知参数太多不能采用直接解法,往往采用迭代法或其他方法。
空间域方法未知参数多、计算时间长、效率低。
Parker (1972,1973,1974)采用了连续模型,得出了频率域重磁位场正反演的理论公式,Oldenburg (1974)把它推广成迭代形式并做了二维计算。
由于引入快速傅立叶变换,在相同精度下,频率域方法比空间域方法反演速度至少要快一个数量级以上。
因此,频率域磁性界面的反演方法成为界面反演的一种常用方法,用于区域磁测资料解释与油气勘探中研究基底构造。
一、磁性界面异常的正演如图7-9-1所示的磁性界面,其上下界面磁化强度差为M ,为简单起见,设M 垂直向下。
图7-9-1 磁性界面示意图若磁化率为常数,考虑n=0时,即泰勒展开式第一项在空间域为常数项,略去n=0项则有: 01()(,)2!n n Hs n s Z u v M e h n μ∞-=⎡⎤-∆=⎢⎥⎣⎦∑ (7-9-13) 上式表示,当给定了平均深度H 及平均深度上的起伏 (),h ξη,取泰勒展开式有限项数n=3~8,就可以计算出 n h和 (,)Z u v ∆,利用快速傅立叶变换即可得到空间域的磁异常值Z ∆(x,y,0)。
二、磁性界面异常的反演式(7-9-13)是磁性界面正演计算公式,稍作一下变化,就可以当作反演迭代公式。
我们把和式中n=1的项写出并移项得202()()!Hs n n n e s h Z h Ms n μ∞=-=-∆-∑ (7-9-14) 若令等式右端的h 值为第i 次近似值,而等式左端的h 为第(i+1)次近似值,即为Parker 提出的迭代反演公式: 1202()[()]!Hs n i n n e s h Z h i Ms n μ∞+=--=∆-∑ (7-9-15) 实践证明,这种迭代方式每次求的都是h 值,而不是h 的修正量h ∆,不容易保证迭代收敛。
磁法数据处理、异常反演与解释的常用方法及常见问题探讨张湖源(安徽省地质矿产勘查局313地质队)磁法勘探是最经典的物探方法,可广泛用于地质调查的各个阶段、工程地质及考古学等众多方面,尤其是在铁矿勘查中,更是必不可少的先行手段。
可以说没有其它的地球物理方法有如此广阔的应用范围,花费少而提供信息的丰富。
因磁参数多为矢量且常见干扰较多(与其它物探方法相比),使数据处理及异常推断解释变得较复杂,在实际应用中,产生不少使人困惑的问题。
结合笔者实际工作经验,对其进行初步探讨,供同行参考。
1磁法数据预处理常用方法对实测数据进行日变、基点、正常场等改正后,应注意消除异常数据的误差与干扰。
误差主要源于仪器的状态和操作及点位误差;干扰主要是指人文或地质因素的干扰。
在严格执行技术规范下,含有人文干扰的数据一般不作为成果;地质干扰通常指与勘探目标无关的地质因素引起的异常。
浅表局部的地质干扰体分为两类:①在空间上有一定分布规律的,如出露的岩石;②孤立的、无规律的,如滚石等。
孤立的地质干扰具有随机性,具有白噪声特征,而一些出露岩石虽然不具备随机特征,但往往具有相同的走向等特征,在某方向上具有有色噪声特点。
由于多数误差和干扰具有随机性特征,其均值为零,因此,可以通过小范围异常进行平均来消除这种误差和干扰。
平均圆滑能够有效地消除随机误差和干扰,但圆滑后有可能改变异常形态特征,给一些利用异常形态特征进行异常解释的工作带来困难,可采用多项式圆滑。
深部大型的地质干扰体引起异常特征为:磁性基底异常强度大,但相对平缓;岩浆岩(强磁性、较强磁性)有相当的强度,但有一定变化;火山岩或火山碎屑岩强度不大,但变化大。
对局部磁异常进行方差统计,方差较大被视为隐伏岩浆岩或火山岩异常,由此推断隐伏岩体的存在。
即可通过场的分离来剥离大型地质干扰体引起的区域异常。
剔处各种干扰后编绘图件时,成图数据位置最好使用xy坐标,以避免用点线号成图造成异常形态扭曲,尤其测线为斜线时更应注意。
利用迭代反演方法解决变磁性磁界面反演问题
张向宇;吴健生
【期刊名称】《物探化探计算技术》
【年(卷),期】2014(036)001
【摘要】南海地层岩石磁性在横向上呈现较大差异,而目前常用的磁界面反演方法多是针对常磁性磁界面进行的,若直接将其应用于南海这样磁性差异较大的地区,则反演结果中必定会带进很大误差.通过引入滑动时窗,对较常用的磁界面迭代反演方法进行研究和改进,获取适用于变磁性磁界面深度反演的方法,提供了一个可用于解决此种情况的方法手段.
【总页数】9页(P35-43)
【作者】张向宇;吴健生
【作者单位】广州海洋地质调查局,广州 510760;同济大学海洋地质国家重点实验室,上海 200092
【正文语种】中文
【中图分类】P318.2
【相关文献】
1.双界面模式Parker算法磁性界面正反演方法 [J], 相鹏;刘展
2.变磁性磁界面迭代反演方法模型计算 [J], 张向宇;吴健生
3.磁性下界面的正则化反演方法 [J], 刘展;潘作枢
4.一种二度体单一界面的空间域重力迭代反演方法 [J], 高艳
5.基于Padé有理展开的磁性界面正反演方法(英文) [J], 张冲;黄大年;张凯;浦义涛;于平
因版权原因,仅展示原文概要,查看原文内容请购买。
浅析磁异常反演计算方法的体会
磁异常反演计算方法是在资源勘查领域中应用比较广泛的一种方法。
它在利用磁场测量空间数据对示性地质单元变化的研究中有着独特的优势。
其基本方法是将地面测量到的垂直磁场信号投射到有效载荷空间,并表述出底部地表的垂直磁感应强度或磁域重建图。
磁异常反演计算方法的计算主要分为三个步骤。
第一步施加磁场载荷,即在磁异常反演计算空间中计算地表垂直磁场。
只要将测量得到的空间数据加在反演模型空间中就可以完成第一步,得到磁场贴图。
第二步计算磁异常,橘黄将在第一步中计算得到的磁贴图和初始模型界面比较,找出磁场偏差的分布以及应用介质的示性地质参数分布,得到地面测量磁异常地形图,从而引导下一步的勘探工作。
第三步骤地心强度计算,即拟合磁异常的地质单元结构,根据地性强度模型估算地心强度空间分布图,以合理反映地表示性单元的结构特征。
综上所述,磁异常反演计算方法具有低成本、有效利用、快速响应、可靠性高等特点,对于资源勘探及开发具有重要意义。
它能够有效有效预测地面负载,根据异常分布特征解离不同单元,改善勘测定位的精度和准确性,是资源勘探中客观贴近岩性结构的有力手段。
《重金属层-铁磁层界面反演对称性破缺诱导的零场SOT-MRAM》篇一重金属层-铁磁层界面反演对称性破缺诱导的零场SOT-MRAM一、引言随着信息技术的飞速发展,磁性随机存取存储器(MRAM)因其非易失性、高速度和低功耗等优点,在下一代存储技术中扮演着重要角色。
其中,基于自旋轨道力矩(SOT)的MRAM技术(SOT-MRAM)因其在零场下的快速切换和低功耗特性,备受关注。
然而,实现高性能的SOT-MRAM仍面临诸多挑战,其中之一便是如何克服重金属层/铁磁层界面反演对称性破缺对SOT-MRAM性能的影响。
本文将深入探讨这一问题的研究进展,并分析其潜在应用。
二、重金属层/铁磁层界面反演对称性破缺的概述在SOT-MRAM中,重金属层与铁磁层的界面是影响自旋轨道力矩性能的关键因素。
在传统结构中,反演对称性的保持对提高磁性材料性能至关重要。
然而,在实际应用中,由于材料生长过程中的缺陷、杂质以及界面结构的不对称性等因素,往往导致界面处反演对称性的破缺。
这种破缺会导致磁化动力学特性的改变,从而影响SOT-MRAM的性能。
三、反演对称性破缺的诱导效应与优化策略针对重金属层/铁磁层界面反演对称性破缺的问题,本文提出了一种新的优化策略。
通过引入特定的材料或结构,如特定形状的纳米颗粒或特殊的合金材料,来调节界面处的反演对称性。
同时,利用第一性原理计算和实验手段相结合的方法,深入研究了这些优化策略对SOT-MRAM性能的影响。
四、零场SOT-MRAM的优化与性能提升在零场下实现SOT-MRAM的高效读写操作是研究的重要目标之一。
通过对重金属层/铁磁层界面的精确控制与优化,我们可以有效地提高SOT-MRAM的切换速度和降低功耗。
本文研究了界面反演对称性破缺对零场下SOT-MRAM性能的影响,并提出了相应的优化方案。
通过实验验证,这些方案显著提高了SOT-MRAM的读写速度和稳定性。
五、潜在应用与展望基于上述研究,我们可以预见SOT-MRAM在高性能计算和嵌入式存储等领域中的潜在应用。
第九节 磁性界面反演方法
具有一定磁性差异的地质界面,如结晶基底面、大岩体的上顶面等,是找矿勘探与基础地质研究中常见的地质现象。
磁性界面反演方法是确定这一类地质模型界面深度的方法。
磁性界面反演方法有空间域和频率域两大类。
在磁性界面反演中,常常把磁性界面划分成大量的离散二度水平棱柱体或三度直立棱柱体组合模型,由于未知参数太多不能采用直接解法,往往采用迭代法或其他方法。
空间域方法未知参数多、计算时间长、效率低。
Parker (1972,1973,1974)采用了连续模型,得出了频率域重磁位场正反演的理论公式,Oldenburg (1974)把它推广成迭代形式并做了二维计算。
由于引入快速傅立叶变换,在相同精度下,频率域方法比空间域方法反演速度至少要快一个数量级以上。
因此,频率域磁性界面的反演方法成为界面反演的一种常用方法,用于区域磁测资料解释与油气勘探中研究基底构造。
一、磁性界面异常的正演
如图7-9-1所示的磁性界面,其上下界面磁化强度差为M ,为简单起见,设M 垂直向下。
图7-9-1 磁性界面示意图
若磁化率为常数,考虑n=0时,即泰勒展开式第一项在空间域为常数项,略去n=0项则有: 01()(,)2!n n Hs n s Z u v M e h n μ∞-=⎡⎤-∆=⎢⎥⎣⎦
∑ (7-9-13) 上式表示,当给定了平均深度H 及平均深度上的起伏 (),h ξη,取泰勒展开式
有限项数n=3~8,就可以计算出 n h
和 (,)Z u v ∆,利用快速傅立叶变换即可得到空间域的磁异常值Z ∆(x,y,0)。
二、磁性界面异常的反演
式(7-9-13)是磁性界面正演计算公式,稍作一下变化,就可以当作反演迭代公式。
我们把和式中n=1的项写出并移项得
2
02()()!Hs n n n e s h Z h Ms n μ∞=-=-∆-∑ (7-9-14) 若令等式右端的h 值为第i 次近似值,而等式左端的h 为第(i+1)次近似值,即为Parker 提出的迭代反演公式: 12
02()[()]!Hs n i n n e s h Z h i Ms n μ∞+=--=∆-∑ (7-9-15) 实践证明,这种迭代方式每次求的都是h 值,而不是h 的修正量h ∆,不容易保证迭代收敛。
可以采用如下的迭代方式:
取式(7-9-15)的一级近似值: (1)02Hs e h Z Ms
μ=-∆ (7-9-16) 其中 (1)h 表示h 的一级近似值(1)h 的频谱。
把一级近似值 (1)
h
代入式(7-9-13)得 (1)(1)102()()!n n Hs n s Z Ms e h n μ∞-=⎡⎤-∆=-⎢⎥⎣⎦
∑ (7-9-17) 把 (1)Z ∆按下式求(2)h : (2)(1)02[]Hs e h
Z Z Ms μ=-∆-∆ (7-9-18) 把 (2)h
反变换求得(2)h ,则(1)(2)h h +为h 的二次近似值,把它带入式(7-9-13)求得 Z ∆的二级近似值 (2)Z ∆。
按下式求 (3)
h
: (3)(2)02(,)[(,)(,)]Hs e h u v Z u v Z u v Ms μ=-∆-∆ (7-9-19) 则(1)(2)(3)h h h ++为h 的三次近似值,再继续求 (3)Z
∆……如此下去,直到h 的第n 次近似值与第(n-1)次近似值之差满足所要求的精度为止。
这时得到的界面起伏为
()1n i i h h ==∑(7-9-20)。