Three-dimensional assembly technologies for the interpretation of gravimetric data

The assembly method of solving a three-dimensional ore-type inverse problem is considered. The method is based on the unity of a grid (finite-element) approximation of a dense medium and a special type of structuring the iteration process of selection which makes it possible to effectively observe all a priori restrictions on the model of field sources. The examples of interpretations of gravity anomalies are presented.     

Tilt梯度在区域重力资料解释中的应用

线性构造识别是区域重力资料解释中的重要一步。本文介绍了用于探测重力异常边界和场源深度的Tilt梯度法;同时,应用Tdr零值线及±π/4曲线距离半值或零值与+π/4(或-π/4)曲线距离对研究区重力异常源的位置和埋深进行了推断,并把推断结果与区内大地电磁测深剖面解释成果进行了对比。结果表明,应用Tilt梯度对区域重力资料进行线性构造识别是可行的。     

Deterministic approach to the problem of reliability in interpretation results of gravimetric data

A method for comparative estimation of reliability in detection of gravity anomaly sources in indicated fragments of studied geological space is given. The possibility of creation of this estimation is connected with disparate ideas on different peculiarities of the structure of studied geological space and noise characteristics in observed fields. The model example and results of practical problem solution are viewed.     

航磁梯度数据解释新方法在迁安铁矿勘探中的应用

为更加清晰和准确地利用磁梯度数据获得地质体的分布特征,笔者采用解析信号总水平导数法和梯度欧拉反褶积技术组合实现磁梯度数据的解释,解析信号总水平导数的极值可有效地圈定出矿体的分布范围,梯度欧拉反褶积法可有效地避免背景异常的干扰,从而准确地给出磁性体的位置和深度信息。理论模型证明解析信号水平导数法和梯度欧拉反褶积法的组合可准确地获得磁性体的水平范围和深度。应用于河北迁安地区磁梯度数据的解释,圈定出5个呈条带状分布的较大型铁矿藏,埋深在200~450 m。     

Relationship between chemical composition of Japanese Island-Arc volcanic rocks and gravimetric data

1670 published chemical analyses from 137 Quaternary eruption centers of Kurile, Japan and Izu—Mariana arcs are collected and K₅₅ and K₆₀ values of individual centers are calculated. Along with the well-established trend of increase of K₅₅ and K₆₀ values across the arc towards the continent, longitudinal change is also recognizable. Interrelation between K₅₅ or K₆₀ value pattern and a map of Bouguer anomaly distribution along the volcanic front is characteristic.To develop one possible simplified model to explain the facts, we assume first that the asthenosphere has a homogeneous thermal state beneath the volcanic zone, second that the thickness of the crust is quite variable areally, and third that the mean density of the crust is constant throughout the area. Thus, the composition of the magma generated along a zone corresponding to equal depths along the Benioff zone should be identical throughout the volcanic arc, and lower values of the Bouguer anomaly should be observed where the crust is thicker. We then assume that magma generated in the asthenosphere ascends into the crust, then cools as fractional crystallization proceeds. Depths of fractional crystallization would be deeper where the crust is thicker. The ratio of the clinopyroxene phase separated from the magma to the plagioclase phase should be larger where the level of fractionation is deeper. Thus, magmatic liquid should be enriched in K where K-poor crystal fractions are subtracted at the deeper level.     

The Angara-Vitim batholite: Data on morphology and sizes based on gravimetric data

Doklady Earth Sciences -     

利用地震资料进行煤层厚度解释预测

利用C++语言基于Windows操作环境开发了适用于煤田地震资料解释的煤层厚度辅助解释系统。该系统采用人工神经网络非线性反演方法对煤层厚度变化进行解释,对煤层厚度和地震属性参数之间的非线性关系给出了定量描述,具有较高的解释精度。理论模型和实际资料试算结果表明:利用该系统进行煤层厚度解释,得到了较好的地质效果。     

突变理论在地震资料解释中的应用

阐述了尖点突变理论,建立了地震信号尖点突变模型并转化为尖点突变的标准形式,在时间域内提取突跳势、突跳时间;在频率域内提取突跳势、突跳次数等参数。运用突变理论对实际地震资料进行构造解释表明:突跳参数曲线上异常的位置与断层的位置有关,为地震资料解释提供了一种新的途径。      

Automatic interpretation of pumping tests data using metaheuristics

Pumping tests data interpretation is of major importance in groundwater engineering. It is traditionally performed in a subjective manner by means of standard type curves. In this paper, an automatic interpretation of time-drawdown data has been proposed based on two algorithms, the real-coded genetic algorithm and differential evolution. The proposed approaches combine metaheuristic algorithms with an appropriate analytical drawdown solution depending upon the nature of the considered aquifer system, leaky or naturally fractured rock aquifers. The standard error of estimate (SEE) was used as a performance criterion to evaluate the discrepancies between predicted and observed drawdown data in different pumping time periods. Both of the proposed metaheuristic algorithms provide accurate aquifer parameters. For all analyzed pumping tests data, the differential evolution yielded the most accurate results with an improvement in SEE values ranged from 0.2 to 50% compared to previously published results, and exhibits speed and robustness. Furthermore, a new full range numerical evaluation of the Hantush well function based on a tanh-sinh quadrature scheme has been proposed. Our results indicate that our method is accurate, since the maximum relative error was found to be equal to 0.00036%, and practical, since it is free from special functions and could be easily incorporated within an optimization technique to analyze transient time-drawdown data.     

起伏地形条件下FPS不同解释方法对比

固定点源激电测深(FPS)是一种有效的激发极化测深勘探方法,但在数据解释中,常规的定量解释法只适应较为简单的地电条件,在起伏地形条件下其解释效果会受到较为明显的影响。针对常规的定量解释方法所存在的不足与缺陷,这里从有限单元法正演模拟入手,采用最小二乘法进行反演计算,实现了在起伏地形条件下的对地电模型进行模拟计算,并将圆弧交汇-相对强度计算结果与最小二乘反演计算结果进行对比。结果表明:最小二乘反演方法能适应起伏地形对模型进行反演计算,其反演效果优于圆弧交汇-相对强度法计算结果。     

Self-organizing maps for geoscientific data analysis: geological interpretation of multidimensional geophysical data

Data interpretation is a common task in geoscientific disciplines. Interpretation difficulties occur especially if the data that have to be interpreted are of arbitrary dimension. This paper describes the application of a statistical method, called self-organizing mapping (SOM), to interpret multidimensional, non-linear, and highly noised geophysical data for purposes of geological prediction. The underlying theory is explained, and the method is applied to a six-dimensional seismic data set. Results of SOM classifications can be represented as two-dimensional images, called feature maps. Feature maps illustrate the complexity and demonstrate interrelations between single features or clusters of the complete feature space. SOM images can be visually described and easily interpreted. The advantage is that the SOM method considers interdependencies between all geophysical features at each instance. An application example of an automated geological interpretation based on the geophysical data is shown.     

重力、航磁资料在花岗岩型铀矿成矿研究中的应用

本文利用重磁场资料对我国南方一些花岗岩体的侵位状态和岩浆动力场进行了分析，并按动力场的强弱对岩体进行了分类。根据重力资料对苗儿山-越城岭、诸广山和贵东花岗岩体的反演计算，讨论了这些岩体的深部分布形态，提出了岩浆流动的3种方式，指出了富大铀矿床、大型铀矿聚集区的赋存部位及其与航磁异常的密切联系。     

Presentation and interpretation of chemical data for igneous rocks

Arguments are made in favor of using variation diagrams to plot analyses of igneous rocks and their derivatives and modeling differentiation processes by least-squares mixing procedures. These methods permit study of magmatic differentiation and related processes in terms ofall of the chemical data available. Data are presented as they are reported by the chemist and specific processes may be modeled and either quantitatively described or rejected as inappropriate or too simple.Examples are given of the differing interpretations that can arise when data are plotted on an AEM ternary vs. the same data on a full set of MgO variation diagrams. Mixing procedures are illustrated with reference to basaltic lavas from the Columbia Plateau.     

平衡剖面技术在地震资料解释中的应用

地震资料解释的正确与否直接关系到储层预测的结果与精度,所以平衡剖面解释作为一种辅助手段,对地震解释的多解性起到很好地制约效果,并为其定量化解释研究上开辟了道路。介绍了平衡剖面技术的基本原理和解释机理;以依奇克里油田阿依库木齐剖面为例,阐述了平衡剖面技术在地震解释中的应用效果。结果显示,该剖面构造复原效果较好。并针对平衡剖面技术在地震资料解释中的应用提出了需注意的问题。     

The application of cluster analysis in geophysical data interpretation

A clustering algorithm that is based on density and is adaptive density-reachable is developed and presented for arbitrary data point distributions in some real-world applications, especially in geophysical data interpretation. Through comparisons of the new algorithm and other algorithms, it is shown that the new algorithm can reduce the dependency of domain knowledge and the sensitivity of abnormal data points, that it can improve the effectiveness of clustering results in which data are distributed in different shapes and different densities, and that it can get a better clustering efficiency. The application of the new clustering algorithm demonstrates that data mining techniques can be used in geophysical data interpretation and can get meaningful and useful results, and that the new clustering algorithm can be used in other real-world applications.     

The geothermal gradient map of Central Tunisia: Comparison with structural,gravimetric and petroleum data

Five hundred and fifty temperature values, initially measured as either bottom-hole temperatures (BHT) or drill-stem tests (DST), from 98 selected petroleum exploration wells form the basis of a geothermal gradient map of central Tunisia. A “global-statistical” method was employed to correct the BHT measurements, using the DST as references. The geothermal gradient ranges from 23° to 49°C/km. Comparison of the geothermal gradient with structural, gravimetric and petroleum data indicates that:

• 1.(1) the general trend of the geothermal gradient curves reflects the main structural directions of the region,
• 2.(2) zones of low and high geothermal gradient are correlated with zones of negative and positive Bouguer anomalies and
• 3.(3) the five most important oil fields of central Tunisia are located near the geothermal gradient curve of 40° C/km.

Such associations could have practical importance in petroleum exploration, but their significance must first be established through further investigation and additional data.     

Comparison and combination of a gravimetric quasigeoid with a levelled GPS data set by statistical analysis

For GPS levelling applications, it is convenient to express the height reference surface in a suitable geodetic reference system. This can be obtained through a set of levelled GPS points. Unfortunately, available data are sparse. A gravimetric geoid is often used to interpolate the height reference surface issued from GPS and levelling. Both surfaces do not coincide exactly with each other. At this point, one must compare two realisations of the geoid, detect outliers, retrieve (if possible) the causes of the discrepancies and finally combine the two kinds of data. The paper presents some practical solutions to these problems.     

准噶尔盆地结构的石油物探综合解析

准噶尔盆地经过50年石油物探工作,取得了丰富资料,对盆地地质结构的认识有了突破性进展.资料显示,二叠系与上下地层之间,存在一个比较明显的密度分界面,对认识和分析盆地重力场特征和异常极有帮助.各时代地层的磁性差异显示,二叠纪及中新生代沉积岩属弱(非)磁性(18×10-6~85×10-6 CGSM),内部磁性差异只能引起小幅度的微弱异常.石炭纪及其以前时代的变质岩、火成岩的磁化率都很大,平均在数百个单位以上,特别是火成岩的磁化率可达数几千至数万单位.这些具有磁性的火成岩,在准噶尔盆地均产生明显局部异常,20世纪70年代以前,石油工作者认为是盆地的磁性基底.与此同时,盆地周围大规模的区域地质调查证实,周边广泛分布的主要是上古生界地槽型沉积建造,而且多是火山碎屑岩沉积,下古生界很少,未发现可靠的前寒武纪露头.地震解释技术从二维到三维、从区块三维到区带三维,不仅从宏观上认识了盆地基底结构、沉积盖层等盆地地质结构,而且能从微观上认识圈闭、储层,在含油气有利构造的发现及石油勘探方面发挥了巨大的作用.     

新场岩体定压力压水试验数据解译方法的对比

以高放废物地质处置北山预选区新场岩体三个深钻孔的定压力压水试验为例,分别从压水阶段数据求解、压力恢复阶段数据求解和特殊试验段参数求解3个方面,对稳定流公式法和非稳定流参数拟合法进行了对比,得到了如下结论：非稳定流参数拟合法较稳定流公式法适用的渗流类型更多,使用数据量更大,结果更为可靠;在压水阶段和压力恢复阶段数据完整时,优先选择压力恢复阶段非稳定流参数拟合结果作为试验段的渗透系数.依据上述筛选原则,得到了新场岩体BS17～BS19三个深钻孔渗透系数的概率分布图,可初步判断新场岩体为低渗透性岩体,适宜于高放废物地质处置.