[發(fā)明專利]適用于TTI介質(zhì)的頻率域有限元全吸收PML方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201610565960.1
• 申請(qǐng)日: 2016-07-18
• 公開（公告）號(hào): CN106094038B
• 公開（公告）日: 2017-11-14
• 發(fā)明（設(shè)計(jì)）人: 王兵;張闊;馬明明
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油大學(xué)(北京)
• 主分類號(hào): G01V1/50
• 分類號(hào): G01V1/50
• 代理公司: 北京三友知識(shí)產(chǎn)權(quán)代理有限公司11127
• 代理人: 李強(qiáng)

1.適用于TTI介質(zhì)的頻率域有限元全吸收PML方法，用于構(gòu)造頻域有限元全吸收PML邊界條件，實(shí)現(xiàn)TTI介質(zhì)復(fù)雜井眼條件的高精度聲波測(cè)井模擬；其特征在于：

所述方法采用復(fù)拉伸坐標(biāo)系的彈性波方程，所述復(fù)拉伸坐標(biāo)系的彈性波方程在所述頻域下表示形式為：

-ρω2u~=∂τ~xx∂x+∂τ~xz∂z-ρω2w~=∂τ~zx∂x+∂τ~zz∂z]]>

在PML區(qū)域內(nèi)，加入拉伸函數(shù)，將上述方程表示為：

-ρω2u~=1sx∂τ~xx∂x+1sz∂τ~xz∂z-ρω2w~=1sx∂τ~zx∂x+1sz∂τ~zz∂z]]>

所述PML是完全匹配層，s_x和s_z為PML吸收邊界內(nèi)的拉伸函數(shù)，s_x和s_z分別構(gòu)造x和z方向的兩個(gè)拉伸函數(shù)分量，s_x和s_z的表達(dá)式如下：

sx=kx+dx+mdzαx+iω]]>

sz=kz+dz+mdxαz+iω]]>

上述式中：

為x方向的實(shí)數(shù)衰減因子，

為x方向復(fù)頻移因子，

為x方向倏逝波抑制因子，

L_x為吸收層x方向的厚度；

m為切向衰減因子；

為z方向的實(shí)數(shù)衰減因子，

為z方向復(fù)頻移因子，

為z方向倏逝波抑制因子，

L_z為吸收層z方向的厚度，

ω為角頻率，為頻率的2π倍，

i為虛數(shù)單位。

2.根據(jù)權(quán)利要求1所述的適用于TTI介質(zhì)的頻率域有限元全吸收PML方法，其特征在于，所述實(shí)數(shù)衰減因子復(fù)頻移因子α₀＝2ω，倏逝波抑制因子中k₀取值范圍為1～20，切向衰減因子m取值范圍為0.005～0.02。