[发明专利]一种考虑非平衡各向异性相对渗透率的油藏状态预测方法

权利要求书

1.一种考虑非平衡各向异性相对渗透率的油藏状态预测方法，其特征在于，包括：

对三维立体地质模型进行网格区域划分，得到井参数和各网格区域的初始属性参数；所述井参数包括井点位置、生产速度、生产压力和生产时间；所述初始属性参数包括初始压力、初始流体流速和初始含水饱和度；相邻两网格区域中的一个网格区域的初始流体流速是根据相邻两网格区域的初始压力差值以及相邻两网格区域中另一网格区域的初始流体流速计算得到的；所述三维立体地质模型是对待预测区域的河流相储层的结构进行地质建模得到的；

根据所述井点位置所在的网格区域确定第一网格区域，并根据所述生产压力与所述初始压力的差值以及所述生产速度更新所述第一网格区域的初始压力和初始流体流速；更新后的第一网格区域的初始压力为所述生产压力与所述初始压力的差值；更新后的第一网格区域的初始流体流速为所述生产速度；

以所述第一网格区域为中心，依次更新与所述第一网格区域相邻的各网格区域的初始压力和初始流体流速，然后以更新后的任意一网格区域为中心，依次更新与所述更新后的任意一网格区域相邻的各网格区域的初始压力和初始流体流速，直至更新完所有网格区域的初始压力和初始流体流速；

根据各网格区域的初始含水饱和度和更新后的各网格区域的初始流体流速确定各网格区域在X、Y、Z三个方向上的相对渗透率，并根据各网格区域在X、Y、Z三个方向上的相对渗透率以及更新后的各网格区域的初始压力和初始流体流速更新各网格区域的初始含水饱和度；

根据各网格区域在X、Y、Z三个方向上的相对渗透率以及更新后的各网格区域的初始压力、初始流体流速和初始含水饱和度得到下一时刻各网格区域的压力、流体流速和含水饱和度；

将下一时刻各网格区域的压力、流体流速和含水饱和度作为初始属性参数循环计算至所述生产时间，得到各个时刻对应的各网格区域的压力、流体流速和含水饱和度，然后根据各个时刻对应的各网格区域的压力、流体流速和含水饱和度绘制河流相储层各区域的含水饱和度和压力随时间变化图；

根据所述含水饱和度和压力随时间变化图得到河流相储层中油藏和压力分布情况。

2.根据权利要求1所述的考虑非平衡各向异性相对渗透率的油藏状态预测方法，其特征在于，所述根据各网格区域的初始含水饱和度和更新后的各网格区域的初始流体流速确定各网格区域在X、Y、Z三个方向上的相对渗透率，具体包括：

根据更新后的各网格区域的初始流体流速确定各网格区域在X、Y、Z三个方向上的相对渗透率曲线；所述X、Y、Z三个方向上的相对渗透率曲线的横坐标为含水饱和度，纵坐标为X、Y、Z三个方向上的相对渗透率；所述X、Y、Z三个方向上的相对渗透率曲线是对X方向岩心柱、Y方向岩心柱、Z方向岩心柱分别进行驱替实验得到的；所述X方向岩心柱、所述Y方向岩心柱和所述Z方向岩心柱是对待预测区域的河流相储层的野外露头中选取的正方体形状的样品分别从X、Y、Z三个方向上钻取得到的；

根据各网格区域的初始含水饱和度以及各网格区域在X、Y、Z三个方向上的相对渗透率曲线，确定各网格区域在X、Y、Z三个方向上的相对渗透率。

3.根据权利要求1所述的考虑非平衡各向异性相对渗透率的油藏状态预测方法，其特征在于，所述根据各网格区域在X、Y、Z三个方向上的相对渗透率以及更新后的各网格区域的初始压力和初始流体流速更新各网格区域的初始含水饱和度，具体包括：

将各网格区域的油水两相模型渗流控制方程中的相对渗透率替换为各网格区域在X、Y、Z三个方向上的相对渗透率，得到相对渗透率为X、Y、Z三个方向形式的各网格区域的油水两相模型渗流控制方程；

将更新后的各网格区域的初始压力和初始流体流速代入所述相对渗透率为X、Y、Z三个方向形式的各网格区域的油水两相模型渗流控制方程中，结合含水饱和度与含油饱和度的约束条件，计算各网格区域的含水饱和度；

以所述各网格区域的含水饱和度更新各网格区域的初始含水饱和度。

4.根据权利要求3所述的考虑非平衡各向异性相对渗透率的油藏状态预测方法，其特征在于，所述根据各网格区域在X、Y、Z三个方向上的相对渗透率以及更新后的各网格区域的初始压力、初始流体流速和初始含水饱和度得到下一时刻各网格区域的压力、流体流速和含水饱和度，具体包括：

将更新后的各网格区域的初始压力、初始流体流速和初始含水饱和度代入所述相对渗透率为X、Y、Z三个方向形式的各网格区域的油水两相模型渗流控制方程中，计算下一时刻各网格区域的压力、流体流速和含水饱和度。