The influence of gas bubbles on sediment acoustic properties: in situ, laboratory, and theoretical results from Eckernförde Bay, Baltic sea

Acoustic turbidity caused by the presence of gas bubbles in seafloor sediments is a common occurrence worldwide,but is as yet poorly understood. The Coastal Benthic Boundary Layer experiment in the Baltic off northern Germany was planned to better characterize the acoustic response of a bubbly sediment horizon. In this context, in situ measurements of compressional wave speed and attenuation were made over the frequency range of 5–400 kHz in gassy sediments of Eckernförde Bay. Dispersion of compressional speed data was used to determine the upper limit of the frequency of methane bubble resonance at between 20 and 25 kHz. These data, combined with bubble size distributions determined from CT scans of sediments in cores retained at ambient pressure, yield estimates of effective bubble sizes of 0.3–5.0 mm equivalent radius. The highly variable spatial distribution of bubble volume and bubble size distribution is used to reconcile the otherwise contradictory frequency-dependent speed and attenuation data with theory. At acoustic frequencies above resonance (>25 kHz) compressional speed is unaffected by bubbles and scattering from bubbles dominates attenuation. At frequencies below resonance (<1 kHz) ‘compressibility effects’ dominate, speed is much lower (250 m s^-1) than bubble-free sediments, and attenuation is dominated by scattering from impedance contrasts. Between 1.5 and 25 kHz bubble resonance greatly affects speed and attenuation. Compressional speed in gassy sediments (1100–1200 m s^-1) determined at 5–15 kHz is variable and higher than predicted by theory (<250 m s^-1). These higher measured speeds result from two factors: speeds are an average of lower speeds in gassy sediments and higher speeds in bubble-free sediments; and the volume of smaller-sized bubbles which contribute to the lower observed speeds is much lower than total gas volume. The frequency-dependent acoustic propagation is further complicated as the mixture of bubble sizes selectively strips energy near bubble resonance frequencies (very high attenuation) allowing lower and higher frequency energy to propagate. It was also demonstrated that acoustic characterization of gassy sediments can be used to define bubble size distribution and fractional volume.     

饱和砂岩的黏弹行为的实验研究

通过Metravib热机械分析仪,用正弦波进行加载,实验时固定静载为100 N,正弦波动载荷恒为60 N,将总载荷控制在屈服点以下;在温度为-50～125℃,升温速率为1℃/min,频率为5～1000 Hz的条件下,对泵油饱和长石砂岩、彭山砂岩样品进行单轴循环加载实验,求取泵油饱和长石砂岩和彭山砂岩的衰减、耗散角、模量...     

裂缝、裂隙介质衰减频散研究

为研究裂缝、裂隙介质中波致流引起的衰减,将裂缝看作背景孔隙岩石中非常薄且孔隙度非常高的层状介质,并等价成White周期层状模型.分别考虑不同类型的裂隙和孔隙之间的挤喷流影响,结合改进的Biot方程,推导得到裂缝裂隙介质的刚度与频率的关系.当缝隙中饱含流体时,介质的衰减和速度频散受裂缝、孔隙之间和裂隙、孔隙之间流体流动的显著影响.在低频极限下,裂缝裂隙介质的性质由各向异性Gassmann理论和挤喷流模型获得;而在非常高的频率时,由于缝隙中的压力来不及达到平衡,波致流的影响可忽略.分析表明,裂隙密度主要影响波的衰减,而裂隙纵横比主要控制优势衰减频率和速度显著变化的频率范围;由于不同裂隙的衰减机制不同,衰减和速度频散大小有所差异,但基本趋势相同.     

Experimental study on frequency-dependent elastic properties of weakly consolidated marine sandstone: effects of partial saturation

Investigating seismic dispersion and attenuation characteristics of loosely compacted marine sandstone is essential in reconciling different geophysical measurements (surface seismic, well logging and ultrasonic) for better characterization of a shallow marine sandstone reservoir. We have experimented with a typical high-porosity and high-permeability sandstone sample, extracted from the Paleogene marine depositional setting in the Gulf of Mexico, at the low-frequency band (2–500 Hz) as well as ultrasonic point (10⁶ Hz), to investigate the effects of varying saturation levels on a rock's elasticity. The results suggest that the Young's modulus of the measured sample with adsorbed moisture at laboratory conditions (room temperature, 60%–90% humidity) exhibits dispersive behaviours. The extensional attenuation can be as high as 0.038, and the peak frequency occurs around 60 Hz. The extensional attenuation due to moisture adsorption can be dramatically mitigated with the increase of confining pressure. For partial saturation status, extensional attenuation increases as increasing water saturation by 79% with respect to the measured frequencies. Additionally, the results show that extensional attenuation at the fully water-saturated situation is even smaller than that at adsorbed moisture conditions. The Gassmann–Wood model can overall capture the P-wave velocity-saturation trend of measured data at seismic frequencies, demonstrating that the partially saturated unconsolidated sandstone at the measured seismic frequency range is prone to be in the relaxed status. Nevertheless, the ultrasonic velocities are significantly higher than the Gassmann–Wood predictions, suggesting that the rocks are in the unrelaxed status at the ultrasonic frequency range. The poroelastic modelling results based on the patchy saturation model also indicate that the characteristic frequency of the partially saturated sample is likely beyond the measured seismic frequency range.     

Effects of petrophysical parameters on attenuation and dispersion of seismic waves in the simplified poroelastic theory: Comparison between the Biot's theory and viscosity‐extended theory

The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly.     

Compressional-wave Q estimation from full-waveform sonic data

There is significant evidence that the anelastic loss of seismic energy is linked to petrophysical properties such as porosity, permeability and clay content. Thus, reliable estimation of anelastic attenuation from seismic data can lead to improved methods for the prediction of petrophysical properties. This paper is concerned with methods for the estimation of attenuation at sonic frequencies (5–30 KHz) from in situ data. Two independent methods have been developed and tested for estimating compressional‐wave attenuation from full‐waveform sonic data. A well‐established technique, the logarithm spectral ratio (LSR) method, is compared with a new technique, the instantaneous frequency (IF) method. The LSR method uses the whole spectrum of the seismic pulse whilst the IF method uses a carefully estimated value of instantaneous frequency which is representative of the centre frequency of the pulse. In the former case, attenuation estimation is based on the relative variation of amplitudes at different frequencies, whilst in the latter case it is based on the shift of the centre frequency of the pulse to lower values during anelastic wave propagation. The IF method does not assume frequency independence of Q which is a necessary assumption for the LSR method, and it provides a stable frequency log, the peak instantaneous frequency (PIF) log, which may be used as an indicator for attenuation under certain limitations. The development and implementation of the two methods is aimed at minimizing the effect of secondary arrivals, such as leaky modes, and involved a series of parameter tests. Testing of the two methods using full‐waveform sonic data of variable quality, obtained from a gas‐bearing sandstone reservoir, showed that the IF method is in general more stable and suitable for full‐waveform sonic data compared with the LSR method. This was evident especially in data sets with high background noise levels and wave‐interference effects. For good quality data, the two methods gave results that showed good agreement, whilst comparison with other log types further increased confidence in the results obtained. A significant decrease (approximately 5 KHz) in the PIF values was observed in the transition from an evaporite/shale sequence to the gas‐bearing sandstone. Average Q values of 54 and 51 were obtained using good quality data from a test region within the gas‐saturated sandstone reservoir, using the LSR and IF methods, respectively.     

周期成层Patchy模型中纵波的频散和衰减研究

本文用孔隙介质力学的方法来研究纵波在介观周期成层Patchy模型的传播问题.相对于White方法,该方法少了不同介质分界面处的流体压力不等的假设,推导过程更为严密;和Dutta的双相方程解耦方法比较,文中方法直接求解双相方程,形式上更为直观.当纵波通过周期成层孔隙模型时,在低频段用孔隙弹性力学得到纵波频散和衰减结果与用White公式得到的结果符合得很好;周期成层的Patchy孔隙模型由于其空间排列的周期性使其成为孔隙弹性声子晶体,在其频散和衰减曲线上的高频段会出现多个禁带和通带,这使得速度不会随频率单调递增,衰减峰也不只一个.对于地震频段,随着渗透率降低模型的衰减峰向低频移动,这和实际观测结果一致,和传统Biot理论预测结果正相反;随着含气量上升纵波衰减峰值先上升后下降,其最大值出现在含气量0.1左右,这和相关实验结果一致.研究结果表明该模型可以定性解释观测到的地震衰减结果.     

基于常规测井数据的纵波本征衰减估计

地震波本征衰减反映了地层及其所含流体的一些特性,对油气勘探开发有重要意义.已有的理论研究与实验发现,地震频带内的衰减主要与中观尺度(波长与颗粒尺度之间)的斑状部分饱和、完全饱和岩石弹性非均匀性情况下波诱导的局部流体流有关.这种衰减与岩石骨架、孔隙度及充填流体的性质密切相关.本文着重讨论均匀流体分布、斑状或非均匀流体分布两种情况下部分饱和岩石的纵波模量差异.以经典岩石物理理论和衰减机制认识为基础,通过分析低频松弛状态、高频非松弛状态岩石的弹性模量,讨论储层参数(如孔隙度、泥质含量以及含水饱和度等)与纵波衰减之间的确定性关系.上述方法与模型在陆相砂泥岩地层与海相碳酸盐岩地层中的适用性通过常规测井资料得到了初步验证.     

An overview of laboratory apparatuses to measure seismic attenuation in reservoir rocks

Intrinsic wave attenuation at seismic frequencies is strongly dependent on rock permeability, fluid properties, and saturation. However, in order to use attenuation as an attribute to extract information on rock/fluid properties from seismic data, experimental studies on attenuation are necessary for a better understanding of physical mechanisms that are dominant at those frequencies. An appropriate laboratory methodology to measure attenuation at seismic frequencies is the forced oscillation method, but technical challenges kept this technique from being widely used. There is a need for the standardization of devices employing this method, and a comparison of existing setups is a step towards it. Here we summarize the apparatuses based on the forced oscillation method that were built in the last 30 years and were used to measure frequency‐dependent attenuation in fluid‐saturated and/or dry reservoir rocks under small strains (10^?8–10^?5). We list and discuss important technical aspects to be taken into account when working with these devices or in the course of designing a new one. We also present a summary of the attenuation measurements in reservoir rock samples performed with these apparatuses so far.     

Forced oscillation measurements of seismic attenuation in fluid saturated sandstone

Adopting the method of forced oscillation, attenuation was studied in Fontainebleau sandstone (porosity 10%, permeability 10 mD) at seismic frequencies (1–100 Hz). Confining pressures of 5, 10, and 15 MPa were chosen to simulate reservoir conditions. First, the strain effect on attenuation was investigated in the dry sample for 11 different strains across the range 1 × 10^?6–8 × 10^?6, at the confining pressure of 5 MPa. The comparison showed that a strain of at least 5 × 10^?6 is necessary to obtain a good signal to noise ratio. These results also indicate that nonlinear effects are absent for strains up to 8 × 10^?6. For all the confining pressures, attenuation in the dry rock was low, while partial (90%) and full (100%) saturation with water yielded a higher magnitude and frequency dependence of attenuation. The observed high and frequency dependent attenuation was interpreted as being caused by squirt flow.     

Wave propagation analysis of porous rocks with the thermal activated relaxation mechanism

Wave attenuation and phase velocity dispersion in the temperature domain are more complicated than those in the frequency domain. To describe wave propagation properties in the temperature domain, a so-called thermal activation mechanism model is built on the experimental result that increasing the temperature or decreasing frequency could obtain similar results on the attenuation. A rheological model (the Zener model) is employed to describe viscoelastic attenuation in saturated porous rocks. The Arrhenius relation is introduced to describe the thermal activation mechanism. The wave propagation model with thermal effects in porous media is then obtained, and 1-D P-wave and S-wave propagation characteristics are analyzed in numeric process, respectively.Two attenuation mechanisms are found in this model, the Biot loss and the thermal activation relaxation. The thermal relaxation attenuation peak and the Biot attenuation peak are observed in both frequency and temperature spectra. These two peaks move towards each other when the temperature increases on frequency spectra. The thermal relaxation peak shifts towards higher frequencies while the Biot peak shifts towards lower frequencies. At some temperature, these two peaks will superpose. The combination of the thermal relaxation and the Biot loss leads to the complexity of wave velocity curves. Similar phenomena could be observed on temperature spectra. The thermal relaxation features may relate to a so-called “local heat transfer” mechanism. These two peaks in the temperature domain have been observed in the experiments by other investigators. The characteristics of velocity and attenuation are more remarkable for high porosity rock samples. The model is helpful for the understanding of wave propagation in the temperature domain.     

Acoustic wave propagation in saturated porous media: reformulation of the Biot/Squirt flow theory

A model of wave propagation in fluid-saturated porous media is developed where the principal fluid/solid interaction mode affecting the propagation of the acoustic wave results from the conjunction of the Biot and the Squirt flow mechanism. The difference between the original Biot/Squirt (BISQ) flow theory and the new theory, which we call the reformulated BISQ, is that the average fluid pressure term appearing in the dynamic equation for a two component solid/fluid continuum is independent of squirt flow length. P-velocity and attenuation relate to measurable rock physical parameters: the Biot's poroelastic constants, porosity, permeability, pore fluid compressibility and viscosity. Modelling shows that velocity and attenuation dispersion obtained using the reformulated BISQ theory are of the same order of magnitude as those obtained using the original BISQ theory. Investigation on permeability effect on velocity and attenuation dispersion indicate that the transition zone in velocity and attenuation peak, occurring both at the relaxation frequency, shifts toward high frequency when permeability decreases. This behaviour agrees with Biot's theory prediction.     

S wave attenuation based on Stokes boundary layer

Biot theory was based on two ideas: the coupling factor to quantify the kinetic energy of fluid and Darcy permeability to quantify the dissipation function. As Biot theory did not well predict attenuation of ultrasonic S wave, we modify the theory to better characterize the S wave attenuation. The range of the coupling factor is at first estimated in view of fluid mechanics. Application of the original theory to water-saturated Boise sandstone and brine-saturated Berea sandstone shows that the model prediction significantly underestimates the S wave attenuation ultrasonically measured. For this reason, we replace Darcy permeability with variable permeability to improve the fluid momentum equation. The new model yields predictions of phase velocity and the quality factor both close to the ultrasonic measurements. The reason why the improved model is superior to Biot theory is that variable permeability is based on the Stokes boundary layer at the fluid–solid interface, thus accurately quantifying the viscous stress between the two phases. Finally, the length scale of the viscous stress is calculated in the mesoscopic sense.     

多孔岩石波传播的热弛豫模型修正

经研究发现热弛豫模型的衰减比BISQ模型大得较多,与地球介质相比衰减量也过大;反演与实验结果相比在虚模量的低频(或低温)端和高频(或高温)端相差较大,仅在峰值附近符合较好.针对上述不足,将Arrhenius关系直接引进Biot模型,替换原模型引进的峰值点的频率对数和温度关系,并重新调整了模型参数.这样既改善了原模型衰减量过大,又克服了模型反演中虚模量峰值曲线两侧差异较大的缺点.进行了P波和S波的波传播分析,仍然在频率谱和温度谱上获得热弛豫峰和Biot峰.分析显示热弛豫峰导致波速随频率升高而上升的普遍规律,Biot峰导致波速随温度升高而上升的异常现象.在相同条件下对Biot模型,BISQ模型和热弛豫模型的P波波速和衰减进行了对比.热弛豫模型得到的速度频散更强,频散范围更宽,所得的衰减峰值频率比BISQ模型要低,衰减幅度比BISQ模型稍大.这些结果与实验结果相近,更符合实际.     

Wave velocities and attenuation of shaley sandstones as a function of pore pressure and partial saturation

We obtain the wave velocities and quality factors of clay‐bearing sandstones as a function of pore pressure, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the coexistence of two solids (sand grains and clay particles) and a fluid mixture. Additional attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/fluid mixing law that satisfies the Wood and Voigt averages at low and high frequencies, respectively. Pressure effects are accounted for by using an effective stress law. By fitting a permeability model of the Kozeny– Carman type to core data, the model is able to predict wave velocity and attenuation from seismic to ultrasonic frequencies, including the effects of partial saturation. Testing of the model with laboratory data shows good agreement between predictions and measurements.     

Laboratory measurements of seismic attenuation and Young's modulus dispersion in a partially and fully water‐saturated porous sample made of sintered borosilicate glass

We measured the extensional‐mode attenuation and Young's modulus in a porous sample made of sintered borosilicate glass at microseismic to seismic frequencies (0.05–50 Hz) using the forced oscillation method. Partial saturation was achieved by water imbibition, varying the water saturation from an initial dry state up to ～99%, and by gas exsolution from an initially fully water‐saturated state down to ～99%. During forced oscillations of the sample effective stresses up to 10 MPa were applied. We observe frequency‐dependent attenuation, with a peak at 1–5 Hz, for ～99% water saturation achieved both by imbibition and by gas exsolution. The magnitude of this attenuation peak is consistently reduced with increasing fluid pressure and is largely insensitive to changes in effective stress. Similar observations have recently been attributed to wave‐induced gas exsolution–dissolution. At full water saturation, the left‐hand side of an attenuation curve, with a peak beyond the highest measured frequency, is observed at 3 MPa effective stress, while at 10 MPa effective stress the measured attenuation is negligible. This observation is consistent with wave‐induced fluid flow associated with mesoscopic compressibility contrasts in the sample's frame. These variations in compressibility could be due to fractures and/or compaction bands that formed between separate sets of forced‐oscillation experiments in response to the applied stresses. The agreement of the measured frequency‐dependent attenuation and Young's modulus with the Kramers–Kronig relations and additional data analyses indicate the good quality of the measurements. Our observations point to the complex interplay between structural and fluid heterogeneities on the measured seismic attenuation and they illustrate how these heterogeneities can facilitate the dominance of one attenuation mechanism over another.     

Velocity and attenuation in partially saturated rocks: poroelastic numerical experiments

We use a poroelastic modelling algorithm to compute numerical experiments on wave propagation in a rock sample with partial saturation using realistic fluid distribution patterns from tomography scans. Frequencies are in the range 10 to 500 kHz. The rock is a homogeneous isotropic sandstone partially filled with gas and water, which are defined by their characteristic values of viscosity, compressibility and density. We assume no mixing and that the two different pore-fills occupy different macroscopic regions. The von Kármán self-similar correlation function is used, employing different fractal parameters to model uniform and patchy fluid distributions, respectively, where effective saturation is varied in steps from full gas to full water saturation. Without resorting to additional matrix–fluid interaction mechanisms, we are able to reproduce the main features of the variation in wave velocity and attenuation with effective saturation and frequency, as those of published laboratory experiments. Furthermore, the behaviour of the attenuation peaks versus water saturation and frequency is similar to that of White's model. The conversion of primary P-wave energy into dissipating slow waves at the heterogeneities is shown to be the main mechanism for attenuating the primary wavefield. Fluid/gas patches are shown to affect attenuation more than equivalent patches in the permeability or solid-grain properties.     

Interaction of multiple courses of wave‐induced fluid flow in layered porous media

Different theoretical and laboratory studies on the propagation of elastic waves in layered hydrocarbon reservoir have shown characteristic velocity dispersion and attenuation of seismic waves. The wave‐induced fluid flow between mesoscopic‐scale heterogeneities (larger than the pore size but smaller than the predominant wavelengths) is the most important cause of attenuation for frequencies below 1 kHz. Most studies on mesoscopic wave‐induced fluid flow in the seismic frequency band are based on the representative elementary volume, which does not consider interaction of fluid flow due to the symmetrical structure of representative elementary volume. However, in strongly heterogeneous media with unsymmetrical structures, different courses of wave‐induced fluid flow may lead to the interaction of the fluid flux in the seismic band; this has not yet been explored. This paper analyses the interaction of different courses of wave‐induced fluid flow in layered porous media. We apply a one‐dimensional finite‐element numerical creep test based on Biot's theory of consolidation to obtain the fluid flux in the frequency domain. The characteristic frequency of the fluid flux and the strain rate tensor are introduced to characterise the interaction of different courses of fluid flux. We also compare the behaviours of characteristic frequencies and the strain rate tensor on two scales: the local scale and the global scale. It is shown that, at the local scale, the interaction between different courses of fluid flux is a dynamic process, and the weak fluid flux and corresponding characteristic frequencies contain detailed information about the interaction of the fluid flux. At the global scale, the averaged strain rate tensor can facilitate the identification of the interaction degree of the fluid flux for the porous medium with a random distribution of mesoscopic heterogeneities, and the characteristic frequency of the fluid flux is potentially related to that of the peak attenuation. The results are helpful for the prediction of the distribution of oil–gas patches based on the statistical properties of phase velocities and attenuation in layered porous media with random disorder.     

Stratigraphic absorption compensation based on multiscale shearlet transform

Seismic waves propagating through viscoelastic media experience stratigraphic absorption and attenuation effects, which directly affect the imaging resolution in seismic exploration. Without stratigraphic absorption, the ratio of deep reflection energy to shallow reflection energy (attenuation ratio) is invariable at different frequencies. If a seismogram is decomposed into different frequency bands, these signals will show similar time–energy distributions. Therefore, the attenuation ratios should be similar across different frequency bands, except for frequency-variable weights. Nevertheless, the frequency-variable weights for different frequency bands can be obtained by benchmarking against the time–energy distributions of low-frequency information because the loss of low-frequency information is relatively insignificant. In this light, we obtained frequency-variable weights for different frequencies and established a stratal absorption compensation (SAC) model. The anisotropic basis of the shearlet enables nearly optimal representation of curved-shape seismic signals, and shearlets at different scales can represent signals for different frequency bands. Then, we combined the SAC model with the shearlet transform and established the new compensation method. As the signal and noise have different distributions in the shearlet domain, we selectively compensated the signals using a thresholding algorithm. Hence, it was possible to avoid noise enhancement. This is the prominent advantage of the proposed method over other compensation methods.     

Relationship between the transition frequency of local fluid flow and the peak frequency of attenuation

Local fluid flow (LFF) at the mesoscopic scale is the main dissipation mechanism of seismic waves in heterogeneous porous media within the seismic frequency band. LFF is easily influenced by the structure and boundary conditions of the porous media, which leads to different behaviors of the peak frequency of attenuation. The associated transition frequency can provide detailed information about the trend of LFF; therefore, research on the transition frequency of LFF and its relationship with the peak frequency of the corresponding attenuation (i.e., inverse of quality factor) facilitates the detailed understanding of the effect of inner structures and boundary conditions in porous media. In this study, we firstly obtain the transition frequency of fluid flux based on Biot’s theory of poroelasticity and the fast Fourier transform algorithm in a sample containing one repeating unit cell (RUC). We then analyze changes of these two frequencies in porous media with different porous properties. Finally, we extend our analysis to the influence of the undrained boundary condition on the transition frequency and peak frequency in porous media with multiple RUCs. This setup can facilitate the understanding of the effect from the undrained boundary condition. Results demonstrate that these two frequencies have the same trend at low water saturation, but amplitude variations differ between the frequencies as the amount of saturation increases. However, for cases of high water saturation, both the trend and the amplitude variation of these two frequencies fit well with each other.