段宏凱 張曉丹 佘翼翀 謝寶林 楊白雪

摘要：在地震波場數(shù)值模擬中，如何使計(jì)算結(jié)果兼?zhèn)涓咝耘c準(zhǔn)確性一直是地震勘探學(xué)者研究的熱點(diǎn)之一. 為保證計(jì)算精度，傳統(tǒng)方法采用統(tǒng)一的較小步長進(jìn)行有限差分計(jì)算，嚴(yán)重降低了正演模擬的計(jì)算效率.提出了一種根據(jù)地質(zhì)模型特點(diǎn)對不同的速度層采用不同尺度網(wǎng)格，并優(yōu)化過渡帶的方法，首先，分析研究對象的速度模型從而確定多尺度網(wǎng)格劃分規(guī)則，其次，確定過渡帶范圍，最后，求取過渡帶內(nèi)外各點(diǎn)波場系數(shù)與差分點(diǎn)數(shù)，最終獲得整個模型網(wǎng)格點(diǎn)的波場值. 通過文中的實(shí)驗(yàn)結(jié)果可見，在保證與常規(guī)網(wǎng)格方法模擬精度不變的前提下，使用多尺度網(wǎng)格方法進(jìn)行數(shù)值模擬的計(jì)算效率明顯提高，就本文的算例而言，平均可以達(dá)到25.16%.

關(guān)鍵詞：多尺度網(wǎng)格；步長；過渡帶；波場模擬；計(jì)算效率

中圖分類號：P631 文獻(xiàn)標(biāo)識碼：A

文章編號：1009-3044（2019）07-0252-03

現(xiàn)代油氣勘探面臨著向地表更復(fù)雜更深層區(qū)域推進(jìn)的現(xiàn)狀，這就對地震資料的采集，處理和解釋提出了更高的要求。數(shù)值模擬方法因其經(jīng)濟(jì)性和準(zhǔn)確性，得到了廣泛的應(yīng)用. 常用的數(shù)值模擬方法有偽譜法，有限差分法，有限元法，有限差分法以其算法易于實(shí)現(xiàn)，運(yùn)算速度較快，模擬結(jié)果分辨率高等優(yōu)點(diǎn)獲得了廣泛應(yīng)用[1]. 傳統(tǒng)有限差分是基于統(tǒng)一尺度的矩形網(wǎng)格進(jìn)行空間劃分，均是以固定的網(wǎng)格步長對整個模型區(qū)域進(jìn)行離散化處理[2]，但是在復(fù)雜的地質(zhì)體和巖性構(gòu)造區(qū)域，介質(zhì)分布是高度不均勻的，而且這些介質(zhì)通常都是帶有低速夾層，這些精細(xì)巖層需要用更小采樣間隔來提高模擬精度，抑制頻散和保障穩(wěn)定性[3]. 但是整個模型區(qū)域剖分越細(xì)，計(jì)算量會成倍增加，從而極大地降低了模擬效率，并且高速層會出現(xiàn)過采樣.為解決以上問題，有的學(xué)者采用提高差分階數(shù)的方法來避免頻散：岳曉鵬[4]使用時間4階、空間2N階差分進(jìn)行正演模擬，李斌[5]在時間域和空間域都使用高精度差分，提高了正演模擬精度；有的學(xué)者則從交錯網(wǎng)格入手，Wei-Zhong Wang[6]提出了多尺度旋轉(zhuǎn)交錯網(wǎng)格，降低了計(jì)算時間，Hongyong結(jié)合極小極大近似和泰勒級數(shù)展開，提出了最優(yōu)變網(wǎng)格有限差分，王建[7]使用余弦函數(shù)修正二項(xiàng)式窗改進(jìn)的交錯網(wǎng)格進(jìn)行數(shù)值模擬，有效控制了數(shù)值頻散；還有的學(xué)者另辟蹊徑，馬吉浩[8]]綜合了不同的地震資料進(jìn)行了正演模擬，梁全文[9]使用改進(jìn)的線性方法確定新模板的有限差分系數(shù)，并提高了計(jì)算精度. 以上學(xué)者都采用固定大小的網(wǎng)格進(jìn)行模擬，因此，一些學(xué)者通過改變網(wǎng)格步長來解決頻散問題：Moczo首次使用網(wǎng)格連續(xù)變化的方法來減少計(jì)算量，提高數(shù)值模擬的效率，1994年，Jastram和Behle[10]提出了二維聲波方程針對某一深度變網(wǎng)格步長的算法，張劍鋒[11]基于彈性波應(yīng)力-速度方程給出了交錯計(jì)算應(yīng)力及速度的非規(guī)則網(wǎng)格差分法，黃超和董良國[12]通過將變化的空間網(wǎng)格與變化的時間步長技術(shù)相結(jié)合，提出了一種空間網(wǎng)格大小與時間步長均可任意變化的高階有限差分模擬方法. 上述的算法，使計(jì)算效率和存儲空間的使用情況都有了較大的改進(jìn)，但在進(jìn)行計(jì)算時，在網(wǎng)格步長突變處的鄰域存在著一個過渡帶，這些算法中一部分是對過渡區(qū)域進(jìn)行插值，這會增加算法的不穩(wěn)定性，另一些是使用平滑函數(shù)對過渡帶區(qū)域進(jìn)行平滑，但這樣使差分算法更難實(shí)現(xiàn)。因此，本文提出了一種無須在過渡帶進(jìn)行插值的多尺度網(wǎng)格算法，能夠有效地解決上述問題.

1 波動方程及其差分格式

4 總結(jié)

通過實(shí)驗(yàn)驗(yàn)證，由于在低速區(qū)域使用小步長網(wǎng)格進(jìn)行采樣，并在過渡帶引入了對稱的有限差分法，較之傳統(tǒng)有限差分法，多尺度網(wǎng)格步長有限差分平均減少計(jì)算時間25.16%，平均節(jié)省儲存空間21.89%。因此，多尺度網(wǎng)格步長有限差分算法具有以下特點(diǎn)：使得使用小步長網(wǎng)格采樣的模型內(nèi)存需求量降低；計(jì)算效率明顯提高；無須在過渡帶進(jìn)行插值，有效地避免了由插值引起的誤差累積. 因而對含有層狀速度地層、階梯狀地層進(jìn)行數(shù)值模擬時有均良好的效果.

參考文獻(xiàn)：

[1]Wang， Fangbo；Sett， Kallol. Time-domain stochastic finite element simulation of uncertain seismic wave propagation through uncertain heterogeneous solids[J].Soil Dynamics and Earthquake Engineering：2016，88：369-385.

[2]張賡，庹先國，張毅. 淺層地震反射波法數(shù)據(jù)采集參數(shù)正演模擬研究[J]. 工程勘察：2015， （1）：93-98.

[3]Jian-Ping Huang；Ying-Ming Qu. Variable-coordinate forward modeling of irregular surface based on dual-variable grid[J]. Applied Geophysics：2015，12（1）：101-110.

[4]李斌，溫明明，牟澤霖. 用于聲波正演模擬的時空域高精度交錯網(wǎng)格有限差分方法[J]. CT理論與應(yīng)用研究：2017，3：259-266.

[5]岳曉鵬，白超英，岳崇旺. 高階交錯網(wǎng)格有限差分彈性波場模擬的精度分析[J]. 煤田地質(zhì)與勘探：2017，1，125-130.

[6]Wei-Zhong Wang，Tian-Yue Hu，x. Variable-order rotated staggered-grid method for elastic-wave forward modeling[J]. Applied Geophysics：2015，12（3）：389-400.

[7]王建，孟小紅，劉洪. 余弦修正二項(xiàng)式窗的交錯網(wǎng)格有限差分地震正演方法[J]. 應(yīng)用地球物理（英文版）：2017，1：115-124.

[8]馮吉浩，曹丹平，秦海旭. 基于波動方程正演模擬的多尺度地震資料反射特征分析[J]. 地球物理學(xué)進(jìn)展 ：2016 （3） ：1058-1065.

[9]梁文全，賴永明，姜琳. 適于聲波方程數(shù)值模擬的時間-空間域有限差分算子改進(jìn)線性算法[J]. 地震工程學(xué)報：2016，5：815-821.

[10]Jastram C，Tessmer E. Elastic modeling on a grid with vertically varying spacing[J].Geophysics Prospecting：1994，42（4）：357-370.

[11]張劍鋒. 彈性波數(shù)值模擬的非規(guī)則網(wǎng)格差分方法[J].地球物理學(xué)報：1998，41（增刊）：357-366.

[12]董良國，李培明. 地震波傳播數(shù)值模擬中的頻散問題[J].天然氣工業(yè)，2004，24（6）：53-56.

【通聯(lián)編輯：梁書】