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

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

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

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

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

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

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

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

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

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

10、喜雙 高級工程師 , 1957 年生 ; 1982 年畢 業(yè)于華東石油 學 院勘探系石油地質(zhì)勘探專 業(yè) , 獲學士 學位 , 1995 年畢 業(yè) 于國家地震局地質(zhì)所 , 獲博士學位。長期從事物探技 術(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