[發明專利]一種基于多GPU計算的三維TTI介質逆時偏移成像方法

有效

• 申請號: 201610059900.2
• 申請日: 2016-01-28
• 公開（公告）號: CN105717539B
• 公開（公告）日: 2018-01-30
• 發明（設計）人: 劉國峰;李春
• 申請（專利權）人: 中國地質大學（北京）
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京中譽威圣知識產權代理有限公司11279
• 代理人: 蔣常雪

1.一種基于多GPU計算的三維TTI介質逆時偏移成像方法，其特征在于：所述方法采用穩定的二階耦合方程組作為計算模型，給定有限SV波分量以增強波場傳播的穩定性，采用各向異性參數匹配法壓制偽橫波干擾，引入隨機速度邊界及ε,δ,θ和四個參數，對這四個模型參數采用對應參數的內邊界進行外擴，均勻各向異性邊界，避免大量的波場數據存儲和頻繁的磁盤讀寫操作，采用基于GPU的CUDA平臺實現高性能并行計算；

所述二階耦合方程組的形式是：

∂2p∂t2=(1+2ϵ)Vpz2H2p+Vpz2H1q+Vsz2H1(p-q)(a)∂2q∂t2=(1+2δ)Vpz2H2p+Vpz2H1q+Vsz2H2(q-p)(b)]]>

H1=sin2θcos2φ∂2x2+sin2θsin2φ∂2y2+cos2θ∂2∂z2+sin2θsin2φ∂2∂x∂y+sin2θsinφ∂2∂y∂z+sin2θcosφ∂2∂x∂z]]>

H2=∂2∂x2+∂2∂y2+∂2∂z2-H1]]>

其中，p是主波場，q是輔助波場；ε和δ是Thomson各向異性參數；θ和分別是對稱軸的傾角和方位角；H₁和H₂是導數算子；V_pz是對稱軸方向縱波速度，V_sz是對稱軸方向橫波速度；V_sz由參數σ和V_pz來控制；

所述該各向異性參數匹配法的過程是：首先確定震源位置，在震源位置周圍設置一個不小于30個采樣點厚度的薄層；然后在該薄層內設置各向異性參數δ等于ε，即改變δ的值，令其隨ε變化；再次，在該薄層的突變帶增加一個大約10個采樣點厚度的薄層的過渡帶，讓δ平滑緩慢的過渡到其真實值，用匹配后的δ參數模型作為該各向異性參數δ的輸入。