MUD BANK SLIDING IN A DAM RESERVOIR

l INTRODUCTIONTo ensure the long-term futhe of dams, the emPtying of storage reservoirs is often realised with tWoobjectives:'-periodic insPection of the sttucthe,-flushing oftalrials Which sedimatd within the reservoir.Both aims are lthely to cause the instability of sedhoent, and these displacements can induceenvirOInnenal distheances that have to be reduced to the Ininimum. This imPlies a good knowledge ofthese maerials and of the processes Which originated the disPlacements.The se…     

On the implementation of Prabhu’s exact solution of the stochastic reservoir equation

After 50 years of Prabhu’s paper on the exact solution of the stochastic reservoir equation for the important class of gamma inflow distributions with an integral shape parameter, a detailed implementation of the exact solution is still lacking, despite its potential usefulness from both theoretical and practical points of view. This paper explores some properties of Prabhu’s exact solution and investigates the numerical difficulties associated with its implementation. The solution is also extended to derive the distributions of deficit, spillage, yield, and actual release from the reservoir. Explicit analytical solutions for three relatively simple cases are given in detail as examples and comparisons with approximate numerical solutions are made, which reveal some shortcomings of approximate methods. The implementation of the solution in the general case reveals some numerical problems associated with large values of the shape parameter of the inflow distribution and large ratios of reservoir size to draft, mainly due to accumulation of round-off errors. A Matlab program has been developed to calculate emptying and filling probabilities over a wide range of reservoir parameters using extended precision. Comparison of Prabhu’s solution with the numerical solution of the reservoir integral equation highlights possible problems with the numerical solution, which may produce inaccurate or even invalid results for large reservoirs, large drift, and large skewness of the inflow distribution. A comparison between gamma and lognormal distributions as models of skew revealed that as the reservoir size, drift, and skewness increase, the probability of emptying of the reservoir becomes smaller for the case of gamma inflow than in the case of lognormal flow having the same skewness coefficient.     

液化型路堤边坡地震安全性分析

地震液化型路堤边坡安全性问题涉及岩土工程与工程地震两个学科领域,是边坡工程与砂土液化的交叉课题。路堤边坡在地震时,因基础之下的砂土发生液化现象,会急剧发生地基承载力丧失、液化沉陷、侧向流动、液化喷砂、整体滑动等问题。应用NCEER液化评估方法(SPT)、Mononobe-Ok-abe方法、Koerner方法和Fellen ius方法,以抗液化安全系数、液化侧向流动安全系数、液化喷砂安全系数和抗滑安全系数为评价指标,系统对地震液化型路堤边坡安全性进行了分析、计算和评价。结果显示,文中方法对地震液化型路堤边坡安全性分析具有适用性。     

江西省东津水库诱发地震可能性研究

本文从库区的地层岩性、断裂构造、水文地质条件和库水能起的作用'判断，认为拟建的东津水库蓄水后发生水库诱发地城的可能性很少。建立在相似地质构造背景和水文地质条件上的黄龙滩水库，至今未发生水库诱发地震就是一个很好的例证。     

SEDIMENTATION IN LOW HEAD RIVER RESERVOIRS

ＳＥＤＩＭＥＮＴＡＴＩＯＮＩＮＬＯＷＨＥＡＤＲＩＶＥＲＲＥＳＥＲＶＯＩＲＳ￥ＷｉｌｈｅｌｍＢｅｃｈｔｅｌｅｒ（Ｐｒｏｆ．Ｄｒ．－Ｉｕｇ．ＩｎｓｔｉｔｕｔｅｏｆＨｙｄｒｏｓｃｉｅｎｃｅｓ，ＦｅｄｅｒａｌＡｒｍｅｄＦｏｒｃｅｓＵｎｉｖｅｒｓｉｔｙＭｕｎｉ...     

Hydrothermal calderas

The standard model of caldera formation is related to the emptying of a magma chamber and ensuing roof collapse during large eruptions or subsurface withdrawal. Although this model works well for numerous volcanoes, it is inappropriate for many basaltic volcanoes (with the notable exception of Hawaii), as these have eruptions that involve volumes of magma that are small compared to the collapse. Many arc volcanoes also have similar oversized depressions, such as Poas (Costa Rica) and Aoba (Vanuatu). In this article, we propose an alternative caldera model based on deep hydrothermal alteration of volcanic rocks in the central part of the edifice. Under certain conditions, the clay-rich altered and pressurized core may flow under its own weight, spread laterally, and trigger very large caldera-like collapse. Several specific mechanisms can generate the formation of such hydrothermal calderas. Among them, we identify two principal modes: mode 1: ripening with summit loading and flank spreading and mode II: unbuttressing with flank subsidence and flank sliding. Processes such as summit loading or flank subsidence may act simultaneously in hybrid mechanisms. Natural examples are shown to illustrate the different modes of formation. For ripening, we give Aoba (Vanuatu) as an example of probable summit loading, while Casita (Nicaragua) is the type example of flank spreading. For unbuttressing, Nuku Hiva Island (Marquesas) is our example for flank subsidence and Piton de la Fournaise (La Réunion) is our example of flank sliding. The whole process is slow and probably needs (a) at least a few tens of thousands of years to deeply alter the edifice and reach conditions suitable for ductile flow and (b) a few hundred years to achieve the caldera collapse. The size and the shape of the caldera strictly mimic that of the underlying weak core. Thus, the size of the caldera is not controlled by the dimensions of the underlying magma reservoir. A collapsing hydrothermal caldera could generate significant phreatic activity and trigger major eruptions from a coexisting magmatic complex. As the buildup to collapse is slow, such caldera-forming events could be detected long before their onset.     

小浪底水库断裂构造分析与诱发地震预测

对小浪底水库区城崖地，塔底，石井河等主要活动断层的基本特征作了 简要介绍，采用断层强度分析法分别估算了小浪底水库蓄水前和一些断层面的抗滑稳定系数，对库区地震活动性，应力状态，岩溶和岩体渗透性等诱震条件进行了分析，对水库诱发地震的最大震级和可能的震中位置进行了初步预测。     

A regional model to assess the hydrological sensitivity of medium size catchments to climate variability

Catchment storage capacity is an important factor in the determination of catchment sensitivity to climate variability. Quantification of catchment sensitivity is in turn important in the regional assessment of the effects of possible climate change. In the present paper, an empirical regional model is proposed that quantifies catchment sensitivity as the ratio of present maximum reservoir storage to catchment storage capacity. Catchment storage capacity is defined theoretically using readily available catchment variables. Present maximum reservoir storage in a catchment, as determined from recession analysis, is expressed as a fraction of catchment storage capacity; the fraction defines catchment sensitivity and depends on storage capacity and annual net precipitation. Average annual conditions for present maximum reservoir storage and average annual net precipitation are used to test the developed model. Although the study used data from only 15 catchments in the Upper Loire region in France, the model proved statistically valid. Storage capacity calculated with the model compares favourably with the baseflow index and a storage index defined in previous research. Values of storage capacity are probable with respect to reported water resources in the area. With the model catchment sensitivity can easily be assessed. Flood or drought prone catchments can be identified as well as a catchment's sensitivity to a catchment-type transition (baseflow versus direct flow dominated catchments). © 1998 John Wiley & Sons, Ltd.     

A survey of geotechnical system identification techniques

Identification and inverse problem techniques play an important role in the characterization and modeling of geotechnical systems. These techniques have been used in estimation of system parameters, model development and calibration, as well as simulation of earthquake ground motions. The recent availability of high quality seismic records of sites equipped with downhole accelerometer arrays led to a burgeoning of identification and inverse problem studies involving geotechnical systems. This paper presents a survey of system identification techniques and analyses of full- and small-scale soil systems, with an emphasis on geotechnical earthquake engineering problems.     

GEOPHYSICAL STUDIES OF WATER LEAKAGES FROM RESERVOIRS*

Considerable water leakages from reservoirs make it difficult to attain the planned storage capacity. In some cases water leakages give rise to suffusion followed by catastrophes. Until recently methods for locating water leakages were extremely imperfect. Geophysical methods offer good prospects in this direction. For solving these problems, it is effective to use streaming potential measurements, water flow rate observations and thermometry. Laboratory experiments were carried out in connection with the fact that water leakages from reservoirs are characterized by negative anomalies of natural currents; the more filtration discharge, the higher these anomalies are. As a result, the relationship governing the intensity of streaming potential and sand granulometric composition, electrolyte concentrations and other factors were revealed. To determine the velocity of water flowing to leakage sites, a special device, based on the relationship between the temperature of a heated body and its resistance, is applied. This device simultaneously makes it possible to measure the water medium temperature. Observations at reservoir sites were effected by moving along the reservoir non-polarizable potential electrodes and water velocity devices. Recordings were carried out automatically by the recording device of the logging apparatus. Under the conditions of ice cover on water surfaces, measurements were made through separate points by digging holes in the ice cover. Practical field observations were carried out at reservoir sites located in regions where fissured massive rocks as well as loose sediments predominate. In the first case field experiments were carried out in alpine reservoirs, in Armenia. The major water leakages were found to be concentrated on the right bank of the reservoir. In this connection it was not only possible to locate water leakage sites, but also to evaluate their relative intensity. These data were used for planning antifiltration measures. In the second case water leakages from a reservoir located in Uzbekistan in the submontane part of the Pamirs were studied. Streaming potential anomalies and high benthonic flow rates made it possible to discover high filtrations in the base and walls of the dam. Further perfection of these methods should not only permit the determination of water leakage sites and their relative intensity, but also filtration discharges in absolute units.     

Optimization and simulation of reservoir operation with sediment evacuation: a case study of the Tarbela Dam,Pakistan

Many reservoirs around the world are being operated based on rule curves developed without considering the evacuation of deposited sediment. Current reservoir simulation and optimization models fall short of incorporating the concept of sustainability because the reservoir storage losses due to sedimentation are not considered. This study develops a new model called Reservoir Optimization‐Simulation with Sediment Evacuation (ROSSE) model. The model utilizes genetic algorithm based optimization capabilities and embeds the sediment evacuation module into the simulation module. The sediment evacuation module is implemented using the Tsinghua university flushing equation. The ROSSE model is applied to optimize the rule curves of Tarbela Reservoir, the largest reservoir in Pakistan with chronic sedimentation problems. In the present study, rule curves are optimized for maximization of net economic benefits from water released. The water released can be used for irrigation, power production, sediment evacuation, and for flood control purposes. Relative weights are used to combine the benefits from these conflicting water uses. Nine sets of rule curves are compared, namely existing rule curves and proposed rule curves for eight scenarios developed for various policy options. These optimized rule curves show an increase of net individual economic benefits ranging from 9 to 248% over the existing rule curves. The shortage of irrigation supply during the simulation period is reduced by 38% and reservoir sustainability is enhanced by 28% through increased sediment evacuation. The study concludes that by modifying the operating policy and rule curves, it is possible to enhance the reservoir's sustainability and maximize the net economic benefits. The developed methodology and the model can be used for optimization of rule curves of other reservoirs with sedimentation problems. Copyright © 2008 John Wiley & Sons, Ltd.     

浅地表地球物理技术在岩土工程中的应用与挑战

地球物理勘探方法在岩土工程的应用已有很长一段时间,但其成效与工程师的期待往往有不小的落差,以致于在一般的工程应用上仍不普遍.近年来浅地表地球物理技术有显着的进展,特别是在走时速度层析成像(Traveltime Tomography)、电阻率层析成像法(Electrical Resistivity Tomography)及多道瞬态面波法(Multichannel Analysis of Surface Wave).本文首先介绍这些方法在台湾岩土工程的应用,主要包括地层土壤液化潜能评估、坝体的安全检测、土壤与地下水污染调查及地基改良的质量管控等,应用案例以台湾常使用的地球物理勘探方法逐一介绍.虽然许多成功案例与新的应用方向对于浅地表地球物理技术在岩土工程应用的推广起了鼓舞作用,本文从工程师的角度提出地球物理勘探工程大量应用的挑战与瓶颈,包括如何提升探测数据的客观性、数据反演非唯一性问题、探测深度与分辨率的限制、实际条件违背反演基本假设的情况、以及地物性质与工程性质链接的不确定性问题,并进一步针对这些问题说明相关研究的进展与实务对策.希冀透过上述探讨,降低物探师与工程师认知上的差距,提升地球物理勘探在工程的应用的合理性与普及性.     

Numerical analysis of coastal cliff failure along the Pembrokeshire Coast National Park,Wales, UK

Cliff stability within the Pembrokeshire Coast National Park was evaluated using a numerical model applied at four sites representative of rock mass failure phenomena and major sedimentary geological sequences. The sites were: Mill Bay, (Old Red Sandstone); St Govan's Head, Carboniferous Limestone (Dinantian); Druidston, Millstone Grit (Namurian) and Lower Coal Measures (Westphalian); and Wiseman's Bridge, Lower Coal Measures (Westphalian). The study integrated a range of geotechnical parameters, measured in the field and laboratory, into a model to predict the likely failure mechanisms. The model is based on the existence of rock prisms delineated by structural parameters, i.e. joints, bedding planes and critical tension fractures behind the cliff face. An iterative approach is used to define the dip of the most probable, stepped failure surface at the base of any potentially unstable multiblock system and to calculate the sliding and toppling forces for each block in the cliff mass. Prediction compared favourably with field observations at three of the four selected sites, i.e. Druidston, St Govan's Head and Wiseman's Bridge. At Druidston sliding is predicted and dominates in the basal blocks, whilst toppling is confined to the upper cliff and is dependent on movement of the lower structural units. St Govan-s Head is shown to have a low risk of toppling and sliding and this was predicted except where basal undercutting reaches a depth of 1·0 m or lateral forces exceed 100 kN m⁻² when failure could occur. Copyright © 1998 John Wiley & Sons, Ltd.     

Separation of time scales in the HCA model for sand

Separation of time scales is used in a high cycle accumulation (HCA) model for sand. An important difficulty of the model is the limited applicability of the Miner’s rule to multiaxial cyclic loadings applied simultaneously or in a combination with monotonic loading. Another problem is the lack of simplified objective HCA formulas for geotechnical settlement problems. Possible solutions of these problems are discussed.     

PHYSICAL MODELING OF HYDRAULIC DESILTATION IN TAPU RESERVOIR

1 INTRODUCTIONIn recent years, the concePt of long-term sustained use of reservoirs has been addressed because areservoir is very much considered to be a nonrenewable resource (Morris and Fan, l998). Technically,many options for reservoir sedAnentation control can be utlized to pursue the sustainable develoPment ofwater resources. In general, reduction of incoming sedimen yields from watersheds is often emPloyedin conjunction with hydraulic methods such as flushing or density currnt vot…     

Recent applications of sliding block theory to geotechnical design

The sliding block theory was proposed by Newmark for determining the permanent displacement of embankments and dams under earthquake loading. This paper highlights recent applications of sliding block theory to different geotechnical structures. The equations to determine seismic factor of safety, yield acceleration and permanent displacement are given for rock block, soil slope, landfill cover, geosynthetic-reinforced soil retaining wall, and composite breakwater. The presented equations for seismic stability degenerate to that of static stability in the absence of earthquake. The permanent displacement for various structures can be obtained from that of a horizontal sliding block through a correction factor. A simplified procedure is included for the permanent displacement under vertical acceleration. The sliding block approach is rational for design under high seismic load.     

A data‐based mechanistic modelling (DBM) approach to understanding dynamic sediment transmission through Wyresdale Park Reservoir,Lancashire, UK

This paper outlines the application of a new data‐based mechanistic (DBM) modelling methodology to the characterization of the sediment transmission dynamics in a small upland reservoir, Wyresdale Park, Lancashire. The DBM modelling strategy exploits advanced statistical procedures to infer the dynamic model structure and its associated parameters directly from the instrumented data, producing a parametrically efficient, continuous time, transfer function model which relates suspended sediment load at the reservoir inflow to the outflow at the event scale. The associated differential equation model parameters have physical attributes which can be interpreted in terms of sediment transmission processes and associated reservoir trap efficiency. Sedigraph analysis suggests that wind‐induced resuspension episodically supplies an additional load to the reservoir outlet. The stochastic nature of the DBM model makes it ideal for evaluating the effects of uncertainty through Monte Carlo simulations (MCS) for discharge and sediment transmission. Copyright © 2000 John Wiley & Sons, Ltd.     

瞬态面波法在滑坡勘察中的应用效果

通过对瞬态面波法的工作原理、野外工作方法、资料整理和成果解释的描述和对应用实例的比较与分析,认为瞬态面波法在滑坡勘察中确定滑坡体的岩土分层、软弱带、滑动面是可行和可靠的。     

EM target location1

We consider the problem of computing the most probable location of a target based on radar measurements of the subsurface. Our algorithm makes use of the maximum likelihood estimator (MLE), which represents a correlation between the measured data and synthetic data generated for the object of interest at different locations. Previous studies assume a plane-wave acquisition geometry and target object(s) embedded in a uniform background. In this paper, a generalization of the MLE method is presented which is valid for discrete point sources (and receivers) and a 2D model (i.e. a 2.5D acquisition geometry). Within this formulation the treatment of a non-uniform background model is also possible. We concentrate on geotechnical ground investigations and assume that the characteristic dimensions of the target object are in the range 1–2λ, (λ being the wavelength). The potential of the method is demonstrated employing cross-hole radar data acquired in a controlled field experiment. The MLE result is also compared with the image obtained employing a full reconstruction method such as diffraction tomography.     

统计预测模型对三峡库区诱发地震预测的应用

通过实地考察并结合前人的研究成果,在地震地质背景研究的基础上,将长江三峡工程坝址至巫山库段划分为5个预测区共12个预测单元,选用8个诱震因素,运用统计预测模型预测了水库诱发地震的可能性及诱发地震震级。结果表明,库水蓄至最高水位,三峡工程坝址至巫山库段中多数预测单元诱发微震和不发生地震的概率值高,无震可能性较大。但九湾溪断层沿线、仙女山断层库段具备诱发3.0～4.5级或4.5～6.0级地震的可能;考虑到龙会观5.1级地震与高桥断裂的关系,预测结果显示高桥断裂沿线具有诱发强烈水库地震(≥6.0级)的可能性,其预测概率值为0.23