本發(fā)明涉及油氣勘探與開發(fā)相關，尤其是涉及一種儲層流體識別方法、系統(tǒng)、裝置和可讀存儲介質。

背景技術：

1、常用的儲層流體識別方法包括地震、測井和電磁法等，地震數(shù)據(jù)通常用于提供高分辨率的構造與層序，用于流體識別的方法包括疊后屬性交匯分析、疊前振幅隨偏移距的變化（amplitude?variation?with?offset，avo）反演方法，但對含油水儲層的區(qū)分力度不夠；測井縱向分辨率高，但僅一孔之見，覆蓋范圍?。浑姶趴碧椒▽λw變化非常靈敏，但分辨率受各種因素影響比較低。因此采用單一方法，多解性強，儲層流體識別精度不夠。

2、目前，現(xiàn)有技術中已存在對地震電阻率、電磁電阻率和測井電阻率進行融合來進行儲層流體識別的方式，但將地震電阻率、電磁電阻率和測井電阻率融合的方法采用的是主成分分析的方法，但主成分分析法會將三個維度的數(shù)據(jù)進行簡化，造成信息的損失，使得數(shù)據(jù)的可信度降低，同時主成分分析只適用于線性數(shù)據(jù)，而三者電阻率之間存在的是非線性關系，融合之后的效果不好，不能完全捕捉到數(shù)據(jù)的特征，使得得到的儲層流體識別準確度不夠。

技術實現(xiàn)思路

1、本發(fā)明旨在至少解決現(xiàn)有技術中存在的技術問題之一。為此，本發(fā)明提出一種儲層流體識別方法，能夠提高儲層流體識別的準確度。

2、本發(fā)明還提供了一種儲層流體識別系統(tǒng)、用于執(zhí)行上述儲層流體識別方法的控制裝置以及計算機可讀存儲介質。

3、根據(jù)本發(fā)明的第一方面實施例的儲層流體識別方法，所述方法包括：

4、獲取地震數(shù)據(jù)、測井數(shù)據(jù)和電磁數(shù)據(jù)；

5、根據(jù)所述地震數(shù)據(jù)、所述測井數(shù)據(jù)和所述電磁數(shù)據(jù)分別計算得到地震視電阻率數(shù)據(jù)、測井視電阻率數(shù)據(jù)和電磁視電阻率數(shù)據(jù)；

6、對所述地震視電阻率數(shù)據(jù)、所述測井視電阻率數(shù)據(jù)和所述電磁視電阻率數(shù)據(jù)進行歸一化處理，得到深度域的地震轉換數(shù)據(jù)、測井轉換數(shù)據(jù)和電磁轉換數(shù)據(jù)；

7、利用小波變換將所述地震轉換數(shù)據(jù)、所述測井轉換數(shù)據(jù)和所述電磁轉換數(shù)據(jù)重構得到復視電阻率數(shù)據(jù)；

8、基于所述復視電阻率數(shù)據(jù)進行儲層流體識別。

9、根據(jù)本發(fā)明實施例的儲層流體識別方法，至少具有如下有益效果：

10、本發(fā)明實施例通過將地震視電阻率數(shù)據(jù)、測井視電阻率數(shù)據(jù)和電磁視電阻率數(shù)據(jù)歸一化為深度域的地震轉換數(shù)據(jù)、測井轉換數(shù)據(jù)和電磁轉換數(shù)據(jù)，并利用小波變換將地震轉換數(shù)據(jù)、測井轉換數(shù)據(jù)和電磁轉換數(shù)據(jù)進行重構得到復視電阻率數(shù)據(jù)。小波變換是一種全局變換，在小波域具有很好的定位能力，可以聚焦到被處理數(shù)據(jù)的任何細節(jié)，因此本發(fā)明實施例運用小波變換融合之后的復視電阻率數(shù)據(jù)分辨率很高，能夠提高儲層流體識別的準確度。

11、根據(jù)本發(fā)明的一些實施例，所述利用小波變換將所述地震轉換數(shù)據(jù)、所述測井轉換數(shù)據(jù)和所述電磁轉換數(shù)據(jù)重構得到復視電阻率數(shù)據(jù)，包括：

12、利用小波變換將所述地震轉換數(shù)據(jù)、所述測井轉換數(shù)據(jù)和所述電磁轉換數(shù)據(jù)變換到相同尺度范圍和相同頻率范圍下，得到地震小波數(shù)據(jù)、測井小波數(shù)據(jù)和電磁小波數(shù)據(jù)；

13、根據(jù)所述地震小波數(shù)據(jù)、所述測井小波數(shù)據(jù)和所述電磁小波數(shù)據(jù)確定不同尺度、不同頻率下對應的多個重構系數(shù)；

14、根據(jù)多個所述重構系數(shù)對所述地震小波數(shù)據(jù)、所述測井小波數(shù)據(jù)和所述電磁小波數(shù)據(jù)進行小波逆變換，重構出所述復視電阻率數(shù)據(jù)。

15、根據(jù)本發(fā)明的一些實施例，所述根據(jù)所述地震小波數(shù)據(jù)、所述測井小波數(shù)據(jù)和所述電磁小波數(shù)據(jù)確定不同尺度、不同頻率下對應的多個重構系數(shù)，包括：

16、計算同一尺度、同一頻率下的所述地震小波數(shù)據(jù)、所述測井小波數(shù)據(jù)和所述電磁小波數(shù)據(jù)的方差數(shù)據(jù)；

17、根據(jù)同一尺度下所有頻率對應的所有所述方差數(shù)據(jù)計算該尺度下對應頻率的各個方差權重，得到不同尺度、不同頻率下對應的多個所述方差權重；

18、將多個所述方差權重作為多個所述重構系數(shù)。

19、根據(jù)本發(fā)明的一些實施例，所述方差權重通過以下步驟得到：

20、計算同一尺度下所有頻率對應的所有方差數(shù)據(jù)之間兩兩乘積之和，得到方差積；

21、將該尺度下除對應頻率之外剩余頻率的所有方差數(shù)據(jù)的乘積，在所述方差積中的占比確定為所述方差權重。

22、根據(jù)本發(fā)明的一些實施例，所述地震視電阻率數(shù)據(jù)通過以下步驟得到：

23、基于地震波阻抗和電阻率之間的巖石物理關系，根據(jù)所述測井數(shù)據(jù)和所述地震數(shù)據(jù)計算得到所述地震視電阻率數(shù)據(jù)，其中，所述巖石物理關系根據(jù)faust公式和gardner公式推導得到。

24、根據(jù)本發(fā)明的一些實施例，所述測井數(shù)據(jù)包括深度數(shù)據(jù)、電阻率曲線、中子曲線、聲波曲線和密度曲線，所述基于地震波阻抗和測井電阻率之間的巖石物理關系，根據(jù)所述測井數(shù)據(jù)和所述地震數(shù)據(jù)計算得到所述地震視電阻率數(shù)據(jù)，包括：

25、根據(jù)所述地震數(shù)據(jù)通過疊后地震反演得到地震波阻抗；

26、根據(jù)所述深度數(shù)據(jù)、所述電阻率曲線、所述中子曲線、所述聲波曲線、所述密度曲線和所述地震波阻抗計算得到所述地震視電阻率數(shù)據(jù)。

27、根據(jù)本發(fā)明的一些實施例，所述對所述地震視電阻率數(shù)據(jù)、所述測井視電阻率數(shù)據(jù)和所述電磁視電阻率數(shù)據(jù)進行歸一化處理，得到深度域的地震轉換數(shù)據(jù)、測井轉換數(shù)據(jù)和電磁轉換數(shù)據(jù)，包括：

28、將所述地震視電阻率數(shù)據(jù)、所述測井視電阻率數(shù)據(jù)和所述電磁視電阻率數(shù)據(jù)配準到平面坐標的校準和縱向上的統(tǒng)一，得到深度域的第一地震數(shù)據(jù)、第一測井數(shù)據(jù)和第一電磁數(shù)據(jù)；

29、統(tǒng)一所述第一地震數(shù)據(jù)、所述第一測井數(shù)據(jù)和所述第一電磁數(shù)據(jù)在所述平面坐標上的處理面元大小和縱向上的采樣間隔大小，得到第二地震數(shù)據(jù)、第二測井數(shù)據(jù)和第二電磁數(shù)據(jù)；

30、將所述第二地震數(shù)據(jù)、所述第二測井數(shù)據(jù)和所述第二電磁數(shù)據(jù)歸一化到統(tǒng)一的數(shù)據(jù)范圍內，得到所述地震轉換數(shù)據(jù)、所述測井轉換數(shù)據(jù)和所述電磁轉換數(shù)據(jù)。

31、根據(jù)本發(fā)明的第二方面實施例的儲層流體識別系統(tǒng)，所述系統(tǒng)包括：

32、數(shù)據(jù)獲取模塊，用于獲取地震數(shù)據(jù)、測井數(shù)據(jù)和電磁數(shù)據(jù)；

33、視電阻率確定模塊，用于根據(jù)所述地震數(shù)據(jù)、所述測井數(shù)據(jù)和所述電磁數(shù)據(jù)分別計算得到地震視電阻率數(shù)據(jù)、測井視電阻率數(shù)據(jù)和電磁視電阻率數(shù)據(jù)；

34、歸一化處理模塊，用于對所述地震視電阻率數(shù)據(jù)、所述測井視電阻率數(shù)據(jù)和所述電磁視電阻率數(shù)據(jù)進行歸一化處理，得到深度域的地震轉換數(shù)據(jù)、測井轉換數(shù)據(jù)和電磁轉換數(shù)據(jù)；

35、小波變換模塊，用于利用小波變換將所述地震轉換數(shù)據(jù)、所述測井轉換數(shù)據(jù)和所述電磁轉換數(shù)據(jù)重構得到復視電阻率數(shù)據(jù)；

36、儲層流體識別模塊，用于基于所述復視電阻率數(shù)據(jù)進行儲層流體識別。

37、根據(jù)本發(fā)明實施例的儲層流體識別系統(tǒng)，至少具有如下有益效果：

38、本發(fā)明實施例通過將地震視電阻率數(shù)據(jù)、測井視電阻率數(shù)據(jù)和電磁視電阻率數(shù)據(jù)歸一化為深度域的地震轉換數(shù)據(jù)、測井轉換數(shù)據(jù)和電磁轉換數(shù)據(jù)，并利用小波變換將地震轉換數(shù)據(jù)、測井轉換數(shù)據(jù)和電磁轉換數(shù)據(jù)進行重構得到復視電阻率數(shù)據(jù)。小波變換是一種全局變換，在小波域具有很好的定位能力，可以聚焦到被處理數(shù)據(jù)的任何細節(jié)，因此本發(fā)明實施例運用小波變換融合之后的復視電阻率數(shù)據(jù)分辨率很高，能夠提高儲層流體識別的準確度。

39、根據(jù)本發(fā)明的第三方面實施例的控制裝置，包括存儲器、處理器及存儲在存儲器上并可在處理器上運行的計算機程序，所述處理器執(zhí)行所述計算機程序時實現(xiàn)如上述第一方面實施例所述的儲層流體識別方法。由于控制裝置采用了上述實施例的儲層流體識別方法的全部技術方案，因此至少具有上述實施例的技術方案所帶來的所有有益效果。

40、根據(jù)本發(fā)明的第四方面實施例的計算機可讀存儲介質，存儲有計算機可執(zhí)行指令，所述計算機可執(zhí)行指令用于執(zhí)行如上述第一方面實施例所述的儲層流體識別方法。由于計算機可讀存儲介質采用了上述實施例的儲層流體識別方法的全部技術方案，因此至少具有上述實施例的技術方案所帶來的所有有益效果。

41、本發(fā)明的其它特征和優(yōu)點將在隨后的說明書中闡述，并且，部分地從說明書中變得顯而易見，或者通過實施本發(fā)明而了解。