專利名稱：一種優(yōu)選風(fēng)化層速度的確定方法
技術(shù)領(lǐng)域：
本發(fā)明涉及石油勘探領(lǐng)域，應(yīng)用于地震數(shù)據(jù)處理過(guò)程中，特別是涉及一種優(yōu)選風(fēng)化層速度的確定方法。
背景技術(shù)：
現(xiàn)有技術(shù)的地震勘探解釋的理論中，一般都假定激發(fā)點(diǎn)與接收點(diǎn)是在一個(gè)水平面上，并且地層速度是均勻的。但實(shí)際上地面常常不平坦，各個(gè)激發(fā)點(diǎn)深度也可能不同，低速帶中的波速與地層中的波速又相差懸殊，所以必將影響實(shí)測(cè)的時(shí)距曲線形狀。為了消除這些影響，對(duì)原始地震數(shù)據(jù)要進(jìn)行地形校正、激發(fā)深度校正、低速帶校正等。由于這些校正對(duì)同一觀測(cè)點(diǎn)的不同地震界面都是不變的，因此統(tǒng)稱靜校正。而在靜校正過(guò)程中，特別是長(zhǎng)波長(zhǎng)靜校正問(wèn)題是復(fù)雜近地表地區(qū)地震資料處理中提高資料處理質(zhì)量、保證獲得地下真實(shí)構(gòu)造形態(tài)需要解決的主要問(wèn)題。目前在地震資料處理領(lǐng)域中有很多關(guān)于折射波靜校正的方法，如GLI，EGRM,加減法，TOMO等等?，F(xiàn)有技術(shù)中折射波靜校正的方法實(shí)現(xiàn)過(guò)程中涉及風(fēng)化層速度填寫方法通常為以下兩種方式第一種方式是當(dāng)研究工區(qū)內(nèi)能預(yù)先提供局部控制點(diǎn)風(fēng)化層速度信息時(shí)，把這些信息分別填入程序各控制點(diǎn)位置對(duì)應(yīng)的風(fēng)化層速度參數(shù)欄。反演程序執(zhí)行時(shí)，其他點(diǎn)的速度參數(shù)通過(guò)控制點(diǎn)信息插值獲得；第二種方式是當(dāng)研究工區(qū)內(nèi)不能預(yù)先提供任何風(fēng)化層速度信息時(shí)。在程序的風(fēng)化層速度參數(shù)欄填入任意一個(gè)常數(shù)，用于全區(qū)反演運(yùn)算，從而由于任意常數(shù)選擇的盲目性，很容易造成選擇的速度和實(shí)際速度有較大偏差，由此引起較大的模型反演錯(cuò)誤和靜校正錯(cuò)誤，超過(guò)了后續(xù)剩余靜校正程序所能修正的限制，不能獲得理想的靜校正結(jié)果。綜上所述，發(fā)明人在實(shí)現(xiàn)本發(fā)明過(guò)程中發(fā)現(xiàn)現(xiàn)有技術(shù)中存在以下缺陷預(yù)獲得可靠的折射波靜校正結(jié)果，必須預(yù)先掌握近地表地層風(fēng)化層速度的橫向變化信息，因此需要預(yù)先進(jìn)行諸如微測(cè)井等的近地表勘測(cè)，從而獲取近地表真實(shí)數(shù)據(jù)，在以上的勘測(cè)過(guò)程中需要巨大的成本投入，同時(shí)由于其勘測(cè)點(diǎn)的分布及數(shù)據(jù)限制，其對(duì)不同勘測(cè)點(diǎn)的選擇與數(shù)據(jù)處理，將直接影響到后期數(shù)據(jù)的精確度。同時(shí)，盲目的任意常數(shù)選擇容易造成選擇的速度和實(shí)際速度有較大偏差，最終無(wú)法達(dá)到預(yù)期效果。

發(fā)明內(nèi)容
本發(fā)明實(shí)施例要解決的問(wèn)題是提供一種優(yōu)選風(fēng)化層速度的確定方法，解決了因復(fù)雜地區(qū)近地表信息無(wú)法預(yù)先精確確定，由此無(wú)法準(zhǔn)獲取風(fēng)化層速度，從而造成后續(xù)靜校正精度偏低的問(wèn)題。為達(dá)到上述目的，本發(fā)明實(shí)施例提出一種優(yōu)選風(fēng)化層速度的確定方法，具體包括預(yù)設(shè)代表風(fēng)化層速度V0 ；根據(jù)所述預(yù)設(shè)代表風(fēng)化層速度Vtl獲取實(shí)際風(fēng)化層速度變化區(qū)間；判斷所述實(shí)際風(fēng)化層速度變化區(qū)間是否介于臨界風(fēng)化層速度區(qū)間中，若是，則將所述預(yù)設(shè)代表風(fēng)化層速度Vtl為處理中所用的代表風(fēng)化層速度，若否，則重新預(yù)設(shè)代表風(fēng)化層速度Vf!。根據(jù)本發(fā)明的實(shí)施例以上技術(shù)方案獲得了以下技術(shù)效果在折射波靜校正過(guò)程中，不再需要預(yù)先進(jìn)行諸如微測(cè)井等的近地表調(diào)查，就可以使復(fù)雜地區(qū)為了獲得地下勘探目的層精確構(gòu)造成像而必需進(jìn)行的短、長(zhǎng)波長(zhǎng)折射波靜校正在精度方面得到提高。


為了更清楚地說(shuō)明本發(fā)明實(shí)施例的技術(shù)方案，下面將對(duì)實(shí)施例描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本發(fā)明的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù)這些附圖獲得其他的附圖。圖I為本發(fā)明實(shí)施例中優(yōu)選風(fēng)化層速度的確定方法的流程示意圖；圖2為本發(fā)明的近地表模型示意圖；圖3為本發(fā)明的采用代表風(fēng)化層速度Vtl為800米/秒、不同的折射層速度和替換速度組合得到的誤差曲線。圖4為本發(fā)明實(shí)施例中實(shí)際得到的近地表模型。圖5為本發(fā)明實(shí)施例中的近地表模型和實(shí)際炮記錄的吻合度檢查圖。圖6為本發(fā)明實(shí)施例中原來(lái)處理的剖面和新處理過(guò)的剖面對(duì)比圖。
具體實(shí)施例方式本發(fā)明實(shí)施例提供一種優(yōu)選風(fēng)化層速度的確定方法，為滿足折射波靜校正過(guò)程中校正精度，通過(guò)對(duì)剩余靜校正值8 t數(shù)據(jù)及截距時(shí)的設(shè)定，通過(guò)誤差曲線公式從而確定其風(fēng)化層速度Vtl值。下面結(jié)合附圖和實(shí)施例，對(duì)本發(fā)明的具體實(shí)施方式
作進(jìn)一步詳細(xì)描述實(shí)施例如圖I所示，為本發(fā)明一種優(yōu)選風(fēng)化層速度的確定方法的流程示意圖，包括以下步驟步驟SlOl :預(yù)設(shè)代表風(fēng)化層速度Vci ;所述預(yù)設(shè)代表風(fēng)化層速度Vci可選取750米/秒 850米/秒之間的任意一個(gè)值,例如800米/秒。步驟S102 :根據(jù)所述預(yù)設(shè)代表風(fēng)化層速度Vtl獲取實(shí)際風(fēng)化層速度變化區(qū)間，具體過(guò)程如下；I)、根據(jù)剩余靜校正值S t及截距時(shí)t的邊界范圍值，確定所述剩余靜校正值5 t、截距時(shí)t，在實(shí)際的地質(zhì)環(huán)境下，為使剩余靜校正量在后續(xù)的剩余靜校正處理中可以被消除掉，通常剩余靜校正值S t具體邊界范圍絕對(duì)值不超過(guò)20毫秒；同時(shí)對(duì)截距時(shí)t的取值范圍，根據(jù)我國(guó)地質(zhì)地貌的實(shí)際情況，可具體選取為截距時(shí)t的邊界范圍值不超過(guò)400毫秒，即風(fēng)化層厚度不超過(guò)300米，因此，所述剩余靜校正值S t的絕對(duì)值為20毫秒，所述截距時(shí)t的取值為400毫秒。2)、將所述剩余靜校正值5 t及截距時(shí)t代入如下誤差曲線公式中獲取誤差指數(shù)、= f-當(dāng)所述剩余靜校正值S t為±20毫秒，及截距時(shí)t為400毫秒時(shí)，所述誤差指數(shù)k等于±0. I。3)、將所述誤差指數(shù)k及分別將所述預(yù)設(shè)代表風(fēng)化層速度Vtl代入如下誤差曲線公式中分別獲取剩余風(fēng)化層速度dV(l ;如將2)中獲取的k等于±0. 1，及在3)中通過(guò)預(yù)設(shè)選取的代表風(fēng)化層速度預(yù)設(shè)值％為800米/秒，代入如下誤差曲線公式中，可獲得對(duì)應(yīng)于預(yù)設(shè)代表風(fēng)化層速度Vtl (800米/秒)的剩余風(fēng)化層速度dV(l為±300米/秒。在以下公式中所示的V1為折射層速度、vr為替換速度，均可通過(guò)設(shè)定獲?。籎c ^ (-1^- I ￥-^----J—)+ dV 0 .
\v r/VkCOS InCOS I /Vr CDS Ini = arcsin - / = arcsin '"~'A 0其中，V1為折射層速度、vr為替換速度；4)、根據(jù)所述預(yù)設(shè)代表風(fēng)化層速度Vtl及所述剩余風(fēng)化層速度dV(l對(duì)應(yīng)獲取實(shí)際風(fēng)化層速度變化區(qū)間。如通過(guò)I)中當(dāng)所述預(yù)設(shè)代表風(fēng)化層速度Vtl為800米/秒，所述剩余風(fēng)化層速度dV(l為±300米/秒時(shí)，實(shí)際風(fēng)化層速度變化區(qū)間為500米/秒 1100米/秒。根據(jù)以上步驟S102，首先，可通過(guò)設(shè)定不同的剩余靜校正值St及截距時(shí)t的邊界范圍值，對(duì)誤差指數(shù)k進(jìn)行設(shè)定，之后，在確定了誤差指數(shù)k的情況下，結(jié)合預(yù)設(shè)的代表風(fēng)化層速度值％，獲得與其對(duì)應(yīng)的不同實(shí)際風(fēng)化層速度變化區(qū)間，如預(yù)設(shè)代表風(fēng)化層速度值V0為800米/秒、750米/秒、850米/秒，或750米/秒 850米/秒中間的任意一個(gè)值時(shí)，將獲得不同的實(shí)際風(fēng)化層速度變化區(qū)間，而獲得的每個(gè)實(shí)際風(fēng)化層速度變化區(qū)間與其預(yù)設(shè)的代表風(fēng)化層速度值％ 一一對(duì)應(yīng)，如當(dāng)實(shí)際風(fēng)化層速度變化區(qū)間為500米/秒 1100米/秒時(shí)，對(duì)應(yīng)的預(yù)設(shè)代表風(fēng)化層速度值％為800米/秒；當(dāng)實(shí)際風(fēng)化層速度變化區(qū)間為450米/秒 1050米/秒時(shí)，對(duì)應(yīng)的預(yù)設(shè)代表風(fēng)化層速度值Vtl為750米/秒；當(dāng)實(shí)際風(fēng)化層速度變化區(qū)間為550米/秒 1150米/秒時(shí)，對(duì)應(yīng)的預(yù)設(shè)代表風(fēng)化層速度值Vtl為850米/秒。步驟S103 :確定為處理中使用的代表風(fēng)化層速度。判斷所述獲得的實(shí)際風(fēng)化層速度變化區(qū)間是否介于臨界風(fēng)化層速度區(qū)間中，若是，則將所述預(yù)設(shè)代表風(fēng)化層速度％確定為代表風(fēng)化層速度，若否，則重新預(yù)設(shè)代表風(fēng)化層速度Vo，直到可獲得代表風(fēng)化層速度。如，當(dāng)步驟S102中，由預(yù)設(shè)代表風(fēng)化層速度為800米/秒獲得的實(shí)際風(fēng)化層速度變化區(qū)間為500米/秒 1100米/秒時(shí)，判斷是否介于臨界風(fēng)化層速度區(qū)間30(Tl500米/秒(其中300米/秒是聲波在空氣中的傳播速度，1500米/秒是聲波在水中的傳播速度)區(qū)間內(nèi)，由于上述獲得的實(shí)際風(fēng)化層速度變化區(qū)間在臨界風(fēng)化層速度區(qū)間中，因此，確定預(yù)設(shè)代表風(fēng)化層速度值Vtl為800米/秒為處理中使用的代表風(fēng)化層速度。需要說(shuō)明的是，通過(guò)以上實(shí)施例中的方法可在臨界風(fēng)化層速度區(qū)間30(T1500米/秒中，通過(guò)設(shè)立多個(gè)預(yù)設(shè)代表風(fēng)化層速度Vtl，對(duì)應(yīng)獲得多個(gè)與之對(duì)應(yīng)的多個(gè)區(qū)間，并可根據(jù)實(shí)際需要通過(guò)調(diào)整誤差指數(shù)k的方式，對(duì)區(qū)間寬度進(jìn)行調(diào)整，同時(shí)，可將獲得值進(jìn)行驗(yàn)證，因此在步驟S103步驟還可包括當(dāng)所述判斷所述實(shí)際風(fēng)化層速度變化區(qū)間是否介于臨界風(fēng)化層速度區(qū)間中，若是，則將所述預(yù)設(shè)代表風(fēng)化層速度Vtl確定為處理中所用的代表風(fēng)化層速度后；根據(jù)所述處理中所用的代表風(fēng)化層速度進(jìn)行靜校正處理，判斷所述剩余靜校正值St是否在剩余靜校正值邊界范圍值取值內(nèi)，若是，則保留所述處理中所用的代表風(fēng)化層速度，若否，則重新預(yù)設(shè)代表風(fēng)化層速度％，對(duì)獲取的處理中所用的代表風(fēng)化層速度在進(jìn)一步給予驗(yàn)證，從而獲得滿足要求的代表風(fēng)化層速度。為進(jìn)一步說(shuō)明實(shí)施例步驟S102中所述誤差曲線公式的獲得，以下將結(jié)合附圖給予詳細(xì)說(shuō)明圖2為本發(fā)明的一個(gè)近地表模型示意圖。其中D為基準(zhǔn)面高程，S為地表面高程，R為假定風(fēng)化層正確時(shí)對(duì)應(yīng)的高速頂界高程，R’為假定風(fēng)化層錯(cuò)誤時(shí)對(duì)應(yīng)的高速頂界高程，h為風(fēng)化層厚度，dh為風(fēng)化層厚度誤差。根據(jù)折射理論，折射層截距時(shí)t與風(fēng)化層厚度h間的關(guān)系由下式表示 其中，i=Vtl是風(fēng)化層速度，V1是折射層速度。由任意物理點(diǎn)的實(shí)際風(fēng)化層速度得到對(duì)應(yīng)的風(fēng)化層厚度為
權(quán)利要求
1.一種優(yōu)選風(fēng)化層速度的確定方法，其特征在于，包括 預(yù)設(shè)代表風(fēng)化層速度Vtl ； 根據(jù)所述預(yù)設(shè)代表風(fēng)化層速度Vtl獲取實(shí)際風(fēng)化層速度變化區(qū)間； 判斷所述實(shí)際風(fēng)化層速度變化區(qū)間是否介于臨界風(fēng)化層速度區(qū)間中，若是，則將所述預(yù)設(shè)代表風(fēng)化層速度Vtl確定為處理中所用的代表風(fēng)化層速度，若否，則重新預(yù)設(shè)代表風(fēng)化層速度Vf!。
2.如權(quán)利要求I所述的方法，其特征在于，所述根據(jù)所述預(yù)設(shè)代表風(fēng)化層速度Vtl獲取實(shí)際風(fēng)化層速度變化區(qū)間步驟還包括 根據(jù)剩余靜校正值S t及截距時(shí)t的邊界范圍值，確定所述剩余靜校正值S t、截距時(shí)t ； 將所述剩余靜校正值S t及截距時(shí)t代入如下誤差曲線公式中獲取誤差指數(shù)k ；
3.如權(quán)利要求I或2所述的方法，其特征在于，所述臨界風(fēng)化層速度區(qū)間為300米/秒 1500米/秒。
4.如權(quán)利要求I或2所述的方法,其特征在于,所述剩余靜校正值St的邊界范圍絕對(duì)值不超過(guò)20毫秒，所述截距時(shí)t的邊界范圍值不超過(guò)400毫秒。
5.如權(quán)利要求4所述的方法，其特征在于，所述將所述剩余靜校正值St及截距時(shí)t代入如下誤差曲線公式中獲取誤差指數(shù)k步驟還包括 當(dāng)所述剩余靜校正值S t為±20毫秒，及截距時(shí)t為400毫秒時(shí)，獲取所述誤差指數(shù)k等于±0. I。
6.如權(quán)利要求I中所述的方法，其特征在于，所述預(yù)設(shè)代表風(fēng)化層速度Vtl包括800米/秒、750米/秒 850米/秒。
7.如權(quán)利要求6所述的方法，其特征在于，所述將所述誤差指數(shù)k及所述預(yù)設(shè)代表風(fēng)化層速度％代入如下誤差曲線公式中分別獲取剩余風(fēng)化層速度dV(l步驟還包括 當(dāng)所述誤差指數(shù)k為±0. I、所述預(yù)設(shè)代表風(fēng)化層速度Vtl為800米/秒時(shí)，獲取剩余風(fēng)化層速度dvQ為±300米/秒。
8.如權(quán)利要求7所述的方法，其特征在于，所述根據(jù)所述預(yù)設(shè)代表風(fēng)化層速度Vtl及所述剩余風(fēng)化層速度dV(l獲取實(shí)際風(fēng)化層速度變化區(qū)間步驟還包括 當(dāng)所述預(yù)設(shè)代表風(fēng)化層速度Vtl為800米/秒，所述剩余風(fēng)化層速度dV(l為±300米/秒時(shí)，獲取所述實(shí)際風(fēng)化層速度變化區(qū)間為500米/秒 1100米/秒；或當(dāng)所述預(yù)設(shè)代表風(fēng)化層速度Vtl為750米/秒，所述剩余風(fēng)化層速度dV(l為±300米/秒時(shí)，獲取所述實(shí)際風(fēng)化層速度變化區(qū)間為450米/秒 1050米/秒； 或當(dāng)所述預(yù)設(shè)代表風(fēng)化層速度％為850米/秒，所述剩余風(fēng)化層速度dV(l為±300米/秒時(shí)，獲取所述實(shí)際風(fēng)化層速度變化區(qū)間為550米/秒 1150米/秒。
9.如權(quán)利要求I所述的方法，其特征在于，所述判斷所述實(shí)際風(fēng)化層速度變化區(qū)間是否介于臨界風(fēng)化層速度區(qū)間中，若是，則將所述預(yù)設(shè)代表風(fēng)化層速度％確定為處理中所用的代表風(fēng)化層速度步驟還包括 判斷所述實(shí)際風(fēng)化層速度變化區(qū)間是否介于臨界風(fēng)化層速度區(qū)間中，若是，則將所述預(yù)設(shè)代表風(fēng)化層速度Vtl確定為處理中所用的代表風(fēng)化層速度； 根據(jù)所述處理中所用的代表風(fēng)化層速度進(jìn)行靜校正處理，判斷所述剩余靜校正值S t是否在剩余靜校正值邊界范圍值取值內(nèi)，若是，則保留所述處理中所用的代表風(fēng)化層速度，若否，則重新預(yù)設(shè)代表風(fēng)化層速度％。
全文摘要
本發(fā)明實(shí)施例公開(kāi)了一種優(yōu)選風(fēng)化層速度的確定方法，具體包括預(yù)設(shè)一個(gè)代表風(fēng)化層速度值，該值在風(fēng)化層厚度不小于300米情況下，與實(shí)際風(fēng)化層速度差引起的靜校正量完全可限制在后續(xù)剩余靜校正程序允許的范圍以內(nèi)，近地表模型的建立就可以獲得穩(wěn)定解。因此，在折射波靜校正過(guò)程中，不再需要預(yù)先進(jìn)行諸如微測(cè)井等的近地表調(diào)查，就可使復(fù)雜地區(qū)構(gòu)造成像而必需進(jìn)行的短、長(zhǎng)波長(zhǎng)折射波靜校正在精度方面得到提高。
文檔編號(hào)G01V1/36GK102749649SQ201210237029
公開(kāi)日2012年10月24日 申請(qǐng)日期2012年7月9日 優(yōu)先權(quán)日2012年4月10日
發(fā)明者吳媚, 孫庚文, 歐陽(yáng)行艷, 王順根, 謝桂生 申請(qǐng)人:恒泰艾普石油天然氣技術(shù)服務(wù)股份有限公司