Use of Seismic Spectral Decomposition Phase and Re
Use of Seismic Spectral Decomposition Phase and Re
Use of Seismic Spectral Decomposition Phase and Re
sciences
Article
Use of Seismic Spectral Decomposition, Phase, and Relative
Geologic Age as Attributes to Improve Quantitative Porosity
Prediction in the Daqing Field, China
David Mora Calderon 1, *, John P. Castagna 1 , Ramses Meza 2 , Shumin Chen 3 and Renqi Jiang 4
1 Department of Earth and Atmospheric Sciences, University of Houston, 4800 Calhoun Rd.,
Houston, TX 77004, USA; [email protected]
2 BHP Billiton, Houston, TX 77056, USA; [email protected]
3 Daqing Exploration and Development Research Institute, PetroChina, 9 Dongzhimen North Street,
Dongcheng District, Beijing 100007, China; [email protected]
4 New Horizon, Ltd., Beijing, China; [email protected]
* Correspondence: [email protected]; Tel.: +1-832-720-0539
Abstract: The high production potential of the Daqing oilfield in China is recognized for seismically
thin sand bodies that usually are not resolved with conventional seismic data. The present study
assesses the usefulness of applying seismic multi-attribute analysis to bandwidth extended data
in resolving and making inferences about these thin layers. In thin layers, tuning can obscure
relationships between seismic amplitude and rock properties. In such cases, the seismic phase varies
with the layer impedance and may hence aid in reservoir characterization. A seismically derived
relative geologic age may also be a useful attribute in predicting rock properties because it helps
Citation: Mora Calderon, D.;
define the stratigraphic position of a layer. When utilized in multi-attribute analysis in the Daqing
Castagna, J.P.; Meza, R.; Chen, S.;
Jiang, R. Use of Seismic Spectral
field, spectral decomposition amplitude, phase, and a relative geological age attribute to improved
Decomposition, Phase, and Relative prediction of well log effective porosity from seismic data and are preferentially selected by stepwise
Geologic Age as Attributes to regression. The study follows standard methodology by implementing seismic multi-attribute
Improve Quantitative Porosity analysis and discusses the improvement of applying it to bandwidth extended data. This will include
Prediction in the Daqing Field, China. a combination of attributes such as relative geologic age, phase, amplitude, and the magnitude
Appl. Sci. 2021, 11, 8034. https:// components of spectrally decomposed data.
doi.org/10.3390/app11178034
Keywords: reservoir characterization; seismic resolution; seismic attributes; thin reservoirs; porosity prediction
Academic Editors: Domenico Patella
and Paolo Mauriello
prediction accuracy and precision when using the spectrally broadened volume as input
data. In that paper and this study, the spectrally broadened volume is obtained by sparse-
layer inversion using basis pursuit decomposition [5]. The bandwidth extension algorithm
is fully described and tested in [5] and theoretical aspects of the validity of the bandwidth
extension are discussed and demonstrated in [6].
In this paper, our objective is to improve the accuracy and resolution of porosity
prediction in the Daqing field, China, utilizing the spectral decomposition of bandwidth
extended seismic data. The use of the seismic phase in rock properties prediction is
investigated and the study introduces a novel seismic attribute termed relative geologic age.
A comparison between predictions with spectral decomposition amplitude, magnitude,
phase, and relative geologic age, in addition to the results obtained with conventional
attributes—such as those utilized by [1] in the multi-attribute analysis—will be discussed.
The Daqing field and, more specifically, the Qingshankou-Putahoa/Shaertu petroleum
system have been broadly explored [7]. Hydrocarbons have been produced from thick
reservoirs and potential has been identified on thinner layers [7] but it has been difficult
to detect and correlate reservoirs that are significantly thinner than conventional seismic
resolution, which is approximately 14 wavelength at the dominant frequency of the original
seismic data (about 33 m for this dataset in the Daqing field). In the study area, well logs
reveal hydrocarbon-bearing intervals with thicknesses ranging from 1 to 15 m that cannot
be seismically resolved.
The spectrally broadened seismic data has been previously shown to yield more
accurate porosity predictions in this field (Mora et al., 2020). Considering the reservoirs
are greatly below the limit of seismic resolution for the original seismic data, bandwidth
extended data are used as inputs to extract seismic attributes and evaluate the predic-
tion improvement when using seismic attributes that are not usually employed for such
purposes [1]. The present study specifically utilizes spectral decomposition amplitude,
magnitude, and phase attributes and employs a relative geologic age volume derived from
picked horizons. Then, constrained least squares spectral analysis [8] is implemented to
perform spectral decomposition. The inputs to the multi-attribute analysis are the am-
plitude, magnitude, and phase attributes extracted at each frequency. The attributes are
divided into four groups: conventional attributes (Group 1), conventional and relative
geologic age (Group 2), conventional and spectral decomposition attributes (Group 3),
and all the above (Group 4). The multi-attribute analysis is performed separately for each
group and the generated porosity volumes are compared. The resulting accuracy, precision,
statistical significance, and resolution in the determination of layer porosity are quantified
and contrasted.
2. Dataset
The dataset consists of a subset of a 3D seismic survey (for details see [2]), covering an
area of about 4 km2 (Figure 1a), consisting of spectrally broadened data (Figure 1b) and
24 wells all including density-derived porosity logs.
Appl. Sci. 2021, 11, 8034 3 of 21
Figure 1. (a) Location of the seismic survey. Black dots represent the relative location of 24 wells.
(b) Inline of spectrally broadened volume. The yellow horizon represents the top of the reservoir interval.
Lian and Castagna [6] shows that spectral broadening accomplished by using sparse-
layer inversion can be valid if the data quality permits, and the subsurface can be repre-
sented as a sparse number of layers. The study concluded that in a worst-case scenario
with a highly stratified random earth reflectivity, no practical spectral broadening could
be achieved outside the band of the original data. However, the researchers reasoned that
this extrapolation could be useful in boosting the high-frequency signal within the original
seismic bandwidth more than the high-frequency noise in spectral bluing, hence improving
the seismic resolution within the original band of the data.
group and the results are compared to determine the accuracy, resolution, and statistical
significance of the predictions. In this study, the spectrally broadened seismic volume
is considered a seismic attribute, without any requirement or conclusion regarding the
accuracy or completeness of the inferred high-frequency content (Figure 2). The spectrally
broadened data has a high-frequency cutoff of about 150 Hz that is more than twice the
original seismic bandwidth. The comparison of the original and spectrally broadened
data is presented in [2]. The rock properties prediction is an empirical statistical approach
that can leverage geological correlations that are not comprehended by a rock properties
inversion based entirely on inverting the physics of wave propagation. For this reason, it is
not required that the spectral broadening be physically correct such that it will reproduce
the seismic data acquired at those frequencies. However, it is required that the spectral
broadening be useful for the research purpose, determined by the well-to-seismic tie and
the validation results of the blind out-of-sample tests [2].
3.1. Well-Tie
Using a statistical wavelet derived from the seismic data spectrum, the correlation
coefficient of synthetic to spectrally broadened seismic data varied between 0.70 and 0.85 in
the 20 wells. No stretch or squeeze was applied in any case.
Figure 3 shows an example of a synthetic tie for one well and the wavelet used.
Figure 3. Well-tie to spectrally broadened data; the correlation coefficient is 0.75. Red trace: composite
seismic trace. Blue trace: synthetic trace. The right side shows the statistical wavelet used. Polarity:
an increase in amplitude corresponds to an increase in acoustic impedance.
Appl. Sci. 2021, 11, 8034 6 of 21
Figure 4. Inline example of relative geologic age volume. Each color represents a different seismically
inferred relative geologic age.
Appl. Sci. 2021, 11, 8034 7 of 21
Figure 5. (a) Amplitude at 100 Hz (polarity: an increase in amplitude corresponds to an increase in acoustic impedance),
Figure 5. (a) Amplitude at 100 Hz (polarity: an increase in amplitude corresponds to an increase in
(b) magnitude at 100 Hz, and (c) phase at 100 Hz.
acoustic impedance), (b) magnitude at 100 Hz, and (c) phase at 100 Hz.
Group 1 Group 2
Cosine Instantaneous Phase Cosine Instantaneous Phase
Integrated Absolute Amplitude Relative Geologic Age
Filter 5 Hz/10 Hz–15 Hz/20Hz Filter 5 Hz/10 Hz–15 Hz/20Hz
X-Coordinate Filter 45 Hz/50 Hz–55 Hz/60Hz
Filter 45 Hz/50 Hz–55 Hz/60Hz Average Frequency
Average Frequency X-Coordinate
Amplitude Weighted Phase Filter 55 Hz/60 Hz–65 Hz/70Hz
Amplitude Weighted Cosine Phase Y-Coordinate
Quadrature Trace Instantaneous Frequency
Raw Seismic Filter 35 Hz/40 Hz–45 Hz/50Hz
Group 3 Group 4
Cosine Instantaneous Phase Cosine Instantaneous Phase
Magnitude at 50 Hz Relative Geologic Age
Phase at 30 Hz Cosine of the Phase at 110 Hz
Cosine of the Phase at 110 Hz Phase at 30 Hz
Magnitude at 30 Hz Filter 45 Hz/50 Hz–55 Hz/60Hz
Amplitude at 70 Hz Semblance
Filter 45 Hz/50 Hz–55 Hz/60Hz Average Frequency
Cosine of the Phase at 70 Hz Cosine of the Phase at 70 Hz
Cosine of the Phase at 30 Hz Phase at 70 Hz
Magnitude at 50 Hz Cosine of the Phase at 30 Hz
Cosine of the instantaneous phase was selected by stepwise regression as the best
attribute in all cases, followed by relative geologic age when available. For Groups 3 and
4, phase attributes represented 50% and 60% of the top ten, respectively. The elimination
of the need for x and y coordinates for both Groups 3 and 4 suggests that the inclusion
of spectral decomposition volumes compensates for any lateral drift in the prediction
calibration and may account for the greater lateral variability in the predicted porosity
resulting from these groups (see below).
Figure 6. (a) Left: inline of synthetic seismic data that represents constant porosity of 20 PU; crossline
variations represent thickness changes. Right: cosine of the phase of the inline on the left. (b) Left:
crossline of synthetic seismic data that represents a constant thickness of 51 m; inline variations
represent porosity changes. Right: cosine of the instantaneous phase of the crossline on the left.
Polarity: an increase in amplitude corresponds to an increase in acoustic impedance.
Figure 6 indicates that amplitude and cosine of the phase vary regularly with changing
porosity or thickness.
To better illustrate such changes, Figure 7 shows a surface plot in which the x-axis
represents porosity, the x-axis represents thickness, and the color is the cosine of the
instantaneous phase from the previously described synthetic model for a constant time
sample at 700 ms from the synthetic cube. The pattern of almost linear and parallel lines
of the constant phase indicates systematic cosine of the instantaneous phase sensitivity
to porosity and thickness, thereby explaining why the seismic phase can be useful in
rock properties prediction. In addition, phase has the advantage of being less affected by
overburden attenuation and transmission losses than amplitude, making it possibly more
robust in relating the seismic response to the layer impedance variations than amplitude.
This a potentially important advantage of phase over amplitude in quantitative analyses of
seismic data for rock properties prediction.
Appl. Sci. 2021, 11, 8034 10 of 21
Figure 7. Cosine of the instantaneous phase variation at 700 ms for porosity versus thickness wedge
models. (a) Crossplot view and (b) 3D view.
Figure 7 shows a plot that represents cosine instantaneous phase variations with porosity.
Figure 8 shows the resulting validation error versus the number attributes used in the
regression.
Figure 8. Average training error versus the number of attributes testing convolutional operator
lengths from 1 to 51 samples every 10 points. The black arrow indicates the point of the lowest
validation error achieved using a 41-point operator. The example shown corresponds to the Group 4
attribute selection process. As the sample rate is 1 ms, the 41-point operator represents 40 ms.
2. Out of the 24 available wells, 21 were used for training the algorithm and 3 were left
as out-of-sample tests for blind validation.
3. The transformation for each group was obtained using a 41-point operator length
(representing 40 ms) and the attributes in Table 1.
Figure 9 shows the measured and predicted porosity logs, filtered to 100 Hz, for one
validation well. Group 4 attributes generally gave the best porosity prediction (quanti-
fied below).
4. Cross-validation among the training wells using the procedure described by Hamp-
son [1] was used to restrict the number of attributes employed to 15.
5. Porosity volumes were built. Figure 10 shows an inline of the predicted porosity
for each volume. The inserted curve is a porosity log filtered to 100 Hz from one
validation well and the color background of the curve represents the filtered porosity
log colored to the same scale as the volume. Notice a tradeoff between matching the
well logs and lateral stability of the predictions.
Average porosity maps over a 20-ms window below the seismic pick for the top of
the most important reservoir layer were extracted to compare to interpolated average
porosities from logs over the same window. Figure 11 shows average porosity maps for
predictions made using the four attribute groups.
Appl. Sci. 2021, 11, 8034 12 of 21
Figure 9. Comparison of resulting predicted log (red) and actual porosity log filtered to 100 Hz
(black) for a well not used in the training process using the different attribute groups. Error is the
mean absolute error over the prediction interval. The blue rectangle highlights the main reservoir
unit. A small time–depth error at the reservoir level is more obvious for Groups 3 and 4.
Figure 10. Inline of the predicted porosity volume obtained using each of the attribute groups. The
inserted curve is a porosity log filtered to 100 Hz from one validation well and the color background
of the curve is the log colored to the same scale as the volume. The dashed black line indicates the
top of the main reservoir. Abbreviation: SD, spectral decomposition.
Appl. Sci. 2021, 11, 8034 13 of 21
Figure 11. (a) Interpolation of the average porosity over a 20-ms window below the main reservoir
top shown in Figure 10. (b) Group 1 map of average porosity over a 20-ms window below the main
reservoir seismic pick in Figure 10. (c) Group 2 map of average porosity over a 20-ms window below
the main reservoir seismic pick shown in Figure 10. (d) Group 3 map of average porosity over a
20-ms window below the main reservoir seismic pick shown in Figure 10. (e) Group 4 map of average
porosity over a 20-ms window below the main reservoir seismic pick shown in Figure 10.
Maps were obtained by horizon-slicing porosity volumes. Black dots are well locations.
Although all maps in Figure 11b–e follow a similar general trend, when compared to the
interpolated porosity map from the logs (Figure 11a), Group 1 fails to accurately predict
average porosity at well locations. The other three maps better match values from logs
but away from well control, they behave differently. The ones from Groups 3 and 4 show
more lateral heterogeneity and apparently sharper detail that could represent geological
variation or could be due to seismic noise or numerical instability. Choosing between
these explanations is an interpretive judgment, keeping in mind that Group 4 has the best
accuracy at validation wells, as will be shown in the next section.
To better observe possibly geological behavior, Figure 12 shows a zoomed-in version
of the map of average porosity obtained by using attributes in Group 4 indicating possible
channels in a lacustrine fan-delta environment; being more channelized to the north and
northeast and more delta front and fan-like in the center and to the south.
Figure 12. Map of average porosity over a 20-ms window below the main reservoir top. The map
was obtained by slicing the porosity volume from Group 4. Black dots are well locations. The black
circle highlights an example of what could be channel-like features in a fan-delta environment.
Figure 13. Cross-plots of the actual versus predicted porosity for predictions made using each of the
attribute groups in three blind validation wells. Abbreviations: CC, correlation coefficient and MAE,
mean absolute error. Data from the entire length of available logs were used in this cross-plot.
Appl. Sci. 2021, 11, 8034 15 of 21
Table 2 gives the resulting values of the correlation coefficient between the mea-
surements and predictions, mean absolute error, standard error, F (Fisher’s F-test), and
significance F as measures of the regression statistical significance for each case.
Table 2. R, mean absolute error (PU), standard error, F, and Significance F for actual versus predicted
porosity logs for only the wells not used in the training process.
Figure 13 shows cross-plots of the predicted versus actual porosity for all attribute
groups’ predictions.
The correlation coefficient (R) is a measure of the degree to which the observation and
prediction vary linearly rather than being a direct measure of accuracy. Mean absolute
error (MAE: the arithmetic average of absolute errors) measures the aver-age magnitude of
discrepancy between pairs of predicted versus observed values. MAE is expected to be as
small as possible. However, as the sign of the error is lost in the magnitude, it does not
reveal bias. It indicates precision better than accuracy.
To evaluate accuracy and bias, the standard error of the predicted mean is used. F is
the ratio of the predicted variance to the error variance adjusted for the degrees of freedom
(one less than the number of data points minus the number of attributes). When the ratio
is less than 1, it can be concluded that the prediction is not useful. Significance F gives
the probability that the correlation of the prediction and observation is a result of random
chance [16]. A value below 0.05 is usually interpreted to mean the prediction is statisti-
cally significant. From Table 2, we can conclude that all the predictions are statistically
significant, with bias and precision both being about two porosity units, approaching the
accuracy of the logs. Compared to the measured log porosity, Group 4 has the best correla-
tion, accuracy, and precision. Adding spectral decomposition to conventional attributes
(Group 3) results in better incremental improvement compared to adding only relative
geologic age (Group 2).
From Table 2, the correlation coefficient increases from 0.64 (Group 1) to 0.79 (Group 4).
F is much greater than 1 and Significance F is virtually zero for all cases, which indicates
that none of the models provides a spurious correlation (a correlation by random chance) of
the predicted versus actual porosity. All the statistics improve systematically from Group 1
to Group 4.
Table 2 was calculated with point-by-point predicted-versus-actual porosity pairs
for about 200 ms of data per well, spanning a larger interval than the main reservoir.
For reservoir volumetrics, the porosity variation with depth (or seismic record time) is less
significant than the average porosity over the reservoir interval. Figure 14 shows actual vs.
predicted average porosity plots over reservoir interval.
Appl. Sci. 2021, 11, 8034 16 of 21
Figure 14. Cross-plots of the actual versus predicted average porosity extracted over a 20-ms window
below the main reservoir. Abbreviations: CC, correlation coefficient and MAE, mean absolute error.
For predicting the average porosity of the interval (Table 3), there is a more significant
systematic improvement in statistics from Group 1 to Group 4, with correlation and
statistical significance increasing and prediction error decreasing. The poor correlations for
Group 1 and Group 2 have a significance of F greater than 0.05, which usually indicates
that the correlation is not statistically significant. Such comparisons suggest that attribute
Groups 3 and 4 are particularly important for accurately predicting average porosity at the
reservoir level.
Table 3. Correlation coefficient and mean absolute error (PU), standard error, F, and Significance F
for the actual versus predicted average porosity over a 20-ms window below the main reservoir.
Average
Group 1 Group 2 Group 3 Group 4
Porosity
Correlation
0.4 0.5 0.65 0.73
Coefficient
Mean Absolute
2.11 1.81 1.7 1.33
Error
Standard Error 1.86 1.83 1.59 1.39
F 0.83 1.55 8.29 17.11
Significance F 3.70 × 10−01 2.20 × 10−01 9.20 × 10−03 5.11 × 10−04
3.7. Resolution
The limit of vertical resolution (λ/4) is 10.04 m for the spectrally broadened seismic
data used as input for this study.
Figure 15 shows thin hydrocarbon-bearing intervals evidenced on gamma-ray and
resistivity logs, as well as associated porosities for a blind validation well that was not
included in the training process.
Appl. Sci. 2021, 11, 8034 17 of 21
Figure 15. Thin hydrocarbon-bearing intervals seen on logs and predictions made using the four
attribute groups. The inserted curve is a porosity log filtered to 100 Hz. The well shown was not used
for training. Black ellipse highlights a reservoir of interest, most detectable using Group 4 prediction.
All highlighted intervals in Figure 15 are at or above the limit of vertical resolution for
the input bandwidth-extended seismic data. The four groups can detect such layers but
only Groups 3 and 4 can accurately predict log porosities in them above 30 PU. The lowest
reservoir just above a record time of 900 ms (black ellipse in Figure 15) is best detected
by Group 4.
The Group 4 porosity volume is also useful to identify geologic features not seen on
the seismic sections alone.
Figure 16 shows examples of channel-like bodies evident on the Group 4 predicted
porosity volume that are not obvious on the original seismic data.
Figure 16. Right: sections from the Group 4 porosity volume showing channel-like features. Left:
same sections from the seismic volume used as input for multi-attribute analysis. Polarity: an increase
in amplitude corresponds to an increase in acoustic impedance.
Appl. Sci. 2021, 11, 8034 18 of 21
High porosity features on the right of Figure 16 resemble incised and stacked-high
porosity channels. This response is expected in the area given that the reservoirs are fluvial
and are lacustrine fan delta deposits. These features are imperceptible or are not clearly
defined on seismic sections on the left.
Another example is shown in Figure 17. In this case, one well passes through the
channel-like feature. Porosity and gamma-ray logs are shown, as well as the corresponding
seismic section.
Figure 17. From left to right: porosity and gamma-ray log. Group 4 predicted porosity volume, with
a red line indicating the well path, corresponding to the spectrally broadened seismic section. The
blue rectangle highlights the channel-like feature. Polarity: an increase in amplitude corresponds to
an increase in acoustic impedance.
The computed porosity section in Figure 17 shows a high porosity possible channel
and this interpretation is supported by well log porosity and gamma-ray responses, which
indicate a very porous clean sand. The blocky shape of the sand layer suggests that the
interval may correspond to a fluvial or deltaic distributary channel. The feature is about 10
m thick at the limit of seismic resolution and it is not obvious on the spectrally broadened
seismic section. Thinner sand intervals are not resolved by the predicted porosity volume.
4. Discussion
Spectral broadening outside the original band of the seismic signal is controversial [6].
In the present study, the spectrally broadened seismic data is considered as an attribute: an
arbitrary mathematical transformation of a seismic trace. Although Liang and Castagna [6]
discussed under what circumstances such broadening can be viewed as physically accurate,
Mora et al. [2] showed—for the same area in the Daqing field—that multi-attribute analysis
of spectrally broadened data resulted in better predictions of well log porosity by all
measures compared with a similar analysis with the original conventional seismic data.
Based on the results, the spectrally broadened data is a significant input for multi-attribute
analysis as indicated by cross-validation and blind well tests. The significance of the
attribute was demonstrated by ties to well data and the performance in predicting porosity
in both the training and validation [2], rather than by theoretical argument.
The attributes were derived from spectrally broadened seismic data and classified
into four groups: conventional attributes (Group 1), conventional attributes and relative
geologic age (Group 2), conventional attributes and components of spectral decomposition
(Group 3), and all attributes combined (Group 4). Then, stepwise regression was performed
for the attributes in each group to create a multi-attribute list that included the ten attributes
within each group that yielded the best prediction error.
From the attributes included in the stepwise regressions, the cosine of the instanta-
neous phase was selected as the most significant attribute in all the cases, followed by
relative geologic age when available. For Groups 3 and 4, more than half of the selections
by stepwise regression were phase attributes. The high ranking of relative geologic age
on the attribute lists is related to its ability to pinpoint where in a stratigraphic sequence
Appl. Sci. 2021, 11, 8034 19 of 21
a particular layer occurs, thereby taking advantage of any geologic succession in rock
properties that may be learned in training. Conversely, the selection of the cosine of the
instantaneous phase first in all the cases was an unexpected outcome, as it was the stepwise
regression process that selected many phase attributes among the top 10 attributes.
The relationship between phase and porosity is explained with a synthetic wedge
model (Figures 6 and 7), showing that phase is sensitive to porosity and thickness in
seismically thin layers. The isophase trends in the cross-plots of phase versus porosity and
thickness suggest that the seismic phase response is related to the phase of the series of
reflection coefficients for a given thickness. In addition, this is directly related to reservoir
porosity through layer impedance. In synthetic data, the cosine of the instantaneous phase
varies systematically with the porosity at a given time thickness. This relationship is
particularly relevant in the present study that focuses on porous, thin, and discontinuous
reservoirs. The cosine of the instantaneous phase removes all amplitude information that
could be ambiguous in predicting porosity in thin layers and with overburden effects. As a
result, the seismic phase potentially is more readily correlated to porosity than amplitude
alone in thin layers. For this reason, it is understandable from a physical perspective that
the attribute is selected first by stepwise regression and that the other phase attributes are
also selected. The agreement between the predicted and actual porosity increases from
Group 1 to Group 4. This is particularly true at the main reservoir level (blue rectangle
in Figure 9 and dashed line in Figure 10), in which the highly porous interval is better
predicted with the 10 attributes selected from all available attributes (Group 4). From
Figure 10, the porosity seismic sections for all groups match the validation porosity log to
some extent. The Groups 1 and 2 sections show continuous layers of the same porosity,
while predicted porosity sections from Groups 3 and 4 show layers with great lateral
variability. This last behavior is typical of a fan delta depositional environment. When
comparing average porosity at the main reservoir level (Figure 11), Groups 2 to 4 seem
to follow the same trend as the interpolated well logs. However, the Group 4 map shows
more apparent detail and features that are possibly small channels (Figure 12). The areal
extent of the predicted porous reservoir appears greatest for Group 4.
When compared to well logs, the porosity predicted with attributes in Group 4 shows
more significant detail and is quantitatively closer to the measured porosity. The resulting
metrics from the blind-well validation data not included in the training process (Table
2) show that Group 4 predictions are more accurate, more precise, and most statistically
significant, with the highest correlation coefficient and F-value, as well as the lowest
absolute and standard errors and lowest Significance F. Adding only relative geologic
age (Group 2) increases the correlation coefficient and decreases error relative to Group 1.
Adding only spectral decomposition components to the conventional attributes (Group 3)
increases the accuracy and statistical significance of the prediction.
The differences in the metrics between the groups are more drastic when predict-
ing average porosity at the reservoir level (Table 3), with a smaller error and both an
increased correlation coefficient and statistical significance for the predictions made with
Groups 3 and 4.
Figure 15 shows a comparison of the measured and predicted porosity at one well
location used in the validation, in which three intervals of thicknesses that are barely above
the λ/4 of 10.04 m are highlighted. The Group 3 and 4 predictions can better resolve these
thin layers and the predicted values are closer to the actual values as compared to the
predictions made with Groups 1 and 2.
The Group 4 predicted porosity volume can also be useful in identifying and delineat-
ing geological features of interest. Figure 16 shows three examples of channel-like features
with high predicted porosity. It also shows a comparison to the corresponding spectrally
broadened seismic section in which the same features are not obvious. Figure 17 shows a
similar predicted high porosity penetrated by a well in which the logs validate the porosity
prediction and channel interpretation.
Appl. Sci. 2021, 11, 8034 20 of 21
5. Conclusions
Resolving thin layers of interest and predicting the rock properties variation in layers
away from the well control are fundamental requirements in seismic reservoir characteriza-
tion. In the present study, both problems were addressed by applying a combination of
bandwidth extension, multi-attribute analysis, and non-conventional attributes to produce
predicted porosity volumes for a portion of the Daqing field in China. The attributes
were derived from spectrally broadened seismic data and multi-attribute analysis was
implemented with four groups of attributes separately. Group 1 is comprised of conven-
tional attributes. Group 2 is comprised of conventional attributes and relative geologic age.
Group 3 is comprised of conventional attributes and spectral decomposition components
including amplitude, magnitude, phase, and the cosine of the phase at 10 Hz increments.
Group 4 is comprised of all the attributes listed prior.
The result presents the comparison of the resulting resolution, accuracy, and statistical
significance of multi-attribute rock properties prediction for each group. It was observed
that the geometries of the predicted porosities was geologically reasonable for the fan delta
environments responsible for depositing the thin reservoir sandstones in the Daqing field.
Group 1 gave the least accurate predictions. In contrast, the predictions with both spectral
decomposition volumes and the relative geologic age gave the best performance by all the
metrics. Phase attributes improved porosity prediction despite their independence from
the seismic amplitude.
The results show that combining spectrally broadened seismic data and multi-attribute
analysis with non-standard seismic attributes can be a powerful approach for resolving,
delineating, and characterizing reservoirs at or below the limit of conventional seismic res-
olution.
Author Contributions: Data preparation and review, R.J. and S.C.; conceptualization and writing,
D.M.C., J.P.C. and R.M.; investigation, D.M.C. and J.P.C.; writing—review and editing, D.M.C., J.P.C.,
R.M., R.J. and S.C. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.
Informed Consent Statement: Not applicable.
Data Availability Statement: Data used in this study is confidential.
Acknowledgments: We thank Brian Russell, Dan Hampson, and CGG for their donation of the
EMERGE software package to the University of Houston. Thanks are also given to dissertation
committee members Evgeni Chesnokov and Lorenzo Colli for their guidance, advice, and support.
Conflicts of Interest: The authors declare no conflict of interest.
References
1. Hampson, D.; Quirein, J.; Schuelke, J. Use of multi-attribute transforms to predict log properties from seismic data. Geophysics
2001, 66, 220–236. [CrossRef]
2. Mora, D.; Castagna, J.; Meza, R.; Chen, S.; Jiang, R. Case study: Seismic resolution and reservoir characterization of thin sands
using multi-attribute analysis and bandwidth extension in the Daqing field, China. Interpretation 2020, 8, T89–T102. [CrossRef]
3. Ray, A.; Biswal, S. An efficient method of effective porosity prediction using an unconventional attribute through multi-attribute
regression and probabilistic neural network: A case study in a deep-water gas field, East Coast of India. SEG Tech. Program.
Expand. Abstr. 2010, 2021, 1413–1417.
4. Abdulaziz, A. The effective seismic attributes in porosity prediction for different rock types: Some implications from four case
studies. Egypt. J. Pet. 2020, 29, 95–104. [CrossRef]
5. Zhang, R.; Castagna, J. Seismic sparse-layer reflectivity inversion using basis pursuit decomposition. Geophysics 2011, 76, 147–158.
[CrossRef]
6. Liang, C.; Castagna, J.; Zavala, R. Tutorial: Spectral bandwidth extension—Invention versus harmonic extrapolation. Geophysics
2017, 82, w1–w16. [CrossRef]
7. Ryder, R.; Quiang, J.; McCabe, P.; Nuccio, V.; Persits, F. Qingshankou-Putaohua/Shaertu and Jurassic Coal—Denglouku/Nongan Total
Petroleum Systems in the Songliao Basin, China; U.S. Geological Survey: Reston, VA, USA, 2003. [CrossRef]
8. Puryear, C.; Castagna, J.; Portniaguine, O.; Cobos, C. Constrained least-squares spectral analysis: Application to seismic data.
Geophysics 2012, 77, V143–V167. [CrossRef]
Appl. Sci. 2021, 11, 8034 21 of 21
9. Chen, S.; Donoho, D. Basis pursuit. In Proceedings of the 1994 28th Asilomar Conference on Signals, Systems and Computers,
Pacific Grove, CA, USA, 2 November 1994; pp. 41–44.
10. Taner, M.; Schuelke, J.; O’Doherty, R.; Baysal, E. Seismic attributes revisited. In SEG Technical Program Expanded Abstracts; Society
of Exploration Geophysicists: Tulsa, OK, USA, 1994; pp. 1104–1106.
11. Barnes, A. The complex seismic trace made simple. Lead. Edge 1998, 17, 473–478. [CrossRef]
12. Chopra, S.; Marfurt, K. Seismic attributes—A historical perspective. Geophysics 2005, 70, 3SO–28SO. [CrossRef]
13. Barnes, A. A tutorial on complex seismic trace analysis. Geophysics 2007, 72, W33–W43. [CrossRef]
14. Stark, T. Relative geologic time (age) volumes—Relating every seismic sample to a geologically reasonable horizon. Lead. Edge
2004, 23, 817–944. [CrossRef]
15. Russell, B.; Hampson, D. Multi-attribute seismic analysis. Lead. Edge 1997, 16, 1439–1444. [CrossRef]
16. Mathews, S. Interpreting Regression Output—Without all the Statistics Theory. 2018. Available online: https://www.
graduatetutor.com/statistics-tutor/interpreting-regression-output (accessed on 5 May 2021).