。
[0124] 通過對原始譜系數(shù)優(yōu)化之后得到了優(yōu)化譜系數(shù)，再通過優(yōu)化譜系數(shù)計算核磁束縛 水飽和度，這是核磁束縛水飽和度獲取單元所要做的工作。具體的，核磁束縛水飽和度獲取 單元404按照上述公式（4)來計算核磁束縛水飽和度，即：
[0125]
[0126] 其中，NSbvi為核磁束縛水飽和度；
[0127] 第j個反演布點對應(yīng)的孔隙度的幅度值。
[0128] 按照公式（4)計算得到的是不同的截止值T2rt下的巖心的核磁束縛水飽和度，即 多個核磁束縛水飽和度，然后需要最優(yōu)核磁束縛水飽和度確定單元405來確定最優(yōu)核磁束 縛水飽和度，這個最優(yōu)核磁束縛水飽和度即最終所需要的。
[0129] 下面通過對比說明本發(fā)明方法和裝置的優(yōu)越性。圖5為使用原始譜系數(shù)計算所得 束縛水飽和度與實驗測得的束縛水飽和度進(jìn)行比較，使用原始譜系數(shù)計算所得束縛水飽和 度與實驗測得的束縛水飽和度兩者之間的相關(guān)性為〇. 661 (1表示相關(guān)性）。圖6為使用優(yōu) 化譜系數(shù)計算所得束縛水飽和度與實驗測得的束縛水飽和度的比較結(jié)果圖，使用優(yōu)化譜系 數(shù)計算所得束縛水飽和度與實驗測得的束縛水飽和度兩者之間的相關(guān)性為〇. 912。由此可 知，相對于傳統(tǒng)的譜系數(shù)計算方法，使用本發(fā)明方法優(yōu)化譜系數(shù)計算束縛水飽和度的計算 精度提高了 25 %左右。
[0130] 下面利用7塊巖心來驗證本發(fā)明方法，也稱為盲點檢測。如表2所示，其中，第二 列為其他核磁方法獲得的束縛水飽和度BVI-NMR，第三列為實驗測量巖心的束縛水飽和度 數(shù)據(jù)BVI-C〇re(Sbvi);第四列為使用本發(fā)明方法獲得的核磁束縛水飽和度NSbvi;第五列為 實驗測量巖心的束縛水飽和度數(shù)據(jù)BVI-Core(Sbvi)與使用本發(fā)明方法獲得的核磁束縛水飽 和度NSbvi之間的誤差。圖7為使用本發(fā)明優(yōu)化譜系數(shù)計算所得束縛水飽和度與實驗測得 的束縛水飽和度進(jìn)行比較，可知，核磁計算束縛水飽和度與實驗測試束縛水飽和度之間的 相關(guān)系數(shù)是0. 895,平均相對誤差為12. 45%，。相對于原始譜系數(shù)的計算，極大的提高了核 磁束縛水飽和度的計算精度。
[0131] 本次發(fā)明的盲點檢測結(jié)果，該發(fā)明方法對于巖石核磁分析數(shù)據(jù)具有較好的適用性 和較好的計算精度。
[0132] 表 2
[0133]

[0134] 顯然，本領(lǐng)項的扠木人員奴該明曰，上還的本友明買施例的谷悮塊現(xiàn)各步驟可以 用通用的計算裝置來實現(xiàn)，它們可以集中在單個的計算裝置上，或者分布在多個計算裝置 所組成的網(wǎng)絡(luò)上，可選地，它們可以用計算裝置可執(zhí)行的程序代碼來實現(xiàn)，從而，可以將它 們存儲在存儲裝置中由計算裝置來執(zhí)行，并且在某些情況下，可以以不同于此處的順序執(zhí) 行所示出或描述的步驟，或者將它們分別制作成各個集成電路模塊，或者將它們中的多個 模塊或步驟制作成單個集成電路模塊來實現(xiàn)。這樣，本發(fā)明實施例不限制于任何特定的硬 件和軟件結(jié)合。
[0135] 以上所述僅為本發(fā)明的優(yōu)選實施例而已，并不用于限制本發(fā)明，對于本領(lǐng)域的技 術(shù)人員來說，本發(fā)明實施例可以有各種更改和變化。凡在本發(fā)明的精神和原則之內(nèi)，所作的 任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。
【主權(quán)項】
1. 一種確定束縛水飽和度的方法，其特征在于，包括： 反演巖心的核磁弛豫回波數(shù)據(jù)，得到核磁共振橫向弛豫時間t2譜； 獲得巖心的束縛水飽和度數(shù)據(jù)； 根據(jù)1~2譜，選擇不同的反演布點的橫向弛豫時間T2作為橫向弛豫時間1~2的截止值T2c;t， 計算不同截止值T2rt下的優(yōu)化譜系數(shù)； 根據(jù)不同截止值T2et下的優(yōu)化譜系數(shù)，計算不同截止值T2et下的巖心的核磁束縛水飽 和度； 分別計算巖心的束縛水飽和度與不同截止值T2rt下的巖心的核磁束縛水飽和度的誤 差，將計算出的最小誤差相對應(yīng)的優(yōu)化譜系數(shù)確定為最優(yōu)譜系數(shù)，最優(yōu)譜系數(shù)相對應(yīng)的核 磁束縛水飽和度確定為最優(yōu)核磁束縛水飽和度。2. 如權(quán)利要求1所述的方法，其特征在于，所述根據(jù)1~2譜，選擇不同的反演布點的橫向 弛豫時間T2作為橫向弛豫時間T2的截止值T2et，計算不同截止值T2rt下的優(yōu)化譜系數(shù)，按照 如下公式進(jìn)行計算：其中，Cbv_為第j個反演布點的橫向弛豫時間對應(yīng)的優(yōu)化譜系數(shù)，無量綱； T2j為第j個反演布點的橫向弛豫時間，ms; T2rt為橫向弛豫時間T2的截止值，ms; a，b是擬合系數(shù)。3. 如權(quán)利要求2所述的方法，其特征在于，所述擬合系數(shù)a，b按照如下方式獲得： 反演多塊巖心的核磁弛豫回波數(shù)據(jù)，得到多塊巖心的核磁共振橫向弛豫時間T2譜； 根據(jù)1~2譜，計算多塊巖心的橫向弛豫時間1~2的幾何平均值； 獲得多塊巖心的束縛水飽和度數(shù)據(jù)； 將所得的多塊巖心的橫向弛豫時間T2的幾何平均值和多塊巖心的束縛水飽和度數(shù)據(jù) 進(jìn)行擬合，獲得擬合系數(shù)a，b。4. 如權(quán)利要求3所述的方法，其特征在于，所述將所得的多塊巖心的橫向弛豫時間T2 的幾何平均值和多塊巖心的束縛水飽和度數(shù)據(jù)進(jìn)行擬合，按照如下公式進(jìn)行計算：其中，Sbvi是巖心的束縛水飽和度； 丁^是巖心的橫向弛豫時間T2的幾何平均值，ms。5. 如權(quán)利要求2所述的方法，其特征在于，所述根據(jù)不同截止值T2c;t下的優(yōu)化譜系數(shù)， 計算不同截止值T2rt下的巖心的核磁束縛水飽和度，按照如下公式計算；其中，NSbvi為核磁束縛水飽和度； Aj為第j個反演布點對應(yīng)的孔隙度的幅度值。6. -種確定束縛水飽和度的裝置，其特征在于，包括： T2譜獲取單元，用于反演巖心的核磁弛豫回波數(shù)據(jù)，得到核磁共振橫向弛豫時間T2譜； 束縛水飽和度獲取單元，用于獲得巖心的束縛水飽和度數(shù)據(jù)； 優(yōu)化譜系數(shù)獲取單元，用于根據(jù)1~2譜，選擇不同的反演布點的橫向弛豫時間T2作為橫 向弛豫時間1~2的截止值T2rt，計算不同截止值T2rt下的優(yōu)化譜系數(shù)； 核磁束縛水飽和度獲取單元，用于根據(jù)不同截止值T2et下的優(yōu)化譜系數(shù)，計算不同截止 值T2rt下的巖心的核磁束縛水飽和度； 最優(yōu)核磁束縛水飽和度確定單元，用于分別計算巖心的束縛水飽和度與不同截止值T2c;t下的巖心的核磁束縛水飽和度的誤差，將計算出的最小誤差相對應(yīng)的優(yōu)化譜系數(shù)確定 為最優(yōu)譜系數(shù)，將最優(yōu)譜系數(shù)相對應(yīng)的核磁束縛水飽和度確定為最優(yōu)核磁束縛水飽和度。7. 如權(quán)利要求6所述的裝置，其特征在于，所述優(yōu)化譜系數(shù)獲取單元具體用于按如下 公式計算不同截止值T2c;t下的優(yōu)化譜系數(shù)：其中，Cbv_為第j個反演布點的橫向弛豫時間對應(yīng)的優(yōu)化譜系數(shù)，無量綱； T2j為第j個反演布點的橫向弛豫時間，ms; T2rt為橫向弛豫時間T2的截止值，ms;a，b是擬合系數(shù)。8. 如權(quán)利要求7所述的裝置，其特征在于，所述優(yōu)化譜系數(shù)獲取單元具體用于按如下 方式獲得所述擬合系數(shù)a，b: 反演多塊巖心的核磁弛豫回波數(shù)據(jù)，得到多塊巖心的核磁共振橫向弛豫時間T2譜； 根據(jù)1~2譜，計算多塊巖心的橫向弛豫時間1~2的幾何平均值； 獲得多塊巖心的束縛水飽和度數(shù)據(jù)； 將所得的多塊巖心的橫向弛豫時間T2的幾何平均值和多塊巖心的束縛水飽和度數(shù)據(jù) 進(jìn)行擬合，獲得擬合系數(shù)a，b。9. 如權(quán)利要求8所述的裝置，其特征在于，所述優(yōu)化譜系數(shù)獲取單元具體用于按如下 公式將所得的多塊巖心的橫向弛豫時間T2的幾何平均值和多塊巖心的束縛水飽和度數(shù)據(jù) 進(jìn)行擬合：其中，sbvi是巖心的束縛水飽和度； 丁^是巖心的橫向弛豫時間T2的幾何平均值，ms。10.如權(quán)利要求7所述的裝置，其特征在于，所述核磁束縛水飽和度獲取單元具體用于 按如下公式計算不同截止值T2c;t下的巖心的核磁束縛水飽和度：其中，NSbvi為核磁束縛水飽和度； Aj為第j個反演布點對應(yīng)的孔隙度的幅度值。
【專利摘要】本發(fā)明提供一種確定束縛水飽和度的方法及裝置，包括：測量并反演巖心的核磁弛豫回波數(shù)據(jù)，得到核磁共振橫向弛豫時間T2譜；測量巖心的束縛水飽和度數(shù)據(jù)；根據(jù)T2譜，選擇不同的反演布點的橫向弛豫時間T2作為T2ct，計算不同T2ct下的優(yōu)化譜系數(shù)；根據(jù)不同T2ct下的優(yōu)化譜系數(shù)，計算不同T2ct下的巖心的核磁束縛水飽和度；分別計算巖心的束縛水飽和度與不同T2ct下的巖心的核磁束縛水飽和度的誤差，當(dāng)誤差取得極小值時，相對應(yīng)的優(yōu)化譜系數(shù)為最優(yōu)譜系數(shù)，相對應(yīng)的核磁束縛水飽和度為最優(yōu)核磁束縛水飽和度。本發(fā)明通過核磁共振方法得到最優(yōu)的譜系數(shù)，根據(jù)最優(yōu)的譜系數(shù)計算束縛水飽和度，顯著提高了束縛水飽和度的計算精度。
【IPC分類】G06F17/11, G01N15/08
【公開號】CN104990854
【申請?zhí)枴緾N201510391658
【發(fā)明人】胡法龍, 劉堂晏, 李潮流, 王飛飛, 王昌學(xué), 李長喜
【申請人】中國石油天然氣股份有限公司
【公開日】2015年10月21日
【申請日】2015年7月6日