Pre-Stack Nonlinear Direct Exact Inversion of Fracture Parameters in Deep Shale Reservoirs

A conventional linear pre-stack inversion method under the conventional stationary convolution model is limited by the assumptions of weak formation contrast change and small angle incidence and fails to take into account the amplitude attenuation of seismic wave propagation. Meanwhile, the resolution and precision of oil and gas evaluation and fracture characterization of shale reservoirs under complex geological conditions are low because the compaction and non-connectivity characteristics of deep shale reservoirs are not fully considered. Therefore, porous rock pores are divided into connected pores and disconnected pores. Combined with the effect of compaction on dry rock skeleton, a petrophysical model considering the compaction and pore dysconnectivity of deep shale reservoir is developed. The quantitative relationship between transverse isotropy with a vertical axis of symmetry (VTI) stiffness matrix, rock physical properties, and fracture parameters is established in this model. It provides a more accurate scheme for the original physical modeling of deep shale. This relationship is incorporated into the exact VTI reflection coefficient equation, and a nonstationary convolution operator is derived by using the attenuation theory of seismic wave propagation. A nonstationary pre-stack nonlinear direct inversion method of fracture parameters of shale reservoirs with horizontal fractures is proposed, which Improves the resolution and accuracy of shale reservoir gas bearing and fracture characteristics prediction. It provides a new way to accurately characterize the fracture development and oil-bearing property of shale reservoirs. A model test and field data test verify the effectiveness of this method.