基于改进的风险矩阵法的生命线网络地震脆弱性分析

生命线网络的脆弱性不单单只表示地震发生后对网络作用而产生的后果,还应该包括网络的连通情况。本文在重新确定生命线网络脆弱性定义的基础上,运用风险评估理论中的风险矩阵方法综合考虑生命线网络的连通性能和失效后果两个方面来评价生命线网络的脆弱性,并以一个供气管网为例说明改进的风险矩阵法评价生命线网络脆弱性的有效性和合理性,找出供气管网中脆弱性等级最高的节点,分析其脆弱性等级最高的原因,以便于重点保护,并降低网络的脆弱性。     

Lava flow hazard and risk at Mt. Cameroon volcano

Mt. Cameroon is one of the most active effusive volcanoes in Africa. About 500,000 people living or working around its fertile flanks are subject to significant threat from lava flow inundation. Lava flow hazard and risk were assessed by simulating probable lava flow paths using the DOWNFLOW code. The vent opening probability density function and lava flow length distribution were determined on the basis of available data from past eruptions at Mt. Cameroon volcano. Code calibration was performed through comparison with real lava flow paths. The topographic basis for simulations was the 90-m resolution SRTM DEM. Simulated lava flows from about 80,000 possible vents were used to produce a detailed lava flow hazard map. The lava flow risk in the area was mapped by combining the hazard map with digitized infrastructures (i.e., human settlements and roads). Results show that the risk of lava flow inundation is greatest in the most inhabited coastal areas comprising the town of Limbe, which constitutes the center of Cameroon’s oil industry and an important commercial port. Buea, the second most important town in the area, has a much lower risk although it is significantly closer to the summit of the volcano. Non-negligible risk characterizes many villages and most roads in the area surrounding the volcano. In addition to the conventional risk mapping described above, we also present (1) two reversed risk maps (one for buildings and one for roads), where each point on the volcano is classified according to the total damage expected as a consequence of vent opening at that point; (2) maps of the lava catchments for the two main towns of Limbe and Buea, illustrating the expected damage upon venting at any point in the catchment basin. The hazard and risk maps provided here represent valuable tools for both medium/long-term land-use planning and real-time volcanic risk management and decision making.     

昆明市区地震震害的不均匀分布特点及其重复性

本文以澜沧耿马地震时昆明震害调查资料为基础,结合地震史料,研究了昆明市区震害的不均匀分布,烈度有两度的差别。相对高烈度区与历史强震破坏的重震害点大部分重合。据此分析了震害重复特点及其与地质构造和沉积盖层的关系。本文可用于研究昆明市区及邻区强震对市区可能造成震害的大致分布,为震害预测、防震抗震和城市规划提供一些依据。     

Road Network and Damaged Buildings in Urban Areas: Short and Long-term Interaction

In this paper, a model able to analyse the seismic behaviour of road network in urban areas, considering interaction between buildings and roads is presented. Damage to buildings and short-term countermeasures, such as propping, can affect roads and even block them, reducing capacity of the road network. Two successive phases are considered. In the first, immediately after the seismic event, emergency services have to reach the relevant buildings. Here, the network topology is of interest. In the second, in the long term, network capacity is not yet completely restored and network demand has to account for displaced people due to unusable buildings. In this case, road serviceability is of interest. In order to consider uncertainties on building behaviour, a probabilistic approach is developed and the network is analysed by means of Monte Carlo simulation. The methodology is applied to the municipality of Potenza (southern Italy) evaluating in the short term, the probability of strategic buildings are not connected and, in the long term, road serviceability.     

Extraction of physically realistic pore network properties from three-dimensional synchrotron X-ray microtomography images of unconsolidated porous media systems

This paper presents application of a series of algorithms used to extract pore network structure from high-resolution three-dimensional synchrotron microtomography images of unconsolidated porous media systems. These algorithms are based on the three-dimensional skeletonization that simplifies the pore space to networks in the form of nodes connected to paths. Dilation algorithms were developed to generate inscribed spheres on the nodes and paths of the medial axis to represent pore-bodies and pore-throats of the network, respectively. The end result is a physically representative pore network structure, i.e. three-dimensional spatial distribution (i.e. x-, y-, and z-coordinates) of pore-bodies and pore-throats, pore-body size distribution, pore-throat size distribution, and the connectivity. Systems analyzed in this study include different glass bead systems and natural marine sand. The media ranged in size from 0.123 to 1.0 mm, while the image volumes ranged between 7.7 and 108.9 mm³. In addition to extracting the pore network structure, the porosity, specific surface area, and representative elementary volume analysis on the porosity were calculated. Spatial correlation between pore-body sizes in the network was investigated using semivariograms and integral scale concepts. The impact of resolution on the calculated property was also investigated.

In this work, we show that microtomography is an effective tool to non-destructively extract the structure of many systems. The quality of the datasets depends on photon energy, photon flux, size of the sample, type of the sample, and size of the sample ‘features’. Results show that the developed method of extracting pore network structure is applicable to ideal and natural porous media systems. The impact of resolution on the quantification of the network structure properties varies in its significance based on feature size of the system and the properties being calculated. Therefore, a thorough resolution sensitivity analysis should be carried out to determine the degree of error associated with a system imaged at a given resolution.     

Distribution of surface imperviousness in small urban catchments predicts runoff peak flows and stream flashiness

Urban growth is a global phenomenon, and the associated impacts on hydrology from land development are expected to increase, especially in peri‐urban catchments. It is well understood that greater peak flows and higher stream flashiness are associated with increased surface imperviousness and storm location. However, the effect of the distribution of impervious areas on runoff peak flow response and stream flashiness of peri‐urban catchments has not been well studied. In this study, a new geometric index, Relative Nearness of Imperviousness to the Catchment Outlet (RNICO), is defined to correlate imperviousness distribution of peri‐urban catchments with runoff peak flows and stream flashiness. Study sites include 21 suburban catchments in New York representing a range of drainage area from 5 to 189 km² and average imperviousness from 10% to 48%. On the basis of RNICO, all development patterns are divided into 3 classes: upstream, centralized, and downstream. Results showed an obvious increase in runoff peak flows and decrease in time to peak when moving from upstream to centralized and downstream urbanization classes. This indicates that RNICO is an effective tool for classifying urban development patterns and for macroscale understanding of the hydrologic behavior of small peri‐urban catchments, despite the complexity of urban drainage systems. We also found that the impact of impervious distribution on runoff peak flows and stream flashiness decreases with catchment scale. For small catchments (A < 40 km²), RNICO was strongly correlated with the average (R² = .95) and maximum (R² = .91) gaged peak flows due to the relatively efficient subsurface routing through stormwater and sewer networks. Furthermore, the Richards–Baker stream flashiness index in small catchments was positively correlated with fractional impervious area (R² = .84) and RNICO (R² = .87). For large catchments (A > 40 km²), the impact of impervious surface distribution on peak flows and stream flashiness was negligible due to the complex drainage network and great variability in travel times. This study emphasizes the need for greater monitoring of discharge in small peri‐urban catchments to support flood prediction at the local scale.     

基于模拟退火算法的供水管网抗震优化设计

本文以管网造价为优化目标,管网拓扑结构为优化参数,管网节点最低可靠度为约束条件,建立了供水管网抗震拓扑优化模型。在供水管网功能可靠性分析方法的基础上,结合单元概率重要度分析方法,利用模拟退火算法提出了供水管网的抗震拓扑优化方法,并对一典型供水管网进行了拓扑优化分析。分析表明,对于管网抗震拓扑优化这样一个组合优化问题,模拟退火算法提供了一类很好的途径,以此为基础进行供水管网抗震可靠性优化设计具有很好的效果。     

Meteorological monitoring of an active volcano: Implications for eruption prediction

Rainfall data collected on and around the Soufriere Hills Volcano, Montserrat between 1998 and 2003 were analysed to assess the impact on primary volcanic activity, defined here as pyroclastic flows, dome collapses, and explosions. Fifteen such rainfall-triggered events were identified. If greater than 20 mm of rain fell on a particular day, the probability of a dome collapse occurring on that day increased by a factor of 6.3% to 9.2%, compared to a randomly chosen day. Similarly, the probability of observing pyroclastic flows and explosions on a day with > 20 mm of rainfall increased by factors of 2.6 and 5.4, respectively. These statistically significant links increased as the rainfall threshold was increased. Seventy percent of these rainfall-induced dome collapse episodes occurred on the same calendar day (most within a few hours) as the onset of intense rainfall, but an extra 3 occurred one or two calendar days later. The state of the volcano was important, with the rainfall–volcanic activity link being strongest during periods of unstable dome growth and weakest during periods of no dome growth or after a recent major collapse.Over 50% of the heavy rain days were associated with large-scale weather systems that can potentially be forecast up to a few days ahead. However, the remaining heavy rain days were associated with small-scale, essentially unpredictable weather systems. There was significant variability in the amount of rainfall recorded by different rain gauges, reflecting topographic variations around the volcano but also the inherent small-scale variability within an individual weather system. Hence, any monitoring/warning program is recommended to use a network, rather than just a single gauge. The seasonal cycle in rainfall was pronounced, with nearly all the heavy rain days occurring in the May–December wet season. Hence, the dome was at its most vulnerable at the beginning of the wet season after a period of uninterrupted growth. Interannual variability in rainfall was related to tropical Pacific and Atlantic sea surface temperature anomalies, and holds out the prospect of some limited skill in volcanic hazard forecasts at even longer lead times.     

油气储运系统钢结构设备地震火灾效应研究

预防和减轻油气储运系统地震火灾危险性是防震减灾中的重要课题。本文论述了油气储运系统的地震火灾特征及防火要求，对钢结构的防火涂料在遭到破坏和平时情况下的抗火性能进行了对比分析，并应用有限元法对油气储运系统中典型钢结构设备在地震火灾荷载作用下的抗火性能进行了数值模拟。研究结果表明，防火涂料遭到强地震荷载破坏后，钢结构的抗火性能大大降低。     

基于微粒群算法的供水管网抗震优化设计

以管网年费用折算值为优化目标、管网拓扑结构与管径为优化参数、管网节点最低可靠度为约束条件,建立了供水管网抗震优化设计模型。利用微粒群算法对这一模型进行了求解,该算法以管网作为微粒个体,通过不断地更新微粒的位置来搜索最优的管网结构,直到最后给出优化的管网结构。利用上述方法对一典型供水管网进行了抗震优化设计分析,给出了3种不同节点最低可靠度约束条件下的优化改造方案。     

Compound and simple lava flows and flood basalts

Compound lava flows, defined as those lavas which are divisible into flowunits, commonly have a shield-like form and are thought to develop when the rate of extrusion of lava is relatively low.Simple lava flows, defined as those lavas which are not divisible into flow-units, are thought to form when the rate of extrusion of lava is relatively high. A logical definition oflava flow must embrace both simple lava sheets and substantial lava shields (compound lava flows) up to 600 m high. Flood basalt piles include both compound and simple flows, but the most extensive and far-reaching flows are simple and they are believed to form when the rate of extrusion of lava is particularly high.     

Genetic algorithm for seismic topology optimization of lifeline network systems

Lifeline systems, such as water distribution and gas supply networks, usually cover large areas. For these systems, seismic design is always a difficult problem because of the complexity of large‐scale networks. In this paper, a topology optimization technology for lifeline networks is established. Firstly, in order to speed up the convergence of optimization process, an element investment importance analysis is carried out to evaluate the importance of components to the lifeline network. Then a topology optimization model is established. The aim of the model is to find the least‐cost network topology while the seismic reliability between the sources and each terminal satisfies prescribed reliability constraints. For this optimization problem, a genetic algorithm, which takes network topologies as the individuals of its population, is used to search for the optimal solutions by suitable operators, including selection, crossover and mutation operators. The capacity of the proposed algorithm is illustrated by its applications to a simple example network consisting of 10 nodes and an actual network with 391 nodes located in a large city of China. Copyright © 2008 John Wiley & Sons, Ltd.     

网络可靠度分析的最小路算法和最小割算法研究

网络可靠度分析是评价城市生命线工程系统整体抗震性能的主要手段。本文分别从最小路和最小割的角度介绍了网络可靠度分析算法,包括:经典不交最小路(割)算法、最小路(割)递推分解算法和改进最小路(割)递推分解算法。在此基础上,通过实例分析,着重进行了改进最小路递推分解算法和改进最小割递推分解算法的对比分析,分析结果表明两种算法在网络单元不同可靠度水平下具有不同的计算效率,并对引起以上区别的三个主要原因进行了分析。     

长白山天池火山千年大喷发火山碎屑流堆积的粒度特征与地质意义

火山碎屑流堆积因其巨大的危害性而成为火山学研究的重要课题之一。文中对长白山区天池火山千年大喷发火山碎屑流堆积进行了粒度分析。分析结果表明,火山碎屑流堆积分选差,伴生的灰云浪堆积分选性好。火山碎屑流堆积中的岩屑和浮岩的平均最大粒度随着离开火山口距离的加大而减小,反映火山碎屑流搬运过程中存在重力分选和机械磨损作用。火山碎屑流在搬运过程中发生了流体化作用,离开火山口越远流体化速度越小,反映火山碎屑流在搬运过程中发生了流体化去气作用,这种作用使火山碎屑流的粘度和屈服强度增加而导致沉积。流体化速度是控制搬运距离的因素之一,远源多渠道火山碎屑流的汇合使流体化速度增大,搬运距离更远,造成的灾害范围也增大     

Coupled 1D–3D hydrodynamic modelling,with application to the Pearl River Delta

Within the hydrodynamic modelling community, it is common practice to apply different modelling systems for coastal waters and river systems. Whereas for coastal waters 3D finite difference or finite element grids are commonly used, river systems are generally modelled using 1D networks. Each of these systems is tailored towards specific applications. Three-dimensional coastal water models are designed to model the horizontal and vertical variability in coastal waters and are less well suited for representing the complex geometry and cross-sectional areas of river networks. On the other hand, 1D river network models are designed to accurately represent complex river network geometries and complex structures like weirs, barrages and dams. A disadvantage, however, is that they are unable to resolve complex spatial flow variability. In real life, however, coastal oceans and rivers interact. In deltaic estuaries, both tidal intrusion of seawater into the upstream river network and river discharge into open waters play a role. This is frequently approached by modelling the systems independently, with off-line coupling of the lateral boundary forcing. This implies that the river and the coastal model run sequentially, providing lateral discharge (1D) and water level (3D) forcing to each other without the possibility of direct feedback or interaction between these processes. An additional disadvantage is that due to the time aggregation usually applied to exchanged quantities, mass conservation is difficult to ensure. In this paper, we propose an approach that couples a 3D hydrodynamic modelling system for coastal waters (Delft3D) with a 1D modelling system for river hydraulics (SOBEK) online. This implies that contrary to off-line coupling, the hydrodynamic quantities are exchanged between the 1D and 3D domains during runtime to resolve the real-time exchange and interaction between the coastal waters and river network. This allows for accurate and mass conserving modelling of complex coastal waters and river network systems, whilst the advantages of both systems are maintained and used in an optimal and computationally efficient way. The coupled 1D–3D system is used to model the flows in the Pearl River Delta (Guangdong, China), which are determined by the interaction of the upstream network of the Pearl River and the open waters of the South China Sea. The highly complex upstream river network is modelled in 1D, simulating river discharges for the dry and wet monsoon periods. The 3D coastal model simulates the flow due to the external (ocean) periodic tidal forcing, the salinity distribution for both dry and wet seasons, as well as residual water levels (sea level anomalies) originating from the South China Sea. The model is calibrated and its performance extensively assessed against field measurements, resulting in a mean root mean square (RMS) error of below 6% for water levels over the entire Pearl River Delta. The model also represents both the discharge distribution over the river network and salinity transport processes with good accuracy, resolving the discharge distribution over the main branches of the river network within 5% of reported annual mean values and RMS errors for salinity in the range of 2 ppt (dry season) to 5 ppt (wet season).     

Disruption index,DI: an approach for assessing seismic risk in urban systems (theoretical aspects)

Urban systems are characterized by very complex interactions. After an earthquake, a wide variety of services, networks and urban facilities may be unavailable to the public during the system failure and recovery processes, thereby causing disruptions in the basic social needs of the affected area. After a disaster, communities face several challenges. For example, the lack of education may impose population migrations, or malfunctions in the electricity distribution system can produce electrical power outages of varying duration with respect to time and space, which generates consequences in the water distribution system, transportation, communications, etc. A methodology called the Disruption index (DI), based on graph theory, includes these multiple interdependencies. It has been developed to estimate the dysfunction of some fundamental dimensions of urban systems on a broad level, starting with the physical damages directly suffered by the exposed assets, proceeding to the impacts that each node has on the functional performance of the nodes depending on them, until reaching the top node. This paper presents the fundamental theory to support the DI concept. The DI provides the likely impacts and consequences of an earthquake in an urban area to fulfill hazard mitigation and provide civil protection agencies and local and state governments with a new decision-making instrument to reduce or prevent severe and recurrent impacts. The DI concept can also be extended to other natural and man-made disasters and may be used as a tool for optimizing the resources of the system components.     

Process oriented use of geostatistics to analyse creeping para‐glacial features

Complex para‐glacial systems may show signs of destabilization leading to frequent and potentially hazardous debris‐flows. Understanding creeping permafrost displacement over a period of time is therefore crucial for hazard management and risks assessment. This paper presents our methodology for estimating creeping permafrost displacement based upon data derived from various survey methods and demonstrates its relevance on the Glacier Bonnard system in southern Switzerland. Geostatistical processing allowed estimation of the displacement intensity over the area of interest, as well as assessment of the interpolation quality. Although the local measurement network needs to be refined, the results largely improve the understanding of the Bonnard para‐glacial system, though highlighting the need to locally refine the measurement network. In the present case, the destabilized front advances at ~1 m/a and the upstream creeping part at ~0.4 m/a. Variance analysis also provide objective thresholds that could be used to distinguish underlying physical processes. Copyright © 2015 John Wiley & Sons, Ltd.     

How wide is a road? The association of roads and mass-wasting in a forested montane environment

A spatial data base of 1609 landslides was analysed using a geographic information system to determine landslide frequency in relation to highways. A 126 km long transportation network in a 201 km² area of humid-tropical, mountainous, forested terrain in Puerto Rico was used in conjunction with a series of 20 buffer (disturbance) zones varying from 5 to 400 m in length, measured perpendicular to the highways. Average landslide frequency in the study area at distances greater than 85 m from roads was about six landslides per square kilometre. At distances of 85 m or less on either side of a highway, landslide frequency was about 30 landslides per square kilometre. On average, this elevated disturbance rate affected 330 m² km⁻² a⁻¹ within the 170 m swath. The mass-wasting rate outside of the disturbance zone affected 40 m² km⁻² a⁻¹. These results indicate that the rate of mass-wasting disturbance is increased from five to eight times in a 170 m wide swath along road corridors. The lateral extent of the environmental impact of roads in the study area is greater than is commonly perceived. The approach described herein demonstrates a simple method to assess the spatial association of mass-wasting with highways. © 1997 by John Wiley & Sons, Ltd.