Comparison of Factorial Kriging Analysis Method and Upward Continuation Filter to Recognize Subsurface Structures — A Case Study: Gravity Data from a Hydrocarbon Field in the Southeast Sedimentary Basins of the East Vietnam Sea

To interpret geophysical anomaly maps, it is necessary to filter out regional and sometimes noise components. Each measured value in a gravity survey consists of different components. Upward continuation (UC) is one of the most widely used filters. The shortcoming of this filter is not to consider the spatial structure of the data, and also the fact that the trial and error approach and expert’s judgment are needed to adjust it. This study aims to compare the factorial kriging analysis (FKA) and UC filters for separation of local and regional anomalies in the gravity data of a hydrocarbon field in the southeast sedimentary basins of the East Vietnam Sea. As shown in this paper, FKA method permits to filter out all of the identified structures, while the UC filter does not possess this capability. Therefore, beside general and classic filtering methods, the FKA method can be used as a strong method in filtering spatial structures and anomaly component.     

基于重力异常分离方法寻找深部隐伏铁矿 ——以安徽泥河铁矿为例

重力异常是地下不同规模、不同形态和不同埋深的不均匀地质体的综合响应,重力勘探主要通过从重力异常中提取感兴趣的局部异常,以便探测深部结构,寻找隐伏矿床.为探讨重力异常分离原则并检验方法效果,本文从各方法原理入手,加上模型试验以及在安徽省泥河矿区深部隐伏铁矿的探寻实践检验,阐明:趋势分析法是整体拟合不同于最小二乘圆滑的局部拟合,由于是多项式拟合区域场,趋势分析法不适宜范围大、地质情况复杂的测区;插值切割法以计算点场值与四点圆周平均值的插值运算为切割算子,通过连续切割,得到重力异常的切割区域场和局部场,插值切割法对于小测区单个异常的分离效果较好,切割次数选择1到2次即可;匹配滤波法通过分析实测异常功率谱曲线、选择合适的滤波段、建造适宜的低通和带通滤波器进行滤波,从而提取不同波数成分的异常场,匹配滤波更适合垂向叠加的异常分离;解析延拓是根据一个面上的一组位场数据确定另一个不同高度面上位场值,应用中要把握延拓高度;垂向二阶导数法可以起到突出浅源异常,区分水平叠加异常,确定异常体的边界,消除或削弱背景场的作用.通过安徽泥河铁矿重力异常分离实验,发现三阶趋势分析、向上延拓以及插值切割法能很好地分离出矿体异常和背景场,同时发现在泥河矿区东南部和东北部还存在剩余重力异常,可为泥河铁矿扩大规模提供新的线索.     

优化滤波方法及其在中国大陆布格重力异常数据处理中的应用

在优选延拓法的理论基础上,研究提出基于格林等效层概念和维纳滤波器的优化滤波法,用于对重力异常数据进行去噪和分离.与传统向上延拓法和优选延拓法相比,优化滤波法分离异常与延拓高度无关,不需要已知延拓高度,具有一定的优势.理论重力模型数据的去噪和异常分离试验表明优化滤波法有效,异常分离效果优于传统向上延拓法和带通滤波法.利用优化滤波法对中国大陆重力异常数据去噪和异常分离,得到有效的布格重力异常和区域重力异常.以中国大陆深地震探测推断的莫霍面深度信息为约束,对区域重力异常数据进行密度界面约束反演,得到中国大陆莫霍面深度分布.本文方法为中国大陆深部探测和区域构造研究提供一定的技术支撑.     

Investigation of Crustal Thickness in Eastern Anatolia Using Gravity,Magnetic and Topographic Data

The tectonic regime of Eastern Anatolia is determined by the Arabia-Eurasia continent-continent collision. Several dynamic models have been proposed to characterize the collision zone and its geodynamic structure. In this study, change in crustal thickness has been investigated using gravity, magnetic and topographic data of the region. In the first stage, two-dimensional low-pass filter and upward analytical continuation techniques were applied to the Bouguer gravity data of the region to investigate the behavior of the regional gravity anomalies. Next the moving window power spectrum method was used, and changes in the probable structural depths from 38 to 52 km were determined. The changes in crustal thickness where free air gravity and magnetic data have inversely correlated and the type of the anomaly resources were investigated applying the Euler deconvolution method to Bouguer gravity data. The obtained depth values are consistent with the results obtained using the power spectrum method. It was determined that the types of anomaly resources are different in the west and east of the 40° E longitude. Finally, using the obtained findings from this study and seismic velocity models proposed for this region by previous studies, a probable two-dimensional crust model was constituted.     

重力异常分离的小波域优化位变滤波方法

在重力异常分离中,频率域滤波分离方法是以全局数据频谱特征设计针对性的滤波器实现的.滤波器参数与空间位置无关,因此无法针对局部数据频谱特征动态调整滤波器参数.针对该局限性,在小波域滤波理论和优化滤波方法的基础上,本文研究提出了小波域优化位变滤波法,该方法具有空间变化滤波能力,在不同空间位置实现不同的滤波器特性,从而能实现局部数据频谱与全局数据频谱存在较大差异的重力异常分离问题.理论模型数据分离实验结果表明,在局部数据频谱与全局数据频谱差异较大的情况下,该方法相对于Butterworth滤波和优化滤波等方法具有优势.最后,用一个实例进行检验计算,体现了所提方法技术的效果和应用前景.     

Geophysical researches (gravity and magnetic) of the Eratosthenes Seamount in the eastern Mediterranean Sea

New free-air gravity and magnetic maps of the Eratosthenes Seamount and its vicinity were regenerated from potential field data. Stages of data processing are power spectrum, upward continuation, filtering on the free-air gravity anomaly data. RTP, pseudo-gravity transformation map, power spectrum, upward continuation, filtering, AS, and HGAS were applied on the magnetic data. A HGAS map shows the images and locations of the Eratosthenes magnetic body. Spectral analysis of the gravity and magnetic anomalies indicates that there is an elliptical elongated structure of the Eratosthenes Seamount in the width of approx. 86 km NW-SE orientation and in the length of 138 km NE-SW orientation, with a strike of N40°E and inclined to NW. It is considered that 22.49 ± 0.08 km obtained from power spectrum of the gravity data may be related to the crust thickness. Also, 15.67 ± 0.02 km obtained from power spectrum of the magnetic data is considered to be related to the magmatic basement of the Eratosthenes Seamount.     

普光气田重力异常的视密度反演

普光气田是我国近年发现的最大气田．根据气田储层和盖层的岩石密度,估算普光气田上约有零点几毫伽量级的重力异常．但由于区域场和其他岩石的影响,气田的负异常受到严重的干扰．我们使用基于位场切割分离和位场大深度向下延拓的视密度反演新方法,在普光气田上得到明显的负密度异常．在川东北其他一些气田上,也得到负密度异常．根据视密度反演结果,对川东北的含气远景区做了预测．     

基于改进的小子域滤波的重力异常及重力梯度张量边缘增强(英文)

针对利用重磁资料增强地质体边界在图像中的视觉效果和提高地质解译的准确性问题,提出应用改进的小子域滤波方法对重力异常及重力梯度张量数据进行增强处理。根据滑动窗口子域平均选择原理,探讨了改进的小子域滤波方法在位场异常数据含有高斯白噪声时,不同窗口大小对地质体边界的识别效果及其在具有不同边界延伸方向地质体中的应用效果。模型试验结果表明,利用改进的小子域滤波对重力梯度张量数据进行增强处理,得到的地质体边界形态失真更小,且受滤波窗口大小、噪声以及地质体边界方向的影响较小;对深部场源体,通过增大滤波窗口的方式,可以较好地反映深部场源体的边界。结合黑龙江省虎林盆地重力异常以及计算的重力梯度张量的处理实例表明改进的小子域滤波法较传统的小子域滤波法增强了对断裂水平位置信息的识别。     

位场向下延拓组合滤波器的设计和应用

一 重、磁异常向下延拓能将水平方向迭加的异常分开,并能突出异常细节,从而丰富异常的解释内容,提高异常解释的可靠性.本文从对引起向下延拓不稳定的各种因素的讨论入手,结合对一些向下延拓方法在频率域的分析,提出了一种与已知资料频谱特征相“匹配”的组合滤波器。经理论模型和实例考核,说明该滤波器既能保证向下延拓有较高的精度,又有较强的抗干扰能力。     

Internal structure of Tenerife (Canary Islands) based on gravity, aeromagnetic and volcanological data

Gravity and magnetic methods have been applied to the Tenerife Island, to provide new information about its internal structure. For this study, 365 gravity stations covering the central part of the island have been selected. The anomalous density maps at different depths were obtained by means of an inversion global adjustment, on fixed density contrast, to describe the three-dimensional (3D) geometry of the anomalous bodies. On the other hand, several analysis techniques, such as reduction to the pole, spectral analysis, low-pass filtering, terrain correction and forward modelling, were applied to process the high-resolution data obtained in an aeromagnetic survey, completed with marine and terrestrial data.The joint analysis of gravity and magnetic anomalies has shown tectonic and volcanic features that define some fundamental aspects of the structural framework and volcanic evolution of the island. A strong gravity anomaly produced by a large and deep source has been associated with an uplifted block of the basement beneath the southern part of Tenerife. The sources of the observed gravity highs from 8 km b.s.l. may be associated with the growth of the submarine shield stage that was clearly controlled by regional tectonic.The long-wavelength magnetic anomalies reveal highly magnetic sources, interpreted as gabbro-ultramafic cumulates associated with the root zone of a large dyke swarm. This intrusive body could be topped by the emplacement zone of magma chambers that correlate with a magnetic horizon at 5.7±0.8 km depth. Rooted in this highly magnetic zone, two dike–like structures can be associated with the magmatic feeding system of large recent basaltic volcanoes. A shallow magnetic horizon (1.4 km a.s.l.) can be correlated with the bottom phonolites of the Las Cañadas Edifice.In the central part of the island the coincidence of some gravity and magnetic lows is consistent with the presence of low-density and low-magnetic materials, that infill a collapsed caldera system. The structures close to the surface are characterised by low-density areas connected with the recent volcanism, in particular the minimum over the Teide volcano. Hydrothermal alteration is assumed to be the cause of a short-wavelength magnetic low over the Teide volcano.     

补偿向下延拓方法研究及应用

位场向下延拓是重、磁处理和解释的常用方法,但其不稳定性限制了其在资料处理及反演中的应用.本文基于补偿圆滑滤波思想以及空间域向下延拓迭代法,通过逐次补偿的办法实现位场的稳定向下延拓.同时,在频率域空间给出了该下延方法的频率域响应因子,并讨论了其低通滤波特性,理论模型和实际位场资料试验表明该方法向下延拓稳定性具有较高的延拓精度.将其应用于重力密度界面反演中,改进反演的稳定性,实际莫霍界面反演表明下延因子具备实用性.     

Interpretation of Gravity Anomalies with the Normalized Full Gradient (NFG) Method and an Example

The Normalized Full Gradient (NFG) method which was put forward about 50 years ago has been used for downward continuation of gravity potential data, especially in the former Union of Soviet Socialist Republics. This method nullifies perturbations due to the passage of mass depth during downward continuation. The method depends on the downwards analytical continuation of normalized full gradient values of gravity data. Analytical continuation discriminates certain structural anomalies which cannot be distinguished in the observed gravity field. This method has been used in various petroleum and tectonic studies. The Trapeze method was used for the determination of Fourier coefficients during the application of this method. No other techniques for calculating these coefficients have been used. However, the Filon method was used for the determination of Fourier coefficients during the application of the NFG method in this work. This method, rather than the Trapeze method, should be preferred for indicating abnormal mass resources at the lower harmonics. In this study, the NFG method using the Filon method has been applied the first time to theoretical models of gravity profiles as example field at the Hasankale-Horasan petroleum exploration province where successful results were achieved. Hydrocarbon presence was shown on the NFG sections by the application of NFG downward continuation operations on theoretical models. Important signs of hydrocarbon structure on the NFG section for field and model data at low harmonics are obtained more effectively using this method.     

火山岩油气藏重磁电震综合预测方法及应用

通过准噶尔盆地陆东地区数十口钻井资料的对比分析、归纳总结,提出了火山岩油气藏重磁电震综合预测方法.将正则化下延与延拓回返垂直二次导数串联形成了一个新的滤波器,该滤波器相当于首先通过正则化下延将位场曲面延拓至地下目的层段,降低火山岩埋深对磁J异常幅值的影响,然后利用延拓回返垂直二次导数对弱信号进行增强,不仅提高了位场异常...     

2D and 3D potential-field upward continuation using splines

The dominant upward‐continuation technique used in the potential‐field geophysics industry is the fast Fourier transform (FFT) technique. However, the spline‐based upward‐continuation technique presented in this paper has some advantages over the FFT technique. The spline technique can be used to carry out level‐to‐uneven surface 2D and 3D potential‐field upward continuation. An example of level‐to‐uneven surface upward continuation of 3D magnetic data using the spline technique is shown, and it is evident that the continued anomalies are very close to the theoretical values. The spacing can be irregular. Synthetic examples using the spline technique to continue noise‐contaminated gravity and magnetic data upward to an altitude of 15 km on irregular grids are shown. Gaussian noise with a zero mean and a standard deviation of 1% does not cause much error and can readily be tolerated. Through comparison with the FFT technique, it is found that for low‐altitude gravity and magnetic upward continuation, both the FFT technique and the spline technique are suitable; for high‐altitude upward continuation, the FFT technique is inaccurate, whereas the spline technique works very well. Also, upward continuation by the spline technique has a smaller edge effect than upward continuation by the FFT technique. The spline‐based upward continuation technique works fairly well even when the periphery of a grid is not quiet: it is rather robust in general. A real example shows that the spline technique can be employed to perform upward continuation of total‐field magnetic data and to de‐emphasize near‐surface noise.     

位场异常分离的递减半径迭代法

本文提出了位场区域-剩余异常分离的空间域递减半径迭代法.由给定半径的圆周上八点位场值的算术平均导出一个新的八点圆周平均公式,它是一个由重磁数据计算区域异常的滤波器.该滤波器的传递函数有一主瓣和多个旁瓣,且半径越大,旁瓣越多,滤波器特性越差.由从大到小不同半径对应的传递函数的连乘积构造了递减半径迭代传递函数,它以大半径为参变量.递减半径迭代传递函数类似低通滤波器,其截止波数与大半径成反比.由递减半径迭代传递函数,给出空间域分离区域异常的递减半径线性迭代法,由重磁数据减去区域异常求得剩余异常.进一步,通过构造非线性修正系数,把空间域递减半径线性迭代法中线性迭代公式变成非线性迭代公式,得到空间域递减半径非线性迭代法.重磁理论模型及安徽泥河铁矿重磁资料的试验表明,空间域递减半径非线性迭代法可有效地压制假异常和高频干扰以及减少异常畸变,并能有效地分离出区域异常和剩余异常.     

AUTOMATIC DECONVOLUTION OF GRAVIMETRIC ANOMALIES*

Existing techniques of deconvolution of gravity anomalies are principally based on upward and downward continuation of measured fields. It can be shown that a unique set of linear filters, depending only on geometrical parameters, relates density distribution at a given depth to gravity measured on the surface. A method to compute the filter coefficients is developed. Very accurate reconstitution of theoretical models of intricate shape, prove the validity of the linear relationship. One of these sets of linear filters is applied to a field case of underground quarries.     

Application of the finite element method to gravity data case study: Western Turkey

Gravity anomalies always include the total effects (combination of the structures which have different densities and depths) of the study area and beyond. And the well-known non-uniqueness of potential field modelling may lead to very different interpretation results.The finite element method (FEM), which has been used in potential field interpretation for decades, makes complex problems to be solved easily and accurately. The first step of FEM is to identify the elements and then to decide on the boundary of the solution space. In this step, the solution space is divided into elements. After determination of the geometrical structure of the solution space, the most suitable elements should be chosen for this geometrical structure. The agreement between the geometry and the elements is quite important for the convergence to the best possible solution.In this work, the methods of trend analysis, filtering, analytical continuation and FEM were applied to a theoretical model and to gravity data from western Turkey to produce the regional and the residual anomalies. The results were compared and it was found that the FEM produced more accurate results than other methods did.     

THE ANALYSIS AND DESIGN OF TWO-DIMENSIONAL FILTERS FOR TWO-DIMENSIONAL DATA *

Some of the methods such as regional removal and second derivative calculations which can be used to outline anomalies on potential data maps can be thought of as a filtering operation. The analysis and design of such two-dimensional filters by means of direct and inverse two-dimensional Fourier transforms have been considered. An analysis of several published sets of second derivative coefficient sets indicates that, in general, they are not a good approximation to the theoretical second derivative filter. Alternate methods of designing regional removal and second derivative filters are discussed. The properties of various two-dimensional filters are further illustrated by means of maps obtained from the convolution of several of these filters with a set of observed field data. These maps show the large changes in anomaly shape which can result from the inclusion or rejection of various wavelength components.     

INTERPRETATION QUANTITATIVE EN GRAVIMETRIE OU MAGNETISME A PARTIR DE CARTES TRANSFORMEES DE GRADIENT VERTICAL*

A discussion is given of the requirements, the advantages and the methods to be considered in attempting the quantitative interpretation of gravity or magnetic fields from computed maps of the vertically derived field. The transform which is used here is the first vertical derivative (or vertical gradient) with or without downward continuation, but the computed maps are in fact obtained by a controllable Fourier method in which two kinds of operations can be simultaneously performed in complete independance: on one hand the separation of any part of the data by frequency cut-off, and on the other hand the transformation by vertical derivation or continuation of the part which is retained. Taking as raw data either actual surveys or artificially constructed maps, it is first shown how separated and transformed maps of this type can be efficiently obtained under quite flexible conditions, using a special computer program. It is further seen that for correctly controlled filterings the accuracy of the computed maps actually permits to take them as the basis for quantitative interpretation. To effect this, any one of the conventional methods which make use of equivalent model computations may in the first place be adapted to the interpretation of gradients, with the benefits, however, of an enhanced lateral separation of the anomalies and of a large attenuation of the regional effects. Particularly, the delineation of horizontal contours for even fairly complex models can often be made directly in a sufficiently safe way on the anomalies as they show on the vertical gradient maps. This greatly accelerates the process of determining equivalent model bodies. More special methods of interpretation can also be designed by taking into account first the fact that the vertical derivation of the field amounts to an operation of separation on the field's sources themselves, and in addition the availability of the frequency form of the information as a result of using a Fourier method of transform computation. Trial utilisations of various such interpretation processes of either the conventional or the less conventional type, are presented especially in connection with an interpretation study on transformed maps of parts of an offshore aeromagnetic survey (English Channel).     

非均匀物性条件下多尺度窗口修正法换算磁源重力异常及在寻找DSO的应用

DSO(Direct Shipping Ore,直运块矿)是一种高品位的富铁矿.在加拿大拉布拉多(Labrador)地槽谢菲尔威利(Schefferville)铁矿成矿带,含铁建造苏克曼(Sokoman)组地层全铁含量低、具高密度、强磁性,能够引起高重力异常与高磁异常;而风化淋滤后富集的赤铁矿和针铁矿等(也称DSO)全铁含量高,具高密度、无磁性,仅能够引起高重力异常.采用一般的滤波方法不能提取DSO的重、磁异常.本文采用基于泊松(Poisson)公式的磁场换算磁源重力异常(pseudo-gravity anomaly)方法,由磁场换算磁源重力异常,再与实测重力异常对比,得到纯粹由高密度、无磁性的DSO产生的剩余重力异常,对剩余重力异常采用密度成像与2.5D反演方法解释DSO.泊松公式虽然提出时间很长,但迄今为止仅仅用在资料解释中的定性分析,本文推导并实现了密度磁性非均匀条件下经典泊松公式的形式与实现过程,提出了多尺度窗口滑动线性回归修正的磁场换算磁源重力异常方法,使该公式的数学原理能够对重、磁异常的反演解释定量化.最后本文将多尺度窗口滑动线性回归修正的换算磁源重力异常方法用于加拿大拉布拉多地槽谢菲尔威利铁矿成矿带铁矿勘探,较好地解决了寻找高密度、无磁性DSO的问题.