
Citation: | Hongtao Zhu, Yuansheng Du, Keyu Liu, Jiaxin Yan, Yajun Xu, Ping Yang, Xinyu Liu. Sequence Stratigraphic Delineation and Correlation of the Dongying Formation in the Nearshore and Adjacent Sea Areas, Bohai Bay Basin. Journal of Earth Science, 2008, 19(1): 54-64. |
Based on the analysis of well logs, seismic data, core studies, the Dongying (东营) Formation along the nearshore and adjacent sea areas, Bohai (渤海) Bay basin has been subdivided into three third-order sequences from base to top: namely, sequences SQ1, SQ2 and SQ3. The three sequences have different wireline (SP) log responses, showing triple-section characteristics with SQ1 being characterized by primarily flat baseline with intercalation of relatively low spontaneous potential, SQ2 generally exhibiting weak or moderate amplitude spontaneous potential with finger-shaped peaks, and SQ3 having relatively high spontaneous potential with funnel-shaped log curves. On the basis of the triple-section characteristics, the stratigraphic sequences can be correlated consistently throughout the entire study area. A stratigraphic and sedimentary model for sequences SQ1 to SQ3 of the Dongying Formation in the study area has been proposed. The accommodation space change in the two sides of the asymmetrical basin was examined and the asymmetrical basin has the feature of the asymmetrical accommodation space change. The asymmetrical physiography is a vital factor to influence the accommodation changes in additional to the lake level change, tectonism and sediment supply. This may have important implications to similar basins in other parts of China or elsewhere.
The precision of the amplitude-coefficientversus-offset(AVO)inversion result depends on the accuracy of reflection amplitude input. However, for field data there are several other important factors affecting AV0, in addition to the variation of reflection coefficient. These factors include the effect of the free surface, geometric spreading, thin-bed effect, medium absorption, multiples, anisotropy and so on (Strudley, 1991; Sheriff, 1975). Geometric spreading and absorption are major factors. Hence, a series of processes of amplitude compensation needs to be performed before AV0 inversion. Huhral and Kev (1 980)deduced an analytic expression for the spreading factor of a compressional wave in plane layered medium and indicated that the expression is not suitable to multiwave.Samec et al.(1991)studied the influence of viscoelasticity and anisotropy on AVO interpretation and pointed out that one cannot ignore these factors when accurate AVO analysis and inversion are desired.Dong(1999)discussed quantitatively the influence of migration stretch and thin-bed tuning on AVo analysis and interpretation, and derived two conditions of determining AVO detectability. Nicolao et al.(1993)and Martinez(1993)studied the reliability of inversion results at a plane reflector in homogenous elastic media. Joakim(1 998)investigated qualitatively and semi-quantitatively the error of inversion results caused by amplitude attenuation in 1D viscoelastic media. These analyses above were limited to the analyses of P-wave. Also the missing from these studies is the expressions of wavefront spreading and absorption for multi-waves in a layered viscoelastic medium. In this paper, we extend the model given by Nicolao et al.(1993) and Joakim (1998)to viscoelastic media and derive expressions of wavefront spreading and absorption for P-P and P-S converted waves in a layered viscoelastic medium. The imaging matrix at a plane reflector between vls coelastic media can be determined in the frequency domain using lineariaed reflection coefficients through Born approximation. We quantitatively analyze the sensitivity by studying eigenvalues and eigenvectors of the imaging matrix. The results show that two linear combinations of petrophysieal parameters can be determined from the multi-wave AVO inversion in the case of amplitude compensation. Muhi-wave AV0 contains the information of attenuation in the media. However, the sensitivity of multi-wave AV0 inversion to attenuation is small. Especially, if we have difficulty in determining accurately quality factor(Q).the average value of Q within layers is the best approximation when the absorption to stratiform media Is compensated in the attenuation of model.
We consider a layered viscoelastic medium with down-going P-wave as the incident wave. The reflected waves shall be up-going P-wave and mode-converted S-wave. Let α, β, ρ and Q respectively be P-wave velocity, S-wave velocity, density and quality factor of the homogenous layers. The reflection amplitude of the p-wave and S-wave are related to source amplitude factor, As; propagation factor, Ap, and reflection response, Ar. It can be expressed in potential function
|
(1) |
where As is the complex amplitude of the source, which can he both frequency dependent and anisotropic; kx and kz are the horizontal and vertlcal wavenumbers of a plane wave, respectively;
|
(2) |
Here θ, Ψare the P-wave incidence angle and the P_S wave reflection angle, respectively.
Ap onsists of wavefront spreading factor Ar, medium absorption As and transmission coefficient At. Since At can be directly calculated from the Zoeppritz equation, we only discuss the calculation of Ar and Aa in a viscoelastic medium.
Based on the three laws of refraction, reflection and transmission (Hubral and Key, 1980), we can derive wavefront spreading factors in a layered or slanted layered medium for both P-P and P-S waves in terms of the radius of curvature matrix. After some algebra, we obtain
|
(3) |
where αkd, βkd, respectively, are the incident and transmission angles of downgoing P-wave or S-wave at the kth reflector; αku, βku are the incident and transmission angle of upgoing P-wave or S-wave at the kth reflector; RE is the radius of curvature matrix for the P-P and P-S wavefronts. v1u, v1d, oare the velocities ofupgoing P or S wave in the first layer. Using the law of wavefront curvature, we have
|
(4) |
where Rx, is the radius of wavefront in the plane perpendicular to the plane of incidence; Ry, is the radius of wavefront in the plane of incidence. Using the three laws of refraction, reflection and transmission, we obtain the radius of curvature of the P-P and P-S waves through recursion. We have
|
(5) |
|
(6) |
where viu, vid, orespectively are the velocities of upgoing and downgoing waves in the ith layer. siu, sid are the ray lengths of upgoing and downgoing waves in the ith layer.
In horizontal layered medium, there exist the following relationships
|
where hi is thickness of the ith layer. After simplifying equation (3), the multiwave wavefront spreadingfactor in nth layer is
|
(7) |
In general, wavenumber in a viscoelastic medium is complex. Let
|
(8) |
where k' is the real part of multiwave wavenumber. The imaginary part
|
(9) |
where k'x is the real part of kx; k'z is the real part of kz. For convenience, from now on we will denote k'x by kx, k'z by kz, , and k' by k.
Using equations (2)-(5) and the relations among
|
(10) |
Consequently, the absorption coefficients of Pwave and S-wave can be written in a single form in horizontally layered medium. That is
|
(11) |
Equations (7) and (11) are the formulae of am-plitude compensation suitable for AVO inversion inlayered viscoelastic media.
In a homogenous model, the relative perturbations of the model parameters can be written as
|
(12) |
Using Born approximation (Nicolao et al., 1993), we obtain the relationship among reflection coefficient, incident angle (θ) and depth(z). It is
|
(13) |
where R(θ) is deterrnined by the equation given by Aki and Richards (1989), which consists of the reflection coefficients of Pwave and S-wave for multi-wave prospecting, i. e.,
|
(14) |
Fourier transforming equation (13) with respect tovariable z yields
|
(15) |
Let the frequency spectrum of the source be 1. The spreading factor and absorption of multiwave in homogenous background model can be obtained by simplifying equations (7) and (11). Further, using equation (15), we can obtain the following matrix equation for the reflection amplitudes and incident angles,
|
(16) |
where d is the vector of reflection amplitudes defined by equation (1), and
|
G is the sensitivity matrix, which will form the base of the inversion algorithm and
|
(17) |
where
|
(18) |
|
(19) |
|
(20) |
The matrix G in equation (16) is the imagingr matrix from model space to data space. The magnitude of the eigenvalues of G determines the stability of model parameters estimated (Gu, 2000), and reveals whether or not the best solution or the best linear combination of solutions can be obtained. The condition number of matrix G can be determined by singular value decomposition (SVD). This enables us to perform quantitative analysis of the sensitivity of imaging matrix in the case of inversion with amplitude compensation.
Assuming that z is equal to 2 000 m, and
Figure 1a and b respectively shows the condition number and the ratio of the two largest eigenvalues for different values of maximum incidence angle and attenuation. The condition number and the ratio in crease with the increasing of incidence angle. This shows that in the case of amplitude compensation the sensitivity of AVO to elastic parameters at large incidence angles is stronger than that at small incidence angles. That is to say, the elastic parameters are more likely to be determined at large incidence angles. Meanwhile, these figures show that the condition number and the eigenvalue ratio increase slightly with the increasing attenuation. Although multiwave AVO contains information of media attenuation, their sensitivity to attenuation is small. In other words, the effect of Q value on elastic parameters is small, which implies that an average Q value for all layers over the reservoir is probably sufficient when it is difficult to determine accurate Q value for each layer in an attenuation model(Clayton and Stolt, 1981).
Figure 2 and 3 show the isochron for each component of the eigenvectors in the model space corresponding to the largest and the 2nd largest eigenvalues, respectively. From these figures, we observe that the eigenvector for the largest eigenvalue points mainly at Δα/α and Δρ/ρ, corresponding to relative P-wave impedance variation. While the eigenvector for the 2nd largest eigenvalue mainly points at Δβ/β, Δρ/ρ, i. e. relative S-wave impedance variation. Therefore, in the case of amplitude compensation, two linear combinations of elastic parameters can be determined from multi-wave AVO inversion using equation (16). Furthermore, since the gradient direction of eigenvectors as a function attenuation is invariable (Fig. 1), variation in attenuation does not change the combination of elastic parameters.
(1) Two linear combinations of elastic parameters can be determined from multi-wave AVO in the case of amplitude compensation. The reliability of the elastic parameters estimated by AVO inversion is higher in large incidence angle.
(2) Although multi-wave AVO contains the information of media attenuation, the sensitivity of multi-wave AVO inversion to attenuation is small. The average Q value of all layers above reservoir is probably the best approximate value when it is difficult to determine an accurate Q value for each layers in the model.
(3) Absorption in viscoelastic media does not change the combination relation of elastic parameters in AVO inversion.
Boyd, R., Diessel, C. F. K., 1994. The Application of Sequence Stratigraphy to Non-marine Clastics and Coal. In: Posamentier, H. W., Mutti, E., eds., The Second High-Resolution Sequence Stratigraphy Conference Proceedings. Tremp, Spain. 13–20 |
Chen, C. M., Huang, J. K., Chen, J. S., et al., 1981. Evolution of Sedimentary Tectonics of Bohai Rift System and Its Bearing on Hydrocarbon Accumulation. Science in China (Series A), 24: 521–529 |
Chi, Y. L., Sun, H. J., Li, W. L., 1998. Cause of Para-sequence and Their Meaning of Dongying Formation of Nearshore Area in the Eastern Parts of Liaohe Depression. Pertroleum Exploration and Development, 25(3): 29–34 (in Chinese with English Abstract) |
Cross, T. A., 1994. High-Resolution Stratigraphic Correlation from the Perspective of Base-Level Cycles and Sediment Accommodation. In: Dolson, J., ed., Unconformity Related Hydrocarbon Exploration and Accumulation in Clastic and Carbonate Settings. Short Course Notes, Rocky Mountain Association of Geologists. 28–41 |
Dong, Y. X., Zhou, H. M., Xia, W, C., 2003. Relationship between Tertiary Sequence Stratigraphy and Oil Reservoiring in Nanpu Depression. Oil & Gas Geology, 24(1): 39–41, 49 (in Chinese with English Abstract) |
Eschard, R., Lemouzy, P., Bacchiana, C., et al., 1998. Combining Sequence Stratigraphy, Geostatistical Simulations and Production Data for Modeling a Fluvial Reservoir in the Chaunoy Field (Triassic, France). AAPG Bulletin, 82(4): 545–568 |
Galloway, W. E., 1989. Genetic Stratigraphic Sequences in Basin Analysis I: Architecture and Genesis of Flooding-Surface Bounded Depositional Units. AAPG Bulletin, 73: 125–142 |
Liu, K. Y., Liang, T. C. K., Paterson, L., et al., 1998. Computer Simulation of the Influence of Basin Physiography on Condensed Section Deposition and Maximum Flooding. Sediment. Geology, 122: 181–191 doi: 10.1016/S0037-0738(98)00105-5 |
Meng, P., Liu, L., Sun, X. M., et al., 2005. The Application of Micropaleontology-Sequence Stratigraphy in Stratigraphic Division. Acta Micropalaeontologica Sinica, 22(4): 417–424 (in Chinese with English Abstract) |
Mitchum, R. M., Van Wagoner, J. C., 1991. High-Frequency Sequences and Their Stacking Patterns: SequenceStratigraphic Evidence of High-Frequency Eustatic Cycles. Sediment. Geology, 70: 131–160 doi: 10.1016/0037-0738(91)90139-5 |
Nadon, G. C., Simo, J. A., 2000. High-Resolution Sequence Stratigraphic Analysis of the St. Peter Sandstone and Glenwood Formation (Middle Ordovician), Michigan Basin, U.S.A. . AAPG Bulletin, 84(7): 975–996 |
Pinous, O. V., Karogodin, Y. N., Ershov, S. V., et al., 1999. Sequence Stratigraphy, Facies, and Sea Level Change of the Hauterivian Productive Complex, Priobskoe Oil Field (West Siberia). AAPG Bulletin, 83(6): 972–989 |
Posamentier, H. W., Allen, G. P., 1993. Variability of the Sequence Stratigraphic Model: Effects of Local Basin Factors. Sediment. Geology, 86: 91–109 doi: 10.1016/0037-0738(93)90135-R |
Posamentier, H. W., Jervey, M. T., Vail, P. R., 1988. Eustatic Controls on Clastic Deposition I—Conceptual Framework. In: Wilgus, C. K., Hastings, B. S., Kendall, C. G. S. C., et al., eds., Sea-Level Changes: An Integrated Approach. SEPM Special Publication, 42: 109–124 |
Posamentier, H. W., Morris, W. R., 2000. Aspects of the Stratal Architecture of Forced Regressive Deposits. In: Hunt, D., Gawthorpe, R. L., eds., Sedimentary Responses to Forced Regressions. Geological Society (London) Special Publication, 172: 19–46 |
Posamentier, H. W., Vail, P. R., 1988. Eustatic Control on Clastic Deposition Ⅱ—Sequence and Systems Tract Models. In: Wilgus, C. K., Hastings, B. S., Kendall, C. G. S. C., et al., eds., Sea-Level Changes: An Integrated Approach. SEPM Special Publication, 42: 125–154 Shanley, K. W., McCabe, P. J., 1994. Perspectives on the Sequence Stratigraphy of Continental Strata. AAPG Bulletin, 78: 544–568 |
Stephen, F., Simon, K., 1998. Application of High-Resolution Sequence Stratigraphy to Northwest Hutton Field, Northern North Sea: Implications for Management of a Mature Brent Group Field. AAPG Bulletin, 82(7): 1416–1436 |
Tian, Z. Y., Zhang, Q. C., 1996. The Theory about the Hydrocarbon-Bearing Basin in China. Petroleum Industry Press, Beijing. 65–71 (in Chinese) |
Thorne, J. A., Swift, D. J. P., 1991. Sedimentation on Continental Margins, Ⅱ: Application of the Regime Concept. In: Swift, D. J. P., Oertel, G. F., Tillman, R. W., et al., eds., Shelf Sand and Sandstone Bodies. Special Publication of the International Association of Sedimentologists, 14: 33–58 |
Vail, P. R., Mitchum, R. M., Thompson, S., 1977. Seismic Stratigraphy and Global Changes of Sea Level: Part 4, Global Cycles of Relative Changes of Sea Level. In: Payton, C. E., ed., Seismic Stratigraphy—Applications to Hydrocarbon Exploration. AAPG Memoir, 26: 83–97 |
Van Wagoner, J. C., Mitchum, R. M., Campion, K. M., 1990. Siliciclastic Sequence Stratigraphy in Well Logs, Cores and Outcrops: Concept for High-Resolution Correlation of Time and Facies. AAPG Methods in Exploration Series, 7: 55 |
Van Wagoner, J. C., Posamentier, H. W., Mitchum, R. M., et al., 1988. An Overview of the Fundamentals of Sequence Stratigraphy and Key Definitions. In: Wilgus, C. K., Hastings, B. S., Kendall, C. G. S. C., et al., eds., Sea-Level Changes: An Integrated Approach. SEPM Special Publication, 42: 39–45 |
Wang, H. L., Deng, H. W., 2000. Tertiary Sequence Stratigraphy and Major Gas Fields in Bohai Bay Basin. China Offshore Oil and Gas (Geology), 14(2): 100–103, 117 (in Chinese with English Abstract) |
Wang, J., Jiang, Z. X., Chen, S. Y., et al., 2005. Sequence Stratigraphy and Subtle Hydrocarbon Reservoir of Paleogene in Chezhen Sag, Bohai Bay Basin. Journal of Jilin University (Earth Science Edition), 35(2): 163–169 (in Chinese with English Abstract) |
Wei, K. S., Shi, S. L., Li, D. L., et al., 2004. High Accuracy Sequence Stratigraphy of Lower Tertiary in Zhanhua Depression, Bohai Bay Basin. Oil & Gas Geology, 25(5): 513–518 (in Chinese with English Abstract) |
Wright, V. P., Marriott, S. B., 1993. The Sequence Stratigraphy of Fluvial Depositional Systems: The Role of Floodplain Sediment Storage. Sediment. Geology, 86: 203–210 doi: 10.1016/0037-0738(93)90022-W |
Wu, X. L., Zou, B. F., 2002. Sedimentary System and Distribution Features of Early Tertiary in Liaohe Tidal Zone. Special Oil & Gas Reservoirs, 9(3): 13–17 (in Chinese with English Abstract) |
Wu, Y. Y., Liu, Z., Wang, W. H., et al., 1996. Sequence Stratigraphy of Shahejie Formation in Qibei Sag. Acta Sedimentologica Sinica, 14(1): 167–175 (in Chinese with English Abstract) |
Xiao, G. L., Chen, J. W., 2003. Oil and Gas in Neogene of the Bohai Sea Area. Marine Geology Letters, 19(8): 1–6 (in Chinese with English Abstract) |
Yan, H., 1990. The Alluvial Fan, Fan-Delta and Sublacustrine Fan of Paleogene Age within Liaohe Rift, Liaoning Province, China. Sediment. Geology, 68: 75–85 doi: 10.1016/0037-0738(90)90120-I |
Yang, S. S., Sheng, L. Y., 2004. The Characteristic and Property of Physical-Mechanics about Deposit of Sea Beach Area of the Bohai Bay West Coast. Hydrogeology & Engineering Geology, 31(3): 93–95 (in Chinese with English Abstract) |
Yang, Y. T., Xu, T. G., 2004. Hydrocarbon Habitat of the Offshore Bohai Basin, China. Marine and Petroleum Geology, 21(6): 691–708 (in Chinese with English Abstract) doi: 10.1016/j.marpetgeo.2004.03.008 |
Zaitlin, B. A., Warren, M. J., Potocki, D. J., et al., 2002. Depositional Styles in a Low Accommodation Foreland Basin Setting: An Example from the Basal Quartz (Lower Cretaceous), Southern Alberta. Bulletin of Canadian Petroleum Geology, 50: 31–72 doi: 10.2113/50.1.31 |
Zhao, X. Z., Zheng, L. H., Li, T. H., 2004. Structural and Depositional Feature of Nearshore Area in Bohai Bay and Favorable Exploration Zones. Oil Geophysical Prospecting, 39(3): 348–353 (in Chinese with English Abstract) |
Zhu, X. M., Wang, G. W., 2000. Sequence Stratigraphy of Lower Tertiary in Nearshore Area of Liaohe Depression. Petroleum Industry Press, Beijing (in Chinese) |