排序方式： 共有7条查询结果，搜索用时 15 毫秒
1
1.
秦皇岛市洋河、戴河滨海平原海水入侵的控制与治理 总被引：5，自引：0，他引：5
秦皇岛市洋河,戴河冲洪积平原,近20年来抽取地下水量迅速增加,80年代初形成了低于海平面的降落漏斗,造成海水从海岸线向内陆和从河流下游河床向两侧同时侵入孔隙含水层。由于水质恶化,北戴河区水源地面临报废的危险。为了保护和开发淡水资源,控制海水入侵,必须对地下水的水质和水量进行综合评价。采用水动力弥散型地下水水质模型来描述滨海含水层中海水和淡水的运动,把氯离子作为水质模型的变量。在水量模型中采用孔隙压力代替水头作为变量。数学模型用有限单元法求解。模拟面积约180km2,剖分为321个单元,184个结点。作为水文地质实体的模型,它被用来模拟和预报不同降水量年份的地下水水位和水质。计算结果表明,如果保持目前的地下水开采量,海水入侵还会深入。为了控制和治理海水入侵,本文探讨了可行措施。根据数学模型模拟的结果,在洋河下游人海口修建防潮闸,同时加强滨海含水层的管理,具有较好的经济效益和环境效益。 相似文献
2.
3.
Dispersive mass transport processes in naturally heterogeneous geological formations (porous media) are investigated based on a particle approach to mass transport and on its numerical implementation using LPT3D, a Lagrangian Particle Tracking 3D code. We are currently using this approach for studying microscale and macroscale space–time behavior (advection, diffusion, dispersion) of tracer plumes, solutes, or miscible fluids, in 1,2,3-dimensional heterogeneous and anisotropic subsurface formations (aquifers, petroleum reservoirs). Our analyses are based on a general advection-diffusion model and numerical scheme where concentrations and fluxes are discretized in terms of particles. The advection-diffusion theory is presented in a probabilistic framework, and in particular, a numerical analysis is developed for the case of advective transport and rotational flows (numerical stability of the explicit Euler scheme). The remainder of the paper is devoted to the behavior of concentration, mass flux density, and statistical moments of the transported tracer plume in the case of heterogeneous steady flow fields, where macroscale dispersion occurs due to geologic heterogeneity and stratification. We focus on the case of perfectly stratified or multilayered media, obtained by generating many horizontal layers with a purely random transverse distribution of permeability and horizontal velocity. In this case, we calculate explicitly the exact mass concentration field C(x,
t), mass flux density field f(x, t), and moments. This includes spatial moments and dispersion variance 2
x
(t) on a finite domain L, and temporal moments on a finite time scale T, e.g., the mass variance of arrival times 2
T
(x). The moments are related to flux concentrations in a way that takes explicitly into account finite space–time scales of analysis (time-dependent tracer mass; spatially variable flow through mass). The multilayered model problem is then used in numerical experiments for testing different ways of recovering information on tracer plume migration, dispersion, concentration and flux fields. Our analyses rely on a probabilistic interpretation that emerges naturally from the particle approach; it is based on spatial moments (particle positions), temporal moments (mass weighted arrival times), and probability densities (both concentrations and fluxes). Finally, as an alternative to direct estimations of the flux and concentration fields, we formulate and study the Moment Inverse Problem. Solving the MIP yields an indirect method for estimating the space–time distribution of flux concentrations based on observed or estimated moments of the plume. The moments may be estimated from field measurements, or numerically computed by particle tracking as we do here. 相似文献
4.
The asymptotic behavior of the solute velocity and dispersivity for a system of parallel fractures with matrix diffusion is
made using numerical modeling and theoretical analyses. The study is limited to linearly sorbing solutes with a constant continuous
source boundary condition. Expressions are provided for solute velocity and effective dispersivity in terms of fracture porosity
during asymptotic stage using spatial moment analyses. The importance of matrix porosity and fracture porosity on solute velocity
as well as the relationship governing effective dispersivity and fracture porosity is discussed for both non-reactive and
linearly sorbing solutes. By using a dimensionless effective dispersivity parameter it is shown that the relationship between
the fracture porosity and dimensionless effective dispersivity is linear for non-reactive solutes. It is also shown that this
holds true for the linearly sorbing solutes with the same proportionality constant. 相似文献
5.
Roughness scale dependence of the relationship between tracer longitudinal dispersion and Peclet number in variable‐aperture fractures 
下载免费PDF全文
下载免费PDF全文 The relationship between the longitudinal dispersion (DL) and Peclet number (Pe) is crucial for predicting and simulating tracer through the variable‐aperture fracture. In this study, the roughness of the self‐affine fracture wall was decomposed into primary roughness (relatively large‐scale waviness) and secondary roughness (relatively small‐scale waviness) by a multiscaled wavelet analysis technique. Based on the complete dispersion mechanism (diffusion, macrodispersion, and Taylor dispersion) in the variable‐aperture fracture, three relationships (second‐order, power‐law, and linear relationships) between the DL and Pe were investigated at large and small scales, respectively. Our results showed that the primary roughness mostly controlled the Taylor dispersion mechanism, whereas the secondary roughness was a dominant factor for the macrodispersion mechanism. Increasing the Hurst exponent and removing the secondary roughness led to the decreasing range of Pe where macrodispersion mechanism dominated the solute transport. It was found that estimating the DL from the power‐law relationship based on Taylor dispersion theory resulted in considerable errors, even in the range of Pe where the Taylor dispersion mechanism dominated. The exponent of the power‐law relationship increased as the secondary roughness was removed. Analysing the linear relationship between the DL and Pe revealed that the longitudinal dispersivity αL increased linearly. However, this linear increase became weak as the Taylor dispersion mechanism dominated. In the range of Pe where the macrodispersion mechanism dominated, increasing the Hurst exponent caused the increase of αL and the secondary roughness played a significant role in enhancing the αL. As the Taylor dispersion mechanism dominated, the αL was insensitive to the influence of multiscale roughness in variable‐aperture fractures. 相似文献
6.
The solutions of advection–dispersion equation in single fractures were carefully reviewed, and their relationships were addressed. The classic solution, which represents the resident or flux concentration within the semi‐infinite fractures under constant concentration or flux boundary conditions, respectively, describes the effluent concentration for a finite fracture. In addition, it also predicts the cumulative distribution of solute particle residence time passing through a single fracture under pulse injection condition, based on which a particle tracking approach was developed to simulate the local advection–dispersion in single fractures. We applied the proposed method to investigate the influence of local dispersion in single fractures on the macrodispersion in different fracture systems with relatively high fracture density. The results show that the effects of local dispersion on macrodispersion are dependent on the heterogeneity of fracture system, but generally the local dispersion plays limited roles on marodispersion at least in dense fracture network. This trend was in agreement with the macrodispersion in heterogeneous porous media. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
7.
Recognizing the spatial heterogeneity of hydraulic parameters, many researchers have studied the solute transport by both groundwater and channel flow in a stochastic framework. One of the methodologies used to up‐scale the stochastic solute transport equation, from a point‐location scale to a grid scale, is the cumulant expansion method combined with the calculus for the time‐ordered exponential and the calculus for the Lie operator. When the point‐location scale transport equation is scaled up to the grid scale, using the cumulant expansion method, a new dispersion coefficient emerges in the dispersive term of the solute transport equation in addition to the molecular dispersion coefficient. This velocity driven dispersion is called ‘macrodispersion’. The macrodispersion coefficient is the integral function of the time‐ordered covariance of the random velocity field. The integral is calculated over a Lagrangian trajectory of the flow. The Lagrangian trajectory depends on the following: (i) the spatial origin of the particle; (ii) the time when the macrodispersion is calculated; and (iii) the mean velocity field along the trajectory itself. The Lagrangian trajectory is a recursive function of time because the location of the particle along the trajectory at a particular time depends on the location of the particle at the previous time. This recursive functional form of the Lagrangian trajectory makes the calculation of the macrodispersion coefficient difficult. Especially for the unsteady, spatially non‐stationary, non‐uniform flow field, the macrodispersion coefficient is a highly complex expression and, so far, calculated using numerical methods in the discrete domains. Here, an analytical method was introduced to calculate the macrodispersion coefficient in the discrete domain for the unsteady and steady, spatially non‐stationary flow cases accurately and efficiently. This study can fill the gap between the theory of the ensemble averaged solute transport model and its numerical implementations. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
1