[发明专利]孔隙喉道网络模型驱替模拟中驱油效率的确定方法及装置

权利要求书

1.一种孔隙喉道网络模型驱替模拟中驱油效率的确定方法，其特征在于，包括：

获取孔隙喉道介质的特征参数，并根据所述特征参数构建二维孔隙喉道网络模型；

驱替模拟所述二维孔隙喉道网络模型的渗流规律；

在模拟所述渗流规律的过程中，当确认所述二维孔隙喉道网络模型内喉道单元的顶点压力达到临界压力值时，使对应的喉道单元张开而形成裂缝，以更新所述二维孔隙喉道网络模型；

在模拟所述渗流规律的过程中，实时获取不同阶段的二维孔隙喉道网络模型的驱油效率；并实时获取所述二维孔隙喉道网络模型的注入PV数；对所述二维孔隙喉道网络模型的注入PV数及驱油效率进行关系拟合，获得注入PV数于驱油效率的关系曲线。

2.如权利要求1所述的孔隙喉道网络模型驱替模拟中驱油效率的确定方法，其特征在于，所述实时获取不同阶段的二维孔隙喉道网络模型的驱油效率，包括：

根据公式计算不同阶段的二维孔隙喉道网络模型的驱油效率；

其中，E_d为二维孔隙喉道网络模型的驱油效率，S_wi为二维孔隙喉道网络模型的束缚水饱和度，S_or为二维孔隙喉道网络模型的残余油饱和度。

3.如权利要求1所述的孔隙喉道网络模型驱替模拟中驱油效率的确定方法，其特征在于，所述实时获取所述二维孔隙喉道网络模型的注入PV数，包括：

根据公式获取所述二维孔隙喉道网络模型在第n驱替步的总注入量；

将第n驱替步的总注入量除以所述二维孔隙喉道网络模型的总孔隙体积，获得所述二维孔隙喉道网络模型在第n驱替步的注入PV数；

其中，V_n为二维孔隙喉道网络模型在第n驱替步的总注入量，V_o为初始状态下二维孔隙喉道网络模型中被侵入的孔隙单元及喉道单元的体积之和，V_i1为二维孔隙喉道网络模型中的喉道单元在第i驱替步后已被侵入的体积，V_i2为二维孔隙喉道网络模型中的孔隙单元在第i驱替步后已被侵入的体积。

4.如权利要求1所述的孔隙喉道网络模型驱替模拟中驱油效率的确定方法，其特征在于，所述根据所述特征参数构建二维孔隙喉道网络模型，包括：

确定用于表征所述特征参数的分布规律的统计分布函数；

按照所述统计分布函数生成孔隙喉道网络模型，并将所述特征参数赋值给所述孔隙喉道网络模型中的喉道单元和孔隙单元，从而形成二维孔隙喉道网络模型。

5.如权利要求4所述的孔隙喉道网络模型驱替模拟中驱油效率的确定方法，其特征在于，所述确定用于表征所述特征参数的分布规律的统计分布函数，包括：

基于公式表征孔隙半径的分布规律；

基于公式表征喉道半径的分布规律；

其中，p(R)为孔隙半径的分布函数，f(x)为喉道半径的分布函数，R为孔隙半径，R_max为孔隙半径的最大值，R_min为孔隙半径的最小 值，σ为分布函数的标准差，e为自然常数，μ为分布函数的期望值，x为喉道半径。

6.如权利要求1所述的孔隙喉道网络模型驱替模拟中驱油效率的确定方法，其特征在于，所述驱替模拟所述二维孔隙喉道网络模型的渗流规律，包括：

根据公式模拟所述二维孔隙喉道网络模型中的喉道单元的渗流规律；

根据公式模拟所述二维孔隙喉道网络模型中的裂缝的渗流规律；

其中，q_ij为二维孔隙喉道网络模型中第i行第j列位置处的流量，r_ij为二维孔隙喉道网络模型中第i行第j列位置处的喉道半径，为流体的动力粘滞系数，L_ij为二维孔隙喉道网络模型中至第i行第j列位置处的裂缝长度，P_j为二维孔隙喉道网络模型中第j列的压力，P_i为二维孔隙喉道网络模型中第i行的压力，θ为接触角，σ^wn为润湿相和非润湿相间的界面张力，r为毛管半径，q为孔隙单元内的流量，e为自然常数，μ为流体的动力粘滞系数，ΔP＝P_j-P_i-P_c。