Recursive integral time extrapolation methods for scalar waves Paul J. Fowler*, Xiang Du, and Robin P. Fletcher, WesternGeco
|
|
- Debra Jackson
- 5 years ago
- Views:
Transcription
1 Reursive integral time extrapolation methods for salar waves Paul J. Fowler*, Xiang Du, and Robin P. Flether, WesternGeo Summary We derive and ompare a variety of algorithms for reursive time extrapolation of salar waves using approximate operators derived from integral solutions of a wave equation. These methods fall into two ategories: those based on ombining or interpolating homogeneous solutions, and those based on series expansions of heterogeneous operators. The former suffer from osillatory noise at large veloity disontinuities unless the time step is small, whereas the latter allow aurate extrapolation at large time steps, but are more ostly. Introdution Reverse time migration requires high-auray time extrapolation of wavefields. Conventional expliit finitedifferene or pseudospetral methods for time extrapolation of waves are usually limited to small time steps by numerial dispersion errors and instability. Reently, several alternative approahes have been proposed for solving salar wave equations based on integral mixeddomain spae / wavenumber time extrapolation methods (Zhang et al., 007; Soubaras and Zhang, 008; Wards et al., 008; Zhang and Zhang, 009; Etgen and Brandsberg- Dahl, 009; Pestana and Stoffa, 009; Stoffa and Pestana, 009; Liu, et al., 009; Du, et al., 010). We refer here to this entire family of approahes generially as reursive integral time extrapolation (RITE), as they are all based on formal integral solutions of a wave equation. These integral solutions are then applied reursively to extrapolate the wavefield forward or bakwards in disrete time steps. These RITE approahes attempt to ahieve stable, dispersion-free extrapolation in heterogeneous media for large time steps, up to the yquist limit or beyond. We ategorize and ompare here a variety of RITE methods within a ommon derivation framework. RITE methods are oneptually based on a heuristi of time extrapolating a wavefield using multipliation by a phase i t shift of the form e ω, where t is time and ω has units of frequeny and is derived from an appropriate dispersion relation. An analogous onept of spatial wavefield extrapolation using multipliation by a phase shift of the ikz z form e, where z is depth and k z is vertial wavenumber, has been widely used previously as a basis for downwardontinuation modeling and migration methods (e.g., Gazdag, 1978). Many RITE algorithms resemble formal translations of existing downward-ontinuation methods. RITE operators for salar aousti waves The onstant-density aousti salar wave equation an be written as: + A u( t) = 0. (1) t The operator A is defined in the spae domain as A[ u( x, t)] = v (x) u( x, t), () or in the wavenumber domain as A[ u( k, t)] = dk v (k k) k u( k, t). (3) A formal integral time extrapolation operator for the salar wave equation 1 an be written as iφ t u( t + t) = e u( t), (4) where Φ A. This form of the extrapolator requires omplex-valued data. Adding solutions from equation 4 for u( t + t) and u( t t) gives a time extrapolator for real-valued data of the form u( t + t) = os( Φ t) u( t) u( t t). (5) We refer to extrapolators of the form of equation 4 as onestep methods, and those from equation 5 as two-step methods. Four-step or higher in time methods an be developed similarly. The square-root, exponential, and osine funtions of operators indiated in equations 4 and 5 an be given formal meaning by deomposition into eigenvetors and eigenvalues of the wavenumber-domain operator A defined in equation 3. This in priniple gives integral solutions for free-spae time extrapolation of the salar wave equation 1 that will be aurate in heterogeneous media for arbitrarily large time steps. This is a time extrapolation analog of the eigen-deomposition downward-ontinuation approah disussed by Jaobs (1980), Pai (1985), and Kosloff and Kessler (1987). However, eigen-deomposition of the operator A is omputationally too expensive for realisti three-dimensional models. The various approahes suggested for pratial approximation of the heterogeneous RITE operator in equation 4 an be divided into two basi families: those based on ombination or interpolation of loally homogeneous solutions, and those based on series expansion approximations of the heterogeneous operator. Methods for ombining homogeneous RITE operators For homogeneous media, the operator A redues to ust A = v k and the one-step extrapolator beomes u k t + t = e u k t iv k t (, ) (, ). (6) 310
2 Reursive integral time extrapolation of waves This operator is now easily implemented in the wavenumber domain, and is essentially free from numerial dispersion, even for arbitrarily large time steps. One simple way to implement a variable-veloity RITE sheme for heterogeneous media would be to design separate onstant-veloity extrapolators for every point in spae, using the partiular veloity appropriate to eah individual loation. We term suh brute fore solutions phase-shift plus seletion (PSPS). Simple PSPS is too expensive to be of muh pratial use. Wards et al. (008) and Etgen and Brandsberg-Dahl (009) suggested less-expensive approahes based on extrapolating a wavefield using only a few disrete onstant veloities, followed by interpolation between these data volumes. We term suh methods phase-shift plus interpolation (PSPI) in analogy with the downward-ontinuation method introdued by Gazdag and Sguazzero (1984). PSPI methods an also be based on the two-step onstantveloity extrapolator u( t + t) = os( v k t) u( t) u( t t). (7) Etgen and Brandsberg-Dahl (009) suggest modifying the osine extrapolator to the form [os( v k t) 1] / ( v t ), alling this normalized operator a pseudo-laplaian. Comparing the Taylor expansions of the operators shows that the pseudo-laplaian normalization anels the leading order veloity dependene. This normalization an thus inrease interpolation auray for small time steps. Zhang and Zhang (009) introdue another sheme for ombining onstant-veloity RITE operators based on earlier downward-ontinuation methods used by Chen and Liu (006). One an build a matrix of one-step onstant-veloity extrapolators mn iv k t m n = e, defined for ordered sets of disrete veloity values { v } and disrete m wavenumber magnitudes { k n }. This matrix an then be deomposed by singular value deomposition into a T produt of three matries of the form =ΞΛΨ. Λ is a diagonal matrix of the orresponding singular values λ. One an then write the elements of matrix as the separable deomposition iv t λ = m n e k λ ξ ( v ) ψ ( k ), (8) = 1 m n where ξ ( v ) = Ξ is a olumn of the matrix Ξ, m m ψ ( k ) = Ψ a olumn of the matrix Ψ, and λ is the n n number of singular values used. In pratie, Zhang and Zhang (009) suggest using the optimal separable approximation (OSA) to reate a parsimonious approximation using only the largest singular values. Given then a method for interpolating the funtions ξ and ψ, this separable deomposition an be extended to represent extrapolators for any value within the desired ranges, yielding an extrapolation sheme of the form { ψ } ( x, ) 1 λ ξ [ ( x)] [ ( k)] ( k, ), (9) u t + t = v FT k u t = 1 where FT indiates a 3-D spatial Fourier transform. Series expansions of heterogeneous RITE operators The approximation methods disussed in the previous setions onstrut a heterogeneous solution using loally homogeneous salar-phase solutions, with no diret use of heterogeneous extrapolators. Series expansions provide one general method for deriving approximations to the wavenumber-onvolutional time extrapolation operators. Series expansions are easier to onstrut for two-step operators, beause only even powers of the phaseφ are needed. Suppose that one has a set of basis polynomials { P ( x) p x } i =. Approximating the two-step osine i = 1 operator as an expansion in these basis polynomials gives where os( Φ t) b P ( Φ t) = a Φ t, (10) a = 1 = 0 = 0 = b p. The series in equation 10 an be i i i evaluated in nested form as a a 1 os( Φ t) a + t + t ( + ). (11) 0 I Φ I Φ I a a 0 1 In this form, evaluation of eah additional term requires multipliation of the data by k in the wavenumber domain, followed by transformation bak to the spae domain and multipliation by the veloity squared, followed by transformation bak to the wavenumber domain, and so forth. ote that trunating suh an evaluation after theφ term yields a onventional seondorder in time, pseudo-spetral in spae algorithm. The exat riterion for operator stability requires omputation of a full eigen-deomposition, whih is usually infeasible. However, a good estimate an be ahieved using the orresponding stability riteria for onstant-veloity operators for the largest values of veloity, wavenumber, and time-step size. For the two-step osine extrapolator, Soubaras and Zhang (008) give the approximate stability riterion that the osine approximation used must satisfy os( φ) 1 on the 311
3 Reursive integral time extrapolation of waves interval [ φ, φ ], where φ = π v t / r and max max max ma x 1/ r = ( x + y + z ). Soubaras and Zhang (008) also present a similar auray riterion, requiring good fitting of the osine up to a value given by φ = π t f. For time stepping at yquist, this redues to ust φ = π. Perhaps the simplest expansion is a Taylor series, but they an onverge very slowly, and so are not omputationally very effiient. Expansions using Legendre or Chebyshev polynomials an provide more effiient and aurate approximations. Very large time steps an be used if enough terms in these expansions are retained. Chebyshev polynomial expansions have been used for downward ontinuation by Kosloff and Kessler (1987) and for time extrapolation by Tal-Ezer et al. (1987), Kosloff et al. (1989), Pestana and Stoffa (009), and Stoffa and Pestana (009). These orthogonal polynomial expansions fit the osine by distributing errors evenly over the entire stability interval. The required auray interval is usually muh smaller. One an instead numerially ompute polynomial oeffiients that optimize the fit over the interval [ φ, φ ], while onstraining the optimization to fore stability over the larger interval. Soubaras and Zhang (008) desribe one suh optimal polynomial two-step RITE approah, using the Remez exhange algorithm to derive optimized oeffiients. This is an extension of the downwardontinuation extrapolators presented by Soubaras (1996). We use here an alternative onstrained weighted leastsquares approah (WLSQ) based on the downwardontinuation method of Thorbeke et al. (004). Examples In Figure 1 we show a series of numerial experiments omparing four different RITE implementations. The model for eah has an upper layer with v=1500 m/s and a lower layer with v=4500 m/s, with an abrupt ump between layers. The spatial sampling is 15 m. A point soure with a maximum frequeny of 50 Hz is loated in the upper layer, and eah image shows a wavefield snapshot at t=0.8s. The olumns show omputations using time-step values t of s, s, 0.01 s, and 0.0 s, equivalent to 0.5, 0.5, 1.0, and.0 times the yquist time step. The first row (a-d) uses a one-step PSPI operator, the seond (e-h) a two-step PSPI operator, the third (i-l) a onestep OSA operator with 6 eigenvalues, and the fourth (m-p) a two-step WLSQ series expansion operator, using 13 terms for the first three images and 18 terms for the last. All the methods give good results for the smallest time step ( t = s ). The one-step PSPI operator in the first row gives good results for time steps up to t = s max but displays osillatory artifats for larger time steps. The two-step PSPI method in the seond row has larger artifats, even for time step t = 0.005s. The OSA results in the third row are nearly idential to the one-step PSPI in the first row. ote here that the OSA and PSPI methods essentially redue to PSPS beause only two veloity values are present, so no atual interpolation is being done. The osillatory artifats seen in the OSA and PSPI results arise from the attempt to ombine homogeneous operators. The two-step WLSQ method in the fourth row avoids these artifats, even for large time steps, by more aurately approximating the orret heterogeneous operator. Further extensions of RITE methods RITE methods an also be implemented in the spae domain only. For PSPS, PSPI, and OSA RITE methods, one an design the extrapolation operators in the wavenumber domain, but transform them bak to beome onvolutional spae-domain operators. For PSPS, this yields an analog of the downward-ontinuation approah of Blaquiere et al. (1988). These spae-domain operators an be preomputed and applied repeatedly at every time step, but the storage requirements an beome very large. Similarly, two-step series expansion methods an be implemented in the spae domain by replaing the wavenumber-domain operator in equation 11 with a finitedifferene spatial Laplaian. Our study here is not exhaustive. Using rational fration expansions instead of series to approximate operators leads to impliit extrapolation shemes (e.g., Claerbout, 1985). For time extrapolation, the resulting sets of impliit equations appear to beome extremely large. Fourier finite-differene (FFD) (Ristow and Ruhl, 1994) and generalized sreen propagator (e.g., Le Rousseau and de Hoop, 1998) methods are based on omputing perturbations from onstant-veloity referene operators (Sava, 000). For time extrapolation, these methods will need to handle a muh larger range of veloities than for downward ontinuation (a problem shared also by PSPI methods), and so an be expeted to be limited to small time steps. However, Song and Fomel (pers. omm.) reently report promising results for a FFD time-stepping method. One an also write a salar equation for P and SV waves in VTI media (Fowler et al., 010), and find formal singlemode exponential operator solutions analogous to equation 4, allowing extension of the various RITE methods disussed here to either P or SV waves in TI media. Etgen and Brandsberg-Dahl (009) disuss TI PSPI methods. Extensions of OSA methods to VTI media were presented by Liu et al. (009) and by Du et al. (010). Two-step osine series expansion extrapolators an also be generalized to VTI or TTI media, but the inrease in omputational ost an beome onsiderable. 31
4 Reursive integral time extrapolation of waves Conlusions The RITE methods we have disussed here allow extension for use in reverse time migration of many of the methods used previously for depth extrapolation. We examined two maor families of mixed-domain RITE methods: ones that ombine homogeneous extrapolators, and ones based on series approximations to heterogeneous extrapolators. The latter allow aurate extrapolation with larger time steps than the former, and exhibit fewer artifats at large abrupt veloity disontinuities. However, this improved auray usually requires muh greater omputational effort per time step, espeially for anisotropi media. a) b) ) d) e) f) g) h) i) ) k) l) m) n) o) p) Figure 1: Wavefield snapshots at time 0.8 s for a two-layer model. The upper layer veloity is 1500 m/s and bottom layer veloity is 4500 m/s. Images from left to right represent time steps of t = s, t = s, t = 0.01 s, and t = 0.0 s. The top row (a-d) uses a one-step PSPI approah. The seond row (e-h) uses a two-step PSPI approah. The third row (i-l) uses a one-step OSA approah. The bottom row (m-p) uses a twostep WLSQ series expansion approah. 313
5 EDITED REFERECES ote: This referene list is a opy-edited version of the referene list submitted by the author. Referene lists for the 010 SEG Tehnial Program Expanded Abstrats have been opy edited so that referenes provided with the online metadata for eah paper will ahieve a high degree of linking to ited soures that appear on the Web. REFERECES Blaquiere, G., H. W. J. Debeye, C. P. A. Wapenaar, and A. J. Berkhout, 1988, 58th Annual International Meeting, SEG, Expanded Abstrats, Chen, J., and H. Liu, 006, Two kinds of separable approximations for the one-way wave operator: Geophysis, 71, no. 1, T1 T5, doi: / Claerbout, J. F., 1985, Imaging the Earth s Interior: Blakwell Sientifi Publiations, In. Du, X., R. P. Flether, and P. J. Fowler, 010, Pure P-wave propagators versus pseudo-aousti propagators for RTM in VTI media: 7nd EAGE Conferene & Exhibition, Extended Abstrats. Etgen, J., and S. Brandsberg-Dahl, 009, The pseudo-analytial method: appliation of pseudo-laplaians to aousti and aousti anisotropi wave propagation: 79th Annual International Meeting, SEG, Expanded Abstrats, Fowler, P. J., X. Du, and R. P. Flether, 010, Coupled equations for reverse time migration in transversely isotropi media : Geophysis, 75, no. 1, S11 S, doi: / Gazdag, J., 1978, Wave equation migration with the phase-shift method: Geophysis, 43, , doi: / Gazdag, J., and P. Sguazzero, 1984, Migration of seismi data by phase shift plus interpolation: Geophysis, 49, , doi: / Jaobs, A., 1980, Wide-angle variable veloity one-way wave-equation modeling: Presented at the 50th Annual International Meeting, SEG. Kosloff, D., and D. Kessler, 1987, Aurate depth migration by a generalized phase-shift method: Geophysis, 5, , doi: / Kosloff, D., A. Q. Filho, E. Tessmer, and A. Behle, 1989, umerial solution of the aousti and elasti wave equation by a new rapid expansion method: Geophysial Prospeting, 37, no. 4, , doi: / tb01.x. Le Rousseau, J. H., and M. V. de Hoop, 1998, Modeling and imaging with the generalized sreen algorithm: 68th Annual International Meeting, SEG, Expanded Abstrats, Liu, F., S. A. Morton, S. Jiang, L. i, and J. P. Leveille, 009, Deoupled wave equations for P and SV waves in an aousti VTI media: 79th Annual International Meeting, SEG, Expanded Abstrats, Pai, D. M., 1985, A new solution method for the wave equation in inhomogeneous media : Geophysis, 50, , doi: / Pestana, R., and P. L. Stoffa, 009, Rapid expansion method (REM) for time-stepping in reverse time migration (RTM): 79th Annual International Meeting, SEG, Expanded Abstrats, Ristow, D., and T. Ruhl, 1994, Fourier finite-differene migration: Geophysis, 59, , doi: /
6 Sava, P. C., 000, A tutorial on mixed-domain wave-equation migration and migration veloity analysis: Stanford Exploration Proet Report, 105, Soubaras, R., 1996, Expliit 3-D migration using equiripple polynomial expansion and Laplaian synthesis: Geophysis, 61, , doi: / Soubaras, R., and Y. Zhang, 008, Two-step expliit marhing method for reverse-time migration: 70th EAGE Conferene & Exhibition, Extended Abstrats. Stoffa, P. L., and R. Pestana, 009, umerial solution of the aousti wave equation by the rapid expansion method (REM) A one-step time evolution algorithm: 79th Annual International Meeting, SEG, Expanded Abstrats, Tal-Ezer, H., D. Kosloff, and Z. Koren, 1987, An aurate sheme for seismi forward modeling: Geophysial Prospeting, 35, no. 5, , doi: / tb00830.x. Thorbeke, J. W., C. P. A. Wapenaar, and G. Swinnen, 004, Design of one-way wavefield extrapolation operators, using smooth funtions in WLSQ optimization: Geophysis, 69, , doi: / Wards, B., G. Margrave, and M. Lamoureux, 008, Phase-shift time-stepping for reverse-time migration: 78th Annual International Meeting, SEG, Expanded Abstrats, Zhang, Y., G. Zhang, D. Yingst, and J. Sun, 007, Expliit marhing method for reverse-time migration: 77th Annual International Meeting, SEG, Expanded Abstrats, Zhang, Y., and G. Zhang, 009, One-step extrapolation method for reverse-time migration: Geophysis, 74, no. 4, A9 A33, doi: /
Fourier finite-difference wave propagation a
Fourier finite-difference wave propagation a a Published in Geophysics, 76, T123-T129, (2011) Xiaolei Song and Sergey Fomel ABSTRACT We introduce a novel technique for seismic wave extrapolation in time.
More informationSUMMARY INTRODUCTION THEORY
Least-squares reverse-time migration using one-step two-way wave extrapolation by non-stationary phase shift Junzhe Sun, Sergey Fomel, and Jingwei Hu, The University of Texas at Austin Downloaded 11/06/14
More informationOne-step extrapolation method for reverse time migration
GEOPHYSICS VOL. 74 NO. 4 JULY-AUGUST 2009; P. A29 A33 5 FIGS. 10.1190/1.3123476 One-step extrapolation method for reverse time migration Yu Zhang 1 and Guanquan Zhang 2 ABSTRACT We have proposed a new
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More informationSeismic dip estimation based on the two-dimensional Hilbert transform and its application in random noise attenuation a
Seismi dip estimation based on the two-dimensional Hilbert transform and its appliation in random noise attenuation a a Published in Applied Geophysis, 1, 55-63 (Marh 015) Cai Liu, Changle Chen, Dian Wang,
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationA Spatiotemporal Approach to Passive Sound Source Localization
A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,
More informationMeasuring & Inducing Neural Activity Using Extracellular Fields I: Inverse systems approach
Measuring & Induing Neural Ativity Using Extraellular Fields I: Inverse systems approah Keith Dillon Department of Eletrial and Computer Engineering University of California San Diego 9500 Gilman Dr. La
More informationSURFACE WAVES OF NON-RAYLEIGH TYPE
SURFACE WAVES OF NON-RAYLEIGH TYPE by SERGEY V. KUZNETSOV Institute for Problems in Mehanis Prosp. Vernadskogo, 0, Mosow, 75 Russia e-mail: sv@kuznetsov.msk.ru Abstrat. Existene of surfae waves of non-rayleigh
More informationAdvances in Radio Science
Advanes in adio Siene 2003) 1: 99 104 Copernius GmbH 2003 Advanes in adio Siene A hybrid method ombining the FDTD and a time domain boundary-integral equation marhing-on-in-time algorithm A Beker and V
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationExtending LMR for anisotropic unconventional reservoirs
Extending LMR for anisotropi unonventional reservoirs Maro A. Perez Apahe Canada Ltd Summary It has beome inreasingly advantageous to haraterize rok in unonventional reservoirs within an anisotropi framework.
More informationThe Laws of Acceleration
The Laws of Aeleration The Relationships between Time, Veloity, and Rate of Aeleration Copyright 2001 Joseph A. Rybzyk Abstrat Presented is a theory in fundamental theoretial physis that establishes the
More informationAn Integrated Architecture of Adaptive Neural Network Control for Dynamic Systems
An Integrated Arhiteture of Adaptive Neural Network Control for Dynami Systems Robert L. Tokar 2 Brian D.MVey2 'Center for Nonlinear Studies, 2Applied Theoretial Physis Division Los Alamos National Laboratory,
More informationLowrank finite-differences and lowrank Fourier finite-differences for seismic wave extrapolation in the acoustic approximation a
Lowrank finite-differences and lowrank Fourier finite-differences for seismic wave extrapolation in the acoustic approximation a a Published in Geophysical Journal International, 193, 960-969, (2013) Xiaolei
More informationAcoustic Waves in a Duct
Aousti Waves in a Dut 1 One-Dimensional Waves The one-dimensional wave approximation is valid when the wavelength λ is muh larger than the diameter of the dut D, λ D. The aousti pressure disturbane p is
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW. P. М. Меdnis
ELECTROMAGNETIC WAVES WITH NONLINEAR DISPERSION LAW P. М. Меdnis Novosibirs State Pedagogial University, Chair of the General and Theoretial Physis, Russia, 636, Novosibirs,Viljujsy, 8 e-mail: pmednis@inbox.ru
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationIDENTIFICATION AND CONTROL OF ACOUSTIC RADIATION MODES
IDENTIFICATION AND CONTROL OF ACOUSTIC RADIATION MODES Arthur P. Berkhoff University of Twente, Faulty of Eletrial Engineering, P.O. Box 217, 7 AE Enshede, The Netherlands email: a.p.berkhoff@el.utwente.nl
More informationSeismic wave extrapolation using lowrank symbol approximation
Geophysical Prospecting doi: 10.1111/j.1365-2478.2012.01064.x Seismic wave extrapolation using lowrank symbol approximation Sergey Fomel 1, Lexing Ying 2 and Xiaolei Song 1 1 Bureau of Economic Geology,
More informationOn Certain Singular Integral Equations Arising in the Analysis of Wellbore Recharge in Anisotropic Formations
On Certain Singular Integral Equations Arising in the Analysis of Wellbore Reharge in Anisotropi Formations C. Atkinson a, E. Sarris b, E. Gravanis b, P. Papanastasiou a Department of Mathematis, Imperial
More informationResolving RIPS Measurement Ambiguity in Maximum Likelihood Estimation
14th International Conferene on Information Fusion Chiago, Illinois, USA, July 5-8, 011 Resolving RIPS Measurement Ambiguity in Maximum Likelihood Estimation Wenhao Li, Xuezhi Wang, and Bill Moran Shool
More informationCharacteristics of P, SV, and SH wave propagation in an anisotropic medium
Wave propagation in an anisotropi medium Charateristis of P, SV, and SH wave propagation in an anisotropi medium Amber C. Kelter and John C. Banroft ABSTRACT Two methods of approximating phase and group
More informationMODELING MATTER AT NANOSCALES. 4. Introduction to quantum treatments Eigenvectors and eigenvalues of a matrix
MODELING MATTER AT NANOSCALES 4 Introdution to quantum treatments 403 Eigenvetors and eigenvalues of a matrix Simultaneous equations in the variational method The problem of simultaneous equations in the
More information10.5 Unsupervised Bayesian Learning
The Bayes Classifier Maximum-likelihood methods: Li Yu Hongda Mao Joan Wang parameter vetor is a fixed but unknown value Bayes methods: parameter vetor is a random variable with known prior distribution
More informationProceedings of Meetings on Acoustics
Proeedings of Meetings on Aoustis Volume 19, 2013 http://aoustialsoiety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Arhitetural Aoustis Session 4aAAa: Room Aoustis Computer Simulation I 4aAAa8.
More informationConformal Mapping among Orthogonal, Symmetric, and Skew-Symmetric Matrices
AAS 03-190 Conformal Mapping among Orthogonal, Symmetri, and Skew-Symmetri Matries Daniele Mortari Department of Aerospae Engineering, Texas A&M University, College Station, TX 77843-3141 Abstrat This
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More informationAnalysis of discretization in the direct simulation Monte Carlo
PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,
More informationWave Propagation through Random Media
Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationGluing Potential Energy Surfaces with Rare Event Simulations
This is an open aess artile published under an ACS AuthorChoie Liense, whih permits opying and redistribution of the artile or any adaptations for non-ommerial purposes. pubs.as.org/jctc Gluing Potential
More informationElectromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.
arxiv:physis/99536v1 [physis.lass-ph] 15 May 1999 Eletromagneti radiation of the travelling spin wave propagating in an antiferromagneti plate. Exat solution. A.A.Zhmudsky November 19, 16 Abstrat The exat
More informationarxiv: v1 [physics.gen-ph] 5 Jan 2018
The Real Quaternion Relativity Viktor Ariel arxiv:1801.03393v1 [physis.gen-ph] 5 Jan 2018 In this work, we use real quaternions and the basi onept of the final speed of light in an attempt to enhane the
More information1 sin 2 r = 1 n 2 sin 2 i
Physis 505 Fall 005 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.5, 7.8, 7.16 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with
More informationA Time-Dependent Model For Predicting The Response Of A Horizontally Loaded Pile Embedded In A Layered Transversely Isotropic Saturated Soil
IOSR Journal of Mehanial and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 232-334X, Volume 16, Issue 2 Ser. I (Mar. - Apr. 219), PP 48-53 www.iosrjournals.org A Time-Dependent Model For Prediting
More informationDetermination of the reaction order
5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order
More informationFrequency Domain Analysis of Concrete Gravity Dam-Reservoir Systems by Wavenumber Approach
Frequeny Domain Analysis of Conrete Gravity Dam-Reservoir Systems by Wavenumber Approah V. Lotfi & A. Samii Department of Civil and Environmental Engineering, Amirkabir University of Tehnology, Tehran,
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationContact Block Reduction Method for Ballistic Quantum Transport with Semi-empirical sp3d5s* Tight Binding band models
Purdue University Purdue e-pubs Other Nanotehnology Publiations Birk Nanotehnology Center -2-28 Contat Redution Method for Ballisti Quantum Transport with Semi-empirial sp3d5s* Tight Binding band models
More informationWavetech, LLC. Ultrafast Pulses and GVD. John O Hara Created: Dec. 6, 2013
Ultrafast Pulses and GVD John O Hara Created: De. 6, 3 Introdution This doument overs the basi onepts of group veloity dispersion (GVD) and ultrafast pulse propagation in an optial fiber. Neessarily, it
More informationSEISMIC ANALYSIS OF SPHERICAL TANKS INCLUDING FLUID-STRUCTURE-SOIL INTERACTION
3 th World Conferene on Earthquake Engineering Vanouver, B.C., Canada August -6, 2004 aper o. 84 SEISMIC AALYSIS OF SHERICAL TAKS ICLUDIG FLUID-STRUCTURE-SOIL ITERACTIO T.L. Karavasilis, D.C. Rizos 2,
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationProcess engineers are often faced with the task of
Fluids and Solids Handling Eliminate Iteration from Flow Problems John D. Barry Middough, In. This artile introdues a novel approah to solving flow and pipe-sizing problems based on two new dimensionless
More informationEinstein s Three Mistakes in Special Relativity Revealed. Copyright Joseph A. Rybczyk
Einstein s Three Mistakes in Speial Relativity Revealed Copyright Joseph A. Rybzyk Abstrat When the evidene supported priniples of eletromagneti propagation are properly applied, the derived theory is
More informationTwo Points Hybrid Block Method for Solving First Order Fuzzy Differential Equations
Journal of Soft Computing and Appliations 2016 No.1 (2016) 43-53 Available online at www.ispas.om/jsa Volume 2016, Issue 1, Year 2016 Artile ID jsa-00083, 11 Pages doi:10.5899/2016/jsa-00083 Researh Artile
More informationMoments and Wavelets in Signal Estimation
Moments and Wavelets in Signal Estimation Edward J. Wegman 1 Center for Computational Statistis George Mason University Hung T. Le 2 International usiness Mahines Abstrat: The problem of generalized nonparametri
More informationStable wide-angle Fourier finite-difference downward extrapolation of 3-D. wavefields
Stable wide-angle Fourier finite-difference downward extrapolation of 3-D wavefields Biondo Biondi Stanford Exploration Project, Mitchell Bldg., Department of Geophysics, Stanford University, Stanford,
More informationSupplementary information for: All-optical signal processing using dynamic Brillouin gratings
Supplementary information for: All-optial signal proessing using dynami Brillouin gratings Maro Santagiustina, Sanghoon Chin 2, Niolay Primerov 2, Leonora Ursini, Lu Thévena 2 Department of Information
More informationA new initial search direction for nonlinear conjugate gradient method
International Journal of Mathematis Researh. ISSN 0976-5840 Volume 6, Number 2 (2014), pp. 183 190 International Researh Publiation House http://www.irphouse.om A new initial searh diretion for nonlinear
More informationA Queueing Model for Call Blending in Call Centers
A Queueing Model for Call Blending in Call Centers Sandjai Bhulai and Ger Koole Vrije Universiteit Amsterdam Faulty of Sienes De Boelelaan 1081a 1081 HV Amsterdam The Netherlands E-mail: {sbhulai, koole}@s.vu.nl
More informationCombined Electric and Magnetic Dipoles for Mesoband Radiation, Part 2
Sensor and Simulation Notes Note 53 3 May 8 Combined Eletri and Magneti Dipoles for Mesoband Radiation, Part Carl E. Baum University of New Mexio Department of Eletrial and Computer Engineering Albuquerque
More informationSingular Event Detection
Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate
More informationLECTURE NOTES FOR , FALL 2004
LECTURE NOTES FOR 18.155, FALL 2004 83 12. Cone support and wavefront set In disussing the singular support of a tempered distibution above, notie that singsupp(u) = only implies that u C (R n ), not as
More informationv v Downloaded 01/11/16 to Redistribution subject to SEG license or copyright; see Terms of Use at
The pseudo-analytical method: application of pseudo-laplacians to acoustic and acoustic anisotropic wae propagation John T. Etgen* and Serre Brandsberg-Dahl Summary We generalize the pseudo-spectral method
More informationFeature Selection by Independent Component Analysis and Mutual Information Maximization in EEG Signal Classification
Feature Seletion by Independent Component Analysis and Mutual Information Maximization in EEG Signal Classifiation Tian Lan, Deniz Erdogmus, Andre Adami, Mihael Pavel BME Department, Oregon Health & Siene
More informationSimplified Buckling Analysis of Skeletal Structures
Simplified Bukling Analysis of Skeletal Strutures B.A. Izzuddin 1 ABSRAC A simplified approah is proposed for bukling analysis of skeletal strutures, whih employs a rotational spring analogy for the formulation
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationLong time stability of regularized PML wave equations
Long time stability of regularized PML wave equations Dojin Kim Email:kimdojin@knu.a.kr Yonghyeon Jeon Email:dydgus@knu.a.kr Philsu Kim Email:kimps@knu.a.kr Abstrat In this paper, we onsider two dimensional
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 4, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume, No 4, 01 Copyright 010 All rights reserved Integrated Publishing servies Researh artile ISSN 0976 4399 Strutural Modelling of Stability
More informationNonreversibility of Multiple Unicast Networks
Nonreversibility of Multiple Uniast Networks Randall Dougherty and Kenneth Zeger September 27, 2005 Abstrat We prove that for any finite direted ayli network, there exists a orresponding multiple uniast
More informationArray Design for Superresolution Direction-Finding Algorithms
Array Design for Superresolution Diretion-Finding Algorithms Naushad Hussein Dowlut BEng, ACGI, AMIEE Athanassios Manikas PhD, DIC, AMIEE, MIEEE Department of Eletrial Eletroni Engineering Imperial College
More informationTime Domain Method of Moments
Time Domain Method of Moments Massahusetts Institute of Tehnology 6.635 leture notes 1 Introdution The Method of Moments (MoM) introdued in the previous leture is widely used for solving integral equations
More informationEXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION
Journal of Mathematial Sienes: Advanes and Appliations Volume 3, 05, Pages -3 EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION JIAN YANG, XIAOJUAN LU and SHENGQIANG TANG
More informationLowrank RTM for converted wave imaging
Lowrank RM for converted wave imaging Lorenzo Casasanta * (CGG), Zhiguang Xue (U Austin), Sam Gray (CGG) Summary his paper is an attempt to fill the technology gap existing between pure P- and PS-wave
More informationA FETI Method for a Class of Indefinite or Complex Second- or Fourth-Order Problems
A FETI Method for a Class of Indefinite or Complex Seond- or Fourth-Order Problems Charbel Farhat 1, Jing Li 2, Mihel Lesoinne 1 and Philippe Avery 1 1 University of Colorado at Boulder, Department of
More informationThree-dimensional morphological modelling in Delft3D-FLOW
Three-dimensional morphologial modelling in Delft3D-FLOW G. R. Lesser, J. van Kester, D.J.R. Walstra and J.A. Roelvink WL delft hydraulis email: giles.lesser@ldelft.nl Abstrat Computer modelling of sediment
More informationAssessing the Performance of a BCI: A Task-Oriented Approach
Assessing the Performane of a BCI: A Task-Oriented Approah B. Dal Seno, L. Mainardi 2, M. Matteui Department of Eletronis and Information, IIT-Unit, Politenio di Milano, Italy 2 Department of Bioengineering,
More informationTHE REFRACTION OF LIGHT IN STATIONARY AND MOVING REFRACTIVE MEDIA
HDRONIC JOURNL 24, 113-129 (2001) THE REFRCTION OF LIGHT IN STTIONRY ND MOVING REFRCTIVE MEDI C. K. Thornhill 39 Crofton Road Orpington, Kent, BR6 8E United Kingdom Reeived Deember 10, 2000 Revised: Marh
More informationLaplace-Transform Finite Element Solution of Nonlocal and Localized Stochastic Moment Equations of Transport
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 6, No., pp. 3-6 Commun. Comput. Phys. July 9 Laplae-Transform Finite Element Solution of Nonloal and Loalized Stohasti Moment Equations of Transport Eri Morales-Casique
More informationSensor Network Localisation with Wrapped Phase Measurements
Sensor Network Loalisation with Wrapped Phase Measurements Wenhao Li #1, Xuezhi Wang 2, Bill Moran 2 # Shool of Automation, Northwestern Polytehnial University, Xian, P.R.China. 1. wenhao23@mail.nwpu.edu.n
More informationQCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines
DOI.56/sensoren6/P3. QLAS Sensor for Purity Monitoring in Medial Gas Supply Lines Henrik Zimmermann, Mathias Wiese, Alessandro Ragnoni neoplas ontrol GmbH, Walther-Rathenau-Str. 49a, 7489 Greifswald, Germany
More informationTransformation to approximate independence for locally stationary Gaussian processes
ransformation to approximate independene for loally stationary Gaussian proesses Joseph Guinness, Mihael L. Stein We provide new approximations for the likelihood of a time series under the loally stationary
More informationLecture 7: Sampling/Projections for Least-squares Approximation, Cont. 7 Sampling/Projections for Least-squares Approximation, Cont.
Stat60/CS94: Randomized Algorithms for Matries and Data Leture 7-09/5/013 Leture 7: Sampling/Projetions for Least-squares Approximation, Cont. Leturer: Mihael Mahoney Sribe: Mihael Mahoney Warning: these
More informationBreakdown of the Slowly Varying Amplitude Approximation: Generation of Backward Traveling Second Harmonic Light
Claremont Colleges Sholarship @ Claremont All HMC Faulty Publiations and Researh HMC Faulty Sholarship 1-1-003 Breakdown of the Slowly Varying Amplitude Approximation: Generation of Bakward Traveling Seond
More informationA Differential Equation for Specific Catchment Area
Proeedings of Geomorphometry 2009. Zurih, Sitzerland, 3 ugust - 2 September, 2009 Differential Equation for Speifi Cathment rea J. C. Gallant, M. F. Huthinson 2 CSIRO Land and Water, GPO Box 666, Canberra
More informationA Combination of BIM and BEM for efficiently analyzing optical. elements
A Combination of BIM and BEM for effiiently analyzing optial elements SUN Fang, LIU Juan, ZHANG Guo-ting, HU Chuan-fei, SU Xiaoxing Institute of Lightwave Tehnology, hool of eletronis and information engineering,
More informationPerturbation Analyses for the Cholesky Factorization with Backward Rounding Errors
Perturbation Analyses for the holesky Fatorization with Bakward Rounding Errors Xiao-Wen hang Shool of omputer Siene, MGill University, Montreal, Quebe, anada, H3A A7 Abstrat. This paper gives perturbation
More informationAn Adaptive Optimization Approach to Active Cancellation of Repeated Transient Vibration Disturbances
An aptive Optimization Approah to Ative Canellation of Repeated Transient Vibration Disturbanes David L. Bowen RH Lyon Corp / Aenteh, 33 Moulton St., Cambridge, MA 138, U.S.A., owen@lyonorp.om J. Gregory
More informationUNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS. V. N. Matveev and O. V. Matvejev
UNCERTAINTY RELATIONS AS A CONSEQUENCE OF THE LORENTZ TRANSFORMATIONS V. N. Matveev and O. V. Matvejev Joint-Stok Company Sinerta Savanoriu pr., 159, Vilnius, LT-315, Lithuania E-mail: matwad@mail.ru Abstrat
More informationSensor management for PRF selection in the track-before-detect context
Sensor management for PRF seletion in the tra-before-detet ontext Fotios Katsilieris, Yvo Boers, and Hans Driessen Thales Nederland B.V. Haasbergerstraat 49, 7554 PA Hengelo, the Netherlands Email: {Fotios.Katsilieris,
More informationRobust Recovery of Signals From a Structured Union of Subspaces
Robust Reovery of Signals From a Strutured Union of Subspaes 1 Yonina C. Eldar, Senior Member, IEEE and Moshe Mishali, Student Member, IEEE arxiv:87.4581v2 [nlin.cg] 3 Mar 29 Abstrat Traditional sampling
More informationGeneration of EM waves
Generation of EM waves Susan Lea Spring 015 1 The Green s funtion In Lorentz gauge, we obtained the wave equation: A 4π J 1 The orresponding Green s funtion for the problem satisfies the simpler differential
More informationPlanning with Uncertainty in Position: an Optimal Planner
Planning with Unertainty in Position: an Optimal Planner Juan Pablo Gonzalez Anthony (Tony) Stentz CMU-RI -TR-04-63 The Robotis Institute Carnegie Mellon University Pittsburgh, Pennsylvania 15213 Otober
More informationCopyright 2018 Society of Photo-Optical Instrumentation Engineers (SPIE). One print or electronic copy may be made for personal use only.
Copyright 018 Soiety of Photo-Optial Instrumentation Engineers (SPIE) One print or eletroni opy may be made for personal use only Systemati reprodution and distribution, dupliation of any material in this
More informationSYNTHETIC APERTURE IMAGING OF DIRECTION AND FREQUENCY DEPENDENT REFLECTIVITIES
SYNTHETIC APERTURE IMAGING OF DIRECTION AND FREQUENCY DEPENDENT REFLECTIVITIES LILIANA BORCEA, MIGUEL MOSCOSO, GEORGE PAPANICOLAOU, AND CHRYSOULA TSOGKA Abstrat. We introdue a syntheti aperture imaging
More informationOptimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach
Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)
More informationMOLECULAR ORBITAL THEORY- PART I
5.6 Physial Chemistry Leture #24-25 MOLECULAR ORBITAL THEORY- PART I At this point, we have nearly ompleted our rash-ourse introdution to quantum mehanis and we re finally ready to deal with moleules.
More informationCOMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and
More informationBottom Shear Stress Formulations to Compute Sediment Fluxes in Accelerated Skewed Waves
Journal of Coastal Researh SI 5 453-457 ICS2009 (Proeedings) Portugal ISSN 0749-0258 Bottom Shear Stress Formulations to Compute Sediment Fluxes in Aelerated Skewed Waves T. Abreu, F. Sanho and P. Silva
More informationThe experimental plan of displacement- and frequency-noise free laser interferometer
7th Edoardo Amaldi Conferene on Gravitational Waves (Amaldi7) Journal of Physis: Conferene Series 122 (2008) 012022 The experimental plan of displaement- and frequeny-noise free laser interferometer K
More informationUTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.
UTC Engineering 329 Proportional Controller Design for Speed System By John Beverly Green Team John Beverly Keith Skiles John Barker 24 Mar 2006 Introdution This experiment is intended test the variable
More informationLagrangian Formulation of the Combined-Field Form of the Maxwell Equations
Physis Notes Note 9 Marh 009 Lagrangian Formulation of the Combined-Field Form of the Maxwell Equations Carl E. Baum University of New Mexio Department of Eletrial and Computer Engineering Albuquerque
More informationGeneralised Differential Quadrature Method in the Study of Free Vibration Analysis of Monoclinic Rectangular Plates
Amerian Journal of Computational and Applied Mathematis 0, (4): 66-7 DOI: 0.59/.aam.0004.05 Generalised Differential Quadrature Method in the Study of Free Vibration Analysis of Monolini Retangular Plates
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationHILLE-KNESER TYPE CRITERIA FOR SECOND-ORDER DYNAMIC EQUATIONS ON TIME SCALES
HILLE-KNESER TYPE CRITERIA FOR SECOND-ORDER DYNAMIC EQUATIONS ON TIME SCALES L ERBE, A PETERSON AND S H SAKER Abstrat In this paper, we onsider the pair of seond-order dynami equations rt)x ) ) + pt)x
More information