Abstract

Revealing the elastic wave properties of carbonate rocks with complex pore structures and improving the reliability of carbonate reservoir descriptions have always been a global challenge in the field of carbonate geophysical exploration. In this study, we established a synthetic borehole model by selecting different particle sizes of cement, carbonate cuttings, and micro-silicon as the matrix, and silicon disks as the pores in carbonate rocks. We conducted four sets of low-porosity (0-3%) borehole models with different pore aspect ratios (ARs) and measured the P- and S-wave velocities ( and ) at the well-logging scale obtained from an acoustic logging system with one source and two receivers. The results indicate that the relationship between velocities and porosities in these borehole models follows a linear relation, with the pore AR significantly influencing the velocities at any given porosity. The velocity variation caused by pore AR reaches 560 m/s and 410 m/s at 3% porosity for the P-wave and S-wave within the AR range of 0.017-0.13. The theoretical DEM models provide a high and broad estimation of and at the well-logging scale in our measurement. They could perform better in fractured formation than in dissolved porous formation in carbonate reservoirs. The linear relation of and is independent of the pore AR and is effective for both fractured and dissolved porous formations. The change of in different pore AR is more responsive to porosity and nonlinear dependent on the pore AR. The relationship between the defined normalized and indicates the pore AR has a more significant effect on than in our model. The constructed borehole models provide a unique opportunity for evaluating the availability of rock physics models at an acoustic logging scale. The study’s findings have significant implications for improving the reliability of carbonate reservoir descriptions and enhancing the accuracy of geophysical exploration in carbonate rocks with complex pore structures.

1. Introduction

Accurate characterization of pore structure in unconventional oil and gas systems has proven to be challenging due to their complex morphology and the strong diagenesis and tectonism [13]. Currently, techniques such as acoustic emission experiments, acoustic logging, and seismic monitoring, which are based on elastic wave theory, have become important methods for detecting underground structures, oil and gas minerals, and geothermal resources. However, the changes in pore structure can significantly impact the elastic frame moduli of naturally occurring rocks, resulting in changes in their elastic wave velocities [49].

The pore aspect ratio (AR; the ratio of the minor axis to the major axis) dependence of velocities of crustal rocks, affecting the energy exchange between the formation and pore fluids, has been widely reported [5, 1013]. The power relationship between elastic wave velocities and pore AR has been proposed in the theoretical model (e.g., [14, 15]) and rock physics experiments at the core scale (e.g., [4, 8]). However, the frequency dependence of velocities can lead to errors and difficulties in the dynamic elastic parameters of rocks when attempting to apply the laboratory results in the ultrasonic frequency band (0.1 MHz-10 MHz) to the acoustic well-logging interpretation in the medium-frequency range (20 Hz-200 kHz) [1619] and seismic exploration (~100 Hz) [20]. It has always been a research hotspot in the field of geophysics to establish the relationship between the elastic wave properties of rocks across frequency bands. Additionally, considering the complex well-logging environment, it is challenging and inaccurate to obtain detailed information about the pore structure beyond porosity using conventional well-logging data [21].

Developing rock physics models with controlled pore structures is an urgent need to quantitatively investigate the effects of pore structures on wave velocity and attenuation (e.g., [4, 2225]). The artificial sandstone cores are relatively mature and simple after decades of development. Usually, quartz sand with different grain sizes and cementing agents (such as Plexiglas and epoxy resin) are mixed to simulate the sandstone. Aluminum, copper, and tin foil are the common materials used to simulate the fractures, interparticle pores, and vugs in the formation [4, 26, 27]. However, for carbonate rocks, the existing methods (e.g., sheet combination method, cutting simulation method, interpolation method, and dissolution method) are more applicable for small-scale artificial cores for the following reasons: (1) the epoxy resin is too expensive to be widely used in well-logging scale models. (2) The consistency after adding epoxy resin is too high to ensure uniform mixing. (3) The final properties of rock physics models are influenced by the times and pressures in the compaction process, while it is difficult to control these parameters in the large-scale models. (4) ccurately controlling pore parameters and distribution in the large-scale models remains difficult with existing methods [28, 29].

To simulate the acoustic well logging in the carbonate reservoir in the laboratory, we have developed a new borehole model with controlled pore structures and examined the influences of pore characteristics on the elastic properties via a homemade acoustic logging system [12]. In this study, we further focus on the pore AR dependence of velocities on a larger set of borehole models, including one blank model without pores and 20 models with different pore ARs and porosities. This work aims to provide a new dataset of measured P- and S-wave velocities ( and ) at the well-logging scale, which can serve as a preliminary reference for acoustic well-logging interpretation in the carbonate reservoir.

2. Methods

2.1. Borehole Model Building

Cement has been widely used in rock mechanics experiments for a long time [30, 31]. Both theoretical and experimental studies have shown that particle gradation determines the cement properties (e.g., density and strength) under the same mass ratio of cement and other matrix materials [32]. To consolidate the model matrix tightly and closely mimic the properties of the target formation, we choose cement slurry, carbonate cuttings, and microsilicon powder with different grain sizes as the solid phases for the borehole model. The carbonate cuttings are clean crystalline grains with a grain size of 150 μm, and the pure calcite content is over 98%. The microsilicon with 1 μm particle size is used to effectively prevent the sedimentation during the curing process. The microsilicon powder has a pure SiO2 content of over 97%.

The particle size distribution of the above three solid-phase materials are measured by the Bettersize2000 laser particle size distribution instrument (Figure 1). The average particle sizes of carbonate cuttings, cement, and microsilicon powder are 142.15 μm, 15.42 μm, and 0.76 μm, respectively. However, the addition of carbonate cuttings and microsilicon powder makes the cement slurry system less fluid and challenging to mix uniformly. The effect of different ratios of cement slurry and carbonate cuttings on the matrix fluidity is shown in Table 1. After considering the matrix density and fluidity of the cement slurry, the mass ratio of carbonate cuttings, cement, and water is determined as 1 : 1 : 0.4. The mass radio of microsilicon is 1% of the matrix material. In addition, graft polymers of sulfonated aldehyde-ketone condensation as the drag reducer and polymers of polyether and organic silicon as the defoamer are added to the matrix to discharge the bubbles formed during agitation. Finally, the required amount of materials for each model is list in Table 2.

In this study, “penny-shaped” silicone disks with different thicknesses and diameters, whose properties are comparable to fluids with relatively low density and propagation velocity, are distributed in the borehole model to simulate the different pore structures in carbonate rocks. A vertical hole with 76 mm width and 600 mm depth was drilled in the center of the model and filled with water to simulate the petroleum drilling procedure. Figure 2 shows the construction process of the borehole model.

The cores were left to cure for 28 days at room temperature (Figure 3(a)) and then examined under a metallurgical microscope (Figures 3(b) and 3(c)). The shiny areas are the carbonate cutting, the dark areas are the cement hydration products, and the black areas are the primary pores resulting from the undischarged bubbles. The porosity and permeability were obtained from the gas measurement. The primary porosity is 24.36%, the permeability is low (<0.029 mD), and the density is 2.12 g/cm3. The density difference is less than 0.05 g/cm3, indicating that the matrix properties are more similar to the actual carbonate rocks than those of the simple cement system. Finally, we constructed four sets of low-porosity (0.6-3%) borehole models with varying pore AR of 0.017, 0.033, 0.067, and 0.13 (Table 3).

2.2. Velocity Measurements

To simulate the acoustic well-logging in the laboratory, we cooperated with Yangzhou Oriental Ultrasound Technology Co. Ltd. to design an acoustic logging system with one ring-emitting piezoelectric source and two piezoelectric receivers (as shown in Figure 4). The dominant frequency of the piezoelectric source is 20-40 kHz. This equipment is placed in the central borehole filled with water. The source-receiver distances are 22.2 cm and 36.5 cm for the two receivers. The typical waveforms received by the receivers are shown in Figure 5. The velocities in the borehole model are calculated from the differential arrival times of the P- and S-waves between the two receivers: where is the velocity of th phases (P-, S-, and Stoneley waves) in km/s, is the length between the two receivers in mm, while and are the picked th arrival times of all phases waves of the two received signals in . The main source of errors in the velocity measurement is associated with the picking of travel times from the acquired raw waveforms.

3. Results

Porosity is the key parameter for characterizing materials microstructures and is the main output from well-log interpretation [33]. A broadly accepted consensus is that the complex pore structures cause the wide scatter in the velocity-porosity relationship for the carbonate rocks [4, 8, 13]. The measurements from our models show that both the and decrease with the increase in porosity at the same pore AR (Figure 6). The relationship can be well-described by a simple negative linear function, which is similar to the results from rock physics experiments [34, 35].

In addition, pore AR has a remarkable effect on elastic wave propagation [3638], affecting the energy exchange between the formation and pore fluids. However, recent researches only focus on building quantitative trends between the velocity and pore type at the core scale [8]. We built four sets of borehole models with different pore AR of 0.017, 0.033, 0.067, and 0.13. The velocity variation caused by pore AR is as high as 560 m/s and 410 m/s at 3% porosity for the P-wave and S-wave within the AR range of 0.017-0.13 in our experiments. It should be noted that more silicon disks are required to achieve the same porosity when pore AR is low in each model set.

The ratio of microstructure size () to wavelength () affects velocity in the porous media. In the past decades, several studies have tested effective medium theories (EMT) in the long-wavelength limit situation (), seeking to relate the rock minerals and pore structures to the effective elastic properties on a macroscopic scale by assuming different pore ARs [35, 3941]. In this study, the dominant frequency of the acoustic logging system is 20-40 kHz, and the wavelength is 50-100 mm. The size of the silicon disks is 30 mm; thus, our model is within the scope of the long-wavelength range. Therefore, we analyze the pore AR dependence of velocities in the following sections and compare it with those of EMT models to gain further insight into the phenomena that control elastic wave propagation in the borehole models. The previous study has indicated that the differential effective medium (DEM) model [42] shows closer outcomes to our results in the borehole models [12]. Therefore, we use the DEM model to calculate the theoretical velocities of the borehole model with different pore ARs. The values of the input parameters in the DEM model are summarized in Table 4. Other parameters used in these theoretical models are the same as those of the borehole model. The results show that the theoretical DEM models provide a high estimation of the and , and the deviations become larger as pore AR increases. It is worth mentioning that relatively small variations of wave velocities are observed at the well-logging scale. Therefore, the theoretical model is more suitable for fractured formation than the dissolved porous formation in carbonate reservoirs.

4. Discussion

4.1. Linear Relation

The relation between the and for sandstone is usually constant. Our newly developed carbonate borehole models follow a linear relation (Figure 7), which has been verified from the test of carbonate samples [43]. The high determinant coefficient ( value) indicates that the linear relation of and is independent of the pore AR.

Wang et al. [24] introduced a general velocity model for carbonate rocks determined by core measurements: where , , , , , , , and are empirical parameters; represents the primary porosity; represents the secondary porosity; and represent the P- and S-wave velocities of the solid matrix, respectively.

Then, Equation (1) can be further rewritten as follows:

Based on the measured data, the fitting parameters and range from -1 to 0; and range from 0.5 to 1. In our study, the secondary porosity is low (<3%). Therefore, the , related to the secondary porosity, is close to 1. Then, Equation (2) could be simplified as follows:

The , , and are related to the solid matrix and constant in the same model. In addition, the fitting parameters in Equation (4) (i.e., , , , and ) are only related to porosity and independent of the pore AR, which means that the linear relation between and is effective for different pore AR conditions. Combined with the results of the same linear relation between and in our borehole models, the method for prediction from well-logging data can be used for carbonate formations with fractures, intraparticle pores, or vugs.

4.2. Ratio

Seismic tomographic results have revealed the importance of as a diagnostic of the physical cause of a seismic anomaly [44, 45]. Although geophysicists have widely recognized the importance of liquid compressibility in changing [46], the effect of pore geometry on has not been commonly acknowledged. In sandstones, the pore structure is simple and is relatively stable. According to classical elastic theory, is approximately equal to 1.732 when the Poisson’s ratio is 0.25. However, caution is needed when applying the suitability of this relation in carbonate rocks due to the variability of pore structures. To quantitatively describe the variation with pore AR, we calculate the change in as

As shown in Figure 8, the is always lower than 0 and has a positive response to porosity, while the influence of pore AR appears to be nonlinear. Due to the low number of models, the accurate relation between pore AR and the cannot be proposed in this study. Therefore, we observe the same negative trend from the DEM modelling. The maximum change in appears at model set C (i.e., ). Overall, the is more responsive to porosity and nonlinear dependent on the pore AR, which indicates a complex effect of pore AR on in carbonate rocks.

4.3. Normalized Velocity

The pore AR influences and differently based on the above analysis. Therefore, we compare the velocities within the same pore AR window. The lower and upper pore AR bounds for the velocities analysis are 0.017 and 0.13, and the normalized velocities ( and ), which quantitatively represent the pore AR dependence of velocities, are defined as

The and usually increase with pore AR and follow a power law; thus, the and are higher than 1 (Figure 9). The and both increase by the porosity. The relationship between and follows a linear function where the value is as high as 0.97. is always higher than with a slope of 1.326, which means the pore AR has a greater effect on than in our model. Interestingly, the calculation from the DEM model shows an astonishing similarity with the trends from our results. Overall, the results suggest that the observed and relationship for carbonate rocks is influenced by the pore structure.

Like any artificial core experiment, the results obtained from our on ideal borehole models have their share of pitfalls in view of the isolated “penny-shaped” pores. First, the P-wave velocity and density of silicon disks are similar to those of water, while the of this material is in low comparability to fluid typically encountered in crustal materials. Second, our borehole model ignores the interaction between the solid matrix and the pores/pore fluids, which is also crucial for their elastic wave velocities [47, 48]. However, in this study, we only focus on the effect of pore structures on velocities. The synthetic borehole models are simplified but similar to the actual case, and the elasticity contrasts between the matrix and pore fluids. Nevertheless, our borehole model is suitable for various petrophysical experiments such as acoustic wave characteristics, electrical characteristics, and seepage mechanisms. The results from our borehole model could be a stepping stone for understanding the velocity variation at lower frequencies compared to the conventional rock physics model.

5. Conclusions

By choosing different particle sizes of cement, carbonate cuttings, and microsilicon as the matrix and silicon disks with different AR as the pores in carbonate rocks, we built four sets of low-porosity borehole models and provided a new dataset of the and at lower frequencies obtained from an acoustic logging system with one source and two receivers. This study focuses on the pore AR dependence of velocities and reaches the following key conclusions: (1)The relation between velocities and porosities measured at the well-logging scale follows a linear function. Meanwhile, the velocity variation caused by pore AR can reach 560 m/s and 410 m/s at 3% porosity for the P-wave and S-wave within the AR range of 0.017-0.13 in our experiments(2)The theoretical DEM models overestimate the and . The deviations between the theoretical velocities at the rock physics scale and the measured results at the well-logging scale gradually increase with AR increase, especially when AR is over 0.05, indicating better performance in fractured formations than dissolved porous formations in carbonate reservoirs(3)The high determinant coefficients indicate that the linear relation of and is independent of pore AR and is effective for both fractured or dissolved porous formation. The change of in different pore AR is more responsive to porosity and nonlinear dependent on the pore AR. The relationship between the defined normalized and indicates that the pore AR has a greater effect on than

Conflicts of Interest

The authors declare that they have no conflicts of interests.

Acknowledgments

This work was jointly funded by the National Natural Science Foundation of China (Grant nos. 42204058 and 42104071), the Science Foundation of Chongqing (Grant no. CSTB2022NSCQ-MSX1660), and the National Nonprofit Institute Research Grant of IGGE (Grant no. AS2022J03).