海湾水库蓄水初期底质与淡水盐分交换的试验研究

在现场调查和钻探的基础上,采集了沐官岛水库库区3种典型底质(粉质粘土、泥质粉沙和中细砂)的原状/扰动土样和地表水样,然后分别采用静水土柱和动水水槽试验测定盐分的时-空变化规律,最后定量计算出不同底质盐分的释放通量,从而为该海湾水库蓄水初期水质的评价和预测提供了科学依据.土柱试验表明,在分子扩散作用下高盐分区主要集中在水-沉积物界面之上7.5 cm的范围内,7.5 cm之上的水体盐分较为均一,底质盐分释放通量按粉质粘土、中细砂、泥质粉沙的顺序递减,盐分释放通量符合负的幂指数形式,而且抽排界面之上高浓度水体对降低水体盐分含量效果显著.根据水槽模拟试验,风的吹拂会影响到界面之上水体盐分的分层,有利于海湾水库中盐分的混合作用.     

海湾水库沉积物对水体咸化影响的时间尺度

为预测海湾水库沉积物盐分释放对库水咸化影响的持续时间,建立了描述沉积物孔隙水盐分剖面分布的非稳态数学模型,并优化选取模型的边界条件与关键参数(沉积物孔隙水盐分的扩散系数)。稳态与非稳态条件下沉积物盐分释放的对比分析表明,在计算污染沉积物对其上覆水水质的影响时间时,按非稳态释放更符合实际情况。非稳态盐分释放的计算结果表明,拟建沐官岛海湾水库沉积物中盐分释放对库水咸化的影响时间将大于600年,该影响时间远长于水库的寿命;表明在水库的整个生命周期内,均需要定期监测海湾水库底层水的盐分浓度,并需采取有效措施预防水体的突然泛咸,而不能仅在建库初期关注咸化问题。     

沐官岛水库蓄水初期盐分运移数值模拟

沐官岛水库是一个拟建的河口海湾水库,库底为富含咸水的潮滩沉积物。这一特殊的地质及水环境条件直接关系到水库未来的水质安全和正常的调度运行,所以必须准确地确定水库水体中盐分空间分布这一关键问题。本文依据现场调查、钻探及室内实验资料,通过垂直方向线性插值技术,建立了沐官岛水库蓄水初期盐分运移准三维模型,模拟了水库蓄水初期在内源盐分释放影响下,库水盐分的演化过程。结果表明:若遇平水年开始蓄水,当蓄水至1.5m时,该层含盐量在平面上的变化范围主要在0.30~1.30g/L之间;当蓄水至4.5m时,该层含盐量在平面上的变化范围主要在0.30~1.55g/L。     

水位变化条件下海湾水库盐分达标探讨

针对中国北方海湾水库间歇来水、连续取水和沉积物动态释盐的特点,建立水量与盐分耦合的数学模型,以青岛拟建的沐官岛水库为例,探讨水位变化条件下混合型海湾水库库水盐分的影响因素、超标风险与达标条件。模拟结果表明,在水库水位连续降低条件下,毫米量级的日蒸发量对库水盐分的累积效应显著。在不利水文条件下,受沉积物释盐、水分蒸发与人工取水的影响,混合型海湾水库长期存在盐分超标的风险。与水位不变时相反,水位降低时库水盐分浓度随着取水量的增大而升高;因此,当库水盐分存在超标风险时,可以通过减少日取水量实现库水盐分达标。为保障安全供水,在海湾水库设计与运营管理阶段,均需要综合考虑水量-盐分因素进行水库的日取水量调算。     

地下咸水与水库水体交换过程中沉积物胶体释放规律

以天津滨海地区北大港水库为研究对象,采用室内柱试验,研究地下咸水与水库水体交换过程中不同位置沉积物胶体释放以及盐分释放/截留的动态特征,同时对沉积物胶体释放、盐分释放/截留机理进行了探讨。研究结果表明:水库不同位置地下咸水与水库水体交换过程中,盐分的归宿不同:接近水库入口处的沉积物能将盐分截留下来,而出水口沉积物却将盐分释放转移到水体。随孔隙体积数的增加,沉积物胶体累计释放量逐渐增加,入库口、库中心、出库口最大累计释放量分别为3.275 mg/g、0.386 mg/g和1.382 mg/g;胶体累计释放量随孔隙体积数的变化曲线符合直线型,胶体释放速率变化很小。盐分的释放或截留是沉积物颗粒的粒径、胶体含量、含盐量等多种因素作用的结果,水库水体与沉积物中的盐分处于动态平衡状态,当沉积物中含盐量高于平衡浓度时,其盐分会向水体中释放,同时吸附在胶体上的盐分也会随着胶体的释放而释放;反之,水体中的盐分会向沉积物中迁移被截留下来,沉积物粒径越小,越易吸附水中的盐分。胶体的释放规律可以用双电层理论得到很好的解释。     

各向异性密度流模型及其验证

针对密度流流动各向异性的特性,采用非线性k-ε模式求取雷诺应力,结合RNG模型及浮力修正建立了三维密度流模型.应用该三维模型计算了水库温差异重流,分析了温差异重流演进过程,并就流速分布与水温下泄过程与实测资料进行了比较.结果表明,所建模型由于引入非线性RNG模式,能有效反映密度流的各向异性特性,较好地模拟了流场和温度场强烈耦合的水库温差异重流,能准确预测温度分层过程与异重流潜行时间,具有较高精度.     

基于微分对策理论的多沙河流水库调度

水库的排沙问题一直是多沙河流水库调度研究的重点之一。由于大部分多沙河流的来水来沙主要集中于汛期,因此水库汛期的蓄水兴利与泄水排沙的矛盾十分突出,具有典型的博弈关系,该关系可以利用博弈的相关理论进行描述。基于微分对策理论,将多沙河流水库汛期的供水兴利与泄水排沙看作博弈的双方,建立起以排沙和供水的综合效益作为性能指标函数的多沙河流水库汛期调度模型,研究了多沙河流水库汛期的调度问题,并通过陕西黑河某水库的实测资料对该模型进行了验证计算。计算结果表明,建立的调度模型能够较好地反映多沙河流水库汛期调度中的主要矛盾,实现对供水与排沙综合效益的优化。     

水库防洪优化调度的恒定出流模型及应用

现行水库洪水优化调度数学模型存在泄流闸门开度不断调整及时段间流量突变的问题。针对这些问题,将水库防洪优化调度数学模型分为时段内出库流量线性变化的瞬时出流模型和时段内出流不变的恒定出流模型。对于1个水库和1个防洪控制点所组成的基本防洪系统,应用矩形入流条件的河道洪水演进方法,以时段内恒定的出库流量为决策变量,构建水库防洪优化补偿调度数学模型。实例计算结果表明:恒定出流模型比瞬时出流模型占用的防洪库容减小0.05%~0.18%、最大下泄流量有增有减,即恒定出流模型既不劣于瞬时出流模型,也能很好地解决瞬时出流模型存在的问题。同时,实例揭示的水库二次补偿调节比等蓄量调度方式减小防洪库容10.6%,为确定水库的防洪库容提供了新的方法。     

水库蓄水对莫洛村滑坡的影响研究

莫洛村滑坡上分布有大量的古藏碉群,猴子岩水库蓄水后,可能影响位于库尾左岸莫洛村滑坡的稳定性,从而对藏碉群文化遗产构成威胁。本文分析了莫洛村滑坡的稳定性现状,预测蓄水可能对莫洛村滑坡稳定性的影响,指出即使采用1 852m初拟最高水位,水库蓄水对堆积体的稳定性影响也是有限的,不会对藏碉群文化遗产构成威胁,今后藏碉群文化遗产的保护应该放到防治莫洛村滑坡的局部解体和藏碉周围的水沟防渗上。     

反调节水库非恒定流数值模拟

采用二维水深平均方程组作为控制方程,利用边界拟合法进行坐标变换,采用有限体积法计算反调节水库的流场。在计算中,由于已知上游泵站抽水流量及下游枢纽供水流量,且上下游均找不到水位流量关系,上下游边界均为流速边界条件,没有压力边界条件,解不惟一,所以采用全水库水量平衡方程作为补充方程来补充定解条件。     

Contained density currents with high volume of release

Contained density currents with high volume of release reflect against the boundaries of the reception environment commonly leading to oscillatory flow. These flows exist in sediment clarifiers, compromising their operations, and deposited signatures of contained turbidity currents are found as part of the infill of sedimentary basins; for operation of the former and interpretation of the latter it is essential to understand the dynamic processes of these flows. Six lock‐exchange experiments with different initial densities were made in a horizontal flume, where the volume of the saline mixture in the lock was equivalent to the volume of the ambient fluid. A further two tests, with a repeated initial density, were made: one with high volume of release and very long duration; and another with low volume of release. Firstly, the movement of the current is discussed, including the oscillations within the experimental tank involving the density current and an upper layer counter‐current. It is shown that the cyclic behaviour is self‐similar with the reduced gravity of the initial density in the lock. Secondly, entrainment and water mixing processes are characterized. The time evolution of mixing is characterized qualitatively by analysing the background potential energy of the density distributions to show that mixing occurs even in the earlier stages of the current, and mainly within the first cycle of the oscillation. Quantified analysis reveals that, in currents with high volume of release, entrainment discharge is one order of magnitude higher, mainly due to the larger interface between the ambient fluid and the current. A model for the evolution of the mixing process is proposed for density currents with high volume of release. Finally, the dynamics of the head of the current is analysed. The entrainment in the head, when compared to the entrainment in the remainder body of the flows, is less important for the configuration with a larger lock.     

水库异重流一维水沙耦合模型

建立了基于库区不规则断面的一维非恒定异重流数学模型,并采用明流与异重流水沙输移模型交替运算的两步模式,即用潜入条件动态判别异重流计算的上游边界位置,将潜入点上游的明流浑水段与下游异重流段计算连接起来。水流运动、泥沙输移与河床变形过程完全耦合,采用TVD(Total Variation Diminishing)形式的MUSCL-Hancock格式进行数值求解。将该模型应用于恒定流量与释放定量悬沙两种条件下的异重流水槽实验模拟,比较了有无水面梯度项对模拟精度的影响,计算结果表明该模型能较为准确地预测异重流的厚度、含沙量分布及传播过程。     

Characteristics of sediment in transport in the swash zone of a steep estuarine foreshore

Cross‐shore grading of sediment has been observed on the surface of estuarine beaches but the swash zone processes responsible for this grading have not been measured. This study was conducted to provide an explanation for the cross‐shore grading of sediment on a predominantly sandy estuarine foreshore. Data on wave and swash characteristics and sediment trapped in the uprush and backwash during 25 swash events were gathered from mid‐rising to mid‐falling tide on a small transgressive barrier in Delaware Bay, New Jersey, USA. Sediment is predominantly quartz and feldspar, medium to coarse sands with a gravel fraction of granules and pebbles. Wave energies increased with tidal rise. The percentage of gravel in transport in the uprush and backwash is similar (11% and 13%) during the rising tide when the swash zone is at mid‐foreshore, decreases in the uprush (9%) and increases in the backwash (18%) when the swash zone is on the upper foreshore. When the swash zone is at mid‐foreshore on the falling tide, the quantity of gravel in the backwash (30%) is greater than in the uprush (24%). The low proportion of gravel within the foreshore prior to trapping, and the increase in the percentage of gravel when the waves and swash are on the upper foreshore, suggests that the step is the primary source of gravel high on the foreshore. The size of the step increases as wave heights increase with tidal rise. The rate of delivery of gravel into the swash is enhanced by sediment entrained during wave breaking and interaction of the uprush with the previous backwash. The lag in the rate of step migration relative to breaker migration during the falling tide increases the likelihood of mining gravel from the step and subsequent transport in the uprush and backwash. These findings are important for low energy estuarine beaches sensitive to small changes in tidal range and wave energy that cause sedimentological change across the foreshore.     

Bedload transport and bed resistance associated with density and turbidity currents

Turbidity currents in the ocean are driven by suspended sediment. Yet results from surveys of the modern sea floor and turbidite outcrops indicate that they are capable of transporting as bedload and depositing particles as coarse as cobble sizes. While bedload cannot drive turbidity currents, it can strongly influence the nature of the deposits they emplace. This paper reports on the first set of experiments which focus on bedload transport of granular material by density underflows. These underflows include saline density flows, hybrid saline/turbidity currents and a pure turbidity current. The use of dissolved salt is a surrogate for suspended mud which is so fine that it does not settle out readily. Thus, all the currents can be considered to be model turbidity currents. The data cover four bed conditions: plane bed, dunes, upstream‐migrating antidunes and downstream‐migrating antidunes. The bedload transport relation obtained from the data is very similar to those obtained for open‐channel flows and, in fact, is fitted well by an existing relation determined for open‐channel flows. In the case of dunes and downstream‐migrating antidunes, for which flow separation on the lee sides was observed, form drag falls in a range that is similar to that due to dunes in sand‐bed rivers. This form drag can be removed from the total bed shear stress using an existing relation developed for rivers. Once this form drag is subtracted, the bedload data for these cases collapse to follow the same relation as for plane beds and upstream‐migrating antidunes, for which no flow separation was observed. A relation for flow resistance developed for open‐channel flows agrees well with the data when adapted to density underflows. Comparison of the data with a regime diagram for field‐scale sand‐bed rivers at bankfull flow and field‐scale measurements of turbidity currents at Monterey Submarine Canyon, together with Shields number and densimetric Froude number similarity analyses, provide strong evidence that the experimental relations apply at field scale as well.     

Sedimentation and welding processes of dilute pyroclastic density currents and fallout during a large-scale silicic eruption, Kikai caldera, Japan

Sedimentation and welding processes of the high temperature dilute pyroclastic density currents and fallout erupted at 7.3 ka from the Kikai caldera are discussed based on the stratigraphy, texture, lithofacies characteristics, and components of the resulting deposits. The welded eruptive deposits, Unit B, were produced during the column collapse phase, following a large plinian eruption and preceding an ignimbrite eruption, and can be divided into two subunits, Units B_l and B_u. Unit B_l is primarily deposited in topographic depressions on proximal islands, and consists of multiple thin (< 1 m) flow units with stratified and cross-stratified facies with various degrees of welding. Each thin unit appears as a single aggradational unit, composed of a lower lithic-rich layer or pod and an upper welded pumice-rich layer. Lithic-rich parts are fines-depleted and are composed of altered country rock, fresh andesite lava, obsidian clasts with chilled margins, and boulders. The overlying Unit B_u shows densely welded stratified facies, composed of alternating lithic-rich and pumice-rich layers. The layers mantle lower units and are sometimes viscously deformed by ballistics. The sedimentary characteristics of Unit B_l such as welded stratified or cross-stratified facies indicate that high temperature dilute pyroclastic density currents were repeatedly generated from limited magma-water interactions. It is thought that dense brittle particles were segregated in a turbulent current and were immediately buried by deposition of hot, lighter pumice-rich particles, and that this process repeated many times. It is also suggested that the depositional temperature of eruptive materials was high and the eruptive style changed from a normal plinian eruption, through surge-generating explosions (Unit B_l), into an agglutinate-dominated fallout eruption (Unit B_u). On the basis of field data, welded pyroclastic surge deposits could be produced only under specific conditions, such as (1) rapid accumulation of pyroclastic particles sufficiently hot to weld instantaneously upon deposition, and (2) elastic particles' interactions with substrate deformation. These physical conditions may be achieved within high temperature and highly energetic pyroclastic density currents produced by large-scale explosive eruptions.     

An experimental investigation of density-stratified inertial gravity currents

Turbidity currents and pyroclastic density currents may originate as stratified flows or develop stratification during propagation. Analogue, density‐stratified laboratory currents are described, using layers of salt solutions with different concentrations and depths to create the initial vertical stratification. The evolving structure of the flow depends on the distribution of the driving buoyancy between the layers, B* (proportional to the layer volumes and densities), and their density ratio, ρ*. When the lower layer contains more salt than the upper layer, and so has a greater proportion of the driving buoyancy (B* < 0·5), this layer can run ahead leading to streamwise or longitudinal stratification (ρ*→0), or the layers can mix to produce a homogeneous current (ρ*→1). If the upper layer contains more salt and thus buoyancy (B* > 0·5), this layer travels to the nose of the current by mixing into the back of the head along the body/wake density interface to produce a homogeneous flow (ρ*→1) or overtaking, leading to streamwise stratification (ρ*→0). Timescales describing the mixing between the layers and the streamwise separation of the layers are used to understand these flow behaviours and are in accordance with the experimental observations. Distance–time measurements of the flow front show that strongly stratified flows initially travel faster than weakly stratified flows but, during their later stages, they travel more slowly. In natural flows that are stratified in concentration and grain size, internal features, such as stepwise grading, gradual upward fining and reverse grading, could be produced depending on B* and ρ*. Stratification may also be expected to affect interactions with topography and overall fan architecture.     

Stacked fluvial and tide-dominated estuarine deposits in high-frequency (fourth-order) sequences of the Eocene Central Basin, Spitsbergen

Eighteen coastal-plain depositional sequences that can be correlated to shallow- to deep-water clinoforms in the Eocene Central Basin of Spitsbergen were studied in 1 × 15 km scale mountainside exposures. The overall mud-prone (>300 m thick) coastal-plain succession is divided by prominent fluvial erosion surfaces into vertically stacked depositional sequences, 7–44 m thick. The erosion surfaces are overlain by fluvial conglomerates and coarse-grained sandstones. The fluvial deposits show tidal influence at their seaward ends. The fluvial deposits pass upwards into macrotidal tide-dominated estuarine deposits, with coarse-grained river-dominated facies followed further seawards by high- and low-sinuosity tidal channels, upper-flow-regime tidal flats, and tidal sand bar facies associations. Laterally, marginal sandy to muddy tidal flat and marsh deposits occur. The fluvial/estuarine sequences are interpreted as having accumulated as a series of incised valley fills because: (i) the basal fluvial erosion surfaces, with at least 16 m of local erosional relief, are regional incisions; (ii) the basal fluvial deposits exhibit a significant basinward facies shift; (iii) the regional erosion surfaces can be correlated with rooted horizons in the interfluve areas; and (iv) the estuarine deposits onlap the valley walls in a landward direction. The coastal-plain deposits represent the topset to clinoforms that formed during progradational infilling of the Eocene Central Basin. Despite large-scale progradation, the sequences are volumetrically dominated by lowstand fluvial deposits and especially by transgressive estuarine deposits. The transgressive deposits are overlain by highstand units in only about 30% of the sequences. The depositional system remained an estuary even during highstand conditions, as evidenced by the continued bedload convergence in the inner-estuarine tidal channels.     

Analogue study of particle segregation in pyroclastic density currents, with implications for the emplacement mechanisms of large ignimbrites

Abstract Analogue flume experiments were conducted to investigate the transport and sedimentation behaviour of turbulent pyroclastic density currents. The experimental currents were scaled approximately to the natural environment in three ways: (1) they were fully turbulent; (2) they had a very wide range of particle sizes and associated Rouse numbers (the ratio of particle settling velocity to effective turbulent eddy velocity in the current); and (3) they contained particles of two different densities. Two sets of surge‐type experiments were conducted in a 5 m long, water‐filled lock‐exchange flume at five different volumetric particle concentrations from 0·6% to 23%. In one set (one‐component experiments), the currents contained just dense particles; in the other set (two‐component experiments), they contained both light and dense particles in equal volume proportions. In both sets of experiments, the population of each component had a log‐normal size distribution. In the two‐component experiments, the size range of the light particle population was selected in order to be in hydrodynamic equivalence with that of the dense particles. Dense particles were normally graded, both vertically and downstream, in the deposits from both sets of experiments. The mass loading (normalized to the initial mass of the suspension) and grain size of the dense component in the deposits decreased with distance from the reservoir and were insensitive to initial total particle concentration in the currents. On the other hand, in the two‐component experiments, the light particles were extremely sensitive to concentration. They were deposited in hydrodynamic equivalence with the dense particles from dilute currents, but were segregated efficiently at concentrations higher than a few per cent. With increasing particle concentration, the large, light particles were carried progressively further down the flume because of buoyancy effects. Deposits from the high‐concentration currents exhibited reverse vertical grading of the large, light particles. Efficient segregation of the light component was observed even if the bulk density of the current was less than that of the light particles. In both sets of experiments, marked inflexions in the rate of downstream decline in mass loading and maximum grain size of the dense component can be attributed to the presence of two different particle settling regimes in the flow: (1) particles with Rouse numbers >2·5, which did not respond to the turbulence and settled rapidly; and (2) particles with Rouse numbers <2·5, which followed the turbulent eddies and settled slowly. The results are applied to the transport and sedimentation dynamics of pyroclastic density currents that generate large, widespread ignimbrites. Field data fail to reveal significant departures from aerodynamic equivalence between pumice and lithic clasts in three such ignimbrites: the particulate loads of some large ignimbrites are transported principally in turbulent suspensions of low concentration. In some ignimbrites, the well‐developed inflexions in curves of maximum lithic (ML) size vs. distance can be attributed to the existence of distinct high and low Rouse number particle settling regimes that mark the transition from an overcharged state to one in which the residual particulate load is transported more effectively by turbulence.     

A GIS-based methodology for hazard mapping of small volume pyroclastic density currents

We present here a methodology implemented within a geographical information system (GIS) for hazard mapping of small volume pyroclastic density currents (PDCs). This technique is implemented as a set of macros written in Visual Basic for Applications (VBA) that run within GIS-software (i.e. ArcGIS). Based on the energy line concept, we calibrated an equation that relates the volume (V) and the mobility (ΔH/L) of single PDCs using data from Soufrière Hills volcano (Montserrat) and Arenal volcano (Costa Rica). Maximum potential run-outs can be predicted with an associated uncertainty of about 30%. Also based on the energy line concept and with data from Soufrière Hills volcano and Mt. St. Helens (USA), we were able to calibrate an equation that predicts the flow velocity as a function of the vertical distance between the energy line and the ground surface (Δh). Velocities derived in this way have an associated uncertainty of 3 m s⁻¹. We wrote code to implement these equations and allow the automatic mapping of run-out and velocity with the inputs being (i) the height and location of the vent (ii) the flow volume and (iii) a digital elevation model (DEM) of the volcano. Dynamic pressure can also be estimated and mapped by incorporating the density of the pyroclastic density current (PDC). This computer application allows the incorporation of uncertainties in the location of the vent and of statistical uncertainties expressed by the 95% confidence limits of the regression model. We were able to verify predictions by the proposed methodology with data from Unzen volcano (Japan) and Mayon volcano (The Philippines). The consistencies observed highlight the applicability of this approach for hazard mitigation and real-time emergency management.