Interaction of multiple courses of wave‐induced fluid flow in layered porous media

Different theoretical and laboratory studies on the propagation of elastic waves in layered hydrocarbon reservoir have shown characteristic velocity dispersion and attenuation of seismic waves. The wave‐induced fluid flow between mesoscopic‐scale heterogeneities (larger than the pore size but smaller than the predominant wavelengths) is the most important cause of attenuation for frequencies below 1 kHz. Most studies on mesoscopic wave‐induced fluid flow in the seismic frequency band are based on the representative elementary volume, which does not consider interaction of fluid flow due to the symmetrical structure of representative elementary volume. However, in strongly heterogeneous media with unsymmetrical structures, different courses of wave‐induced fluid flow may lead to the interaction of the fluid flux in the seismic band; this has not yet been explored. This paper analyses the interaction of different courses of wave‐induced fluid flow in layered porous media. We apply a one‐dimensional finite‐element numerical creep test based on Biot's theory of consolidation to obtain the fluid flux in the frequency domain. The characteristic frequency of the fluid flux and the strain rate tensor are introduced to characterise the interaction of different courses of fluid flux. We also compare the behaviours of characteristic frequencies and the strain rate tensor on two scales: the local scale and the global scale. It is shown that, at the local scale, the interaction between different courses of fluid flux is a dynamic process, and the weak fluid flux and corresponding characteristic frequencies contain detailed information about the interaction of the fluid flux. At the global scale, the averaged strain rate tensor can facilitate the identification of the interaction degree of the fluid flux for the porous medium with a random distribution of mesoscopic heterogeneities, and the characteristic frequency of the fluid flux is potentially related to that of the peak attenuation. The results are helpful for the prediction of the distribution of oil–gas patches based on the statistical properties of phase velocities and attenuation in layered porous media with random disorder.     

Surface‐wave inversion for a P‐velocity profile with a constant depth gradient of the squared slowness

Surface waves are often used to estimate a near‐surface shear‐velocity profile. The inverse problem is solved for the locally one‐dimensional problem of a set of homogeneous horizontal elastic layers. The result is a set of shear velocities, one for each layer. To obtain a P‐wave velocity profile, the P‐guided waves should be included in the inversion scheme. As an alternative to a multi‐layered model, we consider a simple smooth acoustic constant‐density velocity model, which has a negative constant vertical depth gradient of the squared P‐wave slowness and is bounded by a free surface at the top and a homogeneous half‐space at the bottom. The exact solution involves Airy functions and provides an analytical expression for the dispersion equation. If the ratio is sufficiently small, the dispersion curves can be picked from the seismic data and inverted for the continuous P‐wave velocity profile. The potential advantages of our model are its low computational cost and the fact that the result can serve as a smooth starting model for full‐waveform inversion. For the latter, a smooth initial model is often preferred over a rough one. We test the inversion approach on synthetic elastic data computed for a single‐layer P‐wave model and on field data, both with a small ratio. We find that a single‐layer model can recover either the shallow or deeper part of the profile but not both, when compared with the result of a multi‐layer inversion that we use as a reference. An extension of our analytic model to two layers above a homogeneous half‐space, each with a constant vertical gradient of the squared P‐wave slowness and connected in a continuous manner, improves the fit of the picked dispersion curves. The resulting profile resembles a smooth approximation of the multi‐layered one but contains, of course, less detail. As it turns out, our method does not degrade as gracefully as, for instance, diving‐wave tomography, and we can only hope to fit a subset of the dispersion curves. Therefore, the applicability of the method is limited to cases where the ratio is small and the profile is sufficiently simple. A further extension of the two‐layer model to more layers, each with a constant depth gradient of the squared slowness, might improve the fit of the modal structure but at an increased cost.     

The impact of surface‐wave attenuation on 3‐D seismic survey design

Three‐dimensional seismic survey design should provide an acquisition geometry that enables imaging and amplitude‐versus‐offset applications of target reflectors with sufficient data quality under given economical and operational constraints. However, in land or shallow‐water environments, surface waves are often dominant in the seismic data. The effectiveness of surface‐wave separation or attenuation significantly affects the quality of the final result. Therefore, the need for surface‐wave attenuation imposes additional constraints on the acquisition geometry. Recently, we have proposed a method for surface‐wave attenuation that can better deal with aliased seismic data than classic methods such as slowness/velocity‐based filtering. Here, we investigate how surface‐wave attenuation affects the selection of survey parameters and the resulting data quality. To quantify the latter, we introduce a measure that represents the estimated signal‐to‐noise ratio between the desired subsurface signal and the surface waves that are deemed to be noise. In a case study, we applied surface‐wave attenuation and signal‐to‐noise ratio estimation to several data sets with different survey parameters. The spatial sampling intervals of the basic subset are the survey parameters that affect the performance of surface‐wave attenuation methods the most. Finer spatial sampling will reduce aliasing and make surface‐wave attenuation easier, resulting in better data quality until no further improvement is obtained. We observed this behaviour as a main trend that levels off at increasingly denser sampling. With our method, this trend curve lies at a considerably higher signal‐to‐noise ratio than with a classic filtering method. This means that we can obtain a much better data quality for given survey effort or the same data quality as with a conventional method at a lower cost.     

Transient three‐dimensional‐patterned ripples appear during deformation of a two‐dimensional ripple field under wave‐induced oscillatory flow

Wave‐flume experiments on the process of deformation of pre‐existing two‐dimensional ripples revealed that one of two distinct types of three‐dimensional patterned ripple, the rhombic or the barrel type, emerged in the early stages of the experimental runs. Rhombic‐type ripples were characterized by a staggered pattern of similar‐sized mounds, whereas the barrel‐type ripples had a barrel‐like plan shape arranged in echelons. The former ripples were formed under strongly asymmetrical flow conditions with large excursion distances of water particle motion near the pre‐existing rippled bed, whereas the latter were formed under relatively symmetrical flow conditions with small excursion distances. The emergence of both types of ripples was also highly dependent on sediment grain size. Copyright © 2009 John Wiley & Sons, Ltd.     

Kinematic response of single piles to vertically incident P‐waves

Dynamic response of single piles to seismic waves is fundamentally different from the free‐field motion because of the interaction between the pile and the surrounding soil. Considering soil–pile interaction, this paper presents a new displacement model for the steady‐state kinematic response of single piles to vertically incident P‐waves on the basis of a continuum model. The governing equations and boundary conditions of the two undetermined functions in the model are obtained to be coupled by using Hamilton's principle. Then, the two unknown functions are decoupled and solved by an iterative algorithm numerically. A parametric study is performed to investigate the effects of the properties of the soil–pile system on the kinematic response of single piles. It is shown that the effects of the pile–soil modulus ratio, the slenderness ratio of the pile, and the frequency of the incident excitations are very significant. By contrast, the influence of soil damping on the kinematics of the system is slight and can be neglected. Copyright © 2013 John Wiley & Sons, Ltd.     

Evaluation of basin storage–discharge sensitivity in Taiwan using low‐flow recession analysis

Sensitivity analysis of the hydrological behaviour of basins has mainly focused on the correlation between streamflow and climate, ignoring the uncertainty of future climate and not utilizing complex hydrological models. However, groundwater storage is affected by climatic change and human activities. The streamflow of many basins is primarily sourced from the natural discharge of aquifers in upstream regions. The correlation between streamflow and groundwater storage has not been thoroughly discussed. In this study, the storage–discharge sensitivity of 22 basins in Taiwan was investigated by means of daily streamflow and rainfall data obtained over more than 30 years. The relationship between storage and discharge variance was evaluated using low‐flow recession analysis and a water balance equation that ignores the influence of rainfall and evapotranspiration. Based on the obtained storage–discharge sensitivity, this study explored whether the water storage and discharge behaviour of the studied basins is susceptible to climate change or human activities and discusses the regional differences in storage–discharge sensitivity. The results showed that the average storage–discharge sensitivities were 0.056 and 0.162 mm^?1 in the northern and southern regions of Taiwan, respectively. In the central and eastern regions, the values were both 0.020 mm^?1. The storage–discharge sensitivity was very high in the southern region. The regional differences in storage–discharge sensitivity with similar climate conditions are primarily due to differences in aquifer properties. Based on the recession curve, other factors responsible for these differences include land utilization, land coverage, and rainfall patterns during dry and wet seasons. These factors lead to differences in groundwater recharge and thus to regional differences in storage–discharge sensitivity.     

Block‐and‐ash flow deposit of the Narcondam Volcano: Product of dacite–andesite dome collapse in the Burma–Java subduction complex

Narcondam Island in the Andaman Sea represents a dacite–andesite dome volcano in the volcanic chain of the Burma–Java subduction complex. The pyroclasts of andesitic composition are restricted to the periphery of the dome predominantly in the form of block‐and‐ash deposits and minor base surge deposits. Besides pyroclastic deposits, andesitic lava occurs dominantly at the basal part of the dome whereas dacitic lava occupies the central part of the dome. The pyroclasts are represented by non‐vesiculated to poorly vesiculated blocks of andesite, lapilli, and ash. The hot debris derived from dome collapse was deposited initially as massive to reversely‐graded beds with the grain support at the lower part and matrix support at the upper part. This sequence is overlain by repetitive beds of lapilli breccia to tuff breccia. These deposits are recognized as a basal avalanche rather than lahar deposit. This basal avalanche was punctuated by an ash‐cloud surge deposit representing a sequence of thinly bedded units of normal graded unit to parallel laminated beds.     

Impacts of uncertain river flow data on rainfall‐runoff model calibration and discharge predictions

In order to quantify total error affecting hydrological models and predictions, we must explicitly recognize errors in input data, model structure, model parameters and validation data. This paper tackles the last of these: errors in discharge measurements used to calibrate a rainfall‐runoff model, caused by stage–discharge rating‐curve uncertainty. This uncertainty may be due to several combined sources, including errors in stage and velocity measurements during individual gaugings, assumptions regarding a particular form of stage–discharge relationship, extrapolation of the stage–discharge relationship beyond the maximum gauging, and cross‐section change due to vegetation growth and/or bed movement. A methodology is presented to systematically assess and quantify the uncertainty in discharge measurements due to all of these sources. For a given stage measurement, a complete PDF of true discharge is estimated. Consequently, new model calibration techniques can be introduced to explicitly account for the discharge error distribution. The method is demonstrated for a gravel‐bed river in New Zealand, where all the above uncertainty sources can be identified, including significant uncertainty in cross‐section form due to scour and re‐deposition of sediment. Results show that rigorous consideration of uncertainty in flow data results in significant improvement of the model's ability to predict the observed flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparing the performance of empirical black‐box models for river flow forecasting in the Heihe River Basin,Northwestern China

For many practical reasons, the empirical black‐box models have become an increasingly popular modelling tool for river flow forecasting, especially in mountainous areas where very few meteorological observatories exist. In this article, precipitation data are used as the only input to estimate river flow. Using five empirical black‐box models—the simple linear model, the linear perturbation model, the linearly varying gain factor model, the constrained nonlinear system model and the nonlinear perturbation model–antecedent precipitation index—modelling results are compared with actual results in three catchments within the Heihe River Basin. The linearly varying gain factor model and the nonlinear perturbation model yielded excellent predictions. For better simulation accuracy, a commonly used multilayer feed‐forward neural network model (NNM) was applied to incorporate the outputs of the individual models. Comparing the performance of these models, it was found that the best results were obtained from the NNM model. The results also suggest that more reliable and precise predictions of river flow can be obtained by using the NNM model while also incorporating the combined outputs of different empirical black‐box models. Copyright © 2012 John Wiley & Sons, Ltd.     

Impact of lateral flow and spatial scaling on the simulation of semi‐arid urban land surfaces in an integrated hydrologic and land surface model

Understanding and representing hydrologic fluxes in the urban environment is challenging because of fine scale land cover heterogeneity and lack of coherent scaling relationships. Here, the impact of urban land cover heterogeneity, scale, and configuration on the hydrologic and surface energy budget (SEB) is assessed using an integrated, coupled land surface/hydrologic model at high spatial resolutions. Archetypes of urban land cover are simulated at varying resolutions using both the National Land Cover Database (NLCD; 30 m) and an ultra high‐resolution land cover dataset (0.6 m). The analysis shows that the impact of highly organized, yet heterogeneous, land cover typical of the urban domain can cause large variations in hydrologic and energy fluxes within areas of similar land cover. The lateral flow processes that occur within each simulation create variations in overland flow of up to ±200% and ±4% in evapotranspiration. The impact on the SEB is smaller and largely restricted to the wet season for our semi‐arid forcing scenarios. Finally, we find that this seasonal bias, predominantly caused by lateral flow, is displaced by a systematic diurnal bias at coarser resolutions caused by deficiencies in the method used for scaling of land surface and hydrologic parameters. As a result of this research, we have produced land surface parameters for the widely used NLCD urban land cover types. This work illustrates the impact of processes that remain unrepresented in traditional high‐resolutions land surface models and how they may affect results and uncertainty in modeling of local water resources and climate. Copyright © 2015 John Wiley & Sons, Ltd.