用于无网格油藏数值模拟的裂缝网络点云生成方法

1.本发明涉及油藏数值模拟领域。更具体地说，本发明涉及一种用于无网格油藏数值模拟方法的复杂缝网点云生成方法。


背景技术：

2.无网格法一般仅需要点云来离散计算域，而点云的生成技术理论上来说可以直接借用网格生成技术中生成网格节点的部分，无需在点云基础上去确定网格结构中所要求的丰富拓扑结构，从而在计算域离散和自适应分析方面具有更加灵活的优势，这是学者们发展无网格法的最原始的也是最重要的动机之一。当把无网格法应用到油藏数值模拟领域时，具有复杂几何特征的裂缝性油藏或缝洞型油藏是无网格油藏数值模拟方法的良好应用场景。然而，对于计算域中的复杂缝网，由裂缝交叉点导致的裂缝小段/尖端，以及几乎接触的裂缝之间的狭长基质区域等导致了生成匹配性网格的困难。在对带有复杂裂缝网络的计算域进行匹配性网格剖分时，匹配性网格的生成需要在裂缝交点以及裂缝之间的狭窄基质区域处进行局部网格加密。对附图1(a)所示带有复杂缝网的矩形边界油藏计算域，在利用delaunay三角剖分技术获取附图1(b)中该矩形区域内复杂缝网的匹配性网格时，可以看到有很多局部的高密度网格，对于这样一个小模型，计算域的网格节点总数就已达到了18819。嵌入式离散裂缝模型(edfm)可以规避这个困难，因为其仅需对计算计算域基质系统生成网格离散，而将裂缝网络通过几何关系分配给相应的基质网格，可以称之为非匹配性网格。类似得，当我们用点云来离散油藏计算域时，根据对裂缝网络处理方式的不同也可以将计算域点云分为匹配性的和非匹配性的。其中，匹配性点云是指裂缝节点的同一位置也会相应得有一个基质节点，即裂缝节点处有基质和裂缝两个系统，以两相流模拟为例，则该节点处实际上有四个自由度，分别是裂缝系统(节点)压力、裂缝(节点)含水饱和度、基质系统(节点)压力、基质(节点)含水饱和度。所以基于匹配性点云的复杂缝网流动中基质系统与裂缝系统之间的传质则通过裂缝节点和与裂缝节点相同位置的基质节点之间的传质来实现。而非匹配性点云则是通过判断裂缝节点应该在哪一个基质节点的节点控制区域，而建立裂缝节点与该基质节点之间的连接来表征基质系统和裂缝系统之间的传质。


技术实现要素：

3.本发明的目的是提供一种用于无网格油藏数值模拟方法的复杂缝网匹配性和非匹配性点云生成算法。
4.为了实现根据本发明的这些目的和其它优点，提供了一种用于无网格油藏数值模拟的裂缝网络点云生成方法，包括以下步骤：
5.s1、分别建立两条裂缝的参数方程，并计算裂缝交点在所述参数方程中的参数向量；
6.s2、生成裂缝节点时，根据s1中得到的裂缝交点的参数向量计算裂缝的每段裂缝段的真实长度；
7.s3、当对裂缝网络计算域进行非匹配性点云生成时：
8.s3.1、忽略裂缝网络形态生成计算域的基质节点的生成；
9.s3.2、建立基质节点与裂缝节点之间的相邻关系；
10.s4、当对裂缝网络计算域进行匹配性点云生成时：
11.s4.1、忽略裂缝网络形态生成计算域的基质节点的生成；
12.s4.2、设定距离阈值，去掉与裂缝节点同位置的基质节点之间的距离小于所述距离阈值的基质节点。
13.优选的是，所述的用于无网格油藏数值模拟的裂缝网络点云生成方法中，s3.2中建立基质节点与裂缝节点之间的相邻关系时，设定基质节点与裂缝节点的控制区域均为边长为li的正方形区域，裂缝节点f的坐标为(xf,yf)，基质节点i的坐标为(xi,yi)，基质节点与裂缝节点的节点控制体积均记作vi：
14.基质节点与裂缝节点建立相邻关系的准则为：
[0015][0016]
优选的是，所述的用于无网格油藏数值模拟的裂缝网络点云生成方法中，s4.2中与裂缝节点同位置的基质节点的坐标为(x1，y1)，去掉的基质节点的坐标为(x2，y2)，且：
[0017][0018]
其中，dis
tv
为距离阈值。
[0019]
优选的是，所述的用于无网格油藏数值模拟的裂缝网络点云生成方法中，s1中计算裂缝交点在所述参数方程中的参数向量的方法如下：
[0020]
s1.1、设定其中一条裂缝的两端分别为a(a1,b1)和b(a2,b2)，另一条裂缝的两端分别为c(c1,d1)和d(c2,d2)，则两条裂缝的参数方程为
[0021]
(x,y)＝(a1,b1) ξ1(a
2-a1,b
2-b1)
[0022]
(x,y)＝(c1,d1) ξ2(c
2-c1,d
2-d1)
ꢀꢀꢀ
(3)
[0023]
其中，ξ1和ξ2分别为裂缝交点在两条裂缝的参数方程上的参数值，0≤ξ1≤1，0≤ξ2≤1；
[0024]
由上式(3)得到关于ξ1和ξ2的线性方程组：
[0025][0026]
优选的是，所述的用于无网格油藏数值模拟的裂缝网络点云生成方法中，s2中根据s1中得到的裂缝交点的参数向量计算裂缝的每段裂缝段的真实长度的方法如下：
[0027]
s2.1、定义裂缝的总长是lf，其参数向量为ξf，相邻的两个参数值分别是ξ
fi
和ξ
fi 1
，ξ
fi
和ξ
fi 1
对应的裂缝f部分被剖分的段数计算为：
[0028][0029]
其中，floor(x)函数是指取整函数，即当n≤x＜n 1时，n为整数，floor(x)＝n；
[0030]
s2.2、计算裂缝的每段裂缝段的真实长度l
real
:
[0031][0032]
本发明的用于无网格油藏数值模拟的裂缝网络点云生成方法，可以用于复杂缝网匹配性和非匹配性点云生成。
[0033]
本发明的其它优点、目标和特征将部分通过下面的说明体现，部分还将通过对本发明的研究和实践而为本领域的技术人员所理解。
附图说明
[0034]
图1为本发明所述的矩形计算域及三角网格剖分的示意图；
[0035]
图2为本发明一实施例中图1的示意图；
[0036]
图3为本发明一实施例中图1的另一示意图；
[0037]
图4为本发明一实施例中图1的另一示意图；
[0038]
图5为本发明所述的不规则边界计算域及三角网格剖分的示意图；
[0039]
图6为本发明所述的不规则边界复杂缝网计算域的匹配性点云生成的示意图；
[0040]
图7为本发明所述的不规则边界复杂缝网计算域的非匹配性点云生成的示意图。
具体实施方式
[0041]
下面结合附图对本发明做进一步的详细说明，以令本领域技术人员参照说明书文字能够据以实施。
[0042]
需要说明的是，在本发明的描述中，术语“横向”、“纵向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，并不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
[0043]
《实施例1》
[0044]
本发明的实施例提供一种用于无网格油藏数值模拟的裂缝网络点云生成方法，包括以下步骤：
[0045]
s1、分别建立两条裂缝的参数方程，并计算裂缝交点在所述参数方程中的参数向量；
[0046]
计算裂缝交点在所述参数方程中的参数向量的方法如下：
[0047]
s1.1、设定其中一条裂缝的两端分别为a(a1,b1)和b(a2,b2)，另一条裂缝的两端分
别为c(c1,d1)和d(c2,d2)，则两条裂缝的参数方程为
[0048]
(x,y)＝(a1,b1) ξ1(a
2-a1,b
2-b1)
[0049]
(x,y)＝(c1,d1) ξ2(c
2-c1,d
2-d1)
ꢀꢀꢀ
(3)
[0050]
其中，ξ1和ξ2分别为裂缝交点在两条裂缝的参数方程上的参数值，0≤ξ1≤1，0≤ξ2≤1；
[0051]
由上式(3)得到关于ξ1和ξ2的线性方程组：
[0052][0053]
判断上式(4)中的线性方程组系数矩阵是否是奇异的，如果是奇异的，则和是平行的。此时，则没有交点需要被计算。如果非奇异，通过求解上式(4)得到ξ1和ξ2。需要指出的是，两条裂缝的交点需要在这两条裂缝上，即在线段和上，所以要求：0≤ξ1≤1，0≤ξ2≤1，如果计算得到的ξ1和ξ2满足0≤ξ1≤1，0≤ξ2≤1的要求，则ξ1和ξ2分别代表这两条裂缝交点对应参数方程中的参数值。
[0054]
s2、生成裂缝节点时，根据s1中得到的裂缝交点的参数向量计算裂缝的每段裂缝段的真实长度；
[0055]
具体的，在裂缝节点生成时，可以定义裂缝段的特征尺寸(记作l
f,char
)来控制裂缝节点之间裂缝段的长度(也可以说是裂缝段的分辨率)，在第一步中获取了裂缝之间的交点后，每条裂缝会被所得到的裂缝交点分成好几部分，可以用每条裂缝参数方程中的参数ξ构成的向量来表征。例如，对于某条裂缝f，其首尾两个端点对应的参数ξ取值分别是0或者1，通过上式4计算的该裂缝与其他裂缝的裂缝交点对应的参数ξ的取值分别是ξ
fi i＝1,
…
,ηf，其中，ηf是交点的数量。当然，ξ
fi
有可能等于0或者1，或者会存在i≠j时但ξ
fi
＝ξ
fj
，因此需要剔除参数向量中的重复值。在获取了每条裂缝的参数向量后，设某裂缝f的参数向量是ξf，该参数向量中相邻的两个参数值分别是ξ
fi
和ξ
fi 1
，该裂缝总长是lf，则这两个相邻的参数值ξ
fi
和ξ
fi 1
对应的裂缝f部分被剖分的段数计算为
[0056][0057]
其中，floor(x)函数是指取整函数，即当n≤x＜n 1时，n为整数，floor(x)＝n；
[0058]
s2.2、计算裂缝的每段裂缝段的真实长度l
real
:
[0059][0060]
s3、当对裂缝网络计算域进行非匹配性点云生成时：
[0061]
s3.1、忽略裂缝网络形态生成计算域的基质节点的生成；
[0062]
非匹配性点云直接忽略裂缝网络形态生成计算域的基质节点。对于矩形边界计算域，基质节点的生成可以采用笛卡尔式的点云，对于非矩形的具有不规则边界的计算域，可
以首先在y方向按照等间距构建一条条平行于x方向的坐标线，然后在x方向上根据坐标线与边界的交点情况按照等间距配置节点，节点之间的间距可以用x方向的特征长度(记作l
x,char
)和y方向上的特征长度(记作l
y,char
)来控制。可以看到，这种基质节点生成方法相较于计算域的网格剖分的确更简单。当然，由于非匹配式点云无需匹配裂缝网络几何形态，常用的各类网格剖分技术也可以被直接用来生成复杂边界计算域的分布比较均匀的基质节点。
[0063]
s3.2、建立基质节点与裂缝节点之间的相邻关系：设定基质节点与裂缝节点的控制区域均为边长为li的正方形区域，裂缝节点f的坐标为(xf,yf)，基质节点i的坐标为(xi,yi)，基质节点与裂缝节点的节点控制体积均记作vi：
[0064]
有了裂缝节点和基质节点后，本专利采用一个简单的方法来建立裂缝节点与基质节点之间的相邻关系，在定义了节点控制区域的无网格法中，可以计算得到节点控制体积，但一般难以刻画节点控制区域的具体形态。为了能够通过简单的节点间的距离信息来建立基质节点与裂缝节点之间的相邻关系并计算基质点与相邻裂缝点之间的传导率，可以假设节点控制区域是一个正方形区域。设某裂缝节点f，坐标是(xf,yf)，基质节点i的坐标是(xi,yi)，相应的节点控制体积记作vi，则确定哪个基质节点与裂缝节点f建立相邻关系的准则是：
[0065][0066]
式中，li是假设的正方形区域的边长。
[0067]
需要指出的是，上式(1)通过取最小值来确定哪个基质节点与裂缝节点f相邻的做法反映出：在本专利非匹配式复杂缝网点云生成算法中一个裂缝节点只会和一个基质节点间存在相邻关系，这一点类似于edfm中一个裂缝网格仅与其所在的基质网格之间存在相邻关系。
[0068]
s4、当对裂缝网络计算域进行匹配性点云生成时：
[0069]
s4.1、忽略裂缝网络形态生成计算域的基质节点的生成；
[0070]
s4.2、设定距离阈值，去掉与裂缝节点同位置的基质节点之间的距离小于所述距离阈值的基质节点：
[0071]
则去掉节点j，或者说第二类基质节点中的节点j被第一类基质节点中的节点i所取代。
[0072]
s4.2中与裂缝节点同位置的基质节点的坐标为(x1，y1)，去掉的基质节点的坐标为(x2，y2)，且：
[0073][0074]
其中，dis
tv
为距离阈值。
[0075]
离散计算域的点云相较于离散计算域的网格包含的拓扑信息少得多，导致用点云离散计算域比网格剖分简单，相应的点云生成方法比网格剖分也更加多样。对于如图1(a)所示的复杂缝网，图1(b)中复杂缝网的匹配性三角网格之所以会在裂缝交点以及裂缝之间的狭窄基质区域处有大量的小网格，是因为这些区域的几何实体之间的间距太小，在这些区域只能够形成小网格，且由于网格受到更多的拓扑结构限制，这些区域附近的网格只能
逐渐得由小变大，从而导致这些区域附近的网格密度很大。然而，点云不像网格会受到更多拓扑结构的限制，其在计算域上的布置在理论上可以更加灵活随意。本实施例实际上分别提供了复杂缝网匹配性和非匹配性点云生成算法，可以看到：
[0076]
(i)匹配式点云生成算法与非匹配式点云生成算法在裂缝交点计算和裂缝节点生成上是相同的。裂缝节点的分辨率由裂缝段特征长度控制。
[0077]
(ii)在非匹配式点云中，裂缝上的节点仅是裂缝节点，其基质节点的生成与裂缝节点完全独立，因此，非匹配式点云生成技术在生成基质节点离散计算域时类似于edfm那样无需去匹配裂缝形态，从而可以采用各种各样的高效简捷的计算域点云生成技术获取离散计算域基质部分的基质节点构成的点云，此时的基质节点分布一般情况下也是比较均匀的，这样也保证了无网格法对节点影响域内局部广义差分估计的精度。因此，非匹配式点云生成算法在s3中直接忽略裂缝网络生成计算域的基质节点分布。如前所述，因为点云所受拓扑结构的限制很少，所以在生成具有复杂边界计算域的基质节点时相较于网格剖分要更加灵活，在本实施例中的算法中，一般用x方向和y方向节点特征间距来控制基质节点的分辨率。在基质节点生成后，非匹配式算法就基于基质节点的控制体积和裂缝节点到各基质节点之间的距离信息，建立各裂缝节点与某个基质节点之间的非相邻连接。
[0078]
(iii)在匹配式点云中，由于裂缝节点相同位置也有基质节点(为了叙述的清晰，在后续描述中称为第一类基质节点)，因此，匹配式点云中基质节点的生成相较于非匹配式点云中基质节点的生成需要受到s2中生成的裂缝节点(即第一类基质节点)分布的限制，即要求后续生成的基质节点不能与s2中生成的裂缝节点距离过近或者重合。为了能够充分发挥出计算域点云离散相较于网格剖分的优势，在匹配式点云生成算法的s3中，首先采取与非匹配式点云生成算法s3中生成基质节点的相同做法来生成一个初始的基质节点分布(为了叙述的清晰，称为第二类基质节点)，然后在s4中将距离裂缝节点(即第一类基质节点)较近或者重合的第二类基质节点删除，来生成最终的基质节点分布。
[0079]
(iv)值得指出的是，在复杂缝网匹配性和非匹配性点云生成算法的点云生成算法中没有特别指出对井节点(包括直井节点、压裂井射孔点等)的处理方法。此处，我们指出井节点的处理也可以分为匹配性的和非匹配性的，其中非匹配性的处理方法是根据距离最近原则将直井节点取为距该直井节点最近的基质节点，或者将压裂射孔点取为距该压裂射孔点最近的裂缝点。该非匹配性的井节点处理方法因其简单即可以用在非匹配性点云生成中，也可以用在匹配性点云生成中。匹配性的处理方法是将已知坐标的直井节点作为预设的基质节点，或者将已知坐标的压裂射孔点作为预设的裂缝节点，该处理方法可用在匹配性点云生成中，相应得在匹配性点云生成算法的s2中生成裂缝节点前要把已知坐标的压裂射孔点作为已存在的裂缝点(即与s1中计算得到的裂缝交点是同等地位)，且在s4中将已知坐标的直井节点作为预设的基质节点纳入到第一类基质节点之中，即将s3生成的第二类基质节点中离该已知坐标的直井节点距离很近的节点去掉。总体来看，虽然匹配式的井节点处理方法相较于非匹配式的井节点处理方法能够使井节点的位置更符合实际位置，但非匹配式的井节点处理方法比匹配式的井节点处理方法更简单实用。
[0080]
整体上来看，无论是匹配性点云生成算法还是非匹配性点云生成算法，在前三个步骤中，生成裂缝节点和基质节点的步骤基本是相同的，且裂缝节点的生成与基质节点的生成完全独立，只是在s4中对基质节点和裂缝节点之间关系的处理有所不同。整个点云生
成算法流程十分简单，非匹配性点云只需要的l
x,char
、l
y,char
、l
f,char
(后续记作控制参数)来控制点云分辨率，非匹配式点云可以用l
x,char
、l
y,char
、l
f,char
和dis
tv
这四个控制参数来控制点云分辨率。
[0081]
《实施例2》
[0082]
根据本发明提供的匹配性点云生成算法，当令l
x,char
＝10m,l
y,char
＝10m,l
f,char
＝20m,dis
tv
＝5m时，图2(a)(c)(e)分别是图1(a)计算域匹配性点云中的所有节点、基质节点和裂缝节点示意图。可以看到，该匹配性点云相比于图1(b)中的匹配性三角网格，没有在裂缝交点以及裂缝之间的狭窄基质区域处产生大量的节点。当选取不同的控制参数值，例如图3和图4分别展示了l
x,char
＝20m,l
y,char
＝20m,l
f,char
＝20m,dis
tv
＝8m和l
x,char
＝10m,l
y,char
＝10m,l
f,char
＝10m,dis
tv
＝3m时生成的点云，可以看到图3中的点云和图4中的点云分别相较于图2(a)(c)(e)中的点云更粗糙和更精细，说明了本专利给出的复杂缝网点云生成算法中点云分辨率的调整十分灵活，相较于匹配性网格对计算域的加密更简单，便于开展自适应分析。
[0083]
根据本发明提供的非匹配性点云生成算法，当令l
x,char
＝10m,l
y,char
＝10m,l
f,char
＝20m时，图2(b)(d)(f)分别是图1(a)计算域非匹配性点云中的所有节点、基质节点和裂缝节点示意图。可以看到，非匹配性点云中基质节点和裂缝节点的分布完全独立，对于矩形计算域，则可以直接采用笛卡尔式的基质点云。
[0084]
对于图5(a)油藏边界不规则的情况，如前所述，由于本文给出的匹配性点云算法和非匹配性点云算法在前三个步骤中生成基质节点和裂缝节点时是独立的，因此可以采用忽略复杂缝网的油藏计算域的denaulay三角网格的网格节点分别作为匹配性点云生成算法中的第二类基质节点和非匹配性点云生成算法中的基质节点。图5(b)展示了对图5(a)中带复杂缝网油藏计算域的匹配性三角网格，该匹配性网格在裂缝交点和裂缝之间的狭长基质区域有大量的小网格，网格总节点数达到了20458。图6和图7则分别展示了本文算法所给出的匹配性点云和非匹配性点云，其中，图6匹配性点云中基质节点数1751，裂缝节点数1032，总节点数1751 1032*2＝3815。图7中的非匹配性点云中基质节点数2258，裂缝节点数1032，总节点数是2258 1032＝3290。可以看到，生成点云中的总节点数均远小于匹配性网格中的网格节点数。
[0085]
尽管本发明的实施方案已公开如上，但其并不仅仅限于说明书和实施方式中所列运用，它完全可以被适用于各种适合本发明的领域，对于熟悉本领域的人员而言，可容易地实现另外的修改，因此在不背离权利要求及等同范围所限定的一般概念下，本发明并不限于特定的细节和这里示出与描述的实施例。