地震勘探技术在深部找煤中的应用

成矿理论与勘探实践证明,东部一些老矿由于勘探深度、范围的限制以及认识上的局限,在矿区深部、周边仍有一些矿未能发现和查明,其中相当一部分仍具有找矿潜力;而寻找西部新的煤炭资源是国家“稳定东部发展西部”战略需要。针对中国东西部不同的地质条件和勘探目标,其地震勘探野外施工方法及资料处理原则有所差异;而资料解释在东部以查明地层赋存形态、煤层赋存范围、构造发育特征为主,西部则以寻找煤系地层波阻抗差异明显、连续性好的反射波组为解释重点。地震勘探技术在中国东西部的二个典型地震时间剖面佐证了其在深部找矿的应用效果。     

Multi-Scale Texture Modeling

Patterns of crystallographic preferred orientation are referred to as texture. The specific subject of texture analysis is the experimental determination and interpretation of the statistical distribution of orientations of crystals within a specimen of polycrystalline material, which could be metals or rocks. The objective is to relate an observed pattern of preferred orientation to its generating processes and vice versa. In geosciences, texture of minerals in rocks is used to infer constraints on their tectono-metamorphic history. Since most physical properties of crystals, such as elastic moduli, the coefficients of thermal expansion, or chemical resistance to etching depends on crystal symmetry and orientation, the presence of texture imparts directional properties to the polycrystalline material. A major issue of mathematical texture analysis is the resolution of the inverse problem to determine a reasonable orientation density function on SO(3) from measured pole intensities on , which relates to the inverse of the totally geodesic Radon transform. This communication introduces a wavelet approach into mathematical texture analysis. Wavelets on the two-dimensional sphere and on the rotational group SO(3) are discussed, and an algorithms for a wavelet decomposition on both domains following the ideas of Ta-Hsin Li is given. The relationship of these wavelets on both domains with respect to the totally geodesic Radon transform is investigated. In particular, it is shown that the Radon transform of these wavelets on SO(3) are again wavelets on . A novel algorithm for the inversion of experimental pole intensities to an orientation density function based on this relationship is developed.     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

山区煤田地震勘探中静校正存在的问题及其改进方法

在山区或地表复杂地区进行地震勘探,采用基于地表一致性假设的静校正将会严重影响勘探效果,该影响主要源于地表一致性假设存在着不合理因素。如较高的低速带、巨厚的低速带、基岩裸露、地形起伏较大等。为分析一致性假设静校正偏差产生的原因及大小．构建一地形起伏、基岩出露的复杂模型,通过正演其射线路径,对比其时距曲线与理论时距曲线的差异,以及二者静校正量误差大小。模型分析证实该差异与偏移距、地震波穿透深度及基准面高程之间存在直接的联系,据此提出了改进方法,如浮动基准面校正及分块静校正等。理论模型和实际地震资料试算表明,使用改进的方法可有效改善地震时间刮面同相轴聚焦效果及连续性。     

Bayesian estimation of seismic hazard for two sites in Switzerland

Seismic hazard analysis is based on data and models, which both are imprecise and uncertain. Especially the interpretation of historical information into earthquake parameters, e.g. earthquake size and location, yields ambiguous and imprecise data. Models based on probability distributions have been developed in order to quantify and represent these uncertainties. Nevertheless, the majority of the procedures applied in seismic hazard assessment do not take into account these uncertainties, nor do they show the variance of the results. Therefore, a procedure based on Bayesian statistics was developed to estimate return periods for different ground motion intensities (MSK scale).Bayesian techniques provide a mathematical model to estimate the distribution of random variables in presence of uncertainties. The developed method estimates the probability distribution of the number of occurrences in a Poisson process described by the parameter . The input data are the historical occurrences of intensities for a particular site, represented by a discrete probability distribution for each earthquake. The calculation of these historical occurrences requires a careful preparation of all input parameters, i.e. a modelling of their uncertainties. The obtained results show that the variance of the recurrence rate is smaller in regions with higher seismic activity than in less active regions. It can also be demonstrated that long return periods cannot be estimated with confidence, because the time period of observation is too short. This indicates that the long return periods obtained by seismic source methods only reflects the delineated seismic sources and the chosen earthquake size distribution law.     

Some comparisons between mining-induced and laboratory earthquakes

Although laboratory stick-slip friction experiments have long been regarded as analogs to natural crustal earthquakes, the potential use of laboratory results for understanding the earthquake source mechanism has not been fully exploited because of essential difficulties in relating seismographic data to measurements made in the controlled laboratory environment. Mining-induced earthquakes, however, provide a means of calibrating the seismic data in terms of laboratory results because, in contrast to natural earthquakes, the causative forces as well as the hypocentral conditions are known. A comparison of stick-slip friction events in a large granite sample with mining-induced earthquakes in South Africa and Canada indicates both similarities and differences between the two phenomena. The physics of unstable fault slip appears to be largely the same for both types of events. For example, both laboratory and mining-induced earthquakes have very low seismic efficiencies where _a is the apparent stress and is the average stress acting on the fault plane to cause slip; nearly all of the energy released by faulting is consumed in overcoming friction. In more detail, the mining-induced earthquakes differ from the laboratory events in the behavior of as a function of seismic momentM ₀. Whereas for the laboratory events 0.06 independent ofM ₀, depends quite strongly onM ₀ for each set of induced earthquakes, with 0.06 serving, apparently, as an upper bound. It seems most likely that this observed scaling difference is due to variations in slip distribution over the fault plane. In the laboratory, a stick-slip event entails homogeneous slip over a fault of fixed area. For each set of induced earthquakes, the fault area appears to be approximately fixed but the slip is inhomogeneous due presumably to barriers (zones of no slip) distributed over the fault plane; at constant , larger events correspond to larger _a as a consequence of fewer barriers to slip. If the inequality _a / 0.06 has general validity, then measurements of _a =µE _a /M ₀, where is the modulus of rigidity andE _a is the seismically-radiated energy, can be used to infer the absolute level of deviatoric stress at the hypocenter.     

地震综合效应场函数及其在地震预报中的应用──（一）地震综合效应场函数及其计算方法

在系统地分析了目前各种测震学地震预报方法科学思路的基础上,认为测震学地震预报方法基本上可以分为两大类。一类是以已经发生的一些地震作为未来可能发生的地震的“因”,即由于已经发生的地震对区域应力场的影响,导致未来发生较强地震。这一类包括的预报方法较多,如空区、条带、ｂ值、地震迁移、相关地震等等及其由此衍生出来的各种方法。另一类是把已经发生的一些地震作为区域应力场增强的“果”,即已经发生的地震是区域应力场增强过程中的一种反映,而未来地震不一定是已经发生的地震所导致的结果。这一类包括“地震窗口”、小震群活动等方法。针对第一类方法,各种预报方法都是力图从地震三要素中提取未来地震的信息,而具体作法又都是利用地震三要素这个多维空间的某个剖面。为了从地震活动诸要素的多维空间提取综合信息,我们对每个地震加入了破裂面方位,构成了地震第四要素,并依据地震４要素建立了地震综合效应场函数。地震综合效应场函数概括了多种测震学地震预报方法的科学思路和预报经验,从而可以形成测震学的综合预报方法。     

Palaeolimnological studies of the eutrophication of volcanic Lake Albano (Central Italy)

We use palaeolimnological techniques to reconstruct the eutrophication history of a volcanic lake (Lake Albano, central Italy) over the past three centuries. The presence of annual varves down to the bottom of the core (c. 1700 A.D.) indicated the lack of bioturbation and likely long-term meromixis. Sedimentation rates were estimated by varve counts (calcite/diatom couplets), indicating a mean rate of 0.15 cm yr^–1. The reconstruction of eutrophication was traced using past populations of algal and photosynthetic bacteria (through their fossil pigment), and geochemistry, as well as fossil remains of chironomids. Phaeophorbidea and the red carotenoid astaxanthin were used to detect past zooplankton development.The first sign of trophic change related to human activities is datedc. 1870 A.D. From that period onward a sharp increase of authigenic CaCO₃, nitrogen, N:P ratio, and dinoxanthin, a characteristic carotenoid of Chrysophyceae and Dinophyceae, is observed.Chironomid analyses showed the near absence of a deep water fauna throughout the core length. The populations of chironomid larvae are restricted to oxygenated littoral zones. In fact, the few fossil remains found are primarily of littoral origin, representing shallow water midges that were transported to profundal waters. The reduction of total chironomid in the uppermost layers of the core is to be related to human land uses.     

Palaeoclimatic studies in South Shetland Islands,Antarctica, based on numerous stratigraphic variables in lake sediments

The hitherto longest found lake sediment sequence on Byers Peninsula, Livingston Island, South Shetland Islands, was analysed with respect to lithology, chronology, diatoms, Pediastrum, pollen and spores, mosses, mineralogy, and sediment chemistry. During the ca. 5000 year long development the sediments were influenced by frequent tephra fall-outs. This volcanic impact played a major role in the lake's history during two periods, 4700–4600 and 2800–2500 BP, but was of importance during the lake's entire history with considerable influence on many of the palaeoenvironmentally significant indicators. The large and complex data set was analysed and zonated with different types of multivariate analysis. This resulted in a subdivision of the sequence into 8 time periods and 21 variables. Redundancy analysis (RDA) of this data set, both without and with the tephra periods, and with 4–6 of the variables as explanatory environmental variables, reveal that climatic/environmental signals are detectable. The palaeoclimatic picture that emerged out of the tephra noise suggests that the first 100 years were characterized by mild, humid conditions. This was followed by a less mild and humid climate until ca. 4000 BP when a gradual warming seems to have started, coupled with increased humidity. These mild and humid conditions seem to have reached an optimum slightly after 3000 BP. At ca. 2500 BP a distinct climatic deterioration occurred with colder and drier conditions and long seasons with ice cover. This arid, cold phase probably reached its optimum conditions at ca. 1500 BP, when slightly warmer conditions might have prevailed for a while. Except for the modern sample with rather mild climate, the last 1400 years seem to have been fairly arid and cold, and the effects of the frequent volcanic activity during this period is only vaguely seen in the records.     

A comparison between 2D and 3D wavelet analysis method for asteroid family dynamical study

As found in previous studies (Bendjoyaet al. 1991) the wavelet analysis is a reliable tool for the determination of asteroid families. A comparison between the 2D and a newly developed 3D wavelet analysis is performed on 14 fictitious families numerically generated the members of which area priori known. It clearly appears that the 3D method brings improvements to the 2D approach. Indeed the 3D analysis appears well suited for all cases even for those corresponding to the very dense Flora region and for exotic families such as filamentar ones. With the 2D method, satisfying results are obtained for the simulations which correspond to 70% of the real asteroid families previously found. The aim of this paper is to show quantitatively the improvement of the 3D method versus the 2D one.