發(fā)布時間：2021-03-10 19:43

　　為了精確求解零井源距VSP波場記錄,基于體積分方程提出了一種預(yù)條件最小二乘的數(shù)值計算方法。根據(jù)聲波散射理論,將頻率域一維聲波方程通過格林函數(shù)表達為Lippman-Schwinger積分方程;在此基礎(chǔ)上利用矩形積分公式將第二類Fredholm積分方程的求解問題轉(zhuǎn)化為大型線性方程組的求解,為后續(xù)的精確計算奠定了基礎(chǔ);然后,使用最小二乘方法直接求解全波場,而不是通過迭代方法逐步將散射分量加入到全波場中;通過引入預(yù)條件算子,使數(shù)值計算過程更穩(wěn)定;最后,以階梯狀模型和用實際測井?dāng)?shù)據(jù)建立的模型為例,利用提出的方法求解VSP波場,并與常用的有限差分方法結(jié)果對比,驗證了利用所提方法進行零井源距VSP波場模擬的有效性。 

【文章來源】：石油地球物理勘探. 2020,55(03)北大核心

【文章頁數(shù)】：9 頁

【部分圖文】：

圖2 層狀模型不同方法正演的零井源距VSP記錄

表1 不同方法的計算時間對比 數(shù)值模擬方法 時間/s 本文方法 240.3173 Neumann級數(shù)迭代法 350.2856 時間2階、空間2階有限差分法 1.2632 時間2階、空間8階有限差分法 2.0968不同頻率下格林函數(shù)矩陣G的條件數(shù)如圖4a 所示,可見當(dāng)正演頻率較高時,格林函數(shù)矩陣G的條件數(shù)會存在一些異常大的值以及振蕩。圖4b 展示了矩陣A(紅色虛線)和DA(藍色實線)條件數(shù)的對數(shù)及二者之差的對數(shù)(黑色實線)隨頻率的變化情況。由圖可知,在150Hz范圍以內(nèi),隨著頻率的增加,條件數(shù)逐漸增大;在頻率范圍150～200Hz條件數(shù)變化比較劇烈;頻率超過200Hz時,條件數(shù)逐漸增加,同時,二者條件數(shù)之差也有相似的規(guī)律。值得注意的是,頻率較高時,矩陣DA比矩陣A的條件數(shù)大約減小了10%,亦即在理論上,引入預(yù)條件算子D小幅度提高了計算的穩(wěn)定性。另外,在每個頻率下,矩陣DA的條件數(shù)都大于0,根據(jù)條件數(shù)的定義

不同頻率下三個矩陣的條件數(shù)

【參考文獻】：
期刊論文
[1]基于分段快速模擬退火的零偏VSP全波形反演[J]. 韓璇穎,印興耀,曹丹平,梁鍇.  石油物探. 2019(01)
[2]復(fù)雜區(qū)深井VSP層析求速度場的方法[J]. 侯愛源,李慶忠,張文波.  石油地球物理勘探. 2017(06)
[3]井中地震技術(shù)的昨天、今天和明天——井中地震技術(shù)發(fā)展及應(yīng)用展望[J]. 趙邦六,董世泰,曾忠.  石油地球物理勘探. 2017(05)
[4]基于自適應(yīng)優(yōu)化有限差分方法的全波VSP逆時偏移[J]. 蔡曉慧,劉洋,王建民,王維紅,任志明.  地球物理學(xué)報. 2015(09)
[5]基于Backus等效平均的測井資料尺度粗化方法研究[J]. 曹丹平.  石油物探. 2015(01)
[6]基于高斯射線束的斜井VSP正演方法[J]. 楊飛龍,孫淵,李緒宣,王旭,徐長志,俞岱.  地球物理學(xué)進展. 2014(06)
[7]交錯網(wǎng)格任意階導(dǎo)數(shù)有限差分格式及差分系數(shù)推導(dǎo)[J]. 楊慶節(jié),劉財,耿美霞,馮晅,郭智奇,劉洋.  吉林大學(xué)學(xué)報(地球科學(xué)版). 2014(01)
[8]對比源反演算法在二維彈性波成像中的應(yīng)用[J]. 李杰,繆竟鴻.  系統(tǒng)工程與電子技術(shù). 2012(08)
[9]用于零偏移距VSP資料的自適應(yīng)波形反演方法研究[J]. 高靜懷,汪超,趙偉.  地球物理學(xué)報. 2009(12)
[10]標量聲波波動方程高階交錯網(wǎng)格有限差分法[J]. 陳可洋.  中國海上油氣. 2009(04)

碩士論文
[1]多維第二類Fredholm積分方程高精度數(shù)值算法研究[D]. 秦博.電子科技大學(xué) 2015



本文編號：3075179