[發明專利]地層傾角的確定方法和裝置

有效

• 申請號: 201710740537.5
• 申請日: 2017-08-25
• 公開（公告）號: CN107703541B
• 公開（公告）日: 2019-08-06
• 發明（設計）人: 李海山;楊午陽;高建虎;周春雷;杜炳毅;鄢高韓
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/36
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 李輝;徐煥

1.一種地層傾角的確定方法，其特征在于，包括：

計算待測地震數據在多個預設方向中各個預設方向的總能量；

確定多個總能量中的最大值，并確定所述最大值對應的方向；

計算所述最大值對應方向的地層傾角，將確定的地層傾角作為所述待測地震數據相應的地層傾角；

其中，計算所述待測地震數據在多個預設方向中各個預設方向的總能量，包括：

根據所述待測地震數據的曲波變換系數，計算得到所述待測地震數據在多個預設方向在多個預設尺度下的實部分量以及虛部分量；

計算所述多個預設方向中各個預設方向在多個預設尺度下的實部分量之和，得到多個預設方向中各個預設方向的實部和；

計算所述多個預設方向中各個預設方向在多個預設尺度下的虛部分量之和，得到多個預設方向中各個預設方向的虛部和；

根據所述實部分量、所述虛部分量、所述實部和以及所述虛部和，計算得到所述待測地震數據在多個預設方向中各個預設方向的總能量。

2.如權利要求1所述的方法，其特征在于，確定多個總能量中的最大值，包括：

對所述多個預設方向中各個預設方向的總能量進行去噪處理；

根據去噪處理后的各個預設方向的總能量，確定多個總能量中的最大值。

3.如權利要求2所述的方法，其特征在于，按照以下公式對所述多個預設方向中各個預設方向的總能量進行去噪處理：

其中，ThreshTotalE(io,n)表示所述待測地震數據在預設方向為io、下標為n時進行去噪處理后的總能量，1≤io≤No，No表示所述預設方向的總數，n＝(n₁,n₂)，TotalE(io,n)表示所述待測地震數據在預設方向為io、下標為n時的總能量，T表示預設的去噪閾值。

4.如權利要求3所述的方法，其特征在于，所述預設的去噪閾值為TotalE(io,n)×10^-6。

5.如權利要求1所述的方法，其特征在于，根據所述實部分量、所述虛部分量、所述實部和以及所述虛部和，計算得到所述待測地震數據在多個預設方向中各個預設方向的總能量，包括：

對所述實部和進行歸一化處理，得到歸一化后的實部和；

對所述虛部和進行歸一化處理，得到歸一化后的虛部和；

根據所述實部分量、所述虛部分量、所述歸一化后的實部和以及所述歸一化后的虛部和，計算得到所述待測地震數據在多個預設方向中各個預設方向的總能量。

6.如權利要求1所述的方法，其特征在于，按照以下公式計算得到所述待測地震數據在多個預設方向在多個預設尺度下的實部分量以及虛部分量：

其中，RS(is,io,n)表示所述待測地震數據在預設尺度為is、預設方向為io、下標為n時的實部分量，IS(is,io,n)表示所述待測地震數據在預設尺度為is、預設方向為io、下標為n時的虛部分量，n＝(n₁,n₂)，0≤n₁≤Nt,0≤n₂≤Nx表示所述待測地震數據的下標，Nt表示所述待測地震數據在時間方向的采樣點數，Nx表示所述待測地震數據在空間方向的采樣點數，Coef(is,io,k)表示所述曲波變換系數，is表示所述預設尺度，1≤is≤Ns，io表示所述預設方向，1≤io≤No，k＝(k₁,k₂)表示所述曲波變換系數的下標，表示ψ_is,io,k(n₁,n₂)的共軛復數，ψ_is,io,k(n₁,n₂)表示在預設尺度為is、預設方向為io、下標為k時的離散逆曲波，RCoef(is,io,k)表示置Coef(is,io,k)的虛部為零，ICoef(is,io,k)表示置Coef(is,io,k)的實部為零。