50 4 2007 7 CHINESE JOURNAL OF GEOPHYSICS Vol. 50, No. 4 Jul., 2007,,..,2007,50 (4) :1174 1181 Xu T,Ning J R,Liu C C,et al. Influence of the self2organization of the Earth interior upon the traveltime and amplitude of seismic wave. Chinese J. Geophys. (in Chinese),2007,50 (4) :1174 1181 1, 1,2, 1,3, 1,4 1,, 100029 2, 100083 3, 300452 4, 100081 ;,, von Karman ; ;,, ;, ; von Karman,,.,,,,, 0001-5733(2007)04-1174 - 08 P315 2006-06 - 15, 2007-03 - 21 Influence of the self2organization of the Earth interior upon the traveltime and amplitude of seismic wave XU Tao 1,NINGJun2Rui 1,2,LIU Chun2Cheng 1,3,LI Shou2Lin 1,4 1 State Key Laboratory of Lithosphere Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China 2 Exploration and Production Research Institute, SINOPEC, Beijing 100083, China 3 China National Offshore Oil Corporation. Tian Jin, Tianjin 300452, China 4 International Exploration Corporation, China Petroleum & Chemical Corporation, Beijing 100081, China Abstract We utilize random media to describe small2scale inhomogeneity and self2organized structure in the Earth interior. A block model is adopted to describe the self2organized structure instead of traditional grids in 22D. Gaussian, exponential, and von Karman autocorrelation functions are used to represent the isotropic and anisotropic self2organized structure. We use numerical methods to perform ray tracing and rectify amplitudes by dynamic ray tracing, furthermore, analyze the influence upon the kinematic and dynamic characteristics of seismic waves. The results indicate that seismic ray trajectories probably generate significant aberrances and the amplitudes of seismic rays may weaken or enlarge at different offsets for the existence of the self2organized structure in the earth interior ; from Gaussian to von Karman, amplitudes present more obvious statistic characters. Keywords Self2organized structure, Block model, Seismic wave, Ray trajectory, Traveltime, Seismic amplitude (40234044,40374010,40404009) 973 (2002CB412604).,,1978,,2004,,. E2mail : xutao @mail. iggcas. ac. cn
4 : 1175 1 20 60,,,, [1 8 ].,, ( ),. : (1),,,, ; (2) ; (3) (, ).,.,,,,,., [9 11 ].,..,, [12,13 ].,,,.,, [14 21 ] [9 12,22 29 ],,,. ;,, ; ; von Karman. 2 211 [9 14,16,17,22,23,30 ]., v ( x) = v 0 ( x) + v ( x), (1) v 0, ; v,. v ( x) = v 0 ( x) (1 + ), (2), = ( x), ( ) ( ). ( x) = 0, (3) 2 ( x) = 2, (4) ( x) ( x ) = R ( r), (5), 2, R ( r), r = x - x. 212 ( ).,, von Karman ( ) Kummer [9 14,16,17,22,23,30 ], [11,30 ]., ( x) R ( r) ^F ( k) [31, 32 ]. 1 2., K 0, 2 a b x z. ^f ( k),. 1 Table 1 Isotropic autocorrelation functions and associated power spectrums ( ) e r2 Πa 2 e - a2 k 2 Π4 ( ) e - a2 k 2 r Π4 e rπa ( a - 2 + k 2 ) - 1 ( a - 2 + k 2 r) - 3Π2 von Karman K 0 ( rπa) ( a - 2 + k 2 ) - 1Π2 ( a - 2 + k 2 r) - 1
1176 (Chinese J. Geophys. ) 50 2 Table 2 Elliptical anisotropic autocorrelation functions and associated power spectrums ( ) ( ) e x2 Πa 2 + z 2 Πb 2 e - ( a2 k 2 x + b2 k 2 z )Π4 e x 2 Πa 2 + z 2 Πb 2 (1 + a 2 k 2 x + b 2 k 2 z) - 3Π2 von Karman K 0 ( x 2 Πa 2 + z 2 Πb 2 ) (1 + a 2 k 2 x + b 2 k 2 z) - 1 ^F ( k) ^f ( k) : : ^F ( k) = ^f 3 ( k) ^f ( k) ; (6) ^f ( k) = ^f ( k) ; (7), k = ( k T k) 1Π2 ;, k = ( k T Lk) 1Π2, L. [30 ] : ^f ( k) = ^f S ( k) ^f a ( k) ^f ag ( k), (8),^f S ( k), ^f a ( k), ^f S ( k) = k - dπ2- N, (9) ^f a ( k) = (1 + ( ak) - 2 ) - dπ4- NΠ2, (10) ^f ag ( k) ^f ag ( k) = e - a2 G k2 Π8. (11) (9) (11), d ; N Hurst ; a a G von Karman ;.,.,N = - dπ2, ; a G = 0 von Karman, N = 0, N = 1Π2 ; a G = 0, a =, [30]. 1 ^f a ( k) ^f ag ( k), a, k a G, k. 1 ^f a ( k) ^f ag ( k). (10) (11). Fig. 1 Product of high2pass filter ^f a ( k) and low2pass Gaussian filter ^f ag ( k) against a logarithmic wave2number k, : (1) W ( x), ^W ( k) ; (2), ^U ( k) = ^f ( k) ^W ( k) ; (3) ^U ( k), U ( x) ; (4) U ( x) ;. 2 (2a) (2b) von Karman ( ) (2c). 20 km,.,, von Karman,. 213., :,,,,,.,. [15,18 ] ( [ NlgN ], N ).,.,, [19,20 ]. [18,21 ],,,,. ( ) [18 ]. 3, 7, ( 5a), ( 3 4 )., ( 3 4) : (1) 4,, ; (2) ; (3),,, 4,,, ( 5b 5d).
4 : 1177 4 1 new Fig. 4 New incidence angle new updated by an initial incidence angle 1 and an increment 3 2 (a) ; (b) ; (c) von Karman ( ). ( ), ( ). Fig. 2 Three self2organized media (a) Gassian ; (b) Exponential ; (c) Zero von Karman. Blue colors denote fast regions and red colors denote slow regions. 3, 4 Fig. 3 The representation of a self2organized media in a complex block structure. Block # 4 denotes a Gaussian random media 311, [14,16,17 ]., [33 ] : 9 t x = vsin, 9 t z = vcos, 9 t = - cos 9 v 9 x + sin 9 v 9 z. (12), Eluer Runge2Kutta., [21 ] [18 ]. ( 4 ),,,. new 312 = 1 + ( x rcvr - x 1 ) x. (13),. : ; ( ), ; [33,34 ]. 4,.,
1178 (Chinese J. Geophys. ) 50 5 (a) ; (b) ; (c) ; (d) ; (e) c, ; (f) d, ; (g) c,von Karman ; (h) d,von Karman ; (i) (d,f,h) (b) Fig. 5 (a) Initial shooting trajectories based on the background velocity model ; (b) Tracing results based on the background velocity model ; (c) Initial shooting trajectories based on a superimposed Gaussian random perturbation velocity model ; (d) Tracing results based on the velocity model as (c) ; (e) Same as (c) with an exponential random perturbation ; (f) Tracing results based on (e) ; (g) Same as (c) with a von Karman random perturbation ; (h) Tracing results based on (g) ; (i) Traveltime error of tracing trajectories in (d), (f) and (h) to background traveltime (b)
4 : 1179 [35 37 ], 5, 400 km 90 km ; 7., 4,. : ( 5a 5b) ( 5c 5d) ( 5e 5f) von Karman ( 5g 5h)., : 2 km, 8 km,, 1 %( 3 %) ;, ( 5c 5h), von Karman,. 5a 5b ;5c 5h (, ) ; 5i ( 5d 5f 5h) ( 5b). 200 km,,, 100 150 km 250 300 km,,.,, von Karman,, von Karman,. 6a, ; 6b,., 200 250 km, 250 300 km. von Karman,,,,,. 6 (a) (5b) (5d) (5f) von Karman (5h) ; (b) Fig. 6 (a) From the top down, the synthetic seismograms of dynamic tracing based on background (5b), Gaussian (5d), Exponential (5f) and von Karman (5h) velocity models ; (b) Three random amplitude relative errors to the background
1180 (Chinese J. Geophys. ) 50 5,,,.,,.,, ;,.,,,.,,,,,.,. (References) [ 1 ],,.. :,2003 Teng J W, Zhang Z J, Bai W M, et al. Lithosphere Physics (in Chinese). Beijing : Science Press, 2003 [ 2 ],,..,2001, 16(5) : 634 642 Chen Y T, Teng J W, Zhang Z J. Geophysics : The 20 th century in retrospect and the 21 st century in prospect. Advance in Earth Sciences (in Chinese), 2001, 16 (5) : 634 642 [ 3 ] Zhao W, Nelson K D, the Project INDEPTH Team. Deep seismic reflection evidence for continental underthrusting beneath south Tibet. Nature, 1993, 366 : 557 559 [ 4 ]..,1974, 17 (2) : 121 134 Teng J W, et al. Deep reflected waves and the structure of the earth crust of the eastern part of Qaidam Basin. Chinese J. Geophys. (in Chinese), 1974, 17 (2) : 121 134 [ 5 ],,. -.,1995,38 ( ) : 3 14 Gao R, Cheng X Z, Ding Q. Preliminary geodynamic model of Golmud2Ejin Qi geoscience transect. Chinese J. Geophys. ( in Chinese), 1995, 38 (Suppl. ) : 3 14 [ 6 ],,. -.,2003,25(3) : 229 241 Wang C Y, Han W B, Wu J P, et al. Crust structure beneath the Songpan2Garze orogenic belt. Acta Seismologica Sinica (in Chinese), 2003, 25 (3) : 229 241 [ 7 ]. -. :. :, 19901 1 14 Gao E Y. Crustal structure in northern Himalaya2Yarlu Zangbu by artificial deep seismic explosion. In : Tibet Geophysical Corpus. Beijing : Geological Publishing House, 19901 1 14 [ 8 ].. :,1988 Department of Scientific Programming and Earthquake Monitoring, China Earthquake Administration. Developments in the Research of Deep Structures of Chinaπs Continent. Beijing : Seismological Press, 1988 [ 9 ] Levander A R,Holliger K. Small2scale heterogeneity and large2scale velocity structure of the continental crust. J. Geophys. Res., 1992, 97(B6) : 8797 8804 [10 ] Baig A M, Dahlen F A, Hung S H. Traveltimes of waves in three2 dimensional random media. Geophys. J. Int., 2003, 153 : 467 482 [11 ] Ikelle L T, Yung S K, Daubl F. 2 - D random media with ellipsoidal autocorrelation function. Geophysics, 1993, 58(9) : 1359 1372 [12 ],..,2001,36(5) : 546 552 Xi X, Yao Y. 2 - D random media and wave equation forward modeling. OGP (in Chinese), 2001, 36(5) : 546 552 [13 ],..,2005, 40(1) : 71 75 Xi X, Yao Y. Non2stationary random medium model. Chinese), 2005, 40 (1) : 71 75 OGP ( in [14 ] Ojo S B, Mereu R F. The effect of random velocity functions on the travel times and amplitudes of seismic waves. Geophys. J. R. astr. Soc., 1986, 84 : 607 618 [15 ] Moser T J. Shortest path calculation of seismic rays. Geophysics, 1991, 56 : 59 67 [16 ] Crossley D J, Jensen O G. Fractal velocity models in refraction seismology. Pure Appl. Geophys., 1989,131 : 61 76 [17 ] Witte O, Roth M, Muller G. Ray tracing in random media. Geophys. J. Int., 1996, 124 : 159 169 [18 ],,..,2004, 47(6) :1118 1126 Xu T, Xu G M, Gao E G, et al. Block modeling and shooting ray tracing in complex 32D media. Chinese J. Geophys. (in Chinese), 2004, 47 (6) : 1118 1126 [19 ] Zhao A H, Zhang ZJ, Teng J W. Minimum travel time tree algorithm for seismic ray tracing : improvement in efficiency. Journal of Geophysics and Engineering, 2004, 1(4) : 245 251 [20 ],..,2005, 48 (5) : 1141 11471 Zhao A H, Ding Z F. A double2grid algrithm for calculating traveltimes of wide2angle reflection waves. Chinese J. Geophys. (in
4 : 1181 Chinese), 2005, 48 (5) : 1141 1147 [21 ] Xu T, Xu G, Gao E, et al. Block modeling and segmentally iterative ray tracing in complex 3D media. Geophysics, 2006, 71 (3) : T41 T51 [22 ] Frankel A, Clayton R W. Finite difference simulations of seismic scattering : implications for the propagation of short2period seismic waves in the crust and models of crustal heterogeneity. J. Geophys. Res., 1986, 91(B6) : 6465 6489 [23 ] Korn M. Seismic wave in random media. J. Applied Geophys., 1993, 29 : 247 269 [24 ] Zhang Z J, He Q D, Teng J W, et al. Simulation of 32component seismic records in a 22dimensional transversely isotropic media with finite2difference. Can J. Expl. Geophys., 1993, 29 : 51 58 [25 ] Zhang Z J, Wang G J, Harris T M. Multi2component wavefield simulation in viscous extensively dilatancy anisotropic media. Earth Planet Inter., 1999, 114 : 25 38 Phys. [26 ] Yang D H, Liu E R, Zhang Z J, Teng J W. Finite2difference modeling in two2dimensional anisotropic media using a flex2corrected transport technique. Geophys. J. Int., 2002, 148 : 320 328 [27 ]..,2002, 45 (4) : 575 583 Yang D H. Finite element method of the elastic wave equation and wavefield simulation in two2phase anisotropic media. Chinese J. Geophys. (in Chinese), 2002, 45(4) : 575 583 [28 ],,..,2007, 50(1) : 221 232 Liu Y X, Xu T, Zhao B, et al. Seismic sounding of anisotropic self2 similar self2organized medium. Chinese J. Geophys. (in Chinese), 2007, 50 (1) : 221 232 [29 ],,.., 2007, 50(3) : 830 836 Liu Y X, Xu C M, Ning J R. Comparison of acoustic propagation in several different types self2organized media. Chinese J. Geophys. (in Chinese), 2007, 50 (3) :830 836 [30 ] Klimes L. Correlation functions of random media. Pure Appl. Geophys., 2002, 159 : 1811 1831 [31 ],,..,2003 Zhu H, Huang H N, Li Y Q, et al. Random Signal Analysis (in Chinese). Beijing : Beijing Institute of Technology Press, 2003 [32 ] Goff J A, Jordan T H. Stochastic modeling of seafloor morphology : inversion of sea beam data for second2order statistic. J. Geophys. Res., 93B : 13589 13608 [33 ] Cerveny V. Seismic Ray Theory. Cambridge University Press, 2001 [34 ]. [ Xu T. ]. :,2003 Block modeling and fast ray tracing in complex three2 dimensional media [ Ph. D. thesis] (in Chinese). Hefei : University of Science and Technology of China, 2003 [35 ] Zhang Z, Li Y K, Wang GJ, et al. E - W crustal structure under the northern Tibet and down2bowing Moho under the northern Tibet revealed by wide2angle seismic profile. Science in China (Series D). 2002, 31 (11) : 881 888 [36 ] Zhang Z, Teng J W, Li Y K, et al. Crustal structure of seismic velocity in southern Tibet and east2westward escape of the crustal material. Science in China (Series D), 2004, 47 (6) : 500 506 [37 ] Zhang Z, Klemperer S L. West - east variation in crustal thickness in northern Lhasa block, central Tibet, from deep seismic sounding data. J. Geophys. Res., 2005, 110, doi : 1011029Π2004JB003139 ( )