050二维非达西渗流问题


图 1-1-33 非达西渗流速度场云图分布
图 1-1-34 达西渗流速度场云图分布
图 1-1-35 非达西渗流场压力云图分布
图 1-1-36 达西渗流场压力云图分布
图 1-1-37 非达西流右边界中点压力随时间的演进图
图 1-1-38 达西流右边界中点压力随时间的演进图
图 1-1-39 非达西流右边界中点速度随时间的演进图
图 1-1-29 选择渗流速度云图
图 1-1-30 渗流速度场云图分布
(4)点击菜单 【View results】 【Graphs】 【Point evolution】 【unoda0】 , 如图 1-1-31 所示,显示得到右边界中点压力随时间的演进图,如图 1-1-32 所示。
Hale Waihona Puke 图 1-1-31 选择右边界中点压力随时间的演进图
1.1.2 求解步骤
1.1.2.a 选择项目
P x P y
(1)启动 SciFEA,选择“项目”—>“新建项目”菜单或选择新建项目按钮 弹出如图 1-1-2 所示的对话框。
图 1-1-1
选择项目类型对话框
(2)点击“问题类型”栏中的“非达西渗流 Forchheimer 模型”选项。如图 1-1-1 所示。 (3)点击“OK”按钮完成项目类型的选择。
K V (1 V V ) P P V c t
对于二维情况写成分量形式如下:
K K V y (1 V x2 V y2 ) V x V y P c x y t Vx (1 Vx2 V y2 )
图 1-1-27 选择渗流压力水头云图
图 1-1-28 渗流场压力云图分布
(3)点击菜单【View results】-【Contour fill】-【unodb0】-【|unodb0|】 ,如图 1-1-29 所 示 , 显 示 得 到 渗 流 速 度 场 云 图 分 布 , 如 图 1-1-30 所 示 。
1.1.2.c 建模、设置材料属性和施加边界条件
(1) 启动 GID 以创建模型。 点击菜单选择 “前处理” —> “用户自定义模型” , 如图 1-1-6 所示。或者单击工具条中的 窗口。 按钮弹出如图 1-1-7 所示前处理初始化
图 1-1-6
启动前处理
图 1-1-7 前处理初始化窗口
(1)建模。a. 如图 1-1-8 所示,点击【Geometry】—【Create】—【Line】,
图 1-1-10 选择连接
图 1-1-11 形成的封闭线段 如下图 1-1-12 所示,点击【Geometry】—【Creat】—【NURBS surface】—【By contour】菜单项,选择图 1-1-11 中的所有线段,形成如图 1-1-11 所示的几何 面。
图 1-1-12 选择创建面的菜单项
b.选择问题类型。如图 1-1-13 所示,点击菜单【Data】-【Problem Type】【SCI】 ,将弹出如图 1-1-14 所示对话框,点击“确定”按钮即可。
图 1-1-13 选 sci
图 1-1-14
c.定义材料特性、边界条件。如图 1-1-15 点击【Data】—【Conditions】 ,弹 出如图 1-1-16 的 Conditions 对话框。 ① 定义材料。点击其中的 ,加入关于面的条件,如材料号、初值、边值等, 点击如图 1-1-16 所示的 Conditions 对话框中的下拉菜单,选中 Surface-aeq4, 在 mat num 一栏中填入 1,然后点击 Conditions 对话框中的 Assign 按钮,选中 整个区域,然后,点击 Conditions 对话框中的 Finish 按钮。这样,面材料号施 加完毕。点击 Draw-Color 后得到如图 1-1-17 所示的材料属性结果。
不同角度,不同方式来输出计算结果。 (1) 点击“后处理”—>“渗流”—>“二维渗流” ,如图 1-1-26 所示。
图 1-1-26 进入后处理结果分析
(2)点击菜单【View results】-【Contour fill】-【unoda0】 ,如图 1-1-27 所示, 显示得到渗流场压力云图分布,如图 1-1-28 所示。
图 1-1-23
1.1.2.e 有限元计算
转化数据消息框
(1) 点击菜单选择“计算”—>“渗流”—>“二维渗流” ,如图 1-1-24 所示。
或者单击工具条中的
按钮弹出如图 1-1-25 所示计算模拟窗口。
图 1-1-24
启动有限元计算
图 1-1-25
计算模拟窗口
1.1.3 结果分析 对计算结果进行分析属于后处理,GID 提供了丰富的后处理操作,可以从
③定义边界条件。如下图 1-1-19 所示,点击 Conditions 对话框中的 ,加 入关于线的条件,如材料号、初值、边值等。如图 1-1-19 所示,保持 line-nda 不变,保持 u-l 为-1,u-D 改为 1.0e6,然后点击图 1-1-19 中的 Assign 按钮,施 加左边界。 然后, 点击 Conditions 对话框中的 Finish 按钮。 这样, 压力为 P=1.0 6 ×10 Pa 的边界施加完毕。
1.1 非达西渗流 Forchheimer 模型 1.1.1 问题描述 渗流问题是水利工程常见的一种问题,实际工程中,人们常常要对渗流进行 非达西渗流是地下水及石油工程中常见的问题,当孔隙中流体的流动过快 (流动进入湍流状态)或者过慢(受粘性影响)都会出现非达西渗流的情况。以 下算例就是基于 Forchheimer(1901)提出的非达西渗流模型,求解边长为 100m 的正方形区域内非达西渗流的情况。 区域左端作用压力 P=1.0×106Pa, 其余三边界封闭。 其中区域渗透率为 K=1.0 -15 2 ×10 m ;动力黏度μ=0.001Pa.s;渗透容度 K=3.0×10-11Pa;孔隙度φ=0.15; 非线性流动系数β=1.0×108;初始压力为 P=1.0×105Pa。求区域的压力分布 P, 速度 V 分布。 问题分析: 非线性达西流动的偏微分方程的表述如下:
图 1-1-32 右边界中点压力随时间的演进图
(5)将非达西系数β设置为 0.0,问题退化为达西渗流,与上述非达西渗流的计算 结果对比如下图 1-1-33~图 1-1-40 所示,通过比较可以发现,在采用非达西流计 算区域内渗透场演进,在初期渗流速度较达西流情况有明显的减小效果,这跟实 际情况中在到达一定的达西流启动雷诺数之前为非达西流状态是能够符合起来 的。
图 1-1-21
定义单元尺寸大小
此时弹出“Progress in meshing”进度条,网格划分完毕后弹出消息框显示总 的单元数和节点数,如图 1-1-22 所示。
图 1-1-22
消息框
(2).把前处理数据转化为 SciFEA 所需计算格式。点击【Files】-【Save】 菜单, 保存为 GID 文件。 点击菜单 【Calculate】 【Calculate】 , 弹出 “process info”消息框,如图 1-1-23 所示,点击“OK”按钮,然后退出 Gid,至此前 处理工作结束。
图 1-1-8 点击 Line 菜单项
如下图 1-1-9 所示,在 GID 命令栏中先输入 0,0,按 ENTER 键,然后然后依次输 入(100,0)、(100,100)、(0,100)、(0,0)。
图 1-1-9 输入 0,0
输入最后一个 0,0 时, 系统提示是否连接原先节点 , 如下图 1-1-10 。 选 择”JOIN”,形成如图 1-1-11 所示封闭线段。
1.1.2.b 设置材料参数和边界条件
(1)选择“前处理”—>“材料参数”按钮,如图 1-1-2 所示。或者单击工具 条中的 按钮弹出如图 1-1-3 所示材料参数数据输入表格。
图 1-1-2
选择材料参数输入
(2) 按照问题描述中的参数依次填入材料参数数据表格。填写完成后如图 1-1-3 所示。
图 1-1-3
图 1-1-15
图 1-1-16 Conditions
②定义初始条件。参照图图 1-1-16,选择 Surface-D1nda,在 p-D1 栏中填入 1.0e5, 作为计算区域的初始压力 0.1MPa,即一个大气压。得到的赋值结果如图 1-1-18 所示
图 1-1-17 材料属性结果
图 1-1-18 赋初值结果
图 1-1-19
图 1-1-20
点击如图 1-1-19 所示的 Draw->Colors,显示的边界条件的施加情况如图 1-1-20 所示。
1.1.2.d 划分网格和导出数据
(1).划分网格。点击菜单【Meshing】-【Generate】 ,弹出“Enter Value window”
对话框,要求定义单元尺寸大小,输入“2” ,点击“ok”按钮,如图 1-1-21 所 示。
填写完成材料数据输入
(3) 选择“前处理”—>“边界条件”按钮,如图 1-1-4 所示。或者单击工具
条中的
按钮弹如图 1-1-5 所示数据输入表格。
图 1-1-4
选择边界条件输入
(4)按照问题描述中的参数依次填入边界条件数据表格。填写完成后如图 1-1-5 所示。
图 1-1-5
填写完成边界条件输入
图 1-1-40 非达西流右边界中点速度随时间的演进图
合集下载

每天一点土力学(17):二维渗流及应用

每天一点土力学(17):二维渗流及应用

每天一点土力学(17):二维渗流及应用拉普拉斯对于边界条件简单的一维渗流问题,可直接利用达西定律进行渗流计算。

但工程中遇到的渗流问题,常常属于边界条件复杂一些的二维或三维渗流问题。

例如闸坝下透水地基的渗流以及土坝坝身的渗流等(图a),其流线都是弯曲的,不能再视为一维渗流。

为了求解和评价渗流在地基或坝体中是否造成有害的影响,需要知道整个渗流场中各处的测管水头、渗流坡降和渗流速度。

著名的拉普拉斯(Laplace)方程。

它与水力学中的平面势流的拉普拉斯方程一样,该方程描述了渗流场内部的测压管水头h 的分布,是平面稳定渗流的基本方程式。

通过求解一定边界条件下的拉普拉斯方程,即可求得该条件下的渗流场。

流网的特征与绘制上述拉普拉斯方程表明,渗流场内任一点水头是其坐标的函数,知道了水头分布,即可确定渗流场的其他特征。

求解拉普拉斯方程一般有4类方法,即数学解析法、数值解法、电模拟法、图解法。

其中图解法简便、快速,在工程中实用性强,因此,这里简要介绍图解法。

所谓图解法即用绘制流网的方法求解拉普拉斯方程的近似解。

1)流网的特征流网是由流线和等势线所组成的曲线正交网格。

在稳定渗流场中,流线表示水质点的流动路线,流线上任一点的切线方向就是流速矢量的方向。

等势线是渗流场中势能或水头的等值线。

流网的绘制如图所示为板桩墙围堰的流网图。

图中实线为流线,虚线为等势线。

对于各向同性渗流介质,由水力学可知,流网具有下列特征:①流线与等势线互相正交;②流线与等势线构成的各个网格的长宽比为常数,当长宽比为1时,网格为曲线正方形,这也是最常见的一种流网;③相邻等势线之间的水头损失相等;④各个流槽的渗流量相等。

由这些特征可进一步知道,流网中等势线越密的部位,水力梯度越大,流线越密的部位流速越大。

流网的工程应用根据流网,就可以直观地获得渗流特性的总体轮廓,并可定量求得渗流场中各点的测管水头、水力坡降、渗流速度和渗流量。

1)测管水头根据流网的特征可知,任意两相邻等势线间的势能差相等,即水头损失相等,从而相邻两条等势线之间的水头损失Δh 为:2)孔隙水压力渗流场中各点的孔隙水压力,等于该点以上测压管中的水柱高度h u乘以水的容重γ w,故 a点的孔隙水压力为:3)水力坡降流网中任意网格的平均水力坡降,Δl 为该网格处流线的平均长度,可从图中量出。

《非达西渗流实验研究及数学描述》

《非达西渗流实验研究及数学描述》

《非达西渗流实验研究及数学描述》篇一一、引言在多孔介质中,流体流动的研究对于众多领域,如地质工程、石油工程、环境科学等,都具有重要的意义。

传统的达西渗流模型在许多情况下都能很好地描述流体的流动行为,然而,在某些复杂条件下,如高流速、高压力梯度或复杂多孔介质中,达西渗流模型可能无法准确描述流体的实际行为。

为了更深入地理解这些非传统渗流现象,非达西渗流实验研究及其数学描述变得尤为重要。

本文将通过实验研究及数学描述的方式,探讨非达西渗流的特点及规律。

二、非达西渗流实验研究1. 实验装置与材料实验采用多孔介质作为实验材料,并搭建了一套能够进行高精度压力测量和流体流速测量的实验装置。

多孔介质选用不同类型的介质样品,以观察不同介质对渗流行为的影响。

2. 实验过程在实验过程中,我们通过改变流体的流速、压力梯度以及多孔介质的类型等条件,观察并记录了流体的流动行为。

同时,我们还对实验数据进行了实时监测和记录,以便后续的数据分析和处理。

3. 实验结果通过实验观察和记录,我们发现非达西渗流现象在许多情况下都存在。

在特定的条件下,如高流速或高压力梯度时,流体在多孔介质中的流动行为会明显偏离达西渗流模型。

同时,不同类型的多孔介质对渗流行为的影响也不同。

三、数学描述为了更准确地描述非达西渗流现象,我们引入了非达西渗流模型。

该模型基于达西渗流模型的基础上进行改进,考虑了更多的影响因素,如流速、压力梯度以及多孔介质的特性等。

通过数学推导和计算,我们可以得到非达西渗流的数学表达式。

该表达式可以更准确地描述非达西渗流现象的规律和特点。

四、结论通过对非达西渗流实验的研究和数学描述,我们得到了以下结论:1. 非达西渗流现象在许多情况下都存在,尤其是在高流速或高压力梯度等复杂条件下。

2. 不同类型的多孔介质对渗流行为的影响也不同,因此在进行渗流研究时需要考虑多孔介质的特性。

3. 非达西渗流模型能够更准确地描述非达西渗流现象的规律和特点,为相关领域的研究提供了重要的理论依据。

《非达西渗流实验研究及数学描述》范文

《非达西渗流实验研究及数学描述》范文

《非达西渗流实验研究及数学描述》篇一一、引言在地质学、水文学、地下水动力学等多个领域中,渗流实验是非常重要的研究方法之一。

在众多的渗流理论中,达西渗流理论是最为经典和广泛应用的。

然而,在实际的地下环境中,由于多孔介质的复杂性和非均质性,非达西渗流现象时常发生。

非达西渗流研究能够更好地揭示实际地质条件下的地下水流动规律,因此对其实验研究和数学描述具有重要意义。

本文将就非达西渗流实验研究及数学描述进行探讨。

二、非达西渗流实验研究1. 实验原理非达西渗流实验基于地下水的流动特性,在特定的多孔介质条件下,观察和分析其渗流过程和规律。

与达西渗流相比,非达西渗流更多地考虑了多孔介质的复杂性和非均质性,以及不同流速和压力条件下的影响。

2. 实验方法(1)实验材料准备:选择具有代表性的多孔介质材料,如砂土、粘土等。

同时,需要准备实验所需的仪器设备,如压力传感器、流量计、渗透仪等。

(2)实验过程:首先,将多孔介质材料装入渗透仪中,然后施加一定的压力,观察和记录渗透过程中的流量变化和压力变化情况。

在实验过程中,需要注意控制流速和压力的变化范围,以及记录实验数据的时间间隔等。

(3)数据分析:通过分析实验数据,可以得出非达西渗流的特征参数,如渗透系数、流量指数等。

同时,还可以根据实验数据绘制出渗流曲线和压力分布曲线等,以便更好地分析和理解非达西渗流的规律。

三、非达西渗流的数学描述1. 非达西渗流模型非达西渗流模型是一种基于地下水的实际流动特性而建立的数学模型。

该模型能够描述多孔介质中地下水的非线性流动规律,包括流速、压力、渗透系数等因素的影响。

目前,常用的非达西渗流模型包括幂律模型、指数模型等。

2. 数学表达式及解析(1)幂律模型:在幂律模型中,流量与压力梯度之间存在幂律关系。

其数学表达式为:q = -K|▽P|n ,其中q为流量,K为渗透系数,▽P为压力梯度,n为幂律指数。

该模型能够较好地描述低速和高速条件下的非达西渗流现象。

低渗非达西渗流特征及影响因素

低渗非达西渗流特征及影响因素

《高等渗流力学》读书报告----低渗非达西渗流特征及影响因素姓名: 张恒学号:2010050031专业:石油与天然气工程教师:鲁洪江(教授)低渗非达西渗流特征及影响因素1 选题依据及研究现状1.1选题依据随着中国石油工业的发展,低渗透油藏在开发中所占的比例越来越大。

低渗透油藏是我国今后乃至相当长一段时间内增储上产的主要资源基础。

要合理高效地开发这些低渗透油藏,就需要充分合理的认识低渗透油藏本身所具有的特殊规律及其特性参数,并准确地描述低渗透油藏的渗流规律.1.2研究现状国内很多研究人员从实验方面发现了低渗透油藏的启动压力和非线性渗流规律的存在,从理论方面提出了描述启动压力和非线性渗流的模型[1]。

但是,非线性渗流和启动压力梯度的存在并没有得到国内外学术界的普遍认可。

反对者的意见是,引起低渗透油藏非线性达西流和启动压力的原因均为理论推测,而无充分的微观实验科学依据;在流速很低的情况下,受测量手段和如蒸发等现象的影响,对流速和压力的测量误差很大[2]1.3 主要的参考文献[1] 王正波,岳湘安等.影响低渗透油藏低速非线性渗流的实验研究[J].矿物学报,2008,28(1),48-54.[2]王慧明,王恩志等.低渗透岩体饱和渗流研究进展[J].水科学进展, 2003,14(2): 245[3]辛莹娟.低渗透非达西渗流研究[J].西部探矿工程。

2010(10):115-117[4]中国“八五”科技成果.低渗透油层多相渗流机理[M].北京:科学出版社,1996[5]闫庆来,何秋轩,任晓娟,等.低渗透油层中单相液体渗流特征研究[J].西安石油学院学报,1990,5(6):1-6.[6]吴景春,袁满,张继成,等.大庆东部低渗透油藏单相流体低速非达西渗流特征[J].大庆石油学院学报,1999,23(2):82-84[7]阮敏,何秋轩.低渗透多孔介质中新型渗流模型[J].石油勘探与开发,1996[8]程时清,徐论勋,张德超.低速非达西渗流试井典型曲线拟合法[J].石油勘探与开发,1996.[9]宋付权,刘慈群.低渗透多孔介质中新型渗流模型[J].石油勘探与开发,1996.[10]吴景春,袁满,张继成,等.大庆东部低渗透油藏单相流体低速非达西渗流特征[J].大庆石油学院学报,1999.[11]李道品,等低渗透砂岩油田开发[M].北京:石油工业出版社,1997.[12]诺曼,R莫罗.石油开采中的界面现象[M].鄢捷年等译.北京:石油工业出版1992.23~85.[13]邓英尔,闫庆来,马宝歧.表面分子对低渗多孔介质中液体渗流特征的影响[A].渗流力学进展[C].北京:石油工业出版社,1996.9.[14]阮敏.低渗透非达西流临界雷诺数实验研究[J].西安石油学院学报,19992选题研究内容及拟解决的问题达西定律中压力损失完全由粘滞阻力决定,这符合多孔介质比面大这个特点的.而在低渗透岩石中,流体在流动过程中受到岩石孔壁、粘土矿物遇水膨胀以及岩石颗粒的运移等一系列因素的影响而造成附加压力损失,所以流体在低渗透砂岩中的渗流规律不满足达西定律达西定律是渗流的基本规律,但是在低渗透油藏中,渗流表现出对达西定律的偏离,这就使我们有必要对非达西渗流进行深入的研究,从低渗透非达西渗流特征、低渗透非达西渗流模型,非达西渗流过程等几个方面的研究进展进行了总结.为从事相关工作的研究人员提供参考[3]3 方法及路线3.1 低渗透非达西渗流特征同中高渗透率油层相比,低渗透油层具有以下几个特点:低渗透油层一般连续性差、采收率与井网密度关系特别密切;低渗透油层存在“启动生产压差现象”,渗流阻力和压力消耗特别大;低渗透油层见水后,采液和采油指数急剧下降,对油田稳产造成急剧影响;低渗透油田一般裂缝都较发育,注入水沿裂缝窜进十分严重[4].室内实验结果表明,流体在低渗透储层内渗流时,存在非线性段[5,6]压力梯度超过某一定值后,渗流曲线变为直线,见图1.由图1可知,流体通过低渗岩心的渗流特征显示出弹-塑性.低渗透非达西渗流的特征可概括为以下两点.(1)在较宽的渗流速度域内,渗流过程由2个连续过渡而特性各异的渗流曲线段组成,即:低渗流速度下的凹型线性渗流曲线段;较高渗流速度下的直线段;(2)当压力梯度在比较低的范围时,渗流速度是上凹型非线性曲线。

低渗透油藏非达西渗流规律研究.

低渗透油藏非达西渗流规律研究.

dpwD dPD
d t] D
刁rD
了 6︑
. ︑
分
︐ J
a P勺
无限大油藏动边界条件 :
一刁
一 r
一 劝
二一ADe-'
了了
7
︑ ︐
︑
J
rD=R(cp)
封闭油藏外边界条件:
|
aPrD一r。一*(,D)--ADC',I<R ‘: ·,<Re
‘ |
PDIrp > R('o)=0
(8)
十
aPrD一*(ip)一二。一“
本上也是油层实际渗流的速度范围,可见在低渗透油藏级别上可能普遍存在非达西渗流。大 量的实验表明,低渗条件下的油气渗流有如图3所示的特征,图 1,2 中实线是实验所得的
渗流曲线,压力梯度大的部分呈直线段,表现为达西流;压力梯度小的部分表现为非达西 流。为便于工程计算,把上面的直线段延长,用延长的直线段代替渗流曲线,2就是启动压 力梯度,可用公式 (1)描述图3所示非达西渗流规律:
下面以无限大油藏为例说明一下差分方法。 对式 (10)进行数值离散,得:
,‘一,一‘,·箭 ,,、,·,‘一 △一“Der,-, 'OD er22'eO2,t
(15)
(i二0,…,1一1;,=0,’·’,N)。
由式 (11),(12),( 13)分别可得 :
Pi.o=0
(16)
,。,一“+瓮,Pi,十,2,一,△·“+“一‘,一鄂1,j一 (17)
定压外边界条件:
r|
leseseseses之 aPrDl rp=R(' D)一“一‘<R‘!·,<Re
| ..!
PDIrp > R(no)=0
(9)
| PDIR O D)=R .=0

《非达西渗流实验研究及数学描述》

《非达西渗流实验研究及数学描述》

《非达西渗流实验研究及数学描述》篇一一、引言渗流现象在众多领域中具有广泛的应用,如地下水流动、油藏开发等。

达西定律作为描述渗流现象的基本定律,在许多情况下都适用。

然而,在实际的复杂多孔介质中,由于介质结构、压力分布以及流动性质的不同,有时会表现出非达西渗流的特点。

因此,对非达西渗流现象进行实验研究及数学描述显得尤为重要。

本文将介绍非达西渗流实验研究的过程,以及相关数学模型的建立和求解。

二、实验研究1. 实验装置与材料实验采用多孔介质作为研究对象,如砂土、碎石等。

实验装置包括储液器、压力计、测量仪等。

其中,储液器用于储存流体,压力计用于测量不同位置的压力变化,测量仪用于记录实验过程中的数据。

2. 实验步骤(1)将多孔介质装入实验装置中,确保介质紧密且均匀分布。

(2)向储液器中注入流体,并观察多孔介质中的渗流现象。

(3)记录不同位置的压力变化及流体流量数据。

(4)改变流体的性质或施加不同的压力条件,重复上述步骤进行多次实验。

3. 实验结果与分析通过实验数据可以观察到,在一定的压力条件下,多孔介质中的渗流现象表现出非达西特点。

具体表现为流速与压力梯度之间的关系偏离了达西定律的线性关系。

此外,流体的性质、多孔介质的孔隙结构等因素也会对渗流现象产生影响。

三、数学描述为了更好地描述非达西渗流现象,我们引入了非达西渗流模型。

该模型基于达西定律的基本思想,但在描述流速与压力梯度之间的关系时采用了更为复杂的数学表达式。

具体而言,非达西模型将渗流过程分为三个阶段:低速流动阶段、过渡阶段和高速流动阶段。

在每个阶段中,流速与压力梯度之间的关系都有所不同,需要采用不同的数学表达式进行描述。

对于非达西渗流模型的求解,我们可以采用数值计算的方法。

具体而言,可以通过离散化空间和时间,将渗流过程转化为一系列的代数方程组进行求解。

在求解过程中,需要考虑到多孔介质的孔隙结构、流体性质以及边界条件等因素的影响。

通过求解非达西渗流模型,我们可以得到不同位置的压力分布、流体流量等关键参数,从而更好地理解非达西渗流现象的特性和规律。

《非达西渗流实验研究及数学描述》范文

《非达西渗流实验研究及数学描述》篇一一、引言渗流现象在地质学、石油工程、地下水动力学等多个领域具有重要应用。

达西定律作为经典渗流理论的核心,长期以来在描述低速渗流过程中发挥着主导作用。

然而,在高速流动或复杂流体系统内,达西定律可能不再适用,这就需要对非达西渗流进行深入的研究。

本文通过非达西渗流实验研究,探究了高速或复杂条件下流体在多孔介质中的运动规律,并对该过程的数学描述进行了详尽的分析和推导。

二、非达西渗流实验研究1. 实验原理非达西渗流实验的原理主要基于对多孔介质中流体流动的观测和测量。

实验中,我们通过改变流体的速度、压力、介质特性等因素,观察并记录流体在多孔介质中的运动轨迹和流动速度,从而分析其是否符合达西定律。

2. 实验过程(1) 实验设备与材料准备:多孔介质样本、恒速泵、压力传感器、流量计、计算机等。

(2) 实验操作:将多孔介质样本置于实验装置中,通过恒速泵控制流体速度,利用压力传感器和流量计测量并记录数据。

(3) 数据处理与分析:将收集到的数据整理成表格或图形,分析其变化规律,并与达西定律进行对比。

3. 实验结果与讨论通过对非达西渗流实验的观测和数据分析,我们发现:在高速流动或复杂流体系统中,流体在多孔介质中的运动不再遵循达西定律。

非达西渗流具有更复杂的流动形态和速度分布,同时伴随着更多的非线性效应和波动性。

此外,不同因素(如介质结构、流体性质等)对非达西渗流的影响也不同。

三、非达西渗流的数学描述为了更好地描述非达西渗流的特性,我们引入了非线性渗流模型。

该模型基于流体在多孔介质中的实际运动规律,考虑了多种因素的影响,如介质结构、流体性质等。

通过对模型的推导和分析,我们可以得到以下数学表达式:非达西渗流的速度与压力之间的关系可表示为:v = f(p, k, ε, μ),其中v为流体速度,p为压力,k为介质渗透率,ε为孔隙度,μ为流体动力粘度等参数。

根据实际情况,该表达式中的f可根据具体的模型和参数进行调整和优化。

直裂缝井两相非达西椭圆渗流特征线解,差分解及开发

直裂缝井两相非达西椭圆渗流特征线解,差分解及开发
直裂缝井是一种特殊形式的井,在水力良好的场合,可以实现高产量、高效率的可采储层为目标的开发目的。

但是直裂缝井非达西椭圆渗流是较为复杂的,需要进行研究分析以进行更好的井开发。

针对直裂缝井的特殊性质,研究一次渗透率发生变化的非达西椭圆渗流问题。

首先,假设直裂缝井为一椭圆截面,下面根据直裂缝井的特性进行研究。

通过推导和分析,可得到裂缝与截面的夹角为α,截面到上下壁的距离分别为2a和2b,宽度为2c,单位宽度的流体注入压力为pw ,地层压力为p0,流体的动粘度为μ, 不同的参数条件下得到了若干个特征线。

对于直裂缝井两相非达西椭圆渗流的差分解,主要是利用数值计算对其进行仿真来研究。

对于直裂缝井两相渗流的数值计算,采用了单元格有限差分法进行, 数值计算程序采用了Pascal语言开发,数值计算结果可以很好的描述直裂缝井两相非达西椭圆渗流的规律和特性。

为了更好地开发直裂缝井两相非达西椭圆渗流,需要以可靠的数学模型为基础.开发直裂缝井的方法主要包括:
1.利用直裂缝井的特性进行井的设计,优化产量和开采效率。

2.根据差分解所得到的数学模型,在实际计算时进行调整,获得更加真实的数据。

3.考虑直裂缝井对地下水环境的影响,建立监测网站,对井的开采过程进行监测,并做好环境保护工作。

结论
直裂缝井两相非达西椭圆渗流的特征线解,差分解及开发是一项重要的研究,可以为直裂缝井的开发提供理论支持和实践指导。

研究结果可以为井的设计、生产和管理提供可靠的数据,为国家石油钻采行业推进创新提供了新的方向。

《非达西渗流实验研究及数学描述》

《非达西渗流实验研究及数学描述》篇一一、引言渗流现象在多孔介质中广泛存在,如地下水流动、油藏开发等。

传统的达西渗流理论在低流速下具有很高的适用性,然而,随着流速的增加或介质性质的复杂化,非达西渗流现象逐渐显现。

本文将就非达西渗流实验研究及其数学描述进行探讨,旨在为相关领域的研究提供理论依据和实验支持。

二、非达西渗流现象及特点非达西渗流是指在一定条件下,多孔介质中流体的流动不遵循达西定律的渗流现象。

其主要特点包括流速分布的非线性变化、介质内部的非均匀性以及流体与介质之间的相互作用等。

非达西渗流现象在许多工程领域中具有重要影响,如油藏开发、地热利用等。

三、非达西渗流实验研究1. 实验材料与装置本实验主要使用石英砂、活性炭等作为多孔介质材料,利用蠕动泵、压力传感器等设备进行实验装置的搭建。

通过改变流速、压力等参数,观察并记录非达西渗流现象的变化。

2. 实验方法与步骤首先,根据实验目的选择合适的介质材料和装置。

然后,在恒定压力下,逐步提高流速,观察并记录不同流速下的渗流现象。

同时,通过改变介质性质(如孔隙度、渗透率等),分析其对非达西渗流的影响。

最后,根据实验数据绘制相应的图表,分析并总结实验结果。

四、数学描述及模型建立针对非达西渗流现象,我们提出了一种数学描述方法及相应的模型。

该模型基于流体力学原理和实验数据,能够较好地描述非达西渗流现象的规律性。

具体而言,该模型包括以下部分:1. 模型假设与基本原理模型假设多孔介质中的流体流动遵循一定的物理规律,如质量守恒定律、牛顿第二定律等。

在此基础上,通过引入适当的边界条件和初始条件,建立数学模型。

2. 模型建立与求解根据实验数据和理论分析,我们建立了以压力和流速为变量的微分方程组,通过求解该方程组得到不同条件下的非达西渗流规律。

具体而言,我们采用了有限差分法对微分方程进行数值求解,得到了不同条件下的流体速度分布和压力分布。

3. 模型验证与适用性分析为验证模型的正确性和适用性,我们进行了大量的实验数据验证。

非达西流动系数

非达西流动系数一、引言非达西流动系数(Non-Darcy flow coefficient)是描述非达西流动时,流体在孔隙介质中的阻力大小的一个参数。

它是指在非达西流动条件下,单位长度内液体通过孔隙介质所需的压差与单位长度内液体速度平方之比。

非达西流动系数的研究对于理解地下水文学、油藏工程、地热开采等领域具有重要意义。

二、什么是非达西流动在孔隙介质中,当液体速度较小时,液体分子与孔隙壁面之间的摩擦力可以忽略不计,此时称为达西流动(Darcy flow)。

但当液体速度较大时,由于摩擦力增大,导致液体通过孔隙介质的阻力也增大,在这种情况下就称为非达西流动。

三、什么是非达西流动系数在非达西流动条件下,单位长度内液体通过孔隙介质所需的压差与单位长度内液体速度平方之比定义为非达西流动系数。

通常用符号α表示。

α值越小,说明阻力越大。

四、影响因素1. 孔隙结构:孔隙结构的大小、形状、连通性等因素都会影响非达西流动系数的大小。

2. 流体性质:流体的黏度、密度等性质也会影响非达西流动系数。

3. 流速:非达西流动系数随着流速的增加而增加。

五、测量方法1. 实验室试验法:通过实验室试验,测量液体在不同孔隙介质中的流速和压差,计算出非达西流动系数。

2. 数值模拟法:利用计算机模拟液体在孔隙介质中的运动过程,计算出非达西流动系数。

六、应用领域1. 地下水文学:非达西流动系数是描述地下水运移过程中阻力大小的重要参数,对于地下水资源开发和管理具有重要意义。

2. 油藏工程:油藏中存在较多孔隙介质,非达西流动系数对于油藏开采和生产具有重要影响。

3. 地热开采:地热能开采需要通过岩石孔隙介质传递热能,了解非达西流动系数对于优化地热能开采方案具有重要作用。

七、总结非达西流动系数是描述非达西流动时,流体在孔隙介质中的阻力大小的一个重要参数。

它受到孔隙结构、流体性质和流速等因素的影响。

测量非达西流动系数的方法包括实验室试验法和数值模拟法。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
相关文档
最新文档