云南大地热流及地热地质问题

本文报道了云南20多个大地热流测试数据。据热流分布特点并参考前人有关资料,结合地球物理测深资料、岩石热物性与放射性生热率测试数据,对各区地壳深度范围的温度、热流配分组成作了粗略分析,划分出各有特征的若干地质-热流区。云南普遍地热偏高(滇东南小范围除外)的主要原因是新构造运动强烈;怒江以西近代岩浆活动亦强烈,具有典型的板块交汇带的基本地热特点     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.     

A comparative error analysis of manual versus automated methods of data acquisition for algebraic strain estimation

This paper presents a statistical analysis of the algebraic strain estimation algorithm of Shimamoto and Ikeda [Shimamoto, T., Ikeda, Y., 1976. A simple algebraic method for strain estimation from deformed eillipsoidal objects: 1. Basic theory. Tectonophysics 36, 315–337]. It is argued that the error in their strain estimation procedure can be quantified using an expected discrepancy measure. An analysis of this measure demonstrates that the error is inversely proportional to the number of clasts used. The paper also examines the role of measurement error, in particular that incurred under (i) a moment based and (ii) manual data acquisition methods. Detailed analysis of these two acquisition methods shows that in both cases, the effect of measurement error on the expected discrepancy is small relative to the effect of the sample size (number of objects). Given their relative speed advantage, this result favours the use of automated measurement methods even if they incur more measurement error on individual objects. Validation of these results is carried out by means of a simulation study, as well as by reference to studies appearing in previous literature. The results are also applied to obtain an upper bound on the error of strain estimation for various studies published in the literature on strain analysis.     

Improved reliability in the interpretation of geochemical measurements by the quantification of uncertainty from sampling

All geochemical measurements require the taking of field samples, but the uncertainty that this process causes is often ignored when assessing the reliability of the interpretation, of the geochemistry or the health implications. Recently devised methods for the estimation, optimisation and reduction of this uncertainty have been evaluated by their application to the investigation of contaminated land. Uncertainty of measurement caused by primary sampling has been estimated for a range of six different contaminated land site investigations, using an increasingly recognized procedure. These site investigations were selected to reflect a wide range of different sizes, contaminants （organic and metals）, previous land uses （e.g. tin mining, railway sidings and gas works）, intended future use （housing to nature reserves） and routinely applied sampling methods. The results showed that the uncertainty on measurements was substantial, ranging from 25% to 186% of the concentration values at the different sites. Sampling was identified as the dominant source of the uncertainty （〉70% of measurement uncertainty） in most cases. The fitness-for-purpose of the measurements was judged using the optimized contaminated land investigation （OCLI） method. This identifies the optimal level of uncertainty that reduces to overall financial loss caused by the measurement procedures and the misclassification of the contamination, caused by the uncertainty. Generally the uncertainty of the actual measurements made in these different site investigations was found to be sub-optimal, and too large by a factor of approximately two. The uncertainty is usually limited by the sampling, but this can be reduced by increasing the sample mass by a factor of 4 （predicted by sampling theory）. It is concluded that knowing the value of the uncertainty enables the interpretation to be made more reliable, and that sampling is the main factor limiting most investigations. This new approach quantifies this problem for the first time, and allows sampling procedures to be critically evaluated, and modified, to improve the reliability of the geochemical assessment.     

The distribution of radiogenic heat production as a function of depth in the Sierra Nevada Batholith, California

Geochemical analyses and geobarometric determinations have been combined to create a depth vs. radiogenic heat production database for the Sierra Nevada batholith, California. This database shows that mean heat production values first increase, then decrease, with increasing depth. Heat production is 2 μW/m³ within the 3-km-thick volcanic pile at the top of the batholith, below which it increases to an average value of 3.5 μW/m³ at 5.5 km depth, then decreases to 0.5–1 μW/m³ at 15 km depth and remains at these values through the entire crust below 15 km. Below the crust, from depths of 40–125 km, the batholith's root and mantle wedge that coevolved beneath the batholith appears to have an average radiogenic heat production rate of 0.14 μW/m³. This is higher than the rates from most published xenolith studies, but reasonable given the presence of crustal components in the arc root assemblages. The pattern of radiogenic heat production interpreted from the depth vs. heat production database is not consistent with the downward-decreasing exponential distribution predicted from modeling of surface heat flow data. The interpreted distribution predicts a reasonable range of geothermal gradients and shows that essentially all of the present day surface heat flow from the Sierra Nevada could be generated within the 35 km thick crust. This requires a very low heat flux from the mantle, which is consistent with a model of cessation of Sierran magmatism during Laramide flat-slab subduction, followed by conductive cooling of the upper mantle for 70 m.y. The heat production variation with depth is principally due to large variations in uranium and thorium concentration; potassium is less variable in concentration within the Sierran crust, and produces relatively little of the heat in high heat production rocks. Because silica content is relatively constant through the upper 30 km of the Sierran batholith, while U, Th, and K concentrations are highly variable, radiogenic heat production does not vary directly with silica content.     

Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

张家口地下热水资源成因探讨

张家口地热资源丰富，通过对区内地热成矿机制，成矿规律的分析认为，以祁吕系为主的多构造体系复合部位控制了本区地下热水的形成与分布，本区地下热水是大气降水经深循环而成，经计算深循环的深度达为2600m，张家口南部地区有很多个热水异常点，是具有远景意义的地热区。     

Comparison of magnetic parameters with vehicular Br levels in urban roadside soils

Chemical methods are generally chosen to monitor soil pollution but magnetic measurements proved to yield additional information at low cost and less time consumption. In this investigation, the novel use of rapid and non-destructive magnetic measurements to characterize Br levels in Xuzhou (China) urban roadside soils was reported. X-ray fluorescence spectrometry (XRF) was used to quantify Br in the soil samples. Data from 21 roadside soil samples confirm Br contamination, with a mean level of 4.36 mg kg⁻¹ and a range of 2.4–8.7 mg kg⁻¹. These values are higher than that of unpolluted soils in Xuzhou that averages 1.1 mg kg⁻¹. Hierarchical clustering analysis shows Br in Xuzhou roadside soils is mainly from road traffic. Clear correlations between Br levels and simple magnetic parameters [mass specific susceptibility (χ _LF), susceptibility of anhysteretic remanent magnetization (χ _ARM), saturation isothermal remanence (SIRM)] are observed. The present study shows that these three magnetic parameters can be used as a proxy for Br levels in Xuzhou urban roadside soils.