modelling damped acoustic waves by a dissipation preserving conformal symplectic method we propose a novel stable and efficient dissipation preserving method for acoustic wave propagations in attenuating media with both correct phase and amplitude through introducing the conformal multi symplectic structure the intrinsic dissipation law and the conformal symplectic conservation law are revealed for the damped acoustic wave equation the proposed algorithm is exactly designed to preserve a discrete version of the conformal symplectic conservation law more specifically two subsystems in conjunction with the original damped wave equation are derived one is actually the conservative hamiltonian wave equation and the other is a dissipative linear ordinary differential equation ode system standard symplectic method is devoted to the conservative system whereas the analytical solution is obtained for the ode system an explicit conformal symplectic scheme is constructed by concatenating these two parts of solutions by the strang splitting technique stability analysis and convergence tests are given thereafter a benchmark model in homogeneous media is presented to demonstrate the effectiveness and advantage of our method in suppressing numerical dispersion and preserving the energy dissipation further numerical tests show that our proposed method can efficiently capture the dissipation in heterogeneous media

application of the perfectly matched layer pml absorbing boundary condition to elastic wave propagation a method is presented for application of the perfectly matched layer pml absorbing boundary condition abc to the p sv velocity stress finite difference method the pml consists of a nonphysical material containing both passive loss and dependent sources that provides active absorption of fields it has been used in electromagnetic applications where it has provided excellent results for a wide range of angles and frequencies in this work numerical simulations are used to compare the pml and an optimal second order elastic abc peng and toksoz j acoust soc am 95 733 745 1994 reflection factors are used to compare angular performance for continuous wave illumination snapshots of potentials are used to compare performance for broadband illumination these comparisons clearly demonstrate the superiority of the pml formulation within the pml there is a 60 increase in the number of unknowns per grid cell relative to the velocity stress formulation however the high quality of the pml abc allows the use of a smaller grid which can result in a lower overall computational cost

on block triangular preconditioned iteration methods for solving the helmholtz equation abstract to further improve the efficiency of solving the helmholtz equation in heterogeneous media with large wavenumber the krylov subspace methods incorporated with a class of inexact rotated block triangular preconditioners are presented to solve a block two by two linear system derived from the discrete helmholtz equation we further develop the eigenvalue properties of the preconditioned matrices to discuss the convergence of the corresponding preconditioned iteration methods the superiority of such preconditioned iteration methods is prominent according to the numerical results when comparing with other classical iteration methods we also investigate how the wavenumber influences the performance of the corresponding methods and it is shown that the iteration number of our proposed methods linearly increase with the wavenumber roughly furthermore the computational wave fields which conform the real physical law are exhibited to show the correctness of our proposed numerical modeling algorithm

another look at the regression discontinuity design the attractiveness of the regression discontinuity design rdd rests on its close similarity to a formal experimental design on the other hand it is of limited applicability since it is not often the case that units are assigned to the treatment group on the basis of an observable to the analyst pre program measure besides it only allows to identify the mean impact on a very specific sub population in this paper we show that the rdd straightforwardly generalizes to the instances in which the units eligibility is established on an observable pre program measure with eligible units allowed to freely self select into the program this set up also turns out very convenient to build a specification test on conventional non experimental estimators of the program mean impact data requirements are made explicit

p sv wave propagation in heterogeneous media velocity stress finite difference method i present a finite difference method for modeling p sv wave propagation in heterogeneous media this is an extension of the method i previously proposed for modeling sh wave propagation by using velocity and stress in a discrete grid the two components of the velocity cannot be defined at the same node for a complete staggered grid the stability condition and the p wave phase velocity dispersion curve do not depend on the poisson s ratio while the s wave phase velocity dispersion curve behavior is rather insensitive to the poisson s ratio therefore the same code used for elastic media can be used for liquid media where s wave velocity goes to zero and no special treatment is needed for a liquid solid interface typical physical phenomena arising with p sv modeling such as surface waves are in agreement with analytical results the weathered layer and corner edge models show in seismograms the same converted phases obtained by previous authors this method gives stable results for step discontinuities as shown for a liquid layer above an elastic half space the head wave preserves the correct amplitude finally the corner edge model illustrates a more complex geometry for the liquid solid interface as the poisson s ratio v increases from 0 25 to 0 5 the shear converted phases are removed from seismograms and from the time section of the wave field

attenuation compensation in anisotropic least squares reverse time migration abstractanisotropic and attenuating properties of subsurface media cause amplitude loss and waveform distortion in seismic wave propagation resulting in negative influence on seismic imaging to correct the anisotropy effect and compensate amplitude attenuation a compensated amplitude vertical transverse isotropic vti least squares reverse time migration lsrtm method is adopted in this method the attenuation term of an attenuated acoustic wave equation is extended to a vti quasi differential wave equation which takes care of effects from anisotropy and attenuation the finite difference method is used to solve the equation in which attenuation terms are solved in the wavenumber domain and other terms are solved in the space or wavenumber domain stable regularization operators are derived and introduced to the equations to eliminate severe numerical noise in high frequency components during backward propagation meanwhile a demigration operator migration operator and gradient formula for att

convergence analysis of moving finite element methods for space fractional differential equations most existent papers have focused on the fixed mesh methods for solving space fractional differential equations however since some classes of space fractional differential equations may have singular or even finite time blowup solutions it is highly needed to develop adaptive mesh methods to solve these problems in this paper the moving finite element methods are studied for a class of time dependent space fractional differential equations the convergence theories of the methods are derived with l 2 norm and numerical examples are provided to support the theoretical results to simplify the analysis a fractional ritz projection operator is introduced and the error estimation of the projection is derived under the moving mesh setting

accurate simulations of pure quasi p waves in complex anisotropic media abstractreverse time migration rtm in complex anisotropic media requires calculation of the propagation of a single mode wave the quasi p wave this was conventionally realized by solving a 2 2 system of second order partial differential equations the implementation of this 2 2 system required at least twice the computational resources as compared with the acoustic wave equation the s waves an introduced auxiliary function in this system were treated as artifacts in the rtm furthermore the 2 2 system suffered numerical stability problems at the places in which abrupt changes of symmetric axis of anisotropy exist which brings more challenges to real data implementation on the other hand the alkhalifah s equation which governs the pure quasi p wave propagation was hard to solve because it was a pseudodifferential equation we proposed a pure quasi p wave equation that can be easily implemented within current imaging framework our new equation was obtained by decomposing the original pseudodif

efficient wave mode separation in vertical transversely isotropic media abstractwave mode separation can be achieved by projecting elastic wavefields onto mutually orthogonal polarization directions in isotropic media because the p wave s polarization vectors are consistent with wave vectors the isotropic separation operators are represented by divergence and curl operators which are easy to realize in anisotropic media polarization vectors deviate from wave vectors based on local anisotropic strength and separation operators lose their simplicity conventionally anisotropic wave mode separation is implemented either by direct filtering in the wavenumber domain or nonstationary filtering in the space domain which are computationally expensive moreover in conventional anisotropic separation correcting for amplitude and phase changes of waveforms by applying separation operators is also more difficult than in an isotropic case we have developed new operators for efficient wave mode separation in vertical transversely isotropic vti media our separation operators a

multidirectional vector based elastic reverse time migration and angle domain common image gathers with approximate wavefield decomposition of p and s waves abstractelastic reverse time migration e rtm has limitations when the migration velocities contain strong contrasts first the traditional scheme of p s wave mode separation is based on helmholtz s equations which ignore the conversion between p and s waves at the current separation time thus it contains an implicit assumption of the constant shear modulus and requires smoothing the heterogeneous model to approximately satisfy a locally constant condition second the vector based imaging condition needs to use the reflection image normal and it also cannot give the correct polarity of the pp image in all possible conditions third the angle domain common image gathers adcigs calculated using the poynting vectors pvs do not consider the wave interferences that happen at each reflector therefore smooth models are often used for e rtm we relax this condition by proposing an improved data flow that involves three new contributions the first contribution is an improved system of p s wave mode

on simple moving grid methods for one dimensional evolutionary partial differential equations abstract two moving grid algorithms for the numerical integration of evolutionary one dimensional partial differential equations are considered both algorithms discretize the equation on trapezoidal space time elements and combine a prediction step with a remeshing technique in order to determine the orientation of the sides of the trapezoids joining nodes at consecutive time levels the behavior of the discretizations employed is forecast by means of the modified equation technique and then studied in a series of numerical experiments

assessment of rod behaviour theories used in spectral finite element modelling abstract in this work different theories of rods have been discussed and compared the investigated theories are widely used in spectral finite element modelling of rod behaviour associated with propagation of symmetric longitudinal waves these are various single two mode and three mode theories including the elementary classical love and mindlin herrmann approaches as well as new two three and four mode theories proposed by the authors dispersion curves associated with each theory obtained by the use of hamilton s principle have been presented and discussed in the paper the investigation programme carried out by the authors aimed to show major differences and similarities between the rod theories and to discuss certain numerical aspects of their application great attention has been paid on properties limitations as well as difficulties associated with the use of the theories the results obtained from a wide program on numerical tests allowed the authors to draw certain general conclusions that are valid not only in the field of the spectral finite element method but also in the field of dynamics of engineering rod structures

boundary element and finite element coupling for aeroacoustics simulations we consider the scattering of acoustic perturbations in the presence of a flow we suppose that the space can be split into a zone where the flow is uniform and a zone where the flow is potential in the first zone we apply a prandtl glauert transformation to recover the helmholtz equation the well known setting of boundary element method for the helmholtz equation is available in the second zone the flow quantities are space dependent we have to consider a local resolution namely the finite element method herein we carry out the coupling of these two methods and present various applications and validation test cases the source term is given through the decomposition of an incident acoustic field on a section of the computational domain s boundary validations against analytic another numerical method and measurements on different test cases are presented

the pseudospectral method comparisons with finite differences for the elastic wave equation the pseudospectral or fourier method has been used recently by several investigators for forward seismic modeling the method is introduced here in two different ways as a limit of finite differences of increasing orders and by trigonometric interpolation an argument based on spectral analysis of a model equation shows that the pseudospectral method for the accuracies and integration times typical of forward elastic seismic modeling may require in each space dimension as little as a quarter the number of grid points compared to a fourth order finite difference scheme and one sixteenth the number of points as a second order finite difference scheme for the total number of points in two dimensions these factors become 1 16 and 1 256 respectively in three dimensions they become 1 64 and 1 4 096 respectively in a series of test calculations on the two dimensional elastic wave equation only minor degradations are found in cases with variable coefficients and discontinuous interfaces

mimetic seismic wave modeling including topography on deformed staggered grids abstractfinite difference methods for modeling seismic waves are known to be inaccurate when including a realistic topography due to the large dispersion errors that appear in the modelled surface waves and the scattering introduced by the staircase approximation to the topography as a consequence alternatives to finite difference methods have been proposed to circumvent these issues we present a new numerical scheme for 3d elastic wave propagation in the presence of strong topography this finite difference scheme is based upon a staggered grid of the lebedev type or fully staggered grid fsg it uses a grid deformation strategy to make a regular cartesian grid conform to a topographic surface in addition the scheme uses a mimetic approach to accurately solve the free surface condition and hence allows for a less restrictive grid spacing criterion in the computations the scheme can use high order operators for the spatial derivatives and obtain low dispersion results with as few as six points pe

grid dispersion and stability criteria of some common finite element methods for acoustic and elastic wave equations purely numerical methods based on finite element approximation of the acoustic or elastic wave equation are becoming increasingly popular for the generation of synthetic seismograms we present formulas for the grid dispersion and stability criteria for some popular finite element methods fem for wave propagation namely classical and spectral fem we develop an approach based on a generalized eigenvalue formulation to analyze the dispersive behavior of these fems for acoustic or elastic wave propagation that overcomes difficulties caused by irregular node spacing within the element and the use of high order polynomials as is the case for spectral fem analysis reveals that for spectral fem of order four or greater dispersion is less than 0 2 at four to five nodes per wavelength and dispersion is not angle dependent new results can be compared with grid dispersion results of some classical finite difference methods fdm used for acoustic or elastic wave propagation analysis reveals that fdm and classical fem require a larger sampling ratio than a spectral fem to obtain results with the same degree of accuracy the staggered grid fdm is an efficient scheme but the dispersion is angle dependent with larger values along the grid axes on the other hand spectral fem of order four or greater is isotropic with small dispersion making it attractive for simulations with long propagation times

finite difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid we describe the application of the rotated staggered grid rsg finite difference technique to the wave equations for anisotropic and viscoelastic media the rsg uses rotated finite difference operators leading to a distribution of modeling parameters in an elementary cell where all components of one physical property are located only at one single position this can be advantageous for modeling wave propagation in anisotropic media or complex media including high contrast discontinuities because no averaging of elastic moduli is needed the rsg can be applied both to displacement stress and to velocity stress finite difference fd schemes whereby the latter are commonly used to model viscoelastic wave propagation with a von neumann style anlysis we estimate the dispersion error of the rsg scheme in general anisotropic media in three different simulation examples all based on previously published problems we demonstrate the application and the accuracy of the proposed numerical approach

an efficient symplectic reverse time migration method using a local nearly analytic discrete operator in acoustic transversely isotropic media with a vertical symmetry axis abstractreverse time migration rtm has become an extremely important imaging method due to its ability of handling complicated velocity models in recent years there has been an increasing need to perform rtm with anisotropic velocity models such as transversely isotropic media with a vertical symmetry axis vti finite difference fd is one of the most widely used numerical methods in rtm however it often suffers from serious numerical dispersion resulting in significant loss of imaging resolution over the years a nearly analytic discrete operator has been developed to approximate the high order spatial partial differential operators which is effective in suppressing numerical dispersion using this operator we developed a nearly analytic symplectic nas method to numerically solve acoustic vti wave equations which preserves the symplectic structure we have investigated its properties including the stability condition numerical dispersion relation and computational efficiency also we a

maximum norm a posteriori error estimates for a one dimensional convection diffusion problem a singularly perturbed quasi linear two point boundary value problem with an exponential boundary layer is discretized on arbitrary nonuniform meshes using first and second order difference schemes including upwind schemes we give first and second order maximum norm a posteriori error estimates that are based on difference derivatives of the numerical solution and hold true uniformly in the small parameter numerical experiments support the theoretical results

accurate and efficient seismic modeling in random media the optimum method for seismic modeling in random media must 1 be highly accurate to be sensitive to subtle effects of wave propagation 2 allow coarse sampling to model media that are large compared to the scale lengths and wave propagation distances which are long compared to the wavelengths this is necessary to obtain statistically meaningful overall attributes of wavefields high order staggered grid finite difference algorithms and the pseudospectral method combine high accuracy in time and space with coarse sampling investigations for random media reveal that both methods lead to nearly identical wavefields the small differences can be attributed mainly to differences in the numerical dispersion this result is important because it shows that errors of the numerical differentiation which are caused by poor polynomial interpolation near discontinuities do not accumulate but cancel in a random medium where discontinuities are numerous the differentiator can be longer than the medium scale leng

sh wave propagation in heterogeneous media velocity stress finite difference method a new finite difference fd method is presented for modeling sh wave propagation in a generally heterogeneous medium this method uses both velocity and stress in a discrete grid density and shear modulus are similarly discretized avoiding any spatial smoothing therefore boundaries will be correctly modeled under an implicit formulation standard problems quarter plane propagation sedimentary basin propagation are studied to compare this method with other methods finally a more complex example a salt dome inside a two layered medium shows the effect of lateral propagation on seismograms recorded at the surface a corner wave always in phase with the incident wave and a head wave will appear which will pose severe problems of interpretation with the usual vertical migration methods

analysis of moving mesh partial differential equations with spatial smoothing two moving mesh partial differential equations mmpdes with spatial smoothing are derived based upon the equidistribution principle this smoothing technique is motivated by the robust moving mesh method of dorfi and drury j comput phys 69 1987 pp 175 195 it is shown that under weak conditions the basic property of no node crossing is preserved by the spatial smoothing and a local quasi uniformity property of the coordinate transformations determined by these mmpdes is proven it is also shown that discretizing the mmpdes using centered finite differences these basic properties are preserved

source independent full waveform inversion of seismic data a set of seismic trace data is collected in an input data set that is first fourier transformed in its entirety into the frequency domain a normalized wavefield is obtained for each trace of the input data set in the frequency domain normalization is done with respect to the frequency response of a reference trace selected from the set of seismic trace data the normalized wavefield is source independent complex and dimensionless the normalized wavefield is shown to be uniquely defined as the normalized impulse response provided that a certain condition is met for the source this property allows construction of the inversion algorithm disclosed herein without any source or source coupling information the algorithm minimizes the error between data normalized wavefield and the model normalized wavefield the methodology is applicable to any 3 d seismic problem and damping may be easily included in the process

modeling of pseudoacoustic p waves in orthorhombic media with a low rank approximation abstractwavefield extrapolation in pseudoacoustic orthorhombic anisotropic media suffers from wave mode coupling and stability limitations in the parameter range we use the dispersion relation for scalar wave propagation in pseudoacoustic orthorhombic media to model acoustic wavefields the wavenumber domain application of the laplacian operator allows us to propagate the p waves exclusively without imposing any conditions on the parameter range of stability it also allows us to avoid dispersion artifacts commonly associated with evaluating the laplacian operator in space domain using practical finite difference stencils to handle the corresponding space wavenumber mixed domain operator we apply the low rank approximation approach considering the number of parameters necessary to describe orthorhombic anisotropy the low rank approach yields space wavenumber decomposition of the extrapolator operator that is dependent on space location regardless of the parameters a feature necessary for orthorhomb

computationally efficient three dimensional acoustic finite difference frequency domain seismic modeling in vertical transversely isotropic media with sparse direct solver abstractthe computational burden of frequency domain full waveform inversion fwi of wide aperture fixed spread data is conventionally reduced by limiting the inversion to a few discrete frequencies in this framework frequency domain seismic modeling is performed efficiently for multiple sources by solving the linear system resulting from the discretization of the time harmonic wave equation with the massively parallel sparse direct solver frequency domain seismic modeling based on the sparse direct solver dsfdm requires specific design finite difference stencils of compact support to minimize the computational cost of the lower upper decomposition of the impedance matrix in terms of memory demand and floating point operations a straightforward adaptation of such finite difference stencil originally developed for the isotropic acoustic wave equation is proposed to introduce vertical transverse isotropy vti in the modeling without any extra computational cost the method relies on a fourth ord

numerical solution of a singularly perturbed two point boundary value problem using equidistribution analysis of convergence abstract adaptive grid methods are becoming established as valuable computational techniques for the numerical solution of differential equations with near singular solutions adaptive methods are equally effective in approximating solutions of problems with boundary layers or interior layers see for example mulholland et al siam j sci comput 19 4 1998 1261 1289 much is now being done in developing error analyes for methods that are based on adaptivity in this paper we present a rigorous error analysis for the solution of a singularly perturbed two point boundary value problem on a grid that is constructed adaptively from a knowledge of the exact solution the discrete solutions are generated by an upwind finite difference scheme and the grid is formed by equidistributing a monitor function based on arc length an error analysis shows that the discrete solutions are uniformly convergent with respect to the perturbation parameter epsilon the epsilon uniform convergence is confirmed by numerical computations

from association to causation in observational studies the role of tests of strongly ignorable treatment assignment abstract if treatment assignment is strongly ignorable then adjustment for observed covariates is sufficient to produce consistent estimates of treatment effects in observational studies a general approach to testing this critical assumption is developed and applied to a study of the effects of nuclear fallout on the risk of childhood leukemia r a fisher s advice on the interpretation of observational studies was make your theories elaborate formally make causal theories sufficiently detailed that under the theory strongly ignorable assignment has testable consequences

stable p wave modeling for reverse time migration in tilted ti media we present an approach for p wave modeling in inhomogeneous transversely isotropic media with tilted symmetry axis tti media suitable for anisotropic reverse time migration the proposed approach is based on wave equations derived from first principles the equations of motion and hooke s law under the acoustic ti approximation consequently no assumptions are made about the spatial variation of medium parameters a rotation of the stress and strain tensors to a local coordinate system aligned with the ti symmetry axis makes it possible to benefit from the simple and sparse form of the ti elastic tensor in that system the resulting wave equations can be formulated either as a set of five first order or as a set of two second order partial differential equations for the constant density case the second order tti wave equations involve mixed and nonmixed second order spatial derivatives with respect to global nonrotated coordinates we propose a numerical implementation of these equations using high order centered finite differences to minimize modeling artifacts related to the use of centered first derivative operators we use discrete second derivative operators for the nonmixed second order spatial derivatives and repeated discrete first derivative operators for the mixed derivatives such a combination of finite difference operators leads to a stable wave propagator provided that the operators are designed properly in practice stability is achieved by slightly weighting down terms that contain mixed derivatives this has a minor practically negligible effect on the kinematics of wave propagation the stability of the presented scheme in inhomogeneous tti models with rapidly varying anisotropic symmetry axis direction is demonstrated with numerical examples

massively parallel structured direct solver for equations describing time harmonic qp polarized waves in tti media abstractwe considered the discretization and approximate solutions of equations describing time harmonic qp polarized waves in 3d inhomogeneous anisotropic media the anisotropy comprises general tilted transversely isotropic symmetries we are concerned with solving these equations for a large number of different sources we considered higher order partial differential equations and variable order finite difference schemes to accommodate anisotropy on the one hand and allow higher order accuracy to control sampling rates for relatively high frequencies on the other hand we made use of a nested dissection based domain decomposition in a massively parallel multifrontal solver combined with hierarchically semiseparable matrix compression techniques the higher order partial differential operators and the variable order finite difference schemes require the introduction of separators with variable thickness in the nested dissection the development of these and their integration with the multifrontal

multiscale adjoint waveform difference tomography using wavelets abstractfull waveform seismic inversions based on minimizing the distance between observed and predicted seismograms are in principle able to yield better resolved earth models than those minimizing misfits derived from traveltimes alone adjoint based methods provide an efficient way of calculating the gradient of the misfit function via a sequence of forward modeling steps which using spectral element codes can be carried out in realistically complex media convergence and stability of full waveform difference adjoint schemes are greatly improved when data and synthetics are progressively presented to the algorithms in a constructive multiscale approximation using a bi orthogonal wavelet transform wavelets provide the nonredundant spectral decomposition that paves the way for the inversion to proceed successively from long wavelength fitting to detailed exploration of the phases in the seismogram the choice of wavelet class and type the initial depth of the multiscale decomposition and the min

3 d finite difference modeling of elastic waves in borehole environments full wavefield sonic logs may be simulated for arbitrarily complicated three dimensional 3 d borehole environments using 3 d elastic finite differences to ensure reliability through large contrasts in poisson s ratio across mud casing lithology contacts a staggered grid finite difference formulation is used cylindrical symmetry is not assumed so the responses of features such as asymmetrical washouts and dipping structure are easily obtained when features are asymmetrical seismic responses vary significantly at various points around the hole circumference at any depth even for simple hole geometries observed responses are complicated because of coupling between waves inside and outside the hole observations are also sensitive to the source receiver separation output formats include fixed time snapshots of displacement divergence and curl and seismograms for the center and edge of a borehole this allows detailed arrival identification and interpretation to our knowledge this is the most

convergence analysis of moving godunov methods for dynamical boundary layers in this paper a godunov scheme on moving meshes is studied for kinds of time dependent convection dominated equations with dynamical boundary layers the stability and a second order spatial convergence are proved a numerical example is provided to confirm the theoretical results

adaptive mesh movement the mmpde approach and its applications a class of moving mesh algorithms based upon a so called moving mesh partial differential equation mmpde is reviewed various forms for the mmpde are presented for both the simple one and the higher dimensional cases additional practical features such as mesh movement on the boundary selection of the monitor function and smoothing of the monitor function are addressed the overall discretization and solution procedure including for unstructured meshes are briefly covered finally we discuss some physical applications suitable for mmpde techniques and some challenges facing mmpde methods in the future

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

wavefield extrapolation in pseudodepth domain abstractwavefields are commonly computed in the cartesian coordinate frame its efficiency is inherently limited due to spatial oversampling in deep layers where the velocity is high and wavelengths are long to alleviate this computational waste due to uneven wavelength sampling we convert the vertical axis of the conventional domain from depth to vertical time or pseudodepth this creates a nonorthognal riemannian coordinate system isotropic and anisotropic wavefields can be extrapolated in the new coordinate frame with improved efficiency and good consistency with cartesian domain extrapolation results prestack depth migrations are also evaluated based on the wavefield extrapolation in the pseudodepth domain

ineligibles and eligible non participants as a double comparison group in regression discontinuity designs in a sharp regression discontinuity design rdd the participation status deterministically depends on whether a pre programme characteristic is above or below a specified threshold the attractiveness of such a design rests on close similarities with a formal experiment nevertheless it is of limited applicability since participation into a programme is seldom determined according to this rule besides in the presence of heterogeneous effects a sharp rdd only allows identification of mean effects for individuals around the threshold for participation two results are presented in this paper and they both partially overcome the two limitations described above we show that when individuals self select into participation conditional on some eligibility criteria a sharp rdd provides a natural framework to define a specification test for the non experimental estimation of programme effects for participants away from the threshold we also show that in this set up the regularity conditions required for the identification of the mean counterfactual outcome for participants marginally eligible for the programme are essentially the same as in a sharp rdd

the central role of the propensity score in observational studies for causal effects abstract the results of observational studies are often disputed because of nonrandom treatment assignment for example patients at greater risk may be overrepresented in some treatment group this paper discusses the central role of propensity scores and balancing scores in the analysis of observational studies the propensity score is the estimated conditional probability of assignment to a particular treatment given a vector of observed covariates both large and small sample theory show that adjustment for the scalar propensity score is sufficient to remove bias due to all observed covariates applications include matched sampling on the univariate propensity score which is equal percent bias reducing under more general conditions than required for discriminant matching multivariate adjustment by subclassification on balancing scores where the same subclasses are used to estimate treatment effects for all outcome variables and in all subpopulations and visual representation of multivariate adjustment by a two dimensional plot author

two dimensional nonlinear inversion of seismic waveforms numerical results the nonlinear problem of inversion of seismic waveforms can be set up using least squares methods the inverse problem is then reduced to the problem of minimizing a lp nonquadratic function in a space of many 104to106 variables using gradient methods leads to iterative algorithms each iteration implying a forward propagation generated by the actual sources a backward propagation generated by the data residuals acting as if they were sources and a correlation at each point of the space of the two fields thus obtained which gives the updated model the quality of the results of any inverse method depends heavily on the realism of the forward modeling finite difference schemes are a good choice relative to realism because although they are time consuming they give excellent results numerical tests performed with multioffset synthetic data from a two dimensional model prove the feasibility of the approach if only surface recorded reflections are used the high spatial frequency content of the

elastic imaging with exact wavefield extrapolation for application to ocean bottom 4c seismic data abstracta central component of imaging methods is receiver side wavefield backpropagation or extrapolation in which the wavefield from a physical source scattered at any point in the subsurface is estimated from data recorded by receivers located near or at the earth s surface elastic reverse time migration usually accomplishes wavefield extrapolation by simultaneous reversed time injection of the particle displacements or velocities recorded at each receiver location into a wavefield modeling code here we formulate an exact integral expression based on reciprocity theory that uses a combination of velocity stress recordings and quadrupole dipole backpropagating sources rather than the commonly used approximate formula involving only particle velocity data and dipole backpropagating sources the latter approximation results in two types of nonphysical waves in the scattered wavefield estimate first each arrival contained in the data is injected upward and downward rather than unidirectionally as

moving mesh partial differential equations mmpdes based on the equidistribution principle this paper considers several moving mesh partial differential equations that are related to the equidistribution principle several of these are new and some correspond to discrete moving mesh equations that have been used by others their stability is analyzed and it is seen that a key term for most of these moving mesh pdes is a source like term that measures the level of equidistribution it is shown that under weak assumptions mesh crossing cannot occur for most of them finally numerical experiments for these various moving mesh pdes are performed to study their relative properties

a robust adaptive method for a quasi linear one dimensional convection diffusion problem a quasi linear conservative convection diffusion two point boundary value problem is considered to solve it numerically an upwind finite difference scheme is applied the mesh used has a fixed number n 1 of nodes and is initially uniform but its nodes are moved adaptively using a simple algorithm of de boor based on equidistribution of the arc length of the current computed piecewise linear solution it is proved for the first time that a mesh exists that equidistributes the arc length along the polygonal solution curve and that the corresponding computed solution is first order accurate uniformly in varepsilon where varepsilon is the diffusion coefficient in the case when the boundary value problem is linear if n is sufficiently large independently of varepsilon it is shown that after o rm ln 1 varepsilon rm ln n iterations of the algorithm the piecewise linear interpolant of the computed solution achieves first order accuracy in the l infty 0 1 norm uniformly in varepsilon numerical experiments are presented that support our theoretical results

application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media we present and analyze a perfectly matched absorbing layer model for the velocity stress formulation of elastodynamics the principal idea of this method consists of introducing an absorbing layer in which we decompose each component of the unknown into two auxiliary components a component orthogonal to the boundary and a component parallel to it a system of equations governing these new unknowns then is constructed a damping term finally is introduced for the component orthogonal to the boundary this layer model has the property of generating no reflection at the interface between the free medium and the artificial absorbing medium in practice both the boundary condition introduced at the outer boundary of the layer and the dispersion resulting from the numerical scheme produce a small reflection which can be controlled even with very thin layers as we will show with several experiments this model gives very satisfactory results namely the reflection coefficient even in the case of heterogeneous anisotropic media is about 1 for a layer thickness of five space discretization steps

elastic least squares reverse time migration via linearized elastic full waveform inversion with pseudo hessian preconditioning abstracttime domain elastic least squares reverse time migration lsrtm is formulated as a linearized elastic full waveform inversion problem the elastic born approximation and elastic reverse time migration rtm operators are derived from the time domain continuous adjoint state method our approach defines p and s wave impedance perturbations as unknown elastic images our algorithm is obtained using continuous functional analysis in which the problem is discretized at the final stage optimize before discretize approach the discretized numerical versions of the elastic born operator and its adjoint elastic rtm operator pass the dot product test the conjugate gradient least squares method is used to solve the least squares migration quadratic optimization problem in other words the hessian operator for elastic lsrtm is implicitly inverted via a matrix free algorithm that only requires the action of forward and adjoint operators on vectors the diagonal of the pseudo hessian operator is used t

centered and staggered fourier derivatives and hilbert transforms when ringing derivatives are used in seismic modeling migration and other wave propagation problems in heterogeneous media they interact with boundaries of blocks or layers of different physical properties generating noncausal and nonphysical waves these spurious waves can be as strong as reflected or converted waves and other legitimate physical events leading to errors and misinterpretations therefore compact i e short and nonringing derivatives should be used

coupled equations for reverse time migration in transversely isotropic media reverse time migration rtm images reflectors by using time extrapolation modeling codes to synthesize source and receiver wavefields in the subsurface asymptotic analysis of wave propagation in transversely isotropic ti media yields a dispersion relation describing coupled p and sv wave modes this dispersion relation can be converted into a fourth order scalar partial differential equation pde increased computational efficiency can be achieved using equivalent coupled second order pdes analysis of the corresponding dispersion relations as matrix eigenvalue systems allows one to characterize all possible coupled linear second order systems equivalent to a given linear fourth order pde and to determine which ones yield optimally efficient finite difference implementations setting the shear velocity along the axis of symmetry to zero yields a simpler approximate ti wave equation that is more efficient to implement this simpler approximation however can become unstable for some plausible combinations of anisotropic parameters the same eigensystem analysis can be applied using finite vertical shear velocity to obtain solutions that avoid these instability problems

2d isotropic elastic gaussian beam migration for common shot multicomponent records gaussian beam migration is flexible and adaptable for imaging geologically complicated areas it avoids some limitations of traditional ray based migration approaches such as amplitude singularity in caustic zones and inability to image multiple arrivals previous studies on gaussian beam migration mainly focused on acoustic media we develop a 2d isotropic elastic gaussian beam migration scheme for common shot multicomponent records the proposed method extrapolates p and s mode wavefields simultaneously using the kirchhoff helmholtz integral solution of the isotropic elastodynamic equation we separate the extrapolated wavefields into compressional and shear modes using the helmholtz decomposition and then implement a modified dot product imaging condition this approach enables us to produce clear pp images and avoid polarity reversal issues for ps images in addition based on the theory of wavefield approximation in the effective vicinity of central rays we derive a formula to compute the propaga

dispersion relation preserving stereo modeling method beyond nyquist frequency for acoustic wave equation abstractaccording to nyquist shannon sampling theory two sampling points of the displacements per minimal wavelength are required to reconstruct a wavefield however high frequency components above the nyquist frequency can be theoretically recovered by introducing spatial derivatives of the displacements we have developed a new finite difference scheme for solving the acoustic equation named the dispersion relation preserving stereo modeling method which is practical for seismic forward modeling beyond the nyquist frequency the key idea in this model is to approximate the high order spatial derivatives with wavenumber domain optimization in wavefields represented by the wave displacements and their gradients we investigated the theoretical properties of this method including the absolute error of spatial operators numerical error stability criterion and computational efficiency our results indicated that the method can effectively suppress the numerical dispersion on extremely coarse grids wi

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

numerical analysis of a narrow angle one way elastic wave equation and extension to curvilinear coordinates in this paper we review the finite difference implementation of a narrow angle one way vector wave equation for elastic 3d media extrapolation is performed in the frequency domain where the second order accurate lateral spatial difference operators are sufficiently accurate for narrow angle propagation we perform a numerical analysis of the finite difference scheme to highlight the stability and dispersion characteristics the von neumann stability criterion indicates that extracting a reference phase during the extrapolation step noticeably improves the forward marching scheme and dispersion analysis shows that numerical grid anisotropy is minimal for the propagation path lengths source pulse spectral content and angular range of forward propagation of interest although the algorithm is reasonable its computational efficiency is limited by the second order accurate extrapolation step therefore the extrapolation scheme can be improved we extend the cartesian narrow angle formulation to curvilinear coordinates where the computational grid tracks the true wavefront in a reference medium and the wavefield derivative normal to the reference wavefront is evaluated locally using the cartesian propagator an example of curvilinear extrapolation for a simple model consisting of a high velocity sphere within a homogeneous background velocity structure shows that the narrow angle propagator is capable of modeling frequency dependent geometric spreading and diffraction effects in curvilinear coordinates

unstructured grid method for stress wave propagation in elastic media a new numerical method named unstructured grid method is presented to calculate two dimensional stress wave propagation problems in elastic media this method is based on the dynamic equilibrium equations of the investigated lumps that are formed among the auxiliary unstructured grids the approach to the solution is obtained by calculating the node displacements and the central point stresses of the spatial grids alternately the stability criterion of the method is given and tested numerically analysis of dispersion for numerical phase velocity is studied the rule is provided to determine the spatial step size of discretization mesh so as to ensure the numerical accuracy of the unstructured grid method numerical results are compared with the available exact solutions and with the solutions of central difference method for demonstrating the good accuracy of the new method the unstructured grid method possesses higher calculating speed than the central difference method finally the dynamic behavior of a square plate with a hole under an impact loading is researched

on the stability of mesh equidistribution strategies for time dependent partial differential equations abstract we study the stability of several mesh equidistribution schemes for time dependent partial differential equations in one space dimension the schemes move a finite difference or finite element mesh so that a given quantity is uniform over the domain we consider mesh moving methods that are based on solving a system of ordinary differential equations for the mesh velocities and show that some of these methods are unstable with respect to an equidistributing mesh when the partial differential system is dissipative using linear perturbation techniques we are able to develop simple criteria for determining the stability of a particular method and show how to construct stable differential systems for the mesh velocities several examples illustrating stable and unstable mesh motions are presented

3d weak dispersion reverse time migration using a stereo modeling operator reliable 3d imaging is a required tool for developing models of complex geologic structures reverse time migration rtm as the most powerful depth imaging method has become the preferred imaging tool because of its ability to handle complex velocity models including steeply dipping interfaces and large velocity contrasts finite difference methods are among the most popular numerical approaches used for rtm however these methods often encounter a serious issue of numerical dispersion which is typically suppressed by reducing the grid interval of the propagation model resulting in large computation and memory requirements in addition even with small grid spacing numerical anisotropy may degrade images or worse provide images that appear to be focused but position events incorrectly recently stereo operators have been developed to approximate the partial differential operator in space these operators have been used to develop several weak dispersion and efficient stereo modeling methods that have been found to be superior to conventional algorithms in suppressing numerical dispersion and numerical anisotropy we generalized one stereo modeling method fourth order nearly analytic central difference nacd from 2d to 3d and applied it to 3d rtm the rtm results for the 3d seg eage phase a classic data set 1 and the seg advanced modeling project model demonstrated that even when using a large grid size the nacd method can handle very complex velocity models and produced better images than can be obtained using the fourth order and eighth order lax wendroff correction lwc schemes we also applied 3d nacd and fourth order lwc to a field data set and illustrated significant improvements in terms of structure imaging horizon layer continuity and positioning we also investigated numerical dispersion and found that not only does the nacd method have superior dispersion characteristics but also that the angular variation of dispersion is significantly less than for lwc

modeling of pure acoustic wave in tilted transversely isotropic media using optimized pseudo differential operators abstract conventional modeling and imaging for tilted transversely isotropic tti media may suffer from numerical instabilities and shear artifacts due to the coupling of the p and sv wave modes in the coupled equations on the contrary the decoupled equations for tti media provide a more stable solution due to the separated p and s wave modes because the decoupled equations involve complicated pseudo differential operators in space it is more convenient to apply the pseudo spectral method to these equations however the second order time stepping scheme of the pseudo spectral method may suffer from time stepping errors and instabilities for a large time step we have developed an optimized pseudo differential operator opo that incorporates not only the spectral evaluation of the pseudo differential operator but also a temporal correction term that would effectively compensate the time stepping errors associated with the time wavenumber domain extrapolation of the wave equation the opo was const

anisotropic wave propagation through finite difference grids an algorithm is presented to solve the elastic wave equation by replacing the partial differentials with finite differences it enables wave propagation to be simulated in three dimensions through generally anisotropic and heterogeneous models the space derivatives are calculated using discrete convolution sums while the time derivatives are replaced by a truncated taylor expansion a centered finite difference scheme in cartesian coordinates is used for the space derivatives leading to staggered grids the use of finite difference approximations to the partial derivatives results in a frequency dependent error in the group and phase velocities of waves for anisotropic media the use of staggered grids implies that some of the elements of the stress and strain tensors must be interpolated to calculate the hook sum this interpolation induces an additional error in the wave properties the overall error depends on the precision of the derivative and interpolation operators the anisotropic symmetry system its orientation and the degree of anisotropy the dispersion relation for the homogeneous case was derived for the proposed scheme since we use a general description of convolution sums to describe the finite difference operators the numerical wave properties can be calculated for any space operator and an arbitrary homogeneous elastic model in particular phase and group velocities of the three wave types can be determined in any direction we demonstrate that waves can be modeled accurately even through models with strong anisotropy when the operators are properly designed

multiscale modeling of acoustic wave propagation in 2d media abstractconventional finite difference methods produce accurate solutions to the acoustic and elastic wave equation for many applications but they face significant challenges when material properties vary significantly over distances less than the grid size this challenge is likely to occur in reservoir characterization studies because important reservoir heterogeneity can be present on scales of several meters to ten meters here we describe a new multiscale finite element method for simulating acoustic wave propagation in heterogeneous media that addresses this problem by coupling fine and coarse scale grids the wave equation is solved on a coarse grid but it uses basis functions that are generated from the fine grid and allow the representation of the fine scale variation of the wavefield on the coarser grid time stepping also takes place on the coarse grid providing further speed gains another important property of the method is that the basis functions are only computed once and time savin

identification and estimation of treatment effects with a regression discontinuity design the regression discontinuity rd data design is a quasi experimental design with the defining characteristic that the probability of receiving treatment changes discontinuously as a function of one or more underlying variables this data design arises frequently in economic and other applications but is only infrequently exploited as a source of identifying information in evaluating effects of a treatment in the first application and discussion of the rd method thistlethwaite and campbell 1960 study the effect of student scholarships on career aspirations using the fact that awards are only made if a test score exceeds a threshold more recently van der klaauw 1997 estimates the effect of financial aid offers on students decisions to attend a particular college taking into account administrative rules that set the aid amount partly on the basis of a discontinuous function of the students grade point average and sat score angrist and lavy 1999 estimate the effect of class size on student test scores taking advantage of a rule stipulating that another classroom be added when the average class size exceeds a threshold level finally black 1999 uses an rd approach to estimate parents willingness to pay for higher quality schools by comparing housing prices near geographic school attendance boundaries regression discontinuity methods have potentially broad applicability in economic research because geographic boundaries or rules governing programs often create discontinuities in the treatment assignment mechanism that can be exploited under the method although there have been several discussions and applications of rd methods in the literature important questions still remain concerning sources of identification and ways of estimating treatment effects under minimal parametric restrictions here we show that identifying conditions invoked in previous applications of rd methods are often overly strong and that treatment effects can be nonparametrically identified under an rd design by a weak functional form restriction the restriction is unusual in that it requires imposing continuity assumptions in order to take advantage of the known discontinuity in the treatment assignment mechanism we also propose a way of nonparametrically estimating treatment effects and offer an interpretation of the wald estimator as an rd estimator

moving mesh methods with locally varying time steps the time steps associated with moving mesh methods are proportional to the smallest mesh size in space and as a result they are very small at each time level for some practical problems the physical phenomena develop dynamically singular or nearly singular solutions in fairly localized regions and therefore the smallest time step at each time level occurs only in these localized regions in this work we will develop a local time stepping algorithm for the moving mesh methods the principal idea will be demonstrated by investigating the nonlinear hyperbolic conservation laws numerical experiments are carried out to demonstrate the efficiency and robustness of the proposed methods

a new family of finite difference schemes to solve the heterogeneous acoustic wave equation equivalent staggered grid scheme esg is a new family of schemes based on the finite difference method fdm the method is applied to acoustic wave propagation in variable density media and the results are compared with those from some classic fdm approaches the main feature of this new family is that it is designed to generate results numerically equivalent to those using the standard staggered grid formulations ssg but with the same memory requirements of simple grid schemes hence it results in a reduction of memory usage by 33 in 2d and 50 in 3d problems compared to the memory usage of ssg the first order system of equations in terms of pressure and velocity is not used here instead the formulation is based on applying new central difference operators to the secondorder acoustic wave equation in terms of pressure obtaining the same level of accuracy and stability as the ssg schemes the equivalence between the esg and ssg is mathematically demonstrated and issues concerning the application of seismic sources and the boundary conditions are addressed