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

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

多尺度嵌入式离散裂缝模型模拟方法

天然裂缝性油藏和人工压裂油藏内裂缝形态多样,分布复杂,传统的离散裂缝模型将裂缝作为基岩网格的边界,采用非结构化网格进行网格划分,其划分过程复杂,计算量大。嵌入式离散裂缝模型划分网格时不需要考虑油藏内的裂缝形态,只需对基岩系统进行简单的网格剖分,可以大大降低网格划分的复杂度,从而提高计算效率。然而,在油藏级别的数值模拟和人工压裂裂缝下的产能分析中,仍然存在计算量巨大、模拟时间过长的问题。本文提出嵌入式离散裂缝模型的多尺度数值计算格式,使用多尺度模拟有限差分法研究嵌入式离散裂缝模型渗流问题。通过在粗网格上求解局部流动问题计算多尺度基函数,多尺度基函数可以捕捉裂缝与基岩间的相互关系,反映单元内的非均质性,因此该方法既有传统尺度升级法的计算效率,又可以保证计算精度,数值结果表明这是一种有效的裂缝性油藏数值模拟方法。     

裂缝介质渗流多尺度混合有限元数值模拟研究

针对裂缝介质具有多尺度特点,建立了Darcy/Stokes-Brinkman多尺度耦合模型,采用多尺度混合有限元方法,对裂缝介质渗流问题进行了研究.阐述了多尺度混合有限元方法的基本原理,并推导得到Darcy/Stokes-Brinkman方程的多尺度混合有限元计算格式.数值计算结果表明,大尺度Darcy模型能够捕捉到小尺度上裂缝网络渗流特征;与网格粗化、传统有限元方法相比,多尺度混合有限元方法的基函数具有能反映单元内参数变化的优点,在保证计算精度的同时能够减少计算量,对于裂缝油藏具有良好的适用性.     

基于XFEM-MBEM的嵌入式离散裂缝模型流固耦合数值模拟方法

离散缝网的表征与模拟是目前国内外研究的热点. 在非常规油气开发过程中, 由于地应力场的存在会对裂缝的流动属性产生显著影响, 若将裂缝视为静态对象, 与矿场数据会出现极大偏差, 因此要基于动态裂缝做更深入的研究. 本文针对致密油藏应力场?渗流场耦合力学问题, 提出了一种高效的混合数值离散化方法, 其中采用扩展有限元法 (XFEM) 求解岩石的弹性形变, 采用了混合边界元法 (MBEM) 精确计算基岩与裂缝间的非稳态窜流, 这两种数值格式是完全耦合的, 并对整体计算格式的时间项进行了全隐式求解, 可准确表征致密油藏开采过程中的裂缝变形及流体流动机理. 此外, 本文采用了嵌入式离散裂缝前处理算法显式表征大尺度水力压裂缝, 并考虑了支撑剂的作用; 采用了双孔有效应力原理和双重介质隐式裂缝表征方法, 可捕捉基质与小尺度天然裂缝的动态信息; 由此, 本文所提出的混合模型综合表征了基质?天然裂缝?水力压裂缝共同组成的致密油藏复杂渗流环境, 并通过几个实例论证了模型的准确性, 研究表明: 对致密油藏压裂水平井进行产能评价时, 应力场所引起渗流参数的改变及裂缝开度降低的影响不可忽略. 本文研究可为非常规油气资源的开发提供理论指导.      

压裂页岩气藏多尺度耦合流动数值模拟研究

在碳达峰的国策背景之下,页岩气成为传统能源向绿色清洁低碳能源转型的重要过渡和能源支点.压后页岩气藏流体流动力学成为高效开发页岩气的关键力学问题.文章将小尺度低导流天然裂缝等效升级为连续介质,建立有机质-无机质-天然裂缝三重连续介质模型,同时对大尺度高导流裂缝采用离散裂缝模型刻画,嵌入天然裂缝连续介质中,构建多重连续/离散裂缝模型.综合考虑吸附气的非平衡非线性解吸附和表面扩散,自由气的黏性流和克努森扩散,给出页岩气在多尺度复杂介质中的非线性耦合流动数学模型.提出多尺度扩展有限单元法对离散裂缝进行显式求解,创新性构建三类加强形函数捕捉离散裂缝的局部流场特征,解决了压后页岩海量裂缝及多尺度流动通道的流动模拟难题.文章提出的模型和方法既能准确刻画高导流裂缝对渗流的影响,又克服了海量多尺度离散裂缝导致计算量增大的问题.通过算例展示了压后页岩各连续介质的压力衰减规律,发现裂缝中自由气、有机质中自由气、无机质中吸附气依次滞后的压力(浓度)扩散现象,重点分析了吸附气表面扩散系数、自由气克努森扩散系数、天然裂缝连续介质渗透率和吸附气解吸附速率对页岩气产量的影响.文章重点解决压后页岩多尺度流动通道的表征和...     

基于离散裂缝的多段压裂水平井数值试井模型及应用

水平井压裂技术已经成为开发低渗透油气藏、页岩气藏和致密气场等非常规油气藏的关键技术。基于离散裂缝模型,对裂缝进行简化,建立了二维多段压裂水平井有限导流数值试井模型,利用有限元方法求解模型,获得多段压裂水平井试井理论曲线和压力场特征。分析表明:多段压裂水平井的试井理论曲线一共分为七个阶段:井筒储存段、裂缝线性流段、裂缝-地层双线性流段、裂缝干扰段、地层线性流段、系统径向流段和边界作用段,其中裂缝-地层双线性流段和裂缝干扰是其典型特征。分析了裂缝数量、裂缝间距、裂缝不对称、裂缝不等长和裂缝部分缺失等因素对试井理论曲线的影响,结果表明:裂缝数量和裂缝间距对试井理论曲线的影响最大。较多的裂缝、较大裂缝间距、对称的裂缝和等长的裂缝有利于降低压裂水平井井底的流动阻力,提高产能。将建立的数值试井模型应用于四川盆地一口多段压裂水平井的压力恢复测试的数值试井解释,结果表明:本文建立的模型可以较好的拟合压力恢复测试数据,可以获得裂缝的导流能力和裂缝长度,为压裂效果评价和压裂设计提供指导。      

基于离散缝洞网络模型的缝洞型油藏混合多尺度有限元数值模拟

针对缝洞型油藏具有多尺度特征,基于离散缝洞网络模型,建立了Darcy/Stokes-Brinkman多尺度耦合数学模型,采用多尺度混合有限元方法,对缝洞型油藏的流体流动问题进行了研究。阐述了多尺度混合有限元方法的基本原理,推导得到了Darcy/Stokes-Brinkman多尺度模型的多尺度混合有限元计算格式。数值计算结果表明,在大尺度模型上进行计算能够捕捉到小尺度上的流动特征;相对于传统有限元,多尺度混合有限元能够捕捉小尺度上的非均质特征而具有更高的计算精度,在保证计算精度的同时能够减少计算量。     

可压缩流体流动多尺度混合有限元数值方法研究

可压缩流体是天然油藏中广泛存在的一种流体,研究其在多孔介质中的渗流规律对于油藏开发具有重要意义。本文采用多尺度混合有限元方法,对可压缩流体渗流问题进行了研究。考虑流体的可压缩性以及介质形变,推导得到了可压缩流体渗流问题的多尺度计算格式。数值计算结果表明,多尺度混合有限元适于求解非均质性和可压缩流问题,具有节省计算量、计算精度高等优势,对于实际大规模油藏模拟具有重要意义。     

基于连接元法的多尺度裂缝性油藏数值模拟

针对油藏不同尺度复杂几何特征描述和动态连通性识别等难题,近年来发展了一种基于非欧物理连通网络具有无网格特征的油藏数值模拟连接元方法.文章将连接元法推广到裂缝性油藏,从流体流动的角度,利用连接单元将油藏离散为物理连通网络.根据节点物性参数、影响域半径和加权最小二乘法给出了压力扩散项的广义差分近似.结合物质守恒方程计算节点控制体积、基质节点间传导率、裂缝节点间传导率以及基质节点与裂缝节点间传导率,从而构建渗流控制方程组的全隐式离散格式,求解压力、饱和度以及含水率等生产动态参数.引入图论深度优先搜索算法,基于每个时间步求解的节点间压力梯度,计算各时间步注入井的劈分系数,定量表征井节点间的流动关系和连通性.算例验证表明,相较基于网格体系的传统方法,该方法能够自由灵活地刻画包括裂缝复杂分布、不规则油藏边界在内的复杂油藏几何,在粗化模型情况下能够保留更丰富的流动拓扑结构,实现计算精度和计算效率的更优平衡,能更好满足实际大规模裂缝性油藏的生产动态模拟预测需求,同时为具有多尺度几何特征的裂缝性油藏及复杂边界油藏的数值模拟提供了新思路.     

双相介质P波波动方程小波域多尺度模拟

相比于单相介质理论而言,双相介质理论更接近实际地层的真实情况,因此在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用。传统的波动方程数值解法由于本身固有的不足不利于求解诸如双相介质波动方程等复杂的非线性和不规则性问题;而小波方法则由于自身良好的特性可以用来构建解决此类问题的自适应性算法。本文详细推导了双相介质P波波动方程的有限差分矩阵表示形式,利用小波变换将其转移到小波域,设置阈值形成更为稀疏的迭代矩阵以构建自适应算法,从而达到减少计算量,增加地震波场数值模拟灵活性和准确性的目的。地球物理勘探的数值模拟实例验证了方法的有效性。     

基于离散元的杆件断裂全过程数值模拟

杆件的断裂会涉及到大变形、非线性以及不连续等问题,通常的数值计算方法模拟这种复杂力学行为具有局限性。本文基于颗粒离散元法DEM,将接触粘结处的分布式弹簧用梁纤维进行等效,提出了一种适于结构弹塑性分析的DEM纤维梁模型,然后在此基础上构建了构件断裂模拟算法以及纤维破环准则。将该模型应用于悬臂梁结构,模拟了悬臂梁从弹性到弹塑性阶段,再到断裂破坏的全过程,数值模拟得到的结构响应和截面开裂破坏形态均较合理。最后将该方法应用于单层网壳倒塌破坏模拟,并与网壳振动台倒塌试验进行对比,结果表明,数值模拟得到的杆件断裂过程及结构倒塌模式与试验现象一致,验证了该模型的正确性和适用性。     

稀薄流非线性模型方程离散速度坐标法有限差分解

从一般非线性Bo ltzm ann方程出发,发展并实现了一套适于大范围K nudsen数稀薄流问题数值模拟的统一算法。采用BGK模型和Shakov模型近似碰撞项,进而引入两个二速度无量纲简化分布函数,通过关于分子速度第三分量取矩积分,将三速度单一模型方程变换为二速度微分方程组。基于G auss-H erm ite积分公式和正交多项式G auss积分公式,借助离散速度坐标法消除简化模型方程对分子速度空间的连续依赖性,从相空间到物理空间得到一组带源项双曲守恒离散方程,并给出其显式和隐式二阶迎风TVD有限差分解。以二维圆柱A r气体超声速绕流算例,验证了数值算法的有效性,比较分析了漫反射和镜面反射两种气体分子壁面反射模型的计算结果。     

离散元与壳体有限元结合的多尺度方法及其应用

在深入研究复杂结构和非均质材料冲击响应和破坏机理的过程中,往往遇到多尺度计算问题。本文尝试建立三维离散元与壳体有限元结合的多尺度方法用于处理圆柱壳问题,该方法采用三维离散元对感兴趣的局域进行局部模拟,利用平板壳体有限元进行整体模拟,采用一种特殊的过渡层使离散元区和有限元区能很好的衔接。我们将这一方法应用于激光辐照下充压柱壳的热/力耦合冲击破坏响应,得到的模拟结果与文献报道有较好的吻合。     

Extended finite element modeling of production from a reservoir embedded with an arbitrary fracture network

A numerical scheme based on the eXtended Finite Element Method (XFEM) is proposed to simulate complex fluid flow in a fractured porous reservoir. By enriching the elements fully cut by the fracture and the near‐tip region, the flow mechanism including the tip flux singularity can be exactly represented in the XFEM formulation. Fluid transfer between the matrix and the fractures can be easily coupled, and XFEM also overcomes the sensitivity to the mesh used in the traditional unstructured discretizations, regardless of the complexity of the fracture network. The method is validated for a simple case by the exact analytical solution. Results are compared between XFEM and FEM. Case studies are presented to illustrate the power, efficiency, accuracy, and flexibility of the proposed method for simulating transient productive flow in reservoirs with complex fracture networks.     

两个新型的八结点块体单元—三维离散算子差分法

给出了弹性力学三维问题的离散算子差分法 ,讨论离散算子差分法在三维问题中的特点 ,意在为该方法的进一步发展提供依据 ,为应用弱形式进行数值求解的研究提供参考。本文从弹性力学平衡方程更为一般的弱形式出发 ,给出了含边界参数的弱形式方程。由该方程不仅可以得到有限元法 ,还可得到离散算子差分法。给出了两个八结点块体单元 ,虽然单元中位移函数是非协调的 ,不需特殊处理便可保证离散格式收敛 ,并对单元位移有十分好的反映能力。     

Numerical Simulation of Single-Phase and Multiphase Non-Darcy Flow in Porous and Fractured Reservoirs

A numerical method as well as a theoretical study of non-Darcy fluid flow through porous and fractured reservoirs is described. The non-Darcy behavior is handled in a three-dimensional, multiphase flow reservoir simulator, while the model formulation incorporates the Forchheimer equation for describing single-phase or multiphase non-Darcy flow and displacement. The non-Darcy flow through a fractured reservoir is handled using a general dual-continuum approach. The numerical scheme has been verified by comparing its results against those of analytical methods. Numerical solutions are used to obtain some insight into the physics of non-Darcy flow and displacement in reservoirs. In addition, several type curves are provided for well-test analyses of non-Darcy flow to demonstrate a methodology for modeling this type of flow in porous and fractured rocks, including flow in petroleum and geothermal reservoirs.     

缝洞型多孔介质中多相流的有限体积法数值模拟

本文研究的碳酸盐岩油藏储集体属于缝洞型多孔介质.这类缝洞型多孔介质由裂缝、溶蚀孔洞和低孔隙度低渗透率的基岩组成.裂缝是空隙流体流动的主要通道;溶蚀孔洞大小从几厘米到数米不等,渗透率和孔隙度都很高,是流体主要的储集空间.由于缝洞型多孔介质空隙空间的复杂性和强非均质性,数值计算中基本控制方程的空间离散应采用非结构化网格的计算模型.本文采用有限体积法模拟缝洞型多孔介质中多相流体的流动,并给出了相应的单元中心格式有限体积法的计算公式.裂缝介质和溶洞介质中单元间多相流体的流动考虑为高速非达西流,其质量通量采用Forchheimer定律计算.非线性方程的离散选取全隐式格式,并采用Newton-Raphson迭代进行求解.通过两个二维模型注水驱油的数值模拟,验证了本文方法的有效性.     

MODEL AND METHOD OF WELL TEST ANALYSIS FOR WELLS WITH VERTICAL FRACTURE

NomenclatureB———volumecoefficientofformationfluid,dimensionless;C———wellborestoragecoefficient,m3 MPa;Cft ———comprehensivecompressioncoefficientoffluidinfracture ,1 MPa;Ct ———comprehensivecompressioncoefficientofformationfluid,1 MPa ;h———heightofformation ,m ;K———formationpermeability ,μm2 ;Kf ———fracturepermeability ,μm2 ;KfWf ———fractureconductivity ,μm2 ·cm ;Lf ———halfoffracturelength,m ;N———thetotalnumberofselectedtimepoint,dimensionless;P———for…