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，张家口南部地区有很多个热水异常点，是具有远景意义的地热区。     

高速公路桥头跳车病害的分析与治理方法探讨

根据高速公路桥头跳车病害探地雷达检测结果,结合路基实际工程地质资料及沉降观测值,分析出造成此病害的原因是桥与路之间沉降差异所致,故搭板与垫层和垫层与路基之间产生了脱空,造成桥搭板、垫层碎裂而出现了桥头跳车病害.针对病害根源,采用了压密注浆进行搭板加固及重新摊铺沥青混合料进行调坡处理的治理方法.     

Electrical aerosol spectrometer of Tartu University

The electrical aerosol spectrometer (EAS) of the parallel measuring principle at Tartu University is an efficient instrument for rapid measurement of the unstable size spectrum of aerosol particles. The measuring range from 10 nm to 10 μm is achieved by simultaneously using a pair of differential mobility analyzers with two different particle chargers. The particle spectrum is calculated and measurement errors are estimated in real time by using a least-squares method. Experimental calibration ensures reliability of measurement. The instrument is well suited for continuous monitoring of atmospheric aerosol.     

Parameter identification of inelastic structures under dynamic loads

The generalized model of differential hysteresis contains 13 control parameters with which it can curve‐fit practically any hysteretic trace. Three identification algorithms are developed to estimate the control parameters of hysteresis for different classes of inelastic structures. These algorithms are based upon the simplex, extended Kalman filter, and generalized reduced gradient methods. Novel techniques such as global search and internal constraints are incorporated to facilitate convergence and stability. Effectiveness of the devised algorithms is demonstrated through simulations of two inelastic systems with both pinching and degradation characteristics in their hysteretic traces. Owing to very modest computing requirements, these identification algorithms may become acceptable as a design tool for mapping the hysteretic traces of inelastic structures. Copyright © 2002 John Wiley & Sons, Ltd.     

Mantle structure and rifting processes in the Baikal–Mongolia region: geophysical data and evidence from xenoliths in volcanic rocks

The origin of the Baikal rift zone (BRZ) has been debated between the advocates of passive and active rifting since the 1970s. A re-assessment of the relevant geological and geophysical data from Russian and international literature questions the concept of broad asthenospheric upwelling beneath the rift zone that has been the cornerstone of many “active rifting” models. Results of a large number of early and recent studies favour the role of far-field forces in the opening and development of the BRZ. This study emphasises the data obtained through studies of peridotite and pyroxenite xenoliths brought to the surface by alkali basaltic magmas in southern Siberia and central Mongolia. These xenoliths are direct samples of the upper mantle in the vicinity of the BRZ. Of particular importance are suites of garnet-bearing xenoliths that have been used to construct P–T- composition lithospheric cross-sections in the region for the depth range of 35–80 km.Xenolith studies have shown fundamental differences in the composition and thermal regime between the lithospheric mantle beneath the ancient Siberian platform (sampled by kimberlites) and beneath younger mobile belts south of the platform. The uppermost mantle in southern Siberia and central Mongolia is much hotter at similar levels than the mantle in the Siberian craton and also has significantly higher contents of ‘basaltic’ major elements (Ca, Al, Na) and iron, higher Fe/Si and Fe/Mg. The combination of the moderately high geothermal gradient and the fertile compositions in the off-cratonic mantle appears to be a determining factor controlling differences in sub-Moho seismic velocities relative to the Siberian craton. Chemical and isotopic compositions of the off-cratonic xenoliths indicate small-scale and regional mantle heterogeneities attributed to various partial melting and enrichment events, consistent with long-term evolution in the lithospheric mantle. Age estimates of mantle events based on Os–Sr–Nd isotopic data can be correlated with major regional stages of crustal formation and may indicate long-term crust–mantle coupling. The ratios of ^143/144Nd in many LREE-depleted xenoliths are higher than those in MORB or OIB source regions and are not consistent with a recent origin from asthenospheric mantle.Mantle xenoliths nearest to the rift basins (30–50 km south of southern Lake Baikal) show no unequivocal evidence for strong heating, unusual stress and deformation, solid state flow, magmatic activity or partial melting that could be indicative of an asthenospheric intrusion right below the Moho. Comparisons between xenoliths from older and younger volcanic rocks east of Lake Baikal, together with observations on phase transformations and mineral zoning in individual xenoliths, have indicated recent heating in portions of the lithospheric mantle that may be related to localised magmatic activity or small-scale ascent of deep mantle material. Overall, the petrographic, P–T, chemical and isotopic constraints from mantle xenoliths appear to be consistent with recent geophysical studies, which found no evidence for a large-scale asthenospheric upwarp beneath the rift, and lend support to passive rifting mechanism for the BRZ.     

Heat flow and thermal modeling of the Yinggehai Basin, South China Sea

Geothermal gradients are estimated to vary from 31 to 43 °C/km in the Yinggehai Basin based on 99 temperature data sets compiled from oil well data. Thirty-seven thermal conductivity measurements on core samples were made and the effects of porosity and water saturation were corrected. Thermal conductivities of mudstone and sandstone range from 1.2 to 2.7 W/m K, with a mean of 2.0±0.5 W/m K after approximate correction. Heat flow at six sites in the Yinggehai Basin range from 69 to 86 mW/m², with a mean value of 79±7 mW/m². Thick sediments and high sedimentation rates resulted in a considerable radiogenic contribution, but also depressed the heat flow. Measurements indicate the radiogenic heat production in the sediment is 1.28 μW/m³, which contributes 20% to the surface heat flow. After subtracting radiogenic heat contribution of the sediment, and sedimentation correction, the average basal heat flow from basement is about 86 mW/m².Three stages of extension are recognized in the subsidence history, and a kinematic model is used to study the thermal evolution of the basin since the Cenozoic era. Model results show that the peak value of basal heat flow was getting higher and higher through the Cenozoic. The maximum basal heat flow increased from 65 mW/m² in the first stage to 75 mW/m² in the second stage, and then 90 mW/m² in the third stage. The present temperature field of the lithosphere of the Yinggehai Basin, which is still transient, is the result of the multistage extension, but was primarily associated with the Pliocene extension.     

三维地震数据离散光滑插值的共轭梯度法

针对三维地震数据插值,提出采用Laplacian算子进行光滑约束的插值方法,并借鉴Mallet研究的离散光滑插值思路,采用预条件共轭梯度法,直接生成网格节点上的值,从而回避寻求满足插值方程的函数. 为了实现其中Laplacian算子的快速求逆,文中引入Claerbout螺旋坐标系谱因式分解理论. 在螺旋坐标系下,Laplacian算子的表示矩阵具有Toeplitz结构,其快速求逆可由谱法LU分解实现. 基于二维离散光滑插值,文中还给出共轭梯度法与NMO相结合的沿时间切片逐层处理的离散光滑插值流程. 最后,应用该方法对模型数据和实际三维地震数据进行了处理.     

复杂介质结构中折射界面的哈格多恩原理波前成像

在城市活断层探测中 ,浅层结构常常表现为强烈的非均匀性 ,界面横向强烈起伏 ,层内速度变化较大 ,传统的基于平界面均匀层模型的折射资料处理方法不能适用。研究开发能应用于复杂介质结构中折射资料处理的方法就显得十分必要。文中基于惠更斯原理 ,用波前扩张法对波场作正演计算 ,根据哈格多恩折射波前成像原理 ,在lecomte算法和Hole有限差分计算程序的基础上 ,开发出 1种复杂介质结构中折射资料的处理方法与软件 ,并用此方法处理了福州城市活断层折射探测试验中在义序完成的 2条折射剖面资料。结果表明 :探测区浅层为 3层结构 ,分别为盖层、强风化层和基岩。基岩顶界面的埋深约为 5 8～ 5 2m ,盖层P波速度变化较大