Next Article in Journal
The Effects of Isometric Fatigue on Trunk Muscle Stiffness: Implications for Shear-Wave Elastography Measurements
Next Article in Special Issue
Improved Speech Spatial Covariance Matrix Estimation for Online Multi-Microphone Speech Enhancement
Previous Article in Journal
Dynamic Smoothing, Filtering and Differentiation of Signals Defining the Path of the UAV
Previous Article in Special Issue
Performance Analysis of Active Structural Acoustic Control Applied to a Washing Machine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Three-Dimensional Acoustic Device for Testing the All-Directional Anisotropic Characteristics of Rock Samples

1
School of Geosciences and Technology, China University of Petroleum, Qingdao 266555, China
2
Yichang Chengfa Gezhouba Gas Development Corporation, Yichang 443000, China
3
ChangQing Branch of CNPC Logging, Xi’an 710200, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(23), 9473; https://doi.org/10.3390/s22239473
Submission received: 7 November 2022 / Revised: 27 November 2022 / Accepted: 1 December 2022 / Published: 4 December 2022
(This article belongs to the Special Issue Advanced Technology in Acoustic Signal Processing)

Abstract

:
Many oil and gas fields, especially non-conventional shale and compacted sand reservoirs, have formation anisotropy. The acoustic anisotropy measurement of cores in these reservoirs can guide drilling, well logging, and exploitation. However, almost all core holders are designed for cylinder cores, which are not suitable for all-directional measurements. A three-dimensional measurement device was designed on the basis of the cross-hole sonic logging method. This device mainly consisted of two pairs of transducers, a signal generator, an oscillograph, an omnidirectional positioning system, and a computer control system. By adjusting the measurement latitude and longitude circle automatically, this device scanned spherical sample rocks and obtained full-wave waveforms in all directions. Experiments were performed taking granite from the Jiaodong Peninsula, China, as an example, and the arrival times and velocities of the longitudinal and shear waves were calculated based on the full-wave waveforms. Thereafter, anisotropic physical characterizations were carried out on the basis of these velocities. These data play an important role in guiding formation fracturing and analyzing the stability of borehole walls.

1. Introduction

In oil exploration and exploitation, core sampling is often performed on object formations for physical parameter analyses, which can ascertain geological components, rock characteristics, and the processes of geological changes [1,2,3]. Core holders are analysis tools that fix target cores and enable certain kinds of measurements, such as resistivity, permeability, computed tomography, and acoustic velocity measurements. Amongst these measurements, the acoustic parameters of cores can reflect the velocities and amplitudes of longitudinal and shear waves [4]. Then, the lithology, porosity, elasticity, and mechanical characteristics of the cores can be obtained. Since the 1960s, scholars have developed many kinds of devices to obtain the petrophysical parameters of sampling cores through the measurement of acoustic characteristics [5,6]. They discussed the relationships between the acoustic propagation abilities and the physical properties of rocks [7]. Such properties include mechanical parameters, rock structure, mineral composition, and microfractures. These relationships are very important in the exploration and development of oil and gas fields.
Many oil and gas fields, especially non-conventional shale and compacted sand reservoirs, have formation anisotropy [8,9]. The anisotropy measurement of cores in these reservoirs is a challenging issue that can guide drilling, well logging, and exploitation [10]. Given that acoustic signals are highly sensitive to formation anisotropy, acoustic measurement is an effective way to test anisotropy. Almost all core holders are designed for cylinder cores, which are considerably easier to produce and test than other core types. However, cylinder cores are unsuitable for anisotropic measurements [11]. Traditional measurement methods sample and test the same core in different directions to investigate omnidirectional characteristics (Figure 1). However, they cannot obtain the whole three-dimensional (3-D) data, and the cores must be destroyed. In consideration of this issue, a 3-D measurement device is designed to measure the acoustic parameters of anisotropic cores on the basis of the acoustic wave transmission method. By scanning a spherical sample rock in all directions, the full-wave waveforms in all directions are obtained, and the arrival times and velocities of the longitudinal and shear waves are calculated. Thereafter, anisotropic physical characterization can be carried out on the basis of these velocities.

2. Design and Test of the 3-D Measurement Device

2.1. Structural Design of the Device

The pulse wave reflection and cross-hole sonic logging methods are the two main methods for testing rock elastic velocities. However, if object rocks have the characteristic of strong acoustical attenuation, the pulse wave reflection method obtains low signal–noise ratios. In the cross-hole sonic logging method, transmitters and receivers are fixed on the opposite sides of a sample rock, and the propagation times and lengths of the longitude and shear waves are measured to calculate their velocities. The 3-D measurement device was based on the cross-hole sonic logging method. It could realize high-precision positioning and omnidirectional scanning.
Figure 2 shows that the device mainly consisted of two pairs of transducers [12], a signal generator, an oscillograph, an omnidirectional positioning system [13,14], and a computer control system. The transducers were two pairs of high-precision transmitting and receiving transducers for longitudinal and shear waves, and their domain frequencies were 1 MHz. They were used to carry out transformations between excitation electrical signals and acoustic signals [15]. The signal generator produced pulse electrical signals, which were then amplified by a power preamplifier and loaded on the transmitting transducers. The oscillograph collected weak electronic signals from the receiving transducers [16]. The positioning system mainly consisted of two self-designed transducer holders and a rock sample holder. The transducer holders fixed the transducers, carried out their rotational and positioning motions, and controlled the contact pressure between the transducers and the rock sample with an air pump plunger. The sample holder supported the rock sample with a support rod and performed its rotational motion horizontally. The computer control system exerted overall control and recorded and calculated the collected data.

2.2. Workflow of the Measurement Device

A schematic of the measurement conducted by the 3-D measurement device is shown in Figure 3, and the overall workflow was as follows.
  • Measure the radius of the rock sample and enter this data into the computer control software.
  • Push the rock sample into the instrument and lock it vertically and horizontally. Thereafter, automatically lift the rock sample with the support rod to a suitable position.
  • Place the pair of transducers into the holders and fix them to the initial measurement positions.
  • Produce electrical excitation signals with the signal generator that are amplified by the amplifier and loaded on the transmitting transducers to generate acoustic signals.
  • Propagate the acoustic signals into the rock. These reach the receiving transducers and are then collected by the oscilloscope and uploaded to the computer.
  • Rotate the transducers synchronously in the vertical direction by a latitude of ϕ and repeat steps 3–5 until all the measurements in one circle are finished.
  • Loosen the rock sample, rotate it by a longitude of λ in the horizontal direction, and re-lock it again. Then, repeat steps 3–6 until the measurements are completed in all directions.
The device had the function of automatic measurement. By inputting the latitude step ∆ϕ and the longitude step ∆λ into the computer control system, the device performed measurements automatically to reduce manpower consumption and improve efficiency.

3. Examples of Measurements by the Device

The measurement effect was investigated by taking granite from the Jiaodong Peninsula, China, as an example (Figure 4). The produced spherical rock sample had a black-and-white interlaced laminar structure. The white and black components were quartz and biotite, respectively. The sample had a diameter of 0.30 m and a density of 2850 kg/m3. During measurement, the transducers rotated simultaneously to change the latitude angle, and the support rod rotated the rock sample to change the longitude angle.

3.1. Measurement of Longitudinal Waves

Firstly, the rock sample was measured by the transducers for longitudinal waves. Vaseline was smeared over the rock sample and the transducers to couple acoustic signals. The working frequency of the transducers was 1 MHz, and the longitude and latitude steps were set to 15° and 10°, respectively. The full-wave received waveforms in different dimensions measured at the longitude of 150° are shown in Figure 5. The velocity and amplitude of the longitude wave were the strongest when the latitude angle was 90° (parallel to the bedding plane). When the latitude angle was changed to the vertical direction of the bedding plane, the velocities of the waves decreased, and their arrival times followed a sinusoidal trend.
The arrival times of the longitudinal waves were calculated on the basis of the scanning-measurement endings. As the acoustic signals were single and clear, the threshold method was utilized to obtain the arrival times of the longitudinal waves. Figure 6 shows the 3-D imaging results, wherein differences in colors indicate differences in arrival times. The arrival times of the acoustic longitudinal waves could be concluded to conform to the same law for different longitude circles. When the connecting directions of the two transducers were parallel to the bedding direction of the rock sample, the arrival time was the shortest, and the propagation speed was the fastest. When the connection direction was perpendicular to the bedding direction, the arrival time was the longest, and the propagation speed was the slowest.

3.2. Measurement of Shear Waves

The velocities of shear waves were measured with a pair of shear wave transducers. Given that shear waves are excited by shearing motions, the receiving transducer must be in the polarization direction of the firing transducer. The conditions and steps of the measurement of shear waves were the same as those for the measurement of longitudinal waves, except that the excitation frequency was set to 0.5 MHz. The 3-D image of the shear wave arrival times is shown in Figure 7. The arrival times of shear waves were significantly later than those of the longitudinal waves shown in Figure 6 because the propagation velocities of shear waves were considerably lower than those of longitudinal waves. Meanwhile, the anisotropy differences followed the same trend: higher in the latitudinal directions and lower in the longitudinal directions.

4. Discussion

4.1. Calculation of Velocities

Equation (1) could be used to analyze the velocity anisotropy of longitudinal waves in different directions.
k % = v m a x v m i n v m e a n × 10
where k is the degree of anisotropy, v m a x is the maximum velocity, v m i n is the minimum velocity, and v m e a n is the average velocity. In accordance with the arrival times, the anisotropies in the longitudinal direction (vertical direction) and latitudinal direction (parallel bedding direction) were calculated and are shown in Figure 8. The anisotropy differences were approximately 3% in the latitudinal direction and 12% in the longitudinal direction. That is, anisotropy differences existed in the vertical direction, which was in line with the bedding structure of the rock sample.
The velocities of the longitudinal and shear waves of the same testing point were combined and fitted. Their relationships are shown in Figure 9. Their velocities were positively related. Their relationship equation is shown in Equation (2).
V s = 0.3453 V p + 1562.7
The calculated goodness of fit was 0.7416. There were some outside points that were mostly caused by measurement errors, such as the absence of acoustic coupling in these points. Comparison with Ref. [6] shows that the relationship between the compressional and shear wave velocities of the granite samples is similar to those of sandstone and mudstone.

4.2. Calculation of Elastic Physical Parameters

The measurement results could be used to judge the orientation of fractures and to evaluate the brittleness of rocks by inverting their physical parameters because the propagation velocities of acoustic waves in rock are closely related to many physical parameters, such as density, Young’s modulus, Poisson’s ratio, and shear modulus [17]. Therefore, the relationships between the above physical parameters and the velocities of the longitudinal and shear waves were derived.
Young’s modulus:
E = 2 G ( 1 + μ ) = 2 ρ ( 1 + μ ) V s 2 = 9 K ρ V s 2 3 K + ρ V s 2 ,
where ρ is the density of the rock, K is the bulk modulus, and V s is the velocity of the shear waves.
Shear modulus:
G = ρ V s 2 .
Bulk modulus:
K = ρ V p 2 4 3 V s 2 .
Poisson’s ratio:
μ = V p V s 2 2 V p V s 2 1 .
The Poisson’s ratio, bulk modulus, shear modulus, and Young’s modulus values measured in different latitude circles were calculated by substituting the measured velocities of the longitudinal and shear waves into the above formulas and are shown in Figure 10. The above elastic parameters took the maximum values when the latitude angles were approximately from 90° to 100° (parallel to the bedding direction) and the minimum values when the latitude angles were approximately 0° and 160°. This phenomenon also reflected the anisotropic differences in different physical parameters. The data on the sampling core play an important role in guiding formation fracturing and analyzing the stability of borehole walls in the processes of exploration and exploitation.

5. Conclusions

A 3-D acoustic measurement device was designed on the basis of the acoustic wave transmission method to measure the anisotropic physical characteristics of formation cores. Its structural design and workflow were introduced. Its measurement effect was investigated by taking granite from the Jiaodong Peninsula, China, as an example. Longitudinal and shear waves were omnidirectionally excited and collected, and their velocities and acoustic anisotropy characteristics were calculated.
  • The measured granite rock samples exhibited anisotropy characteristics. The velocities of its longitude and shear waves changed with the measured latitude circle, and the velocities were biggest when the circles were parallel to the bedding direction of the rock sample. The anisotropy characteristics were very small in the longitude orientations, and little velocity change happened.
  • The velocities of the longitudinal and shear waves were positively related. The calculated goodness of fit was 0.7416. The relationship between the compressional and shear wave velocities of the granite samples was similar to those of sandstone and mudstone.
  • The anisotropy characteristics of Poisson’s ratio, bulk modulus, shear modulus, and Young’s modulus were calculated based on the measured velocities of the longitudinal and shear waves. These data play an important role in guiding formation fracturing and analyzing the stability of borehole walls in the processes of exploration and exploitation.
  • The design principle of the 3-D acoustic measurement device could also be used in devices for other kinds of measurements, such as resistivity, permeability, and computed tomography measurements.

Author Contributions

Conceptualization, K.Z. and W.W.; validation and formal analysis, S.L.; writing—original draft preparation, K.Z. and S.L.; mechanic design, S.X. and Y.S.; device testing, K.Z. and B.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Major Scientific and Technological Projects of CNPC (Grant No. ZD2019-183-004), the Natural Science Foundation of China (Grant No. 41821002 and 41774138), and the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA14020302 and XDA14020304).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

All authors declare no conflict of interest.

References

  1. Kawamura, T.; Sakamoto, Y.; Ohtake, M.; Yamamoto, Y.; Haneda, H.; Komai, T. Dissociation behavior of hydrate core sample using thermodynamic inhibitor–Part 2: Experimental investigation using long core samples. Int. J. Offshore Polar Eng. 2008, 18, 156–159. [Google Scholar]
  2. Behney, A.C.; O’Shaughnessy, R.; Eichholz, M.W.; Stafford, J.D. Influence of Item Distribution Pattern and Abundance on Efficiency of Benthic Core Sampling. Wetlands 2014, 34, 1109–1121. [Google Scholar] [CrossRef]
  3. Inada, N.; Yamamoto, K. Data report: Hybrid Pressure Coring System tool review and summary of recovery result from gas-hydrate related coring in the Nankai Project. Mar. Pet. Geol. 2015, 66, 323–345. [Google Scholar] [CrossRef]
  4. Jones, L.E.A.; Wang, H.F. Ultrasonic velocities in Cretaceous shales from the Williston Basin. Geophysics 1981, 46, 288–297. [Google Scholar] [CrossRef]
  5. Simmons, G. Velocity of shear waves in rocks to 10 kilobars, 1. J. Geophys. Res. 1964, 69, 1123–1130. [Google Scholar] [CrossRef]
  6. Birch, F. The velocity of compressional waves in rocks to 10 kilobars: 1. J. Geophys. Res. 1960, 65, 1083–1102. [Google Scholar] [CrossRef]
  7. Christensen, N.I. Compressional wave velocities in metamorphic rocks at pressures to 10 kilobars. J. Geophys. Res. 1965, 70, 6147–6164. [Google Scholar] [CrossRef]
  8. Thomsen, L. Weak elastic anisotropy. Geophysics 1986, 51, 1954–1966. [Google Scholar] [CrossRef]
  9. Johnston, J.E.; Christensen, N.I. Seismic anisotropy of shales. J. Geophys. Res. Solid Earth 1995, 100, 5991–6003. [Google Scholar] [CrossRef] [Green Version]
  10. King, M.S. Recent developments in seismic rock physics. Int. J. Rock Mech. Min. Sci. 2009, 46, 1341–1348. [Google Scholar] [CrossRef]
  11. Shakerian, M.; Marica, F.; Afrough, A.; Goora, F.G.; Li, M.; Vashaee, S.; Balcom, B.J. A high-pressure metallic core holder for magnetic resonance based on Hastelloy-C. Rev. Sci. Instrum. 2017, 88, 123703. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, K.; Tan, B.; Liu, X. A Consistency Evaluation and Calibration Method for Piezoelectric Transmitters. Sensors 2017, 17, 985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zhang, K.; Ju, X.D.; Tan, B.H.; Lu, J.Q.; Men, B.Y.; Wu, W.H.; Chen, J.Y. New excitation method for acoustic logging transmitters. J. Geophys. Eng. 2017, 14, 751–757. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, K.; Li, F.; Wu, J.; Tan, B.; Liu, L. Pulse width research on half-sine excitation signal for bending vibrator. J. Geophys. Eng. 2022, 19, 67–74. [Google Scholar] [CrossRef]
  15. Zhu, K.; Ma, J.; Qi, X.; Shen, B.; Liu, Y.; Sun, E.; Zhang, R. Enhancement of Ultrasonic Transducer Bandwidth by Acoustic Impedance Gradient Matching Layer. Sensors 2022, 22, 8025. [Google Scholar] [CrossRef] [PubMed]
  16. Tan, B.; Zhang, K.; Su, Y.; Li, S.; Zhang, L. Research on acoustic logging while drilling transmitting technologies. J. Geophys. Eng. 2022, 19, 511–520. [Google Scholar] [CrossRef]
  17. Sinha, B.K. Rock Static Moduli from Borehole Sonic Data in Stress-Dependent Formations. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2022, 69, 2371–2380. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Traditional measurement methods sample the same core in different directions.
Figure 1. Traditional measurement methods sample the same core in different directions.
Sensors 22 09473 g001
Figure 2. Actual picture of the three-dimensional anisotropic measurement device: (a) omnidirectional positioning system, (b) signal generator and amplifier, and (c) oscillograph.
Figure 2. Actual picture of the three-dimensional anisotropic measurement device: (a) omnidirectional positioning system, (b) signal generator and amplifier, and (c) oscillograph.
Sensors 22 09473 g002
Figure 3. Schematic of measurements conducted by the 3-D measurement device.
Figure 3. Schematic of measurements conducted by the 3-D measurement device.
Sensors 22 09473 g003
Figure 4. Granite rock sample from Jiaodong Peninsula, China.
Figure 4. Granite rock sample from Jiaodong Peninsula, China.
Sensors 22 09473 g004
Figure 5. Full-wave received waveforms in different dimensions measured at the longitude of 150°.
Figure 5. Full-wave received waveforms in different dimensions measured at the longitude of 150°.
Sensors 22 09473 g005
Figure 6. Arrival times of longitudinal waves on the basis of scanning-measurement endings.
Figure 6. Arrival times of longitudinal waves on the basis of scanning-measurement endings.
Sensors 22 09473 g006
Figure 7. Arrival times of shear waves on the basis of scanning-measurement endings.
Figure 7. Arrival times of shear waves on the basis of scanning-measurement endings.
Sensors 22 09473 g007
Figure 8. Velocity anisotropies in the longitudinal direction and latitudinal direction.
Figure 8. Velocity anisotropies in the longitudinal direction and latitudinal direction.
Sensors 22 09473 g008
Figure 9. Measured velocities of the longitudinal and shear waves of the same testing point (blue dots) and their fit curve (red line).
Figure 9. Measured velocities of the longitudinal and shear waves of the same testing point (blue dots) and their fit curve (red line).
Sensors 22 09473 g009
Figure 10. Calculated elastic anisotropic parameters: (a) Poisson’s ratio, (b) bulk modulus, (c) shear modulus, and (d) Young’s modulus.
Figure 10. Calculated elastic anisotropic parameters: (a) Poisson’s ratio, (b) bulk modulus, (c) shear modulus, and (d) Young’s modulus.
Sensors 22 09473 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, K.; Li, S.; Su, Y.; Tan, B.; Wu, W.; Xin, S. Three-Dimensional Acoustic Device for Testing the All-Directional Anisotropic Characteristics of Rock Samples. Sensors 2022, 22, 9473. https://doi.org/10.3390/s22239473

AMA Style

Zhang K, Li S, Su Y, Tan B, Wu W, Xin S. Three-Dimensional Acoustic Device for Testing the All-Directional Anisotropic Characteristics of Rock Samples. Sensors. 2022; 22(23):9473. https://doi.org/10.3390/s22239473

Chicago/Turabian Style

Zhang, Kai, Shengqing Li, Yuanda Su, Baohai Tan, Wenjie Wu, and Shoutao Xin. 2022. "Three-Dimensional Acoustic Device for Testing the All-Directional Anisotropic Characteristics of Rock Samples" Sensors 22, no. 23: 9473. https://doi.org/10.3390/s22239473

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