共查询到20条相似文献，搜索用时 15 毫秒
1.
岩体裂隙中渗流场有限元随机模拟分析 总被引：2，自引：0，他引：2
采用有限元方法数值,模拟了开度随机分布裂隙中的渗流场问题。根据给定的裂隙开度均值和标准差,随机分布生成有限元模型中的单元裂隙开度,模型中的材料参数和单元属性用ANSYS中的APDL参数化语言赋值。根据有限元随机模拟断面的流量分布和稳态渗流问题的达西定律,计算在不同裂隙开度标准差条件下的等效导水系数,研究等效导水系数与裂隙粗糙度之间的关系。数值模拟结果表明,与光滑平直的裂隙相对比,在裂隙平均开度为常数、开度随机正态分布的条件下,表面粗糙裂隙的等效导水系数略有降低;对于裂隙平均开度为0.5 mm的裂隙模型,裂隙开度标准差从0.01 mm变化到 0.2 mm,其等效导水系数变化小于5 %。传统的渗流控制方程无法精确描述地下水在粗糙裂隙中的流动特征,因为在压力水头的定义中忽略了速度项。 相似文献
2.
针对截取侧向径流的抽水井,兰太权提出一个解释单位涌水量的公式,认为“单位涌水量就是导水系数”(q=T),并以砂箱试验的结果作为证据。然而,这个等式的推导过程缺少严谨的水动力学基础。砂箱试验的边界条件等价于含有1个或2个定水头边界的矩形承压含水层,单井抽水的稳定流场可以用镜像井点叠加法求解,从理论上表明单位涌水量与导水系数成正比,而比例系数取决于含水层的形状和抽水井的半径。在特定条件下,含水层和抽水井的尺寸可以造成单位涌水量与导水系数相等的情况,兰太权发表于文献中的2组砂箱试验恰好符合这种条件,但并不能证明q=T具有普适性。新的砂箱试验避免了几何尺寸的这种特殊组合,得到了显著偏离q=T的结果。 相似文献
3.
周训 《水文地质工程地质》2017,(2):184-188
兰太权先生发表3篇文章公布其多年探讨井抽水的单位涌水量和渗透系数或导水系数的关系,认为其提出了一个新的公式和论断“单位涌水量就是导水系数”,并认为传统径向稳定井流理论是错误的,在行业内应改用他提出的新理论、新公式。本文指出其文章中存在的若干问题,说明其公式或理论的不合理和不能用于实际的情形,并阐明抽水井单位涌水量不能等同于含水层导水系数。 相似文献
4.
压剪应力作用下多裂隙岩体的力学特性—断裂损伤演化方程及试验验证 总被引：5，自引：1，他引：5
首先对多裂隙岩体的损伤加以定义，在建立裂纹扩展时产生等效次生裂纹模型的基础上，对多裂隙岩体受荷后发生断裂损伤时其演化方程进行探讨，最后通过大型岩体相似模型试验对理论计算结果分析验证，理论计算与试验结果吻合较好。 相似文献
5.
山东淄博市大武水源地裂隙岩溶水中污染物运移的数值研究 总被引：11，自引：0，他引：11
在分析研究淄博市大武水源地裂隙岩溶含水层的水力性质和污染物运移特征的基础上 ,对裂隙岩溶水的水头和污染物运移进行数值研究。目前国内外对裂隙岩溶水进行数值计算时 ,通常用等价多孔介质模型 ,但裂隙岩溶介质和多孔介质有很大不同。裂隙岩溶介质的储水和导水空间为裂隙网络 ,导水系数大 ,地下水的实际平均流速比孔隙水大得多 ,但给水度和贮水系数小。当用等价多孔介质模型进行模拟时应考虑这些特点。对于污染物运移的模拟 ,要同时求解水头方程和对流弥散方程 ,可采用MODFLOW和MT3D软件进行模拟。研究区裂隙岩溶水水头的数值计算表明 ,等效多孔介质模型水头的拟合误差能满足国标GB/T144 97- 93的要求。各时段地下水水量均衡计算的精度也满足要求。对流弥散方程的数值计算 ,由于Peclet数高达 95 .6 7,对流占绝对优势 ,可能存在数值弥散和数值振荡 ,因而采用多种方法进行了比较。对于同一问题 ,同时采用上游有限差分法 (UFDM) ,混合的欧拉拉格朗日方法 (特征线法MOC、改进特征线法MMOC和混合特征线法HMOC) ,总变异消减法(TVD)进行计算 ,并比较其结果。结果表明 ,混合特征线法 (HMOC)和总变异消减法 (TVD)比较适合于对流占优势的运移问题计算。由于渗透系数K和有效孔隙度θ对溶质运移结果的影响很大 ,? 相似文献
6.
7.
多裂隙岩体的损伤断裂模型及模型试验 总被引：5，自引：0，他引：5
对于由裂隙、岩石组成的多裂隙岩体,本文通过对裂隙岩体内的应力、应变的体积平均,提出了适合于这种岩体的等效连续模型。运用损伤力学以及断裂力学理论,定义了岩体损伤张量、有效应力张量、损伤应变等,建立了损伤演化方程,从而建立了多裂隙岩体的损伤断裂模型。对岩体裂隙的扩展方式进行了模拟试验。 相似文献
8.
Multiple-pumped-well aquifer test to determine the anisotropic properties of a karst limestone aquifer in Pasco County, Florida, USA 总被引：1，自引：0，他引：1
Louis H. Motz 《Hydrogeology Journal》2009,17(4):855-869
Groundwater-level data from an aquifer test utilizing four pumped wells conducted in the South Pasco wellfield in Pasco County, Florida, USA, were analyzed to determine the anisotropic transmissivity tensor, storativity, and leakance in the vicinity of the wellfield. A weighted least-squares procedure was used to analyze drawdowns measured at eight observation wells, and it was determined that the major axis of transmissivity extends approximately from north to south and the minor axis extends approximately from west to east with an angle of anisotropy equal to N4.54°W. The transmissivity along the major axis ${\left( {T_{{\xi \xi }} } \right)}$ is 14,019 m2 day–1, and the transmissivity along the minor axis ${\left( {T_{{\eta \eta }} } \right)}$ is 4,303 m2 day–1. The equivalent transmissivity $T_{e} = {\left( {T_{{\xi \xi }} T_{{\eta \eta }} } \right)}^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}} = 7,767{{\text{m}}^{2} } \mathord{\left/ {\vphantom {{{\text{m}}^{2} } {{\text{day}}^{{ - {\text{1}}}} }}} \right. \kern-0em} {{\text{day}}^{{ - {\text{1}}}} }$ , and the ratio of anisotropy is 3.26. The storativity of the aquifer is 7.52?×?10?4, and the leakance of the overlying confining unit is 1.37?×?10?4 day?1. The anisotropic properties determined for the South Pasco wellfield in this investigation confirm the results of previous aquifer tests conducted in the wellfield and help to quantify the NW–SE to NE–SW trends for regional fracture patterns and inferred solution-enhanced flow zones in west-central Florida. 相似文献
9.
The problem of calculating equivalent grid block permeability tensors for heterogeneous porous media is addressed. The homogenization method used involves solving Darcy's equation subject to linear boundary conditions with flux conservation in subregions of the reservoir and can be readily applied to unstructured grids. The resulting equivalent permeability tensor is stable as defined relative to G-convergence. It is proposed to use both conforming and mixed finite elements to solve the local problems and compute approximations from above and below of the equivalent permeability, respectively. Comparisons with results obtained using periodic, pressure and no-flux boundary conditions and the renormalization method are presented. A series of numerical examples demonstrates the effectiveness of the methodology for two-phase flow in heterogeneous reservoirs. 相似文献
10.
The paper describes a multivariable statistical regression technique for computing the 3D stress tensor applied to the shut-in pressure data obtained on fractures during in-situ stress measurement by hydraulic fracturing in a short borehole. This method computes the mean value of components of the 3D stress tensor and their confidence intervals at a desired confidence level. This has been applied to the data of shut-in pressure from Narwapahar mine, India and Kamaishi mine, Japan. The computed in-situ stress tensor is corroborated with that obtained by other methods. 相似文献
11.
Assessment of the role of facies heterogeneity at the fine scale by numerical transport experiments and connectivity indicators 总被引：2，自引：1，他引：1
The heterogeneity of facies at the scale of individual lithological levels controls, at a macroscopic scale, water flow and contaminant transport in porous sediments. In particular the presence of organized features such as permeable connected levels, has a significant effect on travel times and dispersion. Here, the effects of facies heterogeneity on flow and transport are studied for three blocks, whose volume is of the order of a cubic meter, dug from alluvial sediments from the Ticino valley (Italy). Using the results of numerical tracer experiments on these domains, the longitudinal dispersion coefficient is computed with an Eulerian approach based on the fit of the breakthrough curves with the analytical solution of the convective-dispersive transport equation. Moreover, the dispersion tensor is computed with a Lagrangian approach from the second order moments of particle distributions. Three types of connectivity indicators are tested: (1) connectivity function; (2) flow, transport and statistical connectivity; (3) original (intrinsic, normal and total) indicators of facies connectivity. The connectivity function provides the most complete information. Some of the transport and statistical connectivity indicators are correlated with dispersivity. The simultaneous analysis of the three indicators of facies connectivity emphasizes the fundamental geometrical features that control transport. 相似文献
12.
Zehra Yenihayat 《Mathematical Geology》1996,28(2):203-228
An iterative method for calculating the transmissivity and storage coefficient from pumping test data for a confined aquifer is presented. The method optimizes the fit between the measured and the theoretical data (computed using the Theis equation) in the leastsquare sense. Unlike the existing schemes, this method employs the Levenberg-Marquardt method and the singular value decomposition technique resulting in a stable and rapidly convergent data inversion algorithm. The inverse procedure is initialized by an automatically created starting model derived using a novel technique that operates on the timederivative of the drawdown curve. An important feature of the algorithm is that all the computations are done in logarithmic space which effectively linearizes the pmblem. The proposed method has several advantages over the conventional iterative inversion algorithms because of the linearizing parameterizations at both the forward and inverse stages of the problem. Detailed derivations of the basic equations are provided to guide the potential users as well as applications to field data to demonstrate the usefulness of the proposed algorithm. 相似文献
13.
The problem of calculating equivalent grid block permeability tensors for heterogeneous porous media is addressed. The homogenization method used involves solving Darcy's equation subject to linear boundary conditions with flux conservation in subregions of the reservoir and can be readily applied to unstructured grids. The resulting equivalent permeability tensor is stable as defined relative to G-convergence. It is proposed to use both conforming and mixed finite elements to solve the local problems and compute approximations from above and below of the equivalent permeability, respectively. Comparisons with results obtained using periodic, pressure and no-flux boundary conditions and the renormalization method are presented. A series of numerical examples demonstrates the effectiveness of the methodology for two-phase flow in heterogeneous reservoirs. 相似文献
14.
Hydrogeological framework and estimation of aquifer hydraulic parameters using geoelectrical data: a case study from West Iran 总被引：1，自引：0，他引：1
Abduulah Taheri Tizro Kostantanous S. Voudouris Mortaza Salehzade Hassan Mashayekhi 《Hydrogeology Journal》2010,18(4):917-929
Detailed local geological, geophysical, and hydrogeological investigations were carried out for the alluvial aquifer in the Kangavar basin, West Iran to delineate the architecture of different subsurface geological horizons using lithologs and generated vertical electrical sounding (VES) data. An attempt has also been made to estimate aquifer transmissivity from resistivity data. Forty VESs were recorded with the Schlumberger electrode configuration in the study area; 28 of these were selected for evaluation. The maximum current electrode spacing was 400–500 m. The data obtained were interpreted by computer iterative modeling with curve matching for calibration purposes. In order to ascertain the subsurface geological framework, the general distribution of resistivity responses of the geological formations was obtained and geoelectrical sections along a number of lines were prepared. Probable aquifer horizons from these sections were identified. The transmissivity of the unconfined aquifer was computed by determining the Dar-Zarrouk parameters (longitudinal unit conductance and transverse unit resistance) and were compared with the actual field transmissivity. The results showed a direct relation between aquifer transmissivity and modified transverse resistance. 相似文献
15.
仵彦卿 《水文地质工程地质》1997,24(3):10-14
本文把裂隙岩体看作等效连续介质，用渗透系数张量描述岩体的渗透性能；用应力张量描述岩体的力学性能。介绍岩体渗流场与应力场耦合的等效连续介质模型。 相似文献
16.
Estimation of effective aquifer hydraulic properties from an aquifer test with multi-well observations (Taiwan) 总被引：1，自引：0，他引：1
Jet-Chau Wen Cheng-Mau Wu Tian-Chyi J. Yeh Chung-Ming Tseng 《Hydrogeology Journal》2010,18(5):1143-1155
The representative anisotropic parameters for an aquifer located at the campus of the National Yunlin University of Science and Technology in Taiwan have been acquired. A constant-rate pumping test was carried out, and drawdown-time data were collected from ten observation wells. Applications of the conventional aquifer test analysis, which assumes aquifer homogeneity for each observed well hydrograph, yielded spatially varying transmissivity and storage coefficient estimates, contradicting the homogeneous assumption. A direct approach and a nonlinear-least squares minimization of distance-drawdown data were then employed to analyze anisotropy of the transmissivity of an equivalent homogeneous aquifer. Results show that the direction and values of anisotropic transmissivities vary with time and the estimates depend on the number of observation wells used. These field results are consistent with results from recent theoretical investigations which questioned the suitability of the conventional aquifer test analysis. 相似文献
17.
Andarge Yitbarek Baye Moumtaz Razack Tenalem Ayenew Engida Zemedagegnehu 《Environmental Earth Sciences》2013,69(6):1791-1802
Transmissivity (T) is a basic hydraulic parameter of an aquifer that is utilized in most groundwater flow equations to understand the flow dynamics and is generally estimated from pumping tests. However, the cost of performing a large number of aquifer tests is expensive and time consuming. The fact that specific capacity (S c) is correlated with hydraulic flow properties of aquifers simplifies parameter estimation mainly because specific capacity values are more abundant in groundwater databases than values of transmissivity and they offer another approach to estimate hydraulic parameters of aquifers. In this study, an empirical relation is derived using 214 pairs of transmissivity and specific capacity values that are obtained from pumping tests conducted on water wells penetrating the complex volcanic aquifers of Upper Awash Basin, central Ethiopia. Linear and logarithmic regression functions have been performed and it is found that the logarithmic relationship predicting transmissivity from specific capacity data has a better correlation (R = 0.97) than the linear relationship (R = 0.79). The two parameters are log-normally distributed, in which the logarithmic relation is also better statistically justified than the linear relation. Geostatistical estimations of the transmissivity were made using different inputs and methods. Measured and supplemented transmissivity data obtained from estimates using the derived empirical relation were krigged and cokrigged, spherical and exponential models were fitted to the experimental variograms. The cross-validation results showed that the best estimation is provided using the kriging procedure, the transmissivity field represented by the measured transmissivity data and the experimental variogram fitted with the exponential model. Based on the geostatistical approach, the transmissivity map of the aquifer is produced, which will be used for groundwater flow modeling of the study area that will follow this analysis. 相似文献
18.
The strength of granular materials during triaxial compression is investigated via a grain scale analysis in this paper. A 3D Discrete Element Method (DEM) program provides the triaxial strength data and helps to validate the micromechanical analysis. Some standard methods in statistics are employed first to quantitatively examine the assumptions made when deriving the stress-force-fabric (SFF) equation. After careful validation, a more concise format for the SFF equation is proposed for triaxial compressions. With this SFF equation, the strength is found to be jointly contributed by the magnitudes of the contact force anisotropy and fabric anisotropy. The influence of the initial void ratio, confining pressure and loading direction on the development of contact force anisotropy and fabric anisotropy is examined and presented. With similar techniques, the “force” term in the SFF equation is further decoupled, and an equation is obtained such that it explicitly links the contact force term with the friction coefficient between grains, a tensor defined as a statistic of the normal contact forces and a tensor defined as a statistic of the mobilisation status of contacts. Based on this equation, another equation regarding the stress ratio of granular assembly is obtained, and it clearly indicates two sources that contribute to the phenomenological friction nature of granular assembly. These two sources are caused by the contact force at the grain scale. The first is the anisotropy of the average normal contact forces, and the second is the mobilisation of contacts. 相似文献
19.
We present a new family of flux continuous, locally conservative, finite volume schemes applicable to the diagonal and full tensor pressure equations with generally discontinuous coefficients. For a uniformly constant symmetric elliptic tensor field, the full tensor discretization is second order accurate with a symmetric positive definite matrix. For a full tensor, an M-matrix with diagonal dominance can be obtained subject to a sufficient condition for ellipticity. Positive definiteness of the discrete system is illustrated. Convergence rates for discontinuous coefficients are presented and the importance of modeling the full permeability tensor pressure equation is demonstrated. 相似文献
20.