[發明專利]一種深度域合成地震記錄的制作方法

有效

• 申請號: 202110296191.0
• 申請日: 2021-03-19
• 公開（公告）號: CN113031070B
• 公開（公告）日: 2022-06-28
• 發明（設計）人: 陳樹民;陳可洋;林鐵峰;楊微;姜傳金;喬衛
• 申請（專利權）人: 大慶油田有限責任公司;中國石油天然氣股份有限公司
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/28
• 代理公司: 大慶知文知識產權代理有限公司 23115
• 代理人: 荊曉紅

1.一種深度域合成地震記錄的制作方法，包括以下步驟：

(1)獲得研究區過某井位置的深度域地面地震成像數據道seisdata(z)及對應位置的深度域成像速度模型道vseis(z)、測量的聲波和密度測井曲線、VSP資料提取的深度域層速度曲線；z代表深度；

(2)完成深度域聲波和密度測井曲線的去異常值等常規標準化處理，這里測井曲線共含有n+1個測量點,深度采樣步長為dz；并將聲波和密度測井曲線每一個深度樣點換算得到聲波速度v_i和密度ρ_i，其中聲波速度v_i即為聲波時差的倒數v_i＝1/dt_i，利用公式Z_i＝ρ_i×v_i得到波阻抗曲線，(i＝1，2……，n)，依據具體步驟獲得考慮地層透射能量損失的深度域反射系數曲線；

(3)對步驟(2)標準化后的聲波測井曲線作換算和平滑處理，得到深度域層速度VEL或根據同井測量的VSP資料提取的深度域層速度VEL，并按照步驟(2)中的測井資料深度采樣步長dz對深度域層速度VEL進行重采樣，得到vvel(z)；

(4)利用步驟(1)中的井旁深度域地面地震成像數據道seisdata(z)及對應位置的深度域成像速度模型道vseis(z)提取偽深度域地震子波；

(5)利用步驟(2)中得到的考慮地層透射能量損失的深度域反射系數曲線R'(z)和步驟(3)中的速度vvel(z)換算為偽深度域反射系數曲線dzwell*(z)；

(6)將步驟(4)中得到的偽深度域地震子波seiswavelet(z)和步驟(5)中的偽深度域反射系數曲線R*(z)在偽深度域進行褶積處理，得到偽深度域合成地震記錄syn*(z)，即syn*(z)＝seiswavelet(z)*R*(z)；

(7)利用步驟(6)中偽深度域合成地震記錄syn*(z)和步驟(3)中的速度vvel(z)換算為深度域合成地震記錄syn(z)；

(8)完成深度域合成地震記錄的制作與深度域地震成像剖面的井震標定工作。

2.根據權利要求1所述的深度域合成地震記錄的制作方法，其特征在于：所述的步驟2獲得考慮地層透射能量損失的深度域反射系數曲線依據以下步驟:

1)利用公式R_i＝(Z_i+1-Z_i)/(Z_i+1+Z_i)，(i＝1,2……,n)將深度域波阻抗曲線Z(z)換算得到深度域反射系數曲線R(z)，z代表深度；

2)利用以下公式將深度域反射系數曲線R(z)換算得到考慮地層透射能量損失的反射系數曲線R'(z)：

R'₁＝R₁；

R'₂＝R₂(1-R²₁)；

R'₃＝R₃(1-R²₂)(1-R²₁)；

……

R'_n＝R_n(1-R²_n-1)(1-R²_n-2)……(1-R²₁)。

3.根據權利要求1所述的深度域合成地震記錄的制作方法，其特征在于：所述的步驟4提取偽深度域地震子波具體步驟如下:

1)按照計算公式dzseis*(z)＝dzseis×Vm/vseis(z)將井旁深度域地面地震成像道seisdata(z)轉換為偽深度域地面地震成像道seisdata*(z)；

其中：dzseis*為隨深度z變化的偽深度域采樣步長；dzseis為深度域常數采樣步長；Vm為偽深度域常數層速度；

2)對偽深度域地面地震成像道seisdata*(z)進行重采樣，使得深變的偽深度域采樣步長dzseis*(z)變為常數dzn；

3)采用常規時間域地震子波提取方法獲得偽深度域地震子波seiswavelet(z)。