参考文献
1 Shipley T H，Houston M H，Buller R T，et al.Seismic evidence for wide-spread possible gas hydrate horizons on continental slopes and margins.Am.Assoc.Pet.Geol.Bull.，1979，63，2204～2213
2 宋海斌，耿建华，Wong How-Kin等.南海北部东沙海域天然气水合物的初步研究.地球物理学报，2001，44（5）：684～691
3 宋海斌，松林修，仓本真一.西南海海槽地震资料处理及其BSR特征.地球物理学报，2001，44（6）：785～790
4 Singh S C，Minshull T A，Spence G D.Velocity structure of a gas hydrate reflector.Science，1993，260，204～207
5 MacKay M E，Jarrard R D，Westbrook G K，et al.Origin of bottom simulating reflectors：Geophysical evidence from the Cascadia accretionary prism.Geology，1994，22，459～462
6 Singh S C，Minshull T A.Velocity structure of a gas hydrate reflector at Ocean Drilling Program site 889 from a global seismic waveform inversion.J.Geophys.Res.，1994，99，24221～24233
7 Yuan T，Hyndman R D，Spence G D，et al.Seismic velocity increase and deep-sea hydrate concentration above a bottom simulating reflector on the northern Cascadian slope.J.Geophys.Res.，1996，101：13655～13671
8 Yuan T，Spence G D，Hyndman R D，et al.Seismic velocity studies of a gas hydrate bottom simulating reflector on the northern Cascadia continental margin：Amplitude modeling and full waveform inversion.J.Geophys.Res.，1999，104，1179～1191
9 Minshull T A，Singh S C，Westbrook G K.Seismic velocity structure at a gas hydrate reflector，offshore western Columbia，from full waveform inversion.J.Geophys.Res.，1994，99，4715～4734
10 Pecher I A，Minshull T A，Singh S C，et al.Velocity structure of a bottom simulating reflector offshore Peru：Results from full waveform inversion.Earth Planet.Sci.Lett.，1996，139，459～469
11 Korenga J，Holbrook W S，Singh S C，et al.Natural gas-hydrates on the southeast US margins：Constraints from full waveform and traveltime inversions of wide-angle seismic data.J.Geophys.Res.，1997，102，15345～15365
12 Pecher I A，Ranero C R，von Huene R，et al.The nature and distribution of bottom simulating reflectors at the Costa Rican convergent margin.Geophys.J.Int.，1998，133，219～229
13 Sain K，Minshull T A，Singh S C，et al.Evidence for a thick free gas layer beneath the bottom simulating reflector in the Makran accretionary prism.Marine Geology，2000，164，3～12
14 Kennett B L N，Kerry N J.Seismic waves in a stratified half-space.Geophys.J.R.Astron.Soc.1979，57，557～583(https://www.daowen.com)
15 Kirkpatrick S，Gelatt C D，Vecchi M P.Optimization by simulated annealing.Science，1983，220，671～680
16 Ingber L.Very fast simulated reannealing.Math.Comput.Modeling.，1989，12，967～973
17 Kormendi F，Dietrich M.Non-linear waveform inversion of plane wave seismograms in stratified elastic media.Geophysics，1991，56，664～674
18 Hamilton E L.Vp/Vs and Poisson's ratios in marine sediments and rocks.J.Acou.Soc.Am.，1979，66，1093～1101
19 Posewang J，Mienert J.The engima of double BSRs：Indicators for changes in the hydrate stability field？Geo-Marine Lett.，1999，19，150～156
20 Song Haibin，Ma Zaitian，Zhang Guanquan.Simultaneous inversion of layered compressional velocity and shear velocity by using plane wave seismogram.Science in China，1996，39（6）：618～628
21 宋海斌，马在田，张关泉.层状横向各向同性介质反问题初探.地球物理学报，1997，40（1）：105～119
22 宋海斌.层状介质弹性参数反演研究.北京：高等教育出版社.2002，1～129
FULL WAVEFORM INVERSION OF GAS HYDRATE-RELATED BOTTOM SIMULATING REFLECTORS
SONG Haibin1，2 MATSUBAYASHI Osamu2 KURAMOTO Shin'ichi2
（1.Institute of Geology and Geophysics，CAS，Beijing 100029，China；2.National Institute of Advanced Industrial Science and Technology，Tsukuba 3058567，Japan）
Abstract：A full waveform inversion strategy is proposed for gas hydrate-related bottom simulating reflector.This method is used for stratified elastic media and is implemented in intercept time-slowness domain.It consists of two steps with a combination of a global search method and a nonlinear local search method.In the first step，long wavelength parts of velocity structure are estimated by very fast simulating algorithm，when travelling time information is used.In the second step，short wavelength velocity variations are included into the model by conjugate gradient method，and waveform information in intercept time-slowness domain is used.Multi-grid parameter technique is used in our full waveform inversion.The method is applied to study a double BSR in the eastern Nankai accretionary wedge，and the results show that the full waveform inversion is an important tool for gas hydrate research.
Key words：Full waveform inversion；Gas hydrates；Bottom simulating reflector；Eastern Nankai trough
