Spectral-element simulations of elastic wave propagation in exploration and geotechnical applications
|
|
- Helena Benson
- 5 years ago
- Views:
Transcription
1 Earthq Sci (2014) 27(2): DOI /s RESEARCH PAPER Spectral-element simulations of elastic wave propagation in exploration and geotechnical applications Lin Zheng Qi Zhao Bernd Milkereit Giovanni Grasselli Qinya Liu Received: 18 September 2013 / Accepted: 14 January 2014 / Published online: 22 February 2014 The Seismological Society of China, Institute of Geophysics, China Earthquake Administration and Springer-Verlag Berlin Heidelberg 2014 Abstract We apply the spectral-element method (SEM), a high-order finite-element method (FEM) to simulate seismic wave propagation in complex media for exploration and geotechnical problems. The SEM accurately treats geometrical complexities through its flexible FEM mesh and accurately interpolates wavefields through high-order Lagrange polynomials. It has been a numerical solver used extensively in earthquake seismology. We demonstrate the applicability of SEM for selected 2D exploration and geotechnical velocity models with an open-source SEM software package SPECFEM2D. The first scenario involves a marine survey for a salt dome with the presence of major internal discontinuities, and the second example simulates seismic wave propagation for an open-pit mine with complex surface topography. Wavefield snapshots, synthetic seismograms, and peak particle velocity maps are presented to illustrate the promising use of SEM for industrial problems. Keywords Seismic wave propagation Spectral-element method Exploration seismology L. Zheng Q. Liu (&) Department of Physics, University of Toronto, Toronto, ON, Canada liuqy@physics.utoronto.ca Q. Zhao G. Grasselli Department of Civil Engineering, University of Toronto, Toronto, ON, Canada B. Milkereit Department of Earth Sciences, University of Toronto, Toronto, ON, Canada 1 Introduction Mapping the elastic and anelastic structures of the subsurface, for both global and industrial seismic applications, relies on the ability to solve the seismic wave equations accurately (e.g., Dahlen and Tromp 1998; Aki and Richards 2002). In particular, with the advancing of numerical methods and increasing computational power, full numerical calculations of synthetic seismograms have contributed significantly to the seismic inverse problem in seismological research since the early 80 (e.g., Tarantola 1984a, b; Chaljub et al. 2007; Virieux and Operto 2009; Liu and Gu 2012). The finite-difference (FD) methods are the most widely used among all numerical approaches that model the propagation of seismic waves due to their relatively simple formulation. They employ the strong form of the differential equations of motion and essentially replace derivatives of velocity vectors and stress tensors by their finite difference equivalents between adjacent grid points (e.g., Madariaga 1976; Virieux 1988; Levander 1988; Komatitsch et al. 2005). Other classical gridding methods, such as spectral and pseudo-spectral methods, have also been applied to accurately solve seismic wave propagation in relatively smooth media (e.g., Tessmer et al. 1992; Carcione 1994; Komatitsch et al. 2005). However, these classical methods have inherent limitations for models that exhibit complex geometries, including undulated surface topography and major internal discontinuities, which although could be overcome, are numerically costly (e.g., Robertsson 1996; Ohminato and Chouet 1997; Komatitsch et al. 2005). In comparison, the finite-element method (FEM), another category of numerical techniques in the space-time domain, honors the geometrical model complexity by sophisticated mesh design (Bathe 2008), and has been exploited to predict ground motions for complex media
2 180 Earthq Sci (2014) 27(2): (e.g., Lysmer and Drake 1972; Bao et al. 1998; Bielak et al. 1999; Aagaard et al. 2001). The use of the unstructured meshing technique, which is based upon various types of elements (e.g., triangles, tetrahedra, and hexahedra), provides the FEM with the flexibility to honor strong surface topography and major internal discontinuities, as element boundaries coincide with these interfaces and discontinuities (Komatitsch et al. 2005; Tromp et al. 2008; Casarotti et al. 2008; Peter et al. 2011). Traditionally, lowdegree polynomials are used in the FEM to represent shapes of elements as well as to express functions and integrals on the elements, which may be inadequate to capture large oscillations of wavefield across an element in seismic wave propagation problems. The use of higherorder elements in FEM can increase the accuracy of simulations although at the expense of significantly arising computation cost (Komatitsch et al. 2005). The spectral-element method (SEM) is essential a higherdegree FEM, which combines the boundary composition flexibility of the FEM with the spectral accuracy of the pseudo-spectral method (Komatitsch and Tromp 1999; Komatitsch et al. 2004). Unlike the classic FEM, the use of high-degree polynomials defined over Gauss Lobatto Legendre (GLL) points to discretize functions on elements and as quadrature points for numerical integration, leads to a diagonal mass matrix, which can be easily inverted. It can be also efficiently implemented in parallel (e.g., Komatitsch and Tromp 1999; Komatitsch et al. 2005; Tromp et al. 2008) and makes it possible to solve problems of increasing scales. These favorable features have made SEM an important tool for seismic wave simulations in regional (e.g., Komatitsch and Tromp 1999; Komatitsch et al. 2004; Liu et al. 2004), as well as global (e.g., Komatitsch and Tromp 2002a, b; Chaljub and Valette 2004) earth models. Recently, some efforts have also been made to employ SEM in industrial-scale problems (e.g., Zhu et al. 2009; Luo 2012; Gharti et al. 2012a, b). In this paper, we take advantage of an existing SEM software package SPEC- FEM2D ( to study two cases of exploration and geotechnical problems, including a marine survey for a salt dome, and wave simulations for an open-pit mine scenario. We analyze wavefield snapshots, synthetic seismograms, and peak particle velocity (PPV) maps from these simulations to demonstrate the applicability of SEM in the forward modeling of wave propagation in complex media for industrial problems. 2 Brief introduction to the SEM As the main focus of this paper is on the potential applications of SEM to specific types of industrial scenarios, we only briefly summarize the SEM technique below, with an emphasis on mesh design. More detailed comparisons to other numerical methods, such as FD and low-order FEM, can be found in Komatitsch and Tromp (1999), Komatitsch et al. (2004), Chaljub et al. (2007), and Peter et al. (2011). 2.1 Governing equations Given an Earth model with volume X and a free surface ox, the governing equations for elastodynamic problems can be generally written in the differential form qo 2 t s ¼r Tþf; ð1þ subject to the free-surface boundary condition ^n T ¼ 0 on ox; ð2þ where s(x, t) denotes the displacement motion at position x, q(x) the distribution of mass density, and f(x, t) the external body force. The stress tensor T is related to the displacement gradient rs through the constitutive relation T ¼ c : rs; ð3þ where c is a fourth-order elastic tensor. In addition, the displacement s and traction ^n T should be continuous across solid-solid interfaces. If both solid and fluid regions are present in the earth model, the SEM formulates the elastodynamic problem in terms of velocity potential in the fluid and displacement in solid regions, resulting in separate governing equations that are coupled at fluid-solid boundaries through continuity of traction and the normal component of velocity (Komatitsch et al. 2000a; Komatitsch and Tromp 2002a, b; Chaljub and Valette 2004; Tromp et al. 2008). The integral or weak form of the elastodynamic equation, frequently used in both FEM and SEM, is obtained by dotting the differential equation (1) the strong form for elastostatic problems with an arbitrary test vector function w(x), and integrating over the entire model (Komatitsch and Tromp 1999; Komatitsch et al. 2005; Tromp et al. 2008) Z Z Z qw o 2 t sðdþ3 x ¼ rw : Td 3 x þ f wd 3 x: ð4þ X X Since Eq. (4) is valid for any test vector w, the weak and strong forms are mathematically equivalent. Note that the free-traction surface boundary condition (2) is satisfied naturally through integration by parts, so that the incorporation of surface topography is straightforward. Therefore, surface waves can be simulated more accurately by SEM than those methods based upon the strong formulation (Komatitsch and Tromp 1999). For regional problems, the differential equation (1) needs to be solved subject to absorbing boundary condition in which seismic waves travelling out of the volume X are absorbed on fictitious X
3 Earthq Sci (2014) 27(2): boundaries and do not re-enter the domain (Komatitsch and Tromp 1999; Tromp et al. 2008). For the weak formulation, the absorbing boundary condition can be naturally imposed on Eq. (4) in terms of traction over the fictitious boundaries. 2.2 The spectral-element discretization Similar to FEM, the first step of discretizing the governing equations is to subdivide the model volume X into a number of non-overlapping elements (e.g., Komatitsch and Tromp 1999; Komatitsch et al. 2005; Tromp et al. 2008; Casarotti et al. 2008). Although a variety of elements such as tetrahedra, hexahedra, and prisms may be used in the FEM, the SEM is generally restricted to use meshes of unstructured hexahedral elements (deformed cube) in 3D case or quadrilateral elements (deformed square) in 2D case. Each hexahedral element can be mapped to a reference cube which is defined in terms of natural coordinates n ¼ ðn; g; fþ; 1 n; g; f 1 using the invertible transformation. xðnþ ¼ Xn a x a N a ðnþ; ð5þ a¼1 where x a ðn a ; g a ; f a Þ; a ¼ 1;...; n a denote n a control points within a hexahedron and N a ðnþ corresponding shape functions. The hexahedral shape functions N a are triple products of one-dimensional Lagrange polynomials of degree 1 or 2. Generally, the n? 1 Lagrange polynomials of degree n are defined in terms of n? 1 control points 1 n a 1; a ¼ 0;...; n, by a ¼ ðn n 0Þ ðn n a 1 Þðn n aþ1 Þ ðn n n Þ ðn a n 0 Þ ðn a n a 1 Þðn a n aþ1 Þ ðn a n n Þ : ð6þ To ensure the mapping from the element to the reference cube is unique and invertible, the transformation (5) should have nonzero Jacobian. Similarly, in 2D, each quadrilateral element is mapped to a reference square. In the SEM, a function (or field) f is interpolated over a hexahedral element using the Lagrange polynomials of degree n defined over the Gauss-Lobatto-Legendre (GLL) points as (Komatitsch and Tromp 1999; Komatitsch et al. 2005; Tromp et al. 2008) f ðxðn; g; fþþ ¼ Xn a;b;c¼0 f abc a ðnþ b ðþ g c ðþ; f ð7þ where f abc ¼ f x n a ; g b ; f c denotes the functional value at the GLL point x n a ; g b ; f c and a,b,c the Lagrange polynomials. The n? 1 GLL points in each dimension are determined by the roots of the equation 1 n 2 P 0 n ðnþ ¼ 0, where P 0 n denotes the first derivative of the Legendre polynomial of degree n. SEM simulations typically use polynomials of degree between 4 to 10 to represent functions on the elements, and degree 4 or 5 provides the best trade-off between accuracy and computational cost (Seriani and Priolo 1994; Tromp et al. 2008; Casarotti et al. 2008). In comparison, the classic FEM uses lower-degree polynomials (1 or 2 in general) as basis functions for both the representation of functions and the shape of elements. To evaluate the integrals at the element level, the SEM utilizes the Gauss Labotto Legendre integration rule in conjunction with GLL interpolation points, instead of the Gauss quadrature frequently used in the classical FEM (Komatitsch and Tromp 1999; Komatitsch et al. 2005; Tromp et al. 2008; Casarotti et al. 2008). Although slightly less accurate, the use of the GLL quadrature significantly reduces the complexity and cost of the algorithm, and the mass matrix resulted from the evaluation of the left-hand term in the weak form (4) becomes diagonal (Komatitsch and Tromp 1999). Thus, the SEM can be efficiently run on parallel computers (Komatitsch and Tromp 1999). 2.3 Mesh design for the SEM Generating a high-quality unstructured mesh appropriate for the model of interest, especially in 3D, is the first as well as one of the critical steps of both the FEM and SEM wave modeling processes (Casarotti et al. 2008). The fact that efficient 3D SEM simulations generally require unstructured all-hexahedral mesh significantly complicates the meshing process (Casarotti et al. 2008). As mentioned previously, an unstructured hexahedral mesh not only provides flexibility for the incorporation of topographic variations and major internal discontinuities (Komatitsch et al. 2005), but also leads to a diagonal mass matrix in conjunction with the use of the GLL integration points and quadrature. The effort of 2D/3D hexahedral mesh design is also justified by the fact that once a proper mesh is built for a region, wave simulations can be run repeated without any alteration to the mesh. In order to ensure accuracy and computational stability, the unstructured mesh should be constructed so that the Jacobian of the transformation (5) varies smoothly throughout the model, and the distortion of elements is acceptable (Komatitsch et al. 2005; Casarotti et al. 2008). Generally, the equiangle skewness of the deformed hexahedral elements should be within 0.75 (0.8 is acceptable in some cases) to provide acceptable simulation accuracy (Komatitsch et al. 2005; Casarotti et al. 2008). To resolve a shortest period T 0 of interest, the choice of the element size Dh needs to be subject to the constraint
4 182 Dh\ Earthq Sci (2014) 27(2): nþ1 tmin T0 ; f ðnþ ð8þ where tmin denotes the minimum seismic wave speed in the model, and f(n) is an empirical function related to the degree n of polynomials used to interpolate functions on the elements (Casarotti et al. 2008). For example, if n is chosen as 4, the empirical function f(n) is equal to 5, ensuring that each wavelength k contains at least 5 GLL points (Casarotti et al. 2008). The incremental time stepping duration is in turn constrained by the element size and the maximum seismic wave speed tmax in the element by Dt\C Dh ; tmax (a) ð9þ where C is the Courant constant and can be chosen between 0.3 and 0.4 in the case of a deformed and heterogeneous mesh (Komatitsch et al. 2005; Casarotti et al. 2008). Moreover, the major discontinuities in wave speeds within the model need to be honored by the mesh to avoid numerical diffractions by a staircase sampling of the interfaces (e.g., Zahradnı k et al. 1993; Komatitsch et al. 2004; Casarotti et al. 2008). 3 SEM simulation examples In this paper, a fully-developed 2D SEM community software package SPECFEM2D (downloaded from geodynamics.org) is applied to simulate seismic wave propagation in complex models that include both solid and fluid media, undulation of major discontinuities and strong lateral heterogeneities. In SPECFEM2D, polynomials of degree 4 associated with 5 GLL control points are used as basis functions to interpolate functions on each element. The SPECFEM2D package is an open-source MPI-based package that can be readily used in parallel. Another advantage of the SPECFEM2D package is that it supports not only the built-in mesher which is suitable for models with simple and regular geometries, but also mesh built by external quadrilateral/hexahedral meshing softwares such as CUBIT ( cubit.sandia.gov), which provide more flexibility to the creation of high-quality mesh for complex media. All the quadrilateral meshes used in this paper are created with CUBIT. The SEM simulation examples presented are restricted to isotropic elastic and/or acoustic media. However, features such as full anisotropy, poroelasticity and attenuation are also implemented in the SPECFEM2D package and may be incorporated for future studies. 3.1 Simulation for exploration application In this section, two simulations for a marine survey scenario are presented and analyzed based upon wavefield (b) Fig. 1 Snapshots of velocity vectors at time t = 0.75 s on top of the unstructured quadrilateral meshes for the wave propagation simulations in Sect. 1. The size of modeling domains is 4, ,400 m for both models. The ocean layer with an average thickness of 600 m, the surrounding sediments beneath the ocean floor, and the salt dome are shown in red, blue, and yellow. The arrow indicates the position of the pressure source placed at 2 m below the ocean surface, diamonds indicate the position of the 101 receivers used to record seismograms as presented in Fig. 2, and the white triangle in the sediments indicates the position of the receiver used to record seismic wave propagation in the elastic regions as shown in Fig. 3. One can clearly observe the direct (a) and sea-floor reflected (b) P-waves in the ocean layer; the transmitted P (c) and P-to-S (d) converted waves through sea floor, as well as the reflected P (e) and S (f) waves off the top of the salt dome in the sediments. Phase (g) is the refracted P waves along the sea floor snapshots and seismograms, to illustrate the capability of the SEM, along with the use of unstructured quadrilateral mesh, in accommodating fluid layers (ocean), fluid solid interface topography (bathymetry), and abrupt changes in the seismic velocity model. Both simulations utilize a similar geologic models representing a high-velocity salt dome beneath ocean sediments in a computation domain of 4, ,400 m (Fig. 1), and the only difference is that one includes a flat sea floor and the other one is simulated with an uneven sea floor (Fig. 1). In both models, The first layer (shown in red in Fig. 1) is a water layer with a thickness of approximately 600 m, and the sedimentary
5 Earthq Sci (2014) 27(2): Table 1 Elastic properties of models used in wave propagation simulations in Sect. 1 Density, q (kg/m 3 ) P-wave velocity, V p (m/s) Ocean layer 1,000 1,450 0 Surrounding 2,200 2,500 1,600 sediment Salt dome 2,300 5,200 3,000 S-wave velocity, V s (m/s) layer (shown in blue) extends to *2,400 m beneath the sea floor. A salt dome (shown in yellow) is placed in the sediments with significantly larger P- and S-velocity than the surrounding. The elastic properties of the different materials are listed in Table 1. Taking into account the velocity variations, element sizes need to vary for different blocks to maintain on average a consistent number of grid (e.g., C5) points per S-wavelength for the simulations. The average element sizes for the fluid layer, sedimentary layer, and salt dome, are 12, 18, and 20 m. The CUBIT software automatically generates conforming unstructured quadrilateral meshes that accommodate the transition in element size between different materials (Fig. 1). Absorbing boundary conditions (Stacey 1988) are applied to the fictitious side and bottom boundaries of the model domains to prevent reflections. To investigate the impact of sea floor topography on seismic wave propagation, the same source mechanism and sourcereceiver geometry are used in both simulations. A Ricker pressure source (indicated by vertical arrows at [2,000, -2 m] in Fig. 1) with a dominant frequency of 20 Hz is injected 2 m below the ocean surface, to mimic pulses generated by airguns in marine surveys. A total number of 101 receivers are also evenly placed along a horizontal line 2 meters below the ocean surface with a source offset varying between -1,500 and 1,500 m to synthesize the seismic recordings of a string of hydrophones. Another receiver (at [2200, -800 m]) located in the surrounding sediment is used to record seismic waves motions in the elastic regions. Time steps are estimated for each block based on material wavespeed and mesh size according to Eq. (9), and the final time step employed by the simulation, Dt ¼ s; is chosen to be the smallest of all three estimates. A total duration of 1.5 s is simulated. (a) (b) Fig. 2 Synthetic seismograms of the horizontal (top) and vertical (bottom) component velocity field recorded by the 101 receives denoted by the diamonds in Fig. 1 for the model with a flat ocean floor (a) and varying sea floor topography (b). Data traces are scaled by the maximum vertical velocity values in each scenario. The amplitudes of horizontal records are blown-up by a factor of 10 to better show the much weaker phases. The strong velocity contrast between the water and sedimentary layers results in the large amplitude of the reflected P wave (phase b). Phase (e 0 ) and (f 0 ) are the transmitted P waves in the ocean layer of phase (e) and (f), which are the P and S reflected waves off the top of salt dome
6 184 Earthq Sci (2014) 27(2): (a) (b) Fig. 3 Synthetic seismograms of velocity field recorded by the receiver in the sediment as denoted by the white triangle in Fig. 1. The dashed and solid lines correspond to the synthetics for model with and without sea floor topography. Both traces are scaled by a constant to allow for their proper display. The impact of the sea floor topography on seismic wave propagation is clearly demonstrated by the difference in arrival times of various phases, such as the transmitted P waves (phase c), the P-to-S converted wave at sea floor (phase d), the reflected P (phase e), and S (phase f) waves off the top of salt dome The source receiver geometry as well as snapshots of the velocity vector field at time t = 0.75 s for both models are shown in Fig. 1. The corresponding horizontal- and vertical-component velocity seismogram profiles recorded by the horizontal receiver line are illustrated in Fig. 2. The two-component velocity seismograms recorded by the specific receiver located in the sediment for the models with flat and varying sea-floor topography are shown as the solid and dashed lines in Fig. 3. Several phases can be identified in all three figures, including (a) direct P waves, (b) reflected P waves from seabed, (c) transmitted P waves through seabed, (d) P-to-S converted waves across seabed, (e) reflected P waves of (c) off the top of salt dome, (e 0 ) transmitted waves of (e) through seabed, (f) reflected S waves of (c) off the top of salt dome, (f 0 ) transmitted waves of (f) converted to P waves in the ocean layer, and (g) refracted P waves along seabed. The presence of sea floor topography clearly affects the traveltime and amplitude of various phases, and is nicely captured by the SEM simulations. Both the snapshots and synthetic seismograms illustrate the ability of SEM method in simulating wave propagation in models with both fluid and solid regions as well as sharp velocity contrasts. 3.2 Simulation for geotechnical application In this section, we simulate wave propagation in three homogeneous 2D models with increasing free-surface complexity for an open-pit mine scenario (Fig. 4). Snapshots and synthetic seismograms, as well as the PPV maps for first breaks, are produced to illustrate the ability of the SEM in simulating seismic wave propagation in models that include interfaces with large topographic variations. The first model consists of a flat free-surface and is referred to as the flat-surface model; the second model is referred to as bowl-shape model because of the bowlshaped surface depression; and the last one is inspired by section views of realistic open-pit mine models with staircase free-surface and referred to as the staircase model. In the staircase model, the bench width and height are 15 and m, respectively, and the face angle is 56. The geotechnical ramp width is 40 m and the depth of the pit is m. The simulation domains for all models are approximately equal to 2, , m, and the width of both the bowl-shape and the staircase topographical regions are 1,616 m. To compare the effect of surface topography, the same seismic source mechanism and source-receiver layout, (a) (b) (c) Fig. 4 (a) Geometry of the models described in Sect. 2. The size of all models is approximately equal to 2, , m. The maximum widths and depths of both the bowl-shaped and staircase topographical regions are 1,616 and m, respectively. A Ricker point force (denoted by the black star) with a dominant frequency of 100 Hz is excited vertically at (1108, m). The recordings of 82 ground stations, located unevenly up to 1,010 m offset away from the source, are plotted in Fig. 6. (b) Bench geometry of the staircase model. (c) Unstructured quadrilateral mesh for the portion of the staircase model shown in (b)
7 Earthq Sci (2014) 27(2): =0.05 s 185 =0.25 s =0.15 s (a) Flat-surface model =0.05 s =0.15 s =0.25 s (b) Bowl-shape model =0.05 s =0.15 s =0.25 s (c) Staircase model Fig. 5 Snapshots of vertical displacement fields at t = 0.05, 0.15, and 0.25 s for simulations in Sect. 2. Grey color indicates the modeling domains. Red and blue colors represent positive and negative wavefield values (a) Flat-surface model (b) Bowl-shape model (c) Staircase model Fig. 6 Synthetic seismograms of the vertical velocity fields (black) recorded by the 82 receivers located along the ground surface and derived peak particle velocities (PPVs) for first breaks (red). The lateral variations of PPVs are due to the reflected and scattered waves caused by strong surface topography. The reflected waves (indicated by r) from the left and right absorbing boundaries can be observed, but their amplitudes are generally small compared with other phases as well as elastic properties are used for all three simulations. The seismic source is given vertically at m depth (i.e., 250 m below the bottom of the pit) as a Ricker point source with a dominant frequency of 100 Hz. In order to capture effects of topography on seismic wave propagation and energy radiation, 82 receivers are unevenly located along each ground surface up to a horizontal offset of 1,010 m to the left and right of the source. The homogeneous P- and S-wave velocities are 5,200 and 3,200 m/s, and mass density is 2,800 kg/m3. According to Eq. (8), to resolve waves with a dominant frequency of 100 Hz, the element size should be smaller than 32 m when 5 GLL points are used. Nonetheless, an average element size of 10 m is used to discretize the models so that each bench of the
8 186 Earthq Sci (2014) 27(2): staircase model can be subdivided into two elements (Fig. 4c). The same time stepping Dt ¼ 5: is used for all the simulations, and a total duration of 0.5 s is simulated. Stacey absorbing boundary condition are applied to all the simulations to prevent reflections from the artificial boundaries. The displacement snapshot of seismic wave propagation in the three models are shown in Fig. 5 at t = 0.05 s, t = 0.15 s, and t = 0.25 s. For the homogeneous models, as wavefronts expand from the source, P and S waves clearly separate and show different radiation patterns (t = 0.15 s in Fig. 5a). Both waves are later reflected by the flat surface and travel downward (t = 0.25 s), including P-to-P, P-to-S, S-to-S, and S-to-P reflected and converted waves. In the model with bowl-shaped surface, P and S surface reflections occur much earlier (t = 0.15 s) due to the reduction of topography above the source. In comparison, in the staircase model, strong scattering effects is observed in addition due to abrupt variations of the ground surface. Figure 6 illustrates the vertical-component velocity seismograms recorded by the ground stations (shown in black), along with the derived PPVs for first breaks (in red). The scattered waves off the abrupt changes in surface topography in the staircase model are clearly shown in the synthetic seismograms (Fig. 6c), thanks to the power of SEM in accommodating strong free-surface topographic variations. Some artificial boundary reflected waves can be seen in the seismograms [denoted as (r) in Fig. 6], as Stacey condition works less effectively in this case probably due to the proximity of left and right fictitious boundaries to the surface regions with strong topography (i.e., edge of the model is only 98 m away from the nearest receiver), although amplitudes of these artificial reflections are much smaller than those of the direct P and S waves. The models used in these simulations are fully elastic with no attenuation. Wave amplitudes generally decrease with increasing distance from the source due to geometric spreading effect, as shown by the PPV plots, where smaller PPV value is generally observed for station further away from the source. However, due to the impact of large surface topographic undulation, lateral variations of PPVs can be observed, especially in the staircase simulation (Fig. 6c). These PPV values provide critical input to the analysis of slope stability for mining operations (Dowding 1985). 4 Conclusions In this paper, we presented two simulation examples for exploration and geotechnical problems based upon the SEM. In Sect. 1, complex models that include the fluid layer (ocean), fluid solid interface (sea floor) topography, and the asymmetric salt dome are easily accommodated by the SEM because of the use of the unstructured mesh generated by CUBIT as shown in Fig. 1. The snapshots of the velocity vector field and synthetic seismograms along with the distinct phases of the subsalt reflections illustrate that the SEM can accurately model seismic wave propagation in media with both fluid and solid regions as well as interface variations. The influence of the sea floor topography is demonstrated by differences in arrival times and amplitudes of various phases (see Fig. 3). For geotechnical applications, Sect. 2 presents and analyzes SEM simulation results for an open-pit mine scenario, including highly irregular surfaces representing benches. With the simulations in the flat-surface, bowlshape and staircase models, the capability of the SEM in handling strong surface topography is again proven by the displacement snapshots (Fig. 5) and synthetic seismograms (Fig. 6) recorded by ground stations. Therefore, the SEM can be naturally applied to exploration and geotechnical applications as its weak formulation and flexible mesh design make it possible to deal with model complexities. It shows great potential in simulating seismic wave propagation more accurately and effectively for some industrial problems compared to other numerical methods, such as the FD method, the pseudo-spectral element method, and the traditional low-order FEM. Acknowledgments This research was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) and Center for Excellence in Mining Innovations (CEMI, through SUMIT project). Computations for this study were performed on hardwares purchased through the combined funding of Canada Foundation for Innovation (CFI), Ontario Research Fund (ORF) and University of Toronto Startup Fund. We thank the developers of the SPECFEM2D package for their continued community support. We would also like to thank the two anonymous reviewers for their helpful comments and suggestions. References Aagaard BT, Hall JF, Heaton TH (2001) Characterization of nearsource ground motions with earthquake simulations. Earthq Spectra 17(2): Aki K, Richards PG (2002) Quantitative seismology: theory and methods, 2nd edn. University Science Books, Sausalito Bao H, Bielak J, Ghattas O, Kallivokas LF, O Hallaron DR, Shewchuk JR, Xu J (1998) Large-scale simulation of elastic wave propagation in heterogeneous media on parallel computers. Comput Methods Appl Mech Eng 152(1): Bathe K-J (2008) Finite element method. In: Wiley encyclopedia of computer science and engineering. Wiley, New York, pp Bielak J, Xu J, Ghattas O (1999) Earthquake ground motion and structural response in alluvial valleys. J Geotech Geoenviron Eng 125(5):
9 Earthq Sci (2014) 27(2): Carcione JM (1994) The wave equation in generalized coordinates. Geophysics 59(12): Casarotti E, Stupazzini M, Lee S, Komatitsch D, Piersanti A, Tromp J (2008) CUBIT and seismic wave propagation based upon the spectral-element method: an advanced unstructuredmesher for complex 3D geological media. In: Proceedings of the 16th international meshing roundtable, Seattle, vol 5B.4, pp Chaljub E, Valette B (2004) Spectral element modelling of threedimensional wave propagation in a sel-gravitating Earth with an arbitrarily stratified outer core. Geophys J Int 158(1): Chaljub E, Komatitsch D, Vilotte J-P, Capdeville Y, Valette B, Festa G (2007) Spectral-element analysis in seismology. Adv Geophys 48: Dahlen FA, Tromp J (1998) Theoretical global seismology. Princeton University Press, Princeton Dowding CH (1985) Blast vibration monitoring and control. Prentice Hall, Saddle River Gharti HN, Komatitsch D, Oye V, Martin R, Tromp J (2012) Application of an elastoplastic spectral element method to 3D slope stability analysis. Int J Numer Methods Eng 91(1):1 26 Gharti HN, Oye V, Komatitsch D, Tromp J (2012) Simulation of multistage excavation based on a 3D spectral-element method. Comput Struct :54 69 Komatitsch D, Tromp J (1999) Introduction to the spectral element method for three-dimensional seismic wave propagation. Geophys J Int 139(3): Komatitsch D, Tromp J (2002) Spectral-element simulations of global seismic wave propagation-i. Validation. Geophys J Int 149(2): Komatitsch D, Tromp J (2002) Spectral-element simulations of global seismic wave propagation-ii. 3-D models, oceans, rotation, and self-gravitation. Geophys J Int 150(1): Komatitsch D, Barnes C, Tromp J (2000a) Wave propagation near a fluid-solid interface: a spectral-element approach. Geophysics 65(2): Komatitsch D, Liu Q, Tromp J, Süss P, Stidham C, Shaw JH (2004) Simulations of ground motion in the Los Angeles basin based upon the spectral-element method. Bull Seismol Soc Am 94(1): Komatitsch D, Tsuboi S, Tromp J (2005) The spectral-element method in seismology. In: Levander A, Nolet G (eds) Seismic earth: array analysis of broadband seismograms: geophysical monograph, vol 157. AGU, Washington DC, pp Levander AR (1988) Fourth-order finite-difference p-sv seismograms. Gophysics 53(11): Liu Q, Gu YJ (2012) Seismic imaging: from classical to adjoint tomography. Tectonophysics :31 66 Liu Q, Polet J, Komatitsch D, Tromp J (2004) Spectral-element moment tensor inversions for earthquakes in southern Califonia. Bull Seismol Soc Am 94(5): Luo Y (2012) Seismic imaging and inversion based on spectralelement and adjoint methods. Ph.D. Dissertation, Princeton University Lysmer J, Drake LA (1972) A finite element method for seismology. Methods Comput Phys 11: Madariaga R (1976) Dynamics of an expanding circular fault. Bull Seismol Soc Am 65(3): Ohminato T, Chouet BA (1997) A free-surface boundary condition for including 3d topography in the finite-difference method. Bull Seismol Soc Am 87: Peter D, Komatitsch D, Luo Y, Martin R, Le Goff N, Casarotti E, Loher PL, Magnoni F, Liu Q, Blitz C, Nissen-Meyer T, Basini P, Tromp J (2011) Forward and adjoint simulations of seismic wave propagation on fully unstructured hexahedral meshes. Geophys J Int 186(2): Robertsson JO (1996) A numerical free-surface condition for elastic/ viscoelastic finite-difference modelling in the presence of topography. Geophysics 61(6): Seriani G, Priolo E (1994) Spectral element method for acoustic wave simulation in heterogeneous media. Finite Elem Anal Des 16(3): Stacey R (1988) Improved transparent boundary formulations for the elastic wave equation. Bull Seismol Soc Am 78(6): Tarantola A (1984) Inversion of seismic reflection data in the acoustic approximation. Geophysics 49(8): Tarantola A (1984) Linearized inversion of seismic reflection data. Geophys Prospect 32(6): Tessmer E, Kessler D, Kosloff D, Behle A (1992) Multi-domain Chebyshev Fourier method for the solution of the equations of motion of dynamic elasticity. J Comput Phys 100(2): Tromp J, Komatitsch J, Liu Q (2008) Spectral-element and adjoint methods in seismology. Commun Comput Phys 3(1):1 32 Virieux J (1988) P-SV wave propagation in heterogeneous media: velocity stress finite-difference method. Geophysics 51(4): Virieux J, Operto S (2009) An overview of full-waveform inversion in exploration geophysics. Geoophysics 74(6) WCC1 WCC26 Zahradník J, Maczo P, Hron F (1993) Testing four elastic finitedifference schemes for behavior at discontinuities Bull Seismol Soc Am 83(1): Zhu H, Luo Y, Nissen-Meyer T, Morency C, Tromp J (2009) Elastic imaging and time-lapse migration based on adjoint methods. Geophysics 74(6):WCA167 WCA177
Spectral-element Simulations of Elastic Wave Propagation in Exploration and Geotechnical Applications
Spectral-element Simulations of Elastic Wave Propagation in Exploration and Geotechnical Applications Lin Zheng*, Qi Zhao*, Qinya Liu, Bernd Milkereit, Giovanni Grasselli, University of Toronto, ON, Canada
More informationCO2 sequestration crosswell monitoring based upon spectral-element and adjoint methods
CO2 sequestration crosswell monitoring based upon spectral-element and adjoint methods Christina Morency Department of Geosciences, Princeton University Collaborators: Jeroen Tromp & Yang Luo Computational
More informationSpectral-Element and Adjoint Methods in Seismology
Spectral-Element and Adjoint Methods in Seismology Jeroen Tromp Min Chen, Vala Hjorleifsdottir, Dimitri Komatitsch, Swami Krishnan, Qinya Liu, Alessia Maggi, Anne Sieminski, Carl Tape & Ying Zhou Introduction
More informationMulticomponent land data pre-processing for FWI: a benchmark dataset
Multicomponent land data pre-processing for FWI: a benchmark dataset Raul Cova, Bernie K. Law and Kris Innanen CRWES/University of Calgary Summary Successful full-waveform inversion (FWI) studies using
More informationNUMERICAL SIMULATION OF IRREGULAR SURFACE ACOUSTIC WAVE EQUATION BASED ON ORTHOGONAL BODY-FITTED GRIDS
- 465 - NUMERICAL SIMULATION OF IRREGULAR SURFACE ACOUSTIC WAVE EQUATION BASED ON ORTHOGONAL BODY-FITTED GRIDS LIU, Z. Q. SUN, J. G. * SUN, H. * LIU, M. C. * GAO, Z. H. College for Geoexploration Science
More informationStrategies for elastic full waveform inversion Espen Birger Raknes and Børge Arntsen, Norwegian University of Science and Technology
Strategies for elastic full waveform inversion Espen Birger Raknes and Børge Arntsen, Norwegian University of Science and Technology SUMMARY Ocean-bottom cables (OBC) have become common in reservoir monitoring
More informationNSE 1.7. SEG/Houston 2005 Annual Meeting 1057
Finite-difference Modeling of High-frequency Rayleigh waves Yixian Xu, China University of Geosciences, Wuhan, 430074, China; Jianghai Xia, and Richard D. Miller, Kansas Geological Survey, The University
More informationSEG/San Antonio 2007 Annual Meeting
Yaofeng He* and Ru-Shan Wu Modeling and Imaging Laboratory, Institute of Geophysics and Planetary Physics, University of California, Santa Cruz, CA, 95064, USA Summary A one-way and one-return boundary
More informationSeismic modelling with the reflectivity method
Yongwang Ma, Luiz Loures, and Gary F. Margrave ABSTRACT Seismic modelling with the reflectivity method Numerical seismic modelling is a powerful tool in seismic imaging, interpretation and inversion. Wave
More informationFinite-difference elastic modelling below a structured free surface
FD modelling below a structured free surface Finite-difference elastic modelling below a structured free surface Peter M. Manning ABSTRACT This paper shows experiments using a unique method of implementing
More informationSponge boundary condition for frequency-domain modeling
GEOPHYSIS, VOL. 60, NO. 6 (NOVEMBER-DEEMBER 1995); P. 1870-1874, 6 FIGS. Sponge boundary condition for frequency-domain modeling hangsoo Shin ABSTRAT Several techniques have been developed to get rid of
More informationMIGRATION BY EXTRAPOLATION OF TIME- DEPENDENT BOUNDARY VALUES*
Geophysical Prospecting 31,413-420, 1983. MIGRATION BY EXTRAPOLATION OF TIME- DEPENDENT BOUNDARY VALUES* G.A. McMECHAN** ABSTRACT MCMECHAN, G.A. 1983, Migration by Extrapolation of Time-Dependent Boundary
More informationSurface Wave Propagation Analysis of a Ground with Three-Dimensional Irregularities
Proceedings of the 8th International Conference on Structural Dynamics, EURODYN Leuven, Belgium, 4-6 July G. De Roeck, G. Degrande, G. Lombaert, G. Müller (eds.) ISBN 978-9-76-9-4 574 Surface Wave Propagation
More informationSummary. Introduction
Dmitry Alexandrov, Saint Petersburg State University; Andrey Bakulin, EXPEC Advanced Research Center, Saudi Aramco; Pierre Leger, Saudi Aramco; Boris Kashtan, Saint Petersburg State University Summary
More informationNumerical modelling of seismic waves using imageprocessing
Southern Cross University epublications@scu 23rd Australasian Conference on the Mechanics of Structures and Materials 2014 Numerical modelling of seismic waves using imageprocessing and quadtree meshes
More informationSEG/San Antonio 2007 Annual Meeting
Imaging steep salt flanks by super-wide angle one-way method Xiaofeng Jia* and Ru-Shan Wu, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary The super-wide angle one-way
More informationWe N Converted-phase Seismic Imaging - Amplitudebalancing Source-independent Imaging Conditions
We N106 02 Converted-phase Seismic Imaging - Amplitudebalancing -independent Imaging Conditions A.H. Shabelansky* (Massachusetts Institute of Technology), A.E. Malcolm (Memorial University of Newfoundland)
More informationTravel-Time Sensitivity Kernels in Long-Range Propagation
Travel-Time Sensitivity Kernels in Long-Range Propagation Emmanuel Skarsoulis Foundation for Research and Technology Hellas Institute of Applied and Computational Mathematics P.O. Box 1385, GR-71110 Heraklion,
More informationH003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography
H003 Deriving 3D Q Models from Surface Seismic Data Using Attenuated Traveltime Tomography M. Cavalca* (Schlumberger - Westerngeco) & R.P. Fletcher (Schlumberger - Westerngeco) SUMMARY Estimation of the
More informationAVO Analysis with Multi-Offset VSP Data
AVO Analysis with Multi-Offset VSP Data Z. Li*, W. Qian, B. Milkereit and E. Adam University of Toronto, Dept. of Physics, Toronto, ON, M5S 2J8 zli@physics.utoronto.ca T. Bohlen Kiel University, Geosciences,
More informationImprovements in time domain FWI and its applications Kwangjin Yoon*, Sang Suh, James Cai and Bin Wang, TGS
Downloaded 0/7/13 to 05.196.179.38. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/ Improvements in time domain FWI and its applications Kwangjin Yoon*,
More informationFull waveform inversion of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC)
of physical model data Jian Cai*, Jie Zhang, University of Science and Technology of China (USTC) Summary (FWI) is a promising technology for next generation model building. However, it faces with many
More informationMigration from a non-flat datum via reverse-time extrapolation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 50 Migration from a non-flat datum via reverse-time extrapolation Gopal Palacharla 1 ABSTRACT Land surveys usually have elevation changes,
More information(f) (a) Time (s)
A spectral element method or wave propagation simulation near a uid-solid interace Dimitri Komatitsch (1), Christophe Barnes (2), Jeroen Tromp (1) (1) Harvard University, Mailcode 252-21, Caltech, Pasadena,
More informationReflectivity modeling for stratified anelastic media
Reflectivity modeling for stratified anelastic media Peng Cheng and Gary F. Margrave Reflectivity modeling ABSTRACT The reflectivity method is widely used for the computation of synthetic seismograms for
More informationSUMMARY NOISE REDUCTION BY SPATIAL SLOPE ESTIMATION
: steerable and non-linear filter approach Abdulaziz M. Almuhaidib and M. Nafi Toksöz, Earth Resources Laboratory, MIT SUMMARY We present an approach based on local-slope estimation for the separation
More informationGEOPHYS 242: Near Surface Geophysical Imaging. Class 8: Joint Geophysical Inversions Wed, April 20, 2011
GEOPHYS 4: Near Surface Geophysical Imaging Class 8: Joint Geophysical Inversions Wed, April, 11 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical,
More informationAn imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults
Earth Planets Space, 52, 3 11, 2000 An imaging technique for subsurface faults using Teleseismic-Wave Records II Improvement in the detectability of subsurface faults Takumi Murakoshi 1, Hiroshi Takenaka
More informationModeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics
Modeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics Anthony L. Bonomo *1 and Marcia J. Isakson 1 1 Applied
More informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) doi: 10.1111/j.1365-246X.2012.05676.x Numerical simulation of seismic wave propagation produced by earthquake by using a particle method Junichi
More informationTh C 02 Model-Based Surface Wave Analysis and Attenuation
Th C 02 Model-Based Surface Wave Analysis and Attenuation J. Bai* (Paradigm Geophysical), O. Yilmaz (Paradigm Geophysical) Summary Surface waves can significantly degrade overall data quality in seismic
More informationSEG/New Orleans 2006 Annual Meeting
Accuracy improvement for super-wide angle one-way waves by wavefront reconstruction Ru-Shan Wu* and Xiaofeng Jia, Modeling and Imaging Laboratory, IGPP, University of California, Santa Cruz Summary To
More informationDownloaded 09/16/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Time-domain incomplete Gauss-Newton full-waveform inversion of Gulf of Mexico data Abdullah AlTheyab*, Xin Wang, Gerard T. Schuster, King Abdullah University of Science and Technology Downloaded 9// to
More informationu = v is the Laplacian and represents the sum of the second order derivatives of the wavefield spatially.
Seismic Elastic Modelling Satinder Chopra, Arcis Corporation, Calgary, Canada Summary Seismic modeling experiments were performed to understand the seismic responses associated with a variety of subsurface
More informationOverview of the High-Order ADER-DG Method for Numerical Seismology
CIG/SPICE/IRIS/USAF WORKSHOP JACKSON, NH October 8-11, 2007 Overview of the High-Order ADER-DG Method for Numerical Seismology 1, Michael Dumbser2, Josep de la Puente1, Verena Hermann1, Cristobal Castro1
More informationProgress Report on: Interferometric Interpolation of 3D SSP Data
Progress Report on: Interferometric Interpolation of 3D SSP Data Sherif M. Hanafy ABSTRACT We present the theory and numerical results for interferometrically interpolating and extrapolating 3D marine
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationHTI anisotropy in heterogeneous elastic model and homogeneous equivalent model
HTI anisotropy in heterogeneous elastic model and homogeneous equivalent model Sitamai, W. Ajiduah*, Gary, F. Margrave and Pat, F. Daley CREWES/University of Calgary. Summary In this study, we compared
More informationElectromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah
Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation
More information(x, y, z) m 2. (x, y, z) ...] T. m 2. m = [m 1. m 3. Φ = r T V 1 r + λ 1. m T Wm. m T L T Lm + λ 2. m T Hm + λ 3. t(x, y, z) = m 1
Class 1: Joint Geophysical Inversions Wed, December 1, 29 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical, statistical Deal with data in the same
More informationChapter 5. 3D data examples REALISTICALLY COMPLEX SYNTHETIC INVERSION. Modeling generation and survey design
Chapter 5 3D data examples In this chapter I will demonstrate the e ectiveness of the methodologies developed in the previous chapters using 3D data examples. I will first show inversion results of a realistically
More informationComparison of Numerical Methods for Seismic Wave Propagation and Source Dynamics - the SPICE Code Validation
Third International Symposium on the Effects of Surface Geology on Seismic Motion Grenoble, France, 30 August - 1 September 2006 Paper Number: 62 Comparison of Numerical Methods for Seismic Wave Propagation
More informationMachine-learning Based Automated Fault Detection in Seismic Traces
Machine-learning Based Automated Fault Detection in Seismic Traces Chiyuan Zhang and Charlie Frogner (MIT), Mauricio Araya-Polo and Detlef Hohl (Shell International E & P Inc.) June 9, 24 Introduction
More informationTravel Time Tomography using Neural Networks
Travel Time Tomography using Neural Networks Yoshiya Oda Tokyo Metropolitan University, Japan Tomohisa Ishiyama Tokyo Metropolitan University, Japan Shinya Yokono Tokyo Metropolitan University, Japan SUMMARY:
More informationGeometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT. August Lau and Chuan Yin.
Geometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT August Lau and Chuan Yin January 6, 2017 Abstract The goal of seismic processing is to convert input data
More informationWS10-A03 Alternative Multicomponent Observables for Robust Full-waveform Inversion - Application to a Near-surface Example
1-4 June IFEMA Madrid WS-A3 Alternative Multicomponent Observables for Robust Full-waveform Inversion - Application to a Near-surface Example R. Valensi* (formerly IFSTTAR), R. Brossier (Univ. Grenoble
More informationReverse-time migration imaging with/without multiples
Reverse-time migration imaging with/without multiples Zaiming Jiang, John C. Bancroft, and Laurence R. Lines Imaging with/without multiples ABSTRACT One of the challenges with reverse-time migration based
More informationU043 3D Prestack Time Domain Full Waveform Inversion
U043 3D Prestack Time Domain Full Waveform Inversion D.V. Vigh* (WesternGeco), W.E.S. Starr (WesternGeco) & K.D. Kenneth Dingwall (WesternGeco) SUMMARY Despite the relatively high computational demand,
More informationOn the Scattering Effect of Lateral Discontinuities on AVA Migration
On the Scattering Effect of Lateral Discontinuities on AVA Migration Juefu Wang* and Mauricio D. Sacchi Department of Physics, University of Alberta, 4 Avadh Bhatia Physics Laboratory, Edmonton, AB, T6G
More informationRefraction Full-waveform Inversion in a Shallow Water Environment
Refraction Full-waveform Inversion in a Shallow Water Environment Z. Zou* (PGS), J. Ramos-Martínez (PGS), S. Kelly (PGS), G. Ronholt (PGS), L.T. Langlo (PGS), A. Valenciano Mavilio (PGS), N. Chemingui
More informationModel parametrization strategies for Newton-based acoustic full waveform
Model parametrization strategies for Newton-based acoustic full waveform inversion Amsalu Y. Anagaw, University of Alberta, Edmonton, Canada, aanagaw@ualberta.ca Summary This paper studies the effects
More informationCrosswell Imaging by 2-D Prestack Wavepath Migration
Crosswell Imaging by 2-D Prestack Wavepath Migration Hongchuan Sun ABSTRACT Prestack wavepath migration (WM) is applied to 2-D synthetic crosswell data, and the migrated images are compared to those from
More informationEffects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines
Effects of multi-scale velocity heterogeneities on wave-equation migration Yong Ma and Paul Sava, Center for Wave Phenomena, Colorado School of Mines SUMMARY Velocity models used for wavefield-based seismic
More informationAPPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA.
APPLICATION OF MATLAB IN SEISMIC INTERFEROMETRY FOR SEISMIC SOURCE LOCATION AND INTERPOLATION OF TWO DIMENSIONAL OCEAN BOTTOM SEISMIC DATA. BY: ISAAC KUMA YEBOAH. Department of Engineering, Regent University
More informationVerification and Validation for Seismic Wave Propagation Problems
Chapter 26 Verification and Validation for Seismic Wave Propagation Problems (1989-2-24-25-28-29-21-211-217-) (In collaboration with Dr. Nima Tafazzoli, Dr. Federico Pisanò, Mr. Kohei Watanabe and Mr.
More informationAn Efficient 3D Elastic Full Waveform Inversion of Time-Lapse seismic data using Grid Injection Method
10 th Biennial International Conference & Exposition P 032 Summary An Efficient 3D Elastic Full Waveform Inversion of Time-Lapse seismic data using Grid Injection Method Dmitry Borisov (IPG Paris) and
More informationSUMMARY INTRODUCTION NEW METHOD
Reverse Time Migration in the presence of known sharp interfaces Alan Richardson and Alison E. Malcolm, Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology SUMMARY
More informationNUMERICAL MODELING OF ACOUSTIC WAVES IN 2D-FREQUENCY DOMAINS
Copyright 2013 by ABCM NUMERICAL MODELING OF ACOUSTIC WAVES IN 2D-FREQUENCY DOMAINS Márcio Filipe Ramos e Ramos Fluminense Federal University, Niterói, Brazil mfrruff@hotmail.com Gabriela Guerreiro Ferreira
More informationLARGE-SCALE NORTHRIDGE EARTHQUAKE SIMULATION USING OCTREE-BASED MULTIRESOLUTION MESH METHOD
16th ASCE Engineering Mechanics Conference July 16-18, 3, University of Washington, Seattle LARGE-SCALE NORTHRIDGE EARTHQUAKE SIMULATION USING OCTREE-BASED MULTIRESOLUTION MESH METHOD Eui Joong Kim 1 Jacobo
More informationNumerical simulation of tunnel construction in volcanic rocks
Numerical simulation of tunnel construction in volcanic rocks G. Beer, Ch. Duenser Graz University of Technology, Graz, Austria ABSTRACT: For the design of tunnels, numerical simulation plays an important
More informationSource Estimation for Wavefield Reconstruction Inversion
Source Estimation for Wavefield Reconstruction Inversion Zhilong Fang * and Felix J. Herrmann * * Seismic Laboratory for Imaging and Modeling SLIM), University of British Columbia Abstract Wavefield reconstruction
More informationWave-equation migration from topography: Imaging Husky
Stanford Exploration Project, Report 123, October 31, 2005, pages 49 56 Short Note Wave-equation migration from topography: Imaging Husky Jeff Shragge 1 INTRODUCTION Imaging land seismic data is wrought
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 3, 2017, Lesson 1
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 3, 2017, Lesson 1 1 Politecnico di Milano, February 3, 2017, Lesson 1 2 Outline
More informationRobustness of the scalar elastic imaging condition for converted waves
CWP-830 Robustness of the scalar elastic imaging condition for converted waves Yuting Duan & Paul Sava Center for Wave Phenomena, Colorado School of Mines ABSTRACT For elastic reverse-time migration, one
More informationAdvances in radial trace domain coherent noise attenuation
Advances in radial trace domain coherent noise attenuation ABSTRACT David C. Henley* CREWES, Department of Geology and Geophysics University of Calgary, Calgary, AB henley@crewes.org The radial trace transform,
More informationComputational Seismology: Simulating Seismic Wavefields for AlpArray
Computational Seismology: Simulating Seismic Wavefields for AlpArray Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich 1 Introduction 2 Goals of lecture Introduction
More informationCoherent partial stacking by offset continuation of 2-D prestack data
Stanford Exploration Project, Report 82, May 11, 2001, pages 1 124 Coherent partial stacking by offset continuation of 2-D prestack data Nizar Chemingui and Biondo Biondi 1 ABSTRACT Previously, we introduced
More informationWe B3 12 Full Waveform Inversion for Reservoir Characterization - A Synthetic Study
We B3 12 Full Waveform Inversion for Reservoir Characterization - A Synthetic Study E. Zabihi Naeini* (Ikon Science), N. Kamath (Colorado School of Mines), I. Tsvankin (Colorado School of Mines), T. Alkhalifah
More informationA full-wave equation based seismic illumination analysis method
A full-wave equation based seismic illumination analysis method Xiao-Bi Xie and Hui Yang University of California, Santa Cruz, CA 95064, USA Summary We propose a full-wave equation based method for seismic
More informationFull-waveform inversion for reservoir characterization: A synthetic study
CWP-889 Full-waveform inversion for reservoir characterization: A synthetic study Nishant Kamath, Ilya Tsvankin & Ehsan Zabihi Naeini ABSTRACT Most current reservoir-characterization workflows are based
More informationis decomposed into P and S components localized in both space and slowness where P
eismic wave scattering in D random media: A finite-difference simulation and slowness domain analysis Xiao-Bi Xie* Institute of Geophysics and lanetary hysics, University of California at anta Cruz ummary
More informationFermat s principle and ray tracing in anisotropic layered media
Ray tracing in anisotropic layers Fermat s principle and ray tracing in anisotropic layered media Joe Wong ABSTRACT I consider the path followed by a seismic signal travelling through velocity models consisting
More informationG042 Subsalt Imaging Challenges - A Deepwater Imaging Analysis
G042 Subsalt Imaging Challenges - A Deepwater Imaging Analysis M. Cogan* (WesternGeco), J. Gardner (WesternGeco) & N. Moldoveanu (WesternGeco) SUMMARY Upon completion of the final reverse-time migration
More informationP052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges
P052 3D Modeling of Acoustic Green's Function in Layered Media with Diffracting Edges M. Ayzenberg* (Norwegian University of Science and Technology, A. Aizenberg (Institute of Petroleum Geology and Geophysics,
More information1D internal multiple prediction in a multidimensional world: errors and recommendations
1D internal multiple prediction 1D internal multiple prediction in a multidimensional world: errors and recommendations Pan Pan and Kris Innanen ABSTRACT Internal multiples are more difficult to estimate
More informationMASW Horizontal Resolution in 2D Shear-Velocity (Vs) Mapping
MASW Horizontal Resolution in 2D Shear-Velocity (Vs) Mapping by Choon B. Park Kansas Geological Survey University of Kansas 1930 Constant Avenue, Campus West Lawrence, Kansas 66047-3726 Tel: 785-864-2162
More informationFoolproof AvO. Abstract
Foolproof AvO Dr. Ron Masters, Geoscience Advisor, Headwave, Inc Copyright 2013, The European Association of Geoscientists and Engineers This paper was prepared for presentation during the 75 th EAGE Conference
More information2012 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
EFFECTS OF 3D SURFACE TOPOGRAPHY ON REGIONAL AND TELESEISMIC SIGNALS FROM UNDERGROUND EXPLOSIONS Thorne Lay 1, Megan Avants 1, Xiao-bi Xie 1, and Arthur J. Rodgers 2 University of California Santa Cruz
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Lian Shen Department of Civil Engineering Johns Hopkins University Baltimore, MD 21218 phone: (410) 516-5033 fax: (410) 516-7473 email: LianShen@jhu.edu
More informationInversion after depth imaging
Robin P. Fletcher *, Stewart Archer, Dave Nichols, and Weijian Mao, WesternGeco Summary In many areas, depth imaging of seismic data is required to construct an accurate view of the reservoir structure.
More informationSUMMARY. method to synthetic datasets is discussed in the present paper.
Geophysical modeling through simultaneous Joint Inversion of Seismic, Gravity and Magnetotelluric data Michele De Stefano (1), Daniele Colombo (1) WesternGeco EM - Geosystem, via Clericetti 42/A, 20133
More informationSeismic Modeling, Migration and Velocity Inversion
Seismic Modeling, Migration and Velocity Inversion Aliasing Bee Bednar Panorama Technologies, Inc. 14811 St Marys Lane, Suite 150 Houston TX 77079 May 16, 2014 Bee Bednar (Panorama Technologies) Seismic
More informationNUMERICAL MODELING OF STRONG GROUND MOTION USING 3D GEO-MODELS
th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 24 Paper No. 92 NUMERICAL MODELING OF STRONG GROUND MOTION USING D GEO-MODELS Peter KLINC, Enrico PRIOLO and Alessandro
More informationJournal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp. 59 74 SOLID MECHANICS WAVE PROPAGATION DUE TO AN EMBEDDED SEISMIC SOURCE IN A GRADED HALF-PLANE WITH RELIEF PECULIARITIES.
More informationA hybrid multiresolution meshing technique for finite element three-dimensional earthquake ground motion modelling in basins including topography
Geophys. J. Int. (29) 177, 1221 1232 doi: 1.1111/j.1365-246X.29.4154.x A hybrid multiresolution meshing technique for finite element three-dimensional earthquake ground motion modelling in basins including
More informationA NUMERICAL EXPERIMENT ON NONLINEAR IMAGE RECONSTRUCTION FROM FIRST-ARRIVAL TIMES FOR TWO-DIMENSIONAL ISLAND ARC STRUCTURE
Letter J. Phys. Earth, 34, 195-201, 1986 A NUMERICAL EXPERIMENT ON NONLINEAR IMAGE RECONSTRUCTION FROM FIRST-ARRIVAL TIMES FOR TWO-DIMENSIONAL ISLAND ARC STRUCTURE Ichiro NAKANISHI and Kaoru YAMAGUCHI
More informationElastic full waveform Inversion for land walkaway VSP data
Elastic full waveform Inversion for land walkaway VSP data Olga Podgornova, Scott Leaney, arwan Charara, Schlumberger; Eric von Lunen, Nexen Energy ULC Summary Full waveform inversion (FWI) is capable
More informationContinued Fraction Absorbing Boundary Conditions for Transient Elastic Wave Propagation Modeling
Continued Fraction Absorbing Boundary Conditions for Transient Elastic Wave Propagation Modeling Md Anwar Zahid and Murthy N. Guddati 1 Summary This paper presents a novel absorbing boundary condition
More information3D Finite Difference Time-Domain Modeling of Acoustic Wave Propagation based on Domain Decomposition
3D Finite Difference Time-Domain Modeling of Acoustic Wave Propagation based on Domain Decomposition UMR Géosciences Azur CNRS-IRD-UNSA-OCA Villefranche-sur-mer Supervised by: Dr. Stéphane Operto Jade
More informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationMain Menu. Summary. Introduction
Local-angle domain illumination for full-wave propagators Jun Cao* and Ru-Shan Wu Department of Earth and Planetary Sciences/IGPP, University of California, Santa Cruz Summary We propose an efficient split-step
More informationDetermination of 2D shallow S wave velocity profile using waveform inversion of P-SV refraction data
Determination of 2D shallow S wave velocity profile using waveform inversion of P-SV refraction data Mohamed Amrouche and Hiroaki Yamanaka Dept. of Environmental Science and Technology, Tokyo Institute
More informationAmplitude within the Fresnel zone for the zero-offset case
Amplitude within the Fresnel zone for the zero-offset case Shuang Sun and John C. Bancroft ABSTRACT Reflection energy from a linear reflector comes from an integrant over an aperture often described by
More informationHeadwave Stacking in Terms of Partial Derivative Wavefield
Geosystem Engineering, 7(1), 21-26 (March 2004) Headwave Stacking in Terms of Partial Derivative Wavefield Changsoo Shin School of Civil, Urban and Geosystem Engineering, Seoul National University, San
More informationSUMMARY INTRODUCTION THEORY
Illumination compensation for subsalt image-domain wavefield tomography Tongning Yang, Center for Wave Phenomena, Colorado School of Mines, Jeffrey Shragge, School of Earth and Environment, University
More informationProcessing converted-wave data in the tau-p domain: rotation toward the source and moveout correction
τ-p domain converted-wave processing Processing converted-wave data in the tau-p domain: rotation toward the source and moveout correction Raul Cova and Kris Innanen ABSTRACT The asymmetry of the converted-wave
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Underwater Acoustics Session 2aUW: Wave Propagation in a Random Medium
More informationNumerical investigation of implementation of air-earth boundary by acoustic-elastic boundary approach
GEOPHYSICS, VOL. 72, NO. 5 SEPTEMBER-OCTOBER 27 ; P. SM147SM153, 5 FIGS., 5 TABLES. 1.119/1.2753831 Numerical investigation of implementation of air-earth boundary by acoustic-elastic boundary approach
More informationDownloaded 05/09/13 to Redistribution subject to SEG license or copyright; see Terms of Use at
Elastic converted-wave path migration for subsalt imaging Ru-Shan Wu*, Rui Yan, Xiao-Bi Xie, Modeling and Imaging Laboratory, Earth and Planetary Sciences/IGPP, University of California, Santa Cruz, David
More information