2.1. BENDER Bar and Its Vibration Principles
In this section, we first introduce the transducer vibration principles, and also conducted numerical analyses of the actuator and the whole transducer.
A piezoelectric bender bar is shown in
Figure 1. A metal layer and two piezoelectric ceramic layers are bonded at the two interfaces. The
z-axis is parallel to both the thickness direction of the bender bar and the polarization direction of both piezoelectric ceramic layers. The
x-axis is parallel to the length direction of the bender bar and the bar is symmetric about the
x–
y-plane. The middle plane of the bender bar bisects its thickness. Its vibration principle is that when a voltage is applied to both surfaces of each piezoelectric element, one extends and the other one contracts, which produces the bending motion of the entire bender bar [
13,
14].
In our previous work [
15], we considered the theory of the small deflection of a thin plate. Because the surfaces of the bender bar are traction-free, it is assumed that the stresses in the volume vanish, except for stress
and
, i.e.,
. The displacement in the
z-direction due to bending of the bar is denoted as
. We divided the bender bar into two parts; one is composite three-layer bender with two piezoelectric ceramic layers and one metal layer, and the other one is the isotropic metal layer part. In the middle composite layer part (
), the piezoelectric equations are as follows:
where
are the stress components;
are the strain components;
is the Poisson’s ratio of piezoelectric material when electric displacement is constant;
is the compliance coefficient component;
is the dielectric isolation rate component;
is the piezoelectric constant component; and
are the electric field intensity component and electric displacement component, respectively.
For the isotropic metal layer part (
), the stress and strain equations are shown as follows [
15]:
where
is the Young’s modulus of metal material, and
is the Poisson’s ratio of metal material.
By deriving the formula in our previous work, we obtained the vibration equation of the composite layer:
where
;
is the stiffness coefficient of composite layer;
, with
as the thickness of the piezoelectric material and the metal material, respectively;
is the average density of the composite layer, with
being the density of the piezoelectric material and the metal material, respectively.
For the metal layer at both sides (
), the vibration equation is:
where
is the stiffness coefficient of the metal layer.
Suppose that the usual excitation is symmetric, the symmetrical vibration modes of bender bar are excited. For Equation (3) of the composite layer, the displacement is denoted as:
where
are constants to be determined.
For Equation (4) of the isotropic metal part, the displacement is denoted as:
where
are constants to be determined.
Considering the continuity boundary condition at the junction between the composite layer and the metal layer, the physical function such as the transverse displacement, the shear force, the rotation angle, and the bending moment are continuous [
15]. Moreover, the fixed boundary conditions need to be met at both of the ends of the bender bar. After substituting transverse displacement expressions in Equations (5) and (6) into continuity and fixed boundary conditions, we obtained the following coefficient matrix equations:
According to the coefficient determinant value , we calculated the resonant frequencies and corresponding vibration modes.
For the long and short benders that we used to actuate our cylindrical oscillating transducer, the resonant frequencies were 402 Hz, 3692 Hz, 767 Hz and 4705 Hz, calculated by our theory. The first two frequencies corresponded to the first and third orders of the bending modes of the long bender, and the latter two frequencies corresponded to those of the short bender. We used the finite element method to simulate vibration modes of the bender bar, as shown in
Figure 2.
2.2. Transducer Structures and Its Vibration Principles
For the traditional piezoelectric dipole source in borehole acoustic measurement, due to the limited strain of the piezoelectric ceramic material itself, it is not easy to generate large displacement at low frequencies [
10,
11]. Therefore, it is hard to enhance the transmitting voltage response at low frequencies, especially at frequencies around less than 1 kHz.
In addition, the frequency band of traditional source is not wide enough that it can measure shear wave velocity in soft formations.
Therefore, in order to improve its radiation performance, especially in the low frequency range, and to extend the bandwidth of the dipole transducer, we proposed a new transducer design that we have patented in ref [
16], which is shown in
Figure 3.
The schematic diagram and the sectional view along the vertical axis of the piezoelectric cylindrical transducer with its frame is shown in
Figure 3. The rectangular coordinate system is set up in the acoustic center of the transducer. The
y-axis and
z-axis are shown in
Figure 3, and the
x-axis is in the inward direction perpendicular to the paper.
The dipole transducer consists of an internal piezoelectric actuator and external radiation structures. The internal actuator consists of three parallel piezoelectric benders, one of which is longer in the center, and the other two are shorter and of a same size. Each bender consists of two layers of piezoelectric elements and one layer of metal base element, which is commonly used in practice [
8]. As shown in
Figure 3, three layers are bonded at the two interfaces. The two ends of three parallel positioned benders are fixed to the endcaps bonded with the cylindrical shell. The external radiation structures consist of a cylindrical shell, two endcaps, and two U-shaped springs. The chamber, formed by the cylindrical shell and two endcaps, is filled with air. The cylindrical shell is supported by two U-shaped springs with the end faces fixed to the outside frame of a logging tool. The entire system is immersed in the fluid of the borehole. For numerical modeling, the whole structure is set to be immersed in air or water. During the operation, when a voltage is applied to the electrodes of the two piezoelectric elements, one piezoelectric element elongates and the other one contracts, causing the piezoelectric bender to bend. This makes the cylindrical shell oscillate out of the phase with the inner piezoelectric actuator.
The cylindrical shell, acting as the radiation surface of the transducer, is driven by the piezoelectric actuator to produce a pure dipole signal. Because the vibration of the piezoelectric actuator is confined in the space of cylindrical housing, it will not be harmful to the force imparted to the borehole wall [
17]. Suppose that a force
F is applied to a borehole wall through borehole fluid by a force
F1 from the inner piezoelectric actuator to the cylindrical radiation surface. According to the Newton’s third law, a force
F2, equivalent in magnitude but opposite in direction to the force
F1 is applied to the piezoelectric actuator, and an anti-reaction force harmful to the force
F would not be imparted to the side wall of the borehole because the inner actuator is fixed in the cylindrical chamber and separated from the borehole fluid. Thus, any waves due to the reaction are not radiated toward the borehole wall [
17]. On the other hand, for the conventional bender-type transducers, together with their frame, the opening angle of the sonic window is about 60 degrees due to the limited radiation surface. However, for the proposed transducer, with a full azimuthal angle of radiation surface, the opening angle of the sonic window can be as large as 120 degree or even more, which is efficient to inhibit the hexapole mode in a high frequency range [
10]. Compared with the traditional bender transducer, the radiation structure of the cylindrical shell in the transducer greatly increases the radiation area, adapting to the size of the cylindrical borehole as much as possible, and thus effectively enhances the acoustic radiation resistance and transmitting voltage response.
The transducer vibration modes and responses are simulated by using the FEM. Firstly, the physical model of the transducer is set up; then, the simulation and analyses of electrical conductance, transmitting voltage response, and directivity patterns are conducted, respectively, to illustrate the effectiveness of the transducer for generating acoustic waves in the frequency range of 0.4–6 kHz. Finally, the optimization results are discussed. The finite element simulation model of the piezoelectric cylindrical transducer associated with its frame is shown in
Figure 4.
2.6. Transmitting Voltage Response (TVR)
Based on our design, we fabricated and tested the transducer to compare with simulation results.
As is shown in ref [
18], the transmitting voltage response level is defined with the following formula:
where
is the pressure in the far field, usually far away from the source for about 3–5 wavelengths;
D is the distance from the transducer to the receiver; and
V is the excited voltage. The reference pressure in Equation (9) is taken to be 1 × 10
−6 ·m/V.
The comparison of the simulation and test results of TVR of the proposed transducers in water is shown in
Figure 8. The solid black line and the dashed blue line represent the simulation and the test results of the proposed transducer, respectively.
From the simulation of the transducer in the frequency range of 400–6000 Hz, it is shown that there are multiple resonance peaks in the curve. The four main resonant frequencies are 497 Hz, 980 Hz, 3698 Hz and 5000 Hz, respectively. The TVR values are 99 dB, 121 dB, 123 dB and 127 dB, respectively, corresponding to the four mode modes in the previous modal analysis.
In the test measurement curve, there are mainly four resonance peaks, of which the frequencies are 450 Hz, 925 Hz, 3550 Hz and 4700 Hz, respectively, and TVR values of 101 dB, 121 dB, 133 dB and 127 dB, respectively. For the resonance frequencies, the relative error between the simulation and test results is less than 11%. For the TVR value, the relative error between the simulation and test results is less than 8%. Therefore, the numerical simulation result and the test result are in good agreement.
2.7. Directivity Patterns
The numerical simulation and test measurement results of horizontal directivity of the proposed transducer at 979 Hz, 3690 Hz, and 4816 Hz are compared and shown in
Figure 9. The difference between them is due to the practical measurement error because the minimum angle measurement interval is not as small as the simulation result. The results show that the performance of the developed transducer is a pure dipole source in the operating frequency range.
In theory, the dipole sound field is rather pure due to the cylindrical shell acting as the radiation surface driven by the piezoelectric actuator. Firstly, it is because of an infinite cylinder which oscillates at the velocity of
in the
x-axis direction (
= 0), perpendicular to the vertical axis of the cylinder, the sound pressure (supposing the
z-axis coincides with the vertical axis of the cylinder, then sound pressure
p is unrelated to the
z-coordinate axis) can be expressed as [
19]:
where
represent the radial and angular coordinate, respectively;
are the density and sound velocity of the medium, respectively;
a is the radius of the circular cross section; and
,
k and
t are the angular frequency, wave number and time, respectively [
19]. The normalized form of the horizontal radiation pattern is:
In this work, although the cylindrical shell is finite in length, the horizontal directivity pattern can be considered similar to that of the ideal infinite cylinder because the wavelength is much more than the geometric size of the shell. Hence, the transducer with the cylindrical shell would perform as a pure dipole source. Moreover, with a full azimuthal angle of radiation surface, the vibration of the piezoelectric actuator is confined in the space of cylindrical housing so that the vibration otherwise coupling into the fluid of the borehole is effectively inhibited from interfering with the sound field by the outer radiation surface [
17].