地震灾害预测和应急模拟系统的设计与应用——以永安市城市应急系统为例

城市防震减灾是地震潜在频发区提高自然灾害防治能力建设的重要方面。本文结合福建省永安市城市防震减灾信息管理系统的建设,探讨设计防震减灾数值模型,分析模型集成的关键技术,构建了地震风险评估、建筑物易损性评价、生命与财产损失估算、救援与救灾管理调度等模型和基于GIS开发的震害预测和应急模拟系统。该系统在福建省永安市的实际应用表明,通过对地震灾害预测结果的分析,可加强城市抗震中的薄弱环节,为灾区政府应急响应和制订对策提供决策支持辅助信息,从而显著提高城市防震减灾的综合能力。     

Seismic hazard in central southern Africa

Seismic hazard maps of central-southern Africa where hazard has been expressed in terms of peak ground acceleration for an annual probability in excess of 10^-1 show relatively high values that distinguish the seismic hazard potential of the Deka fault zone, the mid-Zambezi basin-Luangwa rift and western central Mozambique. In areas such as central-southern Africa where little is known about the geology of the region and the fault systems have not been fully mapped, seismic hazard potential may be estimated from seismicity and broad-scale fault features. For this region, such potential is based on earthquake magnitude M_s ≥ 6. Events of such magnitude have recently occurred in the mid-Zambezi basin, southern Zimbabwe and western-central Mozambique. This paper follows the conventional probabilistic hazard analysis procedure, defining seismic source zones from seismicity based on instrumental records from a cataloque that spans a period of 83 years. Geological and geomorphological features in the region are described on the mesoscale and are correlated with the seismicity as broad fault zones. The scarcity of strong-motion accelerogram data necessitated the formulation of attenuation values based on random vibration theory (RVT).     

中国地震灾害人口死亡风险定量评估(英文)

基于地震灾害风险形成机理,在建立人口震害脆弱性曲线与确定地震发生参数的基础上,本文利用评估模型对我国Ⅴ-Ⅺ地震烈度下各县域单元的人口死亡风险进行评估并分析其空间分布格局。主要研究内容有:(1)首次采用基于过去—现在—未来的多方面地震孕灾环境资料来处理地震发生的可能性。具体综合历史地震综合烈度、地震活动断裂带分布、地震动峰值加速度三方面来确定全国2355个县域单元的地震发生参数;(2)利用1990-2009年我国历史地震灾情数据,对地震烈度与人员死亡率之间进行线性拟合,建立适合我国地震灾害风险评估的震害人口死亡脆弱性曲线;(3)利用震害风险评估模型对我国各县域单元的人口死亡风险进行定量评估,并分析风险空间分布格局,彻底摸清Ⅴ-Ⅺ地震烈度下我国各县域单元的地震灾害人口死亡风险。研究表明:在不同地震烈度下,我国广大的东、中部地区面临更高的风险,而西部的人口死亡风险相对较低。高风险区域呈零星状分布于山东与江苏大部、安徽北部、黑龙江与吉林东部等人口分布较密集且孕灾环境发育完备的区域。而无风险区域在全国范围内呈斑块状散布,分布格局基本保持不变。     

Earthquake nests

Summary. Making use of extensive observations of micro-earthquakes available in the Tokyo region, Usami & Watanabe constructed a simple but effective method of exhibiting those regions where small earthquakes have clustered during a 4.5 yr period. Such clusters are defined by contour surfaces, and for convenience are termed earthquake nests .
If we accept the assumption that the magnitude of an earthquake is directly related to the volume previously undergoing intensified strain, and if we regard the nests as measures of such volumes, we can make rough estimates of the magnitudes of potential earthquakes in different parts of the Tokyo region. This method of assigning seismicity can be one way of calculating local risk. Repeated investigations can also detect changes in the tress field.     

Constraints on the frequency–magnitude relation and maximum magnitudes in the UK from observed seismicity and glacio-isostatic recovery rates

Earthquake populations have recently been shown to have many similarities with critical-point phenomena, with fractal scaling of source sizes (energy or seismic moment) corresponding to the observed Gutenberg–Richter (G–R) frequency–magnitude law holding at low magnitudes. At high magnitudes, the form of the distribution depends on the seismic moment release rate M˙ and the maximum magnitude m _max . The G–R law requires a sharp truncation at an absolute maximum magnitude for finite M˙ . In contrast, the gamma distribution has an exponential tail which allows a soft or 'credible' maximum to be determined by negligible contribution to the total seismic moment release. Here we apply both distributions to seismic hazard in the mainland UK and its immediate continental shelf, constrained by a mixture of instrumental, historical and neotectonic data. Tectonic moment release rates for the seismogenic part of the lithosphere are calculated from a flexural-plate model for glacio-isostatic recovery, constrained by vertical deformation rates from tide-gauge and geomorphological data. Earthquake focal mechanisms in the UK show near-vertical strike-slip faulting, with implied directions of maximum compressive stress approximately in the NNW–SSE direction, consistent with the tectonic model. Maximum magnitudes are found to be in the range 6.3–7.5 for the G–R law, or 7.0–8.2 m _L for the gamma distribution, which compare with a maximum observed in the time period of interest of 6.1 m _L . The upper bounds are conservative estimates, based on 100 per cent seismic release of the observed vertical neotectonic deformation. Glacio-isostatic recovery is predominantly an elastic rather than a seismic process, so the true value of m _max is likely to be nearer the lower end of the quoted range.     

Seismic quiescence before the M 7, 1988, Spitak earthquake, Armenia

A detailed analysis of the 35  yr of seismicity between 1962 and 1997 using a gridding technique shows that the M 7, Spitak earthquake of 1988 December 7 was preceded by a quiescence anomaly that started at approximately 1984±0.5, and lasted about 5±0.5  yr, up to the main shock. This quiescence anomaly had a radius of about 20±3  km, estimated from circular areas with 75 per cent rate decrease, centred at the point of maximum significance of the anomaly. The quiescence was clearly present in the aftershock volume during the 5  yr before the 1988 main shock, but its statistically strongest expression was located 30  km NW of the epicentre. This anomaly fulfills the association rules between precursory quiescence anomalies and main shocks, even for a tight definition, and is therefore proposed as a case of precursory quiescence. The largest value of the standard deviate Z , found by random selection of samples by gridding, was Z =14 for a time window of T _w=3  yr, using a sample size of N =300 events. This makes this anomaly the strongest observed so far, and it is the first documented in an environment of continental collision. There are no false alarms exceeding in significance the precursor. The Armenian earthquake catalogue used for this study had 4600 earthquakes with M ≥ M _min=2.2 in the area bounded by 39.5° to 42°N/42.5° to 47°E. From the point of view of homogeneous reporting this is the best catalogue we have analysed so far. The limits of the data used and the density of the grid are dictated by the data, and have no influence on the results. The choice of free parameters does not influence the results significantly within the following limits: 100≤ N ≤500, 2≤ T _w≤7, 2.2≤ M _min≤2.8.     

An improved region–time–length algorithm applied to the 1999 Chi-Chi, Taiwan earthquake

By means of the region–time–length (RTL) algorithm, which is widely used for investigating the precursory seismicity changes in China, Italy, Japan, Russia and Turkey, we examine the precursory seismic activity occurred prior to the 1999, M _w = 7.6, Chi-Chi earthquake around its epicentre. Based on our calculation of the RTL values, the epicentral area has been found to strongly exhibit the signature of anomalous activity, associated with the seismic quiescence and activation, before the main shock. Also proposed in this study is a helpful method for determining two important parameters used in the RTL analysis, the characteristic time and distance. Such method will largely reduce the ambiguity in the original RTL algorithm.     

Statistical analysis of time-dependent earthquake occurrence and its impact on hazard in the low seismicity region Lower Rhine Embayment

The time-dependence of earthquake occurrence is mostly ignored in standard seismic hazard assessment even though earthquake clustering is well known. In this work, we attempt to quantify the impact of more realistic dynamics on the seismic hazard estimations. We include the time and space dependences between earthquakes into the hazard analysis via Monte Carlo simulations. Our target region is the Lower Rhine Embayment, a low seismicity area in Germany. Including aftershock sequences by using the epidemic type aftershock-sequence (ETAS) model, we find that on average the hypothesis of uncorrelated random earthquake activity underestimates the hazard by 5–10 per cent. Furthermore, we show that aftershock activity of past large earthquakes can locally increase the hazard even centuries later. We also analyse the impact of the so-called long-term behaviour, assuming a quasi-periodic occurrence of main events on a major fault in that region. We found that a significant impact on hazard is only expected for the special case of a very regular recurrence of the main shocks.     

Upper-crustal structure of the Benevento area (southern Italy): fault heterogeneities and potential for large earthquakes

The Benevento region is part of the southern Apennines seismogenic belt, which experienced large destructive seismic events both in historical and in recent times. The study area lies at the northern end of the Irpinia fault, which ruptured in 1980 with a M_s = 6.9 normal faulting event, which caused about 3000 casualties. The aims of this paper are to image lateral heterogeneities in the upper crust of the Benevento region, and to try to identify the fault segments that are expected to generate such large earthquakes. This work is motivated by the recognition that lithological heterogeneities along major fault zones, inferred from velocity anomalies, reflect the presence of fault patches that behave differently during large rupture episodes. In this paper, we define the crustal structure of the Benevento region by using the background seismicity recorded during 1991 and 1992 by a local seismic array. These data offer a unique opportunity to investigate the presence of structural discontinuities of a major seismogenic zone before the occurrence of the next large earthquake. The main result that we obtained is the delineation of two NW-trending high-velocity zones (HVZs) in the upper crust beneath the Matese limestone massif. These high velocities are interpreted as high-strength regions that extend for 30-40 km down to at least 12 km depth. The correspondence of these HVZs with the maximum intensity regions of historical earthquakes (1688 AD, 1805 AD) suggests that these anomalies delineate the extent of two fault segments of the southern Apenninic belt capable of generating M = 6.5⁻⁷ earthquakes. The lateral offset observed between the two segments from tomographic results and isoseismal areas is possibly related to transverse right-lateral faults.     

Relocation of M ≥ 2 events of the 1989 Dôbi seismic sequence in Afar: evidence for earthquake migration

The Ethiopian side of central Afar was struck in August 1989 by the largest seismic sequence (three 6.1 ≤ M _s ≤ 6.3 events, 15 with M _s or m _b ≥ 5.0) since that of Serdo in 1969. Using the Djibouti seismological network, we relocated 297 of the events of that sequence. As most of the large events took place outside the network, we assessed the accuracy and stability of earthquake relocations by using three different velocity models and two relocation codes to try to relate individual shocks to distinct faults and surface breaks. A majority of the events apparently occurred underneath the floor of the Dôbi graben, an area about 45  km long and 15  km wide, rupturing boundary and inner floor faults, in agreement with the surface cracks and scarps that we mapped in the area. The relocation shows that the principal events propagated about 50  km northwestwards along the graben in the first 40  hr. A day and a half after the beginning of the sequence, smaller events ( M ≤ 4) started to propagate more than 55  km eastwards, towards Asal Lake. Using two three-component stations installed near the Ethiopian border, we could determine reliable depths for 21 events. The depths are compatible with a seismogenic crust about 14  km thick in the Dôbi and Hanle graben area. Although the Dôbi sequence ruptured about 50  km of the fault array extending from Serdo to Asal, where the regional stress was released by earthquakes in 1969 and 1978, respectively, a seismic gap about 50  km long still subsists along the northern part of the Gaggade region (Der'êla half-graben).     

Seismic risk in southern Europe through to India examined using Gumbel's third distribution of extreme values

Summary. A technique is described for the analysis of seismicity using Gumbel's third asymptotic distribution of extreme values. Seismicity of southern Europe through to India, nominally for the period 1900–74, is subdivided in a cellular manner, without recourse to tectonic discrimination between regions, and a covariance analysis on the three parameters of Gumbel's distribution is performed for each cell of seismicity. The results indicate that the upper bound to the magnitude of earthquake occurrence is often uncertain although it is discernible, while curvature of the earthquake occurrence distributions is usually established. Uncertainties in the forecasts of largest earthquakes, with a return period of 75 yr, are distinctly improved by taking into account the large and negative covariance which is measured between the curvature and upper bound to earthquake magnitude for the observed seismicity. These results are then used to map seismic risk for southern Europe through to India.     

中国地震灾害宏观人口脆弱性评估

地震灾害脆弱性不但是国家和地区中长期规划的基础,也是地震应急时的重要决策依据,因此地震灾害脆弱性评估是当前研究的热点和重点。本文首先介绍了地震灾害人口脆弱性定义,指出了地震灾害人口脆弱性的影响因子和影响模式;然后通过分析传统地震灾害人口脆弱性的宏观和微观评价方法和模型,利用人口总量、儿童人口比例、老年人口比例、人口密度和人均GDP 5 个指标建立了地震灾害宏观人口脆弱性综合评价模型,并对全国各县市区的人口脆弱性进行了评估。结果表明：中国地震灾害人口脆弱性区域差异明显,东部地区人口脆弱性明显高于西部;人口总量是脆弱性的主要因素,宏观人口脆弱性综合评估模型对地震灾害人口脆弱性的修正作用显著。     

Seismicity and seismic gap in the Lesser Antilles arc and earthquake hazard in Guadeloupe

Summary. A seismic study of the Lesser Antilles arc has been carried out, first for the period 1950–1978, for which we can use local seismic networks to draw maps of instrumental seismicity, then for the period 1530–1950, for which we have catalogues of felt earthquakes. The striking feature of the spatial distribution of foci is the cluster of epicentres in the northern half of the arc; all large earthquakes ( M > 7.5) are located north of 14° latitude. Seismicity cross-sections through the arc show a variable dipping subduction zone along the arc; the deep seismic zone is steeper in the centre of the arc than on the extremity.
The time-space diagram for historical seismicity, and the evidence of a seismic gap at the east of Guadeloupe lead us to consider the northern half arc as a likely site for a large earthquake in the near future.
The seismic slip rate calculated from all major earthquakes since 1530 is of much greater value than that obtained from recent plate tectonic models, suggesting that the recurrence rate of earthquakes is more than many hundreds of years with a possible aseismic creep.     

A catalogue of seismicity in Greece and adjacent areas

Summary. A new earthquake catalogue for Greece has been formed to cover the instrumental period 1901–78, in particular 605 earthquakes for the period 1917–63 inclusive are relocated using first arrival data from the International Seismological Summary. These relocations incorporate macro-seismically and other well-controlled master events into an ensuing joint epicentre determination technique. The largest annual average shift is 165 km for 78 earthquakes during the decade after 1917, decreasing to 17 km for the later years 1957–63. Magnitudes are redetermined mainly using readings from the Swedish network and Uppsala since as early as 1908. Catalogue completeness exists for magnitudes around 5.5 for at least the last 60 years, but magnitudes of about 4.7 are completely reported only during the most recent 15 years.
The new epicentral positions lead to a better delineation of the seismic zones than has previously been achieved. The majority of shallow earthquakes are contained in a belt parallel to the Hellenic Arc which extends north into Albania, in the south-east the epicentral locations give a more diffuse extension of this zone into the west coast of Turkey. Depths of intermediate earthquakes along the Hellenic Arc tend to relocate to shallower depths, in the south-west part of Crete no earthquake focus deeper than 100 km is found, although in the south-eastern section of the arc a tendency to increased depths is observed. A second zone starts at Leukas Island in the west and extends through central Greece to near Volos on the east coast, where it divides into two branches, which are less well defined, but which eventually join the seismicity of western Turkey. A third zone follows the Saronikos and Corinth gulfs.     

The first and third asymptotic distributions of extremes as applied to the seismic source regions of India and adjacent areas

The seismic source regions are identified on the basis of spatial and temporal distributions of shocks (1900–1989), recurrence relations and the tectonic architecture of the Indian subcontinent and adjoining areas. The probable occurrence of the maximum magnitude earthquake is estimated using the theory of extreme values of Gumbel. The parameters of the first and third asymptotic distributions of extremes and their uncertainty values are computed for the seven identified seismic source regions of India and adjacent areas. The third-type distribution curve is preferable to the first type in all the regions, as revealed by the χ² test. The results of the third asymptotic distribution indicate the upper bound to earthquake magnitude w is equal to 8.94 ± 0.21 for Assam, 8.56 ± 0.29 for Bihar-Nepal, 8.43 ± 0.10 for Kangra, 8.97 ± 0.27 for Hindukush, 7.61 ± 0.24 for Pakistan-Cutch, 7.34 ± 0.12 for Koyna and 8.98 ± 0.27 for Andaman Sea seismic source regions. The predicted most probable largest earthquake magnitude is computed for return periods of 10, 20, 50, 75 and 100 yr in each source region.     

Large-magnitude Central American earthquakes, 1898–1994

We have collected and re-examined macroseismic information for large Central American earthquakes since the beginning of the period of instrumental recording about one hundred years ago, and combined this with a reassessment of early instrumental information to produce a catalogue of 51 events that, we believe includes ail those with magnitudes ( Ms ) greater than 7.0. We have reassessed surface-wave magnitudes by consulting station bulletins and we have derived a correction that gives an equivalent Ms for events of intermediate depth. We have also developed a regional relationship between Ms and seismic moment, which enables us to estimate the seismic slip rate across the Middle American Trench. Our best estimates give an average slip rate several times smaller than suggested convergence rates, but with the seismic slip in the central segment of the trench almost an order of magnitude smaller than that in the segments on either side. The low seismic slip rate may indicate aseismic crustal deformation     

Station threshold bias in short-period amplitude distance and station terms used to compute body-wave magnitudes mb

A cause of bias in estimates of distance and station terms used to compute seismic magnitudes ( m _b) is demonstrated in a comparative study on the world amplitude–distance curve. This bias results from the presence of station thresholds in the measurement and reporting of seismic amplitudes and is analogous to that known to be present in routinely computed magnitudes. Distance terms corresponding to low amplitudes such as in the core shadow will tend to be positively biased as are station terms for stations with high thresholds.     

Seismic tomographic imaging of P- and S-waves velocity perturbations in the upper mantle beneath Iran

The inverse tomography method has been used to study the P - and S -waves velocity structure of the crust and upper mantle underneath Iran. The method, based on the principle of source–receiver reciprocity, allows for tomographic studies of regions with sparse distribution of seismic stations if the region has sufficient seismicity. The arrival times of body waves from earthquakes in the study area as reported in the ISC catalogue (1964–1996) at all available epicentral distances are used for calculation of residual arrival times. Prior to inversion we have relocated hypocentres based on a 1-D spherical earth's model taking into account variable crustal thickness and surface topography. During the inversion seismic sources are further relocated simultaneously with the calculation of velocity perturbations. With a series of synthetic tests we demonstrate the power of the algorithm and the data to reconstruct introduced anomalies using the ray paths of the real data set and taking into account the measurement errors and outliers. The velocity anomalies show that the crust and upper mantle beneath the Iranian Plateau comprises a low velocity domain between the Arabian Plate and the Caspian Block. This is in agreement with global tomographic models, and also tectonic models, in which active Iranian plateau is trapped between the stable Turan plate in the north and the Arabian shield in the south. Our results show clear evidence of the mainly aseismic subduction of the oceanic crust of the Oman Sea underneath the Iranian Plateau. However, along the Zagros suture zone, the subduction pattern is more complex than at Makran where the collision of the two plates is highly seismic.     

Seismic reflection and transmission coefficients of a self-similar interface

The derivation of seismic reflection and transmission coefficients is generally based on the assumption that the medium parameters behave as step functions of depth, at least in a finite region around the interface. However, outliers observed in well logs generally behave quite differently from step functions. In this paper we represent an interface by a self-similar singularity, embedded between two homogeneous half-spaces, and we derive its frequency-dependent normal-incidence reflection and transmission coefficients. For ω  → 0 the expressions for the coefficients reduce to those for a discrete boundary between two homogeneous half-spaces; for ω → ∞ they become frequency-independent. These asymptotic expressions have a relatively simple form and depend on the singularity exponent α .
  The exact as well as the asymptotic expressions are used to evaluate the time-domain reflection and transmission responses of a self-similar interface. Finally, we use a numerical method to model the response of a smoothed version of a self-similar interface (note that the velocity of a smoothed singularity remains finite). It turns out that smoothing has hardly any effect on the response, provided that the smoothing does not affect the scales corresponding to the seismic frequency range.