复杂边界非均质渗流场流线分布研究

建立了考虑源(汇)影响的含有不渗透区域复杂边界条件下非均质油藏稳定渗流的数学模型。利用扰动边界元方法求解数学模型,获得了地层中任意一点的压力公式.在此基础上,提出了流线场的生成方法。绘制了考虑非均质性、复杂边界和不渗透区域影响的流线分布图,并分析了流线分布的特征。通过分析表明,渗流场的非均质性和不渗透区域的存在都对流线分布存在较大的影响。利用本文方法产生的流线分布图能够较为直观地反映出油藏流体在注采井间的运动轨迹,为优化井网和注入方案提供了重要依据。     

复杂边界稠油油藏蒸汽注采过程自适应网格法的研究

将自适应网格法推广到复杂边界稠油油藏的蒸汽注采过程,针对复杂边界附近的网格提出相应的粗化算法。首先,在实施自适应网格算法之前对边界上最精细网格上的计算参数进行预处理以提高计算精度,然后,再利用同样的预处理方法对自适应网格法中边界处的各层次粗网格的渗透率进行粗化。在建立动态AMR网格系统的网格粗化准则中,仅采用油藏温度和各相饱和度的空间变化作为控制阈值,这样边界区域在相变锋面未到达时将自适应地采用粗网格进行计算。数值算例显示边界附近自适应地采用粗网格进行计算并不影响油藏数值模拟的计算精度,自适应网格法在保持计算精度的同时,大幅度提高了计算速度。     

两套节点格林元嵌入式离散裂缝模型数值模拟方法

对于原始嵌入式离散裂缝模型(EDFM),在计算包含裂缝单元的基质网格内的压力分布时采用了线性分布假设,这导致了油藏开发早期对非稳态窜流量的计算精度不足.因此,本文提出了一种两套节点格林元法的EDFM数值模拟方法.两套节点格林元法的核心思想是将压力节点与流量节点区分开,一套压力节点设置在单元顶点,另一套流量节点设置在网格边的中点,满足局部物质守恒、具有二阶精度的同时,可适用于任意网格类型.本文将两套节点格林元法与EDFM耦合,采用了非稳态渗流控制方程的边界积分形式推导了基质网格与裂缝网格之间传质量的新格式,代替了线性分布假设以提高模拟精度;此外,修正后的EDFM能适应任意形态的基质网格剖分,拓展了原始EDFM仅适用于矩形基质网格、难以考虑复杂油藏边界的局限性.研究表明:通过对比商业模拟软件tNavigator?LGR模块与原始EDFM,验证了本文模型具有较高的早期计算精度;以复杂油藏边界-缝网-SRV分区模型为例,通过对比SFEM-COMSOL商业模拟软件,验证了本文模型处理复杂问题的适应性.本文研究可用于裂缝性油藏开发动态的精确模拟.     

基于反馈力浸入边界法模拟复杂动边界流动

浸入边界法是模拟流固耦合的重要数值方法之一。本文采用反馈力浸入边界方法,对旋转圆柱和水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟。其中,固体边界采用一系列离散的点近似代替,流体为不可压缩牛顿流体,使用笛卡尔自适应加密网格,利用有限差分法进行求解。固体对流场的作用通过构造适宜的反馈力函数实现。本文首先通过旋转圆柱绕流的计算结果同实验结果进行对比,吻合较好,验证了该计算方法的可靠性。然后针对水电站水力过渡过程中水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟,得到导叶动态绕流后的流场分布特性和涡结构的演化特性。     

Intercomparison of boundary schemes in Lattice Boltzmann method for flow simulation in porous media

The Lattice Boltzmann method has been widely adopted to simulate flow in porous media. The choice of appropriate boundary schemes is essential to achieve simulation accuracy; however, the criteria for the most suitable boundary treatment in the simulation of flow in porous media flow remain unresolved. Here, three types of the most commonly used boundary conditions are tested: interpolation bounce back (IBB), partial saturated method (PSM), and immersed boundary method (IBM). The dimensionless drag of face-centered cubic (FCC) sphere array and the dimensionless permeability of a random closely packed (RCP) sphere array are calculated and compared at different viscosities and resolutions. In the FCC sphere array case where spheres are not contacted, the IBB and PSM exhibit the same accuracy and both are of the second-order convergence rate. The IBM is less accurate and is of the first-order convergence rate. In the RCP sphere array case where the spheres are contacted, the IBB shows finer results and a second-order convergence rate. PSM underestimates the dimensionless permeability and increases resolution only slightly improved the results of PSM. The IBM overestimates the dimensionless permeability. These results indicate that among the three methods, the IBB is the most accurate. The PSM has the same accuracy as the IBB when sediments are not contacted; however, it loses its accuracy in the simulation of flow in closely packed porous media. This work could serve as a benchmark for further research in choosing the most appropriate method in the simulation of flow in porous media.     

流体饱和多孔介质黏弹性动力人工边界

基于Biot流体饱和多孔介质本构方程,分别考察具有辐射阻尼性质的外行柱面波和球 面波在圆柱面和球面人工边界上引起的法向、切向应力的表达式. 在应力表达形式上,固相 介质和孔隙流体的法向和切向应力都是由两项组成,它们分别与质点的位移和速度成正比, 因此,可在人工边界的法向和切向设置连续分布的并联弹簧------黏滞阻尼器,用来模拟人工边 界以外的无限域介质对来自有限计算域的外行波动的能量吸收作用,从而形成了流体饱和多 孔介质的黏弹性动力人工边界. 流体饱和多孔介质的黏弹性动力人工边界可方便地与大型通 用软件结合,用于分析饱和土中复杂的结构-地基动力相互作用问题. 算例表明流体饱和多 孔介质黏弹性动力人工边界具有较好的精度和稳定性.     

自由面势流问题的域外奇点边界元法及其数值误差分析

讨论了域外奇点边界元法在自由面势流问题计算中的作用，并以连续及离散Fourier分析对该方法(就m阶面元的一般情况)进行数值误差分析，导出了计及面元阶数、奇点至自由面垂向距离、配置点移动、差分格式等因素影响的数值误差一般表达式。从理论上证明了自由面势流问题计算中采用域外奇点法可改善离散产生的数值色散误差并能结合配置点前移(向上游)等方法以数值满足辐射条件。     

Disturbed boundary layer flow with local time-dependent surface heating

Local flows in a laminar boundary layer flowing over surface heating elements are investigated. Mathematical models of disturbed flows are constructed on the basis of an asymptotic analysis and the similarity parameters are determined. The time-dependent local heating regimes ensuring control of separation and flow stability in the boundary layer are studied. The results of a numerical and analytic analysis are obtained.     

边界元理论在复杂外边界油藏水平井渗流中的应用

受构造作用的影响,实际油藏的外边界往往是复杂多样的.本文从渗流理论出发建立了复杂外边界油藏水平井渗流数学模型,并采用Lord Kelvin点源解、贝塞尔函数积分和泊松叠加公式等方法求解了复杂外边界油藏水平井的边界元基本解,利用边界元的理论建立了复杂外边界油藏水平井井底压力响应数学模型.通过计算得到了无因次压力和压力导数双对数理论图版,并在其基础上分析了复杂外边界油藏水平井渗流特征及其影响因素.     

Dynamic solution of hyperbolic cooling tower with nonsymmetric imperfections by finite element method

In this paper, an approximate method is established to solve the couple problem between different harmonic waves in terms of the needs of nonsymmetric imperfections of cooling tower. By means of this method, the frequencies and response of rotational shell with local geometric imperfections are analysed and calculated.     

Lattice Boltzmann modeling of pore‐scale fluid flow through idealized porous media

In order to capture the hydro‐mechanical impacts on the solid skeleton imposed by the fluid flowing through porous media at the pore‐scale, the flow in the pore space has to be modeled at a resolution finer than the pores, and the no‐slip condition needs to be enforced at the grain–fluid interface. In this paper, the lattice Boltzmann method (LBM), a mesoscopic Navier–Stokes solver, is shown to be an appropriate pore‐scale fluid flow model. The accuracy and lattice sensitivity of LBM as a fluid dynamics solver is demonstrated in the Poiseuille channel flow problem (2‐D) and duct flow problem (3‐D). Well‐studied problems of fluid creeping through idealized 2‐D and 3‐D porous media (J. Fluid Mech. 1959; 5 (2):317–328, J. Fluid Mech. 1982; 115 :13–26, Int. J. Multiphase Flow 1982; 8 (4):343–360, Phys. Fluids A 1989; 1 (1):38–46, Int. J. Numer. Anal. Meth. Geomech. 1999; 23 :881–904, Int. J. Numer. Anal. Meth. Geomech. 2010; DOI: 10.1002/nag.898, Int. J. Multiphase Flow 1982; 8 (3):193–206) are then simulated using LBM to measure the friction coefficient for various pore throats. The simulation results agree well with the data reported in the literature. The lattice sensitivity of the frictional coefficient is also investigated. Copyright © 2010 John Wiley & Sons, Ltd.