NON-LINEAR INVERSION OF SEISMIC DATA FOR THE DETERMINATION OF REFLECTOR GEOMETRY AND VELOCITY STRUCTURE
|
|
- Carmella Golden
- 5 years ago
- Views:
Transcription
1 JOURNAL OF BALKAN GEOPHYSICAL SOCIETY, Vol.8, No 3, Augut 005, p NON-LINEAR INVERSION OF SEISMIC DATA FOR THE DETERMINATION OF REFLECTOR GEOMETRY AND VELOCITY STRUCTURE Panteli M. Soupio (), Contantino B. Papaacho (), Gregory N. Toka (), Antoni Vafidi (3) () Laboratory of Geophyic & Seimology, Department of Natural Reource and Environment, Technological Educational Intitute of Crete, Chania, Greece () Geophyical Laboratory, Aritotle Univerity of Thealoniki, P.O. Box , Macedonia, Greece (3) Applied Geophyic Laboratory, Department of Mineral Reource Engineering, Technical, Univerity of Crete, Chania, Greece Abact: A nonlinear method i ued to compute firt eimic arrival and lection avel time for velocity model with complex lector geomey. The method ue a combination of raction and lection avel-time for imultaneou determination of velocity and interface depth of the model. The elatic wave aumed to be anmitted or lected at interface at which the raypath atify Snell law. The avel time of critically racted wave and derivative maix are computed by applying a revied ray bending method, upplemented by an approximate computation of the firt Frenel one at each point of the ray. Travel-time and raypath of lected wave are computed by applying a fat finite-difference cheme baed on a olution to the eikonal equation and following the avel-time gradient backward from the receiver to the ource, repectively. A damped leat quare inverion cheme i ued to reconuct the velocity model above the lector and the geomey of the interface by minimiing the difference between oberved and calculated (direct, racted and lected) avel-time. In order to reduce inverion artifact both damping and moothing regulariation factor are applied. The applicability of the propoed method i teted uing ynthetic data. Thi imultaneou inverion cheme i appropriate for tatic eimic correction, a well a determination of lateral velocity variation and lector geomey. INTRODUCTION Travel-time tomography i primarily ued by eimologit to tudy the 3-D velocity variation in the upper crut (e.g Papaacho and Nolet, 997). Thi method wa firt inoduced by Aki, Chritofferon and Huebye (977) who divided a horiontally layered media into a Carteian grid of cell with value of lowne aigned at each cell and calculated the lowne diibution by inverting the oberved avel time. The next tep wa the application of the ame method in problem with two type of unknown parameter, uch a, the earthquake parameter and the velocity ucture (e.g Aki and Lee, 976; Thurber, 983). The tomography method of Aki et al. (977) wa adapted to lection data by Bihop et al. (985) and Farra and Madariaga (988). Reflection tomography ha been ued for eimic data analyi primarily by peoleum companie, proceing a large number of recording for different ource - receiver location to etimate the uburface velocity ucture and lector location. Until very recently, application of avel-time inverion method to thee two different type of data (firt break and lection aveltime) had been carried out 74
2 Panteli M. Soupio et al independently. Recently, raction and lection data have been proceed imultaneouly in order to obtain a better knowledge of the tudy area (Zhang et. al., 998; Zelt et. al., 999) by applying imultaneou eimic raction and lection tomography for determination of hallow and/or deep velocity and interface ucture. In lection tomography, we can parameterie the lector and velocity field in ome uitable fahion and imultaneouly invert them. Thi approach i very atactive ince it eat both inverion parameter on the ame bai (Bihop et al., 985; Farra and Madariaga, 988; Hobro et al., 003; Trink et al., 005). However, becaue of the ade-off between media velocity and lector depth (the velocity-depth ambiguity e.g Stork and Clayton, 986; Bickel, 990; Line, 993; Ro, 994; Tieman, 994) the lector interface and velocity cannot be uniquely determined without a priori information on lector location. It i well known that raction tomography provide better longwavelength velocity information, with no or very little conol on the configuration of lecting interface. Theore, if both raction and lection data (wide-angle, normalincidence) are available along a profile, imultaneou inverion of both dataet can provide better conol on both velocity and lector geomey. In the preent tudy we examine poible improvement of the tandard imultaneou tomography uing ray theory. For the avel time of the critically racted wave we follow the approach developed in a previou work (Soupio et al., 00). Specifically, initial ray are calculated uing a fat, approximate ray acing method, known a the ray initialier (Thurber and Ellworth, 980), which i equivalent to olving a one-dimenional approximation to the actual threedimenional problem. Thi method provide an initial etimate of the minimum-time ray path in heterogeneou media for the final ray acing with a revied ray bending technique (Moer et al., 99). Moreover, knowing, that the propertie of the ray are influenced not only by the ucture along the ray, but alo by the ucture in ome "vicinity" of the ray, we incorporate phyical ray by etimating the approximate width of firt Frenel volume which i frequency dependent (Cerveny and Soare, 99, Woodward, 99). The avel-time of the lected wave are computed by a finitedifference algorithm which incorporate propagation of wave in complex velocity model (Hole and Zelt 995). Snell law for lection i ued in the vicinity of the lecting interface. The lector model i allowed to vary moothly with depth, which mean that the depth to the interface d i,k could take any real value and it i not conained to lie at grid node. The advantage of thi method i that avel-time are computed imultaneouly for all receiver. The raypath and the depth derivative are tored in order to conuct the derivative maix. A continuou lector i alo ued and repreented by a erie of upport point whoe depth varie a function of horiontal offet. Ray are aced backward in time from receiver location down to lecting interface and then back up to ource location to determine the expected avel-time (forward modeling). The reult i a ytem of linear equation coniting of two type of parameter, correponding to lowne and lector depth. It i generally aumed (Anderon et. al., 984) that the inverion problem i overdetermined and hence, can be olved in a leatquare ene by minimiing the data mifit. Since there i frequently a 75
3 NON-LINEAR INVERSION OF SEISMIC DATA certain degree of non-uniquene or ill conditioning in the problem, the tandard approach i to conuct a "leat-ucture" olution (Franklin, 970; Tarantola and Nerceian, 984; Contable et al., 987) by conidering additional conaint and uually minimiing the model and/or one of the model derivative norm. The invere problem i olved by either the Levenberg-Marquardt (L-M) method (alo known a the damped leat quare method-levenberg, 944; Marquardt, 963), or by uing both a damping factor and a roughne term, in the form of e.g. a econd derivative moothing filter to the lowne model (Contable et al., 987). To how the efficiency of our algorithm, we conider three problem: firt, we invert avel-time to olve only for the velocity ucture, then we invert firt arrival and lection aveltime to determine only the geomey of the optimum lector and finally we ued the avel-time to imultaneouly olve for depth and velocity ucture. We apply the method uing 3D ynthetic data. The propoed algorithm i ueful for hallow eimic urvey (Soupio, 00) a well a, in the determination of deeper ucture depending on the cale of the problem. FORWARD MODELING Model parameteriation In local inverion tudie, a flat layered Earth model i commonly ued (Aki and Lee, 976). A imilar aumption wa alo adopted here becaue of the mall cale of our tudy area and the convenience of the carteian coordinate ytem, which ignificantly peed up all computation, epecially thoe concerning the 3-D ray acing. Our model i compoed of a threedimenional (3-D) grid of block of ame ie, with value of lowne aigned at the grid point. The lowne at a point which i located between node to 8 i calculated by ilinear interpolation, a given by the following equation, lowne=wf *(wfy *(wfx *n+wfx * n)+wfy *(wfx *n3+wfx *n4))+ wf*(wfy*(wfx*n5+wfx*n6)+wfy *(wfx*n7+wfx*n8)) () where n,n,,n8 denote the lowne value at the node of the grid block and wf(x,y,), are the weight. Either mooth velocity model or model with velocity dicontinuitie can be eaily ued (Hole and Zelt, 995). We alo incorporated a velocity dicontinuity interface, which i defined by a et of point connected by ilinear interpolation. The main advantage of the method i that the lector can be located anywhere in the model (Hole and Zelt, 995) and not necearily lie on the lowne grid node. Ray Tracing The ucceful olution of the forward modelling involve the accurate calculation of the avel-time for the two type of wave (racted and lected) and the etimation of the partial derivative of the avel-time with repect to both velocity and depth at the node (baed on the tarting velocity and lector model). Ray for the firt arrival are calculated uing a revied bending algorithm (Moer et al., 99), which ha uperior convergence propertie. In an effort to alleviate problem uch a convergence and execution peed aociated with the bending method, we incorporated the approximate ray acing cheme of Thurber and Ellworth (980). In thi approach, an equivalent one-dimenional (-D) ucture i calculated from the 3-D model between each ource-receiver pair. The velocitie at the grid point 76
4 Panteli M. Soupio et al that fall within the model are harmonically averaged, which i equivalent to arithmetically averaging the lowne. Thee average value are ued a layer velocitie for the onedimenional model, where the direct and all poible racted wave are conidered. The direct ray i determined uing a combination of the fale poition method and the ecant method (Pre et al., 988), while racted ray are handled uing time term (Officer, 958). The minimum avel-time path from thi ray acing routine i ued a the initial ray path gue for our bending ray acing routine (Moer et al., 99). In thi technique the ray i repreented by a erie of upport point which decribe the ray path with the ue of Beta-pline (Pereyra, 99). The avel time along the Beta-pline curve can be calculated with the apeoidal rule. Finally, we defined the width of the approximate Frenel volume for each point along the raypath from the ource to the receiver, a uggeted by Soupio et al. (00). The firt Frenel i incorporated by applying a imple algebraic manipulation on the element of the raypath maix, M. Thi operation diibute the value of the maix element which correpond to ray point to it Frenel volume. Τhe implementation of the propoed equation (Soupio et al., 00), i much more imple for the calculation of the Frenel volume, ince we only need the raypath connecting the given ource and receiver and the dominant frequency of our ignal for the wavelength etimation. The ray acing for the lected wave i performed uing a fat finite difference cheme (Hole and Zelt, 995). After pecifying the velocity model, the avel time of the incident wave to the lecting interface i firt calculated. Since wave that turn at or below the lector may arrive earlier but urely do not repreent the lected arrival, velocitie below the interface are et to be lower than velocitie above the lector (Podvin and Lecomte, 99). Uing Snell law, lected avel-time are aigned at the node immediately above the lector. The final lected avel-time i calculated by adding the avel-time from the lection point to the receiver, to the one previouly calculated. Raypath are determined and partial derivative of avel-time (maix M) are calculated from the pre-calculated lection aveltime a well a the velocity field and the lector geomey (Williamon, 990). Thi algorithm ha been modified to etimate the derivative moving backward in time. INVERSION PROCEDURE For a two-parameter model including lowne and lector depth, the avel-time deviation δt i given a follow, M 0 δ δt = or, M M δ δt M x = δt (3) where M and M are maice (coniting of lowne partial derivative) for racted and lected wave repectively, and M i a maix of depth partial derivative for lected wave, x i the correction vector which contain information about lowne and lector depth, δ and δ are the velocity and lector correction vector and δt i a vector of the difference between calculated and oberved aveltime. The maix M of equation (3) conit of the partial derivative with repect to both parameter (velocitylector) that have already been calculated. An important iue i the relative weight of the lected and 77
5 NON-LINEAR INVERSION OF SEISMIC DATA racted arrival, a well a the adeoff between the lowne, δ, and lector depth, δ, correction. In the preent work we follow the Bayeian approach (e.g. Tarantola, 987), where equation (3) i written a: lector depth,, and the lowne the lection, σ, avel-time error. Cd M Cm Cm x = Cd δt Cd M CThee m x = error Cd δtconol the relative (4) where C d and C m are the data and model a priori covariance maice and variance, ( = v /v ) where v i the velocity variance. Moreover, the data error can be decribed by the raction σ, and weighting of the two data type. Hence, in order to conol the weight of the two data type (raction-lection), we ue an additional relative-weight factor, x = Cx x are the normalied model r w, which i calculated a: perturbation. The objective function that we are effectively minimiing i raction _ weight σ σ rw = = = given below, lection _ weight σ T T S = ( δt Mx) Cd ( δt Mx) + γ ( x Cm x) σ (7) (5) If we replace C where, γ i a Lagrangian multiplier d and C m in equation (4) by the correponding variance we which conol the ade-off between the derive the following equation: normalied model-norm (meaure of the olution x length) and the normalied fit S δ M + 0 = δt between the oberved and the calculated σ σ data. S δ δ We conider the implet and M + M = δt σ σ σ uually adopted cae, where velocity and lector node have independent (8) and contant a priori error, v and, where δt contain the avel-time repectively, equation (5) can be reidual, M, M and M are anformed to: lowne and lector derivative T T T S = ( δt Mx) Cd ( δt Mx) + γ ( V xv xv + x maice. x ) If we multiply equation (8) (6) with σ and ue equation (7), we rewrite where the perturbation vector x ha equation (8) in maix notation a been plit into two vector, x v and x, follow, each containing the velocity and M x = δt lector depth perturbation to be (9) determined. Thu, the model variance i where, decribed only by the variance of the t t t L 0 0 L 0 n M M O M M M O M t K t K t K L 0 0 L 0 n t t t = M rw rw L r t w t t m rw rw L r w n M M O M t L tl tl M M O M rw rw L r t w L tl tl m rw rw L r w n 78
6 Panteli M. Soupio et al (0) ˆ δt ˆ δt M M ˆ n = xˆ, δtk = δt ˆ δt. rw ˆ δt. rw M M ˆ m δtl. rw () where m i the number of lector node, n i the number of lowne node, K i the number of raction avel-time, L i the number of lection avel-time and ŝ and ẑ denote the normalied (with repect to tandard error) perturbation of lowne and lector depth, repectively. In the current work we preent reult uing the LSQR algorithm where originally propoed by Paige and Saunder (98). The algorithm baed on a bidiagonaliation of the ytem uing Lanco method (Lanco, 950), followed by a QR decompoition to find the olution. In order to peed up convergence and tabilie the inverion everal approache can be adopted, a previouly decribed, uch a damping or moothing. In the preent work we mooth the lowne perturbation uing an appropriate Laplacian operator (Lee and Croon, 989) and the lector geomey by a econd-order difference operator (Ammon and Vidale, 990). The Laplacian operator add the following ytem of equation to the avel time inverion problem, ( 6 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ) i, j, k i, j, k i+, j, k i, j, k ˆ = () ˆ i the normalied lowne i, j, k perturbation at the (i,j,k) node and d i the grid pacing. Similarly the econd-difference operator add the ytem of equation, = 4 ˆ ˆ ˆ ˆ ˆ = ˆ i, j i, j i+, j i, j i, j+ where the (3) ˆ indicate the lector i, j depth perturbation at the (i,j) horiontal poition. In order to ue the above equation in -D problem, we imply ormulate them a, ˆ ˆ ˆ ˆ ˆ 0 and 4 i, k i, k i+, k i, k i, k+ = ˆi ˆ i ˆ i+ = 0, repectively. Moreover, in mot cae the derivative maix ha mall ingular value and i typically highly illconditioned. One approach in dealing with thi i to ue damping factor and augment the ytem of equation by a 0 d i, j+, k i, j, k i, j, k+ = 0 et of additional conaint. Since, we are uing normalied variable, we can imply write thee conain a ε I n m where ε i the damping value, I i the unit maix and the indice or er to lowne or lector depth. Typically, ε hould be et to (Franklin, 970) [or σ in our cae ince the ytem (8) wa multiplied by σ to define ytem (9)], if a priori data and model error are correct. Uually ε i et to value higher than to account for outlier, not expected by the Gauian diibution. The final augmented ytem of equation can be written a, M δtˆ ε I n 0 0 ˆ 0 ε = I m 0 ˆ λs λ D 0 (4) 79
7 NON-LINEAR INVERSION OF SEISMIC DATA where ε and ε are the damping factor, λ and λ are the moothing factor (concerning lowne and lector repectively), S i the Laplacian operator (equation ) and D i the econd difference operator (equation 3). NUMERICAL EXAMPLE In order to demonate the effectivene of the propoed technique we have teted our algorithm uing four model. In thee tet we elected a pecific velocity model for the tudy area and olved the forward problem. Then, the ynthetic racted and lected avel-time are inverted back to a velocity ucture and a pecific lector. Since our ynthetic data are error-free we expect that the input (ynthetic) model hould be recovered. Hence, comparion of the original velocity and what our algorithm recover permit u to analye the efficiency of the inverion cheme. Trying to chooe the proper regulariation parameter we ued the Picard condition and the L-curve criterion a preented by Soupio et al. (00). Thu, for all the tet decribed in thi ection, we conclude that a et of value of λ =5.0 and λ =3.0 (λ and λ are decribed above) hould be elected for the final olution ince our ynthetic data were error free. Alo, a imilar value ε=5.0 wa ued for the damping factor. In the preent work we have teted our method for the following cae: o Auming that lector i completely defined by a priori information, thu we invert only for the velocity field, o Auming that velocity field i completely defined by a priori information (velocity analyi etc.), thu we invert only for lector geomey and, o Simultaneou determination of the velocity model including the lector geomey. o A ynthetic tet with 3-D ource-receiver geomey for a complex geological ucture (anticline). Reconuction of velocity model The tarting model conit of a four-layer background velocity model where the velocity increae linearly with depth according to the following gradient and table (): 0.075Km / Km (5) and a flat lector at a depth of 70 Km (figure ). Model A ha a high velocity block (+0%) at depth from 0 to 40 Km. For the next model (B) we ued a low velocity anomaly in the ame region and with the ame amplitude (- 0%). In the lat model (C), two velocity anomalie (poitive and negative) were ued. A grid with 4 node in the x direction, in the y direction and in the direction (depth), reulted in a total of 86 velocity node (figure ). For the tet we ued 9 ource with 0 Km hot pacing and 0 geophone with geophone pacing 0 Km, reulting in *380 ray (direct, racted and lected wave). For thee ray, ynthetic avel time were calculated. In the inverion we et the r w parameter to a minimum (near ero) value in order to calculate the velocity field mainly from the racted wave. Figure () how the raypath coverage of the racted and lected wave for the tudy area. The ace of the raypath of racted wave are hown by the dotted line avelling through the model. Reflected raypath are preented by the dot on the grid line a the interection of the raypath with the grid. 80
8 Panteli M. Soupio et al Figure. Source-Receiver geomey and the velocity diibution ued for the ynthetic tet. Ray coverage correponding to three (A, B & C) different velocity model. The rectangular block define the region of velocity anomaly, which erred to the background velocity that increae with depth. The racted raypath avelling through the model from and to the urface are alo defined. The black ace on the grid line repreent the interection of the raypath of lected wave with the model grid. It i obviou that the raypath recognie the exitence of the velocity anomalie. Figure () how the final model from the optimiation-invere proce. The contour repreent the percentage (%) of the velocity anomaly compared to the background velocity model. The high velocity rectangular one i reconucted quite well. 8
9 NON-LINEAR INVERSION OF SEISMIC DATA Figure. Cro-ection of the final inverted velocity model a calculated from inverion of firt arrival and lection aveltime. The contour repreent the percent (%) of reconucted velocity model. The conained lector i preented by a olid line at 70 Km. A ditortion of the hape of the anomalie i found and it amplitude ha been underetimated. A) Poitive velocity anomaly i reconucted quite well (>50%), B) The negative velocity anomaly i poorly identified (<50%) ince the ray coverage in the area wa poor and C) The amplitude of the high velocity anomaly i better reconucted than the low velocity anomaly (Wielandt, 987). For model A, the hape of the high velocity anomaly contour i influenced by the ource - receiver geomey. Within the rectangular block, deviation of 0-4% are oberved in the reconucted model compared to the input model. Higher deviation are oberved for model B (-8%) and C (6-0%), due to poor coverage in thoe region (figure ). The high velocity anomaly i better reconucted becaue the ray method work well for poitive anomalie but i clearly inufficient in the preence of negative anomalie (Wielandt, 987). Reconuction of lection interface Uing the ame ource-receiver geomey but a narrower grid ( node in X, node in Y and node in Z) we have inverted the avel-time for 8
10 Panteli M. Soupio et al interface geomey auming that velocity field i a priori known. We ued three ideal lector (figure 3) with different geomey to check the algorithm efficiency. In all cae, a flat lector wa conidered at the depth of 70 Km for the tarting model. Thu, we had to invert for unknown lector node parameter (figure 3). In thi tet, we maximie the effect of lection derivative and minimie the raction, maximiing the r w parameter value a erred in the previou tet. Thu, only the lection avel-time calculated for each ideal lector were the input in the inverion proce. Figure (3) preent the ray coverage from the lected wave a etimated olving the forward problem with finite difference. Low raypath coverage at the edge of the model (lector) i oberved due to the ource/receiver and lector geomey. The ame problem appeared in cae B and C in the area where the lector exhibit, a yncline and a vertical fault, repectively. Thu, low quality of the reconucted image i expected in thoe place. Figure 3. Ray coverage from lection wave a calculated from the forward modelling. The raypath of lected wave are preented by the dot on to the grid line. Low raypath coverage at the edge and in the area of lector anomaly i oberved due to ource-receiver geomey and lector anomaly. 83
11 NON-LINEAR INVERSION OF SEISMIC DATA Figure (4A) how the tomographic reult for the dipping lector. At the end of third iteration, lector i reconucted by 95%. Figure 4. Reconuction of lector, by inverting lection avel time. The ideal lector i decribed by dahed line and the continuou line denote the tarting lector geomey. Solid bullet preent the final lector. A) Dipping layer. The ue model i reconucted by 95%, B) Syncline model, which i finally reconucted by more than 87% and C) Reflector with of a fault-like model, which i poorly reconucted due to low ray coverage cloe to the fault region. 84
12 Panteli M. Soupio et al Figure (4B) preent the tomographic image in cro ection for the yncline model which i ueful in the reconuction of tectonic ucture. The reult how that lector geomey i well reconucted. The reconucted model how dicrepancie from the actual one at low raypath coverage region uch a the edge and the middle of the model. For the lector with fault, ranging between the depth of 60 km and 80 km, the algorithm recover it depth quite well (figure 4C). However, ome region (edge of lector) were reconucted poorly, due to low ray coverage (Figure 3C). The r.m. error wa reduced from 8.96 (initial model coniting of a flat lector model) to.3. Reconuction of the velocity model and the lector geomey The input model i 00 Km long and 00 Km deep, with a background velocity model which conit of four layer, whoe velocity increae linearly with depth (equation 5 and table ) and a flat lector at a depth of 70 Km. The ame velocity anomaly a in figure (A) with a lower amplitude (0%) and the three initial lector of figure (3) are alo ued for the calculation of critically racted and lected wave. We have ued a narrow a well a a dene grid with /4 node in the x direction, in the y direction and / in the direction (depth), reulting in a total of 3/86 velocity node and /4 lector node. 9 ource with 0 geophone along the top of the model were ued, thu reulting in 760 (*380) avel time obervation. Figure (5) preent the raypath coverage from the imultaneou olution of the forward problem uing firt arrival and lection aveltime. Poor ray coverage i alo oberved in region near the lector. Figure (6A), (6B) and (6C) how the inverion reult. We hould note that contour repreent the percentage (%) of the reconucted velocity model. Thu, the high velocity region i reconucted quite well, ince the ray ample better the fater region. Regarding the lector, the inverion recovered it geomey quite well except for the fault model due to low ray coverage. It i important to notice that the main velocity model characteritic, uch a the amplitude and hape of the anomalie were better determined than in the previou tet. 3D ynthetic tet for the modelling in complex geological area The 3D model wa defined a 6 Km and 00 Km deep uing the ame background velocity model a preented in table (). Table. Starting background velocity model Layer Average Velocity (Km/) Depth We ued almot the ame model a in figure () with 3D ource and receiver geomey. A high velocity anomaly (0%) lying on the top of the flat initial lector and an ideal lector like a anticline geological ucture are ued for the calculation of the ynthetic avel time. The initial lector model wa a flat dicontinuity at a depth of 60 Km. We have ued a narrow 3D grid with 6 node in the x direction, 6 85
13 NON-LINEAR INVERSION OF SEISMIC DATA node in the y direction and in the direction (depth), reulting in a total of 5376 velocity node and 56 lector node. A total number of hot point (iangle) diibuted on the top of the model were recorded by 0 urface geophone (quare) a i hown in figure (7). The final data et conited of 840 firt break of P wave and lection aveltime from the initial lector. Figure (7) preent in 3D the tomographic image from the imultaneou inverion of firt arrival and lection aveltime for the lat iteration. The reult how that the initial model i well reconucted. Applying appropriate regulariation parameter ha deflected ome artefact. It wa intereting that in the area of interet both the amplitude and the hape of anomaly (velocity/lector) are ignificantly reconucted for more than 50%. CONCLUSION We preent a tomographic method for joint etimation of the velocity field and lector poition uing firt arrival and lection aveltime. We have alo hown that our algorithm i capable of producing reliable tomographic image uing ynthetic data. Thi method may be alo helpful in conucting complex geological model with local bodie, yncline, anticline and fault. The propoed method i alo able to perform an accurate determination of the hallow velocity model and correct the large time hift (tatic) in lection, reulting in better lector continuity. Moreover, we can ue tomographically determined velocitie to image eimic data applying depth migration, avoiding tacking and normal moveout (pretack migration). The propoed method incorporate the ue of an efficient ray acing technique and a finite difference code for olving the forward problem in conjuction with a combination of inverion tabiliation approache. In pite of the amount of data, the propoed method and algorithm are quite efficient, ince they require relative mall computation time. REFERENCES Aki, K., and W. H. K. Lee, 976, Determination of threedimenional velocity anomalie under a eimic array uing firt p arrival time from local earthquake:, A homogeneou initial model: J. Geophy. Re., 8, Aki, K., Chritofferon, A., and Huebye, E. S., 977, Determination of the three dimenional eimic ucture of the lithophere: J. Geophy. Re., 8, Anderon, D. L., and Diewonki, A. M., 984, Seimic Tomography: Scientific American 5, No 0, Ammon, C. J., and Vidale, J. E., 990, Seimic avel-time tomography uing combinatorial optimiation technique: Sei. Re. Lett., 6, 39. Bickel S. H., 990, Velocity-Depth ambiguity of lection aveltime: Geophyic, 55, Bihop, T. N., Bube, K. P., Cutler, R. T., Langan, R. T., Love, P. L., Renick, J. R., Shuey, J. T., Spindler, D. A. and Wyld, H. W., 985, Tomographic determination 86
14 Panteli M. Soupio et al of velocity and depth in laterally - varying media: Geophyic, 50, Cerveny, V., Soare E. P., 99, Frenel volume ray acing: Geophyic, 57, Contable, S. C., Parker, R. L., and C. Contable, 987, Occam' inverion: A practical algorithm for generating mooth model from elecomagnetic ounding data: Geophyic, 5, Farra, V. and Madariaga, R., 988, Non-linear lection tomography: Geophyic. J., 95, Franklin, J. N., 970, Well-poed tochatic extenion of ill-poed linear problem: J. Math. Anal. Appl., 3, Hobro, J.W.D., Singh, S.C. and Minhull, T.A., 003, Threedimenional tomographic inverion of combined lection and raction eimic aveltime data: Geophyical Journal International, 5, Hole, J. A., and Zelt, B. C., 995, 3-D finite-difference lection aveltime: Geophy. J. Int.,, Lanco, C., 950, An iteration method for the olution of the eigenvalue problem of linear differential and integral operator: J. Re. Nat. Bur. Stand., 45, Lee, J. M., and Croon, R. S., 989, Tomographic inverion for three dimenional velocity ucture at Mount St. Helen uing Earthquake data: JGR, 94, B5, Levenberg, K., 944, A method for the olution of certain nonlinear problem in leat-quare: Quart. Appl. Math.,, Line, L., 993, Ambiguity in analyi of velocity and depth: Geophyic, 58, Marquardt, D. W., 963, An algorithm for leat-quare etimation of nonlinear parameter: Soc. Indu. Appl. Math., J. Appl. Math.,, Moer, T. J., Nolet, G. and Snieder, R., 99, Ray Bending Reviited: Bulletin of the Seimological Society of America, 8, Officer, C. B., 958, Inoduction to the Theory of Sound Tranmiion: McGraw-Hill, New York, pp. 84. Paige, C. C., and Saunder, M. A., 98, LSQR: An algorithm for pare linear equation and pare leat quare: A.C.M. Tran. Math. Software., 8, Papaacho, C.B., and Nolet, G., 997, P and S deep velocity ucture of Hellenic area obtained by robut non-linear inverion of arrival time: J. Geophy. Re., 0, Pereyra, V., 99, Two-point ray acing in general 3D media: Geophyical Propecting, 40, Podvin, P. and Lecomte, I., 99, Finite difference computation of aveltime in very conated velocity model: a maively parallel approach and it aociated tool: Geophy. J. Int., 05, Pre, W. H., B. P. Flannery, S. A. Teukolky, and W. T. Vetterling, 988, Numerical Recipe in Foran - The Art of Scientific Computing: Cambridge Univerity Pre, Cambridge, pp. 9, Ro, W. S., 994, The velocity-depth ambiguity in eimic avel-time data, Geophyic, 59, Soupio, P.M., Papaacho, C.B, Juhlin, C., Toka, G.N., 00, Nonlinear three-dimenional aveltime inverion of crohole data with an application in the area of 87
15 NON-LINEAR INVERSION OF SEISMIC DATA middle Ural, Geophyic, 66, Soupio, P.M., 00, Implemenation of tomographical Method for the Determination of Geotechnical Parameter in the Bridge of Neto River: Technical Report, 45 pp (in Greek). Stork, C. and Clayton, R. W., 986, Analyi of the reolution between ambiguou velocity and lector poition for avel-time tomography: Extended Abact of the 56th SEG International Meeting, Tarantola, A., 987, Invere Problem Theory: Elevier, Amterdam. Tarantola, A., and A. Nerceian, 984, Three-dimenional inverion without block: Geophy. J. R. A. Soc., 79, pp Thurber, C. H., 983, Earthquake location and three-dimenional crutal ucture in the Coyote lake area, Cenal California, J. Geophy. Re., 88, Thurber, C. H., and Ellworth, W. L., 980, Rapid olution of ray acing problem in heterogeneou media: Bulletin of the Seimological Society of America, 70, Tieman, H. J., 994, Invetigating the velocity-depth ambiguity of lection avel-time: Geophyic, 59, Trink, I, Singh, S. C., Chapman, C. H., Barton, P. J., Boch, M. and Cherrett, A., 005, Adaptive aveltime tomography of denely ampled eimic data: Geophyical Journal International, 60, Wielandt, E., 987, On the validity of the ray approximation for interpreting delay time: Seimic Tomography (Guut Nolet), D. Reidel Publihing Company. Williamon, P. R., 990, Tomographic inverion in lection eimology, Geophy. J. Int., 00, Woodward, M. J., 99, Wave-equation tomography: Geophyic, 57, 5-6. Zelt, C. A., A. M. Hojka, E. R. Flueh, and K. D. McIntoh, 999, 3D imultaneou eimic raction and lection tomography of wide-angle data from the Cenal Chilean margin: Geophy. Re. Lett., 6, Zhang, J., U.S. ten Brink, M. N., Toko, 998, Nonlinear raction and lection aveltime tomography: JGR-Solid Earth, 03, B,
Downloaded 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 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 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 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 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 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 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 informationChapter 13 Non Sampling Errors
Chapter 13 Non Sampling Error It i a general aumption in the ampling theory that the true value of each unit in the population can be obtained and tabulated without any error. In practice, thi aumption
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 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 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 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 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 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 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 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 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 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 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 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 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 informationModeling and Analysis of Slow CW Decrease for IEEE WLAN
Modeling and Analyi of Slow CW Decreae for IEEE 82. WLAN Qiang Ni, Imad Aad 2, Chadi Barakat, and Thierry Turletti Planete Group 2 Planete Group INRIA Sophia Antipoli INRIA Rhône-Alpe Sophia Antipoli,
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 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 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 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 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 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 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 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 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 informationThe norm Package. November 15, Title Analysis of multivariate normal datasets with missing values
The norm Package November 15, 2003 Verion 1.0-9 Date 2002/05/06 Title Analyi of multivariate normal dataet with miing value Author Ported to R by Alvaro A. Novo . Original by Joeph
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 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 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 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 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 informationKey Terms - MinMin, MaxMin, Sufferage, Task Scheduling, Standard Deviation, Load Balancing.
Volume 3, Iue 11, November 2013 ISSN: 2277 128X International Journal of Advanced Reearch in Computer Science and Software Engineering Reearch Paper Available online at: www.ijarce.com Tak Aignment in
More informationA Boyer-Moore Approach for. Two-Dimensional Matching. Jorma Tarhio. University of California. Berkeley, CA Abstract
A Boyer-Moore Approach for Two-Dimenional Matching Jorma Tarhio Computer Science Diviion Univerity of California Berkeley, CA 94720 Abtract An imple ublinear algorithm i preented for two-dimenional tring
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 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 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 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 informationMarkov Random Fields in Image Segmentation
Preented at SSIP 2011, Szeged, Hungary Markov Random Field in Image Segmentation Zoltan Kato Image Proceing & Computer Graphic Dept. Univerity of Szeged Hungary Zoltan Kato: Markov Random Field in Image
More informationA Linear Interpolation-Based Algorithm for Path Planning and Replanning on Girds *
Advance in Linear Algebra & Matrix Theory, 2012, 2, 20-24 http://dx.doi.org/10.4236/alamt.2012.22003 Publihed Online June 2012 (http://www.scirp.org/journal/alamt) A Linear Interpolation-Baed Algorithm
More informationComparison of Methods for Horizon Line Detection in Sea Images
Comparion of Method for Horizon Line Detection in Sea Image Tzvika Libe Evgeny Gerhikov and Samuel Koolapov Department of Electrical Engineering Braude Academic College of Engineering Karmiel 2982 Irael
More informationCSE 250B Assignment 4 Report
CSE 250B Aignment 4 Report March 24, 2012 Yuncong Chen yuncong@c.ucd.edu Pengfei Chen pec008@ucd.edu Yang Liu yal060@c.ucd.edu Abtract In thi project, we implemented the recurive autoencoder (RAE) a decribed
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 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 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 informationVariable Resolution Discretization in the Joint Space
Variable Reolution Dicretization in the Joint Space Chritopher K. Monon, David Wingate, and Kevin D. Seppi {c,wingated,keppi}@c.byu.edu Computer Science, Brigham Young Univerity Todd S. Peteron peterto@uvc.edu
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 informationDiffraction imaging and time-migration velocity analysis using oriented velocity continuation a
Diffraction imaging and time-migration velocity analyi uing oriented velocity continuation a a Publihed in Geophyic, 82, U25-U35, (2017) Luke Decker, Dmitrii Merzlikin, and Sergey Fomel ABSTRACT We perform
More informationAnalysis of slope stability
Engineering manual No. 8 Updated: 02/2016 Analyi of lope tability Program: Slope tability File: Demo_manual_08.gt In thi engineering manual, we are going to how you how to verify the lope tability for
More informationP-SV Wave-Field Connection Technique for Regional Wave Propagation Simulation
Bulletin of the Seimological Society of America, Vol. 95, No. 4, pp. 1375 1386, Augut 2005, doi: 10.1785/0120040173 P-SV Wave-Field Connection Technique for Regional Wave Propagation Simulation by Zengxi
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 informationA SIMPLE IMPERATIVE LANGUAGE THE STORE FUNCTION NON-TERMINATING COMMANDS
A SIMPLE IMPERATIVE LANGUAGE Eventually we will preent the emantic of a full-blown language, with declaration, type and looping. However, there are many complication, o we will build up lowly. Our firt
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 informationPerformance Evaluation of an Advanced Local Search Evolutionary Algorithm
Anne Auger and Nikolau Hanen Performance Evaluation of an Advanced Local Search Evolutionary Algorithm Proceeding of the IEEE Congre on Evolutionary Computation, CEC 2005 c IEEE Performance Evaluation
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 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 informationAnalyzing Hydra Historical Statistics Part 2
Analyzing Hydra Hitorical Statitic Part Fabio Maimo Ottaviani EPV Technologie White paper 5 hnode HSM Hitorical Record The hnode i the hierarchical data torage management node and ha to perform all the
More informationMulti-Target Tracking In Clutter
Multi-Target Tracking In Clutter John N. Sander-Reed, Mary Jo Duncan, W.B. Boucher, W. Michael Dimmler, Shawn O Keefe ABSTRACT A high frame rate (0 Hz), multi-target, video tracker ha been developed and
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 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 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 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 informationImage Morphing Using Deformation Techniques
Image Morphing Uing Deformation Technique Seung-Yong Lee, Kyung-Yong Chwa, Jame Hahn, and Sung Yong Shin Department of Computer Science Korea Advanced Intitute of Science and Technology 373-1 Kuong-dong
More informationImage authentication and tamper detection using fragile watermarking in spatial domain
International Journal of Advanced Reearch in Computer Engineering & Technology (IJARCET) Volume 6, Iue 7, July 2017, ISSN: 2278 1323 Image authentication and tamper detection uing fragile watermarking
More informationAnalytical Redundancy and Fuzzy Inference in AUV Fault Detection and Compensation
Analytical Redundancy and Fuzzy Inference in AUV Fault Detection and Compenation A. J. Healey Profeor and Director Center for AUV Reearch Naval Potgraduate School Monterey, CA 93953 healey@me.np.navy.mil
More informationP063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA
1 P063 ADAPTIVE TRAVEL-TIME AND RAY- PARAMETER INVERSION OF DENSELY SAMPLED 2-D SEISMIC DATA I. TRINKS 1, S. SINGH 2, C. CHAPMAN 3, P. BARTON 1, M. BOSCH 4, A. CHERRETT 5 Addresses 1 Bullard Laboratories,
More informationHassan Ghaziri AUB, OSB Beirut, Lebanon Key words Competitive self-organizing maps, Meta-heuristics, Vehicle routing problem,
COMPETITIVE PROBABIISTIC SEF-ORGANIZING MAPS FOR ROUTING PROBEMS Haan Ghaziri AUB, OSB Beirut, ebanon ghaziri@aub.edu.lb Abtract In thi paper, we have applied the concept of the elf-organizing map (SOM)
More informationUsing Mouse Feedback in Computer Assisted Transcription of Handwritten Text Images
2009 10th International Conference on Document Analyi and Recognition Uing Moue Feedback in Computer Aited Trancription of Handwritten Text Image Verónica Romero, Alejandro H. Toelli, Enrique Vidal Intituto
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 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 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 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 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 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 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 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 informationAn Algebraic Approach to Adaptive Scalable Overlay Network Monitoring
An Algebraic Approach to Adaptive Scalable Overlay Network Monitoring ABSTRACT Overlay network monitoring enable ditributed Internet application to detect and recover from path outage and period of degraded
More informationStress-Blended Eddy Simulation (SBES) - A new Paradigm in hybrid RANS-LES Modeling
Stre-Blended Eddy Simulation (SBES) - A new Paradigm in hybrid RANS-LES Modeling Menter F.R. ANSYS Germany GmbH Introduction It i oberved in many CFD imulation that RANS model how inherent technology limitation
More informationWeb Science and additionality
Admin tuff... Lecture 1: EITN01 Web Intelligence and Information Retrieval Meage, lide, handout, lab manual and link: http://www.eit.lth.e/coure/eitn01 Contact: Ander Ardö, Ander.Ardo@eit.lth.e, room:
More informationComputer Arithmetic Homework Solutions. 1 An adder for graphics. 2 Partitioned adder. 3 HDL implementation of a partitioned adder
Computer Arithmetic Homework 3 2016 2017 Solution 1 An adder for graphic In a normal ripple carry addition of two poitive number, the carry i the ignal for a reult exceeding the maximum. We ue thi ignal
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 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 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 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 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 informationInterface Tracking in Eulerian and MMALE Calculations
Interface Tracking in Eulerian and MMALE Calculation Gabi Luttwak Rafael P.O.Box 2250, Haifa 31021,Irael Interface Tracking in Eulerian and MMALE Calculation 3D Volume of Fluid (VOF) baed recontruction
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 information3D MODELLING WITH LINEAR APPROACHES USING GEOMETRIC PRIMITIVES
MAKARA, TEKNOLOGI, VOL. 9, NO., APRIL 5: 3-35 3D MODELLING WITH LINEAR APPROACHES USING GEOMETRIC PRIMITIVES Mochammad Zulianyah Informatic Engineering, Faculty of Engineering, ARS International Univerity,
More informationDistributed Partial Information Management (DPIM) Schemes for Survivable Networks - Part II
IEEE INFOCO 2002 1 Ditributed Partial Information anagement (DPI) Scheme for Survivable Network - Part II Dahai Xu Chunming Qiao Department of Computer Science and Engineering State Univerity of New York
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 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 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 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 informationLecture Outline. Global flow analysis. Global Optimization. Global constant propagation. Liveness analysis. Local Optimization. Global Optimization
Lecture Outline Global flow analyi Global Optimization Global contant propagation Livene analyi Adapted from Lecture by Prof. Alex Aiken and George Necula (UCB) CS781(Praad) L27OP 1 CS781(Praad) L27OP
More informationAccelerating seismic interpolation with a gradient projection method based on tight frame property of curvelet
CSIRO PUBLISHING Exploration Geophyic, 215, 46, 253 26 http://dx.doi.org/1.171/eg1416 Accelerating eimic interpolation with a gradient projection method baed on tight frame property of curvelet Jingjie
More information