THE QUALITY OF SOME TWO-DIMENSIONAL FILTERS IN SEPARATING REGIONAL AND LOCAL GRAVITY ANOMALIES *

The separation of regional from local gravity anomalies by means of the application of two-dimensional linear filters is analyzed. It was found that optimization of the filter in the least squares sense leads to filters that produce strong localized concentrations of the error, which may erroneously be interpreted as anomalies. For this reason the maximum absolute value of the error is a more important criterion for the quality of the filter than the root mean square error. This maximum absolute error is minimized by the minimax filter. Intermediate filters are derived which give a transition zone which comes appreciably closer to that of the optimal filter at only a small price in terms of increase of the maximum absolute error.     

A GENERALIZED METHOD OF COMPUTING SECOND DERIVATIVE OF GRAVITY FIELD*

The unsuitability of a particular method of computation of second derivative to varying types of field problems evinces the need of a method by which it may be possible to obtain a wide range of filters. This requirement is achieved by incorporating a mathematical smoothing operator, e^-λρ2, λ and ρ being the smoothing parameter and radial frequency respectively, in the theoretical frequency response of second derivative operation. A few sets of weights for various smoothing parameters, are presented. The analysis of Bouguer gravity anomaly map of Los Angeles Basin, California, U.S.A., is carried out by Elkins's (1951) and Rosenbach's (1953) methods and the sets of weights given in this paper which clearly establishes their filtering equivalence.     

A NOTE ON CONVOLUTION AND PADDING FOR TWO-DIMENSIONAL DATA*

Gupta , V. K., and Ramani , N., 1978, A Note on Convolution and Padding for Two-Dimensional Data, Geophysical Prospecting 26, 214–217. Filtering of gridded data by means of two-dimensional generalized convolution results in a loss of data along the edges of the map. This loss can be counteracted by ‘padding’– filling the surrounding of the input map with some kind of data. Padding by algebraic extrapolation is assumed to be better than surrounding the maps with zeroes. A Føsrtran IV routine for padding by algebraic extrapolation is available from the authors.     

A CONTROL OF TWO-DIMENSIONAL MAGNETIC INTERPRETATION BY THREE-DIMENSIONAL MODEL BODY ANOMALIES*

In magnetic routine interpretation the comparison of two-dimensional model curves with measured magnetic anomalies is widely used for an approximate evaluation of the position and depth of magnetic models. Before starting an interpretation of a survey by means of two-dimensional models, it is very useful to have an idea of the shape of anomalies caused by extended but finite bodies, taking into account various strike directions: Three sets of anomalies of thin plates (horizontal length 19, downward length 9, width 1) dipping 30°, 60°, and 90° resp. for various strike directions and an inclination of 20° were computed. Some of these anomalies, e.g. those with nearly N-S strike direction look rather complicated, and at the first glance one would not expect that they are caused by such simple bodies. Several profiles crossing the computed anomalies perpendicularly were interpreted two-dimensionally. For less extended anomalies the depths determined for the top of the plates are 10-20% too small, the magnetization amounts to 50–75 % of the value of the finite bodies. The interpretation of the profiles covering more extended anomalies gave very accurately the same values for the position, depth and magnetization for the two-dimensional body as for the original three-dimensional model. Anomalies of vertical prisms with varying extensions in the y-direction were computed. Their differences in amplitude and in the distance maximum-minimum show that interpretation of short anomalies by two-dimensional methods yields depth errors of up to 20 percent. To see the possibilities of the separation of superimposed anomalies dike anomalies were added to the anomaly of a broad body in great depth and several attempts were made to interpret parts of the composite anomalies. The interpreted bodies lie too deep. In complicated cases the depth values can have large errors, but experienced interpreters should be able to keep the errors in the range of one third of the depth values.     

Regional correction of gravity data

Summary Methods for the elimination of the regional effect from gravity data are discussed in this paper. A historical review of the average gradient method used in torsion balance data is given. The connection between the averaging, or grid method, used on the gravity meter data and the average gradient method is pointed out and a criticism of these methods is given. Due to the ambiguity of the interpretation of gravity data no unique determination of the regional effect is possible and no mechanical method for its removal can be found. The determination of a regional gravity effect is an interpretation problem and for this reason it is unavoidably the subject of geological consideration and of the personal judgement of the interpreter. The method of smoothing the isogams and gravity curves is also discussed. A few suggestions are made for breaking the gravity data into two components (regional and residual) by the analysis of the gravity data. Finally the application of the second derivative methods for the elimination of the regional effect is discussed.     

航空重力测量数据的小波滤波处理

构造三类连续小波函数对航空重力测量数据进行小波滤波处理. 三类连续小波函数分别用于对测量数据在某一空间尺度（或时间尺度）上的低通,一阶求导和二阶求导滤波. 着重介绍三类连续小波函数的构造原理与过程,并说明其相对于传统的数字滤波器的优势. 对系统的技术参数（滤波器窗口宽度参数δ 和尺度参数s）进行了调试实验. 实算结果显示了方法的可行性和有效性.     

BASIC PRINCIPLES OF TWO-DIMENSIONAL DIGITAL FILTERING *

The geophysicist involved in the analysis of two-dimensional data should have an understanding of the two-dimensional finite Fourier transform and the mechanics of two-dimensional filtering. Frequency aliasing must be considered when working with sampled data. In two dimensions it is advantageous to consider aliasing in terms of the overlap of the repeating spectra inherent in the finite Fourier transform. Two-dimensional filtering can be performed as a transient convolution in the space domain, as cyclic convolution utilizing the frequency domain or as the multiplication of polynomials using the z-transform. If the “edge” effects are removed, the results of the three methods are identical.     

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).     

GEOELECTRICAL MODEL CALCULATIONS FOR TWO-DIMENSIONAL RESISTIVITY DISTRIBUTIONS*

For the calculation of geoelectrical model curves for a two-dimensional resistivity distribution, the potential equation is transformed by means of a Fourier cosine transform into a two-dimensional Helmholtz equation containing the separation parameter λ. The numerical solution of this equation for different values of λ for an irregular grid is obtained using the method of finite differences combined with the method of overrelaxation. The method by which derivatives are replaced by finite differences turned out to be very important, especially for high resistivity contrasts. After testing several methods designed to deal with any type of resistivity distribution, a method of discretization similar to that used by Brewitt—Taylor and Weaver (1976) for magnetotelluric modeling for H polarization was found the best. Examples are given of model curves for Schlumberger soundings over a vertical fault covered by overburden. The incorrect use of horizontal-layer models leads to erroneous interpretations that are more complex than the real subsurface situations.     

Improved edge detection mapping through stacking and integration: a case study in the Bathurst Mining Camp

Airborne geophysical surveys provide spatially continuous regional data coverage, which directly reflects subsurface petrophysical differences and thus the underlying geology. A modern geologic mapping exercise requires the fusion of this information to complement what is typically limited regional outcrop. Often, interpretation of the geophysical data in a geological context is done qualitatively using total field and derivative maps. With a qualitative approach, the resulting map product may reflect the interpreter's bias. Source edge detection provides a quantitative means to map lateral physical property changes in potential and non‐potential field data. There are a number of Source edge detection algorithms, all of which apply a transformation to convert local signal inflections associated with source edges into local maxima. As a consequence of differences in their computation, the various algorithms generate slightly different results for any given source depth, geometry, contrast, and noise levels. To enhance the viability of any detected edge, it is recommended that one combines the output of several Source edge detection algorithms. Here we introduce a simple data compilation method, deemed edge stacking, which improves the interpretable product of Source edge detection through direct gridding, grid addition, and amplitude thresholding. In two examples, i.e., a synthetic example and a real‐world example from the Bathurst Mining Camp, New Brunswick, Canada, a number of transformation algorithms are applied to gridded geophysical data sets and the resulting Source edge detection solutions combined. Edge stacking combines the benefits and nuances of each Source edge detection algorithm; coincident or overlapping and laterally continuous solutions are considered more indicative of a true edge, whereas isolated points are taken as being indicative of random noise or false solutions. When additional data types are available, as in our example, they may also be integrated to create a more complete geologic model. The effectiveness of this method is limited only by the resolution of each survey data set and the necessity of lateral physical property contrasts. The end product aims at creating a petrophysical contact map, which, when integrated with known lithological outcrop information, can be led to an improved geological map.     

任意曲线上二维重磁位场转换的B样条函数法

重磁位场转换在重磁资料解释中是十分重要和必要的。随着物探仪器的改进和资料解释精度要求的提高,人们对位场转换方法的研究越来越深入。本文提出一种原理简单、计算精度高的新方法--B样条函数法,应用其插值、求导、求积的性质,以解决曲线上二维重磁位场转换的两类问题。一是Z_α分量向H_α分量的转换及向上延拓,二是水平和垂向导数的换算及磁位的计算。 样条函数为分段多项式函数。B样条函数是用σ函数逼近的一类样条函数,其节点     

Corrections for near surface effects: Decomposition of the magnetotelluric impedance tensor and scaling corrections for regional resistivities: A tutorial

This paper primarily examines the effects of small-scale or near-surface conductivity inhomogeneities on the magnetotelluric (MT) impedance tensor. These effects cause three different types of distortion results. (1) The well-knownstatic shifts of sounding curves. (2) When the underlying regional setting is two-dimensional then the two regional impedances are mixed in an arbitrary coordinate system. Thus the level and shape of each sounding curve is distorted as are the phases. (3) At sufficiently high frequencies these effects generate anomalous magnetic fields that in turn alter the background phases.This tutorial first explores the usefulness of various MT tensor analysis techniques to overcome the problem of phase mixing and to recover regional information in the presence of local geological noise. Synthetic and experimental data are considered. A sequence ofa priori models of increasing complexity are described. The use of appropriate decompositions of the MT tensor each with an increasing number of parameters is emphasised. In a second part, phase mixing and static shifts are examined from a synoptic view. Some static shift removal techniques that can be used in conjunction with the decomposition are discussed. This paper is not a review but rather an investigation of a few methods that the authors have found useful with field data.     

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.     

A FREQUENCY DOMAIN APPROACH TO TWO-DIMENSIONAL MIGRATION *

The problem of the propagation of acoustic waves in a two-dimensional layered medium can be easily solved in the frequency domain if the Dix approximation is used, i.e. when only the primary reflections are considered. The migrated data at a depth z are obtained by convolving the time section with a proper two-dimensional operator dependent on z. The same result can be obtained by multiplying their two-dimensional spectra and summing for all the values of the temporal frequency. The aspect of the operator in the time-space domain has the classic hyperbolic structure together with the prescribed temporal and spatial decay. The main advantages of the frequency domain approach consist in the noticeable computer time savings and in the better approximation. On the other hand lateral velocity variations are very difficult to be taken into account. This can be done if a space variant filter is used in the time-space domain. To reduce computer time, this filter has to be recursive; the problem has been solved by Claerbout by transforming the hyperbolic partial differential equation into a parabolic one, and using the latter to generate the recursion operator. In the presentation a method is given for the generation of recursive filters with a better phase characteristics that have a pulse response with the requested hyperbolic shape instead of the parabocli one. This allows a better migration of steeper dips.     

A STUDY OF DIFFERENT METHODS OF WAVEFIELD SEPARATION FOR APPLICATION TO VSP DATA1

During the last couple of years there has been much research in the area of wavefield separation of borehole seismic data, and several articles have been published on various separation techniques. Methods involving the application of two-dimensional Fourier transformation, the Radon transformation, multi-level median filters or optimal filters, are all suggested as possible approaches to the wavefield separation problem. This paper compares some of these methods commonly used in the industry. The theories of the chosen methods are described to see how they are related. Using the different methods on synthetic and real data, we show how this theoretical relation is reflected in the relatively similar results obtained. We also show how the different filters treat coherent and random noise.     

基于磁异常的边界特征增强方法对比研究

基于磁异常进行场源边界识别(如断裂划分)是磁法勘探解释工作的重要内容.由于受到斜磁化、场源埋深、异常叠加以及噪声等多种因素的影响,基于磁异常直接进行场源的边界识别往往难以准确识别出场源体的边界,所以,人们研究了各种场源边界特征增强的方法技术,以便于进行场源边界的识别工作.近年来,国内外出现了很多新的边界特征增强方法,但有些方法的边界特征增强效果与已有方法相近,有的方法效果有待深化检验,存在验证方法的理论模型过于简单,未全面细致地分析出它们的优缺点等问题.针对上述情况,本文运用多种复杂、贴近实际的理论模型,深入对比分析了包括解析信号模法和倾斜角法等在内的15种具有代表性的边界增强方法,分析了它们在不同磁化方向、不同场源埋深、不同场源形态、异常叠加和噪声干扰等多种因素影响下的边界增强效果,重点归纳总结出它们的优缺点和适用条件,以为实际应用时的方法选择提供参考.通过模型实验,我们认为解析信号模垂向导数法、总水平导数法、解析信号模倾斜角法是在多种因素影响下,适用性较强、应用效果较好的三种方法.     

Toward a full multiscale approach to interpret potential fields

The way potential fields convey source information depends on the scale at which the field is analysed. In this sense a multiscale analysis is a useful method to study potential fields particularly when the main field contributions are caused by sources with different depths and extents. Our multiscale approach is built with a stable transformation, such as depth from extreme points. Its stability results from mixing, in a single operator, the wavenumber low‐pass behaviour of the upward continuation transformation of the field with the enhancement high‐pass properties of n‐order derivative transformations. So, the complex reciprocal interference of several field components may be efficiently faced at several scales of the analysis and the depth to the sources may be estimated together with the homogeneity degrees of the field. In order to estimate the source boundaries we use another multiscale method, the multiscale derivative analysis, which utilizes a generalized concept of horizontal derivative and produces a set of boundary maps at different scales. We show through synthetic examples and application to the gravity field of Southern Italy that this multiscale behaviour makes this technique quite different from other source boundary estimators. The main result obtained by integrating multiscale derivative analysis with depth from extreme points is the retrieval of rather effective information of the field sources (horizontal boundaries, depth, structural index). This interpretative approach has been used along a specific transect for the analysis of the Bouguer anomaly field of Southern Apennines. It was set at such scales, so to emphasize either regional or local features along the transect. Two different classes of sources were individuated. The first one includes a broad, deep source with lateral size of 45∼50 km, at a depth of 13 km and having a 0.5 structural index. The second class includes several narrower sources located at shallowest depths, ranging from 3–6 km, with lateral size not larger than 5 km and structural indexes ranging from 1–1.5. Within a large‐scale geological framework, these results could help to outline the mean structural features at crustal depths.     

二维最佳线性数字滤波器的设计原理

针对如何在干扰场的背景上区分出低缓异常,以及在位场的向下延拓一类计算中如何限制因误差的高频放大所导至的解的不稳定性等问题,本文探讨了在“最小二乘”意义下的最佳线性数字滤波器的设计原理,并将它转化为下述数学问题,即在L₂线性赋范函数空间中如何选取最佳滤波函数的问题。在空间域中直接解这个问题是十分复杂和困难的,我们发现在波数域中用变分法中的等周问题的解法直接选取最佳线性滤波器的传输函数（或波数响应）,则在数学方法上既简单又严格。这样选取的最佳线性滤波器的传输函数L（f,k）其表达式也很简单,即L（f,k）=|S_i（f,k）|²/{|S_i（f,k）|²+λ|N_i（f,k）|²}。式中,|S_i（f,k）|²及|N_i（f,k）|²分别代表滤波器输入端讯号和干扰的能谱（或功率谱）,f、k分别代表x、y方向上的波数,λ为大于零的常数。 对上述两类问题以及相关的两种最佳线性滤波器而言,L（f,k）的表达式是相同的,而区别仅在于其参变量λ的选取条件不同而已。 有了最佳线性滤波器的传输函数L（f,k）的理论公式,就可以在最小二乘的意义下分析和评价国内外所发表的解决上述两类问题的各种线性滤波方法,并能指出在不同的讯号与干扰条件下,在理论上线性滤波可能达到的最佳效果,从而为设计二维线性数字滤波器时,提供一个理论上的准则。 对位     

由二维破裂模式导出的地震定标律

使用二维矩形破裂模式,导出了大小不同地震的震源位移谱的表达式.位移谱有三个拐角频率fc1,fc2,fc3.fc1和fc2分别联系着长度方向和宽度方向的破裂时间,fc3联系着震源函数的上升时间.根据三个拐角频率,可以把位移谱u（）分成四个区域,在Ⅰ区,fc1,u（）=u（0）f0;在Ⅱ区,fc1c2,u（）1/f;在Ⅲ区,fc2c3,u（）1/f2;在Ⅵ区,f >fc3,u（）1/f3.由于在四个区域u（）随频率的增加下降的速度不一样,决定了震源参数在不同的震级范围内（也就是不同地震矩范围）有不同的表达式.假定地震满足几何相似、应力环境相似以及动力学相似条件,因而地震矩M0,长度和宽度方向的破裂时间、上升时间都可以用断层长度L来定标。根据Dziewonski and Woodhouse给出的1981——1983年800多个地震的地震矩M0资料,以及BISC给出的面波震级Ms、体波震级mb,来确定定标律中的常数.这样就可以从定标律推导出震源参数之间的统计关系.      

柴达木盆地航磁资料微弱信息增强技术研究及在线性构造识别中的应用

目的利用重磁异常数据微弱信息增强技术识别线性构造的方法,并介绍其实现过程.方法文中提出了一种梯级带滤波增强非线性滤渡技术与Tilt梯度及其水平导数的有效结合,增强放大微弱信息再识别提取线性构造,利用数字图象显示技术成图.结果识别出柴迭木盆地线性构造多条,分析对比柴达木盆地区域地质构造资料及重力解释成果,具有较好的吻合性.为深入研究该地区线性构造、成矿特征、寻找勘探靶区补充了新的证据.结论该方法提取了区域航磁异常资料中的微弱信息,弥补传统方法的不足,对断层边界及异常的边界划定更为准确.梯级带滤波增强技术与Tilt梯度及其水平导数相结合取得良好的解释效果.