[發明專利]一種利用時移地震勘探資料進行差異AVO分析的方法

有效

• 申請號: 201110320920.8
• 申請日: 2011-10-20
• 公開（公告）號: CN102508294A
• 公開（公告）日: 2012-06-20
• 發明（設計）人: 馬勁風
• 申請（專利權）人: 西北大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 西安通大專利代理有限責任公司 61200
• 代理人: 徐文權

1.一種利用時移地震勘探資料進行差異AVO分析的方法，其特征在于：包括如下步驟：

1)根據Zoeppritz方程的近似彈性模量簡化公式推導兩次地震的反射系數差異；Zoeppritz方程的近似彈性模量簡化公式為：

R(p)≈Rf(p)-2Δμρp2=Rf(p)-2Δμρβ2sin2φ---(1)]]>

其中，R_f(p)為流體反射系數，μ為剪切模量，p為地震波傳播的射線參數，α為縱波速度，β為橫波速度，ρ為巖層密度，為P-SV轉換波的轉換角，Δμ為上下巖層切變模量差的平均；

2)由此，計算兩次地震的反射系數差異

ΔRf′(p)=R2′(p)-R1(p)≈Rf2′(p)-Rf1(p)]]>

≈12ln(ρ2′α2′1-(α2′p)2/ρ1α11-(α1p)2)-12ln(ρ2α21-(α2p)2/ρ1α11-(α1p)2)---(2)]]>

=12ln(ρ2′α2′1-(α2′p)2/ρ2α21-(α2p)2)=ρ2′α2′/1-(α2′p)2-ρ2α2/1-(α2p)2ρ2′α2′/1-(α2′p)2+ρ2α2/1-(α2p)2]]>

其中，R’₂(p)和R₁(p)分別為四維地震監測后、前的反射系數；R′_f2(p)和R_f1(p))分別為四維地震監測后、前的流體反射系數；α₁，β₁，ρ₁為蓋層的縱波速度、橫波速度和密度；α₂，β₂，ρ₂為儲層注入流體前或者開發前的縱波速度、橫波速度和密度；α’₂，β’₂，ρ’₂為儲層注入流體后或者開發后的縱波速度、橫波速度和密度；

這個公式也變換為：

ρ2′α2′1-(α2′p)2=ρ2α21-(α2p)21+ΔRf′(p)1-ΔRf′(p)=ρ2α2cosθ21+ΔRf′(p)1-ΔRf′(p)---(3)]]>

其中為Brekhovskikh(1960)定義的聲波阻抗；和ΔRf’(p)都用來分析時移地震觀測中流體的變化；θ₂為四維地震監測前地震波以θ₁入射到儲層時的透射角；

3)差異AVO的正演計算：

結合流體替換模型，計算流體替換或者注入前后彈性參數的變化，再根據上述公式計算差異AVO的變化。

2.根據權利要求1所述的方法，其特征在于，該方法可以直接用于差異AVO的反演計算，利用公式(3)將差異AVO或者反射系數轉換為Brekhovskikh定義的聲波阻抗ρα/cosθ，將聲波阻抗轉化為流體飽和度；或者直接建立流體飽和度變化與公式(2)的關系，從差異反射系數中直接反演流體飽和度。