一種基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法
【專利摘要】本發(fā)明公開了一種基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法，所述方法包括以下步驟：確定預(yù)設(shè)數(shù)量的速度控制層；對(duì)疊后時(shí)間偏移地震數(shù)據(jù)進(jìn)行井震標(biāo)定并對(duì)各速度控制層的頂、底界層位進(jìn)行構(gòu)造解釋，對(duì)比追蹤獲得全區(qū)的各速度控制層的頂、底界層位和斷層；建立全區(qū)的地層框架模型；建立各速度控制層的速度概率密度分布函數(shù)和速度縱向變差函數(shù)；建立各速度控制層的速度橫向變差函數(shù)；在所述全區(qū)的地層框架模型的控制下，利用所述各速度控制層的速度概率密度分布函數(shù)、速度縱向變差函數(shù)和速度橫向變差函數(shù)進(jìn)行地質(zhì)統(tǒng)計(jì)學(xué)模擬，得到多個(gè)速度模型實(shí)現(xiàn)；對(duì)所述多個(gè)速度模型實(shí)現(xiàn)進(jìn)行平均計(jì)算，得到全地層速度模型。
【專利說明】
一種基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法

【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于石油天然氣地震勘探【技術(shù)領(lǐng)域】，更具體地講，涉及一種應(yīng)用于石油地震勘探的野外地震資料處理解釋的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法。

【背景技術(shù)】
[0002]準(zhǔn)確建立全地層速度模型是地震資料處理及解釋工作的核心任務(wù)之一，是獲取合理優(yōu)良地質(zhì)成果的關(guān)鍵工作內(nèi)容和重要前提條件。
[0003]在地震資料處理工作中，由于各類處理方法理論的假設(shè)條件的存在，使得處理人員為了獲得高質(zhì)量的成像而開展速度分析，借此獲得研究區(qū)的最初疊加速度模型。若將該疊加速度模型直接用作后續(xù)疊前時(shí)間或疊前深度偏移的初始速度模型，將帶來較多的偏移速度迭代次數(shù)，這不僅帶來繁復(fù)的工作量，而且容易造成局部速度迭代收斂方向的偏差。如果初始速度模型存在較大局部偏差或者錯(cuò)誤，速度迭代分析也可能無法消除趨勢(shì)誤差，從而降低局部區(qū)域偏移成像的質(zhì)量和精度，甚至在局部區(qū)域無法做到成果資料的保真保幅。
[0004]在地震資料構(gòu)造解釋工作中，全地層速度模型主要用于時(shí)深轉(zhuǎn)換。較為成熟的方法是以處理疊加或偏移速度為背景，利用研究區(qū)的測(cè)井資料進(jìn)行約束校正，并以較為大套的地震層位框架進(jìn)行控制，從而建立全區(qū)的層速度模型。該方法能基本滿足構(gòu)造解釋的需要，而且也不易發(fā)生構(gòu)造畸變。但是在巖性油氣藏及其它各類隱蔽油氣藏的勘探任務(wù)中，需要確定儲(chǔ)層的準(zhǔn)確頂?shù)咨疃?，?duì)于水平井的開發(fā)部署，儲(chǔ)層的準(zhǔn)確位置及深度則顯得更為重要，粗糙的層速度模型的精度遠(yuǎn)遠(yuǎn)無法滿足此類勘探任務(wù)的需要，而在儲(chǔ)層預(yù)測(cè)的地震反演工作中，對(duì)速度模型的精度要求更加不言而喻。隨著疊前深度偏移技術(shù)的進(jìn)步與日趨成熟，時(shí)深轉(zhuǎn)換工作已經(jīng)逐漸淡化，速度研究逐漸集中于偏移速度分析和地震反演的速度建模工作之中。
[0005]地層速度的本質(zhì)包含地層框架的控制因素，在偏移速度模型的生成過程中，多數(shù)情況下未加入地層框架控制，使得速度分析的約束控制條件單一化，僅從成像質(zhì)量和能量聚焦的因素進(jìn)行修整。尤其對(duì)于低信噪比地區(qū)，即使加入早期解釋地層框架，也極易造成縱向速度樣點(diǎn)求取的盲目性和速度橫向變化規(guī)律的凌亂。
[0006]鑒于上述難題，需要從另一個(gè)角度開展速度研究并建立全地層速度模型，以在地層框架控制下模擬地層速度并能較為準(zhǔn)確地反映地層速度的縱橫向變化規(guī)律。


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

[0007]針對(duì)現(xiàn)有技術(shù)中存在的問題，本發(fā)明的目的是解決上述技術(shù)問題中的一個(gè)或多個(gè)。
[0008]本發(fā)明的目的在于提供一種引入地層框架的約束以及各速度控制層地質(zhì)統(tǒng)計(jì)學(xué)特征并更加合理可靠的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法。
[0009]為了實(shí)現(xiàn)上述目的，本發(fā)明提供了一種基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法，所述方法包括以下步驟:利用研究區(qū)內(nèi)的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)確定預(yù)設(shè)數(shù)量的速度控制層；對(duì)疊后時(shí)間偏移地震數(shù)據(jù)進(jìn)行井震標(biāo)定并對(duì)各速度控制層的頂、底界層位進(jìn)行構(gòu)造解釋，對(duì)比追蹤獲得全區(qū)的各速度控制層的頂、底界層位和斷層；利用所述全區(qū)的各速度控制層的頂、底界層位和斷層建立全區(qū)的地層框架模型；利用所述研究區(qū)內(nèi)的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)建立各速度控制層的速度概率密度分布函數(shù)和速度縱向變差函數(shù)；在所述構(gòu)造解釋的基礎(chǔ)上對(duì)各速度控制層進(jìn)行地震屬性研究并獲得各速度控制層的地震屬性，利用所述各速度控制層的地震屬性以及各速度控制層的沉積相研究成果和巖相研究成果建立各速度控制層的速度橫向變差函數(shù)；在所述全區(qū)的地層框架模型的控制下，利用所述各速度控制層的速度概率密度分布函數(shù)、速度縱向變差函數(shù)和速度橫向變差函數(shù)進(jìn)行地質(zhì)統(tǒng)計(jì)學(xué)模擬，得到多個(gè)速度模型實(shí)現(xiàn)；對(duì)所述多個(gè)速度模型實(shí)現(xiàn)進(jìn)行平均計(jì)算，得到全地層速度模型。
[0010]根據(jù)本發(fā)明的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法的一個(gè)實(shí)施例，所述疊后時(shí)間偏移地震數(shù)據(jù)是對(duì)所述研究區(qū)的地震數(shù)據(jù)進(jìn)行常規(guī)預(yù)處理及疊后時(shí)間偏移處理得到的。
[0011]根據(jù)本發(fā)明的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法的一個(gè)實(shí)施例，利用對(duì)所述研究區(qū)的目的層進(jìn)行地震反演獲得的目的層反演速度對(duì)所述全地層速度模型進(jìn)行內(nèi)嵌替換，得到目的層速度更新的全地層速度模型。
[0012]根據(jù)本發(fā)明的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法的一個(gè)實(shí)施例，確定速度控制層的原則是同一個(gè)速度控制層的背景速度基本一致并且同一個(gè)速度控制層具有基本相同的沉積背景。
[0013]本發(fā)明的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法能夠適用于各類復(fù)雜構(gòu)造及低信噪比的地震工區(qū)，所獲得的全地層速度模型較為真實(shí)合理，符合研究區(qū)各地層的地質(zhì)特征，可用于地震時(shí)深轉(zhuǎn)換、確定性反演以及疊前偏移速度迭代的初始速度模型，尤其對(duì)隱蔽油氣藏的精細(xì)刻畫與評(píng)價(jià)，水平井的開發(fā)部署等工作具有重要的應(yīng)用意義。

【專利附圖】

【附圖說明】
[0014]圖1是示出根據(jù)本發(fā)明的示例性實(shí)施例的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法的流程圖。

【具體實(shí)施方式】
[0015]以下，參照附圖來詳細(xì)說明本發(fā)明的實(shí)施例。
[0016]可通過計(jì)算機(jī)軟件模塊實(shí)現(xiàn)本發(fā)明的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法。在本申請(qǐng)中描述的各個(gè)步驟不限制為上述步驟，其中的一些步驟可被進(jìn)一步拆分為更多的步驟，并且一些步驟可合并為較少的步驟。
[0017]對(duì)于某一研究區(qū)的某一地層，在其沉積環(huán)境及巖相變化不大的情況下，地層的彈性性質(zhì)具有較好的地質(zhì)統(tǒng)計(jì)學(xué)規(guī)律。根據(jù)地質(zhì)統(tǒng)計(jì)學(xué)理論，結(jié)合地層速度的概率密度分布函數(shù)和地層框架控制下的縱橫向變差函數(shù)，可以對(duì)地震采樣率下的地層樣點(diǎn)速度進(jìn)行統(tǒng)計(jì)學(xué)模擬，從而建立全區(qū)的全地層速度模型。因此，本發(fā)明基于以上思路提出了一種基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法，對(duì)于地層框架內(nèi)速度樣點(diǎn)的插值采用地質(zhì)統(tǒng)計(jì)學(xué)模擬算法取代以往的線性插值算法，更加符合地質(zhì)統(tǒng)計(jì)學(xué)規(guī)律，其橫向變差函數(shù)將地震數(shù)據(jù)以及其它地質(zhì)資料所表現(xiàn)的同一套地層的沉積學(xué)規(guī)律引入到該速度模型之中，縱向變差函數(shù)則將同一套地層的巖性組合的變化規(guī)律引入到速度模型中，使獲得的速度模型更加符合研究區(qū)的地質(zhì)規(guī)律。
[0018]圖1是示出根據(jù)本發(fā)明的示例性實(shí)施例的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法的流程圖。
[0019]參照?qǐng)D1，在步驟SlOl中，確定速度控制層。根據(jù)本發(fā)明的示例性實(shí)施例，利用研究區(qū)內(nèi)的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)確定預(yù)設(shè)數(shù)量的速度控制層。其中，速度控制層對(duì)應(yīng)的是這樣一套地層，即其地層的背景速度基本一致并且具有基本相同的沉積背景。具體地，利用研究區(qū)內(nèi)鉆井的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)來分析地層速度結(jié)構(gòu)及速度背景，并通過結(jié)合研究區(qū)的區(qū)域地層資料，確定縱向上合適數(shù)量的速度控制層。其中，上述的預(yù)設(shè)數(shù)量實(shí)際上就是所確定的合適數(shù)量。在確定速度控制層時(shí)，其基本原則是同一個(gè)速度控制層的背景速度基本一致并且同一個(gè)速度控制層具有基本相同的沉積背景，從而避免同一個(gè)速度控制層內(nèi)有較大的巖相變化，其中的“基本”的涵義是清楚的，因?yàn)樗俣瓤刂茖拥谋尘八俣群统练e背景不可能是完全相同的，只要基本一致或基本相同就可以劃分為同一個(gè)速度控制層。
[0020]其中，聲波測(cè)井?dāng)?shù)據(jù)是指利用聲學(xué)物理原理制造的測(cè)井儀器，由測(cè)井電纜下入井內(nèi)，使地面電測(cè)儀沿著井筒連續(xù)記錄的隨深度變化的聲學(xué)參數(shù)，聲波測(cè)井能測(cè)量地層的聲波速度。VSP (即垂直地震剖面法)測(cè)井?dāng)?shù)據(jù)(或稱為VSP地震記錄、VSP地震資料、VSP數(shù)據(jù))為地面激發(fā)井中接收的數(shù)據(jù)，即炮點(diǎn)設(shè)置在地面而檢波點(diǎn)設(shè)置在井中，其下行直達(dá)波受噪聲的影響較小，較地面地震記錄能更精確地反映地下介質(zhì)對(duì)振幅的衰減作用，其作用之一是能測(cè)量地層的地震波速度。
[0021]在確定速度控制層之后，在步驟S102中，利用上述聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)對(duì)疊后時(shí)間偏移地震數(shù)據(jù)進(jìn)行井震標(biāo)定，在步驟S103中，對(duì)各速度控制層的頂、底界層位進(jìn)行構(gòu)造解釋，對(duì)比追蹤獲得全區(qū)的各速度控制層的頂、底界層位和斷層。其中，疊后時(shí)間偏移地震數(shù)據(jù)是對(duì)研究區(qū)的地震數(shù)據(jù)進(jìn)行常規(guī)預(yù)處理及疊后時(shí)間偏移處理得到的。預(yù)處理是采用本領(lǐng)域的常規(guī)方法，例如將原始采集的地震數(shù)據(jù)經(jīng)過靜校正、去噪、振幅補(bǔ)償、動(dòng)校正、疊加以及偏移等處理，最終形成疊后或偏移后的地震數(shù)據(jù)，在此不再贅述。
[0022]進(jìn)行井震標(biāo)定是利用聲波測(cè)井制作的合成地震記錄和VSP測(cè)井處理得到的走廊疊加剖面對(duì)地震處理成果剖面進(jìn)行對(duì)比標(biāo)定，確定各速度控制層頂、底界對(duì)應(yīng)的地震層位位置及其地震響應(yīng)特性，進(jìn)行構(gòu)造解釋是根據(jù)井震標(biāo)定的結(jié)果分析并獲得各速度控制層頂、底界的地震層位位置和其地震響應(yīng)特征，結(jié)合區(qū)域地層及構(gòu)造研究成果或認(rèn)識(shí)，對(duì)比追蹤獲得全區(qū)的各速度控制層的頂、底界層位和斷層，井震標(biāo)定和構(gòu)造解釋均可以采用本領(lǐng)域公知的方法進(jìn)行。
[0023]在步驟S104中，利用全區(qū)的各速度控制層的頂、底界層位和斷層建立全區(qū)的地層框架模型。建立地層框架模型是利用構(gòu)造解釋獲得的地震層位和斷層，在地震網(wǎng)格控制下進(jìn)行網(wǎng)格化得到的一種框架模型，其中地震網(wǎng)格是指原始采集的地震數(shù)據(jù)經(jīng)過地震處理之后，由地震測(cè)網(wǎng)和地震采樣率構(gòu)成的三維網(wǎng)格化空間，建立地層框架模型可以采用本領(lǐng)域公知的方法進(jìn)行。
[0024]在步驟S105和S106中，利用研究區(qū)內(nèi)的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)建立各速度控制層的速度概率密度分布函數(shù)和速度縱向變差函數(shù)。
[0025]其中，概率密度分布函數(shù)是一個(gè)描述某個(gè)隨機(jī)變量的輸出值在某個(gè)確定的取值點(diǎn)附近的可能性的函數(shù)，例如，對(duì)于一維實(shí)隨機(jī)變量X，設(shè)它的累積分布函數(shù)是Fx(X)，如果存在可測(cè)函數(shù)fx(x)滿足式1:
[0026]V —oc<6/<.x, Fx (a) = | j\ (x)Jx-式 I
[0027]那么X是一個(gè)連續(xù)型隨機(jī)變量，并且fx(x)就是它的概率密度分布函數(shù)。
[0028]具體地，概率密度分布函數(shù)的含義可以理解為，fx(x)是X落在X處“單位寬度”內(nèi)的概率，速度的概率密度分布函數(shù)描述了速度樣點(diǎn)值落在某一速度值附近的概率，某一速度控制層的速度概率密度分布函數(shù)的建立方法是利用該速度控制層內(nèi)的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)樣點(diǎn)進(jìn)行先驗(yàn)概率密度的計(jì)算，再利用某種理論概率密度分布函數(shù)對(duì)先驗(yàn)概率密度進(jìn)行匹配，并通過調(diào)整理論函數(shù)的參數(shù)使其與先驗(yàn)概率密度樣點(diǎn)達(dá)到較好的擬合度。
[0029]變差函數(shù)是Motheron (1965)年提出的一種矩估計(jì)方法,為區(qū)域化變量的增量平方的數(shù)學(xué)期望，也就是區(qū)域化變量的增量的方差，它反映了區(qū)域化變量在某個(gè)方向上某一距離范圍內(nèi)的變化程度，某一速度控制層的速度縱向變差函數(shù)描述了速度在某一速度控制層內(nèi)縱向上的變化程度，速度縱向變差函數(shù)建立是利用該速度控制層內(nèi)的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)樣點(diǎn)進(jìn)行先驗(yàn)變差值的計(jì)算，再利用某種理論變差函數(shù)對(duì)先驗(yàn)變差值進(jìn)行匹配，通過調(diào)整變差函數(shù)的參數(shù)使其與先驗(yàn)變差值達(dá)到較好的擬合度，其中，縱向變差函數(shù)的參數(shù)包括函數(shù)組分及縱向變程。
[0030]在步驟S107中，在構(gòu)造解釋的基礎(chǔ)上對(duì)各速度控制層進(jìn)行地震屬性研究并獲得各速度控制層的地震屬性，利用所述各速度控制層的地震屬性以及各速度控制層的沉積相研究成果和巖相研究成果建立各速度控制層的速度橫向變差函數(shù)。
[0031]其中，地震屬性主要是以能夠反映巖性變化的屬性為主，例如振幅屬性，它能夠反映某一速度控制層速度的橫向變化，因此某一速度控制層的速度橫向變差函數(shù)可以用該速度控制層的地震屬性的橫向變差函數(shù)來表征，根據(jù)上述變差函數(shù)的涵義，某一速度控制層的速度橫向變差函數(shù)描述了速度在某一速度控制層內(nèi)橫向上的變化程度，其橫向變差函數(shù)的建立是利用該速度控制層的地震屬性數(shù)據(jù)樣點(diǎn)進(jìn)行先驗(yàn)變差函數(shù)值的計(jì)算，再利用某種理論變差函數(shù)對(duì)先驗(yàn)變差值進(jìn)行匹配，通過調(diào)整變差函數(shù)的參數(shù)使其與先驗(yàn)變差值達(dá)到較好的擬合度，其中，橫向變差函數(shù)的參數(shù)包括函數(shù)組分及橫向變程，橫向變程的確定需要兼顧該速度控制層的沉積相研究成果和巖相研究成果。
[0032]其中，在建立各速度控制層的速度橫向變差函數(shù)時(shí)，所利用的地震屬性主要是以能夠反映巖性變化的屬性為主，例如振幅屬性等；而各速度控制層的沉積相研究成果和巖相研究成果則是已有的并且可以采用本領(lǐng)域的現(xiàn)有方法獲得。利用上述地震屬性以及各速度控制層的沉積相研究成果和巖相研究成果就可以確定各速度控制層的速度橫向變差函數(shù)的參數(shù)，如橫向變程、函數(shù)組成類型等，從而建立適用于研究區(qū)各速度控制層的速度橫向變差函數(shù)。
[0033]總之，上述各速度控制層的速度概率密度分布函數(shù)、速度縱向變差函數(shù)的建立是本領(lǐng)域技術(shù)人員可以在各速度控制層的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)的基礎(chǔ)上采用公知的理論和方法獲得；各速度控制層的速度橫向變差函數(shù)的建立亦是本領(lǐng)域技術(shù)人員可以在各速度控制層的地震屬性研究結(jié)果及其沉積相、巖相研究成果的基礎(chǔ)上采用公知的理論和方法獲得。
[0034]在步驟S108中，在全區(qū)的地層框架模型的控制下，利用各速度控制層的速度概率密度分布函數(shù)、速度縱向變差函數(shù)和速度橫向變差函數(shù)進(jìn)行地質(zhì)統(tǒng)計(jì)學(xué)模擬，得到多個(gè)速度模型實(shí)現(xiàn)。其中，多個(gè)速度模型實(shí)現(xiàn)是等概率的并且是通過定義模擬得到。
[0035]其中，地質(zhì)統(tǒng)計(jì)學(xué)模擬主要基于序貫指示模擬類型算法或馬爾科夫鏈-蒙特卡羅算法，后者更為先進(jìn)。對(duì)于某一速度控制層速度的地質(zhì)統(tǒng)計(jì)學(xué)模擬，馬爾科夫鏈-蒙特卡羅算法根據(jù)概率密度分布函數(shù)獲得統(tǒng)計(jì)意義上正確的樣點(diǎn)集，即根據(jù)概率密度分布函數(shù)能夠得到多個(gè)等概率速度組合的結(jié)果，縱向變差函數(shù)約束馬爾科夫鏈-蒙特卡羅模擬算法的縱向分辨率和概率權(quán)重，橫向變差函數(shù)控制橫向未知樣點(diǎn)處速度的模擬結(jié)果，其控制方式可以是依據(jù)已知樣點(diǎn)到未知樣點(diǎn)處速度的橫向變差，對(duì)未知樣點(diǎn)處速度的概率密度函數(shù)和縱向變差函數(shù)的進(jìn)行統(tǒng)計(jì)學(xué)定義。地質(zhì)統(tǒng)計(jì)學(xué)模擬的實(shí)現(xiàn)方式要繁雜許多，但基本原理相同，許多地質(zhì)統(tǒng)計(jì)學(xué)著作和論文中已有詳細(xì)論述，也已產(chǎn)生多個(gè)商業(yè)化計(jì)算機(jī)軟件模塊，在此不再贅述。
[0036]在步驟S109中，對(duì)多個(gè)速度模型實(shí)現(xiàn)進(jìn)行平均計(jì)算，得到全地層速度模型。根據(jù)上述地質(zhì)統(tǒng)計(jì)學(xué)模擬算法的基本原理，模擬算法根據(jù)概率密度分布函數(shù)能夠得到多個(gè)等概率速度組合結(jié)果，一個(gè)等概率速度組合結(jié)果即為一個(gè)速度模型實(shí)現(xiàn)，對(duì)多個(gè)等概率速度模型實(shí)現(xiàn)進(jìn)行平均計(jì)算，可以降低地質(zhì)統(tǒng)計(jì)學(xué)模擬的多解性，更加接近真實(shí)結(jié)果，一般地，定義越多個(gè)速度模型實(shí)現(xiàn)，其平均化結(jié)果越接近真實(shí)結(jié)果。
[0037]如果已對(duì)研究區(qū)的目的層進(jìn)行了地震反演工作，則可以利用對(duì)研究區(qū)的目的層進(jìn)行地震反演獲得的目的層反演速度對(duì)上述全地層速度模型進(jìn)行內(nèi)嵌替換，從而可以得到目的層速度更新的全地層速度模型。
[0038]其中，內(nèi)嵌替換操作的實(shí)現(xiàn)方法可以是，將所述目的層的頂?shù)讓游欢x為時(shí)窗以控制需要進(jìn)行內(nèi)嵌替換的數(shù)據(jù)范圍，在時(shí)窗范圍之內(nèi)VnOT = Vinvwsim，在時(shí)窗范圍之外VnOT=Vmiginal，其中VOTiginal是多個(gè)速度模型實(shí)現(xiàn)進(jìn)行平均計(jì)算得到的原始的全地層速度模型，Vinvers1n是目的層地震反演速度體，Vnew即為目的層速度更新的全地層速度模型。
[0039]以本發(fā)明所闡述的方法或思路建立的全地層速度模型能夠較好地體現(xiàn)研究區(qū)各地層的地質(zhì)學(xué)特征，保證速度模型趨勢(shì)的真實(shí)性，該全地層速度模型不僅能為各類確定性反演提供合理的初始速度模型，也能作為疊前偏移速度分析的初始輸入速度，使偏移速度分析迭代次數(shù)減少并向正確的方向迭代收斂，從而提高偏移處理成果資料的保真度。
[0040]因此，本發(fā)明基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法以地質(zhì)統(tǒng)計(jì)學(xué)模擬算法為基礎(chǔ)建立全地層速度模型，建模過程中地層框架的約束以及各控制層地質(zhì)統(tǒng)計(jì)學(xué)特征的引入，使獲得的速度模型更加符合研究區(qū)的地質(zhì)規(guī)律；目的層反演速度的內(nèi)嵌替換，能夠提高目的層儲(chǔ)層深度及形態(tài)的刻畫精度，本發(fā)明比傳統(tǒng)的基于線性插值的全地層層速度建模方法更加合理可靠。
[0041]盡管已參照優(yōu)選實(shí)施例表示和描述了本發(fā)明，但本領(lǐng)域技術(shù)人員應(yīng)該理解，在不脫離由權(quán)利要求限定的本發(fā)明的精神和范圍的情況下，可以對(duì)這些實(shí)施例進(jìn)行各種修改和變換。
【權(quán)利要求】
1.一種基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法，其特征在于，所述方法包括以下步驟: 利用研究區(qū)內(nèi)的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)確定預(yù)設(shè)數(shù)量的速度控制層； 對(duì)疊后時(shí)間偏移地震數(shù)據(jù)進(jìn)行井震標(biāo)定并對(duì)各速度控制層的頂、底界層位進(jìn)行構(gòu)造解釋，對(duì)比追蹤獲得全區(qū)的各速度控制層的頂、底界層位和斷層； 利用所述全區(qū)的各速度控制層的頂、底界層位和斷層建立全區(qū)的地層框架模型； 利用所述研究區(qū)內(nèi)的聲波測(cè)井及VSP測(cè)井的速度數(shù)據(jù)建立各速度控制層的速度概率密度分布函數(shù)和速度縱向變差函數(shù)； 在所述構(gòu)造解釋的基礎(chǔ)上對(duì)各速度控制層進(jìn)行地震屬性研究并獲得各速度控制層的地震屬性，利用所述各速度控制層的地震屬性以及各速度控制層的沉積相研究成果和巖相研究成果建立各速度控制層的速度橫向變差函數(shù)； 在所述全區(qū)的地層框架模型的控制下，利用所述各速度控制層的速度概率密度分布函數(shù)、速度縱向變差函數(shù)和速度橫向變差函數(shù)進(jìn)行地質(zhì)統(tǒng)計(jì)學(xué)模擬，得到多個(gè)速度模型實(shí)現(xiàn)； 對(duì)所述多個(gè)速度模型實(shí)現(xiàn)進(jìn)行平均計(jì)算，得到全地層速度模型。
2.根據(jù)權(quán)利要求1所述的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法，其特征在于，所述疊后時(shí)間偏移地震數(shù)據(jù)是對(duì)所述研究區(qū)的地震數(shù)據(jù)進(jìn)行預(yù)處理及疊后時(shí)間偏移處理得到的。
3.根據(jù)權(quán)利要求1所述的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法，其特征在于，利用對(duì)所述研究區(qū)的目的層進(jìn)行地震反演獲得的目的層反演速度對(duì)所述全地層速度模型進(jìn)行內(nèi)嵌替換，得到目的層速度更新的全地層速度模型。
4.根據(jù)權(quán)利要求1所述的基于地質(zhì)統(tǒng)計(jì)學(xué)模擬的全地層速度建模方法，其特征在于，確定速度控制層的原則是同一個(gè)速度控制層的背景速度基本一致并且同一個(gè)速度控制層具有基本相同的沉積背景。
【文檔編號(hào)】G06F19/00GK104200115SQ201410465201
【公開日】2014年12月10日 申請(qǐng)日期:2014年9月12日 優(yōu)先權(quán)日:2014年9月12日
【發(fā)明者】章雄, 李樂, 龐崇友, 彭才, 刁永波 申請(qǐng)人:中國石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司