• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.98-104
• /
• 2011

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and threedimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.60-66
• /
• 2006

Three differing sandstones, two synthetic and one field sample, have been tested ultrasonically under a range of confining pressures and pore pressures representative of in-situ reservoir pressures. These sandstones include: a synthetic sandstone with calcite intergranular cement produced using the CSIRO Calcite In-situ Precipitation Process (CIPS); a synthetic sandstone with silica intergranular cement; and a core sample from the Otway Basin Waarre Formation, Boggy Creek 1 well, from the target lithology for a trial $CO_2$ pilot project. Initial testing was carried on the cores at "room-dried" conditions, with confining pressures up to 65 MPa in steps of 5 MPa. All cores were then flooded with $CO_2$, initially in the gas phase at 6 MPa, $22^{\circ}C$, then with liquid-phase $CO_2$ at a temperature of $22^{\circ}C$ and pressures from 7 MPa to 17 MPa in steps of 5 MPa. Confining pressures varied from 10 MPa to 65 MPa. Ultrasonic waveforms for both P- and S-waves were recorded at each effective pressure increment. Velocity versus effective pressure responses were calculated from the experimental data for both P- and S-waves. Attenuations $(1/Q_p)$ were calculated from the waveform data using spectral ratio methods. Theoretical calculations of velocity as a function of effective pressure for each sandstone were made using the $CO_2$ pressure-density and $CO_2$ bulk modulus-pressure phase diagrams and Gassmann effective medium theory. Flooding the cores with gaseous phase $CO_2$ produced negligible change in velocity-effective stress relationships compared to the dry state (air saturated). Flooding with liquid-phase $CO_2$ at various pore pressures lowered velocities by approximately 8% on average compared to the air-saturated state. Attenuations increased with liquid-phase $CO_2$ flooding compared to the air-saturated case. Experimental data agreed with the Gassmann calculations at high effective pressures. The "critical" effective pressure, at which agreement with theory occurred, varied with sandstone type. Discrepancies are thought to be due to differing micro-crack populations in the microstructure of each sandstone type. The agreement with theory at high effective pressures is significant and gives some confidence in predicting seismic behaviour under field conditions when $CO_2$ is injected.