[发明专利]一种提高自动井位优化的优化效率方法

权利要求书

1.一种提高自动井位优化的优化效率方法，其特征在于，包括以下步骤：

1)油藏生产潜力初始化井位：利用生产潜力图初始化井位来代替粒子群算法随机初始化井位，油藏生产潜力图Quality map是用于表征油藏平面生产潜力大小的，生产潜力图可以通过式(1)求得：

J_i,j,k(t)＝[S_o,i,j,k(t)-S_or].[P_o,i,j,k(t)-P_min].LnK_i,j,k.Lnr_i,j,k.h_woc,i,j,k.Φ_i,j,k.h_goc,i,j,k (1)

式中，J_i,j,k(t)为网格(i,j,k)在t时刻的生产潜力；S_o,i,j,k(t)为网格(i,j,k)在t时刻的原始油饱和度；S_or为残余油饱和度；P_o,i,j,k(t)为网格(i,j,k)在t时刻的油相压力；P_min为最小的井底压力；K_i,j,k为网格(i,j,k)的渗透率；r_i,j,k为网格(i,j,k)距最近边界的距离；Φ_i,j,k.为网格(i,j,k)的孔隙度；h_goc,i,j,k为网格(i,j,k)到气水界面的距离；h_woc,i,j,k为网格(i,j,k)到油水界面的距离；

每个网格的潜力值在使用之前需要进行归一化处理，处理方法如式(2)；

式中，为网格(i,j,k)在t时刻的归一化生产潜力，J_max(t)为所有网格潜力值中的最大值；

2)根据步骤1)得到的初始化井位，调用数模软件进行完整油藏数值模拟运算，根据运算结果计算适应度值；

3)根据计算出的适应度值，结合速度和位置更新策略，更新粒子群算法中的速度和位置参数；

4)油藏生产潜力阈值约束：在上述结论中增加一个阈值参数ε，来评价每个粒子所提供的布井，如果该参数大于某一门限值，认为该布井是一个比较良好的选择，之后再调用数值模拟软件进行全部生产时间的模拟计算，否则调用数模软件只进行前半部分生产时间的模拟计算，而后半部分生产时间的累计产量则采用Gompertz产量预测模型完成；阈值参数ε如公式(3)：

其中，ε为潜力阈值；Q_now为目前迭代次数下粒子所提供的布井方案的生产潜力指数；Q_former为上一个迭代次数下粒子所提供的布井方案的生产潜力指数；Q为粒子所提供的布井方案的生产潜力指数；n_z为纵向网格数；n_w为优化的井数；

5)对于那些布井结果评价阈值低于ε的模拟方案，剩余的模拟时间采用数学方法进行预测模拟，对水驱油田的开发过程，油田的产油量由0到可采储量的过程，因此，可以选用Gompertz模型Gompertz model的规律来描述油田产量的变化过程；

y＝e^mn'+c       (5)

其中，y表示产油量，m，n，c可根据产油量的相关数据，回归得到；

对于通过上述(5)式，预测得到的累计产量比较好的布井方案，可以考虑利用数模重启文件继续进行精确模拟未完成的生产年份，定义参数η来评价是否有必要用重启文件继续模拟未完成的生产年份，如果该参数大于某一门限值，那么此时可以认为该布井是一个比较良好的选择，可以进行调用数模进行验证；如公式(6)所示：

其中，η为阈值参数；N_pnow为目前迭代次数下粒子所提供的布井方案的累计产油量；N_pformer为上一个迭代次数下粒子所提供的布井方案的累计产油量；

6)判断是否满足算法迭代停止条件，即最大迭代代数；如果不满足转向步骤3)，否则停止，最终得到最优井位和优化结果。