Stochastic and analytical approaches for sediment accumulation in river reservoirs

ABSTRACT

Sediment accumulation in a river reservoir is studied by stochastic time series models and analytical approach. The first-order moving average process is found the best for the suspended sediment discharge time series of the Juniata River at Newport, Pennsylvania, USA. Synthetic suspended sediment discharges are first generated with the chosen model after which analytical expressions are derived for the expected value and variance of sediment accumulation in the reservoir. The expected value and variance of the volume of sediment accumulation in the reservoir are calculated from a thousand synthetic time series each 38 years long and compared to the analytical approach. Stochastic and analytical approaches perfectly trace the observation in terms of the expected value and variability. Therefore, it is concluded that the expected value and variance of sediment accumulation in a reservoir could be estimated by analytical expressions without the cost of synthetic data generation mechanisms.     

Poynting and polarization vectors based wavefield decomposition and their application on elastic reverse time migration in 2D transversely isotropic media

With the progress in computational power and seismic acquisition, elastic reverse time migration is becoming increasingly feasible and helpful in characterizing the physical properties of subsurface structures. To achieve high-resolution seismic imaging using elastic reverse time migration, it is necessary to separate the compressional (P-wave) and shear (S-wave) waves for both isotropic and anisotropic media. In elastic isotropic media, the conventional method for wave-mode separation is to use the divergence and curl operators. However, in anisotropic media, the polarization direction of P waves is not exactly parallel to the direction of wave propagation. Also, the polarization direction of S-waves is not totally perpendicular to the direction of wave propagation. For this reason, the conventional divergence and curl operators show poor performance in anisotropic media. Moreover, conventional methods only perform well in the space domain of regular grids, and they are not suitable for elastic numerical simulation algorithms based on non-regular grids. Besides, these methods distort the original wavefield by taking spatial derivatives. In this case, a new anisotropic wave-mode separation scheme is developed using Poynting vectors. This scheme can be performed in the angle domain by constructing the relationship between group and polarization angles of different wave modes. Also, it is performed pointwise, independent of adjacent space points, suitable for parallel computing. Moreover, there is no need to correct the changes in phase and amplitude caused by the derivative operators. By using this scheme, the anisotropic elastic reverse time migration is more efficiently performed on the unstructured mesh. The effectiveness of our scheme is verified by several numerical examples.     

Conditional spectrum bidirectional record selection for risk assessment of 3D structures using scalar and vector IMs

Several proposals are explored for the hazard and intensity measure (IM) consistent selection of bidirectional ground motions to assess the performance of 3D structural models. Recent studies have shown the necessity of selecting records that thoroughly represent the seismicity at the site of interest, as well as the usefulness of efficient IMs capable of estimating the response of buildings with low scatter. However, the advances realized are mostly geared towards the structural analysis of 2D models. Few are the combined record, and IM selection approaches suggested expressly for nonlinear dynamic analysis of 3D structural models, especially when plan asymmetry and torsion sensitivity come into play. Conditional spectrum selection is leveraged and expanded here to offer a suite of approaches based on both scalar and vector IMs that convey information from two orthogonal horizontal components of the ground motion. Applications on multiple 3D building models highlight the importance of (a) employing the same IM for both record selection and response assessment and (b) maintaining hazard consistency in both horizontal components, when using either a scalar or a vector IM. All tested approaches that respect these conditions can be viable, yet the one based on the geometric mean of multiple spectral ordinates from both components over a period range seems to hold the most promise for general use.     

秋季太湖水下光场结构及其对水生态系统的影响

水生态系统中光能的分配很大程度上决定了水生态系统的结构和功能,利用2007年11-12月太湖水体光学特性和组分浓度数据,对秋季太湖水下光场结构特征和水体组分光竞争能力的表征光学量(漫衰减系数、平均余弦)和影响因素(吸收系数比重)进行了分析研究.结果表明,秋季太湖水下辐照度呈现单峰分布,最高值为583nm左右:根据Kd可将黄质和非色素物质主导程度的强弱分为弱、较强、强三个等级;Kd(PAR)平均值为4.61±1.54m-1,水体真光层厚度平均值为1.11±0.35m;太湖水下光场的光能主要分布在青光和黄绿光波长范围内,约占总能量的60%,蓝光和红光波长范围内的能量约占30%,这样的光谱结构有利于铜绿微囊藻和斜生栅藻的生长.     

1950～2000年黄河入海水沙的逐日变化及其影响因素

对 195 0～ 2 0 0 0年的 5 1年黄河经利津水文站入海的逐日水沙实测资料进行统计分析 ,展示了黄河入海水沙的极端情况 ,如洪峰和沙峰等逐渐变少到消失的过程及逐日入海水沙量的变化规律 ,并对其影响因素进行了初步分析。结果显示 :在黄河年均入海水沙锐减的背景下 ,2 0世纪 90年代日均流量 <5 0 0m3 /s的天数增加到了是 2 0世纪 5 0年代的 3倍多 ;>2 0 0 0m3 /s的洪季流量由初期频发至 1997年消失 ;>40 0 0m3 /s的洪峰流量自 1989年以后消失 ;>60 0 0m3 /s的洪峰流量从 1986年就不再发生。 >5 0 0t/s的日均输沙率 (4 3 .2× 10 6t/d的沙峰 )都出现在 1989年以前 ;>70 0t/s的日均输沙率 (60 .5× 10 6t/d的沙峰 )出现在 1977年以前 ;含沙量 >80kg/m3 的历时和可能产生异重流的含沙量天数都逐步减少。进入 2 0世纪 90年代以来基本不存在大洪灾及由沙峰造成的河口河道快速堆积而使河流改道现象 ,引起这些变化的因素除降雨外主要为流域耗水量不断增加及大中型水库的建成运用等 ,黄河入海水沙的这种锐减及水沙峰的极端情况都会对河口演变和河口区的生态环境产生重要的影响。     

两个30年气候平均值的差异及其对气候业务的影响

对陕西短期气候预测业务使用的39个气象站气候要素1971～2000年平均值与1961～1990年平均值进行比较,发现大部分地区年平均气温升高,晴天日数增加;年降水量、阴天日数、降水日数、年日照时数减少;年大风日数、沙尘暴日数、冰雹日数、雷暴日数减少等。把两个平均值放在长序列中分析,发现某些气象要素最近的气候平均值的差值在20世纪80～90年代期间的5～6个气候平均值每两个相邻平均值的差值中是最大的,说明陕西20世纪90年代气候较60年代干燥、温暖,且90年代气温高、降水少的特点十分突出。还对某些要素平均值改变对气候影响评价、气候预测业务产生的影响进行简要分析。     

基于GIS和洛伦兹曲线的辽宁省海岛空间分布特征探析

为合理开发利用和保护辽宁海岛资源,优化海岛产业结构和布局,以辽宁纳入中国海域海岛标准名录的633个海岛为研究对象,在对其地理坐标、陆域面积等基础数据进行数理统计的基础上,通过GIS的空间分析和绘制洛伦兹曲线、计算基尼系数,对辽宁海岛的空间分布特征进行了分析与研究。结果表明:辽宁海岛属集聚型分布类型,在空间上呈组团与带状分布,按聚集程度来看,形成了黄海北部长山群岛集聚区、辽东半岛南部浅海集聚带与辽东湾西部浅海集聚带三个区域;约60%的海岛处在离大陆岸线直线距离在10 km以内的范围内,近陆程度高;海岛陆域面积分布不均衡,基尼系数达0.9。     

海原断裂带中东段地貌差异及其成因探讨

以定量化地形因子为切入点的构造地貌学方法已成为活动构造研究的有效手段,被广泛用于定性或半定量解析地貌对新构造运动的响应及其演化过程。针对海原断裂带中东段现今地貌差异,以SRTM 90m分辨率DEM为基础,利用ArcGIS软件和Matlab程序脚本,提取了海原断裂带中东段高程、坡度、地形起伏、地形侵蚀以及河流陡峭系数等地形因子。从空间分布上看,上述各项地形因子沿断裂走向均呈现 “西高东低”的整体分布特征。西段海拔高、坡度陡、起伏大、侵蚀强、抬升快,中段和东段海拔低、坡度缓、起伏小、侵蚀弱、抬升慢,此外,在断裂带的东南尾端呈略微增加趋势,达到小范围内的峰值。在此基础上,通过对比分析地形因子与年降水量、基岩岩性,初步探讨了构造与降水、岩性等因素对地形地貌的控制作用,认为不同降水条件对地貌后期改造起显著作用,基岩岩性与现今地貌之间并无显著关系,该区域地貌类型主要受构造抬升差异所控制。沿断裂带走向上的现今地貌差异表明,西段处于相对快速的构造隆升和强挤压造山构造背景,中段由于受到黄河下切及河流冲积作用影响,地貌参数记录的抬升特征并不显著,而东段则反映出大型断裂带尾端挤压调整效应。     

未打穿散体材料桩复合地基固结简化解研究

利用平均固结度普遍解,将未打穿散体材料桩复合地基转化为等效双层地基,推导出未打穿散体材料桩复合地基固结简化解,给出其平均固结度计算简化公式。根据简化解、现有解和有限单元法编制程序,绘制考虑贯入比、桩体渗透系数、土体扰动区渗透系数和土体扰动区影响因子等因素对固结的影响曲线图,分析了各因素对未打穿散体材料桩复合地基固结性状的影响,并对3种计算方法所求出的平均固结度进行了比较研究,进而探讨了简化解的适用性。研究结果表明,无论是在单面排水还是在双面排水条件下,简化解求出的平均固结度总体上是介于现有解和有限单元法二者之间,尤其在固结过程的“中后期”（即固结度≥50%时）。从实用角度来看,简化解的计算结果具有足够精度,且简便、实用;在各种影响因素中,散体材料桩长度和桩体渗透系数对固结过程的影响最大。     

基于综合变异系数的地基承载力可靠性分析

采用一次二阶矩法,考虑岩土参量变异系数的试验及历史数据,提出了基于综合变异系数的地基承载力简化可靠性与风险分析方法。利用提出的简化可靠性与风险分析方法,可在传统确定性分析的基础上,采用合理的变异系数,分别得到与地基承载力有关的可靠度、破坏概率和平均期望损失的最可能值及其变化范围,为提出优化的地基基础设计方案和工程决策奠定基础;相对于土重度和黏聚力,地基承载力可靠度对于内摩擦角的变化更加敏感;设计时可以综合考虑所需的安全系数、可靠度和破坏概率,确定合适的基础宽度或基底面积。当地基土的场地勘察统计结果的变异系数权重逐渐增加时,综合变异系数不断减小,其相应的地基承载力的可靠度逐渐增加,而相应的破坏概率逐渐减少。