Nested gridding and streamline-based simulation for fast reservoir performance prediction

Detailed reservoir models routinely contain 10⁶–10⁸ grid blocks. These models often cannot be used directly in a reservoir simulation because of the time and memory required for solving the pressure grid on the fine grid. We propose a nested gridding technique that efficiently obtains an approximate solution for the pressure field. The domain is divided into a series of coarse blocks, each containing several fine cells. Effective mobilities are computed for each coarse grid block and the pressure is then found on the coarse scale. The pressure field within each coarse block is computed using flux boundary conditions obtained from the coarse pressure solution. Streamline-based simulation is used to move saturations forward in time. We test the method for a series of example waterflood problems and demonstrate that the method can give accurate estimates of oil production for large 3D models significantly faster than direct simulation using streamlines on the fine grid, making the method overall approximately up to 1,000 times faster than direct conventional simulation.     

EDI electron time-of-flight measurements on Equator-S

We present the first electron time-of-flight measurements obtained with the Electron Drift Instrument (EDI) on Equator-S. These measurements are made possible by amplitude-modulation and coding of the emitted electron beams and correlation with the signal from the returning electrons. The purpose of the time-of-flight measurements is twofold. First, they provide the drift velocity, and thus the electric field, when the distance the electrons drift in a gyro period becomes sufficiently large. Second, they provide the gyro time of the electrons emitted by the instrument, and thus the magnitude of the ambient magnetic field, allowing in-flight calibration of the flux-gate magnetometer with high precision. Results of both applications are discussed.     

The 15 August 2007 Peru Earthquake and Tsunami: Influence of the Source Characteristics on the Tsunami Heights

The tsunami caused by the 2007 Peru earthquake (M_w 8.0) provoked less damage than by the seismic shaking itself (numerous casualties due to the earthquake in the vicinity of Pisco). However, it propagated across the Pacific Ocean and small waves were observed on one tide gauge in Taiohae Bay (Nuku Hiva, Marquesas, French Polynesia). We invert seismological data to recover the rupture pattern in two steps. The first step uses surface waves to find a solution for the moment tensor, and the second step uses body waves to compute the slip distribution in the source area. We find the slip distribution to consist of two main slip patches in the source area. The inversion of surface waves yields a scalar moment of 8.9 10²⁰ Nm, and body-wave inversion gives 1.4 10²¹ Nm. The inversion of tsunami data recorded on a single deep ocean sensor also can be used to compute a fault slip pattern (yielding a scalar moment of 1.1 10²¹ Nm). We then use these different sources to model the tsunami propagation across the Pacific Ocean, especially towards Nuku Hiva. While the source model taken from the body-wave inversion yields computed tsunami waves systematically too low with respect to observations (on the central Pacific Ocean DART buoy as on the Polynesian tide gauge), the source model established from the surface-wave inversion is more efficient to fit the observations, confirming that the tsunami is sensitive to the low frequency component of the source. Finally we also discuss the modeling of the late tsunami arrivals in Taiohae Bay using several friction coefficients for the sea bottom.     

台湾海峡铀的来源及保守性混合

Seawater samples are collected in the spring of 2013 from the Taiwan Strait for the analysis of uranium(U)concentrations and isotopic compositions using MC-ICP-MS, and the geochemical behavior patterns of U in the Taiwan Strait are then investigated. Average concentrations of individual U isotopes are(3.23±0.14) μg/kg for 238 U,(2.34±0.09)×10~(–2) μg/kg for ~(235)U and(2.05±0.07)×10~(–4) μg/kg for 234 U. Correspondingly, the U isotopic compositions are 155±18 for δ234U and 138±2 for 238U:235U. The U concentrations and isotopic ratios in the Taiwan Strait are similar to those of open ocean seawater, suggesting the dominance of the open ocean input to the strait's U pool.However, river input, as suggested by the slightly lower salinity than that of the open ocean, also affected the U concentrations and isotopic compositions in the strait. From a compilation of U concentrations in the Taiwan Strait and adjacent areas, including the Jiulong Estuary and Zhujiang Estuary, the Xiamen Bay and the northern South China Sea, a strong and significant relationship between U concentration and salinity [U:S; U=(0.093 4±0.002 4)S+(0.092 0±0.061 5)] is revealed, suggesting conservative mixing of U in the Taiwan Strait. To better understand the U geochemistry in the Taiwan Strait, a multiple endmembers mixing model is applied to estimate the contributions of potential sources. The open ocean seawater contributed 69%–95% of U in the Taiwan Strait, with river water approximately 2%, and dust deposition only around 0.13%. Therefore, the model results supported the open ocean input source and the conservative mixing behavior of U derived from the observation of U concentrations and isotopic ratios and U:S ratios. The sediment interstitial water may be an important source of U to the Taiwan Strait with a possible contribution of 3%–29%, consistent with previous investigations based on radium isotopes.However, further investigations are warranted to examine the U concentration in the sediment interstitial water and its input to the overlying seawater in the Taiwan Strait.     

Microseismic feasibility studies – assessing the probability of success of monitoring projects

In this paper we present a workflow for microseismic feasibility studies that allows a thorough assessment of the probability of success of the monitoring project to be made. The workflow includes the following elements: assessment of the value of information to establish the business case; identification of hardware options to investigate deployment options; survey of analogue cases to confirm the do‐ability of the project; modelling of event location errors and detectability to establish the most favourable array geometry; generation of full waveform synthetics to anticipate undesired seismic features. A study comprising all these elements provides the reference frame for discussions with service companies, limiting misunderstandings and avoiding missed opportunities.