共查询到20条相似文献，搜索用时 15 毫秒
1.
TAPIO RUOTOISTENMÄKI 《Geophysical Prospecting》1993,41(4):413-433
Magnetic anomalies of complicated 3D sources can be calculated by using a combination of analytical and numerical integration. Two surfaces and the magnetization parameters (the amplitudes of the induced and remanent components and the direction cosines) of the source can be defined by arbitrary functions or by discrete data points in a plane. When combined with a polynomial magnetization function in the direction of the third axis, 3D magnetization distribution can also be modelled. The method gives very general equations for anomaly calculation. It can be used for direct modelling of sources interpreted by seismic or other methods and also for interactive interpretation with fast computers. It is possible to calculate anomalies of, for example, intrusives or folded sedimentary beds whose surfaces are functions of horizontal coordinates and which have polynomial magnetization variations in the vertical direction due to gravitational differentiation and arbitrarily varying magnetization in the horizontal direction due to regional metamorphosis. If the distribution of magnetization parameters in the vertical direction cannot be described satisfactorily by polynomials, models can be used whose surfaces are functions of the vertical coordinate and which can then have any arbitrary magnetization distribution in the vertical direction. 相似文献
2.
Peter Furness 《Geophysical Prospecting》2000,48(3):603-615
Existing techniques for computing the gravitational field due to a homogeneous polyhedron all transform the required volume integral, expressing the field due to a volume distribution of mass, into a surface integral, expressing the potential due to a surface mass distribution over the boundary of the source body. An alternative representation is also possible and results in a surface integral expressing the potential due to a variable-strength double layer located on the polyhedral source boundary. Manipulation of this integral ultimately allows the gravitational field component in an arbitrary direction to be expressed as a weighted sum of the potentials due to two basic source distributions. These are a uniform-strength double layer located on all faces and a uniform-strength line source located along all edges. The derivatives of the gravitational field components can also be expressed in a similar form as can the magnetic field components due to a homogeneous magnetic polyhedron. It follows that the present approach can be used to generate a universal program capable of modelling all the commonly used potential field responses due to 3D bodies of arbitrary shape. 相似文献
3.
在经典位场理论中,许多简单形体位场异常难以通过积分得到全空间的解析式.圆柱体是一类很重要的理论模型体,常用于模拟圆柱状地质体或非地质体(如管线),但目前还不能用解析公式正演有限长圆柱体在三维空间里的磁异常,而多是采用近似简化为有限长磁偶极子或线模型代替.对于有限长圆柱体,特别是半径相对于上顶埋深较大时,这种近似的误差不可忽略.本文利用共轭复数变量替换法,推导出有限长圆柱体在全空间的引力位一阶、二阶导数,利用Poisson关系得到磁异常正演公式,进而利用有限长圆柱体磁异常正演公式求解管状体的磁异常,得到不同磁化方向、不同大小的管线产生的磁场的特征,并将其推广到截面为椭圆的情况.最后通过模拟计算定量给出了将圆柱体近似为线模型的条件. 相似文献
4.
卫星磁异常的理论模型 总被引：3，自引：3，他引：0
本文介绍了计算卫星磁异常理论模型的数学方法 ,即球谐分析方法、冠谐分析方法、矩谐分析方法和等效源方法 .根据相同的 MAGSAT资料 ,计算的卫星磁异常冠谐模型、矩谐模型和等效源模型都能很好地表示卫星磁异常的分布 .由于在整个研究区域 ( 1 0°N～60°N,70°E～ 1 40°E)都有卫星资料 ,所以这些理论模型没有所谓的“边界效应”.这一结论对计算地磁场的区域模型是很有意义的 . 相似文献
5.
We derive expressions for computing the gravitational field (potential and its radial derivative) generated by an arbitrary
homogeneous or laterally varying density contrast layer with a variable depth and thickness based on methods for a spherical
harmonic analysis and synthesis of gravity field. The newly derived expressions are utilised in the gravimetric forward modelling
of major known density structures within the Earth’s crust (excluding the ocean density contrast) beneath the geoid surface.
The gravitational field quantities due to the sediments and crust components density contrasts, shown in numerical examples,
are computed using the 2 × 2 arc-deg discrete data from the global crustal model CRUST2.0. These density contrasts are defined
relative to the adopted value of the reference crustal density of 2670 kgm−3. All computations are realised globally on a 1 × 1 arc-deg geographical grid at the Earth’s surface. The maxima of the gravitational
signal due to the sediments density contrast are mainly along continental shelf regions with the largest sedimentary deposits.
The corresponding maxima due to the consolidated crust components density contrast are over areas of the largest continental
crustal thickness with variable geological structure. 相似文献
6.
通过对欧洲及其邻近地区 MAGSAT 卫星黎明资料进行了处理,得到1°×1°卫星磁异常网格值。本文使用冠谐分析方法,计算该地区卫星矢量磁异常（ΔX,ΔY,ΔZ）冠谐模型。球冠极点位于33°N和26°E,球冠半角为40°.冠谐模型的截断指数为18。根据卫星磁异常冠谐模型和地磁场球谐模型 DGRF1980,计算卫星总强度磁异常（ΔF）的冠谐一球谐模型。根据卫星磁异常的理论模型,计算并绘制不同高度（300,400,500km）的理论卫星磁异常图。对冠谐模型和大磁异常进行了分析和讨论。 相似文献
7.
Summary The solution of the direct magnetic problem is derived for three-dimensional bodies under the assumption that the components of the magnetization vector are analytical functions of the co-ordinates within the body being considered. The solution to the problem is given in the form of Green-type integrals and, much like with gravitational fields of inhomogeneous bodies, it may serve the purpose of solving the problem of the analytical continuation of the external field into the body. 相似文献
8.
M. Ivan 《Geophysical Prospecting》1994,42(7):735-743
Information on the mass and the spatial location of an arbitrary source body can be obtained by performing suitable integrations of 3D gravity and magnetic data along an infinite straight line. No assumptions on the density/magnetization distribution or the shape and location of the source are required. For an oblique borehole, a relationship between the lower limit of the source mass and the distance to the body is obtained. The mass contrast and the magnetic moment of the source can also be estimated. For a vertical borehole, both gravity and vertical magnetic component anomalies have equal areas to the left and right of the depth axis. The particular case of a horizontal gallery not intersecting the body is also studied. If the source is intersected, a lower limit is estimated for the maximum thickness of the body along the gallery. Information on the vertical coordinate of the centre of mass of the source can also be obtained. Numerical tests with synthetic gravity data support the theoretical results. 相似文献
9.
Thc expressions of magnetic gradients due to 3-D homogeneous magnetized polyhedra are systematically derived and presented,
fronr which the forward problem of magnetic gradients of an arbitrary shaped geological body is solved. It is shown that in
the rotation of coordinate systems there is an essential difference between the transformation of magnetic fields and that
of their gradients. In a 2-D coordinate system a unified transformation formula of any order gradients can he derived, hut
cannot in the 3-D case. The calculations of synthetic models show the correctness of the expressioils or magnetic gradients.
Project supported by the National Natural Science Foundation of China. 相似文献
10.
A versatile integral equation technique for magnetic modelling 总被引：1，自引：0，他引：1
Peter Furness 《Journal of Applied Geophysics》1999,41(4):737
A requirement currently exists in both mineral exploration and environmental or engineering geophysics for a technique to model the magnetic fields caused by bodies with large to extreme susceptibilities in which both induced and remanent magnetizations are significant. It is well known that modelling such magnetic fields is not amenable to any known approximation. It is a significantly difficult task that requires the solution of a magnetostatic boundary value problem. Analytical solutions to the problem are extremely useful for providing insight but generally of limited application in practical interpretation due to the geometrical complexity of real situations. Available numerical solutions include both volume and surface integral equation formulations. However neither of these are particularly efficient for the purpose. An alternative surface integral equation formulation is presented here which represents the required magnetic field in terms of a double layer over the surface of the body. The technique accommodates both remanent and induced magnetization and is generally applicable to any 3D body in a magnetic environment for which the Green's function is available. The present technique has significant advantages over other integral equation solutions in the geophysical literature. It is particularly economic in terms of the density of the surface discretization and consequently the computational effort. Moreover, it is extremely robust. It is found to yield accurate solutions for the type of thin bodies that cause numerical instability with other surface integral equation approaches. 相似文献
11.
12.
本文将均质的任意二维、三维物体位场的波谱解析表达式的研究成果推广到变密度、变磁化强度的更一般的情形。对密度差随深度呈指数函数衰减或线性变化的模型,获得了任意倾斜多边形质量面、斜平行六面体以及一般的多面体等形体的重力谱的解析表达式。它们的结构与均质体相应表达式一样简单,易于计算。以上结果表明,在很一般的条件下,位场波谱具有指数函数和的形式。 相似文献
13.
14.
峨眉山大火成岩省位于中国西南部,在晚二叠纪约260 Ma喷发出巨量的大陆溢流型玄武岩.对于大火成岩省的岩浆喷发,在地下必定有一个相应的大规模岩浆聚集和运移系统.地球物理方法是探测岩石圈内部的有效方式.峨眉山大火成岩省为镁铁质岩浆喷发,由于镁铁质-超镁铁质岩石一般具有强磁性,因此,在喷发结束之后,地下岩浆系统如果被镁铁质岩浆填充,冷却固化成为岩石圈的一部分,很有可能会引起磁异常.本文使用区域磁异常数据来对峨眉山大火成岩省的深部构造进行研究.该区域的磁异常由一系列离散的异常组成,通过3D磁化率反演可以得到磁性体的空间分布.由于磁异常中具有明显的剩磁,直接使用经典的反演方法会有较大误差,我们首先将磁异常转换为弱敏感于磁化方向的磁异常模量,再使用模量数据进行3D反演,得到地下空间内磁异常源的分布.经过分析认为这些离散分布的磁异常源反映了岩石圈内部的镁铁质-超镁铁质侵入体.侵入体的位置可能反映了底侵和内侵的镁铁质岩浆固化形成的侵入体,代表镁铁质岩浆房位置或者岩浆运移的主要通道. 相似文献
15.
THE INTERPRETATION OF GRAVITATIONAL AND MAGNETIC DATA FOR THE VOLCANIC ROCKS DISTRIBUTION IN KE-BAI FAULT ZONE,JUNGGAR BASIN 
下载免费PDF全文
下载免费PDF全文 This paper firstly discusses the feasibility of delineating the volcanic rocks distribution by gravitational and magnetic methods on the basis of the statistical results of the magnetic susceptibility and density of the different lithologies in the study area. After the separation of gravitational and magnetic fields by bandpass filter, we determined the residual gravitational and magnetic anomalies caused by volcanic rocks. The results of potential field separation show that the residual anomalies are in beaded NE-directed distribution. In the meantime, the boundary enhancement techniques such as horizontal total gradient and vertical derivative are employed to delineate the distribution of the faults, as the result, fifteen faults including four NE-striking main faults which appear as an arc protruding to the southeast and other secondary faults are distinguished. Furthermore, two fault systems with deep and superficial attribution are revealed from the inversed section of the telluric electromagnetic sounding, and their properties, characteristics and roles in the volcanic activity are fully discussed by combining with the tectonic background. Based on the comprehensive analysis of the correlation between the distribution characteristics of the residual gravitational and magnetic anomalies and the location of the faults, three volcanic activity zones are reasonably delineated. Finally, we carry out the inversion of apparent density by taking advantage of the residual gravitational anomalies and acquire the apparent density anomalies of the top part of Paleozoic. Integrated with the apparent density anomalies, the lithologies exposed from the boreholes, the observed density data of different rocks and the residual gravitational-magnetic anomalies, we elaborately delineate the lithologies of the top of Paleozoic and further define the distribution range and the lithology combination of Kekesayi group of Middle Ordovician, Tailegula group of Lower Carboniferous and Jiamuhe group of Lower Permian. 相似文献
16.
Azimuthally averaged power spectra are widely used in the Curie point depth (CPD) estimation with the implicit assumption that the magnetization distribution is random and uncorrelated. However, the marine magnetic anomalies are caused by bands of normal and reverse magnetization and show obvious trends. To investigate the effects of the anisotropy of marine magnetic anomalies on the CPD estimates, we develop 3D fractal striped magnetization models to produce lineated marine magnetic anomalies for the first time. We analyze the spectra anisotropy of the lineated magnetic anomalies of the synthetic fractal striped magnetization models and investigate its effects on the CPD estimates. The synthetic models and actual data show that the spectra of the lineated marine magnetic anomalies are directionally anisotropic. The amplitude response is strong and the slope of the logarithmic spectrum is large in a direction perpendicular to the stripes of magnetic anomalies, whereas the amplitude response is weak and the slope of the logarithmic spectrum is small in a direction parallel to the stripes of magnetic anomalies. The depth estimates in the perpendicular direction are close to the actual values, whereas the depths estimates in the parallel direction are significantly lower than the actual values. The actual marine magnetic anomalies of the South China Sea exhibit an anisotropic power spectrum that is consistent with the spectral anisotropy of magnetic anomalies of the synthetic fractal striped magnetization models. 相似文献
17.
Detection of potential fields source boundaries by enhanced horizontal derivative method 总被引：1，自引：0，他引：1
A high‐resolution method to image the horizontal boundaries of gravity and magnetic sources is presented (the enhanced horizontal derivative (EHD) method). The EHD is formed by taking the horizontal derivative of a sum of vertical derivatives of increasing order. The location of EHD maxima is used to outline the source boundaries. While for gravity anomalies the method can be applied immediately, magnetic anomalies should be previously reduced to the pole. We found that working on reduced‐to‐the‐pole magnetic anomalies leads to better results than those obtainable by working on magnetic anomalies in dipolar form, even when the magnetization direction parameters are not well estimated. This is confirmed also for other popular methods used to estimate the horizontal location of potential fields source boundaries. The EHD method is highly flexible, and different conditions of signal‐to‐noise ratios and depths‐to‐source can be treated by an appropriate selection of the terms of the summation. A strategy to perform high‐order vertical derivatives is also suggested. This involves both frequency‐ and space‐domain transformations and gives more stable results than the usual Fourier method. The high resolution of the EHD method is demonstrated on a number of synthetic gravity and magnetic fields due to isolated as well as to interfering deep‐seated prismatic sources. The resolving power of this method was tested also by comparing the results with those obtained by another high‐resolution method based on the analytic signal. The success of the EHD method in the definition of the source boundary is due to the fact that it conveys efficiently all the different boundary information contained in any single term of the sum. Application to a magnetic data set of a volcanic area in southern Italy helped to define the probable boundaries of a calderic collapse, marked by a number of magmatic intrusions. Previous interpretations of gravity and magnetic fields suggested a subcircular shape for this caldera, the boundaries of which are imaged with better detail using the EHD method. 相似文献
18.
This paper tries to formulate the C-response of geomagnetic depth sounding(GDS)on an Earth model with finite electrical conductivity. The computation is performed in a spherical coordinate system. The Earth is divided into a series of thin spherical shells. The source is approximated by a single spherical harmonic P10 due to the spatial structure of electrical currents in the magnetosphere. The whole solution space is separated into inner and external parts by the Earth surface. Omitting displacement current, the magnetic field in the external space obeys Laplacian equation, while in the inner part, due to the finite conductivity, the electromagnetic fields obey Helmholtz equation. To connect the magnetic fields in the inner and external space, the continuity condition of magnetic fields is used on the Earth surface. The external magnetic fields are expressed by the inner and external source coefficients, from which a new parameter called C-response is computed from the inner coefficient divided by the external coefficient, thus normalizing the actual source strength. The inner magnetic fields in each layer can be recursively derived by the continuity boundary condition of both normal and tangential components of the magnetic field from the initial boundary condition at core-mantle-boundary. The consistency of our C-responses with that from a typical 1-D global model validates the accuracy of the proposed algorithm. Numerical results also show that the C-response estimated from the geomagnetic transfer function method will deviate exceeding 5%from the actual response at longer periods than about 106s, which means that ignoring the curvature of the Earth at extreme long periods will make inversion result unreliable. Therefore, an accurate C-response should be computed in order to lay a solid foundation for reliable inversion. 相似文献
19.
导了磁化强度矢量层析成像方程,并与磁化强度标量层析成像方程进行了对比. 使该矢量层析成像方程既适用于三维的也适用于二维,既适用于使用磁场垂直分量资料,也适用于使用磁场总强度资料. 本文采用改进的高斯-赛德尔迭代求解磁化强度矢量层析成像方程,在求解方程中引进了与深度有关的权系数. 并对二维模型开展了正反演研究. 当模型层数为二层,且每层51个柱体时,采用零初始模型就能获得较好的磁化强度垂直分量和水平分量反演结果,除了异常体边部外,磁化方向比较可靠. 当模型层数为五层,且每层51个柱体时,采用零初始模型不能得到较好结果,改用接近背景值的均匀初始模型,反演的磁化强度垂直分量和水平分量大致接近真实,但不能刻划某些细节. 相似文献