[發明專利]一種評價儲層孔隙結構的方法及設備

有效

• 申請號: 201910047194.3
• 申請日: 2019-01-18
• 公開（公告）號: CN109736793B
• 公開（公告）日: 2021-01-29
• 發明（設計）人: 謝然紅;吳勃翰;金國文;高倫;谷明宣
• 申請（專利權）人: 中國石油大學（北京）
• 主分類號: E21B49/00
• 分類號: E21B49/00;G06F30/20
• 代理公司: 北京同立鈞成知識產權代理有限公司 11205
• 代理人: 朱穎;黃健

本發明實施例提供一種評價儲層孔隙結構的方法及設備，該方法通過獲取目標儲層的核磁共振回波數據，對核磁共振回波數據進行時域變換，并計算變換后的核磁共振回波數據的振幅的平方及頻率，對振幅的平方和頻率進行線性擬合，得到擬合直線，獲取擬合直線對應的斜率，并根據斜率計算核磁共振回波數據的分形維數，根據分形維數評價目標儲層的孔隙結構，自動根據核磁共振回波數據進行分形分析，得到相應的分形維數，無需進行核磁共振回波數據反演，計算量較小，減少計算誤差。

技術領域

本發明實施例涉及油氣田勘探開發技術領域，尤其涉及一種評價儲層孔隙結構的方法及設備。

背景技術

在油氣開發技術領域，巖石的孔隙結構是指巖石內的孔隙和喉道類型、大小、分布及其相互連通關系。研究油氣儲層的巖石孔隙結構，深入揭示油氣儲層的內部結構，對油氣田勘探和開發有著重要的意義。

目前，為了評價油氣儲層的孔隙結構，獲取油氣儲層的孔隙結構參數，一般采用的是核磁共振測井。當需要評價儲層孔隙結構時，通常對核磁共振測井采集到的核磁共振回波數據進行反演得到核磁共振T₂分布，然后計算核磁共振T₂分布的分形維數來評價儲層孔隙結構。

然而，發明人發現現有技術中至少存在如下問題：對于比常規儲層更加復雜的致密砂巖儲層，進行核磁共振測井時，在低信噪比下核磁共振回波數據的反演會存在較大的誤差，導致計算出來的分形維數誤差較大。

發明內容

本發明實施例提供一種評價儲層孔隙結構的方法及設備，以減少計算分形維數的計算誤差。

第一方面，本發明實施例提供一種評價儲層孔隙結構的方法，包括：

獲取目標儲層的核磁共振回波數據；

對所述核磁共振回波數據進行時域變換，并計算變換后的核磁共振回波數據的振幅的平方及頻率；

對所述振幅的平方和頻率進行線性擬合，得到擬合直線；

獲取所述擬合直線對應的斜率，并根據所述斜率計算所述核磁共振回波數據的分形維數；

根據所述分形維數評價所述目標儲層的孔隙結構。

在一種可能的設計中，所述根據所述斜率計算所述核磁共振回波數據的分形維數的計算公式為：

其中，D_echoe為所述分形維數，k為擬合直線斜率的絕對值。

在一種可能的設計中，所述對所述核磁共振回波數據進行時域變換之前，還包括：

判斷目標儲層核磁共振回波數據對應的概率分布是否符合高斯分布，其中目標儲層核磁共振回波數據為所述核磁共振回波數據中的任一個；

若所述目標儲層核磁共振回波數據的概率分布不符合高斯分布，則對所述目標儲層核磁共振回波數據進行正態變換。

在一種可能的設計中，所述對所述核磁共振回波數據進行時域變換，包括：

按照預設的傅立葉變換公式，對所述核磁共振回波數據進行傅立葉變換；

在一種可能的設計中，所述計算變換后的核磁共振回波數據的振幅的平方及頻率，包括：

獲取核磁共振回波數據的回波個數以及所述傅立葉變換公式中的實部和虛部；

根據所述核磁共振回波個數和所述傅立葉變換公式中的實部與虛部，計算所述變換后的核磁共振回波數據的振幅的平方；

獲取所述核磁共振回波數據的回波間隔；

根據所述回波間隔和所述核磁共振回波個數計算所述核磁共振回波數據中的采樣點的頻率。