共查询到20条相似文献，搜索用时 31 毫秒
1.
含洞穴地层随钻电阻率测井响应自适应hp有限元数值模拟(英文) 总被引：1，自引：1，他引：0
含洞穴的碳酸盐岩地层具有强烈的非均质性及储集空间预测难度大的特点,利用随钻电阻率测井方法对井眼环境含洞穴的储层进行准确识别和划分,是当前研究的一个焦点问题。本文使用一种新型的高效和高精度自适应有限元方法(hp—FEM)模拟和分析了含洞穴地层随钻电阻率测井仪器响应。本文所提的hp-FEM与传统h-FEM相比,其结果具有网格自适应的特点,并且计算能够以指数速率收敛于较高的精度。数值实例使用自适应有限元方法研究地层中洞穴的大小、洞穴距离井眼的远近和仪器发射频率改变对测井响应的影响,并提供了识别含洞穴地层的方法。研究结果可以为实际测井中遇到的各种地层洞穴的准确识别和定量评价提供理论依据。 相似文献
2.
Numerical simulation of resistivity logging-while-drilling (LWD) tool response provides guidance for designing novel logging instruments and interpreting real-time logging data. In this paper, based on self-adaptive hp-finite element method (hp-FEM) algorithm, we analyze LWD tool response against model parameters and briefly illustrate geosteering capabilities of directional resistivity LWD. Numerical simulation results indicate that the change of source spacing is of obvious influence on the investigation depth and detecting precision of resistivity LWD tool; the change of frequency can improve the resolution of low-resistivity formation and high-resistivity formation. The simulation results also indicate that the self-adaptive hp-FEM algorithm has good convergence speed and calculation accuracy to guide the geologic steering drilling and it is suitable to simulate the response of resistivity LWD tools. 相似文献
3.
Many groundwater flow and transport problems, especially those with sharp fronts, narrow transition zones, layers and fingers, require extensive computational resources. In this paper, we present a novel multi-resolution adaptive Fup approach to solve the above mentioned problems. Our numerical procedure is the Adaptive Fup Collocation Method (AFCM), based on Fup basis functions and designed through a method of lines (MOL). Fup basis functions are localized and infinitely differentiable functions with compact support and are related to more standard choices such as splines or wavelets. This method enables the adaptive multi-resolution approach to solve problems with different spatial and temporal scales with a desired level of accuracy using the entire family of Fup basis functions. In addition, the utilized collocation algorithm enables the mesh free approach with consistent velocity approximation and flux continuity due to properties of the Fup basis functions. The introduced numerical procedure was tested and verified by a few characteristic groundwater flow and transport problems, the Buckley–Leverett multiphase flow problem, the 1-D vertical density driven problem and the standard 2-D seawater intrusion benchmark–Henry problem. The results demonstrate that the method is robust and efficient particularly when describing sharp fronts and narrow transition zones changing in space and time. 相似文献
4.
We have developed a new method to analyze the power law based non-Darcian flow toward a well in a confined aquifer with and without wellbore storage. This method is based on a combination of the linearization approximation of the non-Darcian flow equation and the Laplace transform. Analytical solutions of steady-state and late time drawdowns are obtained. Semi-analytical solutions of the drawdowns at any distance and time are computed by using the Stehfest numerical inverse Laplace transform. The results of this study agree perfectly with previous Theis solution for an infinitesimal well and with the Papadopulos and Cooper’s solution for a finite-diameter well under the special case of Darcian flow. The Boltzmann transform, which is commonly employed for solving non-Darcian flow problems before, is problematic for studying radial non-Darcian flow. Comparison of drawdowns obtained by our proposed method and the Boltzmann transform method suggests that the Boltzmann transform method differs from the linearization method at early and moderate times, and it yields similar results as the linearization method at late times. If the power index n and the quasi hydraulic conductivity k get larger, drawdowns at late times will become less, regardless of the wellbore storage. When n is larger, flow approaches steady state earlier. The drawdown at steady state is approximately proportional to r1−n, where r is the radial distance from the pumping well. The late time drawdown is a superposition of the steady-state solution and a negative time-dependent term that is proportional to t(1−n)/(3−n), where t is the time. 相似文献
5.
《Advances in water resources》2005,28(7):729-744
In this paper, we discuss the local discontinuous Galerkin (LDG) method applied to elliptic flow problems and give details on its implementation, focusing specifically on the case of piecewise linear approximating functions. The LDG method is one a family of discontinuous Galerkin (DG) methods proposed for diffusion models. These DG methods allow for very general hp finite element meshes, and produce locally conservative fluxes which can be used in coupling flow with transport. The drawback to DG methods, when compared to their continuous counterparts, is the number of degrees of freedom required to compute the solution. This motivates a coupled approach, discussed herein, where the solution is allowed to be continuous or discontinuous on a node-by-node basis. This coupled approximation is locally conservative in regions where the numerical solution is discontinuous. Numerical results for fully discontinuous, continuous and coupled discontinuous/continuous solutions are given, where we compare solution accuracy, matrix condition numbers and mass balance errors for the various approaches. 相似文献
6.
Efficient, robust simulation of groundwater flow in the unsaturated zone remains computationally expensive, especially for problems characterized by sharp fronts in both space and time. Standard approaches that employ uniform spatial and temporal discretizations for the numerical solution of these problems lead to inefficient and expensive simulations. In this work, we solve Richards’ equation using adaptive methods in both space and time. Spatial adaption is based upon a coarse grid solve and a gradient error indicator using a fixed-order approximation. Temporal adaption is accomplished using variable order, variable step size approximations based upon the backward difference formulas up to fifth order. Since the advantages of similar adaptive methods in time are now established, we evaluate our method by comparison with a uniform spatial discretization that is adaptive in time for four different one-dimensional test problems. The numerical results demonstrate that the proposed method provides a robust and efficient alternative to standard approaches for simulating variably saturated flow in one spatial dimension. 相似文献
7.
Richards Bay Harbour is South Africa’s premier bulk cargo port. It was constructed in the Mhlathuze estuary in 1976 and over the past 34 years has become South Africa’s most modern and largest cargo handling port. Although no official monitoring programme is in progress various studies by different groups have provided relevant data with respect to changing metal levels in brown mussel tissue (Perna perna) over the last 34 years. Eleven elements were analysed in brown mussels from the main channel in Richards Bay Harbour using ICP-MS. The results indicate that the metal concentrations in the mussel tissue remained relatively constant between 1974 and 2005. The mean metal concentrations increased significantly in 2005 possibly due to the construction of the new coal terminal and associated dredging activities. Mean metal concentrations in the 2008 sampling event were also elevated due to increased run off during an above average rainy season. 相似文献
8.
《Advances in water resources》2005,28(8):793-806
In this paper, the numerical errors associated with the finite difference solutions of two-dimensional advection–dispersion equation with linear sorption are obtained from a Taylor analysis and are removed from numerical solution. The error expressions are based on a general form of the corresponding difference equation. The variation of these numerical truncation errors is presented as a function of Peclet and Courant numbers in X and Y direction, a Sink/Source dimensionless number and new form of Peclet and Courant numbers in X–Y plane. It is shown that the Crank–Nicolson method is the most accurate scheme based on the truncation error analysis. The effects of these truncation errors on the numerical solution of a two-dimensional advection–dispersion equation with a first-order reaction or degradation are demonstrated by comparison with an analytical solution for predicting contaminant plume distribution in uniform flow field. Considering computational efficiency, an alternating direction implicit method is used for the numerical solution of governing equation. The results show that removing these errors improves numerical result and reduces differences between numerical and analytical solution. 相似文献
9.
10.
After 50 years of Prabhu’s paper on the exact solution of the stochastic reservoir equation for the important class of gamma inflow distributions with an integral shape parameter, a detailed implementation of the exact solution is still lacking, despite its potential usefulness from both theoretical and practical points of view. This paper explores some properties of Prabhu’s exact solution and investigates the numerical difficulties associated with its implementation. The solution is also extended to derive the distributions of deficit, spillage, yield, and actual release from the reservoir. Explicit analytical solutions for three relatively simple cases are given in detail as examples and comparisons with approximate numerical solutions are made, which reveal some shortcomings of approximate methods. The implementation of the solution in the general case reveals some numerical problems associated with large values of the shape parameter of the inflow distribution and large ratios of reservoir size to draft, mainly due to accumulation of round-off errors. A Matlab program has been developed to calculate emptying and filling probabilities over a wide range of reservoir parameters using extended precision. Comparison of Prabhu’s solution with the numerical solution of the reservoir integral equation highlights possible problems with the numerical solution, which may produce inaccurate or even invalid results for large reservoirs, large drift, and large skewness of the inflow distribution. A comparison between gamma and lognormal distributions as models of skew revealed that as the reservoir size, drift, and skewness increase, the probability of emptying of the reservoir becomes smaller for the case of gamma inflow than in the case of lognormal flow having the same skewness coefficient. 相似文献
11.
《Advances in water resources》1998,21(4):315-324
The paper deals with numerical solutions to the Richards equation to simulate one-dimensional flow processes in the unsaturated zone of layered soil profiles. The equation is expressed in the pressure-based form and a finite-difference algorithm is developed for accurately estimating the values of the hydraulic conductivity between two neighboring nodes positioned in different soil layers, often referred to as the interlayer hydraulic conductivity. The algorithm is based upon flux conservation and continuity of pressure potential at the interface between two consecutive layers, and does not add significantly to simulation run time. The validity of the model is established for a number of test problems by comparing numerical results with the analytical solutions developed by Srivastava and Yeh29 which hold for vertical infiltration towards the water table in a two-layer soil profile. The results show a significant reduction in relative mass balance errors when using the proposed model. Some specific insights into its numerical performance are also gained by comparisons with a numerical model in which the more common geometric averaging operator acts on the interlayer conductivities. 相似文献
12.
《Advances in water resources》2003,26(4):373-394
Richards’ equation (RE) is commonly used to model flow in variably saturated porous media. However, its solution continues to be difficult for many conditions of practical interest. Among the various time discretizations applied to RE, the method of lines (MOL) has been used successfully to introduce robust, accurate, and efficient temporal approximations. At the same time, a mixed-hybrid finite element method combined with an adaptive, higher order time discretization has shown benefits over traditional, lower order temporal approximations for modeling single-phase groundwater flow in heterogeneous porous media. Here, we extend earlier work for single-phase flow and consider two mixed finite element methods that have been used previously to solve RE using lower order time discretizations with either fixed time steps or empirically based adaption. We formulate the two spatial discretizations within a MOL context for the pressure head form of RE as well as a fully mass-conservative version. We conduct several numerical experiments for both spatial discretizations with each formulation, and we compare the higher order, adaptive time discretization to a first-order approximation with formal error control and adaptive time step selection. Based on the numerical results, we evaluate the performance of the methods for robustness and efficiency. 相似文献
13.
We propose a spatially and temporally adaptive solution to Richards’ equation based upon a local discontinuous Galerkin approximation in space and a high-order, backward difference method in time. We cast our approach in terms of a general, decoupled adaption algorithm based upon operators. We define non-unique instances of all operators to result in an adaption method from within the general class of methods that is defined. We formally decouple the spatial adaption from the temporal adaption using a method of lines approach and limit the temporal truncation error so that the total error is dominated by the spatial component. We use a multiple grid approach to guide adaption and support the data structures. Spatial adaption decisions are based upon error and regularity indicators, which are economical to compute. The resultant methods are compared for two test problems. The results show that the proposed adaption methods are superior to methods that adapt only in time and that in cases in which the problem has sufficient smoothness, adapting the order of the elements in addition to the grid spacing can further improve the efficiency of this robust solution approach. 相似文献
14.
Developing robust and efficient numerical solution methods for Richards' equation (RE) continues to be a challenge for certain problems. We consider such a problem here: infiltration into unsaturated porous media initially at static conditions for uniform and non-uniform pore size media. For ponded boundary conditions, a sharp infiltration front results, which propagates through the media. We evaluate the resultant solution method for robustness and efficiency using combinations of variable transformation and adaptive time-stepping methods. Transformation methods introduce a change of variable that results in a smoother solution, which is more amenable to efficient numerical solution. We use adaptive time-stepping methods to adjust the time-step size, and in some cases the order of the solution method, to meet a constraint on nonlinear solution convergence properties or a solution error criterion. Results for three test problems showed that adaptive time-stepping methods provided robust solutions; in most cases transforming the dependent variable led to more efficient solutions than untransformed approaches, especially as the pore-size uniformity increased; and the higher-order adaptive time integration method was robust and the most efficient method evaluated. 相似文献
15.
Analytical solutions for the water flow and solute transport equations in the unsaturated zone are presented. We use the Broadbridge and White nonlinear model to solve the Richards’ equation for vertical flow under a constant infiltration rate. Then we extend the water flow solution and develop an exact parametric solution for the advection-dispersion equation. The method of characteristics is adopted to determine the location of a solute front in the unsaturated zone. The dispersion component is incorporated into the final solution using a singular perturbation method. The formulation of the analytical solutions is simple, and a complete solution is generated without resorting to computationally demanding numerical schemes. Indeed, the simple analytical solutions can be used as tools to verify the accuracy of numerical models of water flow and solute transport. Comparison with a finite-element numerical solution indicates that a good match for the predicted water content is achieved when the mesh grid is one-fourth the capillary length scale of the porous medium. However, when numerically solving the solute transport equation at this level of discretization, numerical dispersion and spatial oscillations were significant. 相似文献
16.
N. I. Berezina V. I. Dmitriev N. A. Merschikova 《Izvestiya Physics of the Solid Earth》2013,49(3):350-355
We consider the iterative numerical method for solving two-dimensional (2D) inverse problems of magnetotelluric sounding, which significantly reduces the computational burden of the inverse problem solution in the class of quasi-layered models. The idea of the method is to replace the operator of the direct 2D problem of calculating the low-frequency electromagnetic field in a quasi-layered medium by a quasi-one dimensional operator at each observation point. The method is applicable for solving the inverse problems of magnetotellurics with either the E- and H-polarized fields and in the case when the inverse problem is simultaneously solved using the impedance values for the fields with both polarizations. We describe the numerical method and present the examples of its application to the numerical solution of a number of model inverse problems of magnetotelluric sounding. 相似文献
17.
A generalized transformation approach for simulating steady-state variably-saturated subsurface flow
A numerical method is proposed to accurately and efficiently compute a direct steady-state solution of the nonlinear Richards equation. In the proposed method, the Kirchhoff integral transformation and a complementary transformation are applied to the governing equation in order to separate the nonlinear hyperbolic characteristic from the linear parabolic part. The separation allows the transformed governing equation to be applied to partially- to fully-saturated systems with arbitrary constitutive relations between primary (pressure head) and secondary variables (relative permeability). The transformed governing equation is then discretized with control volume finite difference/finite element approximations, followed by inverse transformation. The approach is compared to analytical and other numerical approaches for variably-saturated flow in 1-D and 3-D domains. The results clearly demonstrate that the approach is not only more computationally efficient but also more accurate than traditional numerical solutions. The approach is also applied to an example flow problem involving a regional-scale variably-saturated heterogeneous system, where the vadose zone is up to 1 km thick. The performance, stability, and effectiveness of the transform approach is exemplified for this complex heterogeneous example, which is typical of many problems encountered in the field. It is shown that computational performance can be enhanced by several orders of magnitude with the described integral transformation approach. 相似文献
18.
《Advances in water resources》1996,19(4):207-213
It is shown that from any solution of the linear diffusion equation, we may construct a solution of a realistic form of the Richards equation for unsaturated flow. Compared to the usual direct linearization method, our inverse approach involves a quite different sequence of transformations. This opens the possibility of exact solutions with a wider variety of continuously varying flux boundary conditions. Closed-form solutions are presented for two examples. In these, the varying water flux boundary conditions resemble (i) the passage of a peaking storm and (ii) the continuous opening of a valve preceding a steady water supply. Unlike earlier more systematic approaches to this problem, our method does not require the numerical solution of an integral equation. 相似文献
19.
20.
《Advances in water resources》1999,23(3):271-301
Primary variable switching appears as a promising numerical technique for variably saturated flows. While the standard pressure-based form of the Richards equation can suffer from poor mass balance accuracy, the mixed form with its improved conservative properties can possess convergence difficulties for dry initial conditions. On the other hand, variable switching can overcome most of the stated numerical problems. The paper deals with variable switching for finite elements in two and three dimensions. The technique is incorporated in both an adaptive error-controlled predictor–corrector one-step Newton (PCOSN) iteration strategy and a target-based full Newton (TBFN) iteration scheme. Both schemes provide different behaviors with respect to accuracy and solution effort. Additionally, a simplified upstream weighting technique is used. Compared with conventional approaches the primary variable switching technique represents a fast and robust strategy for unsaturated problems with dry initial conditions. The impact of the primary variable switching technique is studied over a wide range of mostly 2D and partly difficult-to-solve problems (infiltration, drainage, perched water table, capillary barrier), where comparable results are available. It is shown that the TBFN iteration is an effective but error-prone procedure. TBFN sacrifices temporal accuracy in favor of accelerated convergence if aggressive time step sizes are chosen. 相似文献