首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 33 毫秒
1.
A simple and fast treatment of hydrogeologic features with irregularly shaped boundaries in two‐dimensional analytic element groundwater flow models is presented. The star domain shapes of the features are restricted to closed shapes represented as smooth and continuous single‐valued functions of distance from a focus point, . The element can be used to treat a variety of boundary and continuity conditions, including those of irregularly shaped lakes or heterogeneities in hydraulic conductivity. The new element is demonstrated via some simple illustrative test cases and shown to be efficient, accurate, and much simpler to implement than existing solutions for irregular shapes.  相似文献   

2.
The conventional velocity scan can be computationally expensive for large‐scale seismic data sets, particularly when the presence of anisotropy requires multiparameter scanning. We introduce a fast algorithm for 3D azimuthally anisotropic velocity scan by generalizing the previously proposed 2D butterfly algorithm for hyperbolic Radon transforms. To compute semblance in a two‐parameter residual moveout domain, the numerical complexity of our algorithm is roughly as opposed to of the straightforward velocity scan, with N being the representative of the number of points in a particular dimension of either data space or parameter space. Synthetic and field data examples demonstrate the superior efficiency of the proposed algorithm.  相似文献   

3.
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.  相似文献   

4.
To more accurately predict the migration behavior of pollutants in porous media, we conduct laboratory scale experiments and model simulation. Aniline (AN) is used in one-dimensional soil column experiments designed under various media and hydrodynamic conditions. The advection-dispersion equation (ADE) and the continuous-time random walk (CTRW) were used to simulate the breakthrough curves (BTCs) of the solute transport. The results show that the media and hydrodynamic conditions are two important factors affecting solute transport and are related to the degree of non-Fickian transport. The simulation results show that CTRW can more effectively describe the non-Fickian phenomenon in the solute transport process than ADE. The sensitive parameter in the CTRW simulation process is , which can reflect the degree of non-Fickian diffusion in the solute transport. Understanding the relationship of with velocity and media particle size is conducive to improving the reactive solute transport model. The results of this study provide a theoretical basis for better prediction of pollutant transport in groundwater.  相似文献   

5.
In this methods note, we present a simple analytical formula to quantify the steady‐state leakage flux (Q) over a perforated aquitard. The flux depends on the aquitard thickness (D), the radius of the perforation (R), the hydraulic conductivity of the material inside the perforation (kfill), the conductivities of the overlying and underlying aquifers (k1 and k2, respectively), and the head difference between the two aquifers (ΔH): This equation assumes an aquitard separating two homogeneous and infinite aquifers (R ? aquifer thickness) in which radial flow to and from the perforation occurs, with no other recharge or discharge boundaries near the perforation. The flux through a perforation in a hypothetical case study with D = 10 m, k1 = 10 m/d, k2 = 20 m/d, R = 0.072 m, and ΔH = 1 m ranges between less than 1 mL/d if the hole is backfilled with bentonite (kfill = 10?4 m/d), to several liters per day if the perforation is backfilled with sand from the overlying aquifer (kfill = 10 m/d), to several m3/d if the perforation forms an open conduit (kfill = 105 m/d). The leakage fluxes calculated with this model agree well with those calculated using a numerical model (MODFLOW).  相似文献   

6.
  相似文献   

7.
We introduce a simple correction to coastal heads for constant‐density groundwater flow models that contain a coastal boundary, based on previous analytical solutions for interface flow. The results demonstrate that accurate discharge to the sea in confined aquifers can be obtained by direct application of Darcy's law (for constant‐density flow) if the coastal heads are corrected to ((α + 1)/α)hs ? B/2α, in which hs is the mean sea level above the aquifer base, B is the aquifer thickness, and α is the density factor. For unconfined aquifers, the coastal head should be assigned the value . The accuracy of using these corrections is demonstrated by consistency between constant‐density Darcy's solution and variable‐density flow numerical simulations. The errors introduced by adopting two previous approaches (i.e., no correction and using the equivalent fresh water head at the middle position of the aquifer to represent the hydraulic head at the coastal boundary) are evaluated. Sensitivity analysis shows that errors in discharge to the sea could be larger than 100% for typical coastal aquifer parameter ranges. The location of observation wells relative to the toe is a key factor controlling the estimation error, as it determines the relative aquifer length of constant‐density flow relative to variable‐density flow. The coastal head correction method introduced in this study facilitates the rapid and accurate estimation of the fresh water flux from a given hydraulic head measurement and allows for an improved representation of the coastal boundary condition in regional constant‐density groundwater flow models.  相似文献   

8.
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.  相似文献   

9.
Under certain circumstances, seismic propagation within porous media may be associated to the conversion of mechanical energy to electromagnetic energy, which is known as a seismo‐electromagnetic phenomenon. The propagation of fast compressional P‐waves is more specifically associated to the manifestations of a seismoelectric field linked to the fluid flows within the pores. The analysis of seismoelectric phenomena, which requires the combination of the theory of electrokinetics and Biot's theory of poroelasticity, provides us with transfer function that links the coseismic seismoelectric field E to the seismic acceleration . To measure the transfer function, we have developed an experimental setup enabling seismoelectric laboratory observation in unconsolidated quartz sand within the kilohertz range. The investigation focused on the impact of fluid conductivity and water saturation over the coseismic seismoelectric field. During the experiment, special attention was given to the accuracy of electric field measurements. We concluded that, to obtain a reliable estimate of the electric field amplitude, the dipole from which the potential differences are measured should be of much smaller length than the wavelength of the propagating seismic field. Time‐lapse monitoring of the seismic velocities and seismoelectric transfer functions were performed during imbibition and drainage experiments. In all cases, the quantitative analysis of the seismoelectric transfer function was in good agreement with theoretical predictions. While investigating saturation variations from full to residual water saturation, we showed that the ratio undergoes a switch in polarity at a particular saturation , which also implies a sign change of the filtration, traducing a reversal of the relative fluid displacement with respect to the frame. This sign change at critical saturation stresses a particular behaviour of the poroelastic medium: the dropping of the coseismic electric field to zero traduces the absence of relative pore/fluid displacements representative of a Biot dynamically compatible medium. We concluded from our experimental study in loose sand that the measurements of the coseismic seismoelectric coupling may provide information on fluid distribution within the pores and that the reversal of the seismoelectric field may be used as an indicator of the dynamically compatible state of the medium.  相似文献   

10.
Release of nitrogen compounds into groundwater, particularly those compounds from excessive agricultural fertilization, is a major concern in an aquifer recharge. Among the nitrogen compounds, ammonium ( ) is a common one. In order to assess the risk of agricultural fertilizer contamination to an aquifer through infiltration, adsorption onto a loamy agricultural soil profile (0–0.60 m depth) was studied using a soil column experiment and modelling simulation. The soil used in the experiment was drawn from an agricultural field in Xinzhen, Fangshan district, Beijing, China, and reconstituted in laboratory soil columns. Column experiments were conducted using bromide (conservative tracer) and ‐bearing aqueous solutions. The ammonium concentrations in the soil water samples were measured, and their values were plotted as the breakthrough curves. The chemical's soil–water distribution coefficients (Kd) were calculated using breakthrough curves. Then the retardation factor (R) in saturated soil was calculated. For the ‐bearing aqueous solutions, the strongest adsorption occurred at the soil depth of 0.30–0.45 m. The convection–dispersion equation model and chemical non‐equilibrium model in Hydrus‐1D were used to simulate transport in the loamy soil. The two‐site chemical non‐equilibrium model in Hydrus‐1D was best to simulate transport through the soil column. Parameter sensitivity study was conducted to investigate the influences of solute transport by Kd, the fraction of exchange sites assuming to be in equilibrium with the solution phase (f), the longitudinal dispersivity (λ), and the first‐order rate coefficients (ω). The sensitivity analysis results indicate Kd is the most critical parameter.  相似文献   

11.
D. Markovic  M. Koch 《水文研究》2015,29(7):1806-1816
Hydrological processes commonly exhibit long‐term persistence, also known as the ‘Hurst phenomenon’. Here, we examine long‐term precipitation and streamflow time series from the Elbe River Basin to quantify differences in the spectral properties and in the Hurst parameter estimates () of the individual hydrological cycle components. Precipitation‐runoff modelling is performed for the Elbe River sub‐catchment Striegis using the Soil and Water Assessment Tool (SWAT). For 38 daily 50 years long streamflow time series from the Elbe River Basin, baseflow separation and spectral analysis is performed. The results show a spectral shift towards low‐frequency scales (>2 years) from precipitation to baseflow, with a parallel increase of from 0.52 (precipitation) to 0.65 (baseflow). The SWAT model is able to reproduce both, the main low‐frequency mode (≈7 yr.) and the (0.62) of the observed Striegis River flow time series. The baseflow appears to be the main component which shapes the low‐frequency response and of streamflow in the Elbe River Basin to the input precipitation. This conclusion is further confirmed through PMWIN‐MODFLOW groundwater modelling of a hypothetic phreatic stream‐connected aquifer system that consists of various soils (sand, loamy sand and silt). A power shift towards lower frequencies and an increase of for the hydraulic heads is obtained, as the aquifer's lateral dimensions increase and its hydraulic conductivity decreases. The average of the groundwater heads is 0.80, 0.90 and 1.0 for sand, loamy sand and silt aquifers, respectively. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

12.
Quantifying biogeochemical cycles of nitrogen (N) and the associated fluxes to surface waters remains challenging, given the need to deal with spatial and temporal variability and to characterize complex and heterogeneous landscapes. We focused our study on catchments S14 and S15 located in the Adirondack Mountains of New York, USA, which have similar topographic and hydrologic characteristics but contrasting stream nitrate ( ) concentrations. We characterized the mechanisms by which reaches the streams during hydrological events in these catchments, aiming to reconcile our field data with our conceptual model of factors that regulate nutrient exports from forested catchments. Combined hydrometric, chemical and isotopic (δ ) data showed that the relative contributions of both soil and ground water sources were similar between the two catchments. Temporal patterns of stream chemistry were markedly different between S14 and S15, however, because the water sources in the two catchments have different solute concentrations. During late summer/fall, the largest source of in S14 was till groundwater, whereas shallow soil was the largest source in S15. concentrations in surface water decreased in S14, whereas they increased in S15 because an increasing proportion of stream flow was derived from shallow soil sources. During snowmelt, the largest sources of were in the near‐surface soil in both catchments. Concentrations of increased as stream discharge increased and usually peaked before peak discharge, when shallow soil water sources made the largest contribution to stream discharge. The timing of peaks in stream concentrations was affected by antecedent moisture conditions. By elucidating the factors that affect sources and transport of N, including differences in the soil nutrient cycling and hydrological characteristics of S14 and S15, this study contributes to the overall conceptualization of release from temperate forested catchments. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

13.
Presence of noise in the acquisition of surface nuclear magnetic resonance data is inevitable. There are various types of noise, including Gaussian noise, spiky events, and harmonic noise that affect the signal quality of surface nuclear magnetic resonance measurements. In this paper, we describe an application of a two‐step noise suppression approach based on a non‐linear adaptive decomposition technique called complete ensemble empirical mode decomposition in conjunction with a statistical optimization process for enhancing the signal‐to‐noise ratio of the surface nuclear magnetic resonance signal. The filtering procedure starts with applying the complete ensemble empirical mode decomposition method to decompose the noisy surface nuclear magnetic resonance signal into a finite number of intrinsic mode functions. Afterwards, a threshold region based on de‐trended fluctuation analysis is defined to identify the noisy intrinsic mode functions, and then the no‐noise intrinsic mode functions are used to recover the partially de‐noised signal. In the second stage, we applied a statistical method based on the variance criterion to the signal obtained from the initial phase to mitigate the remaining noise. To demonstrate the functionality of the proposed strategy, the method was evaluated on an added‐noise synthetic surface nuclear magnetic resonance signal and on field data. The results show that the proposed procedure allows us to improve the signal‐to‐noise ratio significantly and, consequently, extract the signal parameters (i.e., and V0) from noisy surface nuclear magnetic resonance data efficiently.  相似文献   

14.
A simple relationship is proposed in this paper to construct damage‐based inelastic response spectra including the effect of ground motion duration that it can be used for damage control in seismic design of structures. This relation is established for three groups of ground motions with short‐duration, moderate‐duration, and long‐duration ranges. To develop the model, the duration effect is included in the cyclic ductility of structures by an energy‐based method, and then strength reduction factors are computed based on this modified ductility (named ). The strength reduction factors were calculated for 44 stiffness‐degrading oscillators having vibration periods between 0.05 and 4.0 s, four ultimate ductility capacities, and five damage levels subjected to 296 earthquake records. The results showed that ductility capacity, damage level, and ground motion duration are effective parameters in the energy dissipation of structures, which affect the spectra. The values of short‐period oscillators (e.g., low‐rise structures) under short‐duration records are generally greater than those under moderate‐duration and long‐duration records. Residual analysis has been made in terms of magnitude and distance to examine the validity of the proposed simple expression. Finally, the introduced spectra were compared with three previously published proposals. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

15.
Pre‐stack seismic data are indicative of subsurface elastic properties within the amplitude versus offset characteristic and can be used to detect elastic rock property changes caused by injection. We perform time‐lapse pre‐stack 3‐D seismic data analysis for monitoring sequestration at Cranfield. The time‐lapse amplitude differences of Cranfield datasets are found entangled with time‐shifts. To disentangle these two characters, we apply a local‐correlation‐based warping method to register the time‐lapse pre‐stack datasets, which can effectively separate the time‐shift from the time‐lapse seismic amplitude difference without changing the original amplitudes. We demonstrate the effectiveness of our registration method by evaluating the inverted elastic properties. These inverted time‐lapse elastic properties can be reliably used for monitoring plumes.  相似文献   

16.
The atmospheric chloride mass balance (CMB) method allows spatial evaluations of the average diffuse aquifer recharge by rainfall () in large and varied territories when long‐term steady conditions can be assumed. Often, the distributed average CMB variables necessary to calculate have to be estimated from the available variable‐length data series, which may be of suboptimal quality and spatial coverage. This paper explains the use of these data and the reliability of the results in continental Spain, chosen as a large and varied territory. The CMB variables have been regionalized by ordinary kriging at the same 4976 nodes of a 10 km × 10 km grid. Nodal values vary from 14 to 810 mm year–1, 90% ranging from 30 to 300 mm year–1. The recharge‐to‐precipitation ratios vary from 0.03 in low‐permeability formations and semiarid areas to 0.65 in some carbonate massifs. Integrated average results for the whole of continental Spain yield a potential aquifer recharge of 64 km3 year?1, the net recharge over permeable formations (40% of the territory) being 32 km3 year?1. Two main sources of uncertainty affecting (given by the coefficient of variation, CV), induced by the inherent natural variability of the variables (CVR) and from mapping (), have been segregated. The average CVR is 0.13 and could be improved with longer data series. The average is 0.07 and may be decreased with better data coverage. The estimates were compared with other regional and local recharge estimates, being 4% and 1% higher, respectively. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

17.
The nonlinearity of the seismic amplitude‐variation‐with‐offset response is investigated with physical modelling data. Nonlinearity in amplitude‐variation‐with‐offset becomes important in the presence of large relative changes in acoustic and elastic medium properties. A procedure for pre‐processing physical modelling reflection data is enacted on the reflection from a water‐plexiglas boundary. The resulting picked and processed amplitudes are compared with the exact solutions of the plane‐wave Zoeppritz equations, as well as approximations that are first, second, and third order in , , and . In the low angle range of 0°–20°, the third‐order plane‐wave approximation is sufficient to capture the nonlinearity of the amplitude‐variation‐with‐offset response of a liquid‐solid boundary with , , and ρ contrasts of 1485–2745 m/s, 0–1380 m/s, and 1.00–1.19 gm/cc respectively, to an accuracy value of roughly 1%. This is in contrast to the linear Aki–Richards approximation, which is in error by as much as 25% in the same angle range. Even‐order nonlinear corrective terms are observed to be primarily involved in correcting the angle dependence of , whereas the odd‐order nonlinear terms are involved in determining the absolute amplitude‐variation‐with‐offset magnitudes.  相似文献   

18.
Surface water quality can vary a lot with fluctuating discharge during a Rainfall – runoff event. This paper uses a set of hydrological and hydrochemical variables to explain concentration–discharge loops and hysteresis of ${\rm NO}_{3}^{- } $ , ${\rm NH}_{4}^{ + } $ and total suspended solids in a brook dewatering a small upland agricultural catchment in the Czech Republic. Our study is based on data collected by a continuous monitoring approach provided by an automatic ISCO sampler both from snow thawing and rainfall – runoff events. Methods of correlation, regression and principal component analysis (PCA) were employed to reveal possible relationships among the variables. For ${\rm NO}_{3}^{- } $ and ${\rm NH}_{4}^{ + } $ , we found several types of concentration–discharge loops due to the loop rotation direction and also the loop curvature shape, in mutual combinations, no matter which type of a hydrological event it was related to. PCA indicated that ${\rm NO}_{3}^{- } $ loops correlated mostly with the length of a rising hydrograph limb and with the slope of the initial phase of a falling hydrograph limb, 5‐day amount of precipitation and runoff coefficient. In case of ${\rm NH}_{4}^{ + } $ , the concentrations usually increased with elevated discharge, whereas PCA did not detect any closer linkages. For suspended solids, an unambiguous positive monotonic relationship was discovered. Although no definite pattern was found, this study showed the necessity of a continuous water quality monitoring system as an approach for capturing and understanding relationships between solute concentrations and runoff formation for tracing and modelling catchment pollution sources and describing transport processes.  相似文献   

19.
We present an approach based on local‐slope estimation for the separation of scattered surface waves from reflected body waves. The direct and scattered surface waves contain a significant amount of seismic energy. They present great challenges in land seismic data acquisition and processing, particularly in arid regions with complex near‐surface heterogeneities (e.g., dry river beds, wadis/large escarpments, and karst features). The near‐surface scattered body‐to‐surface waves, which have comparable amplitudes to reflections, can mask the seismic reflections. These difficulties, added to large amplitude direct and back‐scattered surface (Rayleigh) waves, create a major reduction in signal‐to‐noise ratio and degrade the final sub‐surface image quality. Removal of these waves can be difficult using conventional filtering methods, such as an filter, without distorting the reflected signal. The filtering algorithm we present is based on predicting the spatially varying slope of the noise, using steerable filters, and separating the signal and noise components by applying a directional nonlinear filter oriented toward the noise direction to predict the noise and then subtract it from the data. The slope estimation step using steerable filters is very efficient. It requires only a linear combination of a set of basis filters at fixed orientation to synthesize an image filtered at an arbitrary orientation. We apply our filtering approach to simulated data as well as to seismic data recorded in the field to suppress the scattered surface waves from reflected body waves, and we demonstrate its superiority over conventional techniques in signal preservation and noise suppression.  相似文献   

20.
Soils in post‐wildfire environments are often characterized by a low infiltration capacity with a high degree of spatial heterogeneity relative to unburned areas. Debris flows are frequently initiated by run‐off in recently burned steeplands, making it critical to develop and test methods for incorporating spatial variability in infiltration capacity into hydrologic models. We use Monte Carlo simulations of run‐off generation over a soil with a spatially heterogenous saturated hydraulic conductivity (Ks) to derive an expression for an aerially averaged saturated hydraulic conductivity ( ) that depends on the rainfall rate, the statistical properties of Ks, and the spatial correlation length scale associated with Ks. The proposed method for determining is tested by simulating run‐off on synthetic topography over a wide range of spatial scales. Results provide a simplified expression for an effective saturated hydraulic conductivity that can be used to relate a distribution of small‐scale Ks measurements to infiltration and run‐off generation over larger spatial scales. Finally, we use a hydrologic model based on to simulate run‐off and debris flow initiation at a recently burned catchment in the Santa Ana Mountains, CA, USA, and compare results to those obtained using an infiltration model based on the Soil Conservation Service Curve Number.  相似文献   

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

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