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

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

Long-term seismic noise investigations on Shikotan Island: First results

In 2003–2004, long-term seismic noise observations were launched on Shikotan Island (Lesser Kuril Range) based on the “Shikotan” dormant regional seismic station. The geological and geophysical data on the registration area are reported. Information about the equipment and its technical specifications is given. The precursors to the strongest local earthquakesthat occurred in the Shikotan Island region in January 2005–March 2007 are identified.     

(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

北京及邻区微震活动的初步分析

通过对北京遥测地震台网近年来记录到的北京及邻区地震的震中分布,地震活动频度及能量释放强度的分析,得到本区地震活动在时间分布上具有“聚堆性”。在年发震频度,地震强度和能释放方面均具有双峰值特征,并且具有较好的一致性和同步性,在空间分布上具有条带特征,且形成北东～南西和北西～南东的两条相互交汇的条带。又通过统计分析得到本区发震概率最大的时间段是每年的１０月前后,而地震主要发生在北西～南东带上。     

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.     

Formation of the mid-fifteenth century Kuwae caldera (Vanuatu) by an initial hydroclastic and subsequent ignimbritic eruption

In the mid-fifteenth century, one of the largest eruptions of the last 10 000 years occurred in the Central New Hebrides arc, forming the Kuwae caldera (12x6 km). This eruption followed a late maar phase in the pre-caldera edifice, responsible for a series of alternating hydromagmatic deposits and airfall lapilli layers. Tuffs related to caldera formation ( 120 m of deposits on a composite section from the caldera wall) were emitted during two main ignimbritic phases associated with two additional hydromagmatic episodes. The lower hydromagmatic tuffs from the precaldera maar phase are mainly basaltic andesite in composition, but clasts show compositions ranging from 48 to 60% SiO₂. The unwelded and welded ashflow deposits from the ignimbritic phases and the associated intermediate and upper hydromagmatic deposits also show a wide compositional range (60–73% SiO₂), but are dominantly dacitic. This broad compositional range is thought to be due to crystal fractionation. The striking evolution from one eruptive style (hydromagmatic) to the other (magmatic with emission of a large volume of ignimbrites) which occurred either over the tuff series as a whole, or at the beginning of each ignimbritic phase, is the most impressive characteristic of the caldera-forming event. This strongly suggests triggering of the main eruptive phases by magma-water interaction. A three-step model of caldera formation is presented: (1) moderate hydromagmatic (sequences HD 1–4) and magmatic (fallout deposits) activity from a central vent, probably over a period of months or years, affected an area slightly wider than the present caldera. At the end of this stage, intense seismic activity and extrusion of differentiated magma outside the caldera area occurred; (2) unhomogenized dacite was released during a hydromagmatic episode (HD 5). This was immediately followed by two major pyroclastic flows (PFD 1 and 2). The vents spread and intense magma-water interaction at the beginning of this stage decreased rapidly as magma discharge increased. Subsequent collapse of the caldera probably commenced in the southeastern sector of the caldera; (3) dacitic welded tuffs were emplaced during a second main phase (WFD 1–5). At the beginning of this phase, magma-water interaction continued, producing typical hydromagmatic deposits (HD 6). Caldera collapse extended to the northern part of the caldera. Previous C¹⁴ dates and records of explosive volcanism in ice from the south Pole show that the climactic phase of this event occurred in 1452 A.D.     

The upper mantle structure of the yakutian kimberlite province,eastern Russia

Deep seismic sounding in the region of the Mirnyi kimberlite field indicates that the boundary velocity of the uppermost mantle is elevated (v _b=8.6–8.8 km/sec) and extremely variable near the Mir kimberlite pipe. These velocity heterogeneities are probably associated with the kimberlite magmatism and may be useful in the identification of other kimberlite fields.     

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

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

Estimation of statistical parameters for censored lognormal hydraulic conductivity measurements

Field tests of hydraulic conductivity (e.g., injection test, pumping test, etc.) in low permeability formations are subject to censoring due to the detection limit of the instruments used. An iterative method of estimating the mean and variance of hydraulic conductivity data with a presumed log-normal distribution function is presented. This method accounts for the data that are actually below the lower detection limit (called truncated data) and thus gives distribution parameters that are more representative for the underlying distribution. The proposed method is then tested on two simulated normally distributed random datasets having different variances. The results show that the means and variances estimated by the proposed method are very accurate. Finally, the method is used to estimate the mean and variance of hydraulic conductivity data from single hole water injection tests in a fractured geological formation.