地震偏移成像的方法、裝置、計算機設備及存儲介質制造方法及圖紙

本申請提供了一種地震偏移成像的方法、裝置、計算機設備及存儲介質，屬于地球物理勘探技術領域。該方法包括：獲取待測工區內的油氣儲層對應的地震波數據；確定油氣儲層對應的正演模擬參數、共軛轉置偏移算子參數、敏感核函數參數；基于正演模擬算子參數、共軛轉置偏移算子參數和敏感核函數參數，確定油氣儲層對應的反射率向量；通過預設迭代次數的迭代計算，確定每次迭代的反射率向量對應的正演地震波數據與地震波數據之間的偏移參數，確定偏移參數最小時對應的反射率向量；確定與反射率向量對應的橫波疊加剖面和縱波疊加剖面。由于通過最小二乘法確定了與實際地震波數據最接近的正演數據，進而提高了確定的疊加剖面的準確性。性。性。

【技術實現步驟摘要】
地震偏移成像的方法、裝置、計算機設備及存儲介質


[0001]本申請涉及地球物理勘探
，特別涉及一種地震偏移成像的方法、裝置、計算機設備及存儲介質。

技術介紹

[0002]目前，隨著常規油氣儲層越來越少，深海油氣儲層、煤層氣儲層和頁巖氣儲層等非常規油氣儲層成為勘探重點。在對非常規油氣儲層進行勘探時，因其構造復雜，常規偏移成像方法無法滿足對非常規油氣儲層進行高精度成像的要求。
[0003]相關技術中，逆時偏移成像方法可以對非常規油氣儲層進行成像；逆時偏移成像方法是通過求解油氣儲層內縱波和橫波的雙程波動方程，來模擬縱波和橫波的傳播過程，然后應用成像條件得到彈性波成像剖面，進而實現對復雜油氣儲層進行成像。
[0004]但是，高精度逆時偏移成像方法在應用成像條件前，需要對縱波和橫波進行波場分離，這個過程容易產生相位畸變、縱波和橫波振幅不準確的問題，最終降低了成像精度。

技術實現思路

[0005]本申請實施例提供了一種地震偏移成像的方法、裝置、計算機設備及存儲介質，可以提高成像的準確性。所述技術方案如下：
[0006]一方面，本申請提供了一種地震偏移成像的方法，所述方法包括：
[0007]獲取待測工區內的油氣儲層對應的地震波數據；
[0008]確定所述油氣儲層對應的正演模擬參數，以及，確定所述油氣儲層對應的共軛轉置偏移算子參數，以及確定所述油氣儲層對應的敏感核函數參數；
[0009]基于所述正演模擬算子參數、所述共軛轉置偏移算子參數和所述敏感核函數參數，確定所述油氣儲層對應的反射率向量；
[0010]對所述反射率向量進行預設迭代次數的迭代計算，確定每次迭代的反射率向量對應的正演地震波數據與地震波數據之間的偏移參數，確定所述偏移參數最小時對應的反射率向量，所述偏移參數最小時對應的反射率向量包括轉換橫波反射率參數和縱波反射率參數；
[0011]確定與所述轉換橫波反射率參數對應的橫波疊加剖面以及與所述縱波反射率參數對應的縱波疊加剖面，所述橫波疊加剖面和所述縱波疊加剖面用于對所述油氣儲層進行偏移成像。
[0012]在一種可能的實現方式中，所述基于所述正演模擬算子參數、所述共軛轉置偏移算子參數和所述敏感核函數參數，確定所述油氣儲層對應的反射率向量，包括：
[0013]基于所述正演模擬算子參數，確定波恩算子參數；以及，基于所述正演模擬算子參數、所述共軛轉置偏移算子參數和所述敏感核函數參數，確定伴隨偏移算子參數；
[0014]確定初始模型反射率向量對應的初始正演地震波數據，根據所述波恩算子參數、所述伴隨偏移算子參數、所述初始正演地震波數據和所述油氣儲層對應的地震波數據，確
定所述初始正演地震波數據對應的初始變化梯度；
[0015]根據所述初始模型反射率向量和所述初始變化梯度，確定更新迭代的反射率向量，將更新迭代的反射率向量作為所述油氣儲層對應的反射率向量。
[0016]在另一種可能的實現方式中，所述確定每次迭代的反射率向量對應的正演地震波數據與地震波數據之間的偏移參數，包括：
[0017]獲取所述油氣儲層對應的調解參數；
[0018]根據所述波恩算子參數、所述調解參數、每次迭代的所述反射率向量對應的正演地震波數據與所述地震波數據，通過以下公式一，確定每次迭代的所述反射率向量對應的正演地震波數據與所述地震波數據之間的偏移參數；
[0019]公式一：J(m)＝||Lm
?
d
obs
||2+λ||m||
TV
[0020]其中，J(m)表示所述偏移參數，m表示所述反射率向量，d
obs
表示所述地震波數據，L表示所述波恩算子參數，λ表示所述調解參數。
[0021]在另一種可能的實現方式中，所述地震波數據包括橫波數據和縱波數據；
[0022]所述確定所述油氣儲層對應的正演模擬參數，包括：
[0023]確定所述橫波數據對應的橫波速度、所述縱波數據對應的縱波速度、所述油氣儲層的密度、表征垂向和橫向縱波速度差異的各向異性參數、表征垂向入射縱波速度兩次導的各向異性參數和理論震源數據，所述理論震源數據包括理論震源發射的震源水平分量和震源垂直分量；
[0024]根據所述橫波速度、所述縱波速度、所述密度、所述表征垂向和橫向縱波速度差異的各向異性參數和所述表征垂向入射縱波速度兩次導的各向異性參數，確定剛度矩陣參數；
[0025]根據所述密度、所述剛度矩陣分量參數所述理論震源數據，確定所述油氣儲層對應的水平速度、垂直速度、x方向的正應力、z方向的正應力和所述切應力，確定所述油氣儲層對應的正演模擬參數。
[0026]在另一種可能的實現方式中，所述確定所述油氣儲層對應的共軛轉置偏移算子參數，包括：
[0027]確定所述油氣儲層內的目標點的水平速度和垂直速度以及所述目標點所在位置的儲層的密度、表征垂向和橫向縱波速度差異的各向異性參數、表征垂向入射縱波速度兩次導的各向異性參數、x方向的正應力、z方向的正應力和切應力，所述x方向的正應力和所述z方向的正應力相互垂直；
[0028]確定所述水平速度對應的水平伴隨速度和所述垂直速度對應的垂直伴隨速度，以及，確定所述x方向的正應力對應的x方向的伴隨正應力、所述z方向的正應力對應的z方向的伴隨正應力和所述切應力對應的伴隨切應力；
[0029]所述根據所述水平伴隨速度、所述垂直伴隨速度、所述x方向的伴隨正應力、所述z方向的伴隨正應力和所述伴隨切應力，確定所述油氣儲層對應的共軛轉置偏移算子參數。
[0030]在另一種可能的實現方式中，所述地震波數據包括橫波數據和縱波數據；
[0031]所述根據所述水平伴隨速度、所述垂直伴隨速度、所述x方向的伴隨正應力、所述z方向的伴隨正應力和所述伴隨切應力，確定所述油氣儲層對應的共軛轉置偏移算子參數，包括：
[0032]根據所述水平伴隨速度、所述垂直伴隨速度、所述x方向的伴隨正應力、所述z方向的伴隨正應力和所述伴隨切應力，通過以下公式七，得到水平分量的數據殘差和垂直分量的數據殘差，確定所述水平分量的數據殘差和所述垂直分量的數據殘差為所述油氣儲層對應的共軛轉置偏移算子參數；
[0033]公式七：
[0034]其中，Δd
x
表示所述水平分量的數據殘差，Δd
z
表示所述垂直分量的數據殘差，C
33
表示所述第一剛度矩陣分量，C
55
表示所述第二剛度矩陣分量，C
11
表示所述第三剛度矩陣分量，C
13
表示所述第四剛度矩陣分量，表示所述水平伴隨速度，氛示所述垂直伴隨速度，表示所述x方向的伴隨正應力，表示所述z方向的伴隨正應力，表示所述伴隨切應力，ρ表示所述密度。
[0035]在另一種可能的實現方式中，所述確定所述油氣儲層對應的敏感核函數參數，包括：
[0036]根據所述目標點對應的表征垂向和橫向縱本文檔來自技高網...

【技術保護點】

【技術特征摘要】
1.一種地震偏移成像的方法，其特征在于，所述方法包括：獲取待測工區內的油氣儲層對應的地震波數據；確定所述油氣儲層對應的正演模擬參數，以及，確定所述油氣儲層對應的共軛轉置偏移算子參數，以及確定所述油氣儲層對應的敏感核函數參數；基于所述正演模擬算子參數、所述共軛轉置偏移算子參數和所述敏感核函數參數，確定所述油氣儲層對應的反射率向量；對所述反射率向量進行預設迭代次數的迭代計算，確定每次迭代的反射率向量對應的正演地震波數據與地震波數據之間的偏移參數，確定所述偏移參數最小時對應的反射率向量，所述偏移參數最小時對應的反射率向量包括轉換橫波反射率參數和縱波反射率參數；確定與所述轉換橫波反射率參數對應的橫波疊加剖面以及與所述縱波反射率參數對應的縱波疊加剖面，所述橫波疊加剖面和所述縱波疊加剖面用于對所述油氣儲層進行偏移成像。2.根據權利要求1所述的方法，其特征在于，所述基于所述正演模擬算子參數、所述共軛轉置偏移算子參數和所述敏感核函數參數，確定所述油氣儲層對應的反射率向量，包括：基于所述正演模擬算子參數，確定波恩算子參數；以及，基于所述正演模擬算子參數、所述共軛轉置偏移算子參數和所述敏感核函數參數，確定伴隨偏移算子參數；確定初始模型反射率向量對應的初始正演地震波數據，根據所述波恩算子參數、所述伴隨偏移算子參數、所述初始正演地震波數據和所述油氣儲層對應的地震波數據，確定所述初始正演地震波數據對應的初始變化梯度；根據所述初始模型反射率向量和所述初始變化梯度，確定更新迭代的反射率向量，將所述更新迭代的反射率向量作為所述油氣儲層對應的反射率向量。3.根據權利要求1所述的方法，其特征在于，所述確定每次迭代的反射率向量對應的正演地震波數據與地震波數據之間的偏移參數，包括：獲取所述油氣儲層對應的調解參數；根據所述波恩算子參數、所述調解參數、每次迭代的所述反射率向量對應的正演地震波數據與所述地震波數據，通過以下公式一，確定所述反射率向量對應的正演地震波數據與所述地震波數據之間的偏移參數；公式一：J(m)＝‖Lm
?
d
obs
‖2+λ‖m‖
TV
其中，J(m)表示所述偏移參數，m表示所述反射率向量，d
obs
表示所述地震波數據，L表示所述波恩算子參數，λ表示所述調解參數。4.根據權利要求1所述的方法，其特征在于，所述地震波數據包括橫波數據和縱波數據；所述確定所述油氣儲層對應的正演模擬參數，包括：確定所述橫波數據對應的橫波速度、所述縱波數據對應的縱波速度、所述油氣儲層的密度、表征垂向和橫向縱波速度差異的各向異性參數、表征垂向入射縱波速度兩次導的各向異性參數和理論震源數據，所述理論震源數據包括理論震源發射的震源水平分量和震源垂直分量；根據所述橫波速度、所述縱波速度、所述密度、所述表征垂向和橫向縱波速度差異的各向異性參數和所述表征垂向入射縱波速度兩次導的各向異性參數，確定剛度矩陣參數；
根據所述密度、所述剛度矩陣分量參數所述理論震源數據，確定所述油氣儲層對應的水平速度、垂直速度、x方向的正應力、z方向的正應力和所述切應力，確定所述油氣儲層對應的正演模擬參數。5.根據權利要求1所述的方法，其特征在于，所述確定所述油氣儲層對應的共軛轉置偏移算子參數，包括：確定所述油氣儲層內的目標點的水平速度和垂直速度以及所述目標點所在位置的儲層的密度、表征垂向和橫向縱波速度差異的各向異性參數、表征垂向入射縱波速度兩次導的各向異性參數、x方向的正應力、z方向的正應力和切應力，所述x方向的正應力和所述z方向的正應力相互垂直；確定所述水平速度對應的水平伴隨速度和所述垂直速度對應的垂直伴隨速度，以及，確定所述x方向的正應力對應的x方向的伴隨正應力、所述z方向的正應力對應的z方向的伴隨正應力和所述切應力對應的伴隨切應力；所述根據所述水平伴隨速度...

【專利技術屬性】
技術研發人員：劉文平，趙楊，左佳卉，雍銳，楊學鋒，吳建發，文山師，石學文，楊繼東，王星皓，
申請(專利權)人：中國石油天然氣股份有限公司，
類型：發明
國別省市：