利用热红外温差识别沙尘

沙尘识别是沙尘灾害监测和沙尘气溶胶特性研究的首要工作。利用辐射传输方程进行了沙尘气溶胶的辐射计算 ,对不同沙尘气溶胶光学厚度下的热红外通道温差ΔT(T11μm-T12 μm)的变化进行了分析。理论分析表明 ,利用热红外通道的温度差ΔT进行陆地沙尘识别是可行的。并利用NOAA AVHRR热红外通道的温度差ΔT进行了陆地沙尘识别的实验 ,经与地面气象站实测的沙尘天气现象相比较 ,结果一致。     

Splash-saltation trajectories of soil particles under wind-driven rain

It is usually recognized that relatively large amounts of soil particles cannot be transported by raindrop splashes under windless rain. However, the splash-saltation process can cause net transportation in the prevailing wind direction since variations in splash-saltation trajectory due to the wind are expected in wind-driven rain. Therefore, determining the combined effect of rain and wind on the process should enable improvement of the estimation of erosion for any given prediction technique. This paper presents experimental data on the effects of slope aspect, slope gradient, and horizontal wind velocity on the splash-saltation trajectories of soil particles under wind-driven rain. In a wind tunnel facility equipped with a rainfall simulator, the rains driven by horizontal wind velocities of 6, 10, and 14 m s⁻¹ were allowed to impact three agricultural soils packed into 20×55 cm soil pans placed at both windward and leeward slopes of 7%, 15%, and 20%. Splash-saltation trajectories were measured by trapping the splashed particles at distances downwind on a 7-m uniform slope segment in the upslope and downslope directions, respectively, for windward and leeward slopes. Exponential decay curves were fitted for the mass distribution of splash-saltation sediment as a function of travel distance, and the average splash-saltation trajectory was derived from the average value of the fitted functions. The results demonstrated that the average trajectory of a raindrop-induced and wind-driven soil particle was substantially affected by the wind shear velocity, and it had the greatest correlation (r=0.96 for all data) with the shear velocity; however, neither slope aspect nor slope gradient significantly predicted the splash-saltation trajectory. More significantly, a statistical analysis conducted with nonlinear regression model of C₁(u_*²/g) showed that average trajectory of splash saltation was approximately three times greater than that of typical saltating sand grain.     

A parameterisation for the threshold shear velocity to initiate deflation of dry and wet sediment

A new parameterisation for the threshold shear velocity to initiate deflation of dry and wet particles is presented. It is based on the balance of moments acting on particles at the instant of particle motion. The model hence includes a term for the aerodynamic forces, including the drag force, the lift force and the aerodynamic-moment force, and a term for the interparticle forces. The effect of gravitation is incorporated in both terms. Rather than using an implicit function for the effect of the aerodynamic forces as reported earlier in literature, a constant aerodynamic coefficient was introduced. From consideration of the van der Waals force between two particles, it was further shown that the effect of the interparticle cohesion force between two dry particles on the deflation threshold should be inversely proportional to the particle diameter squared. The interparticle force was further extended to include wet bonding forces. The latter were considered as the sum of capillary forces and adhesive forces. A model that expresses the capillary force as a function of particle diameter squared and the inverse of capillary potential was deduced from consideration of the well-known model of Fisher and the Young–Laplace equation. The adhesive force was assumed to be equal to tensile strength, and a function which is proportional to particle diameter squared and the inverse of the potential due to adhesive forces was derived. By combining the capillary-force model and the adhesive force model, the interparticle force due to wet bonding was simplified and written as a function of particle diameter squared and the inverse of matric potential. The latter was loglinearly related to the gravimetric moisture content, a relationship that is valid in the low-moisture content range that is important in the light of deflation of sediment by wind. By introducing a correction to force the relationship to converge to zero moisture content at oven dryness, the matric potential–moisture content relationship contained only one unknown model parameter, viz. moisture content at −1.5 MPa. Working out the model led to a rather simple parameterisation containing only three coefficients. Two parameters were incorporated in the term that applies to dry sediment and were determined by using experimental data as reported by Iversen and White [Sedimentology 29 (1982) 111]. The third parameter for the wet-sediment part of the model was determined from wind-tunnel experiments on prewetted sand and sandy loam aggregates. The model was validated using data from wind-tunnel experiments on the same but dry sediment, and on data obtained from simulations with the model of Chepil [Soil Sci. Soc. Am. Proc. 20 (1956) 288]. The experiments showed that soil aggregates should be treated as individual particles with a density equal to their bulk density. Furthermore, it was shown that the surface had to dry to a moisture content of about 75% of the moisture content at −1.5 MPa before deflation became sustained. The threshold shear velocities simulated with our model were found to be in good agreement with own observations and with simulations using Chepil's model.     

Fractally Invariant Distributions and an Application in Geochemical Exploration

Fractal modelling has been applied extensively as a means of characterizing the spatial distribution of geological phenomena that display self-similarity at differing scales of measurement. A fractal distribution exists where the number of objects exhibiting values larger than a specified magnitude displays a power-law dependence on that magnitude, and where this relationship is scale-invariant. This paper shows that a number of distributions, including power-function, Pareto, log-normal and Zipf, display fractal properties under certain conditions and that this may be used as the mathematical basis for developing fractal models for data exhibiting such distributions. Population limits, derived from fractal modelling using a summation method, are compared with those derived from more conventional probability plot modelling of stream sediment geochemical data from north-eastern New South Wales. Despite some degree of subjectivity in determining the number of populations to use in the models, both the fractal and probability plot modelling have assisted in isolating anomalous observations in the geochemical data related to the occurrence of mineralisation or lithological differences between sub-catchments. Thresholds for the main background populations determined by the fractal model are similar to those established using probability plot modelling, however the summation method displays less capacity to separate out anomalous populations, especially where such populations display extensive overlap. This suggests, in the geochemical data example provided, that subtle differences in the population parameters may not significantly alter the fractal dimension.     

Detection and 2d modelling of cavities using pole–dipole array

The resistivity method is often used in cave prospecting. In this paper the pole–dipole array ability to detect cavities at different depths and with different water contents is investigated. The research was performed using analogical and numerical modelling. According to the results, empty caves at a depth less than four times its diameter can be easily detected. The ability of the pole–dipole array to detect water-filled caves reaches a minimum at filling percentages of 30–50%. Overburden effects research shows that low resistivity overburden reduces the resolution capability of the array. This study shows that equivalent results can be obtained by modelling the empty caves as infinite resistivity bodies or, alternatively, as very high resistivity ones. The analysis of field data acquired in the Maciço Calcário Estremenho (Portugal) shows the practical importance of the pole–dipole array in cavities prospecting.     

复杂条件下的地下管线探测技术

以深圳市某大道顶管施工地下管线的探测成果为例,研究了在复杂场地条件下,探测地下管线的地球物理方法和技术,着重研究了电磁感应探测法.认为在复杂条件下,对钢质煤气管,应采用电磁感应探测法;对钢质或铁质给水管,应采用直接法;对电力管和电信管,应尽量采用夹钳法;对排水(管)渠等非金属材质管线,应采用地质雷达探测的方法.     

Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of the response of a reinforced wall to a high speed train passage

The dynamic response of a mechanically stabilized earth wall to the passing of a high‐speed train is modelled using the finite element method. A three‐dimensional analysis is carried out, using a specific framework that allows performing the analysis with a moderate computational effort. In the first place, a so‐called multiphase approach is used to take into account the reinforcing strips. The moving load is taken into account by performing the calculation in a mobile referential using the properties of symmetry of the train cars and a simplifying assumption of periodicity for the whole train. We also assume a steady state. A partial validation of the approach is obtained by means of a comparison with an analytical solution. The quick increase in displacements induced by the train passing when the speed comes close to the celerity of Rayleigh waves clearly appears in the results. The vertical displacements, vertical stresses in the backfill, tensile forces in the strips and the influence of the stiffness of the soil are discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Characterisation of hydroclimatological trends and variability in the Lake Naivasha basin,Kenya

Recent hydro‐climatological trends and variability characteristics were investigated for the Lake Naivasha basin with the aim of understanding the changes in water balance components and their evolution over the past 50 years. Using a Bayesian change point analysis and modified Mann–Kendall tests, time series of annual mean, maximum, minimum, and seasonal precipitation and flow, as well as annual mean lake volumes, were analysed for the period 1960–2010 to uncover possible abrupt shifts and gradual trends. Double cumulative curve analysis was used to investigate the changes in hydrological response attributable to either human influence or climatic variability. The results indicate a significant decline in lake volumes at a mean rate of 9.35 × 10⁶ m³ year^?1. Most of the river gauging stations showed no evidence of trends in the annual mean and maximum flows as well as seasonal flows. Annual minimum flows, however, showed abrupt shifts and significant (upward/downward) trends at the main outlet stations. Precipitation in the basin showed no evidence of abrupt shifts, but a few stations showed gradual decline. The observed changes in precipitation could not explain the decline in both minimum flows and lake volumes. The findings show no evidence of any impact of climate change for the Lake Naivasha basin over the past 50 years. This implies that other factors, such as changes in land cover and infrastructure development, have been responsible for the observed changes in streamflow and lake volumes. Copyright © 2015 John Wiley & Sons, Ltd.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.