高精度傅里葉有限差分法模型正演 - 圖文-

1、 2007 年 12 月 石油地球物理勘探 第 42 卷 第6期 作 者 介 紹 侯成福 高級(jí)工程師 , 1966 年生 ; 1989 年畢業(yè)于長(zhǎng)春 地質(zhì)學(xué) 院地球物理專業(yè) , 獲工學(xué)學(xué) 士學(xué)位 , 2002 年獲石油 大學(xué) ( 華東 礦產(chǎn)普查與勘探工學(xué)碩士學(xué)位 , 正在攻讀 中科院 地化專業(yè)博士學(xué)位。現(xiàn)在 東方地 球物 理公司 國(guó)際 部從 事野外地震采集工作。 丁仁偉 博士研究生 , 1980 年生 ; 2005 年畢業(yè)于中國(guó) 石油大 學(xué) ( 華東 , 獲計(jì)算機(jī) 應(yīng)用技 術(shù)專業(yè) 工學(xué) 碩士學(xué) 位 , 同年 到中國(guó)石油大學(xué) ( 華 東 地質(zhì)資 源與地 質(zhì)工程 專業(yè) 攻讀 博士學(xué)位 ; 主要從事

2、地震成像與偏移速度分析理 論與方 法、 高性能計(jì)算技術(shù)等方面的研究工作。 葉月明 博士研究生 , 1982 年生 ; 2004 年畢業(yè)于中國(guó) 石油大 學(xué) ( 華東 , 獲勘查技 術(shù)與工 程專業(yè) 學(xué)士 學(xué)位 , 同年 到中 國(guó)石油大學(xué) ( 華東 地球探測(cè)與信息專業(yè)攻讀碩士學(xué) 位 , 2006 年免試錄取為中國(guó)石油大學(xué) ( 華 東 地質(zhì)資源 與地 質(zhì)工程專業(yè)博 士研 究生 ; 主 要 從事 地震 資 料真 振 幅成 像、 地震照明度分析理論與方法以及復(fù)雜地質(zhì)條 件下的 地震偏移成像方法研究。 劉文革 博士研究 生 , 1972 年 生 ; 現(xiàn)主 要 從事 地震 正、 反演 方法研究及地震資料處理。

3、 朱生旺 教授級(jí)高級(jí)工程師 , 1962 年生 ; 1983 年畢業(yè) 于華東 石油學(xué)院物探專業(yè) , 現(xiàn)在中石化石油勘探開發(fā)研 究院從 事地球物理方法研究工作。 吳國(guó)忱 教授 , 1965 年生 ; 1987 年獲長(zhǎng)春地質(zhì)學(xué)院石 油物探 專業(yè)學(xué)士學(xué)位 ; 1996 年 獲石 油大 學(xué) ( 華 東 應(yīng) 用地 球物 理專業(yè)碩士學(xué)位 ; 2005 年 獲同 濟(jì)大 學(xué)固 體地 球物 理專 業(yè)博士學(xué)位 ; 現(xiàn)在中國(guó)石油大學(xué) ( 華東 從事地球 物理方 面的教學(xué)和科研工作。 張海燕 副教授 , 1968 年生 ; 1990 年畢業(yè)于清華大學(xué) 精密儀 器專業(yè) , 獲工學(xué)學(xué)士學(xué)位 ; 1993 年 畢業(yè)于 中國(guó)

4、海洋 大學(xué) 海洋物理專業(yè) , 獲理 學(xué)碩士 學(xué)位 ; 目前 正在職 攻讀 中國(guó) 海洋大學(xué)海洋地質(zhì)專業(yè)博士學(xué)位 , 師從李 慶忠院士。研 究方向?yàn)?A VO 疊前反演、 計(jì)算 智能。已 經(jīng)在國(guó)內(nèi) 外核 心學(xué)術(shù)期刊上發(fā)表文章十余篇 , 其中有多 篇被 EI 收錄。 楊旭明 高級(jí)工程師 , 1963 年生 ; 2002 年畢業(yè)于成都 理工大 學(xué) , 獲博士學(xué)位 , 主要從事地震波場(chǎng)正反演、 靜校 正及地 震散射降噪等方面的研究工作。 裴正林 博士 , 副研 究員 , 1962 年生 ; 2000 年畢業(yè) 于中 國(guó)地 質(zhì)大學(xué) ( 北京 , 獲博士學(xué)位 , 2003 年于石 油大學(xué) ( 北京 博士后出站

5、。主要從事地震波傳播理論及數(shù)值模擬、 井 間地震層析成像、 小波域地震屬性提取和油氣預(yù) 測(cè)以及 小波變換應(yīng)用研究等。 董 瑤 碩士研究生 , 1982 年生 ; 2005 年畢業(yè)于河北 工業(yè)大 學(xué)電子信息工程專業(yè) , 獲學(xué)士學(xué)位?，F(xiàn)在河北工 業(yè)大學(xué) 攻讀通信與信息系統(tǒng)專業(yè)碩士學(xué)位 , 主要從事計(jì) 算智能 與數(shù)據(jù)挖掘的研究 , 發(fā)表論文多篇。 嚴(yán)建文 高級(jí)工程師 , 1966 年生 ; 1989 年畢業(yè)于中國(guó) 地質(zhì)大 學(xué)( 武漢 應(yīng)用地球物理專業(yè) , 獲學(xué)士學(xué)位。一直 在中國(guó) 石化石油勘探開發(fā) 研究院南 京石油 物探 研究所 從事 專 業(yè)學(xué)術(shù)期刊的編輯 出版工 作和 專業(yè)情 報(bào)研 究工作。 現(xiàn)

6、為中國(guó)地質(zhì)大學(xué) ( 武漢 資源學(xué)院在職博士研 究生。 劉玉柱 講師 , 1979 年 生 ; 2001 年和 2004 年先后獲 得同 濟(jì) 大學(xué)地球物理專業(yè) 學(xué)士和 碩士 學(xué)位?，F(xiàn) 為同濟(jì) 大學(xué) 地 球物理專業(yè)博士研究生 , 研究方向?yàn)榈卣饘游龀上衽c 反 演研究。 李紅星 博士研究生 , 1981 年生 ; 現(xiàn)為吉 林大學(xué) 地球探 測(cè)科 學(xué)與技術(shù)學(xué)院固體地球物理專 業(yè)博士研究 生 , 主要 從事 地震波及聲波傳播特性及圖像處理研究。 黃捍東 副教授 , 碩 士生 導(dǎo)師 , 1964 年 生 ; 主要 從事 復(fù)雜 油 氣藏儲(chǔ)層地球物理方法的研究工作 , 完成多項(xiàng)中國(guó)石 油 和中國(guó)石化集團(tuán)公司 生

7、產(chǎn)、 科研 項(xiàng)目 , 獲 國(guó)家科 技進(jìn) 步 二等獎(jiǎng) 1 項(xiàng) , 省部級(jí)科技進(jìn)步一等獎(jiǎng) 1 項(xiàng) , 二 等獎(jiǎng) 2 項(xiàng) , 授權(quán)發(fā)明專利 1 項(xiàng) , 發(fā) 表論 文 20 余篇。 現(xiàn)在中 國(guó)石 油 大學(xué) ( 北京 資源與信息學(xué)院從事教學(xué)和科研工作。 張明振 高級(jí)工程師 , 1967 年生 ; 1991 年畢 業(yè)于成都地質(zhì) 學(xué) 院石油地質(zhì)勘 查專 業(yè) , 1999 年獲 同濟(jì) 大學(xué) 固體 地球 物 理專業(yè)理學(xué)碩士學(xué)位 , 現(xiàn)為中國(guó)石油大學(xué) ( 華 東 地質(zhì) 資 源與地質(zhì)工程專業(yè)博士研究生 , 主要從事地震資料的 綜 合研究工作。 路順行 高 級(jí)工程師 , 碩士 , 1961 年 生 ; 1985 年

8、畢業(yè)于山 東 海洋學(xué) 院海 洋 地 質(zhì) 地球 物 理 專 業(yè) , 獲 理 學(xué)學(xué) 士 學(xué) 位 ; 2003 年畢業(yè)于中國(guó) 海洋 大學(xué)地 質(zhì) 工程 專業(yè) , 獲 工程 碩 士學(xué)位?，F(xiàn)在中國(guó)石 化勝利 油田分 公司 地質(zhì)科 學(xué)研 究 院從事石油地質(zhì)勘探研究工作 , 同時(shí)為中國(guó)海洋大學(xué) 海 洋地球科學(xué)學(xué)院在 讀博士 研究 生。一直 在勝利 油田 地 質(zhì)科學(xué)研究院從事石油地質(zhì)勘探研究、 勘探生產(chǎn)技術(shù) 管 理等工作。 薛國(guó)強(qiáng) 工程師 , 1966 年生 ; 2002 年在長(zhǎng)安 大學(xué)獲得碩士 學(xué) 位 , 2005 年獲西安交通大學(xué)博士 學(xué)位 , 現(xiàn)在中 國(guó)科學(xué) 院 地質(zhì)與地球物理研究所博士后流動(dòng) 站做研究工

9、作。 駱 遙 碩士研究生 , 1982 年生 ; 2005 年畢 業(yè)于中國(guó)地質(zhì) 大 學(xué) ( 北京 地球物 理學(xué)專業(yè)。 現(xiàn)在中 國(guó)科 學(xué)院地 質(zhì)與 地 球物理研究所空間物理學(xué)專業(yè)攻讀 碩士學(xué)位。 王衛(wèi)華 教授級(jí)高級(jí)工程師 , 中國(guó)石油集團(tuán)公司高級(jí)技術(shù) 專 家 , 1952 年生 ; 1990 年 畢業(yè) 于成都 理工 大學(xué) 研究 生院 , 獲碩士學(xué)位。早年從事海上地震勘探工作 , 近年主要 負(fù) 責(zé)我國(guó)西部沙漠和山 地山前 帶等復(fù) 雜地 區(qū)的地 震勘 探 技術(shù)工作 , 現(xiàn)任東方地球物理公司副總工程師。先后 從 事地震資料解釋方法、 數(shù)據(jù)處理方法和采集方法研究 等 工作 , 發(fā)表論文 20 多篇。 王

10、喜雙 高級(jí)工程師 , 1957 年生 ; 1982 年畢 業(yè)于華東石油 學(xué) 院勘探系石油地質(zhì)勘探專 業(yè) , 獲學(xué)士 學(xué)位 , 1995 年畢 業(yè) 于國(guó)家地震局地質(zhì)所 , 獲博士學(xué)位。長(zhǎng)期從事物探技 術(shù) 管理工作。 O il Geophysical Pr ospecting 2007 by t he t est of M armousi model. Key words: cont ro lled illuminat ion, finite! diff erence tr av el!t im e com put at ion, synt het ic operat or, w av e equ

11、at ion, prestack dept h m igrat io n Ding Ren!wei, Co lleg e of Eart h Resource and Inf or! m at ion, China U niv ersit y of P et ro leum, Dongying Cit y, Shandong Pro vince, 257061, China Synthetic plane wave migration in relief surface con! dition and implementation in paral lel way. Ye Yue! ming,

12、 Li Zhen!chun, Tong Zhao !qi and An Qi. OGP , 2007, 42( 6 : 622 628 Wav e equat ion!based co nt inuat io n approach keeps t he f eat ure of dynam ics and kinem at ics o f w av ef ield, can bett er overco me t he short com ing s o f tim e!shift st at ic correct ion approach that dist ort s w av ef ie

13、ld and elim inat e t he inf luence of relief sur! f ace, and ref orm s records in non!horizo nt al survey! ing int erf ace t o t hat in horizo nt al sur veying inter! f ace ( new dat um . Synt hesis of plane w ave source and it s r ecords in r econst ruct ive horizo nt al dat um can overcome the iss

14、ues t hat is incapable to synt he! size t he plane w ave source and it s records because of relief g eo morpholog y, finishing plane w av e mi! g rat io n in r elief surf ace condit ion. In o rder t o fur! t her improv e operation eff iciency, and acco rding t o the principl e dist ribut ing t he pl

15、ane w av es in dif fer! ent incident angl es and f requency range t o separat e no des, t he M PI par al lel computat ion w as im ple! m ent ed, w hich do ubled and redo ubled t he computa! t ional ef ficiency. T he co ncret e steps are as f ollo w s: continuing the CSP g athers to new datum that is

16、 below t he low est relief plane; ! so rt ing CRP g ath! ers f rom cont inued CSP g athers dat a; continuing the CRP gat hers to new dat um t hat is below t he low est relief plane; . sor ting CSP g athers fr om CRP g at her s cont inued t o new dat um . T he to tal abov em ent io ned pro cess corr

17、esponds t o dow nshift the g eo metr y t o new dat um. T he results o f m odel tests and pract ical dat a processing show ed t hat t he result s o f parallel plane w ave dept h mig ration af ter reconst ruct ion of dat um are same as t hat of dept h migrat ion dat a by directly dow n!co nt inuing w

18、av e! f ield based on raw shot records, but the for mer comput at ional eff iciency doubled and r edoubled a! long w it h incr easing t he number of applied no des. Key words: relief surf ace, prest ack dept h m ig ra! t ion, dat um r econst ruct ion, plane w ave migrat ion, synthet ic shot reco rd,

19、 M PI parallel Ye Yue!ming, College of Eart h Resource and Infor! mat ion, China Universit y of Pet roleum, Dong ying Cit y, Shandong P rovince, 257061, China Forward modeling by Fourier finite !difference ap! proach. Liu Wen!ge, He Zhen!hua, Huang De!ji and Du Zeng ! li. OGP, 2007, 42( 6 : 629 633

20、T here are tow kinds of w ave equat ion num eric m odeling met hods: f init e!dif ference appr oach and approach in f requency!w avenumber dom ain, amo ng w hich t he comput at ion precisio n o f f init e!dif f erence m et hod depends on t he approx imat e degree of w avefield ex tr apolat ion opera

21、t or, discret e g rid in! t er val and o rder of finite! diff erence equat ion, it is a! dapted fo r any lat er al variat ions o f v elo cit y, but the dispersion and backg round no ise are easily presen! t ed in larg e dip. T he alg orit hm in f requency!w ave! num ber dom ain is charact er ized by

22、 sim ple, hig h! precision, sm all noise and adapt ive t o any dips of str at a, but inadequat e f or any lateral variat ions of w avefields. Combining t he advant ages of fr equency! w avenumber domain w it h F init e!Diff erence ap! proach, the paper used Fo urier F init e!Dif f erence ( F FD appr

23、oach t o finish the numeric sim ulat ion of seismic reco rds of m odel by using hig h!precision continuat ion operato r and in multi!dom ain. T he ex! t rapo lat ion oper at or of w av ef ield is com po sed of phase!shift it em, ref ract ion it em ( t ime!shif t it em and finite!diff erence com pens

24、ation it em. T he t heo! ret ical and er ror analyses are carried out f or F FD m et hod, and one!w ay acoust ic w av e equat ion is used t o do num eric simulation t est for lay er ed mod! el and SEG/ EAGE salt m odel respect ively. T he correlat io n of numeric tests show ed t he FF D met h! od is

25、 adapt ive f or sev ere lat er al change of velocity f ield and characterist ics of high pr ecision, no!dis! persion and w eak backg round noise, and t he ref lec! V ol. 42 N o. 6 A bstracts / t ion f eat ure is distinguishable in simulated result s, w hich can carry o ut accur at e numeric seismic

26、simu! l at ion of com plex geolo gic st ruct ure. Key words: w ave equat ion, numeric seismic sim ula! t ion, f init e!dif ference met hod, approach in f requen! cy!w av enum ber dom ain, F FD Liu Wen!ge, Col leg e of Inf ormat ion Engineering, Cheng du U niversit y of Science and Engineering, Cheng

27、 du Cit y, Sichuan Pro vince, 610059, China Numeric simulation by grid!various finite!di fference elastic wave equation. Zhu Sheng !wang, Qu Shou!li, Wei Xiu!cheng and Liu Chun!yuan. OGP , 2007, 42( 6 : 634 639 St udy ing seism ic w av e propag at ing rule and seismic respo nse in com plex medium ne

28、eds carrying out numeric sim ulat ion of elast ic w ave equat ion un! der tiny gr id separat ion, w hich brings an issue o f giant comput at ional w orks, and using grid!various com put at ion is ef fective w ay t o reduce computa! t ional w or ks. T he paper gives a m et hod for im ple! m ent ing g

29、 rid !v ar io us f init e!dif ference com put ation: the computat ion in t iny grid is used in local com! plex medium ar ea, t he comput at ion in gross g rid is carried o ut in r est area. In t ransit io nal area betw een tw o g rids, the t ransit ional link is finished by chan! g ing diff er ent i

30、al oper at or and w av ef ield int er pola! t ion. T he t heor et ical analysis and numeric sim ula! t ion show ed t hat it do es not aff ect t he sim ulat ion o f seismic w ave propag ation in t he tr ansit ional area bet w een t he gross and tiny g rid, w hich at t ained t he goal bo th reducing c

31、o mput ational w orks and g uar! ant eeing t he precision. Key words: elastic w av e equat io n, num er ic sim ula! t ion, g rid ! various f init e dif f erence, dispersio n Zhu Sheng !wang, Instit ute of P et ro leum Explora! t ion and Depart ment o f SINOPEC, No. 31, Xueyuan Road, H aidian Dist ri

32、ct , Beijing Cit y, 100083, China High!order one !way generalized!screen propagation operator of qP!wave in VTI medium. Wu Guo !chen, Liang Kai and Wang Hua ! zhong. OGP, 2007, 42( 6 : 640 650 St art ing f rom elastic w av e equatio n in VT I m edium and t hrough eig envalue equatio n and a! co usti

33、c w ave t heo ry, t he paper approx imat ely crea! t ed dispersion equat ion o f qP!w av e in VT I medium and deeply st udied high!o rder o ne!w ay generalized! screen propagat io n operato r of qP!w ave in VT I me! dium based on dist ur bed t heory of medium . Disso l! ving the par am et er f ield

34、in VT I medium int o uni! f orm anisot ro pic velocity f ield, v elo cit y dist urbed f ield and anisot ro pic parameter dist ur bed f ield, the high!o rder one!w ay g eneralized !scr een pr opag at ion operato r of qP!w ave is composed of tw o part s: free propag ation it em of w avefield in backg

35、round m edi! um and anisot ropic dist urbed item. T he co mput a! t ion of backg round f ield is implement ed in w av enu! mber domain and no n!unifo rm anisot ropic dis! t ur bed item ( including slow ness dist urbed and ani! so tr opic par am et er dist urbed items is co mput ed in space domain. When vio lent lateral var