共查询到20条相似文献，搜索用时 15 毫秒
1.
The spatial distribution of hydraulic properties in the subsurface controls groundwater flow and solute transport. However, many approaches to modeling these distributions do not produce geologically realistic results and/or do not model the anisotropy of hydraulic conductivity caused by bedding structures in sedimentary deposits. We have developed a flexible object-based package for simulating hydraulic properties in the subsurface—the Hydrogeological Virtual Realities (HyVR) simulation package. This implements a hierarchical modeling framework that takes into account geological rules about stratigraphic bounding surfaces and the geometry of specific sedimentary structures to generate realistic aquifer models, including full hydraulic-conductivity tensors. The HyVR simulation package can create outputs suitable for standard groundwater modeling tools (e.g., MODFLOW), is written in Python, an open-source programming language, and is openly available at an online repository. This paper presents an overview of the underlying modeling principles and computational methods, as well as an example simulation based on the Macrodispersion Experiment site in Columbus, Mississippi. Our simulation package can currently simulate porous media that mimic geological conceptual models in fluvial depositional environments, and that include fine-scale heterogeneity in distributed hydraulic parameter fields. The simulation results allow qualitative geological conceptual models to be converted into digital subsurface models that can be used in quantitative numerical flow-and-transport simulations, with the aim of improving our understanding of the influence of geological realism on groundwater flow and solute transport. 相似文献
2.
现有海洋可控源电磁三维数值模拟方法大多基于电导率各向同性介质理论,不能模拟海底地层电导率各向异性的实际情况.本文给出了电导率各向异性三维介质中电性源海洋可控源电磁二次电场的边值问题以及相应的变分问题,采用长方体单元对研究区域剖分,将场分量定义在剖分单元的边上,利用矢量有限单元法求解变分问题,实现了电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟.这个新的正演方法可以计算电导率任意各向异性三维地电模型的海洋可控源电磁响应,基于二次场矢量有限元法直接求解电磁场,避免了传统有限元方法可能遇到的伪解问题和难于处理电场法向分量不连续的问题,提高了数值模拟计算精度.一维电导率各向异性模型电磁场数值解与解析解吻合得相当好,无论在源附近还是远离源处相对误差均不超过1%.电导率各向异性二维模型的计算结果与已有文献采用的非结构有限元模拟结果十分吻合.三维地电模型数值模拟结果显示,电导率各向异性张量电导率主轴分量和欧拉角对不同装置海洋可控源电磁响应均有着明显的影响. 相似文献
3.
4.
本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究. 相似文献
5.
Influence of bedding orientation on the anisotropy of hydraulic conductivity in a well-sorted fluvial sediment 总被引：2，自引：0，他引：2
The paper describes a permeameter test method for determination of the hydraulic conductivity(AT) along multi-directions in fluvial sediments with cross beddings.Unlike existing in-situ permeameter methods that determine hydraulic conductivity for submerged streambeds,our method was intended to measure hydraulic conductivity of exposed streambeds or fluvial sediments.The method was applied to the Wei River,Shaanxi Province,Central China for characterization of the anisotropy of K in a well-sorted fluvial sediment.The results illustrated that even in well-sorted sediments,cross-bedding and sediment fabrication(or texture) can lead to varied K values along different measurement directions.The K value was the largest along the dip direction(or the major direction) that is parallel to the orientation of cross bedding and the smallest in the direction perpendicular to the bedding(or the minor direction). The K value in a given direction between the major and minor direction often fell in the range bounded by the K values in the major and minor directions.The anisotropy ratio of K(the ratio of K value between the major and minor directions) in two trenches for this well-sorted fluvial sediment was up to 1.14 to 1.23,respectively.Our results also demonstrated that even for well-sorted sediments,the K values between two sampling points only about 10 cm apart can differ.It is clear that the K distribution strongly correlates to the bedding orientation. 相似文献
6.
On the relation between steep monoclinal flexure zones and steep hydraulic gradients 总被引：1，自引：0，他引：1
Steep hydraulic gradients are found in association with steep monoclinal flexures. However, the physics of the reduction of the hydraulic conductivity, which is responsible for the steep gradients, has seldom been studied. We present results of hydrological and mechanical modeling aiming to study the effect of such steep hydraulic gradients demonstrated in the Judea Group Aquifer system, Israel. The hydrological configuration of steep dips and anisotropy between flows parallel and perpendicular to the bedding planes was simulated using the FEFLOW code. It exhibited a situation whereby part of the flow is oblique to the bedding planes and therefore some steepening of the hydraulic gradients occurred due to actual conductivity reduction. However, this reduction is not enough to account for the steeper gradients observed. The effect of a deep-seated reverse fault under the monocline on the permeability distribution within the structure was examined by numerical mechanical simulations. It exhibited a compressional stress distribution in the steep part of the monocline, which, due to shortening and closure of joints and voids, is presumably responsible for a significant pressure-induced permeability reduction. This process by itself in a layered structure, including interlayering of thin marl layers, could be responsible for the steep hydraulic gradients in the steep part of the monocline. 相似文献
7.
Electrical anisotropy occurs when electric current flow varies with azimuth. In porous media, this may correspond to anisotropy in the hydraulic conductivity resulting from sedimentary fabric, fractures, or dissolution. In this study, a 28‐electrode resistivity imaging system was used to investigate electrical anisotropy at 13 sites in the Biscayne Aquifer of SE Florida using the rotated square array method. The measured coefficient of electrical anisotropy generally ranged from 1.01 to 1.12 with values as high as 1.36 found at one site. The observed electrical anisotropy was used to estimate hydraulic anisotropy (ratio of maximum to minimum hydraulic conductivity) which ranged from 1.18 to 2.83. The largest values generally were located on the Atlantic Coastal Ridge while the lowest values were in low elevation areas on the margin of the Everglades to the west. The higher values of anisotropy found on the ridge may be due to increased dissolution rates of the oolitic facies of the Miami formation limestone compared with the bryozoan facies to the west. The predominate trend of minimum resistivity and maximum hydraulic conductivity was E‐W/SE‐NW beneath the ridge and E‐W/SW‐NE farther west. The anisotropy directions are similar to the predevelopment groundwater flow direction as indicated in published studies. This suggests that the observed anisotropy is related to the paleo‐groundwater flow in the Biscayne Aquifer. 相似文献
8.
This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow‐through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream‐aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow‐through streams. 相似文献
9.
Hydraulic Gradient Comparison Method to Estimate Aquifer Hydraulic Parameters Under Steady-State Conditions 总被引：1，自引：0，他引：1
The hydraulic gradient comparison method is an inverse method for estimation of aquifer hydraulic conductivity (or trans-missivity) and boundary conductance for a ground water flow model under steady-state conditions. This method, following formal optimization techniques, defines its objective function to minimize differences between interpreted (observed) and simulated hydraulic gradients, which results in minimization of differences between observed and simulated hydraulic heads. The key features of this method are that (1) the derived optimality conditions have an explicit form with a clear hydrology concept that is con-sistent with Darcy's law, and (2) the derived optimality conditions are spatially independent as they are a function of only local hydraulic conductivity and local hydraulic gradient. This second feature allows a multidimensional optimization problem to be solved by many one-dimensional optimization procedures simultaneously, which results in a substantial reduction in computation time. The results of the numerical performance testing on a heterogeneous hypothetical case confirm that minimizing gradient residuals in the entire model domain leads to minimizing head residuals. Application of the method in real-world projects requires rigorous conceptual model development, use of a global calibration target, and an iterative calibration proess. The conceptual model development includes interpretation of a potentiometric surface and estimation of other hydrologic parameters. This method has been applied to a wide range of real-world modeling projects, including the Rocky Mountain Arsenal and Rocky Flats sites in Colorado, which demonstrates that the method is efficient and practical. 相似文献
10.
The airborne electromagnetic (AEM) method has a high sampling rate and survey flexibility. However, traditional numerical modeling approaches must use high-resolution physical grids to guarantee modeling accuracy, especially for complex geological structures such as anisotropic earth. This can lead to huge computational costs. To solve this problem, we propose a spectral-element (SE) method for 3D AEM anisotropic modeling, which combines the advantages of spectral and finite-element methods. Thus, the SE method has accuracy as high as that of the spectral method and the ability to model complex geology inherited from the finite-element method. The SE method can improve the modeling accuracy within discrete grids and reduce the dependence of modeling results on the grids. This helps achieve high-accuracy anisotropic AEM modeling. We first introduced a rotating tensor of anisotropic conductivity to Maxwell’s equations and described the electrical field via SE basis functions based on GLL interpolation polynomials. We used the Galerkin weighted residual method to establish the linear equation system for the SE method, and we took a vertical magnetic dipole as the transmission source for our AEM modeling. We then applied fourth-order SE calculations with coarse physical grids to check the accuracy of our modeling results against a 1D semi-analytical solution for an anisotropic half-space model and verified the high accuracy of the SE. Moreover, we conducted AEM modeling for different anisotropic 3D abnormal bodies using two physical grid scales and three orders of SE to obtain the convergence conditions for different anisotropic abnormal bodies. Finally, we studied the identification of anisotropy for single anisotropic abnormal bodies, anisotropic surrounding rock, and single anisotropic abnormal body embedded in an anisotropic surrounding rock. This approach will play a key role in the inversion and interpretation of AEM data collected in regions with anisotropic geology. 相似文献
11.
为了实现交错沉积等复杂环境中的电磁场数值模拟,本文在常规横向同性模型的基础上引入了电导率主轴坐标系相对地层坐标系的层理方位角和倾角,建立了交错地层模型.并利用传播矩阵法建立了一维层状交错地层模型中的多分量感应测井仪器响应的正演模拟算法.首先将频率-波数域中的电磁场分解为上行和下行模式波,给出了任意朝向的磁偶极子在无限大地层中模式波的解析解.进一步通过引入地层界面上的透射、局部反射以及广义反射系数矩阵,推导了一维层状地层中的模式波表达式.在此基础上,利用二维Gauss-Legendre积分实现了Fourier逆变换,得到了可用于多分量感应测井模拟的频率-空间域磁场并矢格林函数.最后,通过多个数值模拟结果考察了井眼倾角、层理方位角和倾角变化对多分量感应测井响应的影响. 相似文献
12.
在一些地层层理发育的地区,地下介质存在显著的电各向异性,此时基于各向同性模型解释含各向异性效应的可控源音频大地电磁（CSAMT）测深观测数据会导致错误的结果.本文通过引入3×3的对称正定张量表征电导率各向异性,采用非结构四面体网格和矢量有限元方法离散电场满足的矢量Helmholtz方程,并将电磁场源等效为系列电偶极子,实现任意各向异性介质中CSAMT高效数值模拟.本文首先通过层状各向异性模型检验三维有限元算法的精度和有效性,进一步建立三维地电模型研究异常体各向异性和围岩各向异性对CSAMT响应的影响,最后使用视电阻率极性图来识别各向异性电导率主轴方向.数值模拟结果表明,各向异性电导率对CSAMT视电阻率幅值及分布规律都有很大影响,视电阻率极性图能够很好地识别各向异性主轴方向. 相似文献
13.
The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution 
下载免费PDF全文
下载免费PDF全文 The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
14.
WANG Xiaoqiong XU Jianguo ZHAO Chenxu LIU Tongyuan GE Hongkui SHEN Yinghao WU Shan YU Jiayao HUANG Rongyan 《中国地震研究》2020,34(2):167-186
The premise of hydraulic fracturing is to have an accurate and detailed understanding of the rock mechanical properties and fracture propagation law of shale reservoirs. In this paper, a comprehensive evaluation of the mechanical properties of the shale oil reservoir in the south of Songliao Basin is carried out. Based on the experiments and the in-situ stress analysis, the fracture propagation law of three types of shale reservoirs is obtained, and the suggestions for fracturing are put forward. The results have shown that the fracture propagation of pure shale and low mature reservoir is easy to open along the bedding plane under compression loading, which is greatly influenced by the bedding. Sand-bearing shale is slightly better, the fractures of which are not easy to open along the bedding plane. The mechanical experimental results show that all the samples have the characteristics of low compressive strength, low Young''s modulus and strong anisotropy, indicating that the shale oil reservoir is certain plastic, which is related to its high clay mineral content and controlled by the bedding development. Compared with pure shale and low mature shale, the sand-bearing shale has less clay content and less developed bedding, which maybe the main reason for its slightly better brittleness. Overall, the expansion of hydraulic fracture is controlled by in-situ stress and bedding. Because of the development of bedding, it is easy to form horizontal fractures. Thus it is not suitable for horizontal well fracturing. Because of the high content of clay minerals, the applicability of conventional slick hydraulic fracturing fluid is poor. It is suggested to use vertical well or directional well to carry out volume fracturing. In this way, the effect of bedding can be effectively used to open and connect the bedding and form a larger fracture network. 相似文献
15.
目前,瞬变电磁法(TEM)数据基本都是基于各向同性模型进行反演解释,这对于存在明显电性各向异性的勘探区域会产生较大的反演解释误差.为分析电各向异性对回线源瞬变电磁信号的影响方式与程度,本文通过求解离散化的全张量电导率时间域Helmholtz方程,实现了基于有限体积法的TEM任意各向异性的三维正演算法.该算法采用基于交错网格的拟态有限体积法(MFV)对时域Maxwell方程组进行空间域离散,并利用后退欧拉算法(Backward Euler Method)进行时间域离散.为提高时域电磁场的求解精度与效率,该算法将时间分段等步长算法与方程直接求解法相结合.通过对一维各向异性模型以及三维复杂各向同性模型进行测试,验证了本算法对于回线源瞬变电磁响应计算的正确性及有效性.最后,通过对几类典型电各向异性介质中大回线源瞬变电磁信号响应的分析,总结了不同电各向异性类型对TEM电磁信号的影响模式,结果表明,主轴各向异性情况下TEM信号主要受水平方向电导率的影响,倾斜各向异性对TEM信号的影响程度远大于水平各向异性,而通过水平各向异性信号能较清晰判断出各向异性主轴方向. 相似文献
16.
Vertically Integrated Hydraulic Conductivity: A New Parameter for Groundwater-Surface Water Analysis
Numerical modeling of groundwater-surface water interactions provides vital information necessary for determining the extent of nutrient transport, quantifying water budgets, and delineating zones of ecological support. The hydrologic data that drive these models are often collected at disparate scales and subsequently incorporated into numerical models through upscaling techniques such as piecewise constancy or geostatistical methods. However, these techniques either use basic interpolation methods, which often simplifies the system of interest, or utilize complex statistical methods that are computationally expensive, time consuming, and generate complex subsurface configurations. We propose a bulk parameter termed “vertically integrated hydraulic conductivity” (KV), and defined as the depth-integrated resistance to fluid flow sensed at the groundwater-surface water interface, as an alternative to hydraulic conductivity when investigating vertical fluxes across the groundwater-surface water interface. This bulk parameter replaces complex subsurface configurations in situations dominated by vertical fluxes and where heterogeneity is not of primary importance. To demonstrate the utility of KV, we extracted synthetic temperature time series data from a forward numerical model under a variety of scenarios and used those data to quantify vertical fluxes using the amplitude ratio method. These quantified vertical fluxes and the applied hydraulic head gradient were subsequently input into Darcy's Law and used to quantify KV. This KV was then directly compared to the equivalent hydraulic conductivity (KT) assuming an infinitely extending layer. Vertically integrated hydraulic conductivity allows for more accurate and robust flow modeling across the groundwater-surface water interface in instances where complex heterogeneities are not of primary concern. 相似文献
17.
当地表存在三维非均匀电导率分布时,区域大地电磁响应发生畸变. 以往对这种畸变研究多假设近地表为三维,区域构造为一维或二维. 对于更一般的三维/三维构造,为了分析并消除这种畸变影响,真实反映地下三维区域构造信息,本文实现了三维大地电磁相位张量积分方程数值算法,并研究在不同地质模型下相位张量响应. 结果表明,相位张量不仅可以反映一般三维构造信息,亦可有效反映复杂近地表构造下三维区域构造信息,而无须假设区域构造为一维或二维,证明相位张量具有较强抗近地表局部非均匀构造干扰能力,能够保持更为一般的三维区域构造信息. 为了加快正演计算,同时保持一定精度,算法采用了积分方程多网格法. 相似文献
18.
实际地层可能出现 2种电阻率的宏观各向异性 :微观的统计平均效应和裂缝走向等因素的构造效应。三分量感应测井方法可能是确定地层电阻率各向异性的最好方法 ,但是 ,该仪器目前在中国还没有应用。已有的数据一般是梯度和双感应曲线 ,梯度和双感应曲线单独使用都难以揭示各向异性 ,但是 ,将二者联合使用则有可能。在反演模型中考虑了电阻率各向异性的存在 ,给出了梯度和双感应曲线联合求地层电阻率各向异性的方法。通过直井水平层的二维人工模型表明 ,用感应方法可确定地层的水平电阻率和厚度。而梯度方法的视电阻率可以近似看作水平电阻率和纵向电阻率的几何平均值 ,梯度方法的视厚度是各向异性参数和真实厚度的乘积。因而不能用梯度方法单独确定各向异性参数和真实厚度 2个参数中的任何 1个。然而 ,联合梯度和双感应的方法则可能确定 3个参数 :各向异性参数 ,水平电阻率和地层厚度。人工数据的例子表明了这种联合的可行性 相似文献
19.
Fractures in porous media have been documented extensively. However, they are often omitted from groundwater flow and mass transport models due to a lack of data on fracture hydraulic properties and the computational burden of simulating fractures explicitly in large model domains. We present a MATLAB toolbox, FracKfinder, that automates HydroGeoSphere (HGS), a variably saturated, control volume finite-element model, to simulate an ensemble of discrete fracture network (DFN) flow experiments on a single cubic model mesh containing a stochastically generated fracture network. Because DFN simulations in HGS can simulate flow in both a porous media and a fracture domain, this toolbox computes tensors for both the matrix and fractures of a porous medium. Each model in the ensemble represents a different orientation of the hydraulic gradient, thus minimizing the likelihood that a single hydraulic gradient orientation will dominate the tensor computation. Linear regression on matrices containing the computed three-dimensional hydraulic conductivity (K) values from each rotation of the hydraulic gradient is used to compute the K tensors. This approach shows that the hydraulic behavior of fracture networks can be simulated where fracture hydraulic data are limited. Simulation of a bromide tracer experiment using K tensors computed with FracKfinder in HGS demonstrates good agreement with a previous large-column, laboratory study. The toolbox provides a potential pathway to upscale groundwater flow and mass transport processes in fractured media to larger scales. 相似文献
20.
实际地层可能出现2种电阻率的宏观各向异性:微观的统计平均效应和裂缝走向等因素的构造效应.三分量感应测井方法可能是确定地层电阻率各向异性的最好方法,但是,该仪器目前在中国还没有应用.已有的数据一般是梯度和双感应曲线,梯度和双感应曲线单独使用都难以揭示各向异性,但是,将二者联合使用则有可能.在反演模型中考虑了电阻率各向异性的存在,给出了梯度和双感应曲线联合求地层电阻率各向异性的方法.通过直井水平层的二维人工模型表明,用感应方法可确定地层的水平电阻率和厚度.而梯度方法的视电阻率可以近似看作水平电阻率和纵向电阻率的几何平均值,梯度方法的视厚度是各向异性参数和真实厚度的乘积.因而不能用梯度方法单独确定各向异性参数和真实厚度2个参数中的任何1个.然而,联合梯度和双感应的方法则可能确定3个参数:各向异性参数,水平电阻率和地层厚度.人工数据的例子表明了这种联合的可行性. 相似文献