Estimation of fracture parameters from P-wave AVOA data based on equivalent media theory

Xue Jiao1,2 Gu Hanming1,2 Cai Chengguo1,2 Li Zongjie3 Zhu Ding3

(1.Institute of Geophysics and Geomatics, China University of Geosciences (Wuhan), Wuhan, Hubei, China 430074)
(2.Hubei Subsurface Multi-scale Imaging Key Laboratory, Wuhan, Hubei, China 430074)
(3.Research Institute of Exploration and Development, Northwest Oilfield Branch Co., SINOPEC, Urumqi, Xinjiang, China 830011)

【Abstract】The contrasts of the fracture weakness may be positive or negative while the fracture weakness is nonnegative. Layer stripping is normally used to invert the fracture parameters of individual layers in sections. Otherwise, the upper layer of each interface is assumed to be isotropic to simplify inversion algorithms. Large errors in the final estimations of the fracture parameters may be introduced by layer stripping or inversion methods assuming that upper layer of each interface is isotropic. Thus, we propose in this paper a method for estimation of fracture parameters from P-wave AVOA data based on equivalent media theory. We express the contrasts of the fracture weakness as multiplication of the first-order difference matrix with the fracture weakness itself, and express the reflected amplitude as a function of the fracture weakness, which can be used to estimate normal and tangential weakness. A priori information about the model can be added to the objective function as penalty term according to Bayesian theory. The AVOA inversion results for a synthetic model show that the fracture weakness can be reliably estimated. Fracture weakness and fracture density are estimated from real seismic data in Makit area in the Tarim Basin. The AVOA inversion results are consistent with well and logging information, which means they can be used to determine fractured reservoir locations.

【Keywords】 equivalent media theory for fractures; fracture weakness; fracture density; P-wave AVOA;


【Funds】 National Science and Technology Major Project (2011ZX05000-023-004-001) Special Fund for Research in the Public Interest by Ministry of Land and Resources (201211082-03) Open Fund Project of Hubei Subsurface Multi-scale Imaging Key Laboratory (SMIL-2016-02)

Download this article


    [1]Ruger A, Tsvankin I. Using AVO for fracture detection: Analytic basis and practical solutions. The Leading Edge, 1997, 16 (10): 1429–1434.

    [2]Tsvankin I. Reflection moveout and parameter estimation for horizontal transverse isotropy. Geophysics, 1997, 62 (2): 614–629.

    [3]Qu S L, Ji Y X, Wang X et al. Seismic method for using full-azimuth P wave attribution to detect fracture. OGP, 2001, 36 (4): 390–397 (in Chinese).

    [4]Zhao S W, He Z H, Xiong X J et al. Application of azimuthal P-wave fracture detection technique in fracture reservoir in Northeast Sichuan, China. OGP, 2010, 45 (2): 260–264 (in Chinese).

    [5]Wang H Q, Yang W Y, Xie C H et al. Azimuthal anisotropy analysis of different seismic attributes and fracture prediction. OGP, 2014, 49 (5): 925–931 (in Chinese).

    [6]Zhen J J, Liu Y. Review over physical model of fractured rock medium. Progress in Geophysics, 2011, 26 (5): 1708–1716 (in Chinese).

    [7] Gui Z X. Seismic wave field feature theory for fractured media and its application. Wuhan, Hubei Province, Hubei Science & Technology Press, 2005 (in Chinese).

    [8]Hsu C J, Schoenberg M. Elastic waves through a simulated fractured medium. Geophysics, 1993, 58 (7): 964–977.

    [9]Schoenberg M, Sayers C M. Seismic anisotropy of fractured rock. Geophysics, 1995, 60 (1): 204–211.

    [10]Shaw R K, Sen M K. Use of AVOA data to estimate fluid indicator in a vertically fractured medium. Geophysics, 2006, 71 (3): 15–24.

    [11]Alhussain M, Sen M K. Quantitative estimation of fracture parameters after removing anisotropic overburden effect. SEG Technical Program Expanded Abstracts, 2013, 32:2443–2447.

    [12]Chen H Z, Yin X Y, Gao C G et al. AVAZ inversion for fluid factor based on fracture anisotropic rock physics theory. Chinese Journal of Geophysics, 2014, 57 (3): 968–978 (in Chinese).

    [13]Chen H, Zhang G, Chen J et al. Fracture filling fluids identification using azimuthally elastic impedance based on rock physics. Journal of Applied Geophysics, 2014, 110 (11): 98–105.

    [14]Ruger A. Variation of P-wave reflectivity with offset and azimuth in anisotropic media. Geophysics, 1998, 63 (3): 935–947.

    [15]Bakulin A, Grechka V, Tsvankin I. Estimation of fracture parameters from reflection seismic data. Part 1:HTI model due to a single fracture set. Geophysics, 2000, 65 (6): 1788–1802.

    [16]Far M E, Hardage B, Wagner D. Fracture parameter inversion for Marcellus Shale. Geophysics, 2014, 79 (3): C55-C63.

This Article


CN: 13-1095/TE

Vol 51, No. 06, Pages 1171-1179+1051-1052

December 2016


Article Outline


  • 1 Introduction
  • 2 Inversion theory using fracture parameters
  • 3 Tentative calculation of fracture parameter inversion
  • 4 Conclusions
  • References