[發(fā)明專利]儲層裂縫預測方法及系統(tǒng)

在審

• 申請?zhí)枺?/td>: 201610758001.1
• 申請日: 2016-08-29
• 公開（公告）號: CN107783204A
• 公開（公告）日: 2018-03-09
• 發(fā)明（設計）人: 劉玄;計秉玉;呂成遠;倫增珉;趙淑霞;王友啟;何應付;潘偉義;周銀邦;張慶新
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發(fā)研究院
• 主分類號: G01V9/00
• 分類號: G01V9/00
• 代理公司: 北京思創(chuàng)畢升專利事務所11218
• 代理人: 孫向民

技術領域

本發(fā)明涉及地震勘探領域，更具體地，涉及一種儲層裂縫預測方法及系統(tǒng)。

背景技術

在石油開采技術領域中，儲層注水、氣、聚合物開發(fā)等提高采收率過程中，儲層裂縫對流體流動影響巨大，并直接影響開發(fā)效果：儲層天然主裂縫往往成為水竄、氣竄的重要通道，儲層天然裂縫預測位置不準確會導致注水、氣效果下降，甚至使生產(chǎn)井被迫關井。

目前對儲層天然裂縫位置、走向、導流能力等的識別方法主要有三種：首先是根據(jù)成像測井的結果進行判斷，該方法識別裂縫準確，但僅能發(fā)現(xiàn)與井筒相交的裂縫面，無法獲得未穿過井筒的裂縫信息；其次是通過地震數(shù)據(jù)進行判斷，但由于大部分儲層裂縫的的長度均較小，因此無法通過地震做到準確識別；第三種方法是根據(jù)注采井的生產(chǎn)動態(tài)數(shù)據(jù)進行反推，從而獲得天然主裂縫的信息，該方法在數(shù)學上是一個典型的反問題求解過程，難以實現(xiàn)。另外，集合卡曼濾波方法(ENKF)是處理反問題的方法之一，目前已在儲層物性反演、儲層裂縫預測方面取得了一些應用。

發(fā)明人發(fā)現(xiàn)，現(xiàn)有的基于集合卡曼濾波方法的裂縫預測技術均建立在常規(guī)的結構化網(wǎng)格(正交網(wǎng)格或角點網(wǎng)格等)基礎之上，此類網(wǎng)格在處理裂縫型儲層時都只能做簡化處理，必須將裂縫的導流能力等效進基質(zhì)網(wǎng)格中，因此無法精確且簡便的模擬裂縫型儲層，相應的裂縫反演結果也就不夠準確。因此，有必要開發(fā)一種能夠精確預測儲層裂縫位置的方法及系統(tǒng)。

公開于本發(fā)明背景技術部分的信息僅僅旨在加深對本發(fā)明的一般背景技術的理解，而不應當被視為承認或以任何形式暗示該信息構成已為本領域技術人員所公知的現(xiàn)有技術。

發(fā)明內(nèi)容

本發(fā)明提出了一種儲層裂縫預測方法及系統(tǒng)，所述方法能夠基于離散裂縫網(wǎng)格模型及卡曼增益進行儲層裂縫預測，從而更精確地預測儲層裂縫的信息。

根據(jù)本發(fā)明的一方面，提出了一種儲層裂縫預測方法，所述方法可以包括：a)、在工區(qū)范圍內(nèi)，隨機生成多個裂縫，獲得多個裂縫的特征參數(shù)；b)、基于所述多個裂縫的特征參數(shù)，建立離散裂縫網(wǎng)格模型；c)、基于所述離散裂縫網(wǎng)格模型，進行生產(chǎn)歷史數(shù)值模擬，獲得生產(chǎn)動態(tài)參數(shù)模擬值；d)、基于所述生產(chǎn)動態(tài)參數(shù)模擬值，獲取所述生產(chǎn)動態(tài)參數(shù)模擬值和實際的生產(chǎn)動態(tài)參數(shù)的協(xié)方差矩陣，進而獲得卡曼增益；e)、基于所述卡曼增益，對所述多個裂縫的特征參數(shù)進行調(diào)整并判斷收斂性；f)、在所述多個裂縫的特征參數(shù)不收斂的情況下，重復步驟b)～e)；以及g)、在所述多個裂縫的特征參數(shù)收斂的情況下，將收斂后的特征參數(shù)作為儲層裂縫的特征參數(shù)。

根據(jù)本發(fā)明的另一方面，提出了一種儲層裂縫預測系統(tǒng)，所述系統(tǒng)可以包括：用于在工區(qū)范圍內(nèi)，隨機生成多個裂縫，獲得所述多個裂縫的特征參數(shù)的單元A；用于基于所述多個裂縫的特征參數(shù)，建立離散裂縫網(wǎng)格模型的單元B；用于基于所述離散裂縫網(wǎng)格模型，進行生產(chǎn)歷史數(shù)值模擬，獲得生產(chǎn)動態(tài)參數(shù)模擬值的單元C；用于基于所述生產(chǎn)動態(tài)參數(shù)模擬值，獲取所述生產(chǎn)動態(tài)參數(shù)模擬值和實際的生產(chǎn)動態(tài)參數(shù)的協(xié)方差矩陣，進而獲得卡曼增益的單元D；用于基于所述卡曼增益，對所述多個裂縫的特征參數(shù)進行調(diào)整并判斷收斂性的單元E；以及用于在所述多個裂縫的特征參數(shù)收斂的情況下，將收斂后的特征參數(shù)作為儲層裂縫的特征參數(shù)的單元G。

本發(fā)明的方法和裝置具有其它的特性和優(yōu)點，這些特性和優(yōu)點從并入本文中的附圖和隨后的具體實施方式中將是顯而易見的，或者將在并入本文中的附圖和隨后的具體實施方式中進行詳細陳述，這些附圖和具體實施方式共同用于解釋本發(fā)明的特定原理。

附圖說明

通過結合附圖對本發(fā)明示例性實施方式進行更詳細的描述，本發(fā)明的上述以及其它目的、特征和優(yōu)勢將變得更加明顯，其中，在本發(fā)明示例性實施方式中，相同的參考標號通常代表相同部件。

圖1示出了根據(jù)本發(fā)明的儲層裂縫預測方法的步驟的流程圖。

圖2示出了根據(jù)本發(fā)明的一個實施方式的工區(qū)范圍及注采井位示意圖。

圖3a至圖3i分別示出了根據(jù)本發(fā)明的一個實施方式的各井的生產(chǎn)數(shù)據(jù)的示意圖。

圖4示出了根據(jù)本發(fā)明的一個實施方式的所有初始實現(xiàn)中的裂縫的示意圖。

圖5a至圖5d分別示出了根據(jù)本發(fā)明的一個實施方式的卡曼增益調(diào)整后裂縫特征參數(shù)的收斂情況的示意圖。

圖6a至圖6f分別示出了根據(jù)本發(fā)明的一個實施方式的卡曼增益調(diào)整后裂縫位置的收斂情況的示意圖。

圖7a至圖7j分別示出了根據(jù)本發(fā)明的一個實施方式的通過儲層裂縫預測方法預測的生產(chǎn)動態(tài)的示意圖。

具體實施方式