[發(fā)明專利]一種大尺度裂縫性油藏數值模擬方法

有效

• 申請?zhí)枺?/td>: 201310138269.1
• 申請日: 2013-04-19
• 公開（公告）號: CN104112057B
• 公開（公告）日: 2018-02-13
• 發(fā)明（設計）人: 呂心瑞;張宏方;劉中春
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發(fā)研究院
• 主分類號: G06F19/00
• 分類號: G06F19/00
• 代理公司: 北京思創(chuàng)畢升專利事務所11218
• 代理人: 郭韞

1.一種大尺度裂縫性油藏數值模擬方法，其特征在于：所述方法將油藏中的大尺度裂縫作為油藏的內邊界進行降維處理，建立油藏幾何模型，然后采用三角形網格對油藏幾何模型進行幾何剖分形成離散單元，同時形成油藏幾何信息描述文件、幾何離散信息文件及井信息描述文件，在離散單元基礎上構建數值計算單元，建立表征此類油藏流體流動的數學模型，采用積分法在每個數值計算單元上建立其數值計算格式，形成線性代數方程組，并進行求解，從而實現大尺度裂縫性油藏數值模擬，

所述方法具體包括以下步驟：

(1)根據實際油藏特征及巖心、測井及流體的測試資料確定油藏儲層物性參數，建立表征大尺度裂縫性油藏流體流動的數學模型；

(2)依據油藏實際情況及油藏描述成果，確定油藏幾何邊界范圍，建立油藏幾何模型，形成油藏幾何信息描述文件；

(3)對步驟(2)建立的油藏模型進行幾何離散，生成油藏模型的幾何離散信息文件：平面上基巖部分采用具有厚度的Delaunay三角形單元進行劃分，縱向上所有三角形厚度相加即為儲層厚度，大尺度裂縫平面上用線單元代表，所述線單元具有代表裂縫直實高度的厚度，生成油藏模型的幾何離散信息文件；

(4)基于步驟(3)的幾何離散信息文件，連接所有三角形的重心與三邊的中點形成聯絡線，以三角形的頂點為中心，其周圍的所有聯絡線組成的多邊形即為數值計算單元；

(5)計算數值計算單元節(jié)點參數值：設定飽和度變量(S_w，S_o)在每個數值計算單元內恒定不變，其數值通過測井數據取得，而數值計算單元的流動勢變量(Φ_w，Φ_o，Φ_c)由以該數值計算單元的節(jié)點為頂點的相鄰三角網格的值通過線性逼近計算得到：

Φ(x)=Σi=1nNi(x)Φi---(7)]]>

其中，w,o分別表示水相和油相；S_w為油藏含水飽和度；S_o為油藏含油飽和度；Φ_w為油藏水流動勢；Φ_o為油藏原油流動勢；Φ_c為毛管壓力勢；n為三角形單元的頂點數；Φ_i代表節(jié)點i的任意流動勢變量；x為節(jié)點i的坐標；N_i(x)為形函數，為坐標的獨立函數；

由式(7)可得在一個三角形內任意變量的梯度恒為：

▿Φ=Σi=1nΦi▿Ni(x)---(8)]]>

(6)對步驟(1)建立的數學模型在步驟(4)建立的數值計算單元上進行積分，建立流動控制方程的數值計算方程，數值計算單元節(jié)點參數值是由步驟(5)計算得到的；

(7)確定油藏的邊界條件與初始條件，然后根據步驟(6)建立的每一個數值計算單元上的數值計算方程，構建線性代數方程組，并進行求解，完成大尺度裂縫性油藏的數值模擬，

所述步驟(2)包括：

(21)建立外邊界軌跡文件，該文件的內容包括：描述油藏邊界軌跡的點的X、Y坐標值及儲層深度值，其中，所述外邊界軌跡文件為油藏邊界范圍文件；

(22)將大尺度裂縫作為油藏內邊界進行簡化降維處理，根據油藏描述得到的大尺度裂縫具體位置，建立大尺度裂縫的具體信息文件，具體信息文件的內容包括；裂縫軌跡點對應的裂縫編號、X、Y坐標值、裂縫高度、裂縫開度；

(23)根據油藏中井的情況，建立井信息文件，該井信息文件的內容包括井號、X、Y坐標值、井所在的組別及井類別。