DOE OSTI2012
The Source Physics Experiment (SPE) is designed to improve ground motion prediction capabilities for underground nuclear testing through the development of new physics-based models together with validation via new chemical explosion data. The SPE shots consist of different size and depth chemical explosions in a granite borehole at the Climax Stock on the Nevada National Security Site (NNSS). The third shot in the series (SPE3) will be the same size (1000 kg TNT equivalent) and depth (45 m canister centroid) as the previous shot (SPE2). The main purpose of the SPE3 is to investigate the effects of the damage zone caused by the SPE2 chemical explosion on the SPE3 explosion, keeping all other parameters (e.g. size, depth, explosive type) the same. This SPE3 pre-shot report addresses progress made in simulation of the expected wavefield in the both the near- and the far-field. In this report we consider the far-field to be where seismic wave motion can be characterized by elastic or visco-elastic approximations, in contrast to the near-field where other, often non-linear effects, play a major role. The current simulations build on previous work from the SPE1 and SPE2 shots. The simulations are intended to highlight differences between SPE2 and SPE3 using a variety to approximations intended to mimic the expected effect of the damage zone. The simplest physics of damage zone effects suggest changes in the source region material properties (density, velocity), and in the yield strength (coupling) of the rock, should both impact the far-field signal. We can illustrate this in a simplistic way for P-waves in Figure A using the empirical Mueller-Muphy (1971) far-field explosion model (MM). First if the SPE2 explosion-induced damage reduces velocities and densities without affecting the yield strength or coupling, then by energy conservation the amplitudes for SPE 3 would increase as shown in the top of Figure A. Similarly if the explosion induced damage increased density (say by compacting pore space) the amplitude would decrease. This impedance effect, where the product of the density, wave-velocity, and ground-velocity squared is conserved, is well known. The second and harder effect to quantify, are changes in the energy coupling or yield strength. In the simple MM model this is represented by the coupling factor between cavity radius and yield. If we assume the granite was slightly weakened by SPE2 such that the cavity radius increases slightly then the overall amplitudes of SPE3 decrease relative to SPE2 as shown in the bottom of Figure A. In reality both effects, may exist due to the damage zone, and furthermore we know they are interrelated in complex ways. For S-waves we expect even more complex effects due to motion on joints and faults. For this reason we are developing physics-based near field simulation models coupled to far-field propagation codes to predict the overall effects of changes in the source emplacement region on observed signals. Near-field simulations have been performed using the Eulerian code (GEODYN). An ensemble model, calibrated to match records from past large-scale nuclear events is used. Both 2D axisymmetric (no faults) and fully three-dimensional GEODYN calculations (with /without faults) were conducted in which the damage zones sustained by SPE2 are retained as the initial condition for the SPE3 calculations. The results indicate that the waveform peak amplitudes are slightly decayed since the ensemble model includes no softening mechanism. The pulse widths are reduced, indicating the energy is dissipated more by the damage zone as expected. The dominant motions calculated in all the scenarios are spherically symmetric and non-spherical components are small as seen for the SPE2 calculations. The arrival time of the SPE3 waveform at the boreholes is shorter than SPE2 due to compaction during the previous explosion. As shown in the SPE 1 and 2 observations shear waves are generated, indicating the potential contributions from the induced motion on distributed joints. The effects of joints are evaluated using the Lagrangian hydrocode (GEODYN-L), which was benchmarked against the GEODYN calculations for the ensemble model. A new granite model calibrated for samples from the SPE site is used. The new model includes effects of poroelasticity and water saturation and may be scaled allowing application to both small and large scales. Parametric tests were performed to evaluate the effect of joints within the granite and the results demonstrate that in situ joints introduce anisotropy and significant tangential motions.