獲取風(fēng)化層的地質(zhì)參數(shù)的方法
【專利摘要】本發(fā)明提供一種獲取風(fēng)化層的地質(zhì)參數(shù)的方法，所述方法包括：(a)獲取目標(biāo)區(qū)域的標(biāo)準(zhǔn)的初至數(shù)據(jù)；(b)確定最小偏移距；(c)確定每個(gè)點(diǎn)的風(fēng)化層速度值；(d)計(jì)算選擇的m個(gè)點(diǎn)中的任意兩點(diǎn)之間的風(fēng)化層速度值的變差函數(shù)值；(e)計(jì)算選擇的n個(gè)點(diǎn)中的任意兩點(diǎn)之間的微測(cè)井資料的速度值的變差函數(shù)值；(f)計(jì)算m個(gè)點(diǎn)中的任意一點(diǎn)與n個(gè)點(diǎn)中的任意一點(diǎn)之間的關(guān)于速度的交叉編譯函數(shù)；(g)獲得與m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值對(duì)應(yīng)的m個(gè)加權(quán)系數(shù)和與n個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的速度值對(duì)應(yīng)的n個(gè)加權(quán)系數(shù)；(h)計(jì)算風(fēng)化層速度的估計(jì)值。采用本發(fā)明獲取風(fēng)化層的地質(zhì)參數(shù)的方法可得到更符合地質(zhì)分布規(guī)律的淺層速度結(jié)構(gòu)。
【專利說(shuō)明】獲取風(fēng)化層的地質(zhì)參數(shù)的方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明總體來(lái)說(shuō)涉及石油地震資料分析處理領(lǐng)域，更具體地講，涉及一種基于協(xié) 同克里金技術(shù)融合微測(cè)井資料的獲取風(fēng)化層的地質(zhì)參數(shù)的方法。

【背景技術(shù)】
[0002] 在地震勘探中，近地表的低降速帶（S卩，風(fēng)化層）的速度和厚度不均勻變化會(huì)引起 地震波傳播旅行時(shí)的不均勻延遲，致使實(shí)際觀測(cè)得到的反射波到達(dá)時(shí)間發(fā)生畸變，極大地 影響著地震剖面的成像和構(gòu)造形態(tài)。靜校正就是消除地表低降速帶對(duì)勘探資料的影響，對(duì) 于復(fù)雜地表?xiàng)l件下地震資料的處理，求取正確的靜校正量是非常關(guān)鍵和必要的。如果不能 很好的解決嚴(yán)重的靜校正問(wèn)題，就會(huì)在地震剖面上產(chǎn)生假象或不成像，就會(huì)給地質(zhì)解釋帶 來(lái)嚴(yán)重困難，甚至錯(cuò)誤的地質(zhì)構(gòu)造認(rèn)識(shí)。
[0003] 目前，在解決靜校正問(wèn)題過(guò)程中，對(duì)于存在穩(wěn)定折射層的區(qū)域，一般采用折射靜校 正方法（例如，基于迭代法的折射靜校正方法和廣義互換靜校正方法）進(jìn)行計(jì)算。但這兩 類方法在使用時(shí)，對(duì)于風(fēng)化層速度，一般參照工作經(jīng)驗(yàn)，給定常速度值進(jìn)行計(jì)算。然而，很多 地區(qū)的風(fēng)化層速度橫向變化劇烈，實(shí)際風(fēng)化層速度與前述參照工作經(jīng)驗(yàn)給定的常速度值存 在較大差距。
[0004] 此外，有的方法還會(huì)對(duì)實(shí)際風(fēng)化層速度進(jìn)行擬合計(jì)算或者結(jié)合微測(cè)井資料進(jìn)行速 度反演，其效果優(yōu)于給定常速度值計(jì)算的結(jié)果。但采用速度擬合的方法仍存在地震剖面的 成像橫向不連續(xù)現(xiàn)象，需要進(jìn)一步處理；而進(jìn)行速度反演，則計(jì)算時(shí)間長(zhǎng)，微測(cè)井資料的好 壞將直接影響速度反演的結(jié)果，不確定程度高，從而影響處理后疊加剖面成像效果，使得構(gòu) 造形態(tài)不準(zhǔn)確，信噪比大大降低，不能達(dá)到資料精細(xì)處理的目的。


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

[0005] 本發(fā)明的目的是提出一種基于協(xié)同克里金技術(shù)融合微測(cè)井資料的獲取風(fēng)化層的 地質(zhì)參數(shù)的方法，通過(guò)變異函數(shù)這一數(shù)學(xué)工具，分別刻畫(huà)出微測(cè)井資料的速度值、厚度值和 由初至數(shù)據(jù)擬合出的風(fēng)化層速度、厚度的變化趨勢(shì)，再運(yùn)用交叉編譯函數(shù)得到兩類數(shù)據(jù)的 整體變化趨勢(shì)，并加以它們的相關(guān)系數(shù)約束，最終得到符合地質(zhì)分布規(guī)律的淺層速度結(jié)構(gòu)。
[0006] 本發(fā)明的一方面提供一種獲取風(fēng)化層的地質(zhì)參數(shù)的方法，所述方法包括：(a)獲 取目標(biāo)區(qū)域的標(biāo)準(zhǔn)的初至數(shù)據(jù)；(b)確定最小偏移距；（c)在所述目標(biāo)區(qū)域內(nèi)選擇m個(gè)點(diǎn)， 從獲取的初至數(shù)據(jù)中選擇m個(gè)點(diǎn)中的每個(gè)點(diǎn)在最小偏移距范圍內(nèi)的初至數(shù)據(jù)，并對(duì)選擇的 m個(gè)點(diǎn)中的每個(gè)點(diǎn)在最小偏移距范圍內(nèi)對(duì)應(yīng)的初至數(shù)據(jù)進(jìn)行擬合，確定每個(gè)點(diǎn)的風(fēng)化層速 度值，m為大于等于1的自然數(shù)；（d)根據(jù)確定的每個(gè)點(diǎn)的風(fēng)化層速度值，計(jì)算選擇的m個(gè) 點(diǎn)中的任意兩點(diǎn)之間的風(fēng)化層速度值的變差函數(shù)值；(e)在所述目標(biāo)區(qū)域內(nèi)選擇η個(gè)點(diǎn)，從 所述目標(biāo)區(qū)域的微測(cè)井資料中提取η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井資料的速度值，并根據(jù)提取 的每個(gè)點(diǎn)的微測(cè)井資料的速度值，計(jì)算選擇的η個(gè)點(diǎn)中的任意兩點(diǎn)之間的微測(cè)井資料的速 度值的變差函數(shù)值，η為大于等于1的自然數(shù)；（f)根據(jù)確定的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層 速度值和η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井資料的速度值，計(jì)算m個(gè)點(diǎn)中的任意一點(diǎn)與η個(gè)點(diǎn)中 的任意一點(diǎn)之間的關(guān)于速度的交叉編譯函數(shù)；（g)根據(jù)步驟（d)?步驟（f)的結(jié)果獲得與 m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值對(duì)應(yīng)的m個(gè)加權(quán)系數(shù)和與η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù) 據(jù)的速度值對(duì)應(yīng)的η個(gè)加權(quán)系數(shù)；(h)根據(jù)獲得的與m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值對(duì) 應(yīng)的m個(gè)加權(quán)系數(shù)、與η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的速度值對(duì)應(yīng)的η個(gè)加權(quán)系數(shù)、確定 的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值和η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井資料的速度值，計(jì)算風(fēng) 化層速度的估計(jì)值。
[0007] 可選地，步驟（g)可包括：根據(jù)步驟（d)?步驟（f)的結(jié)果利用協(xié)同克里金方程組 獲得與m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值對(duì)應(yīng)的m個(gè)加權(quán)系數(shù)和與η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微 測(cè)井?dāng)?shù)據(jù)的速度值對(duì)應(yīng)的η個(gè)加權(quán)系數(shù)。
[0008] 可選地，步驟（g)可包括：利用下面的協(xié)同克里金方程組計(jì)算與m個(gè)點(diǎn)中的每個(gè)點(diǎn) 的風(fēng)化層速度值對(duì)應(yīng)的m個(gè)加權(quán)系數(shù)和與η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的速度值對(duì)應(yīng)的 η個(gè)加權(quán)系數(shù)，

【權(quán)利要求】
1. 一種獲取風(fēng)化層的地質(zhì)參數(shù)的方法，所述方法包括： (a) 獲取目標(biāo)區(qū)域的標(biāo)準(zhǔn)的初至數(shù)據(jù)； (b) 確定最小偏移距； (c) 在所述目標(biāo)區(qū)域內(nèi)選擇m個(gè)點(diǎn)，從獲取的初至數(shù)據(jù)中選擇m個(gè)點(diǎn)中的每個(gè)點(diǎn)在最小 偏移距范圍內(nèi)的初至數(shù)據(jù)，并對(duì)選擇的m個(gè)點(diǎn)中的每個(gè)點(diǎn)在最小偏移距范圍內(nèi)對(duì)應(yīng)的初至 數(shù)據(jù)進(jìn)行擬合，確定每個(gè)點(diǎn)的風(fēng)化層速度值，m為大于等于1的自然數(shù)； (d) 根據(jù)確定的每個(gè)點(diǎn)的風(fēng)化層速度值，計(jì)算選擇的m個(gè)點(diǎn)中的任意兩點(diǎn)之間的風(fēng)化 層速度值的變差函數(shù)值； (e) 在所述目標(biāo)區(qū)域內(nèi)選擇η個(gè)點(diǎn)，從所述目標(biāo)區(qū)域的微測(cè)井資料中提取η個(gè)點(diǎn)中的 每個(gè)點(diǎn)的微測(cè)井資料的速度值，并根據(jù)提取的每個(gè)點(diǎn)的微測(cè)井資料的速度值，計(jì)算選擇的η 個(gè)點(diǎn)中的任意兩點(diǎn)之間的微測(cè)井資料的速度值的變差函數(shù)值，η為大于等于1的自然數(shù)； (f) 根據(jù)確定的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值和η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井資料 的速度值，計(jì)算m個(gè)點(diǎn)中的任意一點(diǎn)與η個(gè)點(diǎn)中的任意一點(diǎn)之間的關(guān)于速度的交叉編譯函 數(shù)； (g) 根據(jù)步驟（d)?步驟（f)的結(jié)果獲得與m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值對(duì)應(yīng)的 m個(gè)加權(quán)系數(shù)和與η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的速度值對(duì)應(yīng)的η個(gè)加權(quán)系數(shù)； (h) 根據(jù)獲得的與m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值對(duì)應(yīng)的m個(gè)加權(quán)系數(shù)、與η個(gè)點(diǎn)中 的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的速度值對(duì)應(yīng)的η個(gè)加權(quán)系數(shù)、確定的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層 速度值和η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井資料的速度值，計(jì)算風(fēng)化層速度的估計(jì)值。
2. 根據(jù)權(quán)利要求1所述的方法，其中，步驟（g)包括：根據(jù)步驟（d)?步驟（f)的結(jié)果 利用協(xié)同克里金方程組獲得與m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值對(duì)應(yīng)的m個(gè)加權(quán)系數(shù)和與 η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的速度值對(duì)應(yīng)的η個(gè)加權(quán)系數(shù)。
3. 根據(jù)權(quán)利要求2所述的方法，其中，步驟（g)包括： 利用下面的協(xié)同克里金方程組計(jì)算與m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值對(duì)應(yīng)的m個(gè)加 權(quán)系數(shù)和與η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的速度值對(duì)應(yīng)的η個(gè)加權(quán)系數(shù)，
其中,Cx (Xi,χρ)為η個(gè)點(diǎn)中的第i點(diǎn)的坐標(biāo)位置Xi與η個(gè)點(diǎn)中的第ρ點(diǎn)的坐標(biāo)位置χρ 對(duì)應(yīng)的微測(cè)井資料的速度值的變差函數(shù)值，Cx (0, Xi)為待求位置O與η個(gè)點(diǎn)中的第i點(diǎn)的坐 標(biāo)位置Xi對(duì)應(yīng)的微測(cè)井資料的速度值的變差函數(shù)值，Cy (y,，0為m個(gè)點(diǎn)中的第q點(diǎn)的坐標(biāo) 位置!與m個(gè)點(diǎn)中的第j點(diǎn)的坐標(biāo)位置yj對(duì)應(yīng)的風(fēng)化層速度值的變差函數(shù)值,Cy (0,yj為 待求位置0與m個(gè)點(diǎn)中的第j點(diǎn)的坐標(biāo)位置t對(duì)應(yīng)的風(fēng)化層速度值的變差函數(shù)值，CJbyj) 和Cc (yj，Xi)為關(guān)于η個(gè)點(diǎn)中的第i點(diǎn)的坐標(biāo)位置Xi對(duì)應(yīng)的微測(cè)井資料的速度值與m個(gè)點(diǎn) 中的第j點(diǎn)的坐標(biāo)位置&對(duì)應(yīng)的風(fēng)化層速度值的交叉編譯函數(shù)，為協(xié)同克里金方程組中 與η個(gè)點(diǎn)中的第i點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的速度值對(duì)應(yīng)的加權(quán)系數(shù)，/>\,為協(xié)同克里金方程組中與 m個(gè)點(diǎn)中的第j點(diǎn)的風(fēng)化層速度值對(duì)應(yīng)的加權(quán)系數(shù)，μχ、μy為計(jì)算方程組的解與原始值的 誤差的參數(shù)。
4. 根據(jù)權(quán)利要求1所述的方法，其中，步驟（h)包括： 利用下面的公式計(jì)算風(fēng)化層速度的估計(jì)值Z(O)，
其中，Zi(O)為風(fēng)化層速度在待求位置0處的估計(jì)值，Zx(Xi)為η個(gè)點(diǎn)中的第i點(diǎn)的坐 標(biāo)位置Xi對(duì)應(yīng)的微測(cè)井資料的速度值，A為與η個(gè)點(diǎn)中的第i點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的速度值對(duì) 應(yīng)的加權(quán)系數(shù)，Zy (yp為m個(gè)點(diǎn)中的第j點(diǎn)的坐標(biāo)位置&對(duì)應(yīng)的風(fēng)化層速度值，Py,.為與m 個(gè)點(diǎn)中的第j點(diǎn)的風(fēng)化層速度值對(duì)應(yīng)的加權(quán)系數(shù)。
5. 根據(jù)權(quán)利要求1所述的方法，其中，步驟（c)還包括：對(duì)選擇的m個(gè)點(diǎn)中的每個(gè)點(diǎn)在 最小偏移距范圍內(nèi)對(duì)應(yīng)的初至數(shù)據(jù)進(jìn)行擬合，確定每個(gè)點(diǎn)的炮點(diǎn)延遲時(shí)的初始值， 其中，所述方法還包括： (〇)確定最大偏移距； (P)從獲取的初至數(shù)據(jù)中選擇m個(gè)點(diǎn)中的每個(gè)點(diǎn)在最小偏移距至最大偏移距范圍內(nèi)的 初至數(shù)據(jù)，對(duì)選擇的m個(gè)點(diǎn)中的每個(gè)點(diǎn)在最小偏移距至最大偏移距范圍內(nèi)對(duì)應(yīng)的初至數(shù)據(jù) 進(jìn)行擬合，確定每個(gè)點(diǎn)的高速層速度的初始值； (q)根據(jù)確定的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的高速層速度的初始值和m個(gè)點(diǎn)中的每個(gè)點(diǎn)的炮點(diǎn) 延遲時(shí)的初始值，利用折射初至?xí)r間公式計(jì)算m個(gè)點(diǎn)中的每個(gè)點(diǎn)的檢波點(diǎn)延遲時(shí)、炮點(diǎn)延 遲時(shí)和高速層速度。
6. 根據(jù)權(quán)利要求5所述的方法，其中，步驟（q)包括： (ql)根據(jù)m個(gè)點(diǎn)中的每個(gè)點(diǎn)的高速層速度、炮點(diǎn)延遲時(shí)、初至數(shù)據(jù)中每個(gè)點(diǎn)對(duì)應(yīng)的炮 檢對(duì)從激發(fā)到接收產(chǎn)生的時(shí)間、每個(gè)點(diǎn)對(duì)應(yīng)的炮點(diǎn)與檢波點(diǎn)之間的距離，利用折射初至?xí)r 間公式計(jì)算m個(gè)點(diǎn)中的每個(gè)點(diǎn)的檢波點(diǎn)延遲時(shí)，其中，將m個(gè)點(diǎn)中的每個(gè)點(diǎn)的高速層速度的 初始值和m個(gè)點(diǎn)中的每個(gè)點(diǎn)的炮點(diǎn)延遲時(shí)的初始值作為初值代入； (q2)根據(jù)m個(gè)點(diǎn)中的每個(gè)點(diǎn)的高速層速度、步驟（ql)得到的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的檢波 點(diǎn)延遲時(shí)、初至數(shù)據(jù)中每個(gè)點(diǎn)對(duì)應(yīng)的炮檢對(duì)從激發(fā)到接收產(chǎn)生的時(shí)間、每個(gè)點(diǎn)對(duì)應(yīng)的炮點(diǎn) 與檢波點(diǎn)之間的距離，利用折射初至?xí)r間公式計(jì)算m個(gè)點(diǎn)中的每個(gè)點(diǎn)的炮點(diǎn)延遲時(shí)； (q3)根據(jù)步驟（ql)得到的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的檢波點(diǎn)延遲時(shí)、步驟（q2)得到的m個(gè) 點(diǎn)中的每個(gè)點(diǎn)的炮點(diǎn)延遲時(shí)、初至數(shù)據(jù)中每個(gè)點(diǎn)對(duì)應(yīng)的炮檢對(duì)從激發(fā)到接收產(chǎn)生的時(shí)間、 每個(gè)點(diǎn)對(duì)應(yīng)的炮點(diǎn)與檢波點(diǎn)之間的距離，利用折射初至?xí)r間公式計(jì)算m個(gè)點(diǎn)中的每個(gè)點(diǎn)的 高速層速度； (q4)根據(jù)步驟（ql)得到的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的檢波點(diǎn)延遲時(shí)、步驟（q2)得到的m個(gè) 點(diǎn)中的每個(gè)點(diǎn)的炮點(diǎn)延遲時(shí)、步驟（q3)得到的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的高速層速度、每個(gè)點(diǎn)對(duì) 應(yīng)的炮點(diǎn)與檢波點(diǎn)之間的距離，利用折射初至?xí)r間公式計(jì)算對(duì)應(yīng)的初至數(shù)據(jù)中每個(gè)點(diǎn)對(duì)應(yīng) 的炮檢對(duì)從激發(fā)到接收產(chǎn)生的時(shí)間； (q5)當(dāng)所述時(shí)間小于等于設(shè)定的時(shí)間范圍時(shí)，將步驟（ql)?步驟（q3)的結(jié)果作為m個(gè)點(diǎn)中的每個(gè)點(diǎn)的炮點(diǎn)延遲時(shí)、檢波點(diǎn)延遲時(shí)和高速層速度，當(dāng)所述時(shí)間大于設(shè)定的時(shí)間 范圍時(shí)，則返回執(zhí)行步驟（ql)。
7. 根據(jù)權(quán)利要求6所述的方法，其中，折射初至?xí)r間公式為，
其中，為初至數(shù)據(jù)中m個(gè)點(diǎn)中的第j點(diǎn)對(duì)應(yīng)的炮檢對(duì)從激發(fā)到接收產(chǎn)生的時(shí)間，TA(yj)為m個(gè)點(diǎn)中的第j點(diǎn)的坐標(biāo)位置y』對(duì)應(yīng)的炮點(diǎn)延遲時(shí)，TB(yj)為m個(gè)點(diǎn)中的第j點(diǎn) 的坐標(biāo)位置yj對(duì)應(yīng)的檢波點(diǎn)延遲時(shí)，/iZ?/為m個(gè)點(diǎn)中的第j點(diǎn)對(duì)應(yīng)的炮點(diǎn)與檢波點(diǎn)之間的 距離，S(yj)為m個(gè)點(diǎn)中的第j點(diǎn)的坐標(biāo)位置yj對(duì)應(yīng)的高速層速度。
8. 根據(jù)權(quán)利要求5所述的方法，所述方法還包括： (k)根據(jù)m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層速度值、m個(gè)點(diǎn)中的每個(gè)點(diǎn)的高速層速度、m個(gè)點(diǎn) 中的每個(gè)點(diǎn)的炮點(diǎn)延遲時(shí)或者檢波點(diǎn)延遲時(shí)，計(jì)算m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層厚度值； 其中，步驟（d)還包括：根據(jù)確定的每個(gè)點(diǎn)的風(fēng)化層厚度值，計(jì)算選擇的m個(gè)點(diǎn)中的任 意兩點(diǎn)之間的風(fēng)化層厚度值的變差函數(shù)值； 步驟（e)還包括：從所述目標(biāo)區(qū)域的微測(cè)井資料中提取η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井資 料的厚度值，并根據(jù)提取的每個(gè)點(diǎn)的微測(cè)井資料的厚度值，計(jì)算選擇的η個(gè)點(diǎn)中的任意兩 點(diǎn)之間的微測(cè)井資料的厚度值的變差函數(shù)值； 步驟（f)還包括：根據(jù)確定的m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層厚度值和η個(gè)點(diǎn)中的每個(gè)點(diǎn) 的微測(cè)井資料的厚度值，計(jì)算m個(gè)點(diǎn)中的任意一點(diǎn)與η個(gè)點(diǎn)中的任意一點(diǎn)之間的關(guān)于厚度 的交叉編譯函數(shù)； 步驟（g)還包括：根據(jù)步驟（d)?步驟（f)的結(jié)果獲得與m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層 厚度值對(duì)應(yīng)的m個(gè)加權(quán)系數(shù)和與η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的厚度值對(duì)應(yīng)的η個(gè)加權(quán) 系數(shù)； 步驟（h)還包括：根據(jù)獲得的與m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層厚度值對(duì)應(yīng)的m個(gè)加權(quán)系 數(shù)、與η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井?dāng)?shù)據(jù)的厚度值對(duì)應(yīng)的η個(gè)加權(quán)系數(shù)、確定的m個(gè)點(diǎn)中的每 個(gè)點(diǎn)的風(fēng)化層厚度值和η個(gè)點(diǎn)中的每個(gè)點(diǎn)的微測(cè)井資料的厚度值，計(jì)算風(fēng)化層厚度的估計(jì) 值。
9. 根據(jù)權(quán)利要求8所述的方法，其中，步驟（k)包括： 利用下面的公式計(jì)算m個(gè)點(diǎn)中的每個(gè)點(diǎn)的風(fēng)化層厚度值，
其中，Dy (yp為m個(gè)點(diǎn)中的第j點(diǎn)的坐標(biāo)位置t對(duì)應(yīng)的風(fēng)化層厚度值，T(yp為m個(gè)點(diǎn) 中的第j點(diǎn)的坐標(biāo)位置yj對(duì)應(yīng)的炮點(diǎn)延遲時(shí)或檢波點(diǎn)延遲時(shí)，Zy (y」）為m個(gè)點(diǎn)中的第j點(diǎn) 的坐標(biāo)位置y」對(duì)應(yīng)的風(fēng)化層速度值，S(y」）為m個(gè)點(diǎn)中的第j點(diǎn)的坐標(biāo)位置y」對(duì)應(yīng)的高速 層速度。
10.根據(jù)權(quán)利要求8所述的方法，所述方法還包括：(m)根據(jù)風(fēng)化層速度的估計(jì)值和風(fēng) 化層厚度的估計(jì)值，計(jì)算折射靜校正量。
【文檔編號(hào)】G01V1/36GK104316961SQ201410612194
【公開(kāi)日】2015年1月28日 申請(qǐng)日期:2014年11月4日 優(yōu)先權(quán)日:2014年11月4日
【發(fā)明者】張亨, 陳愛(ài)萍, 何光明, 劉奇琳, 曹中林, 周強(qiáng), 漆韜, 李若禹 申請(qǐng)人:中國(guó)石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司