Three‐dimensional finite‐difference (FD) modeling of seismic scattering from free surface topography has been pursued. We have developed exact 3-D free surface topography boundary conditions for the particle velocities. A velocity‐stress formulation of the full elastic wave equations together with the boundary conditions has been numerically modeled by an eighth‐order FD method on a staggered grid. We give a numerical stability criterion for combining the boundary conditions with curved‐grid wave equations, where a curved grid represents the physical medium with topography. Implementation of this stability criterion stops instabilities from arising in areas of steep and rough topographies. We have simulated scattering from teleseismic P-waves using a plane, vertically incident wavefront and real topography from a 40 × 40 km area centered at the NORESS array of seismic receiver stations in southeastern Norway. Synthetic snapshots and seismograms of the wavefield show clear conversion from P-waves to Rg (short period fundamental mode Rayleigh) waves in an area of rough topography approximately 10 km east of NORESS. This result is consistent with numerous observations.

No abstract

Summary A method based on the ensemble Kalman filter (EnKF) for continuous model updating with respect to the combination of production data and 4D seismic data is presented. When the seismic data are given as a difference between two surveys, a combination of the ensemble Kalman filter and the ensemble Kalman smoother has to be applied. Also, special care has to be taken because of the large amount of data assimilated. Still, the method is completely recursive, with little additional cost compared to the traditional EnKF. The model system consists of a commercial reservoir simulator coupled with a rock physics and seismic modeling software. Both static variables (porosity, permeability, and rock physic parameters) and dynamic variables (saturations and pressures) may be updated continuously with time based on the information contained in the assimilated measurements. The method is applied to a synthetic model and a real field case from the North Sea. In both cases, the 4D seismic data are different variations of inverted seismic. For the synthetic case, it is shown that the introduction of seismic data gives a much better estimate of reservoir permeability. For the field case, the introduction of seismic data gives a very different permeability field than using only production data, while retaining the production match. Introduction The Kalman filter was originally developed to update the states of linear systems (Kalman 1960). For a presentation of this method in a probabilistic, linear least-squares setting, see Tarantola (2005). However, this method is not suitable for nonlinear models, and the ensemble Kalman filter (EnKF) method was introduced in 1994 by Geir Evensen for updating nonlinear ocean models (Evensen 1994). The method may also be applied to a combined state and parameter estimation problem (Evensen 2006; Lorentzen 2001; Anderson 1998). Several recent investigations have shown the potential of the EnKF for continuous updating of reservoir simulation models, as an alternative to traditional history matching (Nævdal et al. 2002a, b; Nævdal et al. 2005; Gu and Oliver 2004; Gao and Reynolds 2005; Wen and Chen 2005). The EnKF method is a Monte Carlo type sequential Bayesian inversion, and provides an approximate solution to the combined parameter and state-estimation problem. The result is an ensemble of solutions approximating the posterior probability density function for the model input parameters (e.g., permeability and porosity), state variables (pressures and saturations), and other output data (e.g., well production history) conditioned to measured, dynamic data. Conditioning reservoir simulation models to seismic data is a difficult task (Gosselin et al. 2003). In this paper, we show how the ensemble Kalman filter method can be used to update a combined reservoir simulation/seismic model using the combination of production data and inverted 4D seismic data. There are special challenges involved in the assimilation of the large amount of data available with 4D seismic, and the present work is based on the work presented by Evensen (2006, 2004) and Evensen and van Leeuwen (2000). In the following, the combined state and parameter estimation problem is described in a Bayesian framework, and it is shown how this problem is solved using the EnKF method, with emphasis on the application to 4D seismic data. When the seismic data are given as a difference between two surveys, a combination of the ensemble Kalman filter and the ensemble Kalman smoother has to be applied. Special challenges involved when the amount of data is very large are discussed. The validity of the method is examined using a synthetic model, and finally, a real case from the North Sea is presented.

Summary A method based on the ensemble Kalman filter (EnKF) for continuous model updating with respect to the combination of production data and 4D seismic data is presented. When the seismic data are given as a difference between two surveys, a combination of the ensemble Kalman filter and the ensemble Kalman smoother has to be applied. Also, special care has to be taken because of the large amount of data assimilated. Still, the method is completely recursive, with little additional cost compared to the traditional EnKF. The model system consists of a commercial reservoir simulator coupled to a rock physics and seismic modelling software. Both static variables (porosity, permeability, rock physic parameters, etc.) and dynamic variables (saturations and pressures) may be updated continuously with time based on the information contained in the assimilated measurements. The method is applied to a synthetic model and a real field case from the North Sea. In both cases, the 4D seismic data are different variations of inverted seismic. For the synthetic case, it is shown that the introduction of seismic data gives a much better estimate of reservoir permeability. For the field case, the introduction of seismic data gives a very different permeability field than using only production data, while retaining the production match. Introduction The Kalman filter was originally developed to update the states of linear systems.[1] For a presentation of this method in a probabilistic, linear least-squares setting, see e.g., Tarantola.[2] However, this method is not suitable for non-linear models, and the ensemble Kalman filter (EnKF) method was introduced in 1994 by Geir Evensen for updating non-linear ocean models.[3] It may also be applied to a combined state and parameter estimation problem.[4] Several recent investigations have shown the potential of the EnKF for continuous updating of reservoir simulation models, as an alternative to traditional history-matching.[5–10] The EnKF method is a Monte Carlo type sequential Bayesian inversion, and provides an approximate solution to the combined parameter and state estimation problem. The result is an ensemble of solutions approximating the posterior probability density function for the model input parameters (e.g., permeability and porosity), state variables (pressures and saturations), and other output data (e.g., well production history) conditioned to measured, dynamic data. Conditioning reservoir simulation models to seismic data is a difficult task.[11] In this paper we show how the ensemble Kalman filter method can be used to update a combined reservoir simulation/seismic model using the combination of production data and inverted 4D seismic data. Special challenges are involved in the assimilation of the large amount of data available with 4D seismic, and the present work is based on the work presented by Evensen,[4,12] and Evensen and van Leeuwen.[13] In the following, the combined state and parameter estimation problem is described in a Bayesian framework, and it is shown how this problem is solved using the EnKF method, with emphasis on the application to 4D seismic data. When the seismic data are given as a difference between two surveys, a combination of the ensemble Kalman filter and the ensemble Kalman smoother has to be applied. Special challenges involved when the amount of data is very large are discussed. The validity of the method is examined using a synthetic model, and finally a real case from the North Sea is presented.

We have pursued two‐dimensional (2D) finite‐difference (FD) modelling of seismic scattering from free‐surface topography. Exact free‐surface boundary conditions for the particle velocities have been derived for arbitrary 2D topographies. The boundary conditions are combined with a velocity–stress formulation of the full viscoelastic wave equations. A curved grid represents the physical medium and its upper boundary represents the free‐surface topography. The wave equations are numerically discretized by an eighth‐order FD method on a staggered grid in space, and a leap‐frog technique and the Crank–Nicholson method in time. In order to demonstrate the capabilities of the surface topography modelling technique, we simulate incident point sources with a sinusoidal topography in seismic media of increasing complexities. We present results using parameters typical of exploration surveys with topography and heterogeneous media. Topography on homogeneous media is shown to generate significant scattering. We show additional effects of layering in the medium, with and without randomization, using a von Kármán realization of apparent anisotropy. Synthetic snapshots and seismograms indicate that prominent surface topography can cause back‐scattering, wave conversions and complex wave patterns which are usually discussed in terms of inter‐crust heterogeneities.