首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 550 毫秒
1.
《地下水》2020,(2)
目前,水库调洪计算有试算法和半图解法,由于两种方法迭代试算或辅助线重复插补求解非线性超越方程的速度较慢,工作量大,故本文将水平衡方程和水库流量曲线归纳为一阶常微分方程,根据边界条件,采用四阶龙格库塔法求解水库洪水调度过程,该方法是一种求解微分方程的数值方法,便于程序设计,精度较高,该方法在水库洪水调度中具有一定的适用性。  相似文献   

2.
采用颗粒离散单元法进行动力计算时,人工截断边界上需设置吸收边界条件,以防止波的反射。鉴于颗粒离散单元数值计算模型的人工边界上颗粒单元半径大小不一、边界面凸凹不平,在连续介质的黏性、黏弹性、自由场边界条件方程基础之上,推导出适用于离散介质的等效方程。在离散介质的黏性边界条件等效方程中引入微调系数,提出比值迭代法以快速确定其最优值,以实现对波的最佳吸收。采用二维颗粒离散单元计算软件PFC2D,分别建立黏性、黏弹性、自由场边界条件相关数值分析模型,探讨颗粒分布模式对黏性边界上颗粒单元半径、速度分布及比值迭代过程的影响;采用外源波动算例及经典Lamb问题算例验证黏弹性边界设置方法的正确性;通过隧洞算例检验提出的自由场边界条件设置方法的正确性。  相似文献   

3.
研究了二阶声波方程频域PML边界条件和频域变网格并行计算技术。PML边界是一种较为理想的吸收边界方法,多用在求解时域应力速度方程中,但对于频域声波正演,二阶位移方程更常用。从一阶声波方程PML吸收边界条件导出频域二阶位移方程PML边界条件,模拟算例得到的频率切片、时间切片和地震记录对比都说明该边界条件吸收效果很好。频域单炮正演不同频率间是独立的,据此低频部分采用大网格计算,高频采用小网格,实现变网格步长计算技术,这是较时间域正演的一个优势,在保证模拟质量的同时,减少计算量和内存消耗。  相似文献   

4.
在水文地质计算中,正确认识和处理边界条件是至关重要的。文献《一维优化搜索联合有限元方法计算渗流二类边界》(以下简称《边界》)利用一维搜索结合有限单元方法计算流量和水位随时间和随空间变化的边界条件(以下简称变边值问题)。笔者认为该文所述方进在反映变边值问题的特点方面有值得商榷的地方。本文分析了《边界》方法的不足。并提出了适合于该文所述变边值问题的计算方法。  相似文献   

5.
弹性波数值模拟的混合边界与频散抑制   总被引:1,自引:0,他引:1  
针对刚性边界条件、衰减边界条件与透射边界条件对弹性波场的影响及计算中产生的数值频散现象,将衰减边界与透射边界双重作用于人工边界,引入了混合边界条件,并在计算中采用了微量阻尼的方法,结果表明,混合边界条件能更好地吸收人工边界的反射;数值频散的阻尼算法能大大减弱由网格剖分、计算误差造成的频散现象。该方法可有效提高波场模拟精度。   相似文献   

6.
EFG法在土体固结中的应用   总被引:4,自引:0,他引:4  
EFGM作为一种新的数值计算方法,具有只需节点信息而无段单元的特性,故在解固结方程方面有很大的优势。在计算分析中,此法容易构造固结方程的EFGM刚度矩阵和处理不同边界条件。对单面排水等条件的计算结果表明,EFGM在解决固结变形问题上,精度较高,处理边界准确。  相似文献   

7.
在降雨入渗模拟中,现有研究无法准确实现复杂降雨条件下降雨边界在流量与压力两类边界之间的双向动态转换。基于饱和?非饱和渗流理论,将降雨补给流量与实际入渗流量的差值作为两类边界动态转换的判别条件,改进了单重渗透介质降雨边界动态转换控制方程,并进一步提出了双重渗透介质降雨边界处理方法及控制方程,在此基础上进行了数值试验及工程实例验证。结果表明:改进的降雨边界动态转换控制方程克服了传统降雨边界只能实现流量到压力边界单向转换的局限性,可准确实现流量边界与压力边界的双向动态实时转换;案例边坡的降雨入渗模拟结果说明了改进后降雨边界的正确性和工程适用性。  相似文献   

8.
地下水数值模拟中人为边界的处理方法研究   总被引:10,自引:1,他引:9  
边界条件预测是地下水数值模拟中的关键问题。本文针对人为边界受自然和人为因素影响,边界流量变化复杂难于预测的问题,讨论了人为边界条件的处理方式和实现方法,并将其应用到具有人为边界的二维潜水含水层系统。预测结果表明,受区内人工开采影响,边界流量呈逐年上升趋势,部分边界由流出边界转化为流入边界,结果符合地下水系统变化趋势,提高了预测精度。  相似文献   

9.
利用格林公式及点源场的边界条件,将地下半空间位函数的体积积分变为地面位函教的积分方程以后,用边界元方法解此积分方程便可得到地面各节点的电位值。其计算精度和速度取决于剖分方式和节点密度。本文采用变网格的剖分方式,用较少的点实现了计算,并兼顾了精度和速度。  相似文献   

10.
闫富有  吴义章  郭院成 《岩土力学》2011,32(Z2):604-609
基于有限压缩层地基模型和Reissner板的边界积分方程,考虑柱荷载作用区域的形状和筏板的横向剪切变形效应,建立了有限压缩层地基上厚筏基础与地基相互作用分析的边界元方程和系统的数值方法。对于弹性半空间模型,其柔度方程可视为有限压缩层模型的特殊情况。把筏板作为自由边界条件处理,被剖分为一系列三角形或矩形网格,假设基底反力在网格内均匀分布,以便与现有的地基沉降计算模式相一致。计算表明,虽然基底反力在内部网格相接处不连续,但并不影响计算结果,反而消除了边界基底压力计算值过大的现象。将该方法与其他方法的计算结果进行比较,显示了该方法的有效性。计算结果表明,对于实际复杂的筏板基础,无需划分太多单元即可得到较高的计算精度  相似文献   

11.
油井开采过程中油层变形的流固耦合分析   总被引:7,自引:0,他引:7  
在油气开采过程中,随着油气的不断采出,必然造成孔隙流体压力的逐渐降低,由此导致储层岩石骨架的有效应力增大,使得油层产生变形或压实。当油层产生变莆或压实时,对油气生产将造成不利影响。比如:使得油藏的渗透率降低,继而使油井的产能降低,同时,油层的变形直接影响着油井和套管的变形与破坏等等。敢开采过程中油层的变形可以描述为三维变形与三维流体流动场的耦合问题,利用可变形多孔介质中流体渗流的流固耦合有限元数值  相似文献   

12.
储层应力场的数值模拟   总被引:5,自引:0,他引:5  
该文给出了油田储层应力场的两种数值模拟方法。一种是利用现场实测地应力资料,借助空间膜单元模型,反演远场应力边界条件,从而计算出储层应力场。另一种是借助于活动的最小二乘法直接插值得到应力场。针对鄯善油田三间房油藏,用所提出的方法计算模拟了储层应力场,得到了比较满意的结果,同时论述了这两种数值模拟方法的可行性。   相似文献   

13.
林永生  陈胜宏 《岩土力学》2013,34(4):1191-1196
提出一种基于有限元计算的边坡三维滑裂面搜索的方法,根据有限元计算的应力结果,计算可能滑裂面的安全系数,然后用遗传算法搜索出最危险滑裂面,并得出相应的安全系数作为边坡的安全度评价指标。该方法的优越性是只需进行一次有限元计算,不需要假定滑裂面,可以用于具有复杂地貌和地质构造的三维边坡。通过算例,验证该方法的可行性和程序的正确性,以大岗山水电站库区边坡工程实例验证该方法的实用性。  相似文献   

14.
油藏多相渗流的面向对象有限元程序设计   总被引:3,自引:0,他引:3  
董平川  韩德金  牛彦良  李莉 《岩土力学》2009,30(4):1115-1121
近年来,面向对象的编程技术越来越受到人们的重视,并用来解决和模拟各种物理问题和工程问题,特别是在有限元方法编程和模拟中得到了广泛地应用。介绍了面向对象的程序设计方法和它与有限元程序相结合的原理和优点。在此基础上提出了面向对象有限元方法在油藏多相渗流中的应用实例,并按照面向对象的程序设计方法,基于Windows XP平台采用Visual C++.NET语言,建立了有关描述多相渗流的有限元模型的类,给出了类的描述和它的实现方法。建立的相关类及其实现方法主要包括处理矩阵计算的类、节点类、单元类、材料类、形函数类等。据此编制了油藏多相渗流的有限元分析的数值计算程序,并给出了一个数值实例。  相似文献   

15.
堆积层滑坡的岩土体渗透系数具有一定的不确定性,且渗透系数是饱和-非饱和渗流分析的重要参数,开展考虑其空间变异性的库岸堆积层滑坡渗流变形分析具有重要意义。以三峡库区中的白水河滑坡为研究对象,基于地面核磁共振技术获取的岩土体渗透系数,分析滑坡体渗透系数的空间变异特征,采用半变异函数方法求得滑坡体渗透系数的竖直波动范围,在此基础上建立渗透系数的非平稳随机场模型。以非侵入式随机有限元的方式开展库水升降两种工况下不确定模型与确定模型的流固耦合模拟,分析两种模型的渗流场、位移变形特征及其差异。结果表明:相比于确定模型,不确定模型孔压改变的滞后性更为明显,且库水下降工况下整体的变形更大,若忽略滑体渗透系数的非平稳空间变异特征将会低估滑坡的实际变形。  相似文献   

16.
Based on the theory of double-porosity, a novel mathematical model for multiphase fluid flow in a deforming fractured reservoir is developed. The present formulation, consisting of both the equilibrium and continuity equations, accounts for the significant influence of coupling between fluid flow and solid deformation, usually ignored in the reservoir simulation literature. A Galerkin-based finite element method is applied to discretize the governing equations both in the space and time domain. Throughout the derived set of equations the solid displacements as well as the fluid pressure values are considered as the primary unknowns and may be used to determine other reservoir parameters such as stresses, saturations, etc. The final set of equations represents a highly non-linear system as the elements of the coefficient matrices are updated during each iteration in terms of the independent variables. The model is employed to solve a field scale example where the results are compared to those of ten other uncoupled models. The results illustrate a significantly different behaviour for the case of a reservoir where the impact of coupling is also considered. © 1997 by John Wiley & Sons, Ltd.  相似文献   

17.
含天然气水合物沉积物分解过程的有限元模拟   总被引:1,自引:0,他引:1  
温度和压力的变化会引起含天然气水合物沉积物的分解,其过程伴随着相态转换、孔隙水压力和气压力耗散、热传导、骨架变形等过程的相互耦合作用。基于多孔介质理论建立了描述含天然气水合物沉积物分解过程的数学模型,考虑了水合物分解产生的水、气流动、水合物相变和分解动力学过程、热传导、骨架变形等过程的耦合作用。基于有限元法,建立了模拟水合物分解过程的数值模型,并编制了计算机分析程序。通过对降压法和升温法开采过程的数值模拟,揭示了在水合物分解过程中沉积物储层的变形、压力、温度等因素的变化规律。结果表明:降压法和升温法都会导致储层变形以及产生超孔隙压力,但两种方法作用效果不同;同时,水合物分解过程包含渗流及热传导作用。  相似文献   

18.
Reservoir simulators typically use cell‐centered finite volume schemes and do not model directly the coupling of the flow processes with the geomechanics. Coupling of geomechanics with fluid flow can be important in many cases, but introducing fully coupled geomechanical effects in those simulators is not a trivial issue, because the geomechanics is better done by using the Galerkin vertex‐centered finite element methods by which the solid displacements are computed at the vertices of the cells. This creates difficulties in interfacing cell variables with nodal variables. Uncoupled or loosely coupled models are used by many researchers/practitioners by which a reservoir model is coupled to a geomechanical model by staggering in‐time flow and deformation via a sophisticated interface that repeatedly calls first flow and then mechanics. The method therefore requires projection of the reservoir cell variables onto the nodes of the geomechanics Galerkin finite element mesh. In this note, we attempt to quantify the errors associated with cell to node projection operations. For that purpose, we use a simple model of the pressure equation for a heterogeneous medium in one dimension. We are able to derive the exact analytical solution for this problem for both nodal and cell pressures. This allows us to compute the errors due to projection analytically, function of meshing refinement and permeability field variations. We compute upper and lower bounds for the errors, and analyze their magnitude for a variety of cases. We conclude that, in general, cell to node projection operations lead to substantial errors. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

19.
Thermal oil recovery processes involve high pressures and temperatures, leading to large volume changes and induced stresses. These cannot be handled by traditional reservoir simulation because it does not consider coupled geomechanics effects. In this paper we present a fully coupled, thermal half‐space model using a hybrid DDFEM method. A finite element method (FEM) solution is adopted for the reservoir and the surrounding thermally affected zone, and a displacement discontinuity method is used for the surrounding elastic, non‐thermal zone. This approach analyzes stress, pressure, temperature and volume change in the reservoir; it also provides stresses and displacements around the reservoir (including transient ground surface movements) in a natural manner without introducing extra spatial discretization outside the FEM zone. To overcome spurious spatial temperature oscillations in the convection‐dominated thermal advection–diffusion problem, we place the transient problem into an advection–diffusion–reaction problem framework, which is then efficiently addressed by a stabilized finite element approach, the subgrid‐scale/gradient subgrid‐scale method. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

20.
土体固结弹塑性分析的参数二次规划理论及有限元解   总被引:1,自引:0,他引:1  
基于广义Biot理论对土体弹塑性固结过程进行求解,建立了问题对应的参变量变分理论,并给出了数学证明,对此基础上推导了有限元分析列式,问题的求解最终化为参数二次规划问题,本文提出的方法适用于固结弹塑性分析的关联与非关流流动问题,也可处理各类软化问题。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号