基于擴展彈性阻抗的地質儲層檢測方法
【技術領域】
[0001] 本發(fā)明總體來說涉及油氣田勘探技術領域，屬于地震資料解釋范疇，更具體地講， 涉及一種基于擴展彈性阻抗的地質儲層檢測方法。
【背景技術】
[0002] 地球物理勘探利用常規(guī)的縱波疊后偏移地震資料處理中損失的信息來解決構造 油氣藏問題不大，但是要解決巖性復雜的油氣藏就顯得力不從心。隨著疊前偏移技術的發(fā) 展，疊前道集中蘊藏的不同方位角、入射角地震信息可用來解決裂縫、巖性和流體識別技術 問題。
[0003] 彈性阻抗是聲阻抗的概念延伸和推廣，它建立在非零偏移距的基礎上，是縱波速 度、橫波速度、密度以及入射角的函數。彈性波阻抗反演屬于疊前反演技術，它包含了豐富 的巖性以及流體信息，但是由于AVO效應的影響，常規(guī)彈性阻抗值隨入射角發(fā)生劇烈變化 (即，彈性阻抗值隨入射角的增大而急劇減?。?，可能會掩蓋遠、近角彈性阻抗值分析的某 些信息（例如，包含流體或巖性變化差異的信息）。

【發(fā)明內容】

[0004] 本發(fā)明的目的是提出一種基于擴展彈性阻抗的地質儲層檢測方法，以解決彈性阻 抗值隨入射角的增大而急劇減小，無法有效識別地下地質儲層的巖性和流體性質的技術問 題。
[0005] 本發(fā)明的一方面提供一種基于擴展彈性阻抗的地質儲層檢測方法，所述方法包 括：(a)獲取目的層段的測井數據，在所述目的層段中選取多個采樣點，并從所述測井數據 中提取與所述多個采樣點中的每個采樣點對應的縱波速度、橫波速度和密度；(b)選取多 個實際入射角，并根據提取的縱波速度、橫波速度和密度，來計算所述多個實際入射角中的 每個實際入射角在所述每個采樣點處的反射系數；（c)選取多個理論入射角，并根據提取 的縱波速度、橫波速度和密度，來計算所述多個理論入射角中的每個理論入射角在所述每 個采樣點處的反射系數；（d)根據所述多個理論入射角的反射系數和所述多個實際入射角 的反射系數，分別確定與所述多個實際入射角中的每個實際入射角對應的等效理論入射 角；（e)利用基于所述等效理論入射角的擴展彈性阻抗來進行巖性和流體檢測，以識別地 下地質儲層的狀態(tài)。
[0006] 可選地，在步驟（b)中，選取多個實際入射角的步驟可包括：獲取實際入射角的初 始值，并將獲取的實際入射角的初始值作為最大實際入射角，在零度至最大實際入射角的 范圍之內以第一增量步長為間隔，選取所述多個實際入射角。
[0007] 可選地，在步驟（C)中，選取所述多個理論入射角的步驟可包括：在預定角度范圍 內以第二增量步長為間隔，選取所述多個理論入射角。
[0008] 可選地，第二增量步長的數值可小于第一增量步長的數值。
[0009] 可選地，在步驟（C)中，根據提取的縱波速度、橫波速度和密度，來計算所述多個 理論入射角中的每個理論入射角在所述每個采樣點處的反射系數的步驟可包括：（Cl)根 據提取的縱波速度、橫波速度和密度，計算所述多個理論入射角中的所述每個理論入射角 在所述每個采樣點處的擴展彈性阻抗；(C2)根據所述每個理論入射角在所述每個采樣點 處的擴展彈性阻抗，來計算所述每個理論入射角在所述每個采樣點處的反射系數。
[0010] 可選地，在步驟（Cl)中，可利用下面的公式的來計算所述每個理論入射角在所述 每個采樣點處的擴展彈性阻抗，
[0012] 其中，
[0013] Pj= cos X j+sin X j
[0014] qj= _8Ksin x j
[0015] rj= cos x j-4Ksin x s
[0016] 其中，EEIk ( x ,)為所述多個理論入射角中的第j個理論入射角x ,在所述多個采 樣點中的第k個采樣點處的擴展彈性阻抗，I < j <M，M為選取的理論入射角的總個數， I < k < N，N為選取的采樣點的總個數，a k為與第k個采樣點對應的縱波速度，β k為與第 k個采樣點對應的橫波速度，Pk為與第k個采樣點對應的密度，α。為所述目的層段中的縱 波速度的初始值，為所述目的層段中的橫波速度的初始值，P ^為所述目的層段中的密 度的初始值，K為縱橫波速度比的平均值，K= [(0k+1/ak+1)2+(0k/a k)2]/2，ak+1為與第 k+1個采樣點對應的縱波速度，β k+1為與第k+1個采樣點對應的橫波速度。
[0017] 可選地，在步驟（c2)中，可利用下面的公式的來計算所述每個理論入射角在所述 每個采樣點處的反射系數，
[0019] 其中，Rk ( X ,)為所述多個理論入射角中的第j個理論入射角X ,在所述多個采樣 點中的第k個采樣點處的反射系數，I < j <M，M為選取的理論入射角的總個數，I < N， N為選取的采樣點的總個數，EEIk+1( X ,)為所述多個理論入射角中的第j個理論入射角X , 在第k+Ι個采樣點處的擴展彈性阻抗，EEIk ( X ,)為所述多個理論入射角中的第j個理論入 射角 '在第k個采樣點處的擴展彈性阻抗。
[0020] 可選地，在步驟（b)中，根據提取的縱波速度、橫波速度和密度，來計算所述多個 實際入射角中的每個實際入射角在所述每個采樣點處的反射系數的步驟可包括：(bl)根 據提取的縱波速度、橫波速度和密度，計算所述多個實際入射角中的每個實際入射角在所 述每個采樣點處的彈性阻抗；(b2)根據所述每個實際入射角在所述每個采樣點處的彈性 阻抗，來計算所述每個實際入射角在所述每個采樣點處的反射系數。
[0021] 可選地，在步驟（d)中，根據所述多個理論入射角的反射系數和所述多個實際入 射角的反射系數，分別確定與所述多個實際入射角中的每個實際入射角對應的等效理論入 射角的步驟可包括：（dl)針對所述多個實際入射角中的任一實際入射角，分別計算所述任 一實際入射角的反射系數與所述多個理論入射角中的每個理論入射角的反射系數的相關 系數；(d2)確定所述相關系數中的最大值，并將所述相關系數中的最大值所對應的理論入 射角作為所述任一實際入射角對應的等效理論入射角。
[0022] 可選地，在步驟（dl)中，可利用下面的公式的來計算所述任一實際入射角的反射 系數與所述多個理論入射角中的每個理論入射角的反射系數的相關系數，
[0024] 其中，Corr[R( Θ丄R( X )]為所述多個實際入射角中的第i個實際入射角Θ ^勺 反射系數R(S1)與所述多個理論入射角中的第j個理論入射角X,的反射系數R(x ,)的 相關系數，I < i < η，η為選取的實際入射角的總個數，I < j < M，M為選取的理論入射角 的總個數，I < k < N，N為選取的采樣點的總個數。
[0025] 采用本發(fā)明示例性實施例的上述基于擴展彈性阻抗的地質儲層檢測方法，可準確 地預測地下地質儲層的巖性和流體性質，以真實反映地下地質儲層的真實狀態(tài)，有效降低 地質勘探風險，提尚鉆探成功率。
【附圖說明】
[0026] 通過下面結合附圖進行的詳細描述，本發(fā)明的上述和其它目的、特點和優(yōu)點將會 變得更加清楚，其中：
[0027] 圖1示出根據本發(fā)明示例性實施例的基于擴展彈性阻抗的地質儲層檢測方法的 流程圖；
[0028] 圖2示出根據本發(fā)明示例性實施例的計算每個實際入射角在每個采樣點處的反 射系數的步驟的流程圖；
[0029] 圖3示出根據本發(fā)明示例性實施例的計算每個理論入射角在每個采樣點處的反 射系數的步驟的流程圖；
[0030] 圖4示出根據本發(fā)明示例性實施例的確定與每個實際入射角對應的等效理論入 射角的步驟的流程圖；
[0031] 圖5示出根據本發(fā)明示例性實施例的彈性阻抗與擴展彈性阻抗分別隨入射角變 化的關系圖；
[0032] 圖6示出根據本發(fā)明示例性實施例的實際入射角與理論入射角的等效換算圖。
【具體實施方式】
[0033] 下面，將參照附圖詳細描述本發(fā)明的示例性實施例。
[0034] 提供參照附圖的以下描述以幫助對由權利要求及其等同物限定的本發(fā)明的示例 性實施例的全面理解。包括各種特定細節(jié)以幫助理解，但這些細節(jié)僅被視為是示例性的。因 此，本領域的普通技術人員將認識到在不脫離本發(fā)明的范圍和精神的情況下，可對描述于 此的示例性實施例進行各種改變和修改。此外，為了清楚和簡潔，省略對公知的功能和結構 的描述。
[0035] 圖1示出根據本發(fā)明示例性實施例的基于擴展彈性阻抗的地質儲層檢測方法的 流程圖。
[0036] 在步驟Sl