黃土層靜校正方法及裝置的制造方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及石油勘探技術(shù)領(lǐng)域，特別涉及一種黃土層靜校正方法及裝置。
【背景技術(shù)】
[0002] 靜校正問題作為地質(zhì)勘探中經(jīng)常出現(xiàn)的問題，一直是地質(zhì)勘探的難點，目前常有 的靜校正方法主要有：小折射法、微測井法和折射靜校正方法。
[0003] 然而，上述幾種靜校正方法都存在不同層面的問題，例如，對于小折射法，在低速 帶調(diào)查中，由于排列的長度有限，導致觀測到的初至折射波的偏移距較小，難以全面反映低 速帶底界的形態(tài)，因此，對地形起伏劇烈或者是表層速度變化較大的地區(qū)并不適用；對于微 測井方法，能夠較準確地確定測點處表層速度隨深度變化的規(guī)律，但對于低速帶厚度較大 或者是速度變化較大的工區(qū)，必須增加測點數(shù)量，測點數(shù)量的增多必然導致成本升高；對于 折射靜校正方法，僅適用于地表起伏不大、表層速度橫向均質(zhì)性較好、存在明顯折射界面且 界面比較平滑的地區(qū)。
[0004] 此外，還有一些其它的靜校正方法，例如：基于淺層折射和統(tǒng)計分析的靜校正方 法，具體地，包括：交互迭代靜校正、多域迭代靜校正、模擬退火靜校正等。然而，這些靜校正 方法對野外數(shù)據(jù)的精度有很高的要求，因此也無法廣泛應用。
[0005] 針對如何進行低成本準確性高的靜校正，目前尚未提出有效的解決方案。

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

[0006] 本發(fā)明實施例提供了一種黃土層靜校正方法及裝置，以達到在保證靜校正結(jié)果準 確性的同時降低成本的目的。
[0007] 本發(fā)明實施例提供了一種黃土層靜校正方法，包括：采集工區(qū)的電磁數(shù)據(jù)，其中， 所述工區(qū)中有Μ個電磁測點，所述Μ個電磁測點中有N個電磁測點為微測井點，有Z個電磁 測點為深微測井點，Μ、Ν、Ζ為正整數(shù)；對所述電磁數(shù)據(jù)進行反演，得到各個電磁測點的黃土 層的反演厚度和反演電阻率；根據(jù)Ν個微測井點的微測井數(shù)據(jù)，得到各個微測井點的黃土 層厚度和黃土層速度；對Ζ個深微測井點進行黃土層電阻率標定；根據(jù)所述黃土層電阻率 的標定結(jié)果，對各個電磁測點的反演電阻率進行追蹤，得到所述各個電磁測點的黃土層厚 度；檢測各個微測井點的電阻率，并根據(jù)各個微測井點的黃土層速度和電阻率，建立速度與 電阻率擬合關(guān)系式，并求取各個微測井點不同階的多項式擬合系數(shù)；將各個微測井點的黃 土層的反演厚度、反演電阻率和不同階的多項式擬合系數(shù)代入所述速度與電阻率擬合關(guān)系 式，得到各個微測井點不同階的多項式擬合系數(shù)所對應的擬合速度，并根據(jù)所述擬合速度 確定統(tǒng)一多項式擬合階數(shù)，以求取各個微測井點在該統(tǒng)一多項式擬合階數(shù)下的擬合系數(shù)； 根據(jù)所述各個微測井點在該統(tǒng)一多項式擬合階數(shù)下的擬合系數(shù)，進行線性插值得到所述Μ 個電磁測點中各個電磁測點的擬合系數(shù)；將各個電磁測點的擬合系數(shù)、黃土層的反演厚度 和反演電阻率，代入所述擬合關(guān)系式，得到各個電磁測點的速度；根據(jù)所述各個電磁測點的 黃土層厚度和各個電磁測點的速度，進行黃土層靜校正。
[0008] 在一個實施例中，所述速度與電阻率擬合關(guān)系式為：
[0010]其中，Vlgl表示第i個微測井點在深度Dig處的速度，a1表示擬合系數(shù)，Rigi表示 第i個微測井點在深度Dig處的電阻率，ηι表示第i個微測井點的多項式擬合總階數(shù)，1^表 示微測井點的個數(shù)。
[0011] 在一個實施例中，將各個微測井點的黃土層的反演厚度、反演電阻率和不同階的 多項式擬合系數(shù)代入所述速度與電阻率擬合關(guān)系式，得到各個微測井點不同階的多項式擬 合系數(shù)所對應的擬合速度，并根據(jù)所述擬合速度確定統(tǒng)一多項式擬合階數(shù)，以求取各個微 測井點在該統(tǒng)一多項式擬合階數(shù)下的擬合系數(shù)，包括：將各個微測井點的黃土層的反演厚 度作為微測井點的深度、反演電阻率作為在微測井點的深度處的電阻率、不同階的多項式 擬合系數(shù)作為擬合系數(shù)，代入所述速度與電阻率擬合關(guān)系式，得到各個微測井點不同階的 多項式擬合系數(shù)所對應的擬合速度；根據(jù)各個微測井點不同階的多項式擬合系數(shù)所對應的 擬合速度，計算各個微測井點不同多項式擬合階數(shù)下的相關(guān)系數(shù)；對每個微測井，選取最大 相關(guān)系數(shù)所對應的多項式擬合階數(shù)作為該微測井的最大多項式擬合階數(shù)，以得到N個最大 多項式擬合階數(shù)；根據(jù)所述N個最大多項式擬合階數(shù)計算得到統(tǒng)一多項式擬合階數(shù)：將所 述統(tǒng)一多項式擬合系數(shù)作為所述N個微測井點的多項式擬合階數(shù)，代入所述速度與電阻率 擬合關(guān)系式，得到所述各個微測井點在該多項式擬合階數(shù)下的擬合系數(shù)。
[0012] 在一個實施例中，按照以下公式計算各個微測井點不同多項式擬合階數(shù)下的相關(guān) 系數(shù)：
[0014] 其中，Q為第i個微測井點的相關(guān)系數(shù)，Nk為每個微測井點的深度點數(shù)，U為第 i個微測井點第k個深度點擬合后的速度，歹為第i個微測井點所有深度點擬合后的速度 平均值，F(xiàn)/g為第i個微測井點第k個深度點的速度，％，為第i個微測井點所有深度點的 速度平均值。
[0015] 在一個實施例中，按照以下公式計算得到統(tǒng)一多項式擬合階數(shù)：
[0017] 其中，η表示所述統(tǒng)一多項式擬合階數(shù)，Ns表示微測井點的個數(shù)，n i表示第i個微 測井點的最大多項式擬合階數(shù)。
[0018] 在一個實施例中，所述電磁數(shù)據(jù)包括：二維電磁數(shù)據(jù)，或者三維電磁數(shù)據(jù)；相應 地，當所述電磁數(shù)據(jù)為二維電磁數(shù)據(jù)時，進行一維線性插值得到所述Μ個電磁測點中各個 電磁測點的擬合系數(shù)；當所述電磁數(shù)據(jù)為三維電磁數(shù)據(jù)時，進行二維線性插值得到所述Μ 個電磁測點中各個電磁測點的擬合系數(shù)。
[0019] 本發(fā)明實施例還提供了一種黃土層靜校正裝置，包括：采集單元，采集工區(qū)的電磁 數(shù)據(jù)，其中，所述工區(qū)中有Μ個電磁測點，所述Μ個電磁測點中有N個電磁測點為微測井點， 有Ζ個電磁測點為深微測井點，Μ、Ν、Ζ為正整數(shù)；反演單元，對所述電磁數(shù)據(jù)進行反演，得到 各個電磁測點的黃土層的反演厚度和反演電阻率；黃土層數(shù)據(jù)獲取單元，根據(jù)Ν個微測井 點的微測井數(shù)據(jù)，得到各個微測井點的黃土層厚度和黃土層速度；標定單元，對Ζ個深微測 井點進行黃土層速度標定和電阻率標定；追蹤單元，根據(jù)所述黃土層電阻率的標定結(jié)果，對 各個電磁測點的反演電阻率進行追蹤，得到所述各個電磁測點的黃土層厚度；擬合關(guān)系式 建立單元，檢測各個微測井點的電阻率，并根據(jù)各個微測井點的黃土層速度和電阻率，建立 速度與電阻率擬合關(guān)系式，并求取各個微測井點不同階的多項式擬合系數(shù)；第一擬合系數(shù) 獲取單元，將各個微測井點的黃土層的反演厚度、反演電阻率和不同階的多項式擬合系數(shù) 代入所述速度與電阻率擬合關(guān)系式，得到各個微測井點不同階的多項式擬合系數(shù)所對應的 擬合速度，并根據(jù)所述擬合速度確定統(tǒng)一多項式擬合階數(shù)，以求取各個微測井點在該統(tǒng)一 多項式擬合階數(shù)下的擬合系數(shù)；第二擬合系數(shù)獲取單元，根據(jù)所述各個微測井點在該統(tǒng)一 多項式擬合階數(shù)下的擬合系數(shù)，進行線性插值得到所述Μ個電磁測點中各個電磁測點的擬 合系數(shù)；速度獲取單元，將各個電磁測點的擬合系數(shù)、黃土層的反演厚度和反演電阻率，代 入所述擬合關(guān)系式，得到各個電磁測點的速度；靜校正單元，根據(jù)所述各個電磁測點的黃土 層厚度和各個電磁測點的速度，進行黃土層靜校正。
[0020] 在一個實施例中，所述速度與電阻率擬合關(guān)系式為：
[0022] 其中，Vlgl表示第i個微測井點在深度Dig處的速度，a1表示擬合系數(shù)，Rigi表示 第i個微測井點在深度Dig處的電阻率，ηι表示第i個微測井點的多項式擬合總階數(shù)，1^表 示微測井點的個數(shù)。
[0023] 在一個實施例中，第一擬合系數(shù)獲取單元包括：擬合速度獲取模塊，將各個微測井 點的黃土層的反演厚度作為微測井點的深度、反演電阻率作為在微測井點的深度處的電阻 率、不同階的多項式擬合系數(shù)作為擬合系數(shù)，代入所述速度與電阻率擬合關(guān)系式，得到各個 微測井點不同階的多項式擬合系數(shù)所對應的擬合速度；相關(guān)系數(shù)計算模塊，根據(jù)各個微測 井點不同階的多項式擬合系數(shù)所對應的擬合速度，計算各個微測井點不同多項式擬合階數(shù) 下的相關(guān)系數(shù)；最大多項式擬合階數(shù)獲取模塊，對每個微測井，選取最大相關(guān)系數(shù)所對應的 多項式擬合階數(shù)作為該微測井的最大多項式擬合階數(shù)，以得到N個最大多項式擬合階數(shù)； 統(tǒng)一多項式擬合階數(shù)獲取模塊，根據(jù)所述N個最大多項式擬合階數(shù)計算得到統(tǒng)一多項式擬 合階數(shù)：擬合系數(shù)獲取模塊，將所述統(tǒng)一多項式擬合系數(shù)作為所述N個微測井點的多項式 擬合階數(shù)，代入所述速度與電阻率擬合關(guān)系式，得到所述各個微測井點在該多項式擬合階 數(shù)下的擬合系數(shù)。
[0024] 在一個實施例中，所述相關(guān)系數(shù)計算模塊具體用于按照以下公式計算各個微測井 點不同多項式擬合階數(shù)下的相關(guān)系數(shù)：
[0026] 其中，Q為第i個微測井點的相關(guān)系數(shù)，Nk為每個微測井點的深度點數(shù)，U為第 i個微測井點第k個深度點擬合后的速度，歹為第i個微測井點所有深度點擬合后的速度 平均值，F(xiàn)/g,4為第i個微測井點第k個深度點的速度，.為第i個微測井點所有深度點的 速度平均值。
[0027] 在上述實施例中，在計算不同深度的黃土塬對應地層速度的過程中，在原有的電 磁測點的基礎上布設了微測井點和深微測井點，然后基于電磁測點的電磁數(shù)據(jù)和微測井點 的微測井數(shù)據(jù)，確定各個微測井點在統(tǒng)一多項式擬合階數(shù)下的擬合系數(shù)，然后基于這些微 測井點在統(tǒng)一多項式擬合階數(shù)下的