In this paper, we record P- and S-waveforms in three sandstones used in the Scottish heritage and obtain the P- and S-velocities as well as the calculation of wavelength of the P-waves in both fresh and salt weathered samples. We use an automatic methodology for the calculation of onset time that includes the signal pre-processing and the analysis, in the time-domain, of the first pulse symmetry, amplitude and duration criteria. This triple check provides greater confidence in the obtained results. This method limits human error, while improving the accuracy and reproducibility of P and S wave onset time estimations.
Materials and Methods
Stones
In this paper three sandstones used in the Scottish heritage are investigated:
Doddington sandstone (D) is a common building stone used throughout Scotland and is currently used as a replacement stone at Jedburgh Abbey, Scotland. The sandstone is quarried at Doddington quarry near Wooler in the Scottish Borders region; and forms part of the Fell formation, deposited during the early-mid Carboniferous. The sandstone is mineralogically mature and textural immature; with rounded-sub rounded grain shape and an average grain size of 0.25 mm. Mineralogy consists of: quartz (90 %); feldspar (mainly orthoclase) (10 %); and minor amounts of muscovite, lithic fragments and clay within the matrix. The sandstone is classified as a quartz-arenite.
St. Bees sandstone (BC) is currently used as a replacement stone at Arbroath Abbey and is quarried on the west coast of Cumbria, England. The red-dark red sandstone is part of the Chester sandstone formation, deposited during the early Triassic. The sandstone is texturally and mineralogically immature; with (sub)rounded – (sub)angular grain shape and an average grain size of 0.12 mm. The sandstone is classified as a lithic arkose. Mineralogy consists of: quartz (60 %); feldspar (orthoclase and plagioclase) (20 %); lithic fragments (10 %); mica (biotite and muscovite) (10 %); and clay-rich matrix. BC exhibits planar 193bedding on a scale of 1–2 mm. Therefore, during this study BC samples were measured in two directions: parallel and perpendicular to bedding, to determine the degree of anisotropy and potential effects on output signals.
Forest of Dean (F) is a distinctive grey-green sandstone used throughout the UK, for example in the restoration of Dunkeld Cathedral, Scotland. The sandstone forms part of the Pennant sandstone formation deposited during the mid-late Carboniferous; and is currently quarried at Barnhill quarry in the county of Gloucestershire, England. The sandstone is mineralogically and texturally immature; with (sub) angular grain shape and an average grain size of 0.19 mm. Mineralogy consists of: quartz (70 %); feldspar (orthoclase and plagioclase) (10 %); lithic fragments (15 %); muscovite (5 %); and a clay-rich matrix. The sandstone is classified as a (sub)litharenite.
The ultrasonic measurements were carried out by means of the transmission method, which consists of two piezoelectric sensors coupled to the sample at constant pressure. Compressive (P) and shear (S) waves were measured using polarised Panametric transducers (1 MHz). Emitting-receiving equipment (Panametrics-NDT 5058PR) and an oscilloscope (TDS 3012B-Tektronix) were used to acquire and digitalize the waveforms to be displayed, manipulated and stored. Two types of visco-elastic couplants were used to achieve good coupling between the transducer and the sample: one of a fluid consistency for the P-wave transducers (eco-gel), and another of a more viscous consistency for the S-waves (SWC, shear wave couplant, GE Panametrics).
In this work, we have applied a new method for automatic characterisation of the first pulse in both compressional and shear waves described in Benavente et al., (2020). Firstly, the recorded signal is pre-processed in the wavelet domain. This step filters noise and removes low-frequency disruption. Secondly, the recorded signal is analysed in the time-domain. This step identifies and characterises the first pulse, allowing for an estimation of onset time. Using the automatic approach allows for all detected pulses in the output signal to be analysed, where the first pulse of the P or S wave is determined based on criteria related to symmetry, amplitude and duration. The method outlined above is carried out using Matlab based on Galiana-Merino et al. (2013), where parts of the pre-processing stage are carried out using the Wavelet Toolbox. This triple check provides greater confidence in the results. We obtain the automatic onset time from recorded P- and S- waveforms and they are compared to manual picking, which is considered as a true or reference value.
Finally, the wavelength, λ, was calculated using the velocity, v, and frequency, f, of the first pulse as λ=v/f.
Salt crystallisation test
Salt crystallisation tests were in concordance with EN-12370 (1999) recommendations, respectively. 4 cm cubic samples of each lithotype were tested and underwent cycles of saline immersion (14 % w/w Na 2SO 4solution, at 20 °C for 4 hours), drying (at 60 °C for 16 hours) and cooling at room conditions (20 °C, for 4 hours). Samples were exclusively cleaned at the end of the test (after the 15 cycles). P and S waves were measured before and after the durability test. Every measurement of the P and S waves was repeated three times in order to test the reproducibility of the experiments and the corresponding results. Figures 1 and 2 display, respectively, P and S signals for fresh and weathered samples. They aim to show the evolution of waveform rather than wave velocities because the sample length may vary from samples to sample. Tested samples have a cubic shape and the edge length range from 40–50 mm. Moreover, weathered samples may also suffer superficial weight loss, resulting in a decrease in edge length.
For the studied sandstones, the arrival time of P-wave signals are easily identifiable with a high SNR, although their waveforms and wavelength values are dependent on the type of rock being 194analysed. The frequency of the samples decreases with particle size (Fig. 1). Inversely, coarse-grained rocks have the highest wavelength values for the studied rocks ( Table 1) and also the degree of salt crystallization. The wavelength of output signals depends on grain size, while microstructural components of the rocks such as pores, fractures, grains and the presence of salts operate as a wavelength filter. The frequency associated with the first pulse of the output signal differs from the central frequency of the input elastic waves which is fixed at 1 MHz. This is the result of the interaction between elastic wave and microstructural rock components. Manual picking of the onset time for elastic waves becomes more difficult with increase in wavelength, although this observation is less important for P-waves than S-waves (Figs. 1 and 2).
Table 1: Automatic mean values of grain size, P- (VP) and S-wave (VS) and wavelength (λ) for fresh and weathered sandstones.
The recorded signals highlight that the microstructural components of rocks and their modification by salt crystallization affect the output signal. Values for manual and calculated P-wave are almost equal (Fig. 3). Without concluding which results are closer to reality, one may compare the results of a proposed method with the results obtained by a human analyst (Sarout et al., 2009). If we assume that manual measuring offers a true or reference value (e. g.: Siegesmund and Dürrast, 2011), then it can be concluded that the proposed method accurately calculates P-waves velocities for a range of studied rocks.
Читать дальше