Next Article in Journal
Analysis of Physiological, Physical, and Tactical Responses in Small-Sided Games in Women’s Soccer: The Effect of Numerical Superiority
Previous Article in Journal
A Review of Magnetoelectric Composites Based on ZnO Nanostructures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Background Porosity on Seismic Anisotropy in Fractured Rocks: An Experimental Study

1
State Key Laboratory of Petroleum Resource and Prospecting, China University of Petroleum (Beijing), Beijing 102249, China
2
CNPC Key Laboratory of Geophysical Exploration, China University of Petroleum (Beijing), Beijing 102249, China
3
Chinese Academy of Geological Sciences, Beijing 100037, China
4
Bureau of Geophysical Prospecting Inc., China National Petroleum Corporation, Zhuozhou 072750, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(14), 8379; https://doi.org/10.3390/app13148379
Submission received: 6 June 2023 / Revised: 14 July 2023 / Accepted: 18 July 2023 / Published: 20 July 2023

Abstract

:
Fractures are widely distributed in the subsurface and are crucial for hydrocarbon, CCS, offshore infrastructure (windfarms), and geothermal seismic surveys. Seismic anisotropy has been widely used to characterize fractures and has been shown to be sensitive to background matrix porosities in theoretical studies. An understanding of the effects of background porosity on seismic anisotropy could improve seismic characterization in different fractured reservoirs. Based on synthetic rocks with controlled fractures, we conducted laboratory experiments to investigate the influence that background porosity has on P-wave anisotropy and shear wave splitting. A set of rocks containing the same fracture density (0.06) with varying porosities of 15.3%, 22.1%, 26.1% and 30.8% were constructed. The P- and S-wave velocities were measured at 0.5 MHz as the rocks were water saturated. The results show that when porosity increased from 15.3% to 22.1%, P-wave anisotropy and shear wave splitting exhibited slight fluctuations. However, when porosity continued to increase to 30.8%, P-wave anisotropy declined sharply, whereas shear wave splitting stayed nearly constant. The measured results were compared with predictions from equivalent medium theories. Qualitative agreements were found between the theoretical predictions and the measured results. In the Eshelby–Cheng model, an increase in porosity reduces fracture-induced perturbation in the normal direction of the fracture, resulting in lower P-wave anisotropy. In the Gurevich model, an increase in porosity can reduce the compressional stiffness in parallel directions to a larger extent than that in perpendicular directions, thus leading to lower P-wave anisotropy.

1. Introduction

Fractures are found throughout crustal rocks, and fracture detection is of critical importance for hydrocarbon exploration, water resource management, CO2 capture and storage, and geo-thermal exploitation. Aligned fractures can cause seismic anisotropy, making it possible to obtain fracture information from seismic data. Fracture-induced anisotropy is an integrated result of various fracture parameters (fracture density, scale, aspect ratio, etc.), fluid properties (fluid modulus and viscosities), and background properties (modulus, porosity, and permeability). Background matrix porosity has been thought to greatly influence P-wave velocity anisotropy [1,2]. Because fractures are present in low, mid, and high porosity hydrocarbon reservoirs, understanding the effect of background porosity on seismic anisotropy could further improve seismic fracture characterization in reservoirs of different porosities.
To date, theoretical and experimental studies have been conducted to relate the physical properties of rock to seismic velocity and anisotropy, and various equivalent medium theories have been proposed [1,3,4,5,6]. These theories were based on different assumptions and were derived using different methods, resulting in different predictions. Thus, there is a desire to verify and calibrate these theoretical models through laboratory experiments. Laboratory experiments typically use synthetic samples with controlled and known rock physical properties and fracture parameters. A number of experimental studies have been conducted to observe elastic wave velocity and anisotropy in fractured rocks [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. These studies have investigated how parameters, such as fracture density, scale, and aspect ratio, can influence elastic wave velocity and anisotropy, providing significant insights into the physical mechanisms of elastic wave propagation in fractured rocks. However, natural rocks have a complex depositional and diagenetic/tectonic history, and fracture parameters (fracture geometry and density) cannot be quantitatively controlled, making them too variable for experimental investigations concerning the effect of a single variable on seismic wave propagation. Thus, synthetic samples must be used in laboratory experiments to observe seismic wave propagation in cracked rocks. Despite variations in the background porosity of fractured reservoirs, few experimental studies, both in laboratory and field observations, have focused on the effect of background porosity on seismic anisotropy in these reservoirs.
Gurevich (2003) discussed P-wave anisotropy variation with background porosity in water-saturated rocks with aligned fractures in a low-frequency range using numerical simulations [1]. When there is no background porosity, P-wave anisotropy (ε) is very small. In the absence of background porosity, fluid is ‘isolated’ in fractures. The ‘isolated’ fluid acts to ‘stiffen’ the fractures and reduce fracture compliance. This results in a relatively low P-wave anisotropy for the overall fractured rock. As the characteristic porosity increases, there is a sharp increase in P-wave anisotropy. When surrounded by pores, the fluid in these fractures can escape during compression; therefore, the fracture is more compliant than it would be in ‘isolated’ conditions. As porosity continues to increase, P-wave anisotropy flattens and decreases gradually. Other theoretical studies have also discussed the influence of porosity on anisotropic parameters [2,23]. However, no experimental studies have been conducted to test and verify this, owing to the difficulties faced in manufacturing artificial fractured rocks with low porosity.
This study experimentally examined the influence of background porosity on P-wave anisotropy in fractured rocks under water saturation conditions. A set of artificial fractured sandstone samples with the same fracture density (0.06), but different porosities (ranging from 15.3% to 30.8%), were built. Laboratory ultrasonic measurements were conducted to obtain the P- and S-wave velocities and anisotropy. We analysed the variations in P-wave anisotropy and shear wave splitting (SWS) with background porosity. The results were compared with theoretical models and the causes of the variations in P-wave anisotropy and shear wave splitting with porosity were discussed.

2. Equivalent Medium Theories

2.1. High Frequency (Unrelaxed) Moduli

High frequency theories have modelled fractures as “isolated” inclusions. Thus, there was no fluid exchange between inclusions and the background medium. They assumed that if the frequency was higher than the squirt characteristic frequency, f f c , the fluid pressure does not have time to equilibrate between stiff and compliant pores during the half-wave cycle, or the so-called “unrelaxed” state. Cheng (1978, 1993) derived the elastic stiffness tensor of a fractured medium under high frequency (unrelaxed) conditions based on the Eshelby inclusion theory. Mathematically, the Eshelby–Cheng model is given by the following [5,6]:
C i j e f f = C i j ( 0 ) φ f C i j ( 1 ) ,
where C i j ( 0 ) is the background matrix stiffness tensor and φ f is the fracture porosity, and
C i j ( 1 ) = C i k ( 0 ) A k j ,
where A i j relates the applied strain e i A to the “stress-free” strain e i T of the inclusion:
e i T = A i j e j A

2.2. Low-Frequency (Relaxed) Moduli

In contrast to high-frequency theories, low-frequency theories consider the fluid exchange between fractures and the surrounding pores and micro-cracks. They assume that if the frequency is lower than the squirt characteristic frequency, the fluid pressure can equilibrate locally between the cracks and adjacent pores. By combining linear slip theory and the anisotropic Gassmann equation [24], Gurevich (2003) derived expressions for a fluid-saturated cracked porous medium under low frequency conditions in the following form:
C i j s a t = C i j 0 + α i α j M i , j = 1 , 6
where C i j 0 is the dry modulus of the fractured medium and α m is Biot’s coefficient, which is as follows:
α m = 1 Σ n = 1 3 C m n 0 3 K g
For i = 1, 2, and 3, α 4 = α 5 = α 6 = 0 . Scalar M is the direct analog of Gassmann’s pore space modulus, which is as follows:
M = K g ( 1 K K g ) + φ ( K g K f 1 )
where K g and K f is the bulk modulus of mineral grains and pore fluid, respectively. φ is the overall porosity of the fractured rock. K denotes the so-called generalized drained bulk modulus, which is defined as:
K = 1 9 i = 1 3 j = 1 3 C i j 0

3. Sample Preparation

The method of Ding was adapted to build artificial porous sandstones with similar mineral components [17,18], pore structures, and cementation to natural rocks. In this method, powders of silica sand, feldspar, and kaolinite, were mixed in a ball mill for 24 h to ensure homogeneity. Then, these mineral powders were mixed with sodium silicate solution. The mixture was then paved into a mould layer by layer. In order to make the “meso-scale” fractures, touch paper discs were spread on the surface of the mineral powder mixtures when layering the mixture in the mould. Touch paper is a special kind of paper which has been soaked in saltpetre. At high temperatures, touch paper decomposes into gas and leaves nearly no remains.
In order to pave the mixture in the mould in layers, a special tool was made (Figure 1). The tool consisted of a frame, a glass plate, and a threaded rod. By rotating the threaded rod, the glass plate can be moved vertically. The procedure of paving the mixture is shown in Figure 2. Afterwards, a predetermined number of discs were carefully placed on top of each mixture layer. Then, the next mixture layer was paved and the same number of discs were placed on it (except for the top layer).
After this, a uniaxial compressing pressure was applied to the mould for 10 min in order to ensure that the mineral grains were compact. The compressing direction was perpendicular to the fractures and layers. During this process, the pressure was adjusted to control the porosity of the samples. During the next 10 h, the sodium silicate in the mixture consolidated. In this way, a block was formed and the consolidated sodium silicate gave the block initial mechanical strength. Next, the block was demoulded and sintered in a muffle oven at 900 °C. During sintering, the touch paper discs decomposed, leaving penny-shaped voids as fractures. The final step was to polish the sample surface for the convenience of obtaining ultrasonic velocity measurements. The size of all the cuboid-shaped samples was 69.9 ± 0.1 mm in the two directions parallel to the fractures and 49.9 ± 0.1 mm in the direction perpendicular to the fractures.
In order to produce artificial rocks with the same fracture densities but different porosities, we made the following modifications to the method of Ding (2014) [17,18]: (1) Different amounts of mineral material mixture were layered in the mould (70 mm × 70 mm) in 37 layers when building the different porosity samples. (2) The mixture of different samples was compressed into the same thickness (50 ± 0.3 mm) by adjusting the applied pressure to form blocks of the same volume. Therefore, the blocks that were made using more mixture would experience higher loading pressure and would have lower porosity. (3) A total of 126 touch paper discs were spread on the surface of each layer (except the top layer) in all the blocks. Fracture density is unitless and is calculated using the following formula:
e = N a 3 V
where N is the number of fractures, a is the radius of the fractures, and V is the total volume of the sample. Thus, all the blocks had the same volume (70 mm × 70 mm × 50 mm) and contained the same number of fractures (36 layers × 126 fractures/layer), making the fracture density practically identical for the different porosities.
For each fractured sample, a blank sample (containing no penny-shaped fractures) was made as a reference sample to provide background properties. The blank samples were created using the same method as the fractured samples, excluding the introduction of fractures.
Four groups of artificial rocks with different porosities were constructed (Figure 3a). Each group consisted of two samples (one fractured sample containing penny shaped fractures with a fracture density of 0.06 and a corresponding unfractured sample). The fracture diameter was set at 3 mm and the fracture thickness was 0.06 mm. Table 1 lists the main parameters of all the samples. The porosity was measured using helium porosimetry. In the fractured samples, the fracture porosity was calculated to be 0.75%, but the porosity difference between the fractured and blank (unfractured) samples within each group varied from 0.4% to 1.5% larger or smaller than the fracture porosity, indicating a small difference in background porosity between the fractured and blank samples. This was caused by tiny errors in the manufacturing process. A section obtained by cutting a fractured sample is shown in Figure 3b, where parallelly distributed fractures can be observed. A scanning electron microscope (SEM) image of the section is shown in Figure 3c,d, where pore and fracture structures with controlled geometry are clearly shown.

4. Ultrasonic Measurement

The P- and S-wave velocity measurements were conducted on an ultrasonic bench-top pulse transmission system at room temperature and an atmosphere pressure suitable for water-saturated conditions. The ultrasonic transducers, which were heavily damped with a central frequency of 500 kHz, were excited by a voltage spike, giving rise to broadband ultrasonic pulses. The transmitted P- and S-wave signals were recorded with a digital oscilloscope and a desktop computer. The time sampling interval for all the experiments was 0.04 μs for both P- and S-wave signals. The measurement error was about 0.5% for P-wave velocity and 0.8% for S-wave velocity. Water saturation was obtained by immersing the samples in a water-filled container. The container was placed in a sealed bin to extract the air. To ensure full water saturation, the saturation rate was calculated using the following formula:
S w ( w a t e r ) = m s a t m d r y V φ ρ f × 100 %
where m s a t and m d r y is the wet and dry mass of the sample, V and φ represent the volume and helium porosimetry measured porosity of the sample, respectively, and ρ f is the bulk density of water. The P-wave velocities were measured both parallel (90°) and perpendicular (0°) to the fractures. The S-wave velocities were measured parallel to the fractures, and polarization was performed in two directions (parallel (S1 wave) and perpendicular (S2 wave) to the fractures) by rotating the transducers. Then, the Thomsen parameter (ε), which represents the difference in P-wave velocity between the parallel and perpendicular directions, was calculated using the following formula:
ε ( % ) = 100 × ( V P V P V P )
where V P and V P represent the P-wave velocities in parallel and perpendicular directions, respectively. Shear wave splitting (SWS) in the parallel direction can be calculated as follows:
S W S ( % ) = 100 × ( V s 1 V s 2 V s 2 )
where V s 1 and V s 2 are fast and slow shear wave velocities. The numerical modelling experiment of wave propagation in anisotropic media by Dellinger and Vernik (1994) showed that if the wave front was propagating parallel or perpendicular to the layering (or in this case, to the fractures) [25], a true phase velocity was measured in laboratory ultrasonic experiments. Therefore, in this study, P- and S-wave phase velocities were measured in both parallel and perpendicular directions. The transmitted signals of the blank samples were analysed using Fourier transform to obtain the dominant frequency and to calculate the wavelength (Figure 4). The ratio of wavelength to fracture diameter (λ/d ratio) was between 5.2 and 10.8 for the P wave and between 3.9 and 12 for the S wave. Based on experimental and numerical studies [21], for this λ/d ratio, both the P and S wave may exhibit Rayleigh scattering. This is undesirable because the equivalent medium criteria are violated. This problem has been frequently encountered in previous experimental studies [7,11]. However, similar to previous studies, we also found some interesting correlations between the experimental results and theoretical predictions.

5. Results and Discussions

5.1. Laboratory Results

Table 2 shows the P- and S-wave velocities of the blank samples. The P-wave velocity in the parallel direction was slightly higher than that in the perpendicular direction, and the S1 velocity (polarized parallel to the fractures) was slightly higher than the S2 velocity (polarized perpendicular to the fractures). This was due to layering and uniaxial loading during the manufacturing process.
Figure 5 shows the measured P- and S-wave velocities in the water-saturated fractured samples. Significant velocity anisotropy was evident in all the fractured samples, with the faster P-wave being measured parallel and the slower one perpendicular. The measured velocity anisotropy is discussed further in the following paragraphs. When the porosity increased from 15.3% to 22.1%, the P-wave velocities (both parallel and perpendicular) dropped sharply. However, as the porosity continued to increase, the decreasing trend in the perpendicular direction was slower than that in the parallel direction, narrowing the difference between the parallel and perpendicular directions and indicating a decrease in anisotropy. However, when porosity increased from 15.1% to 30.8%, the decrease observed in both S1 and S2 wave velocities was generally similar.
Figure 6 shows the P-wave anisotropy parameter (ε) and SWS in all the blank samples. Both the ε parameter and SWS decreased with increasing porosity in the blank samples. Because the higher porosity samples suffered lower uniaxial loading (Table 1), they exhibited weaker anisotropy.
We aimed to investigate the effect of background porosity on fracture-induced P-wave anisotropy and SWS. Because the blank samples exhibited transverse isotropy, we inferred that the background matrix in the fractured samples had the same anisotropy as the blank samples. In fractured samples, the measured anisotropy was a superposition of background matrix anisotropy and fracture-induced anisotropy. To model the anisotropy caused by fractures, the anisotropy of the background matrix had to be considered. Therefore, we approximately represented the fracture-induced P-wave anisotropy and SWS by subtracting the measured ε parameter and SWS of the blank samples from that of the fractured ones. Figure 7 depicts the measured results of the fracture-induced ε parameter and SWS, as well as the theoretical predictions. When porosity increased from 15.3% to 22.1%, both the measured ε parameter and SWS exhibited small fluctuations. In this porosity range, the ε parameter and SWS are not so sensitive to porosity changes. However, when porosity increased to 26.1%, P-wave anisotropy dramatically declined from 4.5% to 3.5%, whereas SWS remained fairly constant. Afterwards, as porosity increased to 30.8%, P-wave anisotropy continued to decline from 3.5%to 2.2%, but SWS fluctuated very slightly.

5.2. Modelling Insights and Discussions

The inputs for equivalent medium theories are given in Table 1 and Table 3. The Eshelby–Cheng theory (high frequency) predicted P-wave anisotropy to be close to the measured results [3,4]. Although the expressions of the Eshelby–Cheng theory do not contain background porosity, its predicted P-wave anisotropy showed a decreasing trend with increasing background porosity. The Eshelby–Cheng theory is expressed by Equations (1) and (2). Note that in Equation (1), φ f refers to the fracture-occupied porosity, not the background porosity, and in all the fractured samples, the fracture porosities are practically identical. The predicted P-wave anisotropy fluctuations were mainly dominated by parameter A i j in Equation (2), which relates the applied strain to the “stress-free” strain of the inclusion. A i j is expressed as follows:
A i j = B i j 1 L k j
where:
B i j = G S j + 2 H S i j + λ b + 2 μ b δ i j
L i j = G 2 H δ i j
H = μ μ b
G = λ λ b
S j = i = 1 3 S i j
where λ b and μ b are the lame constants of the background medium, and λ and μ are the lame constants of the inclusion material (water). S i j is the Eshelby S tensor, which relates the strain field (due to an ellipsoidal inclusion) to the “stress-free” strain in the inclusion. For very thin fractures (with a small aspect ratio), S i j can be simplified as follows:
S 11 = S 22 = π α ( 9 5 R ) / 16
S 33 = 1 π α R / 2
S 12 = S 21 = π α ( 3 7 R ) / 16
S 13 = S 23 = π α R / 4
S 31 = S 32 = 1 2 R π α ( 3 5 R ) / 4
S 55 = S 44 = 2 π α ( 3 2 R ) / 2
S 66 = π α ( 3 + R )
where α is the fracture aspect ratio, and
R = μ b λ b + 2 μ b
In the normal fracture direction (direction three), since α is so small ( α = 0.018 ), S 33 1 and S 3 1 . Because the shear modulus of the inclusion material (water) μ is zero, B 33 λ and L 33 λ b + 2 μ b λ . With increasing porosity, L 33 decreases, whereas B 33 stays nearly constant. This leads to a decrease in A i j . As a result, the fracture-induced perturbation in the normal fracture direction decreases with increasing porosity. However, as the aspect ratios of the fractures were extremely small in this study, the fracture-induced perturbations were nearly negligible in the parallel direction of the fractured samples. As the fracture-induced perturbations decreased with increasing porosity in the perpendicular direction, the stiffness difference between the parallel and perpendicular directions also decreased. This led to the declining trend observed for the predicted P-wave anisotropy with increasing porosity.
In Figure 7b, we can see that the Eshelby–Cheng theory predicted that SWS would be higher than the measured results. The shear moduli of S1 ( C 44 ) and S2 ( C 66 ) in the Eshelby–Cheng theory are as follows:
C 44 = μ [ 1 φ f ( μ μ ) μ + 2 ( μ μ ) S 44 ]
C 66 = μ [ 1 φ f ( μ μ ) μ + 2 ( μ μ ) S 66 ]
As μ = 0 in our study, SWS can be approximately expressed as follows:
S W S = C 66 C 44 2 C 44 4 e 3 2 R ε = C 11 s a t C 33 s a t 2 C 33 s a t = 1 2 ( C 11 d r y + α 1 2 M C 33 d r y + α 3 2 M 1 ) α 3 γ sat 4 3 e S W S
where e is the fracture density and φ f = 2 π α e . As the value of R in the four rocks was close (between 0.25 and 0.3), the SWS predictions were similar and not sensitive to changes in porosity.
In the theoretical predictions of [1], with increasing porosity, P-wave anisotropy showed a decreasing trend, whereas SWS remained practically constant. This trend was similar to the measured results. However, the theoretical predictions were higher than the measured results. For water-saturated conditions, the theoretically predicted P-wave anisotropy of [1] was as follows:
ε = C 11 s a t C 33 s a t 2 C 33 s a t = 1 2 ( C 11 d r y + α 1 2 M C 33 d r y + α 3 2 M 1 )
where direction three is the normal fracture direction. C 11 s a t and C 33 s a t are the P-wave moduli in the parallel and perpendicular directions, respectively, in water-saturated conditions. C 11 d r y and C 33 d r y are the P-wave moduli in the parallel and perpendicular directions, respectively, in dry conditions, and
C 33 d r y = 3 μ ( λ + μ ) 3 μ ( λ + μ ) + 4 ( 4 μ 2 + 4 λ μ ) e C 11 d r y
where λ and μ are Lame parameters of the dry background matrix. Equation (29) can be rewritten as follows:
ε = 1 2 ( C 11 d r y + α 1 2 M 3 μ ( λ + μ ) 3 μ ( λ + μ ) + 4 ( 4 μ 2 + 4 λ μ ) e C 11 d r y + α 3 2 M 1 )
Firstly, with increasing porosity, the value of 3 μ ( λ + μ ) 3 μ ( λ + μ ) + 4 ( 4 μ 2 + 4 λ μ ) e increases from 0.62 to 0.75. According to Equation (6), when the porosity increases, scalar M decreases greatly. Since α 1 is larger than α 3 (Equation (5)), α 1 2 M decreases more than α 3 2 M . As a result, with increasing porosity, these two factors cause the C 11 s a t to decrease more than C 33 s a t , resulting in a decrease in P-wave anisotropy. The predicted P-wave anisotropy was higher than the measured results. Potential explanations for this phenomenon include: (i) Under an experimental frequency of 0.5 MHz, the fluid did not have enough time to flow between the pores and fractures, thus the system did not reach pressure equilibrium in the pore space. This equilibrium is assumed by the Gurevich theory. This can cause a “stiffening” effect on the fractures, reducing P-wave anisotropy and causing the measured results to be lower than the theoretical predictions. (ii) Interactions between the fractures were not included in the Gurevich theory, which has been shown to reduce the overall anisotropy in a fractured medium in theoretical studies [5,6].
In the Gurevich theory, the expression of SWS is as follows:
S W S = C 55 s a t C 44 s a t 2 C 44 s a t
where:
C 55 s a t = μ
C 44 s a t = μ ( 1 Δ T )
Δ T = μ Z T 1 + μ Z T
Z T = 16 ( λ + 2 μ ) e 3 μ ( 3 λ + 4 μ )
In combining Equations (32) and (36), the following can be obtained:
S W S = 4 e 3 ( 1 λ 3 λ + 4 μ )
where e is the fracture density. The term λ 3 λ + 4 μ was extremely small in the four fractured samples (less than 0.077); therefore, γ sat 4 3 e . Compared with P-wave anisotropy, the S W S prediction was relatively close to the measured results.
Since the fracture aspect ratio is 0.02 in the fractured samples, we assumed that no partial or complete closure of the fractures occurred during the passage of ultrasonic waves. The above analysis is based on this assumption. For cases with smaller aspect ratios, we believe the partial or complete closure of fractures should be considered.

6. Conclusions

In this study, we investigated the effects of background porosity on seismic anisotropy in fractured rocks by conducting laboratory experiments. P- and S-wave velocities were measured on a set of synthetic samples with the same fracture density but different porosities (ranging from 15.1% to 30.8%) under water-saturated conditions. The results exhibited a significant effect of background porosity on P-wave anisotropy for the measured porosity range. The P-wave anisotropy fluctuated slightly when porosity increased from 15.1% to 22.1%, and then declined sharply when porosity continued to increase to 30.8%. However, SWS stayed fairly constant for the entire porosity range. The results were compared with predictions from equivalent medium theories. In both the Eshelby–Cheng model and the Gurevich model, increasing porosity reduces P-wave anisotropy, but has little influence on SWS. The experimental results show good qualitative agreement with the theoretical predictions.
Understanding the influence of background porosity on seismic anisotropy could improve the characterization of different fractured reservoirs. The findings of this study have potential applications for fracture detection in different areas with different reservoir physical properties. However, considering the inevitable differences between synthetic samples and natural rocks, further studies concerning this topic using field observations are required, such as post-survey and post-drilling calibrations of the seismic anisotropy after VSP.

Author Contributions

Sample manufacturing and ultrasonic measurement, Y.Z.; data analysis, Y.Z., F.G. and L.L.; supervision, B.D.; writing, Y.Z. and L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that supports the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gurevich, B. Elastic properties of saturated porous rocks with aligned fractures. J. Appl. Geophys. 2003, 54, 203–218. [Google Scholar] [CrossRef]
  2. Sil, S.; Sen, M.K.; Gurevich, B. Analysis of fluid substitution in a porous and fractured medium. Geophysics 2011, 76, WA157–WA166. [Google Scholar] [CrossRef]
  3. Eshelby, J.D. The determination of the elastic field of an ellipsoidal inclusion, and related problems. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 1957, 241, 376–396. [Google Scholar]
  4. Schoenberg, M.; Sayers, C.M. Seismic anisotropy of fractured rock. Geophysics 1995, 60, 204–211. [Google Scholar] [CrossRef] [Green Version]
  5. Cheng, C.H. Seismic Velocities in Porous Rocks: Direct and In-Verse Problems. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 1978. [Google Scholar]
  6. Cheng, C.H. Crack models for a transversely isotropic medium. J. Geophys. Res. Solid Earth 1993, 98, 675–684. [Google Scholar] [CrossRef]
  7. Rathore, J.S.; Fjaer, E.; Holt, R.M.; Renlie, L. P- and S-wave anisotropy of a synthetic sandstone with controlled crack geometry. Geophys. Prospect. 1995, 43, 711–728. [Google Scholar] [CrossRef]
  8. Ass’ad, J.M.; Tatham, R.H.; McDonald, J.A. A physical model study of microcrack-induced anisotropy. Geophysics 1992, 57, 1562–1570. [Google Scholar] [CrossRef]
  9. Wei, J.; Di, B. A physical model study of effect of fracture aperture on seismic wave. Sci. China Ser. D Earth Sci. 2008, 38, 211–218. [Google Scholar] [CrossRef]
  10. Tillotson, P.; Sothcott, J.; Best, A.I.; Chapman, M.; Li, X.-Y. Experimental verification of the fracture density and shear-wave splitting relationship using synthetic silica cemented sandstones with a controlled fracture geometry. Geophys. Prospect. 2011, 60, 516–525. [Google Scholar] [CrossRef]
  11. Tillotson, P.; Chapman, M.; Sothcott, J.; Best, A.I.; Li, X.-Y. Pore fluid viscosity effects on P- and S-wave anisotropy in synthetic silica-cemented sandstone with aligned fractures. Geophys. Prospect. 2014, 62, 1238–1252. [Google Scholar] [CrossRef]
  12. Amalokwu, K.; Best, A.I.; Sothcott, J.; Chapman, M.; Minshull, T.; Li, X.-Y. Water saturation effects on elastic wave attenuation in porous rocks with aligned fractures. Geophys. J. Int. 2014, 197, 943–947. [Google Scholar] [CrossRef] [Green Version]
  13. Amalokwu, K.; Chapman, M.; Best, A.I.; Minshull, T.A.; Li, X.-Y. Water saturation effects on P-wave anisotropy in synthetic sandstone with aligned fractures. Geophys. J. Int. 2015, 202, 1088–1095. [Google Scholar] [CrossRef] [Green Version]
  14. Amalokwu, K.; Chapman, M.; Best, A.I.; Sothcott, J.; Minshull, T.A.; Li, X.-Y. Experimental observation of water saturation effects on shear wave splitting in synthetic rock with fractures aligned at oblique angles. Geophys. J. Int. 2014, 200, 17–24. [Google Scholar] [CrossRef] [Green Version]
  15. Amalokwu, K.; Best, A.I.; Chapman, M. Effects of aligned fractures on the response of velocity and attenuation ratios to water saturation variation: A laboratory study using synthetic sandstones. Geophys. Prospect. 2016, 64, 942–957. [Google Scholar] [CrossRef] [Green Version]
  16. Amalokwu, K.; Papageorgiou, G.; Chapman, M.K.; Best, A.I. Modelling ultrasonic laboratory measurements of the saturation dependence of elastic modulus: New insights and implications for wave propagation mechanisms. Int. J. Greenh. Gas Control 2007, 59, 148–159. [Google Scholar] [CrossRef] [Green Version]
  17. Ding, P.; Di, B.; Wang, D.; Wei, J.; Li, X. P and S wave anisotropy in fractured media: Experimental research using synthetic samples. J. Appl. Geophys. 2014, 109, 1–6. [Google Scholar] [CrossRef]
  18. Ding, P.; Di, B.; Wei, J.; Li, X.; Deng, Y. Fluid-dependent anisotropy and experimental measurements in synthetic porous rocks with controlled fracture parameters. J. Geophys. Eng. 2014, 11, 015002. [Google Scholar] [CrossRef]
  19. Ding, P.; Di, B.; Wang, D.; Wei, J.; Li, X. Measurements of Seismic Anisotropy in Synthetic Rocks with Controlled Crack Geometry and Different Crack Densities. Pure Appl. Geophys. 2017, 174, 1907–1922. [Google Scholar] [CrossRef]
  20. Ding, P.; Wang, D.; Di, G.; Li, X. Investigation of the Effects of Fracture Orientation and Saturation on the Vp/Vs Ratio and Their Implications. Rock Mech. Rock Eng. 2019, 52, 3293–3304. [Google Scholar] [CrossRef]
  21. Ding, P.; Wang, D.; Li, X.-Y. An Experimental Study on Scale-Dependent Velocity and Anisotropy in Fractured Media Based on Artificial Rocks with Controlled Fracture Geometries. Rock Mech. Rock Eng. 2020, 53, 3149–3159. [Google Scholar] [CrossRef]
  22. Ding, P.; Wei, J.; Di, B.; Li, X.; Zeng, L. Observation and theoretical calibration of the fluid flow mechanism of artificial porous rocks with various size fractures. Geophys. Prospect. 2021, 69, 1235–1247. [Google Scholar] [CrossRef]
  23. Huang, L.; Stewart, R.R.; Dyaur, N.; Sil, S. Fluid substitution effects on seismic anisotropy. J. Geophys. Res. Solid Earth 2015, 120, 850–863. [Google Scholar] [CrossRef]
  24. Gassmann, F. Elastic Waves through a packing of spheres. Geophysics 1951, 16, 673–685. [Google Scholar] [CrossRef]
  25. Dellinger, J.; Vernik, L. Do travel times in pulse-transmission experiments yield anisotropic group or phase velocities. Geophysics 1994, 59, 1774–1779. [Google Scholar] [CrossRef]
Figure 1. (a) A tool used to strickle mineral mixture in the mould; (b) a structural representation.
Figure 1. (a) A tool used to strickle mineral mixture in the mould; (b) a structural representation.
Applsci 13 08379 g001
Figure 2. (a) Pouring the mixture into the mould; (b) setting the tool on the mould, vertically adjusting the glass plate by rotating the threaded rod and horizontally removing the tool, so that the glass plate could ‘flatten’ the mixture; (c) the mixture was paved flatly in the mould.
Figure 2. (a) Pouring the mixture into the mould; (b) setting the tool on the mould, vertically adjusting the glass plate by rotating the threaded rod and horizontally removing the tool, so that the glass plate could ‘flatten’ the mixture; (c) the mixture was paved flatly in the mould.
Applsci 13 08379 g002
Figure 3. (a) Four groups of samples with different porosities, (b) a section of the fractured sample, (c) an SEM image of the pore structure in the section, and (d) an SEM image of the fracture structure in the section.
Figure 3. (a) Four groups of samples with different porosities, (b) a section of the fractured sample, (c) an SEM image of the pore structure in the section, and (d) an SEM image of the fracture structure in the section.
Applsci 13 08379 g003aApplsci 13 08379 g003bApplsci 13 08379 g003c
Figure 4. P and S waveforms and the amplitude spectra of blank samples. (a) P waveform, (b) P-wave amplitude spectrum, (c) S waveform, and (d) S-wave amplitude spectrum.
Figure 4. P and S waveforms and the amplitude spectra of blank samples. (a) P waveform, (b) P-wave amplitude spectrum, (c) S waveform, and (d) S-wave amplitude spectrum.
Applsci 13 08379 g004aApplsci 13 08379 g004b
Figure 5. P-wave (a) and shear wave (b) velocities in water-saturated fractured samples (the P-wave velocity measurement error is smaller than the size of the labels).
Figure 5. P-wave (a) and shear wave (b) velocities in water-saturated fractured samples (the P-wave velocity measurement error is smaller than the size of the labels).
Applsci 13 08379 g005
Figure 6. P-wave anisotropy and shear wave splitting (SWS) in water-saturated blank samples.
Figure 6. P-wave anisotropy and shear wave splitting (SWS) in water-saturated blank samples.
Applsci 13 08379 g006
Figure 7. Fracture-induced P-wave anisotropy (a) and SWS (b) variation with porosity.
Figure 7. Fracture-induced P-wave anisotropy (a) and SWS (b) variation with porosity.
Applsci 13 08379 g007aApplsci 13 08379 g007b
Table 1. Physical parameters of the three sample groups (water-saturated).
Table 1. Physical parameters of the three sample groups (water-saturated).
Group 1Group 2Group 3Group 4
1-B1-F2-B2-F3-B3-F4-B4-F
Fracture density00.0600.0600.0600.06
Fracture aspect ratio/0.018/0.018/0.018/0.018
Porosity (%)14.115.320.622.125.326.130.430.8
Fracture porosity (%)/0.75/0.75/0.75/0.75
Bulk density
(g/cm3)
2.222.212.162.152.072.052.032.03
Compression pressure (MPa)186186535321211515
Table 2. P- and S-wave velocities in the blank samples (water-saturated).
Table 2. P- and S-wave velocities in the blank samples (water-saturated).
PorosityP-Wave Velocity (m/s)S-Wave Velocity (m/s)
(%)ParallelPerpendicularS1S2
1-B14.13930389621702120
2-B20.63442342018571820
3-B25.33128300515831495
4-B30.42817281113971390
Table 3. Main inputs for equivalent medium theories.
Table 3. Main inputs for equivalent medium theories.
Sample
Number
Background Matrix Modulus (GPa)Grain Modulus (GPa)Fluid Modulus (GPa)
K d r y μ d r y K s a t μ s a t K s μ s
1-F7.057.4120.3610.4530.624.52.15
2-F5.525.8215.667.4530.624.52.15
3-F4.694.1413.335.1930.624.52.15
4-F3.123.1110.833.9630.624.52.15
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Di, B.; Gao, F.; Li, L. Effects of Background Porosity on Seismic Anisotropy in Fractured Rocks: An Experimental Study. Appl. Sci. 2023, 13, 8379. https://doi.org/10.3390/app13148379

AMA Style

Zhang Y, Di B, Gao F, Li L. Effects of Background Porosity on Seismic Anisotropy in Fractured Rocks: An Experimental Study. Applied Sciences. 2023; 13(14):8379. https://doi.org/10.3390/app13148379

Chicago/Turabian Style

Zhang, Yuangui, Bangrang Di, Feng Gao, and Lei Li. 2023. "Effects of Background Porosity on Seismic Anisotropy in Fractured Rocks: An Experimental Study" Applied Sciences 13, no. 14: 8379. https://doi.org/10.3390/app13148379

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop