福州地区的Q结构和它的横向不均匀性

本文利用人工爆破地震资料,在频率域中,用波组频谱比方法计算了福州地区的Q结构,结果发现:该区中地壳低速层Q值明显偏低,仅为40——80,且存在横向变化,地壳上部和下部Q均为200左右;福州盆地与永泰——梧桐地段地壳平均Q也明显偏低,且福州盆地中部Q仅120—150,为全区最低值,低Q区与温泉分布和断裂位置有良好的对应关系。     

非渠化海相碳酸盐岩滩储集地质体的构型研究

１９８５年Ｍｉａｌｌ创立于河流沉积体系的构型单元分析方法可通过旋回性分级特征运用于非渠化海相碳酸盐岩滩沉积体中，并且在构型界面和构型 与河流沉积体的划分可进行对比，所获得的非均质性储层地质模型表明，沉积体最小的均一性单元是由二级界面所限定的较小规模中层型，该构型规模以上的非均质性特征受构型界面级次和构型单元多样性的控制。     

中国东部新生代玄武岩的地球化学（Ⅱ）Sr、Nd、Ce同位素组成

对有代表性的中国东部新生代玄武岩的Ｓｒ、Ｎｄ和Ｃｅ同位素组成进行了系统的研究，结合系统的微量元素研究结果，进一步阐明其地幔源区同位素组成的不均一性；利用Ｃｅ、Ｎｄ同位素地球化学给出的信息，探讨富集Ⅰ型（ＥＭⅠ）地幔端元的成因。     

利用地震层析成象技术确定横向非均匀介质中的速度分布和界面位置

本文提出一种横向非均匀介质中的地震层析成象方法，由地震走时资料确定介质的速度分布和界面位置。计算中采用了参数分离、变阻尼、奇异值分解等技术，方程的解较稳定，且可得到解的分办。除反射波外，该方法还可同时利用分析波的资料。因此，适宜于折射及宽角度反射地震测深剖面的资料处理。数字试验结果表明，在模型内部，该方法能得到可靠的反演结果。在边界附近效果较差，受初始值影响较大。     

Micro-scale modelling in the study of plume evolution in heterogeneous media

The migration of contaminants in heterogeneous aquifers involves dispersive processes that act at different scales. The interaction of these processes as a plume evolves can be studied by micro-scale modelling whereby two scales, a local- or micro-scale and an aquifer- or macro-scale, are covered simultaneously. Local-scale dispersive processes are represented through the local dispersion coefficient in the transport equation, while large-scale dispersion due to heterogeneities is represented through the resolution of the flow field and the diffusive exchange between streamtubes. The micro-scale model provides both the high degree of resolution compatible with local-scale processes, and the extent required for the approach to asymptotic conditions, using grids of up to a million nodal points. The model is based on the dual potential-streamfunction formulation for flow, and the transport problem is formulated in a natural coordinate system provided by the flownet. Simulations can be used to verify stochastic theories of dispersion, without the restrictive assumptions inherent in the theory. For the two-dimensional case, results indicate convergence of the effective dispersivity to the theoretical macrodispersivity value. Convergence takes place within a travel distance of about 50 correlation lengths of the hydraulic conductivity field. However, the approach taken to asymptotic conditions, as well as the macrodispersivity value, may differ for different realizations of the same medium. The influence of early-time events such as plume splitting on the asymptotic convergence remains to be investigated.     

Coda waves: A review

The analysis and interpretation of coda waves have received increasing attention since the early seventies. In the past few years interest in this subject has spread worldwide, and the study of high-frequency seismic coda waves has become a very important seismological topic. As a conclusion of the studies accomplished in this time, coda waves are considered the result of scattering processes caused by heterogeneities acting on seismic waves.P andS waves play a particularly important role in this interaction. The process introduces an attenuation which, added to the intrinsic absorption, gives the observed apparent attenuation. Therefore, coda waves constitute a thumbprint left by the heterogeneities on the seismograms. Coda waves offer decisive information about the mechanism of how scattering and attenuation take place. This review describes coda waves in detail, and summarizes the work done in this subject to 1986. The relation between coda waves and attenuation in the context of research on seismic scattering problems is stressed. Particular attention has been given to the application of coda waves to estimate source and medium parameters. The state-of-the-art of the temporal variations of coda wave shape, and the possible use of these variations as an earthquake precursor also are presented. Care has been taken to introduce the statistical models used to deal with the heterogeneities responsible for scattering.     

Three-dimensional elastic wave scattering and diffraction due to a rigid cylinder embedded in an elastic medium by a point source

The formal solutions of displacement field to the problem of elastic wave scattering and diffraction due to an infinitely long rigid cylinder embedded in an infinite elastic medium by an impulsive point source have been obtained in the integral form. The integrals for the reflected and the diffracted waves both in the shadow zone and in the illuminated zone are evaluated asymptotically for the early time motion by the Reisdue-Cagniard method and the Saddle-point-Cagniard method.Numerical results of the diffractedP, S andPS waves at a fixed circum-distance from the surface of the rigid cylinder show noticeably that (1) the energy partition for the diffractedS wave is small in comparison with that for the diffractedP wave, (2) the wave form of the diffractedS wave is broader and more diffused than that of the diffractedP wave, (3) the direction of the radial motions of the diffractedP andS waves varies as a function of the observational point, and (4) the energy partition for the diffractedP wave is much smaller than that for the direct or the reflectedP waves.This paper has been presented at the 46th Annual International Meeting of Society of Exploration Geophysicists in Houston, Texas, Oct. 28, 1976.     

A Combined Three-Dimensional Geological-Geostatistical-Numerical Model of Underground Excavations in Rock

Summary. This paper exploits geological and borehole geotechnical data obtained in the exploratory phase of a tunneling project to investigate in a first place if the kriging interpolation scheme may effectively reproduce the spatial variability of rock mass quality (Rock Mass Rating, RMR) in the vicinity of tunnels. For this purpose a quick solver in Fortran has been developed that performs variography analysis of 3D spatial data, fast kriging estimations of RMR between borehole sampling locations at the centroids of the elements of the numerical model, and model validation. For the purpose of an integrated underground excavation design, a step further is made by incorporating into the 3D mechanical numerical model of the rock mass, the three-dimensional (3D) solid geological model, thus coupling the geology with the ground (geotechnical) model (i.e. each element of the numerical model is assigned a geological material). The mechanical properties of each finite difference cell (or Representative Elementary Volume) of the ground model were then prescribed according to its geological type, the spatial heterogeneity of the rock mass expressed quantitatively with the kriging model, and the upscaling calculations of the mechanical properties of the intact rocks determined in the laboratory, based on the size-effect (strength dependence on size) and Damage Theory. Furthermore, a preliminary numerical simulation of the advance of unsupported tunnels in the model of the heterogeneous rock mass was performed for illustration purposes.