亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        矩形旋轉(zhuǎn)交錯網(wǎng)格PML邊界條件的研究與實現(xiàn)

        2017-03-01 10:56:24孫盡堯
        電子設(shè)計工程 2017年2期
        關(guān)鍵詞:模型

        張 亮,孫盡堯

        (武漢大學(xué) 物理科學(xué)與技術(shù)學(xué)院,湖北 武漢430072)

        矩形旋轉(zhuǎn)交錯網(wǎng)格PML邊界條件的研究與實現(xiàn)

        張 亮,孫盡堯

        (武漢大學(xué) 物理科學(xué)與技術(shù)學(xué)院,湖北 武漢430072)

        現(xiàn)有的頻率域波全形反演多是默認(rèn)基于正方形旋轉(zhuǎn)交錯網(wǎng)格,而基于矩形旋轉(zhuǎn)交錯網(wǎng)格的全波形反演有更大的靈活性,能更好地適應(yīng)實際地下介質(zhì)結(jié)構(gòu),有更好的內(nèi)存使用率和計算效率。為了解決原始的基于正方形網(wǎng)格的PML(完全匹配層)吸收邊界形式不再適用的問題,文中將通過聲波方程完全匹配層吸收邊界的經(jīng)典方法,構(gòu)建并推導(dǎo)出PML吸收邊界條件在矩形旋轉(zhuǎn)交錯網(wǎng)格中的具體形式,通過Marmousi模型正演和反演實驗,得出本文所構(gòu)建的PML邊界獲得良好的吸收效果,對反演結(jié)果有很好的改善。

        頻率域全波形反演;矩形旋轉(zhuǎn)交錯網(wǎng)格;正演;完全匹配層

        頻率域全波形反演是在時域波形反演的基礎(chǔ)上發(fā)展起來的,經(jīng)過Marfurt[1]、Shin[2]、Worthington[3]、Pratt[4-5]等進一步地深入研究與完善,利用頻率域全波形反演的方法來重建地下介質(zhì)結(jié)構(gòu),越來越受到人們的關(guān)注。與時域全波形反演相比,頻域全波形反演計算效率更高,能處理更大的地下介質(zhì)模型,有更好的實用性。

        交錯網(wǎng)格最早由Madariage[6]提出,Pratt改進了經(jīng)典五點有限差分法求解波動方程,Jo[7]等人提出了九點旋轉(zhuǎn)交錯網(wǎng)格有限差分法來提高拉普拉斯域和質(zhì)量加速項的精度[8-9]。而這些方法是基于正方形網(wǎng)格的,這樣造成的影響是在一些領(lǐng)域無法適應(yīng)地下結(jié)構(gòu),同時可能增加不必要的內(nèi)存和計算量,這對于大模型全波形反演是一個難點。而實際地下介質(zhì)結(jié)構(gòu),橫向速度變化較小,縱向速度變化較大,基于矩形網(wǎng)格的全波形反演有更大的靈活性,能更好地適應(yīng)實際地質(zhì)結(jié)構(gòu),同時可減少網(wǎng)格點數(shù)目,有更好的內(nèi)存使用率和計算效率。

        在實際地震勘探中,地震波在地下是無限傳播的,而在我們波場模擬的過程中,由于計算機存儲的限制,在我們所需研究的區(qū)域外,加入了人為的截斷邊界。地震波會在人為邊界上產(chǎn)生反射,反射回的波會產(chǎn)生新的波場,對我們預(yù)期的波場產(chǎn)生干擾。因此,需要引入吸收邊界條件,消除人為截斷邊界的影響。現(xiàn)在應(yīng)用最廣泛的人工邊界條件是Berenger提出的完全匹配層法[10],有著比較理想的吸收效果。其基本思想[11]是,在模擬區(qū)域周圍添加完全匹配層,在完全匹配層中加入吸收衰減因子,吸收人為截斷邊界差生的反射波。

        在矩形旋轉(zhuǎn)交錯網(wǎng)格頻率域全波形反演中,隨著矩形網(wǎng)格寬長比的變化,原始的基于正方形網(wǎng)格的PML邊界形式已經(jīng)不再適用。因此,文中從聲波方程完全匹配層吸收邊界的經(jīng)典方法出發(fā),構(gòu)建并推導(dǎo)出PML邊界條件在矩形旋轉(zhuǎn)交錯網(wǎng)格中的形式。在原矩形旋轉(zhuǎn)交錯網(wǎng)格頻域全波形反演算法的基礎(chǔ)上,添加文中所構(gòu)建的矩形旋轉(zhuǎn)交錯網(wǎng)格中的PML邊界條件,以Marmousi模型為目標(biāo)模型,通過正演和反演實驗,表明文中所提方法的有效性。

        1 方法原理

        1.1 一般網(wǎng)格PML層聲波方程的離散形式

        二維線性標(biāo)量聲波方程形式如下[12]:

        式中P(x,z,t)為模型中(x,z)點處的聲波波場,v(x,z)為模型中(x,z)點處的介質(zhì)速度,S(x,z,t)為時域震源,?2為拉普拉斯算子。

        為簡單起見,不考慮PML層中震源S的作用。按照經(jīng)典PML邊界的理論[13],將波場P分裂為兩項,即P=Px+Pz,通過引入中間變量、加入衰減因子,可構(gòu)造如下PML吸收邊界條件:

        式中Q為中間變量,γ為PML層衰減系數(shù)。再將將式(2)變換回頻域、差分離散、消去中間變量、合并,可推出一般網(wǎng)格PML層聲波方程的離散形式為:

        式中Δx為橫向網(wǎng)格間距,Δz為縱向網(wǎng)格間距,ξ=1+iγ/ω。

        1.2 矩形旋轉(zhuǎn)交錯網(wǎng)格中PML邊界的構(gòu)建

        矩形旋轉(zhuǎn)交錯網(wǎng)格的基本思路如圖1所示,對任意點和其8鄰域點進行有限差分,將微分形式近似轉(zhuǎn)換為差分形式,對聲波方程進行離散,然后求解離散后的聲波方程。如圖2所示,里面包含兩種坐標(biāo)系,一種是直角坐標(biāo)系,另外一種是旋轉(zhuǎn)后的坐標(biāo)系。

        圖1 矩形網(wǎng)格示意圖

        圖2 旋轉(zhuǎn)坐標(biāo)系示意圖

        矩形旋轉(zhuǎn)交錯網(wǎng)格PML邊界的構(gòu)建分為3個部分:一是旋轉(zhuǎn)坐標(biāo)系,推導(dǎo)出聲波方程在新旋轉(zhuǎn)坐標(biāo)系下的具體形式;二是參照經(jīng)典PML邊界的構(gòu)建方法,推導(dǎo)出新旋轉(zhuǎn)坐標(biāo)系下的離散形式;三是將原正交坐標(biāo)系和新旋轉(zhuǎn)坐標(biāo)系下的PML邊界混合加權(quán),融合生成矩形旋轉(zhuǎn)交錯網(wǎng)格下的PML邊界。

        旋轉(zhuǎn)后坐標(biāo)系x′z′下的坐標(biāo)與原正交坐標(biāo)系xz下的坐標(biāo)關(guān)系為:

        對(4)式兩邊求二階微分并整理可得:

        由(1)(4)式可以得出旋轉(zhuǎn)后新坐標(biāo)系x′z′下的聲波方程為:

        與構(gòu)建一般網(wǎng)格PML邊界同樣的過程,可推導(dǎo)出旋轉(zhuǎn)新坐標(biāo)系x′z′下PML層聲波方程的離散形式為:

        設(shè)原正交坐標(biāo)系xz下的PML層所占比重為η,新旋轉(zhuǎn)坐標(biāo)系x′z′下的PML層所占比重為1-η,構(gòu)建權(quán)重因子如下:

        將式(3)、(6)分別乘以η和1-η,即可將兩個坐標(biāo)系下PML層聲波方程按權(quán)重融合即可得到矩形旋轉(zhuǎn)交錯網(wǎng)格中PML層的離散形式。至此,矩形旋轉(zhuǎn)交錯網(wǎng)格中PML邊界的離散形式構(gòu)建完畢。

        2 算法數(shù)值實驗

        為了驗證文中提出的基于矩形旋轉(zhuǎn)交錯網(wǎng)格PML邊界的有效性,文中將以Marmousi模型進行基于矩形旋轉(zhuǎn)交錯網(wǎng)格的頻域全波形反演實驗,從正演和反演結(jié)果兩個方面進行說明。在算法實現(xiàn)過程中,軟件平臺使用Visual Studio 2010,運行環(huán)境為64位windows7操作系統(tǒng),編程語言是C++,在處理線性方程組時調(diào)用MUMPS[14]庫來進行求解。

        Marmousi模型橫向距離范圍為0~9 200米,橫坐標(biāo)網(wǎng)格間距為12.5米,共737個網(wǎng)格;縱向深度0~3 000米,縱坐標(biāo)網(wǎng)格間距12.5米,共240個網(wǎng)格。為將Marmousi按矩形網(wǎng)格劃分,文中對Marmousi模型的縱向距離進行線性插值,在不改變整個模型大小的情況下,將模型縱坐標(biāo)網(wǎng)格間距減少到10米,縱向網(wǎng)格數(shù)目變?yōu)?60。炮點共240個,分布于地表,炮點覆蓋范圍分布于模型的橫坐標(biāo)從3 000米到8 975米處,炮間距25米。每個炮點的接收點均分布于該炮點左側(cè),最大炮檢距(炮點和接收點之間的距離)為2 575米,最小炮檢距200米,道間距25米,每個炮點有96個接收點,即96道。

        2.1 正演實驗

        正演是利用初始速度模型得到模擬波場記錄,是反演的前提和基礎(chǔ)。文中初始速度模型選取勻速遞增模型。正演的采樣時間為3 s,采樣間隔為4 ms,共750個采樣點。震源為Ricker子波,主頻為8 Hz。

        以Marmousi模型第一炮的數(shù)據(jù)為例,第一個炮點位于地表3 000米處,96個接收點分布在425~2 800米間,即接收道的第1個接收點位于425米處,第96個接收點位于2 800米,中間94個檢波點平均分布于此區(qū)域內(nèi)。分別記錄無PML層和添加本文所構(gòu)建的PML層(20層網(wǎng)格厚度)時,第一個炮點正演的波形記錄,其中橫坐標(biāo)表示96個接收點的道號,縱坐標(biāo)表示采樣時間。

        圖3(a)顯示,當(dāng)原原矩形旋轉(zhuǎn)交錯網(wǎng)格頻域全波形反演算法不添加PML吸收層時,在邊界處反射現(xiàn)象明顯。由于正演是反演的前提和基礎(chǔ),正演結(jié)果的精度將在很大程度上影響反演結(jié)果,所以邊界處產(chǎn)生的反射必將對反演過程帶來很大的干擾。圖3(b)顯示,加入文中所構(gòu)建的PML層后反射現(xiàn)象得到明顯抑制。由此可見,文中所構(gòu)建的PML邊界取得了良好的吸收效果,這也將會對反演結(jié)果帶來積極影響。

        2.2 反演實驗

        反演是建立理論波場與模擬波場之間殘差的目標(biāo)函數(shù),通過使目標(biāo)函數(shù)極小來更新速度模型,是一個對目標(biāo)函數(shù)最優(yōu)化的問題。文中使用的是預(yù)處理梯度算法[15]。

        圖3 Marmousi模型第一個炮點正演模擬的波形記錄

        圖4 Marmous模型的頻率域全波形反演

        圖4 (c)和4(d)分別為不添加PML層和添加文中所構(gòu)建PML層反演結(jié)果的反演結(jié)果,可以看出添加本文所構(gòu)建的PML層后,反演結(jié)果整體上結(jié)構(gòu)更清晰。圖4(e)給出了在水平距離左邊界2.50 km處,兩組結(jié)果的對比,可以看出添加本文構(gòu)建PML層的反演結(jié)果,速度分布更接近理論模型。這是由于PML層很好地吸收了地震波在人工邊界處的反射,減少了反射波的干擾,實驗結(jié)果更精確,從而更接近理論模型,這與正演的實驗結(jié)果也相符合,由此可以驗證文中所提算法的有效性。

        3 結(jié)束語

        文中構(gòu)建并推導(dǎo)出PML邊界在矩形旋轉(zhuǎn)交錯網(wǎng)格中的形式,通過正演和反演實驗,驗證其對人為邊界處的反射波有很好的吸收效果,對最終的反演結(jié)果也有很好的提升。另外,由于本文在選取PML層衰減因子時,沿用了正方形網(wǎng)格PML層衰減因子的形式。而隨著矩形網(wǎng)格寬長比的變化,構(gòu)建新的衰減因子也是一個值得研究的方向。

        [1]Marfurt K J.Accuracy of finite difference and finite element modeling of the scalar and elastic waveequations[J].Geophysics,1984,49(5):533-549.

        [2]Shin C S.Nonlinear elastic wave inversion by blocky parameterization [D].Tulsa,Oklahoma,United States:University of Tulsa,1988.

        [3]Pratt R G,Worthington M H.Inverse theory applied to multi-source cross-hole tomography.Part I:Acoustic wave-equation method[J].Geophys Prosp,1990,38(3):287-310.

        [4]Pratt R G.Inverse theory applied to multi-source cross-hole tomography, PartII:Elastic waveequation method[J].Geophys Prosp,1990,38(3): 287-310.

        [5]Pratt R G.Frequency-domain elastic wave modeling by finite differences:A tool for crosshole seismic imaging[J].Geophysics,1990,55(5):626-632.

        [6]Madariaga R.Dynamics of expanding circular fault [J].Bulletin of the Seismological Society of America,1976,66(3):639-666.

        [7]Jo CH,Shin C,Suh JH.An optimal 9-point,finite-difference,frequency-space,2-D scalar wave extrapolator[J].Geophysics,1996(61):529-537.

        [8]Erik H S,Norbert G,Serge A S.Modeling the propagation of elastic waves using a moodified finitedifference grid[J].Wave Motion,2000(31):7-792.

        [9]Reeshidev B,Mrinal K S.Finite-difference modeling of S-wave splitting in anisotropic media[J]. Geophysical Prospecting,2008(56):293-312.

        [10]Berenger J A.perfectly matched layer for the absorption of electromagnetic waves[J].Journal of Computational Physics,1994(114):185-200.

        [11]Hastings F D,Schneider J B,Broschat S L. Application of the perfectly matched layer(PML)absorbing boundary condition to elastic wave propagation[J].Journal of the Acoustical Society of America,1996,100(100):3061-3069.

        [12]Bjorn Engquist.Computationalhigh frequency wave propagation[J].Acta Numerica,2003:181-266.

        [13]Collino F,Tsogka C.Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media[J].Geophysics,2012,66(1):294-307.

        [14]Solver.MUltifrontal Massively ParallelSolver(MUMPS 4.10.0)Users guide[R].French:MUMPS,2011.

        [15]Operto S,Virieux J,Dessa J X,et al.Crustal imaging from multifold ocean bottom seismometers data by frequency-domain fullwaveform tomography: Application to the eastern ankaiTrough[J].Journal of Geophysical Research,2006,111(B9):171.

        Research and implementation of the PML boundary conditions based on rotated rectangular-grid modeling method

        ZHANG Liang,SUN Jin-yao
        (School of Physics and Technology,Wuhan University,Wuhan 430072,China)

        The present Frequency-domain full waveform inversion(FWI)is mostly default to a rotated square-grid,but the FWI based on rotated rectangular-grid is more flexible,and can better adapt to the actualstructure ofunderground media,and have greateradventage in memory utilization and computational efficiency.To overcome the problem that the original square-grid PML(Perfectly Latched Layer)absorbing boundary conditions no longer apply,the paper give the new form of PML absorbing boundary conditions based on rotated rectangular-grid through the classic building strategy.Using the forward modeling and inversion experiment with Marmousi model,we show that the proposed PML boundary in this article can obtain satisfactory absorbing effect,and improve the inversion result.

        frequency-domain full waveform inversion;rotate rectangular-grid;forward modeling;PML

        TN91

        :A

        :1674-6236(2017)02-0163-05

        2016-04-12稿件編號:201604110

        大型油氣田及煤層氣開發(fā)國家科技重大專項課題(2011ZX05003-003)

        張 亮(1991—),男,湖北棗陽人,碩士。研究方向:二維地震波頻率域全波形反演。

        猜你喜歡
        模型
        一半模型
        一種去中心化的域名服務(wù)本地化模型
        適用于BDS-3 PPP的隨機模型
        提煉模型 突破難點
        函數(shù)模型及應(yīng)用
        p150Glued在帕金森病模型中的表達(dá)及分布
        函數(shù)模型及應(yīng)用
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
        3D打印中的模型分割與打包
        久久天天躁狠狠躁夜夜av浪潮| 国产av一区二区三区丝袜| 亚洲一二三四五区中文字幕| 日韩有码在线免费视频| 国产日产欧产精品精品蜜芽| 国产精品免费看久久久8 | 九九99国产精品视频| 好看的国内自拍三级网站| 亚洲中文久久精品字幕| 亚洲av永久精品爱情岛论坛| 欧美精品中文| 亚洲天堂av另类在线播放| 自拍偷自拍亚洲精品第按摩 | 男女野外做爰电影免费| 精品国产一区二区三区亚洲人| 白色白色白色在线观看视频 | 中文字幕国产精品专区| 国产综合精品久久99之一| v一区无码内射国产| 亚洲线精品一区二区三区八戒| 国产又色又爽的视频在线观看91| 亚洲视频在线一区二区| 国产69精品久久久久999小说| 妺妺窝人体色www聚色窝韩国| av免费在线播放一区二区| 欧美xxxxx在线观看| 青草国产精品久久久久久| 国产成人亚洲综合无码品善网| 久久久久麻豆v国产精华液好用吗 欧美性猛交xxxx乱大交丰满 | 国产毛片三区二区一区| 亚洲成人免费av影院| 亚洲国产精品久久久久婷婷老年| AV无码中文字幕不卡一二三区| 色综合久久五十路人妻| 波多野结衣中文字幕一区二区三区| 亚洲五月天综合| 国产一极毛片| 日本一区二区日韩在线| 久久精品99国产精品日本 | 国产精品亚洲av无人区一区蜜桃| 日韩经典午夜福利发布|