# Seismic rock physics

Rock physics studies the relationship between physical and elastic properties of rocks and is the basis of quantitative seismic interpretation. It has arguably given exploration geophysicists a solid quantitative basis to their interpretation of seismic data. It can be tackled at different scales of investigation by people with various backgrounds — from geophysicists to civil engineers, from mathematicians to petrophysicists.

In this tutorial, and the Python code in the accompanying Jupyter notebook (https://github.com/seg/tutorials), I offer a practical introduction to rock physics, showing how to code simple models and apply them to analyze well-log data and predict their seismic responses. Here is the workflow:

1. data exploration — inspecting what we have;
2. rock-physics modeling — two simple models;
3. rock-physics templates — modeling more parameters; and
4. seismic response modeling — predicting seismic.

The well data used in this tutorial is one of the data sets available on the Quantitative Seismic Interpretation book website.

## Data exploration

We get our bearings by plotting the logs together. Figure 1a shows $V_{sh}$ (shale volume) log in the first track, $I_{P}$ (P-wave impedance) in the second track, and $V_{P}/V_{S}$ (the ratio of P-wave to S-wave velocities) in the third. I highlight sands in green and shales in red. Next to these three tracks are crossplots of $V_{P}$ against effective porosity $\phi _{e}$ (Figure 1b), and of $V_{P}/V_{S}$ against $I_{P}$ (Figure 1c). Figure 1. (a) Shale volume ($V_{sh}$ ), P-wave impedance ($I_{P}$ in 103 kg $s^{-1}$ $m^{-2}$ ), and $V_{P}/V_{S}$ for the interval 2100 m to 2250 m. (b) P-wave velocity ($V_{P}$ in m/s) versus effective porosity ($\phi _{e}$ ). (c) $V_{P}/V_{S}$ versus $I_{P}$ .

In the $V_{P}-\phi _{e}$ crossplot, I have decided to avoid plotting shales. Since effective porosity is zero for pure shales, this plot will not help in the analysis of nonreservoir rocks.

## Rock-physics modeling

Rock-physics models output dry-rock frame moduli: bulk modulus $K_{d}$ and shear modulus $\mu _{d}$ . These moduli represent the resistance of the rock’s dry frame to normal and tangential stress, respectively. They are used together with Gassmann’s equation to calculate the saturated rock moduli ($K$ and $\mu$ ), which are, together with density, the building blocks for P-wave and S-wave velocities, $V_{P}$ and $V_{S}$ :

$V_{P}={\sqrt {\frac {K+4\mu /3}{\rho }}}\ \ \ V_{S}={\sqrt {\frac {\mu }{\rho }}}.$ The product of velocity and density is called acoustic impedance. Impedance contrasts are ultimately what generate seismic reflections; this is why rock physics, by linking rock and fluid properties to velocities, allows us to understand and explain seismic reflections.

In many cases, we know the velocity and density of a nonreservoir sandstone, and we would like to know how these quantities will change if the porosity were only higher, the net:gross (N:G) more favorable, or the fluid more valuable. We need a model!

In the notebook, I introduce two simple but useful models, the Soft Sand model and the Stiff Sand model. The notebook provides the theory and the code to help you understand these models. Here, I only show what they produce.

We would like to know how well our data fit the models we make. The easiest way is to look at them both in one plot. Figure 2a shows the same interval as before — the reservoir sand from 2100 m to 2250 m — with the Soft Sand model and the Stiff Sand model superimposed. At first sight, both models may offer a good starting point to investigate these sands.

The sands are not hydrocarbon-bearing, so I can dismiss fluid effects and start leveraging the rock-physics model by displaying how mineralogy affects the velocities. To do that, I will simplify the problem assuming variations in N:G. When N:G is 1, then the mineral moduli will be equal to those of quartz; when it is lower than 1, I will use composite mineral moduli given by the Voigt–Reuss–Hill average of quartz and a generic shale. So to compute the bulk modulus given the constituent volume fractions f:

$K_{Voigt}=(f_{qz}\cdot K_{qz})+(f_{sh}\cdot K_{sh})$ $K_{Reuss}={\frac {1}{(f_{qz}/K_{qz})+(f_{sh}/K_{sh})}}$ $K_{Hill}={\frac {K_{Voigt}+K_{Reuss}}{2}}$ Now let’s put it all together and show how the models change for different mineralogies. There’s just one catch: to have the models spread out to cover the entire range of the well data, I had to code them with different parameters. In the notebook, these are the lines that call the two models:

 Kdry, MUdry = softsand(K0, MU0, phi, phic=0.5, Cn=12, P=45) 

Kdry, MUdry = stiffsand(K0, MU0, phi, phic=0.4, Cn=8, P=45) 

The Soft Sand model uses a higher critical porosity and coordination number than the Stiff Sand model, which means that the reservoir sands can be described by either a stiffer “soft” sand or a softer “stiff” sand (the coordination number Cn suggests a closer packing of the individual grains that make up our sand, which in turn means a more compact and harder material).

Figures 2b and 2c show the result of plotting the data against the two models. So far, both models fit the data points. It is up to the interpreter to decide which model is more appropriate. Adding the effect of hydrocarbons and changing other parameters — especially porosity — might help us figure it out. Figure 2. $V_{P}$ versus $\mu _{e}$ for the reservoir interval shown in Figure 1. (a) The Soft Sand model is shown as a black line; the Stiff Sand model is a discontinuous line. (b) The Soft Sand model and (c) the Stiff Sand model with varying N:G.

## Rock-physics templates

The set of models for varying degrees of N:G shown above are a basic example of rock-physics templates (RPTs), as promoted by Avseth et al. (2005). From that starting point, it is fairly easy to build a function to create RPTs that show not only variations in porosity and mineralogy but also fluid content.

In the notebook, I show how to combine the functions to compute $K_{dry}$ and $\mu _{dry}$ for the Soft Sand and Stiff Sand models with the velocities given by Gassmann’s equation, to produce arrays of P-wave impedance $I_{P}$ and $V_{P}\ V_{S}$ that we can study individually — as in Figure 3a — or overlay on top of well data as an aid to interpretation — as in Figures 3b and 3c.

We can go further. For example, one can check the validity of the elastic moduli assigned to the second mineral (the one I have referred to as “shale”) by varying them and displaying different RPTs. In Figure 3c, I have changed the bulk and shear modulus values to get a better fit against the data. In the absence of other information, this is perfectly reasonable — and valuable — knowledge to add to the exploration puzzle. Figure 3. Rock-physics templates for the Soft Sand model. (a) The template by itself, showing how we expect $V_{P}/V_{S}$ and $I_{P}$ to vary with water saturation ($S_{W}$ ) and porosity ($\mu _{e}$ ). (b) The template (gray) plotted with the data (green). (c) The template (blue) has been recalculated by changing the rock matrix moduli, assuming higher proportion of quartz (i.e., cleaner sands).

## Seismic response modeling

Once I have a good grasp of the RPM and the way it reacts to different input parameters, I can start using it to explore seismic data too. It is very easy to use the velocities and densities output by our rock-physics model to build simple two-layer models and analyze seismic and AVO signatures of the modeled top reservoir. One example is shown in Figure 4; the idea is fully developed in the accompanying Jupyter notebook. Figure 4. The expected seismic response of a gas-saturated sand with 15% porosity. From left to right: P-wave impedance $I_{P}$ , $V_{P}/V_{S}$ , the response at 0 degrees and 30 degrees of offset, and the amplitude-versus-angle curve.

## Conclusions

I have shown here how I implement rock-physics models and templates using standard Python libraries. The functions introduced here can be extended fairly easily for additional rock-physics models, some of which are shown in the accompanying notebook (https://github.com/seg/tutorials). The workflow can be used to study real data sets, analyze the elastic behavior of individual reservoirs, and build a collection of seismic and AVO signatures to be associated to particular combinations of cap rocks and reservoirs.