確定束縛水飽和度的方法及裝置的制造方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及巖石滲流特征研宄技術(shù)領(lǐng)域，特別涉及一種確定束縛水飽和度的方法 及裝置。
【背景技術(shù)】
[0002] 束縛水飽和度是表征儲(chǔ)層物性最重要的參數(shù)之一，在砂泥巖儲(chǔ)層評(píng)價(jià)中起著關(guān)鍵 作用。隨著儲(chǔ)層泥質(zhì)含量越高，滲透率變化越大，孔隙結(jié)構(gòu)越復(fù)雜，其作用越顯重要，確定難 度越大。譜系數(shù)方法是計(jì)算巖石束縛水飽和度的一種方法，用于研宄巖石的滲流特征。現(xiàn) 有的譜系數(shù)方法，只是簡(jiǎn)單利用數(shù)據(jù)擬合的方法，得到譜系數(shù)，然后，根據(jù)譜系數(shù)計(jì)算巖石 的束縛水飽和度。這個(gè)方法的問題是并不能得到具有較好精度的譜系數(shù)，即無法保證得到 最優(yōu)的譜系數(shù)來提高束縛水飽和度的計(jì)算精度。

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

[0003] 本發(fā)明實(shí)施例提供了一種確定束縛水飽和度的方法，通過優(yōu)化方法獲得最優(yōu)譜系 數(shù)，從而得到最優(yōu)核磁束縛水飽和度，提高了核磁束縛水飽和度的計(jì)算精度。該方法包括：
[0004] 反演巖心的核磁弛豫回波數(shù)據(jù)，得到核磁共振橫向弛豫時(shí)間t2譜；
[0005] 獲得巖心的束縛水飽和度數(shù)據(jù)；
[0006] 根據(jù)1~2譜，選擇不同的反演布點(diǎn)的橫向弛豫時(shí)間T2作為橫向弛豫時(shí)間T2的截止 值T2rt，計(jì)算不同截止值T2c;t下的優(yōu)化譜系數(shù)；
[0007] 根據(jù)不同截止值T2et下的優(yōu)化譜系數(shù)，計(jì)算不同截止值T2et下的巖心的核磁束縛 水飽和度；
[0008] 分別計(jì)算巖心的束縛水飽和度與不同截止值T2et下的巖心的核磁束縛水飽和度的 誤差，將計(jì)算出的最小誤差相對(duì)應(yīng)的優(yōu)化譜系數(shù)確定為最優(yōu)譜系數(shù)，最優(yōu)譜系數(shù)相對(duì)應(yīng)的 核磁束縛水飽和度確定為最優(yōu)核磁束縛水飽和度。
[0009] 在一個(gè)實(shí)施例中，所述根據(jù)1譜，選擇不同的反演布點(diǎn)的橫向弛豫時(shí)間T2作為橫 向弛豫時(shí)間1~2的截止值T2rt，計(jì)算不同截止值T2et下的優(yōu)化譜系數(shù)，按照如下公式進(jìn)行計(jì) 算：
[0010]
[0011] 其中，Cbvi_w_為第j個(gè)反演布點(diǎn)的橫向弛豫時(shí)間對(duì)應(yīng)的優(yōu)化譜系數(shù)，無量綱；
[0012] T2j為第j個(gè)反演布點(diǎn)的橫向弛豫時(shí)間，ms;
[0013] T2ct為橫向弛豫時(shí)間T2的截止值，ms;
[0014] a，b是擬合系數(shù)。
[0015] 在一個(gè)實(shí)施例中，所述擬合系數(shù)a，b按照如下方式獲得：
[0016] 反演多塊巖心的核磁弛豫回波數(shù)據(jù)，得到多塊巖心的核磁共振橫向弛豫時(shí)間T2 譜；
[0017] 根據(jù)1譜，計(jì)算多塊巖心的橫向弛豫時(shí)間丁2的幾何平均值；
[0018] 獲得多塊巖心的束縛水飽和度數(shù)據(jù)；
[0019] 將所得的多塊巖心的橫向弛豫時(shí)間T2的幾何平均值和多塊巖心的束縛水飽和度 數(shù)據(jù)進(jìn)行擬合，獲得擬合系數(shù)a，b。
[0020] 在一個(gè)實(shí)施例中，所述將所得的多塊巖心的橫向弛豫時(shí)間t2的幾何平均值和多塊 巖心的束縛水飽和度數(shù)據(jù)進(jìn)行擬合，按照如下公式進(jìn)行計(jì)算：
[0021]
[0022] 其中，Sbvi是巖心的束縛水飽和度；
[0023] T21n^巖心的橫向弛豫時(shí)間T2的幾何平均值，ms。
[0024] 在一個(gè)實(shí)施例中，所述根據(jù)不同截止值T2et下的優(yōu)化譜系數(shù)，計(jì)算不同截止值T2c;t 下的巖心的核磁束縛水飽和度，按照如下公式計(jì)算；
[0025]
[0026] 其中，NSbvi為核磁束縛水飽和度；
[0027] A」為第j個(gè)反演布點(diǎn)對(duì)應(yīng)的孔隙度的幅度值。
[0028] 本發(fā)明實(shí)施例還提供了一種確定束縛水飽和度的裝置，通過優(yōu)化方法獲得最優(yōu)譜 系數(shù)，從而得到最優(yōu)核磁束縛水飽和度，提高了核磁束縛水飽和度的計(jì)算精度。該裝置包 括：
[0029] T2譜獲取單元，用于反演巖心的核磁弛豫回波數(shù)據(jù)，得到核磁共振橫向弛豫時(shí)間 丁 2譜；
[0030] 束縛水飽和度獲取單元，用于獲得巖心的束縛水飽和度數(shù)據(jù)；
[0031] 優(yōu)化譜系數(shù)獲取單元，用于根據(jù)1~2譜，選擇不同的反演布點(diǎn)的橫向弛豫時(shí)間1~ 2作 為橫向弛豫時(shí)間1~2的截止值T2rt，計(jì)算不同截止值T2rt下的優(yōu)化譜系數(shù)；
[0032] 核磁束縛水飽和度獲取單元，用于根據(jù)不同截止值T2rt下的優(yōu)化譜系數(shù)，計(jì)算不同 截止值T2rt下的巖心的核磁束縛水飽和度；
[0033] 最優(yōu)核磁束縛水飽和度確定單元，用于分別計(jì)算巖心的束縛水飽和度與不同截止 值T2c;t下的巖心的核磁束縛水飽和度的誤差，將計(jì)算出的最小誤差相對(duì)應(yīng)的優(yōu)化譜系數(shù)確 定為最優(yōu)譜系數(shù)，將最優(yōu)譜系數(shù)相對(duì)應(yīng)的核磁束縛水飽和度確定為最優(yōu)核磁束縛水飽和 度。
[0034] 在一個(gè)實(shí)施例中，所述優(yōu)化譜系數(shù)獲取單元具體用于按如下公式計(jì)算不同截止值 T2c;t下的優(yōu)化譜系數(shù)：
[0035]
[0036] 其中，Cbvi_w_為第j個(gè)反演布點(diǎn)的橫向弛豫時(shí)間對(duì)應(yīng)的優(yōu)化譜系數(shù)，無量綱；
[0037] T2j為第j個(gè)反演布點(diǎn)的橫向弛豫時(shí)間，ms;
[0038] T2c;t為橫向弛豫時(shí)間T2的截止值，ms;
[0039] a，b是擬合系數(shù)。
[0040] 在一個(gè)實(shí)施例中，所述優(yōu)化譜系數(shù)獲取單元具體用于按如下方式獲得所述擬合系 數(shù)a，b:
[0041] 反演多塊巖心的核磁弛豫回波數(shù)據(jù)，得到多塊巖心的核磁共振橫向弛豫時(shí)間t2 譜；
[0042] 根據(jù)1~2譜，計(jì)算多塊巖心的橫向弛豫時(shí)間T2的幾何平均值；
[0043] 獲得多塊巖心的束縛水飽和度數(shù)據(jù)；
[0044] 將所得的多塊巖心的橫向弛豫時(shí)間T2的幾何平均值和多塊巖心的束縛水飽和度 數(shù)據(jù)進(jìn)行擬合，獲得擬合系數(shù)a，b。
[0045] 在一個(gè)實(shí)施例中，所述優(yōu)化譜系數(shù)獲取單元具體用于按如下公式將所得的多塊巖 心的橫向弛豫時(shí)間T2的幾何平均值和多塊巖心的束縛水飽和度數(shù)據(jù)進(jìn)行擬合：
[0046]
[0047] 其中，Sbvi是巖心的束縛水飽和度；
[0048] 丁^是巖心的橫向弛豫時(shí)間T2的幾何平均值，ms。
[0049] 在一個(gè)實(shí)施例中，所述核磁束縛水飽和度獲取單元具體用于按如下公式計(jì)算不同 截止值T2rt下的巖心的核磁束縛水飽和度：
[0050]
[0051] 其中，NSbvi為核磁束縛水飽和度；
[0052] A#第j個(gè)反演布點(diǎn)對(duì)應(yīng)的孔隙度的幅度值。
[0053] 在本發(fā)明實(shí)施例中，根據(jù)1~2譜，選擇不同的反演布點(diǎn)的橫向弛豫時(shí)間1~ 2作為截止 值T2rt，對(duì)譜系數(shù)進(jìn)行優(yōu)化，得到不同T2rt下的優(yōu)化譜系數(shù)，根據(jù)優(yōu)化的譜系數(shù)計(jì)算得到核磁 束縛水飽和度，計(jì)算巖心的束縛水飽和度數(shù)據(jù)和核磁束縛水飽和度的誤差，將計(jì)算出的最 小誤差相對(duì)應(yīng)的優(yōu)化譜系數(shù)確定為最優(yōu)譜系數(shù)，將最優(yōu)譜系數(shù)相對(duì)應(yīng)的核磁束縛水飽和度 確定為最優(yōu)核磁束縛水飽和度，與現(xiàn)有技術(shù)中通過簡(jiǎn)單的數(shù)據(jù)擬合獲得譜系數(shù)，再根據(jù)譜 系數(shù)計(jì)算巖石的束縛水飽和度的方法相比較，采用本發(fā)明方法可以獲得具有較好精度的譜 系數(shù)，從而提高束縛水飽和度的計(jì)算精度。
【附圖說明】
[0054] 此處所說明的附圖用來提供對(duì)本發(fā)明的進(jìn)一步理解，構(gòu)成本申請(qǐng)的一部分，并不 構(gòu)成對(duì)本發(fā)明的限定。在附圖中：
[0055] 圖1是本發(fā)明實(shí)施例提供的一種確定束縛水飽和度的方法流程圖；
[0056] 圖2是本發(fā)明實(shí)施例提供的一種核磁共振橫向弛豫時(shí)間T2譜；
[0057] 圖3是本發(fā)明實(shí)施例提供的一種選擇T2c;t改造原始譜系數(shù)的方法示意圖；
[0058] 圖4是本發(fā)明實(shí)施例提供的一種確定束縛水飽和度的裝置結(jié)構(gòu)圖；
[0059] 圖5是本發(fā)明實(shí)施例提供的一種使用原始譜系數(shù)計(jì)算的束縛水飽和度和實(shí)驗(yàn)測(cè) 得的束縛水飽和度之間的相關(guān)性示意圖；
[0060] 圖6是本發(fā)明實(shí)施例提供的一種使用改進(jìn)譜系數(shù)計(jì)算的束縛水飽和度和實(shí)驗(yàn)測(cè) 得的束縛水的比較結(jié)果圖；
[0061] 圖7是本發(fā)明實(shí)施例提供的一種改進(jìn)譜系數(shù)計(jì)算束縛水飽和度的盲點(diǎn)檢查交會(huì) 圖。
【具體實(shí)施方式】
[0062] 為使本發(fā)明的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚明白，下面結(jié)合實(shí)施方式和附圖，對(duì) 本發(fā)明做進(jìn)一步詳細(xì)說明。在此，本發(fā)明的示意性實(shí)施方式及其說明用于解釋本發(fā)明，但并 不作為對(duì)本發(fā)明的限定。
[0063] 在根據(jù)核磁實(shí)驗(yàn)數(shù)據(jù)計(jì)算束縛水飽和度的算法中，關(guān)鍵是要得到合理的譜系數(shù)， 而現(xiàn)有的只是簡(jiǎn)單利用數(shù)據(jù)擬合的方法獲得譜系數(shù)，這種方法不能獲得較好精度的譜系 數(shù)，同樣根據(jù)譜系數(shù)計(jì)算所得的束縛水飽和度的精度也是不高的。如果可以提出一種方法， 通過對(duì)譜系數(shù)進(jìn)行優(yōu)化獲得最