This website uses cookies to improve your experience. If you continue without changing your settings, you consent to our use of cookies in accordance with our cookie policy. You can disable cookies at any time.


Theoretical study on inversing the anisotropic permeability of VTI two-phase media by flexural waves in a fluid-filled borehole



We developed a theoretical method of estimating the vertical and horizontal permeability of a transversely isotropic two-phase formation with a vertical symmetry principle axis (VTI) by using the attenuation characteristic of the dipole-flexural waves in a fluid-filled borehole. The attenuation value of flexural waves was extracted from synthetic data. The inversion was operated by the least-squares method. Two main conclusions were obtained based on the analysis of the results of inversion implemented on multiple sets of synthetic waveform data. First, the simultaneous estimation of vertical and horizontal permeability was feasible for the slow VTI formation, which means its shear-wave velocity is slower than the acoustic velocity of fluid in a borehole that is different from the previous understanding in some studies. Second, the fast formation corresponded to the opposite case from the slow formation, the precision of the estimated horizontal permeability relates to the accuracy of the vertical permeability, but the vertical permeability is difficult to estimate by the attenuation of the dipole-flexural waves. We added random noise with different intensities to synthetic data in order to simulate the errors that might exist in actual logging data. We adopted Prony’s method to extract attenuation to eliminate the effect of noise to a certain degree.


Permeability is one of the most important parameters in geology to describe the production of an oil and gas reservoir. Accurately estimating permeability is one of the most basic tasks in oil and gas exploration and development, as well as one of the most important subjects of logging interpretation. Permeability has a relatively wide range of possible values and is difficult to measure. There is yet no single way with enough credibility to replace all other methods. The most prevalent current approach for determining permeability is via comprehensive interpretation and evaluation by several methods. So continuing to seek an effective new method to determine permeability is still a practically significant subject and a worthy direction for researchers to pursue.

It is generally believed that, compared to other measurement methods, the permeability value extracted from logging corresponds to the reservoir description scale and reflects in situ permeability better. Permeability is a mechanical parameter to describe fluid flow in two-phase media. The permeability measurement resulting from acoustic logging based on the principle of elastodynamics can be considered as a direct method. Based on the theoretical and field data analysis, the inference and prediction about the phase velocity and attenuation of the Stoneley waves in a borehole have correlation with formation permeability has been proposed by several researchers (Rosenbaum, 1974; White, 1983; Williams et al., 1984). Through a series of convincing laboratory experiments, and comparison with theory, Winkler et al. (1989) verify the validation of Biot’s model used to describe the dispersion and attenuation of the Stoneley waves in a borehole surrounded by permeable formation. The studies on estimating permeability quantitatively by applying Stoneley waves of acoustic logging were performed successively. Tang et al. (1991) realize rapid permeability estimation by the Stoneley wave in a simplified Biot-Rosenbaum model (Biot, 1956a; Rosenbaum, 1974) and took into consideration the influence of the logging tool size (Tang and Cheng, 1996). Wu et al. (1995) study the permeability inversion by the attenuation of the Stoneley wave based on the complete Biot-Rosenbaum model. Furthermore, the correct program for the impermeable intrinsic attenuation was proposed using the analysis of argillaceous by Wu et al. (1996). Liu and Johnson (1997) and Ticchelaar (1999) consider the effect of mud-cake on the Stoneley wave and achieved the quantitative inversion of permeability. Brie et al. (1998) realize the quantitative evaluation of formation permeability from Stoneley waves, taking into account the modification of the effect of mud and mud-cake properties. This methodology of estimating permeability by Stoneley waves was also used to process the actual field data (Pampuri et al., 1998).

In the last twenty years, the technology of acoustic logging has been developing rapidly that greatly expands the function of acoustic detection and provides a vast amount of measured data. Wang (1999) discovers that there is better correlation between the attenuation of the flexural wave and the permeability of two-phase medium formation, which is similar to the situation for Stoneley waves, and proposes the program of permeability inversion using the amplitude ratio of the flexural wave based on the study of forward modeling. The feasibility of their method has been verified by processing the synthetic signal and using practical examples, which offer a different way to get the permeability of formation from the Stoneley wave permeability. Modern devices used in acoustic logging are all combinations of monopole and cross-dipole array and other functional tools. So applying two sets of field data gotten from different acoustic logging methods to estimate identical reservoir parameters, such as shear velocity and permeability, should always be useful because they can confirm and complement each other.

However, the above results are all acquired in the case of isotropic formations. In fact, reservoirs and permeability are usually anisotropic. Therefore, it would be necessary in further study to estimate anisotropic permeability by acoustic logging. A transversely isotropic (TI) formation is one of the most common anisotropic formations because of the formation deposition and the directional distribution of the cracking. There are five independent elastic constants for a TI elastic solid, and there are three more constants that can be calculated for TI two-phase media. The TI media with the vertical symmetry axis (VTI) was studied earlier than other kinds of TI formations. Schmitt (1989) studies the problem of multipole acoustic fields in a fluid-filled borehole surrounded by a VTI two-phase formation.

Because of the demands of acoustic logging in a horizontal or a high-angle borehole surrounded by a TI formation, research on the numerical simulation of a 3D dipole transient acoustic field in TI media with arbitrary symmetry axis has been conducted (e.g., Liu and Sinha, 2000; He and Hu, 2009; Yan, 2011). The study on the forward model of acoustic logging in anisotropic media has been promoted, but there were few works on inversion parameters of anisotropic media. For example, there is still no way to quantitatively invert all the elastic modulus in either TI medium with the horizontal symmetry axis (HTI) or VTI medium. In addition, the detection and inversion of anisotropic permeability have barely been developed. One possible reason is that the dipole-flexural wave is not sensitive to vertical permeability (Schmitt, 1989). Recently, He et al. (2010) come to a similar conclusion through investigating the sensitivity coefficients of the flexural wave attenuation, and estimated the single horizontal permeability by simulated data based on a given vertical permeability. We also found that there is a big disparity of the sensitivity degree of vertical permeability to the flexural wave attenuation between the fast and slow formations (the slow formation means the velocity of a shear wave in the skeleton is slower than the acoustic velocity in the fluid of borehole VS<Vf the fast formation means the velocity of a shear wave in the skeleton is faster than the acoustic velocity in the fluid of the borehole VS>Vf). The sensitivity of flexural wave attenuation in a slow VTI two-phase formation to vertical permeability is significantly greater than it is in fast VTI two-phase formation. Furthermore, the slow formation is actually a kind of important common reservoir also. Therefore, the synchronous inversion of vertical and horizontal permeability for VTI two-phase formation is still possible and should be explored continuously, at least for the slow formation.

Detecting formation anisotropy is an important task of cross-dipole acoustic logging. Developing the functions of the cross-dipole acoustic logging from detecting only velocity anisotropy to determining permeability anisotropy is a natural and logical advance in methodology. In continuously cross-dipole acoustic logging, many kinds of anisotropic formations, including HTI, VTI, or even more complicated anisotropic formations, can all be encountered. Using cross-dipole acoustic logging to realize the inversion of anisotropic permeability is a more reasonable choice in practical application also. This does not exclude the effect of the monopole acoustic logging. In fact, the completion of the inversion of anisotropic formation parameters usually needs the combination of two methods.

In accordance with the discussion concerning the possibility of the simultaneous inversion of vertical and horizontal permeability as outlined above, we focus our study on the inversion of anisotropic permeability of VTI two-phase media by dipole-flexural waves, which is also a beginning for solving more complicated problems of anisotropic permeability.

As an important stage before practical application, this paper only describes theoretical methods and attempts to invert vertical and horizontal permeability simultaneously and emphasizes for slow formation. Several works realize permeability inversion by Stoneley waves in a borehole (Tang et al., 1991; Wu et al., 1995, 1996; Tang and Cheng, 1996; Brie, 1998), which were all based on the traditional model studied to a higher frequency range (Biot, 1956b; Johnson, 1987), and improve it being applicable to anisotropic (VTI) two-phase media (Schmitt, 1989; Zhang 1995). To simulate the influence of possible errors in practical logging data, we added random noise of different intensities to them. Finally, the vertical and horizontal permeability inversion for fast formations using synthetic data is also conducted and discussed.


Borehole configuration

The borehole is assumed to be a fluid-filled cylinder with the radius a. The density and acoustic speed of the fluid in the borehole are ρf and Vf, respectively. The medium outside the borehole is the transversely isotropic two-phase medium, whose symmetric axis is parallel to the borehole axis. The dipole source and receivers, which are used to excite and receive the acoustic field, are all located on the borehole axis (Figure 1).

Figure 1.

Figure 1. The configuration of the VTI two-phase formation and borehole.

The displacement potentials

We adopt the cylindrical coordinate system (r, θ, and z) centered on the dipole source and oriented along the borehole axis. Because the excitation frequency of the acoustic logging tool is obviously higher than seismic exploration frequency and covers a wider frequency range, Biot’s (1956b) high-frequency theory is adopted (Johnson, 1987). The dynamic differential equations in frequency domain are as follows:

where σ is the stress tensor for the skeleton and σ is the effective fluid stress for saturated fluid (Schmitt, 1989; Zhang, 1995). Here, ρ11, ρ12, and ρ22 are the mass coefficients tensor, and these coefficients are considered when the borehole is not uniform:
where I is the unit tensor and ρs and ρf represent the skeleton density and fluid density. Where ϕ is the porosity of media and E is the intrinsic tortuosity tensor, which is defined by
Symbol b is the second-order tensor, which is defined by
where η is the intrinsic viscosity and k is the intrinsic anisotropic permeability tensor defined by
where FJ(ω) is the viscous corrected vector factor, which should be considered when the viscous skin depth becomes comparable to the relevant pore size (Biot, 1956b; Johnson, 1987) that is taken for the expression by Johnson (1987):
where m is a nondimensional parameter, and usually takes a value from 8 to 12 (Johnson, 1987).

With the viscous correction factor, the viscosity is frequency dependent and so is permeability. Here, U and V represent displacements of the skeleton and fluid, respectively. Then, we define γ as follows:

where the superscripts H and V correspond to the horizontal and vertical components of the parameters.

With the time dependence factor eiωt, equation 1 can be expressed in the frequency domain as


The relationships of stress and strain in the two-phase VTI formation are as follows:

where A, N, C, F, and L are the elastic modulus of the TI formation, M and Q are for solid-fluid coupling, and R is for the fluid-phase. Here, eij and ε represent the strain components of the solid and fluid, respectively. The M, Q, and R are calculated by
where Ks and Kf are bulk compressional modulus of grains and fluid and Kf=Vf2ρf, respectively.

The displacement potentials are defined by the formulas below:

where φ and φ are compressional potentials, χ and χ are horizontal polarized shear potentials (SH-wave), and ψ and ψ are vertical polarized shear potentials (SV-wave).

Wave potentials

Based on the analysis above, we get the solutions of displacement potential equations by a set of derivations, (given in Appendix A), as


Equations 1213 are the solutions of SH-wave, and equations from 14 to 17 are the solutions of P–SV-wave (P- and SV-wave are coupled strongly in the transversely isotropic media). The In and Kn are the nth modified Bessel functions of the first and second kind. We use the symbol ν to represent the radial wavenumber in fluid (ν=kz2kf2), kz represents the axial wavenumber, and kf=ω/Vf. For the SH-wave (equations 12 and 13), χ0 is the weighting coefficient, αSH is the radial wavenumber, and a0=γ12H/γ22H (the meaning of γ is given in Appendix A). For the P–SV-wave (equations 1417), gi and αi(i=1,2,3) express weighting coefficients and radial wavenumber of quasi-P1(i=1), quasi-P2(i=2), and quasi-SV wave (i=3), respectively. The ai, bi, and ci(i=1,2,3) are unknown coefficients, which can be solved by substituting equations from 14 to 17 to the boundary conditions. The order of the multipole source is denoted by n, the flexural wave is excited by the dipole source, which corresponds to the condition of n=1, which is discussed below, so the subscript n is omitted. Based on equations 1217, the field equations A-12, A-23A-28 can be derived, and the derivation is given in Appendix A.

Then we can get the representations of displacement and stress, which are shown in Appendix A.

The representation of displacement potential of the fluid inside the borehole is given by

Here, r0 and θ0 are the separation and azimuth angle of the dipole source. Based on equation 18, the associated radial displacement and the stress of fluid-phase can be solved (equations A-29 and A-30), the specific process is given in Appendix A.

Wave propagation

In this paper, we consider the interface between the fluid-filled borehole and the VTI two-phase formation as a permeable wall only, so the boundary conditions at the position r=a of the cylindrical coordinate system are

where the superscript I and II represent the media inside and outside borehole, respectively. Here, T represents the total stress, which is defined by Tij=σij+σδij. We substitute the acoustic-field equations (from A-10A-12 and A-23A-30) into equation 19, and we obtain the following equation:
where the elements mij and b1,b2,b3 of the matrix above are given in Appendix B. The reflection coefficient A1 can be represented as

Then, the acoustic field excited by a dipole source in the borehole in frequency–wavenumber domain can be written as

here, F(ω) is the frequency spectrum of the source excitation pulse.

Inversion method

The flexural wave is a type of guided wave whose information is contained in the reflection coefficient. The poles of equation 21 correspond to all the modes of flexural waves. The amplitude of the flexural wave in the spatial-frequency domain can be calculated by the poles’ residue and can be written as


The flexural wave has attenuation characteristics because the VTI two-phase formation outside the borehole is dissipative (White, 1983; Tang and Cheng, 2004). Thus, the wavenumber of flexural wave is complex, which is expressed by kzfl.


where kR and α represent the propagation and attenuation of the flexural wave, respectively. It is conceivable that the attenuation can be extracted by amplitude ratio of the flexural wave in theoretical study. However, the practical logging data are more complicated; this method seems too be simple to extract accurate attenuation. To approach a practical application, we use Prony’s method to extract attenuation of the flexural waves. Prony’s method has become an effective way to estimate power spectral density (Kay and Marple, 1981). Moreover, many researchers applied Prony’s method to estimate slowness dispersion and attenuation from arrays of sonic logging waveforms (Parks et al., 1982, 1983; Lang, 1987; Ma et al., 2010).

First, we analyze the effect of the intrinsic anisotropic permeability of VTI two-phase media on the flexural wave attenuation by numerical simulation.

The elastic modulus of the VTI two-phase media are given in Table 1, these parameters are quoted from Schmitt (1989). We take three and four for the intrinsic horizontal and vertical tortuosity, and nine for m. Theoretically, the value of tortuosity relates to the value of permeability, but that effect is so little that can be ignored. Therefore, we always take these values in the calculation of different permeabilities. For reflecting the features of slow and fast formation and their anisotropy directly, the parts of characteristic wave velocities in the VTI two-phase medium are shown in Table 2; they include fast and slow P-waves, SV-, and SH-waves, respectively. These velocities are in the condition of 6-kHz frequency with the intrinsic vertical and horizontal permeability — both 300 mD. Both the fluid-phase in the two-phase formation and the fluid in the borehole are assumed to be water. The borehole radius and the dipole source separation were 0.1 and 0.01 m, respectively. In VTI two-phase formation, the permeability is related to the directions, so we make kH and kV represent the intrinsic horizontal and vertical permeability, respectively.

Table 1. The elastic modulus and characteristic wave velocities of the VTI two-phase medium (elastic modulus unit: 109N/m2, density unit: kg/m3, these parameters are quoted from Schmitt, 1989).

 Slow formationFast formation

Table 2. The vertical and horizontal velocity of compressional waves and shear waves (unit: m/s).

FormationDirection of propagationFast P-waveSlow P-waveSV-waveSH-wave
Slow formationVertical222011711481148
Fast formationVertical3676199.221032103

Figure 2 displays the relationship between the attenuation α(ω) of the flexural wave and different permeability values in the frequency domain. The intrinsic permeability values kH and kV are changed from 10 to 1000 mD in the slow and fast VTI two-phase formations. It shows from Figure 2 that the attenuation α(ω) increases along with the frequency. The α(ω) relates to intrinsic permeability values kH and kV. The effect of the permeability on α(ω) in the high-frequency area is greater than that in the low-frequency area. It could be found that the effect of kH on α(ω) is significantly greater than that of kV on α(ω), which is consistent to the conclusion that the attenuation is not sensitive to kV. Figure 2 also shows that there is a significant difference in the sensitivity between the slow and fast formations. The attenuation in the slow formation is much more sensitive to the vertical permeability than that in the fast formation. This is a favorable conclusion for estimating the vertical and horizontal permeability simultaneously in the slow formation. In this paper, we focus our attention on the intrinsic anisotropic permeability inversion in the slow formation that covers some of the known oil reservoirs. Meanwhile, we also study the intrinsic anisotropic permeability inversion by the dipole-flexural wave in fast VTI two-phase media as a comparative work.

Figure 2.

Figure 2. The attenuation of the flexural wave in the slow (a and b) and fast (c and d) VTI two-phase media. (a and c) The curves when the horizontal permeability kH is changed from 10 to 1000 mD and vertical permeability kV=300mD. (b and d) The curves when kV is changed from 10 to 1000 mD and kH=300mD.

The inversion analysis is then carried out. The attenuation of the flexural wave, which is extracted from the time-domain full waveforms by Prony’s method, is named measured value and denoted as D(ω,kH,kV), where kH and kV are the real intrinsic horizontal and vertical permeability values in the formation. The attenuation α(ω) calculated by equation 25 expressed by G(ω,kH,kV), in which kH and kV are the estimated intrinsic permeability values. Other parameters of the formation are given values. Then, D and G are substituted in the variance function:

as the target function, where N is the number of the frequency points in selection. The least-squares method of double unknown quantities is used to invert kH and kV.

Feasibility investigation

To ensure the smoothness of inversion and the adaptability of the method, the variation function Q was studied for different values of estimated intrinsic permeability (kH and kV). The vertical and horizontal permeability values of slow VTI two-phase formation are investigated, and they are assumed to be 400 and 500 mD, respectively. We selected 80 points in the frequency domain from 4 to 8 kHz (N=80), the range of kH is from 200 to 600 mD and that of kV is from 200 to 800 mD. The variation function Q is calculated when kH and kV are changed, and the result is shown in Figure 3.

Figure 3.

Figure 3. (a) Three-dimensional image of variance square error function Q when the estimated permeability kH and kV changes and (b) projection image of (a) in the plane of kH and kV.

From Figure 3a, we can see that the value of Q has many poles when kH and kV changes, and reaches the minimum value when kH=kH and kV=kV. It is clearer in its plane projection image in Figure 3b, the values of kH and kV in the minimum pole correspond to 400 and 500 mD, respectively. This is equal to the true value in the formation, which means the variation function Q reaches the global minimum when the estimated permeability is close to the true value. Through changing initial values of kH and kV to test the convergence of the inversion method, for different true values kH and kV, trying many times, it always converges smoothly to the true values, even though there are some local minimum in the neighborhood of the global minimum as shown in Figure 3. A possible reason is that the variation function Q of local minimums is too large to satisfy the least-squares method. In conclusion, the convergence of the inversion method is smooth and not sensitive to the initial values of kH and kV in quite a wide range. It is feasible to estimate permeability for the theoretical model by the least-squares method.

Therefore, the least-squares method is suitably used to invert the vertical and horizontal permeability simultaneously in slow and fast VTI two-phase media. The measured value of the attenuation can be extracted from the flexural wave of dipole full waveforms by Prony’s method. We observe that the flexural wave calculated by pole residue is basically the same as that extracted from the full waveforms, so the flexural wave calculated by equation 24 can be used as actual waveform as a theory study. The attenuation extracted from the full waveforms above could be seen as measured data in inversion study. Our aim is to numerically investigate the feasibility, adaptability, and accuracy of estimating the anisotropic permeability by the flexural wave attenuation in VTI two-phase formation.


In this section, the simultaneous inversion of vertical and horizontal permeabilities are performed on a series of synthetic waveforms data, and the effect of noise on the inversion results is also discussed numerically.

First, we estimate permeability by attenuations extracted from synthetic data of dipole-flexural waves. The inversion is implemented by the least-squares method, and its convergence is very smooth for any initial value.

Table 3 shows the inversion results of different vertical and horizontal permeability of the slow VTI two-phase formation whose elastic modulus are given in Table 1. Here, kH and kV are the true intrinsic permeability values of the simulation reservoir, and kH and kV are the inversion results. Meanwhile the relative errors of them are also shown in the third and sixth rows of Table 3. We infer that horizontal and vertical permeabilities of slow porous formations can be estimated simultaneously based on our previous investigation. The curves of different horizontal permeabilities clearly separate from each other (Figure 2a) and that of the different vertical permeabilities separate a little (Figure 2b). It can be seen that the results coincident with our deduction, the relative errors of inversed horizontal permeability are all less than 0.3%, whereas that of inversed vertical permeability are under 2% too. The inversing precision is very high, and the convergence of inversion is quite good. We notice that when the horizontal permeability is smaller than the vertical permeability, the relative error of the inversion horizontal permeability is usually large, such as the sixth row and eighth row of Table 3. Therefore, we preliminarily infer that the error of the inversion result may increase when the vertical permeability is much larger than the horizontal permeability. That means the inversion results of vertical and horizontal permeability are influencing each other; if one of them is inaccurate, the other would be far away from the true value, too.

Table 3. The inversion of vertical and horizontal permeability of the slow VTI two-phase formation using the attenuation extracted from waveforms received by five receivers.

kH (mD)kH (mD)Relative errorkV (mD)kV (mD)Relative error

The inversion discussed above is idealized and indicates that it is possible to invert vertical and horizontal permeability simultaneously in theory. However, at this point, we cannot affirm yet that this method is feasible in practical application. The actual logging data are much more complicated than the synthetic data because the theoretical model is always more simple than the actual formation, and there is a higher probability of systematic and accidental errors. To reflect the error of the actual data, we investigate the inversion results on the synthetic data with random noise.

We add 5% random noise to waveforms of the five receivers with equal spacing (0.5 m) in the same formation. We adopt Prony’s method to extract attenuation of flexural waves to eliminate the effect of noise. Generally, it is believed that, compared to other acoustical characteristics, velocity, for example, the accuracy of extracted attenuation is much more affected by noise. The extracted attenuation is unstable in low- and high-frequency areas, except around the central frequency range. We abandon these unstable data and only select the middle frequency area (5–6 kHz). It definitely promotes the inversion precision. We then invert several groups of vertical and horizontal permeability in the slow two-phase formation by the method discussed above. The results are shown in Table 4. We can see that when the true horizontal and vertical permeability is relatively small (50 mD), the relative error of the inversion horizontal permeability result is a little large, but the absolute error of which is acceptable. On the contrary, the inversion result of vertical permeability on the low permeability condition is far from true value. Low permeability causes small attenuation, which means the attenuation barely changes, whereas vertical permeability changes. As the true permeability increases, the relative errors of inversion results decrease rapidly. For horizontal permeability, the relative errors are all smaller than 10%, most of them are less than 4%. The precision can be accepted in practical applications. For vertical permeability, except the low permeable situation, the relative errors are acceptable too. Most of them are less than 10%. We also notice that, corresponding to Table 3, the inversion results’ error of the sixth and eighth rows are relatively large as well. This is an interesting phenomenon, and it partly proves our conclusion that the precision of inversion results of vertical and horizontal permeability interact somehow. That illustrates the importance of estimating horizontal and vertical permeability simultaneously.

Table 4. The permeability inversion results of slow VTI two-phase formation of five receivers adding 5% noise.

kH (mD)kH (mD)Relative errorkV (mD)kV (mD)Relative error

Furthermore, the receiver array with eight or 12 receivers is usually applied in a practical logging system, and the spacing is 0.15 m, as usual. We adopt a receiver array with eight receivers to invert the permeability, to approach the actual application. With spacing of 0.15 m instead of 0.5 m previously, the attenuation will be smaller and more difficult to extract. The inversion results are shown in Table 5. Compared with the results of five receivers, the accuracy of inversions are not as good as Table 3, especially for vertical permeability. In addition, as in Table 3, the relative error of inversion horizontal permeability is large, whereas that of inversion vertical permeability is large. We infer that the shorter spacing of receivers causes the larger inversion error of the results.

Table 5. The permeability inversion results of slow VTI two-phase formation of eight receivers adding 5% noise.

kH (mD)kH (mD)Relative errorkV (mD)kV (mD)Relative error

Finally, the vertical and horizontal permeability of a fast VTI two-phase formation are also inversed by theoretical simulation data, and the results are shown in Table 6. In Figure 2c, we can see that the attenuation curves of different horizontal permeabilities separate clearly, but that curves of different vertical permeabilities almost superpose in Figure 2d. The results are in accordance with the expectation that the inversion precision of horizontal permeability is high. It is notable that most of the vertical permeability results are credible also, these results greatly exceed our expectation. Though the method still needs improvement before practical application because the practical waveforms are more complicated than the synthetic waveforms, we could still get enlightening results. From Table 6, we can see that the errors of some horizontal permeability inversion are larger than others’, whereas errors of the corresponding vertical permeability are large too, and the relative error is high when vertical permeability value is larger than horizontal permeability (kV>kH). This illustrates that the precision of the inversion of horizontal permeability strongly relates to the vertical permeability. This conclusion is suitable for slow VTI two-phase formation too. In practical inversion of VTI two-phase formations, a priori value of vertical permeability would affect directly on the precision of inversion results of horizontal permeability when the vertical and horizontal permeability could not invert simultaneously.

Table 6. The permeability inversion results of fast VTI two-phase formation.

kH (mD)kH (mD)Relative errorkV (mD)kV (mD)Relative error


The estimation of anisotropic permeability of a reservoir is an interesting and challenging subject. Because the inversion of vertical permeability of the VTI two-phase formation using the dispersion and attenuation of dipole-flexural waves is considered impossible according to previous understandings (Schmitt, 1989; He and Hu 2009; He et al., 2010), the study in this paper is still focused on the simultaneous estimation of vertical and horizontal permeability of VTI two-phase formations in theory.

In the next step of the actual data processing, we have to solve the influence of a variety of factors on the permeability inversion. We think that the influence factors mainly include two aspects. One is the value assignment of other parameters except unknown permeabilities; their value can be obtained by petrophysics measurements and other wireline logs as well as both monopole and dipole acoustic logging, and the analysis of the sensitive coefficient to some quantity can usually simplify treatment process. The other one is the change or damage of the physical environment, among which the influence of the mud cake and impairment layer of a borehole is a particularly important and perplexing problem. The flexural waves with horizontal polarization propagating along the skeleton of a poroelastic formation are different from the character of Stoneley waves; therefore, a new study on this topic should be developed. Actually, these problems can only be solved combining with the actual data processing, so making a feasible strategy (technical route) of the inversion would be necessary.


In this paper, we carry out a theoretical and numerical investigation of the simultaneous estimation of the intrinsic vertical and horizontal permeability for VTI two-phase (Biot) formations using the attenuation characteristic of the dipole-flexural waves in a fluid-filled borehole. The attenuation value of flexural waves is extracted by Prony’s method from synthetic waveforms data of dipole-flexural waves in frequency domain. Based on the analysis of the extreme value distribution of variance functions structured by model function and measuring data, the inversion of vertical and horizontal permeability is implemented using the least-squares method; the convergence of the inversion for multiple sets of parameters are all smooth.

For the slow VTI two-phase formation, the precision of inversion results for multiple sets of permeability are all higher, which means it is feasible to estimate vertical and horizontal permeability simultaneously for a slow VTI two-phase formation, at least in theory. The influence of error that is possibly presented in field data on the inversion precision is investigated numerically by adding random noise to synthetic data. Noise largely effects the accuracy of extracted attenuation, especially on the two poles of the frequency band. Thus, we abandon the data of the low- and high-frequency areas to eliminate the noise effect.

For the fast VTI two-phase formation, the inversion of vertical and horizontal permeability is tentatively carried out also for synthetic data. It is found that the inversion precision of horizontal permeability strongly relates to the precision of vertical permeability. This is a new understanding, which means that although the attenuation of dipole-flexural waves is not very sensitive to vertical permeability of a fast formation, it is not yet expected that a precise result of horizontal permeability can be obtained by single parameter inversion, no matter how much the true value of vertical permeability is.


This work is supported by the National Nature Science Foundation of China (grant nos. 11134011 and 11174322). We thank G. Harrison for his earnest help in corrections for language usage.


Six displacement potential equations can be obtained by equations 811:


The displacement potential equations of the SH-wave are equations A-1 and A-2, so we can easily get equations 12 and 13 from equations A-1 and A-2.

The P–SV-wave displacement potential function are equations A-3A-6; we can take the trial solution of them in the form below:


Equations 1114 can be solved as follows:


The coefficient determinant of equations A-8A-11 must be zero to ensure R1, R2, R3, and R4 are nonzero. Then we get

where di(i=1,,10) are in the form below:
There are always three solutions of α in equation A-12 corresponding to α1,α2 and α3, respectively. For a given αi(i=1,2,3), only one independent value of R1,R2,R3,R4 exists. Meanwhile α1,α2 and α3 correspond to quasi-P1, quasi-P2, and quasi-SV, respectively. Because P1 and P2 waves are coupled with SV-waves strongly in TI media, we can gain the equations from 13 to 17.

Then, based on the six potential equations (equations 1217), the displacement and stress of the formation outside the borehole are in the following form:


Here, T represents the total stress tensor, which is defined by Tij=σij+σδij.

For the fluid field in the borehole, the associated radial displacement Urf and the stress Pf can be derived by the potential function (equation 18). Given that Pf=ω2ρfϕf and Urf=ϕf, we determine the representations of fluid displacement and stress inside the borehole as follows:





  • Biot, M. A., 1956a, Theory of propagation of elastic waves in a fluid-saturated porous solid. I: Lower frequency range: Journal of the Acoustical Society of America, 28, 168–178, doi: 10.1121/1.1908239.JASMAN0001-4966CrossrefWeb of ScienceGoogle Scholar
  • Biot, M. A., 1956b, Theory of propagation of elastic waves in a fluid-saturated porous solid. II: Higher frequency range: Journal of the Acoustical Society of America, 28, 179–191, doi: 10.1121/1.1908241.JASMAN0001-4966CrossrefWeb of ScienceGoogle Scholar
  • Brie, A., T. Endo, D. L. Johnson, and F. Pampuri, 1998, Quantitative formation permeability evaluation from Stoneley waves: SPE Reservoir Evaluation and Engineering, 3, 109–117.SREEFG1094-6470Web of ScienceGoogle Scholar
  • He, X., and H. S. Hu, 2009, Borehole flexural modes in anisotropic formations: The low-frequency asymptotic velocity: Geophysics, 74, no. 4, E149–E158, doi: 10.1190/1.3141442.GPYSA70016-8033AbstractWeb of ScienceGoogle Scholar
  • He, X., H. S. Hu, and W. Guan, 2010, Fast and slow flexural waves in a deviated borehole in a homogeneous or layered anisotropic formation: Geophysical Journal International, 181, 417–426, doi: 10.1111/j.1365-246X.2010.04503.x.GJINEA0956-540XCrossrefWeb of ScienceGoogle Scholar
  • Johnson, D. L., J. Koplik, and R. Dashen, 1987, Theory of dynamic permeability and tortuosity in fluid-saturated porous media: Journal of Fluid Mechanics, 176, 379–402, doi: 10.1017/S0022112087000727.JFLSA70022-1120CrossrefWeb of ScienceGoogle Scholar
  • Kay, S. M., and S. L. MarpleJr., 1981, Spectrum analysis — A modern perspective: Proceedings of the IEEE, 69, 1380–1419, doi: 10.1109/PROC.1981.12184.IEEPAD0018-9219CrossrefWeb of ScienceGoogle Scholar
  • Lang, S. W., A. L. Kurkjian, J. H. McClellan, C. F. Morris, and T. W. Parks, 1987, Estimating slowness dispersion from arrays of sonic logging data: Geophysics, 52, 530–544, doi: 10.1190/1.1442322.GPYSA70016-8033AbstractWeb of ScienceGoogle Scholar
  • Liu, H., and D. L. Johnson, 1997, Effects of an elastic membrane on tube waves in permeable formations: Journal of the Acoustical Society of America, 101, 3322–3329, doi: 10.1121/1.418347.JASMAN0001-4966CrossrefWeb of ScienceGoogle Scholar
  • Liu, Q. H., and B. K. Sinha, 2000, Multipole acoustic waveforms in fluid-filled boreholes in biaxially stressed formations: A finite-difference method: Geophysics, 65, 190–201, doi: 10.1190/1.1444710.GPYSA70016-8033AbstractWeb of ScienceGoogle Scholar
  • Ma, J., P. J. Matuszyk, R. K. Mallan, C. Torres-Veredin, and B. C. Voss, 2010, Joint processing of forward and backward extended Prony and weighted spectral semblance methods for robust extraction of velocity dispersion data: Presented at SPWLA 51st Annual Logging Symposium.Google Scholar
  • Pampuri, F., M. Rovellini, A. Brie, and T. Fukushima, 1998, Effective evaluation of fluid mobility from Stoneley waves using full Biot model inversion: Two case histories: Presented at SPE Annual Technical Conference and Exhibition, SPE 49132.CrossrefGoogle Scholar
  • Parks, T. W., J. H. McClellan, and C. F. Morris, 1983, Algorithms for full-waveform sonic logging, in Proceedings of the Second Acoustics, Speech, and Signal Processing Workshop on Spectral Estimation: IEEE, 186–191.Google Scholar
  • Parks, T. W., C. F. Morris, and J. D. Ingram, 1982, Velocity estimation from short-time temporal and spatial frequency estimates, in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing, ICASSP’82, vol. 7, IEEE, 399–402.CrossrefGoogle Scholar
  • Rosenbaum, J. H., 1974, Synthetic microseismograms: Logging in porous formations: Geophysics, 39, 14–32, doi: 10.1190/1.1440407.GPYSA70016-8033AbstractWeb of ScienceGoogle Scholar
  • Schmitt, D. P., 1989, Acoustic multipole logging in transversely isotropic poroelastic formations: Journal of the Acoustical Society of America, 86, 2397–2421, doi: 10.1121/1.398448.JASMAN0001-4966CrossrefWeb of ScienceGoogle Scholar
  • Tang, X. M., and A. Cheng, 2004, Quantitative borehole acoustic methods: Elsevier Science Publishing Co., Inc.Google Scholar
  • Tang, X. M., and C. H. Cheng, 1996, Fast inversion of formation permeability from Stoneley wave logs using a simplified Biot-Rosenbaum model: Geophysics, 61, 639–645, doi: 10.1190/1.1443993.GPYSA70016-8033AbstractWeb of ScienceGoogle Scholar
  • Tang, X. M., C. H. Cheng, and M. N. Toksöz, 1991, Dynamic permeability and borehole Stoneley waves: A simplified Biot-Rosenbaum model: Journal of the Acoustical Society of America, 90, 1632–1646, doi: 10.1121/1.401904.JASMAN0001-4966CrossrefWeb of ScienceGoogle Scholar
  • Tichelaar, B. W., H. Liu, and D. L. Johnson, 1999, Modeling of borehole Stoneley waves in the presence of skin effect: Journal of the Acoustical Society of America, 105, 601–609, doi: 10.1121/1.426250.JASMAN0001-4966CrossrefWeb of ScienceGoogle Scholar
  • Wang, K. X., J. Ma, X.-Y. Wu, and B.-X. Zhang, 1999, Determination of permeability from flexural waves in dipole acoustic logging: 69th Annual International Meeting, SEG, Expanded Abstracts, 33–36.AbstractGoogle Scholar
  • White, J. E., 1983, Underground sound: Application of seismic waves: Elsevier Publishing Company.Google Scholar
  • Williams, D. M., J. Zemanik, F. A. Angona, C. L. Denis, and R. L. Caldwell, 1984, The long spaced acoustic logging tool: Presented at 25th Annual Logging Symposium, SPWLA.Google Scholar
  • Winkler, K. W., H.-L. Liu, and D. L. Johnson, 1989, Permeability and borehole Stoneley waves: Comparison between experiment and theory: Geophysics, 54, 66–75, doi: 10.1190/1.1442578.GPYSA70016-8033AbstractWeb of ScienceGoogle Scholar
  • Wu, X. Y., and K. Kang, 1996, Estimation of permeability from attenuation of the Stoneley wave in a borehole: 66th Annual International Meeting, SEG, Expanded Abstracts, 174–177.AbstractGoogle Scholar
  • Wu, X. Y., K. X. Wang, L. Guo, S.-C. Yu, and Q.-D. Dong, 1995, Inversion of permeability from the full waveform acoustic logging data: Chinese Journal of Geophysics, 38, 224–231.0898-9591Google Scholar
  • Yan, S. G., R. L. Song, W. G. Lv, K. X. Wang, and J. Ma, 2011, The finite-difference simulation of the porous cylindrical pipe simulating source and formation of axis-asymmetry transverse isotropy medium: Chinese Journal of Geophysics, 54, 2412–2418.0898-9591Google Scholar
  • Zhang, B. X., K. X. Wang, and Q. D. Dong, 1995, Acoustic multipole logging in transversely isotropic two-phase medium: Journal of the Acoustical Society of America, 97, 3462–3472, doi: 10.1121/1.412432.JASMAN0001-4966CrossrefWeb of ScienceGoogle Scholar