Constitutive and numerical framework for modeling joints and faults in rock

A three‐dimensional constitutive model for joints is described that incorporates nonlinear elasticity based on volumetric elastic strain, and plasticity for both compaction and shear with emphasis on compaction. The formulation is general in the sense that alternative specific functional forms and evolution equations can be easily incorporated. A corresponding numerical structure based on finite elements is provided so that a joint width can vary from a fraction of an element size to a width that occupies several elements. The latter case is particularly appropriate for modeling a fault, which is considered simply to be a joint with large width. For small joint widths, the requisite equilibrium and kinematic requirements within an element are satisfied numerically. The result is that if the constitutive equation for either the joint or the rock is changed, the numerical framework remains unchanged. A unique aspect of the general formulation is the capability to handle either pre‐existing gaps or the formation of gaps. Representative stress–strain plots are given to illustrate both the features of the model and the effects of changes in values of material parameters. Copyright © 2015 John Wiley & Sons, Ltd.     

Diagnosing snow accumulation errors in a rain‐snow transitional environment with snow board observations

Diagnosing the source of errors in snow models requires intensive observations, a flexible model framework to test competing hypotheses, and a methodology to systematically test the dominant snow processes. We present a novel process‐based approach to diagnose model errors through an example that focuses on snow accumulation processes (precipitation partitioning, new snow density, and snow compaction). Twelve years of meteorological and snow board measurements were used to identify the main source of model error on each snow accumulation day. Results show that modeled values of new snow density were outside observational uncertainties in 52% of days available for evaluation, while precipitation partitioning and compaction were in error 45% and 16% of the time, respectively. Precipitation partitioning errors mattered more for total winter accumulation during the anomalously warm winter of 2014–2015, when a higher fraction of precipitation fell within the temperature range where partition methods had the largest error. These results demonstrate how isolating individual model processes can identify the primary source(s) of model error, which helps prioritize future research.     

From soil aggregates to riverine flocs: a laboratory experiment assessing the respective effects of soil type and flow shear stress on particles characteristics

Particles eroded from hillslopes and exported to rivers are recognized to be composite particles of high internal complexity. Their architecture and composition are known to influence their transport behaviour within the water column relative to discrete particles. To‐date, hillslope erosion studies consider aggregates to be stable once they are detached from the soil matrix. However, lowland rivers and estuaries studies often suggest that particle structure and dynamics are controlled by flocculation within the water column. In order to improve the understanding of particle dynamics along the continuum from hillslopes to the lowland river environment, soil particle behaviour was tested under controlled laboratory conditions. Seven flume erosion and deposition experiments, designed to simulate a natural erosive event, and five shear cell experiments were performed using three contrasting materials: two of them were poorly developed and as such can not be considered as soils, whilst the third one was a calcareous brown soil. These experiments revealed that soil aggregates were prone to disaggregation within the water column and that flocculation may affect their size distribution during transport. Large differences in effective particle size were found between soil types during the rising limb of the bed shear stress sequence. Indeed, at the maximum applied bed shear stress, the aggregated particles median diameter was found to be three times larger for the well‐developed soil than for the two others. Differences were smaller in the falling limb, suggesting that soil aggregates underwent structural changes. However, characterization of particles strength parameters showed that these changes did not fully turn soil aggregates into flocs, but rather into hybrid soil aggregate–floc particles. Copyright © 2013 John Wiley & Sons, Ltd.     

Variational analysis of dynamics of fissured poroelastic rocks

Rocks can be modeled in a continuum framework as fissured, poroelastic materials, i.e., materials with two degrees of porosity, one due to the fissures and another one due to the pores. The governing equations of motion of fissured poroelastic rocks established by Beskos are rederived here by establishing a variational statement which also provides the boundary conditions of the problem. This is accomplished by considering strain, dissipation and kinetic energies as well as the work of external forces. The above statement is also derived here by employing the method of weighted residuals.     

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.     

Retrieval of effective leaf area index (LAIe) and leaf area density (LAD) profile at individual tree level using high density multi-return airborne LiDAR

As an important canopy structure indicator, leaf area index (LAI) proved to be of considerable implications for forest ecosystem and ecological studies, and efficient techniques for accurate LAI acquisitions have long been highlighted. Airborne light detection and ranging (LiDAR), often termed as airborne laser scanning (ALS), once was extensively investigated for this task but showed limited performance due to its low sampling density. Now, ALS systems exhibit more competing capacities such as high density and multi-return sampling, and hence, people began to ask the questions like—“can ALS now work better on the task of LAI prediction?” As a re-examination, this study investigated the feasibility of LAI retrievals at the individual tree level based on high density and multi-return ALS, by directly considering the vertical distributions of laser points lying within each tree crown instead of by proposing feature variables such as quantiles involving laser point distribution modes at the plot level. The examination was operated in the case of four tree species (i.e. Picea abies, Pinus sylvestris, Populus tremula and Quercus robur) in a mixed forest, with their LAI-related reference data collected by using static terrestrial laser scanning (TLS). In light of the differences between ALS- and TLS-based LAI characterizations, the methods of voxelization of 3D scattered laser points, effective LAI (LAIe) that does not distinguish branches from canopies and unified cumulative LAI (ucLAI) that is often used to characterize the vertical profiles of crown leaf area densities (LADs) was used; then, the relationships between the ALS- and TLS-derived LAIes were determined, and so did ucLAIs. Tests indicated that the tree-level LAIes for the four tree species can be estimated based on the used airborne LiDAR (R² = 0.07, 0.26, 0.43 and 0.21, respectively) and their ucLAIs can also be derived. Overall, this study has validated the usage of the contemporary high density multi-return airborne LiDARs for LAIe and LAD profile retrievals at the individual tree level, and the contribution are of high potential for advancing forest ecosystem modeling and ecological understanding.     

基于规则格网的DEM插值实验

通过地形建模,将6个地形隶属函数按照101×101,112×112,126×126,144×144,168×168,201×201,257×257 7种格网密度生成6种局部地形单元的规则格网DEM;使用反距离加权(IDW)等7种插值算法,将前6种格网密度下的DEM插值成257×257规格;从原始257×257DEM中随机抽取检查点计算残差,并对残差中误差进行分析。通过分组插值实验,运用控制变量法、方差分析等方法研究地貌类型、采样密度和插值算法对DEM插值精度的影响。     

水库群水沙调控的单-多目标调度模型及其应用

为缓解内蒙古河段"二级悬河"形势,以黄河上游沙漠宽谷河段为研究对象,以龙羊峡水库、刘家峡水库为调控主体,开展黄河上游水沙调控研究。建立了输沙量、发电量最大的单目标模型以及多目标模型;分别采用自迭代模拟算法、逐次逼近动态规划算法(DPSA)和改进的非支配排序遗传优化算法(NSGA-Ⅱ)求解模型;设置了初始、常规、优化和联合优化4种方案。通过实例计算,联合优化调度方案的区间总冲刷量达到了0.38亿t,梯级发电量148.22亿kW·h。该方案以较小的电量损失换来了输沙量的大幅度增加,水沙调控效果显著,推荐为最优方案。研究成果量化了水沙调控效果和各目标间的转化规律,为开展黄河上游水沙调控提供了决策依据,具有重要的应用价值和实际指导意义。     

基于重力全张量数据的Parker-Oldenburg算法研究

相比重力数据,重力张量数据通常包含更多的异常信息。本文根据重力数据与重力张量数据的关系,利用位场转化技术,将重力张量数据应用于传统的Parker-Oldenburg密度界面反演算法中。通过模型试验,证明了在网格间距较大或者数据存在一定噪音时,使用本文算法进行反演能得到更好的效果。实验结果说明利用重力张量数据可以有效地提高密度界面反演的分辨率。     

基于峰值因子最大化的鬼波压制方法

针对传统拖缆采集得到的海上地震资料受鬼波影响而导致频带宽度变窄,笔者提出一种去鬼波方法,可以压制鬼波效应,拓宽频带。首先在频率-空间域推导出镜像数据生成公式,并基于峰值因子最大化的参数搜索方法来估计最优的鬼波延迟时间,将最优的鬼波延迟时代入镜像记录生成公式,获得最优的镜像记录,利用联合反褶积算法得到最优的去鬼波结果。将本文中的方法应用于模拟的变深度采集数据和海上实测的变深度采集数据,证实该方法可以有效地压制鬼波,消除频谱中的陷波点,得到宽频数据。