[發明專利]一種油藏水驅油開發動態分析方法、裝置及可讀存儲介質

有效

• 申請號: 201910451637.5
• 申請日: 2019-05-28
• 公開（公告）號: CN112100788B
• 公開（公告）日: 2023-06-27
• 發明（設計）人: 雷剛;曹耐;吳亞紅;李虹
• 申請（專利權）人: 中國石油大學（北京）
• 主分類號: G06F30/20
• 分類號: G06F30/20;E21B43/20;G06F119/14
• 代理公司: 北京集佳知識產權代理有限公司 11227
• 代理人: 羅滿

1.一種油藏水驅油開發動態分析方法，其特征在于，包括：

依據預先采集的油藏數據信息及動態參數模型得到與油藏對應的各個生產指標參數；

依據各個所述生產指標參數建立油藏開發方案；其中，所述動態參數模型為預先基于儲層介質形變建立的；其中：

所述動態參數模型的建立過程為：

依據預設條件建立基于儲層介質形變的油水兩相滲流數學模型；

依據所述油水兩相滲流數學模型建立壓力模型、流模型、有效動用系數模型及含水飽和度模型；

所述油水兩相數學模型包括水相速度方程、油相速度方程和調和方程，其中：

所述水相速度方程為所述油相速度方程為所述調和方程為▽·(▽U)＝0，其中，k(p)＝k₀exp[α(p-p₀)]，U＝exp[α(p-p₀)]，為水相滲流速度，k(p)為巖石絕對滲透率，k_rw為水相相對滲透率，p為地層壓力，▽為Laplace算子，μ_w為水相流體粘度，為油相滲流速度，k_ro為油相相對滲透率，μ_o為油相流體粘度，k₀為巖石初始滲透率，α為介質變形系數，p₀為初始地層壓力，U為無因次壓力；

所述壓力模型包括無窮井網在地層中任意一點M的擬勢函數方程和地層壓力分布函數方程，其中：

所述無窮井網在地層中任意一點M的擬勢函數為：

其中，ψ_M為擬勢函數，Q為油井產量，h為有效油層厚度，n為油井或注入井數量，m為井排排數，x為M點的橫坐標，y為M點的縱坐標，a為y軸方向井間距，d為x軸方向井間距，C_mn為第m排井排中第n口油井所在位置處的常數，C′_mn為第m排井排中第n口注入井所在位置處的常數；

所述地層壓力分布函數方程為：

其中，τ＝Qα/(4πhλ_t)，C₁＝C/λ_t，μ_rw為水相相對流體粘度，μ_ro為油相相對流體粘度；

所述流模型為基于流函數建立的，所述流函數為：

其中：

U(1,0)-U(0,0)表示無因次注采壓差，Γ為流函數，C₂為常數；

所述有效動用系數模型為基于井網單元有效動用系數函數建立的，所述井網單元有效動用系數函數為：

其中，ξ為井網單元有效動用系數，T為見水時間，為地層孔隙度，A為井網單元面積，其中：

T＝min(T_x,T_y)，

其中，

所述含水飽和度模型為基于等飽和度運行特征方程建立的，其中所述等飽和度運行特征方程為：

其中，t為生產時間，f′_w(S_w)為分流量函數對含水飽和度S_w的導數，S_w為含水飽和度，v_tx(ω,y_D)為點(ω,y_D)處在x方向上的滲流速度，ω為積分變量，v_ty(x_D,ω)為點(x_D,ω)處在y方向上的滲流速度。

2.根據權利要求1所述的油藏水驅油開發動態分析方法，其特征在于，所述預設條件包括：

水平均質等厚的無限大地層上布有無限延伸的五點井網；地層只有油水兩相，并且流體在多孔介質中的流動均為二維流動；流體在多孔介質中流動不受毛管壓力影響；油井以定液量生產，注入井注入量與油井產液量相等；多孔介質發生形變時，儲層滲流隨壓力發生變化，孔隙度不隨壓力變化。