Shock and Vibration
Volume 2017 (2017), Article ID 2670218, 12 pages
https://doi.org/10.1155/2017/2670218

Research Article

Dynamic Characteristics Analysis of a Seismic Vibrator-Ground Coupling System

Liu Jun, Huang Zhi-Qiang, Li Gang, Li Shi-De, Liu Xian, and Zhang Kai

School of Mechatronic Engineering, Southwest Petroleum University, No. 8 Xindu Avenue, Xindu District, Chengdu, Sichuan 610500, China

Correspondence should be addressed to Liu Jun; 99803392@qq.com

Received 13 March 2017; Accepted 22 May 2017; Published 4 July 2017

Academic Editor: Ivo Caliò

Copyright © 2017 Liu Jun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A mathematical model to describe the seismic vibrator and ground coupling is proposed based on dynamic analyses of the vibrator baseplate structure and properties of soil. This mathematical model is solved using Newmark’s method. It produces reasonable results and the trend of amplitude envelope follows the experimental data. With this model, characterization of the vibrator-ground coupling is performed. Modal analysis study shows that the motions between the vibrator reaction mass and baseplate experience a combination of a first-order vibration mode and a second-order vibration mode. The corresponding frequencies of two vibration modes are located at 2.77 Hz and 488.7 Hz, respectively. As the frequency goes up, the motion behavior of the reaction mass and baseplate is dominated by the second-order vibration mode. It is realized that the second-order vibration mode is responsible for the phase difference between the input sweep signal and the output force signal. Furthermore, the impact of the coupling system parameters on the vibrator output force is also investigated. It is observed that the vibrator output force decreases as the sweep frequency increases. The weight ratio of the reaction mass and the baseplate has an impact on the vibrator output force.

1. Introduction

Due to efficient operation and being environment friendly, seismic vibrators (Figure 1) have become a main energy source on land in oil and gas exploration. With extension of seismic acquisition frequency bandwidth towards low and high ends, seismic vibrators face more challenges to produce a predictable and invariant output force signal. To achieve this goal, it is very necessary to understand the impact of key parameters in the vibrator-ground system on the vibrator output force [1].

Figure 1: Seismic vibration system and vibrator.

Chapman et al. [2] presented a new design of the vibrator actuator to improve the vibrator output force in the frequency range of 200 Hz. They stressed that the weight and stiffness of the vibrator baseplate played a big role in extension of the vibrator output force bandwidth towards 200 Hz. Sallas and Weber [3] used a mass-spring-damper model to describe the vibrator system and the coupled ground system. Under rigid-body assumptions of reaction mass and baseplate, they theoretically demonstrated that a vibrator output force could be estimated using a weighted-sum method that is summing reaction mass and baseplate acceleration weighted by their respective masses. This weighted-sum ground force estimate provided a mean to approach the vibrator-ground force. However, this weighted-sum ground force estimate was proved to be valid in a narrow frequency bandwidth by many researchers [4–8].

Wei [9, 10] demonstrated that the frequency bandwidth of the vibrator output force was constrained by hydraulic and mechanical limits of the vibrator system. Based on experimental testing results, he derived formulas that can be used to estimate the maximum fundamental force from the vibrator at low frequencies. Subsequently, Wei [11] simulated the dynamic behavior of the vibrator baseplate as well as its coupled ground using finite element analysis tools. Then, he [11] concluded that the vibrator performance was significantly affected by the baseplate stiffness.

To improve the vibrator force output towards a wide frequency range, one way is to design broadband sweeps using nonlinear sweep design techniques [12, 13]. The other way is to improve the performance of fundamental components in the vibrator mechanical and hydraulic system [14, 15]. Although improvements of individual components in the vibrator system are seen [5–16], the theoretical progress on the seismic vibrator design moves very slowly. This is because the dynamics between the vibrator and ground coupling are unclear.

The purpose of this paper is to present a mathematical model that is used to simulate the vibrator-ground coupling. Next, an impact of vibrator system parameters on the vibrator-ground coupling is investigated. Last, an impact of these parameters on the quality of the vibrator force output is performed.

2. Ground Stiffness and Damping

2.1. Surface Displacement Field of a Half-Infinite Ground Space under Concentrated Vertical Harmonic Disturbances

The displacement field of a half-infinite ground space under a concentrated vertical harmonic disturbance force (Figure 2) can be written as follows [17, 18]:

with boundary conditions

where , , volume strain , and rotational component refer to Table 3 and, here, , , and are the density, elastic modulus, and Poisson’s ratio of the ground, respectively. and are displacement components and and are stress components.

Figure 2: Half-infinite space under a concentrated vertical harmonic disturbance force.

The surface displacement field of the half-infinite space can be derived as follows:

where and are, respectively, zeroth-order and first-order Bessel functions of the first kind and , , , , , , and refer to Table 3.

The second formula in (3) yields the surface displacement component of the half-infinite space under a vertical disturbance as

where is the dimensionless frequency and and are the functions of . The specific expression refers to Table 3.

2.2. Wave Displacement Field of a Circular Baseplate under a Simple Harmonic Load

Using the wave model for a dynamic half space presented in Section 2.1, the displacement field of an infinite half space with a circular surface under a simple harmonic load can be derived as shown in Figure 3.

Figure 3: Half-infinite space model for circular baseplate vibration.

A vibrational coupling model of a circular baseplate and the ground is shown in Figure 3, where , , , and denote, respectively, the displacement components along the horizontal axis , the vertical axis , and the direction . The displacements produced by a body wave are given as follows [19, 20]:

where and refer to Table 3.

The stress field of the circular baseplate-ground system under the action of a harmonic disturbance can be obtained from (5). To determine the vertical displacement component of the circular baseplate with , the distribution of the reaction force imposed on the bottom of the baseplate needs to be determined.

The force beneath the baseplate is determined by the stress distribution pattern beneath the baseplate. Three stress distribution patterns, which are (a) a static stiffness, (b) a uniform distribution, and (c) a parabolic distribution, are assumed (Figure 4). The stress boundary based on this three stress distribution patterns is obtained. The specific expression of refers to Table 3.

(a)

(b)

(c)

Figure 4: Distribution of reaction force imposed on the bottom of the baseplate.

The displacement at the center of the baseplate will be equal to the displacement of the ground, :

where and are the dimensionless frequencies. For each assumption of the stress distribution and refer to Table 3.

Redefining and , (6) can be written as

That is,

where is the driving signal, is the response, and is the impedance function of the vibrational system.

Equation (7) is the dynamic model of the infinite half space where circular surface is under a uniform stress .

2.3. Dynamic Model of Baseplate-Ground Vibrational System

For the baseplate-ground system, if nonlinear factors are neglected, the system vibration can be treated as a superposition vibration under a vertical disturbance force (Figure 5).

Figure 5: Decomposition and superposition of the baseplate-ground vibrational system.

According to (8), the reaction force on the bottom of the baseplate can be written as

Therefore, the differential equation describing the dynamics of the baseplate can be written as [19]

That is,

where and are, respectively, the mass and the displacement response of the baseplate, is the amplitude of the external force, and is the phase difference.

For harmonic vibrations, and and (11) can be written as

Under the assumptions

(12) can be simplified to

where and are the stiffness and damping . They are functions of the disturbance frequency . Expressions for and are given in Table 1 for the different modes of vibration.

Table 1: and for different modes of vibration.

Mode of vibration
Poisson’s ratio

0
0.25
0.5

Vertical vibration ()

Horizontal vibration ()

Swing vibration ()

—
—

Twisting vibration ()

3. Dynamic Model of Seismic Vibrator-Ground Coupling

3.1. Vibrational Coupling Model

A simple seismic vibrator-ground coupling system is shown in Figure 6 [3].

Figure 6: Mechanical calculation diagram for seismic vibrator-ground coupling system.

The dynamic equations for the coupling system can be written as

where and are, respectively, the masses of the reaction mass and the baseplate, and are, respectively, the stiffness and damping of the ground, and and are, respectively, the stiffness and damping of the hydraulic oil. is the sweep signal and can be written as follows for a linear signal:

where is the amplitude, is the sweep length, is the time, and are the start frequency and end frequency, and is the sweep bandwidth.

3.2. Solution of the Model

Equation (15) is solved using Newmark’s integral method, with the displacement and velocity at instant being expressed, respectively, as

where is the time interval. and are parameters related to the precision and stability, and the values are 0.5 and 0.25, respectively.

Using (17a) and (17b), the velocity and acceleration at instant can be expressed approximately by the responses at instant and the displacement at instant :

Substitution of (18a) and (18b) into (15) yields the following equation:

The system responses, including displacements, velocities, and acceleration, at instant can be obtained by substituting the solution of (19) into (18a) and (18b).

3.3. Verification of the Model

To verify the validity of the dynamic model, the force output of a real seismic vibrator-ground system is calculated and compared with the experiment data. The experiment was taken with an AHV-IV model-362 vibrator positioned on load cells using 5 Hz to 105 Hz in a 20-second linear sweep. A linear sweep was used to run a vibrator equipped with an advanced actuator which carries a modified hydraulic system and a new stiffer baseplate that is 2.5 times stiffer than the standard vibrator baseplate [21, 22]. This vibrator in Figure 7 was placed on load cells from which the true vibrator-ground force can be measured. All load cells were firmly mounted on a concrete pad and the target force was set at 42000 lbs or 186816 N [14]. In the experiment, it is assumed that the ground is homogeneous. Key parameters in the model and the experiment are listed in Table 2, and the parameters of the sweep signal in (16) are  kN,  Hz,  Hz, and  s.

Table 2: Parameters of seismic vibrator-ground system.

Hydraulic oil
Air spring
Vibrator vehicle
Reaction mass
Baseplate
Ground

Stiffness (N/m)

—
—
—

Damping (N⋅s/m)

—
—
—

Mass (kg)
—
—
28,000
2600
1380
—

Density (kg/m³)
—
—
—
—
—
2600

Area (m²)
—
—
—
—
2.64
—

Young’s modulus (N/m²)
—
—
—
—
—

Table 3

,

,

Volume strain

Rotational component

,

(a) Stress distribution which is a static stiffness

(b) Stress distribution which is a uniform distribution

(c) Stress distribution which is a parabolic distribution

(a) The stress distribution and which is a static stiffness
,

;

(b) The stress distribution and which is a uniform distribution
,

;

(c) The stress distribution and which is a parabolic distribution
,

Figure 7: The vibrator equipped with an advanced actuator located on load cells [14].

The time-history response of the output force from the model is compared with experimental results in Figure 8. Both model data and experimental results show that the vibrator output force decreases as the sweep frequency goes up. At very high frequencies, the torque motor current on the pilot servo-valve reaches its limit of 40 mA (milliampere). Often times, the instability of the servo-valve system occurs once the torque motor current is saturated and causes the poor vibrator performance [23]. The amplitude envelope in Figure 9 produced by model data shows a similarity to the amplitude envelope of experimental data although some discrepancy is still seen. This shows that the model is reasonable and can be used to estimate the impact of parameters in vibrator-ground system on the vibrator output force.

(a)

(b)

Figure 8: Comparison between theoretical and measured values of the output force: (a) output force time-history response from the present model; (b) output force time-history response from the weighted-sum method and experimental results [14].

Figure 9: Comparison of output force amplitude time-history responses.

4. Dynamic Characteristics Analysis

4.1. Modal Analysis

The characteristic equation for the dynamic equations (15) of the coupling system can be written as

The natural frequency of the system can be obtained by solving this characteristic equation and is found to be

The vibrational modal matrix of the coupling system can be written as

If a first-order vibration system is used to describe the vibrator actuator motion at low frequencies, the reaction mass and baseplate motions can be depicted as a synchronous vibration where the natural frequency is at 2.77 Hz (Figure 10(a)). The displacement amplitude of the reaction mass is much greater than the baseplate. This means that contribution to the vibrator output force from the baseplate motion can be neglected and decoupling hardly occurs.

(a)

(b)

Figure 10: First two vibrational modes of the coupled system: (a) first-order mode and (b) second-order mode.

As the frequency goes up, a second-order vibration system is capable of describing the reaction mass and baseplate motions where the natural frequency is at 488.7 Hz (Figure 10(b)). The displacement amplitude of the baseplate becomes greater than the reaction mass. As the vibration frequency approaches the natural frequency, the phase difference of the reaction mass and baseplate becomes large and goes to 180 degrees apart. This out-phase relationship will possibly cause the decoupling. It is important that the sweep frequency should not be close to this natural frequency.

4.2. Effects of System Parameters on Vibrational Output

The effects of reaction (inertial) mass and baseplate mass on the vibrational output are shown in Figures 11 and 12 which illustrate the effects of the masses of the reaction mass and baseplate on the vibrator output force, respectively. It can be seen from Figures 11 and 12 that an increase of the reaction mass results in an increase of the vibrator force amplitude while an increase of the mass of the baseplate results in a decrease of the vibrator force amplitude. These findings provide guidelines for the design of the vibrator actuator. It is authors’ opinion that the mass ratio of the reaction mass and the baseplate should be kept at 3 : 1.

Figure 11: Effect of reaction (inertial) mass on vibrational output.

Figure 12: Effect of baseplate mass on vibrational output.

Figure 13 displays how the elastic modulus of the ground affects the vibrator output force. It can be observed that increasing the elastic modulus of the ground will improve the vibrator output force. Increasing the elastic modulus means that the ground becomes stiff. Although it helps improve the amplitude of the vibrator output force, it reduces the force-energy that is transmitted to the deep ground. For example, a vibrator vibrating on a piece of limestone will produce an improved force amplitude. However, the force-energy radiated into deep ground will become poor.

Figure 13: Effect of elastic modulus on vibrational output.

4.3. Effects of System Parameters on Phase Difference

To obtain a vibrational output force of high quality, it is necessary that the phase difference between the sweep signal and the vibrational output be as small as possible. To investigate the influence of the system parameters on the phase difference, the following harmonic signal is used to replace the actual linear sweep signal shown in (16):

Figure 14 shows the vibrational output response under a harmonic sweep signal. It can be seen that, with increasing frequency of the sweep signal, the phase difference becomes larger.

(a)

(b)

Figure 14: Phase comparison between input signal and vibrational output: (a) sweep signal frequency  Hz and (b) sweep signal frequency  Hz.

To further investigate the influence of the frequency of the sweep signal on the phase difference, the load vector of the system dynamic (15) is replaced by the following expression:

The steady-state solution of (15) can be assumed to be

where and are the displacement state solutions for the vibrating plate and hammer and , , , and are undetermined coefficients.

Substituting , , and and , , and into (15) and separating the sine and cosine functions yield the following equations:

The undetermined coefficients , , , and , and thus and , can be determined by solving (26a), (26b), (26c), and (26d).

The force transmitted to the ground through the baseplate can be written as

Substituting into this equation then gives

where and are, respectively, the amplitude of the vibrational output and the phase difference between the sweep signal and the vibrational output.

Figure 15 shows that the effect of the reaction mass on the phase difference is significant, with greater reaction mass leading to a bigger phase difference, but, as shown in Figure 16, there is no obvious effect of baseplate mass on the phase difference.

Figure 15: Effect of reaction mass on phase difference for a baseplate mass of 2600 kg.

Figure 16: Effect of baseplate mass on phase difference for a reaction mass of 1380 kg.

Figure 17 shows that there is an obvious influence of baseplate area on phase difference, with greater baseplate area leading to a smaller phase difference. It can be seen from Figure 18 that the greater the ground stiffness, the greater the phase difference.

Figure 17: Effect of baseplate area on phase difference for a baseplate mass of 1380 kg and a reaction mass of 2600 kg.

Figure 18: Effect of ground stiffness on phase difference for a baseplate mass of 1380 kg and a reaction mass of 2600 kg.

5. Conclusions

The dynamic characteristics of a seismic vibrator-ground coupling system have been investigated by analyzing the system’s vibrational model and the influence of the main system parameters on the vibrational output and the phase difference between the input sweep signal and the vibrational output. The following conclusions can be drawn from the results obtained:

(1)The main vibrational modes of the coupling system are the first- and second-order modes with frequencies 2.77 Hz and 488.7 Hz, the latter being reverse movement of the reaction mass and baseplate will lead to an output phase difference of 180 degrees. This out-phase relationship will possibly cause the decoupling. Therefore, in an actual vibrator design, the frequency of the sweep signal should be lower than the second-order frequency of the system.(2)The design of the two most important components of the coupling system, namely, the baseplate and the reaction mass, should aim at improving the amplitude of the vibrational output and reducing the phase difference. It is found in this paper that an increase of the reaction mass results in an increase of the vibrator force amplitude while an increase of the mass of the baseplate results in a decrease of the vibrator force amplitude. These findings provide guidelines for the design of the vibrator actuator. The mass ratio of the reaction mass and the baseplate should be kept at 3 : 1.(3)The physical properties of the ground will affect the natural frequency of the coupling system and ultimately the quality of the output signal. Therefore, the sweep signal should be determined according to the physical properties of the ground and the modal analysis of the coupled system. That is, the frequency of the sweeping signal should be kept away from the frequency of second-order vibration of the vibrator-ground system.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

Z. Q. Huang, Y. P. Ding, Z. F. Tao, L. Hao, and G. Li, “Research status and development direction of vibroseis vibrator baseplate at home and abroad,” Chinese Journal of Oil Field Equipment, vol. 44, no. 6, pp. 1–5, 2015.

W. L. Chapman, G. L. Brown, and D. W. Fair, “The Vibroseis system: a high‐frequency tool,” Journal of Geophysical, vol. 46, no. 12, pp. 1657–1666, 1981.

J. J. Sallas and R. M. Weber, “The amplitude and phase response of a seismic vibrator’ by W. E. Lerwill,” Geophysical Prospecting, vol. 30, no. 6, pp. 935–938, 1982.

K. P. Allen, M. L. Johnson, and J. S. May, “High fidelity vibratory seismic (HFVS) method for acquiring seismic data,” in Proceedings of the SEG Technical Program Expanded Abstracts 1998, pp. 140–143.

A. V. Lebedev and I. A. Beresnev, “Nonlinear distortion of signals radiated by vibroseis sources,” Journal of Geophysical, vol. 69, no. 4, pp. 968–977, 2004.

N. Nagarajappa and D. Wilkinson, “Source measurement effect on high fidelity vibratory seismic separation,” Geophysical Prospecting, vol. 58, no. 1, pp. 55–68, 2010.

C. Saragiotis, P. Scholtz, and C. Bagaini, “On the accuracy of the ground force estimated in vibroseis acquisition,” Geophysical Prospecting, vol. 58, no. 1, pp. 69–80, 2010.

R. Zhang, Y. M. Chen, and S. L. Xie, “Analysis of dynamic characteristics of micro-vibroseis,” Equipment for Geophysical Prospecting, vol. 24, no. 2, pp. 91–95, 2014 (Chinese).

Z. Wei, “Pushing the vibrator envelope: extending low and high frequency limits,” First Break, vol. 26, no. 3, pp. 37–43, 2008.

Z. Wei, “Pushing the vibrator ground-force envelope towards low frequencies,” Geophysical Prospecting, vol. 57, no. 1, pp. 151–161, 2009.

Z. Wei, “Modelling and modal analysis of seismic vibrator baseplate,” Geophysical Prospecting, vol. 58, no. 1, pp. 19–32, 2010.

G. M. Zheng, X. M. Lei, H. Zhang, and Z. F. Chen, “Researches and simulation of vibroseis signal generator based on MATLAB GUI,” Geophysical and Geochemical Exploration, vol. 37, no. 1, pp. 160–164, 2013 (Chinese).

W. Sun, “Vibration output simulation of controlled source,” Chinese Journal of Geophysical Prospecting Equipment, vol. 22, no. 2, pp. 79–82, 2012.

Z. Wei, T. F. Phillips, and M. A. Hall, “Fundamental discussions on seismic vibrators,” Journal of Geophysical, vol. 75, no. 6, pp. W13–W25, 2010.

M. Hall, “Analysis of field tests with an improved hydraulic vibrator,” SEG Technical Program Expanded Abstracts, vol. 28, no. 1, pp. 104–108, 2009.

W. T. van Horssen and N. V. Gaiko, “On the transverse, low frequency vibrations of a traveling string with boundary damping,” Journal of Vibration and Acoustics, Transactions of the ASME, vol. 137, no. 4, Article ID 041004, 2015.

J. M. Tian and C. Z. Liang, “Foundation design of dynamic machine foundation on natural base,” Cement Technology, vol. 6, pp. 85–89, 2014 (Chinese).

R. J. Yan, Y. S. Wang, and Q. Y. Han, Introduction to Semi Space Theory of Dynamic Foundation, China Building Industry Press, Beijing, China, 1981.

H. J. Qian, D. M. Zhang, and J. X. Wang, Dynamic Machine Foundation Design, China Building Industry Press, Beijing, China, 1980.

B. D. Wu and L. B. Wu, “Analog method for the equivalent lumped system of foundation half-space theories and associated field measurements,” China Earthquake Engineering Journal, vol. 37, no. 4, pp. 1029–1036, 2015 (Chinese).

Z. Wei, “Reducing harmonic distortion on vibrators—stiffening the vibrator baseplate,” in Proceedings of the 70th EAGE Conference and Exhibition incorporating SPE EUROPEC 2008, Rome, Italy, June 2008.

Z. Wei, “How good is the weighted-sum estimate of the vibrator ground force?” The Leading Edge, vol. 28, no. 8, pp. 960–965, 2009.

Z. Wei, “Pushing the vibrator envelope—how high can we go?” in Proceedings of the 69th EAGE Conference and Exhibition incorporating SPE EUROPEC '07, London, UK, June 2007.