三维波动方程有限差分正演方法

三维地震资料的处理和解释都需要有效的三维正演模型予以验证．本文提出一种在ｘ－ｔ域实现快速、高精度的有限差分正演方法，采用了独特的“平行四边形网格”，并用Ｐ－Ｒ交替差分格式使三维波动方程可以局部地分裂成二维求解方程，从而有效地减少运算量．     

黑水─泉州地学断面的重磁解释

讨论华南黑水─泉州地学断面的重磁解释。在解释中除应用常规的２．５维重磁异常反演外，还采用等密度线算法构制梯度层密度模型，用于分辨地壳密度细结构，以及计算自由地幔而深度用于分析上地幔密度的横向不均匀性。重磁模拟结果揭示了１０ｋｍ深度内的上地壳构造以及地壳与上地幔密度的分布。提出华南造山带以低密度上地幔为特征，它可能与上地幔的改造有关；四川盆地具有较高的上地幔密度，为未经改造的原始地幔。扬子克拉通与华南造山带的分界线与上地幔密度的分界线一致。根据布格异常、地表岩石密度和地形资料的综合分析，圈定出反映内生成矿作用深部标志的密度倒转区，可作为进一步找矿的远景区。     

Simple stochastic modelling of the eruption history of a basaltic volcano: Nyamuragira,Zaire

Three simple models of the behaviour of a series of basaltic eruptions have been tested against the eruptive history of Nyamuragira. The data set contains the repose periods and the volumes of lava emitted in 22 eruptions since 1901. Model 1 is fully stochastic and eruptions of any volume with random repose intervals are possible. Models 2 and 3 are constrained by deterministic limits on the maximum capacity of the magma reservoir and on the lowest drainage level of the reservoir respectively. The method of testing these models involves (1) seeking change points in the time series to determine regimes of uniform magma supply rate, and (2) applying linear regression to these regimes, which for models 2 and 3 are the determinsstic limits to those models. Two change points in the time series for Nyamuragira, in 1958 and 1980, were determined using a Kolmogorov-Smirnov technique. The latter change involved an increase in the magma supply rate by a factor of 2.5, from 0.55 to 1.37 m³s^-1. Model 2 provides the best fit to the behavior of Nyamuragira with the ratio of variation explained by the model to total variation. R², being greater than 0.9 for all three regimes. This fit can be interpreted to mean that there is a determinstic limit to the elastic strength of the magma reservoir 4–8 km below the summit of the volcano.     

Application of numerical wave-propagation techniques to study local soil effects: The case of Benevento (Italy)

Many of the numerical techniques used for seismic zonation studies treat one-dimensional structural models and/or the incidence of plane polarized body waves. These techniques are often not adequate for laterally heterogeneous structures and for sources that are not located beneath the site of interest. In such cases a more rigorous treatment of the combined effects of the source, the path and the site response is needed. This can be accomplished with a hybrid approach combining modal summation and the finite-difference technique. To demonstrate the differences between these techniques, the ground motion in the city of Benevento (Italy) is modelled. We first compare the results obtained with one-and two-dimensional structural models for vertical incidence of plane polarized body waves. These results are then compared with those obtained with the hybrid approach for two-dimensional structural models.The comparisons have allowed us to find important differences in the response obtained with the different modelling techniques. For the same site, these differences consist of strong variations in amplitude and in the shape of the spectral amplifications. For a seismic source which is not located beneath the site, vertical incidence of waves significantly overestimates the local hazard in a laterally homogeneous structure. For a laterally heterogeneous area, we can conclude that one-dimensional modelling fails to estimate the seismic hazard, whereas for a seismic source which is not located beneath the site of interest, two-dimensional modelling with vertical incidence of plane polarized body waves may not allow reliable estimates to be made of the frequency bands at which amplifications occur. The results obtained for two-dimensional structural models are used for a zonation of the city of Benevento.     

Geomagnetic depth sounding over the Singhbhum and the surrounding regions of eastern India

Magnetovariational studies have been carried out in Singhbhum and surrounding regions during 1987 and 1989. Three deep-seated linear conductors have been identified. One of them is located to the north of Ranchi, Bokaro and Purulia extending in E-W direction coinciding with high heat flow region and Gondwana sediments. The trend of anomaly at Ranchi and Purulia at longer periods suggests a conductivity anomaly due to the mafic and ultramafic intrusions, considered to be responsible for the uplift of Chhotanagpur plateau. The second conductor is associated with the basin margin fault that separates the Singhbhum craton and Chhotanagpur plateau from the West Bengal basin. This conductive zone appears to extend further south and join the high heat flow region of Attri-Tarabalo. This conductor could be isolated only after eliminating the coast effect from the observed induction vectors. The third conductive zone follows the trend of Mahanadi valley located south of the Sukinda thrust. Conductive anomaly associated with the Sukinda and Singhbhum thrust zones could not be resolved due to the interference from neighbouring conductive structures. These two thrusts may not be very deep-seated structures. The Singhbhum granite batholith is found to be highly resistive and seems to extend to greater depths.     

Semi-empirical method for calculation of maximum subsidence in coal mines

Summary Most empirical methods of subsidence prediction are based on the prior knowledge of the maximum possible subsidence, S_max, for a particular coalfield. The subsidence profile and other associated parameters along a desired line may be derived, in most cases, as a function of S _max. These methods are site specific and are not based upon the rational concepts of mechanics. Purely mechanistic methods, on the other hand, could not find wide application because of their limitations in representing the complex behaviour of a rockmass.Therefore a semi-empirical method of calculation of S _max or S, which is an improvement over the empirical method reported elsewhere (Bahuguna et al., 1991a), has been derived from combining the mechanistic and empirical approach. The effects of various parameters are studied by numerical modelling. These qualitative results are then used to influence the development of a more general semiempirical method. The method has been tested for 125 coal mine workings in India and 22 mines of the North Appalachian basin. The method is easily adaptable to other countries.     

长江三峡工程构造应力场数值模拟

本文分别进行了长江三峡工程区域1/100万现今构造应力场有限元数值模拟和长江三峡工程库首区1/50万燕山期、新构造期和现今构造应力场演化规律模拟和分析。区域模拟结果表明:坝区不存在发生大于5级地震的构造应力场背景,库首区模拟结果表明:坝区不存在发生大于3级地震的构造应力场背景。