真實(shí)的反應(yīng)了地下油藏在開發(fā)過程中的變化，因此，以巖石物理模型為橋梁將地下油藏的 變化轉(zhuǎn)換為聲學(xué)上的變化，即得到地下油藏在開發(fā)中縱波速度、橫波速度、密度的變化，將 這種真實(shí)的變化再應(yīng)用到實(shí)測的縱波速度、橫波速度、密度中，則可以將實(shí)測的測井?dāng)?shù)據(jù)回 歸到地震數(shù)據(jù)采集的時間點(diǎn)上，這樣在進(jìn)行測井、地震綜合研究時，二者處于同一測量時間 點(diǎn)，為準(zhǔn)確預(yù)測地下油藏位置及油藏變化建立良好的基礎(chǔ)。
【附圖說明】
[0044] 圖1為基于油藏模型的測井?dāng)?shù)據(jù)歷史回歸處理方法流程；
[0045] 圖2為根據(jù)油藏模型和井點(diǎn)解釋結(jié)果，調(diào)整為最佳的巖石物理模型參數(shù)；
[0046] 圖3為巖石物理模型模板，縱波速度與孔隙度之間的對應(yīng)關(guān)系；
[0047] 圖4為在井點(diǎn)位置通過油藏模型合成的縱波速度、橫波速度；
[0048] 圖5為在井點(diǎn)位置分別計(jì)算的在兩個地震數(shù)據(jù)采集時間點(diǎn)，縱波速度、橫波速度 的變化率；
[0049] 圖6為在井點(diǎn)位置分別回歸到與兩個地震數(shù)據(jù)采集時間點(diǎn)的實(shí)測的縱波速度、橫 波速度；
[0050] 圖7為在井點(diǎn)位置在第一個地震數(shù)據(jù)采集時間點(diǎn)的測井正演結(jié)果與實(shí)際采集地 震數(shù)據(jù)對比；
[0051] 圖8為在井點(diǎn)位置在第二個地震數(shù)據(jù)采集時間點(diǎn)的測井正演結(jié)果與實(shí)際采集地 震數(shù)據(jù)對比。
【具體實(shí)施方式】
[0052] 下面將結(jié)合附圖，對本發(fā)明的技術(shù)方案作進(jìn)一步說明。
[0053] 如圖1所示，測井?dāng)?shù)據(jù)歷史回歸處理方法流程。
[0054] 首先，根據(jù)地震解釋數(shù)據(jù)、測井解釋數(shù)據(jù)、生產(chǎn)數(shù)據(jù)以及巖心數(shù)據(jù)，通過歷史擬合， 建立油藏模型。地質(zhì)學(xué)家對油田進(jìn)行詳細(xì)的油藏描述工作，基于地震解釋數(shù)據(jù)，測井解釋數(shù) 據(jù)，巖芯數(shù)據(jù)以及結(jié)合地質(zhì)家對油田的認(rèn)識建立三維地質(zhì)模型。流程進(jìn)入到步驟2。
[0055] 其二，利用工區(qū)測井解釋數(shù)據(jù)和油藏模型數(shù)據(jù)對巖石物理模型標(biāo)定，建立地震響 應(yīng)參數(shù)與巖石參數(shù)之間的非線性關(guān)系。油藏模型數(shù)據(jù)里沒有速度和密度這些值，因此需要 利用測井?dāng)?shù)據(jù)解釋結(jié)果與油藏模型對應(yīng)的測井位置的一些動態(tài)參數(shù)進(jìn)行擬合，這樣一個過 程稱之為標(biāo)定。即：巖石物理模型標(biāo)定是將通過對各種巖心資料、測井資料和地震資料進(jìn) 行綜合分析，研究巖石巖性、孔隙度、孔隙類型、孔隙流體、流體飽和度和頻率等對巖石中彈 性性質(zhì)如縱波速度和橫波速度及衰減的影響，并提出利用地震響應(yīng)預(yù)測巖石物理性質(zhì)的理 論和方法，是地震響應(yīng)與儲層巖石參數(shù)之間聯(lián)系的橋梁，是進(jìn)行定量儲層預(yù)測的基本前提。 圖2為標(biāo)定后巖石物理模型參數(shù)，圖3為巖石物理模型模板，縱波速度與孔隙度之間的關(guān) 系。
[0056] 其三，合成油藏模型的動態(tài)縱波速度、橫波速度及密度，即根據(jù)油藏模型的各個時 間點(diǎn)，合成各個時間的縱波速度、橫波速度、密度。圖4為在井點(diǎn)位置在某一時間點(diǎn)合成的 縱波速度及橫波速度。
[0057] 其四，根據(jù)合成的動態(tài)縱波速度、橫波速度及密度，以地震數(shù)據(jù)采集時間點(diǎn)作為基 準(zhǔn)，分別計(jì)算縱波速度、橫波速度及密度的變化率，圖5為測井?dāng)?shù)據(jù)采集時間點(diǎn)相對地震數(shù) 據(jù)采集時間點(diǎn)的縱波速度、橫波速度、密度的變化率，可以看到在本實(shí)例中有兩次采集的地 震數(shù)據(jù)分別為1991年和2010年。
[0058] 其五，利用實(shí)測的縱波速度、橫波速度、密度以及步驟四中計(jì)算的縱波速度、橫波 速度、密度的變化率，計(jì)算對應(yīng)地震數(shù)據(jù)采集時間點(diǎn)的實(shí)測縱波速度、橫波速度、密度。圖6 為本實(shí)例中選取的井計(jì)算的經(jīng)過回歸后的縱波速度、橫波速度、密度。
[0059] 其六，利用回歸后的縱波速度、橫波速度、密度進(jìn)行正演，得到的合成記錄與采集 的地震記錄進(jìn)行一致性對比分析，以判別回歸的測井?dāng)?shù)據(jù)是否正確，圖7和圖8分別為向兩 個時間點(diǎn)（1991年和2010年）回歸的測井?dāng)?shù)據(jù)正演的結(jié)果與采集的地震記錄的對比分析。
[0060] 下面通過典型實(shí)施例對本發(fā)明作更加細(xì)化的說明。
[0061] 典型本實(shí)施例的包括下列步驟：
[0062] 1)根據(jù)勝利油田單56工區(qū)地震解釋數(shù)據(jù)、測井解釋數(shù)據(jù)、生產(chǎn)數(shù)據(jù)以及巖心數(shù) 據(jù)，通過歷史擬合，建立油藏模型；
[0063] 所述地震解釋數(shù)據(jù)包括：由地震數(shù)據(jù)解釋的層位，斷層數(shù)據(jù)；
[0064] 所述測井解釋數(shù)據(jù)包括：地質(zhì)分層、單井相、砂泥巖分布、滲透率、含油飽和度、含 水飽和度等；
[0065] 所述生產(chǎn)數(shù)據(jù)包括：實(shí)時的流體的流量、含水飽和度、溫度、壓力等；
[0066] 所述巖心數(shù)據(jù)是指鉆井時從地下取出的巖心，在實(shí)驗(yàn)室進(jìn)行測量，主要包括：滲透 率、孔隙度、電阻率等；
[0067] 所述歷史擬合是指由油藏模型被建立起來以后，是否完全反映油氣藏實(shí)際，并未 經(jīng)過檢驗(yàn)。只有利用將生產(chǎn)和注入的歷史數(shù)據(jù)輸入模型并運(yùn)行模擬器，再將計(jì)算結(jié)果與油 氣藏實(shí)際動態(tài)相比，才能確定模型中采用的油氣藏描述是否有效。若計(jì)算獲得的動態(tài)數(shù)據(jù) 與油藏實(shí)際動態(tài)數(shù)據(jù)差別甚遠(yuǎn)，我們就必須不斷調(diào)整輸入模型基本數(shù)據(jù)，直到由模擬器計(jì) 算得到的動態(tài)與油藏生產(chǎn)的實(shí)際動態(tài)達(dá)到滿意的擬合為止。
[0068] 所述的建立油藏模型：第一步是利用地震解釋的層位和斷層數(shù)據(jù)，通過空間上的 網(wǎng)格剖分，建立空間上的格架模型，并將油藏剖分為一系列獨(dú)立的單元體；第二步是結(jié)合測 井?dāng)?shù)據(jù)以及巖心數(shù)據(jù)，通過空間上的插值，在空間上的每一個單元中賦值，賦予不同的物理 屬性；第三步：通過模擬器對歷史的油氣藏動態(tài)數(shù)據(jù)進(jìn)行模擬，并且與實(shí)際的井點(diǎn)上的油 氣藏動態(tài)數(shù)據(jù)進(jìn)行對比，以確定模擬的油氣藏動態(tài)數(shù)據(jù)是否準(zhǔn)確。如果模擬的動態(tài)數(shù)據(jù)與 油藏實(shí)際動態(tài)數(shù)據(jù)差別較大，則調(diào)整輸入的模型參數(shù)，重新進(jìn)行擬合，最終得到與實(shí)際油氣 藏動態(tài)數(shù)據(jù)一致的油藏模型。
[0069] 2)利用工區(qū)測井解釋數(shù)據(jù)和油藏模型數(shù)據(jù)對巖石物理模型標(biāo)定，建立地震響應(yīng)參 數(shù)與巖石參數(shù)之間的非線性關(guān)系；
[0070] 所述的巖石物理模型是Gassmann理論模型。
[0071] 所述標(biāo)定是指確定巖石物理模型中的各項(xiàng)參數(shù)。
[0072] 所述的建立地震響應(yīng)參數(shù)與巖石參數(shù)之間的非線性關(guān)系是建立聲波參數(shù)、彈性參 數(shù)與巖石物性、流體參數(shù)之間的關(guān)系。
[0073] 3)根據(jù)油藏模型的流體參數(shù)及彈性參數(shù)，利用標(biāo)定后的巖石物理模型，計(jì)算油藏 模型的縱波速度、橫波速度以及密度；再將縱波速度與密度相乘得到合成的阻抗；
[0074] 所述的巖石物理模型是Gassmann理論模型。
[0075] 所述計(jì)算采用gassmann模型，巖石體積彈性模量K計(jì)算：
[0076]
[0077] 其中，Ks為巖石顆粒彈性模量；Kd為干巖石彈性模量；〇為孔隙度；
[0078] Kf為流體彈性模量，由下式計(jì)算：
[0079]
;其中，Kw，K。，Kg分別是水、油、氣的彈性模量，S w， S。分別是水、油的飽和度，1-SW-S。為含氣飽和度；
[0080] 聲波在巖石中的縱波速度Vp由下式計(jì)算：
[0081]
〔中，y為巖石剪切模量；P為巖石密度；Vs為橫 波速度^
[0082]
其中，P w、P。、P g分別是水、 油、氣密度，Pm為巖石骨架密度。
[0083] 4)利用步驟3)中計(jì)算的各個時間點(diǎn)的縱波速度、橫波速度、密度，與地震數(shù)據(jù)采 集時間點(diǎn)時的縱波速度、橫波速度、密度，計(jì)算油藏模型在各個開發(fā)時間點(diǎn)的變化率；
[0084] 所述油藏模型的各個開發(fā)時間點(diǎn)是在建立油藏模型時