共查询到20条相似文献，搜索用时 218 毫秒
1.
A local artificial boundary for transient seepage problems with unsteady boundary conditions in unbounded domains 
下载免费PDF全文
下载免费PDF全文 Numerical analysis of transient seepage in unbounded domains with unsteady boundary conditions requires a more sophisticated artificial boundary approach to deal with the infinite character of the domain. To that end, a local artificial boundary is established by simplifying a global artificial boundary. The global artificial boundary conditions (ABCs) at the truncated boundary are derived from analytical solutions for one‐dimensional axisymmetric diffusion problems. By applying Laplace transforms and introducing some specially defined auxiliary variables, the global ABCs are simplified to local ABCs to significantly enhance the computational efficiency. The proposed local ABCs are implemented in a finite element computer program so that the solutions to various seepage problems can be calculated. The proposed approach is first verified by the computation of a one‐dimensional radial flow problem and then tentatively applied to more general two‐dimensional cylindrical problems and planar problems. The solutions obtained using the local ABCs are compared with those obtained using a large element mesh and using a previously proposed local boundary. This comparison demonstrates the satisfactory performance and obvious superiority of the newly established boundary to the other local boundary. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
2.
3.
局域地应力场获取的插值平衡方法 总被引：4，自引：1，他引：3
针对由整体地应场插值计算局域地应力场的问题计算量大、受地质结构的影响大、不过分追求高精度等特点,改进了反距离加权插值法,考虑了单元体积的影响,提高了该方法在此问题中的适用性。基于此,提出了平面问题、三维地面和地下工程局域地应力场的求取策略,建议了边界全平动约束和法向约束两种边界条件,并验证了其合理性和适用场合。最后,在大型水电工程深切河谷坝肩边坡分析中进行了应用分析。结果表明：该插值方法简单实用,能够满足工程分析的精度要求;局域地应力场的求取策略计算效率高,计算精度能够满足工程分析的要求。 相似文献
4.
本文利用面积广、勘探深度大的重磁数据通过增强型局部相位边界识别方法和优化欧拉反褶积法获得东海陆架盆地的构造特征,为油气勘探提供基础资料。增强型局部相位边界识别法是重磁数据一阶和二阶水平导数的比值函数,能有效地均衡不同深度地质体的边界响应,从而同时显示不同深度构造的边界,提高对较深层构造的分辨率;优化欧拉反褶积法有效地降低了背景异常的干扰,从而可获得更加准确的构造反演结果。理论模型研究表明,增强型局部相位边界识别法可更加准确和清晰地划分边界信息,优化欧拉反褶积法可获得更加准确和收敛的反演结果。综合利用增强型局部相位边界识别和优化欧拉反褶积法对东海陆架断裂构造分布进行划分,获得了分析区域构造背景和评估油气潜力的基础资料。 相似文献
5.
6.
We describe a fast geomechanically-based paleostress inversion technique that uses observed fracture data to constrain stress through multiple simulations. The method assumes that the local stress field around individual fractures is heterogeneous and derives the far field tectonic stress, that we also call the far field boundary conditions. We show how such far field tectonic stress can be recovered through a mechanical stress inversion technique using local observations of natural fractures (i.e. mechanical type, orientation and location). We test the paleostress inversion against outcrop analogues of fractured carbonates from both Nash Point, U.K., where there are well exposed faults and joints and the Matelles, France, where there are well exposed faults, veins and stylolites. We demonstrate through these case studies how the method can be efficiently applied to natural examples and we highlight its advantages and limitations. We discuss how such method could be applied to subsurface problems and how it can provide complementary constraints to drive discrete fracture models for better fractured reservoir characterization and modelling. 相似文献
7.
We present a variational multiscale mixed finite element method for the solution of Darcy flow in porous media, in which both
the permeability field and the source term display a multiscale character. The formulation is based on a multiscale split
of the solution into coarse and subgrid scales. This decomposition is invoked in a variational setting that leads to a rigorous
definition of a (global) coarse problem and a set of (local) subgrid problems. One of the key issues for the success of the
method is the proper definition of the boundary conditions for the localization of the subgrid problems. We identify a weak
compatibility condition that allows for subgrid communication across element interfaces, a feature that turns out to be essential
for obtaining high-quality solutions. We also remove the singularities due to concentrated sources from the coarse-scale problem
by introducing additional multiscale basis functions, based on a decomposition of fine-scale source terms into coarse and
deviatoric components. The method is locally conservative and employs a low-order approximation of pressure and velocity at
both scales. We illustrate the performance of the method on several synthetic cases and conclude that the method is able to
capture the global and local flow patterns accurately. 相似文献
8.
Knut-Andreas Lie Jostein R. Natvig Stein Krogstad Yahan Yang Xiao-Hui Wu 《Computational Geosciences》2014,18(3-4):357-372
A Dirichlet–Neumann representation method was recently proposed for upscaling and simulating flow in reservoirs. The DNR method expresses coarse fluxes as linear functions of multiple pressure values along the boundary and at the center of each coarse block. The number of flux and pressure values at the boundary can be adjusted to improve the accuracy of simulation results and, in particular, to resolve important fine-scale details. Improvement over existing approaches is substantial especially for reservoirs that contain high-permeability streaks or channels. As an alternative, the multiscale mixed finite-element (MsMFE) method was designed to obtain fine-scale fluxes at the cost of solving a coarsened problem, but can also be used as upscaling methods that are flexible with respect to geometry and topology of the coarsened grid. Both methods can be expressed in mixed-hybrid form, with local stiffness matrices obtained as “inner products” of numerically computed basis functions with fine-scale sub-resolution. These basis functions are determined by solving local flow problems with piecewise linear Dirichlet boundary conditions for the DNR method and piecewise constant Neumann conditions for MsMFE. Adding discrete pressure points in the DNR method corresponds to subdividing faces in the coarse grid and hence increasing the number of basis functions in the MsMFE method. The methods show similar accuracy for 2D Cartesian cases, but the MsMFE method is more straightforward to formulate in 3D and implement for general grids. 相似文献
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 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. 相似文献
11.
Large-scale flow models constructed using standard coarsening procedures may not accurately resolve detailed near-well effects.
Such effects are often important to capture, however, as the interaction of the well with the formation can have a dominant
impact on process performance. In this work, a near-well upscaling procedure, which provides three-phase well-block properties,
is developed and tested. The overall approach represents an extension of a recently developed oil–gas upscaling procedure
and entails the use of local well computations (over a region referred to as the local well model (LWM)) along with a gradient-based
optimization procedure to minimize the mismatch between fine and coarse-scale well rates, for oil, gas, and water, over the
LWM. The gradients required for the minimization are computed efficiently through solution of adjoint equations. The LWM boundary
conditions are determined using an iterative local-global procedure. With this approach, pressures and saturations computed
during a global coarse-scale simulation are interpolated onto LWM boundaries and then used as boundary conditions for the
fine-scale LWM computations. In addition to extending the overall approach to the three-phase case, this work also introduces
new treatments that provide improved accuracy in cases with significant flux from the gas cap into the well block. The near-well
multiphase upscaling method is applied to heterogeneous reservoir models, with production from vertical and horizontal wells.
Simulation results illustrate that the method is able to accurately capture key near-well effects and to provide predictions
for component production rates that are in close agreement with reference fine-scale results. The level of accuracy of the
procedure is shown to be significantly higher than that of a standard approach which uses only upscaled single-phase flow
parameters. 相似文献
12.
Greco提出的临界滑裂面搜寻的Monte Carlo搜索技术具有原理简单、适应性强、搜索性能好及容易编程实现等优点,目前该类方法已在边坡稳定的优化计算中有许多成功的应用。为能够在多块体上限法求解地基极限承载力的优化计算中应用Monte Carlo搜索技术,基于上限法相容速度场的要求及地基承载力问题的特点,对Monte Carlo搜索技术实现中的目标函数、几何约束条件、初始破坏面的产生以及收敛准则等重新进行了设置。由于Monte Carlo搜索技术随机性的特点,某一次的搜索优化往往存在陷入局部极值的危险,为解决这一不足,采取了随机设置多次初始破坏面分别进行搜索优化的办法,计算表明,优化效果很好。通过对实际问题的计算以及对比发现,Monte Carlo搜索技术在此处多块体上限法求解地基承载力问题中的优化应用是相当成功的。此外,由于此处的优化是以破坏面上的节点为对象的,因此,不但可以方便地考虑三角形块体的优化,而且也可以方便地考虑四边形块体的优化问题,其适用性更强。 相似文献
13.
紧支径向基函数能使支配方程中的刚度矩阵具有稀疏性,很适合应用于无网格方法中,其缺点是在插值计算时精度不高.点插值方法的插值函数具有Delta函数性质,可以很方便的施加本质边界条件,但在计算插值函数时矩阵易出现奇异.为了提高计算精度并避免点插值法的局限性,首先对紧支径向基函数进行完备性修正,然后用完备性修正的紧支径向基函数代替多项式来形成插值函数,建立了紧支径向基函数点插值方法.由于该方法中的形函数满足Delta函数性质,因此本质边界条件可以像传统的有限元方法一样很容易施加.然后将该方法用于二维弹性静力问题的求解,导出了其相应的离散方程.最后将该方法应用于一个悬臂梁的分析中,初步验证了该方法的有效性与合理性. 相似文献
14.
W. A. Mulder 《Computational Geosciences》1997,1(1):85-108
Simulation of wave propagation for seismic purposes is usually restricted to a small portion of the earth. Artificial boundary
conditions are required where the subsurface model is truncated. Absorbing boundaries should ensure that waves hitting the
artificial boundaries are not reflected. The vast amount of literature on the subject suggests that “good” conditions have
not been found, and only “reasonable” solutions exist. A cursory overview of existing and a few new ideas is presented that
may guide the construction of suitable boundary conditions. Because the intended application of the boundary conditions was
a high-order finite-difference code that runs on a parallel computer, we have restricted our attention to local boundary conditions.
A fundamental problem in the design of accurate local boundary conditions is pointed out: accuracy is required to keep the
amount of reflected energy small, but at the same time allows for growing low-frequency modes. We have settled for Higdon’s
boundary conditions. Higdon proposes to include some damping to suppress the growing low-frequency modes. We show that third-order
conditions provide acceptable results for the simple scalar wave equation and the acoustic equation. In the elastic case,
an additional low-frequency growing mode may occur. This mode can be suppressed by using a dissipative boundary scheme and
by increasing the amount of damping. The increase in damping results in an increase in the amount of reflected energy, which
is larger than in the scalar case. Numerical experiments exhibit a reasonable performance, although some improvement would
be useful, particularly in the anisotropic elastic case.
This revised version was published online in July 2006 with corrections to the Cover Date. 相似文献
15.
16.
The paper deals with the application of the Boundary Element Method (BEM) to problems where the boundary geometry and conditions change with time. An example of this is the simulation of sequential excavation in underground construction. A novel approach is presented that is more efficient than the currently used methods. The method involves only one region, where the boundary conditions and the geometry are changed as time progresses. Test examples are presented and results are compared with the ones obtained from other methods. The examples show that the new method is more efficient than currently used methods. 相似文献
17.
18.
A new artificial boundary approach for transient seepage problems in unbounded domain is presented. The artificial boundary condition at the truncated boundary is derived from the analytical solutions for transient seepage problems in one dimension, including solutions, respectively, for flow in one‐dimensional infinite space and for radial flow in an infinite layer, and then it is tentatively applied for some two dimensional problems in addition to the one‐dimensional problems mentioned above. The boundary conditions derived relate the time‐dependent boundary flux with the time derivative of the hydraulic head at the truncated boundary, which makes the implementation much easier compared with the infinite element method. The accuracy and efficiency of the artificial boundary are validated by several numerical examples, which shows that the proposed boundary can give very good results for one‐dimensional transient seepage problems, as expected, whereas reasonable results can be also obtained for two‐dimensional problems, such as two‐dimensional axisymmetric flow and flow in an infinite plane. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
19.
Xiao Ma Setare Hajarolasvadi Gabriele Albertini David S. Kammer Ahmed E. Elbanna 《国际地质力学数值与分析法杂志》2019,43(1):317-338
The finite element method (FEM) and the spectral boundary integral method (SBI) have both been widely used in the study of dynamic rupture simulations along a weak interface. In this paper, we present a hybrid method that combines FEM and SBI through the consistent exchange of displacement and traction boundary conditions, thereby benefiting from the flexibility of FEM in handling problems with nonlinearities or small-scale heterogeneities and from the superior performance and accuracy of SBI. We validate the hybrid method using a benchmark problem from the Southern California Earthquake Center's dynamic rupture simulation validation exercises.We further demonstrate the capability and computational efficiency of the hybrid scheme for resolving off-fault heterogeneities by studying a 2D in-plane shear crack in two different settings: one where the crack is embedded in a high-velocity zone and another where it is embedded in a low-velocity zone. Finally, we discuss the potential of the hybrid method for addressing a wide range of problems in geophysics and engineering. 相似文献
20.
A new finite element scheme is proposed, in this paper, for solving two-dimensional wave propagation problems in multilayered soils resting on a rigid base. The multilayered soils are treated as multiple horizontal layers of lateral infinite extension in geometry. Since these horizontal layers can be truncated by two artificially truncated vertical boundaries, two high-order artificial boundary conditions are applied for propagating the incoming waves from the interior domain into the far field of the system. Both the semi-analytical method and the truncated boundary migration procedure are used to derive the high-order artificial boundary conditions, which are comprised of a physically meaningful dashpot and a generalized energy absorber. The main advantage of using the proposed finite element scheme is that the derived artificial boundary condition can be straightforwardly implemented in the finite element analysis, without violating the band/sparse structure of the conventional finite element equation. The related numerical examples have demonstrated that the proposed finite element scheme is of high accuracy in dealing with wave propagation problems in multiple horizontal layers. 相似文献