分形渗流分析研究

wyt2010

分形渗流是指分形孔隙介质和分形裂缝介质中流体的流动，关于分形渗流的研究，近年来有了很大发展。详见合肥辰工科技有限公司网站。
上世纪90年代，曾有人提出过一种分形渗流的数学模型，导出了分形油藏的压力扩散方程。但是他们的模型中，有的分形参数物理意义不太明确，使用不便，受到了质疑。
近10年国内外在研究分形多孔介质渗透率方面做了大量工作，这些文章的应用背景大多偏向于塑料成型的模拟，基本不涉及分形油藏。研究的参数重点是针对一定的单元体的渗透率。对孔隙度研究较少，大多是给出一些半经验性质的表达式。
针对分形油藏，对分形孔隙介质和分形裂缝介质给出了渗流速度、渗透率和孔隙度三个基本公式，建立起分形油藏中渗流的压力扩散方程，并求得了相应的解析解，绘制出典型的压力曲线。
着重对分形油藏的一些特性参数作进一步的分析研究，讨论了分形渗流与传统达西渗流的某种相关性，建立起分形双重介质（裂缝和孔隙）中渗流的数学模型，对致密性油藏渗流的分析研究提供一种分析计算方法。
第一节分型油藏的数学描述
为区别起见，将描述裂缝和孔隙介质的参数分别用下标1和2表示。众所周知，将Navier-Stokes方程用于求解单位深度裂缝和圆截面毛管中流动，可求得体积流量Q与压力梯度之间的关系式分别为：
                        
其中和分别表示裂缝宽度和毛管直径，是裂缝的单位深度，是流体粘度，是压力梯度。考虑流动路径迂曲的分形特性，有：               
对于径向流，上式中Ｌ应改为ｒ，于是式和可分别写成：               
其中称为分形曲线的迂曲分维，L是沿流动方向的外观长度，是弯曲流线的有效长度。
考虑通过垂直于流动方向单位截面积有大量的毛管，其最大直径和最小直径分别用和表示。管径尺寸分布具有分形特性。即毛管累积数与管径之间遵从以下标度关系：
               
其中是参考长度，称为管径分维。因而通过单位截面毛管的累积数的微分和毛管总数分别为：               
               
对于裂缝分形介质，首先需对通过截面的裂缝进行量化。将裂缝沿流动方向的尺度称为裂缝长度，将垂直于流动方向的尺度称为裂缝深度（无论是铅垂或倾斜方向）。裂缝深度ｄ有大有小，设单位截面上，同一宽度，不同深度的裂缝有ｎ条，将它们的深度相加，给出：
               
其中为单位深度，是上宽度为的诸裂缝深度之和折算为单位深度的条数，它是正数，但一般不是整数。类似地有：        
               
其中称为缝宽分维。这样量化的涵义是将裂缝介质化为二维分形体。和分别为通过单位截面上的最大和最小缝宽。
由Q对从到积分，可给出渗流速度V（或比流量q）从而可定义出分形渗透率。裂缝或毛管截面积乘以对积分，可给出孔隙度。这样，可以得出分形油藏渗流的三个基本公式（以径向流为例）写成：
分形渗流速度：
        
其中：，j=1，2分别对应于裂缝和孔隙介质。和分别为渗透率常数和孔隙度常数。它们只与介质本身的结构特性有关。
        
其中，常数m对平行流、平面径向流和球心向心流分别为0、1和2。若考虑表皮S和储集常数，则无量纲的压力扩散方程和内边界条件可分别写成：
                        
铅直井的压力分析的典型曲线如图1所示。虽然分形渗流与经典达西渗流的理论基础不同，经典渗流是建立在达西定律的基础上，而分形渗流是基于介质的分形理论和流体力学的Navier-Stokes方程，但是其数学描述普遍存在明显的相关性。我们可以看到对的极限情形，分形渗流的扩散方程、内边界条件以及其它各式在形式上都退化为经典渗流的情形。