Optimal Design of a Gravity Survey Network and its Application to Delineate the Jabera-Damoh Structure in the Vindhyan Basin,Central India

Fractal dimension analysis was carried out for optimal designing of 2-D gravity survey network and to determine an optimum range of gridding interval to generate least aliased Bouguer anomaly maps. As a test case, this method has been successfully applied to the Jabera-Damoh region of the Vindhyan Basin, which is considered as a potential hydrocarbon bearing area. In particular, we aim to delineate accurately the lateral extent of a possible hydrocarbon bearing structure. To achieve this aim, fractal dimension of survey network was computed using 2-D distributions of observation points in the planning phase of the survey so that the optimum station spacing for gravity survey can be obtained. A range of optimum gridding interval for the gravity data set was suggested using the box-counting method of fractal dimension determination. Bouguer anomaly maps of the region are prepared utilizing the optimum gridding interval. For the first time, these anomaly maps clearly outline the gravity evidence of an anomalous rifted structure, which is bounded by parallel faults on either side. This structure is interpreted as a favorable basin for the occurrence of hydrocarbons. Another finding of this study has been the delineation of an apparently small ridge-like structure running east-west, dividing this basin in two parts. A subsurface geological model along a profile across the Jabera structure has also been presented.     

Error estimation and optimization of gravity surveys1

Gravity data are often acquired over long periods of time using different instruments and various survey techniques, resulting in data sets of non-uniform accuracy. As station locations are inhomogeneously distributed, gravity values are interpolated on to a regular grid to allow further processing, such as computing horizontal or vertical gradients. Some interpolation techniques can estimate the interpolation error. Although estimation of the error due to interpolation is of importance, it is more useful to estimate the maximum gravity anomaly that may have gone undetected by a survey. This is equivalent to the determination of the maximum mass whose gravity anomaly will be undetected at any station location, given the data accuracy at each station. Assuming that the maximum density contrast present in the survey area is known or can be reasonably assumed from a knowledge of the geology, the proposed procedure is as follows: at every grid node, the maximum mass whose gravity anomaly does not disturb any of the surrounding observed gravity values by more than their accuracies is determined. A finite vertical cylinder is used as the mass model in the computations. The resulting map gives the maximum detection error and, as such, it is a worst-case scenario. Moreover, the map can be used to optimize future gravity surveys: new stations should be located at, or near, map maxima. The technique is applied to a set of gravity observations obtained from different surveys made over a period of more than 40 years in the Abitibi Greenstone Belt in eastern Canada.     

Mean Free-Air Gravity Anomalies in the Mountains

Mean free-air gravity anomalies are often needed in geodesy for gravity field modelling. Two possible ways of compiling the mean free-air gravity anomalies are discussed. One way is via simple Bouguer gravity anomalies and the second and more time consuming way is via refined Bouguer gravity anomalies. In flat areas the differences between using any of the two ways should not be significant. In the mountains however, every effect introducing a high dependency, such as e.g. terrain effect, can negatively affect the interpolation process. In fact, a numerical experiment conducted in one part of Rocky Mountains revealed large and systematic differences. The effect of these differences on the geoid model is more then two meters in the test area. Our investigation shows that this bias is caused by the location of gravity measurement points, chosen mostly on hill-tops. At such points, the terrain correction to gravity is systematically larger than the mean value of the correction. Therefore, it is not possible to prevent the mean free-air gravity anomalies obtained from simple Bouguer gravity anomalies from having a systematic bias. One can see this bias as a result of the aliasing effect because the simple Bouguer gravity anomalies in the mountains contain a higher frequency signal (terrain effect) that is, according to sampling theorem, impossible to reconstruct by sparse measured gravity data, see e.g. (Goos et al., 2003). Therefore, the more rigorous way of computing the mean free-air gravity anomalies is via refined Bouguer gravity anomalies.     

分形插值对共炮集记录叠前逆时偏移的影响

利用分形插值方法对稀疏采样的地震记录进行了插值加密处理,并对原始记录和分形插值重建记录进行了叠前逆时偏移处理．数值实验表明,分形插值重建的地震图是原始地震图的良好近似,而分形插值叠前偏移剖面是原始记录偏移剖面的良好近似,二者的分辨率几乎一样．利用分形插值重建技术可以在节约经济成本的情况下,得到地下结构的高精度成像结果．      

Monte Carlo Studies of Relations between Fractal Dimensions in Monofractal Data Sets

—Within the fractal approach to studying the distribution of seismic event locations, different fractal dimension definitions and estimation algorithms are in use. Although one expects that for the same data set, values of different dimensions will be different, it is usually anticipated that the direction of fractal dimension changes among different data sets will be the same for every fractal dimension.¶Mutual relations between the three most popular fractal dimensions, namely the capacity, cluster and correlation dimensions, have been investigated in the present work. The studies were performed on the Monte Carlo generated data sets. The analysis has shown that dependence of the fractal dimensions on epicenter distribution, and relations among the fractal dimensions, are complex and variable. Neither values nor even inequalities among dimension estimates are preserved when different fractal dimensions are used. The correlation and the capacity dimensions seem to be good tools to trace collinear tendencies of eipicenters while the cluster dimension is more appropriate to studying uniform clustering of points.     

地表负载及GPS测站分布对参考框架转换的影响分析

GPS数据处理通常采用Helmert七参数转换将瞬时站坐标转换到指定的框架下,但瞬时站坐标中尚未模型化的季节性地表负载会影响平移参数(地心运动)和尺度参数的估值,进而影响测站坐标;不均匀的测站分布会加剧这一影响.本文利用GRACE重力场系数仿真GPS地表负载的实验表明,基于网平移法采用实际的IGS站至少能够恢复90%的地心运动信号.地表负载及GPS实际测站的不均匀分布可以解释大约30%的GPS尺度的周年变化.相对于IGb08的所有框架站,目前采用91个全球均匀分布的核心站作为框架转换的基准是合理的.采用IGb08的所有框架站进行转换会导致U方向误差增加,特别是对框架站密集的欧洲区域(误差均值约为1mm).因此框架转换时,应尽量选取均匀分布的测站,同时不估计尺度参数.     

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.     

RADON TRANSFORM APPLICATION TO THE IMPROVED GRIDDING OF AIRBORNE GEOPHYSICAL SURVEY DATA1

The Radon transform is applied to airborne geophysical data, which consist of parallel profiles, analogous to a seismic record. The plane-wave decomposition (PWD) thus becomes the strike-direction decomposition (SDD) since the observed spatially distributed information is represented by its strike directions in a domain achieved by the transformation. It is important that, after the SDD, we can identify anomalies and work on them according to their strikes. In particular, for gridding purposes, we may guide the second interpolation of the bi-directional gridding approach along the strike directions. In principle, the proposed Radon transform gridding method (RTGM) transforms the observed parallel profiles into a domain where information is mapped as its strike-direction ‘traces’ against its wavelengths. The number of strike directions into which the data are decomposed is equal to the number of lines to be interpolated. The Fourier spectrum of the grid is reconstructed from the strike-wavenumber domain by using the projection-slice theorem and the final square grid is obtained by performing an inverse Fourier transformation on the spectrum. The SDD is restricted to the Nyquist wavenumber bandwidth imposed by the survey line-spacing, so that there is no addition of ambiguous short wavelengths in the gridded data. A tapering window is employed to prevent any Gibb's oscillation in the final grid because of the sharp Nyquist cut-off in the reconstructed spectrum due to the survey line-spacing. The RTGM is first tested on a set of synthetic line-based data. It is also applied to aeromagnetic profile data from northern Botswana as a practical example.     

重震反演中国东北地壳上地幔三维密度结构

本文利用重力和地震P波到时数据反演得到了中国东北地区地壳上地幔三维密度结构.与单一的重力或地震反演相比,重震反演一方面有效地克服了重力反演结果垂向分辨率低的问题,另一方面也提高了地震反演结果的可靠性.结果显示:中国东北地区的地壳及上地幔剩余密度异常分布与构造单元具有明显的相关性,造山带对应低密度异常,盆地对应高密度异常;区域内火山下方有明显的低密度体存在,可能是由于太平洋板块俯冲进入上地幔并部分滞留,在滞留板块深部脱水和软流圈热物质共同作用下产生了上涌岩浆,喷发后形成了火山.     

离散时变重力数据的可视化、指标量定义与解释

针对流动重力观测获得的数据成果特点,提出一种基于流动重力段差变化的时变重力数据可视化方法,并定义了两个指标量G值和C值,用以评价区域性重力场变化的显著性程度.在此基础上,应用该方法处理和分析了首都圈地区的流动重力数据.结果表明:与传统采用等值线方式来刻画时变重力场特征不同,该方法更能突出发生变化的重力测点位置、测量误差及其变化的显著性程度,可为研究与地震孕育、发生过程有关的重力场变化提供更多的定量依据.      

地壳垂直形变的分维特征─以京、津、唐复测网为例

本文应用分数维的概念及方法，对京、津、唐地区１９６８─１９７５年的大面积地壳垂直形变复测成果数据处理后的结果进行了分维特征分析。结果表明，地震前该地区地壳垂直形变空间分布的分数维明显降低，地壳的有序性增强。本文认为，地壳垂直形变空间分布分数维能够明确地揭示地壳的有序程度，是一种较好地利用形变测量资料提取地震前兆异常的方法。     

A practical review of geostatistical processing applied to geophysical data: methods and applications

Nowadays, geostatistics is commonly applied for numerous gridding or modelling tasks. However, it is still under used and unknown for classical geophysical applications. This paper highlights the main geostatistical methods relevant for geophysical issues, for instance to improve the quality of seismic data such as velocity cubes or interpreted horizons. These methods are then illustrated through four examples. The first example, based on a gravity survey presents how a geostatistical interpolation can be used to filter out a global trend, in order to better define real anomalies. In the second case study, dedicated to refraction surveying, geostatistical filtering is used to filter out acquisition artefacts and identify the main geological structures. The third one is an example of porosity being integrated geostatistically with a seismic acoustic impedance map. The last example deals with geostatistical time to depth conversion; the interest of performing geostatistical simulations is finally discussed.     

陆态网络重力测网的分形特征与地震监测能力分析

介绍了“中国大陆构造环境监测网络”（简称“陆态网络”）重力观测体系。从分形几何学角度研究了陆态网络重力测网的分形特征,计算得到测网分维数D_f=1.5598、最佳网格化间距r=139km。综合测网分维数与重力时变距、重力异常范围、测网范围等量化指标,对陆态网络重力测网的地震监测能力进行分析,结果表明,测网最佳网格化间距等特征值与M_S6.0地震引发的重力场变化区域的特征异常区半径等统计量化指标接近,具备监测我国除边界及藏北无人区以外大陆区域的M_S6.0及以上地震的能力。     

Reduction of distortion and improvement of efficiency for gridding of scattered gravity and magnetic data

This paper presents a reasonable gridding-parameters extraction method for setting the optimal interpolation nodes in the gridding of scattered observed data. The method can extract optimized gridding parameters based on the distribution of features in raw data. Modeling analysis proves that distortion caused by gridding can be greatly reduced when using such parameters. We also present some improved technical measures that use human-machine interaction and multi-thread parallel technology to solve inadequacies in traditional gridding software. On the basis of these methods, we have developed software that can be used to grid scattered data using a graphic interface. Finally, a comparison of different gridding parameters on field magnetic data from Ji Lin Province, North China demonstrates the superiority of the proposed method in eliminating the distortions and enhancing gridding efficiency.     

应用雷达拼图数据估测降水试验

利用2008年6月5～7日广东省新一代天气雷达网(广州、梅州、韶关、阳江、深圳及汕头6部雷达)的雷达原始体扫资料及自动雨量站资料,对广东省3 km高度上的雷达网系统观测值差异进行了分析,发现广州雷达的观测值比周围雷达偏高1～3 dBz,梅州雷达比其周围雷达的观测值偏低1～2 dBz。用Z-R关系和最优插值校准法分别进行6 min和1 h估测降水,并用面雨量偏差和均方根误差对该次降水估测试验做了简单评估。结果显示:6 min定量估测降水,两种方法都会低估,面雨量越大,估测效果越好;对于1 h定量估测降水,各种估测方法都有偏高或偏低情况,但普遍偏高,其中用先累加再用最优插值校准法校准的雷达-雨量计联合估测方法效果最好;短时降水估测可以很好地反映降水过程变化,而长时间降水估测可以较准确估测降水大小。     

Gravity Delineation of a Buried Valley in Quartzite

A gravity survey was conducted in a 94-square-mile area of northeastern Hanson County, South Dakota. The 340 measured gravity values, together with test-hole data, were used to approximately delineate a buried valley eroded into the Sioux Quartzite of Precambrian age. This valley contains, in places, an aquifer composed of quartzose sand of pre-Cretaceous age derived from the Sioux Quartzite. The bottom of the valley is approximately 450–600 ft below land surface. Simple Bouguer values were determined from measured gravity data, referenced to a local base station, and interpolated to a 0.5- by 0.5-mile grid. The interpolated simple Bouguer values of residual gravity were determined using a five-ring, inverse-weighted filtering method. The second derivative of the interpolated gravity values, as well as their downward continuation, did not delineate the buried valley as well. Subsequent drilling of nine test holes showed that the gravity method can be used for approximately delineating subsurface features     

考虑航磁水平梯度变化的ΔT网格化方法研究

实测航磁横向水平梯度反映垂直于测线的磁场梯度,比传统航磁ΔT数据包含测线之间更多的磁场信息.针对航磁数据网格化问题,采用Hardwick提出的方法,利用航磁水平梯度与ΔT数据构建拟测线,并结合Akima插值法,开展了双方向测线型ΔT网格化方法研究,最终实现了考虑航磁水平梯度变化的ΔT网格化;针对网格化结果中的虚假异常采取了有效滤波方法.通过理论模型数据和实际数据网格化处理,表明该方法可以突出航磁测线之间的异常细节、更清晰地反映线性构造或磁性体走向,提高了网格化的精度和分辨率.     

青藏高原东部玛多-沙马地区的重力场与深部构造

根据青藏高原东部玛多-沙马（下察隅）重力剖面的重力数据资料，对该地区的重力场和深部地壳构造特征作了分析研究，提出青藏高原东部的布格重力异常是高原边缘高，内部低，地壳厚度是边缘薄，内部厚，平均地壳厚度为60km左右，在察隅-沙马地区，为负均衡异常区，因此，该地区是属于地壳上升的地区，此项结果，填补了察隅-沙马地区的均衡重力异常的空白。     

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

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

Crustal structure in Dabieshan UHP metamorphic belt and its tectonic implication

IntroductionGeologistsfirstlyfoundcoesite-bearingecologitesattheendof1980'sandthenthemicrodiamond(Xu,elal,1992)inDabieshanarea.Theultra-highpressure(UHP)metamorphismandthegeodynamicprocessesofDabieorogenhaveattractedmanygeoscientists(Wang,etal,1989,O...