此，本申請(qǐng)實(shí)施例的示意性實(shí)施例及其說(shuō)明用于 解釋本申請(qǐng)實(shí)施例，但并不作為對(duì)本申請(qǐng)實(shí)施例的限定。
[0044]下面結(jié)合附圖，對(duì)本申請(qǐng)實(shí)施例的【具體實(shí)施方式】作進(jìn)一步的詳細(xì)說(shuō)明。
[0045] 參考圖1所示，本申請(qǐng)一實(shí)施例基于流度的可采稠油油藏的監(jiān)測(cè)方法，包括以下步 驟：
[0046] 步驟Sl，確定儲(chǔ)層中飽和巖石的巖石密度、滲透率和流體粘度。
[0047] 儲(chǔ)層中飽和巖石的滲透率、流體粘度可通過(guò)地質(zhì)報(bào)告、錄井報(bào)告、巖心數(shù)據(jù)和測(cè)井 曲線(xiàn)等得到。而儲(chǔ)層中飽和巖石的巖石密度則可通過(guò)以下方式得到：
[0048] 首先，根據(jù)公式
獲取所述儲(chǔ)層中飽和巖石的基質(zhì)密度口8，其中P 1表示 第i種組成礦物的密度，fi表示第i種礦物的體積分?jǐn)?shù)；
[0049] 然后，根據(jù)公式Pb=(I-Ci) )pg+(i)pf獲取所述儲(chǔ)層中飽和巖石的巖石密度，其中，Φ 為所述飽和巖石的孔隙度，Pf所述飽和巖石內(nèi)的流體密度。
[0050] 步驟S2,在開(kāi)采所述儲(chǔ)層的過(guò)程中，獲取所述儲(chǔ)層在至少兩次地震勘探下采集的 地震數(shù)據(jù)，并對(duì)各次地震勘探下采集的地震數(shù)據(jù)進(jìn)行時(shí)移互均化處理，各次地震勘探的時(shí) 間間隔為預(yù)設(shè)值。本文的時(shí)移互均化處理指的是對(duì)疊前時(shí)移地震數(shù)據(jù)或疊后時(shí)移地震數(shù)據(jù) 進(jìn)行互均化處理。
[0051] -般地，在通過(guò)蒸汽輔助重力泄油等方式開(kāi)采儲(chǔ)層的過(guò)程中，往往需要每隔一段 時(shí)間對(duì)儲(chǔ)層進(jìn)行一次地震勘探并采集地震數(shù)據(jù)。根據(jù)實(shí)際需要，各次地震勘探的時(shí)間間隔 可以預(yù)先設(shè)定，比如設(shè)定為等間隔，例如每年一次、每月一次、每周一次、甚至是每天一次。 當(dāng)然，根據(jù)需要也可是不等間隔。
[0052]步驟S3,對(duì)各次經(jīng)時(shí)移互均化處理后的地震數(shù)據(jù)進(jìn)行頻譜分析，并從各次的頻譜 分析結(jié)果中選取相同頻帶的頻譜能量，所述相同頻帶為低頻頻帶，因?yàn)榈皖l頻帶最為相關(guān)。 其中，可通過(guò)連續(xù)小波變換(例如Morlet小波變換)對(duì)各次經(jīng)時(shí)移互均化處理后的地震數(shù)據(jù) 進(jìn)行頻譜分析。合成信號(hào)如圖2a所示，合成信號(hào)經(jīng)過(guò)Morlet小波變換得到的頻譜如圖2b所 示，如圖3所示，為本申請(qǐng)實(shí)施例的一個(gè)示例性的頻帶選取示意圖。
[0053] Si Iin等基于Biot彈性波傳播理論和Barenblatt方法得到了Biot-Barenblatt孔 彈性模型，在引入一個(gè)無(wú)量綱參數(shù)對(duì)彈性介質(zhì)與飽含流體多孔介質(zhì)分界面的地震反射進(jìn)行 了低頻漸近分析，當(dāng)流體和巖石的特性在合理范圍時(shí)，無(wú)量綱參數(shù)在低頻(一般低于IKHz) 時(shí)很小，則在反射界面上，可得到飽含流體多孔介質(zhì)中低頻地震反射系數(shù)R的漸近表示。具 體的：
[0059]其中，V1是聲速，即介質(zhì)的一個(gè)特征參數(shù)，該參數(shù)不依賴(lài)于頻率；
為飽和巖石內(nèi)流體相對(duì)于骨架的平均速度，Pb為飽和巖石的巖石密 度，Pf為飽和巖石內(nèi)流體密度，PgS飽和巖石的基質(zhì)密度，其關(guān)系如下：

4為干骨架的壓縮系數(shù)。
[0073] βι，β2，β?分別為介質(zhì)1，介質(zhì)2和流體的體積模量的倒數(shù)，顯然Ro，Ri是與頻率無(wú)關(guān)的 量，反射系數(shù)對(duì)頻率求導(dǎo)公式為：
[0075]進(jìn)而得到流度m:
[0077]其中在實(shí)際應(yīng)用中，一般將C常數(shù)化，ω為優(yōu)勢(shì)頻率。
[0078]步驟S4，確定各次選取的頻譜能量隨頻率變化的導(dǎo)數(shù)。
[0079] 在本申請(qǐng)的一實(shí)施例中，具體可根據(jù)公式（
確定各次選取的頻 譜能量隨頻率變化的導(dǎo)數(shù)，其中，
為第i次的頻譜能量隨頻率變化的導(dǎo)數(shù)，Ph1為第i次 的低頻頻帶的高頻段對(duì)應(yīng)的頻譜能量，PU為第i次的低頻頻帶的低頻段對(duì)應(yīng)的頻譜能量，η 謂低頻頻帶的高、低頻段所取的樣點(diǎn)個(gè)數(shù)，R為飽和巖石的，w為頻率。
[0080] 步驟S5,基于所述飽和巖石的滲透率、流體粘度以及所述各次選取的頻譜能量隨 頻率變化的導(dǎo)數(shù)對(duì)應(yīng)確定各次的流度比例系數(shù)。
[0081] 在本申請(qǐng)的一實(shí)施例中，具體可根據(jù)公式 計(jì)算各次的流度比例系 數(shù)，其中，C1為第i次的流度比例系數(shù)，K為所述飽和巖石的滲透率，η所述飽和巖石的流體粘 度，
為第i次的頻譜能量隨頻率變化的導(dǎo)數(shù)。
[0082] 步驟S6,基于所述各次的流度比例系數(shù)及所述各次選取的頻譜能量隨頻率變化的 導(dǎo)數(shù)對(duì)應(yīng)確定各次的流度。
[0083] 在本申請(qǐng)的一實(shí)施例中，具體可根據(jù)公式
計(jì)算各次的流度，其 中，nu第i次的流度，Ci為第i次的流度比例系1
為第i次的頻譜能量隨頻率變化的導(dǎo) 數(shù)，w為頻率。
[0084] 步驟S7,根據(jù)各次的流度確定所述儲(chǔ)層的流度變化，并基于所述儲(chǔ)層的流度變化 確定所述儲(chǔ)層的可采稠油油藏的變化剖面。
[0085] 其中，各次的流度確定所述儲(chǔ)層的流度變化可通過(guò)計(jì)算所述各次的流度中相鄰次 流度間的差值，獲得流度差值集合，所述流度差值集合反應(yīng)所述儲(chǔ)層的流度變化。例如各次 的流度分別為，a、b、c、d、e、f、g，那么各次的流度中相鄰次流度間的差值則為:b-a、c-b、(Ι- ο'θ-?Ι,?'-θΑ-?·。
[0086] 稠油油藏在受熱狀態(tài)下其流動(dòng)性會(huì)相應(yīng)的增強(qiáng)，而利用蒸汽輔助重力泄油開(kāi)采稠 油的過(guò)程中，稠油油藏中流動(dòng)性強(qiáng)的部分會(huì)隨著隨著蒸汽腔的擴(kuò)展在原始稠油油藏邊界范 圍內(nèi)擴(kuò)展。例如某一儲(chǔ)層在09年和11年分別進(jìn)行了一次地震勘探，獲得的地震數(shù)據(jù)在經(jīng)過(guò) 本申請(qǐng)實(shí)施例的方法處理后，得到如圖4所示的可采稠油油藏的變化剖面。在該圖中可以看 出，隨著時(shí)間推移，蒸汽腔在不斷的擴(kuò)展，蒸汽腔的邊界變化則反映了可采稠油油藏的邊界 變化。
[0087] 由此可見(jiàn)，本申請(qǐng)實(shí)施例充分利用稠油油藏粘度變化大，引起流度變化大，利用多 次采集的地震資料預(yù)測(cè)儲(chǔ)層的流度變化，從而根據(jù)流度變化獲得樂(lè)可采稠油油藏的邊界變 化。
[0088] 參考圖5所示，本申請(qǐng)一實(shí)施例基于流度的可采稠油油藏的監(jiān)測(cè)裝置，包括：
[0089] 參數(shù)獲取模塊51，用于確定儲(chǔ)層中飽和巖石的巖石密度、滲透率和流體粘度；
[0090] 地震數(shù)據(jù)獲取模塊52,用于在開(kāi)采所述儲(chǔ)層的過(guò)程中，獲取所述儲(chǔ)層在至少兩次 地震勘探下采集的地震數(shù)據(jù)，并對(duì)各次地震勘探下采集的地震數(shù)據(jù)進(jìn)行時(shí)移互均化處理， 各次地震勘探的時(shí)間間隔為預(yù)設(shè)值；
[0091] 頻譜處理模塊53,用于對(duì)各次經(jīng)時(shí)移互均化處理后的地震數(shù)據(jù)進(jìn)行頻譜分析，并 從各次的頻譜分析結(jié)果中選取相同頻帶的頻譜能量，所述相同頻帶為低頻頻帶；
[0092] 導(dǎo)數(shù)獲取模塊54,用于確定各次選取的頻譜能量隨頻率變化的導(dǎo)數(shù)；
[0093] 比例系數(shù)獲取模塊55,用于基于所述飽和巖石的滲透率、流體粘度以及所述各次 選取的頻譜能量隨頻率變化的導(dǎo)數(shù)對(duì)應(yīng)確定各次的流度比例系數(shù)；
[0094] 流度獲取模塊56,用于基于所述各次的流度比例系數(shù)及所述各次選取的頻譜能量 隨頻率變化的導(dǎo)數(shù)對(duì)應(yīng)確定各次的流度；
[0095] 油藏剖面獲取模塊57,用于根據(jù)各次的流度確定所述儲(chǔ)層的流度變化，并基于所 述儲(chǔ)層的流度變化確定所述儲(chǔ)層的可采稠油油藏的變化剖面。
[0096] 由于該基于流度的可采稠油油藏的監(jiān)測(cè)裝置與上述實(shí)施例中，基于流度的可采稠 油油藏的監(jiān)測(cè)方法對(duì)應(yīng)，其各模塊的具體情況請(qǐng)參見(jiàn)上述方法實(shí)施例中對(duì)應(yīng)的步驟，在此 不再贅述。
[0097] 由此可見(jiàn)，本申請(qǐng)實(shí)施例充分利用稠油油藏粘度變化大，引起流度變化大，利用多 次采集的地震資料預(yù)測(cè)儲(chǔ)層的流度變化，從而根據(jù)流度變化獲得樂(lè)可采稠油油藏的邊界變 化。
[0098] 本領(lǐng)域技術(shù)人員還可以了解到本申請(qǐng)實(shí)施例列出的各種說(shuō)明性邏輯塊、單元和步 驟可以通過(guò)硬件、軟件或兩者的結(jié)合來(lái)實(shí)現(xiàn)。至于是通過(guò)硬件還是軟件來(lái)實(shí)現(xiàn)取決于特定 的應(yīng)用和整個(gè)系統(tǒng)的設(shè)計(jì)要求。本領(lǐng)域技術(shù)人員可以對(duì)于每種特定的應(yīng)用，可以使用各種 方法實(shí)現(xiàn)所述的功能，但這種實(shí)現(xiàn)不應(yīng)被理解為超出本申請(qǐng)實(shí)施例保