1. SWPI,China University of Petroleum(East China),Qingdao 266580,China 2. Ocean Oil Production Plant,Shengli Oilfield Company,China Petrochemical Group Co. Ltd.,Dongying 257237,China
As a method for effectively reducing the velocity model of the subsurface media,seismic tomography provides a reliable initial velocity model for full waveform inversion.The finite frequency characteristics of seismic wave propagation are realized from primitive ray toe to phase shift travel time tomography and instantaneous travel time tomography.From the acoustic wave equation to the elastic wave equation and from the isotropic medium to the VTI,TTI media,the simulation of the real underground medium is realized.The morbidity of mitigating tomographic inversion has also been a research hotspot.The commonly used methods have regularization,and the sensitive nucleus of Gaussian beam tomography has replaced the traditional ray-sensitive nucleus.Furthermore,in order to avoid the dependence of the accuracy of the imaging results on the true depth of the reflection bits on the common imaging gather,the angular domain double differential reflection tomography can converge stably and efficiently to the accurate migration velocity model.At present,tomography is gradually transitional to anisotropic media,data used are transitional from VSP to WVSP,and a single waveform is developed into multiple waveforms combined inversion.However,problems related to resolution and computational efficiency still require attention.
Langan R T, Lerche I, Cutler R T . Tracing of rays through heterogeneous media:An accurate and efficient procedure[J]. Geophysics, 1985,50(9):1456-1465.
[2]
Virieux J, Farra V . Ray tracing in 3-D comples isotropic media:An analysis of the problem[J]. Geophysics, 1991,56(12):2057-2069.
[3]
Rawlinson N, Hauser J, Sambridge M . Seismic ray tracing and wavefront tracking in laterally heterogeneous media[J].Advances in Geophysics, 2008(49):203-273.
[4]
Um J, Thurber C . A fast algorithm for two-point ray tracing[J]. Bulletin of the Seismological Society of America, 1987,77(3):972-986.
[5]
Vinje V, Iversen E, Gjoystdal H . Traveltime and amplitude estimation using wavefront construction[J]. Geophysics, 1993,58(8):1157-1166.
[6]
Vinje V, Iversen E, Astebol K , et al. Estimation of multivalued arrivals in 3D models using wavefront construction-Part I[J]. Geophysical Processing, 1996,44(5):819-842.
[7]
Vinje V, Iversen E, Astebol K , et al. Estimation of multivalued arrivals in 3D models using wavefront construction-Part II[J]. Geophysical Processing, 1996,44(5):843-858.
[8]
Bai C Y, Greenhalgh S A . 3-D local earthquake hypocenter determination with an “irregular” shortest-path method[J]. Bull. Seism. Soc. Am., 2006,96(6):2257-2268.
[9]
Zhou B, Greenhalgh S A . “Shortest path” ray tracing for most general 2D/3D anisotropic media[J]. Journal of Geophysics & Engineering, 2005,2(1):54-63.
[10]
Bona A, Slawinski M A, Smith P . Ray tracing by simulated annealing: Bending method[J]. Geophysics, 2009,74(2):25-32.
[11]
Kirkpatrick S, Jr Gelatt C D, Vecchi M P .Optimization by simulated annealing[J]. Science, 1983,220(11):650-671.
[12]
Vidale J E . Finite difference traveltime calculation[J].Seismics, 1988(78):2062-2076.
[13]
Qin F H, Luo Y, Olsen K B , et al. Finite difference solution of the eikonal equation along expanding wavefronts[J]. Geophysics, 1992,57(3):478-487.
[14]
Vidale J E . Finite difference calculation of traveltimes in three dimensions[J]. Geophysics, 1990,55(3):521-526.
[15]
Hole J A, Zelt B C . 3-D finite difference reflection traveltimes[J]. Geophysical Journal International, 1995,121(2):427-434.
[16]
Afnimar, Koketsu K . Finite difference traveltime calculation for head waves travelling along an irregular interface[J].Geophysics, 2000(143):729-734.
[17]
Podvin P, Lecomte I . Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools[J]. Geophysical Journal International , 1991,105(1):271-284.
[18]
Van T J, Symes W W . Upwind finite difference calculation of traveltimes[J]. Geophysics, 1991,56(6):812-821.
[19]
Kim S, Cook R . 3-D traveltime computation using second order ENO scheme[J].Geophysics, 1999(64):1867-1876.
[20]
William A, Schneider J . Robust and efficient upwind finite difference traveltime calculatioins in three dimensions[J]. Geophysics, 1995,60(4):1108-1117.
[21]
Qian J L, Symes W W . A adaptive finite difference method for traveltimes and amplitudes[J]. Geophysics, 1999,64(1):167-176.
[22]
Popovici A M, Sethian J A . Three-dimensional travel-time computation using the Fast marching method[J]. Geophysics, 1999,64(3):516-523.
[23]
Popovici A M, Sethian J A . 3-D imaging using higher order fast marching traveltimes[J]. Geophysics, 2002,67(10):604-609.
[24]
Alkhalifah T, Fomel S . Implementing the fast marching eikonal solver:spherical Versus Cartesian coordinates[J]. Geophysical Prospecting, 2001,49(2):165-178.
Huang X G, Sun J G, Sun Z Q , et al. A method for calculating complex travel time of seismic waves based on L-BFGS theory[J]. Geophysical and Geochemical Exploration, 2016,40(5):961-967.
[26]
Chavent G . Identification of functional parameters in partial differential equations[J].Joint Automatic Control Conference, 1974(12):155-156.
[27]
Taillandier C, Noble M, Chauris H , et al. First-arrival traveltime tomography based on the adjoint-state method[J]. Geophysics, 2009,74(6):1-10.
[28]
Waheed U, Flagg G, Yarman C E . First-arrival traveltime tomography for anisotropic media using the adjoint-state method[J]. Geophysics, 2016,81(4):147-155.
Li Y D . A new preconditioned adjoint state method for travel time chromatography of first arrival waves[J]. Chinese Journal of Geophysics, 2017,60(10):3934-3941.
Gao F, Wang J Y, Sun C . Wireless tomography algorithm based on exponential moving average[J]. Journal of Beijing Institute of Technology, 2018,38(11):1141-1148.
[31]
Dahlen F A, Hung S H, Nolet G . Fréchet kernels for finite-frequency traveltimes—I.Theory[J]. Geophysical Journal International, 2000,141(1):157-174.
[32]
Spetzler J, Snieder R . The formation of caustics in two and three dimensional media[J]. Geophys. J. Int., 2001,144(1):175-182.
[33]
Van der Hilst R D, de Hoop M V . Banana_doughnut kernels and mantle tomography[J]. Geophys. J. Int., 2005,148(6):956-961.
[34]
Baig A, Dahlen F A . Statistics of traveltimes and amplitudes in random media[J]. Geophys. J. Int., 2015,158(1):187-210.
[35]
Zhao L, Jordan T H . Structural sensitivities of finite-frequency seismic waves:a full-wave approach[J]. Geophys. J. Int., 2006,165(3):981-990.
[36]
Zhang Z G, Shen Y, Zhao L . Finite-frequency sensitivity kernels for head waves[J]. Geophys. J. Int., 2007,171(2):847-856.
[37]
Sieminski A, Liu Q Y, Tranmpert J , et al. Finite-frequency sensitivity of surface waves to anisotropy based upon adjoint methods[J]. Geophys. J. Int., 2007,168(3):1153-1174.
Xu X M, Shi D N, Li X F . Research progress of finite frequency tomography[J]. Progress in Geophysics, 2009,24(2):432-438.
[39]
Liu Y Z, Wu Z, Geng Z C. First-arrival phase-traveltime tomography [C]//SEG 2017 Workshop:Full-waveform Inversion and Beyond, 2017: 83-86.
[40]
Liu Y Z, Wu Z. Which traveltime tomography should be used in near-surface velocity inversion [C]//SEG International Exposition and 88 th Annual Meeting , 2018: 2617-2620.
[41]
张晓辉 . 直达波与反射波联合层析成像方法研究[D]. 秦皇岛:燕山大学, 2017.
[41]
Zhang X H . Study on direct wave and reflected wave tomography[D]. Qinhuangdao:Yanshan University, 2017.
[42]
付翠 . 井间地震初至波与反射波旅行时联合层析反演[D]. 西安:长安大学, 2014.
[42]
Fu C . Combined tomographic inversion of the first and reflected waves traveling between Wells[D]. Xi’an:Chang’an University, 2014.
Yu D, Sun Y, Lu J , et al. Parallel tomographic algorithm and its application in ground fracture investigation during shallow initial wave travel[J]. Geophysical and Geochemical Exploration, 2017,41(5):977-985.
[44]
黄光南 . 基于射线理论的地震速度层析成像方法研究[D]. 北京:中国石油大学, 2013.
[44]
Huang G N . Research on seismic velocity tomography based on ray theory[D]. Peking:China University of Petroleum, 2013.
Zhang X M Translate . The p-wave and s-wave tomography of the subduction zone in Japan is inverted by using local seismic,remote seismic travel time and surface wave data[J]. World Earthquake Translation, 2019,50(1):35-63.
[46]
Bishop T N, Bube K P, Cutler R T , et al. Tomographic determination of velocity and depth in laterally varying media[J]. Geophysics, 2012,50(6):903-923.
[47]
查树贵 . 约束条件下的反射层析成像方法[J]. 江汉石油科技, 1994,4(3):39-44.
[47]
Cha S G . Reflection tomography under constraint conditions[J]. Jianghan petroleum technology, 1994,4(3):39-44.
[48]
Wang Y, Houseman G A . Tomographic inversion of reflection seismic amplitude data for velocity variation[J]. Geophysical Journal International, 1995,123(2):355-372.
[49]
Kosloff D D, Sudman Y . Uncertainty in determining interval velocities from surface reflection seismic data[J]. Geophysics, 2002,67(3):952-963.
Bai X, Li Z C, Zhang K , et al. High precision tomographic velocity inversion based on well data constraints[J]. Geophysical and Geochemical Exploration, 2015,39(4):805-811.
Ren F, Li Z C, Zhang M , et al. Migration velocity analysis by differential similarity optimization method based on common image point trace set in Angle domain[J]. Geophysical and Geochemical Exploration, 2017,41(5):872-880.
Zhang R Z, Gao C X . Robust time domain migration velocity analysis for establishing initial model in full waveform tomography[J]. Petroleum Geophysics, 2018,16(1):63-66.
Liang S X, Wang Q, Jiao Y J , et al. Analysis and Evaluation of LSQR Method in Position Field Inversion[J]. Geophysical and Geochemical Exploration, 2019,43(2):359-366.
[55]
Sava P, Fomel S. Coordinate-independent angle-gathers for wave equation migration [C]//SEG Expanded Abstracts, 2005: 1850-1853.
[56]
Tang Y, Biondi B . Target oriented wavefield tomography using synthesized Born data[J]. Geophysics, 2011,76(5):191-207.
Liu R H, Yin X Y, Pu Y T . Analysis of factors affecting PML absorption by wave equation[J]. Geophysical and Geochemical Exploration, 2016,40(4):757-762.
Xiao Y L, Fan X, Wang X T , et al. Application of grid tomography velocity inversion technique in depth migration of Qigu anticline prestack[J]. Oil Geophysical Prospecting, 2017,52(2):98-105.
Xue H, Du M, Wen P F , et al. Application of grid tomography velocity inversion method in quasi-three-dimensional Xisha hydrate[J]. Geophysical and Geochemical Exploration, 2017,41(5):846-851.
[60]
Nolan C J, Symes W W. Imaging and coherency in complex structure [C]//SEG Expanded Abstracts, 1996: 359-363.
[61]
Nolan C J, Symes W W . Global solution of a linearized inverse problem for the wave equation[J]. Communications in Partial Differential Equations, 1997,22(1):127-149.
[62]
Kroode A P E, Smit D J, Verdel A R . Amicrolocal analysis of migration[J]. Wave Motion, 1998,28(1):149-172.
[63]
Xu S, Chauris H, Lambare G, et al. Common angle image gather:A new strategy for imaging complex media [C]//SEG Expanded Abstracts, 1998: 1538-1541.
[64]
Prucha M L, Biondi B, Symes W. Angle domain common image gathers by Wave equation migration [C]//SEG Expanded Abstracts, 1999: 824-827.
[65]
Lehmann I . Publications du Bureau central seismologique international[M]. Bureau, 1936.
[66]
Lailly P. The seismic inverse problem as a sequence of before stack migrations:Conference on Inverse Scattering,Theory and Application [C]//Society for Industrial and Applied Mathematics,Expanded Abstracts, 1983: 206-220.
[67]
Gauthier O, Virieux J, Tarantola A . Two dimensional nonlinear inversion of seismic waveforms:Numerical results[J]. Geophysics, 1986,51(2):1387-1403.
[68]
Beylkin G . Imaging of discontinuities in the inverse scattering problem by inversion of a causal generalized Radon transform[J]. Journal of Mathematical Physics, 1985,26(1):99-108.
[69]
Bleistein N . On the imaging of reflectors in the earth[J]. Geophysics, 1987,52(4):931-942.
[70]
Tarantola A . Inverse problem theory:Method for data fitting and model parameter estimation[M]. Elsevier Science Publ. Co. Inc., 1987.
[71]
Mora P . Nonlinear 2-D elastic inversion of multi-offset seismic data[J]. Geophysics, 1987,52(7):1200-1228.
[72]
Cerveny V . Direct and inverse kinematic problems for inhomegeneous anisotropic media—linearization approach[J].Contrib. Geophys. Inst. Acad. Sci, 1982(13):127-133.
[73]
Hiraharak Ishikaway . Travel time inversion for three-dimensional P-wave velocity anisotropy[J]. Journal of Physics of the Earth, 1984,32(3):197-218.
[74]
Jech J., Pšenĉìk I . Three dimensional inversion problem for inhomogeneous transversely Isotropic media[J]. Studia Geophysica et Geodetica, 1988,32(3):136-143.
[75]
Jech J, Psencik I . First order perturbation method for anisotropic media[J]. Geophysical J. Int., 1989,99(2):369-376.
[76]
Chapman C H, Pratt R G . Traveltime tomography in anisotropic media-I. Theory[J]. Geophysical J. Int., 1992,109(1):1-19.
[77]
Zhou B, Greenhalgh S . Nonlinear traveltime inversion for 3D seismic tomography in strongly anisotropic media[J]. Geophysical J. Int., 2008,172(1):383-394.
[78]
黄国娇, 白超英, 钱卫 .一般 TI 介质中多震相走时层析成像——以井间成像为例[J]. 石油地球物理勘探, 2016,51(1):115-126.
[78]
Huang G J, Bai C Y, Qian W . Multi-seismic traveltime tomography in general TI media — taking interwell imaging as an example[J]. Oil Geophysical Prospecting, 2016,51(1):115-126.
He L Y . Cross-well anisotropic multi-wave traveltime tomography — taking wave field simulation picking and taking time as an example[J]. Progress in Geophysics, 2018,33(5):1796-1806.
Huang J, Yang G Q, Li Z C , et al. Numerical simulation of pseudo-acoustic wave equation in TTI medium[J]. Geophysical and Geochemical Exploration, 2018,42(1):134-143.
Jia Y Y, Dong M . Finite element model simulation of seismic wave anisotropic traveltime tomography[J]. Computer Simulation, 2019,36(3):216-219.
[82]
Hagedoorn J G . A process ofseismic reflection interpretation[J].Geophysical Processing, 1954(2):85-127.
[83]
Tikhonov A N, Arsenin V Y . Solution of Ⅲ-posed Problems[J].New York:John Wiley, 1977(25):11-23.
[84]
王宏伟 . 菲涅尔带地震层析[J]. 石油地球物理勘探, 1996,31(3):309-326.
[84]
Wang H W . Fresnel zone seismic tomography[J]. Oil Geophysical Prospecting, 1996,31(3):309-326.
[85]
Marquering H, Dahlen F A, Nolet G . Three dimensional sensitivity kernels for finite frequency traveltimes: the banana doughnut paradox[J]. Geophysical J. I., 1999,137(3):805-815.
[86]
Dahlen F A . Finite frequency sensitivity kernels for boundary topography Perturbations[J]. Geophysical J. Int., 2005,162(2):525-540.
[87]
Gelis C, Virieux J, Grandjean G . 2D elastic waveform inversion using Born and Rytov formulations in the frequency domain[J]. Geophysical J. Int., 2007,168(2):605-633.
[88]
Spetzler G, Snieder R . The effect of small scale heterogeneity on the arrival time of waves[J]. Geophysical J. Int., 2001,145(3):768-796.
[89]
Zhang Z G, Shen Y, Zhao L . Finite frequency sensitivity kernels for head waves[J]. Geophysical J. Int., 2007,171(3):847-856.
[90]
杨国辉 . 菲涅尔带旅行时层析成像方法及应用研究[M]. 厦门:厦门大学, 2009.
[90]
Yang G H . Study on the method and application of tone imaging in Fresnel zone[M]. Xiamen:Xiamen University, 2009.
[91]
Xia F, Jin S. Angle domain wave equation tomography using RTM image gathers [C]//SEG Technical Program Expanded Abstracts, 2013: 4822-4826.
[92]
Zhang K, Yin Z, Li Z C , et al. Wave equation tomographic velocity Inversion method based on the Born/Rytov approximation[J]. Apply Geophysics, 2013,10(3):314-322.
[93]
Xie X B, Geng Y . Gaussian beam based finite frequency turing wave tomography[J]. Journal of Applied Geophysics, 2014,109:71-79.
Wang G C . Application of multi-component tomographic velocity inversion method based on three-dimensional angular gathers in Shengli area[J]. Progress in Geophysics, 2017,32(4):1659-1664.
Yuan M L, Jiang F Y, Yang H F , et al. Research on Gaussian beam linear forward modeling method[J]. Geophysical and Geochemical Exploration, 2017,41(5):881-889.
Song Z D, Jiang Z H, Wei C W , et al. Research on forward modeling of seismic tomography based on three-dimensional layered medium model[J]. Computerized Tomography Theory and Applications, 2018,27(1):1-8.
Cheng G . Large-scale sparse matrix compression storage and solution in reflection seismic traveltime tomography[J]. Progress in Geophysics, 2008,23(3):676-678.