[发明专利]一种TTI介质中纯qP波地震数据模拟方法

摘要

本说明书实施例公开了一种TTI介质中纯qP波地震数据模拟方法。本说明书实施例所提供的方案，通过将TTI介质中耦合在一起的qP波与qSV波解耦，相比于各向异性全弹性波波动方程正演模拟，解耦的纯qP波的模拟方法天然地解决了波场分离困难的问题，实际应用也证明了纯qP波模拟很好地应用于纵波勘探；模拟过程总不存在伪横波干扰且对于各向异性参数εδ仍稳定的优点；在各向异性参数剧烈变化时，模拟过程中波场仍可稳定传播，稳定性更好。

主权项

1.一种TTI介质中纯qP波地震数据模拟方法，包括：获取初始地震地质各向异性参数场，各向异性参数包括ε、δ、φ和v_p，其中，ε和δ为介质的Thomsen参数，φ为各向异性倾角，v_p表示介质纵波速度；将所述参数场划分为多个网格；针对每个网格，获取与所述各向异性参数相关的模拟方程系数λ₁,λ₂,λ₃；采用预设的波场传播算子，基于所述与该网格位置相关的模拟方程系数和所述各向异性参数场，确定炮点和检波点位置后，进行数值模拟，得到检波器接收到的应力场值，其中，所述波场传播算子包括：c₁＝(1+λ₃+λ₁cos⁴φ+2εcos²φ+λ₂cos²φ)，c₂＝(1+λ₃+λ₁sin⁴φ+2εsin²φ+λ₂sin²φ)，c₃＝(2+2λ₃+λ₂+2ε+1.5λ₁sin²2φ)，c₄＝(‑2εsin2φ‑λ₂sin2φ‑2λ₁sin2φcos²φ)，c₅＝(‑2εsin2φ‑λ₂sin2φ‑2λ₁sin2φsin²φ)；其中，u_p表示应力场值，u_q表示辅助应力场值，φ是各向异性对称轴与垂直方向夹角，t表示传播时间，i、j分别表示纵横网格点位置，k表示离散的时间值，△x表示横向网格间距，△z表示纵向网格间距，△t表示时间采样间隔，a₀和a_nn表示差分系数，FFT表示快速傅里叶变换，FFT^‑1表示快速傅里叶变换的反变换，k_x和k_z分别表示纵横向的波数值，N表示差分阶数，x和z分别表示横坐标和纵坐标。