非侵入式模拟复杂裂缝的EDFM方法及装置与流程

技术特征：

1.一种非侵入式模拟复杂裂缝的EDFM方法，其特征在于，包括：

获取裂缝属性参数、油藏地质属性参数和井筒属性参数；所述裂缝属性参数包括：裂缝条数、每条裂缝的位置坐标、每条裂缝的宽度、每条裂缝的方位角和每条裂缝的渗透率；所述油藏地质属性参数包括：油藏基质渗透率、油藏大小、基质网格、油藏深度；所述井筒属性参数包括：井的数量、井的位置坐标和井的半径；

确定每一条所述裂缝与所述基质网格之间的几何相交关系，并基于所述基质网格的边缘离散所述裂缝，将所述裂缝分割为一系列的裂缝段，每个裂缝段对应一个裂缝网格；

确定所述裂缝网格的数量，并增加所述裂缝网格的属性参数；

计算所述裂缝网格的连接关系和相应连接的传导率，所述连接关系包括：裂缝网格与基质网格之间的NNC连接关系、裂缝网格与裂缝网格之间的NNC连接关系、裂缝网格与井筒之间的连接关系；所述传导率包括：裂缝网格与基质网格之间的传导率、裂缝网格与裂缝网格之间的传导率、裂缝网格与井筒之间的传导率；所述连接关系和传导率用于确定流体在裂缝中和油藏地质中的流动和传质。

2.根据权利要求1所述的方法，其特征在于，还包括：

根据所述连接关系和传导率修改第三方油藏数值模拟软件的相对应的关键字；

基于所述第三方油藏数值模拟软件耦合求解流体在裂缝中和油藏地质中的流动和传质。

3.根据权利要求1所述的方法，其特征在于，计算裂缝网格与基质网格之间的传导率包括：

与裂缝网格对应的裂缝段穿过基质网格，并将基质网格分成第一单元和第二单元，且所述第一单元与裂缝段之间的面为第一裂缝面，所述第二单元与裂缝段之间的面为第二裂缝面；

确定所述第一单元到所述第一裂缝面的垂直距离df-A和所述第二单元到所述第二裂缝面的垂直距离df-B；且其中VA表示所述第一单元的体积，xnA表示第一单元的单元体到第一裂缝面的垂直距离；

确定流体从所述第一裂缝面到所述第一单元的传质速度Qf-A和流体从所述第二裂缝面到所述第二单元的传质速度Qf-B；其中，Qf-A＝Tf-Aλi(pf-pA)；pf代表裂缝段中的平均压力，Tf-A代表从所述第一裂缝面到所述第一单元的传导率，λi代表流体的相对流动性；其中，Af代表第一裂缝面的面积，K代表基质网格的渗透率，代表第一裂缝面的法向矢量；

确定裂缝段与所述裂缝段穿过的基质网格之间的传导率：

其中，df-m是基质网格到裂缝面的垂直距离，且xn是基质网格的单元体到裂缝面的垂直距离，V表示基质网格的体积。

4.根据权利要求1所述的方法，其特征在于，计算裂缝网格与裂缝网格之间的传导率包括：

计算不同裂缝间的第一裂缝网格和第二裂缝网格之间的传导率Tint：

其中，Lint是两条不同裂缝相交线的长度，df1和df2分别是两条不同的裂缝到相交线的加权平均距离：

式中，dSi是第i个裂缝面的裂缝积分面单元，Si代表不同裂缝被离散成的第i个裂缝面的面积，i＝1，2，3，4，且第1个和第3个裂缝面为与第一裂缝网格对应的裂缝的面，第2个和第4个裂缝面为与第二裂缝网格对应裂缝的面；xni代表第i个裂缝面的裂缝面单元到相交线的距离；

计算同一条裂缝内的第三裂缝网格和第四裂缝网格之间的传导率Tseg：

其中，

式中，kf为裂缝段的渗透率，Ac是两个裂缝片段共有接触面的面积，dseg3和dseg4分别是与第三裂缝网格和第四裂缝网格对应的裂缝段的中心到共有接触面的距离。

5.根据权利要求1所述的方法，其特征在于，计算裂缝网格与井筒之间的传导率包括：

计算裂缝网格与穿过与所述裂缝网格对应的裂缝段的井筒之间的传导率：

其中，kf是裂缝的渗透率，wf是裂缝宽度，L和W分别表示裂缝段的长度和高度，rw是井筒半径。

6.一种非侵入式模拟复杂裂缝的EDFM装置，其特征在于，包括：

获取模块，用于获取裂缝属性参数、油藏地质属性参数和井筒属性参数；所述裂缝属性参数包括：裂缝条数、每条裂缝的位置坐标、每条裂缝的宽度、每条裂缝的方位角和每条裂缝的渗透率；所述油藏地质属性参数包括：油藏基质渗透率、油藏大小、基质网格、油藏深度；所述井筒属性参数包括：井的数量、井的位置坐标和井的半径；

离散模块，用于确定每一条所述裂缝与所述基质网格之间的几何相交关系，并基于所述基质网格的边缘离散所述裂缝，将所述裂缝分割为一系列的裂缝段，每个裂缝段对应一个裂缝网格；

确定模块，用于确定所述裂缝网格的数量，并增加所述裂缝网格的属性参数；

计算模块，用于计算所述裂缝网格的连接关系和相应连接的传导率，所述连接关系包括：裂缝网格与基质网格之间的NNC连接关系、裂缝网格与裂缝网格之间的NNC连接关系、裂缝网格与井筒之间的连接关系；所述传导率包括：裂缝网格与基质网格之间的传导率、裂缝网格与裂缝网格之间的传导率、裂缝网格与井筒之间的传导率；所述连接关系和传导率用于确定流体在裂缝中和油藏地质中的流动和传质。

7.根据权利要求6所述的装置，其特征在于，还包括：

修改模块，用于根据所述连接关系和传导率修改第三方油藏数值模拟软件的相对应的关键字；

处理模块，用于基于所述第三方油藏数值模拟软件耦合求解流体在裂缝中和油藏地质中的流动和传质。

8.根据权利要求6所述的装置，其特征在于，所述计算模块用于计算裂缝网格与基质网格之间的传导率，包括：

与裂缝网格对应的裂缝段穿过基质网格，并将基质网格分成第一单元和第二单元，且所述第一单元与裂缝段之间的面为第一裂缝面，所述第二单元与裂缝段之间的面为第二裂缝面；

确定所述第一单元到所述第一裂缝面的垂直距离df-A和所述第二单元到所述第二裂缝面的垂直距离df-B；且其中VA表示所述第一单元的体积，xnA表示第一单元的单元体到第一裂缝面的垂直距离；

确定流体从所述第一裂缝面到所述第一单元的传质速度Qf-A和流体从所述第二裂缝面到所述第二单元的传质速度Qf-B；其中，Qf-A＝Tf-Aλi(pf-pA)；pf代表裂缝段中的平均压力，Tf-A代表从所述第一裂缝面到所述第一单元的传导率，λi代表流体的相对流动性；其中，Af代表第一裂缝面的面积，K代表基质网格的渗透率，代表第一裂缝面的法向矢量；

确定裂缝段与所述裂缝段穿过的基质网格之间的传导率：

其中，df-m是基质网格到裂缝面的垂直距离，且xn是基质网格的单元体到裂缝面的垂直距离，V表示基质网格的体积。

9.根据权利要求6所述的装置，其特征在于，所述计算模块用于计算裂缝网格与裂缝网格之间的传导率，包括：

计算不同裂缝间的第一裂缝网格和第二裂缝网格之间的传导率Tint：

其中，Lint是两条不同裂缝相交线的长度，df1和df2分别是两条不同的裂缝到相交线的加权平均距离：

式中，dSi是第i个裂缝面的裂缝积分面单元，Si代表不同裂缝被离散成的第i个裂缝面的面积，i＝1，2，3，4，且第1个和第3个裂缝面为与第一裂缝网格对应的裂缝的面，第2个和第4个裂缝面为与第二裂缝网格对应裂缝的面；xni代表第i个裂缝面的裂缝面单元到相交线的距离；

计算同一条裂缝内的第三裂缝网格和第四裂缝网格之间的传导率Tseg：

其中，

式中，kf为裂缝段的渗透率，Ac是两个裂缝片段共有接触面的面积，dseg3和dseg4分别是与第三裂缝网格和第四裂缝网格对应的裂缝段的中心到共有接触面的距离。

10.根据权利要求6所述的装置，其特征在于，所述计算模块用于计算裂缝网格与井筒之间的传导率，包括：

计算裂缝网格与穿过与所述裂缝网格对应的裂缝段的井筒之间的传导率：

其中，kf是裂缝的渗透率，wf是裂缝宽度，L和W分别表示裂缝段的长度和高度，rw是井筒半径。