A gradient‐ratio method for the semi‐automatic interpretation of gravity map data sets

The application of semi‐automatic interpretation techniques to potential field data can be of significant assistance to a geophysicist. This paper generalizes the magnetic vertical contact model tilt‐depth method to gravity data using a vertical cylinder and buried sphere models. The method computes the ratio of the vertical to the total horizontal derivative of data and then identifies circular contours within it. Given the radius of the contour and the contour value itself, the depth to the source can be determined. The method is applied both to synthetic and gravity data from South Africa. The Matlab source code can be obtained from the author upon request.     

川滇地区次级地块及其主要断裂带现今活动研究

以位于川滇区域的陆态网和四川省CORS站速度场为约束,结合研究区已知的断裂分布情况,通过使用欧拉矢量计算方法和反欧拉矢量计算方法求出每个次级块体的欧拉矢量,进而对川滇区域次级块体作进一步的调整。同时对该区域进行1°×1°的网格插值,并利用各次级地块的欧拉矢量反演出每个次级块体的长期运动方向及其速度;对相邻块体的边界处,结合各自块体的欧拉矢量求出同一点在2个块体的速度并作差,进而研究分析各主要断裂带的现今活动速率。     

Multiscale anisotropic diffusion for ringing artifact suppression in geophysical deconvolution data

Ringing artifact degradations always appear in the deconvolution of geophysical data. To address this problem, we propose a postprocessing approach to suppress ringing artifacts that uses a novel anisotropic diffusion based on a stationary wavelet transform (SWT) algorithm. In this paper, we discuss the ringing artifact suppression problem and analyze the characteristics of the deconvolution ringing artifact. The deconvolution data containing ringing artifacts are decomposed into different SWT subbands for analysis, and a new multiscale adaptive anisotropic filter is developed to suppress these degradations. Finally, we demonstrate the performance of the proposed method and describe the experiments in detail.     

A Euler number-based topological computation model for land parcel database updating

Intersection relations are important topological considerations in database update processes. The differentiation and identification of non-empty intersection relations between new updates and existing objects is one of the first steps in the automatic incremental update process for a land parcel database. The basic non-empty intersection relations are meet, overlap, cover, equal and inside, but these basic relationships cannot reflect the complex and detailed non-empty relations between a new update and the existing objects. It is therefore necessary to refine the basic non-empty topological relations to support and trigger the relevant update operations. Such relations have been refined by several researchers using topological invariants (e.g., dimension, type and sequence) to represent the intersection components. However, the intersection components often include only points and lines, and the refined types of 2-dimensional intersection components that occur between land parcels have not been defined. This study examines the refinement of non-empty relations among 2-dimensional land parcels and proposes a computation model. In this model, an entire spatial object is directly used as the operand, and two set operations (i.e., intersection (∩) and difference (\)) are applied to form the basic topological computation model. The Euler number is introduced to refine the relations with a single 2-dimensional intersection (i.e., cover, inside and overlap) and to distinguish the refined types of 2-dimensional intersection components for the relations with multiple intersections. In this study, the cover and overlap relations with single intersections between regions are refined into seven cases, and nine basic types of 2-dimensional intersection components are distinguished. A composite computation model is formed with both Euler number values and dimensional differences. In this model, the topological relations with single intersections are differentiated by the value of the dimension and the Euler number of the resulting set of the whole-object intersection and differences, whereas the relations with multiple intersections are discriminated by the value of the resulting set at a coarse level and are further differentiated by the type and sequence of the whole-object intersection component in a hierarchical manner. Based on the refined topological relations, an improved method for automatic and incremental updating of the land parcel database is presented. The effectiveness of the models and algorithms was verified by the incremental update of a land cover database. The results of this study represent a new avenue for automatic spatial data handling in incremental update processes.     

求解模型对欧拉矢量的影响

利用GNSS水平运动场求解欧拉矢量时,待估参数的个数会对欧拉矢量结果产生影响,针对该问题本文分别从理论推导和具体算例角度进行了分析.首先从理论上推导出不同的求解模型得到的欧拉矢量的差异;然后以2004 ～2007年中国大陆GNSS水平运动场为基础,选用两种常用的求解模型(块体整体旋转模型和块体的整体旋转与均匀应变模型),讨论了求解模型对欧拉矢量及后续研究的影响.结果表明,两种模型下块体整体旋转的差异最大可达2.60mm·a-1,是不能忽略的.因此认为选用不同的块体运动模型会得到不同的地壳水平运动图像,在地壳水平运动分析中对此需加以重视.     

Fast mapping of magnetic basement depth,structure and nature using aeromagnetic and gravity data: combined methods and their application in the Paris Basin

Assessment of deep buried basin/basement relationships using geophysical data is a challenge for the energy and mining industries as well as for geothermal or CO₂ storage purposes. In deep environments, few methods can provide geological information; magnetic and gravity data remain among the most informative and cost‐effective methods. Here, in order to derive fast first‐order information on the basement/basin interface, we propose a combination of existing and original approaches devoted to potential field data analysis. Namely, we investigate the geometry (i.e., depth and structure) and the nature of a deep buried basement through a case study SW of the Paris Basin. Joint processing of new high‐resolution magnetic data and up‐to‐date gravity data provides an updated overview of the deep basin. First, the main structures of the magnetic basement are highlighted using Euler deconvolution and are interpreted in a structural sketch map. The new high‐resolution aeromagnetic map actually offers a continuous view of regional basement structures and reveals poorly known and complex deformation at the junction between major domains of the Variscan collision belt. Second, Werner deconvolution and an ad hoc post‐processing analysis allow the extraction of a set of magnetic sources at (or close to) the basin/basement interface. Interpolation of these sources together with the magnetic structural sketch provides a Werner magnetic basement map displaying realistic 3D patterns and basement depths consistent with data available in deep petroleum boreholes. The last step of processing was designed as a way to quickly combine gravity and magnetic information and to simply visualize first‐order petrophysical patterns of the basement lithology. This is achieved through unsupervised classification of suitably selected gravity and magnetic maps and, as compared to previous work, provides a realistic and updated overview of the cartographic distribution of density/magnetization of basement rocks. Altogether, the three steps of processing proposed in this paper quickly provide relevant information on a deep buried basement in terms of structure, geometry and nature (through petrophysics). Notwithstanding, limitations of the proposed procedure are raised: in the case of the Paris Basin for instance, this study does not provide proper information on Pre‐Mesozoic basins, some of which have been sampled in deep boreholes.     

应用改进欧拉算法解算磁性目标空间位置参数

指出目前磁性目标探测中常规欧拉方法存在的局限性,提出在欧拉窗口内视地磁场和其他目标干扰异常联合影响为线性变化,将构造指数作为动态变化参数代入欧拉方程一并求解,并利用多元线性回归方法解决了欧拉方程非线性问题.在质量控制方案中,根据构造指数和目标深度变化规律,给出离散欧拉解滤波措施;采用层次聚类分析方法对滤波后欧拉解进行分...     

欧拉梁动力反应初边值问题的微分求积解

动荷载作用下欧拉梁动响应的计算是一个初边值问题,通常很难得到解析解,传统数值方法一般是把空间和时间分别离散进行求解,计算相对复杂,效率也不高.针对分布动荷载作用下欧拉梁的振动偏微分方程,采用传统微分求积法,在空间和时间上同时进行离散;对于所有非0阶的初/边值条件,采用嵌入法在权系数计算中予以考虑.算例的数值结果与精确解的对比证明采用传统微分求积法处理此问题是可行的,而且是高效的.对于实际工程中的其他类似问题,该方法同样适用.     

对古地理再造中古地磁方法的回顾与探讨

古地磁数据的获取和应用是进行古地理再造的重要基础.同时,计算机技术的发展以及GIS技术的推广,极大地推动了古地理再造方法的革新,为古地理再造提供了更加方便、高效的软件辅助平台.本文主要围绕古地磁数据资料、应用方法及软件开发等方面,对其发展历史与研究现状作一详细的回顾与探讨.     

γ测井资料反褶积分层解释法在8411矿床的试验对比

γ测井资料的反褶积分层解释法具有方法简单,精度高,速度快,分辨能力强,用微机进行计算等优点。作者通过用反褶积法对8411矿床γ测井资料的解释,对视特征参数的选择条件,γ强度底数的确定和划分矿层边界的原则,提出了初步看法。将该法的试验结果与传统的平均含量法和矿心分析结果进行了对比,还分析了视特征参数的取值大小对矿层线储量及矿化分层的影响,并讨论了引起线储量相对误差的原因。