擬合系數(shù)進行插值以得到各個電磁測點的擬合系數(shù)，從 而計算出各個電測測點的厚度和速度信息，以實現(xiàn)對黃土層的靜校正。通過上述方式解決 了如有想要保證靜校正的精度，就必然要增加大量成本的技術(shù)問題，達到了在保證靜校正 結(jié)果準確性的同時降低成本的目的。
【附圖說明】
[0028] 此處所說明的附圖用來提供對本發(fā)明的進一步理解，構(gòu)成本申請的一部分，并不 構(gòu)成對本發(fā)明的限定。在附圖中：
[0029] 圖1是本發(fā)明實施例的黃土層靜校正的方法流程圖；
[0030] 圖2是本發(fā)明實施例的二維音頻大地電磁勘探數(shù)據(jù)和微測井數(shù)據(jù)結(jié)合進行二維 地震剖面的黃土層靜校正的方法流程圖；
[0031] 圖3是本發(fā)明實施例的黃土層靜校正的裝置結(jié)構(gòu)框圖。
【具體實施方式】
[0032] 為使本發(fā)明的目的、技術(shù)方案和優(yōu)點更加清楚明白，下面結(jié)合實施方式和附圖，對 本發(fā)明做進一步詳細說明。在此，本發(fā)明的示意性實施方式及其說明用于解釋本發(fā)明，但并 不作為對本發(fā)明的限定。
[0033] 考慮到現(xiàn)有的基于微測井方法進行靜校正時，如果需要較準確地確定測點處表層 速度隨深度變化的規(guī)律，就需要增加大量測點，如果測點太少，則無法為地震資料提供可靠 的表層模型。在本例中，主要是考慮到可以按照一定的規(guī)則增加測點，從而使得在保證準確 度的情況下，盡可能少的增加測點，從而降低成本。具體地，通過在已有的微測井點的基礎(chǔ) 上，引入了電磁測點，并進行一系列處理后，可以對工區(qū)內(nèi)的地震資料進行更加準確的黃土 層靜校正。具體地，在本例中，提供了一種黃土層靜校正方法，如圖1所示，可以包括以下步 驟：
[0034] 步驟101 :采集工區(qū)的電磁數(shù)據(jù)，其中，所述工區(qū)中有Μ個電磁測點，所述Μ個電磁 測點中有Ν個電磁測點為微測井點，有Ζ個電磁測點為深微測井點，Μ、Ν、Ζ為正整數(shù)；
[0035] 在實際應(yīng)用中，可以根據(jù)地震工區(qū)的勘探范圍和地震資料對黃土層靜校正數(shù)據(jù)密 度的要求，設(shè)計并采集相應(yīng)密度的電磁數(shù)據(jù)。其中，相應(yīng)密度的電磁數(shù)據(jù)可以包括：音頻大 地電磁、瞬變電磁、大地電磁、連續(xù)電磁剖面，這些數(shù)據(jù)可以是二維電磁數(shù)據(jù)，也可以是三維 電磁數(shù)據(jù)。
[0036] 上述工區(qū)中的電磁測點可以按照以下規(guī)則布設(shè)：根據(jù)電磁測點的十分之一到二十 分之一的密度，在地震工區(qū)內(nèi)布設(shè)常規(guī)微測井，在布設(shè)的過程中，各個微測井的深度可以根 據(jù)實際要求確定。布設(shè)的這些常規(guī)微測井在工區(qū)內(nèi)相對均勻分布。具體地，當電磁數(shù)據(jù)為 二維測線時，可以在測線的起點和終點處各布設(shè)一口微測井，當電磁數(shù)據(jù)為三維測網(wǎng)時，可 以在測網(wǎng)的拐角和邊界處布設(shè)一定數(shù)量的微測井。進一步的，還可以在工區(qū)內(nèi)布設(shè)2 口鉆 穿黃土層的深微測井。
[0037] 步驟102:對所述電磁數(shù)據(jù)進行反演，得到各個電磁測點的黃土層的反演厚度和 反演電阻率；
[0038] 具體地，當電磁數(shù)據(jù)為二維電磁數(shù)據(jù)時，對該數(shù)據(jù)進行二維反演處理，當電磁數(shù)據(jù) 為三維電磁數(shù)據(jù)時，對該數(shù)據(jù)進行三維反演處理，從而獲得地震工區(qū)內(nèi)各個電磁測點的黃 土層的反演厚度和反演電阻率。
[0039] 步驟103:根據(jù)N個微測井點的微測井數(shù)據(jù)，得到各個微測井點的黃土層厚度和黃 土層速度；
[0040] 具體地，可以采用時距圖法求取各個微測井點的黃土層厚度和黃土層速度。
[0041] 時距圖表示地震波到達各觀測點的時間與震源至觀測點間距離關(guān)系的圖形。它是 運動地震學(xué)研究的基本內(nèi)容。對時距圖的分析和處理，可以對地質(zhì)得出定性定量的解釋。
[0042] 步驟104:對Z個深微測井點進行黃土層電阻率標定；其中，這Z個深微測井點是 鉆穿黃土層的深微測井。
[0043] 步驟105 :根據(jù)所述黃土層電阻率的標定結(jié)果，對各個電磁測點的反演電阻率進 行追蹤，得到所述各個電磁測點的黃土層厚度；
[0044] 步驟106 :檢測各個微測井點的電阻率，并根據(jù)各個微測井點的黃土層速度和電 阻率，建立速度與電阻率擬合關(guān)系式，并求取各個微測井點不同階的多項式擬合系數(shù)；
[0045] 具體的，根據(jù)各個微測井點的黃土層速度和電阻率，所建立的速度與電阻率擬合 關(guān)系式如下：
[0047] 其中，Vlgl表示第i個微測井點在深度Dig處的速度，a1表示擬合系數(shù)，Rigi表示 第i個微測井點在深度Dig處的電阻率，ηι表示第i個微測井點的多項式擬合總階數(shù)，1^表 示微測井點的個數(shù)。
[0048] 將步驟103所得到的各個微測井點的黃土層速度和黃土層厚度，以及步驟104中 深微測井點所得到的各個微測井點的黃土層的電阻率代入所建立的速度與電阻率擬合關(guān) 系式，可以得到地震工區(qū)內(nèi)每個微測井點所對應(yīng)的1至5階多項式擬合系數(shù) &1。
[0049] 步驟107 :將各個微測井點的黃土層的反演厚度、反演電阻率和不同階的多項式 擬合系數(shù)代入所述速度與電阻率擬合關(guān)系式，得到各個微測井點不同階的多項式擬合系數(shù) 所對應(yīng)的擬合速度，并根據(jù)所述擬合速度確定統(tǒng)一多項式擬合階數(shù)，以求取各個微測井點 在該統(tǒng)一多項式擬合階數(shù)下的擬合系數(shù)；
[0050] 具體地，可以按照以下方式求取各個微測井點在統(tǒng)一多項式擬合階數(shù)下的擬合系 數(shù)，包括：
[0051] 步驟1 :將各個微測井點的黃土層的反演厚度作為微測井點的深度、反演電阻率 作為在微測井點的深度處的電阻率、不同階的多項式擬合系數(shù)作為擬合系數(shù)，代入速度與 電阻率擬合關(guān)系式（即上述公式1)，得到各個微測井點不同階的多項式擬合系數(shù)所對應(yīng)的 擬合速度；
[0052] 通過將N個微測井點在步驟102中所得到的黃土層反演厚度和反演電阻率以及步 驟106中所得到的地震工區(qū)內(nèi)各個微測井點所對應(yīng)的1至5階多項式擬合系數(shù)代入公 式1，可以得到各個微測井點1至5階多項式擬合系數(shù)所對應(yīng)的擬合速度。
[0053] 步驟2 :根據(jù)各個微測井點不同階的多項式擬合系數(shù)所對應(yīng)的擬合速度，計算各 個微測井點不同多項式擬合階數(shù)下的相關(guān)系數(shù)；
[0054] 具體地，可以按照以下公式計算各個微測井點不同多項式擬合階數(shù)下的相關(guān)系 數(shù)：
[0056] 其中，Q為第i個微測井點的相關(guān)系數(shù)，Nk為每個微測井點的深度點數(shù)，：^為第 i個微測井點第k個深度點擬合后的速度，為第i個微測井點所有深度點擬合后的速度 平均值，F(xiàn)/g為第i個微測井點第k個深度點的速度，％為第i個微測井點所有深度點的 速度平均值。
[0057] 通過將步驟1所得到的擬合速度代入公式2中，可以得到各個微測井點1至5階 多項式擬合階數(shù)下的相關(guān)系數(shù)。
[0058] 步驟3 :對每個微測井，選取最大相關(guān)系數(shù)所對應(yīng)的多項式擬合階數(shù)作為該微測 井的最大多項式擬合階數(shù)，以得到N個最大多項式擬合階數(shù)；
[0059] 在本實施例中，對于N個微測井中的每個微測井均有1至5階共5個多項式擬合 階數(shù)，對于每個微測井都選取其最大相關(guān)系數(shù)所對應(yīng)的多項式的擬合階數(shù)，作為該微測井 的最大多項式擬合階數(shù)，也就相當于在該步驟中會得到N個最大多項式擬合階數(shù)，每個微 測井有一個最大多項式擬合階數(shù)。
[0060] 步驟4 :根據(jù)所述N個最大多項式擬合階數(shù)計算得到統(tǒng)一多項式擬合階數(shù)；
[0061] 具體地，可以按照以下公式計算統(tǒng)一多項式擬合階數(shù)：
[0063] 其中，η表示所述統(tǒng)一多項式擬合階數(shù)，Ns表示微測井點的個數(shù)，ni表示第i個微 測井點的最大多項式擬合階數(shù)。
[0064] 步驟5 :將所述統(tǒng)一多項式擬合階數(shù)作為所述N個微測井點的多項式擬合階數(shù)，代 入所述速度與電阻率擬合關(guān)系式，得到所述各個微測井點在該多項式擬合階數(shù)下的擬合系 數(shù)。
[0065] 具體的，將步驟4中所得到的統(tǒng)一多項式擬合階數(shù)作為一個固定值代入公式1，并 分別將步驟103所得到的各個微測井點的黃土層速度和黃土層厚度，以及步驟104中深微 測井點所得到的各個微測井點的黃土層的電阻率代入公式1即速度與電阻率擬合關(guān)系式， 可以得到地震工區(qū)內(nèi)每個微測井點所對應(yīng)的在該統(tǒng)一多項式擬合階數(shù)下的擬合系數(shù)％。
[0066] 步驟108 :根據(jù)所述各個微測井點在該統(tǒng)一多項式擬合階數(shù)下的擬合系數(shù)，進行 線性插值得到所述Μ個電磁測點中各個電磁測點的擬合系數(shù)；
[0067] 例如，Ν取值為10,Μ取值為100,Ζ的取值為2,就相當于通過2個深微測井點，以 及微測井數(shù)據(jù)，計算得到10個微測井點的擬合系數(shù)，然后基于這10個微測井點在統(tǒng)一的多 項式階數(shù)下的擬合系數(shù)進行插值，得到100個電磁測點的擬合系數(shù)。之所采用這種方式，主 要是考慮到僅采用微測井方法進行靜校正時，由于測點較少，因而無法為地震資料提供可 靠的表層模型。本實施例在已有的Ν個微測井點的基礎(chǔ)上，引入了Ζ個深微測井點，根據(jù)已 有的Ν個微測井點在統(tǒng)一多項式擬合階數(shù)下的擬合系數(shù)，對其進行線性插值，從而可以得 到本實施例中的Μ個電磁測點中各個電磁測點的擬合系數(shù)。
[0068] 具體地，在插值的過程中，如果電磁數(shù)據(jù)為二維電磁數(shù)據(jù)，那么就可以進行一維線 性插值得到Μ個電磁測點中各個電磁測點的擬合系數(shù)，如果電磁數(shù)據(jù)為三維電磁數(shù)據(jù)，那 么就可以進行二維線性插值得到Μ個電磁測點中各個電磁測點的擬合系數(shù)。