This is a continutation of Rock Physics of Geopressure and Prediction of Abnormal Pore Fluid Pressures Using Seismic Data (Part 1).

4. Seismic-to-pressure conversion

The sand/clay acoustic model for shaley sandstones, developed by Carcione et al. (2000) yields the seismic velocities as a function of clay (shale) content, porosity, saturation, dry-rock moduli, and fluid and solid-grain properties. As stated before, the large change in seismic velocity is mainly due to the fact that the dry-rock moduli are sensitive functions of the effective pressure, with the largest changes occurring at low differential pressures. The major effect of porosity changes is implicit in the dry-rock moduli. Explicit changes in porosity and saturation are important but have a lesser influence on wave velocities than changes in the moduli. This is due to the fact that the moduli are highly affected by the contact stiffnesses between grains. In this sense, porosity-based methods can be highly unreliable.

In order to use the theory to predict pore pressure, we need to obtain the expression of the dry-rock moduli versus effective pressure. The calibration process should be based on well, geolvogical and laboratory data, mainly sonic and density data, and porosity and clay content inferred from logging profiles.

Let us assume a rock at depth z. The lithostatic or confining pressure can be obtained by integrating the density log. We have that

Equation 4-1

where ρ is the density and g is the acceleration of gravity. Furthermore, the hydrostatic pore pressure is approximately given by ph = gρwz, where ρw is the density of water. As a good approximation (Prasad and Manghnani, 1997), compressional and shear-wave velocities and bulk and shear moduli depend on effective pressure pe = pnp, where p is the pore pressure and n is the effective stress coefficient, which can be different for velocities and moduli. Note that the effective pressure equals the confining pressure at zero pore pressure.

Fig. 4-1
Figure 4.1: Dry-rock bulk modulus of several reservoir rocks for different confining pressures, compared to the HS upper bounds. The solid line represents the analytical curve (4.4). The data (z1-z23) is taken from table 3.1 of Zimmerman (1991).

In general, n is approximately linearly dependent on the differential pressure pd = pcp in dynamic experiments (Gangi and Carlson, 1996; Prasad and Manghnani, 1997). Therefore, we assume n = n1 – n2pd = n1 – n2(pc – p). This dependence of n versus differential pressure is in good agreement with the experimental values, corresponding to the compressional velocity obtained by Christensen and Wang (1985) and Prasad and Manghnani (1997). It is clear that to obtain n1 and n2, we need two evaluations of n at different pore pressures, preferably a normally pressured well and an overpressured well. If one well or equally pressured wells are available, the algorithm provides an average value for n.

Calibration of the model

Ideally, a precise determination of n requires laboratory experiments on saturated samples for different confining and pore pressures. However, even this “laboratory” n does not reflect the behaviour of the rock at the in-situ conditions, due to two main reasons. First, laboratory measurements of wave velocity are performed at ultrasonic frequencies, and second, the in-situ stress distribution is different from the stress applied in the experiments. In the absence of laboratory data, or for shales, we perform the following steps with the data available from a calibration well (see Figure 4.2):

  1. We use Krief et al.’s model (Krief et al., 1990) modified by Goldberg and Gurevich (1998) to obtain the dry-rock moduli as a function of depth, using the clay-content and porosity profiles.
  2. We calculate the upper limit of the dry-rock moduli — at infinite effective pressure — using the Hashin-Shtrikman (HS) upper bounds (Hashin and Shtrikman, 1963), and assume an exponential law for the bulk moduli versus effective pressure.
  3. We compute the exponential coefficients (see equation (4.3)) using the values of the moduli obtained in step 1, the confining and (measured) pore pressures, and the effective stress coefficients predicted by Biot’s theory.
  4. We obtain the in-situ effective stress coefficients by fitting the theoretical wave velocities (Carcione et al., 2000) to the soniclog wave velocities, using the dry-rock moduli versus effective pressure obtained in steps 2 and 3, and n as a fitting parameter. The effective stress coefficient versus pore pressure, corresponding to the same geological unit, is obtained by using the linear law n = n1 – n2pd.

It is important to point out here that the HS bounds and Biot’s effective stress coefficient do not depend on the size and shape of the grain and pores. In this sense, the model has a general character. The only conditions are linearity, isotropy and the low-frequency approximation.

Fig. 4-2
Figure 4.2: Flow chart illustrating the pore pressure prediction method from seismic data.

If sandstone cores are available, we proceed as follows:

  • The upper limits (infinite confining pressure) and exponential coefficients of the moduli are obtained by fitting the dry-rock moduli, which are calculated from the dry-rock wave velocities, while n is obtained from experiments on saturated samples for different confining and pore pressures (Carcione and Gangi, 2000a,b).
  • Continue with step 4 of the preceding list. This step should improve the determination of n, estimated with laboratory experiments in step 1.

More precisely, we use the following data of the study area, to calibrate the model and obtain the effective-stress-coefficient profile for the formations under consideration:

  • An estimation of the porosity profile, ϕ(z) , to use in Krief et al’s model (see below) and in the sand/clay acoustic model (from a series of logs using artificial neural networks (Helle et al, 2001)).
  • An estimation of the clay-content profile C(z), to use in Krief et al’s model (see below) and in the sand/clay acoustic model (shale volumes obtained from well logs using neural networks (Helle et al., 2001)).
  • Direct measurements of pore pressure, p(z) from repeat formation tests (RFT) and/or mud weights provided by the mudlogging operator).
  • Sonic-log information, that is, the P-wave and S-wave velocity profiles, VP(z) and VS(z), used to obtain nK and nμ for the whole range of effective pressures by fitting the theoretical wave velocities to VP and VS. Here, nK and nμ are the effective stress coefficients corresponding to the dry-rock bulk and shear moduli, respectively.

No laboratory experiments available

Firstly, we consider the model of Krief et al. (1990), to obtain an estimation of the dry-rock moduli, Ksm , μsm (sand matrix), and Kcm, μcm (clay matrix) versus porosity and clay content. The porosity dependence of the sand and clay matrices should be consistent with the concept of critical porosity (Mavko et al., 1998, p. 244) since the moduli should vanish above a certain value of the porosity (usually from 0.4 to 0.5). This dependence is determined by the empirical coefficient A (see equation (4.2)). This relation was suggested by Krief et al. (1990) and applied to sand/clay mixtures by Goldberg and Gurevich (1998). The bulk and shear moduli of the sand and clay matrices are respectively given by

Equation 4-2

where Ks and μs are the bulk and shear moduli of the sand grains, and Kc and μc those of the clay particles. Krief et al. (1990) set the A parameter to 3 regardless of the lithology, and Goldberg and Gurevich (1998) obtain values between 2 and 4, while Carcione et al. (2000) use A = 2. Alternatively, the value of A can be estimated by using regional data from the study area. We use a general form of Goldberg and Gurevich’s equation. Experimental data is fitted in Carcione et al. (2000), showing that the model has been successfully tested. The model is not based on a dual porosity theory, but there is only one (connected) porosity. The clay moduli are taken from fit to experimental data in Goldberg and Gurevich’s paper.

Secondly, we assume the following functional form for the dry-rock moduli as a function of effective pressure:

Equation 4-3

where α(z), and p*(z) should be obtained (for each moduli) by fitting Krief et al’s expressions (4.2). The effective pressure at depth z is assumed to be pe = pc – n0p, where pc is given by equation (4.1), p is given, and n0 is a first estimation of the effective stress coefficient, based on Biot’s theory (Todd and Simmons, 1972; Zimmerman, 1991, p. 33).

Since there are two unknown parameters (α(z), and p*(z)) and one value of M for each depth, α(z) is assumed to be equal to the Hashin-Shtrikman (HS) upper bounds (Hashin and Shtrikman, 1963, Mavko et al., 1998, p.106). Note that the HS lower bounds are zero, and that the Voigt bounds are (1 – ϕ) Ks and (1–ϕ)μs, respectively. For quartz grains with clay, Ks = 39 GPa and μs = 33 GPa (Mavko et al., 1998, p. 307), and if the limit porosity is 0.2, the HS upper bounds for the bulk and shear moduli are 26 GPa and 22 GPa, compared to the Voigt upper bounds 31 GPa and 26 GPa, respectively. However, the HS bounds are still too large to model the moduli of in-situ rocks. These contain clay and residual water saturation, inducing a chemical weakening of the contacts between grains (Knight and Dvorkin, 1992; Mavko et al., 1998, p. 203).

Figure 4.1 shows the dry-rock bulk modulus of several reservoir rocks for different confining pressures (Zimmerman, 1991, p. 29, Table 3.1), compared to the HS upper bounds. The solid line represents the analytical curve (4.4). On the basis of these data we apply a constant weight factor β = 0.8 to the HS bounds — due to the softening effects.

Equation (4.3) for the sand matrix can be written as

Equation 4-4


Equation 4-5

where we assume that the effective stress coefficients are given by Biot’s expressions (Todd and Simmons, 1972),

Equation 4-6

and pe depends on the specific coefficient. The last equality results from the third equation(4.2), and the second equation is merely an extension of the first to the case of shear deformations.

Laboratory experiments available

The evaluation can be improved if laboratory data of dry-rock P-wave and S-wave velocity are available, and serves to constrain the values of α and p* in equation (4.3). These sparse calibration points are based on sandstone or shaley sandstone cores, since dry measurements in shales are practically impossible to perform. The seismic bulk moduli Ksm and μsm versus confining pressure can be obtained from laboratory measurements in dry samples. If VP(dry) and VS(dry) are the experimental compressional and shear velocities, the moduli are given approximately by

Equation 4-7

where ρs is the grain density. We recall that Ksm is the rock modulus at constant pore pressure, i.e. the case when the bulk modulus of the pore fluid is negligible compared with the dry-rock bulk modulus, as for example air at room conditions. Then, we perform experiments on saturated samples for different confining and pore pressures, to obtain the effective stress coefficient n. Because these experiments yield the P-wave and S-wave velocities, and the effective stress coefficients of wave velocity and wave moduli may differ from each other, we obtain n for

Equation 4-8

where K and μ are the undrained moduli.

Calculation of the effective stress coefficients

The last step of the calibration process is to consider equation pe = pc – np and obtain the effective stress coefficients nK(z) and nμ(z) by fitting the theoretical velocities to the corresponding sonic-log P-wave and S-wave velocities, by using expressions (4.3). First, we obtain nμ by fitting the S-wave velocity, because this velocity only depends on μsm, and then, we obtain nK by fitting the P-wave velocity. If shear-wave velocity data is not available, we assume nμ = n, i.e., the Biot estimate. From the values of n obtained at the two wells, we obtain the linear law for the geological unit under investigation. The values of clay content and porosity away from the wells are assumed to be equal to those of the nearest well. Interpretation is required to follow the geological units laterally, as a function of depth, so that the n profiles can be properly extrapolated. In this study, the clay-matrix modulus is simply given by Krief et al’s expressions, respectively, with no explicit dependence on pressure.

Pore-pressure calculation

Finally, the seismic velocity, derived from velocity analysis and inversion techniques, can be fitted with the theoretical velocities by using the pore pressure as fitting parameter. The theory (Carcione et al., 2000) allows us to introduce different kinds of information explicitly, such as composition (clay content), fluid saturation, porosity, permeability and viscosity (Pham et al., 2002). Before dealing with the seismic data we should test the above procedure in a nearby over-pressured well. The pore pressure prediction flow chart is shown in Figure 4.2 (Carcione et al., 2001, 2002).

AVO-based verification

In some cases, velocity information alone is not enough to distinguish between a velocity inversion due to overpressure and a velocity inversion due to pore fluid and lithology (e.g., base-of-salt reflections (Miley, 1999; Miley and Kessinger, 1999)). There are cases, where overpressuring is not associated with large velocity variations, as in smectite/illite transformations. Best et al. (1990) use AVO analysis to treat these cases. Modeling analysis of AVO signatures of pressure transition zones are given in Miley (1999), Miley and Kessinger (1999), Carcione (2000a,c) and Tinivella et al. (2001). This type of analysis should complement the present prediction method on the basis of geological information of the study area.

5. Application to a North Sea gas field

We consider the Tune Field area in the Viking Graben of the North Sea. This basin is 170-200 km wide, and represents a fault-bounded north-trending zone of extended crust, flanked by the mainland of western Norway and the Shetland platform. The area is characterized by large normal faults with north, northeast and northwest orientations which define tilted blocks. Such blocks contain the sequences present within the well used for this study. The main motivation for selecting this area is the fact that highly overpressured compartments were identified by drilling, and that higher overpressure is expected in future wells down the flank side towards the central Viking Graben. A detailed analysis of the fault sealing and pressure distribution in Tune Field is given by Childs et al. (2002).

Figure 5.1 displays a structural time map of Top Ness, showing the pressure compartments and the locations of three wells. Well 2 and 3 are overpressured (by about 15 MPa) and well 1 has almost normal (hydrostatic) pore pressure. The dashed line indicates the location of the seismic section shown in Figure 5.2. The calibration well (well 1) is an exploration well drilled to a depth of 3720 m (driller’s depth) to test the hydrocarbon potential of the Jurassic Brent Group. The well includes reservoir rocks of the Tarbert and Ness formations. The Tarbert sands are the target units considered in the present study.

Fig. 5-1
Figure 5.1: Structural time map of Top Ness (base reservoir) showing the pressure compartments in the study area. Well 2 and 3 are overpressured (by about 15 MPa). In well 1 there is almost normal (hydrostatic) pore pressure. The dashed line indicates the location of the seismic section shown in Figure 5.2. Notice that well 3 is highly deviated as indicated by the well location at surface (east) and at depth (west).

The 3-D marine seismic data was acquired by using a system of 6 streamers of 3 km length with a group interval of 12.5 m and cross-line separation of 100 m. The shot spacing was 25 m and the sampling rate 2 ms. The conventional stacked section is displayed in Figure 5.2, where the location of the wells is shown.

Fig. 5-2
Figure 5.2: Seismic section through Tune wells (Figure 5.3) showing the location of the Top Tarbert -Top Ness interval. The mean reservoir fluid pressures are indicated. The depths of interest are between Top Tarbert (green) and Top Ness (pink).

Figure 5.3 shows pressure and formation data (porosity ϕ, clay content C, density ρ, water saturation Sw and sonic-log velocities VP and VS) for the Tune wells. Note that well 1 is water bearing with moderate pore pressures while wells 2 and 3 are gas bearing and overpressured. Reservoir properties and fluid saturation are derived from wireline data using the neural net approach of Helle et al. (2001) and Helle and Bhatt (2002).

Fig. 5-3
Figure 5.3: Pressure and formation data (porosity ϕ clay content C, density ρ, water saturation Sw and soniclog velocities VP and VS for the Tune wells.

Velocity determination by tomography of depth migrated gathers

Recent advances in depth migration have improved subsurface model determination based on reflection seismology. Subsurface imaging is linked to velocity, and an acceptable image can be obtained only with a highly accurate velocity field. It has been recognized that prestack migration is a powerful velocity analysis tool that yields better imaging results than poststack migration in complicated structures. The basic assumption underlying the velocity determination methods based on prestack migration is that when the velocity is correct, all the migrations with data in different domains (e.g. different offset, different shots, different migration angles etc.) must yield a consistent image. In order to obtain the velocity field, we use the seismic inversion algorithm described by Koren et al. (1998).

We start with an initial model based on the depth converted time model, using a layer velocity cube based on conventional stacking velocities, and the interpreted timehorizons from the Tune project. Line by line, we perform the 3-D prestack depth migration using the initial velocity model and an appropriate aperture (3 km x 3 km at 3 km depth) in the 3-D cube. Through several iteration loops the model is gradually refined in velocity and hence depth. Each loop includes re-interpretation of the horizons in the depth domain, residual moveout analysis and residual moveout picks in the semblance volume. This is performed for each reflector of significance, starting at the seabed and successively stripping the layers down to the target.

The tomography considers, 1) an initial velocity model and, 2) the errors as expressed by the depth gather residual moveout and the associated 3-D residual maps. From these two inputs a new velocity model is derived where the layer depths and layer velocities are updated iteratively in order to yield flat gathers. The refined model is derived using a tomographic algorithm that establishes a link between perturbation in velocity and interface location, and traveltime errors along the common reflection point (CRP) rays traced across the model. CRP rays are ray pairs that obey Snell’s law and emanate from points along the reflecting horizon, arriving at the surface with predefined offsets, corresponding to the offset locations for the migrated gathers. Each pair establishes a relationship between the CRP and the midpoint of the rays at the surface. Depth errors, indicating the difference in depth of layer images and reference depth, are picked on the migrated gather along the horizon and converted to time errors along the CRP rays. The equations relating the time errors to changes in the model are solved by a weighted least squares technique.

Fig. 5-4
Figure 5.4: Seismic velocity map (left), overpressure prediction (center) and difference in overpressure due to gasbearing Tarbert and water-bearing Tarbert (right).

The final model consists of seven layers, i.e., the sea water layer, seabed-Top Diapir (clay diapirism is a characteristic feature of Tertiary throughout the area), Top Diapir-Top Balder, Top Balder-Top Cretaceous, Top Cretaceous-Base Cretaceous, Base Cretaceous-Top Tarbert, and the target layer, Top Tarbert-Top Ness. The velocity maps for Top Tarbert -Top Ness is shown in Figure 5.4 (left panel), where the well locations are indicated. The Cretaceous layer velocity and the depth to Base Cretaceous reveal a remarkable similarity, i.e., where the Cretaceous is deep the velocity is high and where the Cretaceous is shallow the velocity is low, indicating that the velocity of Cretaceous is essentially governed by the overburden (e.g. compaction). Base Cretaceous, Top Tarbert and Top Ness depth maps display the more dramatic geometry of Upper Jurassic, where the Tarbert Formation is completely eroded in the northwest and where patches of missing Ness are present along the major faults. These erosion features are, of course, well reflected in the layer velocities since here the layer thickness tends to zero and layer velocity is hence not defined. Structural features are well displayed in the velocity maps of the Top Tarbert and Top Ness. The Base Cretaceous- Top Tarbert velocity map reveals, however, a fairly scattered distribution, with small patches of highs and lows within the main fault blocks.

Table 5-1
Table 5.1: Statistics of Top Tarbert-Top Ness tomography velocities (in m/s) at the well locations.

For the reservoir itself, represented by the Top Ness velocity map, the distribution is far more coherent. In the Tarbert formation at wells 2 and 3 in the North fault block, there are consistently lower velocities than at well 1 in the East block. This feature is fairly constant for several independent velocity analysis, with a velocity increase of about 200 m/s across the fault separating the gas-bearing reservoir in the North block from the water-bearing reservoir in the East block. A high-velocity ridge separates the lows at well 2 and 3. Distinct low-velocity zones are also seen to the south and southeast that are not correlated with the depth variations. On the other hand, the high-velocity zones in the southwest may be related to the Tarbert dipping down at the western flank. Table 5.1 shows the results of seven independent velocity analysis obtained in the three well locations, revealing that the differences in mean velocity between the “normal” well 1 and the two overpressured wells 2 and 3 are far beyond the estimated error.

Application of the velocity model for pressure prediction

In order to estimate the pressure map in the Tarbert formation, we follow the procedure outlined above. Table 5.2 shows the values of the basic physical quantities used to compute the theoretical velocities. Gas density and gas bulk modulus are computed by using the van der Waals equation, as described in Carcione and Gangi (2000b). (The values indicated in the table are for hydrostatic pressure.) We assume, according to Biot’s theory, that n1 = 1, i.e., that at zero differential pressure the frame bulk modulus vanishes. The same assumption has been used for the effective stress coefficient related to the frame rigidity modulus. Figure 5.4 shows the velocity map (left panel) and the overpressure map assuming Sw = 0.35 and a gas saturation Sg = 0.65 (center panel). The picture at the right represents the difference in pore pressure by assuming gas-bearing Tarbert (the center picture) and water-bearing Tarbert (Sw= 0.94 and Sg = 0.06). An overpressure of about 15 MPa is predicted for well 2, while slightly higher overpressure (18 MPa) is predicted for well 3. The direct measurements indicate overpressures of about 15 MPa (see Figure 5.3). Figure 5.4 (right panel) shows that the sensitivity of the model to fluid saturation is about 8 MPa, indicating that significant uncertainity can be attributed to assumption of uniform saturation.

Table 5-2
Table 5.2: Material properties of the single constituents.

The velocity obtained by careful analysis of pre-stack 3-D data from the deep and complex Tarbert reservoir in the Tune gas field is sufficiently sensitive to pressure and pore fluid to perform a meaningful analysis. The velocity and pressure distribution complies well with the structural features of the target and the general geological understanding of the pressure compartments in the Tune field. The partial saturation model used for pressure prediction can conveniently be calibrated against well data, provided that a complete set of logging data are available for the zone of interest. The most important part of the prediction process is the determination of the effective stress coefficients and dry-rock moduli versus effective pressure, since these properties characterise the acoustic behaviour of the rock. The inversion method based on the shaley sandstone model must fix some parameters while inverting the others. For instance, assuming the reservoir and fluid properties (mainly, the saturation values), formation pressure can be inverted. Conversely, assuming the pore pressure, the saturations can be obtained. The latter implies that this method may be used in reservoir monitoring where the pressure distribution is known while saturation, i.e., the remaining hydrocarbon reserves, are uncertain. We have neglected velocity dispersion, which is not easy to take into account, since Q factor measurements are rare and difficult to obtain with enough reliability. When using laboratory data for the calibration, the effect of velocity dispersion can be significant (Pham et al., 2002).

6. Conclusions

Abnormal pressure, or pressures above or below hydrostatic pressure, occurs on all continents in a wide range of geological conditions. In young (Tertiary) sequences, compaction disequilibrium is the dominant cause of abnormal pressure. In older (pre-Tertiary) rocks, hydrocarbon generation and tectonics are most often cited as the causes of overpressure. Hydrocarbon accumulations are frequently found in close association with abnormal pressure. Thus, in exploration for hydrocarbons and exploitation of the reserves, knowledge of the pressure distribution is of vital importance for prediction of migration routes and location reserves, for the safety of the drilling and for optimising the recovery rate in production.

In this study, we have quantified the effect on seismic properties caused by the common mechanisms of overpressure generation such as hydrocarbon generation, oil-to-gas cracking and disequilibrium compaction. Fluid pressure due to hydrocarbon generation in source-rock shale significantly reduces the seismic velocities and enhances the anisotropy. Attenuation and attenuation anisotropy are also strongly effected, thus further enhancing the seismic visibility of overpressure in the shale. In deeply buried reservoirs, oil-to-gas cracking may increase the fluid pressure to reach lithostatic pressure, and the seismic velocities are shown to drop significantly when only a small fraction of oil in a closed reservoir is converted to gas. Non-equilibrium compaction generates abnormal pressures that, under certain conditions, can be detected with seismic methods. In this case, the fluid mixture filling the pore-space has a major influence on P-wave velocity and my cause under- or over-pressures depending on its compressibility and thermal expansion coefficient. Rocks saturated with fluids of low compressibility and high thermal expansion are generally overpressured and can be seismically visible.

The feasibility of pressure prediction from seismic data rely upon two main factors:

  • Accuracy of the relationships linking the rock physics parameters to pore pressure and properties of the seismic wavefield.
  • The quality of the seismic data and the analysis of seismically derived parameters such as the wave velocity and amplitude (AVO).

We have attempted to link the rock physics parameters, based on first-principles, directly to pore pressure whereas the common conversion methods are based on empirical porosity-stress models valid only for shales saturated with water. Our method, on the other hand, applies to shaley reservoir rocks as well as pure sandstone and shale. Since drilling hazard more often is associated with overpressure in reservoir rocks containing gas rather than with low-permeability water saturated shale, a more general method that accounts for variation in lithology and pore-fill can be justified. However, no reliable method of pore pressure estimation from seismic data can be applied without model calibration against rock physics parameters, as obtained from well logs or laboratory experiments. In absence of local wells a calibration based on regional data is normally in favour of importing empirical relationships developed for a different basin.

Recent developments in 3-D seismic data acquisition and processing techniques have provided an appropriate platform for accurate estimation of seismic velocity and amplitude for input to the pressure analysis. This is further fertilised by the recent developments in techniques and availability of software for 3-D prestack depth migration, velocity tomography and AVO, which greatly enhance the feasibility of pore pressure estimation from seismic data.

The velocity obtained from the Tune gas field is sufficiently sensitive to pressure and pore fluid to perform a meaningful analysis. A critical step in the prediction procedure is the determination of the effective stress coefficients and dry-rock moduli versus effective pressure, since these properties characterise the acoustic behaviour of the rock. The inversion method based on the shaley sandstone model must fix some parameters while inverting the others. For instance, assuming the reservoir and fluid properties, formation pressure can be inverted. Conversely, assuming the pore pressure, the saturations can be obtained. The latter also illustrates the ambiguity imbedded in the problem that may be resolved by AVO.



The European Union under the project “Detection of overpressure zones by seismic and well data” has supported a substantial part of this work.


About the Author(s)

José M. Carcione was born in Buenos Aires, Argentina, in 1953. He received the degree “Licenciado; in Fisicas’ from Buenos Aires University in 1978, the “Dottore in Fisica” from Milan University in 1984, and Ph.D. in Geophysics from Tel Aviv University in 1987. In 1987 he was awarded the Alexander von Humboldt scholarship for a position at the Geophysical Institute of Hamburg University, where he stayed from 1987 to 1989.

From 1978 to 1980 he worked at the “ Comisión Nacional de Energfa Atómica” at Buenous Aires. From 1981 to 1987 he worked as a research geophysicist at “Yacimientos Petrolfferos Fiscales” the national oil company of Argentina. Presently, he is a senior geophysicist at the “Instituto Nazionale de Oceanografia et di Geofisica Sperimentale (OGS)” (former Osservatorio Geofisico Sperimentale) in Trieste, where he was Head of the Department of Geophysics from 1996 to 2000. He has been editor of GEOPHYSICS since 1999. His current research deals with numerical modeling, the theory of wave propagation in acoustic and electromagnetic media, and their application to geophysical problems.

Hans B. Helle, received his MSc (Oceanography) from University of Bergen (Norway) in 1973. From 1970 to 1977 he worked as a research fellow and assistant professor at Geophysical Institute (UoB), and in 1978-1980 as a research engineer at SINTEF, Trondheim. From 1980 to 1983 he was an assistant professor at Seismological Observatory (UoB). From 1984 to the present he has been a research geophysicist at Norsk Hydro, E & P Research Center, Bergen. He has written a number of papers on oceanography, seismology and seismic. His current interests are numerical modeling of seismic wave propagation in realistic media and its application to rock and fluid characterization from seismic data. One of his recent sidetracks is computational neural network applied to reservoir characterization from well logs. He is a member of SEG, EAGE and Association of Marine Researchers (Norway).


Audet, D. M., 1996, Compaction and overpressuring in Pleistocene sediments on the Louisiana shelf, Gulf of Mexico, Marine and Petroleum Geology, 13, 467-474.

Barker, C., 1990, Calculated volume and pressure changes during the thermal cracking of oil to gas in reservoirs: AAPG Bull., 74, 1254-1261.

Belotti, P. and Giacca, D., 1978, Seismic data can detect overpressure in deep drilling, Oil and Gas Journal, August issue.

Berg, R. R., and Gangi, A. F., 1999, Primary migration by oil-generation microfracturing in low-permeability source rocks: Application to the Austin chalk, Texas: AAPG Bull., 83 (5), 727-756.

Berryman, J. G., 1992, Effective stress for transport properties of inhomogeneous porous rock, J. Geophys. Res.,97, 17409-17424.

Best, M. E., Cant, D. J., Mudford, B. S., and Rees, J. L., 1990, Can velocity inversion help map overpressure zones? Examples from an offshore margin, 60th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts, 763-765.

Bilgeri, D., and Ademeno, E. B., 1982, Predicting abnormally pressured sedimentary rocks, Geophys. Prosp., 30, 608-621.

Biot, M. A., 1962, Mechanics of deformation and acoustic propagation in porous media: J. Appl. Phys., 33, 1482-1498.

Bowers, G. L., 1995, Pore pressure estimation from velocity data, Accounting for overpressure mechanisms besides undercompaction, IADC/SPE Drilling Conf., #27488, 515- 530.

Bradley, J. S., and Powley, D. E., 1994, Pressure compartments in sedimentary basins: A review: AAPG Memoir 61, 3- 26. The American Association of Petroleum Geologists.

Bredehoeft, J. D., and Hanshaw, B. B., 1968, On the maintenance of anomalous fluid pressures: I. Thick sedimentary sequences, Geol. Soc. Am. Bull., 79, 1097-1106.

Brie, A., Pampuri, F., Marsala, A. F., and Meazza, O., 1995, Shear sonic interpretation in gas-bearing sands, SPE Annual Technical Conf.,#30595, 701-710.

Brown, R. J. S. and Korringa, J., 1975, On the dependence of the elastic properties of a porous rock on the compressibility of the pore fluid: Geophysics, 40, 608-616.

Bryant, T. M., 1989, A dual pore pressure detection technique, SPE #18714.

Cadoret, T., Marion, D., and Zinszner, B., 1995, Influence of frequency and fluid distribution on elastic wavevelocities in partially saturated limestones, J. Geophys. Res., 100, 9789- 9803.

Carcione, J. M., 2000a, Amplitude variations with offset of pressure-seal reflections, Geophysics, 66, 283-293.

Carcione, J. M., 2000b, A model for seismic velocity and attenuation in petroleum source rocks, Geophysics, 66, 1080-1092.

Carcione, J. M., 2000c, AVO effects of a hydrocarbon sourcerock layer, Geophysics, 66, 419-427.

Carcione, J. M., 2001, Wave fields in real media: Wave propagation in anisotropic, anelastic and porous media: Handbook of Geophysical Exploration, vol. 31, Pergamon Press Inc.

Carcione, J. M. and Gangi, A., 2000a, Non-equilibrium compaction and abnormal pore-fluid pressures: effects on seismic attributes, Geophys. Prosp., 48, 521-537.

Carcione, J. M. and Gangi, A., 2000b, Gas generation and overpressure: effects on seismic attributes, Geophysics, 65, 1769-1769.

Carcione, J. M., Gurevich, B. and Cavallini, F., 2000, A generalized Biot-Gassmann model for the acoustic properties of clayley sandstones, Geophys. Prosp., 48, 539-557.

Carcione, J. M., Helle, H. B., Pham, N. H., and Toverud, T., 2001, Pore pressure estimation from seismic reflection data, EAGE/SEG Research Workshop on Reservoir Rocks, Pau, France.

Carcione, J. M., Helle, H. B., Pham, N. H., and Toverud, T., 2002, Pore pressure estimation in reservoir rocks from seismic reflection data, Geophysics (submitted)

Carcione, J. M. and Seriani, G. 2001, Wave simulation in frozen sediments: J. Comput. Phys., 170, 1-20.

Carcione, J. M. and Tinivella, U., 2001, The seismic response to overpressure: a modeling methodology based on laboratory, well and seismic data, Geophys. Prosp., 49, 523-539.

Carcione, J. M. and Cavallini, F., 2002, Poisson’s ratio at high pore pressure, Geophys. Prosp., 50, 97-106.

Chaney, P. E., 1950, Abnormal pressure and lost circulation: World Oil, 130, 122-126.

Chiarelli, A. and Duffaud, F., 1980, Pressure origin and distribution in Jurassic of Viking basin (United Kingdom-Norway): AAPG Bulletin, 64, 1245-1266.

Childs, C., Manzocchi, T., Nell, P. A. R., Walsh, J. J., Strand, J. A., Heath, A. E., and Lygren, T. H., 2002, Geological implications of a large pressure difference across a small fault in the Viking Graben, in Hydrocarbon Seal Quantification. (Koestler, A. G. and Hunsdale, R., Eds.), Norwegian Petroleum Society (NPF), Special Publication, 11, 127-139.

Christensen, N. I., and Wang, H. F., 1985, The influence of pore pressure and confining pressure on dynamic elastic properties of Berea sandstone, Geophysics, 50, 207-213.

Dutta, N. C., 1983, Shale compaction and abnormal pore pressures: A model of geopressures in the Gulf of Mexico Basin, 53rd Ann, SEG Meeting.

Dutta, N. C., and Levin, F. K., 1990, Geopressure, Geophysical Reprint Series No. 7, Society of Exploration Geophysicists.

Domenico, S. N., 1977, Elastic properties of unconsolidated porous sand reservoirs, Geophysics, 42, 1339-1368.

Dvorkin, J. and Walls, J. 2000, Detecting overpressure from seismic velocity calibrated to log and core measurements, 32nd Annual Offshore technology conference, Paper OTC 11912.

Eaton, B. A. 1976, Graphical method predicts geopressure worldwide, World Oil, July issue, 100-104.

Eaton, B. A., and Eaton, T. L., 1997, Fracture gradient prediction for the new generation, World Oil, October issue, 93-100.

Foster, J. B., and Whalen, J. E., 1966, Estimation of formation pressures from electrical surveys Offshore Louisiana, Journal of Petroleum Technology, February issue.

Gangi, A. F., and Carlson, R. L., 1996, An asperity-deformation model for effective pressure, Tectonophysics, 256, 241-251.

Gist, G. A., 1994, Fluid effects on velocity and attenuation in sandstones, J. Acoust. Soc. Am., 96, 1158-1173.

Goldberg, I., and Gurevich, B., 1998, A semi-empirical velocity-porosity-clay model for petrophysical interpretation of P- and S- velocities, Geophys. Prosp., 46, 271-285.

Harrold, T. W. D., Swarbick, R. E., and Goulty, N. R., 1999, Pore pressure estimation from mudrock porosities in tertiary basins, Southeast Asia, AAPG Bulletin, 83, 1057- 1067.

Han, D., Nur, A., and Morgan, D., 1986, Effects of porosity and clay content on wave velocities in sandstones, Geophysics, 51, 2093-2107.

Hart, B. S., Flemings, P. B., and Desphande, A., 1995, Porosity and pressure: role of compaction disequilibrium in the development of geopressures in a Gulf Coast Pleistocene basin, Geology, 23.

Hashin, Z., and Shtrikman, S., 1963, A variational approach to the elastic behavior of multiphase materials, J. Mech. Phys. Solids, 11, 127-140.

Helle, H. B., Bhatt, A., and Ursin, B., 2001, Porosity and permeability prediction from wireline logs using artificial neural networks, A North Sea case study, Geophys. Prosp., 49, 431-444.

Helle, H.B. and Bhatt, A., 2002, Fluid saturation from well logs using committee neural networks: Petroleum Geoscience, 8, 109-118.

Holbrook, P., Maggiori, D. A., and Hensley, R., 1995, Real time pore pressure and fracture pressure determination in all sedimentary lithologies, SPE Formation Evaluation, December issue, 215-222.

Iverson, W. P., Martinsen, R. S., and Surdam, R. C., 1994, Pressure seal permeability and two-phase flow: AAPG Memoir 61, 313-319. The American Association of Petroleum Geologists.

Johnson, D. L., Koplik, J., and Dashen, R., 1987, Theory of dynamic permeability and tortuosity in fluid-saturated porous media: J. Fluid Mech., 176, 379-402.

Kan, T. K., and Sicking, C. J., 1994, Pre-drill geophysical methods for geopressure detection and evaluation in abnormal formation pressures, Chapman, R. Ed., Elsevier, 155- 186.

Keyser, W. Johnston, L. Reeses, R., and Rodriguez, G., 1991, Pore pressure prediction from surface seismic, World Oil, 212, September issue, 115-125.

King, M. S., Marsden, J. R., and Dennis, J. W., 2000, Biot dispersion for P- and S-waves velocities in partially and fully saturated sandstones, Geophys. Prosp., 48, 1075-1089.

Kjartansson, E., 1979, Constant Q-wave propagation and attenuation, J. Geophys. Res., 84, 4737-4748.

Knight, R., and Dvorkin, J., 1992, Seismic and electrical properties of sandstones at low saturations, J. Geophys. Res.,97, 17425-17432.

Koren, Z., Kosloff, D. Zackhem, U., and Fagin, S., 1998, Velocity model determination by tomography of depth migrated gathers, in S. Fagin, ed., Model-Based Depth Imaging, Society of Exploration Geophysicist, Course Note Series Nr. 10, 119-130.

Krief, M., Garat, J., Stellingwerff, J., and Ventre, J., 1990, A petrophysical interpretation using the velocities of P and S waves (full waveform sonic), The Log Analyst, 31, 355-369.

Kuster, G. T., and Toksöz, M. N., 1974, Velocity and attenuation of seismic waves in two-phase media: Part I. Theoretical formulations, Geophysics, 39, 587-606.

Law, B. E., Ulmishek, G. F., and Slavin, V. I. (eds.), 1998, Abnormal pressures in hydrocarbon environments: AAPG Memoir 70. The American Association of Petroleum Geologists.

Leclaire, Ph., Cohen-Tenoudji, F., and Aguirre-Puente, J., 1994, Extension of Biot’s theory of wave propagation to frozen porous media, J. Acoust. Soc. Am., 96, 3753-3768.

Lesso, W. G. Jr., and Burgess, T. M., 1986, Pore pressure and porosity from MWD measurements, IADC/SPE Drilling Conf., #14801.

Louis, J. N., and Asad, A. M., 1994, Seismic amplitude versus offset (AVO) character of geopressured transition zones: AAPG Memoir 61, 131-137. The American Association of Petroleum Geologists.

Luo, X., and Vasseur, G., 1996, Geopressuring mechanism of organic matter cracking: numerical modeling, AAPG Bull., 80, 856-874.

Mann, D. M., and Mackenzie, A. S., 1990, Prediction of pore fluid pressures in sedimentary basins: Marine and Petroleum Geology, 7, 55-65.

Mavko, G., Mukerji, T., and Dvorkin, J., 1998, The rock physics handbook: tools for seismic analysis in porous media, Cambridge Univ. Press.

Matusevich, V. M., Myasnikova, G. P., Maximov, A. M., Volkov, A. M., Christiakova, N. F., Kanalin, V.G. and Pupilli, M., 1997, Abnormal formation pressure in the West Siberian Mega-basin, Russia: Petroleum Geoscience, 3, 269-283.

Meissner, F. F., 1978, Petroleum geology of the Bakken Formation, Williston Basin, North Dakota and Montana: Montana Geol. Soc. Symp., 207-227.

Miley, M. P., 1999, Converted modes in subsalt seismic exploration, M.Sc., Thesis, Rice University.

Miley, M. P., and Kessinger, W. P., 1999, Overpressure prediction using converted mode reflections from base of salt, 60th Ann. Internat. Mtg., Soc.Expl. Geophys., Expanded Abstracts, 880-883.

Murphy, W. F., 1982, Effect of partial water saturation on attenuation of Massilon sandstone and Vycor porous glass, J. Acoust. Soc. Am., 71, 1458-1468.

Pennebaker, E., 1968, Seismic data indicate depth and magnitude of abnormal pressures, World Oil, 166, 73-78.

Pham, N. H., Carcione, J. M., Helle, H. B., and Ursin, B., 2002, Wave velocities and attenuation of shaley sandstones as a function of pore pressure and partial saturation. Geophys. Prosp.(in press).

Prasad, M., and Manghnani, M. H., 1997, Effects of pore and differential pressure on compressional wave velocity and quality factor in Berea and Michigan sandstones, Geophysics, 62, 1163-1176.

Rubey, W. W., and Hubbert, M. K., 1959, Role of fluid pressure mechanics of overthrust faulting, II. Overthrust belt in geosynclinal area of Western Wyoming in light of fluid pressure hypothesis, Geol. Soc. Am., 70, 167-205.

Sayers, C. M., Johnson, G. M., and Denyer, G., 2000, Predrill pore pressure prediction using seismic data, IADC/SPE Drilling Conf., #59122.

Schoenberg, M., and Muir, F., 1989, A calculus for finely layered media: Geophysics, 54, 582-590.

Smith, J. E., 1971, The dynamics of shale compaction and evolution of pore-fluid pressure, Math. Geol., 3, 239-263.

Terzaghi, K., 1936, The shearing resistance of saturated soils and the angle between the planes of shear, in: Proceedings of the International Conference on Soil Mechanics and Foundation Engineering, Harvard University Press, Cambridge, 1, 54-56.

Todd, M. P., and Simmons, W. S., 1972, Effect of pore pressure on the velocity of compressional waves in low-porosity rocks, J. Geophys. Res., 70 (20), 3731-3743.

Tinivella, U., Carcione, J. M., and Helle, H. B., 2001, Estimation of pore pressure by AVO inversion, 71th Ann. Internat. Mtg., Soc. Expl. Geophys., Expanded Abstracts (CD-rom).

Traugott, M.,1997, Pore/fracture pressure determinations in deep water, Deepwater Technology supplement, World Oil, 218 August issue.

Vernik, L., 1994, Hydrocarbon-generation-induced microcraking of source rocks, Geophysics, 59, 555-563.

Vernik, L., 1995, Petrophysics of the Kimmeridge shale, North Sea, Stanford Rock Physics Laboratory.

Vernik, L. and Nur A., 1992, Ultrasonic velocity and anisotropy of hydrocarbon source rocks, Geophysics, 57, 727-735.

Vernik, L. and Landis, C., 1996, Elastic anisotropy of source rocks: implications for hydrocarbon generation and primary migration: AAPG Bulletin, 80, 531-544.

Winkler, K. W., 1985, Dispersion analysis of velocity and attenuation in Berea sandstone, J. Geophys. Res., 90, 6793-6800.

Yin, C. S., Batzle, M. L., and Smith, B. J., 1992, Effects of partial liquid/gas saturation on extensional wave attenuation in Berea sandstone, Geophys. Res. Letters, 19, 1399-1402.

Zimmerman, R. W., 1991, Compressibility of sandstones, Elsevier, New York


Join the Conversation

Interested in starting, or contributing to a conversation about an article or issue of the RECORDER? Join our CSEG LinkedIn Group.

Share This Article