Premium
Lithology and hydrocarbon mapping from multicomponent seismic data
Author(s) -
Özdemir Hüseyin,
Flanagan Kevin,
Tyler Emma
Publication year - 2010
Publication title -
geophysical prospecting
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.735
H-Index - 79
eISSN - 1365-2478
pISSN - 0016-8025
DOI - 10.1111/j.1365-2478.2009.00821.x
Subject(s) - amplitude versus offset , amplitude , geology , offset (computer science) , lithology , shear waves , regional geology , prestack , poisson's ratio , shear modulus , shear (geology) , gemology , poisson distribution , economic geology , petrophysics , frequency band , reservoir modeling , mineralogy , hydrogeology , seismology , engineering geology , optics , petrology , physics , geotechnical engineering , mathematics , computer network , metamorphic petrology , bandwidth (computing) , computer science , tectonics , telmatology , thermodynamics , programming language , statistics , volcanism , porosity
Elastic rock properties can be estimated from prestack seismic data using amplitude variation with offset analysis. P‐wave, S‐wave and density ‘reflectivities’, or contrasts, can be inverted from angle‐band stacks. The ‘reflectivities’ are then inverted to absolute acoustic impedance, shear impedance and density. These rock properties can be used to map reservoir parameters through all stages of field development and production. When P‐wave contrast is small, or gas clouds obscure reservoir zones, multicomponent ocean‐bottom recording of converted‐waves (P to S or Ps) data provides reliable mapping of reservoir boundaries. Angle‐band stacks of multicomponent P‐wave (Pz) and Ps data can also be inverted jointly. In this paper Aki‐Richards equations are used without simplifications to invert angle‐band stacks to ‘reflectivities’. This enables the use of reflection seismic data beyond 30° of incident angles compared to the conventional amplitude variation with offset analysis. It, in turn, provides better shear impedance and density estimates. An important input to amplitude variation with offset analysis is the V s / V p ratio. Conventional methods use a constant or a time‐varying V s / V p model. Here, a time‐ and space‐varying model is used during the computation of the ‘reflectivities’. The V s / V p model is generated using well log data and picked horizons. For multicomponent data applications, the latter model can also be generated from processing V s / V p models and available well data. Reservoir rock properties such as λρ, μρ, Poisson's ratio and bulk modulus can be computed from acoustic impedance, shear impedance and density for pore fill and lithology identification. λ and μ are the Lamé constants and ρ is density. These estimations can also be used for a more efficient log property mapping. V p / V s ratio or Poisson's ratio, λρ and weighted stacks, such as the one computed from λρ and λ/μ, are good gas/oil and oil/water contact indicators, i.e., pore fill indicators, while μρ mainly indicates lithology. μρ is also affected by pressure changes. Results from a multicomponent data set are used to illustrate mapping of gas, oil and water saturation and lithology in a Tertiary sand/shale setting. Whilst initial log crossplot analysis suggested that pore fill discrimination may be possible, the inversion was not successful in revealing fluid effects. However, rock properties computed from acoustic impedance, shear impedance and density estimates provided good lithology indicators; pore fill identification was less successful. Neural network analysis using computed rock properties provided good indication of sand/shale distribution away from the existing wells and complemented the results depicted from individual rock property inversions.