A numerical study of effective velocities in fractured media using the rotated staggered finite difference grid

2000 ◽  
Author(s):  
Erik H. Saenger ◽  
Serge A. Shapiro
Keyword(s):  
Finite Difference ◽  
Numerical Study ◽  
Fractured Media ◽  
Difference Grid

Finite difference code for velocity and surface traction of a Fluid between Two Eccentric Spheres

2015 ◽  
Vol 11 (1) ◽  
pp. 2960-2971
Author(s):  
M.Abdel Wahab
Keyword(s):  
Velocity Field ◽  
Angular Velocity ◽  
Finite Difference ◽  
Numerical Study ◽  
Outer Sphere ◽  
Equation Of Motion ◽  
Annular Region ◽  
Surface Traction ◽  
Angular Velocities ◽  
Axis Passing

The Numerical study of the flow of a fluid in the annular region between two eccentric sphere susing PHP Code isinvestigated. This flow is created by considering the inner sphere to rotate with angular velocity 1  and the outer sphererotate with angular velocity 2  about the axis passing through their centers, the z-axis, using the three dimensionalBispherical coordinates (, ,) .The velocity field of fluid is determined by solving equation of motion using PHP Codeat different cases of angular velocities of inner and outer sphere. Also Finite difference code is used to calculate surfacetractions at outer sphere.

Simulation of a seismic survey with combined surface and in-mine data acquisition for imaging steeply dipping ore veins

2021 ◽  
Author(s):  
Olaf Hellwig ◽  
Stefan Buske
Keyword(s):  
Finite Difference ◽  
Message Passing Interface ◽  
Seismic Survey ◽  
Difference Operators ◽  
Mining District ◽  
Prestack Depth Migration ◽  
Data Set ◽  
Depth Migration ◽  
Difference Grid ◽  
Triangulated Surfaces

<p>The polymetallic, hydrothermal deposit of the Freiberg mining district in the southeastern part of Germany is characterised by ore veins that are framed by Proterozoic orthogneiss. The ore veins consist mainly of quarz, sulfides, carbonates, barite and flourite, which are associated with silver, lead and tin. Today the Freiberg University of Mining and Technology is operating the shafts Reiche Zeche and Alte Elisabeth for research and teaching purposes with altogether 14 km of accessible underground galleries. The mine together with the most prominent geological structures of the central mining district are included in a 3D digital model, which is used in this study to study seismic acquisition geometries that can help to image the shallow as well as the deeper parts of the ore-bearing veins. These veins with dip angles between 40° and 85° are represented by triangulated surfaces in the digital geological model. In order to import these surfaces into our seismic finite-difference simulation code, they have to be converted into bodies with a certain thickness and specific elastic properties in a first step. In a second step, these bodies with their properties have to be discretized on a hexahedral finite-difference grid with dimensions of 1000 m by 1000 m in the horizontal direction and 500 m in the vertical direction. Sources and receiver lines are placed on the surface along roads near the mine. A Ricker wavelet with a central frequency of 50 Hz is used as the source signature at all excitation points. Beside the surface receivers, additional receivers are situated in accessible galleries of the mine at three different depth levels of 100 m, 150 m and 220 m below the surface. Since previous mining activities followed primarily the ore veins, there are only few pilot-headings that cut through longer gneiss sections. Only these positions surrounded by gneiss are suitable for imaging the ore veins. Based on this geometry, a synthetic seismic data set is generated with our explicit finite-difference time-stepping scheme, which solves the acoustic wave equation with second order accurate finite-difference operators in space and time. The scheme is parallelised using a decomposition of the spatial finite-difference grid into subdomains and Message Passing Interface for the exchange of the wavefields between neighbouring subdomains. The resulting synthetic seismic shot gathers are used as input for Kirchhoff prestack depth migration as well as Fresnel volume migration in order to image the ore veins. Only a top mute to remove the direct waves and a time-dependent gain to correct the amplitude decay due to the geometrical spreading are applied to the data before the migration. The combination of surface and in-mine acquisition helps to improve the image of the deeper parts of the dipping ore veins. Considering the limitations for placing receivers in the mine, Fresnel volume migration as a focusing version of Kirchhoff prestack depth migration helps to avoid migration artefacts caused by this sparse and limited acquisition geometry.</p>

Finite‐difference modeling of faults and fractures

Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1514-1526 ◽  
Author(s):  
Richard T. Coates ◽  
Michael Schoenberg
Keyword(s):  
Finite Difference ◽  
Seismic Wave ◽  
Computation Time ◽  
Anisotropic Materials ◽  
Model Material ◽  
Material Behavior ◽  
Interface Condition ◽  
Medium Theory ◽  
Difference Grid ◽  
Receiver Array

For the purposes of seismic propagation, a slip fault may be regarded as a surface across which the displacement caused by a seismic wave is discontinuous while the stress traction remains continuous. The simplest assumption is that this slip and the stress traction are linearly related. Such a linear slip interface condition is easily modeled when the fault is parallel to the finite‐difference grid, but is more difficult to do for arbitrary nonplanar fault surfaces. To handle such situations we introduce equivalent medium theory to model material behavior in the cells of the finite‐ difference grid intersected by the fault. Virtually identical results were obtained from modeling the fault by (1) an explicit slip interface condition (fault parallel to the grid) and (2) using the equivalent medium theory when the finite‐difference grid was rotated relative to the fault and receiver array. No additional computation time is needed except for the preprocessing required to find the relevant cells and their associated moduli. The formulation is sufficiently general to include faults in and between arbitrary anisotropic materials with slip properties that vary as a function of position.

scholarly journals Density Difference Grid Design in a Point-Mass Filter

Energies ◽  
2020 ◽  
Vol 13 (16) ◽  
pp. 4080
Author(s):  
Jakub Matoušek ◽  
Jindřich Duník ◽  
Ondřej Straka
Keyword(s):  
Numerical Study ◽  
The State ◽  
Point Mass ◽  
Density Difference ◽  
Conditional Density ◽  
Stochastic Dynamic ◽  
Mass Filter ◽  
Difference Grid ◽  
Grid Design ◽  
Stochastic Dynamic Systems

The paper deals with the Bayesian state estimation of nonlinear stochastic dynamic systems. The stress is laid on the point-mass filter, solving the Bayesian recursive relations for the state estimate conditional density computation using the deterministic grid-based numerical integration method. In particular, the grid design is discussed and the novel density difference grid is proposed. The proposed grid design covers such regions of the state-space where the conditional density is significantly spatially varying, by the dense grid. In other regions, a sparse grid is used to keep the computational complexity low. The proposed grid design is thoroughly discussed, analyzed, and illustrated in a numerical study.

Sign in / Sign up

Export Citation Format

Share Document