MVO曲线方法在海洋可控源电磁数值模拟中的应用研究

采用有限差分方法,对水平电偶极子源在海洋二维介质情况下电磁响应进行数值模拟。利用油气高阻异常体模型的数值模拟结果,研究MVO(Magnitude Versus Offset)曲线解释方法的有效性,进而分析高阻油气异常体的水平位置的变化,对电场和磁场幅值的MVO曲线产生的影响,得出了通过MVO曲线定性分析高阻油气异常体相对水平位置的方法。本研究可为海洋可控源电磁法数据的定性分析提供部份理论依据。     

基于重磁场特征的敦化盆地构造格架研究

依据重力、磁力异常数据及其处理结果( 水平梯度模和斜导数) ,对敦化盆地边界、基底起伏、断裂位置及以火成岩为代表的磁性体分布进行了研究。盆地重力异常的分析和水平梯度模及斜导数的计算结果表明,盆地基底具南部凹陷、中央凸起和北部凹陷的“两凹一凸”的起伏形态特征,盆地内断裂以SW--NE 向为主,盆地为单断半地堑式盆地。依据航磁异常,将盆地划分为4 个异常区: 东北部磁异常区、中部低磁异常区、西南高磁异常区和西南边部相对低磁异常区。结合磁异常水平梯度模和斜导数的计算结果显示,以火成岩为代表的磁性体受SW--NE 向构造控制。     

估计磁场源位置的倾斜角梯度算法改进

倾斜角梯度法在未知场源物性先验知识和构造指数情况下,可快速估计出场源的边界位置和深度分布,对于大面积实测磁异常的快速解释具有应用前景。笔者在Salem的倾斜角梯度法基础上,提出了一些算法改进。在导数换算上采用空间域有限差分法换算磁异常的水平导数,采用ISVD算法换算磁异常的高阶垂直导数。在剔除坏解上,采用磁异常广义倾斜角梯度的总水平梯度滤波法和场源深度有效性筛选法相结合。理论模型数据试验表明本方法简单快速,效果明显。     

A gradient method for interpreting magnetic anomalies due to horizontal circular cylinders,infinite dykes and vertical steps

A new method of interpreting magnetic anomalies of arbitrarily-magnetised horizontal circular cylinders, dipping dykes and vertical steps is presented. The method makes use of both horizontal and vertical gradients of the magnetic field of the model under consideration, rather than the observed magnetic anomaly. Vertical and horizontal gradients are calculated from the observed anomalies, and plotted one against the other to find out the locus of tip of the resultant gradient vector. This locus is a symmetrical curve for each of the three models mentioned above. The properties of these curves are used to deduce the various parameters of these models and the direction of magnetisation.     

Tomographic Imaging of the India-Asia Plate Collisional Tectonics and Mantle Upwelling Beneath Western Tibet

To better understand the lithosphere mantle collision tectonics between the India plate and Asia plate, we determine three dimensional P wave velocity structure beneath western Tibet using 27,439 arrival times from 2,174 teleseismic events recorded by 182 stations of Hi-CLIMB Project and 16 stations in the north of Hi-CLMB. Our tomographic images show the velocity structure significantly difference beneath northern and southern Qiangtang, which can further prove that the Longmu Co-Shuanghu ophiolitic belt is a significant tectonic boundary fault zone. There are two prominent high velocity anomalies and two prominent low velocity anomalies in our images. One obvious high velocity anomalies subduct beneath the Tibet at the long distance near 34°N, whereas it is broke off by an obvious low velocity anomaly under the IYS. We interpret them as northward subducting Indian lithosphere mantle and the low velocity anomanly under IYS likely reflects mantle material upwelling triggered by tearing of the northward subduction Indian lithosphere. The other prominent high velocity anomaly was imaged at a depth from 50 km to 200 km horizontal and up to the northern Qiangtang with its southern edge extending to about 34°N through Hoh Xil block. We infer it as the southward subducting Asia lithosphere mantle. The other widely low velocity anomaly beneath the Qiangtang block lies in the gap between the frontier of India plate and Asia plate, where is the channel of mantle material upwelling.     

Geomagnetic field anomalies over the Lombok Island region: an attempt to understand the local tectonic changes

During the last years, several investigations on the earthquakes and related tectonic structures along the Java Trench have been conducted. In this study, we focus on the Lombok Island—West Nusa Tenggara (Indonesia), which lies between the centres of the highest intensity of contiguous negative–positive geomagnetic anomalies in this region. The geological and tectonic structures, however, are not known in detail for this island, whereas a better knowledge of these conditions could lead to better estimate the potential earthquake risks and thus a possible improvement of the existing early warning system. We have performed a ground-based geomagnetic survey at 56 stations in the Lombok Island during October–November 2006 for a detailed mapping and interpretation of geomagnetic anomalies related to the geological and tectonic characteristics. The 2D and 3D magnetic maps show a general geomagnetic anomaly pattern in the Lombok Island which consists of repeated contiguous negative–positive anomalies. Two forward models have been proposed for a profile connecting a strongest apparent dipolar structure. The first model assumes a uniformly magnetized sphere as the source of the anomaly, and could be interpreted as a specific local structure composed by a quite large magnetic body. The second model considers several rocks with different susceptibilities and magnetizations, and could be interpreted as a discontinuity in the geological structures. This model agrees with the local geological surface conditions and the known large scale regional tectonic structure. Therefore, it is used to interpret our results in terms of tectonic characteristics, which suggests the potential existence of a new tectonic element (e.g. a local normal fault) in this region.     

测量磁场水平分量H y的电性源广域电磁测深法

基于层状大地表面水平电偶极源电磁场公式系统,提出了比值法获取H y分量全区视电阻率计算公式。计算了层状地电模型多个收发距的电磁场各场分量E x、H y、H z,并用这三个场分量计算了全区视电阻率,以及E x/H y卡尼亚电阻率,并与模型的MT卡尼亚电阻率频谱曲线进行对比。发现H y全区视电阻率在赤道装置和轴向装置的广大区域对地电模型均有较好的响应,可进行单分量广域电磁测深。轴向装置H y全区视电阻率与赤道装置垂直磁场H z全区视电阻率相似,在低频段与层状模型大地电磁测深卡尼亚电阻率有相似的频率响应特征,在较小收发距条件下对大埋深基底就能有较好的响应。赤道装置和轴向装置测量的水平电场E x分量全区视电阻率在低频区均会进入“饱和”区。对于各电阻率定义方法,合适的收发距是较好地反映出地电特征深度变化的前提,多收发距测量有利于揭示深部电性特征。     

Modeling and Interpreting CHAMP Magnetic Anomaly Field over China Continent Using Spherical Cap Harmonic Analysis

Based on the CHAMP Magsat data set, spherical cap harmonic analysis was used to model the magnetic fields over China continent. The data set used in the analysis includes the 15′×15′ gridded values of the CHAMP anomaly fields (latitude φ=25°N to 50°N and longitude λ=78°E to 135°E). The pole of the cap is located at φ=35°N and λ=110°E with half-angle of 30°. The maximum index (Kmax) of the model is 30 and the total number of model coefficients is 961, which corresponds to the minimum wavelength at the earth's surface about 400 km. The root mean square (RMS) deviations between the calculated and observed values are ～ 4 nT for ΔX, ～ 3 nT for ΔY and ～ 3.5 nT for ΔZ, respectively. Results show that positive anomalies are found mainly at the Tarim basin with ～6- 8 nT, the Yangtze platform and North China platform with ～4 nT, and the Songliao basin with ～4-6 nT. In contrast, negative anomaly is mainly located in the Tibet orogenic belt with the amplitude ～ (-6)-(-8) nT. Upward continuation of magnetic anomalies was used to semi-quantitatively separate the magnetic anomalies in different depths of crust. The magnetic anomalies at the earth's surface are from -6 to 10 nT for upper crust, middle crust -27 to 42 nT and lower crust -12 to 18 nT, respectively. The strikes of the magnetic anomalies for the upper crust are consistent with those for the middle crust, but not for the lower crust. The high positive magnetic anomalies mainly result from the old continental nucleus and diastrophic block (e.g. middle Sichuan continental nucleus, middle Tarim basin continental nucleus, Junggar diastrophic block and Qaidam diastrophic block). The amplitudes of the magnetic anomalies of the old continental nucleus and diastrophic block are related to evolution of deep crust. These results improve our understanding of the crustal structure over China continent.     

高精度磁测技术在西成铅锌矿勘探中的应用

近年来,甘肃西成铅锌矿的勘探向周边地区扩展,高精度磁法测量发挥了一定作用。利用化极磁异常垂向一阶导数和归一化总水平导数垂向导数,推断出5处地质体边界和4条断裂。根据该区地质和矿产分布特征,推测矿体主要分布于黑云母石英片岩向大理岩的接触或过渡地带,对应中泥盆统安家岔组厂坝层D2a1条带状大理岩。根据化极磁异常和该区构造,结合化探圈定出与成矿有关的4个靶区,这些靶区与Pb、Zn化探异常对应,可见若干矿化点。选取靶区内一磁测剖面进行反演,反演结果说明矿体向北西倾斜。经钻孔验证,已在不同深度见矿,证明了靶区的找矿潜力。     

Seismic ray direction anomalies and relative residuals of earthquakes recorded at Gauribidanur array

Slowness and azimuthal anomalies provide valuable information about lateral inhomogeneities within the crust and mantle of the earth. Over 300 earthquakes (distance range 14°–36° and azimuth 0°–360°) recorded at Gauribidanur seismic array (GBA) in southern India, were analysed using adaptive processing techniques. Slowness anomalies upto 1·3 sec/deg and azimuthal anomalies upto 8° have been observed in the present analysis. Slowness anomaly patterns for Java trench, Mid-Indian oceanic ridge earthquakes are more consistent as compared to the events originating in the Himalayan and Hindukush regions. A significant feature of the azimuthal anomaly pattern was the distinct absence of any positive anomalies from earthquakes occurring in mid-oceanic ridge. These anomalies have also been analysed as a function of epicentral distance and are mainly attributed to the transition zones occurring between 400–700 km depth ranges in the Indian upper mantle regions. Relative residuals between the stations of GBA have very little dependence on azimuth and distance. An anomalous structure beneath the array in the direction of the Java trench region (azimuth 116–126°) has been postulated on the basis of large systematic slowness vectors observed.     

地磁异常数据处理技术在低纬度地区的应用研究——以柬埔寨王国Mesam金矿为例北大核心CSCD

柬埔寨王国位于低纬度地区,地质勘查程度较低,地面高精度磁测技术能够快速获取磁异常信息,进行找矿评价。在低纬度地区磁化方向主要以水平磁化为主,磁性体产生的磁异常以负磁异常为主。本文通过对柬埔寨王国Mesam金矿磁测数据进行倒相180°、化赤和低纬度化极处理对比,认为三种处理方法得到的结果都能够完整地反映磁场分布特点,在不能忽略剩磁影响的低纬度地区对实测数据进行低纬度化极后再进行解释更加贴近真实且细节更加丰富。通过对低纬度化极后的磁异常进行水平总梯度模和小波多尺度分解计算等位场异常处理,认为磁异常水平总梯度模计算能够较好地反映深部地质体的边界范围,小波多尺度分解计算能够提取特定深度地质体的空间信息,在该区选取的磁异常处理技术是有效的,为进一步扩大该金矿区远景储量提供了新的思路。     

Aeromagnetic study of buried basement structures and lineaments of Sahel region (Eastern Tunisia,North Africa)

This study presents original results regarding the use of aeromagnetic to explore deep subsurface structuring in southern part of Tunisian Sahel petroleum province (Eastern Tunisia, North Africa). Several filters and techniques were applied to the total magnetic intensity (TMI) grid. First, an adequate reduce to the pole (RTP) grid was generated. The RTP map shows 17 positive and negative anomalies associated with short- and long-wavelength amplitude anomalies. Positive anomalies correspond to high magnetic basement structures estimated from seismic lines in the area. Circular anomalies are also distinguished locally and can be explained by Cretaceous magmatic rocks recognized in numerous petroleum wells. Second, specific qualitative and quantitative filters (e.g., residual-regional separation, horizontal tilt angle (TDX), tilt angle (TILT), total horizontal derivative of tilt angle (THDTILT)) were applied to elucidate the form and the extent of buried magnetic anomalies and lineaments. The resulting deep structural map revealed the presence of NW-SE, N-S, and E-W regional magnetic basement structures and lineaments and a regional tectonic node surrounding Henchir Keskes-Agareb-Mahres-Hachichina area. Two magnetic inversion models calculated across the study area highlight west to east crustal thinning trend and permit depth to basement estimation. These results are valuable for future conventional and unconventional petroleum exploration in this underexplored southern part of Sahel plain.     

位场反演计算的场比值法及在宽城地区磁异常解释中的应用

在磁异常的反演计算中, 为了避免由于选错了磁性体的形状而产生很大的误差, 前人已研究了多种所谓自动反演的计算方法。与传统的方法不同, 用这些方法, 可以同时确定磁性体的几何形状和埋藏深度等参数。但是前人研究的方法需要计算磁异常的三阶导数, 因此所受高频噪声干扰十分严重, 以至于影响到计算结果的可靠性。因此很难用到实际资料处理中。 本文提出的场比值法可以同时确定磁性体的形状和埋藏深度。场比值法通过场比值field-ratio来确定磁性体的形状, field-ratio的物理含义与Euler反褶积公式中的形状因子类似。场比值法的优点是在反演计算中只需要计算磁异常的一、二阶导数, 因此比前人的方法受高频随机干扰小, 可以用于实际资料的处理中。模型实验证明了场比值法的正确性, 在河北省宽城地区用场比值法对磁异常反演计算, 展示了方法的实用性。 关键词: 磁异常; 自动反演; 磁性体形状的判别     

不同产状磁性体的磁异常特征

目前井中磁测解释大都采用水平矢量图指向确定矿体的方法,但由于磁场水平分量的定向问题,致使仅使用水平分量矢量对磁性体解释常会出现错误结果。通过描述水平产状与倾斜产状磁性体的磁场特征,对磁性体表面不同位置磁异常垂直分量及水平分量差异进行了说明,通过三个实例对反映在磁性体上、下位置的垂直分量与水平分量磁异常的空间分布特征进行了具体描述,并利用井中磁测曲线有效地解决了井底及井旁磁性体的赋存情况。多年来的工作实践证明,使用磁异常场垂直分量及水平分量特征,能准确的判定磁性矿体的主体方向及产状,该解释方法的推广有利于寻找深部磁性矿体。     

重磁异常的小波模极大值分解及其在北衙矿区中的应用

位场数据异常存在的线性特征往往对应着地下断裂构造、不同岩性地质体的边界接触带或者其他具有一定密度或磁性差异的构造特征。对这些线性特征进行增强、提取并进行半定量地解释是重磁资料处理的主要内容。然而,位场数据中若混入噪声干扰,再利用总水平导数、总水平梯度倾斜角等方法进行高次求导运算会把噪声放大,导致提取的线性构造位置发生偏离甚至出现错误。为了获取较为准确的线性构造及边界位置,选取对噪声干扰不敏感的小波模极大值方法,将该方法用于模型试验和实际矿区数据处理中,都能较好地定位出异常体模型边界和矿体的投影边界,表明小波模极大值是一种有效的重磁异常线性特征增强与提取方法。     

利用磁异常模量进行磁性体边界检测

地下磁性体具有较强的剩磁或退磁会给磁异常的解释带来困难。由于磁异常模量对磁化方向不敏感,将磁异常转换为磁异常模量并用来分析和处理,可很大程度上提高磁异常解释的准确性和可靠性。本文介绍了五种磁异常模量,并且利用改进的Tilt梯度方法对各磁异常模量数据做边界检测,确定地下磁性体的水平位置。理论模型试验表明,在存在强剩磁的情况下,不同模量数据对磁性体的反应程度也不一样,运用改进的Tilt梯度方法对磁异常模量进行处理可不同程度地加强对异常体的识别,其中模量E最为准确,可较清晰地识别磁性体边界。笔者还将该方法运用到坦桑尼亚基戈马地区的航磁异常处理,并推测了该地区的构造带及磁性体边界,为深入了解基戈马成矿区带的地质和成矿背景及可能矿床分布提供了重要信息。     

航空物探综合测量在森林和草原区寻找铜多金属矿的试验研究

通过航空物探(电、磁)综合站在已知铜钼多金属矿床上的飞行测量结果，根据从已知到未知的原则，建立航空物探(电、磁)综合测量在草原、森林、沼泽区寻找斑岩型铜多金属矿的找矿模式。根据此试验模式圈定的航空物探异常，在进行地面查证时取得很好的地质效果。     

Velocity of small-scale auroral ionospheric current systems over Indian Antarctic station Maitri

The Indian Antarctic station Maitri (geog. 70‡45’S, 11‡45’E, geom. 66‡.03S, 53‡.21E) occupies a sub-auroral location during magnetically quiet conditions(Σκρ < 10), but attains an auroral position when the auroral oval shifts equatorwards with increasing strength of magnetic disturbance. At the latter times, triangulation with 3 fluxgate magnetometers located at the vertices of a suitable triangle provides a means of monitoring mobile auroral ionospheric current systems over Maitri. The spacing between the magnetometers is typically kept at 75–200 km, keeping in mind the scale-sizes of ∼ 100 km for these mobile current systems. This work reports the results of two triangulation experiments carried out around Maitri in January 1992 and January 1995, both during Antarctic summer. The velocities estimated for pulsations of the Pc4 and Pc5 type were about 0.59 km/sec in the direction 102‡.7 east of due north, in the first case, and about 1–3 km/sec in the second case in the east-west direction. While several magnetometer arrays exist in the northern auroral regions (e.g., the Alberta array in Canada, the Alaskan array in the U.S. and the IMS Scandinavian array), there is no report in literature of triangulation through arrays in Antarctica, except for a one-day study by Neudegget al 1995 for ULF pulsations of the Pc1 and Pc2 type. The velocities obtained for the Pi3 type of irregular pulsations over Antarctica in the present study tally well with those obtained for northern auroral locations.     

Magnetic anomalies over thin plates and their analysis

The magnetic anomaly over a two-dimensional thin horizontal plate is similar to the first horizontal derivative of the magnetic anomaly over a thick dipping dike of infinite depth extent but with a different direction of magnetization. Hence, the magnetic anomalies of thin plates may be integrated along the profile and the pseudomagnetic potential anomaly thus obtained may be interpreted using any standard method of interpreting dike anomalies. Expressions for the Fourier amplitude and phase spectra of the magnetic anomaly over a thin plate are also derived and procedures to evaluate the parameters of the plate from the spectra are formulated.