一種基于油藏模型的測井?dāng)?shù)據(jù)歷史回歸處理方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于油田開發(fā)的測井、地震、油藏綜合數(shù)據(jù)處理技術(shù)領(lǐng)域，特別是涉及到一 種基于油藏模型的測井?dāng)?shù)據(jù)歷史回歸處理方法。
【背景技術(shù)】
[0002] 石油測井是石油勘探、開發(fā)的"眼睛"。在油氣田勘探、開發(fā)的不同階段，石油測井 的目的和任務(wù)是不同的。一般來說，裸眼井測井（下套管之前的井稱裸眼井，因此下套管之 前進(jìn)行的測井稱為裸眼井測井）的主要目的和任務(wù)是發(fā)現(xiàn)和評價(jià)油氣層的儲集性能及生 產(chǎn)能力；而生產(chǎn)測井（油水井投入生產(chǎn)以后進(jìn)行的測井稱為生產(chǎn)測井）的主要目的是見識 和分析油氣層的開發(fā)動(dòng)態(tài)及生產(chǎn)情況。開發(fā)測井是指在油氣田整個(gè)開發(fā)期間進(jìn)行的所有測 井項(xiàng)目。開發(fā)測井的主要對象為裸眼完成的生產(chǎn)井和下套管的生產(chǎn)井，用于分析目前的生 產(chǎn)動(dòng)態(tài)及井內(nèi)技術(shù)狀況。
[0003] 測井?dāng)?shù)據(jù)詳細(xì)的記錄了地層巖性、物性在縱向上的細(xì)微變化，具有較高的縱向分 辨率，在井位置處，能準(zhǔn)確的計(jì)算儲層參數(shù)，正確評價(jià)儲層。地震資料連續(xù)地記錄了地層在 橫向上的變化，具有較好的橫向連續(xù)性。因此，在勘探和開發(fā)階段往往是結(jié)合測井和地震數(shù) 據(jù)進(jìn)行綜合研究。通常測井?dāng)?shù)據(jù)的采集時(shí)間點(diǎn)與地震數(shù)據(jù)的時(shí)間點(diǎn)很少一致，因此，當(dāng)?shù)?下儲層物性變化較大時(shí)，如果直接拿來進(jìn)行綜合研究，則會(huì)因采集時(shí)間點(diǎn)的不一致給綜合 研究帶來很大誤差，因此，需要對測井?dāng)?shù)據(jù)進(jìn)行歷史回歸處理，將實(shí)際測井?dāng)?shù)據(jù)校正到地震 數(shù)據(jù)的采集時(shí)間點(diǎn)，為更好、更準(zhǔn)確的研究儲層建立基礎(chǔ)。

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

[0004] 本發(fā)明的目的是針對現(xiàn)有測井技術(shù)在開發(fā)應(yīng)用中存在的問題，提出基于油藏模型 的測井?dāng)?shù)據(jù)歷史回歸處理方法，即測井與地震進(jìn)行綜合研究時(shí)，二者的采集時(shí)間點(diǎn)往往不 一致所帶來的系統(tǒng)誤差。
[0005] 本發(fā)明總體技術(shù)方案是
[0006] -種基于油藏模型的測井?dāng)?shù)據(jù)歷史回歸方法，包括：
[0007] 1)根據(jù)地震解釋數(shù)據(jù)、測井解釋數(shù)據(jù)、生產(chǎn)數(shù)據(jù)以及巖心數(shù)據(jù)，通過歷史擬合，建 立油藏模型；
[0008] 2)利用工區(qū)測井解釋數(shù)據(jù)和油藏模型數(shù)據(jù)對巖石物理模型標(biāo)定，建立地震響應(yīng)參 數(shù)與巖石參數(shù)之間的非線性關(guān)系；
[0009] 3)合成油藏模型的動(dòng)態(tài)縱波速度、橫波速度及密度；
[0010] 4)根據(jù)合成的動(dòng)態(tài)縱波速度、橫波速度及密度，以地震數(shù)據(jù)采集時(shí)間點(diǎn)作為基準(zhǔn)， 分別計(jì)算縱波速度、橫波速度及密度的變化率；
[0011] 5)利用實(shí)測的縱波速度、橫波速度、密度以及步驟四中計(jì)算的縱波速度、橫波速 度、密度的變化率，計(jì)算對應(yīng)地震數(shù)據(jù)采集時(shí)間點(diǎn)的實(shí)測縱波速度、橫波速度、密度；
[0012] 6)利用回歸后的縱波速度、橫波速度、密度進(jìn)行正演，得到的合成記錄與采集的地 震記錄進(jìn)行一致性對比分析，以判別回歸的測井?dāng)?shù)據(jù)是否正確。
[0013] 上述方案的細(xì)化方案包括：
[0014] 1)根據(jù)工區(qū)地震解釋數(shù)據(jù)、測井解釋數(shù)據(jù)、生產(chǎn)數(shù)據(jù)以及巖心數(shù)據(jù)，通過歷史擬 合，建立油藏模型；
[0015] 所述地震解釋數(shù)據(jù)包括：由地震數(shù)據(jù)解釋的層位，斷層數(shù)據(jù)；
[0016] 所述測井解釋數(shù)據(jù)包括：地質(zhì)分層、單井相、砂泥巖分布、滲透率、含油飽和度、含 水飽和度；
[0017] 所述生產(chǎn)數(shù)據(jù)包括：實(shí)時(shí)的流體的流量、含水飽和度、溫度、壓力；
[0018] 所述巖心數(shù)據(jù)是指鉆井時(shí)從地下取出的巖心，在實(shí)驗(yàn)室進(jìn)行測量，主要包括：滲透 率、孔隙度、電阻率；
[0019] 所述的建立油藏模型：第一步是利用地震解釋的層位和斷層數(shù)據(jù)，通過空間上的 網(wǎng)格剖分，建立空間上的格架模型，并將油藏剖分為一系列獨(dú)立的單元體；第二步是結(jié)合測 井?dāng)?shù)據(jù)以及巖心數(shù)據(jù)，通過空間上的插值，在空間上的每一個(gè)單元中賦值，賦予不同的物理 屬性；第三步：通過模擬器對歷史的油氣藏動(dòng)態(tài)數(shù)據(jù)進(jìn)行模擬，并且與實(shí)際的井點(diǎn)上的油 氣藏動(dòng)態(tài)數(shù)據(jù)進(jìn)行對比，以確定模擬的油氣藏動(dòng)態(tài)數(shù)據(jù)是否準(zhǔn)確，如果模擬的動(dòng)態(tài)數(shù)據(jù)與 油藏實(shí)際動(dòng)態(tài)數(shù)據(jù)差別較大，則調(diào)整輸入的模型參數(shù)，重新進(jìn)行擬合，最終得到與實(shí)際油氣 藏動(dòng)態(tài)數(shù)據(jù)一致的油藏模型；
[0020] 2)利用工區(qū)測井解釋數(shù)據(jù)和油藏模型數(shù)據(jù)對巖石物理模型標(biāo)定，建立地震響應(yīng)參 數(shù)與巖石參數(shù)之間的非線性關(guān)系；
[0021] 所述的巖石物理模型是Gassmann理論模型；
[0022] 所述標(biāo)定是指確定巖石物理模型中的各項(xiàng)參數(shù)；
[0023] 所述的建立地震響應(yīng)參數(shù)與巖石參數(shù)之間的非線性關(guān)系是建立聲波參數(shù)、彈性參 數(shù)與巖石物性、流體參數(shù)之間的關(guān)系；
[0024] 3)根據(jù)油藏模型的流體參數(shù)及彈性參數(shù)，利用標(biāo)定后的巖石物理模型，計(jì)算油藏 模型的縱波速度、橫波速度以及密度；再將縱波速度與密度相乘得到合成的阻抗；
[0025] 所述計(jì)算采用gassmann模型，巖石體積彈性模量K計(jì)算：
[0026]
[0027] 其中，Ks為巖石顆粒彈性模量；Kd為干巖石彈性模量；〇為孔隙度；
[0028] Kf為流體彈性模量，由下式計(jì)算：
[0029]
;其中，Kw，K。，Kg分別是水、油、氣的彈性模量，S w， S。分別是水、油的飽和度，1-SW-S。為含氣飽和度；
[0030] 聲波在巖石中的縱波速度Vp由下式計(jì)算：
[0031]
其中，y為巖石剪切模量；P為巖石密度；Vs為橫 波速度。
[0032]
?，其中，Pw、P。、Pg*別是水、 油、氣密度，Pm為巖石骨架密度。
[0033] 4)利用步驟3)中計(jì)算的各個(gè)時(shí)間點(diǎn)的縱波速度、橫波速度、密度，與地震數(shù)據(jù)采 集時(shí)間點(diǎn)時(shí)的縱波速度、橫波速度、密度，計(jì)算油藏模型在各個(gè)開發(fā)時(shí)間點(diǎn)的變化率；
[0034] 所述油藏模型的各個(gè)開發(fā)時(shí)間點(diǎn)是在建立油藏模型時(shí)，人為給定的具有相同時(shí)間 間隔的各個(gè)時(shí)間狀態(tài)；
[0035] 所述縱波為傳播方向與振動(dòng)防線一致的波，橫波為振動(dòng)方向與傳播方向垂直的 波；
[0036] 5)對比實(shí)測測井?dāng)?shù)據(jù)時(shí)間點(diǎn)的縱波速度、橫波速度，密度與油藏模型對應(yīng)時(shí)間點(diǎn) 的縱波速度、橫波速度，密度，判斷二者是否一致，如不一致，則返回步驟2)，調(diào)整巖石物理 模型參數(shù)，如一致，則可進(jìn)入步驟6);
[0037] 所述實(shí)測的測井?dāng)?shù)據(jù)是指利用測井儀器在井筒中實(shí)際測量的曲線數(shù)據(jù)，包括聲波 時(shí)差（DT)、橫波聲波時(shí)差（DTS)、密度（RH0B)，其中，縱波速度為DT的倒數(shù)乘以106,橫波速 度為DTS的倒數(shù)乘以10 6 ;
[0038] 所述二者是否一致是指同一時(shí)間狀態(tài)標(biāo)準(zhǔn)下合成縱波速度、橫波速度、密度與實(shí) 際測量的縱波速度、橫波速度、密度在趨勢以及數(shù)值范圍上是否一致；
[0039] 6)結(jié)合實(shí)測測井?dāng)?shù)據(jù)與油藏模型計(jì)算的各個(gè)開發(fā)時(shí)間段的變化率，將實(shí)測測井?dāng)?shù) 據(jù)回歸到地震數(shù)據(jù)采集時(shí)間點(diǎn)；
[0040] 所述變化率是指以地震數(shù)據(jù)采集的時(shí)間點(diǎn)作為基礎(chǔ)，計(jì)算別的時(shí)間點(diǎn)的油藏動(dòng)態(tài) 變化百分比；
[0041] 所述回歸是指以計(jì)算的變化率，對實(shí)測測井?dāng)?shù)據(jù)進(jìn)行基于時(shí)間點(diǎn)的校正。
[0042] 其中：步驟2)中標(biāo)定巖石物理模型中的各項(xiàng)參數(shù)包括巖石巖性、孔隙度、孔隙類 型、孔隙流體、流體飽和度和頻率，以及反應(yīng)巖石中彈性性質(zhì)的包括縱波速度和橫波速度及 衰減的參數(shù)。
[0043] 本發(fā)明的技術(shù)效果是：由于油藏模型數(shù)據(jù)是一個(gè)隨時(shí)間變化的動(dòng)態(tài)數(shù)據(jù)，其客觀、