Re-Look to Conventional Techniques for Trapping Efficiency Estimation of a Reservoir

All reservoirs are subjected to sediment inflow and deposition up to a certain extent leading to reduction in their capacity. Thus, the important practical problem related to the life of reservoir is the estimation of sedimentation quantity in the reservoirs. Large number of methods and models are available for estimation of reservoir sedimentation process. However, each model differs greatly in terms of their complexity, inputs and other requirements. In the simplest way, the fraction of sediment deposit in the reservoir can be determined through the knowledge of its trap efficiency. Trap efficiency （Te） is the proportion of the incoming sediment that is deposited or trapped in a reservoir. Most of the Te estimation methods define a relationship of the T, of the reservoir to their capacity and annual inflow, generally through curves. In this study, the empirical relationships given by Brune and Brown were used and compared for estimating the trap efficiency of Gobindsagar Reservoir （Bhakra Dam） on Satluj River in Bilaspur district of Himachal Pradesh, in the Himalayan region of India. A new set of regression equations has been developed for Brune＇s method and compared with Brown and other available Brune＇s equations. It has been found that Brune＇s equations developed in the present study estimated better than the other Brune＇s equations reported in literature. Later, in the present study it was found that Brown＇s approach was over estimating the T,. Hence it was again modified for Gobindsagar reservoir. It was also identified that sediments coming to this particular reservoir were mainly of coarse nature.     

NUMERICAL SIMULATION OF TURBIDITY CURRENT IN RESERVOIR

I. INTRODUCTIONWhen a sediment--laden flow reaches the backwater zone of a reservoir, the suddenreduction of flow velocity causes sediment particles to settle toward the river bed. Undercertain circumstsnces, it will plunge and form a layer of sediment--water mixture flowingbeneath the water surface. This flowing layer is called the turbidity current. A turbiditycurrent is relatively stable and has important impacts on reservoir sedimentation.In the case of deep reservoirs, due to temper…     

One-dimensional morphodynamic model for retrogressive erosion based on a sediment entrainment theory at high flow velocity

Retrogressive erosion is a high-speed erosion process that usually occurs during the rapid release of stored water in reservoirs built on sandy rivers.Retrogressive erosion has been utilized in the practice of reservoir sedimentation control,but accurate prediction of the bed deformation process by numerical models has rarely been reported.The current study presents a one-dimensional morphodynamic model for simulating the evolution process of retrogressive erosion induced by high-velocity flows on steep slopes.The governing equations apply a Cartesian coordinate system with a vertically oriented z axis.The bed surface gradient and friction terms in the flow equations include correction factors to take account of the effects of high slope on flow movement.The net vertical sediment flux term in the sediment transport and bed deformation equations is calculated using an equation of erosion velocity.Particularly,this equation is based on an empirical relation between the sediment entrainment rate and the Shields parameter in contrast to the traditional sediment transport capacity,and the critical Shields parameter is modified by taking into account the permeability of the sediment layer and the stability of particles on a slope.The feedback of scoured sediment on the flow movement is considered by additional terms in the governing equations.Flume experiments of retrogressive erosion in literature were simulated to validate the model.The temporal variations of the longitudinal profiles of the free surface and channel bed and the sediment transport rate were well predicted.The algorithm calculating sediment entrainment in the proposed model also was validated for an experiment measuring entrainment rate from the literature.More importantly,it was found that the morphodynamic model using the sediment transport capacity equation predicts the trend of cumulative erosion contrary to the measurements,while results of the proposed model can follow a similar trend with the observed data in the retrogressive erosion process.     

Induced seismicity at wujiangdu reservoir,China: a case induced in the karst area

To date 19 cases of reservoir-induced seismicity have been acknowledged in China and 15 of them are associated with karst. The Wujiangdu case is a typical one induced in the karst area. The dam with a height of 165 m is the highest built in a karst area in China. Seismic activity has been successively induced in five reservoir segments seven months after the impoundment in 1979. A temporary seismic network consisting of 8 stations was set up in one of the segments some 40 km upstream from the dam. The results indicate that epicenters were distributed along the immediate banks, composed of karstified carbonate, and focal depths were only several hundred meters. Most of the focal mechanisms were of thrust and normal faulting. It is suggested that karst may be an important factor in inducing seismicity. It can provide an hydraulic connection to change the saturation and pressure and also weak planes for dislocation to induce seismicity.     

Stochastic partial differential equations in groundwater hydrology

Many problems in hydraulics and hydrology are described by linear, time dependent partial differential equations, linearity being, of course, an assumption based on necessity.Solutions to such equations have been obtained in the past based purely on deterministic consideration. The derivation of such a solution requires that the initial conditions, the boundary conditions, and the parameters contained within the equations be stipulated in exact terms. It is obvious that the solution so derived is a function of these specified, values.There are at least four ways in which randomness enters the problem. i) the random initial value problem; ii) the random boundary value problem; iii) the random forcing problem when the non-homogeneous part becomes random and iv) the random parameter problem.Such randomness is inherent in the environment surrounding the system, the environment being endowed with a large number of degrees of freedom.This paper considers the problem of groundwater flow in a phreatic aquifer fed by rainfall. The goveming equations are linear second order partial differential equations. Explicit form solutions to this randomly forced equation have been derived in well defined regular boundaries. The paper also provides a derivation of low order moment equations. It contains a discussion on the parameter estimation problem for stochastic partial differential equations.     

于桥水库风季流场的计算与分析

本文依据水流运动的基本理论,并考虑风对水体的剪切作用,应用改进的移步ADI法,建立了风作用下流场的数值模似系统,其结果与连续四年的实测资料验证相吻合。计算预报了于桥水库三个典型水文年定常风不同风速、风向的流场,并由计算机绘制出彩色流场图。文中还针对计算应用成果进行了分析。     

于桥水库风季流埸的计算与分析

本文依据水流运动的基本理论,并考虑风对水体的剪切作用,应用改进的移步ADI法,建立了风作用下流场的数值模似系统,其结果与连续四年的实测资料验证相吻合。计算预报了于桥水库三个典型水文年定常风不同风速、风向的流场,并由计算机绘制出彩色流场图。文中还针对计算应用成果进行了分析。     

龙羊峡水库地震分形特征

着重研究了龙羊峡水库地震的时间频次和时间强度分形的ｆ（ａ）－α谱，以及Ｄ＋∞、Ｄ－∞的变化特点。结果认为龙羊峡水库地震对周围强震的震前异常反应很好。需要指出的是，在周围强震发生前存在较为明显的ｆ（α）－α谱的收缩现象Ｄｑ－ｑ曲线也趋于平缓。这种由于水库蓄水而诱发地震并对周围强震有良好反应的现象应引起一定的重视。     

储集层构造裂缝定量预测方法研究

选用裂缝体积密度作为构造裂缝定量化参数,对基于钻井岩心的裂缝体积密度的计算方法作了改进。统计分析得出,破裂试件裂缝体积密度与试件破裂应变能密度的关系受试件所受的围压变化及试件物性的控制,钻井岩心构造裂缝体积密度与数值模拟应变能密度的关系同试件研究结果有可比之处。在上述研究基础上,建立了利用钻井岩心裂缝参数和数值模拟应变能密度来定量预测碳酸盐岩储层构造裂缝体积密度的方法,并已应用于油田生     

Influence of reservoir geometry and conditions on the seismic response of arch dams

An analysis of the influence that reservoir levels and bottom sediment properties (especially on the degree of saturation) have on the dynamic response of arch dams is carried out. For this purpose, a Boundary Element Model developed by the authors that allows the direct dynamic study of problems that incorporate scalar (dammed up water), viscoelastic (dam and soil site) and poroelastic media (bottom sediments in the reservoir) is used. All of the regions are discretized using boundary elements, later formulating equations of compatibility and equilibrium that allow their interaction to be rigorously established. The seismic excitation consists in plane longitudinal waves (P waves) and shear waves (S waves) impinging the dam site with an angle of vertical incidence. The analysis is carried out in the frequency domain, and the time response is obtained, for synthesized artificial accelerograms defined in terms of the elastic response spectrum taken from Eurocode 8, using a FFT algorithm. The variables used to characterize the response are: Amplitude of the complex-valued frequency-response function, acceleration response spectra and the integral of velocity of points located at the structure. These variables clearly indicate the importance that the factors analyzed have on the dynamic response.     

Linear mechanism by which internal gravity waves are generated and intensified in the ionosphere when interacting with an inhomogeneous zonal wind: 2. Internal gravity wave generation and intensification during the linear evolution stage

The linear mechanism by which internal gravity waves (IGWs) are generated and subsequently intensified in a stably stratified dissipative ionosphere in the presence of an inhomogeneous zonal wind (shear flow) has been studied. In the case of shear flows, the operators of linear problems are nonself-adjoint and the corresponding eigenfunctions are nonorthogonal; a canonical approach can hardly be used to study such motions. It is more adequate to apply the so-called nonmodal calculation. Dynamic equations and equations of energy transfer of IGW disturbances in the ionosphere with a shear flow have been obtained based on a nonmodal approach. Exact analytical solutions for the constructed dynamic equations have been found. The growth rate of the IGW shear instability has been determined. It has been established that IGW disturbances are intensified in an algebraically power manner rather than exponentially in the course of time. The effectiveness of the linear mechanism by which IGWs are intensified when interacting with an inhomogeneous zonal wind is analyzed. It has been indicated that IGWs effectively obtain the shear flow energy during the linear evolution stage and substantially increase (by an order of magnitude) their amplitude and energy. The frequency and the wave vector of generated IGW modes depend on time; therefore, a wide spectrum of wavelike disturbances, depending on the linear, rather than nonlinear, turbulent effects, is formed in the ionosphere with a shear flow. Thereby, a new degree of freedom appears, and the turbulent state of atmospheric—ionospheric layers can be formed on IGW disturbances.     

基于逆建模理论的储层特征定量预测方法

定量描述储层特征的物理量称为储层参数,在储层描述工作中储层参数预测是一个重要环节.本文以岩石物理理论为基础研究了利用岩石物理逆建模理论进行储层参数预测的方法.在合理、有效的岩石物理模型的基础上,逆建模方法通过弹性参数的等值面的空间交会可以预测出储层参数.在此基础上提出了一种扩展的岩石物理逆建模方法,在两个不同的储层参数域进行三维岩石物理逆建模,通过逆建模结果的参数域转换可以预测更多的储层参数,为储层特征描述提供了更加充分的数据支持.输入数据的种类对于预测结果的准确性有较大影响,利用本文建立的弹性参数适用性分析方法可以选择出适用性好的弹性参数组合作为输入数据.将本方法应用于模型数据和实际工区数据,取得了较好的预测效果,证明方法具有一定的实用性.     

Markov chain Monte Carlo algorithms for target-oriented and interval-oriented amplitude versus angle inversions with non-parametric priors and non-linear forward modellings

In geophysical inverse problems, the posterior model can be analytically assessed only in case of linear forward operators, Gaussian, Gaussian mixture, or generalized Gaussian prior models, continuous model properties, and Gaussian-distributed noise contaminating the observed data. For this reason, one of the major challenges of seismic inversion is to derive reliable uncertainty appraisals in cases of complex prior models, non-linear forward operators and mixed discrete-continuous model parameters. We present two amplitude versus angle inversion strategies for the joint estimation of elastic properties and litho-fluid facies from pre-stack seismic data in case of non-parametric mixture prior distributions and non-linear forward modellings. The first strategy is a two-dimensional target-oriented inversion that inverts the amplitude versus angle responses of the target reflections by adopting the single-interface full Zoeppritz equations. The second is an interval-oriented approach that inverts the pre-stack seismic responses along a given time interval using a one-dimensional convolutional forward modelling still based on the Zoeppritz equations. In both approaches, the model vector includes the facies sequence and the elastic properties of P-wave velocity, S-wave velocity and density. The distribution of the elastic properties at each common-mid-point location (for the target-oriented approach) or at each time-sample position (for the time-interval approach) is assumed to be multimodal with as many modes as the number of litho-fluid facies considered. In this context, an analytical expression of the posterior model is no more available. For this reason, we adopt a Markov chain Monte Carlo algorithm to numerically evaluate the posterior uncertainties. With the aim of speeding up the convergence of the probabilistic sampling, we adopt a specific recipe that includes multiple chains, a parallel tempering strategy, a delayed rejection updating scheme and hybridizes the standard Metropolis–Hasting algorithm with the more advanced differential evolution Markov chain method. For the lack of available field seismic data, we validate the two implemented algorithms by inverting synthetic seismic data derived on the basis of realistic subsurface models and actual well log data. The two approaches are also benchmarked against two analytical inversion approaches that assume Gaussian-mixture-distributed elastic parameters. The final predictions and the convergence analysis of the two implemented methods proved that our approaches retrieve reliable estimations and accurate uncertainties quantifications with a reasonable computational effort.     

Modeling uncertainties of reservoir flushing simulations

Every dam or barrage construction affects the watercourse and the retention of sediment that previously was carried by the river, which can lead to siltation of the reservoir and obstruction of water intakes over time, reducing their capacities. However, the information available regarding the effect of sediment and drawdown parameters, sediment management at reservoirs, as well as different equational approaches, is scarce. The current research aims to evaluate the effect of parameters associated with numerical modeling of sediment management in reservoirs considering scenarios with different drawdowns, transport equations, sediment size distributions, and thickness of the initial sediment layer. The case study of the Aimorés Hydropower Plant (HPP) is used, applying the Delft3D-FLOW model for two-dimensional modeling. All parameters influenced the volume of mobilized sediment, among which the initial layer thickness was the parameter that resulted in the greatest changes in simulated results. In general, the results show that the uncertainties in the input parameters outweigh the uncertainties between the techniques, which found large variations in results when evaluating the use of different transport equations. These results indicate the importance of proper estimation of model parameters for predicting effects with accuracy and the need for such studies before planning and management operations are evaluated to avoid environmental harm and energy waste.     

Effect of surface stress on magneto-elastic surface waves in finitely conducting media

This paper deals with the investigation of the effect of surface stress and conductivity on the propagation of surface wave in isotropic, homogeneous, elastic media under the action of a primary magnetic field. Formulation of the general surface wave propagation problem has been made, and the corresponding frequency equation has been derived. Frequency equations for Rayleigh wave, surface shear wave and Stoneley wave have been deduced from that of general surface wave as special cases. The effects of surface stress and magnetic field on the wave velocities and attenuation factors of Rayleigh wave and surface shear wave are shown by numerical calculations and graphs. Some important wave velocity equations, as obtained by other authors, have been deduced as special cases from the wave velocity equation for Stoneley wave. It is found that the combined effect of surface stress and magnetic field modulates the wave velocity ratios and attenuation factors of Rayleigh wave and surface shear wave to a considerable extent.     

储层弹性与物性参数地震叠前同步反演的确定性优化方法

储层弹性与物性参数可直接应用于储层岩性预测和流体识别,是储层综合评价和油气藏精细描述的基本要素之一.现有的储层弹性与物性参数地震同步反演方法大都基于Gassmann方程,使用地震叠前数据,通过随机优化方法反演储层弹性与物性参数;或基于Wyllie方程,使用地震叠后数据,通过确定性优化方法反演储层弹性与物性参数.本文提出一种基于Gassmann方程、通过确定性优化方法开展储层弹性和物性参数地震叠前反演的方法,该方法利用Gassmann方程建立储层物性参数与叠前地震观测数据之间的联系,在贝叶斯反演框架下以储层弹性与物性参数的联合后验概率为目标函数,通过将目标函数的梯度用泰勒公式展开得到储层弹性与物性参数联合的方程组,其中储层弹性参数对物性参数的梯度用差分形式表示,最后通过共轭梯度算法迭代求解得到储层弹性与物性参数的最优解.理论试算与实际资料反演结果证明了方法的可行性.     

Radio signal propagation effects in a nonstationary ionosphere based on the numerical solution of the one-dimensional wave equation

The goal of this work is to solve Maxwell equations analytically and numerically in a one-dimensional case under the conditions of a nonstationary medium. Analytical solutions to the Maxwell equations have been obtained in two partial cases of the linear and quadratic time dependence of medium permittivity. Since the number of models for which the wave equation can be solved analytically is limited, it becomes also necessary to apply numerical methods, specifically the method of finite differences, in a time domain Finite Difference Time Domain method. The effects of the decameter wave dynamic reflection from structures with considerable spatial gradients (the scales of which are comparable with the sounding pulse wavelength) have been studied based on this method. It has been shown that the spectrum can broaden and a Doppler frequency shift of a reflected signal can originate can take place.     

黄河尾闾段一维耦合水沙数学模型研究及其应用

小浪底水库于1999年运用以后,该河道经历了长时间持续冲刷过程.为掌握小浪底水库运用后黄河尾闾段洪水演进特点及河床冲淤规律,采用一维水沙数学模型研究是一条重要的途径.本研究首先采用浑水控制方程,建立了一维耦合水沙数学模型,并利用2003年利津-西河口段汛期实测水沙及汛前断面地形资料对该模型进行率定,计算的流量、水位及含沙量等过程与实测值吻合较好;然后采用2015年利津—汊3段汛期实测资料对该模型进行验证,结果显示水位与冲淤量计算值与实测值较为符合;最后基于2015年实测洪水过程,计算了若干组不同断面间距下的洪水演进及冲淤过程,分析了不同断面间距对沿程水位及河段冲淤量等计算结果的影响,结果表明:采用不同断面间距对水位计算结果影响较小,而对冲淤量计算结果会产生一定影响;在河段水沙及冲淤特性复杂的情况下,采用一维数学水沙模型计算时应考虑断面间距的选择.     

The historically adjusted range and the historically rescaled adjusted range

It was remarked by Hurst in 1951 that the adjusted range gives the size of the smallest reservoir capable of providing a constant discharge equal to the mean inflow. Since that time this range and its rescaled modification, the Hurst range, have been widely discussed, not however primarily with a view to applying them to reservoir design problems, but rather on account of their possible relevance to the simulation of geophysical time series.Acknowledging the well-known conceptual weaknesses of adjusted ranges and the theoretical difficulties that inhibit their direct utilisation in the design and operation of real reservoirs, the authors argue that the interest displayed on ranges during the past few decades justifies the effort of eliminating one in particular of these weakness, namely their non-implementability as operating policies, a consequence of the fact that they can only be retrospectively evaluated. The paper proposes modifications in which the unknowable mean and standard deviation of future samples are replaced by the known mean and sample standard deviation of historical data, leading to the historically adjusted range and the historically rescaled and adjusted range. The latter is produced as an implementable approximation to Hurst's (1951) solution to the optimal reservoir problem.The expected values of the new ranges are evaluated and numerically tabulated.     

CHARACTERISTICS OF SEDIMENT TRANSPORT ALONG A RIVER REACH WITH A RESERVOIR

Based on long-term measurements at three gauging stations, Toudaoguai, Fugu and Hequ, and one meteorological station, this article discusses the features of discharge (Q) and sediment concentration (Cs) of a river reach of the Yellow River with a reservoir located in the Loess Plateau. The impacts of the local sub-watershed between Toudaoguai and Fugu gauging stations on sediment budget to the Yellow River have been analyzed. In addition, the deposition processes in the Tianqiao Reservoir have been investigated. Results show over 80% of the precipitation that falls in the local subwatershed is unable to contribute to the Yellow River runoff process. It is found that the annualmaximum sediment concentration is usually less than 30 kg/m^3 during flood seasons at Toudaoguai Gauging Station, but the sediment concentration varies dramatically at Fugu Gauging Station. About 35% of the sediment eroded in the sub-watersheds between Toudaoguai and Fugu gauging stationswas produced from the Huangfuchuan sub-watershed which has a drainage area accounting only for 10% of the drainage area between Toudaoguai and Fugu gauging stations. The Tianqiao Reservoir generally has deposition during the summer flood season, and scouring during the non-flood season.On average, over 85% of deposited sediment in the reservoir occurs in the 12 km long lower reservoir reach. The volume of annual deposition in the reservoir mainly depends on the volume of water from the local region between Hequ and Fugu gauging stations.