Simultaneous inversion application for characterizing Hamra Quartzite tight sand reservoir: a case study from Hassi Messaoud (Algeria)

摘要

Channel sand acts as a stratigraphic trap for hydrocarbon accumulation in many parts of the world. Delineation of this type of reservoir is crucial as channel sand may be scarce, and inaccurate location of the drilling wells could lose a huge currency. The Hassi Messaoud (HMD) field was subjected to multi-phase tectonic events, where deep-seated structures were rejuvenated leading to intensive fault complexity. The main effective tectonic events upon the studied area are the Hercynian compression and deep erosion till the OrdovicianHamra Quartzite (HQZ) oil reservoir, followed by active Triassic rifting and filling the deeply eroded areas or the graben areas by eruptive volcanic rocks at Triassic time. Hercynian erosion and volcanic rocks distribution introduce a big uncertainty to the reservoir structural model. Amplitude versus offset (AVO) method is used as a helpful tool to differentiate channel sand from surrounding formations. Several attributes (P-impedance, S-impedance, longitudinal velocity Vp, shear velocity Vs and density.) are estimated from pre-stack seismic inversion. They have different sensitivity to the reservoir properties. Derived attributes such as Lame parameters, incompressibility x density (lambda rho) and rigidity x density (mu rho) can provide key lithology and fluid indicators (Goodway et al. 1997, Goodway CSEG Rec 26(6): 39-60 2001). Petrophysically relating AVO attributes both to lambda rho and mu rho and to each other in Lambda-Mu-Rho (LMR) cross-plot space can be a good tool for AVO interpretation (Rutherford and Williams Geophysics 54: 680-688 1989 and Castagna and Swan Lead Edge 16(4): 337-342 1997). After proper data conditioning, simultaneous inversion of pre-stack angle gathers is performed to get acoustic wave impedance (P-impedance), elastic wave impedances (S-impedance) and density rho, then to calculate lambda rho and mu rho volumes. In the studied area, lambda rho and mu rho are used as a very important key to separate reservoir sands. The lambda rho and mu rho curves are generated at each well location. Cross plots showed a fair separation of sand in the formation, i.e. higher mu rho and lower lambda rho can detect sand. The output lambda rho and mu rho volumes after simultaneous inversion follow the distribution of the sand which is consistent with the wells penetrating the target reservoir. This finding on the extension of the sand reservoir in terms of lambda rho and mu rho 3D cross-plot zonations are used for lithology discrimination. In this study, well logs were used to constrain lithology and to control the zonation filters by reducing the limits ambiguity. Other types of advanced attributes are calculated and tested. The obtained (mu rho-lambda rho) volume acts as a good indicator for the sand distribution. It was finally used as sand presence index in the area. Also mu rho. has shown a good linear relationship with porosity. To note that the porosity volume is created based on the linear relationship with mu rho. A product of derived porosity and the sand presence index (mu rho-lambda rho) provides a good tool for reservoir characterization and lead to reservoir management, future planning of the field, and setting location for new wells.