潍坊市浅层地热能资源调查及开发利用研究

潍坊市浅层地热能资源丰富,调查区内浅层地热热容量可折合标准煤1355．28万t,利用范围较广。该文对资源开发利用中存在的勘查论证不充分,开发过程和开发方式不合理等问题提出了合理建议,为综合开发利用该资源,改善环境,建设可持续发展社会提供理论参考。     

降雨诱发浅层滑坡渐进破坏分析研究

为分析降雨诱发浅层滑坡的演变过程。本文以湖南湘西古丈滑坡为例,基于Green-Ampt入渗模型,进行了降雨诱发浅层滑坡渐进破坏分析。研究结果表明:在强降雨作用下,滑坡的失稳破坏主要是由于前缘土体以及中前部土体的局部破坏,而逐渐发展为整体破坏。并且,受滑坡地形影响,地形平缓的区域虽然湿润锋下渗较快,土体抗剪强度较低,但由于土体饱和带的渗流作用较小,而重力提供垂直于滑面的分力较大,该部分稳定性较为良好,故湿润锋对于滑坡稳定性的影响还应该根据不同地形条件加以分析。渐进式滑坡破坏分析方法对滑坡的监测和防治具有重要的指导意义。     

Identification of sources and groundwater recharge zones from hydrochemistry and stable isotopes of an agriculture-based paleo-lacustrine basin of drought-prone northeast Mexico

Over exploitation for agricultural activities and consumption has depleted the groundwater resources of drought-prone northeast Mexico. Major ion concentrations along with δ¹⁸O_H2O, δ²H_H2O and d-excess values of shallow groundwater from the Cieneguilla Basin (near Tula) located at a distance of ∼200 km from coast of the Gulf of Mexico helped to contribute new data about drought vulnerability in this region through identification of the moisture source and groundwater recharge zone. Different degrees of rock-water interaction through gypsum, anhydrite and halite dissolutions and minor silicate weathering controlled the hydrochemistry. Stable isotopes yielded a least square regression and slope similar to the local as well as global meteoric water lines, indicating minimal effect of evaporation during the recharge as well as in the subsoil. Isotopic fractionations along with a digital elevation model demarcated the recharge zones at north and east of the basin, with altitudinal difference of ≥1000 m, and indicated that the recharge occurred through warm season moisture sourced from the Gulf of Mexico. Less frequent landfalling of tropical storms caused by warmer sea surface temperature, however, has reduced this rainfall over the last few decades. If the trend of global warming continues unabated, the depleted groundwater resources would trigger reduction in agricultural activities in this drought-prone region and lead to enhanced socio-economic challenges.     

Nearshore sticky waters

Wind- and current-driven flotsam, oil spills, pollutants, and nutrients, approaching the nearshore will frequently appear to slow down/park just beyond the break zone, where waves break. Moreover, the portion of these tracers that beach will do so only after a long time. Explaining why these tracers park and at what rate they reach the shore has important implications on a variety of different nearshore environmental issues, including the determination of what subscale processes are essential in computer models for the simulation of pollutant transport in the nearshore. Using a simple model we provide an explanation for the underlying mechanism responsible for the parking of tracers, not subject to inertial effects, the role played by the bottom topography, and the non-uniform dispersion which leads, in some circumstances, to the eventual landing of all or a portion of the tracers. We refer to the parking phenomenon in this environment as nearshore sticky waters.     

Optical remote sensing of submerged aquatic vegetation: Opportunities for shallow clearwater streams

Remote sensing has rarely been used as a tool to map and monitor submerged aquatic vegetation (SAV) in rivers, due to a combination of insufficient spatial resolution of available image data and strong attenuation of light in water through absorption and scattering. The latter process reduces the possibility to use spectral reflectance information to accurately classify submerged species. However, increasing availability of very high resolution (VHR) image data may enable the use of shape and texture features to help discriminate between species by taking an object based image analysis (OBIA) approach, and overcome some of the present limitations.This study aimed to investigate the possibility of using optical remote sensing for the detection and mapping of SAV. It firstly looked at the possibilities to discriminate submerged macrophyte species based on spectral information only. Reflectance spectra of three macrophyte species were measured in situ across a range of submergence depths. The results showed that water depth will be a limiting factor for the classification of species from remote sensing images. Only Spiked Water Milfoil (Myriophyllum spicatum) was indicated as spectrally distinct through ANOVA analysis, but subsequent Jeffries–Matusita distance analysis did not confirm this. In particular Water Crowfoot (Ranunculus fluitans) and Pondweed (Potamogeton pectinatus) could not be discriminated at 95% significance level. Spectral separability of these two species was also not possible without the effect of an overlying water column.Secondly, the possibility to improve species discrimination, using spatial and textural information was investigated for the same SAV species. VHR image data was acquired with a Near Infrared (NIR) sensitive DSLR camera from four different heights including a telescopic pole and a Helikite UAS. The results show that shape and texture information can improve the detection of the spectrally similar Pondweed and Water Crowfoot from VHR image data. The best performing feature ‘length/width ratio of sub-objects’ was obtained through expert knowledge. All of the shape and texture based features performed better at species differentiation than the spectrally based features.In conclusion this study has shown that there is considerable potential for the combination of VHR data and OBIA to map SAV in shallow stream environments, which can benefit species monitoring and management.     

Balancing the source terms in a SPH model for solving the shallow water equations

A shallow flow generally features complex hydrodynamics induced by complicated domain topography and geometry. A numerical scheme with well-balanced flux and source term gradients is therefore essential before a shallow flow model can be applied to simulate real-world problems. The issue of source term balancing has been exhaustively investigated in grid-based numerical approaches, e.g. discontinuous Galerkin finite element methods and finite volume Godunov-type methods. In recent years, a relatively new computational method, smooth particle hydrodynamics (SPH), has started to gain popularity in solving the shallow water equations (SWEs). However, the well-balanced problem has not been fully investigated and resolved in the context of SPH. This work aims to discuss the well-balanced problem caused by a standard SPH discretization to the SWEs with slope source terms and derive a corrected SPH algorithm that is able to preserve the solution of lake at rest. In order to enhance the shock capturing capability of the resulting SPH model, the Monotone Upwind-centered Scheme for Conservation Laws (MUSCL) is also explored and applied to enable Riemann solver based artificial viscosity. The new SPH model is validated against several idealized benchmark tests and a real-world dam-break case and promising results are obtained.     

Simulation of dam-break waves on movable beds using a multi-stage centered scheme

This paper presents the application of the multi-stage first-order centered scheme GMUSTA to solve a two-phase flow model with four equations for simulating dam-break floods without and with sediment transport.Computation of generalized Riemann invariants can be particularly complex and costly in simulating dam-break floods with sediment transport.GMUSTA numerical scheme,which does not require complete knowledge of the eigenstructure of the hyperbolic mathematical model,offers a suitable and attractive option.The quality of the dam-break flood simulations with GMUSTA scheme is verified by comparing the results against laboratory tests and some experimental data available in the literature,on fixed and mobile bed conditions,with different bed materials and flush or stepped bottoms.The numerical results reproduce quite well the experimental evidence,proving that the model is capable of predicting the temporal evolution of the free-surface and the bed.The GMUSTA scheme,which is not only simple to implement but also both accurate and computationally efficient,is proposed as an appropriate tool for integrating non-equilibrium sediment-transport models.     

Evaluation of undrained failure envelopes of caisson foundations under combined loading

In this paper, results of a three-dimensional finite element study addressing the effect of embedment ratio (L/D) of caisson foundations on the undrained bearing capacity under uniaxial and combined loadings are discussed. The undrained response of caisson foundations under uniaxial vertical (V), horizontal (H) and moment (M) loading are investigated. A series of equations are proposed to predict the ultimate vertical, moment and maximum horizontal bearing capacity factors. The undrained response of caisson foundations under combined V-H and V-M load space is studied and presented using failure envelopes generated with side-swipe method. The kinematic mechanism accompanying failure under uniaxial loading is addressed and presented for different embedment ratios. Predictions of the uniaxial bearing capacities are compared with other models and it is confirmed that the proposed equations appropriately describe the capacity of caisson foundations under uniaxial vertical, horizontal and moment loading in homogenous undrained soils. The results of this paper can be used as a basis for standard design codes of off-shore skirted shallow foundations which will be the first of its kind.     

Refined seaward boundary conditions for modelling shallow water waves in semi-enclosed water bodies

Based on the theory of characteristics, this research elaborates on the numerical treatment of two types of seaward boundary conditions for modelling long-wave dynamics in truncated estuarine and coastal domains. These seaward boundary conditions are devised for the solution of the fully non-linear shallow water equations in the time domain. The first type is the clamped boundary, at which the water level variation is given and the velocity is computed along the characteristic line going out of the domain. The second type is the non-reflecting boundary, where the incident wave information is introduced and the reflected waves from inside the computational domain are allowed to escape at the same time. The essence of its numerical implementation is to distinguish the inward and outward characteristics and to disconnect the incoming characteristic relation from the actual flow inside the domain. Compared with previous techniques, the present method includes extra terms in the derivation to account for the effects of the uneven bed, bottom friction and shape of the characteristic lines. A shock-capturing finite difference method is used to solve the shallow water equations in the deviatoric format, but the seaward boundary algorithms constructed herein are generic and applicable to other solvers. The necessity of these refinements is highlighted by simulating the tidal oscillation in the Persian/Arabian Gulf, periodic wave runup on the coastline and the wave resonance in a narrow harbour. It is found that neglecting the bed slope at the boundary may result in biased mean water levels in the prediction.     

Interpretation of gravity profiles across the northern Oaxaca terrane,its boundaries and the Tehuacán Valley,southern Mexico

A gravity study was conducted across the northern Oaxaca terrane and its bounding faults: the Caltepec and Oaxaca Faults to the west and east, respectively. These faults juxtapose the Oaxaca terrane against the Mixteca and Juarez terranes, respectively. The Oaxaca Fault also forms the eastern boundary of the Cenozoic Tehuacán depression. On the west, at depth, the Tehuacán valley is limited by the normal buried Tehuacán Fault. This gravity study reveals that the Oaxaca Fault system gives rise to a series of east tilted basamental blocks (Oaxaca Complex). The tectonic depression is filled with Phanerozoic rocks and has a deeper depocenter to the west. The gravity data also indicate that on the west, the Oaxaca Complex, the Caltepec and Santa Lucia faults continue northwestwards beneath Phanerozoic rocks. A major E–W to NE–SW discontinuity is inferred to exist between profiles 1 and 2.