generalized bremmer series with rational approximation for the scattering of waves in inhomogeneous media the bremmer series solution of the wave equation in generally inhomogeneous media requires the introduction of pseudodifferential operators in this paper sparse matrix representations of these pseudodifferential operators are derived the authors focus on designing sparse matrices keeping the accuracy high at the cost of ignoring any critical scattering angle phenomena such matrix representations follow from rational approximations of the vertical slowness and the transverse laplace operator symbols and of the vertical derivative as they appear in the parabolic equation method sparse matrix representations lead to a fast algorithm an optimization procedure is followed to minimize the errors in the high frequency limit for a given discretization rate the bremmer series solver consists of three steps directional decomposition into up and downgoing waves one way propagation and interaction of the counterpropagating constituents each of these steps is represented by a sparse matrix equation the resulting algorithm provides an improvement of the parabolic equation method in particular for transient wave phenomena and extends the latter method systematically for backscattered waves

interaction cohomology of forward or backward self similar systems abstract we investigate the dynamics of forward or backward self similar systems iterated function systems and the topological structure of their invariant sets we define a new cohomology theory interaction cohomology for forward or backward self similar systems we show that under certain conditions the space of connected components of the invariant set is isomorphic to the inverse limit of the spaces of connected components of the realizations of the nerves of finite coverings u of the invariant set where each u consists of backward images of the invariant set under elements of finite word length we give a criterion for the invariant set to be connected moreover we give a sufficient condition for the first cohomology group to have infinite rank as an application we obtain many results on the dynamics of semigroups of polynomials moreover we define postunbranched systems and we investigate the interaction cohomology groups of such systems many examples are given

flux normalized wavefield decomposition and migration of seismic data separation of wavefields into directional components can be accomplished by an eigenvalue decomposition of the accompanying system matrix in conventional pressurenormalized wavefield decomposition the resulting one way wave equations contain an interaction term which depends on the reflectivity function applying directional wavefield decomposition using flux normalized eigenvalue decomposition and disregarding interaction between up and downgoing wavefields these interaction terms were absent by also applying a correction term for transmission loss the result was an improved estimate of the up and downgoing wavefields in the wave equation angle transform a crosscorrelation function in local offset coordinates was fourier transformed to produce an estimate of reflectivity as a function of slowness or angle we normalized this wave equation angle transform with an estimate of the planewave reflection coefficient the flux normalized oneway wave propagation scheme was applied to imaging and to the normalized wave equation angle transform on synthetic and field data this proved the effectiveness of the new methods

application of unsplit convolutional perfectly matched layer for scalar arbitrarily wide angle wave equation abstracta classical split perfectly matched layer pml method has recently been applied to the scalar arbitrarily wide angle wave equation awwe in terms of displacement however the classical split pml obviously increases computational cost and cannot efficiently absorb waves propagating into the absorbing layer at grazing incidence our goal was to improve the computational efficiency of awwe and to enhance the suppression of edge reflections by applying a convolutional pml cpml we reformulated the original awwe as a first order formulation and incorporated the cpml with a general complex frequency shifted stretching operator into the renewed formulation a staggered grid finite difference fd method was adopted to discretize the first order equation system for wavefield depth continuation the first order awwe with the cpml saved memory compared with the original second order awwe with the conventional split pml with the help of numerical examples we verified the correctness of the staggered

target oriented wave equation least squares migration inversion with phase encoded hessian prestack depth migration produces blurred images resulting from limited acquisition apertures complexities in the velocity model and band limited characteristics of seismic waves this distortion can be partially corrected using the model space least squares migration inversion approach where a target oriented wave equation hessian operator is computed explicitly and then inverse filtering is applied iteratively to deblur or invert for the reflectivity however one difficulty is the cost of computing the explicit hessian operator which requires storing a large number of green s functions making it challenging for large scale applications a new method to compute the hessian operator for the wave equation based least squares migration inversion problem modifies the original explicit hessian formula enabling efficient computation of this operator an advantage is that the method eliminates disk storage of green s functions the modifications however also introduce undesired crosstalk artifacts two different phase encoding schemes plane wave phase encoding and random phase encoding suppress the crosstalk when the randomly phase encoded hessian operator is applied to the sigsbee2a synthetic data set an improved subsalt image with more balanced amplitudes is obtained

reducing two way splitting error of ffd method in dual domains the fourier finite difference ffd method is very popular in seismic depth migration but its straightforward 3d extension creates two way splitting error due to ignoring the cross terms of spatial partial derivatives traditional correction schemes either in the spatial domain by the implicit finite difference method or in the wavenumber domain by phase compensation lead to substantially increased computational costs or numerical difficulties for strong velocity contrasts we propose compensating the two way splitting error in dual domains alternately in the spatial and wavenumber domains via fourier transform first we organize the expanded square root operator in terms of two way splitting ffd plus the usually ignored cross terms second we select a group of optimized coefficients to maximize the accuracy of propagation in both inline and crossline directions without yet considering the diagonal directions finally we further optimize the constant coefficient of the compensation part to further improve the overall accuracy of the operator in implementation the compensation terms are similar to the high order corrections of the generalized screen method but their functions are to compensate the two way splitting error rather than the expansion error numerical experiments show that optimized one term compensation can achieve nearly perfect circular impulse responses and the propagation angle with less than 1 error for all azimuths is improved up to 60 from 35 compared with traditional single domain methods our scheme can handle lateral velocity variations even for strong velocity contrasts much more easily with only one additional fourier transform based on the two way splitting ffd method which helps retain the computational efficiency

exact and operator rational approximate solutions of the helmholtz weyl composition equation in underwater acoustics the quadratic profile the application of phase space and path integral methods to both mathematical and computational direct and inverse wave propagation modeling at the level of the scalar one way helmholtz equation is briefly reviewed the role of operator symbols is stressed and their properties briefly discussed the construction of the operator symbol requires the exact or approximate solution of the helmholtz weyl composition equation in the weyl pseudodifferential operator calculus the exact symbols for several quadratic profiles are presented and briefly analyzed the symbols exactly corresponding to a family of operator rational approximations are also presented for the quadratic case these results are used to illustrate several points pertinent to wide angle propagation modeling and the refractive index profile reconstruction problem in underwater acoustics

discovering fuzzy exception and anomalous rules nowadays searching for specific kind of knowledge that deviates from the usual standards is very useful in several domains such as network traffic anomalies fraud detection economic analysis or medical diagnosis fuzzy association rules have been developed as a powerful tool for dealing with imprecision in databases that may come from the source i e imprecise measures taken by the machine or from the human understanding of a concept and offering a comprehensive representation of found knowledge in this paper we introduce the notion of fuzzy exception and fuzzy anomalous rule for the recognition of these types of deviations the deviations are associated with the common patterns which usually are hidden in data affected by some kind of fuzziness we present a new approach for mining such rules based on a recently proposed model for representing and evaluating fuzzy rules important advantages are to obtain more understandable results and that the mining process can be parallelized an algorithm following the proposed model is developed and some experiments are performed in data where some numerical attributes have been fuzzified and also in some real fuzzy transactional datasets for testing the proposed algorithm from experimentation we have seen that the proposed fuzzy rules give some insights on the exception and anomaly detection in credit payments

examples counterexamples and measuring fuzzy associations this paper examines the measurement of the degree to which tuples in a database support a relation among attributes based on a comparison of the number of examples and counterexamples of the relation in particular we are concerned with associations that represent imprecise constraints placed upon the value of one attribute by those of other attributes associations of this form may be described by fuzzy rules and the analysis requires an assessment of the degree to which a tuple satisfies the imprecise constraint specified by the rule standard measures of rule validity are extended to fuzzy associations based on the degree that the tuples are examples counterexamples or irrelevant to imprecise relations a scaling of the relevance of a tuple is proposed to minimize the impact of the accumulation of small membership values on the confidence based validity measures

derivation and application of extended parabolic wave theories i the factorized helmholtz equation the reduced scalar helmholtz equation for a transversely inhomogeneous half space supplemented with an outgoing radiation condition and an appropriate boundary condition on the initial value plane defines a direct acoustic propagation model this elliptic formulation admits a factorization and is subsequently equivalent to a first order weyl pseudodifferential equation which is recognized as an extended parabolic propagation model perturbation treatments of the appropriate weyl composition equation result in a systematic development of approximate wave theories which extend the narrow angle weak inhomogeneity and weak gradient ordinary parabolic schrodinger approximation the analysis further provides for the formulation and exact solution of a multidimensional nonlinear inverse problem appropriate for ocean acoustic and seismic studies the wave theories foreshadow computational algorithms the inclusion of range dependent effects and the extension to 1 the vector formulation appropriate for elas

riemann map and holomorphic dynamics f f r a f r a and f u c c l a f r a i e fra repels to the side of a n o assume that u is always sufficiently small so that f has no singularities in u c a let r d z a be a riemann map a conformal homeomorphism of the unit disc onto a then there exists a holomorphic extension g of r l o f o r to a neighbourhood of s 1 and gls happens to be an expanding map i e there exists n 0 such that for every z s l l g z l although the proofs are straightforward we include them for sceptics in the last section sect 7 together with a remark about examples of a and f satisfying the above assumptions denote by gthe non tangential limit of r and by k s 1 the domain where it exists denote by ijl g the space of all borel probability g invariant ergodic positive entropy measures on s 1 it is known see p that for every egji g k exists almost everywhere so the f invariant measure r on fra can be considered the functions log g log lf are z respectively r 0 integrable see 1 5 denote these integrals by u g g f in this paper we consider derivatives with respect to the riemannian metric on 2 define the hausdorff dimension of any probability measure v

fuzzy association rules general model and applications the theory of fuzzy sets has been recognized as a suitable tool to model several kinds of patterns that can hold in data in this paper we are concerned with the development of a general model to discover association rules among items in a crisp set of fuzzy transactions this general model can be particularized in several ways each particular instance corresponds to a certain kind of pattern and or repository of data we describe some applications of this scheme paying special attention to the discovery of fuzzy association rules in relational databases

wide angle fd and ffd migration using complex pade approximations seismic migration by downward continuation using the one way wave equation approximations has two shortcomings imaging steep dip reflectors and handling evanescent waves complex pade approximations allow a better treatment of evanescent modes stabilizing finite difference migration without requiring special treatment for the migration domain boundaries imaging of steep dip reflectors can be improved using several terms in the pade expansion we discuss the implementation and evaluation of wide angle complex pade approximations for finite difference and fourier finite difference migration methods the dispersion relation and the impulsive response of the migration operator provide criteria to select the number of terms and coefficients in the pade expansion this ensures stability for a prescribed maximum propagation direction the implementations are validated on the marmousi model data set and seg eage salt model data

friwe fuzzy rule identification with exceptions in this paper the friwe method is proposed to identify fuzzy models from examples such a method has been developed trying to achieve a double goal accuracy and interpretability in order to do that maximal structure fuzzy rules are firstly obtained based on a method proposed by castro et al in a second stage the conflicts generated by the maximal rules are solved thus increasing the model accuracy the resolution of conflicts are carried out by including exceptions in the rules this strategy has been identified by psychologists with the learning mechanism employed by the human being thus improving the model interpretability besides in order to improve the interpretability even more several methods are presented based on reducing and merging rules and exceptions in the model the exhaustive use of the training examples gives the method a special suitability for problems with small training sets or high dimensionality finally the method is applied to an example in order to analyze the achievement of the goals

target oriented beamlet migration based on gabor daubechies frame decomposition we develop beamlet propagation and imaging using gabor daubechies g d frame decomposition based on local perturbation theory and apply it to target oriented prestack depth migration the method is formulated with local background velocities and local perturbations in wavefield extrapolation the localized propagators and phase correction operators are obtained analytically or semianalytically by one way operator decomposition and screen approximation in the beamlet and space beamlet mixed domain beamlet wavefields have superior localization properties in both local space and direction wavenumber over gaussian beams in the sense that localizations are not limited within short propagation distances in either homogeneous or heterogeneous media comparisons of the prestack depth migrated images for the 2d seg eage salt model and the marmousi model indicate that for seismic wave propagation and imaging in complex structural environments the g d beamlet propagator has higher accuracy and better wide angle properties than do global propagators target oriented prestack g d beamlet migration is performed by means of local angledomain imaging and controlled superposition of commondip angle images based on the local directivity features of the target structures this considers the spatial and direction localizations of beamlets as a numerical example we process the 2d seg eage salt model prestack data the results show that the proposed migration method has considerable advantages in suppressing noise and enhancing structural features image quality for subsalt structures especially for steep faults is improved through structure based superposition of common angle images this demonstrates the potential and capability of beamlet migration in target oriented seismic imaging

globally optimized finite difference extrapolator for strongly vti media implicit finite difference fd migration is unconditionally stable and is popular in handling strong velocity variations but its extension to strongly transversely anisotropic media with vertical symmetric axis media is difficult traditional local optimizations generate the optimized coefficients for each pair of thomsen anisotropy parameters independently which can degrade results substantially for large anisotropy variations and lead to a huge table we developed an implicit fd method using the analytic taylor series expansion and used a global optimization scheme to improve its accuracy for wide phase angles we first extended the number of the constant coefficients then we relaxed the coefficient of the time delay extrapolation term by tuning a small factor such that error is less than 0 1 finally we optimized the constant coefficients using a simulated annealing algorithm by constraining that all the error functions on a fine grid of the whole anisotropic region did not exceed 0 5 simultaneously the extended number of the constant coefficients and the relaxed coefficient greatly enhanced the flexibility of matching the dispersion relation and significantly improved the ability of handling strong anisotropy over a much wider range compared with traditional local optimization our scheme does not need any table and table lookup for each order of the fd method only one group of optimized coefficients is enough to handle strong variations in velocity and anisotropy more importantly our global optimization scheme guarantees the accuracy for various possible ranges of anisotropy parameters no matter how strong the anisotropy is for the globally optimized second order fd method the accurate phase angle is up to 58 and the increase is about 18 22 for the globally optimized fourth order fd method the accurate phase angle is up to 77 and the increase is about 22 27

fourier finite difference migration many existing migration schemes cannot simultaneously handle the two most important problems of migration imaging of steep dips and imaging in media with arbitrary velocity variations in all directions for example phase shift k migration is accurate for nearly all dips but it is limited to very simple velocity functions on the other hand finite difference schemes based on one way wave equations consider arbitrary velocity functions but they attenuate steeply dipping events we propose a new hybrid migration method named fourier finite difference migration wherein the downward continuation operator is split into two downward continuation operators one operator is a phase shift operator for a chosen constant background velocity and the other operator is an optimized finite difference operator for the varying component of the velocity function if there is no variation of velocity then only a phase shift operator will be applied automatically on the other hand if there is a strong variatio

self affine tiles in ℝn abstract a self affine tile in r n is a set t of positive measure with a t d d t d where a is an expanding n n real matrix with det a m an integer and d d d 2 d m r n is a set of m digits it is known that self affine tiles always give tilings of r n by translation this paper extends known characterizations of digit sets d yielding self affine tiles it proves several results about the structure of tilings of r n possible using such tiles and gives examples showing the possible relations between self replicating tilings and general tilings which clarify results of kenyon on self replicating tilings

rational square root approximations for parabolic equation algorithms in this article stable pade approximations to the function 1 z are derived by choosing a branch cut in the negative half plane the pade coefficients are complex and may be derived analytically to arbitrary order from the corresponding real coefficients associated with the principal branch defined by z 1 i z 0 i z denotes the imaginary part of z the characteristics of the corresponding square root approximation are illustrated for various segments of the complex plane in particular for waveguide problems it is shown that an increasingly accurate representation may be obtained of both the evanescent part of the mode spectrum for the acoustic case and the complex mode spectrum for the elastic case an elastic parabolic equation algorithm is used to illustrate the application of the new pade approximations to a realistic ocean environment including elasticity in the ocean bottom

wave equation based seismic illumination analysis we present a wave equation based method for seismic illumination analysis a one way wave equation based generalized screen propagator is used to extrapolate the wavefields from sources and receivers to the subsurface target a local plane wave analysis is used at the target to calculate localized directionalenergyfluxesforbothsourceandreceiver wavefields we construct an illumination matrix using these energy fluxes to quantify the target illumination conditions the target geometry information is used to manipulate the illumination matrix and generate different types of illumination measures the wave equation based approach can properly handle forward multiple scattering phenomena including focusing defocusing diffraction and interference effects itcanbedirectlyappliedtocomplexvelocitymodels velocity model smoothing and fresnel zone smoothing are not required different illumination measurements derived from this method can be applied to target oriented or volumetric illumination analyses this new method is flexible and practical for illumination analysis in complex 2d and 3d velocity models with nontrivial acquisition and target

elicitation of fuzzy association rules from positive and negative examples the aim of this paper is to provide a crystal clear insight into the true semantics of the measures of support and confidence that are used to assess rule quality in fuzzy association rule mining to achieve this we rely on two important pillars the identification of transactions in a database as positive or negative examples of a given association between attributes and the correspondence between measures of support and confidence on one hand and measures of compatibility and inclusion on the other hand in this way we remove the mystery from recently suggested quality measures for fuzzy association rules

an acoustic wave equation for anisotropic media a wave equation derived using the acoustic medium assumption for p waves in transversely isotropic ti media with a vertical symmetry axis vti media yields a good kinematic approximation to the familiar elastic wave equation for vti media the wavefield solutions obtained using this vti acoustic wave equation are free of shear waves which significantly reduces the computation time compared to the elastic wavefield solutions for exploding reflector type applications from this vti acoustic wave equation the eikonal and transport equations that describe the ray theoretical aspects of wave propagation in a ti medium are derived these equations based on the acoustic assumption shear wave velocity 0 are much simpler than their elastic counterparts yet they yield an accurate description of traveltimes and geometrical amplitudes numerical examples prove the usefulness of this acoustic equation in simulating the kinematic aspects of wave propagation in complex ti models

a formal model for mining fuzzy rules using the rl representation theory data mining techniques managing imprecision are very useful to obtain meaningful and interesting information for the user among some other techniques fuzzy association rules have been developed as a powerful tool for dealing with imprecision in databases and offering a good representation of found knowledge in this paper we introduce a formal model for managing the imprecision in fuzzy transactional databases using the restriction level representation theory a recent representation of imprecision that extends that of fuzzy sets this theory introduces some new operators keeping the usual crisp properties even when negation is involved the model allows us to mine fuzzy association rules in a straightforward way extending the accuracy measures from the crisp case in addition we introduce several ways of representing and summarizing the obtained results in order to offer new and very interesting semantics as an application we present how to extract fuzzy association rules involving both the presence and the absence of items using the proposed model and we also perform some experiments with real fuzzy transactional datasets

coding and tiling of julia sets for subhyperbolic rational maps abstract let f c c be a subhyperbolic rational map of degree d we construct a set of proper coding maps cod f r j r of the julia set j by geometric coding trees where the parameter r ranges over mappings from a certain tree to the riemann sphere using the universal covering space s s for the corresponding orbifold we lift the inverse of f to an iterated function system i g i i 1 2 d for the purpose of studying the structure of cod f we generalize kenyon and lagarias wang s results if the attractor k of i has positive measure then k tiles 1 j and the multiplicity of r is well defined moreover we see that the equivalence relation induced by r is described by a finite directed graph and give a necessary and sufficient condition for two coding maps r and r to be equal

factorization and path integration of the helmholtz equation numerical algorithms the propagator for the reduced scalar helmholtz equation plays a significant role in both analytical and computational studies of acoustic direct wave propagation path functional integrals are taken to provide the principal representation of the propagator and are computed directly the path integral is the primary tool in extending the classical fourier methods so appropriate for wave propagation in homogeneous media to inhomogeneous media for transversely inhomogeneous environments the n dimensional helmholtz equation can be exactly factored into separate forward and backward one way wave equations a parabolic based one way phase space path integral construction provides the generalization of the tappert hardin split step fft algorithm to the full one way factored helmholtz wave equation these extended marching algorithms can readily accommodate density profiles and range updating and further in conjunction with imbedding methods provide the basis for incorporating backscatter effects in a complementary manner for general range dependent environments elliptic based two way path integral constructions lead to an approximate representation of the propagator feynman garrod and a natural statistical monte carlo means of evaluation taken together the path integrals provide the basis for a global analysis in addition to providing a unifying framework for dynamical approximations resolution of the square root operator and the concept of an underlying stochastic process the one way marching algorithms are applied to ocean acoustic environments seismological environments and extreme model environments designed to establish their range of validity and manner of breakdown

arbitrarily wide angle wave equations for complex media by combining various ideas related to one way wave equations owwes half space stiffness relation special finite element discretization and complex coordinate stretching a systematic procedure is developed for deriving a series of highly accurate space domain versions of owwes the resulting procedure is applicable to complex media where the governing equation full wave equation is a second order differential system making the procedure applicable for general heterogeneous anisotropic porous viscoelastic media owing to their high accuracy in representing waves propagating in an arbitrarily wide range of angles the resulting equations are named arbitrarily wide angle wave equations awwes in order to illustrate the proposed procedure awwes are derived for one way propagation in acoustic as well as elastic media while acoustic awwes can be considered as modified versions of well known space domain owwes based on rational approximations of the square root operator the elastic awwes are significantly different from the existing elastic owwes unlike the existing elastic owwes elastic awwes are displacement based and are applicable to general anisotropic media furthermore awwes are simple in their form and appear amenable to easy numerical implementation

a stable one way wave propagator for vti media for the purpose of one way wave equation imaging a pseudoscreen propagator is developed for transversely isotropic media with vertical axes of symmetry the phase shift for propagation through a depth slice is decomposed into three terms a gazdag phase shift for propagation in a laterally homogeneous reference medium a correction for lateral variability of vertical propagation and a remaining wide angle term for oblique directions of propagation based on rational function approximation for this remaining wide angle term a fourier finite difference ffd approach with four way splitting is applied fourth order pade approximation is unsatisfactory in anelliptic media for large propagation angles with respect to the vertical direction therefore a method of coefficient optimization is developed in conjunction with a method of choosing an adequate homogeneous reference medium in a depth slice by symmetrizing the finite difference operators and because of the choice of the optimized coefficients the propagator is stable in the sense that the least squares norm of the wavefield measured for a frequency depth slice does not grow with increasing depth of propagation a small amount of artificial damping is applied to suppress artifacts that appear at the critical angle defined by the velocities in the reference medium and the actual medium synthetic examples confirm that good kinematic accuracy is achieved for a wide range of propagation angles typically up to 60

scalar generalized screen algorithms in transversely isotropic media with a vertical symmetry axis the scalar generalized screen method in isotropic media is extended here to transversely isotropic media with a vertical symmetry axis vti although wave propagation in a transversely isotropic medium is essentially elastic we use an equivalent acoustic system of equations for the qp waves which we prove to be accurate for both the dispersion relation and the polarization angle in the case of mild anisotropy the enhanced accuracy of the generalized screen method as compared to the split step fourier methods allows the extension to vti media the generalized screen expansion of the one way propagator follows closely the method used in the isotropic case the medium is defined in terms of a background and a perturbation the generalized screen expansion of the vertical slowness is based upon an expansion of the medium parameters simultaneously into magnitude and smoothness of variation we cast the theory into numerical algorithms and assess the accuracy of the generalized screen method in a partic

optimized chebyshev fourier migration a wide angle dual domain method for media with strong velocity contrasts a wide angle propagator is essential when imaging complex media with strong lateral velocity contrasts in one way waveequation migration we have developed a dual domain one way propagator using truncated chebyshev polynomials and a globally optimized scheme our method increases the accuracy of the expanded square root operator by adding two high order terms to the traditional split step fourier propagator first we approximate the square root operator using taylor expansion around the reference background velocity then we apply the first kindchebyshevpolynomialstoeconomizetheresultsofthe taylor expansion finally we optimize the constant coefficients using the globally optimized scheme which are fixed and feasible for arbitrary velocity models theoretical analysis and numericalexperimentshavedemonstratedthatthemethodhasvery highaccuracyandexceedstheunoptimizedfourierfinite difference propagator for the entire range of practical velocity contrasts the accurate propagation angle of the method is always about 60 under the relative error of 1 for complex media with weak moderate and even strong lateral velocity contrasts the method allows us to handle wide angle propagations and strong lateral velocity contrast simultaneously by using fourier transform alone only four 2d fourier transforms are required for each step of 3d wavefield extrapolation and the computing cost is similar to that of the fourier finite difference method compared with the finite difference method our method has no twowaysplittingerrori e azimuthal anisotropyerrorfor3dcases and almost no numerical dispersion for coarse grids in addition it has strong potential to be accelerated when an enhanced fast fouriertransformalgorithmemerges

anisotropic complex pade hybrid finite difference depth migration standard real valued finite difference fd and fourier finite difference ffd migrations cannot handle evanescent waves correctly which can lead to numerical instabilities in the presence of strong velocity variations a possible solution to these problems is the complex pade approximation which avoids problems with evanescent waves by rotating the branch cut of the complex square root we have applied this approximation to the acoustic wave equation for vertical transversely isotropic media to derive more stable fd and hybrid fd ffd migrations for such media our analysis of the dispersion relation of the new method indicates that it should provide more stable migration results with fewer artifacts and higher accuracy at steep dips our studies lead to the conclusion that the rotation angle of the branch cut that should yield the most stable image is 60 for fd migration as confirmed by numerical impulse responses and work with synthetic data

stable wide angle fourier finite difference downward extrapolation of 3 d wavefields i present an unconditionally stable implicit finite difference operator that corrects the constant velocity phase shift operator for lateral velocity variations the method is based on the fourier finite difference ffd method contrary to previous results my correction operator is stable even when the medium velocity has sharp discontinuities and the reference velocity is higher than the medium velocity the stability of the new correction enables the definition of a new downward continuation method based on the interpolation of two wavefields the first wavefield is obtained by applying the ffd correction starting from a reference velocity lower than the medium velocity the second wavefield is obtained by applying the ffd correction starting from a reference velocity higher than the medium velocity the proposed fourier finite difference plus interpolation ffdpi method combines the advantages of the ffd technique with the advantages of interpolation a simple and economical procedure for defining

directional illumination analysis using the local exponential frame we have developed an efficient method of directional illuminationanalysisinthelocalangledomainusinglocalexponential frame beamlets the space domain wavefields with different shot receiver geometries are decomposed into the local angle domain by using the local exponential beamlets which form a tight frame with the redundancy ratio two and are implemented by a linear combination of local cosine and localsinetransforms becauseofthefastalgorithmsofthelocal cosine sine transforms this method is much more efficient than the previously used decomposition methods in directional illumination analysis such as the local slant stackingmethodandthegabor daubechiesframemethod theresults of directional illumination di maps and the acquisitiondipresponsesadrforthe2dseg eagesaltmodel and the 45 shot 3d seg eage model demonstrated the validity and feasibility of our method compared with the illumination results using local slant stacking decomposition thenewmethodproducesilluminationmapsofsimilarquality but it does so a few times faster furthermore because of its high computational efficiency and saving in memory usage the new method makes the 3d directional illumination analysisreadilyapplicableintheindustry

derivation and application of extended parabolic wave theories ii path integral representations the n dimensional reduced scalar helmholtz equation for a transversely inhomogeneous medium is naturally related to parabolic propagation models through 1 the n dimensional extended parabolic weyl pseudodifferential equation and 2 an imbedding in an n 1 dimensional parabolic schrodinger equation the first relationship provides the basis for the parabolic based hamiltonian phase space path integral representation of the half space propagator the second relationship provides the basis for the elliptic based path integral representations associated with feynman and fradkin feynman and garrod and feynman and dewitt morette exact and approximate path integral constructions are derived for the homogeneous and transversely inhomogeneous cases corresponding to both narrow and wide angle extended parabolic wave theories the path integrals allow for a global perspective of the transition from elliptic to parabolic wave theory in addition to providing a unifying framework for dynamical approximations resolution of the square root operator and the concept of an underlying stochastic process

fourier finite difference wave propagation we introduce a novel technique for seismic wave extrapolation in time the technique involves cascading a fourier transform operator and a finite difference operator to form a chain operator fourier finite differences ffd we derive the ffd operator from a pseudoanalytical solution of the acoustic wave equation two dimensional synthetic examples demonstrate that the ffd operator can have high accuracy and stability in complex velocity media applying the ffd method to the anisotropic case overcomes some disadvantages of other methods such as the coupling of qp waves and qsv waves the ffd method can be applied to enhance accuracy and stability of seismic imaging by reverse time migration

prestack migration by split step dsr the double square root dsr prestack migration equation though defined for depth variable velocity can be used to image media with strong velocity variations using a phase shift plus interpolation pspi or split step correction the split step method is based on applying a phase shift correction to the extrapolated wavefield a correction that attempts to compensate for the lateral velocity variations i show how to extend dsr prestack migration to lateral velocity media and exemplify the method by applying the new algorithm to the marmousi data set the split step dsr migration is very fast and offers excellent imaging results

generalization of the bremmer coupling series an operator formalism is developed to expand the acoustic wave field in a multi dimensionally smoothly varying medium generated by a source localized in space and time into a sum of constituents each of which can be interpreted as a wave that has traveled up and down with respect to a direction of preference a definite number of times this expansion is a generalization of the bremmer coupling series the condition of smoothness of the medium relates to the width of the signature of the source in the configuration both the existence and the convergence in the weak sense of the expansion are discussed the operator calculus involved leads to a natural generalization of the concept of slowness surface to multi dimensionally smoothly varying media the operator associated with the corresponding generalized vertical slowness induces the full one way wave operator in the type of media under consideration in addition a wavefield decomposition operator as well as an interaction operator that couples the decomposed constituents are derived

wide angle one way wave equations a one way wave equation also known as a paraxial or parabolic wave equation is a differential equation that permits wave propagation in certain directions only such equations are used regularly in underwater acoustics in geophysics and as energy absorbing numerical boundary conditions the design of a one way wave equation is connected with the approximation of 1 s2 1 2 on 1 1 by a rational function which has usually been carried out by pade approximation this article presents coefficients for l2 l and other alternative classes of approximants that have better wide angle behavior for theoretical results establishing the well posedness of these wide angle equations see the work of trefethen and halpern well posedness of one way wave equations and absorbing boundary conditions math comput 47 421 435 1986

recursive integral time extrapolation methods for waves a comparative review abstractwe compared several families of algorithms for recursive integral time extrapolation rite algorithms for waves in isotropic and anisotropic media these methods allow simulating accurate wave extrapolation with little numerical dispersion even when using larger time steps than are usually possible for conventional finite difference methods these various rite algorithms all share the use of mixed space wavenumber domain operators derived from fourier integral solutions of single mode wave equations we evaluated a taxonomy for rite methods based on how they approximated the influence of medium heterogeneity one family of methods uses mixed domain series expansions to provide accurate approximations to heterogeneous extrapolators even for large time steps we compared several methods for deriving coefficients for such series approximations another family of methods uses interpolation between different homogeneous extrapolations to approximate heterogeneous time extrapolation such methods can b

representation by levels an alternative to fuzzy sets for fuzzy data mining abstract in this paper we describe and discuss the main contributions of the representation by levels approach to fuzzy data mining representation by levels is an alternative representation of fuzziness in information and data which is complementary to fuzzy sets in the sense that it provides tools and algebraic structures beyond the capabilities of fuzzy set theories based on t norms t conorms and fuzzy negations our approach allows to extend any crisp mining technique to the fuzzy case in a simple way keeping all the properties of the crisp technique we illustrate our discussion with examples and existing approaches based on representation by levels to fuzzy association rules and the related issues of mining exception anomalous rules and mining fuzzy bag databases

optimization of the parameters in complex pade fourier finite difference migration abstractcomplex pade fourier finite difference migration is a stable one way wave equation technique that allows for better treatment of evanescent modes than its real counterpart in this way producing fewer artifacts as for real fourier finite difference ffd migration its parameters can be optimized to improve the imaging of steeply dipping reflectors the dip limitation of the ffd operator depends on the variation of the velocity field we have developed a wide angle approximation for the one way continuation operator by means of optimization of the pade coefficients and the most important velocity dependent parameter we have evaluated the achieved quality of the approximate dispersion relation in dependence on the chosen function of the ratio between the model and reference velocities under consideration of the number of terms in the pade approximation and the branch cut rotation angle the optimized parameters are chosen based on the migration results and the computational cost we found that by

subsurface reflectivity estimation from imaging of primaries and multiples using amplitude normalized separated wavefields abstractseveral studies have shown the benefits of including multiple reflections together with primaries in the structural imaging of subsurface reflectors however to characterize the reflector properties there is a need to compensate for propagation effects due to multiple scattering and to properly combine the information from primaries and all orders of multiples from this perspective and based on the wave equation and rayleigh s reciprocity theorem recent works have suggested computing the subsurface image from the green s function reflection response or reflectivity by inverting a fredholm integral equation in the frequency space domain by following claerbout s imaging principle and assuming locally reacting media the integral equation may be reduced to a trace by trace deconvolution imaging condition for a complex overburden and considering that the structure of the subsurface is angle dependent this trace by trace deconvolution does not properly solve the fredholm integral equation we

one way wave equation migration in log polar coordinates we obtained acoustic wave and wavefield extrapolation equations in log polar coordinates lpcs and tried to enhance the imaging to achieve this goal it was necessary to decrease the angle between the wavefield extrapolation axis and wave propagation direction in the one way wave equation migration wem if we were unable to carry it out more reflection wave energy would be lost in the migration process it was concluded that the wavefield extrapolation operator in lpcs at low frequencies has a large wavelike region and at high frequencies it can mute the evanescent energy in these coordinate systems an extrapolation operator can readily lend itself to high order finitedifference schemes therefore even with the use of inexpensive operators wem in lpcs can clearly image varied horizontal and vertical events in complex geologic structures using wide angle and turning waves in these coordinates we did not encounter any problems with reflections from opposing dips dispersion played important roles not only as a filter operator but also as a gain function prestack and poststack migration results were obtained with extrapolation methods in different coordinate systems and it was concluded that migration in lpcs can image steeply dipping events in a much better way when compared with other methods

finite difference modeling with adaptive variable length spatial operators most finite difference simulation algorithms use fixed length spatial operators to compute spatial derivatives the choice of length is dictated by computing cost stability and dispersion criteria that are satisfied globally we propose finite difference schemes with adaptive variable length spatial operators to decrease computing costs significantly without reducing accuracy these schemes adopt long operators in regions of low velocity and short operators in regions of high velocity two methods automatically determine variable operator lengths dispersion analysis along with 1d and 2d modeling demonstrates the validity and efficiency of our schemes in addition a hybrid absorbing boundary condition helps reduce unwanted reflections from model boundaries our scheme is more efficient than those based on variable grid methods for modeling migration and inversion of models with complex velocity structures because the latter require local grid refinement which usually increases memory requirements

explicit time dependent schrodinger propagators the authors compute explicitly the time dependent schrodinger and heat propagator for the potentials lambda 2m m 1 cosh2 lambda x a delta function potential several cases of periodic delta function potentials and a delta function potential in a square well they also discuss the semiclassical approximation for the delta function potential

migration of seismic data by phase shift plus interpolation under the horizontally layered velocity assumption migration is defined by a set of independent ordinary differential equations in the wavenumber frequency domain the wave components are extrapolated downward by rotating their phases this paper shows that one can generalize the concepts of the phase shift method to media having lateral velocity variations the wave extrapolation procedure consists of two steps in the first step the wave field is extrapolated by the phase shift method using l laterally uniform velocity fields the intermediate result is l reference wave fields in the second step the actual wave field is computed by interpolation from the reference wave fields the phase shift plus interpolation pspi method is unconditionally stable and lends itself conveniently to migration of three dimensional data the performance of the methods is demonstrated on synthetic examples the pspi migration results are then compared with those obtained from a finite difference method

full wave directional illumination analysis in the frequency domain directional illumination analysis based on one way wave equations has been studied extensively however its inherent limitations e g one way propagation wide angle error and amplitude inaccuracy can severely hinder its applications for accurate survey design and true reflection imaging corrections in complex media we have analyzed the illumination in the frequency domain using full two way wave propagators considering the extensive computation and huge storage required for time domain methods and the fact that the illumination is frequency dependent this full wave analysis can provide frequency dependent full angle true amplitude illumination not only for the downgoing waves but also for the upgoing waves including turning waves and reflected waves two methods were considered to decompose the full wavefield into the local angle domain a direct full dimensional decomposition and more efficient split step decomposition composed of three lower dimensional decompositions the results of illumination analysis demonstrated the advantages of this method the two decomposition methods produced similar results