地震学报
地震學報
지진학보
ACTA SEISMOLOGICA SINICA
2015年
4期
648-660
,共13页
杨富森%李振春%王小丹%张凯%刘成斋
楊富森%李振春%王小丹%張凱%劉成齋
양부삼%리진춘%왕소단%장개%류성재
各向异性介质%一阶拟声波%PML 边界条件%交错网格%正演模拟%qSV 人为干扰波
各嚮異性介質%一階擬聲波%PML 邊界條件%交錯網格%正縯模擬%qSV 人為榦擾波
각향이성개질%일계의성파%PML 변계조건%교착망격%정연모의%qSV 인위간우파
anisotropic media%first-order pseudo-acoustic wave%PML boundary condition%staggered-grid%modeling%qSV-wave artifact
为克服各向异性弹性波动方程正演模拟的局限,本文研究了各向异性介质拟声波方程的交错网格有限差分数值解法。首先,从 VTI 介质胡克定律和 qP-qSV 波频散关系两种思路出发,通过声假设近似,给出了两种不同形式的 VTI 介质一阶拟声波方程,并通过引入波场的伪速度分量,推导了一种新的 VTI 介质一阶应力-速度方程,并通过旋转坐标系将其推广到 TTI 介质中;其次,构造了一阶拟声波方程的交错网格高阶有限差分格式,并推导了相应的 PML 边界条件;最后,对本文方法中固有的 qSV 人为干扰波的产生机制和压制方法进行了简单讨论。数值结果表明:3种一阶拟声波方程在运动学和动力学上是等价的,相对于各向异性弹性波正演模拟,其节省了内存,提高了计算效率;各向异性因素会影响反射波旅行时和振幅等波场特征,在后续的处理、反演和解释中不可忽略;VTI 介质 HESS 模型的逆时偏移结果也验证了本文方法的合理性。
為剋服各嚮異性彈性波動方程正縯模擬的跼限,本文研究瞭各嚮異性介質擬聲波方程的交錯網格有限差分數值解法。首先,從 VTI 介質鬍剋定律和 qP-qSV 波頻散關繫兩種思路齣髮,通過聲假設近似,給齣瞭兩種不同形式的 VTI 介質一階擬聲波方程,併通過引入波場的偽速度分量,推導瞭一種新的 VTI 介質一階應力-速度方程,併通過鏇轉坐標繫將其推廣到 TTI 介質中;其次,構造瞭一階擬聲波方程的交錯網格高階有限差分格式,併推導瞭相應的 PML 邊界條件;最後,對本文方法中固有的 qSV 人為榦擾波的產生機製和壓製方法進行瞭簡單討論。數值結果錶明:3種一階擬聲波方程在運動學和動力學上是等價的,相對于各嚮異性彈性波正縯模擬,其節省瞭內存,提高瞭計算效率;各嚮異性因素會影響反射波旅行時和振幅等波場特徵,在後續的處理、反縯和解釋中不可忽略;VTI 介質 HESS 模型的逆時偏移結果也驗證瞭本文方法的閤理性。
위극복각향이성탄성파동방정정연모의적국한,본문연구료각향이성개질의성파방정적교착망격유한차분수치해법。수선,종 VTI 개질호극정률화 qP-qSV 파빈산관계량충사로출발,통과성가설근사,급출료량충불동형식적 VTI 개질일계의성파방정,병통과인입파장적위속도분량,추도료일충신적 VTI 개질일계응력-속도방정,병통과선전좌표계장기추엄도 TTI 개질중;기차,구조료일계의성파방정적교착망격고계유한차분격식,병추도료상응적 PML 변계조건;최후,대본문방법중고유적 qSV 인위간우파적산생궤제화압제방법진행료간단토론。수치결과표명:3충일계의성파방정재운동학화동역학상시등개적,상대우각향이성탄성파정연모의,기절성료내존,제고료계산효솔;각향이성인소회영향반사파여행시화진폭등파장특정,재후속적처리、반연화해석중불가홀략;VTI 개질 HESS 모형적역시편이결과야험증료본문방법적합이성。
In order to overcome the limitations of anisotropic elastic wave equa-tions modeling,this paper studied anisotropic pseudo-acoustic wave modeling by using staggered-grid finite-difference method.Starting from two different ideas, i.e .,Hooke’s law and qP-qSV dispersion relation for VTI media,two first-order pseudo-acoustic wave equations in two forms for VTI media are given.