[發(fā)明專利]一種油藏的動態(tài)離散裂縫模擬方法和裝置

在審

• 申請?zhí)枺?/td>: 201410612927.0
• 申請日: 2014-11-04
• 公開（公告）號: CN104392109A
• 公開（公告）日: 2015-03-04
• 發(fā)明（設(shè)計）人: 雷征東;田昌炳;石成方;王文環(huán);王方;彭緩緩;侯建鋒;李軍詩;王錦芳
• 申請（專利權(quán)）人: 中國石油天然氣股份有限公司
• 主分類號: G06F19/00
• 分類號: G06F19/00;G06F17/50
• 代理公司: 北京三友知識產(chǎn)權(quán)代理有限公司 11127
• 代理人: 黨曉林;李永強

1.一種油藏的動態(tài)離散裂縫模擬方法，其特征在于，包括：

對油藏進行網(wǎng)格剖分，得到油藏的網(wǎng)格；

對所述油藏的網(wǎng)格進一步劃分為多重次級網(wǎng)格；

基于所述油藏的網(wǎng)格和多重次級網(wǎng)格，計算基質(zhì)與基質(zhì)、基質(zhì)與裂縫、裂縫與裂縫之間的傳導(dǎo)率，以及計算儲層與井筒的流體交換系數(shù)，并進一步得到每一時間步長下的壓力和飽和度；

建立注水過程中裂縫形成和動態(tài)延伸的判別準則，并建立天然裂縫與注水動態(tài)裂縫、壓裂縫之間的耦合作用關(guān)系，判斷不同時間點的裂縫延伸情況；

以延伸的新裂縫為邊界條件動態(tài)劃分網(wǎng)格，同時將所述新裂縫的屬性信息賦給網(wǎng)格邊界，以更新所述網(wǎng)格邊界；

以所述更新后的網(wǎng)格邊界重新計算壓力場和飽和度場；

修正油藏與巖石力學參數(shù)，進行歷史數(shù)據(jù)擬合，計算當前時間步的飽和度和壓力情況。

2.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述對油藏進行網(wǎng)格剖分包括：

采用以裂縫面為幾何約束的Delaunary剖分方法對所述油藏進行網(wǎng)格剖分，將裂縫處理為曲面，將基質(zhì)處理為四面體網(wǎng)格。

3.根據(jù)權(quán)利要求2所述的方法，其特征在于，所述對油藏進行網(wǎng)格剖分還包括：

將基質(zhì)網(wǎng)格體積減去相鄰的裂縫網(wǎng)格的體積的一半。

4.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述方法還包括：

判斷是否出現(xiàn)新裂縫。

5.根據(jù)權(quán)利要求4所述的方法，其特征在于，所述判斷是否出現(xiàn)新裂縫包括：

以原生裂縫的尖端坐標作為新裂縫的起點坐標，在每個時間步判斷與新裂縫相連的原生裂縫尖端網(wǎng)格的壓力是否大于其臨界壓力，若大于其臨界壓力，則判斷為出現(xiàn)新裂縫。

6.根據(jù)權(quán)利要求1所述的方法，其特征在于，所述新裂縫的屬性信息賦給網(wǎng)格邊界包括：

如果新裂縫產(chǎn)生后的系統(tǒng)的裂縫網(wǎng)格中心與原系統(tǒng)的裂縫網(wǎng)格中心重合或者在原系統(tǒng)的裂縫網(wǎng)格上，則將原系統(tǒng)的網(wǎng)格壓力和飽和度直接賦值給新裂縫產(chǎn)生后的系統(tǒng)的裂縫網(wǎng)格上；

如果新裂縫產(chǎn)生后的系統(tǒng)的裂縫網(wǎng)格中心在原系統(tǒng)的基質(zhì)網(wǎng)格中，則將原系統(tǒng)的裂縫網(wǎng)格壓力和飽和度賦值到新裂縫產(chǎn)生后的系統(tǒng)的裂縫網(wǎng)格上；

如果新裂縫產(chǎn)生后的系統(tǒng)的基質(zhì)網(wǎng)格中心與原系統(tǒng)的基質(zhì)網(wǎng)格中心重合或在原系統(tǒng)的基質(zhì)網(wǎng)格中，則將原系統(tǒng)的基質(zhì)網(wǎng)格的壓力和飽和度直接賦值到新裂縫產(chǎn)生后的系統(tǒng)的基質(zhì)網(wǎng)格上。

7.一種油藏的動態(tài)離散裂縫模擬裝置，其特征在于，包括：

油藏網(wǎng)格劃分單元，用于對油藏進行網(wǎng)格劃分，得到油藏的網(wǎng)格；

次級網(wǎng)格劃分單元，用于對所述油藏的網(wǎng)格進一步劃分為多重次級網(wǎng)格；

流體交換系統(tǒng)計算單元，用于基于所述油藏的網(wǎng)格和多重次級網(wǎng)格，計算基質(zhì)與基質(zhì)、基質(zhì)與裂縫、裂縫與裂縫之間的傳導(dǎo)率，以及計算儲層與井筒的流體交換系數(shù)，并進一步得到每一時間步長下的壓力和飽和度；

裂縫延伸情況判斷單元，用于建立注水過程中裂縫形成和動態(tài)延伸的判別準則，并建立天然裂縫與注水動態(tài)裂縫、壓裂縫之間的耦合作用關(guān)系，判斷不同時間點的裂縫延伸情況；

屬性信息賦值單元，用于以延伸的新裂縫為邊界條件動態(tài)劃分網(wǎng)格，同時將所述新裂縫的屬性信息賦給網(wǎng)格邊界，以更新所述網(wǎng)格邊界；

壓力飽和度修正單元，用于以所述更新后的網(wǎng)格邊界重新計算壓力場和飽和度場；

當前時間步計算單元，用于修正油藏與巖石力學參數(shù)，進行歷史數(shù)據(jù)擬合，計算當前時間步的飽和度和壓力情況。

8.根據(jù)權(quán)利要求7所述的裝置，其特征在于，所述油藏網(wǎng)格劃分單元具體用于采用以裂縫面為幾何約束的Delaunary剖分方法對所述油藏進行網(wǎng)格剖分，將裂縫處理為曲面，將基質(zhì)處理為四面體網(wǎng)格。

9.根據(jù)權(quán)利要求8所述的裝置，其特征在于，所述油藏網(wǎng)格劃分單元還用于將基質(zhì)網(wǎng)格體積減去相鄰的裂縫網(wǎng)格的體積的一半。

10.根據(jù)權(quán)利要求7所述的裝置，其特征在于，所述裝置還包括：新裂縫出現(xiàn)判斷單元，用于判斷是否出現(xiàn)新裂縫。