城市地下孤石的地震波场响应特征及其在地震层析成像中的应用研究

城市浅层地质条件复杂,地下孤石是常见的地质灾害体.为降低城市地下空间开发过程中孤石造成的安全隐患,本文通过构建城市地下孤石的数值模型,分析地震波在含孤石介质中的传播特征,根据地震波的传播特征利用层析成像技术对孤石进行速度成像,并通过工程实例进行验证.研究表明:(1)地震波场的数值模拟是研究地球物理勘探的有效方法,能高效...     

Flash flood mitigation as a positive consequence of anthropogenic forcing on the groundwater resource in a karst catchment

The Mediterranean coastal region is prone to high-intensity rainfall events that are frequently associated with devastating flash floods. This paper discusses the role of a karst aquifer system in the flash floods of a Mediterranean river, the Lez river. Most of the Lez river watershed is located on karst terrains where interactions between surface water and groundwater take place. During extreme rainfall events, the presence of fractures and well-developed karst features in carbonate terrains enhances the infiltration processes and involves the concentration of the recharge into highly organized and permeable flow paths. The groundwater, therefore, quickly moves towards the natural outlets of the karst system. The influence of the Lez karst aquifer system on the associated river floods dynamics is analysed while considering the spatially distributed rainfall, as well as the time series of the groundwater level within the aquifer and of the Lez river discharge measured at various gauging stations. Special attention is given to the relative importance of the surface and underground processes involved in flash flood genesis. It is shown that the karst groundwater contributes to flash floods under certain conditions, while high-rate pumping within the karst aquifer, which generates significant drawdown, may mitigate flash floods under other conditions.     

Geometry and kinematics of dunes during steady and unsteady flows in the Calamus River, Nebraska, USA

The geometry and kinematics of river dunes were studied in a reach of the Calamus River, Nebraska. During day-long surveys, dune height, length, steepness, migration rate, creation and destruction were measured concurrently with bedload transport rate, flow depth, flow velocity and bed shear stress. Within a survey, individual dune heights, lengths and migration rates were highly variable, associated with their three-dimensional geometry and changes in their shape through time. Notwithstanding this variability, there were discernible changes in mean dune height, length and migration rate in response to changing discharge over several days. Changes in mean dune height and length lagged only slightly behind changes in discharge. Therefore, during periods of both steady and unsteady flow, mean dune lengths were quite close to equilibrium values predicted by theoretical models. Mean dune steepnesses were also similar to predicted equilibrium values, except during high, falling flows when the steepness was above that predicted. Variations in mean dune height and length with discharge are similar to those predicted by theory under conditions of low mean dune excursion and discharge variation with a short high water period and long low water period. However, the calculated rates of change of height of individual dunes vary considerably from those measured. Rates of dune creation and destruction were unrelated to discharge variations, contrary to previous results. Instead, creations and destructions were apparently the result of local variations in bed shear stress and sediment transport rate. Observed changes in dune height during unsteady flows agree with theory fairly well at low bed shear stresses, but not at higher bed shear stresses when suspended sediment transport is significant.     

Inverse modeling of velocities and inferred cause of overwash that emplaced inland fields of boulders at Anegada, British Virgin Islands

A combination of numeric hydrodynamic models, a large-clast inverse sediment-transport model, and extensive field measurements were used to discriminate between a tsunami and a storm striking Anegada, BVI a few centuries ago. In total, 161 cobbles and boulders were measured ranging from 1.5 to 830?kg at distances of up to 1?km from the shoreline and 2?km from the crest of a fringing coral reef. Transported clasts are composed of low porosity limestone and were derived from outcrops in the low lying interior of Anegada. Estimates of the near-bed flow velocities required to transport the observed boulders were calculated using a simple sediment-transport model, which accounts for fluid drag, inertia, buoyancy, and lift forces on boulders and includes both sliding and overturning transport mechanisms. Estimated near-bed flow velocities are converted to depth-averaged velocities using a linear eddy viscosity model and compared with water level and depth-averaged velocity time series from high-resolution coastal inundation models. Coastal inundation models simulate overwash by the storm surge and waves of a category 5 hurricane and tsunamis from a Lisbon earthquake of M 9.0 and two hypothetical earthquakes along the North America Caribbean Plate boundary. A modeled category 5 hurricane and three simulated tsunamis were all capable of inundating the boulder fields and transporting a portion of the observed clasts, but only an earthquake of M 8.0 on a normal fault of the outer rise along the Puerto Rico Trench was found to be capable of transporting the largest clasts at their current locations. Model results show that while both storm waves and tsunamis are capable of generating velocities and temporal acceleration necessary to transport large boulders near the reef crest, attenuation of wave energy due to wave breaking and bottom friction limits the capacity of storm waves to transport large clast at great inland distances. Through sensitivity analysis, we show that even when using coefficients in the sediment-transport model which yield the lowest estimated minimum velocities for boulder transport, storm waves from a category 5 hurricane are not capable of transporting the largest boulders in the interior of Anegada. Because of the uncertainties in the modeling approach, extensive sensitivity analyses are included and limitations are discussed.     

Use of LSPIV in assessing urban flash flood vulnerability

The probability of the occurrence of urban flash floods has increased appreciably in recent years. Scientists have published various articles related to the estimation of the vulnerability of people and vehicles in urban areas resulting from flash floods. However, most published works are based on research performed using numerical models and laboratory experiments. This paper presents a novel approach that combines the implementation of image velocimetry technique (large-scale particle image velocimetry—LSPIV) using a flash flood video recorded by the public locally and the estimation of the vulnerability of people and vehicles to high water velocities in urban areas. A numerical one‐dimensional hydrodynamic model has also been used in this approach for water velocity characterization. The results presented in this paper correspond to a flash flood resulting on November 29, 2012, in the city of Asunción in Paraguay. During this flash flood, people and vehicles were observed being carried away because of high water velocities. Various sequences of the recorded flash flood video were characterized using LSPIV. The results obtained in this work validate the existing vulnerability criterion based on the effect of the flash flood and resulting high water velocities on people and vehicles.     

Experimental and theoretical study on the effect of unsteady flow on the fracturing pressure in hydraulic fracturing test

Reasonable determination of formation fracturing pressure concerns the stable operation of underground fluid injection projects. In this work, we studied the effect of unsteady flow on fracturing pressure. Hydraulic fracturing tests on low permeable sandstone were conducted with the injection rate between 0.1 and 2.0 ml/min. Then, the fracturing pressure prediction models for hollow cylinder under both unsteady flow and steady flow conditions were deduced. Finally, the effect of unsteady flow on the fracturing pressure was studied based on the experimental result and several influence factors. It was shown that fracturing pressure increased with the elevated pressurization rate in the tests, while the slope of the variation curve decreases. The model considering unsteady flow can reflect the variation tendency of fracturing pressures in experiments, while fracturing pressures from the model considering steady flow are invariant with different pressurization rates. Fracturing pressure decreases with the elevated rock permeability and increases with the elevated fluid viscosity, and these two effects are actually generated by the unsteady flow. Whether to consider the unsteady flow has no significant influence on the effect of rock tensile strength on fracturing pressure when the tensile strength is very low. However, when the tensile strength is high, the effect of unsteady flow cannot be neglected.     

MEASUREMENT OF SAND TRANSPORT IN RIVERS WITH SPECIAL REFERENCE TO TRACER METHODS

The paper describes tests made on a small river to examine the practicability of measuring quantitatively, by tracer methods, sand transport occurring over a rippled bed in a natural channel. Three principal tracer methods, all previously verified under steady flow conditions in a laboratory channel, and identified as the spatial integration, time integration and steady dilution methods, were examined by introducing radioactive and fluorescent tracers to represent a narrow size fraction of the bed sand. The serious limitations of the time integration and steady dilution methods under conditions of unsteady solids discharge were exposed by a sudden flood discharge occurring during the course of the study. An approximate method, based on the time integration approach, but also requiring some spatial tracer data, is suggested, whereby an estimate of the transport is possible in the case of highly variable sand discharge. However, it is seen that in the sand case, the spatial integration method offers the most workable and informative tracer method for use in natural rivers. The results, by this method, were generally consistent with those given by the summation of the measured suspended-solids flux and the transport occurring either by surface creep or in close vicinity to the bed. The latter quantity is determined from the downstream advance of the ripple bed-forms, and methods are described for its derivation from ultrasonic echo sounder records of bed profiles. Finally sand transport values are calculated for a number of river discharges and for two laboratory channel tests using various well-known computational procedures—Einstein, Bagnold, modified Einstein, Meyer-Peter and Muller, and Laursen procedures. Significant differences in the computed values for individual flow conditions are observed.     

渠道非恒定流水力学响应研究

输水渠道中的水位将随流量的改变而产生波动,且由一种恒定流状态过渡到另一种恒定流状态需要经过一段稳定的时间。这一水力变化过程与渠道自身的结构特性、渠道运行方式的选择和流量变化的幅度、速率有关。通过对渠道中的非恒定流进行数值模拟,对不同情况下渠道的水力学响应进行了研究。     

River Warren boulders, Minnesota, USA: catastrophic paleoflow indicators in the southern spillway of glacial Lake Agassiz

Boulders resting on meltwater-sculpted and striated-granite bedrock near the head of the southern outlet spillway of glacial Lake Agassiz are used to generate paleodischarge calculations. The rounded nature of many boulders suggests fluvial transport and a corestone origin. The distribution of boulders in clusters and linear trains records the interaction of clasts during transport and deposition. The geomorphology of the spillway with streamlined erosional remnant hills is characteristic of other large flood spillways. Using the Manning equation and a variety of empirical equations to determine paleovelocity, preferred discharges between 0.364 and 0.102 Sv are calculated. These discharges agree well with flood discharges using modeling methodologies, and most likely represent ephemeral and catastrophic flood events linked to either episodic incision at the outlet or the result of rapid inputs of meltwater to Lake Agassiz.     

Early warning system for detection of urban pluvial flooding hazard levels in an ungauged basin

Prolonged and high intensity rainfall often saturates urban drainage systems and generates urban pluvial flooding, resulting in hazardous flash floods. The city most affected by urban flooding in Colombia (South America) is Barranquilla since lacks a proper storm water drainage system. Heavy rainfall produces flash floods to quickly become torrential streams that flow down the streets endangering pedestrians. This research describes a low-cost early warning system (EWS) to detect in real time the hazard level of a stream in an ungauged basin. The EWS indicates whether it is safe or not for pedestrians to cross a flooded street, based on certain criteria used to assess the hazard level of the torrent. A hydrological and hydraulic model calculates the flow, velocity and water level in all cross sections along the stream. The model uses only real-time measurements of rain gauges and topographic survey data to determine the hazard level. Finally, a wireless sensor network sends the alert to a web platform and warns the community in real time.

     

Laboratory sustained turbidity currents form elongate ridges at channel mouths

The turbulent flow structure, suspended sediment dynamics and deposits of experimental sustained turbidity currents exiting a channel across a break in slope into a wide tank are documented. The data shed light on the flow evolution and deposit geometry of analogous natural channel‐fed submarine fans. Flows generated in a 0·3 m wide, sloping (0°, 6°, 9° or 20°) channel crossed an angular slope break and spread onto a horizontal tank floor. Flow development comprised: (i) channelized phase (unsteady channel flow developing into steady channel flow); (ii) initial lateral expansion phase (unsteady‐spreading wall jet phase); (iii) constant lateral expansion phase (steady‐spreading wall jet phase); and (iv) rapid waning phase. Phases (i) and (iv) are similar to laterally constrained turbidity currents, but phases (ii) and (iii) are considerably different from such two‐dimensional currents. Steeper channel slopes produced greater flow velocities and turbulence intensities, but these effects diminished markedly with distance from the channel mouth. Flow velocity vectors in the tank had similar patterns for all channel slopes, with a central core of faster velocity and narrow vector dispersion and slower flow with larger dispersion at the jet margins. Suspended sediment concentrations were higher within flow heads and dense basal layers in flow bodies. Time‐averaged acoustic backscatter data showed vertical concentration gradients, confirmed by siphon samples. The deposits comprised a thick central ridge, of similar order width to the channel mouth, with abrupt margins and a surrounding, very thin, fan‐like sheet. The ridge was coarser grained and better sorted than the original sediment, with grain‐size fining downstream, particularly over the fan‐like sheet. The formation of a central ridge suggests that, in the tank, vertical turbulent momentum exchange is more significant for sediment dynamics than spanwise momentum exchange due to lateral expansion. The streamwise elongate geometry of the ridge contrasts with conventional fan‐like geometry found with surge‐type turbidity flows, a result that has widespread stratigraphic and economic implications.     

Multi-hazard susceptibility mapping based on Convolutional Neural Networks

Multi-hazard susceptibility prediction is an important component of disasters risk management plan. An effective multi-hazard risk mitigation strategy includes assessing individual hazards as well as their interactions. However, with the rapid development of artificial intelligence technology, multi-hazard susceptibility prediction techniques based on machine learning has encountered a huge bottleneck. In order to effectively solve this problem, this study proposes a multi-hazard susceptibility mapping framework using the classical deep learning algorithm of Convolutional Neural Networks (CNN). First, we use historical flash flood, debris flow and landslide locations based on Google Earth images, extensive field surveys, topography, hydrology, and environmental data sets to train and validate the proposed CNN method. Next, the proposed CNN method is assessed in comparison to conventional logistic regression and k-nearest neighbor methods using several objective criteria, i.e., coefficient of determination, overall accuracy, mean absolute error and the root mean square error. Experimental results show that the CNN method outperforms the conventional machine learning algorithms in predicting probability of flash floods, debris flows and landslides. Finally, the susceptibility maps of the three hazards based on CNN are combined to create a multi-hazard susceptibility map. It can be observed from the map that 62.43% of the study area are prone to hazards, while 37.57% of the study area are harmless. In hazard-prone areas, 16.14%, 4.94% and 30.66% of the study area are susceptible to flash floods, debris flows and landslides, respectively. In terms of concurrent hazards, 0.28%, 7.11% and 3.13% of the study area are susceptible to the joint occurrence of flash floods and debris flow, debris flow and landslides, and flash floods and landslides, respectively, whereas, 0.18% of the study area is subject to all the three hazards. The results of this study can benefit engineers, disaster managers and local government officials involved in sustainable land management and disaster risk mitigation.     

Influence of temporal fluctuations and spatial heterogeneity on pollution transport in porous media

The combined influence of temporal fluctuations and spatial heterogeneity on non-reactive solute transport mechanisms in porous media can be understood by performing simulations of steady and unsteady flow and transport in heterogeneous media. The study focuses on issues such as the degree of heterogeneity, correlation length, separation of the combined effects of temporal and spatial variations, and ergodicity conditions under unsteady flow conditions. It is shown that the effect of temporal variations on solute transport is masked by the strong effect of spatial heterogeneity. There is no obvious difference in plume shape between steady and unsteady flow conditions; the first and the second spatial moments of the plume of the unsteady-state flow condition fluctuate around the steady-state flow condition with the same period of oscillations as the input signal at small storage coefficient (S????0.001). At a relatively high standard deviation in hydraulic conductivity and a small storage coefficient, the unsteady flow condition sharpens the temporal variations in macrodispersion coefficients. The magnitude of the longitudinal macrodispersion coefficient under unsteady flow condition is almost doubled at the maximum values. However, the transverse macrodispersion coefficient fluctuates around zero. The Kubo number and Peclet number ranges are 1.2?C64 and 10?C250, respectively.     

Modeling flood dynamics in a temporary river draining to an eutrophic reservoir in southeast Portugal

Enxoé is a small temporary river with a flushy regime, which flash floods carry significant loads to the reservoir. As a result, the reservoir, which supplies 25,000 inhabitants, exhibits a high trophic state and cyanobacteria blooms since its construction in 1998, with water abstractions requiring extensive treatment. This study aimed to understand the contribution of flash floods to the Enxoé’s reservoir high trophic state using a modeling approach. This was the first time the river was monitored and that a modeling study was carried out. The MOHID-Land model was implemented to assess the water path in the catchment, and was integrated with field data to compute river loads. Results confirmed the importance of flash events. During flash floods, water properties were determined by soil surface and river bottom wash out, and depended mostly on the flush sequence and intensity. Model simulations showed that soil surface permeability reduction was an important factor regulating surface runoff while soil moisture was low. The first flood after the dry period contributed to 2% of the yearly discharge, 3% of yearly N load, and 7% of the yearly P loads. Winter floods contribution differed, producing 10% of both yearly discharge and loads. However, concentration of particulate matter and organic compounds in the first flood were one order of magnitude higher than in winter floods. This was due to river bottom resuspension and erosion of riparian areas, representative dynamics of a flushy regime. During subsequent winter floods, nutrient concentrations tended to remain constant as the watershed surface and respective soils were washed. Further work should link a watershed model to a reservoir model to depict the flood impact in the reservoir, and test management strategies to reduce the reservoir trophic state.     

Finite element analysis of axisymmetric flow into partially penetrating wells

The finite element method was used to solve the problem of flow to a partially penetrating strainer well in a leaky confined aquifer under steady and unsteady conditions. Predictions of the model were compared with reported pumping test data. A good agreement was obtained between them. The effects of penetration and screen position on discharge through imperfect well were studied. Results were compared with analytical solutions and presented graphically.     

The numerical modelling of advective transport in the presence of fluid pressure transients

Conventional modelling of transport problems for porous media usually assumes that the Darcy flow velocities are steady. In certain practical situations, the flow velocity can exhibit time‐dependency, either due to the transient character of the flow process or time dependency in the boundary conditions associated with potential flow. In this paper, we consider certain one‐ and three‐dimensional problems of the advective transport of a chemical species in a fluid‐saturated porous region. In particular, the advective flow velocity is governed by the piezo‐conduction equation that takes into account the compressibilities of the pore fluid and the porous skeleton. Time‐ and/or mesh‐refining adaptive schemes used in the computational modelling are developed on the basis of a Fourier analysis, which can lead to accurate and optimal solutions for the advective transport problem with time‐ and space‐dependent advective flow velocity distributions. Copyright © 2006 John Wiley & Sons, Ltd.     

Distribution, origin and transport process of boulders deposited by the 2004 Indian Ocean tsunami at Pakarang Cape, Thailand

The tsunami of 2004 in the Indian Ocean transported thousands of meters-long boulders shoreward at Pakarang Cape, Thailand. We investigated size, position and long axis orientation of 467 boulders at the cape. Most of boulders found at the cape are well rounded, ellipsoid in shape, without sharp broken edges. They were fragments of reef rocks and their sizes were estimated to be < 14m³ (22.7t). The distribution pattern and orientation of long axis of boulders reflect the inundation pattern and behavior of the tsunami waves. It was found that there is no clear evidence indicating monotonous fine/coarse shoreward trends of these boulders along each transect line. On the other hand, the large boulders were deposited repeatedly along the three arcuate lines at the intertidal zone with a spacing of approximately 136m interval. This distribution pattern may suggest that long-lasting oscillatory flows might have repositioned the boulders and separated the big ones from small. No boulders were found on land, indicating that the hydraulic force of the tsunami wave rapidly dissipated on reaching the land due to the higher bottom friction and the presence of a steep slope. We further conducted numerical calculation of tsunami inundation at Pakarang Cape. According to the calculation, the sea receded and the major part of the tidal bench (area with boulders at present) was exposed above the sea surface before the arrival of the first tsunami wave. The first tsunami wave arrived at the cape from west to east at approximately 130min after the tsunami generation, and then inundated inlands. Our calculation shows that tsunami wave was focused around the offshore by a small cove at the reef edge and spread afterwards in a fan-like shape on the tidal bench. The critical wave velocities necessary to move the largest and average-size boulders by sliding can be estimated to be approximately 3.2 and 2.0m/s, respectively. The numerical result indicates that the maximum current velocity of the first tsunami wave was estimated to be from 8 to 15m/s between the reef edge and approximately 500m further offshore. This range is large enough for moving even the largest boulder shoreward. These suggest that the tsunami waves that were directed eastward, struck the reef rocks and coral colonies, originally located on the shallow sea bottom near the reef edge, and detached and transported the boulders shoreward.     

Potential hazard analysis from the viewpoint of flow measurement in large open-channel junctions

A robust and reliable risk assessment procedure for hydrological hazards deserves attention. The transport of mixed material in the types of hazard, such as flash floods and debris flows, plays an important role and should be taken into consideration. At present, however, the mixed material transport phenomenon is not systematically and simply included in the procedures for the elaboration of hazard analysis. The consequence is the risk of losing prediction accuracy and of underestimating hazard impacts. The phenomenon of open-channel junction flow such as in irrigation ditches and at wastewater treatment facilities is an important aspect of hydrological engineering. Extensive experimental studies of the flow characteristics have been conducted. In this study, an Acoustic Doppler Velocimeter is used to monitor the flow in large open-channel junctions at the Chihtan purification plant in Taipei. The condition of the experiment is made up of different discharge. The result of the measurement data shows that the mean velocity profile, aspect ratio, measured vertical line location and the ratio of the distance of the vertical line maximum velocity from the bed to water depth. An artificial neural network (ANN) is used to simulate the mean velocity profiles in this study. The results show that the ANN can accurately and reliably simulate the mean velocity profiles. The result is believed to provide hints for the consideration of the effects, induced by mixed material transport, in the elaboration of hazard analysis.     

用水位恢复数据反演越流承压含水层参数

水文地质勘测中,由于水位恢复阶段没有人力和机械因素干扰,由测量数据可以画出平滑的曲线,更适宜于分析水文地质参数,但经典非稳定流Theis公式和Theis水位恢复法所依据的若干假设条件在实际中难以满足。针对越流承压含水层具有补给条件而很快稳定以及非稳定抽水阶段流量常“抖动”的特点,论证了利用稳定流量和Hantush-Jacob越流模型计算的可行性,并用扩展的卡尔曼滤波器反演水文地质参数,经文献中历史数据和实地抽水试验检验结果表明,由该方法反演的预测值和实测值吻合良好。     

人工阶梯-深潭破坏案例与稳定性分析

通过2006—2010年4个修建人工阶梯-深潭系统的治理山区河流案例,总结其治理效果和最终破坏原因.以单个阶梯为分析对象,给出其受力表达式,建立单个阶梯-深潭的简化稳定性模型,进而分析来流量和冲刷角变化对其稳定性的影响.单个阶梯的稳定性取决于关键石块粒径、河道坡降、流量和冲刷角.洪水期的洪峰流量和阶梯下游冲刷是阶梯破坏的主要原因,上游来流量增加和冲刷角越大,阶梯越易发生破坏.人工阶梯-深潭系统在洪水期的稳定性是其发挥长期治理效果的关键.