海洋可控源電磁法二維有限元正演及反演.pdf

1、進入21世紀,石油和天然氣勘探的重點已由陸地向海洋轉(zhuǎn)移。三維海洋地震勘探一直是為海洋鉆探提供海底地質(zhì)構(gòu)造特征的重要手段,但石油公司面對海洋油氣鉆探的巨額風(fēng)險,并不滿足于僅獲知海底地質(zhì)構(gòu)造的特征,還希望在鉆探前得知構(gòu)造內(nèi)流體的性質(zhì),以便提高鉆探成功率。對巖性變化十分敏感的電阻率一直是油層解釋(如電阻率測井資料解釋)的重要標(biāo)志性參數(shù)之一。海洋可控源電磁法是勘探海底油氣資源的一種地球物理方法,通過觀測油氣儲層和低阻沉積圍巖之間的電阻率差異所產(chǎn)

2、生的電磁異常,從而可以確定油氣儲層的分布范圍。在國外,海洋可控源電磁裝備研制和方法研究已較為成熟,正反演研究成果較多,而我國海洋可控源電磁探測技術(shù)的研究相對滯后。因此,我們需要加強對海洋可控源電磁場的數(shù)值模擬和反演方法的研究工作。本論文在總結(jié)前人研究成果的基礎(chǔ)上,深入探討了海洋可控源電磁一維和二維正反演方法和關(guān)鍵技術(shù),主要研究內(nèi)容和取得的研究成果如下:
　　(1)從麥克斯韋方程組出發(fā),利用矢量位函數(shù),詳細推導(dǎo)了直角坐標(biāo)系下三個方向

3、上電偶源分別在層狀各向同性和垂直各向異性介質(zhì)中的電磁場解析表達式。利用水平旋轉(zhuǎn)角與傾角來分解任意取向的電偶源,求得三個方向上的電偶極距,采用漢克爾變換計算含虛宗量貝塞爾函數(shù)的無窮積分,得到三個方向電偶源產(chǎn)生的電磁場,經(jīng)矢量疊加,即可求得任意取向電偶源的電磁場。采用積分法或者剖分法計算有限長電偶源的電磁場,并分析了這兩種方法的計算效率及精度。在收發(fā)距較小的情況下,不能將有限長電偶源看作為電偶極子計算電磁場響應(yīng),而應(yīng)采用積分法來計算。用Fo

4、rtran計算機語言編寫了計算任意取向有限長電偶源電磁響應(yīng)的計算機程序,并驗證了算法的正確性和精度。
　　(2)以空氣-海水-巖層為背景模型,分析了電磁場響應(yīng)的振幅與海底巖層電阻率之間的關(guān)系,提出了一種計算等效視電阻率的方法,以海洋油氣模型為例證明了等效視電阻率方法的有效性。
　　(3)從層狀各向同性介質(zhì)的電磁場解析表達式出發(fā),詳細推導(dǎo)了電磁場對各層電導(dǎo)率的偏導(dǎo)數(shù),采用高斯牛頓方法,實現(xiàn)了海洋可控源電磁場一維反演。
　

5、　(4)海洋可控源電磁場2D正演模擬中采用自適應(yīng)有限元法分別求解關(guān)于總場和異常場的邊值問題。異常場法中,通過水平層狀模型的一次場解析解和二次場有限元數(shù)值解的矢量疊加得到電磁場響應(yīng)。
　　對海洋可控源電磁2D正演模擬中的難點問題進行了詳細的測試和分析:
　　1)采用偽delta函數(shù)對源進行離散,分析了不同離散方法的精度;
　　2)針對海洋可控源電磁場2D有限元模擬的線性方程組,采用不同的求解器求解,確定了最優(yōu)的求解器;<

6、br>　　3)有限元求解過程中,采用非結(jié)構(gòu)化三角單元網(wǎng)格模擬復(fù)雜地電模型,利用后驗誤差估計方法指導(dǎo)網(wǎng)格加密細化,這樣在同樣網(wǎng)格數(shù)量下能夠得到更高精度的有限元解,該方法具有精度高、易于模擬復(fù)雜地電模型等優(yōu)點;
　　4)由于海洋可控源電磁法正演計算耗時較長,計算時間難以接受,正演模擬的并行計算必不可少。
　　(5)在正演研究基礎(chǔ)上,開展了海洋可控源電磁場2D反演方法研究。詳細分析了海洋可控源電磁場2D反演中的關(guān)鍵問題,如正則化因