專利名稱：測(cè)井聲波時(shí)差信號(hào)校正方法
技術(shù)領(lǐng)域：
本發(fā)明涉及地質(zhì)和石油開發(fā)，具體涉及一種測(cè)井聲波時(shí)差信號(hào)校正方法。
背景技術(shù)：
對(duì)地震資料解釋者而言，聲波時(shí)差曲線是最重要的測(cè)井曲線，它可以用來建立時(shí)間一深度關(guān)系，進(jìn)而聯(lián)通時(shí)間域地震資料和深度域測(cè)井、地質(zhì)資料；它是構(gòu)成波阻抗的重要因子，可以用于測(cè)井資料約束的測(cè)井、地震、地質(zhì)聯(lián)合反演工作。聲波時(shí)差曲線的質(zhì)量直接影響地震資料標(biāo)定、解釋質(zhì)量和測(cè)井、地震聯(lián)合反演精度。塔河油田的鉆井在鉆探石炭系、三疊系砂泥巖地層序列時(shí)，普遍遇到泥巖段井徑擴(kuò)大引起的聲波時(shí)差測(cè)井曲線嚴(yán)重失真的問題；在鉆探奧陶系碳酸鹽巖地層時(shí)，也遇到放空、漏失段聲波時(shí)差嚴(yán)重失真問題，主要表現(xiàn)在聲波時(shí)差曲線受周波跳躍、時(shí)差展寬影響， 數(shù)值增大畸變。受環(huán)境影響聲波時(shí)差測(cè)井曲線失真是一個(gè)較普遍問題，不同地區(qū)聲波時(shí)差曲線受環(huán)境影響程度是不同的，因而，工業(yè)化測(cè)井解釋軟件包內(nèi)沒有適用的環(huán)境校正方法或程序?qū)Υ诉M(jìn)行校正。地震資料解釋者只能應(yīng)用手工編輯的方法對(duì)聲波時(shí)差進(jìn)行主觀拉伸一壓縮、方波化或削截修改。針對(duì)聲波時(shí)差野值段進(jìn)行削截編輯可以消除不合理的大反射系數(shù)， 拉伸-壓縮手段可以改變合成地震道與井旁地震道的相似度。但表面問題的解決引發(fā)了更深層次的問題編輯后聲波時(shí)差曲線對(duì)應(yīng)的速度成果(特別是層速度)存在突變、不能用于下步變速成圖工序；方波化降低了測(cè)井資料的垂向分辨率，方波化改變了反射系數(shù)的相對(duì)變化、降低了測(cè)井資料與地震資料的匹配關(guān)系，將引起測(cè)井、地震約束反演多解性的增加。同時(shí)，塔河油田碳酸鹽巖縫洞型油藏具有儲(chǔ)集空間類型多樣化、儲(chǔ)集體連通性差、 流體流動(dòng)狀態(tài)及油水關(guān)系復(fù)雜的特點(diǎn)，在國(guó)外很少見，這類油藏采用衰竭式開發(fā)采收率很低，采收率只有10%，且產(chǎn)量遞減快；由于對(duì)縫洞體發(fā)育規(guī)律和油藏流體流動(dòng)特征的認(rèn)識(shí)程度低，國(guó)內(nèi)外在縫洞型油藏開發(fā)方面還沒有形成相應(yīng)的開發(fā)理論和開發(fā)方法。做好聲波時(shí)差環(huán)境校正將產(chǎn)生準(zhǔn)確的速度模型、將為測(cè)井地震聯(lián)合反演、儲(chǔ)層描述提供最基礎(chǔ)資料， 將為碳酸鹽巖孔隙度計(jì)算模型、油藏模型的建立提供最準(zhǔn)確的素材。本發(fā)明要解決的問題是提升對(duì)縫洞型油藏的認(rèn)識(shí)、提高油藏采收率、實(shí)現(xiàn)油田快速開發(fā)迫切要解決的關(guān)鍵問題。

發(fā)明內(nèi)容
本發(fā)明需要解決的技術(shù)問題是，如何提供一種測(cè)井聲波時(shí)差信號(hào)校正方法，實(shí)現(xiàn)對(duì)聲波時(shí)差畸變段的自動(dòng)修復(fù)校正，提高測(cè)井聲波時(shí)差資料的質(zhì)量。本發(fā)明的技術(shù)問題這樣解決構(gòu)建一種測(cè)井聲波時(shí)差信號(hào)校正方法，其特征在于， 包括以下步驟1. 1)獲取測(cè)井聲波時(shí)差信號(hào)并按實(shí)際情況設(shè)置當(dāng)前環(huán)境；1. 2)根據(jù)當(dāng)前環(huán)境的巖性為砂巖(SP < 90)，自動(dòng)根據(jù)遞歸公式ACRD = ACRDkl^log (RD)+ACRDkO R2 = 29. 57，24410 點(diǎn)[1]
ACRS = ACRSkl^log(RS)+ACRSkO R2 = 28. 18,23867 點(diǎn)[2]ACDEN = ACDENkl*DEN+ACDENkO R2 = 29. 57，16301 點(diǎn)[3]ACCNL = ACCNLkl*CNL+ACCNLkO R2 = 20. 01 18583 點(diǎn) W]AC 平均=Prd*ACRD+Prs*ACRS+Pden*ACDEN+Pcnl*ACCNL [5]識(shí)別畸變點(diǎn)，實(shí)測(cè)聲波時(shí)差與AC平均差值大于標(biāo)準(zhǔn)偏差為畸變點(diǎn)，畸變點(diǎn)按公式 [5]對(duì)其校正。式中SP為自然電位，RD為深電阻率，RS為淺電阻率，DEN為密度，CNL為中子，ACRD、ACRS, ACDEN、ACCNL 分別為據(jù) RD、RS、DEN、CNL 遞歸的聲波時(shí)差，Prd, Prs, Pden, Pcnl分別為RD、RS、DEN、CNL的權(quán)重，R2為遞歸方程的相關(guān)性指標(biāo)，AC平均為遞歸平均值， ACRDkl、ACRDkO等后兩個(gè)字節(jié)為kl、k0的參數(shù)分別代表遞歸方程的一次系數(shù)和常數(shù)，其含義見參數(shù)文件，下同。按照本發(fā)明提供的校正方法，該校正方法還包括步驟1. 3)根據(jù)當(dāng)前環(huán)境的巖性為泥巖段(SP > 80)，按以下四個(gè)層次進(jìn)行校正1. 3. 1)當(dāng)井徑資料大于鉆頭尺寸0. 5英寸時(shí)按以下公式ACCAL = AC-1. 63* (CAL-CAL 鉆)R2 = 11. 93 36440 點(diǎn) B]計(jì)算井徑校正后的聲波時(shí)差值，當(dāng)井徑資料大于鉆頭尺寸不足0. 5英寸時(shí)不進(jìn)行井徑校正。ACCAL為井徑校正后的聲波時(shí)差，AC為原始聲波時(shí)差，CAL為井徑，CAL鉆為鉆頭直徑，下同。1. 3. 2)當(dāng)DEN、CNL曲線存在且正常時(shí)，根據(jù)同類曲線(AC、DEN、CNL均屬孔隙度類)遞歸公式ACDEN = ACDENkl*DEN+ACDENkOR2 = 31. 60,45580 點(diǎn)[7]ACCNL = ACCNLkl*CNL+ACCNLkOR2 = 22. 51,38540 點(diǎn)[8]AC 平均=0. 8*ACDEN+0. 2*ACCNL[9]識(shí)別畸變點(diǎn)，當(dāng)AC平均落在正常值范圍內(nèi)，聲波時(shí)差與AC平均差值大于標(biāo)準(zhǔn)偏差為畸變點(diǎn)，參照公式[7-9]對(duì)聲波時(shí)差進(jìn)行校正。1. 3. 3)當(dāng)DEN、CNL曲線沒測(cè)或不正常時(shí)，根據(jù)井徑穩(wěn)定泥巖段(SP > 80，CAL < 13. 5)聲波時(shí)差與電阻率遞歸公式ACRD = ACRDkl*log(RD)+ACRDkOR2 = 39. 71 53366 點(diǎn)[10]ACRS = ACRSkl*log(RS)+ACRSkOR2 = 52. 47 52832 點(diǎn)[11]AC 平均=0. 8*ACRD+0. 2*ACRS[12]識(shí)別畸變點(diǎn)，當(dāng)AC平均落在正常值范圍內(nèi)，聲波時(shí)差與AC平均差值大于標(biāo)準(zhǔn)偏差為畸變點(diǎn)，參照公式[10-12]對(duì)聲波時(shí)差進(jìn)行校正。1. 3. 4)當(dāng)AC、DEN、CNL、RD、RS均不正常時(shí)，參照井徑穩(wěn)定段聲波時(shí)差與深度的關(guān)系A(chǔ)CTVD = ACTVDk2*TVD2+ACTVDkl^TVD+ACTVDkO R2 = 57. 94 661841 點(diǎn)[13]對(duì)聲波時(shí)差曲線進(jìn)行重建或校正。按照本發(fā)明提供的校正方法，所述校正方法綜合利用了自然電位、中子、密度、深電阻率、淺電阻率、井徑資料，分巖性進(jìn)行多層次、綜合校正。
按照本發(fā)明提供的校正方法，以塔河油田測(cè)井曲線為樣本，將公式遞歸公式 [1-13]糅合成一體識(shí)別畸變點(diǎn)、建立目標(biāo)函數(shù)、對(duì)聲波時(shí)差曲線進(jìn)行校正。按照本發(fā)明提供的校正方法，以塔河油田測(cè)井曲線為樣本，產(chǎn)生密度、中子、深電阻率、淺電阻率資料的統(tǒng)計(jì)分析數(shù)據(jù)，包括最小值、下限值、平均值、中值、標(biāo)準(zhǔn)偏差等，根據(jù)公式[1-4]、W-8]、[10-11]進(jìn)行遞歸值計(jì)算時(shí)，可以根據(jù)自變量的質(zhì)量產(chǎn)生相應(yīng)的浮動(dòng)權(quán)重，自變量質(zhì)量越高對(duì)應(yīng)的權(quán)重越大。按照本發(fā)明提供的校正方法，該校正方法采用基于數(shù)據(jù)庫的參數(shù)文件，上述所有遞歸公式系數(shù)、權(quán)重、最小值、下限值、平均值、中值、最大值、標(biāo)準(zhǔn)偏差等參數(shù)均放置在數(shù)據(jù)庫參數(shù)區(qū)內(nèi)，不同的目的層段調(diào)用不同的參數(shù)，主處理流程、模塊不變，因此可以很方便將該技術(shù)推廣至其它盆地或油田。本發(fā)明提供的測(cè)井聲波時(shí)差信號(hào)校正方法，通過大量統(tǒng)計(jì)、分析建立聲波時(shí)差曲線畸變的識(shí)別標(biāo)準(zhǔn)或標(biāo)志，充分利用聲波時(shí)差資料與自然電位、中子、密度、電阻率、井徑、 鉆時(shí)、埋深資料的遞歸關(guān)系和統(tǒng)計(jì)規(guī)律，編寫自動(dòng)批量處理程序，實(shí)現(xiàn)對(duì)聲波時(shí)差畸變段的修復(fù)校正，提高測(cè)井聲波時(shí)差資料的質(zhì)量。


圖1是本發(fā)明與傳統(tǒng)技術(shù)效果對(duì)比示意圖。
具體實(shí)施例方式下面通過附圖和實(shí)施例，對(duì)本發(fā)明實(shí)施例的技術(shù)方案做進(jìn)一步的詳細(xì)描述。首先，說明本發(fā)明關(guān)鍵建立了一套聲波時(shí)差曲線環(huán)境校正技術(shù)，包括畸變點(diǎn)、段識(shí)別和畸變點(diǎn)、段校正兩個(gè)主要環(huán)節(jié)，畸變點(diǎn)、段的識(shí)別和畸變量的提取建立在大量統(tǒng)計(jì)、分析工作基礎(chǔ)上其中 為使統(tǒng)計(jì)分析工作準(zhǔn)確、真實(shí)，首先要考慮參與統(tǒng)計(jì)工作的資料歸一化問題井身結(jié)構(gòu)的差異導(dǎo)致不能簡(jiǎn)單地根據(jù)井徑值大小判斷井徑擴(kuò)大程度，需要根據(jù)井身結(jié)構(gòu)和鉆頭尺寸對(duì)其進(jìn)行歸一化處理；泥漿類型和地層水礦化度的變化導(dǎo)致自然電位曲線的基準(zhǔn)線不一致，為使得全區(qū)能用統(tǒng)一的自然電位資料指示巖性、計(jì)算泥質(zhì)含量和有效孔隙度，需要對(duì)自然電位曲線進(jìn)行了基準(zhǔn)線校準(zhǔn)處理；為消除測(cè)井儀器類型不同引起的基準(zhǔn)線(刻度器)差異，我們選擇標(biāo)志層為樣本進(jìn)行了儀器因素的歸一化處理。第二步，說明本發(fā)明的核心基于歸一化處理后的大量統(tǒng)計(jì)分析資料表明自然電位資料能夠很好的識(shí)別巖性，砂巖段井徑正常，聲波、密度曲線畸變較輕；泥巖段井徑擴(kuò)大，聲波時(shí)差曲線畸變嚴(yán)重。 因此，校正思路按巖性分別進(jìn)行。砂巖段聲波時(shí)差與中子、密度、電阻率存在很好的遞歸關(guān)系[1-4]，偏離遞歸平均值上限的個(gè)別離散點(diǎn)(段)認(rèn)為其發(fā)生了畸變，根據(jù)公式[1-5]加以校正。泥巖段聲波時(shí)差畸變嚴(yán)重，校正工作分四個(gè)層次首先，分析表明聲波時(shí)差與井徑有較明顯的遞歸關(guān)系聲波時(shí)差隨井徑的擴(kuò)大而往大畸變。遞歸關(guān)系的一次系數(shù)被用來提取井徑擴(kuò)大引起的系統(tǒng)畸變量，根據(jù)公式[6]將畸變量從原始聲波時(shí)差曲線中減除完成井徑校正；其次，當(dāng)DEN、CNL曲線存在且正常時(shí)，根據(jù)聲波時(shí)差與中子、密度同類曲線間的遞歸關(guān)系[7-8]識(shí)別畸變點(diǎn)，當(dāng)AC平均落在正常值范圍內(nèi)，聲波時(shí)差與AC平均差值大于標(biāo)準(zhǔn)偏差為畸變點(diǎn)，根據(jù)公式[9]對(duì)畸變點(diǎn)進(jìn)行校正。第三，考慮到感應(yīng)電阻率曲線橫向探測(cè)半徑相對(duì)較大，受井眼環(huán)境影響相對(duì)較小，在中子、密度資料缺失或不正常的情況下，根據(jù)井徑穩(wěn)定泥巖段聲波時(shí)差與電阻率間的遞歸關(guān)系[10-11]識(shí)別畸變點(diǎn)，當(dāng)AC平均落在正常值范圍內(nèi)，聲波時(shí)差與AC平均差值大于標(biāo)準(zhǔn)偏差為畸變點(diǎn)，依據(jù)公式[12]對(duì)畸變點(diǎn)進(jìn)行校正。第四，當(dāng)AC、DEN、CNL、RD、RS均不正常時(shí)，根據(jù)井徑穩(wěn)定段聲波時(shí)差與深度的關(guān)系[13] 對(duì)聲波時(shí)差曲線進(jìn)行重建或校正。第三步，說明本發(fā)明的具體實(shí)現(xiàn)為使校正工作規(guī)范化，為使本技術(shù)推廣工作簡(jiǎn)單，我們將上述校正工作的核心思想編寫為基于Oracle數(shù)據(jù)庫的實(shí)用程序或主程序，而將所有統(tǒng)計(jì)分析得到的遞歸公式、平均值、標(biāo)準(zhǔn)偏差和流程指示參數(shù)放在一個(gè)單獨(dú)的參數(shù)文件內(nèi)。主程序分為兩部分，前半部分對(duì)砂巖進(jìn)行校正，后半部分對(duì)泥巖進(jìn)行校正，校正結(jié)果在數(shù)據(jù)庫內(nèi)新產(chǎn)生DT曲線，源聲波時(shí)差保留不變，主程序結(jié)構(gòu)如下(一為注釋標(biāo)志)—定義變量VARIABLE ACRD；ACRS ；ACDEN ；ACCNL ；ACCAL ；ACVSH ；VARIABLE krd；krs ；kden ；kcnl ；kcal ；acmix ；acmid ；acmin；index；i ；-注釋部分—if index = 1 for RD RS DEN CNL CAL ；—if index = 2 for RD RS DEN CNL ；—if index = 3 for RD RS DEN ；—if index = 4 for RD RS ；；校正選項(xiàng)index = 1 ；IF(〃 SP" < 90. 0)THEN —For Sand Lithology (選擇砂巖段)；i = 1 ；krd = 0 ；krs = 0 ；kden = 0 ；kcnl = 0 ；kcal = 0 ；—初女臺(tái)權(quán)系數(shù)；while ((i < 6) and ((krd+krs+kden+kcnl+kcal) < 0. 05)) loop—砂巖遞歸分析循環(huán)；if ((SACRDk ！ = 0) and ( “ RD_M “ > = (srdmin/i)) and ( “ RD_M “ < =
(srdmax+(i_l)^srddev))) then...... ；else...... ；end if ；―聲波時(shí)差 深電阻率遞歸計(jì)算；if ((SACRSkl ！ = 0) and ( “ RS_M “ > = (srsmin/i)) and ( “ RS_M “ < = (srsmax+ (i_l)^srsdev))) then...... ；else...... ；end if ；―聲波時(shí)差 淺電阻率遞歸計(jì)算；if(〃 DEN" ！ = null) then—判斷密度曲線存在情況；if ((SACDENkl ！ = 0) and ( “ DEN “ > = (sdenmin+(i_3)/4*sdendev)) and(" DEN"<= (sdenmax+(i_l)*sdendev)))then...... ；else...... ；end if;—密度存在
時(shí)聲波時(shí)差 密度遞歸處理；
else...... ；end if ；―密度不存在情況；if(" CNL" ！ = null) then —判斷中子曲線存在情況；if ((SACCNLkl ！ = 0) and ( “ CNL “ > = (scnlmin+(i_3)/4*scnldev)) and(" CNL“ <= (scnlmax+(i-l) ^scnldev))) then...... ；else...... ； end if;—中子存在聲
波時(shí)差 中子遞歸處理；else...... ；end if ；—中子不存在情況；i = i+1 ；end loop ；—砂巖遞歸分析循環(huán)結(jié)束；if ((krd+krs+kden+kcnl) > 0) then —計(jì)算權(quán)系數(shù)；if ((index > 0. 5) and (index < 2. 5)) then......；—校正選項(xiàng)為 1 或 2 時(shí)遞歸平均值；else if ((index > 2. 5) and (index < 3. 5)) then......；—校正選項(xiàng)為 3 時(shí)遞歸平均值；else if ((krd+krs) > 0) then......； —校正選項(xiàng)為 4 時(shí)遞歸平均值；else......end if ；—聲波時(shí)差 深度遞歸平均值；end if ；end if ；else acmix = 〃 AC" ；end if ；―遞歸平均值計(jì)算結(jié)束；if ((" AC" >= (acmix+sacdev))and(〃 AC" <= (acmix+2*sacdev))) then............-聲波時(shí)差上偏遞歸值一倍標(biāo)準(zhǔn)偏差時(shí)時(shí)；else if(〃 AC" > (acmix+2*sacdev)) then......—聲波時(shí)差上偏遞歸值兩倍標(biāo)準(zhǔn)偏差時(shí)；else if(〃 AC" >= (acmix-sacdev)) then......—聲波時(shí)差下偏遞歸值一倍標(biāo)準(zhǔn)偏差時(shí)；else if (“ AC" < sacminavg) then............—聲波時(shí)差低于聲波時(shí)差下限；else......end if ；end if ；end if ；—聲波時(shí)差正常時(shí)；end if ；—砂巖段校正結(jié)束；ELSE—開始泥巖段校正；i = 1 ；krd = 0 ；krs = 0 ；kden = 0 ；kcnl = 0 ；kcal = 0 ；—定義初始權(quán)系數(shù)；while ((i < 6) and ((krd+krs+kden+kcnl+kcal) < 0. 05)) loop—泥巖遞歸分析循環(huán)；if ((ACRDkl ！ = 0) and ( “ RD_M “ > = (rdmin/i)) and ( “ RD_M “ < = (rdmax+(i-1)*rddev
))) then......else......end if ；—聲波時(shí)差 深電阻率遞歸計(jì)算；if ((ACRSkl ！ = 0) and ( “ RS_M “ > = (rsmin/i)) and ( “ RS_M “ < = (rsmax+ (i_l)*rsdev)))then......else......end if ； —聲波時(shí)差 淺電阻率遞歸計(jì)算；if(〃 DEN" ！ = null) then—判斷密度曲線是否存在；if((ACDENkl ！ =0)and(〃 DEN" >= (denmin+(i_3)/4*dendev))and(〃 DEN" < =(denmax+(i_l)*dendev)))then......Else......end if ；―密度存在時(shí)聲波時(shí)
差 密度遞歸處理；Else......end if； —密度不存在情況；if(" CNL" ！ = null) then —判斷中子曲線是否存在；if((ACCNLkl ！ =0)and(" CNL" >= (cnlmin+(i_3)/4*cnldev))and(〃 CNL“ <=(cnlmax+(i-l)*cnldev)))then......Else......end if ；―中子存在時(shí)聲波時(shí)
差 中子遞歸處理；Else......end if；—中子不存在情況；if(〃 CAL_BIT〃 ！ = null) then—判斷井徑曲線是否存在；if ((ACCALkl ！ = 0) and (“ CAL_BIT〃 >= (calmin/i)) and (〃 CAL_BIT〃 < = (calmax+(i_l)*caldev)))then......else......end if ；―井徑存在聲波時(shí)差 井徑遞歸處理；else......end if ； —井徑不存在或不正常情況；i = i+1 ；end loop ； —泥巖遞歸分析循環(huán)結(jié)束；if ((krd+krs+kden+kcnl) > 0) then—計(jì)算砂巖遞歸平均值；if ((index > 0. 5) and (index < 1. 5)) then......—校正選項(xiàng)為 1 時(shí)遞歸平均值；else if ((index > 1. 5) and (index < 2. 5)) then —校正選項(xiàng)為 2時(shí)遞歸平均值；else if ((index > 2· 5) and (index < 3· 5)) then—校正選項(xiàng)為 3 時(shí)遞歸平均值；else if((krd+krs) > 0) then —校正選項(xiàng)為 4 時(shí)遞歸平均值；else. . . end if ；—聲波時(shí)差 深度遞歸平均值；
end if ；end if ；end if ；else...—遞歸分析無效情況；end if ；—遞歸平均值計(jì)算結(jié)束；if ((" AC" >= (acmix-acdev))and(〃 AC" <= (acmix+acdev))) then—聲波時(shí)差正常時(shí)；else if((" AC" > (acmix+acdev))and(" AC" <= (acmix+2*acdev)))then—聲波時(shí)差上偏遞歸值1_2倍標(biāo)準(zhǔn)偏差時(shí)；else if(〃 AC" >= (acmix+2*acdev)) then—聲波時(shí)差上偏遞歸值兩倍標(biāo)準(zhǔn)偏差以上時(shí)；else if(" AC" < acminavg) then —聲波時(shí)差低于聲波時(shí)差下限時(shí)；else—聲波時(shí)差下偏遞歸值一倍標(biāo)準(zhǔn)偏差 下限；end if ；end if ；end if ；end if;—泥巖段校正結(jié)束；end if； ‘‘ DT" = acmid ；―校正后聲波時(shí)差曲線命名為“DT，，；參數(shù)文件具有缺省值一某井適用值一某井某地層單元適用值參數(shù)遺傳功能?？梢圆捎靡粋€(gè)復(fù)雜參數(shù)文件的方式(將所有地層單元參數(shù)按地層框架表父子、兄弟關(guān)系放置就位)，也可以采用多個(gè)簡(jiǎn)單參數(shù)文件的方式(一個(gè)參數(shù)文件放置一個(gè)地層單元參數(shù))；下面以塔里木盆地石炭系卡拉沙依組地層單元為例，介紹參數(shù)文件內(nèi)的變量和含義(一為注釋標(biāo)志)Prd =1.0；—深電阻率權(quán)重；Prs = 0. 5 ；—淺電阻率權(quán)重；Pden = 0. 8 ；—密度權(quán)重；Pcnl = O. 5;—中子權(quán)重；Pcal = 0. 8 ；—井徑權(quán)重ACTVDkO = 128. 37 ；—聲波時(shí)差 垂深遞歸系數(shù)0 ；ACTVDkl = -0. 0157 ；—聲波時(shí)差 垂深遞歸系數(shù)1 ；ACTVDk2 = 7*10-7 ；—聲波時(shí)差 垂深遞歸系數(shù)2 ；—For Clkl2 shale —以下為泥巖段參數(shù)；ACRDkl = -8. 30 ；—聲波時(shí)差 深電阻率遞歸系數(shù)1 ；ACRDkO = 70. 27 ；—聲波時(shí)差 深電阻率遞歸系數(shù)0 ；ACRSkl = -8. 82 ；—聲波時(shí)差 淺電阻率遞歸系數(shù)1 ；ACRSkO = 71. 16 ；—聲波時(shí)差 淺電阻率遞歸系數(shù)0 ；-ACCALkl = 1. 63 ；—聲波時(shí)差 井徑遞歸系數(shù)1 ；-ACCALkO = 47. 78 ；—聲波時(shí)差 井徑遞歸系數(shù)0 ；ACDENkl = -5. 17 ；—聲波時(shí)差 密度遞歸系數(shù)1 ；ACDENkO = 76. 33 ；—聲波時(shí)差 密度遞歸系數(shù)0 ；ACCNLkl = 0. 31 ；—聲波時(shí)差 中子遞歸系數(shù)1 ；ACCNLkO = 58. 24 ；—聲波時(shí)差 中子遞歸系數(shù)0 ；
rdmin = 1. 64 ；—深電阻率最小值；rdmid = 4. 22 ；—深電阻率中值；rdmax = 9. 25 ；—深電阻率最大值；rddev = 0. 28 ；—深電阻率標(biāo)準(zhǔn)偏差；rsmin = 0. 80 ；—淺電阻率最小值；rsmid = 3. 69 ；—淺電阻率中值；rsmax = 9. 20 ；—淺電阻率最大值；rsdev = 0. 51 ；—淺電阻率標(biāo)準(zhǔn)偏差；denmin = 2. 01 ；—密度最小值；denmid = 2. 43 ；—密度中值；denmax = 2. 72 ；—密度最大值；dendev = 0. 05 ；—密度標(biāo)準(zhǔn)偏差cnlmin = 7. 50 ；__ 中子最々、值；cnlmid = 25. 06 ；—中子中值；cnlmax = 48. 80 ；—中子最大值；cnldev = 25. 31 ；—中子標(biāo)準(zhǔn)偏差；calmin = 11. 97 ；—井徑最小值；calmid = 14. 89 ；—井徑中值；calmax = 16. 94 ；—井徑最大值；acavg = 71. 63 ；—聲波時(shí)差平均值；acdev = 4. 74 ；—聲波時(shí)差標(biāo)準(zhǔn)偏差；acminavg = 61. 95 ；—聲波時(shí)差下限；acmindev = 1. 15 ；—聲波時(shí)差下限標(biāo)準(zhǔn)偏差；—For Clkl2 Sand—以下為砂巖段參數(shù)；SACRDkl = -8. 09 ；—聲波時(shí)差 深電阻率遞歸系數(shù)1 ；SACRDkO = 72. 01 ；—聲波時(shí)差 深電阻率遞歸系數(shù)0 ；SACRSkl = -9. 48 ；—聲波時(shí)差 淺電阻率遞歸系數(shù)1 ；SACRSkO = 75. 73 ；—聲波時(shí)差 淺電阻率遞歸系數(shù)0 ；SACDENkl = -31. 84；—聲波時(shí)差 密度遞歸系數(shù)1 ；SACDENkO = 147. 81；—聲波時(shí)差 密度遞歸系數(shù)0 ；SACCNLkl = 0. 71 ；—聲波時(shí)差 中子遞歸系數(shù)1 ；SACCNLkO = 59. 43 ；—聲波時(shí)差 中子遞歸系數(shù)0
srdmin = 0. 27 ；—深電阻率最小值；srdmid = 2. 80 ；—深電阻率中值；srdmax = 9. 11 ；—深電阻率最大值；srddev = 0. 71 ；—深電阻率標(biāo)準(zhǔn)偏差；srsmin = 0. 85 ；—淺電阻率最小值；srsmid = 4. 85 ；—淺電阻率中值；srsmax = 9. 53 ；—淺電阻率最大值；
srsdev = 1. 49 ；—淺電阻率標(biāo)準(zhǔn)偏差；sdenmin = 2. 05 ； —密度最小值；sdenmid = 2. 46 ； —密度中值；sdenmax = 2. 68 ； —密度最大值；sdendev = 0. 05 ； —密度標(biāo)準(zhǔn)偏差；scnlmin = 4. 67 ； —中子最小值；scnlmid = 14. 01 ； —中子中值；scnlmax = 31. 53 ； —中子最大值；scnldev = 4. 33 ； —中子標(biāo)準(zhǔn)偏差；sacavg = 69. 75 ； —聲波時(shí)差平均值；sacdev = 5. 18 ；—聲波時(shí)差標(biāo)準(zhǔn)偏差；sacminavg = 62. 19 ；—聲波時(shí)差下限平均值；sacmindev = 1. 12 ；—聲波時(shí)差下限標(biāo)準(zhǔn)偏差；最后，說明本發(fā)明效果校正工作按采樣間隔逐井、逐點(diǎn)進(jìn)行，保留了測(cè)井資料的垂向高分辨率，校正效果如圖1所示，校正前聲波時(shí)差(虛線)高頻干擾成分多，校正后(實(shí)線)聲波時(shí)差高頻干擾被修正、曲線穩(wěn)定可靠。直觀校正效果圖1可以看出砂巖段校正量小、泥巖段校正量大，聲波時(shí)差欄右側(cè)小值部分基本沒變，聲波時(shí)差欄左側(cè)由周波跳躍、時(shí)差展寬引起的大值畸變部分被往右校正歸位，校正后聲波時(shí)差的動(dòng)態(tài)范圍減小，說明校正尊重了聲波時(shí)差的畸變機(jī)理。量化校正效果衡量聲波時(shí)差曲線校正效果的最好手段是合成地震記錄，最可信、 可靠的參數(shù)是合成道與井旁道的互相關(guān)系數(shù)和測(cè)井約束的從井旁地震道中提取實(shí)際地震子波的視信噪比。根據(jù)校正后聲波時(shí)差曲線從井旁地震道中提取的實(shí)際地震子波的信噪比、地震道與合成地震道之間的互相關(guān)系數(shù)明顯較依據(jù)校正前曲線的高，以塔河油田艾協(xié)克三維工區(qū)內(nèi)的96 口井為例校正前提取的子波信噪比=0.81，合成道與井旁道的相關(guān)系數(shù)=0. 58 ；校正后提取的子波信噪比=1. 82，合成道與井旁道的相關(guān)系數(shù)=0. 78。96 口井的統(tǒng)計(jì)資料表明砂巖地層校正前后聲波時(shí)差平均值沒有變化，標(biāo)準(zhǔn)偏差由2. Ius/ ft — 1.6us/ft ；泥巖地層聲波時(shí)差平均值由70. 9us/ft —69. lus/ft，校正量約為2. 64%, 標(biāo)準(zhǔn)偏差由2. 2us/ft — 1. ;3us/ft。砂、泥巖校正量不同，說明校正尊重了聲波時(shí)差的畸變機(jī)理；校正后聲波時(shí)差的標(biāo)準(zhǔn)偏差變小，從測(cè)量學(xué)的角度看意味著校正后的測(cè)量精度提高。最后所應(yīng)說明的是，以上實(shí)施例僅用以說明本發(fā)明實(shí)施例的技術(shù)方案而非限制， 盡管參照較佳實(shí)施例對(duì)本發(fā)明實(shí)施例進(jìn)行了詳細(xì)說明，本領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解， 可以對(duì)本發(fā)明實(shí)施例的技術(shù)方案進(jìn)行修改或者等同替換，而不脫離本發(fā)明實(shí)施例技術(shù)方案和權(quán)利要求的精神和范圍。
權(quán)利要求
1.一種測(cè)井聲波時(shí)差信號(hào)校正方法，其特征在于，包括以下步驟 獲取測(cè)井聲波時(shí)差信號(hào)并按實(shí)際情況設(shè)置當(dāng)前環(huán)境；根據(jù)當(dāng)前環(huán)境的巖性為砂巖(SP < 90)，自動(dòng)根據(jù)遞歸公式 ACRD = ACRD*kl*log(RD)+ACRD*kO R2 = 29.57,24410 點(diǎn) ACRS = ACRS*kl*log(RS)+ACRS*kO R2 = 28.18，23867 點(diǎn) ACDEN = ACDEN*kl*DEN+ACDEN*kO R2 = 29.57，16301 點(diǎn) ACCNL = ACCNL*kl*CNL+ACCNL*kO R2 = 20. 01 18583 點(diǎn) AC 平均=Prd*ACRD+Prs*ACRS+Pden*ACDEN+Pcnl*ACCNL計(jì)算AC平均，與實(shí)測(cè)聲波時(shí)差比較差值大于標(biāo)準(zhǔn)偏差的為畸變點(diǎn)，再對(duì)畸變點(diǎn)按AC平均對(duì)其校正；式中SP為自然電位，RD為深電阻率，RS為淺電阻率，DEN為密度，CNL為中子，ACRD、ACRS, ACDEN、ACCNL 分別為據(jù) RD、RS、DEN、CNL 遞歸的聲波時(shí)差，Prd, Prs, Pden, Pcnl分別為RD、RS、DEN、CNL的權(quán)重，R2為遞歸方程的相關(guān)性指標(biāo)，AC平均為遞歸平均值， kl、k0是遞歸方程的一次系數(shù)和常數(shù)。
2.根據(jù)權(quán)利要求1所述校正方法，其特征在于，該校正方法還包括步驟1.3)根據(jù)當(dāng)前環(huán)境的巖性為泥巖段(SP > 80)，自動(dòng)按以下步驟進(jìn)行校正當(dāng)井徑資料大于鉆頭尺寸0. 5英寸時(shí)按公式 ACCAL = AC-1. 63*(CAL-CAL 鉆)R2 = 11. 93 36440 點(diǎn)計(jì)算井徑校正后的聲波時(shí)差值，當(dāng)井徑資料大于鉆頭尺寸不足0. 5英寸時(shí)不進(jìn)行井徑校正，式中ACCAL為井徑校正后的聲波時(shí)差，AC為原始聲波時(shí)差，CAL為井徑，CAL鉆為鉆頭直徑；當(dāng)DEN、CNL曲線存在且正常時(shí)，根據(jù)同類曲線遞歸公式 ACDEN = ACDEN*kl*DEN+ACDEN*kO R2 = 31. 60，45580 點(diǎn) ACCNL = ACCNL*kl*CNL+ACCNL*kO R2 = 22.51，38540 點(diǎn) AC 平均=0. 8*ACDEN+0. 2*ACCNL識(shí)別畸變點(diǎn)，當(dāng)AC平均落在正常值范圍內(nèi)，聲波時(shí)差與AC平均差值大于標(biāo)準(zhǔn)偏差為畸變點(diǎn)，并按該同類曲線遞歸公式對(duì)聲波時(shí)差進(jìn)行校正；當(dāng)DEN、CNL曲線沒測(cè)或不正常時(shí)，根據(jù)井徑穩(wěn)定泥巖段(SP > 80，CAL < 13. 5)聲波時(shí)差與電阻率遞歸公式ACRD = ACRD*kl*log(RD)+ACRD*kO R2 = 39. 71 53366 點(diǎn) ACRS = ACRS*kl*log(RS)+ACRS*kO R2 = 52. 47 52832 點(diǎn) AC 平均=0. 8*ACRD+0. 2*ACRS識(shí)別畸變點(diǎn)，當(dāng)AC平均落在正常值范圍內(nèi)，聲波時(shí)差與AC平均差值大于標(biāo)準(zhǔn)偏差為畸變點(diǎn)，并按該聲波時(shí)差與電阻率公式對(duì)聲波時(shí)差進(jìn)行校正；當(dāng)AC、DEN、CNL、RD、RS均不正常時(shí)，按井徑穩(wěn)定段聲波時(shí)差與深度的關(guān)系 ACTVD = ACTVD*k2*TVD2+ACTVD*kl*TVD+ACTVD*kO R2 = 57. 94 661841 點(diǎn)對(duì)聲波時(shí)差曲線進(jìn)行重建或校正。
3.根據(jù)權(quán)利要求1或2所述校正方法，其特征在于，所述校正方法綜合利用自然電位、 中子、密度、深電阻率、淺電阻率和井徑資料分巖性進(jìn)行多層次、綜合校正。
4.根據(jù)權(quán)利要求1或2所述校正方法，其特征在于，將所有所述遞歸公式糅合成一體識(shí)別畸變點(diǎn)、建立目標(biāo)函數(shù)、對(duì)聲波時(shí)差曲線進(jìn)行校正。
5.根據(jù)權(quán)利要求4所述校正方法，其特征在于，該校正方法應(yīng)用于塔河油田以塔河油田測(cè)井曲線為樣本，將遞歸公式糅合成一體識(shí)別畸變點(diǎn)、建立目標(biāo)函數(shù)、對(duì)聲波時(shí)差曲線進(jìn)行校正。
6.根據(jù)權(quán)利要求1或2所述校正方法，其特征在于，根據(jù)所述遞歸公式進(jìn)行遞歸值計(jì)算時(shí)，可以根據(jù)自變量RD、RS、DEN、CNL的質(zhì)量產(chǎn)生相應(yīng)的浮動(dòng)權(quán)重，自變量質(zhì)量越高對(duì)應(yīng)的權(quán)重越大。
7.根據(jù)權(quán)利要求6所述校正方法，其特征在于，該校正方法應(yīng)用于塔河油田以塔河油田測(cè)井曲線為樣本，產(chǎn)生RD、RS、DEN、CNL資料的統(tǒng)計(jì)分析數(shù)據(jù)，包括最小值、下限值、平均值、中值、標(biāo)準(zhǔn)偏差。
8.根據(jù)權(quán)利要求1或2所述校正方法，其特征在于，該校正方法是采用基于數(shù)據(jù)庫的參數(shù)文件的處理流程。
9.根據(jù)權(quán)利要求8所述校正方法，其特征在于，所述數(shù)據(jù)庫參數(shù)文件包括所有所述遞歸公式的系數(shù)參數(shù)、權(quán)重參數(shù)、最小值參數(shù)、下限值參數(shù)、平均值參數(shù)、中值參數(shù)、最大值參數(shù)、標(biāo)準(zhǔn)偏差參數(shù)。
10.根據(jù)權(quán)利要求9所述校正方法，其特征在于，該校正方法應(yīng)用于塔河油田不同的目的層段的測(cè)井聲波對(duì)應(yīng)的所述參數(shù)不同，所述處理流程及其包含的模塊相同。
全文摘要
本發(fā)明涉及一種測(cè)井聲波時(shí)差信號(hào)校正方法，獲取測(cè)井聲波時(shí)差信號(hào)并按實(shí)際情況設(shè)置當(dāng)前環(huán)境；分巖性進(jìn)行測(cè)井聲波時(shí)差信號(hào)校正。這種校正方法通過大量統(tǒng)計(jì)、分析建立聲波時(shí)差曲線畸變的識(shí)別標(biāo)準(zhǔn)或標(biāo)志，充分利用聲波時(shí)差資料與自然電位、中子、密度、電阻率、井徑、鉆時(shí)、埋深資料的遞歸關(guān)系和統(tǒng)計(jì)規(guī)律，編寫自動(dòng)批量處理程序，實(shí)現(xiàn)對(duì)聲波時(shí)差畸變段的修復(fù)校正，提高測(cè)井聲波時(shí)差資料的質(zhì)量。
文檔編號(hào)E21B49/00GK102454399SQ20101051971
公開日2012年5月16日 申請(qǐng)日期2010年10月26日 優(yōu)先權(quán)日2010年10月26日
發(fā)明者劉群, 安精文, 李家蓉, 李愛紅, 蔣進(jìn)勇 申請(qǐng)人:中國(guó)石油化工股份有限公司