[發明專利]基于時間推移的微元化導電模型求水淹層相關參數的方法

在審

• 申請號: 201710704339.3
• 申請日: 2017-08-17
• 公開（公告）號: CN107741605A
• 公開（公告）日: 2018-02-27
• 發明（設計）人: 丁磊;雷霄;譚偉;張恒榮;王利娟
• 申請（專利權）人: 中國海洋石油總公司;中海石油（中國）有限公司湛江分公司
• 主分類號: G01V3/18
• 分類號: G01V3/18;G01N15/08
• 代理公司: 廣州市南鋒專利事務所有限公司44228
• 代理人: 李慧

技術領域

本發明涉及地球物理測井技術領域，尤其涉及基于時間推移的微元化導電模型求水淹層相關參數的方法。

背景技術

目前水淹層的水淹級別主要是根據產水率來劃分的，要想求準產水率，準確計算出水淹層含水飽和度是關鍵。油層水淹后由于注入水礦化度與原生水礦化度不同，導致混合水電阻率不斷變化，使得含水飽和度難以求準。目前水淹層混合水電阻率求取方法主要有根據自然電位曲線計算以及通過理論推導或者數值模擬方法建模兩種，但是這些方法尚存在如下技術問題：

1、受限于測井系列的影響，大量未測量自然電位曲線的水淹井無法應用自然電位法計算混合水電阻率；

2、絕大部分的理論推導或者數值模擬方法都是以未水淹時為起始點，而未考慮到水淹其實是一個長期的過程，各相關參數在水淹過程中會發生實時的變化。

發明內容

本發明提供了一種基于時間推移的微元化導電模型求水淹層相關參數的方法，首先通過巖心分析相關實驗建立束縛水、殘余油、相對滲透率及產水率相關模型，然后采用基于時間推移的微元化導電模型迭代聯合求解出水淹層各關鍵參數，提高了水淹層解釋的精度。

為解決上述技術問題，本發明提供了一種基于時間推移的微元化導電模型求水淹層相關參數的方法，包括以下步驟：

S1：通過巖電實驗數據，計算出區塊巖電參數a、b、m、n；

S2：通過毛管壓力試驗、離心機實驗及油水相滲實驗建立區塊束縛水、殘余油和油水相對滲透率計算模型；

S3：通過巖性曲線(GR、SP曲線)計算出泥質含量Vsh，通過物性曲線：DEN、CNL和AC曲線計算出孔隙度φ，通過電性曲線讀出泥巖電阻率Rsh；

S4：令含水飽和度由Swi到1-Sor變化，步長ΔSw＝0.1％；

S5：含水飽和度取某一個定值時，令未參與混合的原生水飽和度Swb由Swi到0變化，步長ΔSwb＝0.1％；

S6：計算出此時的產水率Fw，再根據產水率計算出產出水礦化度為：

根據Pz₂計算出混合水電阻率Rw2，然后再根據并聯導電原理，將混合水電阻率與未參與混合的原生地層水并聯，得到視地層水電阻率：

再由印尼公式：

計算出新的Sw2’，算出含水飽和度誤差Δ＝|Sw₂'-Sw₂|。

S7：循環Swb、Sw，采用時間推移微元化思想，以上一時刻的混合液電阻率和含水飽和度作為起點，計算新的混合液電阻率和含水飽和度，找到令Δ最小的Sw以及Swb作為真實值，同時還可以得到混合水電阻率Rwz以及產水率Fw。

本申請實施例中提供的一個或多個技術方案，至少具有如下技術效果或優點：

本發明通過時間推移的思想，建立了能更好的反映水淹層客觀規律的混合液礦化度變化規律模型，經循環一次性求解出混合水電阻率Rwz、含水飽和度Sw以及產水率Fw等多個水淹層關鍵參數。解決了水淹層含水飽和度難以求準的問題，提高了水淹油田開發經濟效益。

附圖說明

圖1方法處理效果圖。

具體實施方式

為了更好的理解上述技術方案，下面以具體的實施方式對上述技術方案進行詳細的說明。

本實施例所述的一種基于時間推移的微元化導電模型求水淹層相關參數的方法，S1：通過巖電實驗數據，計算出區塊巖電參數a、b、m、n；

S2：通過毛管壓力試驗、離心機實驗及油水相滲實驗建立區塊束縛水、殘余油和油水相對滲透率計算模型；

S3：通過巖性曲線(GR、SP曲線)計算出泥質含量Vsh，通過物性曲線：DEN、CNL和AC曲線計算出孔隙度φ，通過電性曲線讀出泥巖電阻率Rsh；

S4：令含水飽和度由Swi到1-Sor變化，步長ΔSw＝0.1％；

S5：含水飽和度取某一個定值時，令未參與混合的原生水飽和度Swb由Swi到0變化，步長ΔSwb＝0.1％；

S6：計算出此時的產水率Fw，再根據產水率計算出產出水礦化度為：

根據Pz₂計算出混合水電阻率Rw2，然后再根據并聯導電原理，將混合水電阻率與未參與混合的原生地層水并聯，得到視地層水電阻率：

再由印尼公式：

計算出新的Sw2’，算出含水飽和度誤差Δ＝|Sw₂'-Sw₂|。