首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 156 毫秒
1.
各向异性普遍存在于强磁性体矿物中,为研究各向异性强磁性体磁场响应特征,本文提出一种空间波数混合域三维各向异性磁场数值模拟方法.该方法首先将各向异性强磁性体磁位满足的三维偏微分方程进行水平方向二维傅里叶变换,将其降为不同波数之间相互独立的一维常微分方程;然后加载准确的上下边界条件,采用二次插值有限单元法计算一维常微分方程,得到五对角方程,采用追赶法进行高效求解;最后采用迭代法求解场分量,引入紧算子保证迭代稳定收敛;综合傅里叶变换的高效性、一维方程求解的快速性和迭代算法的稳定性,实现各向异性强磁性体磁场的三维高效、高精度数值模拟.设计各向异性椭球模型验证算法的正确性,并分析紧算子对不同各向异性磁化率模型的迭代收敛性;与COMSOL Multiphysics软件对比计算效率,表明相同节点下本文算法效率优于常规三维有限元方法,且计算节点总数越多优势越明显.重点研究各向异性参数改变对VTI、HTI、TTI强磁介质异常场响应的影响.最后采用某磁铁矿DEM高程数据模拟起伏地形对各向异性强磁性体磁异常场幅值和形态的影响,体现出本文算法对各向异性强磁性体大规模复杂地形的适应性.  相似文献   

2.
强剩磁强退磁条件下的二维井中磁测反演   总被引:5,自引:4,他引:1       下载免费PDF全文
强剩磁、强退磁改变了总磁化强度的大小和方向,给磁测资料解释带来困难.为此,本文利用二维井中磁测数据反演磁化强度矢量的二维分布.首先利用井中磁测的磁异常模量反演磁化强度大小的分布.然后,在已知磁化强度大小分布的前提下,拟合磁场分量,反演磁化强度方向的分布.其中,磁化强度大小和方向均用共轭梯度法求解,并通过预优矩阵改善磁化强度大小的反演效果.理论模拟说明,该方法能准确获得磁化强度矢量分布.磁化强度矢量反演结果包括感磁、剩磁及退磁的影响,这为研究强剩磁、高磁化率矿床提供了一种有效方法.  相似文献   

3.
忽视自退磁效应的影响,会导致磁异常反演结果的误差,影响地质解释的效果,这使得人们在实际应用中不得不考虑退磁的影响.本文对考虑自退磁影响的三维复杂形体磁场正演和反演方法进行了系统总结,论述了各种方法的特点和实际应用中的困难,提出了考虑自退磁影响的三维复杂形体磁场正反演研究的发展方向.  相似文献   

4.
基于电场Helmholtz方程的回线源瞬变电磁法三维正演   总被引:5,自引:5,他引:0       下载免费PDF全文
正演是电磁法勘探野外工作参数选取、室内资料处理与解释的基础,精确、稳定、高效的三维正演算法尤为重要.本文采取先求解拉普拉斯域电场、再由Gaver-Stehfest算法获得时间域磁场的思路,基于电场异常场Helmholtz方程实现了交错网格有限差分法和有限体积法对回线源瞬变电磁法的三维正演.通过对比低阻块状体的积分方程法、时域有限差分法、矢量有限单元法和SLDM法的数值解,验证了交错网格有限差分法和有限体积法的正确性.由于交错网格有限差分法、有限体积法和基于矩形块单元的矢量有限单元法将待求电场均定义在矩形块单元棱边上,因此三种数值算法可采用相同方法进行电场待求量编码、计算背景场和后处理.然而,与矢量有限单元法相比,交错网格有限差分法和有限体积法的系数矩阵更加稀疏,求解效率更高.通过对水平低阻板状体三维模型的数值模拟,我们发现本研究中交错网格有限差分法比有限体积法精度更高;再利用一维解析法求解相应三层层状地电模型的感应电动势,我们还发现两种数值算法和一维解析法计算的感应电动势等值线形状吻合程度高,只是数值范围略有差异.  相似文献   

5.
重磁异常相关成像法是一种快速确定场源位置的有效方法,其相关系数极大值表征地质体中心位置.现有方法需已知地质体构造指数,本文提出重磁不同阶梯度比值的相关成像方法,可有效去除构造指数的影响,并讨论不同组合梯度比值的应用效果.对于磁异常数据,考虑到剩磁的影响,将采用解析信号及其梯度的比值来获取地质体的分布.通过理论模型试验,证明梯度比值相关成像法可以确定地质体中心位置,也具备良好的抗噪性.此外,解析信号二阶垂直梯度与解析信号比值的相关成像结果最稳定,精度和分辨率最高,为了降低噪声的干扰,在二阶及以上导数计算时采用Laplace方程来完成.将本文方法应用于埃及Hamrawien地区的实测磁数据的解释,反演获得地下异常体的深度在680 m和808 m.  相似文献   

6.
本文讨论稳定磁化场中强磁性三度体内部磁化强度的数值解法。展示了由于退磁作用而造成的不均匀性,并以长方体内一些典型截面为例,较详细地讨论了磁化强度的特征,作了分布图;推导了较均匀磁化体磁化强度分布的近似公式,并计算了一组长方体的视退磁系数。  相似文献   

7.
球形磁性地质体是地质勘探中广泛遇到的基本地质体之一,选择合适的磁异常分析方法对其磁场进行准确正演计算具有重要的意义.为了准确快速地求解球形磁性地质体磁场特征,本文提出了一种基于磁偶极子构造原理求解球形磁性地质体磁场的数值计算方法,对球形磁性地质体按半径由小增大进行多层单元分割,再对每层单元进行块单元分割,将每个块单元视为"磁偶极子",利用基本磁偶极子公式计算了所有块单元在探测平面(线)的磁场大小,同时分析了球体磁场的剖面特征,并将数值解与解析解结果进行了对比和误差分析,最后建立双球组合模型结合Comsol多物理场仿真工具对其仿真对比,验证了磁偶极子构造法求解球体磁场的准确性.研究结果表明:数值解法结果与解析解结果磁异常波动趋势完全一致;不同测点处绝对误差有所差别,但磁异常值均在1 nT以下.本文提出的数值解法无需复杂的数学推导,计算结果稳定可靠.  相似文献   

8.
为了提高海洋可控源电磁的地形响应模拟精度,探讨了三维交错网格剖分有限体积法求解频率域电磁场解的数值方法.在似稳场近似条件下,推导了有限差分法和有限体积法低频电磁场控制方程的离散表达式.通过对离散表达式的分析比较表明,有限体积法较之有限差分具有更高的计算精度,同时也具有与有限差分相当的计算效率.对含地形界面网格单元电导率采用加权平均处理,通过与二维有限元程序的计算结果对比,验证了该方法可以有效提高有限体积法对地形变化的模拟精度.应用有限体积法计算了一个三维带地形储层模型的可控源电磁响应,分析表明地形变化对电场分量影响明显,磁场分量对地形变化不敏感.  相似文献   

9.
含菱铁矿砂岩磁性的初步研究   总被引:12,自引:0,他引:12       下载免费PDF全文
作为煤田盖层产出的中国内蒙古平庄J3砂岩,其平均磁化率为60×10-5(SI);平均剩磁为5×10-3A/m.经热退磁后,在磁化率、剩余磁化、等温剩磁曲线、粘滞剩磁及磁组构等方面都发生了明显的变化,出现一系列奇异磁特征.岩石透射光与反射光鉴定、物相热红外分析、差热分析、X射线衍射分析表明,该砂岩的胶结物为弱磁性的菱铁矿,它在逐步热退磁过程中逐渐变为磁铁矿,部分变为赤铁矿.并大致确定,该砂岩胶结物含量近50%.此外,对这种砂岩样品仔细地进行了常规方法的退磁试验,得到其特征剩磁方向,确定其J3古地磁极位置为198.8°E,84.6°N(α95=4.8,k=94.5).  相似文献   

10.
在地质体具有较强剩磁或自退磁效应的情况下,获得地质体总磁化方向对磁法勘探的数据处理、反演及解释具有重要意义.与传统航磁总场测量相比,航磁矢量(三分量)测量能够获得地磁场的矢量信息.基于三分量磁异常一阶矩与磁源磁矩的积分关系,我们实现了基于航磁矢量数据的磁源总磁化方向估算.针对该方法易受相邻异常影响,从而使多异常的估算难以确定磁源中心位置的问题,分析了在不同磁化倾角下,估算的磁化偏角、倾角与实际磁源中心位置的关系,并根据该关系采取了限定估算区域的方法,获得了多异常的估算结果.将该方法应用于东天山启鑫地区实测航磁三分量数据,估算了具有强剩磁的启鑫岩体总磁化方向,并将估算的方向应用于化极计算.化极结果显示,减少了启鑫岩体磁异常受斜磁化的影响,且与磁异常模量计算结果类似.该方法假设磁源均匀磁化,估算结果可视为地质体总磁化方向的整体体现,对强剩磁地质体磁异常的数据处理、反演及解释具有实际意义.  相似文献   

11.
Summary Conditions are given under which two thick plates, differing in dip, apparent susceptibility, and remanence, will produce similar magnetic anomalies. From these conditions correction formulae are developed. Using these formulae the dip and susceptibility of a plate with remanent magnetization can be obtained from those of non-remanent plate. An interpretation procedure is suggested where the magnetic anomaly is first interpreted by means of a plate without remanence, dip and apparent susceptibility are then estimated by using the correction formulae developed. Thickness, position and depth of the plate are unaffected by the remanence correction procedure. The procedure is independent of the field component measured.  相似文献   

12.
To further evaluate the potential of magnetic anisotropy techniques for determining the origin of the natural remanent magnetization (NRM) in sedimentary rocks, several new remanence anisotropy measurement techniques were explored. An accurate separation of the remanence anisotropy of magnetite and hematite in the same sedimentary rock sample was the goal.In one technique, Tertiary red and grey sedimentary rock samples from the Orera section (Spain) were exposed to 13 T fields in 9 different orientations. In each orientation, alternating field (af) demagnetization was used to separate the magnetite and hematite contributions of the high field isothermal remanent magnetization (IRM). Tensor subtraction was used to calculate the magnetite and hematite anisotropy tensors. Geologically interpretable fabrics did not result, probably because of the presence of goethite which contributes to the IRM. In the second technique, also applied to samples from Orera, an anisotropy of anhysteretic remanence (AAR) was applied in af fields up to 240 mT to directly measure the fabric of the magnetite in the sample. IRMs applied in 2 T fields followed by 240 mT af demagnetization, and thermal demagnetization at 90°C to remove the goethite contribution, were used to independently measure the hematite fabric in the same samples. This approach gave geologically interpretable results with minimum principal axes perpendicular to bedding, suggesting that the hematite and magnetite grains in the Orera samples both carry a depositional remanent magnetization (DRM). In a third experiment, IRMs applied in 13 T fields were used to measure the magnetic fabric of samples from the Dome de Barrot area (France). These samples had been demonstrated to have hematite as their only magnetic mineral. The fabrics that resulted were geologically interpretable, showing a strong NW-SE horizontal lineation consistent with AMS fabrics measured in the same samples. These fabrics suggest that the rock's remanence may have been affected by strain and could have originated as a DRM or a CRM.Our work shows that it is important to account for the presence of goethite when using high field IRMs to measure the remanence anisotropy of hematite-bearing sedimentary rocks. It also shows that very high magnetic fields (>10 T) may be used to measure the magnetic fabric of sedimentary rocks with highly coercive magnetic minerals without complete demagnetization between each position, provided that the field magnetically saturates the rock.  相似文献   

13.
A comprehensive rock magnetic, magnetic anisotropy and paleomagnetic study has been undertaken in the brecciated LL6 Bensour ordinary chondrite, a few months only after its fall on Earth. Microscopic observations and electronic microprobe analyses indicate the presence of Ni-rich taenite, tetrataenite and rare Co-rich kamacite. Tetrataenite is the main carrier of remanence. Magnetization and anisotropy measurements were performed on mutually oriented 125 mm3 sub-samples. A very strong coherent susceptibility and remanence anisotropy is evidenced and interpreted as due to the large impact responsible for the post-metamorphic compaction of this brecciated material and disruption of the parent body. We show that the acquisition of remanent magnetization postdates metamorphism on the parent body and predates the entering of the meteorite in Earth’s atmosphere. Three components of magnetization could be isolated. A soft coherent component is closely related to the anisotropy of the meteorite and is interpreted as a shock remanent magnetization acquired during the same large impact on the parent body. Two harder components show random directions at a few mm scale. This randomness is attributed either to the formation mechanism of tetrataenite or to post-metamorphic brecciation. All components are likely acquired in very low (≈μT) to null ambient magnetic field, as demonstrated by comparison with demagnetization behavior of isothermal remanent magnetization. Two other LL6 meteorites, Kilabo and St-Mesmin, have also been studied for comparison with Bensour.  相似文献   

14.
Magnetic carriers in remagnetized Cretaceous granitic rocks of northeast Japan were studied using paleomagnetism, rock magnetism, optical microscopy and scanning electron microscopy (SEM) by comparison with unremagnetized granitic rocks. The natural remanent magnetization (NRM) of the remagnetized rocks is strong (0.3–1.7 A/m) and shows a northwesterly direction with moderate inclination (NW remanence), whereas the unremagnetized rocks preserve weak NRM (<0.5 A/m) with westerly and shallow direction (W remanence). Although thermal demagnetization shows that both NRMs are carried by magnetite, the remagnetized rocks reveal a higher coercivity with respect to alternating field demagnetization (20 mT相似文献   

15.
本文报道利用岩石剩磁组构对华北下三叠统红层进行磁倾角浅化效应的进一步识别与校正研究结果.首先,采用45°等温剩磁各向异性方法,即通过沿与样品原始水平面(即层面)呈45°夹角方向施加磁化场获得等温剩磁,并进行逐步热退磁,获得平行于层面和垂直于层面的等温剩磁分量随外加磁场和热退磁温度的变化趋势,计算获得浅化因子f=0.70.其次,应用高场等温剩磁各向异性方法,结合峰值为100 mT的交变退磁和120℃热退磁处理,分离获得碎屑赤铁矿对剩磁各向异性的贡献;由直接测量获得的单颗粒碎屑赤铁矿的各向异性度(a=1.35),计算获得f=0.59.该结果与前人对刘家沟组红层进行E/I法磁倾角浅化校正的结果(f=0.60)具有很好的一致性;表明华北下三叠统刘家沟组红层磁倾角浅化效应显著,其浅化因子为f=0.59;高场等温剩磁各向异性方法是红层磁倾角浅化校正的最有效方法.同时,如果有足够的独立样品数,且特征剩磁来自于单一剖面或有证据表明多条采样剖面之间未发生显著的相对运动,E/I法对红层磁倾角浅化因子的估计也是可信的.  相似文献   

16.
This paper reports the alternating field demagnetization characteristics of glass–ceramic magnetite assemblages carrying weak-field thermoremanent magnetization (TRM), weak-field anhysteretic remanent magnetization (ARM), and saturation remanence (Jrs). Average grain sizes vary from less than 0.1 μm to approximately 100 μm, and hysteresis parameters indicate that these assemblages encompass single-domain (SD) through truly multidomain (MD) behavior. In all assemblages, weak-field TRM and weak-field ARM are more stable to alternating field demagnetization than is (Jrs). This response is especially remarkable in the 100 μm assemblage, which otherwise displays truly MD behavior. Although the SD samples pass the Lowrie–Fuller test for SD behavior, calculations presented here show that populations of noninteracting, uniaxial SD grains should behave in just the opposite sense to that reported originally by Lowrie and Fuller. This discrepancy could indicate that SD, glass–ceramic magnetite populations are more affected by magnetic interactions than would be expected for magnetite crystals that nucleated individually from a silicate matrix. This interpretation is supported by the SD assemblages failing the ‘Cisowski' test: that is, the curves for acquisition and AF demagnetization of (Jrs) intersect well below the 50% mark. However, a second and intriguing explanation of the SD-like results obtained from all samples is that alternating field demagnetization characteristics reflect a strong dependence of local energy minimum domain state, and its associated stability, on the state of magnetization.  相似文献   

17.
To identify the stable remanence carrier in rock samples, we conducted magnetic microscopic observations combined with conventional stepwise demagnetization experiments. The instrument, which employs an amorphous wire-based magneto-impedance sensor (30 μm diameter, 5 mm length), can document magnetic anomalies (vertical component) of the millimeter to sub-millimeter-thick rock samples with a resolution of 500 μm. Our new technique allows identification of the sources of both stable and unstable remanence components in meteorite and shocked granite samples. However, stray magnetic fields from the sensor magnetize the magnetic minerals in the sample and makes serious artifacts on the magnetic images. Although the artifacts of the induced magnetization should be solved, this new corroborative technique leads to a microscopic discrimination of stable paleomagnetic records from terrestrial and extraterrestrial materials.  相似文献   

18.
Mars Global Surveyor (MGS) observations of crustal magnetic fields over Tharsis provide new constraints on models for the thermal and magmatic evolution of this region. We analyze the distribution of magnetic field anomalies over Tharsis surface units of Noachian, Hesperian and Amazonian age. These data suggest that early Noachian crust underlies the Tharsis province, and formed contemporaneously with the existence of a martian dynamo. This crust either pre-dates the formation of Tharsis, or formed during the earlier phases of Tharsis volcanism. The preservation of strong magnetic field anomalies over some of the earliest Noachian and topographically high units, together with the observation of magnetic field anomalies over Hesperian- and Amazonian-age surface units, indicate that a large fraction of the magnetized crust has remained cool (below the blocking temperature of the magnetic carrier) throughout the construction of Tharsis. Moreover, the distributions of magnetic anomaly amplitudes over Noachian, Hesperian, and Amazonian surface units suggest that the youngest units overlie sites of prolonged intrusion and have undergone a greater extent of thermal demagnetization. The absence of magnetic anomalies around the Tharsis Montes and Olympus Mons argues for strong, localized heating, as would be expected at volcanic centers. We show that end-member models for progressive thermal demagnetization of a Noachian magnetized crustal layer are consistent with the anomaly amplitude distributions. We integrate the magnetic field observations with constraints from tectonics, gravity, and topography, and present a revised scenario for the evolution of the Tharsis region.  相似文献   

19.
The magnetic and gravity field produced by a given homogeneous source are related through Poisson's equation. Starting from this consideration, it is shown that some 2D interpretation tools, widely applied in the analysis of aeromagnetic data, can also be used for the interpretation of gravity gradiometric data (vertical gradient). This paper deals specifically with the Werner deconvolution, analytic signal and Euler's equation methods. After a short outline of the mathematical development, synthesized examples have been used to discuss the efficiency and limits of these interpretation methods. These tools could be applied directly to airborne gravity gradiometric data as well as ground gravity surveys after transformation of the Bouguer anomalies into vertical gradient anomalies. An example is given of the application of the Werner deconvolution and Euler's equation methods to a microgravity survey.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号