P-SV Wave-Field Connection Technique for Regional Wave Propagation Simulation
|
|
- Sydney James
- 5 years ago
- Views:
Transcription
1 Bulletin of the Seimological Society of America, Vol. 95, No. 4, pp , Augut 2005, doi: / P-SV Wave-Field Connection Technique for Regional Wave Propagation Simulation by Zengxi Ge, Li-Yun Fu,* and Ru-Shan Wu Abtract A boundary element (BE) method i developed to calculate the twodimenional P-SV elatic repone for crutal wave guide with irregular topographic feature. To imulate long-range propagation of regional wave, a connection technique i propoed to avoid large matrix inverion that become formidable for longrange, high-frequency problem. By uing thi technique, a long crutal wave guide can be divided into relatively horter ection, and the BE method can be ued ection by ection to model the effect of rough topography on wave propagation at extended regional ditance. The validity of the technique i teted by comparion with a direct calculation. Numerical imulation with thi cheme how that rough topography can catter the P and Rayleigh wave and attenuate the energy propagating in the wave guide. Thi method can be ued in computing the ite effect on ite uch a canyon, mountain, and valley. The connection technique expand thi method to deal with large earth model with irregular topography. Introduction Regional phae have long been recognized a an important iue in the tudy of large-cale crutal tructure, mall-cale crutal heterogeneitie, eimic ource, and analyi of underground exploion and earthquake. To improve path correction and verify empirical obervation, numerical imulation are needed to etimate the path effect on regional wave propagation. A number of regional numerical modeling method have been developed to model Lg propagation behavior and path effect, depending on the complexity of the crut heterogeneity in the region conidered. Finite-difference (FD) and finite-element (FE) method are univeral numerical technique for imulation in general heterogeneou media. Becaue of the computational intenity, thee full-waveform method are often limited to low frequencie for Lg imulation at regional ditance. To reduce computation time, ome alternative and flexible approache, for example, ray diagram (Botock and Kennett, 1990; Keer et al., 1996) and dynamic raytracing method (Gibon and Campillo, 1994), have been ued to obtain intuitive undertanding of path effect on Lg propagation. For long-range Lg-wave propagation, Wu et al. (1996, 2000) introduced a half-pace GSP (generalized creen propagator) for modeling the main characteritic of Lg (2D SH cae) in heterogeneou crutal wave guide. The method wa ex- *Preent addre: Intitute of Geology and Geophyic, Chinee Academy of Science, P.O. Box 9825, Beijing , People Republic of China. lfu@mail.igca.ac.cn tended to include the cae of irregular urface topographie (Wu and Wu, 2001) Baed on obervation and numerical experiment, large-cale crutal tructure with variation in free urface and Moho topography control the principal characteritic of Lg propagation. In contrat to the FD and FE method, the boundary integral method are more uitable for modeling complex reflection and tranmiion effect acro largecale crutal tructure with rugged topography. For intance, the boundary element (BE) method provide a geometrically accurate decription of irregular interface. Becaue the BE method i formulated in term of integral along boundarie, the traction-free condition for a rugged free urface i eaily and naturally treated. However, for regional wave guide up to thouand of kilometer, the BE method in it original form lead to very large matrice to be inverted and it implementation i prohibitively expenive. In a previou paper (Fu and Wu, 2001), a SH-wavefield connection technique wa developed by which the BE method can be ued ection by ection for an event at farregional ditance up to everal thouand kilometer; the method take the output of the previou ection a the input of the next ection to complete the entire crutal wave-guide computation. Full-wave BE modeling i implemented within each ection. The wave-field connection technique i ued to couple the field calculated in two adjoining ection. The diviion of ection i baed on the complexity of crutal tructure with a criterion of minimizing the poible multiple backcattering between ection. Thi approach for 1375
2 1376 Z. Ge, L.-Y. Fu, and R.-S. Wu long regional waveguide lead to ignificant computational aving in time and memory compared with the whole waveguide BE method. The ection-by-ection approach alo lead to a hybrid modeling cheme of BE and GSP. The BE method i implemented in the frequency domain and ha a kernel function compatible with the GSP method. In the hybrid cheme, the time-conuming BE method can be ued to handle the ection with complicated boundary tructure and evere urface topographie. Subequently, the output will be ued a the input to the GSP method for modeling ection with a large volume of moderately heterogeneou media and mild topographie. The hybrid method ha been applied to two crutal wave-guide model from the Tibet region, one with Lg blockage and another without blockage (Fu and Wu, 2001). The object of thi article i to develop a P-SV wavefield connection technique for imulating elatic wave propagation in regional crutal waveguide. We firt introduce an integral equation repreentation for the crutal waveguide problem. We then briefly decribe the principle of the elatic BE method and validate the computation program uing previouly publihed reult. We develop the P-SV wave-field connection technique and tet it by uing numerical experiment. Finally, we apply the ection-by-ection approach to Lg-propagation imulation in long regional wave guide. Boundary Integral Equation for Crutal Wave Guide Conider 2D teady-tate elatic wave propagation in a implified crutal wave guide X 1 bounded by a rough free urface C 1 and an irregular Moho interface C 2. Figure 1 depict the geometry of the problem. The wave-guide medium i iotropic and homogeneou, decribed by the Lame contant (k 1 and l 1 ) and denity (q 1 ). The diplacement vector u(r) at a location r(x, z) atifie the following elatic wave equation: 2 2 l u(r) (k l) u(r) qx u(r) f(r, x) (1) where f(r, x) i the body force occupying a region X. u(r) alo atifie the traction-free boundary condition on C 1 and the continuity condition of diplacement and traction acro the Moho. The medium in the mantle X 2 i decribed by k 2, l 2, and q 2. We add two artificial boundarie C at the two truncated edge of the wave guide to form a cloed olution domain C C 1 C 2 C. Baed on the repreentation theorem (Aki and Richard, 1980), equation (1) can be tranformed into the following integral formulation: C(r)u(r) [u(r)r(r, r) t(r)g(r, r)]dc(r) C f(r, x)g(r, r)dx(r) (2) X Figure 1. Geometry of a implified crutal wave guide with irregular topography and Moho interface. where t(r) i the traction vector, the coefficient C(r) generally depend on the local geometry at r, and G(r, r) and R(r, r) are the fundamental olution (Green tenor) for diplacement and traction, repectively. C(r), R(r, r) and G(r, r) are 2 2 matrice for 2D problem. In thi equation, r i the poition of obervation point and r i the poition of cattering point. For implicity, the ource ditribution conit of a point ource at r 0 located inide X 1. The ource integral over X in the right ide of equation (2) can be reduced to G(r,r)f(r,x)dr G(r,r )f(x) (3) 0 X where f (x) i the ource pectrum vector. The diplacement Green tenor G(r, r) atifie: 2 [(k l) l ]G(r,r) 2 qx G(r,r) d(r r) (4) with the olution given by More and Fehbach (1953). The traction Green tenor R(r, r) can be derived from G(r, r) by uing Hooke law, R(r,r) ki( G(r,r)) l(g(r,r) G(r,r)), (5) where I i the unit dyadic and G(r, r)d i the tranpoe of DG(r, r) with repect to the correponding coordinate, ee Wu (1989, Appendix). For analytic expreion of G and R in an iotropic homogeneou elatic medium, thee Green tenor for 2D problem can be expreed a (Fu, 1996): i Gij (wdij ci c) j 4l i w ij (dijck nk cj n i) 4 r r 2 (c n 2c ccn ) 2 cccn r r k w c i n l r r r j i j i j k k i y k k
3 P-SV Wave-Field Connection Technique for Regional Wave Propagation Simulation 1377 where i, j 1, 2, n i the outward normal of the boundary C, and (1) 1 (1) b (1) w H 0 (kr) H 1 (kr) H 1 (kpr) kr 2 b (1) (1) 2 2 p 2 H (k r) H (k r) with k x/b, k p x/, r r r and c j (r j r j )/r. The boundary integral repreentation for wave propagation naturally atifie the Sommerfeld radiation boundary condition that are impoed on the far-field behavior at infinity. No wave come back to X 1 through C, that i, the following integral on C for the interior problem vanihe: C [R(r,r)u(r) G(r,r)t(r,r)]dr 0, r X (6) 1 For numerical calculation, truncating the wave guide i neceary and an infinite boundary element aborbing boundary technique (Fu and Wu, 2000) ha to be applied to the element at the truncating point on C 1 and C 2. Conidering equation (3) and (6), and applying the traction-free condition to the free urface C 1, equation (2) for the interior problem in the crut (region X 1 ) i implified to C(r)u(r) C1 R(r, r)u(r)dr [R(r, r)u(r) G(r, r)t(r)]dr (7) C2 G(r, r )f(x) 0 Equation (7) i a tarting point for numerical implementation for wave-propagation imulation. To olve u(r) and t(r) on C 2 we mut build the correponding boundary integral equation in X 2. The following integral formulation can be etablihed for r X 2 bounded by a cloed urface C C 2 C : C(r)u(r) [R(r,r)u(r) G(r,r)t(r)]dr 0 (8) C Boundary Element Method for Elatic Wave Simulation The boundary integral repreentation decribed previouly can be ued to calculate the field at any point inide X 1 once u(r) and t(r) are known on the boundarie. We ue the BE method to olve the value of u(r) and t(r) onthe boundarie. We dicretize C 1 into L 1 element and N 1 node, and C 2 into L 2 element and N 2 node. The total node number i N. By uing the linear interpolation hape function (n) in an element between the node I 1 and I 2, the variable u(r) and t(r) are approximated by the linear combination of their node value over the element, for example (ee Fig 1): I 2 l l li 1 u(n) u(r ) (n), (10) where n and l denote the local coordinate and local node index of an element. Uing the following Kronecker delta function notation relating the local node code l of an element to the global node code j we have 0, l j d, (11) ij 1, l j N I 2 l l ij j1 li 1 u(n) u(r ) (n)d, (12) Letting C i C(r i ) and u i u(r i ), equation (7) for i 1 to N i dicretized into Cu N L I 1 2 i i (r i,r(n)) l(n)dr(n) dijuj j1 e1 li1 C e L2 I2 (r i,r(n)) l(n)dr(n) dljuj e1 li1 C e L2 I2 G(r i,r(n)) l(n)dr(n) dljtj e1 li1 C e Similarly, the integration over the tranparent artificial boundary C vanihe for the interior problem, implifying equation (8) a G(r,r )f(x) i 0 which can be further rewritten a (13) C(r)u(r) [R(r,r)u(r) G(r,r)t(r)]dr 0. (9) C2 Equation (8) and (9) provide a decription of the field through the crutal wave guide, making poible the imultaneou evaluation of the unknown [u(r) onc 1 and u(r) and t(r) onc 2 ] by uing the continuity of diplacement and traction acro C 2. N (1) (1) (2) (2) (2) (2) ij j ij j ij j i 0 j1 with h u h u g t G(r,r )f(x), (14) L1 I2 (1) (1) ij j l ij j ij e1 li1 Ce h [R(r,r(n)) (n)dr(n)]d C d, (15)
4 1378 Z. Ge, L.-Y. Fu, and R.-S. Wu L2 I2 (2) (2) ij i l ij j ij e1 li1 Ce h [R(r,r(n)) (n)dr(n)]d C d, (16) L2 I2 (2) gij [G(r i,r(n)) l(n)dr(n)]d ij, (17) e1 li1 Ce (1) In thi expreion u j i the diplacement vector on the free urface C 1 and u (2) j and tj (2) are the diplacement and traction vector on the Moho C 2. The coefficient matrice h (1) ij, hij (2) (2) and g ij, obtained by numerically integrating the product of the Green tenor with interpolation hape function over element, denote a concentrated force generated at the jth cattering point on C and applied at the ith obervation point. For i 1toN, equation (14) can be further compacted a a matrix equation: (1) (1) (2) (2) (2) (2) H u H u G t F, (18) where F qg (r i, r 0 )f(x). Similarly, equation (9) for i 1 to N can be dicretized and compacted a (2) (2) (2) (2) H u G t 0, (19) where the coefficient matrice, H (2) and G (2), are calculated by uing the medium propertie of the mantle X 2. A imultaneou ytem of matrix equation (18) and (19) can be aembled by uing the continuity of diplacement and traction acro C 2, The computation program of the elatic wave BE method decribed previouly i teted by dimenionle frequency repone of a emicircular canyon of radiu a. Previouly publihed reult (Sánchez-Sema et al., 1985; Sánchez-Sema and Campillo, 1991) for thi typical topographic tructure are ued for comparion. The two harp edge at x a provide a crucial target to validate variou numerical method. Figure 2 how the comparion with good agreement in both horizontal and vertical amplitude between our reult (olid line) and the Sánchez-Sema and Campillo olution for vertically incident plane P wave and for variou normalized frequencie g defined a g 2a / l, where l i the wavelength. Poion ratio i aumed to be 1/3. We can ee ome minor departure becaue of different element approximation ued by thee two numerical method. Fewer element per wavelength will reduce the ize of the reultant coefficient matrice. To determine an applicable element number per wavelength, the comparion for 0 Incidence and g 1.0 are given in Figure 3 for different dicretization rate of point per wavelength to dicretize the arc of the emicircular canyon. We ee that a ampling at three point per wavelength could be ued for general application. The preceding comparion confirm the validity of our formulation and computation code. Later we will ue thi method to validate the connection technique for a larger model. We now demontrate the applicability of our program by yntheizing wave propagation through a ingle- u (2) (2) ũ (2) (2) t t (20) by which the unknown [u(r)onc 1 and u(r) and t(r)onc 2 ] can be found. Thee matrice in equation (18) and (19) are full with complex coefficient that are function of frequency, material property, and geometry. The BE method decribed previouly can be directly extended to complex geological tructure with multiple region (Fu, 1996) for explorationoriented eimic modeling. Alvarez-Rubio et al. (2004) applied the BE method to ite effect aement of laterally varying layered media. Becaue a large number of matrix operation are involved and the matrix for each frequency component mut be inverted, the BE method i computationally intenive at high frequencie for far-regional wave guide. To improve computation peed, a frequency-dependant element dimenion technique can be adopted in the program implementation (Fu and Mu, 1994). Bouchon et al. (1995) uggeted a threhold criterion approach to make the coefficient matrice parer by removing very mall entrie in the BE coefficient matrice. An efficient modification to the BE method can be made uing the ection-by-ection approach with the wave-field connection technique developed in thi tudy. Figure 2. The horizontal and vertical amplitude repone by our method (dotted line) and Sanchez- Sema and Campillo (1991) (olid line) of a emicircular canyon topography to vertically incident P wave for variou dimenionle frequencie g.
5 P-SV Wave-Field Connection Technique for Regional Wave Propagation Simulation 1379 Figure 3. The horizontal and vertical amplitude repone of a emicircular canyon topography to vertically incident P wave for g 1.0 with the olid line denoting the olution with 21 point per wave length and the dot denoting the pecified ampling rate. Figure 4. A ingle-layer flat crutal wave guide with an exploive ource at 2 km depth. The wave generated in the wave guide are grouped into three ytem: the firt directly from the ource to the receiver (olid line), the econd off the free urface (dahed line), and the third off the Moho (dotted line). layer crutal wave guide with flat free urface. The homogeneou wave guide hown in Figure 4 i 100 km long and 32 km thick, overlaying a flat mantle half-pace. The point ource (P-wave exploive ource) i located on the left boundary at 2.0 km depth. Figure 5 how the ynthetic eimogram calculated in the frequency range of Hz, with receiver at 1-km pacing along a vertical profile at the ditance of 80 km from the ource. We ee that for elatic wave propagation thi imple wave guide with both the flat topography and flat Moho lead to complex uperpoition of variou wave, demontrating the development of converted wave and the formation of guided wave a repetitive reflection at both the free urface and Moho. In term of propagation path (ee Fig. 4), we can identify three major ytem of wave generated in the wave guide. A hown in Figure 5, the direct P wave from the ource to receiver carrie a major part of the ource energy in thi flat boundary wave guide. The P-wave incident on the free urface and the Moho at oblique angle generate two ytem of wave, one et off the free urface and another et off the Moho. Thee two ytem of wave and their converted wave bounce back and forth between the free urface and the Moho. Becaue no cattering mechanim i preent in the wave guide, the contructive interference of the repeatedly reflected wave preent a checkerboard-like pattern (Jih, 1996) that adequately explain the formation of crutal guided wave either a multiple reflection or a higher mode.
6 1380 Z. Ge, L.-Y. Fu, and R.-S. Wu Figure 5. Synthetic eimogram for receiver at 1-km pacing along a vertical profile at 80 km for a point ource. The horizontal (left panel) and vertical (right panel) component are computed in the frequency range 0 3.5Hz. P-SV Wave-Field Connection Technique We aim to develop a ection-by-ection approach for imulating wave propagation in regional wave guide to reduce the computation cot of extremely large matrix operation. The wave-field connection technique couple the field between two adjacent ection. The connection configuration i illutrated in Figure 6a. An artificial boundary C AB i introduced a a wave-field connection boundary to a crutal wave guide coniting of an irregular free urface C 1 and an interface C 2. The wave guide i divided into four ubdomain, X 1 and X 2 in the crut, and Xl and X2 in the mantle. The left boundarie of the firt ection X 1 and X l, and the right boundarie of the econd ection X 2 and X 2 are aumed to extend to infinity at left and right, repectively. The field in X 1 and X l are calculated from the ource uing boundary element method and the output wave field u 0 (r) on the connection boundary C AB will be ued to atify the boundary condition acro C AB when the BE method i ued to calculate wave propagation in X 2 and X 2. The output field are received along the next connection interface C CD and will be ued a the input to the next propagation. With the initial field u 0 (r) known on the connection boundary C AB, we analyze wave propagation in X 2 and X 2 by building boundary integral equation for wave propagation in thee ection. Apparently, all the pace point in X 2 and X 2 except thoe on C AB only receive the cattered field. Therefore, the total field u(r) on the boundarie of X 2 and X 2 i compoed of u 0(r) u (r) r CAB u(r). (21) u (r) r C Applying the traction-free condition to the free urface C 1, the boundary integral equation can be built for r X 2 and r C AB AB C(r)u (r) C 1 C 2 CAB R(r,r)u (r)dr (r,r)[u (r) u (r)]dr (22) 0 CCD (r,r)u (r)dr The artificial boundarie C AB and C CD are aumed to be tranparent, implying the outward radiation of energy acro C AB and C CD alway hould be in the outward direction with no reflection returning to X 2. Therefore, there i no energy contribution cattering from C AB and C CD, that i, [R(r,r)u (r) G(r,r)t (r)]dr 0 r X (23) 2 CAB and [R(r,r)u (r) G(r,r)t (r)]dr 0 r X (24) 2 CCD Scattering from the artificial truncated point, A, B, C, and D, can be handled by uing an infinite-element aborbing boundary technique (Fu and Wu, 2000). Therefore, equation (22) i reduced to C(r)u (r) C1 R(r,r)u (r)dr [(r,r)[u (r) G(r,r)t (r)]dr (25) C2 CAB [(G(r,r)t (r) (r,r)u (r)]dr 0 0
7 P-SV Wave-Field Connection Technique for Regional Wave Propagation Simulation 1381 Figure 6. The boundary connection technique. (a) The diagram how the connection formulation. (b) The vertical incident wave field along the connecting boundary ABB. (c) Comparion of the horizontal component of the wave field along CDD calculated directly from the ource (thick line) and thoe calculated by uing the connection technique (thin line). (d) Comparion of the vertical component. In equation (25), both the initial diplacement and traction on the connection boundary are aumed to be known. u 0 (r) can be calculated from the previou ection, and incident traction t 0 (r) can alo be calculated from the previou ection. Alternatively, we can ue the elatic Rayleigh integral (Wu, 1989), which contain only either the diplacement or the traction field. We can reduce the urface integral in the right ide of equation (25) to the elatic wave Rayleigh integral that only contain the initial diplacement u 0 (r): C(r)u (r) C1 R(r,r)u (r)dr C2 [(r,r)u (r) G(r,r)t (r)]dr (26) 2 CAB [(r,r)u (r)]dr. 0 We ee that the unknown in equation (26) are u (r) onc 1 and u (r) and t (r) onc 2. To olve for u (r) and t (r), we mut build the correponding boundary integral equation in ubdomain X 2 : C(r)u (r) R(r,r)u (r)dr (r,r)[(u (r) 0 C2 CBB u (r)]dr (r,r)u 0(r)dr CDD (27) G(r,r)t (r)dr G(r,r)[(t (r) 0 C2 CBB G(r,r)t 0(r)dr, CDD t (r)]dr where a ufficiently long boundary C BB i ued with it end et a an infinite element. Similarly, becaue there i no di-
8 1382 Z. Ge, L.-Y. Fu, and R.-S. Wu continuity acro C BB and C DD, we can aume that C BB and C DD are tranparent and reduce equation (27) to C(r)u (r) [R(r,r)u (r) G(r,r)t (r)]dr C2 (28) [G(r,r)t (r) (r,r)u (r)]dr. 0 0 CBB The integration term containing the initial traction t 0 (r) in the right ide of equation (28) can be removed by uing the elatic wave Rayleigh integral repreentation, yielding C(r)u (r) [R(r,r)u (r) G(r,r)t (r)]dr C2 (29) 2 CBB (r,r)u (r)]dr. 0 The continuity of the diplacement and traction acro interface C BD i ued when equation (26) and (29) are combined to olve the problem. By olving the joint boundary integral equation of X 2 and X 2, we can obtain the wave field u (r) onc 1, u (r) and t (r) onc 2. The oberved field along C CD i calculated explicitly from the field on the boundarie. Similar to equation (18) and (19), equation (26) and (29) can be expreed in the matrix form and (1) (1) (2) (2) (2) (2) H u H u G t F, (30) (2) (2) (2) (2) H ũ G t F, (31) with the continuou condition acro C 2 u (2) (2) ũ (2) (2) t t (32) By uing thi technique, wave propagation in a long regional wave guide can be partitioned into everal ection contribution for great aving in both CPU time and memory. For example, if we divide a large model into N ection, we can reduce the ize of matrix by N time. Thu the total compute time will be reduced by N * (1/N) 3 1/N 2 time. Validation for P-SV Connection Technique To validate the connection technique, we preent a comparion between the wave field obtained by uing the BE method to directly calculate wave propagation from the ource to the obervation urface C CD and the wave field calculated by the connection cheme (Fig. 6). In both cae the ource time function are the ame, with a dominant frequency of 1 Hz. Firt the intermediate wave field u 0 (r) on C AB (hown in Fig. 6b) calculated from the ource i ued a the incident field for wave propagation in X 2 and X 2. The dominant arrival for the incident field at C AB conit of direct P wave, ps, pp, and multireflection between the two layer. Then the BE method i ued to calculate wave propagation from C AB to C CD. More multiply reflected wave between the free urface and the interface can be clearly een in the eimogram at C CD. The excellent agreement between the wave field calculated by the two method hown in Figure 6c and 6d confirm the validity of the connection technique. Although thi model i only deigned to tet the validity of the connection technique, we can ee computation aving. For a ingle-proceor Pentium IV 2.0 Hz computer, it take about 2 hr to calculate the wave field along C CD directly from the ource, wherea it take only 45 min to compute the wave field along C CD with the connection technique. The memory requirement i alo reduced by four time. Numerical Example and Application In thi ection, we ue the BE connection technique given previouly to imulate the wave field for everal model with rough topographie. The firt model hown in the top panel of Figure 7 i a two-layer crutal model with an irregular topography. There are two dicontinuou interface located at the depth of 15 km and 37 km, repectively, with the propertie lited in Table 1. The correlation length of the topographic fluctuation i 5 km and the root-mean-quare amplitude i 0.5 km. The receiver are along the urface. An exploive ource i located at the depth of 2 km, and the ource time function i a Ricker wavelet with a dominant frequency of 1 Hz. Figure 7 how the ynthetic eimogram received along the rough urface. Comparing with the wave field (Fig. 8) for the model with a flat urface, we can ee that the amplitude of the direct P and Rayleigh wave are diminihed becaue of the cattering effect of the rough topography. Both the forward-cattering and backwardcattering wave can be een in Figure 7. The top panel of Figure 9 how four different topographic curve ued with the crutal model of Figure 7. The ynthetic eimogram are calculated for the exponential and Gauian random topographie, repectively, with each including two rm amplitude of 0.5 km and 0.25 km. The energy attenuation curve correponding to thee topographie are calculated with increaing propagation ditance and are hown in the lower panel of Figure 9. Thee energy Table 1 Propertie of Crut Model Depth (km) V p (km/ec) V (km/ec) P (g/cm 3 ) Infinity
9 P-SV Wave-Field Connection Technique for Regional Wave Propagation Simulation 1383 Figure 7. Synthetic eimogram for a two-layered model with random topography. The ource i located at the depth of 2 km. The topographic fluctuation ha an exponential correlation function with a rm of 0.5 km and a correlation length of 5 km. We can ee both forward-cattering and backcattering in the eimogram. Figure 8. Synthetic eimogram for the model hown in Figure 7 but with flat urface. Multiply-reflected wave and Rayleigh wave can be clearly een in the figure.
10 1384 Z. Ge, L.-Y. Fu, and R.-S. Wu Figure 9. Energy attenuation for earth model with different topographie. On the top are the topographie ued in our calculation, the thick olid line i a Gauian random topography with a rm of 0.5 km, the thin olid line i a Gauian random topography with a rm of 0.25 km, the thick dahed line i exponential random topography with rm 0.5 km, and the thin dahed line i exponential random topography with rm 0.25 km, while the reference model (flat urface) i hown by the dotted line. On the bottom are the energy attenuation curve for the whole time window, the line have the ame meaning a thoe in the top figure. attenuation curve manifet different degree of topographycattering attenuation. We can ee that the exponential topographie lead to more effective energy attenuation than Gauian topographie, and alo the larger the topographic fluctuation i, the tronger i the topographic cattering. Next we apply the connection technique to regional wave propagation imulation in another ynthetic model. A hown in Figure 10a, the model i a laterally varying crutal wave guide, 600 km long and 32 km thick, overlaying a mantle half-pace. The model i divided into three 200- km-long egment. The firt egment contain a Gauian hill which i given by h(x) h 2 (x x 0) 0 exp 2 where x 0 62 km, h 0 4 km, r km. The econd egment preent a thinned wave guide with the 7-km-high tep. The third egment contain an exponential random topography with correlation length of 5 km and rm height of 0.6 km. The point ource i located at 8 km depth. The receiver are along both the free urface and vertical profile that lice thee wave guide. The ynthetic eimogram are calculated for the frequency range of Hz. We firt compute wave propagation from the ource to produce the 2r
11 P-SV Wave-Field Connection Technique for Regional Wave Propagation Simulation 1385 Figure 10. Application to a large model. (a) The model. (b) Horizontal wave field on C 1. (c) Horizontal wave field on C 2. (d) Horizontal wave field on C 3. (e) Horizontal wave field along the urface. (f) Energy ditribution along the urface. incident wave field on the firt connection boundary C 1 hown in Figure 10b, where the multiply reflected wave, converted wave, head wave, and Rayleigh wave can all be clearly een. Subequently the BE method i ued to calculate wave propagation in the econd ection and obtain the wave field on the econd connection boundary C 2 hown in Figure 10c. Finally, we ue the wave field on C 2 a the incident field and calculate the wave field in the third ection. Figure 10d how the horizontal wave field on C 3, where the cattering effect of the random topography i obviou in that the Rayleigh wave are motly cattered and ome coda wave appear. The wave field received, repectively, along the free urface of the three ection are put together and hown in Figure 10e. Figure 10f i the energy ditribution curve along the entire free urface, where a relatively low energy appear around the hill and the attenuation curve varie dramatically along the random topography. The topography ha a great effect on the energy attenuation in the wave guide.
12 1386 Z. Ge, L.-Y. Fu, and R.-S. Wu Concluion and Dicuion We preent the boundary element method for the 2D P- SV elatic wave problem. To model the effect of irregular topography on regional wave propagation, we further developed a connection technique to imulate long-range wave propagation ection by ection. Numerical comparion with independent method indicated that the preented method including the connection technique i accurate for regional wave imulation. Numerical reult how that irregular topography can attenuate the total wave energy propagating in the crutal waveguide. The boundary element method can be ued in computing the ite effect on ite uch a canyon, mountain, and valley. The connection technique expand thi method to deal with large earth model with irregular topography. To extend the boundary element method to 3D cae, everal change have to be made. Firt, the Green function of diplacement and traction of the 3D cae hould be ued. Second, becaue the ingularity of the 3D Green function i 1/r and 1/r 2 while in 2D cae thoe are 1/r 1/2 and 1/r, the integral hould be treated more carefully. Third, the ize of matrice to olve the problem will be quared to the 2D cae, becaue thi method include the matrix inverion, the requirement of memory and computation time will be enlarged. Extenion to 3D cae will be conducted in the future. Acknowledgment Thi work wa done while we were working in the Intitute of Geophyic and Planetary Phyic, Univerity of California, Santa Cruz. The reearch i upported by the Air Force Reearch Laboratory under Grant no. DTRA01-01-C The facility upport from the W. M. Keck Foundation i alo acknowledged. Reference Alvarez-Rubioa, S., F. J. Sanchez-Sema, J. J. Benitoc, and E. Alarcond (2004). The direct boundary element method: 2D ite effect aement on laterally varying layered media (methodology), Soil Dyn. Earthquake Eng. 24, Aki, K., and P. G. Richard (1980). Quantitative Seimology. Theory and Method. W. H. Freeman, San Francico. Botock, M. G., and B. L. N. Kennett (1990). The effect of three-dimenional tructure on Lg propagation pattern, Geophy. J. Int. 101, Bouchon, M., C. A. Schultz, and M. N. Toköz (1995). A fat implementation of boundary integral equation method to calculate the propagation of eimic wave in laterally varying layered media, Bull. Seim. Soc. Am. 85, Fu, L. Y. (1996). 3-D boundary element eimic modeling in complex geology, in 66th Ann. Internat. Mtg., Soc. Expl. Geophy., Expanded Abtract, Fu, L. Y., and Y. G. Mu (1994). Boundary element method for elatic wave forward modeling, Acta Geophy. Sinica, 37, Fu, L. Y., and R. S. Wu (2000). Infinite boundary element aborbing boundary for wave propagation imulation, Geophyic 65, Fu, L. Y., and R. S. Wu (2001). A hybrid BE-GS method for modeling regional wave propagation, Pure Appl. Geophy. 158, Gibon, R. L., Jr., and M. Campillo (1994). Numerical imulation of highand low-frequency Lg-wave propagation, Geophy. J. Int. 118, Jih, R. S. (1996). Waveguide effect of large-cale tructural variation, anelatic attenuation, and random heterogeneity on SV Lg propagation: a finite-difference modeling tudy, in Proc. 18th Annual Seimic Reearch Smpoium on Monitoring a Comprehenive Tet Ban Treaty, Keer, H., G. Nolet, and F. A. Dahlen (1996). Ray theoretical analyi of Lg, Bull. Seim. Soc. Am. 86, More, P. M., and H. Fehbach (1953). Method of Theoretical Phyic, McGraw-Hill, New York. Sánchez-Sema, F. J., M. A. Bravo, and I. Herrera (1985). Surface motion of topographical irregularitie for incident P, SV, and Rayleigh wave, Bull. Seim. Soc. Am. 75, Sánchez-Sema, F. J., and M. Campillo (1991). Diffraction of P, SV and Rayleigh wave by topographic feature: a boundary integral formulation, Bull. Seim. Soc. Am. 81, Wu, R. S. (1989). Repreentation integral for elatic wave propagation containing either the diplacement term or the tre term alone, Phy. Rev. Lett. 62, Wu, R. S., S. Jin, and X. B. Xie (1996). Synthetic eimogram in heterogeneou crutal waveguide uing creen propagator, in Proc. 18th Annual Seimic Reearch Sympoium on Monitoring a Comprehenive Tet Ban Treaty, Wu, R. S., S. Jin, and X. B. Xie (2000). Seimic wave propagation and cattering in heterogeneou crutal waveguide uing creen propagator. I: SH wave, Bull. Seim. Soc. Am. 90, Wu, X. Y., and R. S. Wu (2001). Lg wave imulation in heterogeneou crut with irregular topography uing halfpace creen propagator, Geophy. J. Int. 146, Intitute of Geophyic and Planetary Phyic Univerity of California 1156 High Street Santa Cruz, California (Z.G., L.-Y.F., R.-S.W.) Laboratory of Computational Geodynamic Department of Geophyic, School of Earth and Space Science Peking Univerity Beijing , People Republic of China zxgai@pku.edu.cn (Z.G.) Manucript received 20 Augut 2004.
SEG/San Antonio 2007 Annual Meeting
Yaofeng He* and Ru-Shan Wu Modeling and Imaging Laboratory, Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, CA, 95064, USA Summary A one-way and one-return boundary
More informationUniversität Augsburg. Institut für Informatik. Approximating Optimal Visual Sensor Placement. E. Hörster, R. Lienhart.
Univerität Augburg à ÊÇÅÍÆ ËÀǼ Approximating Optimal Viual Senor Placement E. Hörter, R. Lienhart Report 2006-01 Januar 2006 Intitut für Informatik D-86135 Augburg Copyright c E. Hörter, R. Lienhart Intitut
More information( ) subject to m. e (2) L are 2L+1. = s SEG SEG Las Vegas 2012 Annual Meeting Page 1
A new imultaneou ource eparation algorithm uing frequency-divere filtering Ying Ji*, Ed Kragh, and Phil Chritie, Schlumberger Cambridge Reearch Summary We decribe a new imultaneou ource eparation algorithm
More informationMirror shape recovery from image curves and intrinsic parameters: Rotationally symmetric and conic mirrors. Abstract. 2. Mirror shape recovery
Mirror hape recovery from image curve and intrinic parameter: Rotationally ymmetric and conic mirror Nuno Gonçalve and Helder Araújo Λ Intitute of Sytem and Robotic Univerity of Coimbra Pinhal de Marroco
More informationDownloaded 11/08/14 to Redistribution subject to SEG license or copyright; see Terms of Use at
Downloaded 11/08/14 to 50.131.5.33. Reditribution ubject to SEG licene or copyright; ee Term of Ue at http://library.eg.org/ Gauian beam baed enitivity kernel calculation and it application in turning
More informationPrestack Migration/Imaging Using Synthetic Beam-Sources and Plane-Sources
Pretack Migration/Imaging Uing Synthetic Beam-Source and Plane-Source Ru-Shan Wu, Ling Chen and Yongzhong Wang Modeling and Imaging Laboratory, IGPP, Univerity of California, Santa Cruz, CA 95064, USA
More information3D SMAP Algorithm. April 11, 2012
3D SMAP Algorithm April 11, 2012 Baed on the original SMAP paper [1]. Thi report extend the tructure of MSRF into 3D. The prior ditribution i modified to atify the MRF property. In addition, an iterative
More informationPerformance of a Robust Filter-based Approach for Contour Detection in Wireless Sensor Networks
Performance of a Robut Filter-baed Approach for Contour Detection in Wirele Senor Network Hadi Alati, William A. Armtrong, Jr., and Ai Naipuri Department of Electrical and Computer Engineering The Univerity
More informationPlanning of scooping position and approach path for loading operation by wheel loader
22 nd International Sympoium on Automation and Robotic in Contruction ISARC 25 - September 11-14, 25, Ferrara (Italy) 1 Planning of cooping poition and approach path for loading operation by wheel loader
More informationNumerical Modeling of Material Discontinuity Using Mixed MLPG Collocation Method
umerical odeling of aterial Dicontinuit Uing ied LPG Collocation ethod B. alušić 1,. Sorić 1 and T. arak 1 Abtract A mied LPG collocation method i applied for the modeling of material dicontinuit in heterogeneou
More informationKINEMATIC BENDING OF FIXED-HEAD PILES IN NON- HOMOGENEOUS SOIL
KINEMATIC BENDING OF FIXED-HEAD PILES IN NON- HOMOGENEOUS SOIL Raffaele Di Laora Univerity of Napoli Parthenope raffaele.dilaora@uniparthenope.it Emmanouil Rovithi Intitute of Engineering Seimology and
More informationES205 Analysis and Design of Engineering Systems: Lab 1: An Introductory Tutorial: Getting Started with SIMULINK
ES05 Analyi and Deign of Engineering Sytem: Lab : An Introductory Tutorial: Getting Started with SIMULINK What i SIMULINK? SIMULINK i a oftware package for modeling, imulating, and analyzing dynamic ytem.
More information[N309] Feedforward Active Noise Control Systems with Online Secondary Path Modeling. Muhammad Tahir Akhtar, Masahide Abe, and Masayuki Kawamata
he 32nd International Congre and Expoition on Noie Control Engineering Jeju International Convention Center, Seogwipo, Korea, Augut 25-28, 2003 [N309] Feedforward Active Noie Control Sytem with Online
More informationRepresentations and Transformations. Objectives
Repreentation and Tranformation Objective Derive homogeneou coordinate tranformation matrice Introduce tandard tranformation - Rotation - Tranlation - Scaling - Shear Scalar, Point, Vector Three baic element
More informationMotion Control (wheeled robots)
3 Motion Control (wheeled robot) Requirement for Motion Control Kinematic / dynamic model of the robot Model of the interaction between the wheel and the ground Definition of required motion -> peed control,
More informationTexture-Constrained Active Shape Models
107 Texture-Contrained Active Shape Model Shuicheng Yan, Ce Liu Stan Z. Li Hongjiang Zhang Heung-Yeung Shum Qianheng Cheng Microoft Reearch Aia, Beijing Sigma Center, Beijing 100080, China Dept. of Info.
More informationCompressed Sensing Image Processing Based on Stagewise Orthogonal Matching Pursuit
Senor & randucer, Vol. 8, Iue 0, October 204, pp. 34-40 Senor & randucer 204 by IFSA Publihing, S. L. http://www.enorportal.com Compreed Sening Image Proceing Baed on Stagewie Orthogonal Matching Puruit
More informationAnalysis of Surface Wave Propagation Based on the Thin Layered Element Method
ABSTACT : Analyi of Surface Wave Propagation Baed on the Thin ayered Element Method H. AKAGAWA 1 and S. AKAI 1 Engineer, Technical Centre, Oyo Corporation, Ibaraki, Japan Profeor, Graduate School of Engineering,
More informationxy-monotone path existence queries in a rectilinear environment
CCCG 2012, Charlottetown, P.E.I., Augut 8 10, 2012 xy-monotone path exitence querie in a rectilinear environment Gregory Bint Anil Mahehwari Michiel Smid Abtract Given a planar environment coniting of
More informationANALYSIS OF PILE DRIVING IN VERTICAL AND HORIZONTAL DIRECTIONS USING A HYBRID MODEL
ANALYSIS OF PILE DRIVING IN VERTICAL AND HORIZONTAL DIRECTIONS USING A HYBRID MODEL TATSUNORI MATSUMOTO i) PASTSAKORN KITIYODOM ii) EIJI KOJIMA iii) HIROMICHI KUMAGAI iv) SATOSHI NISHIMOTO v) and KOUICHI
More informationA METHOD OF REAL-TIME NURBS INTERPOLATION WITH CONFINED CHORD ERROR FOR CNC SYSTEMS
Vietnam Journal of Science and Technology 55 (5) (017) 650-657 DOI: 10.1565/55-518/55/5/906 A METHOD OF REAL-TIME NURBS INTERPOLATION WITH CONFINED CHORD ERROR FOR CNC SYSTEMS Nguyen Huu Quang *, Banh
More informationDistributed Packet Processing Architecture with Reconfigurable Hardware Accelerators for 100Gbps Forwarding Performance on Virtualized Edge Router
Ditributed Packet Proceing Architecture with Reconfigurable Hardware Accelerator for 100Gbp Forwarding Performance on Virtualized Edge Router Satohi Nihiyama, Hitohi Kaneko, and Ichiro Kudo Abtract To
More informationA note on degenerate and spectrally degenerate graphs
A note on degenerate and pectrally degenerate graph Noga Alon Abtract A graph G i called pectrally d-degenerate if the larget eigenvalue of each ubgraph of it with maximum degree D i at mot dd. We prove
More informationFinite Elements Method in Split Hopkinson Pressure Bar developing process
6th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING, Venice, Italy, November 1-3, 007 63 Finite Element Method in Split Hopkinon Preure Bar developing proce ROTARIU ADRIAN
More informationAnisotropic filtering on normal field and curvature tensor field using optimal estimation theory
Aniotropic filtering on normal field and curvature tenor field uing optimal etimation theory Min Liu Yuhen Liu and Karthik Ramani Purdue Univerity, Wet Lafayette, Indiana, USA Email: {liu66 liu28 ramani}@purdue.edu
More informationNew DSP to measure acoustic efficiency of road barriers. Part 2: Sound Insulation Index
New DSP to meaure acoutic efficiency of road barrier. Part 2: Sound Inulation Index LAMBERTO TRONCHIN 1, KRISTIAN FABBRI 1, JELENA VASILJEVIC 2 1 DIENCA CIARM, Univerity of Bologna, Italy 2 Univerity of
More informationDAROS: Distributed User-Server Assignment And Replication For Online Social Networking Applications
DAROS: Ditributed Uer-Server Aignment And Replication For Online Social Networking Application Thuan Duong-Ba School of EECS Oregon State Univerity Corvalli, OR 97330, USA Email: duongba@eec.oregontate.edu
More informationDark-field Tomography: Modeling and Reconstruction
Dark-field Tomography: Modeling and Recontruction W. Cong, F. Pfeiffer, M. Bech, O. Bunk, C. David, and G. Wang Abtract: Dark-field image are formed by mall-angle cattering of x-ray photon. The mall-angle
More informationMAT 155: Describing, Exploring, and Comparing Data Page 1 of NotesCh2-3.doc
MAT 155: Decribing, Exploring, and Comparing Data Page 1 of 8 001-oteCh-3.doc ote for Chapter Summarizing and Graphing Data Chapter 3 Decribing, Exploring, and Comparing Data Frequency Ditribution, Graphic
More informationLaboratory Exercise 6
Laboratory Exercie 6 Adder, Subtractor, and Multiplier The purpoe of thi exercie i to examine arithmetic circuit that add, ubtract, and multiply number. Each type of circuit will be implemented in two
More informationRouting Definition 4.1
4 Routing So far, we have only looked at network without dealing with the iue of how to end information in them from one node to another The problem of ending information in a network i known a routing
More informationKS3 Maths Assessment Objectives
KS3 Math Aement Objective Tranition Stage 9 Ratio & Proportion Probabilit y & Statitic Appreciate the infinite nature of the et of integer, real and rational number Can interpret fraction and percentage
More informationThe Association of System Performance Professionals
The Aociation of Sytem Performance Profeional The Computer Meaurement Group, commonly called CMG, i a not for profit, worldwide organization of data proceing profeional committed to the meaurement and
More informationDevelopment of an atmospheric climate model with self-adapting grid and physics
Intitute of Phyic Publihing Journal of Phyic: Conference Serie 16 (2005) 353 357 doi:10.1088/1742-6596/16/1/049 SciDAC 2005 Development of an atmopheric climate model with elf-adapting grid and phyic Joyce
More informationCENTER-POINT MODEL OF DEFORMABLE SURFACE
CENTER-POINT MODEL OF DEFORMABLE SURFACE Piotr M. Szczypinki Iintitute of Electronic, Technical Univerity of Lodz, Poland Abtract: Key word: Center-point model of deformable urface for egmentation of 3D
More informationA System Dynamics Model for Transient Availability Modeling of Repairable Redundant Systems
International Journal of Performability Engineering Vol., No. 3, May 05, pp. 03-. RAMS Conultant Printed in India A Sytem Dynamic Model for Tranient Availability Modeling of Repairable Redundant Sytem
More informationOn successive packing approach to multidimensional (M-D) interleaving
On ucceive packing approach to multidimenional (M-D) interleaving Xi Min Zhang Yun Q. hi ankar Bau Abtract We propoe an interleaving cheme for multidimenional (M-D) interleaving. To achieved by uing a
More informationOperational Semantics Class notes for a lecture given by Mooly Sagiv Tel Aviv University 24/5/2007 By Roy Ganor and Uri Juhasz
Operational emantic Page Operational emantic Cla note for a lecture given by Mooly agiv Tel Aviv Univerity 4/5/7 By Roy Ganor and Uri Juhaz Reference emantic with Application, H. Nielon and F. Nielon,
More informationIMPLEMENTATION OF AREA, VOLUME AND LINE SOURCES
December 01 ADMS 5 P503I1 IMPEMENTATION OF AREA, VOUME AND INE SOURCES The Met. Office (D J Thomon) and CERC 1. INTRODUCTION ADMS model line ource, and area and volume ource with conve polgon bae area.
More informationAdvanced Encryption Standard and Modes of Operation
Advanced Encryption Standard and Mode of Operation G. Bertoni L. Breveglieri Foundation of Cryptography - AES pp. 1 / 50 AES Advanced Encryption Standard (AES) i a ymmetric cryptographic algorithm AES
More informationLaboratory Exercise 6
Laboratory Exercie 6 Adder, Subtractor, and Multiplier The purpoe of thi exercie i to examine arithmetic circuit that add, ubtract, and multiply number. Each circuit will be decribed in Verilog and implemented
More informationResearch note: Calculating spectral irradiance indoors
Lighting Re. Technol. 217; Vol. 49: 122 127 Reearch note: Calculating pectral irradiance indoor S Bará PhD a and J Ecofet PhD b a Área de Óptica, Facultade de Óptica e Optometría, Univeridade de Santiago
More informationDomain-Specific Modeling for Rapid System-Wide Energy Estimation of Reconfigurable Architectures
Domain-Specific Modeling for Rapid Sytem-Wide Energy Etimation of Reconfigurable Architecture Seonil Choi 1,Ju-wookJang 2, Sumit Mohanty 1, Viktor K. Praanna 1 1 Dept. of Electrical Engg. 2 Dept. of Electronic
More informationMotivation: Level Sets. Input Data Noisy. Easy Case Use Marching Cubes. Intensity Varies. Non-uniform Exposure. Roger Crawfis
Level Set Motivation: Roger Crawfi Slide collected from: Fan Ding, Charle Dyer, Donald Tanguay and Roger Crawfi 4/24/2003 R. Crawfi, Ohio State Univ. 109 Eay Cae Ue Marching Cube Input Data Noiy 4/24/2003
More informationMinimum congestion spanning trees in bipartite and random graphs
Minimum congetion panning tree in bipartite and random graph M.I. Otrovkii Department of Mathematic and Computer Science St. John Univerity 8000 Utopia Parkway Queen, NY 11439, USA e-mail: otrovm@tjohn.edu
More informationQuadrilaterals. Learning Objectives. Pre-Activity
Section 3.4 Pre-Activity Preparation Quadrilateral Intereting geometric hape and pattern are all around u when we tart looking for them. Examine a row of fencing or the tiling deign at the wimming pool.
More informationNUMERICAL MODELING ON THE DAMPING CONTROL OF TLD STRUCTURE
4th International Conference on Earthquake Engineering Taipei, Taiwan October 12-13, 2006 Paper No. 183 NUMERICAL MODELING ON THE DAMPING CONTROL OF TLD STRUCTURE Han jun 1, Li Yingmin 2, Liu Liping 3,
More informationLecture 14: Minimum Spanning Tree I
COMPSCI 0: Deign and Analyi of Algorithm October 4, 07 Lecture 4: Minimum Spanning Tree I Lecturer: Rong Ge Scribe: Fred Zhang Overview Thi lecture we finih our dicuion of the hortet path problem and introduce
More informationFocused Video Estimation from Defocused Video Sequences
Focued Video Etimation from Defocued Video Sequence Junlan Yang a, Dan Schonfeld a and Magdi Mohamed b a Multimedia Communication Lab, ECE Dept., Univerity of Illinoi, Chicago, IL b Phyical Realization
More informationIntroduction to PET Image Reconstruction. Tomographic Imaging. Projection Imaging. PET Image Reconstruction 11/6/07
Introduction to PET Image Recontruction Adam Aleio Nuclear Medicine Lecture Imaging Reearch Laboratory Diviion of Nuclear Medicine Univerity of Wahington Fall 2007 http://dept.wahington.edu/nucmed/irl/education.html
More informationIncreasing Throughput and Reducing Delay in Wireless Sensor Networks Using Interference Alignment
Int. J. Communication, Network and Sytem Science, 0, 5, 90-97 http://dx.doi.org/0.436/ijcn.0.50 Publihed Online February 0 (http://www.scirp.org/journal/ijcn) Increaing Throughput and Reducing Delay in
More informationDrawing Lines in 2 Dimensions
Drawing Line in 2 Dimenion Drawing a traight line (or an arc) between two end point when one i limited to dicrete pixel require a bit of thought. Conider the following line uperimpoed on a 2 dimenional
More informationBuilding a Compact On-line MRF Recognizer for Large Character Set using Structured Dictionary Representation and Vector Quantization Technique
202 International Conference on Frontier in Handwriting Recognition Building a Compact On-line MRF Recognizer for Large Character Set uing Structured Dictionary Repreentation and Vector Quantization Technique
More informationEdits in Xylia Validity Preserving Editing of XML Documents
dit in Xylia Validity Preerving diting of XML Document Pouria Shaker, Theodore S. Norvell, and Denni K. Peter Faculty of ngineering and Applied Science, Memorial Univerity of Newfoundland, St. John, NFLD,
More informationModeling of underwater vehicle s dynamics
Proceeding of the 11th WEA International Conference on YTEM, Agio Nikolao, Crete Iland, Greece, July 23-25, 2007 44 Modeling of underwater vehicle dynamic ANDRZEJ ZAK Department of Radiolocation and Hydrolocation
More informationKinematics Programming for Cooperating Robotic Systems
Kinematic Programming for Cooperating Robotic Sytem Critiane P. Tonetto, Carlo R. Rocha, Henrique Sima, Altamir Dia Federal Univerity of Santa Catarina, Mechanical Engineering Department, P.O. Box 476,
More informationA NEW APPROACH IN MEASURING OF THE ROUGHNESS FOR SURFACE CONSTITUTED WITH MACHINING PROCESS BY MATERIAL REMOVAL
International Journal of Mechanical and Production Engineering Reearch and Development (IJMPERD) ISSN 49-689 Vol. 3, Iue, Mar 3, 4-5 TJPRC Pvt. Ltd. A NEW APPROACH IN MEASURING OF THE ROUGHNESS FOR SURFACE
More informationx y z Design variable positions A
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2001 00:1{7 [Verion: 2000/03/22 v1.0] A tabilied peudo-hell approach for urface parametriation in CFD deign problem O. Soto,R.Lohner
More information3D High-Resolution Imaging Using Separated Wavefields
3D High-Reolution Imaging Uing Separated Wavefield S. Lu (Petroleum Geo-Service), N.D. Whitmore (Petroleum Geo- Service), H. LeGleut (Petroleum Geo-Service) & A. Long (Petroleum Geo- Service) SUMMARY Conventional
More informationANALYSIS OF THE FIRST LAYER IN WEIGHTLESS NEURAL NETWORKS FOR 3_DIMENSIONAL PATTERN RECOGNITION
ANALYSIS OF THE FIRST LAYER IN WEIGHTLESS NEURAL NETWORKS FOR 3_DIMENSIONAL PATTERN RECOGNITION A. Váque-Nava * Ecuela de Ingeniería. CENTRO UNIVERSITARIO MEXICO. DIVISION DE ESTUDIOS SUPERIORES J. Figueroa
More informationUC Berkeley International Conference on GIScience Short Paper Proceedings
UC Berkeley International Conference on GIScience Short Paper Proceeding Title A novel method for probabilitic coverage etimation of enor network baed on 3D vector repreentation in complex urban environment
More informationImplementation of a momentum-based distance metric for motion graphs. Student: Alessandro Di Domenico (st.no ), Supervisor: Nicolas Pronost
Implementation of a momentum-baed ditance metric for motion graph Student: Aleandro Di Domenico (t.no 3775682), Supervior: Nicola Pronot April 3, 2014 Abtract Thi report preent the procedure and reult
More informationGeneration of nearly nondiffracting Bessel beams with a Fabry Perot interferometer
Horváth et al. Vol. 14, No. 11/November 1997/J. Opt. Soc. Am. A 3009 Generation of nearly nondiffracting Beel beam with a Fabry Perot interferometer Z. L. Horváth, M. Erdélyi, G. Szabó, and Z. Bor Department
More informationAdaptive Beam Director for a Tiled Fiber Array
Adaptive Beam Director for a Tiled Fiber Array Mikhail A. Vorontov Intelligent Optic Laboratory, Computational and Information Science Directorate, U.S. Army Reearch Laboratory, 2800 Powder Mill Road,
More informationGray-level histogram. Intensity (grey-level) transformation, or mapping. Use of intensity transformations:
Faculty of Informatic Eötvö Loránd Univerity Budapet, Hungary Lecture : Intenity Tranformation Image enhancement by point proceing Spatial domain and frequency domain method Baic Algorithm for Digital
More informationA Practical Model for Minimizing Waiting Time in a Transit Network
A Practical Model for Minimizing Waiting Time in a Tranit Network Leila Dianat, MASc, Department of Civil Engineering, Sharif Univerity of Technology, Tehran, Iran Youef Shafahi, Ph.D. Aociate Profeor,
More informationCutting Stock by Iterated Matching. Andreas Fritsch, Oliver Vornberger. University of Osnabruck. D Osnabruck.
Cutting Stock by Iterated Matching Andrea Fritch, Oliver Vornberger Univerity of Onabruck Dept of Math/Computer Science D-4909 Onabruck andy@informatikuni-onabrueckde Abtract The combinatorial optimization
More informationLinkGuide: Towards a Better Collection of Hyperlinks in a Website Homepage
Proceeding of the World Congre on Engineering 2007 Vol I LinkGuide: Toward a Better Collection of Hyperlink in a Webite Homepage A. Ammari and V. Zharkova chool of Informatic, Univerity of Bradford anammari@bradford.ac.uk,
More informationAalborg Universitet. Published in: Proceedings of the Working Conference on Advanced Visual Interfaces
Aalborg Univeritet Software-Baed Adjutment of Mobile Autotereocopic Graphic Uing Static Parallax Barrier Paprocki, Martin Marko; Krog, Kim Srirat; Kritofferen, Morten Bak; Krau, Martin Publihed in: Proceeding
More information/06/$ IEEE 364
006 IEEE International ympoium on ignal Proceing and Information Technology oie Variance Etimation In ignal Proceing David Makovoz IPAC, California Intitute of Technology, MC-0, Paadena, CA, 95 davidm@ipac.caltech.edu;
More informationAspects of Formal and Graphical Design of a Bus System
Apect of Formal and Graphical Deign of a Bu Sytem Tiberiu Seceleanu Univerity of Turku, Dpt. of Information Technology Turku, Finland tiberiu.eceleanu@utu.fi Tomi Weterlund Turku Centre for Computer Science
More informationUNIAXIAL COMPRESSION TEST AND STRESS WAVE PROPAGATION MODELLING USING SPH
Fifth International Conference on CFD in the Proce Indutrie CSIRO, Melbourne, Autralia - December UNIAXIAL COMPRESSION TEST AND STRESS WAVE PROPAGATION MODELLING USING SPH Rajarhi DAS and Paul W. CLEARY
More informationA CLUSTERING-BASED HYBRID REPLICA CONTROL PROTOCOL FOR HIGH AVAILABILITY IN GRID ENVIRONMENT
Journal of Computer Science 10 (12): 2442-2449, 2014 ISSN: 1549-3636 2014 R. Latip et al., Thi open acce article i ditributed under a Creative Common Attribution (CC-BY) 3.0 licene doi:10.3844/jcp.2014.2442.2449
More informationComputer Aided Drafting, Design and Manufacturing Volume 25, Number 3, September 2015, Page 10
Computer Aided Drafting, Deign and Manufacturing Volume 5, umber 3, September 015, Page 10 CADDM Reearch of atural Geture Recognition and Interactive Technology Compatible with YCbCr and SV Color Space
More informationA Load Balancing Model based on Load-aware for Distributed Controllers. Fengjun Shang, Wenjuan Gong
4th International Conference on Machinery, Material and Computing Technology (ICMMCT 2016) A Load Balancing Model baed on Load-aware for Ditributed Controller Fengjun Shang, Wenjuan Gong College of Compute
More informationHow to. write a paper. The basics writing a solid paper Different communities/different standards Common errors
How to write a paper The baic writing a olid paper Different communitie/different tandard Common error Reource Raibert eay My grammar point Article on a v. the Bug in writing Clarity Goal Conciene Calling
More informationA Study of a Variable Compression Ratio and Displacement Mechanism Using Design of Experiments Methodology
A Study of a Variable Compreion Ratio and Diplacement Mechanim Uing Deign of Experiment Methodology Shugang Jiang, Michael H. Smith, Maanobu Takekohi Abtract Due to the ever increaing requirement for engine
More informationLaboratory Exercise 6
Laboratory Exercie 6 Adder, Subtractor, and Multiplier The purpoe of thi exercie i to examine arithmetic circuit that add, ubtract, and multiply number. Each circuit will be decribed in VHL and implemented
More informationKeywords: Defect detection, linear phased array transducer, parameter optimization, phased array ultrasonic B-mode imaging testing.
Send Order for Reprint to reprint@benthamcience.ae 488 The Open Automation and Control Sytem Journal, 2014, 6, 488-492 Open Acce Parameter Optimization of Linear Phaed Array Tranducer for Defect Detection
More informationHow to Select Measurement Points in Access Point Localization
Proceeding of the International MultiConference of Engineer and Computer Scientit 205 Vol II, IMECS 205, March 8-20, 205, Hong Kong How to Select Meaurement Point in Acce Point Localization Xiaoling Yang,
More informationShortest Paths with Single-Point Visibility Constraint
Shortet Path with Single-Point Viibility Contraint Ramtin Khoravi Mohammad Ghodi Department of Computer Engineering Sharif Univerity of Technology Abtract Thi paper tudie the problem of finding a hortet
More informationIMPLEMENTATION OF CHORD LENGTH SAMPLING FOR TRANSPORT THROUGH A BINARY STOCHASTIC MIXTURE
Nuclear Mathematical and Computational Science: A Century in Review, A Century Anew Gatlinburg, Tenneee, April 6-, 003, on CD-ROM, American Nuclear Society, LaGrange Park, IL (003) IMPLEMENTATION OF CHORD
More informationTrainable Context Model for Multiscale Segmentation
Trainable Context Model for Multicale Segmentation Hui Cheng and Charle A. Bouman School of Electrical and Computer Engineering Purdue Univerity Wet Lafayette, IN 47907-1285 {hui, bouman}@ ecn.purdue.edu
More informationRefining SIRAP with a Dedicated Resource Ceiling for Self-Blocking
Refining SIRAP with a Dedicated Reource Ceiling for Self-Blocking Mori Behnam, Thoma Nolte Mälardalen Real-Time Reearch Centre P.O. Box 883, SE-721 23 Väterå, Sweden {mori.behnam,thoma.nolte}@mdh.e ABSTRACT
More informationSLA Adaptation for Service Overlay Networks
SLA Adaptation for Service Overlay Network Con Tran 1, Zbigniew Dziong 1, and Michal Pióro 2 1 Department of Electrical Engineering, École de Technologie Supérieure, Univerity of Quebec, Montréal, Canada
More informationA Novel Feature Line Segment Approach for Pattern Classification
12th International Conference on Information Fuion Seattle, WA, USA, July 6-9, 2009 A Novel Feature Line Segment Approach for Pattern Claification Yi Yang Intitute of Integrated Automation Xi an Jiaotong
More informationA Multi-objective Genetic Algorithm for Reliability Optimization Problem
International Journal of Performability Engineering, Vol. 5, No. 3, April 2009, pp. 227-234. RAMS Conultant Printed in India A Multi-objective Genetic Algorithm for Reliability Optimization Problem AMAR
More informationA Simplified Method for Calculating the Effective Solar Optical Properties of a Venetian Blind Layer for Building Energy Simulation
A Simplified Method for Calculating the Effective Solar Optical Propertie of a Venetian Blind Layer for Building Energy Simulation N. A. Kotey, Department of Mechanical Engineering, Univerity of Waterloo,
More informationLaboratory Exercise 6
Laboratory Exercie 6 Adder, Subtractor, and Multiplier a a The purpoe of thi exercie i to examine arithmetic circuit that add, ubtract, and multiply number. Each b c circuit will be decribed in Verilog
More informationDynamically Reconfigurable Neuron Architecture for the Implementation of Self- Organizing Learning Array
Dynamically Reconfigurable Neuron Architecture for the Implementation of Self- Organizing Learning Array Januz A. Starzyk,Yongtao Guo, and Zhineng Zhu School of Electrical Engineering & Computer Science
More informationLocalized Minimum Spanning Tree Based Multicast Routing with Energy-Efficient Guaranteed Delivery in Ad Hoc and Sensor Networks
Localized Minimum Spanning Tree Baed Multicat Routing with Energy-Efficient Guaranteed Delivery in Ad Hoc and Senor Network Hanne Frey Univerity of Paderborn D-3398 Paderborn hanne.frey@uni-paderborn.de
More information3-D Visualization of a Gene Regulatory Network: Stochastic Search for Layouts
3-D Viualization of a Gene Regulatory Network: Stochatic Search for Layout Naoki Hooyama Department of Electronic Engineering, Univerity of Tokyo, Japan hooyama@iba.k.u-tokyo.ac.jp Abtract- In recent year,
More informationEnhanced subsurface illumination from separated wavefield imaging
firt break volume 32, November 2014 pecial topic Marine Seimic Enhanced uburface illumination from eparated wavefield imaging Shaoping Lu 1*, Dan Whitmore 1, Alejandro Valenciano 1 and Nizar Chemingui
More informationTopics. Lecture 37: Global Optimization. Issues. A Simple Example: Copy Propagation X := 3 B > 0 Y := 0 X := 4 Y := Z + W A := 2 * 3X
Lecture 37: Global Optimization [Adapted from note by R. Bodik and G. Necula] Topic Global optimization refer to program optimization that encompa multiple baic block in a function. (I have ued the term
More informationUSING ARTIFICIAL NEURAL NETWORKS TO APPROXIMATE A DISCRETE EVENT STOCHASTIC SIMULATION MODEL
USING ARTIFICIAL NEURAL NETWORKS TO APPROXIMATE A DISCRETE EVENT STOCHASTIC SIMULATION MODEL Robert A. Kilmer Department of Sytem Engineering Unite State Military Acaemy Wet Point, NY 1996 Alice E. Smith
More informationCERIAS Tech Report EFFICIENT PARALLEL ALGORITHMS FOR PLANAR st-graphs. by Mikhail J. Atallah, Danny Z. Chen, and Ovidiu Daescu
CERIAS Tech Report 2003-15 EFFICIENT PARALLEL ALGORITHMS FOR PLANAR t-graphs by Mikhail J. Atallah, Danny Z. Chen, and Ovidiu Daecu Center for Education and Reearch in Information Aurance and Security,
More informationAn Intro to LP and the Simplex Algorithm. Primal Simplex
An Intro to LP and the Simplex Algorithm Primal Simplex Linear programming i contrained minimization of a linear objective over a olution pace defined by linear contraint: min cx Ax b l x u A i an m n
More information1-D N umerical S imulation for Comprehensive G ophysical e P rospecting of S eismic R fraction e and H h igdensity R sistivity e
International Conference on Mechatronic, Electronic, Indutrial and Control Engineering (MEIC 04) -D N umerical S imulation for Comprehenive G ophyical e P ropecting of S eimic R fraction e and H h igdenity
More informationAN ALGORITHM FOR RESTRICTED NORMAL FORM TO SOLVE DUAL TYPE NON-CANONICAL LINEAR FRACTIONAL PROGRAMMING PROBLEM
RAC Univerity Journal, Vol IV, No, 7, pp 87-9 AN ALGORITHM FOR RESTRICTED NORMAL FORM TO SOLVE DUAL TYPE NON-CANONICAL LINEAR FRACTIONAL PROGRAMMING PROLEM Mozzem Hoain Department of Mathematic Ghior Govt
More informationA PROBABILISTIC NOTION OF CAMERA GEOMETRY: CALIBRATED VS. UNCALIBRATED
A PROBABILISTIC NOTION OF CAMERA GEOMETRY: CALIBRATED VS. UNCALIBRATED Jutin Domke and Yianni Aloimono Computational Viion Laboratory, Center for Automation Reearch Univerity of Maryland College Park,
More information