Assessment of gas hydrate reservoir from inverted seismic impedance and porosity in the northern Hikurangi margin, New Zealand

Abstract

Hikurangi margin, New Zealand, the shallowest subduction zone on Earth is the unique place to study the slow slip creep-like deformation and seafloor failure, and their probable link to the presence of gas hydrate, cold seeps and fluid migration. International Ocean Discovery Program expedition 372 discovered gas hydrate within thin intervals (ranging from ~5 to 25 m) of the deformed sediments in the Tuaheni landslide complex area of the northern Hikurangi margin. Seismic profiles show typical BSRs (bottom simulating reflectors) cutting across the dipping strata throughout the sections. To characterize the reservoir, here, we invert both acoustic impedance (product of velocity and density) and porosity along two perpendicular seismic profiles crossing the well using a model-based post-stack seismic inversion technique. Results show the layer structure of alternate high-impedance with low-porosity and low-impedance with high-porosity sediments. Interestingly, each layer is consisting of thin and relatively low impedance intra-layers, which clearly indicates the possible pathways for fluid and gas migration, and one of the reasons for seafloor collapse in this region. Very low impedance observed below the BSR is due to the gas-charged sediments at the dipping strata, and is the probable source of free-gas migrating upward. Next, we apply rock physics theory to know the distribution of gas hydrate at Site U1517 as well as along two perpendicular seismic profiles. Rock physics modeling using sonic velocity at well location shows that gas hydrate is distributed mainly within the depth intervals of 5–50 m, 90–115 m and 130–155 m with an average saturation of about 10–15% and the maximum concentration of about 35% of the pore space at 100 m depth. Prediction of gas hydrate saturation from the resistivity method using neutron porosity and NMR (nuclear magnetic resonance) logs in Archie's empirical formula shows less saturation compared to that from the sonic log. Whereas, estimated gas hydrate saturation from the chloride concentrations ranges from 0 to 45% of the pore space. The spatial distributions of gas hydrate along seismic profiles are predicted by rock physics theory using the inverted impedance and porosity. Saturation of gas hydrate along the seismic profiles varies from 0 to 40% with an average of ~10% of the pore space. The correlation between inverted impedance and porosity indicates that the high impedance layers are due to low porosity consolidated sediment without significant gas hydrate concentration. Gas hydrate is distributed in relatively high porosity and low impedance layers along the seismic profiles. Our result shows that the amount of gas hydrate concentration is low in this region and sediments are highly deformed/disturbed in the presence of many fluid/gas migration pathways.