[发明专利]一种复杂裂缝网络性油藏流动模拟方法

权利要求书

1.一种复杂裂缝网络性油藏流动模拟方法，其特征在于，包括如下步骤：

步骤1：根据压裂水平井控制区域复杂裂缝网络分布规律将油藏离散成一系列节点，采用节点间流动路径表征流体在不同尺度介质中流动特征；

步骤2：利用节点物性、影响域和加权最小二乘算法估计压力拉普拉斯算子，给出反映连接单元渗流能力的传导率的计算方法以及物质基础的属性参数的连通体积的计算方法；

步骤3：以连接单元为对象，通过节点压力方程求解获取单元内的流动通量，沿着连接单元进行饱和度方程的半解析追踪计算产出动态，得到不同射孔点生产动态，实现对复杂裂缝网络性油藏流动的模拟。

2.根据权利要求1所述的复杂裂缝网络性油藏流动模拟方法，其特征在于，所述步骤1中针对具有复杂裂缝网络分布特征的压裂水平井，采用节点体系对储层油藏进行等效刻画，再针对节点体系的节点分布情况，根据储层流体运移特征连接节点表征流体流动路径，基于节点体系和节点间流动路径连接原则，对非常规水平井压裂改造储层进行储层表征和复杂裂缝形态刻画。

3.根据权利要求2所述的复杂裂缝网络性油藏流动模拟方法，其特征在于，根据裂缝分布特征将储层划分为三个区域，分别为区域1、区域2和区域3；

所述区域1是压裂次生裂缝分布位置，在次生裂缝根部、前端以及各拐点处分别布置③类节点，根据裂缝扩展规律将各所述③类节点依次连接，精确刻画裂缝形态，其中，所述③类节点处的物性参数采用裂缝物性参数；

所述区域2为压裂次生裂缝周围区域，在所述区域2中划分②类节点，以精确描述压裂次生裂缝周围区域内基质向裂缝的复杂流动特征；

所述区域3为储层未改造区域，在所述区域3中划分①类节点，所述①类节点用以表征流体从未改造区域向改造区域内的流体；

其中，所述③类节点的分布密度大于所述①类节点的分布密度，所述②类节点的分布密度大于所述③类节点的分布密度；

其中，所述①类节点和②类节点处的物性参数采用基质物性参数，所述基质物性参数可根据精细地质模型进行非均质赋值。

4.根据权利要求3所述的复杂裂缝网络性油藏流动模拟方法，其特征在于，所述步骤1中的节点间流动路径包括基质内部流动路径、基质-压裂裂缝流动路径、压裂裂缝内部流动路径；

所述基质内部流动路径为①类节点与其周边的①类节点全部连接，且任意相连接的两个①类节点之间的流体的流动方向为双向流动；

所述基质-压裂裂缝流动路径为压裂裂缝周围的②类节点与周围的①类节点以及距其最近的③类节点均连接；

所述压裂裂缝内部流动路径为裂缝内流体沿着裂缝形态流向井筒处。

5.根据权利要求4所述的复杂裂缝网络性油藏流动模拟方法，其特征在于，所述步骤2中所述传导率T_ij的计算公式如式(4)：

其中，V_i是节点i的控制体积，单位为m³；

λ_ij为流体总流度，单位为mD·(mPa·s)^-1，

其中，

k_ro和k_rw分别是油相和水相相对渗透率，单位为mD；

μ_o和μ_w分别是油相和水相粘度，单位为mPa·s；

是Laplace算子广义有限差分估计的系数，单位为1；

α是单位换算系数，取值0.0864；

k_ij是连接单元(i·j)的平均渗透率，单位为mD，取连接单元i节点处的渗透率k_i和j节点处的渗透率k_j的调和平均，如式(5)：