專利名稱：用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法及裝置的制作方法
技術領域：
本發(fā)明涉及地震資料層位精細解釋與地震反演所用地質(zhì)小層數(shù)據(jù)的處理技術，是一種用于地震資料解釋的基于聲波測井曲線的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法。
背景技術：
據(jù)統(tǒng)計，我國油田絕大多數(shù)儲層為陸相碎屑巖沉積，其中內(nèi)陸河流-三角洲和沖積扇-扇三角洲沉積的儲量占總量的90%以上。這類沉積呈多旋回性，使得我國陸上油田縱向上地質(zhì)小層多，有的多達數(shù)十層甚至百余層。這樣的地質(zhì)條件決定我國油田開發(fā)通常以地質(zhì)小層為單元，長期油田開發(fā)生產(chǎn)實踐形成了一系列以小層為單元的技術和方法，地層劃分對比、沉積相分析、堵水調(diào)剖、分層注水和油水井注采關系完善等油田開發(fā)措施，都需要依托對各個地質(zhì)小層斷層、構(gòu)造和砂體分布情況的翔實認識。 盡管在油藏開發(fā)階段，特別是高含水期，井網(wǎng)密度較大，但井間仍然有很大的不確定性，僅靠井數(shù)據(jù)分析推測和插值、模擬，無法解決井間的地質(zhì)認識問題，必須結(jié)合具有明確井間信息的地震資料加以研究。在油田開發(fā)階段應用的地震技術被稱為開發(fā)地震。基于上述油田開發(fā)以地質(zhì)小層為單元的實際情況，開發(fā)地震技術在油田開發(fā)領域服務于油田開發(fā)開展油藏描述研究，同樣需要以地質(zhì)小層為單元進行精細成圖。實現(xiàn)開發(fā)地震技術以地質(zhì)小層為單元進行成圖，前提條件是必須做到在時間域?qū)Φ卣鹳Y料能夠進行有效的地質(zhì)小層解釋。然而，我國陸相儲層多數(shù)是典型的砂泥巖薄互層儲層，這類儲層砂泥巖均較薄，單層厚度遠小于四分之一地震波長。地震反射波組是若干個地質(zhì)小層反射相互干涉的結(jié)果，在地震剖面上，就某一具體地質(zhì)小層而言，一般不具有單獨的連續(xù)可追蹤的地震反射特征。正是由于這個原因，在開展地震資料層位解釋時，人們通常只能解釋幾個大套層系的頂?shù)?，即所謂地震標志層，而無法對地質(zhì)小層進行連續(xù)追蹤解釋。解決上述問題，最簡單及最有效的方法是通過深時轉(zhuǎn)換，將地質(zhì)小層數(shù)據(jù)由深度域轉(zhuǎn)換到時間域，在時間域進行標志層約束插值，生成時間域地質(zhì)小層層面數(shù)據(jù)，再將時間域地質(zhì)小層層面數(shù)據(jù)加載到地震解釋系統(tǒng)，實現(xiàn)對地質(zhì)小層的有效解釋。而實現(xiàn)該方法的關鍵是深時轉(zhuǎn)換。馬勁鳳(馬勁鳳，聲波測井與地震數(shù)據(jù)匹配的新方法.西北大學學報(自然科學版)，1998，28 (3) :229 232)于1998年給出了深時轉(zhuǎn)換公式tm= Σ (2 Δ h/Vj)，其中tm為聲波到達第η層時的旅行時，Δ h為深度域采樣間隔，Vj為第j小層的速度。但是該公式假設起始時間和深度值均為0，這與實際情況不符，通常測井的有效值深度和地震基準面海拔并不一致。另外一個問題是這種計算方法沒有考慮測井和地震兩種技術存在的速度誤差，以及誤差的積累問題。單剛義、韓立國、張麗華等(單剛義、韓立國、張麗華等.子波整形提高合成地震記錄質(zhì)量.吉林大學學報(地球科學版)，2006，36增刊182 184)于2006年給出了如下深時轉(zhuǎn)換計算公式
HFa = T0 + 2 j Δ/(/ ν/ X 101，其中Htl為起始深度、H為終止深度、Ttl為起始深度對應的時間、Th為深度h所對應的時間、AtOO為深度h處的聲波時差值。盡管該方法考慮了起始深度和時間，克服了馬勁鳳方法的部分問題，但是仍然沒有考慮到測井和地震兩種技術存在的速度誤差，以及誤差的積累問題。經(jīng)過四、五十年的開采，我國已開發(fā)油田總體上進入了高含水后期。高含水油田開發(fā)地震技術需要以地質(zhì)小層為成圖單元。實現(xiàn)開發(fā)地震技術以小層為單元進行成圖，前提條件是必須做到在時間域?qū)Φ卣鹳Y料能夠進行有效的地質(zhì)小層解釋。但是我國陸相儲層多為砂泥巖薄互層，地震反射波組是若干個地質(zhì)小層反射相互干涉的結(jié)果。在地震剖面上，通常僅能解釋少數(shù)大的標志層，絕大多數(shù)地質(zhì)小層無法進行連續(xù)追蹤解釋。然而，盡管地質(zhì)小層在地震剖面上無法直接解釋，但是客觀上地質(zhì)小層深度值畢竟對應一個地震波旅行時間，實現(xiàn)地震資料地質(zhì)小層的有效解釋，研發(fā)地質(zhì)小層深時轉(zhuǎn)換方法是關鍵。需要指出的是，目前人們提出的深時轉(zhuǎn)換方法，都是將測井曲線由深度域轉(zhuǎn)換到 時間域的方法，主要是用于解決井震匹配的問題，并沒有涉及地質(zhì)小層深時轉(zhuǎn)換的方法。如何將地質(zhì)小層由深度域轉(zhuǎn)換到時間域，需要研發(fā)精度更高的用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換(由深度域轉(zhuǎn)換到時間域)方法。

發(fā)明內(nèi)容
本發(fā)明的目的就是提供一種用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法，針對現(xiàn)有技術開展深時轉(zhuǎn)換未考慮測井和地震兩種技術存在速度誤差，以及誤差的積累等問題，采用上下地震標志層約束和采樣點誤差分擔的方法，以聲波時差測井曲線提供的深度采樣點間隔和聲波時差平均速度逐點遞推計算，實現(xiàn)高精度地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換。該方法的建立能夠?qū)⒌刭|(zhì)小層數(shù)據(jù)由深度域轉(zhuǎn)換到時間域，在時間域?qū)Ω骶c的地質(zhì)小層數(shù)據(jù)進行插值，生成時間域地質(zhì)小層層面，即可快捷高效地實現(xiàn)時間域地質(zhì)小層的有效解釋。本發(fā)明實施例提供一種用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法，包括提取初始值，將位于目標層系上部最近的地震標志層位作為計算的起始位置，獲得該地震解釋層位與地質(zhì)層位的對應關系，如地震層位Tl對應地質(zhì)層位S2-1，將Tl地震解釋結(jié)果插值，形成網(wǎng)格數(shù)據(jù)，求取Tl網(wǎng)格數(shù)據(jù)中與待求井位最近點的時間值，與S2-1地質(zhì)層位深度值一起，形成起始位置的時間值和深度值對應的數(shù)據(jù)，即Ttl和Htl ;獲取地質(zhì)小層深度值，從油田地質(zhì)分層數(shù)據(jù)表中獲取要計算的地質(zhì)小層深度值Hn，作為要進行深時轉(zhuǎn)換的目標數(shù)值；深時轉(zhuǎn)換計算，以聲波時差測井曲線提供的深度采樣點間隔和聲波時差平均速度，從起始深度Htl向下逐點遞推計算，直至深度值扎處，求取目標地質(zhì)小層對應的時間值Tn;誤差校正，考慮到測井與地震兩種技術存在速度誤差和誤差的積累問題，采取上述方法計算至目標層系下方最近的地震解釋層位，比較采取測井速度計算的時間值和實際地震波旅行時間值存在的誤差，采取采樣點誤差平均分擔的方法，開展誤差校正；將校正后的結(jié)果輸出，獲得最終的地質(zhì)小層時間值。本發(fā)明實施例還提供一種用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法裝置，其特征在于，所述裝置包括提取初始值單元，用于提取時間和深度的初始值數(shù)據(jù)；獲取地質(zhì)小層深度值單元，用于獲取待深時轉(zhuǎn)換的地質(zhì)小層深度值；深時轉(zhuǎn)換計算單元，用于將要深時轉(zhuǎn)換的地質(zhì)小層的深度值轉(zhuǎn)換為時間值；
誤差校正單元，用于對深時轉(zhuǎn)換計算單元計算的時間值進行測井和地震誤差校正；結(jié)果輸出單元，用于將校正后的結(jié)果輸出，獲得最終的地質(zhì)小層時間值。本發(fā)明實施例用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法及裝置，能夠自動計算，將地質(zhì)小層數(shù)據(jù)由深度域轉(zhuǎn)換到時間域，在時間域?qū)Ω骶c時深轉(zhuǎn)換后的地質(zhì)小層數(shù)據(jù)進行插值，便可快捷高效地實現(xiàn)時間域地質(zhì)小層的有效解釋。該方法具有較強的準確性，而且實現(xiàn)了多口井的自動批量化處理，較采取通過制作合成地震記錄逐口井進行時深轉(zhuǎn)換的方法極大地提高了工作效率，節(jié)省了科技人員寶貴的人力。


此處所說明的附圖用來提供對本發(fā)明的進一步理解，構(gòu)成本申請的一部分，并不構(gòu)成對本發(fā)明的限定。在附圖中 圖I為本發(fā)明實施例的用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法的流程圖；圖2為本發(fā)明實施例的基于測井曲線的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換裝置的結(jié)構(gòu)示意圖；圖3為本發(fā)明實施例的用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法針對問題的示意圖；圖4為本發(fā)明實施例的用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法的效果圖；圖5為基于本發(fā)明實施例的用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法繪制的地質(zhì)小層砂體分布圖。
具體實施例方式為使本發(fā)明的目的、技術方案和優(yōu)點更加清楚明白，下面結(jié)合附圖對本發(fā)明實施例做進一步詳細說明。在此，本發(fā)明的示意性實施例及其說明用于解釋本發(fā)明，但并不作為對本發(fā)明的限定。圖I為本發(fā)明實施例的用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法的流程圖,如圖I所示,所述方法包括步驟S101，提取時間和深度的初始值。本實施例中，將位于目標層系上部最近的地震標志層位作為計算的起始位置，獲得該地震解釋層位與地質(zhì)層位的對應關系，如圖3中，地震層位Tl對應地質(zhì)層位S2-1，將Tl地震解釋結(jié)果插值，形成網(wǎng)格數(shù)據(jù)，求取Tl網(wǎng)格數(shù)據(jù)中與待求井位最近點的時間值，與S2-1地質(zhì)層位深度值一起，形成起始位置的時間值和深度值對應的數(shù)據(jù)，即Ttl和Η”步驟S102，獲取待深時轉(zhuǎn)換的地質(zhì)小層深度值。本實施例中，從油田地質(zhì)分層數(shù)據(jù)表中獲取要計算的地質(zhì)小層深度值Hn，作為要進行深時轉(zhuǎn)換的目標數(shù)值。地質(zhì)分層數(shù)據(jù)表為一個數(shù)據(jù)庫文件，其結(jié)構(gòu)主要包括jh表示井號、X表示當前層井位橫坐標、y表示當前層井位縱坐標、xcfzmc表示小層分組名稱、xcfzds表示小層分組頂深,如下表一所示。表一地質(zhì)分層數(shù)據(jù)庫
權利要求
1.一種用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法，其特征在于，所述方法包括 將位于目標層系上部最近的地震標志層位作為計算的起始位置，提取時間和深度初始值； 從油田地質(zhì)分層數(shù)據(jù)表中獲取要計算的地質(zhì)小層深度值，作為要進行深時轉(zhuǎn)換的目標深度數(shù)值； 采用聲波時差曲線，從起始深度向下逐點遞推進行深時轉(zhuǎn)換計算，直至目標深度，求取目標地質(zhì)小層對應的時間值； 以目標層系下方地震解釋層位為參考，比較采取測井速度計算的時間值和實際地震波旅行時間值存在的誤差，采取誤差平均方法開展誤差校正，將校正后的結(jié)果輸出，獲得最終的地質(zhì)小層時間值。
2.如權利要求I所述的方法，其特征在于，所述提取時間和深度初始值，包括 將位于目標層系上部最近的地震標志層位作為計算的起始位置，將起始層位地震解釋結(jié)果插值，形成網(wǎng)格數(shù)據(jù)，求取該網(wǎng)格數(shù)據(jù)中與待求井位最近點的時間值，與對應地質(zhì)層位深度值一起，構(gòu)成時間Ttl和深度初始值
3.如權利要求I所述的方法，其特征在于，所述獲取要計算的地質(zhì)小層深度值，包括 從油田地質(zhì)分層數(shù)據(jù)表中獲取要計算的地質(zhì)小層深度值Hn，作為要進行深時轉(zhuǎn)換的目標數(shù)值。
4.如權利要求I所述的方法，其特征在于，所述采用聲波時差曲線從起始深度向下逐點遞推進行深時轉(zhuǎn)換計算，包括 以聲波時差測井曲線提供的深度采樣點間隔和聲波時差平均速度，從起始深度Htl采用如下公式向下逐點遞推計算，直至深度值Hn處，求取目標地質(zhì)小層對應的時間值Tn ； 深時轉(zhuǎn)換的公式為
5.如權利要求I所述的方法，其特征在于，所述采取誤差平均方法開展誤差校正，包括 采取上述深時轉(zhuǎn)換公式計算至目標層系下方的地震解釋層位，比較采取測井速度計算的時間值和實際地震波旅行時間值存在的誤差，采取采樣點誤差平均分擔的方法，開展誤差校正，將校正后的結(jié)果輸出，獲得最終的地質(zhì)小層時間值； 采樣點誤差平均分擔校正公式為
6.一種用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法裝置，其特征在于，所述裝置包括 提取初始值單元，用于提取時間和深度的初始值數(shù)據(jù)； 獲取地質(zhì)小層深度值單元，用于獲取待深時轉(zhuǎn)換的地質(zhì)小層深度值； 深時轉(zhuǎn)換計算單元，用于將要深時轉(zhuǎn)換的地質(zhì)小層的深度值轉(zhuǎn)換為時間值； 誤差校正單元，用于對深時轉(zhuǎn)換計算單元計算的時間值進行測井和地震誤差校正； 結(jié)果輸出單元，用于將校正后的結(jié)果輸出，獲得最終的地質(zhì)小層時間值。
7.如權利要求6所述的裝置，其特征在于，所述的提取初始值單元用于提取時間和深度的初始值數(shù)據(jù)，包括 將位于目標層系上部最近的地震標志層位作為計算的起始位置，將起始層位地震解釋結(jié)果插值，形成網(wǎng)格數(shù)據(jù)，求取該網(wǎng)格數(shù)據(jù)中與待求井位最近點的時間值，與對應地質(zhì)層位深度值一起，構(gòu)成時間和深度初始值，即Ttl和
8.如權利要求6所述的裝置，其特征在于，所述的獲取地質(zhì)小層深度值單元用于獲取待深時轉(zhuǎn)換的地質(zhì)小層深度值，包括 從油田地質(zhì)分層數(shù)據(jù)表中獲取要計算的地質(zhì)小層深度值Hn，作為要進行深時轉(zhuǎn)換的目標數(shù)值。
9.如權利要求6所述的裝置，其特征在于，所述的深時轉(zhuǎn)換計算單元用于將要深時轉(zhuǎn)換的地質(zhì)小層的深度值轉(zhuǎn)換為時間值，包括 以聲波時差測井曲線提供的深度采樣點間隔和聲波時差平均速度，從起始深度Htl采用如下公式向下逐點遞推計算，直至深度值Hn處，求取目標地質(zhì)小層對應的時間值Tn ； 深時轉(zhuǎn)換的公式為
10.如權利要求6所述的裝置，其特征在于，所述的誤差校正單元用于計算的時間值進行誤差校正，包括 采取上述深時轉(zhuǎn)換公式計算至目標層系下方的地震解釋層位，比較采取測井速度計算的時間值和實際地震波旅行時間值存在的誤差，采取采樣點誤差平均分擔的方法，開展誤差校正； 采樣點誤差平均分擔校正公式為
全文摘要
本發(fā)明提供一種用于地震資料解釋的地質(zhì)小層數(shù)據(jù)深時轉(zhuǎn)換方法及裝置，所述方法包括將位于目標層系上部最近的地震標志層位作為計算的起始位置，提取時間和深度初始值；從油田地質(zhì)分層數(shù)據(jù)表中獲取要計算的地質(zhì)小層深度值，作為要進行深時轉(zhuǎn)換的目標深度數(shù)值；采用聲波時差曲線，從起始深度向下逐點遞推進行深時轉(zhuǎn)換計算，直至目標深度，求取目標地質(zhì)小層對應的時間值；以目標層系下方地震解釋層位為參考，比較采取測井速度計算的時間值和實際地震波旅行時間值存在的誤差，采取誤差平均方法開展誤差校正，將校正后的結(jié)果輸出，獲得最終的地質(zhì)小層時間值。
文檔編號G01V1/28GK102967881SQ201210450988
公開日2013年3月13日 申請日期2012年11月12日 優(yōu)先權日2012年11月12日
發(fā)明者劉文嶺 申請人:中國石油天然氣股份有限公司