首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 893 毫秒
1.
2.
3.
4.
5.
6.
Concentration–discharge (CQ) relationships are widely used to assess the link between hydrological and biogeochemical processes at the catchment scale. CQ relationships are mainly calibrated using mono-objective methods to represent, either concentrations or discharge-weighted concentrations (i.e., load). Based on its wide use in hydrological modelling, we test a multi-objective calibration for the CQ relationship parameters, using both concentration and load, and compare it to a mono-objective calibration applied on either concentrations or load. This work is carried out on a high-frequency dataset (ORACLE-Orgeval Observatory, France). Our findings show that the multi-objective calibration yield a better representation of CQ relationships parameters during the high and low-flow events. The multi-objective calibration can be used for all forms of CQ relationships and avoids issues of under-representation of dilution processes characterized by high-discharge, low-concentration periods.  相似文献   

7.
Nearshore sandbars are important features in the surf zone of many beaches because they strongly influence the mean circulation and evolving morphology. Due to variations in wave conditions, sandbars can experience cross-shore migration and vary in shape from alongshore uniform (shore-parallel) to alongshore rhythmic (crescentic). Sandbar dynamics have been studied extensively, but existing observational studies usually do not quantify the processes leading to crescentic bar formation and straightening. This study analyses the dynamics of crescentic bar events at the fetch-limited beach of Castelldefels (northwestern Mediterranean Sea, Spain) using 7.5 years of hourly time-exposure video images and detailed wave conditions. The results show that, despite the generally calm wave conditions, the sandbars were very dynamic in the cross-shore and longshore directions. They often migrated rapidly offshore during storms (up to 70 m in one day) and more slowly onshore during post-storm conditions. Crescentic bars were often present at the study site (48% of the time), but only when the sandbar was at least 10 m from the shoreline. They displayed a large variability in wavelengths (100–700 m), alongshore migration speeds (0–50 m/day) and cross-shore amplitudes (5–20 m). Wavelengths increased for larger bar–shoreline distances and the alongshore migration speeds were strongly correlated with the alongshore component of the radiation stresses. Crescentic patterns typically developed during low–medium energetic waves with limited obliquity ( θ20° at 10 m depth), while bar straightening occurred during medium–high energetic waves with strong oblique angles of incidence ( θ15°). Overall, this study provides further proof for the important role of wave direction in crescentic bar dynamics and highlights the strong dependence of crescentic bar development on the initial bathymetric configuration.  相似文献   

8.
9.
To this day, field capacity (FC) is rarely defined in the context of soil properties, and the use of non-physical simplistic models is the common way to normalize water content at FC. In this study, the problem of water drainage redistribution in a soil column with and without the presence of evaporation (EV) was extensively studied. Analytical solutions for the Richards equation were established for the case of water drainage redistribution through a deeply wetted soil water column with and without EV at FC conditions. Water retention and depth evolution curves were plotted first, using different EV values of (2 mmday, 5 mmday and 8 mmday) and second, for different drainage redistribution durations of (1 day, 4 days and 6 days) where EV was set to zero for the case with no EV or to a fixed value of 5 mmday for the case with EV. The results suggest that EV plays a significant role in soil water drainage suggesting that, in the presence of EV, the FC drying front reaches much higher depths in the soil water profile than if EV is turned off. It was also concluded that FC reaches deeper depths faster the stronger EV is acting at the surface of a soil water column. Additionally, the results suggest that the texture of the soil receiving drainage controls the amount of water available for EV and as a result, EV was found to play a stronger role the smaller the hydraulic conductivity of the soil is.  相似文献   

10.
11.
Forest fires and post-fire practices influence sediment connectivity (SC). In this study, we use the ‘aggregated index of connectivity’ (AIC) to assess SC in five Mediterranean catchments (198–1090 ha) affected by a wildfire in 2012 in south-eastern Spain. Two temporal scenarios were considered, immediately after the fire and before post-fire management, and 2 years after the fire including all practices (hillslope barriers, check-dams, afforestation, salvage logging and skid trails). One LiDAR (light detection and ranging)-derived digital elevation model (DEM, 2 m × 2 m resolution) was generated, per scenario. The five catchment outlets were established as the computation target (AICOUT), and structural and functional SC were calculated. Index outputs were normalized to make the results of the non-nested catchments comparable (AICN-OUT). The output analysis includes the SC distribution along the catchments and at local scale (929 sub-catchments, 677 in the burned area), the hillslope and channel measures' effect on SC, and a sedimentological analysis using observed area-specific sediment yield (SSY) at 10 new (built after post-fire practices) concrete check-dams located in the catchments (SSY = 1.94 Mg ha−1 yr−1; σ = 1.22). The catchments with more circular shapes and steeper slopes were those with higher AICN-OUT. The structural SC maps – removing the rainfall erosivity influence – allowed evaluating the actual role played by the post-fire practices that reduced SC ( x¯= − 1.19%; σ = 0.41); while functional SC was linked to the actual change of SC ( x¯= + 5.32%; σ = 0.62). Hillslope treatments resulted in significant changes on AICN-OUT at sub-catchment scale with certain disconnectivity. A good and positive correlation was found between the SSY and the changes of AICN-OUT. However, the coarse DEM resolution explained the lack of effect of the rock check-dams – located on the secondary channels – on AICN-OUT. AICN-OUT proved to be a useful tool for decision making in post-fire restoration, but an optimal input data is still necessary to refine calculations.  相似文献   

12.
We investigate the dependence of quasi P-wave phase velocity propagating in orthotropic media on particular elasticity parameters. Specifically, due to mathematical facilitation, we consider the squared-velocity difference, , resulted from propagation in two mutually perpendicular symmetry planes. In the context of the effective medium theory, may be viewed as a parameter evaluating the influence of cracks – embedded in the background medium – parallel to one or both aforementioned planes. Our investigation is both theoretical and numerical. Based on Christoffel's equations, we propose two accurate approximations of . Due to them, we interpret the aforementioned squared-velocity difference as being twice more dependent on , than on . To describe the magnitude of the dependence, we consider the proportions between the partial derivatives of . Further, it occurs that is influenced by the ratio of vertically propagating quasi P-wave to vertically propagating quasi S-wave. Anomalously high might be caused by the low P/S ratio, which in turn can be an indicator of the presence of gas in natural fractures or aligned porosity. Also, we carry out numerical sensitivity study, according to which is approximately twice more dependent on than on , twice more sensitive to than to , and equally dependent on as on . The dependence on and can be neglected, especially for small phase angles. We verify the approximations and perform the sensitivity study, using eight examples of the elasticity tensors.  相似文献   

13.
Logistic regression, neural networks and support vector machines are tested for their effectiveness in isolating surface waves in seismic shot records. To distinguish surface waves from other arrivals, we train the algorithms on three distinguishing features of surface-wave dispersion curves in the domain: spectrum coherency of the trace's magnitude spectrum, local dip and the frequency range for a fixed wavenumber k in the spectrum. Numerical tests on synthetic data show that the kernel-based support vector machines algorithm gives the highest accuracy in predicting the surface-wave window in the domain compared to neural networks and logistic regression. This window is also used to automatically pick the fundamental dispersion curve. The other two methods correctly pick the low-frequency part of the dispersion curve but fail at higher frequencies where there is interference with higher-order modes.  相似文献   

14.
Full waveform inversion in transversely isotropic media with a vertical symmetry axis provides an opportunity to better match the data at the near and far offsets. However, multi-parameter full waveform inversion, in general, suffers from serious cycle-skipping and trade-off problems. Reflection waveform inversion can help us recover a background model by projecting the residuals of the reflected wavefield along the reflection wavepath. Thus, we extend reflection waveform inversion to acoustic transversely isotropic media with a vertical symmetry axis utilizing the proper parameterization for reduced parameter trade-off. From a radiation patterns analysis, an acoustic transversely isotropic media with a vertical symmetry axis is better described by a combination of the normal-moveout velocity and the anisotropic parameters η and δ for reflection waveform inversion applications. We design a three-stage inversion strategy to construct the optimal resulting model. In the first stage, we only invert for the background by matching the simulated reflected wavefield from the perturbations of and δ with the observed reflected wavefield. In the second stage, the background and η are optimized simultaneously and the far-offset reflected wavefield mainly contribute to their updates. We perform Born modelling to compute the reflected wavefield for the two stages of reflection waveform inversion. In the third stage, we perform full waveform inversion for the acoustic transversely isotropic media with a vertical symmetry axis to delineate the high-wavenumber structures. For this stage, the medium is described by a combination of the horizontal velocity , η and ε instead of , η and δ. The acoustic multi-parameter full waveform inversion utilizes the diving waves to improve the background as well as utilizes reflection for high-resolution information. Finally, we test our inversion algorithm on the modified Sigsbee 2A model (a salt free part) and a two-dimensional line from a three-dimensional ocean bottom cable dataset. The results demonstrate that the proposed reflection waveform inversion approach can recover the background model for acoustic transversely isotropic media with a vertical symmetry axis starting from an isotropic model. This recovered background model can mitigate the cycle skipping of full waveform inversion and help the inversion recover higher resolution structures.  相似文献   

15.
Image-based methods have compelling, demonstrated potential for characterizing flow fields in rivers, but algorithms like particle image velocimetry (PIV) must be further tested and improved to enable more effective use of these techniques. This paper presents a framework designed for this exact purpose: Simulating Hydraulics and Images for Velocimetry Evaluation and Refinement (SHIVER). The approach involves coupling a hydrodynamic model with a synthetic particle generator to advect particles between frames, as dictated by local velocity vectors and thus construct a plausible image sequence specific to the reach of interest. The resulting time series can then be used as input to a velocimetry algorithm to compare image-derived estimates with known (modelled) velocities to perform an exhaustive, spatially distributed accuracy assessment. As an example application of SHIVER, we examined the effects of interrogation area (IA) size, frame rate, flow velocity, and image sequence duration on the performance of a standard PIV algorithm. This analysis indicated that image-derived velocities were generally in close agreement with those from the flow model (root mean square error <10% and mean bias <3%), except when small IAs were coupled with low frame rates. Velocity estimates were most accurate for the lowest modelled discharge ( at baseflow) and became less reliable as the mean flow velocity increased ( for an intermediate discharge and at bankfull). Accuracy was essentially independent of image sequence duration, implying that long occupations might not be necessary. Errors were concentrated along channel margins, where PIV-based velocities tended to be greater than those from the flow model. Small IAs led to underpredictions of velocity, while larger IAs led to overpredictions. SHIVER is highly modular and could be updated to make use of different hydrodynamic models or image simulators. The framework could also facilitate more thorough sensitivity analyses and comparison of various velocimetry algorithms.  相似文献   

16.
Ill-posedness is one of the most common and intractable issues that arise when solving geophysical inverse problems. Ill-posedness could be induced by various factors such as noise, band-limited intrinsic property of seismic data and inappropriate forward operators. Regularization has been proven to be an effective method widely accepted for mitigating the adverse effects of ill-posedness. Aiming to improve the stability and fidelity of the pre-stack seismic inversion process, we implement the inversion in a Bayesian framework, with a logarithmic absolute criterion taken as a likelihood function, and an -norm metric as a priori constraint. Here, we exploit the linear approximation as the forward operator, and optimize the regularized misfit function by the alternating direction method of multipliers. Applications of the method to synthetic and real data sets yielded improved inversion results in terms of accuracy and resolution, and demonstrated the robustness of the method to noise.  相似文献   

17.
Tensor algebra provides a robust framework for multi-dimensional seismic data processing. A low-rank tensor can represent a noise-free seismic data volume. Additive random noise will increase the rank of the tensor. Hence, tensor rank-reduction techniques can be used to filter random noise. Our filtering method adopts the Candecomp/Parafac decomposition to approximates a N-dimensional seismic data volume via the superposition of rank-one tensors. Similar to the singular value decomposition for matrices, a low-rank Candecomp/Parafac decomposition can capture the signal and exclude random noise in situations where a low-rank tensor can represent the ideal noise-free seismic volume. The alternating least squares method is adopted to compute the Candecomp/Parafac decomposition with a provided target rank. This method involves solving a series of highly over-determined linear least-squares subproblems. To improve the efficiency of the alternating least squares algorithm, we uniformly randomly sample equations of the linear least-squares subproblems to reduce the size of the problem significantly. The computational overhead is further reduced by avoiding unfolding and folding large dense tensors. We investigate the applicability of the randomized Candecomp/Parafac decomposition for incoherent noise attenuation via experiments conducted on a synthetic dataset and field data seismic volumes. We also compare the proposed algorithm (randomized Candecomp/Parafac decomposition) against multi-dimensional singular spectrum analysis and classical prediction filtering. We conclude the proposed approach can achieve slightly better denoising performance in terms of signal-to-noise ratio enhancement than traditional methods, but with a less computational cost.  相似文献   

18.
We have derived a convergent scattering series solution for the frequency-domain wave equation in acoustic media with variable density and velocity. The convergent scattering series solution is based on the homotopy analysis of a vectorial integral equation of the Lippmann–Schwinger type. By using the Green's function and partial integration, we have derived the vectorial integral equation of the Lippmann–Schwinger type that involves the pressure gradient field as well as the pressure field from the wave equation. The vectorial Lippmann–Schwinger equation can in principle be solved via matrix inversion, but the computational cost of matrix inversion scales like N 3 , where N is the number of grid blocks. The computational cost can be significantly reduced if one solves the vectorial Lippmann–Schwinger equation iteratively. A simple iterative solution is the Born series, but it is only convergent when the scattering potential is sufficiently small. In this study, we have used the so-called homotopy analysis method to derive an iterative solution for the vectorial Lippmann–Schwinger equation which can be made convergent even in strongly scattering media. The computational cost of our convergent scattering series scales as N 2 . Our algorithm, which is based on the homotopy analysis method, involves a convergence control operator that we select using hierarchical matrices. We use a three-layer model and a resampled version of the SEG/EAGE salt model to show the performance of the developed convergent scattering series.  相似文献   

19.
Our objective was to discover the effect of variations in fluid properties and fracture geometry on the velocity of seismic wave propagation in fluid‐saturated media with parallel planar fractures. We used numerical models calculated by analytical solutions to examine the behaviour of P‐wave phase velocity dispersion in the normal direction to layering, in non‐porous and porous media with planar fractures. We also examined the anisotropy of low frequency phase and group velocities of fast and slow P‐waves and angular‐dependent reflection coefficients in media with planar fractures, under conditions of saturation by fluids with varying bulk moduli, densities, and fracture apertures. We defined several parametre , , and characterising dispersion, characterising anisotropy, characterising the difference between fast and slow modes, and R0 and characterising reflection. Our results show that the behaviour of dispersion shows wider stopbands in the case of gas saturation. Concavity indicator of dispersion for gas saturation was greater than that for liquid saturation and is usually greater than one. Anisotropy is more sensitive to bulk modulus contrast than to density contrast between the solid and the fluid, and is more sensitive to density contrast than to bulk modulus contrast. The case of gas saturation usually had a greater negative R0 and a greater value of compared with those of brine and heavy and light oil saturations. Our results are helpful in distinguishing fluid types saturating geophysical fractures and estimating the aperture and spacing of planar fractures. In seismic exploration, bulk modulus and fluid density can provide useful information in distinguishing among brine, oil, and gas; fracture geometry is important to estimate the permeability of reservoirs.  相似文献   

20.
While seismic reflection amplitudes are generally determined by real acoustical impedance contrasts, there has been recent interest in reflections due to contrasts in seismic‐Q. Herein we compare theoretical and modelled seismic reflection amplitudes for two different cases of material contrasts. In case A, we examine reflections from material interfaces that have a large contrast in real‐valued impedance () with virtually no contrast in seismic‐Q. In case B, we examine reflections from material interfaces that have virtually no contrast in but that have very large seismic‐Q contrasts. The complex‐valued reflection coefficient formula predicts non‐zero seismic reflection amplitudes for both cases. We choose physical materials that typify the physics of both case A and case B. Physical modelling experiments show significantly large reflections for both cases – with the reflections in the two cases being phase shifted with respect to each other, as predicted theoretically. While these modelling experiments show the existence of reflections that are predicted by theory, there are still intriguing questions regarding the size of the Q‐contrast reflections, the existence of large Q‐contrast reflections in reservoir rocks and the possible application of Q‐reflection analysis to viscosity estimation in heavy oilfields.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号