Ultrasonic Transmission Measurements Research Articles

Poroelastic Response of a Fractured Rock to Hydrostatic Pressure Oscillations

AbstractPoroelastic coupling between fractures and the surrounding rock is important to numerous applications in geosciences. We measure the in‐situ fluid pressure and local strain response of a fractured carbonate sample to hydrostatic pressure oscillations. A linear poroelastic model that represents the rock sample is parameterized using X‐ray imaging and ultrasonic wave transmission measurements. The numerical solution, based on Biot's quasistatic equations, is consistent with the measured frequency dependent dispersion of the apparent bulk modulus of the background matrix and the in‐situ pore pressure response, which is caused by fluid pressure diffusion from the compliant fractures into the stiffer matrix. The observed fluid pressure diffusion is causally related to the numerically quantified intrinsic attenuation at seismic frequencies, which is a major contributor to the dissipation of seismic waves. Our analysis supports the use of a simple approximation of fractures as compliant and planar inclusions in numerical simulations based on linear poroelasticity.

Ultrasonic Methods for Characterizing Layered Electrodes of Lithium-Ion Batteries

Lithium-ion battery electrodes have a multilayered structure consisting of active porous electrode materials coated on one or both sides of a thin metal current collector. The porosity and tortuosity of the active materials play a vital role in determining the battery's performance metrics, such as energy density and power density. However, existing methods for quantifying these parameters are often limited and cannot be easily applied in a production environment. In this work, we propose two ultrasonic methods based on a newly developed stiffness matrix theory [1] to achieve in-production quantification. The first method establishes a quantitative relationship to link angled ultrasound transmission to the porous parameters of electrodes. This relationship is then used to inversely determine porosity and tortuosity from ultrasonic transmission measurements using a non-contact, air-coupled setup. From the first method, we observe that ultrasound is highly transmissible as guided wave modes are excited in electrodes. Based on this phenomenon, our second method directly excites guided waves in the layers using electromagnetic acoustic transducers (EMATs), enabling us to obtain more accurate measurements with a better signal-to-noise ratio. Both methods have demonstrated good sensitivity to porosity and tortuosity, thus providing accurate quantification using ultrasonic transmission and guided wave results. We anticipate that our methods have the potential to provide real-time online characterization of electrodes during production, which would enable closed-loop control of the manufacturing processes and thus dramatically reduce wastage and improve product quality. Reference [1] M. Huang, F. Cegla, and B. Lan, Stiffness matrix method for modelling wave propagation in arbitrary multilayers, http://arxiv.org/abs/2302.12868.

Performance Evaluation of AE Sensors Installed Like Hydrophones in Adaptive Monitoring Networks During a Decametre-Scale Hydraulic Stimulation Experiment

AbstractIn the framework of the STIMTEC and STIMTEC-X hydraulic stimulation experiments at the Reiche Zeche mine, Freiberg (Germany), we installed acoustic emission (AE) sensors for the recording of picoseismicity both conventionally using pneumatic coupling and experimentally like a hydrophone, i.e. the sensors were placed in the borehole without a further coupling system or cementing. We investigate performance measures of the hydrophone-like acoustic emission (HAE) sensors such as frequency bandwidth, sensitivity, first motion polarity, coupling and placement quality to assess the sensor’s applicability in adaptive monitoring networks. HAE sensors can be paired with hydraulic equipment, especially with the double packer probe used for stimulation at the decametre scale because the monitored frequency content differs from injection-related noise. This offers a unique opportunity to improve the network geometry and consequently the quality of a seismic catalogue. We analyse the sensor characteristics using active ultrasonic transmission measurements from boreholes with different orientations in the rock volume, noise measurements preceding active centre punch hits in the access galleries and passive recordings of induced acoustic emission events. HAE sensors placed in water-filled boreholes show good sensitivity performance even without optimal coupling to the crystalline rock for recording distances up to 17 m. The HAE sensors record the wavefield adequately for first-arrival identification, polarity picking and amplitude characteristics but are less suitable for detecting S-waves. Due to the borehole geometry HAE sensors record waves with incidence angles from the side, resulting in opposite polarity compared to side-view AE sensors as observed in the field and lab. We discuss the advantages of adaptive monitoring networks with HAE sensors being optimally placed for each stimulation interval configuration anew to improve seismic event detection and quality of event hypocentre locations during hydraulic stimulations. We show that we are able to significantly reduce the azimuthal gap, halve the location uncertainties and improve the network coverage for the purpose of focal mechanism estimations.

Metre-scale damage zone characterization usingS-coda waves from active ultrasonic transmission measurements in the STIMTEC project, URL Reiche Zeche, Germany

SUMMARYStudies of controlled hydraulic stimulation experiments with active and passive seismic monitoring conducted in Underground Research Laboratories (URLs) benefit from specific knowledge of hydraulic parameters, close by microseismic monitoring revealing structural details of the rock mass, and detailed evolution of seismicity in response to injection operations. Microseismic monitoring is commonly used to characterize a stimulated reservoir volume, for example, in terms of damage evolution of the rock mass. Since seismic attenuation is affected by damage of the rock volume, active seismic sources covering sizes from the centimetre to decimetre scale may help us to investigate space–time varying attenuation properties in a reservoir. This may allow us to monitor damage evolution of the stimulated rock volume in more detail, also since active seismic sources produce stronger signals leading to a broader frequency range that can be analysed compared to passive seismic signals. Within the STIMTEC project in the URL Reiche Zeche (URL-RZ) in Freiberg (Germany), more than 300 active Ultrasonic Transmission (UT) measurements were performed before and after hydraulic stimulations in two boreholes in the targeted rock volume, an anisotropic metamorphic gneiss. The signal-frequency content ranges between 1 and 60 kHz. Assuming scattering attenuation to dominate over intrinsic attenuation, we here apply the single isotropic scattering model. S-coda waves of 88 spatially representative UT measurements are used to estimate the coda quality factor (QC). We obtain stable QC estimates for centre frequencies of octave-width frequency bands between 3 and 21 kHz. We group neighbouring UT measurements to stabilize the observations and form eight UT groups in total, covering different depth intervals in three boreholes and four different time periods to investigate scattering attenuation changes in a spatiotemporal manner. Our final mean QC ($\overline {{Q}_C} $) estimates show characteristic frequency-dependence as observed at the field scale in geological reservoirs. We find temporal variations of QC are strongly connected to hydraulic stimulation, and these variations are more significant than those resolved from velocity changes. $\overline {{Q}_C} $ estimates at frequencies above 15 kHz indicate healing of injection-induced small-scale fractures during a two-months post-stimulation phase. Larger fractures, mostly sampled by lower frequencies (<15 kHz), seem to be more persistent with time (over 15 months). We observe spatial differences of $\overline {{Q}_C} $ values near the mine galleries (driftway and vein drift) and relate these observations to different extents and characteristics of the galleries’ excavation damage zones. Our results further support previous assumptions based on borehole televiewer logs and mapped structures of an existing fault with larger damage zone that crosses the stimulated rock volume NW-SE between the galleries. We conclude that the coda analysis of active UT measurements complements established imaging methods used during experiments in URLs. In particular, coda analysis is a powerful tool for the detection of damage zones and for monitoring local fracture networks with immediate application for imaging georeservoirs considered for exploitation or underground storage of gases and liquids.

Ultrasonic Through-Transmission Measurements of Human Musculoskeletal and Fat Properties.

The study described here was aimed at investigating the feasibility of using the ultrasonic through-transmission technique to estimate human musculoskeletal and fat properties. Five hundred eighty-two volunteers were assessed by dual-energy X-ray absorptiometry (DXA) and ultrasonic transmission techniques. Bone mineral density (BMD), muscle and fat mass were measured for both legs and the whole body. Hip BMD and spine BMD were also measured. Ultrasonic transmission measurements were performed on the heel, and the measured parameters were broadband ultrasound attenuation (BUA), speed of sound (SOS), ultrasonic stiffness index (SI), T-score and Z-score, which were significantly correlated with all measured BMDs. The optimal correlation was observed between SI and left-leg BMD (p < 0.001) before and after adjustment for age, sex and body mass index (BMI). The linear and partial correlation analyses revealed that BUA and SOS were closely associated with muscle and fat mass, respectively. Multiple regressions revealed that muscle and fat mass significantly contributed to the prediction of transmission parameters, explaining up to 17.83% (p < 0.001) variance independently of BMD. The results suggest that the ultrasonic through-transmission technique could help in the clinical diagnosis of skeletal and muscular system diseases.

Seismic monitoring of the STIMTEC hydraulic stimulation experiment in anisotropic metamorphic gneiss

Abstract. In 2018 and 2019, we performed STIMulation tests with characterising periodic pumping tests and high-resolution seismic monitoring for improving prognosis models and real-time monitoring TEChnologies for the creation of hydraulic conduits in crystalline rocks (STIMTEC). The STIMTEC underground research laboratory is located at 130 m depth in the Reiche Zeche mine in Freiberg, Germany. The experiment was designed to investigate the rock damage resulting from hydraulic stimulation and to link seismic activity and enhancement of hydraulic properties in strongly foliated metamorphic gneiss. We present results from active and passive seismic monitoring prior to and during hydraulic stimulations. We characterise the structural anisotropy and heterogeneity of the reservoir rocks at the STIMTEC site and the induced high-frequency (&gt;1 kHz) acoustic emission (AE) activity, associated with brittle deformation at the centimetre-to-decimetre scale. We derived the best velocity model per recording station from over 300 active ultrasonic transmission measurements for high-accuracy AE event location. The average P-wave anisotropy is 12 %, in agreement with values derived from laboratory tests on core material. We use a 16-station seismic monitoring network comprising AE sensors, accelerometers, one broadband sensor and one AE hydrophone. All instrumentation was removable, providing us with the flexibility to use existing boreholes for multiple purposes. This approach also allowed for optimising the (near)-real-time passive monitoring system during the experiment. To locate AE events, we tested the effect of different velocity models and inferred their location accuracy. Based on the known active ultrasonic transmission measurement points, we obtained an average relocation error of 0.26±0.06 m where the AE events occurred using a transverse isotropic velocity model per station. The uncertainty resulting from using a simplified velocity model increased to 0.5–2.6 m, depending on whether anisotropy was considered or not. Structural heterogeneity overprints anisotropy of the host rock and has a significant influence on velocity and attenuation, with up to 4 % and up to 50 % decrease on velocity and wave amplitude, respectively. Significant variations in seismic responses to stimulation were observed ranging from abundant AE events (several thousand per stimulated interval) to no activity with breakdown pressure values ranging between 6.4 and 15.6 MPa. Low-frequency seismic signals with varying amplitudes were observed for all stimulated intervals that are more correlated with the injection flow rate rather than the pressure curve. We discuss the observations from STIMTEC in context of similar experiments performed in underground research facilities to highlight the effect of small-scale rock, stress and structural heterogeneity and/or anisotropy observed at the decametre scale. The reservoir complexity at this scale supports our conclusion that field-scale experiments benefit from high-sensitivity, wide-bandwidth instrumentation and flexible monitoring approaches to adapt to unexpected challenges during all stages of the experiment.

Laboratory seismic measurements for layer-specific description of fluid mud and for linking seismic velocities to rheological properties.

The velocities of the seismic waves propagating in the fluid-mud layer are governed by the rheological properties and density of the fluid mud. Performing seismic transmission measurements inside the fluid mud can give good estimates of the seismic velocities and, thus, of the rheological properties and density. Laboratory ultrasonic transmission measurements of the wave velocities in the fluid-mud layer and their temporal evolution are shown. It is found that the shear-wave velocity and yield stress are positively correlated. Performing a seismic reflection survey for characterization of the fluid-mud layers could be more practical because it allows towing the sources and receivers above the top of fluid-mud layer. Interpretation of the results from a reflection survey, though, is influenced by the water layer above the fluid mud. Applying seismic interferometry to reflection measurements can eliminate the influence of the water layer and retrieve a reflection response from inside the fluid-mud layer. This eliminates the influence of the temperature and salinity of the water layer to obtain information about the seismic properties of the fluid-mud layer. To introduce the approach of retrieving and extracting the reflection response from inside the fluid-mud layer, data from laboratory measurements are used. The obtained compressional- and shear-wave velocities are validated by comparing them with values from current transmission measurements.

Analysis of Eshelby-Cheng’s model in anisotropic porous cracked medium: An ultrasonic physical modeling approach

Many effective medium theories are designed to describe the macroscopic properties of a medium (the rock, or reservoir in this case) in terms of the properties of its constituents (the background matrix of the rock and the inclusions, for our scenario). A very well known effective medium theory is the Eshelby-Cheng model, which was studied by us in previous work, being tested for the case where the background medium was weakly-anisotropic and porous. The analysis was done testing elastic velocities and Thomsen parameters - as a function of crack density for fixed values of aspect ratio - predicted by the model with data acquired from synthetic rock samples. In this work, we aim to complete the analysis of the Eshelby-Cheng model capabilities when applied to rocks with porous and vertical transversely isotropic (VTI) backgrounds, testing the model for the elastic velocities as functions of aspect ratio - for fixed values of crack density - against experimental data. The data used to test the model were obtained from 17 synthetic rock samples, one uncracked and 16 cracked, the latter divided into four groups of four samples each, each group with cracks having the same aspect ratio, but with the samples having different crack densities. In these samples, ultrasonic pulse transmission measurements were performed to obtain the experimental velocities used to test the model. As was not possible to acquire data for velocity as a function of aspect ratio for fixed values of crack density, we performed interpolations of the experimental data to estimate these velocities. Eshelby-Cheng model effective velocities and Thomsen parameters were calculated using three formulations proposed for the crack porosity: one proposed by Thomsen, the second one proposed in our previous work (which depends only on the crack density) and the third one proposed in this work (which depends on the crack porosity and the aspect ratio, just like Thomsen’s proposal). The comparisons between elastic velocities and Thomsen parameters - as function of crack aspect ratio, for fixed values of crack density - predicted by the model and estimated from the data via interpolation showed that the third formulation produced better fittings (lower root-mean-square errors) between model and experimental data for all ranges of aspect ratio and crack density.

Detection of structural changes in concrete using embedded ultrasonic sensors based on autoregressive model

Embedded ultrasonic transmission measurements can be a cost effective and more user-friendly alternative in comparison to commonly used structural health monitoring systems used in civil engineering to detect operational or environmental changes in structure. They can be used to detect small structural changes in large concrete structures without necessity of placing a sensor on the spot where the changing is taking place. This paper presents the investigations on the possibility of utilising autoregressive model, where the velocity of ultrasonic wave in a medium is dependent on the operational state. The goal is to use the model for localization of operational changes in the large concrete structure by means of embedded ultrasonic transducer networks. In this study, several static load tests and dynamic test on large reinforced concrete beams have been performed using embedded ultrasonic sensors. Using the autoregressive model it is possible to localize operational changes in the concrete structure. The proposed approach of diagnostic signal processing allows for precise evaluation of structural changes in concrete.

Processing Ultrasonic Data by Coda Wave Interferometry to Monitor Load Tests of Concrete Beams.

Ultrasonic transmission measurements have been used for decades to monitor concrete elements, mostly on a laboratory scale. Recently, coda wave interferometry (CWI), a technique adapted from seismology, was introduced to civil engineering experiments. It can be used to reveal subtle changes in concrete laboratory samples and even large structural elements without having a transducer directly at the place where the change is taking place. Here, several load tests until failure on large posttensioned concrete beams have been monitored using networks of embedded transducers. To detect subtle effects at the beginning of the experiments and cope with severe changes due to cracking close to failure, the coda wave interferometry procedures had to be modified to an adapted step-wise approach. Using this methodology, we were able to monitor stress distribution and localize large cracks by a relatively simple technique. Implementation of this approach on selected real structures might help to make decisions in infrastructure asset management.

High-Resolution Ultrasonic Imaging of Artworks with Seismic Interferometry for Their Conservation and Restoration

ABSTRACTArtworks are an inseparable part of the cultural heritage of societies and provide us with a unique look at cultural developments through time and space. For the best possible conservation, it is paramount to know the constituent materials, condition, and construction techniques of the objects (e.g. painting on wood, fresco, sculpture). Such information is required not only for the surfaces of the objects, but also for the interiors; in the imaging discipline, this is known as depth imaging. Here, we introduce a new method for non-invasive depth imaging as an alternative to traditional non-invasive methods when the latter cannot be used to obtain the required information. We use ultrasonic transverse-wave transmission measurements and turn them into virtual reflection measurements. We achieve this by applying seismic interferometry with active sources. Obtaining reflection measurements by seismic interferometry allows us to apply an advanced imaging technique – prestack depth migration, as used in seismic exploration – to produce a high-resolution depth image of an object. We apply our method to ultrasonic data recorded on a mockup of a painting on a wooden support. We validate our method by comparing our results with an image from X-ray computed tomography.

A program to calculate pulse transmission responses through transversely isotropic media

We provide a program (AOTI2D) to model responses of ultrasonic pulse transmission measurements through arbitrarily oriented transversely isotropic rocks. The program is built with the distributed point source method that treats the transducers as a series of point sources. The response of each point source is calculated according to the ray-tracing theory of elastic plane waves. The program could offer basic wave parameters including phase and group velocities, polarization, anisotropic reflection coefficients and directivity patterns, and model the wave fields, static wave beam, and the observed signals for pulse transmission measurements considering the material's elastic stiffnesses and orientations, sample dimensions, and the size and positions of the transmitters and the receivers. The program could be applied to exhibit the ultrasonic beam behaviors in anisotropic media, such as the skew and diffraction of ultrasonic beams, and analyze its effect on pulse transmission measurements. The program would be a useful tool to help design the experimental configuration and interpret the results of ultrasonic pulse transmission measurements through either isotropic or transversely isotropic rock samples.

Uncertainty estimation for linearised inverse problems comparing Bayesian inference and a pseudoinverse approach for acoustic transmission measurements

Abstract The determination of acoustic material parameters using ultrasonic transmission measurements can mathematically be described as an inverse problem. The question concerning the influence of uncertainties on the problem's solution can be answered using a statistical approach. Therefore, the sources of uncertainty have to be identified statistically. A method for linearising the model function using the Guide to the expression of uncertainty in measurement is introduced for inverse problems. The statistical inversion of the model is executed by means of two different approaches, that are compared and interpreted exemplarily.

Ultrasonic transmission measurements in the characterization of viscoelasticity utilizing polymeric waveguides

For the numerical simulation of acoustic wave propagation in (measurement) systems and their design, the use of reliable material models and material parameters is a central issue. Especially in polymers, acoustic material parameters cannot be evaluated based on quasistatically measured parameters, as are specified in data sheets by the manufacturers.In this work, a measurement method is presented which quantifies, for a given polymeric material sample, a complex-valued and frequency-dependent material model. A novel three-dimensional approach for modeling viscoelasticity is introduced. The material samples are designed as hollow cylindrical waveguides to account for the high damping characteristics of the polymers under test and to provide an axisymmetric structure for good performance of waveguide modeling and reproducible coupling conditions arising from the smaller coupling area in the experiment. Ultrasonic transmission measurements are carried out between the parallel faces of the sample. To account for the frequency dependency of the material properties, five different transducer pairs with ascending central frequency from to are used. After passing through the sample, each of the five received signals contains information on the material parameters which are determined in an inverse procedure. The solution of the inverse problem is carried out by iterative comparison of an innovative forward SBFEM-based simulations of the entire measurement system with the experimentally determined measurement data. For a given solution of the inverse problem, an estimate of the measurement uncertainty of each identified material parameter is calculated.Moreover, a second measurement setup, based on laser-acoustic excitation of Lamb modes in plate-shaped specimens, is presented. Using this setup, the identified material properties can be verified on samples with a varied geometry, but made from the same material.

How diffraction limits ultrasonic screening in phononic plate composed of a periodic array of resonant slits

We explore experimentally the role played by diffraction in the phenomenon of acoustic shielding provided by a plate that is periodically perforated with subwavelength slits and immersed in water. We carried out ultrasonic transmission measurements for all directions of propagation in order to check the omnidirectionality of acoustic shielding. While a single slit acts as a Fabry–Perot resonator in the frequency range of interest, the coupling between adjacent slits provides an attenuation frequency band centered around the resonant frequency that is mostly independent of the angle of incidence. Beyond the incident angle of 45 degrees, however, we observe the appearance of scattered radiation that limits the attenuation of ultrasound. This spurious scattering is shown to arise from diffraction by the grating of slits.

Ultrasonic measurements of attenuation and velocity of compressional and shear waves in partially frozen unconsolidated sediment and synthetic porous rock

The presence of partially frozen liquid in the pore spaces of porous materials has significant effects on elastic wave propagation. Although the characterization of partially frozen systems using velocity information has been well developed, application of the attenuation information is limited because the attenuation mechanisms in partially frozen systems are poorly understood. We have conducted ultrasonic wave transmission measurements with changing temperatures from 0°C to [Formula: see text] to estimate the effect of partially frozen liquids grown in unconsolidated (unconsolidated sediment) and consolidated (synthetic porous rock) materials on the velocity and attenuation of P- and S-waves. Our experimental results determined that the existence of partially frozen liquid in the unconsolidated and consolidated materials increases the velocity and attenuation for temperatures of 0°C to around the freezing point (i.e., [Formula: see text]), thus experimentally validating the unintuitive observations of high velocity and high attenuation. We interpreted the differences in velocity-versus-temperature curves as both the difference in inherent stiffness between the matrix of the consolidated material and the ice frame of the partially frozen unconsolidated material and the microscale ice distribution in pore spaces. We have also attributed the difference in the attenuation-versus-temperature curves in the unconsolidated and consolidated materials between the P- and S-waves to the difference of attenuation mechanism between the P- and S-waves. Our findings can be used for interpreting the velocity and attenuation results from the sonic logging measurements.

Setting Control of Completely Recyclable Concrete with Slag and Aluminate Cements

A completely recyclable concrete (CRC) is designed to have a chemical composition equivalent to the one of general raw materials for cement production. By doing so, this CRC can be used at the end of its service life in cement manufacturing without the need for ingredient adjustments. In one of the designed CRC compositions, blast-furnace slag cement (BFSC) was combined with calcium aluminate cement (CAC), which resulted in fast setting. In an attempt to control this fast setting, different retarders and/or the combination of lime and calcium sulfate were added to the system. The workability (slump and flow), setting time (ultrasonic transmission measurements and Vicat), strength development (compressive strength tests), and hydration behavior (isothermal calorimetry) were studied. It was found that the combined addition of lime and calcium sulfate results in a workable mixture that becomes even more workable if a retarder is also added to the system.

Wood elastic characterization from a single sample by resonant ultrasound spectroscopy

The goal of this paper is to propose an experimental method allowing the identification of the complete elastic tensor of anisotropic biological materials such as wood using only one sample. To do so, two complementary methods are used. First, the wood eigen-directions are defined from a sample of spherical shape that is then cut into a cube in a way to perform resonant ultrasound spectroscopy (RUS). The method is successfully applied on a reference beech sample with known orthotropic directions. A comparison of the identified elastic constants with those from the literature and some inferred from ultrasonic transmission measurements is given.

Calculation of Particle Size Distributions from Ultrasonic Measurements: A Mineral Slurry Case Study

AbstractUltrasonic transmission measurements have been used extensively to determine the particle size of solids in slurries. This case study examines the application of mathematical inversion techniques to the determination of the particle size distribution of a mineral slurry from data collected at a minerals processing plant. A new mathematical inversion technique, based on an extension of modified Chahine iteration combined with the principle of maximum entropy has been developed. Four algorithms were constructed and used to calculate particle size distributions from synthetic and raw plant data. These incorporated modified Chahine iteration and its extension, together with two different approaches to applying a density measurement constraint on the particle size distribution. In general the algorithms performed well with regression errors below 3 %. The correlation coefficients and slopes for this technique were 0.86 and 1.35 for the weight fraction of particles less than 75 microns when compared with the laser diffraction analysis. A better match was obtained for the plant data by using the new inversion technique, into which the principle of maximum entropy has been incorporated whereas this was not the case with the synthetic data, illustrating the need to match the inversion technique to the problem.

A poroelastic model for ultrasonic wave attenuation in partially frozen brines

Although there are many possible mechanisms for the intrinsic seismic attenuation in composite materials that include fluids, relative motion between solids and fluids during seismic wave propagation is one of the most important attenuation mechanisms. In our previous study, we conducted ultrasonic wave transmission measurements on an ice-brine coexisting system to examine the influence on ultrasonic waves of the unfrozen brine in the pore microstructure of ice. In order to elucidate the physical mechanism responsible for ultrasonic wave attenuation in the frequency range of 350–600 kHz, measured at different temperatures in partially frozen brines, we employed a poroelastic model based on the Biot theory to describe the propagation of ultrasonic waves through partially frozen brines. By assuming that the solid phase is ice and the liquid phase is the unfrozen brine, fluid properties measured by a pulsed nuclear magnetic resonance technique were used to calculate porosities at different temperatures. The computed intrinsic attenuation at 500 kHz cannot completely predict the measured attenuation results from the experimental study in an ice-brine coexisting system, which suggests that other attenuation mechanisms such as the squirt-flow mechanism and wave scattering effect should be taken into account.