任意各向异性介质三维非结构谱元法直流电阻率正演模拟研究

各向异性是地电异常解释中不可忽视的因素,广泛存在于裂隙或层理发育的地质环境中.本文针对任意各向异性条件下直流电阻率法三维正演问题进行研究,结合非结构谱元法建立模拟算法,充分利用谱方法的指数收敛性以及非结构有限元对地形和复杂异常体刻画能力,提高计算精度和效率.通过灵活的四面体网格剖分和高阶谱插值,实现了复杂介质任意各向异性模型电阻率响应的高精度数值模拟.我们首先通过层状各向异性模型验证本文非结构谱元法的计算精度,进而我们以半空间中立方体模型为例分析各向异性对电阻率响应的影响特征,并通过计算针对不同各向异性参数的视电阻率极性图,探究地下介质各向异性特征识别方法.最后,我们针对典型的山脊模型计算和分析存在地形效应条件下各向异性电流场分布及视电阻率特征.模型计算结果表明基于四面体网格的谱元法模拟带有复杂地形和异常体的任意各向异性模型具有很高的计算精度.本文的研究成果将在推进电阻率方法用于解决裂隙及层理等环境和工程地质问题中发挥积极作用.

     

Climate-change driven increased flood magnitudes and frequency in the British uplands: Geomorphologically informed scientific underpinning for upland flood-risk management

Upland river systems in the UK are predicted to be prone to the effects of increased flood magnitudes and frequency, driven by climate change. It is clear from recent events that some headwater catchments can be very sensitive to large floods, activating the full sediment system, with implications for flood risk management further down the catchment. We provide a 15-year record of detailed morphological change on a 500-m reach of upland gravel-bed river, focusing upon the geomorphic response to an extreme event in 2007, and the recovery in the decade following. Through novel application of two-dimensional (2D) hydrodynamic modelling we evaluate the different energy states of pre- and post-flood morphologies of the river reach, exploring how energy state adjusts with recovery following the event. Following the 2007 flood, morphological adjustments resulted in changes to the shear stress population over the reach, resulting in higher shear stresses. Although the proportion of shear stresses in excess of those experienced using the pre-flood digital elevation model (DEM) varied over the recovery period, they remained substantially in excess of those experienced pre-2007, suggesting that there is still potential for enhanced bedload transport and morphological adjustment within the reach. Although volumetric change calculated from DEM differencing does indicate a reduction in erosion and deposition volumes in the decade following the flood, we argue that the system still has not fully recovered to the pre-flood state. We further argue that Thinhope Burn, and other similarly impacted catchments in upland environments, may not recover under the wet climatic phase currently being experienced. Hence systems like Thinhope Burn will continue to deliver large volumes of sediment further down river catchments, providing new challenges for flood risk management into the future.     

Evaluating soil water routing approaches in watershed-scale,ecohydrologic modelling

Soil water dynamics are central in linking and regulating natural cycles in ecohydrology, however, mathematical representation of soil water processes in models is challenging given the complexity of these interactions. To assess the impacts of soil water simulation approaches on various model outputs, the Soil and Water Assessment Tool was modified to accommodate an alternative soil water percolation method and tested at two geographically and climatically distinct, instrumented watersheds in the United States. Soil water was evaluated at the site scale via measured observations, and hydrologic and biophysical outputs were analysed at the watershed scale. Results demonstrated an improved Kling–Gupta Efficiency of up to 0.3 and a reduction in percent bias from 5 to 25% at the site scale, when soil water percolation was changed from a threshold, bucket-based approach to an alternative approach based on variable hydraulic conductivity. The primary difference between the approaches was attributed to the ability to simulate soil water content above field capacity for successive days; however, regardless of the approach, a lack of site-specific characterization of soil properties by the soils database at the site scale was found to severely limit the analysis. Differences in approach led to a regime shift in percolation from a few, high magnitude events to frequent, low magnitude events. At the watershed scale, the variable hydraulic conductivity-based approach reduced average annual percolation by 20–50 mm, directly impacting the water balance and subsequently biophysical predictions. For instance, annual denitrification increased by 14–24 kg/ha for the new approach. Overall, the study demonstrates the need for continued efforts to enhance soil water model representation for improving biophysical process simulations.     

An overview of hydrometeorological datasets from a small agricultural catchment (Nučice) in the Czech Republic

We introduce the freely available web-based Water in an Agricultural Landscape—NUčice Database (WALNUD) dataset that includes both hydrological and meteorological records at the Nučice experimental catchment (0.53 km²), which is representative of an intensively farmed landscape in the Czech Republic. The Nučice experimental catchment was established in 2011 for the observation of rainfall–runoff processes, soil erosion processes, and water balance of a cultivated landscape. The average altitude is 401 m a.s.l., the mean land slope is 3.9%, and the climate is humid continental (mean annual temperature 7.9°C, annual precipitation 630 mm). The catchment is drained by an artificially straightened stream and consists of three fields covering over 95% of the area which are managed by two different farmers. The typical crops are winter wheat, rapeseed, and alfalfa. The installed equipment includes a standard meteorological station, several rain gauges distributed across the basin, and a flume with an H-type facing that is used to monitor stream discharge, water turbidity, and basic water quality indicators. Additionally, the groundwater level and soil water content at various depths near the stream are recorded. Recently, large-scale soil moisture monitoring efforts have been introduced with the installation of two cosmic-ray neutron sensors for soil moisture monitoring. The datasets consist of observed variables (e.g. measured precipitation, air temperature, stream discharge, and soil moisture) and are available online for public use. The cross-seasonal, open access datasets at this small-scale agricultural catchment will benefit not only hydrologists but also local farmers.     

Dynamic soil–structure interaction analysis of a telescope at the Javalambre Astrophysical Observatory

This paper presents the dynamic soil–structure analysis of the main telescope T250 of the Observatorio Astrofísico de Javalambre (OAJ, Javalambre Astrophysical Observatory) on the Pico del Buitre. Vibration control has been of prime concern in the design, since astrophysical observations may be hindered by mechanical vibration of optical equipment due to wind loading. The telescope manufacturer therefore has imposed a minimal natural frequency of 10 Hz for the supporting telescope pier. Dynamic soil–structure interaction may significantly influence the lowest natural frequency of a massive construction as a telescope pier. The structure clamped at its base has a resonance frequency of 14.3 Hz. A coupled finite element–boundary element (FE–BE) model of the telescope pier that accounts for the dynamic interaction of the piled foundation and the soil predicts a resonance frequency of 11.2 Hz, demonstrating the significant effect of dynamic soil–structure interaction. It is further investigated to what extent the coupled FE–BE model can be simplified in order to reduce computation time. The assumption of a rigid pile cap allows us to account for dynamic soil–structure interaction in a simplified way. A coupled FE–BE analysis with a rigid pile cap predicts a resonance frequency of 11.7 Hz, demonstrating a minor effect of the pile cap flexibility on the resonance frequency of the telescope pier. The use of an analytical model for the pile group results in an overestimation of the dynamic soil stiffness. This error is due to the large difference between the actual geometry and the square pile cap model for which the parameters have been tuned.     

Structure of residual oil as a function of wettability using pore-network modelling

In the water flooding of mixed-wet porous media, oil may drain down to relatively low residual oil saturations (S_or). Various studies have indicated that such low saturations can only be reached when oil layers in pore corners are included in the pore-scale modelling. These processes within a macroscopic porous medium can be modelled at the pore-scale by incorporating the fundamental physics of capillary dominated displacement within idealised pore network models. Recently, the authors have developed thermodynamic criteria for oil layer existence in pores with non-uniform wettability which takes as input geometrically and topologically representative networks, to calculate realistic S_or values for mixed-wet and oil-wet sandstones [16, 21]. This previous work is developed in this paper to include (i) the visualisation of the 3D structure of this residual oil, and (ii) a statistical analysis of this “residual/remaining” oil. Both the visualisation and the statistical analysis are done under a wide range of wettability conditions, which is reported for the first time in this paper.The structure of residual oil for strongly water wet systems is well known (where residual = remaining oil) and our model agrees with this but this structure changes radically for mixed wet systems (where residual ≠ remaining) and this has not yet been visualised experimentally. We find that for more water-wet systems high final residual oil saturations are reached at relatively small amounts of water injected and this oil is present in the pores as bulk oil. On the other hand, for more oil-wet systems we find a slow decrease of the amount of remaining oil with increasing amounts of injected water. During the process, the remaining connectivity of the oil phase is increasingly provided by oil layers only, hence the slow drainage. The final residual oil saturation, only reached in the theoretical limit of an infinite amount of injected water, is almost entirely contained in large number of (relatively low volume) oil layers, which are present in pores of most radius sizes.     

基于二进小波的数字形变资料趋势异常分析探讨

利用二进正交离散小波，通过Mallat快速算法对形变倾斜资料进行了处理，并将计算结果与该地区的地震活动进行了对比分析。综合分析认为，小波分解的本质是窗口加权滑动平均的结果，能有效地对各种不同频率成分的形变序列按不同的频段分离和重构，可以起到明显放大低频段趋势成份的作用，为有效识别形变趋势异常提供一种可行方法。     

全新世初期气候的数值模拟研究

本文尝试用一个完全耦合的大气-海洋-植被模式(Atmosphere-Ocean-Vegetation General Circulation Model, AOVGCM)模拟全新世初期的地球气候.模拟结果表明,耦合模式成功地再现了全新世初期复杂气候条件下的基本特征,为研究这一时期的气候状况提供了重要的模型支持.就全球平均而言,11ka BP冬季地表气温比现代约低1.6 K,夏季比现代低约0.3 K;大气温度从低层到高层有不同的表现,高层反映了太阳辐射的重要作用,而低层气候对下垫面（如冰川、植被和海洋等）的影响比较敏感.从区域分布来看,11 ka BP冬季大部地区比现在偏干,但热带太平洋和南半球少部地区降水偏多;夏季大部地区比现在偏湿,亚洲和非洲季风偏强,主要的季风区降水偏多.     

地下RC结构物地震响应特征土工离心试验的模拟

地下结构物抗震稳定性的研究具有重要实际意义。基于动态土工离心试验,应用有效应力法分析了不同水平地震作用下土层和地下RC结构物的最大水平位移、不同埋深的土层加速度和超孔隙水压力变化规律,以及结构物的破坏发展特征,并进一步与试验结果进行了相互验证和对比,取得了较好效果,为重要地下结构物的抗震设计提供了依据和参考。     

Seismic performance of two-story subway station structures with different isolating systems

When the sliding bearing is fixed only at the top of the middle column of the underground structure, the cracks at the side end of the middle plate should be aggravated while the seismic damage of the mid-column should be alleviated. To enhance the seismic performance of the mid-plate, a new isolation design method has been mentioned while the elastic sliding bearings are set at the top of the mid-columns and between the side end of the mid-plate and the side wall at the same time. By establishing a nonlinear finite element analysis model for the static-dynamic coupling interaction system, the seismic response characteristics of the cast-in-place station structure without a sliding bearing have been analyzed and compared with those of the station structure with the sliding bearing fixed only at the top of the middle columns, and those of the station structure with sliding bearing be fixed between the mid-plate and the sidewall at the same time. The results show that the new isolation station structures suffer fewer earthquake damages at the mid-plate and mid-columns at the same time, which can improve the overall seismic performance of the subway station structure.